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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

SIMULATED FARADAY ROTATION MEASURES TOWARD HIGH GALACTIC LATITUDES

We study the Faraday rotation measure (RM) due to the Galactic magnetic field (GMF) toward high Galactic latitudes. The RM arises from the global, regular component as well as from the turbulent, random component of the GMF. We model the former based on observations and the latter using the data of magnetohydrodynamic turbulence simulations. For a large number of different GMF models, we produce mock RM maps around the Galactic poles and calculate various statistical quantities with the RM maps. We find that the observed medians of RMs toward the north and south Galactic poles, {approx}0.0 {+-} 0.5 rad m{sup -2} and {approx} + 6.3 {+-} 0.5 rad m{sup -2}, are difficult to explain with any of our many alternate GMF models. The standard deviation of observed RMs, {approx}9 rad m{sup -2}, is clearly larger than that of simulated RMs. The second-order structure function of observed RMs is substantially larger than that of simulated RMs, especially at small angular scales. We discuss other possible contributions to RM toward high Galactic latitudes. Besides observational errors and the intrinsic RM of background radio sources against which RM is observed, we suggest that the RM due to the intergalactic magnetic field may account for a substantial fraction of the observed RM. Finally, we note that reproducing the observed medians may require additional components or/and structures of the GMF that are not present in our models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

IMPROVING THE PRECISION OF TIME-DELAY COSMOGRAPHY WITH OBSERVATIONS OF GALAXIES ALONG THE LINE OF SIGHT

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line-of-sight convergence {kappa}{sub ext}. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with {kappa}{sub ext} obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of {kappa}{sub ext} from {sigma}{sub {kappa}} {approx} 0.06 to {approx}0.04 for very overdense LOSs like that of the system B1608+656. For more common LOSs, {sigma}{sub {kappa}} is reduced to {approx}<0.03, corresponding to an uncertainty of {approx}< 3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain {kappa}{sub ext} almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable {kappa}{sub ext} estimator gives {sigma}{sub {kappa}} = 0.047 down from 0.065 using only galaxy counts. Although deeper multiband observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

IDENTIFICATION OF MEMBERS IN THE CENTRAL AND OUTER REGIONS OF GALAXY CLUSTERS

The caustic technique measures the mass of galaxy clusters in both their virial and infall regions and, as a byproduct, yields the list of cluster galaxy members. Here we use 100 galaxy clusters with mass M{sub 200} {>=} 10{sup 14} h {sup -1} M{sub Sun} extracted from a cosmological N-body simulation of a {Lambda}CDM universe to test the ability of the caustic technique to identify the cluster galaxy members. We identify the true three-dimensional members as the gravitationally bound galaxies. The caustic technique uses the caustic location in the redshift diagram to separate the cluster members from the interlopers. We apply the technique to mock catalogs containing 1000 galaxies in the field of view of 12 h {sup -1} Mpc on a side at the cluster location. On average, this sample size roughly corresponds to 180 real galaxy members within 3r{sub 200}, similar to recent redshift surveys of cluster regions. The caustic technique yields a completeness, the fraction of identified true members, f{sub c} = 0.95 {+-} 0.03, within 3r{sub 200}. The contamination, the fraction of interlopers in the observed catalog of members, increases from f{sub i}=0.020{sup +0.046}{sub -0.015} at r{sub 200} to f{sub i}=0.08{sup +0.11}{sub -0.05} at 3r{sub 200}. No other technique for the identification of the members of a galaxy cluster provides such large completeness and small contamination at these large radii. The caustic technique assumes spherical symmetry and the asphericity of the cluster is responsible for most of the spread of the completeness and the contamination. By applying the technique to an approximately spherical system obtained by stacking the individual clusters, the spreads decrease by at least a factor of two. We finally estimate the cluster mass within 3r{sub 200} after removing the interlopers: for individual clusters, the mass estimated with the virial theorem is unbiased and within 30% of the actual mass; this spread decreases to less than 10% for the spherically symmetric stacked cluster.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

