Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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.

At least 451 records · Page 25

Chandra measurements of gas homogeneity and turbulence at intermediate radii in the Perseus Cluster

ABSTRACT We present a Chandra study of surface brightness fluctuations in the diffuse intracluster medium of the Perseus Cluster. Our study utilizes deep, archival imaging of the cluster core as well as a new mosaic of 29 short $\, 5\,$-ks observations extending in eight different directions out to radii of $r_{500} \sim 2.2\, r_{2500}$. Under the assumption that the distribution of densities at a given radius is log-normally distributed, two important quantities can be derived from the width of the log-normal density distribution on a given spatial scale: the density bias, that is equal to the square root of the clumping factor C; and the one-component turbulent velocity, vk,1D. We forward-model all contributions to the measured surface brightness, including astrophysical and particle background components, and account for the Poisson nature of the measured signal. Measuring the distribution of surface brightness fluctuations in 1 arcmin2 regions, spanning the radial range $0.3\rm {-}2.2\, r_{2500}$ ($7.8\rm {-}57.3\,$ arcmin), we find a small to moderate average density bias of around $3{{\ \rm per\ cent}}$ at radii below $1.6\, r_{2500}$. We also infer an average turbulent velocity at these radii of V1D < 400 km s−1. Direct confirmation of our results on turbulent velocities inferred from surface brightness fluctuations should be possible using the X-ray calorimeter spectrometers to be flown aboard the XRISM and Athena observatories.

79 ASTRONOMY AND ASTROPHYSICS↗

The satellite population around luminous red galaxies in the 25 deg2 DESI Legacy Imaging Surveys Early Data Release

ABSTRACT Luminous red galaxies, or LRGs, are representative of the most massive galaxies and were originally selected in the Sloan Digital Sky Survey as good tracers of large-scale structure. They are dominated by by uniformly old stellar populations, have low star formation rates, early type morphologies, and little cold gas. Despite having old stellar populations and little in situ star formation, studies have shown that they have grown their stellar mass since z = 1, implying that they grow predominantly via the accretion of satellites. Tests of this picture have been limited because of the lack of deep imaging data sets that both covers a large enough area of the sky to contain substantial numbers of LRGs and that also is deep enough to detect faint satellites. We use the 25 deg2 Early Data Release (EDR) of the DESI Legacy Imaging Surveys to characterize the satellite galaxy population of LRGs out to z = 0.65. The DESI Legacy Imaging Surveys are comprised of grz imaging to 2–2.5 mag deeper than SDSS and with better image quality. We use a new statistical background technique to identify excess populations of putative satellite galaxies around 1823 LRGs at 0.2 < z < 0.65. In three redshift and luminosity bins we measure the numbers of satellite galaxies and their r−z colour distribution down to rest-frame g-band luminosity limits at least 3.6 times fainter than L*. In addition, we develop a forward modeling technique and apply it to constrain the mean number of satellites in each of our redshift and luminosity bins. Finally, we use these estimates to determine the amount of stellar mass growth in LRGs down to the local Universe.

Astronomy & Astrophysics↗

Cosmology from weak lensing peaks and minima with Subaru Hyper Suprime-Cam Survey first-year data

We present cosmological constraints derived from peak counts, minimum counts, and the angular power spectrum of the Subaru Hyper Suprime-Cam first-year (HSC Y1) weak lensing shear catalogue. Weak lensing peak and minimum counts contain non-Gaussian information and hence are complementary to the conventional two-point statistics in constraining cosmology. In this work, we forward-model the three summary statistics and their dependence on cosmology, using a suite of N-body simulations tailored to the HSC Y1 data. We investigate systematic and astrophysical effects including intrinsic alignments, baryon feedback, multiplicative bias, and photometric redshift uncertainties. We mitigate the impact of these systematics by applying cuts on angular scales, smoothing scales, signal-to-noise ratio bins, and tomographic redshift bins. By combining peaks, minima, and the power spectrum, assuming a flat-ΛCDM model, we obtain S 8 ≡ σ 8 $\sqrt{Ω_m/0.3}$ = 0.810$^{+0.022}_{-0.026}$⁠, a 35 per cent tighter constraint than that obtained from the angular power spectrum alone. Our results are in agreement with other studies using HSC weak lensing shear data, as well as with Planck 2018 cosmology and recent CMB lensing constraints from the Atacama Cosmology Telescope and the South Pole Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

