Engineering Papers⌕ Search

Engineering topics

Evrard, August E.

Publications and source records attributed to Evrard, August E..

HSC-XXL: Baryon budget of the 136 XXL groups and clusters

Abstract We present our determination of the baryon budget for an X-ray-selected XXL sample of 136 galaxy groups and clusters spanning nearly two orders of magnitude in mass (M500 ∼ 1013–1015 M⊙) and the redshift range 0 ≲ z ≲ 1. Our joint analysis is based on the combination of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) weak-lensing mass measurements, XXL X-ray gas mass measurements, and HSC and Sloan Digital Sky Survey multiband photometry. We carry out a Bayesian analysis of multivariate mass-scaling relations of gas mass, galaxy stellar mass, stellar mass of brightest cluster galaxies (BCGs), and soft-band X-ray luminosity, by taking into account the intrinsic covariance between cluster properties, selection effect, weak-lensing mass calibration, and observational error covariance matrix. The mass-dependent slope of the gas mass–total mass (M500) relation is found to be $1.29_{-0.10}^{+0.16}$, which is steeper than the self-similar prediction of unity, whereas the slope of the stellar mass–total mass relation is shallower than unity; $0.85_{-0.09}^{+0.12}$. The BCG stellar mass weakly depends on cluster mass with a slope of $0.49_{-0.10}^{+0.11}$. The baryon, gas mass, and stellar mass fractions as a function of M500 agree with the results from numerical simulations and previous observations. We successfully constrain the full intrinsic covariance of the baryonic contents. The BCG stellar mass shows the larger intrinsic scatter at a given halo total mass, followed in order by stellar mass and gas mass. We find a significant positive intrinsic correlation coefficient between total (and satellite) stellar mass and BCG stellar mass and no evidence for intrinsic correlation between gas mass and stellar mass. All the baryonic components show no redshift evolution.

Akino, Daichi↗

Probing Galaxy Evolution in Massive Clusters Using ACT and DES: Splashback as a Cosmic Clock

In this work, we measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev–Zel'dovich effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius is consistent with CDM-only simulations and is located at 2.4$_{-0.4}^{+0.3}$ Mpch -1 . We split the galaxies on color and find significant differences in their profile shapes. Red and green-valley galaxies show a splashback-like minimum in their slope profile consistent with theory, while the bluest galaxies show a weak feature at a smaller radius. We develop a mapping of galaxies to subhalos in simulations and assign colors based on infall time onto their hosts. We find that the shift in location of the steepest slope and different profile shapes can be mapped to the average time of infall of galaxies of different colors. The steepest slope traces a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in clusters over 3.2 Gyr, green galaxies about 2.2 Gyr, while blue galaxies have been accreted most recently and have not reached apocenter. Using the full radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching occurs after a delay of about a gigayear and that galaxies quench rapidly thereafter with an exponential timescale of 0.6 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗

Rhapsody-G simulations I: the cool cores, hot gas and stellar content of massive galaxy clusters

We present the rhapsody-g suite of cosmological hydrodynamic zoom simulations of 10 massive galaxy clusters at the M vir ~10 15 M ⊙ scale. These simulations include cooling and subresolution models for star formation and stellar and supermassive black hole feedback. The sample is selected to capture the whole gamut of assembly histories that produce clusters of similar final mass. We present an overview of the successes and shortcomings of such simulations in reproducing both the stellar properties of galaxies as well as properties of the hot plasma in clusters. In our simulations, a long-lived cool-core/non-cool-core dichotomy arises naturally, and the emergence of non-cool cores is related to low angular momentum major mergers. Nevertheless, the cool-core clusters exhibit a low central entropy compared to observations, which cannot be alleviated by thermal active galactic nuclei feedback. For cluster scaling relations, we find that the simulations match well the M 500 –Y 500 scaling of Planck Sunyaev–Zeldovich clusters but deviate somewhat from the observed X-ray luminosity and temperature scaling relations in the sense of being slightly too bright and too cool at fixed mass, respectively. Stars are produced at an efficiency consistent with abundance-matching constraints and central galaxies have star formation rates consistent with recent observations. In conclusion, while our simulations thus match various key properties remarkably well, we conclude that the shortcomings strongly suggest an important role for non-thermal processes (through feedback or otherwise) or thermal conduction in shaping the intracluster medium.

79 ASTRONOMY AND ASTROPHYSICS↗

