Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mock”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

The Lyman β forest as a cosmic thermometer

We present a comprehensive analysis of high resolution hydrodynamic simulations in terms of Lyman α and Lyman β one dimensional flux power spectra P{sub αα} and P{sub ββ}). In particular, we focus on the behaviour that the flux auto-power spectra and cross-power spectra P{sub αβ}) display when the intergalactic medium (IGM) thermal history is changed in a range of values that bracket a reference model, while cosmological parameters are kept fixed to best fit the cosmic microwave background data. We present empirical fits that describe at the sub-percent level the dependence of the power spectra on the thermal parameters. At the largest scales, the power spectra show a constant bias between each other that is set by the parameters describing the IGM thermal state. The cross-power spectrum has an oscillatory pattern and crosses zero at a scale which depends on T{sub 0}, the IGM temperature at the mean density, for reasonable values of the power-law index γ of the IGM temperature-density relation T=T{sub 0}(1+δ){sup γ-1}). By performing a Fisher matrix analysis, we find that the power spectrum P{sub ββ} is more sensitive to the thermal history than P{sub αα} alone, due to the fact that it probes denser regions than Lyman α. When we combine the power and cross spectra the constraints on γ can be improved by a factor ∼ 4, while the constraints on T{sub 0} improve by a factor of ∼ 2. We address the role of signal-to-noise and resolution by mocking realistic observations and we conclude that the framework presented in this work can significantly improve the knowledge of the IGM thermal state, which will in turn guarantee better constraints on IGM-derived cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Accurate shear measurement with faint sources

For cosmic shear to become an accurate cosmological probe, systematic errors in the shear measurement method must be unambiguously identified and corrected for. Previous work of this series has demonstrated that cosmic shears can be measured accurately in Fourier space in the presence of background noise and finite pixel size, without assumptions on the morphologies of galaxy and PSF. The remaining major source of error is source Poisson noise, due to the finiteness of source photon number. This problem is particularly important for faint galaxies in space-based weak lensing measurements, and for ground-based images of short exposure times. In this work, we propose a simple and rigorous way of removing the shear bias from the source Poisson noise. Our noise treatment can be generalized for images made of multiple exposures through MultiDrizzle. This is demonstrated with the SDSS and COSMOS/ACS data. With a large ensemble of mock galaxy images of unrestricted morphologies, we show that our shear measurement method can achieve sub-percent level accuracy even for images of signal-to-noise ratio less than 5 in general, making it the most promising technique for cosmic shear measurement in the ongoing and upcoming large scale galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Optimization of spectroscopic surveys for testing non-Gaussianity

We investigate optimization strategies to measure primordial non-Gaussianity with future spectroscopic surveys. We forecast measurements coming from the 3D galaxy power spectrum and compute constraints on primordial non-Gaussianity parameters f{sub NL} and n{sub NG}. After studying the dependence on those parameters upon survey specifications such as redshift range, area, number density, we assume a reference mock survey and investigate the trade-off between number density and area surveyed. We then define the observational requirements to reach the detection of f{sub NL} of order 1. Our results show that power spectrum constraints on non-Gaussianity from future spectroscopic surveys can improve on current CMB limits, but the multi-tracer technique and higher order correlations will be needed in order to reach an even better precision in the measurements of the non-Gaussianity parameter f{sub NL}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Approximate Bayesian computation for forward modeling in cosmology

Bayesian inference is often used in cosmology and astrophysics to derive constraints on model parameters from observations. This approach relies on the ability to compute the likelihood of the data given a choice of model parameters. In many practical situations, the likelihood function may however be unavailable or intractable due to non-gaussian errors, non-linear measurements processes, or complex data formats such as catalogs and maps. In these cases, the simulation of mock data sets can often be made through forward modeling. We discuss how Approximate Bayesian Computation (ABC) can be used in these cases to derive an approximation to the posterior constraints using simulated data sets. This technique relies on the sampling of the parameter set, a distance metric to quantify the difference between the observation and the simulations and summary statistics to compress the information in the data. We first review the principles of ABC and discuss its implementation using a Population Monte-Carlo (PMC) algorithm and the Mahalanobis distance metric. We test the performance of the implementation using a Gaussian toy model. We then apply the ABC technique to the practical case of the calibration of image simulations for wide field cosmological surveys. We find that the ABC analysis is able to provide reliable parameter constraints for this problem and is therefore a promising technique for other applications in cosmology and astrophysics. Our implementation of the ABC PMC method is made available via a public code release.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Position-dependent correlation function from the SDSS-III Baryon Oscillation Spectroscopic Survey Data Release 10 CMASS sample

