Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained variational analysis”

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 235 records · Page 13

THE MASS OF THE BLACK HOLE IN XTE J1118+480

We present contemporaneous, broadband, near-infrared spectroscopy (0.9-2.45 {mu}m) and H-band photometry of the black hole X-ray binary, XTE J1118+480. We determined the fractional dilution of the NIR ellipsoidal light curves of the donor star from other emission sources in the system by comparing the absorption features in the spectrum with field stars of known spectral type. We constrained the donor star spectral type to K7 V-M1 V and determined that the donor star contributed 54% {+-} 27% of the H-band flux at the epoch of our observations. This result underscores the conclusion that the donor star cannot be assumed to be the only NIR emission source in quiescent X-ray binaries. The H-band light curve shows a double-humped asymmetric modulation with extra flux at orbital phase 0.75. The light curve was fitted with a donor star model light curve, taking into account a constant second flux component based on the dilution analysis. We also fitted models that included emission from the donor star, a constant component from the accretion disk, and a phase-variable component from the bright spot where the mass accretion stream impacts the disk. These simple models with reasonable estimates for the component physical parameters can fully account for the observed light curve, including the extra emission at phase 0.75. From our fits, we constrained the binary inclination to 68 Degree-Sign {<=} i {<=} 79 Degree-Sign . This leads to a black hole mass of 6.9 M{sub Sun} {<=} M{sub BH} {<=} 8.2 M{sub Sun }. Long-term variations in the NIR light curve shape in XTE J1118+480 are similar to those seen in other X-ray binaries and demonstrate the presence of continued activity and variability in these systems even when in full quiescence.

79 ASTRONOMY AND ASTROPHYSICS↗

BROADBAND SPECTROSCOPY USING TWO SUZAKU OBSERVATIONS OF THE HMXB GX 301-2

We present the analysis of two Suzaku observations of GX 301-2 at two orbital phases after the periastron passage. Variations in the column density of the line-of-sight absorber are observed, consistent with accretion from a clumpy wind. In addition to a cyclotron resonance scattering feature (CRSF), multiple fluorescence emission lines were detected in both observations. The variations in the pulse profiles and the CRSF throughout the pulse phase have a signature of a magnetic dipole field. Using a simple dipole model we calculated the expected magnetic field values for different pulse phases and were able to extract a set of geometrical angles, loosely constraining the dipole geometry in the neutron star. From the variation of the CRSF width and energy, we found a geometrical solution for the dipole, making the inclination consistent with previously published values.

79 ASTRONOMY AND ASTROPHYSICS↗

Spitzer/MIPS 24 {mu}m OBSERVATIONS OF HD 209458b: THREE ECLIPSES, TWO AND A HALF TRANSITS, AND A PHASE CURVE CORRUPTED BY INSTRUMENTAL SENSITIVITY VARIATIONS

We report the results of an analysis of all Spitzer/MIPS 24 {mu}m observations of HD 209458b, one of the touchstone objects in the study of irradiated giant planet atmospheres. Altogether, we analyze two and a half transits, three eclipses, and a 58 hr near-continuous observation designed to detect the planet's thermal phase curve. The results of our analysis are: (1) a mean transit depth of 1.484% {+-} 0.033%, consistent with previous measurements and showing no evidence of variability in transit depth at the 3% level. (2) A mean eclipse depth of 0.338% {+-} 0.026%, somewhat higher than that previously reported for this system; this new value brings observations into better agreement with models. From this eclipse depth we estimate an average dayside brightness temperature of 1320 {+-} 80 K; the dayside flux shows no evidence of variability at the 12% level. (3) Eclipses in the system occur 32 {+-} 129 s earlier than would be expected from a circular orbit, which constrains the orbital quantity ecos {omega} to be 0.00004 {+-} 0.00033. This result is fully consistent with a circular orbit and sets an upper limit of 140 m s{sup -1} (3{sigma}) on any eccentricity-induced velocity offset during transit. The phase curve observations (including one of the transits) exhibit an anomalous trend similar to the detector ramp seen in previous Spitzer/IRAC observations; by modeling this ramp we recover the system parameters for this transit. The long-duration photometry which follows the ramp and transit exhibits a gradual {approx}0.2% decrease in flux over {approx}30 hr. This effect is similar to that seen in pre-launch calibration data taken with the 24 {mu}m array and is better fit by an instrumental model than a model invoking planetary emission. The large uncertainties associated with this poorly understood, likely instrumental effect prevent us from usefully constraining the planet's thermal phase curve. Our observations highlight the need for a thorough understanding of detector-related instrumental effects on long timescales when making the high-precision mid-infrared measurements planned for future missions such as EChO, SPICA, and the James Webb Space Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