HOST GALAXY IDENTIFICATION FOR SUPERNOVA SURVEYS

Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey and the Large Synoptic Survey Telescope, which will discover SNe by the thousands. Spectroscopic resources are limited, and so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve photometric classification of SNe. In addition, SN luminosities are known to correlate with host-galaxy properties. Therefore, reliable identification of host galaxies is essential for cosmology and SN science. We simulate SN events and their locations within their host galaxies to develop and test methods for matching SNe to their hosts. We use both real and simulated galaxy catalog data from the Advanced Camera for Surveys General Catalog and MICECATv2.0, respectively. We also incorporate “hostless” SNe residing in undetected faint hosts into our analysis, with an assumed hostless rate of 5%. Our fully automated algorithm is run on catalog data and matches SNe to their hosts with 91% accuracy. We find that including a machine learning component, run after the initial matching algorithm, improves the accuracy (purity) of the matching to 97% with a 2% cost in efficiency (true positive rate). Although the exact results are dependent on the details of the survey and the galaxy catalogs used, the method of identifying host galaxies we outline here can be applied to any transient survey.

79 ASTRONOMY AND ASTROPHYSICS↗

VDES J2325-5229 a z = 2.7 gravitationally lensed quasar discovered using morphology-independent supervised machine learning

In this paper, we present the discovery and preliminary characterization of a gravitationally lensed quasar with a source redshift z s = 2.74 and image separation of 2.9 arcsec lensed by a foreground z l = 0.40 elliptical galaxy. Since optical observations of gravitationally lensed quasars show the lens system as a superposition of multiple point sources and a foreground lensing galaxy, we have developed a morphology-independent multi-wavelength approach to the photometric selection of lensed quasar candidates based on Gaussian Mixture Models (GMM) supervised machine learning. Using this technique and gi multicolour photometric observations from the Dark Energy Survey (DES), near-IR JK photometry from the VISTA Hemisphere Survey (VHS) and WISE mid-IR photometry, we have identified a candidate system with two catalogue components with i AB = 18.61 and i AB = 20.44 comprising an elliptical galaxy and two blue point sources. Spectroscopic follow-up with NTT and the use of an archival AAT spectrum show that the point sources can be identified as a lensed quasar with an emission line redshift of z = 2.739 ± 0.003 and a foreground early-type galaxy with z = 0.400 ± 0.002. We model the system as a single isothermal ellipsoid and find the Einstein radius θ E ~ 1.47 arcsec, enclosed mass M enc ~ 4 × 10 11 M ⊙ and a time delay of ~52 d. Finally, the relatively wide separation, month scale time delay duration and high redshift make this an ideal system for constraining the expansion rate beyond a redshift of 1.

79 ASTRONOMY AND ASTROPHYSICS↗

A Multiphase Model for the Intracluster Medium

Constraints on the clustered mass density of the universe derived from the observed population mean intracluster gas fraction of x-ray clusters may be biased by reliance on a single-phase assumption for the thermodynamic structure of the intracluster medium (ICM). We propose a descriptive model for multiphase structure in which a spherically symmetric ICM contains isobaric density perturbations with a radially dependent variance. Fixing the x-ray emission and emission weighted temperature, we explore two independently observable signatures of the model in the parameter space. For bremsstrahlung dominated emission, the central Sunyaev-Zel'dovich (SZ) decrement in the multiphase case is increased over the single-phase case and multiphase x-ray spectra in the range 0.1-20 keV are flatter in the continuum and exhibit stronger low energy emission lines than their single-phase counterpart. We quantify these effects for a fiducial 10e8 K cluster and demonstrate how the combination of SZ and x-ray spectroscopy can be used to identify a preferred location in the plane of the model parameter space. From these parameters the correct value of mean intracluster gas fraction in the multiphase model results, allowing an unbiased estimate of clustered mass density to he recovered.

Nagai, Daisuke↗

ROSAT observations of compact groups of galaxies

We have systematically analyzed a sample of 13 new and archival ROSAT Position Sensitive Proportional Counter (PSPC) observations of compact groups of galaxies: 12 Hickson compact groups plus the NCG 2300 group. We find that approximately two-thirds of the groups have extended X-ray emission and, in four of these, the emission is resolved into diffuse emission from gas at a temperature of kT approximately 1 keV in the group potential. All but one of the groups with extended emission have a spiral fraction of less than 50%. The baryon fraction of groups with diffuse emission is 5%-19%, similar to the values in clusters of galaxies. However, with a single exception (HCG 62), the gas-to-stellar mass ratio in our groups has a median value near 5%, somewhat greater than the values for individual early-type galaxies and two orders of magnitude than in clusters of galaxies. The X-ray luminosities of individual group galaxies are comparable to those of similar field galaxies, although the L(sub X)-L(sub B) relation for early-type galaxies may be flatter in compact groups than in the field.

Pildis, Rachel A.↗

A simulation of the intracluster medium with feedback from cluster galaxies