A PLATEAU IN THE PLANET POPULATION BELOW TWICE THE SIZE OF EARTH

We carry out an independent search of Kepler photometry for small transiting planets with sizes 0.5-8.0 times that of Earth and orbital periods between 5 and 50 days, with the goal of measuring the fraction of stars harboring such planets. We use a new transit search algorithm, TERRA, optimized to detect small planets around photometrically quiet stars. We restrict our stellar sample to include the 12,000 stars having the lowest photometric noise in the Kepler survey, thereby maximizing the detectability of Earth-size planets. We report 129 planet candidates having radii less than 6 R{sub E} found in three years of Kepler photometry (quarters 1-12). Forty-seven of these candidates are not in Batalha et al., which only analyzed photometry from quarters 1-6. We gather Keck HIRES spectra for the majority of these targets leading to precise stellar radii and hence precise planet radii. We make a detailed measurement of the completeness of our planet search. We inject synthetic dimmings from mock transiting planets into the actual Kepler photometry. We then analyze that injected photometry with our TERRA pipeline to assess our detection completeness for planets of different sizes and orbital periods. We compute the occurrence of planets as a function of planet radius and period, correcting for the detection completeness as well as the geometric probability of transit, R{sub *}/a. The resulting distribution of planet sizes exhibits a power law rise in occurrence from 5.7 R{sub E} down to 2 R{sub E} , as found in Howard et al. That rise clearly ends at 2 R{sub E} . The occurrence of planets is consistent with constant from 2 R{sub E} toward 1 R{sub E} . This unexpected plateau in planet occurrence at 2 R{sub E} suggests distinct planet formation processes for planets above and below 2 R{sub E} . We find that 15.1{sup +1.8}{sub -2.7}% of solar type stars-roughly one in six-has a 1-2 R{sub E} planet with P = 5-50 days.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

A SYSTEMATIC ANALYSIS OF CAUSTIC METHODS FOR GALAXY CLUSTER MASSES

We quantify the expected observed statistical and systematic uncertainties of the escape velocity as a measure of the gravitational potential and total mass of galaxy clusters. We focus our attention on low redshift (z {<=}0.15) clusters, where large and deep spectroscopic datasets currently exist. Utilizing a suite of Millennium Simulation semi-analytic galaxy catalogs, we find that the dynamical mass, as traced by either the virial relation or the escape velocity, is robust to variations in how dynamical friction is applied to ''orphan'' galaxies in the mock catalogs (i.e., those galaxies whose dark matter halos have fallen below the resolution limit). We find that the caustic technique recovers the known halo masses (M{sub 200}) with a third less scatter compared to the virial masses. The bias we measure increases quickly as the number of galaxies used decreases. For N{sub gal} > 25, the scatter in the escape velocity mass is dominated by projections along the line-of-sight. Algorithmic uncertainties from the determination of the projected escape velocity profile are negligible. We quantify how target selection based on magnitude, color, and projected radial separation can induce small additional biases into the escape velocity masses. Using N{sub gal} = 150 (25), the caustic technique has a per cluster scatter in ln (M|M{sub 200}) of 0.3 (0.5) and bias 1% {+-} 3{r_brace} (16% {+-} 5{r_brace}) for clusters with masses >10{sup 14} M{sub Sun} at z < 0.15.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE MASSIVE SATELLITE POPULATION OF MILKY-WAY-SIZED GALAXIES