The Variable Hard X-Ray Emission of NGC4945 as Observed by NuSTAR

We present a broadband (approx. 0.5 - 79 keV) spectral and temporal analysis of multiple NuSTAR observations combined with archival Suzaku and Chandra data of NGC4945, the brightest extragalactic source at 100 keV. We observe hard X-ray (> 10 keV) flux and spectral variability, with flux variations of a factor 2 on timescales of 20 ksec. A variable primary continuum dominates the high energy spectrum (> 10 keV) in all the states, while the reflected/scattered flux which dominates at E< 10 keV stays approximately constant. From modelling the complex reflection/transmission spectrum we derive a Compton depth along the line of sight of Thomson approx.2.9, and a global covering factor for the circumnuclear gas of approx. 0.15. This agrees with the constraints derived from the high energy variability, which implies that most of the high energy flux is transmitted, rather that Compton-scattered. This demonstrates the effectiveness of spectral analysis in constraining the geometric properties of the circumnuclear gas, and validates similar methods used for analyzing the spectra of other bright, Compton-thick AGN. The lower limits on the e-folding energy are between 200 − 300 keV, consistent with previous BeppoSAX, Suzaku and Swift BAT observations. The accretion rate, estimated from the X-ray luminosity and assuming a bolometric correction typical of type 2 AGN, is in the range approx. 0.1 − 0.3 lambda(sub Edd) depending on the flux state. The substantial observed X-ray luminosity variability of NGC4945 implies that large errors can arise from using single-epoch X-ray data to derive L/L(sub Edd) values for obscured AGNs.

galaxies: active aEuro" galaxies: individual (NGC4↗

Subgrid-scale horizontal and vertical variation of cloud water in stratocumulus clouds: a case study based on LES and comparisons with in situ observations

Abstract. Stratocumulus clouds in the marine boundary layer cover a large fraction of ocean surface and play an important role in the radiative energy balance of the Earth system. Simulating these clouds in Earth system models (ESMs) has proven to be extremely challenging, in part because cloud microphysical processes such as the autoconversion of cloud water into precipitation occur at scales much smaller than typical ESM grid sizes. An accurate autoconversion parameterization needs to account for not only the local microphysical process (e.g., the dependence on cloud water content qc and cloud droplet number concentration Nc) but also the subgrid-scale variability of the cloud properties that determine the process rate. Accounting for subgrid-scale variability is often achieved by the introduction of a so-called enhancement factor E. Previous studies of E for autoconversion have focused more on its dependence on cloud regime and ESM grid size, but they have largely overlooked the vertical dependence of E within the cloud. In this study, we use a large-eddy simulation (LES) model, initialized and constrained with in situ and surface-based measurements from a recent airborne field campaign, to characterize the vertical dependence of the horizontal variation of qc in stratocumulus clouds and the implications for E. Similar to our recent observational study (Zhang et al., 2021), we found that the inverse relative variance of qc, an index of horizontal homogeneity, generally increases from cloud base upward through the lower two-thirds of the cloud and then decreases in the uppermost one-third of the cloud. As a result, E decreases from cloud base upward and then increases towards the cloud top. We apply a decomposition analysis to the LES cloud water field to understand the relative roles of the mean and variances of qc in determining the vertical dependence of E. Our analysis reveals that the vertical dependence of the horizontal qc variability and enhancement factor E is a combined result of condensational growth throughout the lower portion of the cloud and entrainment mixing at cloud top. The findings of this study indicate that a vertically dependent E should be used in ESM autoconversion parameterizations.

54 ENVIRONMENTAL SCIENCES↗

Comet P/Halley's periodic brightness variations in 1910

Periodic variations in Comet Halley's production of gas and dust during its 1910 apparition were searched for using measurements from a homogeneous set of photographic plates obtained at the National Observatory of Argentina at Cordoba. After applying corrections to the observations to compensate for geometric and heliocentric distance effects, the data, spanning the two months immediately following perihelion, were split into two sections for analysis. Of the possible periods indicated in either half of the dataset, only those in the range of 7.3-7.8 days were present in both halves. The additional constraint of similarity in the shape of the resulting, phased light curves further constrained the period from about 7.35 to 7.60 days. Allowing for the inherent limitations of the 1910 data, comparison of these light curves to the equivalent 1986 light curve indicates that a match in features between the two apparitions is made possible by utilizng one of two phase shifts. While not conclusive, this ability to match features suggests long-term stability in many aspects of Halley's intrinsic activity.