We report on the first measurement of the three-point function with the position-dependent correlation function from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 10 CMASS sample. This new observable measures the correlation between two-point functions of galaxy pairs within different subvolumes, ξ-circumflex (ř,ř{sub L}), where ř{sub L} is the location of a subvolume, and the corresponding mean overdensities, δ-bar (ř{sub L}). This correlation, which we call the 'integrated three-point function', iζ(r)≡(ξ-circumflex (ř,ř{sub L})δ-bar (ř{sub L})), measures a three-point function of two short- and one long-wavelength modes, and is generated by nonlinear gravitational evolution and possibly also by the physics of inflation. The iζ(r) measured from the BOSS data lies within the scatter of those from the mock galaxy catalogs in redshift space, yielding a ten-percent-level determination of the amplitude of iζ(r). The tree-level perturbation theory in redshift space predicts how this amplitude depends on the linear and quadratic nonlinear galaxy bias parameters (b{sub 1} and b{sub 2}), as well as on the amplitude and linear growth rate of matter fluctuations (σ{sub 8} and f). Combining iζ(r) with the constraints on b{sub 1σ{sub 8}} and fσ{sub 8} from the global two-point correlation function and that on σ{sub 8} from the weak lensing signal of BOSS galaxies, we measure b{sub 2}=0.41±0.41 (68% C.L.) assuming standard perturbation theory at the tree level and the local bias model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A general reconstruction of the recent expansion history of the universe

Distance measurements are currently the most powerful tool to study the expansion history of the universe without specifying its matter content nor any theory of gravitation. Assuming only an isotropic, homogeneous and flat universe, in this work we introduce a model-independent method to reconstruct directly the deceleration function via a piecewise function. Including a penalty factor, we are able to vary continuously the complexity of the deceleration function from a linear case to an arbitrary (n+1)-knots spline interpolation. We carry out a Monte Carlo (MC) analysis to determine the best penalty factor, evaluating the bias-variance trade-off, given the uncertainties of the SDSS-II and SNLS supernova combined sample (JLA), compilations of baryon acoustic oscillation (BAO) and H(z) data. The bias-variance analysis is done for three fiducial models with different features in the deceleration curve. We perform the MC analysis generating mock catalogs and computing their best-fit. For each fiducial model, we test different reconstructions using, in each case, more than 10{sup 4} catalogs in a total of about 5× 10{sup 5}. This investigation proved to be essential in determining the best reconstruction to study these data. We show that, evaluating a single fiducial model, the conclusions about the bias-variance ratio are misleading. We determine the reconstruction method in which the bias represents at most 10% of the total uncertainty. In all statistical analyses, we fit the coefficients of the deceleration function along with four nuisance parameters of the supernova astrophysical model. For the full sample, we also fit H{sub 0} and the sound horizon r{sub s}(z{sub d}) at the drag redshift. The bias-variance trade-off analysis shows that, apart from the deceleration function, all other estimators are unbiased. Finally, we apply the Ensemble Sampler Markov Chain Monte Carlo (ESMCMC) method to explore the posterior of the deceleration function up to redshift 1.3 (using only JLA) and 2.3 (JLA+BAO+H(z)). We obtain that the standard cosmological model agrees within 3σ level with the reconstructed results in the whole studied redshift intervals. Since our method is calibrated to minimize the bias, the error bars of the reconstructed functions are a good approximation for the total uncertainty.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Initial conditions for cosmological N-body simulations of the scalar sector of theories of Newtonian, Relativistic and Modified Gravity

We present a description for setting initial particle displacements and field values for simulations of arbitrary metric theories of gravity, for perfect and imperfect fluids with arbitrary characteristics. We extend the Zel'dovich Approximation to nontrivial theories of gravity, and show how scale dependence implies curved particle paths, even in the entirely linear regime of perturbations. For a viable choice of Effective Field Theory of Modified Gravity, initial conditions set at high redshifts are affected at the level of up to 5% at Mpc scales, which exemplifies the importance of going beyond Λ-Cold Dark Matter initial conditions for modifications of gravity outside of the quasi-static approximation. In addition, we show initial conditions for a simulation where a scalar modification of gravity is modelled in a Lagrangian particle-like description. Our description paves the way for simulations and mock galaxy catalogs under theories of gravity beyond the standard model, crucial for progress towards precision tests of gravity and cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Weak lensing of large scale structure in the presence of screening

A number of alternatives to general relativity exhibit gravitational screening in the non-linear regime of structure formation. We describe a set of algorithms that can produce weak lensing maps of large scale structure in such theories and can be used to generate mock surveys for cosmological analysis. By analysing a few basic statistics we indicate how these alternatives can be distinguished from general relativity with future weak lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Probing cosmology and gravity with redshift-space distortions around voids

