Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Power Law”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modified gravity constraints from the full shape modeling of clustering measurements from DESI 2024

We present cosmological constraints on deviations from general relativity (GR) from the first-year of clustering observations from the Dark Energy Spectroscopic Instrument (DESI) in combination with other available datasets including the CMB data from Planck with CMB-lensing from Planck and ACT, BBN constraints on the physical baryon density, the galaxy weak lensing and clustering from DESY3 and supernova data from DESY5. We first consider the μ(a,k)–Σ(a,k) modified gravity (MG) parameterization (as well as η(a,k)) in a ΛCDM and a w 0 w a CDM cosmological backgrounds. Using a functional form for time-only evolution gives μ 0 = 0.11 +0.44 -0.54 from DESI(FS+BAO)+BBN and a wide prior on n s . Using DESI(FS+BAO)+CMB+DESY3+DESY5-SN, we obtain μ 0 = 0.05 ± 0.22 and Σ 0 = 0.008 ± 0.045 and similarly μ 0 = 0.02 +0.19 -0.24 and η 0 = 0.09 +0.36 -0.60 , in an ΛCDM background. In w 0 w a CDM we obtain μ 0 = -0.24 +0.32 -0.28 and Σ 0 = 0.006 ± 0.043, consistent with GR, and we still find a preference of the data for a dynamical dark energy with w 0 > -1 and w a < 0. Using functional dependencies in both time and scale gives μ 0 and Σ 0 with a same level of precision as above but other scale MG parameters remain hard to constrain. We then move to binned parameterizations in a ΛCDM background starting with two bins in redshift and obtain, μ 1 = 1.02 ± 0.13, μ 2 = 1.04 ± 0.11, Σ 1 = 1.021 ± 0.029 and Σ 2 = 1.022 +0.027 -0.023 , all consistent with the unity value of GR in the binning formalism. We then extend the analysis to combine two bins in redshift and two in scale giving 8 MG parameters that we find all consistent with GR. We note that we find here that the tension reported in previous studies about Σ 0 being inconsistent with GR when using Planck PR3 data goes away when we use the recent LoLLiPoP+HiLLiPoP likelihoods. As noted in previous studies, this seems to indicate that the tension is indeed related to the CMB lensing anomaly in PR3 which is also resolved when using the recent likelihoods. We then constrain the class of Horndeski theory in the effective field theory of dark energy approach. We consider both EFT-basis and α-basis in the analysis. Assuming a power law parameterization for the EFT function Ω, which controls non-minimal coupling, we obtain Ω 0 = 0.012 +0.001 -0.012 and s 0 = 0.996 +0.54 -0.20 from the combination of DESI(FS+BAO)+DESY5SN+CMB in a ΛCDM background, which are consistent with GR. Similar results are obtained when using the α-basis and assuming no-braiding (α B = 0) giving c M < 1.14 at 95% CL in a ΛCDM background, also in agreement with GR. However, we see a mild yet consistent indication for c B > 0 when α B is allowed to vary which will require further study to determine whether this is due to systematics or new physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Weak Lensing Mass Calibration of the ACT DR5 Galaxy Clusters with the DES Year 3 Weak Lensing Data

We use weak gravitational lensing measurements from Year 3 Dark Energy Survey data to calibrate the masses of 443 galaxy clusters selected via the Sunyaev-Zel'dovich effect from Atacama Cosmology Telescope Data Release 5 maps of the cosmic microwave background. We incorporate redshift and SZ measurements for individual clusters into a hierarchical model for the stacked lensing signals and perform Bayesian analyses to constrain the hydrostatic mass bias of the clusters. Our treatment of systematic uncertainties includes a prescription for measuring and accounting for the weak lensing boost factor, consideration of a miscentering effect, as well as marginalization over uncertainties in the source galaxy photometric redshift distributions and shear calibration. The resultant constraints on the normalization of the mass-observable relation have a precision of approximately 7%, with the mean WL halo mass of M $_{500c}$ = 5.4 × 10$^{14}$ M $_{⊙}$. We measure the bias between the true cluster mass and the mass estimated from the SZ signal based on an X-ray-calibrated scaling relation assuming hydrostatic equilibrium, to be 1 - b = 0.74$^{+0.06}$ $_{-0.05}$ over the full sample. When splitting the clusters into high (z = 0.43-0.70) and low (z = 0.15-0.43) redshift bins, we measure 1 - b = 0.58$^{+0.06}$ $_{-0.05}$ and 0.81$^{+0.08}$ $_{-0.06}$, respectively. When introducing additional freedom in redshift and mass to the hydrostatic bias model, we find that 1 - b decreases with redshift (with the power law of -1.8$^{+0.5}$ $_{-0.6}$, 99.95% confidence), consistent with findings from other recent studies, while we do not find any significant trend in mass. We also demonstrate that our result is robust against various systematics such as a scale cut, priors on baryonic and miscentering parameters, and degree of scatter in mass-observable relation. The weak-lensing mass calibration presented in this study will be a useful tool for using the ACT clusters as probes of astrophysics, and as a step towards using their abundance as a cosmological probe.