Several occupational distributions for satellite galaxies more massive than m{sub *} Almost-Equal-To 4 Multiplication-Sign 10{sup 7} M{sub Sun} around Milky-Way (MW)-sized hosts are presented and used to predict the internal dynamics of these satellites as a function of m{sub *}. For the analysis, a large galaxy group mock catalog is constructed on the basis of (sub)halo-to-stellar mass relations fully constrained with currently available observations, namely the galaxy stellar mass function decomposed into centrals and satellites, and the two-point correlation functions at different masses. We find that 6.6% of MW-sized galaxies host two satellites in the mass range of the Small and Large Magellanic Clouds (SMC and LMC, respectively). The probabilities of the MW-sized galaxies having one satellite equal to or larger than the LMC, two satellites equal to or larger than the SMC, or three satellites equal to or larger than Sagittarius (Sgr) are Almost-Equal-To 0.26, 0.14, and 0.14, respectively. The cumulative satellite mass function of the MW, N{sub s} ({>=}m{sub *}) , down to the mass of the Fornax dwarf is within the 1{sigma} distribution of all the MW-sized galaxies. We find that MW-sized hosts with three satellites more massive than Sgr (as the MW) are among the most common cases. However, the most and second most massive satellites in these systems are smaller than the LMC and SMC by roughly 0.7 and 0.8 dex, respectively. We conclude that the distribution N{sub s} ({>=}m{sub *}) for MW-sized galaxies is quite broad, the particular case of the MW being of low frequency but not an outlier. The halo mass of MW-sized galaxies correlates only weakly with N{sub s} ({>=}m{sub *}). Then, it is not possible to accurately determine the MW halo mass by means of its N{sub s} ({>=}m{sub *}); from our catalog, we constrain a lower limit of 1.38 Multiplication-Sign 10{sup 12} M{sub Sun} at the 1{sigma} level. Our analysis strongly suggests that the abundance of massive subhalos should agree with the abundance of massive satellites in all MW-sized hosts, i.e., there is not a missing (massive) satellite problem for the {Lambda}CDM cosmology. However, we confirm that the maximum circular velocity, v{sub max}, of the subhalos of satellites smaller than m{sub *} {approx} 10{sup 8} M{sub Sun} is systematically larger than the v{sub max} inferred from current observational studies of the MW bright dwarf satellites; different from previous works, this conclusion is based on an analysis of the overall population of MW-sized galaxies. Some pieces of evidence suggest that the issue could refer only to satellite dwarfs but not to central dwarfs, then environmental processes associated with dwarfs inside host halos combined with supernova-driven core expansion should be on the basis of the lowering of v{sub max}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

AN IMPRINT OF MOLECULAR CLOUD MAGNETIZATION IN THE MORPHOLOGY OF THE DUST POLARIZED EMISSION

We describe a morphological imprint of magnetization found when considering the relative orientation of the magnetic field direction with respect to the density structures in simulated turbulent molecular clouds. This imprint was found using the Histogram of Relative Orientations (HRO), a new technique that utilizes the gradient to characterize the directionality of density and column density structures on multiple scales. We present results of the HRO analysis in three models of molecular clouds in which the initial magnetic field strength is varied, but an identical initial turbulent velocity field is introduced, which subsequently decays. The HRO analysis was applied to the simulated data cubes and mock-observations of the simulations produced by integrating the data cube along particular lines of sight. In the three-dimensional analysis we describe the relative orientation of the magnetic field B with respect to the density structures, showing that: (1) the magnetic field shows a preferential orientation parallel to most of the density structures in the three simulated cubes, (2) the relative orientation changes from parallel to perpendicular in regions with density over a critical density n{sub T} in the highest magnetization case, and (3) the change of relative orientation is largest for the highest magnetization and decreases in lower magnetization cases. This change in the relative orientation is also present in the projected maps. In conjunction with simulations, HROs can be used to establish a link between the observed morphology in polarization maps and the physics included in simulations of molecular clouds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

REIONIZATION ON LARGE SCALES. I. A PARAMETRIC MODEL CONSTRUCTED FROM RADIATION-HYDRODYNAMIC SIMULATIONS