Cosmological constraints from density-split clustering in the BOSS CMASS galaxy sample

We present a clustering analysis of the BOSS DR12 CMASS galaxy sample, combining measurements of the galaxy two-point correlation function and density-split clustering down to a scale of $1 \, h^{-1}\, \text{Mpc}$. Our theoretical framework is based on emulators trained on high-fidelity mock galaxy catalogues that forward model the cosmological dependence of the clustering statistics within an extended-ΛCDM framework, including redshift-space and Alcock–Paczynski distortions. Our base-ΛCDM analysis finds ω cdm = 0.1201 ± 0.0022, σ 8 = 0.792 ± 0.034, and n s = 0.970 ± 0.018, corresponding to fσ 8 = 0.462 ± 0.020 at z ≈ 0.525, which is in agreement with Planck 2018 predictions and various clustering studies in the literature. We test single-parameter extensions to base-ΛCDM, varying the running of the spectral index, the dark energy equation of state, and the density of mass-less relic neutrinos, finding no compelling evidence for deviations from the base model. We model the galaxy–halo connection using a halo occupation distribution framework, finding signatures of environment-based assembly bias in the data. We validate our pipeline against mock catalogues that match the clustering and selection properties of CMASS, showing that we can recover unbiased cosmological constraints even with a volume 84 times larger than the one used in this study.

79 ASTRONOMY AND ASTROPHYSICS↗

The mass profiles of dwarf galaxies from Dark Energy Survey lensing

We present a novel approach to extracting dwarf galaxies from photometric data to measure their average halo mass profile with weak lensing. We characterize their stellar mass and redshift distributions with a spectroscopic calibration sample. By combining the ${\sim} 5000\,\mathrm{deg}^2$ multiband photometry from the Dark Energy Survey and redshifts from the Satellites Around Galactic Analogs Survey with an unsupervised machine learning method, we select a low-mass galaxy sample spanning redshifts $z\lt 0.3$ and divide it into three mass bins. From low to high median mass, the bins contain [146 420, 330 146, 275 028] galaxies and have median stellar masses of $\log _{10}(M_*/\text{M}_\odot)=\left[8.52\substack{+0.57 -0.76},\, 9.02\substack{+0.50 -0.64},\, 9.49\substack{+0.50 -0.58}\right]$ . We measure the stacked excess surface mass density profiles, $\Delta \Sigma (R)$, of these galaxies using galaxy–galaxy lensing with a signal-to-noise ratio of [14, 23, 28]. Through a simulation-based forward-modelling approach, we fit the measurements to constrain the stellar-to-halo mass relation and find the median halo mass of these samples to be $\log _{10}(M_{\rm halo}/\text{M}_\odot)$ = [$10.67\substack{+0.2 -0.4}$, $11.01\substack{+0.14 -0.27}$, $11.40\substack{+0.08 -0.15}$]. The cold dark matter profiles are consistent with NFW (Navarro, Frenk, and White) profiles over scales ${\lesssim} 0.15 \, {h}^{-1}$ Mpc. We find that ${\sim} 20$ per cent of the dwarf galaxy sample are satellites. This is the first measurement of the halo profiles and masses of such a comprehensive, low-mass galaxy sample. The techniques presented here pave the way for extracting and analysing even lower mass dwarf galaxies and for more finely splitting galaxies by their properties with future photometric and spectroscopic survey data.

dark matter↗

Joint cosmic density reconstruction from photometric and spectroscopic samples

ABSTRACT We reconstruct the dark matter density field from spatially overlapping spectroscopic and photometric redshift catalogues through a field-level forward modelling approach. Instead of directly inferring the underlying density field, we find the best-fitting initial Gaussian fluctuations that will evolve into the observed cosmic volume. To account for the substantial uncertainty of photometric redshifts we employ a differentiable continuous Poisson process. As an initial test, we construct a mock based on the upcoming Prime Focus Spectrograph combined with photometric sample modelled on the Subaru Hyper Suprime-Cam. Depending on the statistic of interest, we find improvements in cosmic structure classification equivalent to 50–100 per cent more spectroscopic targets by combining relatively sparse spectroscopic with dense photometric samples.