We detail method and report first results from a three-dimensional hydrodynamical and N-body simulation of the formation and evolution of a Coma-sized cluster of galaxies, with the intent of studying the history of the hot, X-ray emitting intracluster medium. Cluster gas, galaxies, and dark matter are included in the model. The galaxies and dark matter fell gravitational forces; the cluster gas also undergoes hydrodynamical effects such as shock heating and PdV work. For the first time in three dimensions, we include modeling of ejection of processed gas from the simulated galaxies by winds, including heating and heavy element enrichment. For comparison, we employ a `pure infall' simulation using the same initial conditions but with no galaxies or winds. We employ an extreme ejection history for galactic feedback in order to define the boundary of likely models. As expected, feedback raises the entropy of the intracluster gas, preventing it from collapsing to densities as high as those attained in the infall model. The effect is more pronounced in subclusters formed at high redshift. The cluster with feedback is always less X-ray luminous, but experiences more rapid luminosity evolution, than the pure infall cluster. Even employing an extreme ejection model, the final gas temperature is only approximately 15% larger than in the infall model. The radial temperature profile is very nearly isothermal within 1.5 Mpc. The cluster galaxies in the feedback model have a velocity dispersion approximately 15% lower than the dark matter. This results in the true ratio of specific energies in galaxies to gas being less than one, beta(sub spec) approximately 0.7. The infall model predicts beta(sub spec) approximately 1.2. Large excursions in these values occur over time, following the complex dynamical history of the cluster. The morphology of the X-ray emission is little affected by feedback. The emission profiles of both clusters are well described by the standard beta-model with beta(sub fit) approximately equal to 0.7 - 0.9. X-ray mass estimates based on the assumptions of hydrostatic equilibrium and the applicability of the beta-model are quite accurate in both cases. A strong, radial iron abundance gradient is present, which develops as a consequence of the steepening of the galaxy density profile over time. Spectroscopic observations using nonimaging detectors with wide (approximately 45 min) fields of view dramatically smear the gradient. Observations with arcminute resolution, made available with the ASCA satellite, would readily resolve the gradient.

Metzler, Christopher A.↗

The cosmological dependence of cluster density profiles

We use N-body simulations to study the shape of mean cluster density and velocity profiles in the nonlinear regime formed via gravitational instability. The dependence of the final structure on both cosmology and initial density field is examined, using a grid of cosmologies and scale-free initial power spectra P(k) varies as k(exp n). Einstein-de Sitter, open (Omega(sub 0) = 0.2 and 0.1) and flat, low density (Omega(sub 0) = 0.2 lambda(sub 0) = 0.8) models are examined, with initial spectral indices n = -2, -1 and 0. For each model, we stack clusters in an appropriately scaled manner to define an average density profile in the nonlinear regime. The profiles are well fit by a power law rho(r) varies as r(exp -alpha) for radii whereat the local density contrast is between 100 and 3000. This covers 99% of the cluster volume. We find a clear trend toward steeper slopes (larger alphas) with both increasing n and decreasing Omega(sub 0). The Omega(sub 0) dependence is partially masked by the n dependence; there is degeneracy in the values of alpha between the Einstein-de Sitter and flat, low-density cosmologies. However, the profile slopes in the open models are consistently higher than the Omega = 1 values for the range of n examined. Cluster density profiles are thus potentially useful cosmological diagnostics. We find no evidence for a constant density core in any of the models, although the density profiles do tend to flatten at small radii. Much of the flattening is due to the force softening required by the simulations. An attempt is made to recover the unsoftened profiles assuming angular momentum invariance. The recovered profiles in Einstein-de Sitter cosmologies are consistent with a pure power law up to the highest density contrasts (10(exp 6)) accessible with our resolution. The low-density models show significant deviation from a power law above density contrasts approximately 10(exp 5). We interpret this curvature as reflecting the non-scale-invariant nature of the background cosmology in these models. These results are at the limit of our resolution and so should be tested in the future using simulations with larger numbers of particles. Such simulations will also provide insight on the broader problem of understanding, in a statistical sense, the full phase space structure of collapsed, cosmological halos.

Crone, Mary M.↗

A comparison of cosmological hydrodynamic codes

We present a detailed comparison of the simulation results of various hydrodynamic codes. Starting with identical initial conditions based on the cold dark matter scenario for the growth of structure, with parameters h = 0.5 Omega = Omega(sub b) = 1, and sigma(sub 8) = 1, we integrate from redshift z = 20 to z = O to determine the physical state within a representative volume of size L(exp 3) where L = 64 h(exp -1) Mpc. Five indenpendent codes are compared: three of them Eulerian mesh-based and two variants of the smooth particle hydrodynamics 'SPH' Lagrangian approach. The Eulerian codes were run at N(exp 3) = (32(exp 3), 64(exp 3), 128(exp 3), and 256(exp 3)) cells, the SPH codes at N(exp 3) = 32(exp 3) and 64(exp 3) particles. Results were then rebinned to a 16(exp 3) grid with the exception that the rebinned data should converge, by all techniques, to a common and correct result as N approaches infinity. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainites in even primitive quantities such as (T), (rho(exp 2))(exp 1/2) persists at the 3%-17% level achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high-density, high-temeprature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by rho(exp 2) is only 4.5%. Examined at high resolution we suspect that the density resolution is better in the SPH codes and the thermal accuracy in low-density regions better in the Eulerian codes. In the low-density, low-temperature regions the SPH codes have poor accuracy due to statiscal effects, and the Jameson code gives the temperatures which are too high, due to overuse of artificial viscosity in these high Mach number regions. Overall the comparison allows us to better estimate errors; it points to ways of improving this current generation ofhydrodynamic codes and of suiting their use to problems which exploit their best individual features.

Kang, Hyesung↗

Optical signatures of high-redshift galaxy clusters

We combine an N-body and gasdynamic simulation of structure formation with an updated population synthesis code to explore the expected optical characteristics of a high-redshift cluster of galaxies. We examine a poor (2 keV) cluster formed in a biased, cold dark matter cosmology and employ simple, but plausible, threshold criteria to convert gas into stars. At z = 2, the forming cluster appears as a linear chain of very blue (g-r approximately equals 0) galaxies, with 15 objects brighter than r = 25 within a 1 square arcmin field of view. After 2 Gyr of evolution, the cluster viewed at z = 1 displays both freshly infalling blue galaxies and red galaxies robbed of recent accretion by interaction with the hot intracluster medium. The range in G-R colors is approximately 3 mag at z = 1, with the reddest objects lying at sites of highest galaxy density. We suggest that red, high-redshift galaxies lie in the cores of forming clusters and that their existence indicates the presence of a hot intracluster medium at redshifts z approximately equals 2. The simulated cluster viewed at z = 2 has several characteristics similar to the collection of faint, blue objects identified by Dressler et al. in a deep Hubble Space Telescope observation. The similarities provide some support for the interpretation of this collection as a high-redshift cluster of galaxies.

Evrard, August E.↗

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗

Galaxy tracers in N-body simulations

Using the method of smoothed particle hydrodynamics, we have modeled the formation of a compact group of galaxies with sufficient resolution to trace galaxies. Radiative cooling allows the baryons to dissipate their thermal energy and collapse to overdensities characteristic of real galaxies. With their cross section greatly reduced, these galaxy tracers remain distinct during cluster formation while their dark matter halos merge. In addition, the number density, the mass distribution function, and even the morphology of these objects are similar to those of observed galaxies. A viable population of galaxy tracers can be unambiguously defined.

Summers, F. J.↗

Gas dynamic simulations of galaxy formation

Results are presented from a simulation modeling the formation of a group of galaxies in a 'standard' cold, dark matter universe with delta = 1, h sub 0 = 50 km/(s(Mpc)), baryon fraction omega sub b = 0.1 and spectrum normalization sigma sub 8 = 0.6 (bias parameter b = 1.7). Initial conditions are generated within a periodic box with comoving length 16 Mpc in a manner constrained to produce a small cluster of total mass approximately 10 exp 14 solar mass. Two sets of 643 particles are used to model the dark matter and baryon fluids. Each gas particle represents 1.08 x 10 exp -8 solar mass, implying an L* galaxy is resolved by approximately 1000 particles. The system is evolved self-consistently in three dimensions using the combined N-body/hydrodynamic scheme P3MSPH up to a final redshift z = 1. Evolving to the present is prohibited by the fact that the mean density in the simulated volume is above critical and the entire volume would be going nonlinear beyond this point, We are currently analyzing another run with somewhat poorer mass resolution which was evolved to the present.

Evrard, August E.↗

Previrialization

The results of N-body experiments employing Gaussian random initial conditions designed to address the effect of previrialization recently examined by Peebles are reported. The central issue is whether the development of small-scale structure within a collapsing protocluster can significantly retard its collapse. A series of runs with progressively diminished small-scale power indicates that the collapse of objects is insensitive to the degree of clustering on small scales. It is concluded that previrialization does not have significant impact on the formation epoch of galaxies and clusters arising from gravitational clustering of initially Gaussian random density fields.

Evrard, August E.↗

Clues to galaxy activity from rich cluster simulations

New simulations of rich cluster evolution are used to evaluate the first infall hypothesis of Gunn and Dressler - the idea that the enhanced fraction of active galaxies seen in high redshift clusters is due to a one-time burst of star formation triggered by the rapid rise in external pressure as a galaxy plows into the hot intracluster medium (ICM). Using three-dimensional simulations which contain both baryonic gas and collisionless dark material, local static pressure histories for test orbits of galaxies are generated and a simple trigger threshold based on dP/dt/P sub ISM is applied to define an active fraction of the population. The results lend qualitative and quantitative support to the first infall interpretation.

Evrard, August E.↗