Cosmic voids in the large-scale structure of the Universe affect the peculiar motions of objects in their vicinity. Although these motions are difficult to observe directly, the clustering pattern of their surrounding tracers in redshift space is influenced in a unique way. This allows to investigate the interplay between densities and velocities around voids, which is solely dictated by the laws of gravity. With the help of N-body simulations and derived mock-galaxy catalogs we calculate the average density fluctuations around voids identified with a watershed algorithm in redshift space and compare the results with the expectation from general relativity and the ΛCDM model. We find linear theory to work remarkably well in describing the dynamics of voids. Adopting a Bayesian inference framework, we explore the full posterior of our model parameters and forecast the achievable accuracy on measurements of the growth rate of structure and the geometric distortion through the Alcock-Paczyński effect. Systematic errors in the latter are reduced from ∼15% to ∼5% when peculiar velocities are taken into account. The relative parameter uncertainties in galaxy surveys with number densities comparable to the SDSS MAIN (CMASS) sample probing a volume of 1h{sup −3}Gpc{sup 3} yield σ{sub f/b}/(f/b).∼2% (20%) and σ{sub D{sub A{sub H}}}/D{sub AH∼0}.2% (2%), respectively. At this level of precision the linear-theory model becomes systematics dominated, with parameter biases that fall beyond these values. Nevertheless, the presented method is highly model independent; its viability lies in the underlying assumption of statistical isotropy of the Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Simulated Milky Way analogues: implications for dark matter indirect searches

We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the 'Evolution and Assembly of GaLaxies and their Environments' (EAGLE) project, and identify those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

The local value of H{sub 0} in an inhomogeneous universe

The effects of local inhomogeneities on low redshift H{sub 0} determinations are studied by estimating the redshift-distance relation of mock sources in N-body simulations. The results are compared to those obtained using the standard approach based on Hubble's law. The comparison shows a clear tendency for the standard approach to yield lower values of H{sub 0} than the approach based on the scheme using light rays. The difference is, however, small. More precisely, it is found that the overall effect of inhomogeneities on the determination of H{sub 0} is a small increase in the local estimates of about 0.3% compared to the results obtained with Hubble's law, when based on a typical distribution of supernovae in the redshift range 0.01 < z < 0.1. The overall conclusion of the study is a verification of the results that have earlier been obtained by using Hubble's law: the effects of inhomogeneities on local H{sub 0} estimates are not significant enough to make it plausible that differences in high- and low-redshift estimates of H{sub 0} are due to small inhomogeneities within the setting of standard cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Cosmic reionization study: principle component analysis after Planck

The study of reionization history plays an important role in understanding the evolution of our universe. It is commonly believed that the intergalactic medium (IGM) in our universe are fully ionized today, however the reionizing process remains to be mysterious. A simple instantaneous reionization process is usually adopted in modern cosmology without direct observational evidence. However, the history of ionization fraction, x{sub e}(z) will influence CMB observables and constraints on optical depth τ. With the mocked future data sets based on featured reionization model, we find the bias on τ introduced by instantaneous model can not be neglected. In this paper, we study the cosmic reionization history in a model independent way, the so called principle component analysis (PCA) method, and reconstruct x{sub e} (z) at different redshift z with the data sets of Planck, WMAP 9 years temperature and polarization power spectra, combining with the baryon acoustic oscillation (BAO) from galaxy survey and type Ia supernovae (SN) Union 2.1 sample respectively. The results show that reconstructed x{sub e}(z) is consistent with instantaneous behavior, however, there exists slight deviation from this behavior at some epoch. With PCA method, after abandoning the noisy modes, we get stronger constraints, and the hints for featured x{sub e}(z) evolution could become a little more obvious.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Joint constraints on galaxy bias and σ{sub 8} through the N-pdf of the galaxy number density

We present a full description of the N-probability density function of the galaxy number density fluctuations. This N-pdf is given in terms, on the one hand, of the cold dark matter correlations and, on the other hand, of the galaxy bias parameter. The method relies on the assumption commonly adopted that the dark matter density fluctuations follow a local non-linear transformation of the initial energy density perturbations. The N-pdf of the galaxy number density fluctuations allows for an optimal estimation of the bias parameter (e.g., via maximum-likelihood estimation, or Bayesian inference if there exists any a priori information on the bias parameter), and of those parameters defining the dark matter correlations, in particular its amplitude (σ{sub 8}). It also provides the proper framework to perform model selection between two competitive hypotheses. The parameters estimation capabilities of the N-pdf are proved by SDSS-like simulations (both, ideal log-normal simulations and mocks obtained from Las Damas simulations), showing that our estimator is unbiased. We apply our formalism to the 7th release of the SDSS main sample (for a volume-limited subset with absolute magnitudes M{sub r} ≤ −20). We obtain b-circumflex = 1.193 ± 0.074 and σ-bar{sub 8} = 0.862 ± 0.080, for galaxy number density fluctuations in cells of the size of 30h{sup −1}Mpc. Different model selection criteria show that galaxy biasing is clearly favoured.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Wavelet-based techniques for the gamma-ray sky