We present a new method for modeling inhomogeneous cosmic reionization on large scales. Utilizing high-resolution radiation-hydrodynamic simulations with 2048{sup 3} dark matter particles, 2048{sup 3} gas cells, and 17 billion adaptive rays in a L = 100 Mpc h {sup –1} box, we show that the density and reionization redshift fields are highly correlated on large scales (∼> 1 Mpc h {sup –1}). This correlation can be statistically represented by a scale-dependent linear bias. We construct a parametric function for the bias, which is then used to filter any large-scale density field to derive the corresponding spatially varying reionization redshift field. The parametric model has three free parameters that can be reduced to one free parameter when we fit the two bias parameters to simulation results. We can differentiate degenerate combinations of the bias parameters by combining results for the global ionization histories and correlation length between ionized regions. Unlike previous semi-analytic models, the evolution of the reionization redshift field in our model is directly compared cell by cell against simulations and performs well in all tests. Our model maps the high-resolution, intermediate-volume radiation-hydrodynamic simulations onto lower-resolution, larger-volume N-body simulations (∼> 2 Gpc h {sup –1}) in order to make mock observations and theoretical predictions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

PREDICTING MERGER-INDUCED GAS MOTIONS IN ΛCDM GALAXY CLUSTERS

In the hierarchical structure formation model, clusters of galaxies form through a sequence of mergers and continuous mass accretion, which generate significant random gas motions especially in their outskirts where material is actively accreting. Non-thermal pressure provided by the internal gas motions affects the thermodynamic structure of the X-ray emitting intracluster plasma and introduces biases in the physical interpretation of X-ray and Sunyaev-Zeldovich effect observations. However, we know very little about the nature of gas motions in galaxy clusters. The ASTRO-H X-ray mission, scheduled to launch in 2015, will have a calorimeter capable of measuring gas motions in galaxy clusters at the level of ∼< 100 km s{sup –1}. In this work, we predict the level of merger-induced gas motions expected in the ΛCDM model using hydrodynamical simulations of galaxy cluster formation. We show that the gas velocity dispersion is larger in more massive clusters, but exhibits a large scatter. We show that systems with large gas motions are morphologically disturbed, while early forming, relaxed groups show a smaller level of gas motions. By analyzing mock ASTRO-H observations of simulated clusters, we show that such observations can accurately measure the gas velocity dispersion out to the outskirts of nearby relaxed galaxy clusters. ASTRO-H analysis of merging clusters, on the other hand, requires multi-component spectral fitting and enables unique studies of substructures in galaxy clusters by measuring both the peculiar velocities and the velocity dispersion of gas within individual sub-clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Mergers in galaxy groups. I. Structure and properties of elliptical remnants

We present collisionless simulations of dry mergers in groups of 3 to 25 galaxies to test the hypothesis that elliptical galaxies form at the centers of such groups. Mock observations of the central remnants confirm their similarity to ellipticals, despite having no dissipational component. We vary the profile of the original spiral's bulge and find that ellipticals formed from spirals with exponential bulges have too low Sersic indices. Mergers of spirals with de Vaucouleurs (classical) bulges produce remnants with larger Sersic indices correlated with luminosity, as with Sloan Digital Sky Survey ellipticals. Exponential bulge mergers are better fits to faint ellipticals, whereas classical bulge mergers better match luminous ellipticals. Similarly, luminous ellipticals are better reproduced by remnants undergoing many (>5) mergers, and fainter ellipticals by those with fewer mergers. The remnants follow tight size-luminosity and velocity dispersion-luminosity (Faber-Jackson) relations (<0.12 dex scatter), demonstrating that stochastic merging can produce tight scaling relations if the merging galaxies also follow tight scaling relations. The slopes of the size-luminosity and Faber-Jackson relations are close to observations but slightly shallower in the former case. Both relations' intercepts are offset—remnants are too large but have too low dispersions at fixed luminosity. Some remnants show substantial (v/σ > 0.1) rotational support, although most are slow rotators and few are very fast rotators (v/σ > 0.5). These findings contrast with previous studies concluding that dissipation is necessary to produce ellipticals from binary mergers of spirals. Multiple, mostly minor and dry mergers can produce bright ellipticals, whereas significant dissipation could be required to produce faint, rapidly rotating ellipticals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Cross-correlation of SDSS DR7 quasars and DR10 BOSS galaxies: The weak luminosity dependence of quasar clustering at z ∼ 0.5