Horowitz, B.↗

Validation of the ESPRESSO wavelength calibration using iodine absorption cell spectra

ABSTRACT High-quality wavelength calibration is crucial for science cases like radial-velocity studies of exoplanets, the search for a possible variation of fundamental constants, and the redshift drift experiment. However, for state-of-the-art spectrographs, it has become difficult to verify the wavelength calibration on sky because no astrophysical source provides spectra with sufficiently stable or accurate wavelength information. We therefore propose to use iodine absorption cells to validate the wavelength calibration. Observing a bright and featureless star through the iodine cell emulates an astrophysical target with exactly known spectral features that can be analysed like any other science target, allowing to verify the wavelength calibration derived from the internal calibration sources and to identify systematics in the data processing. As demonstration, we temporarily installed an $\mathrm{I_2}$ absorption cell at ESPRESSO. Employing a full forward modelling approach of the $\mathrm{I_2}$ spectrum, including the instrumental line-spread function, we demonstrate wavelength calibration accuracy at the level of a few $\mathrm{m\,s^{-1}}$. We also show that wavelength measurements do depend on the geometry of the light-injection into the spectrograph fibers. This highlights the importance of probing exactly the same light path as science targets, something not possible with internal calibration sources alone. We also demonstrate excellent radial-velocity stability at the ${{< 20\, \mathrm{cm\,s^{-1}}}}$ level in a full end-to-end fashion, from sky to data product. Our study therefore showcases the great potential of absorption cells for the verification and long-term monitoring of the wavelength calibration as well as the unique insights they can provide.

Schmidt, Tobias M. (ORCID:0000000248337273)↗

Mass dependence of halo baryon fractions from the kinetic Sunyaev–Zeldovich effect

ABSTRACT We detect the kinetic Sunyaev–Zeldovich imprint of peculiar motions of galaxy groups and clusters, using the photometric Dark Energy Spectroscopic Instrument Legacy Surveys together with cosmic microwave background (CMB) maps from the Atacama Cosmology Telescope (ACT). We develop a comprehensive forward model based on the AbacusSummit cosmological simulations: Mock galaxy group catalogues and synthetic kinetic Sunyaev–Zeldovich (kSZ) maps are generated, together with a reconstructed peculiar velocity field that allows for photo-z errors, redshift-space distortions, and survey masks. We investigate possible contamination from the cosmic infrared background (CIB), finding that CIB effects are subdominant to the kSZ signal in the relevant ACT frequency channel. We then predict the kSZ signal expected when stacking CMB temperature maps around groups, taking account of their estimated radial velocities. Comparing the model with observations, we are able to constrain the total baryon fraction within haloes, as well as their internal gas profiles. We find evidence for mass dependence of the halo baryon fraction within the virial radius. The gas fraction in massive groups is consistent with the universal baryon fraction, but low-mass groups ($10^{12.5} \lesssim M\, /\, h^{-1}\, \mathrm{M}_\odot \lesssim 10^{14}$) are depleted to $0.38 \pm 0.11$ times the universal baryon fraction. We find this low virial baryon fraction to be consistent with an extended gas profile, for which the total baryon content reaches the universal value well beyond the virial radius. This conclusion is consistent with previous analyses using X-ray, kSZ, and weak lensing, and plausibly reflects energetic feedback processes from the galaxies in these haloes.

