SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The files are final forward model simulations for four geophysical array configurations. Protocol files for each of the geophysical array and the forward model code (ResIPy) is included.
Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modeling and real-time control applications. More accurate simulations with fast electron diffusion are even slower, requiring multiple hours of run time with parallel processing. The machine learning models use a database of 16,000+ GEN-RAY/CQL3D simulations for training, validation, and testing. Latin hypercube sampling methods implemented in πScope ensure that the database covers the range of 9 input parameters (n e0 , T e0 , I p , B t , R 0 , n ∥︀ , Z e f f , V loop , P LHCD ) with sufficient density in all regions of parameter space. The surrogate models reduce the computation time from minutes-hours to ms with high accuracy across the input parameter space. Data-driven surrogate models also allow for solving inverse and “lateral” problems. A surrogate model for the inverse problem maps from a desired current drive or power deposition profile to a set of input parameters that would result in such a profile, while a surrogate model for the lateral problem maps from a measured experimental quantity such as hard x-ray emission to a current drive or power deposition profile. In conclusion, the πScope database creation workflow is flexible and applicable to other RF simulation codes such as TORIC.
Like most magnetic confined fusion experiments, the ST40 tokamak started off with a small subset of diagnostics and gradually increased the diagnostic set to include more complex and comprehensive systems. To make the most of each operational phase, forward models of various diagnostics are used and developed to aid design, provide consistency-checks during commissioning, test analysis methods, and build workflows to constrain high-level parameters to inform interpretation, theory and modelling. For new models and new analysis workflows, minimum-viable-products are released early, and their complexity is increased in a step-wise manner, facilitating the support of all programme phases on multiple parallel applications, while enabling learning opportunities and feedback loops. In this contribution we review the philosophy, scope and architecture of the framework under development. We discuss the details of some forward models, with examples on how they are used to aid diagnostic design, to investigate analysis methodologies through synthetic data, and how they are embedded in experimental analysis workflows. We compare previously published experimental results with new, more advanced analysis workflows employing more recent, detailed models and new diagnostic data, providing confirmation of the published material from the 2021–22 experimental campaign.
The development of nuclear reaction models for the production of evaluated nuclear data has traditionally been performed by comparing measured cross sections with predictions from reaction model codes whose physical input parameters are adjusted to obtain the best agreement between measured and modeled results. To more directly probe reaction model inputs, this work introduces a forward modeling approach to experimental reaction cross-section determination, where the most important physical input parameters to reaction model calculations are obtained via 𝜒 2 minimization between measured and calculated observables. This was demonstrated using data collected by the Gamma Energy Neutron Energy Spectrometer for Inelastic Scattering (GENESIS) at the 88-inch cyclotron at Lawrence Berkeley National Laboratory, a detection array consisting of organic liquid scintillators and high-purity germanium (HPGe) detectors. Using a broad-spectrum neutron beam and a 99.98%-enriched 56 Fe target, GENESIS was used to perform a simultaneous measurement of 56 Fe 𝛾-ray production cross sections and secondary neutron energy and angle distributions. The results of the forward modeling approach to the determination of energy-differential 𝛾-ray production cross sections for the yrast 4 + → 2 + and 6 + → 4 + transitions, as well as eight other off-yrast transitions, were compared against those obtained using conventional techniques, and the results are in good agreement. In addition to discrete 𝛾-ray yield total scattered neutron energy-angular distributions as a function of incident neutron energy were also obtained using forward modeling and found to agree with evaluated data, with the exception of elastic scattering at small angles. The fitted reaction model parameters obtained through forward modeling were also used to calculate the cross section for the unobserved (𝑛, 2𝑛) reaction; excellent agreement with the current evaluation was obtained, providing a validation of the predictive capabilities of the forward model approach. This work bridges the gap between nuclear data experiment and evaluation by providing a new means for extracting inelastic neutron-scattering cross sections and neutron-induced 𝛾-ray production data while directly probing reaction model physics.
The Jovian Auroral Distributions Experiment Ion sensor (JADE-I) on NASA’s Juno mission provides in-situ measurements of ions from 0.1 to 46.2 keV/q inside Jupiter’s magnetosphere. JADE-I is used to study the plasma with two types of datasets from the same measurement: Time-of-flight (TOF) and SPECIES. The TOF dataset provides mass-per-charge measurements with a range of 1–64 amu/q but oversamples particles over 6π steradian viewing per spacecraft spin and has little directional information. On the other hand, the SPECIES dataset can provide a good measurement of the flow direction but does not provide mass-per-charge information due to the telemetry limit. In this study, we developed a 2-step forward modeling method that combines the advantages and avoids the disadvantages of TOF and SPECIES data to derive the 3-D properties of heavy ions. Assuming that the ion velocity distribution can be described with the kappa distribution, we first perform the forward model fit of the TOF data to calculate the relative abundance of heavy ion species. Then we fix the relative abundance and perform the second forward model fit on the SPECIES data. Here, using this method, we obtain the densities of different heavy ions, the shared temperature and kappa value, and the 3-D flow velocity vector. Some data examples of the equatorial plasma disk before Perijove 24 are included to demonstrate the method. Plasma properties can then be mapped to explore spatial and temporal variabilities in Jupiter’s magnetosphere.
To study nanostructures on substrates, surface-sensitive reflection-geometry scattering techniques such as grazing incident small angle X-ray scattering are commonly used to yield an averaged statistical structural information of the surface sample. Grazing incidence geometry can probe the absolute three-dimensional structural morphology of the sample if a highly coherent beam is used. Coherent surface scattering imaging (CSSI) is a powerful yet non-invasive technique similar to coherent X-ray diffractive imaging (CDI) but performed at small angles and grazing-incidence reflection geometry. A challenge with CSSI is that conventional CDI reconstruction techniques cannot be directly applied to CSSI because the Fourier-transform-based forward models cannot reproduce the dynamical scattering phenomenon near the critical angle of total external reflection of the substrate-supported samples. To overcome this challenge, we have developed a multislice forward model which can successfully simulate the dynamical or multi-beam scattering generated from surface structures and the underlying substrate. The forward model is also demonstrated to be able to reconstruct an elongated 3D pattern from a single shot scattering image in the CSSI geometry through fast-performing CUDA-assisted PyTorch optimization with automatic differentiation.
The Nancy Grace Roman Space Telescope (Roman) is an observatory for both wide-field observations and coronagraphy that is scheduled for launch in the mid-2020s. Part of the planned survey is a deep, cadenced field or fields that enable cosmological measurements with type Ia supernovae (SNe Ia). With a pixel scale of 0$_.^"$11, the Wide Field Instrument will be undersampled, presenting a difficulty for precisely subtracting the galaxy light underneath the SNe. We use simulated data to validate the ability of a forward-model code (such codes are frequently also called “scene-modeling” codes) to perform precision supernova photometry for the Roman SN survey. Our simulation includes over 760,000 image cutouts around SNe Ia or host galaxies (~10% of a full-scale survey). To have a realistic 2D distribution of underlying galaxy light, we use the VELA simulated high-resolution images of galaxies. We run each set of cutouts through our forward-modeling code which automatically measures time-dependent SN fluxes. Given our assumed inputs of a perfect model of the instrument point-spread functions and calibration, we find biases at the millimagnitude level from this method in four red filters (Y106, J129, H158, and F184), easily meeting the 0.5% Roman inter-filter calibration requirement for a cutting-edge measurement of cosmological parameters using SNe Ia. Simulated data in the bluer Z087 filter shows larger ~ 2–3 mmg biases, also meeting this requirement, but with more room for improvement. Our forward-model code has been released on Zenodo.
ABSTRACT Among the properties shaping the light of a galaxy, the star formation history (SFH) is one of the most challenging to model due to the variety of correlated physical processes regulating star formation. In this work, we leverage the stellar population synthesis model fsps, together with SFHs predicted by the hydrodynamical simulation IllustrisTNG and the empirical model universemachine, to study the impact of star formation variability on galaxy colours. We start by introducing a model-independent metric to quantify the burstiness of a galaxy formation model, and we use this metric to demonstrate that universemachine predicts SFHs with more burstiness relative to IllustrisTNG. Using this metric and principal component analysis, we construct families of SFH models with adjustable variability, and we show that the precision of broad-band optical and near-infrared colours degrades as the level of unresolved short-term variability increases. We use the same technique to demonstrate that variability in metallicity and dust attenuation presents a practically negligible impact on colours relative to star formation variability. We additionally provide a model-independent fitting function capturing how the level of unresolved star formation variability translates into imprecision in predictions for galaxy colours; our fitting function can be used to determine the minimal SFH model that reproduces colours with some target precision. Finally, we show that modelling the colours of individual galaxies with per cent-level precision demands resorting to complex SFH models, while producing precise colours for galaxy populations can be achieved using models with just a few degrees of freedom.
We present cosmological constraints from a simulation-based inference (SBI) analysis of galaxy clustering from the SimBIG forward modeling framework. SimBIG leverages the predictive power of high-fidelity simulations and provides an inference framework that can extract cosmological information on small nonlinear scales. In this work, we apply SimBIG to the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy sample and analyze the power spectrum, P ℓ (k), to k max = 0.5 h/Mpc. We construct 20,000 simulated galaxy samples using our forward model, which is based on 2,000 high-resolution Quijote N -body simulations and includes detailed survey realism for a more complete treatment of observational systematics. We then conduct SBI by training normalizing flows using the simulated samples and infer the posterior distribution of ΛCDM cosmological parameters: Ω m , Ω b , h, n s , σ 8 . We derive significant constraints on Ω m and σ 8 , which are consistent with previous works. Our constraint on σ 8 is 27% more precise than standard P ℓ analyses because we exploit additional cosmological information on nonlinear scales beyond the limit of current analytic models, k > 0.25 h/Mpc. This improvement is equivalent to the statistical gain expected from a standard P ℓ analysis of galaxy sample ~ 60% larger than CMASS. While we focus on P ℓ in this work for validation and comparison to the literature, SimBIG provides a framework for analyzing galaxy clustering using any summary statistic. We expect further improvements on cosmological constraints from subsequent SimBIG analyses of summary statistics beyond P ℓ .
ABSTRACT Novel summary statistics beyond the standard 2-point correlation function (2PCF) are necessary to capture the full astrophysical and cosmological information from the small-scale (r < 30h−1Mpc) galaxy clustering. However, the analysis of beyond-2PCF statistics on small scales is challenging because we lack the appropriate treatment of observational systematics for arbitrary summary statistics of the galaxy field. In this paper, we develop a full forward modelling pipeline for a wide range of summary statistics using the large high-fidelity AbacusSummit light cones that account for many systematic effects as well as remain flexible and computationally efficient to enable posterior sampling. We apply our forward model approach to a fully realistic mock galaxy catalog and demonstrate that we can recover unbiased constraints on the underlying galaxy–halo connection model using two separate summary statistics: the standard 2PCF and the novel k-th nearest neighbour (kNN) statistics, which are sensitive to correlation functions of all orders. We will demonstrate its strong constraining power on extended galaxy–halo connection models and cosmology in follow up papers. We expect this to become a powerful approach when applying to upcoming surveys such as DESI where we can leverage a multitude of summary statistics across a wide redshift range to maximally extract information from the non-linear scales.
We forward-model mass-weighted stellar ages (MWAs) and quiescent fractions (f Q ) in projected phase space (PPS), using data from the Sloan Digital Sky Survey, to jointly constrain an infall quenching model for galaxies in log (M vir /M ⊙ ) > 14 galaxy clusters at z ∼ 0. We find the average deviation in MWA from the MWA–M⋆ relation depends on position in PPS, with a maximum difference between the inner cluster and infalling interloper galaxies of ∼1 Gyr. Our model employs infall information from N-body simulations and stochastic star-formation histories from the universemachine model. We find total quenching times of t Q = 3.7 ± 0.4 Gyr and t Q = 4.0 ± 0.2 Gyr after first pericentre, for 9 < log (M ⋆ /M ⊙ ) < 10 and 10 < log (M ⋆ /M ⊙ ) < 10.5 galaxies, respectively. By using MWAs, we break the degeneracy in time of quenching onset and time-scale of star formation rate (SFR) decline. We find that time of quenching onset relative to pericentre is $t_{\mathrm{delay}}=3.5^{+0.6}_{-0.9}$ Gyr and $t_{\mathrm{delay}}=-0.3^{+0.8}_{-1.0}$ Gyr for 9 < log (M ⋆ /M ⊙ ) < 10 and 10 < log (M ⋆ /M ⊙ ) < 10.5 galaxies, respectively, and exponential SFR suppression time-scales are τenv ≤ 1.0 Gyr for 9 < log (M ⋆ /M ⊙ ) < 10 galaxies and τenv ∼ 2.3 Gyr for 10 < log (M ⋆ /M ⊙ ) < 10.5 galaxies. Stochastic star formation histories remove the need for rapid infall quenching to maintain the bimodality in the SFR of cluster galaxies; the depth of the green valley prefers quenching onsets close to first pericentre and a longer quenching envelope, in slight tension with the MWA-driven results. Taken together these results suggest that quenching begins close to, or just after pericentre, but the time-scale for quenching to be fully complete is much longer and therefore ram-pressure stripping is not complete on first pericentric passage.
We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). Imaging systematics refers to the false fluctuation of galaxy densities due to varying observing conditions and astrophysical foregrounds corresponding to the imaging surveys from which DESI LRG target galaxies are selected. We update the Obiwan pipeline, which we previously developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition allows simulating the DESI LRGs sample, which utilizes WISE data in the target selection. Deep DESI imaging data combined with a method to account for biases in their shapes is used to define a truth sample of potential LRG targets. We inject these data evenly throughout the DESI Legacy Imaging Survey footprint at declinations between -30 and 32.375 degrees. We simulate a total of 15 million galaxies to obtain a simulated LRG sample (Obiwan LRGs) that predicts the variations in target density due to imaging properties. We find that the simulations predict the trends with depth observed in the data, including how they depend on the intrinsic brightness of the galaxies. We observe that faint LRGs are the main contributing source of the imaging systematics trend induced by depth. We also find significant trends in the data against Galactic extinction that are not predicted by Obiwan. These trends depend strongly on the particular map of Galactic extinction chosen to test against, implying systematic contamination in the Galactic extinction maps is a likely root cause (e.g., Cosmic-Infrared Background, dust temperature correction). We additionally observe a morphological change of the DESI LRGs population evidenced by a correlation between OII emission line average intensity and the size of the z-band PSF. This effect most likely results from uncertainties in background subtraction. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRGs sample.
We present a novel analysis for cluster cosmology that fully forward models the abundances, weak lensing, and the clustering of galaxy clusters. Our analysis notably includes an empirical model for the anisotropic boosts impacting the lensing and clustering signals of optical clusters. These boosts arise from a preferential selection of clusters surrounded by anisotropic large-scale structure, a consequence of the limited discrimination between line-of-sight interlopers and true cluster members offered by photometric surveys. We validate our analysis via a blind cosmology challenge on mocks, and find that we can obtain tight and unbiased cosmological constraints without informative priors or external calibrations on any of our model parameters. We then apply our analysis on the SDSS redMaPPer clusters, and find results favoring low Ωm and high σ8, combining to yield the lensing strength constraint $S_8 = 0.715_{-0.021}^{+0.024}$. We investigate potential drivers behind these results through a series of post-unblinding tests, noting that our results are consistent with existing cluster cosmology constraints but clearly inconsistent with other CMB/LSS based cosmology results. From these tests, we find hints that a suppression in the cluster lensing signal may be driving our results.
Reconstructing the initial conditions of the Universe from late-time observations has the potential to optimally extract cosmological information. Due to the high dimensionality of the parameter space, a differentiable forward model is needed for convergence, and recent advances have made it possible to perform reconstruction with nonlinear models based on galaxy (or halo) positions. In addition to positions, future surveys will provide measurements of galaxies' peculiar velocities through the kinematic Sunyaev-Zel'dovich effect (kSZ), type Ia supernovae, the fundamental plane relation, and the Tully-Fisher relation. Here we develop the formalism for including halo velocities, in addition to halo positions, to enhance the reconstruction of the initial conditions. We show that using velocity information can significantly improve the reconstruction accuracy compared to using only the halo density field. Here we study this improvement as a function of shot noise, velocity measurement noise, and angle to the line of sight. We also show how halo velocity data can be used to improve the reconstruction of the final nonlinear matter overdensity and velocity fields. We have built our pipeline into the differentiable Particle-Mesh FlowPM package, paving the way to perform field-level cosmological inference with joint velocity and density reconstruction. This is especially useful given the increased ability to measure peculiar velocities in the near future.
In the Double Shell Inertial Confinement Fusion concept, characterizing the shape asymmetry of imploding metal shells is vital for understanding energy-efficient compression and radiative losses of the thermonuclear fuel. The Monte Carlo N-Particle MCNP® code forward models radiography of Double Shell capsule implosions using the Advanced Radiographic Capability at the National Ignition Facility. A procedure is developed for using MCNP to reconstruct density profiles from the radiograph image intensity. For a given Double Shell imploding target geometry, MCNP radiographs predict image contrast, which can help guide experimental design. Finally, in future work, the calculated MCNP synthetic radiographs will be compared with experimental radiographs to determine the radial and azimuthal density profiles of the Double Shell capsules.
Abstract Recent studies have shown that the observed color distributions of Type Ia supernovae (SNe Ia) can be well described by a combination of a dust distribution and an intrinsic color distribution. Using the Pantheon+ sample of 1701 SN Ia, we apply a new forward-modeling fitting method (Dust2Dust) to measure the parent dust and color distributions, including their dependence on host-galaxy mass. At each fit step, the SN Ia selection efficiency is determined from a large simulated sample that is reweighted to reflect the proposed distributions. We use five separate metrics to describe the goodness of fit: distribution of fitted light-curve color c , cosmological residual trends with c , cosmological residual scatter with c , fitted color–luminosity relationship β SALT2 , and intrinsic scatter σ int . We present the results and the uncertainty in 12-dimensional space. Furthermore, we measure that the uncertainty on this modeling propagates to an upper threshold uncertainty in the equation of state of dark energy w of 0.014(1) for the Pantheon+ cosmology analysis and contributes negligible uncertainty to the Hubble constant H 0 . The Dust2Dust code is made publicly available at https://github.com/djbrout/dustdriver .
Abstract Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo- z ) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo- z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce PZFlow, a Python package for the probabilistic forward modeling of galaxy catalogs with normalizing flows. For catalogs simulated with PZFlow, there is a natural notion of “true” redshift posteriors that can be used for photo- z validation. We use PZFlow to simulate a photometric galaxy catalog where each galaxy has a redshift, noisy photometry, shape information, and a true redshift posterior. We also demonstrate the use of an ensemble of normalizing flows for photo- z estimation. We discuss how PZFlow will be used to validate the photo- z estimation pipeline of the Dark Energy Science Collaboration, and the wider applicability of PZFlow for statistical modeling of any tabular data.
Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.