We present the measurement of the two-point cross-correlation function (CCF) of 8198 Sloan Digital Sky Survey Data Release 7 quasars and 349,608 Data Release 10 CMASS galaxies from the Baryonic Oscillation Spectroscopic Survey at 0.3 < z < 0.9. The CCF can be reasonably well fit by a power-law model ξ{sub QG}(r) = (r/r {sub 0}){sup –γ} on projected scales of r{sub p} = 2-25 h {sup –1} Mpc with r {sub 0} = 6.61 ± 0.25 h {sup –1} Mpc and γ = 1.69 ± 0.07. We estimate a quasar linear bias of b{sub Q} = 1.38 ± 0.10 at (z) = 0.53 from the CCF measurements, which corresponds to a characteristic host halo mass of ∼4 × 10{sup 12} h {sup –1} M {sub ☉}, compared with a ∼10{sup 13} h {sup –1} M {sub ☉} characteristic host halo mass for CMASS galaxies. Based on the clustering measurements, most quasars at z-bar ∼0.5 are not the descendants of their higher luminosity counterparts at higher redshift, which would have evolved into more massive and more biased systems at low redshift. We divide the quasar sample in luminosity and constrain the luminosity dependence of quasar bias to be db{sub Q} /dlog L = 0.20 ± 0.34 or 0.11 ± 0.32 (depending on different luminosity divisions) for quasar luminosities –23.5 > M{sub i} (z = 2) > –25.5, implying a weak luminosity dependence of clustering for luminous quasars at z-bar ∼0.5. We compare our measurements with theoretical predictions, halo occupation distribution (HOD) models, and mock catalogs. These comparisons suggest that quasars reside in a broad range of host halos. The host halo mass distributions significantly overlap with each other for quasars at different luminosities, implying a poor correlation between halo mass and instantaneous quasar luminosity. We also find that the quasar HOD parameterization is largely degenerate such that different HODs can reproduce the CCF equally well, but with different satellite fractions and host halo mass distributions. These results highlight the limitations and ambiguities in modeling the distribution of quasars with the standard HOD approach.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Methods for Bayesian power spectrum inference with galaxy surveys

We derive and implement a full Bayesian large scale structure inference method aiming at precision recovery of the cosmological power spectrum from galaxy redshift surveys. Our approach improves upon previous Bayesian methods by performing a joint inference of the three-dimensional density field, the cosmological power spectrum, luminosity dependent galaxy biases, and corresponding normalizations. We account for all joint and correlated uncertainties between all inferred quantities. Classes of galaxies with different biases are treated as separate subsamples. This method therefore also allows the combined analysis of more than one galaxy survey. In particular, it solves the problem of inferring the power spectrum from galaxy surveys with non-trivial survey geometries by exploring the joint posterior distribution with efficient implementations of multiple block Markov chain and Hybrid Monte Carlo methods. Our Markov sampler achieves high statistical efficiency in low signal-to-noise regimes by using a deterministic reversible jump algorithm. This approach reduces the correlation length of the sampler by several orders of magnitude, turning the otherwise numerically unfeasible problem of joint parameter exploration into a numerically manageable task. We test our method on an artificial mock galaxy survey, emulating characteristic features of the Sloan Digital Sky Survey data release 7, such as its survey geometry and luminosity-dependent biases. These tests demonstrate the numerical feasibility of our large scale Bayesian inference frame work when the parameter space has millions of dimensions. This method reveals and correctly treats the anti-correlation between bias amplitudes and power spectrum, which are not taken into account in current approaches to power spectrum estimation, a 20% effect across large ranges in k space. In addition, this method results in constrained realizations of density fields obtained without assuming the power spectrum or bias parameters in advance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

CO component estimation based on the independent component analysis