Schleicher, David G.↗

How reliable are coronal hole temperatures deduced from observations?

Given the importance of the temperature at the base of the corona for the modeling of the solar wind, we investigate the range of temperatures which have been deduced from remote measurements in coronal holes, within a heliocentric distance of 1.6 solar radii, and the accuracy to which these temperatures have been inferred. Results are presented from an analysis of EUV observations made simultaneously in three spectral lines at the limb in a polar coronal hole, with little contamination from quiet region emission. A temperature range of 7.8-9.3 x 10 exp 5 K is obtained, between 1.02 to 1.07 solar radii for the coronal hole, with a very different temperature range of 9.4 x 10 exp 5 - 1.2 x 10 exp 6 K for the quiet regions bordering it. Inhomogeneities within the coronal hole contribute to a 14 percent variation in inferred temperature. The elemental abundance, which is one of the parameters that influence the temperature inference, can in turn be significantly constrained when intensity ratios from three spectral lines are used.

Habbal, S. R.↗

Domain Decomposition for Integer Optimal Control with Total Variation Regularization

Total variation integer optimal control problems admit solutions and necessary optimality conditions via geometric variational analysis. In spite of the existence of said solutions, algorithms which solve the discretized objective suffer from high numerical cost associated with the combinatorial nature of integer programming. Hence, such methods are often limited to small and medium-sized problems. We propose a globally convergent, coordinate descent–inspired algorithm that allows tractable subproblem solutions restricted to a partition of the domain. Our decomposition method solves relatively small trust-region subproblems that modify the control variable on a subdomain only. Given nontrivial subdomain overlap, we prove that a global first-order necessary optimality condition is equivalent to a first-order necessary optimality condition per subdomain. We additionally show that a sufficient decrease is achieved on a single subdomain by way of a trust-region subproblem solver using geometric measure–theoretic arguments, which we integrate with a greedy patch selection to prove convergence of our algorithm. In conclusion, we demonstrate the practicality of our algorithm on a benchmark large-scale, PDE-constrained integer optimal control problem and find that our method is faster than the state of the art.

domain decomposition↗

Sampling strategies for subsampled segmented EPI PRF thermometry in MR guided high intensity focused ultrasound

Purpose: To investigate k-space subsampling strategies to achieve fast, large field-of-view (FOV) temperature monitoring using segmented echo planar imaging (EPI) proton resonance frequency shift thermometry for MR guided high intensity focused ultrasound (MRgHIFU) applications. Methods: Five different k-space sampling approaches were investigated, varying sample spacing (equally vs nonequally spaced within the echo train), sampling density (variable sampling density in zero, one, and two dimensions), and utilizing sequential or centric sampling. Three of the schemes utilized sequential sampling with the sampling density varied in zero, one, and two dimensions, to investigate sampling the k-space center more frequently. Two of the schemes utilized centric sampling to acquire the k-space center with a longer echo time for improved phase measurements, and vary the sampling density in zero and two dimensions, respectively. Phantom experiments and a theoretical point spread function analysis were performed to investigate their performance. Variable density sampling in zero and two dimensions was also implemented in a non-EPI GRE pulse sequence for comparison. All subsampled data were reconstructed with a previously described temporally constrained reconstruction (TCR) algorithm. Results: The accuracy of each sampling strategy in measuring the temperature rise in the HIFU focal spot was measured in terms of the root-mean-square-error (RMSE) compared to fully sampled “truth.” For the schemes utilizing sequential sampling, the accuracy was found to improve with the dimensionality of the variable density sampling, giving values of 0.65 °C, 0.49 °C, and 0.35 °C for density variation in zero, one, and two dimensions, respectively. The schemes utilizing centric sampling were found to underestimate the temperature rise, with RMSE values of 1.05 °C and 1.31 °C, for variable density sampling in zero and two dimensions, respectively. Similar subsampling schemes with variable density sampling implemented in zero and two dimensions in a non-EPI GRE pulse sequence both resulted in accurate temperature measurements (RMSE of 0.70 °C and 0.63 °C, respectively). With sequential sampling in the described EPI implementation, temperature monitoring over a 192 × 144 × 135 mm{sup 3} FOV with a temporal resolution of 3.6 s was achieved, while keeping the RMSE compared to fully sampled “truth” below 0.35 °C. Conclusions: When segmented EPI readouts are used in conjunction with k-space subsampling for MR thermometry applications, sampling schemes with sequential sampling, with or without variable density sampling, obtain accurate phase and temperature measurements when using a TCR reconstruction algorithm. Improved temperature measurement accuracy can be achieved with variable density sampling. Centric sampling leads to phase bias, resulting in temperature underestimations.