Roper, Finn A [Institute for Astronomy, University↗

The effect of selection – a tale of cluster mass measurement bias induced by correlation and projection

ABSTRACT Cosmology analyses using galaxy clusters by the Dark Energy Survey have recently uncovered an issue of a previously unknown selection effect affecting weak-lensing mass estimates. In this letter, we use the Illustris-TNG simulation to demonstrate that selecting on galaxy counts induces a selection effect because of projection and correlation between different observables. We compute the weak-lensing-like projected mass estimations of dark matter haloes and examine their projected subhalo counts. In the two-dimensional (2D) projected space, haloes that are measured as more massive than truth have higher subhalo counts. Thus, projection along the line of sight creates cluster observables that are correlated with cluster mass measurement deviations, which in turn creates a mass measurement bias when the clusters are selected by this correlated observable. We demonstrate that the bias is predicted in a forward model using the observable-mass measurement correlation.

79 ASTRONOMY AND ASTROPHYSICS↗

Bayesian quantification of observability and equation of state of twin stars

The possibility of discovering twin stars, two neutron stars (NSs) with the same mass but different radii, is usually studied in forward modelings by using a restricted number of NS matter equations of state (EOSs) encapsulating a first-order phase transition from hadronic to quark matter (QM). Informing our likelihood function with the NS radius data from GW170817 and using a metamodel with nine parameters capable of mimicking most NS EOSs available in the literature, we conduct a Bayesian quantification of the observability and underlying EOSs of twin stars. Of the accepted EOSs, between 12 and 18% yield twin stars, depending on the restrictions we place on the second branch. The possibility of twin stars remains robust even under recent observational constraints. Here, we show that many of these twin star scenarios are observable with currently available levels of accuracy in measuring NS radii. We also present the marginalized posterior probability density functions (PDFs) of every EOS parameter for each of four mass-radius correlation topologies. We find that the inferred EOS depends sensitively on not only whether twin stars are present, but also the category of twin stars, indicating that the observation of twin stars would provide a strong constraint on the underlying EOS. In particular, for two coexisting hybrid stars having QM cores at different densities, the PDF for QM speed of sound squared 𝑐$^2_{qm}$ has two peaks, one below and another above the conformal limit 𝑐$^2_{qm}$ = 1/3 predicted by perturbative QCD.

QCD in nuclear reactions↗

Monte Carlo control loops for cosmic shear cosmology with DES Year 1 data

Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. Monte Carlo control loops (MCCL) is a forward modeling method designed to tackle this problem. It relies on the ultra fast image generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. Here, we present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_ℓ$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, point spread function (PSF), an approximate Bayesian computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and a fast method to estimate the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a nontomographic setting, we derive a set of $C_ℓ$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $Ω_m$, $σ_8$, and intrinsic alignment amplitude $A_{IA}$, defined as $S_8D_{IA}=σ_8(Ω_m/0.3)^{0.5}D_{IA}$, where $D_{IA}=1-0.11(A_{IA}-1)$. We find $S_8D_{IA}=0.8954_{-0.039}^{+0.054}$, where systematics are at the level of roughly 60% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Masses of Sunyaev-Zel’dovich galaxy clusters detected by the Atacama Cosmology Telescope: Stacked lensing measurements with Subaru HSC year 3 data

We present a stacked lensing analysis of 96 galaxy clusters selected by the thermal Sunyaev-Zel’dovich (SZ) effect in maps of the cosmic microwave background (CMB). We select foreground galaxy clusters with a 5 σ -level SZ threshold in CMB observations from the Atacama Cosmology Telescope, while we define background source galaxies for the lensing analysis with secure photometric redshift cuts in Year 3 data of the Subaru Hyper Suprime Cam survey. We detect the stacked lensing signal in the range of 0.1 < R [ h - 1 Mpc ] < 100 in each of three cluster redshift bins, 0.092 < z ≤ 0.445 , 0.445 < z ≤ 0.695 , and 0.695 < z ≤ 1.180 , with 32 galaxy clusters in each bin. The cumulative signal-to-noise ratios of the lensing signal are 14.6, 12.0, and 6.6, respectively. Using a halo-based forward model, we then constrain statistical relationships between the mass inferred from the SZ observation (i.e. SZ mass) and the total mass derived from our stacked lensing measurements. At the average SZ mass in the cluster sample ( 2.1 - 2.4 × 10 14 h - 1 M ⊙ ), our likelihood analysis shows that the average total mass differs from the SZ counterpart by a factor of 1.3 ± 0.2 , 1.6 ± 0.2 , and 1.6 ± 0.3 (68%) in the aforementioned redshift ranges, respectively. Our limits are consistent with previous lensing measurements, and we find that the cluster modeling choices can introduce a 1 σ -level difference in our parameter inferences.

79 ASTRONOMY AND ASTROPHYSICS↗

Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations

Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale mol­ecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.

X-ray imaging↗

Data-Driven Optimal Control of Traffic Signals for Urban Road Networks

This paper studies the issue of data-driven optimal control design for traffic signals of oversaturated urban road networks. The signal control system based on the store and forward model is generally uncontrollable for which the controllable decomposition is needed. Instead of identifying the unknown parameters like saturation rates and turning ratios, a finite number of measured trajectories can be used to parametrize the system and help directly construct a transformation matrix for Kalman controllable decomposition through the fundamental lemma of J. C. Willems. On top of that, an infinite-horizon linear quadratic regulator (LQR) problem is formulated considering the constraints of green times for traffic signals. The problem can be solved through a two-phase data-driven learning process, where one solves an infinite-horizon unconstrained LQR problem and the other solves a finite-horizon constrained LQR problem. The simulation result shows the theoretical analysis is effective and the proposed data-driven controller can yield desired performance for reducing traffic congestion.

Liu, Tong↗

TRINIDI: Time-of-Flight Resonance Imaging With Neutrons for Isotopic Density Inference

Accurate reconstruction of 2D and 3D isotope densities is a desired capability with great potential impact in applications such as evaluation and development of next-generation nuclear fuels. Neutron time-of-flight (TOF) resonance imaging offers a potential approach by exploiting the characteristic neutron absorption spectra of each isotope. However, it is a major challenge to compute quantitatively accurate images due to a variety of confounding effects such as severe Poisson noise, background scatter, beam non-uniformity, absorption non-linearity, and extended source pulse duration. We present the TRINIDI algorithm which is based on a two-step process in which we first estimate the neutron flux and background counts, and then reconstruct the areal densities of each isotope and pixel. Both components are based on the inversion of a forward model that accounts for the highly non-linear absorption, energy-dependent emission profile, and Poisson noise, while also modeling the substantial spatio-temporal variation of the background and flux. Further, to do this, we formulate the non-linear inverse problem as two optimization problems that are solved in sequence. We demonstrate on both synthetic and measured data that TRINIDI can reconstruct quantitatively accurate 2D views of isotopic areal density that can then be reconstructed into quantitatively accurate 3D volumes of isotopic volumetric density.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Conceptual Design of a Scintillator-Based Fast-Ion Loss Detector for the Wendelstein 7-X Stellarator

Here, a conceptual design of a scintillator-based fast-ion loss detector (FILD) has been developed for the Wendelstein 7-X stellarator (W7-X). Simulations using the Monte Carlo codes ASCOT5 and BEAMS3D have been performed to calculate the expected flux of neutral beam injection (NBI)-generated fast hydrogen ions onto the conceptual detector probe head. These fast-ion loss fluxes have been calculated for several magnetic field configurations as well as probe insertion positions. At the maximum insertion position, both co-and counter-going losses with high incident pitch angles are observed; however, at retracted positions, only co-going fast ions reach the probe head. The FILDSIM code has been used to optimize the geometry of the detector entrance and collimating elements to achieve a wide velocity space coverage as well as a high velocity–space resolution. A synthetic FILD signal is calculated for the expected loss distribution via forward modeling using the instrument response function. The synthetic signal is found to largely retain the velocity space features of the loss distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

danish

Danish is a forward modeling package for curvature wavefront sensing in the geometric limit.

Meyers, JoshuaE.↗

sf-23-021 SONICC (SCALE-DEPENDENT OMEGA FOR NEUTRINO INITIAL CONDITION CODES)

SONICC computes density and velocity transfer functions for use in the initial conditions of massive neutrino cosmological simulations. The transfer functions are backscaled using a scale-dependent growth function that produces self-consistent evolution for all three matter species (cold dark matter, baryons, massive neutrinos) in the Newtonian forward model of cosmological simulations.

EMBERSON, JEFFREYD↗