Fast Independent Component Analysis (FastICA) is a component separation algorithm based on the levels of non-Gaussianity. Here we apply FastICA to the component separation problem of the microwave background, including carbon monoxide (CO) line emissions that are found to contaminate the PLANCK High Frequency Instrument (HFI) data. Specifically, we prepare 100 GHz, 143 GHz, and 217 GHz mock microwave sky maps, which include galactic thermal dust, NANTEN CO line, and the cosmic microwave background (CMB) emissions, and then estimate the independent components based on the kurtosis. We find that FastICA can successfully estimate the CO component as the first independent component in our deflection algorithm because its distribution has the largest degree of non-Gaussianity among the components. Thus, FastICA can be a promising technique to extract CO-like components without prior assumptions about their distributions and frequency dependences.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Measuring distances and reddenings for a billion stars: Toward a 3D dust map from Pan-STARRS 1

We present a method to infer reddenings and distances to stars based only on their broad-band photometry, and show how this method can be used to produce a three-dimensional (3D) dust map of the Galaxy. Our method samples from the full probability density function of distance, reddening, and stellar type for individual stars, as well as the full uncertainty in reddening as a function of distance in the 3D dust map. We incorporate prior knowledge of the distribution of stars in the Galaxy and the detection limits of the survey. For stars in the Pan-STARRS 1 (PS1) 3π survey, we demonstrate that our reddening estimates are unbiased and accurate to ∼0.13 mag in E(B – V) for the typical star. Based on comparisons with mock catalogs, we expect distances for main-sequence stars to be constrained to within ∼20%-60%, although this range can vary, depending on the reddening of the star, the precise stellar type, and its position on the sky. A later paper will present a 3D map of dust over the three quarters of the sky surveyed by PS1. Both the individual stellar inferences and the 3D dust map will enable a wealth of Galactic science in the plane. The method we present is not limited to the passbands of the PS1 survey but may be extended to incorporate photometry from other surveys, such as the Two Micron All Sky Survey, the Sloan Digital Sky Survey (where available), and in the future, LSST and Gaia.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Properties of the remnant clockwise disk of young stars in the galactic center

We present new kinematic measurements and modeling of a sample of 116 young stars in the central parsec of the Galaxy in order to investigate the properties of the young stellar disk. The measurements were derived from a combination of speckle and laser guide star adaptive optics imaging and integral field spectroscopy from the Keck telescopes. Compared to earlier disk studies, the most important kinematic measurement improvement is in the precision of the accelerations in the plane of the sky, which have a factor of six smaller uncertainties (σ ∼ 10 μas yr{sup –2}). We have also added the first radial velocity measurements for eight young stars, increasing the sample at the largest radii (6''-12'') by 25%. We derive the ensemble properties of the observed stars using Monte Carlo simulations of mock data. There is one highly significant kinematic feature (∼20σ), corresponding to the well-known clockwise disk, and no significant feature is detected at the location of the previously claimed counterclockwise disk. The true disk fraction is estimated to be ∼20%, a factor of ∼2.5 lower than previous claims, suggesting that we may be observing the remnant of what used to be a more densely populated stellar disk. The similarity in the kinematic properties of the B stars and the O/WR stars suggests a common star formation event. The intrinsic eccentricity distribution of the disk stars is unimodal, with an average value of (e) = 0.27 ± 0.07, which we show can be achieved through dynamical relaxation in an initially circular disk with a moderately top-heavy mass function.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Comparing simulated emission from molecular clouds using experimental design

We propose a new approach to comparing simulated observations that enables us to determine the significance of the underlying physical effects. We utilize the methodology of experimental design, a subfield of statistical analysis, to establish a framework for comparing simulated position-position-velocity data cubes to each other. We propose three similarity metrics based on methods described in the literature: principal component analysis, the spectral correlation function, and the Cramer multi-variate two-sample similarity statistic. Using these metrics, we intercompare a suite of mock observational data of molecular clouds generated from magnetohydrodynamic simulations with varying physical conditions. Using this framework, we show that all three metrics are sensitive to changing Mach number and temperature in the simulation sets, but cannot detect changes in magnetic field strength and initial velocity spectrum. We highlight the shortcomings of one-factor-at-a-time designs commonly used in astrophysics and propose fractional factorial designs as a means to rigorously examine the effects of changing physical properties while minimizing the investment of computational resources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