42 ENGINEERING↗

Aerosol optical hygroscopicity measurements during the 2010 CARES campaign

Measurements of the effect of water uptake on particulate light extinction or scattering made at two locations during the 2010 Carbonaceous Aerosols and Radiative Effects Study (CARES) study around Sacramento, CA are reported. The observed influence of water uptake, characterized through the dimensionless optical hygroscopicity parameter γ, is compared with calculations constrained by observed particle size distributions and size-dependent particle composition. A closure assessment has been carried out that allowed for determination of the average hygroscopic growth factors (GFs) at 85% relative humidity and the dimensionless hygroscopicity parameter &kappa; for oxygenated organic aerosol (OA) and for supermicron particles (defined here as particles with aerodynamic diameters between 1 and 2.5 microns), yielding κ = 0.1–0.15 and 0.9–1.0, respectively. The derived range of oxygenated OA &kappa; values are in line with previous observations. The relatively large values for supermicron particles is consistent with substantial contributions of sea-salt-containing particles in this size range. Analysis of time-dependent variations in the supermicron particle hygroscopicity suggest that atmospheric processing, specifically chloride displacement by nitrate and the accumulation of secondary organics on supermicron particles, can lead to substantial depression of the observed GF.

54 ENVIRONMENTAL SCIENCES↗

Planck intermediate results: XLII. Large-scale Galactic magnetic fields

Recent models for the large-scale Galactic magnetic fields in the literature have been largely constrained by synchrotron emission and Faraday rotation measures. In this paper, we use three different but representative models to compare their predicted polarized synchrotron and dust emission with that measured by the Planck satellite. We first update these models to match the Planck synchrotron products using a common model for the cosmic-ray leptons. We discuss the impact on this analysis of the ongoing problems of component separation in the Planck microwave bands and of the uncertain cosmic-ray spectrum. In particular, the inferred degree of ordering in the magnetic fields is sensitive to these systematic uncertainties, and we further show the importance of considering the expected variations in the observables in addition to their mean morphology. We then compare the resulting simulated emission to the observed dust polarization and find that the dust predictions do not match the morphology in the Planck data but underpredict the dust polarization away from the plane. We modify one of the models to roughly match both observables at high latitudes by increasing the field ordering in the thin disc near the observer. Finally, though this specific analysis is dependent on the component separation issues, we present the improved model as a proof of concept for how these studies can be advanced in future using complementary information from ongoing and planned observational projects.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of Terrestrial Water Storage Changes from GRACE and GLDAS

Since March 2002, the Gravity Recovery and Climate Experiment (GRACE) has provided first estimates of land water storage variations by monitoring the time-variable component of Earth's gravity field. Here we characterize spatial-temporal variations in terrestrial water storage changes (TWSC) from GRACE and compare them to those simulated with the Global Land Data Assimilation System (GLDAS). Additionally, we use GLDAS simulations to infer how TWSC is partitioned into snow, canopy water and soil water components, and to understand how variations in the hydrologic fluxes act to enhance or dissipate the stores. Results quantify the range of GRACE-derived storage changes during the studied period and place them in the context of seasonal variations in global climate and hydrologic extremes including drought and flood, by impacting land memory processes. The role of the largest continental river basins as major locations for freshwater redistribution is highlighted. GRACE-based storage changes are in good agreement with those obtained from GLDAS simulations. Analysis of GLDAS-simulated TWSC illustrates several key characteristics of spatial and temporal land water storage variations. Global averages of TWSC were partitioned nearly equally between soil moisture and snow water equivalent, while zonal averages of TWSC revealed the importance of soil moisture storage at low latitudes and snow storage at high latitudes. Evapotranspiration plays a key role in dissipating globally averaged terrestrial water storage. Latitudinal averages showed how precipitation dominates TWSC variations in the tropics, evapotranspiration is most effective in the midlatitudes, and snowmelt runoff is a key dissipating flux at high latitudes. Results have implications for monitoring water storage response to climate variability and change, and for constraining land model hydrology simulations.