Shin, T. [Carnegie Mellon U.] (ORCID:0000000263895↗
DOE OSTI2022

Measurement of muon-induced neutron yield at the China Jinping Underground Laboratory *

Abstract Solar, terrestrial, and supernova neutrino experiments are subject to muon-induced radioactive background. The China Jinping Underground Laboratory (CJPL), with its unique advantage of a 2400 m rock coverage and long distance from nuclear power plants, is ideal for MeV-scale neutrino experiments. Using a 1-ton prototype detector of the Jinping Neutrino Experiment (JNE), we detected 343 high-energy cosmic-ray muons and (7.86 3.97) muon-induced neutrons from an 820.28-day dataset at the first phase of CJPL (CJPL-I). Based on the muon-induced neutrons, we measured the corresponding muon-induced neutron yield in a liquid scintillator to be μ −1 g −1 cm 2 at an average muon energy of 340 GeV. We provided the first study for such neutron background at CJPL. A global fit including this measurement shows a power-law coefficient of (0.75 0.02) for the dependence of the neutron yield at the liquid scintillator on muon energy.

Zhao, Lin↗
DOE OSTI2023

The dependence of tungsten fuzz layer thickness and porosity on tungsten deposition rate and helium ion fluence

Abstract Fuzz formation on a heated tungsten surface in the presence of a helium-containing plasma and tungsten deposition source was investigated. Tungsten samples were exposed at 1123 K to pure helium plasma with ion incident energy of 76 eV, W/He ion flux ratio of ∼ 0.4 × 10 − 4 , and varied helium ion fluence from 0.18 to 3.4 × 10 26 m −2 . Fuzz thickness was measured by cross-sectional scanning electron microscopy to increase from 0.22 to 15 µ m with increasing helium ion fluence. No indication of saturation in fuzz thickness at high fluence was observed, in contrast to fuzz produced on a tungsten surface without tungsten deposition. Additional tungsten samples were exposed at 1123 K to pure helium plasma with ion incident energy of 76 eV, helium ion fluence of ∼ 3.4 × 10 26 m −2 , and varied W/He ion flux ratio from 0.26 to 3.0 × 10 − 4 . Fuzz thickness increased from 7.5 to 120 µ m with increasing W/He ion ratio. A final sample exposed at 1123 K to a mixed helium-deuterium plasma with ion incident energy of 76 eV, helium ion fluence of 0.18 × 10 26 m −2 , and W/He ion flux ratio of 2.2 × 10 − 4 developed nearly identical fuzz structures to that developed in a pure He plasma. As a function of deposited tungsten fluence, all results were found to trace out a single layer-growth curve given by a power law relation, indicating that fuzz thickness is independent of the W/He ion flux ratio in the range investigated and independent of any deuterium present in the plasma. As a result, for tungsten plasma facing walls in magnetic fusion devices at 1000–2000 K with 10 −4 W/He ion flux ratio, fuzz with thicknesses greater than hundreds of microns may form in as little as 10 4 s (in the absence of ELM-induced erosion or annealing), and may more significantly affect its thermophysical properties than fuzz generated without a tungsten deposition source.

Physics↗
DOE OSTI2025

Metrics and extrapolation of resonant magnetic perturbation thresholds for ELM suppression

This large database study of resonant magnetic perturbation (RMP) edge localized mode (ELM) suppression thresholds in the AUG, DIII-D, EAST, and KSTAR tokamaks details the key strengths and weaknesses of RMP metrics. The RMP ELM suppression database used for this work contains plasma information at the time of transition from ELMing to ELM suppressed states where a clear experimental threshold is identified. The experimental threshold distributions are compared for five metrics: (1) the island overlap width, (2) pedestal top Chirikov overlap, (3) peeling edge displacement, (4) pedestal top resonant drive, and (5) edge dominant mode overlap. The distributions, the regularity of the dependence on RMP coil currents, and the sensitivities of a given metric to equilibrium reconstruction details are compared. The overlap metric proves to be a good compromise between including the appropriate plasma response physics and maintaining a numerical robustness. This quantity does not exhibit clear power-law scalings for projection, but machine learning can assist in predicting thresholds within the existing parameter ranges and providing uncertainty quantification of those predictions. Two new first-principles models, one utilizing a threshold from the non-linear Modified Rutherford equation evaluated at the pedestal top and one utilizing the SLAYER code to calculate the linear tearing threshold from torque balance, offer possible paths to extrapolation beyond the existing database parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

HIPED: Machine learning framework for spherical tokamak pedestal prediction and optimization

We introduce a Machine Learning framework, HIPED (HeIght and width Predictor for Edge Dynamics), for predicting and optimizing pedestal and core performance in spherical tokamak plasmas. Trained on pedestal and core datasets from the third MAST-U campaign, HIPED provides accurate estimates of pedestal height and width. The results reveal notable differences compared with conventional aspect-ratio studies; for instance, a simple power-law relation between pedestal width and height has very low accuracy. Instead, additional parameters such as normalized plasma pressure, elongation, and Greenwald fraction significantly improve accuracy. HIPED can also be trained only on `control room parameters' to inform experimentalists of which controllable parameters to adjust for improving core-integrated performance. The framework further includes a multi-objective optimization scheme that helps guide experimental planning and optimization. We find Pareto-optimal discharges with respect to various features, including distance from edge-localized modes and normalized plasma pressure, track their parameter trajectories over time, and identify the control room parameters required for these Pareto-optimal discharges. This provides a framework for systematically optimizing core and edge performance according to different experimental priorities.

Parisi, Jason F. [Princeton Plasma Physics Laborat↗
DOE OSTI2026

On how avalanches penetrate the SOL and broaden heat loads

Recent experiments reported a correlation between power law core temperature spectra and Dα emission, suggesting that heat avalanches penetrate the SOL. This paper derives a threshold criterion for avalanche penetration using a reduced model. Avalanches with ($\bigtriangledown$$\tilde{T}$) rms > $\bigtriangledown$$\tilde{T}$ crit at the separatrix are predicted to penetrate, and so broaden the SOL and heat load distribution. $\bigtriangledown$$\tilde{T}$ crit is ~ 1/$τ$ ∥ , where $τ$ ∥ is the parallel heat flow time through the SOL. Penetration occurs when avalanches are strong enough to steepen sufficiently to shock at the separatrix. A positive correlation is found between the nonlinear drive for steepening and the penetration depth. In particular, penetration depth exceeds that of the heuristic drift limit when shocks form. Implications for numerical and physical experiments are also discussed.

SOL heat load↗
DOE OSTI2024

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth-order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Correlated single- and few-electron backgrounds milliseconds after interactions in dual-phase liquid xenon time projection chambers

In this work, we characterize single- and few-electron backgrounds that are observed in dual-phase liquid xenon time projection chambers at timescales greatly exceeding a maximum drift time after an interaction. These instrumental backgrounds limit a detector's sensitivity to dark matter and cosmogenic neutrinos. Using the ~150g liquid xenon detector at Purdue University, we investigate how these backgrounds, produced after 122keV 57 Co Compton interactions, behave under different detector conditions. We find that the rates of single- and few-electron signals follow power-laws with time after the interaction. We observe linearly increasing rates with increased extraction field, and increased rates in the single-electron background with increased drift field. Normalizing the rates to the primary interaction's measured ionization signal, the rates increase linearly with the depth of the interaction. We test the hypothesis that infrared photons (1550nm) would stimulate and accelerate electron emission via photodetachment from impurities, but find that even 1 Watt of infrared light fails to reduce these backgrounds. We thus provide a characterization that can inform background models for low-energy rare event searches.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Precise relative magnitude measurement improves fracture characterization during hydraulic fracturing

SUMMARY Microseismic monitoring is an important technique to obtain detailed knowledge of in-situ fracture size and orientation during stimulation to maximize fluid flow throughout the rock volume and optimize production. Furthermore, considering that the frequency of earthquake magnitudes empirically follows a power law (i.e. Gutenberg–Richter), the accuracy of microseismic event magnitude distributions is potentially crucial for seismic risk management. In this study, we analyse microseismicity observed during four hydraulic fracture treatments of the legacy Cotton Valley experiment in 1997 at the Carthage gas field of East Texas, where fractures were activated at the base of the sand-shale Upper Cotton Valley formation. We perform waveform cross-correlation to detect similar event clusters, measure relative amplitude from aligned waveform pairs with a principal component analysis, then measure precise relative magnitudes. The new magnitudes significantly reduce the deviations between magnitude differences and relative amplitudes of event pairs. This subsequently reduces the magnitude differences between clusters located at different depths. Reduction in magnitude differences between clusters suggests that some attenuation-related biases could be effectively mitigated with relative magnitude measurements. The maximum likelihood method is applied to understand the magnitude frequency distributions and quantify the seismogenic index of the clusters. Statistical analyses with new magnitudes suggest that fractures that are more favourably oriented for shear failure have lower b-value and higher seismogenic index, suggesting higher potential for relatively larger earthquakes, rather than fractures subparallel to maximum horizontal principal stress orientation.

58 GEOSCIENCES↗
DOE OSTI2020

Upper boundaries of active galactic nucleus regions in optical diagnostic diagrams

The distribution of galaxies in optical diagnostic diagrams can provide information about their physical parameters when compared with ionization models under proper assumptions. By using a sample of central emitting regions from the Mapping Nearby Galaxies at Apache Point Observatory survey (MaNGA), we find evidence of the existence of upper boundaries for narrow-line regions (NLRs) of active galactic nuclei (AGNs) in optical Baldwin, Phillips & Terlevich (BPT) diagrams, especially in diagrams involving [Sii]λλ6716, 6731/Hα. Photoionization models can reproduce the boundaries well, as a consequence of the decrease of [Sii]λλ6716, 6731/Hα and [Oiii]λ5007/Hβ ratios at very high metallicity. Whilst the exact location of the upper boundary in the [Sii] BPT diagram depends only weakly on the electron density of the ionized cloud and the secondary nitrogen prescription, its dependence on the shapes of the input spectral energy distributions (SEDs) is much stronger. This allows us to constrain the power-law index of the AGN SED between 1 Ryd and ~100 Ryd to be less than or equal to -1.40 ± 0.05. The coverage of photoionization models in the [Nii] BPT diagram has a stronger dependence on the electron density and the secondary nitrogen prescription. With the density constrained by the [Sii] doublet ratio and the input SED constrained by the [Sii] BPT diagram, we find that the extent of the data in the [Nii] BPT diagram favours those prescriptions with high N/O ratios. Furthermore, although shock-ionized clouds can produce line ratios similar to those from photoionization, the resulting shapes of the upper boundaries, if they exist, would likely be different from those of photoionizing origin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Two-dimensional simulations of internal gravity waves in the radiation zones of intermediate-mass stars

Intermediate-mass main-sequence stars have large radiative envelopes overlying convective cores. This configuration allows internal gravity waves (IGWs) generated at the convective–radiative interface to propagate towards the stellar surface. The signatures of these waves can be observed in the photometric and spectroscopic data from stars. Here, we have studied the propagation of these IGWs using two-dimensional (2D) fully non-linear hydrodynamical simulations with realistic stellar reference states from the 1D stellar evolution code, Modules for Stellar Astrophysics (MESA). When a single wave is forced, we observe wave self-interaction. When two waves are forced, we observe non-linear interaction (i.e. triadic interaction) between these waves forming waves at different wavelengths and frequencies. When a spectrum of waves similar to that found in numerical simulations is forced, we find that the surface IGW frequency slope is consistent with recent observations. This power law is similar to that predicted by linear theory for the wave propagation, with small deviations that can be an effect of non-linearities. When the same generation spectrum is applied to 3 M ⊙ models at different stellar rotation and ages, the surface IGW spectrum slope is very similar to the generation spectrum slope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Discovery and follow-up of ASASSN-19dj: an X-ray and UV luminous TDE in an extreme post-starburst galaxy

Here, we present observations of ASASSN-19dj, a nearby tidal disruption event (TDE) discovered in the post-starburst galaxy KUG 0810+227 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of d ≃ 98 Mpc. We observed ASASSN-19dj from -21 to 392 d relative to peak ultraviolet (UV)/optical emission using high-cadence, multiwavelength spectroscopy and photometry. From the ASAS-SN g-band data, we determine that the TDE began to brighten on 2019 February 6.8 and for the first 16 d the rise was consistent with a flux ∝t 2 power law. ASASSN-19dj peaked in the UV/optical on 2019 March 6.5 (MJD = 58548.5) at a bolometric luminosity of L = (6.2 ± 0.2) × 10 44 erg s -1 . Initially remaining roughly constant in X-rays and slowly fading in the UV/optical, the X-ray flux increased by over an order of magnitude ~225 d after peak, resulting from the expansion of the X-ray emitting region. The late-time X-ray emission is well fitted by a blackbody with an effective radius of ~1 × 10 12 cm and a temperature of ~6 × 10 5 K. The X-ray hardness ratio becomes softer after brightening and then returns to a harder state as the X-rays fade. Analysis of Catalina Real-Time Transient Survey images reveals a nuclear outburst roughly 14.5 yr earlier with a smooth decline and a luminosity of L V ≥ 1.4 × 10 43 erg s -1 , although the nature of the flare is unknown. ASASSN-19dj occurred in the most extreme post-starburst galaxy yet to host a TDE, with Lick Hδ A = 7.67 ± 0.17 Å.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

To TDE or not to TDE: the luminous transient ASASSN-18jd with TDE-like and AGN-like qualities

ABSTRACT We present the discovery of ASASSN-18jd (AT 2018bcb), a luminous optical/ultraviolet(UV)/X-ray transient located in the nucleus of the galaxy 2MASX J22434289–1659083 at z = 0.1192. Over the year after discovery, Swift UltraViolet and Optical Telescope (UVOT) photometry shows the UV spectral energy distribution of the transient to be well modelled by a slowly shrinking blackbody with temperature $T \sim 2.5 \times 10^{4} \, {\rm K}$, a maximum observed luminosity of $L_{\rm max} = 4.5^{+0.6}_{-0.3}\times 10^{44} \, {\rm erg \,s}^{-1}$, and a radiated energy of $E = 9.6^{+1.1}_{-0.6} \times 10^{51} \, {\rm erg}$. X-ray data from Swift X-Ray Telescope (XRT) and XMM–Newton show a transient, variable X-ray flux with blackbody and power-law components that fade by nearly an order of magnitude over the following year. Optical spectra show strong, roughly constant broad Balmer emission and transient features attributable to He ii, N iii–v, O iii, and coronal Fe. While ASASSN-18jd shares similarities with tidal disruption events (TDEs), it is also similar to the newly discovered nuclear transients seen in quiescent galaxies and faint active galactic nuclei (AGNs).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A cautionary tale of attenuation in star-forming regions

ABSTRACT The attenuation of light from star-forming galaxies is correlated with a multitude of physical parameters including star formation rate, metallicity and total dust content. This variation in attenuation is even more evident on kiloparsec scales, which is the relevant size for many current spectroscopic integral field unit surveys. To understand the cause of this variation, we present and analyse Swift/UVOT near-UV (NUV) images and SDSS/MaNGA emission-line maps of 29 nearby (z < 0.084) star-forming galaxies. We resolve kiloparsec-sized star-forming regions within the galaxies and compare their optical nebular attenuation (i.e. the Balmer emission line optical depth, $\tau ^{l}_{B}\equiv \tau _{\textrm {H}\beta }-\tau _{\textrm {H}\alpha }$) and NUV stellar continuum attenuation (via the NUV power-law index, β) to the attenuation law described by Battisti et al. We show the data agree with that model, albeit with significant scatter. We explore the dependence of the scatter of the β–$\tau ^{l}_{B}$ measurements from the star-forming regions on different physical parameters, including distance from the nucleus, star formation rate and total dust content. Finally, we compare the measured $\tau ^{l}_{B}$ and β values for the individual star-forming regions with those of the integrated galaxy light. We find a strong variation in β between the kiloparsec scale and the larger galaxy scale that is not seen in $\tau ^{l}_{B}$. We conclude that the sightline dependence of UV attenuation and the reddening of β due to the light from older stellar populations could contribute to the scatter in the β–$\tau ^{l}_{B}$ relation.

Molina, Mallory↗
DOE OSTI2021

The rapid transition from star formation to AGN-dominated rest-frame ultraviolet light at z ≃ 4

ABSTRACT With the advent of deep optical-to-near-infrared extragalactic imaging on the degree scale, samples of high-redshift sources are being selected that contain both bright star-forming (SF) galaxies and faint active galactic nuclei (AGN). In this study, we investigate the transition between SF- and AGN-dominated systems at z ≃ 4 in the rest-frame ultraviolet (UV). We find a rapid transition to AGN-dominated sources brightward of MUV ≃ −23.2. The effect is observed in the rest-frame UV morphology and size–luminosity relation, where extended clumpy systems become point-source-dominated, and also in the available spectra for the sample. These results allow us to derive the rest-frame UV luminosity function (LF) for the SF- and AGN-dominated subsamples. We find the SF-dominated LF is best fit with a double power law, with a lensed Schechter function being unable to explain the existence of extremely luminous SF galaxies at MUV ≃ −23.5. If we identify AGN-dominated sources according to a point-source morphology criterion, we recover the relatively flat faint-end slope of the AGN LF determined in previous studies. If we instead separate the LF according to the current spectroscopic AGN fraction, we find a steeper faint-end slope of α = −1.83 ± 0.11. Using a simple model to predict the rest-frame AGN LF from the z = 4 galaxy LF, we find that the increasing impact of host galaxy light on the measured morphology of faint AGN can explain our observations.

Bowler, R. A. A.↗
DOE OSTI2021

Dissipative dark matter on FIRE – I. Structural and kinematic properties of dwarf galaxies

We present the first set of cosmological baryonic zoom-in simulations of galaxies including dissipative self-interacting dark matter (dSIDM). These simulations utilize the Feedback In Realistic Environments galaxy formation physics, but allow the dark matter to have dissipative self-interactions analogous to standard model forces, parametrized by the self-interaction cross-section per unit mass, (σ/m), and the dimensionless degree of dissipation, 0 < f diss < 1. We survey this parameter space, including constant and velocity-dependent cross-sections, and focus on structural and kinematic properties of dwarf galaxies with $M_{\rm halo} \sim 10^{10-11}{\, \rm M_\odot }$ and $M_{\ast } \sim 10^{5-8}{\, \rm M_\odot }$. Central density profiles (parametrized as ρ ∝ rα) of simulated dwarfs become cuspy when $(\sigma /m)_{\rm eff} \gtrsim 0.1\, {\rm cm^{2}\, g^{-1}}$ (and f diss = 0.5 as fiducial). The power-law slopes asymptote to α ≈ -1.5 in low-mass dwarfs independent of cross-section, which arises from a dark matter ‘cooling flow’. Through comparisons with dark matter only simulations, we find the profile in this regime is insensitive to the inclusion of baryons. However, when $(\sigma /m)_{\rm eff} \ll 0.1\, {\rm cm^{2}\, g^{-1}}$, baryonic effects can produce cored density profiles comparable to non-dissipative cold dark matter (CDM) runs but at smaller radii. Simulated galaxies with $(\sigma /m) \gtrsim 10\, {\rm cm^{2}\, g^{-1}}$ and the fiducial fdiss develop significant coherent rotation of dark matter, accompanied by halo deformation, but this is unlike the well-defined thin ‘dark discs’ often attributed to baryon-like dSIDM. The density profiles in this high cross-section model exhibit lower normalizations given the onset of halo deformation. For our surveyed dSIDM parameters, halo masses and galaxy stellar masses do not show appreciable difference from CDM, but dark matter kinematics and halo concentrations/shapes can differ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The history of metal enrichment traced by X-ray observations of high-redshift galaxy clusters

ABSTRACT We present the analysis of deep X-ray observations of 10 massive galaxy clusters at redshifts 1.05 < z < 1.71, with the primary goal of measuring the metallicity of the intracluster medium (ICM) at intermediate radii, to better constrain models of the metal enrichment of the intergalactic medium. The targets were selected from X-ray and Sunyaev–Zel’dovich effect surveys, and observed with both the XMM–Newton and Chandra satellites. For each cluster, a precise gas mass profile was extracted, from which the value of r500 could be estimated. This allows us to define consistent radial ranges over which the metallicity measurements can be compared. In general, the data are of sufficient quality to extract meaningful metallicity measurements in two radial bins, r < 0.3r500 and 0.3 < r/r500 < 1.0. For the outer bin, the combined measurement for all 10 clusters, Z/Z⊙ = 0.21 ± 0.09, represents a substantial improvement in precision over previous results. This measurement is consistent with, but slightly lower than, the average metallicity of 0.315 solar measured at intermediate-to-large radii in low-redshift clusters. Combining our new high-redshift data with the previous low-redshift results allows us to place the tightest constraints to date on models of the evolution of cluster metallicity at intermediate radii. Adopting a power-law model of the form Z ∝ (1 + z)γ, we measure a slope $\gamma = -0.5^{+0.4}_{-0.3}$, consistent with the majority of the enrichment of the ICM having occurred at very early times and before massive clusters formed, but leaving open the possibility that some additional enrichment in these regions may have occurred since a redshift of 2.

79 ASTRONOMY AND ASTROPHYSICS↗