LyMAS: Predicting large-scale Lyα forest statistics from the dark matter density field

We describe Lyα Mass Association Scheme (LyMAS), a method of predicting clustering statistics in the Lyα forest on large scales from moderate-resolution simulations of the dark matter (DM) distribution, with calibration from high-resolution hydrodynamic simulations of smaller volumes. We use the 'Horizon-MareNostrum' simulation, a 50 h {sup –1} Mpc comoving volume evolved with the adaptive mesh hydrodynamic code RAMSES, to compute the conditional probability distribution P(F{sub s} |δ {sub s}) of the transmitted flux F{sub s} , smoothed (one-dimensionally, 1D) over the spectral resolution scale, on the DM density contrast δ {sub s}, smoothed (three-dimensionally, 3D) over a similar scale. In this study we adopt the spectral resolution of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) at z = 2.5, and we find optimal results for a DM smoothing length σ = 0.3 h {sup –1} Mpc (comoving). In its simplest form, LyMAS draws randomly from the hydro-calibrated P(F{sub s} |δ {sub s}) to convert DM skewers into Lyα forest pseudo-spectra, which are then used to compute cross-sightline flux statistics. In extended form, LyMAS exactly reproduces both the 1D power spectrum and one-point flux distribution of the hydro simulation spectra. Applied to the MareNostrum DM field, LyMAS accurately predicts the two-point conditional flux distribution and flux correlation function of the full hydro simulation for transverse sightline separations as small as 1 h {sup –1} Mpc, including redshift-space distortion effects. It is substantially more accurate than a deterministic density-flux mapping ({sup F}luctuating Gunn-Peterson Approximation{sup )}, often used for large-volume simulations of the forest. With the MareNostrum calibration, we apply LyMAS to 1024{sup 3} N-body simulations of a 300 h {sup –1} Mpc and 1.0 h {sup –1} Gpc cube to produce large, publicly available catalogs of mock BOSS spectra that probe a large comoving volume. LyMAS will be a powerful tool for interpreting 3D Lyα forest data, thereby transforming measurements from BOSS and other massive quasar absorption surveys into constraints on dark energy, DM, space geometry, and intergalactic medium physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

PRIMUS: Galaxy clustering as a function of luminosity and color at 0.2 < z < 1

We present measurements of the luminosity and color-dependence of galaxy clustering at 0.2 < z < 1.0 in the Prism Multi-object Survey. We quantify the clustering with the redshift-space and projected two-point correlation functions, ξ(r{sub p} , π) and w{sub p} (r{sub p} ), using volume-limited samples constructed from a parent sample of over ∼130, 000 galaxies with robust redshifts in seven independent fields covering 9 deg{sup 2} of sky. We quantify how the scale-dependent clustering amplitude increases with increasing luminosity and redder color, with relatively small errors over large volumes. We find that red galaxies have stronger small-scale (0.1 Mpc h {sup –1} < r{sub p} < 1 Mpc h {sup –1}) clustering and steeper correlation functions compared to blue galaxies, as well as a strong color dependent clustering within the red sequence alone. We interpret our measured clustering trends in terms of galaxy bias and obtain values of b {sub gal} ≈ 0.9-2.5, quantifying how galaxies are biased tracers of dark matter depending on their luminosity and color. We also interpret the color dependence with mock catalogs, and find that the clustering of blue galaxies is nearly constant with color, while redder galaxies have stronger clustering in the one-halo term due to a higher satellite galaxy fraction. In addition, we measure the evolution of the clustering strength and bias, and we do not detect statistically significant departures from passive evolution. We argue that the luminosity- and color-environment (or halo mass) relations of galaxies have not significantly evolved since z ∼ 1. Finally, using jackknife subsampling methods, we find that sampling fluctuations are important and that the COSMOS field is generally an outlier, due to having more overdense structures than other fields; we find that 'cosmic variance' can be a significant source of uncertainty for high-redshift clustering measurements.

79 ASTRONOMY AND ASTROPHYSICS↗