Syed, Tajdarul H.↗

Variability in Ice Nucleating Particles Across Greater Houston Texas

The concentration and cloud-forming potential of a region's ice nucleating particle (INP) population have uncertain impacts on deep convective clouds. Specifically, ice nucleating particles (INPs) may affect various cloud properties related to the formation, lifetime, and precipitation of deep convective clouds. As part of the U.S. Department of Energy's TRacking Aerosol and Convection interaction ExpeRiment (TRACER) campaign, researchers from Texas A&M University deployed three Davis Rotating-drum Universal-size-cut Monitoring (DRUM) samplers throughout Greater Houston, Texas from June through September 2022. Ambient particles, collected at the surface with the DRUM samplers in four aerodynamic diameter size ranges (>3, 3–1.2, 1.2–0.34, and 0.34–0.15 μm), were analyzed in offline cold-stage ice nucleation experiments. The INP population in Greater Houston is complex, varying by site and day, but can be generalized by a weak to moderately efficient mode of INPs at −24°C and an efficient mode at −15°C. Analysis reveals that supermicron particles are largely responsible for ice nucleation warmer than −20°C across the region while submicron particles dominate at temperatures colder than −20°C. Additionally, significant spatial diversity in the INP population was observed, with differences in mean nucleation temperature between sites for nearly every size cut. Although INP concentrations were typically ∼0.08 L −1 at −20°C throughout the campaign, a notable region-wide increase in INP concentration for particles freezing at temperatures warmer than −20°C occurred from mid-August to mid-September. This comprehensive characterization of Greater Houston's INP population, including spatial, temporal, and particle size variations, can help constrain ice microphysics parameterizations in weather and climate models.