Here, we demonstrate how the image analysis technique of wavelet decomposition can be applied to the gamma-ray sky to separate emission on different angular scales. New structures on scales that differ from the scales of the conventional astrophysical foreground and background uncertainties can be robustly extracted, allowing a model-independent characterization with no presumption of exact signal morphology. As a test case, we generate mock gamma-ray data to demonstrate our ability to extract extended signals without assuming a fixed spatial template. For some point source luminosity functions, our technique also allows us to differentiate a diffuse signal in gamma-rays from dark matter annihilation and extended gamma-ray point source populations in a data-driven way.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Cosmology and the neutrino mass ordering

We propose a simple method to quantify a possible exclusion of the inverted neutrino mass ordering from cosmological bounds on the sum of the neutrino masses. The method is based on Bayesian inference and allows for a calculation of the posterior odds of normal versus inverted ordering. We apply the method for a specific set of current data from Planck CMB data and large-scale structure surveys, providing an upper bound on the sum of neutrino masses of 0.14 eV at 95% CL. With this analysis we obtain posterior odds for normal versus inverted ordering of about 2:1. If cosmological data is combined with data from oscillation experiments the odds reduce to about 3:2. For an exclusion of the inverted ordering from cosmology at more than 95% CL, an accuracy of better than 0.02 eV is needed for the sum. We demonstrate that such a value could be reached with planned observations of large scale structure by analysing artificial mock data for a EUCLID-like survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Measuring the power spectrum of dark matter substructure using strong gravitational lensing

In recent years, it has become possible to detect individual dark matter subhalos near images of strongly lensed extended background galaxies. Typically, only the most massive subhalos in the strong lensing region may be detected this way. In this work, we show that strong lenses may also be used to constrain the much more numerous population of lower mass subhalos that are too small to be detected individually. In particular, we show that the power spectrum of projected density fluctuations in galaxy halos can be measured using strong gravitational lensing. We develop the mathematical framework of power spectrum estimation, and test our method on mock observations. We use our results to determine the types of observations required to measure the substructure power spectrum with high significance. We predict that deep observations (∼10 hours on a single target) with current facilities can measure this power spectrum at the 3σ level, with no apparent degeneracy with unknown clumpiness in the background source structure or fluctuations from detector noise. Upcoming ALMA measurements of strong lenses are capable of placing strong constraints on the abundance of dark matter subhalos and the underlying particle nature of dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Measuring the velocity field from type Ia supernovae in an LSST-like sky survey

In a few years, the Large Synoptic Survey Telescope will vastly increase the number of type Ia supernovae observed in the local universe. This will allow for a precise mapping of the velocity field and, since the source of peculiar velocities is variations in the density field, cosmological parameters related to the matter distribution can subsequently be extracted from the velocity power spectrum. One way to quantify this is through the angular power spectrum of radial peculiar velocities on spheres at different redshifts. We investigate how well this observable can be measured, despite the problems caused by areas with no information. To obtain a realistic distribution of supernovae, we create mock supernova catalogs by using a semi-analytical code for galaxy formation on the merger trees extracted from N-body simulations. We measure the cosmic variance in the velocity power spectrum by repeating the procedure many times for differently located observers, and vary several aspects of the analysis, such as the observer environment, to see how this affects the measurements. Our results confirm the findings from earlier studies regarding the precision with which the angular velocity power spectrum can be determined in the near future. This level of precision has been found to imply, that the angular velocity power spectrum from type Ia supernovae is competitive in its potential to measure parameters such as σ{sub 8}. This makes the peculiar velocity power spectrum from type Ia supernovae a promising new observable, which deserves further attention.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The variance of the locally measured Hubble parameter explained with different estimators

We study the expected variance of measurements of the Hubble constant, H {sub 0}, as calculated in either linear perturbation theory or using non-linear velocity power spectra derived from N -body simulations. We compare the variance with that obtained by carrying out mock observations in the N-body simulations, and show that the estimator typically used for the local Hubble constant in studies based on perturbation theory is different from the one used in studies based on N-body simulations. The latter gives larger weight to distant sources, which explains why studies based on N-body simulations tend to obtain a smaller variance than that found from studies based on the power spectrum. Although both approaches result in a variance too small to explain the discrepancy between the value of H {sub 0} from CMB measurements and the value measured in the local universe, these considerations are important in light of the percent determination of the Hubble constant in the local universe.

79 ASTRONOMY AND ASTROPHYSICS↗