Thompson, Seth A. [Texas A & M Univ., College Stat↗

Study of GRB light-curve decay indices in the afterglow phase

In this work, we study the distribution of temporal power-law decay indices, α, in the gamma-ray burst (GRB) afterglow phase, fitted for 176 GRBs (139 long GRBs, 12 short GRBs with extended emission, and 25 X-ray flashes) with known redshifts. These indices are compared with the temporal decay index, α {sub W}, derived with the light-curve fitting using the Willingale et al. model. This model fitting yields similar distributions of α {sub W} to the fitted α, but for individual bursts a difference can be significant. Analysis of (α, L {sub a}) distribution, where L {sub a} is the characteristic luminosity at the end of the plateau, reveals only a weak correlation of these quantities. However, we discovered a significant regular trend when studying GRB α values along the Dainotti et al. correlation between L {sub a} and the end time of the plateau emission in the rest frame, T{sub a}{sup ∗}, hereafter LT correlation. We note a systematic variation of the α parameter distribution with luminosity for any selected T{sub a}{sup ∗}. We analyze this systematics with respect to the fitted LT correlation line, expecting that the presented trend may allow us to constrain the GRB physical models. We also attempted to use the derived correlation of α(T{sub a}) versus L{sub a}(T{sub a}) to diminish the luminosity scatter related to the variations of α along the LT distribution, a step forward in the effort of standardizing GRBs. A proposed toy model accounting for this systematics applied to the analyzed GRB distribution results in a slight increase of the LT correlation coefficient.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding PDF uncertainty in W boson mass measurements*

Abstract We study the dependence of the transverse mass distribution of charged leptons and the missing energy on parton distributions (PDFs) adapted to W boson mass measurements at the CDF and ATLAS experiments. We compare the shape variations of the distribution induced by different PDFs and find that the spread of predictions from different PDF sets can be significantly larger than the PDF uncertainty predicted by a specific PDF set. We suggest analyzing the experimental data using up-to-date PDFs to gain a better understanding of the PDF uncertainties in W boson mass measurements. We also perform a series of Lagrange multiplier scans to identify the constraints on the transverse mass distribution imposed by individual data sets in the CT18 global analysis. In the case of the CDF measurement, the distribution is mostly sensitive to d -quark PDFs in the intermediate x region, which are largely constrained by DIS and Drell-Yan data on deuteron targets and Tevatron lepton charge asymmetry data.

79 ASTRONOMY AND ASTROPHYSICS↗

Description and Evaluation of the Specified-Dynamics Experiment in the Chemistry-Climate Model Initiative

We provide an overview of the REF-C1SD specified-dynamics experiment that was conducted as part of phase 1 of the Chemistry-Climate Model Initiative (CCMI). The REF-C1SD experiment, which consisted of mainly nudged general circulation models (GCMs) constrained with (re)analysis fields, was designed to examine the influence of the large-scale circulation on past trends in atmospheric composition. The REF-C1SD simulations were produced across various model frameworks and are evaluated in terms of how well they represent different measures of the dynamical and transport circulations. In the troposphere there are large (~40 %) differences in the climatological mean distributions, seasonal cycle amplitude, and trends of the meridional and vertical winds. In the stratosphere there are similarly large (~50 %) differences in the magnitude, trends and seasonal cycle amplitude of the transformed Eulerian mean circulation and among various chemical and idealized tracers. At the same time, interannual variations in nearly all quantities are very well represented, compared to the underlying reanalyses. We show that the differences in magnitude, trends and seasonal cycle are not related to the use of different reanalysis products; rather, we show they are associated with how the simulations were implemented, by which we refer both to how the large-scale flow was prescribed and to biases in the underlying free-running models. In most cases these differences are shown to be as large or even larger than the differences exhibited by free-running simulations produced using the exact same models, which are also shown to be more dynamically consistent. Overall, our results suggest that care must be taken when using specified-dynamics simulations to examine the influence of large-scale dynamics on composition.

Clara Orbe↗

A Study of the Kepler K2 Variable EPIC 211957146 Exhibiting a Variable O’ Connell Effect

We present the multi-band photometric and spectroscopic study of an over-contact binary system, EPIC 211957146. The light curves exhibit a variable O’ Connell effect, confirmed from our observational data and the Kepler K2 data. The best photometric solution incorporating a dark spot over the primary component unveils that the system has a low-mass ratio ( q ∼ 0.17) and a high inclination ( i ∼ 85°). To confirm the solution and constrain the uncertainty, Monte-Carlo simulations are performed and the results are reported. Based on the O–C diagram analysis, we see that the variable shows a period increase at the rate of dP / dt ∼ 1.06 × 10{sup −6} days yr{sup −1}, which is higher than the theoretically predicted value. Presence of a third body having a period of ∼16.23 years is evident from the O–C diagram. No filled-in effect is observed in the H α line, while the effect is vividly present in the Na line. From the Kepler K2 data, we found that the primary and secondary minima exhibit an anti-correlated O–C variation followed by an erratic behavior. This is possibly caused by the longitudinal motion of the spot, and hence, we set a lower limit of ∼40 days for the spot modulation. We also observe a possibly associated photometric difference in the primary depth by comparing our light curves with Kepler K2 normalized light curves. This system has a low-mass ratio and a high fill-out factor, and, theoretically, such a physical configuration would lead to a merger.

79 ASTRONOMY AND ASTROPHYSICS↗

The Outer Limits of Galaxy Clusters: Observations to the Virial Radius with Suzaku, XMM,and Chandra

The outskirts of galaxy clusters, near the virial radius, remain relatively unexplored territory and yet are vital to our understanding of cluster growth, structure, and mass. In this presentation, we show the first results from a program to constrain the sate of the outer intra-cluster medium (ICM) in a large sample of galaxy clusters, exploiting the strengths of three complementary X-ray observatories: Suzaku (low, stable background), XMM-Newton (high sensitivity),and Chandra (good spatial resolution). By carefully combining observations from the cluster core to beyond r200, we are able to identify and reduce systematic uncertainties that would impede our spatial and spectral analysis using a single telescope. Our sample comprises nine clusters at z is approximately 0.1-0.2 fully covered in azimuth to beyond r200, and our analysis indicates that the ICM is not in hydrostatic equilibrium in the cluster outskirts, where we see clear azimuthal variations in temperature and surface brightness. In one of the clusters, we are able to measure the diffuse X-ray emission well beyond r200, and we find that the entropy profile and the gas fraction are consistent with expectations from theory and numerical simulations. These results stand in contrast to recent studies which point to gas clumping in the outskirts; the extent to which differences of cluster environment or instrumental effects factor in this difference remains unclear. From a broader perspective, this project will produce a sizeable fiducial data set for detailed comparison with high-resolution numerical simulations.

Miller, Eric D.↗