Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Unity”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Feedback in the dark: a critical examination of CMB bounds on primordial black holes

Abstract If present in the early universe, primordial black holes (PBHs) would have accreted matter and emitted high-energy photons, altering the statistical properties of the Cosmic Microwave Background (CMB). This mechanism has been used to constrain the fraction of dark matter that is in the form of PBHs to be much smaller than unity for PBH masses well above one solar mass. Moreover, the presence of dense dark matter mini-halos around the PBHs has been used to set even more stringent constraints, as these would boost the accretion rates.In this work, we critically revisit CMB constraints on PBHs taking into account the role of the local ionization of the gas around them. We discuss how the local increase in temperature around PBHs can prevent the dark matter mini-halos from strongly enhancing the accretion process, in some cases significantly weakening previously derived CMB constraints. We explore in detail the key ingredients of the CMB bound and derive a conservative limit on the cosmological abundance of massive PBHs.

Astronomy & Astrophysics↗
DOE OSTI2025

Neutrino and pair creation in reconnection-powered coronae of accreting black holes

A ubiquitous feature of accreting black hole systems is their hard X-ray emission which is thought to be produced through Comptonization of soft photons by electrons and positrons in the vicinity of the black hole, in a region with optical depth of order unity. The origin and composition of this Comptonizing region, known as the corona, is a matter open for debate. In this paper we investigate the role of relativistic protons accelerated in black-hole magnetospheric current sheets for the pair enrichment and neutrino emission of AGN coronae. Our model has two free parameters, namely the proton plasma magnetization σ p , which controls the peak energy of the neutrino spectrum, and the Eddington ratio λ X,Edd (defined as the ratio between X-ray luminosity L X and Eddington luminosity L Edd ), which controls the amount of energy transferred to secondary particles. For sources with λ X,Edd ≳ λ Edd,crit (where λ Edd,crit ∼ 10 -1 for σ p = 10 5 or ∼ 10 -2 for σ p = 10 7 ), proton-photon interactions and γγ annihilation produce enough secondary pairs to achieve Thomson optical depths τ T ∼ 0.1-10. In the opposite case of λ X,Edd ≲ λ Edd,crit , the coronal pairs cannot originate only from hadronic interactions. Additionally, we find that the neutrino luminosity scales as L 2 X /L Edd for λ X,Edd ≲ λ Edd,crit , while it is proportional to L X for higher λ X,Edd values. We apply our model to four Seyfert galaxies, including NGC 1068, and discuss our results in light of recent IceCube observations.

active galactic nuclei↗
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 OSTI2020

Magnetic shear effect on plasma transport at T e / T i ~ 1 through electron cyclotron heating in DIII-D plasmas

The effect of magnetic shear on plasma transport for an electron to ion temperature ratio ( T e / T i ) near unity has been explored in DIII-D utilizing electron cyclotron heating (ECH). Previous reports showed that significant confinement degradation occurred at T e / T i ~ 1 in positive shear (PS) plasmas in DIII-D, whereas reduced confinement degradation was observed in negative central shear (NCS) plasmas. In this study, plasma transport in weak magnetic shear (WS) plasmas with ECH is investigated and compared with that in NCS and PS plasmas. Here the magnetic shears ($\hat s$) are $\hat s$ > 0.5, ~0 and <-0.1 in the core region (ρ~ 0.3–0.4) of PS, WS and NCS plasmas, respectively, and flat or negative inside ρ~ 0.4 in the WS and NCS plasmas. Weak magnetic shear is found to be effective in minimizing degradation of ion thermal confinement as T e / T i increases through ECH application, and an improved confinement factor of H 98y2 ~ 1.2 is maintained, similar to NCS plasmas. At T e / T i ~ 1, the ion thermal diffusivity around an internal transport barrier decreases when changing the magnetic shear from positive to weak or negative shear. Also, reduced local particle and momentum transport was indicated by steeper density and toroidal rotation profiles in the weak and negative shear regimes. Linear gyrokinetic simulations predict little change in growth rates of low-k turbulence with ECH application in the WS and NCS plasmas, which is consistent with the transport and profile analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Plasma radiation behavior approaching high-radiation scenarios in W7-X

The W7-X stellarator has so far performed experiments under both limiter and divertor conditions. The plasma is mostly generated by ECR-heating with powers up to 6.5 MW, and the plasma density is usually limited by the radiation losses from low-Z impurities (such as carbon and oxygen) released mainly from the graphite targets. The present work first summarizes the radiation loss fractions f rad achieved in quasi-stationary hydrogen plasmas in both operational phases, and then shows how impurity radiation behaves differently with the two different boundary conditions as the plasma density increases. The divertor operation is emphasized and some beneficial effects (with respect to impurity radiation) are highlighted: (1) intensive radiation is located at the edge (r/a > 0.8) even at high radiation loss fractions, (2) the plasma remains stable up to f rad approaching unity, (3) the reduction in the stored energy is about 10% for high f rad scenarios. Moreover, effects of wall boronisation on impurity radiation profiles are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Tungsten–carbon surface evolution and erosion modeling for a small angle slot divertor in DIII-D

Abstract We modeled tungsten–carbon mixed surface evolution, sputtering erosion, and transport for the tungsten coated region of a small angle slot (SAS) divertor design for the DIII-D tokamak. This divertor concept aims to achieve a closed slot dissipative plasma to minimize heat load and surface erosion, and to study high- Z material performance. Our advanced simulations use coupled ITMC-DYN material mixing/response and 3D full kinetic REDEP/WBC erosion/redeposition code packages, with divertor plasma solution from the SOLPS-ITER package with 4 MW power input. The SAS design geometry and resulting in-slot plasma parameters cause significant differences in predicted sputter/transport from a conventional divertor. For 2% C/D incident plasma ratio, an equilibrium mixed C/W surface is attained at ∼30 s of discharge, from wall sputtered carbon transported to the 10 cm long tungsten divertor region. Tungsten remains exposed to the plasma, but the evolved surface composition varies with different C/D ratios. Tungsten is primarily sputtered from the mixed surface by impinging carbon ions in the +1 to +4 charge states, with some self-sputtering. Redeposition of sputtered tungsten to the divertor is significant, ∼80% along the higher plasma temperature attached plasma SAS entrance region, but this is less than the typically near-unity values for a conventional divertor. Plasma-incident carbon is highly backscattered (∼50%) from the mixed surface, with little redeposition (<10%); this helps maintain tungsten in the surface sputter zone. Carbon is mainly sputtered from the mixed surface by D + ions, also with low redeposition (∼10%–30%). Finally, the modeling shows non-zero but low sputtered tungsten current from the divertor to the core plasma direction. These results appear favorable for effective testing of a tungsten-containing SAS divertor in DIII-D, and extrapolation of mixed-material evolution/response findings to the analogous low- Z /high- Z , Be/W, ITER plasma facing system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modelling of sawtooth-induced fast ion transport in positive and negative triangularity in TCV

Internal kinks are a common magneto hydro-dynamic (MHD) instability observed in tokamak operation when the q profile in the plasma core is close to unity. This MHD instability impacts both the transport of the bulk plasma (current, particle and energy transport) and minority species, such as fast ions. In TCV (R 0 /a = 0.88 m/0.25 m) the fast ion population is generated in the plasma by neutral beam tangential injection of energies up to 28 keV. TCV features 16 active shaping coils permitting a great flexibility in plasma shape, including negative triangularity (δ) configurations that show surprisingly high confinement. This study focuses on the transport of fast ions induced by sawteeth, by comparing two triangularity cases and simulation results with experimental data. Comparison of two equilibria with opposite δ shows that the fast ion drifts are larger for δ < 0. Furthermore, the sawtooth-induced transport in this case is larger than δ > 0 in similar conditions. Additionally, comparison with experimental data confirms the dominance of the modification of thermal kinetic profiles following the sawtooth crash in explaining drops in the neutron rates and FIDA signals. Additional fast ion diffusion, however, improves the interpretation of the experimental data. For δ < 0, the amplitude of the perturbation better representing the experimental data is larger. Finaly, an exploratory study for 50 keV particles (soon available in TCV) shows that the situation does not worsen for such particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analysis of hydrogen fueling, recycling, and confinement at Wendelstein 7-X via a single-reservoir particle balance

Abstract A single-reservoir particle balance for the main plasma species hydrogen has been established for W7-X. This has enabled the quantitative characterization of the particle sources in the standard island divertor configuration for the first time. Findings from attached scenarios with two different island sizes with a boronized wall and turbo molecular pumping are presented. Fueling efficiencies, particle flows and source locations were measured and used to infer the total particle confinement time τ p . Perturbative gas injection experiments served to measure the effective particle confinement time τ p * . Combining both confinement times provides access to the global recycling coefficient R ¯ . Hydrogen particle inventories have been addressed and the knowledge of particle sources and sinks reveals the core fueling distribution and provides insight into the capability of the magnetic islands to control exhaust features. Measurements of hydrogen fueling efficiencies were sensitive to the precise fueling location and measured between 12% and 31% with the recycling fueling at the strike line modeled at only 6%, due to much higher densities. 15% of the total 5.2 × 10 22 a/s recycling flow ionizes far away from the recycling surfaces in the main chamber. It was shown that 60% of recycled particles ionize above the horizontal and 18% above the vertical divertor target, while the remainder of the recycling flow ionizes above the baffle (7%). Combining these source terms with their individual fueling efficiencies resolves the core fueling distribution. Due to the higher fueling efficiency in the main chamber, up to 51% of the total 5.1 × 10 21 s −1 core fueling particles are entering the confined plasma from the main chamber. τ p values in the range of 260 ms were extracted for these discharges. Together with τ p , the global recycling coefficient R ¯ was resolved for every τ p * measurement and a typical value close to unity was obtained. An increase of the island size, resulted in no change of τ p , but doubled τ p * , indicating the feasibility of the control coils as an actuator to control exhaust features without affecting core confinement properties.

recycling↗
DOE OSTI2024

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Multi-X-point radiation and its dynamics in up/down asymmetry in plasma detachment regimes in Wendelstein 7-X

In Wendelstein 7-X (W7-X), stable plasma detachment in ECRH hydrogen plasmas is routinely achieved with the standard magnetic configuration, containing five magnetic islands in the scrape-off layer (SOL). Plasma detachment induced by intrinsic carbon impurities is characterized by a high radiation fraction (ƒ rad ) with significantly reduced divertor heat load and particle flux. In this work, we present the radiation dynamics during the detached plasma (DP) phase, such as the shift of radiation zones in radial and poloidal directions towards the X-points as ƒ rad increases (∼0.6–0.9), as well as the penetration of the radiation layer into the confinement region in the deep DP phase with ƒ rad approaching unity. In particular, a structure of multi-X-point radiation (multi-XPR) with an up/down asymmetry in the DP phase is highlighted, which is revealed by bolometer tomography and is further confirmed by video diagnostics. The multi-XPR structure forms helical 3D bands aligned with W7-X field periodicity. Field reversal experiments show that the brightest XPR displaces between the upper and lower SOL regions, suggesting the E× B drift effects. This paper presents the multi-XPR structure observed in the W7-X plasma in detail for the first time. A simplified model considering the influence of the poloidal E× B drift ( V d ) on the impurity flow in the SOL shows that the downstream drift toward the target or target shadow region ( V d < 0) decreases the impurity content, while the upstream drift toward the LCFS ( V d > 0) increases the impurity content. On this basis, the poloidal drift potentially leads to an up/down asymmetry of impurity density in the SOL despite the symmetry magnetic topology. The dynamics of the up/down asymmetry in the multi-XPR structure is also related to the magnitude V d /D (normalized to the impurity diffusivity), with an additional effect owing to the radial inward shift of the emission zone. These results provide new insights into impurity-induced detachment dynamics and provide a basis for improving 3D modeling of impurity transport, which typically does not consider drifts.

E × B drift↗
DOE OSTI2024

A numerical Poisson solver with improved radial solutions for a self-consistent locally scaled self-interaction correction method

Abstract The universal applicability of density functional approximations is limited by self-interaction error made by these functionals. Recently, a novel one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities was proposed. The locally scaled SIC (LSIC) method is exact for the one-electron densities, and unlike the well-known Perdew–Zunger SIC (PZSIC) method recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to PZSIC method as a special case when isoorbital indicator is set to the unity. Here, we present a numerical scheme that we have adopted to evaluate the Coulomb potential of the electron density scaled by the iso-orbital indicator required for the self-consistent LSIC calculations. After analyzing the behavior of the finite difference method (FDM) and the green function solution to the radial part of the Poisson equation, we adopt a hybrid approach that uses the FDM for the Coulomb potential due to the monopole and the GF for all higher-order terms. The performance of the resultant hybrid method is assessed using a variety of systems. The results show improved accuracy than earlier numerical schemes. We also find that, even with a generic set of radial grid parameters, accurate energy differences can be obtained using a numerical Coulomb solver in standard density functional studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Modular categories of dimension $p^3m$ with $m$ square-free

We give a complete classification of modular categories of dimension $p^3m$ where $p$ is prime and $m$ is a square-free integer. When $p$ is odd, all such categories are pointed. For $p$ = 2 one encounters modular categories with the same fusion ring as orthogonal quantum groups at certain roots of unity, namely SO (2$m$) 2 . As a result, we also classify the more general class of modular categories with the same fusion rules as SO (2$N$) 2 with $N$ odd.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Gravito-turbulence and dynamo in poorly ionized protostellar discs – I. Zero-net-flux case

In their early stages, protoplanetary discs are sufficiently massive to undergo gravitational instability (GI). This instability is thought to be involved in mass accretion, planet formation via gas fragmentation, the generation of spiral density waves, and outbursts. A key and very recent area of research is the interaction between the GI and magnetic fields in young protoplanetary discs, in particular whether this instability is able to sustain a magnetic field via a dynamo. We conduct 3D, stratified shearing-box simulations using two independent codes, PLUTO and Athena++, to characterize the GI dynamo in poorly ionized protostellar discs subject to ambipolar diffusion. We find that the dynamo operates across a large range of ambipolar Elssaser number Am (which characterizes the strength of ambipolar diffusion) and is particularly strong in the regime Am = 10–100, with typical magnetic to thermal energy ratios of order unity. The dynamo is only weakly dependent on resolution (at least for Am ≲ 100), box size, and cooling law. The magnetic field is produced by the combination of differential rotation and large-scale vertical roll motions associated with spiral density waves. Our results have direct implications for the dynamo process in young protoplanetary discs and possibly some regions of active galactic nucleus discs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Ionizing radiation from AGNs at z > 3.3 with the Subaru Hyper Suprime-Cam Survey and the CFHT Large Area U-band Deep Survey (CLAUDS)

ABSTRACT We use deep and wide imaging data from the CFHT Large Area U-band Deep Survey (CLAUDS) and the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to constrain the ionizing radiation (Lyman Continuum, LyC) escape fraction from active galactic nuclei (AGNs) at z ∼ 3–4. For 94 AGNs with spectroscopic redshifts at 3.3 < z < 4.0, we use their U/i-band flux ratios to estimate LyC transmission of individual AGNs. The distribution of their LyC transmission shows values lower than the range of LyC transmission values for IGM of the same redshift range, which suggests that LyC escape fraction of AGNs at z > 3.3 is considerably lower than unity in most cases. We do not find any trend in LyC transmission values depending on their UV luminosities. Based on the photometry of stacked images we find the average flux ratio of LyC and non-ionizing UV photons escaping from the objects (fLyC/fUV)out = 0.182 ± 0.043 for AGNs at 3.3 < z < 3.6, which corresponds to LyC escape fraction fesc = 0.303 ± 0.072 if we assume a fiducial intrinsic SED of AGN. Based on the estimated LyC escape fraction and the UV luminosity function (UVLF) of AGNs, we argue that UV-selected AGNs’ contribution to the LyC emissivity at the epoch is minor, although the size of their contribution largely depends on the shape of the UVLF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Estimating the local dark matter density in a non-axisymmetric wobbling disc

The density of dark matter near the Sun, ρDM, ⊙, is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. Here, we assess the impact of these assumptions on determinations of ρDM, ⊙ by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia–Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρDM, ⊙ from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2 kpc above the disc plane, implying a fractional error on ρDM, ⊙ of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρDM, ⊙, consistent with the true value within our 95 percent confidence intervals for realistic 20 percent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 percent to ρDM, ⊙ (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The half-mass radius of MaNGA galaxies: effect of IMF gradients

Gradients in the stellar populations (SP) of galaxies – e.g. in age, metallicity, stellar initial mass function (IMF) – can result in gradients in the stellar-mass-to-light ratio, M * / L . Such gradients imply that the distribution of the stellar mass and light is different. For old SPs, e.g. in early-type galaxies at z ∼ 0, the M * / L gradients are weak if driven by variations in age and metallicity, but significantly larger if driven by the IMF. A gradient which has larger M * / L in the centre increases the estimated total stellar mass ( M * ) and reduces the scale which contains half this mass ( R e,* ), compared to when the gradient is ignored. For the IMF gradients inferred from fitting MILES simple SP models to the H β, 〈Fe〉, [MgFe], and TiO 2SDSS absorption lines measured in spatially resolved spectra of early-type galaxies in the MaNGA survey, the fractional change in R e,* can be significantly larger than that in M * , especially when the light is more centrally concentrated. The R e,* – M * correlation which results from accounting for IMF gradients is offset to smaller sizes by 0.3 dex compared to when these gradients are ignored. Comparisons with ‘quiescent’ galaxies at higher z must account for evolution in SP gradients (especially age and IMF) and in the light profile before drawing conclusions about how R e,* and M * evolve. Here, the implied merging between higher z and the present is less contrived if R e,* / R e at z ∼ 0 is closer to our IMF-driven gradient calibration than to unity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Magnetic fields in the formation of the first stars – II. Results

ABSTRACT Beginning with cosmological initial conditions at z = 100, we simulate the effects of magnetic fields on the formation of Population III stars and compare our results with the predictions of Paper I. We use gadget-2 to follow the evolution of the system while the field is weak. We introduce a new method for treating kinematic fields by tracking the evolution of the deformation tensor. The growth rate in this stage of the simulation is lower than expected for diffuse astrophysical plasmas, which have a very low resistivity (high magnetic Prandtl number); we attribute this to the large numerical resistivity in simulations, corresponding to a magnetic Prandtl number of order unity. When the magnetic field begins to be dynamically significant in the core of the minihalo at z = 27, we map it on to a uniform grid and follow the evolution in an adaptive mesh refinement, MHD simulation in orion2. The non-linear evolution of the field in the orion2 simulation violates flux-freezing and is consistent with the theory proposed by Xu & Lazarian. The fields approach equipartition with kinetic energy at densities ∼1010–1012 cm−3. When the same calculation is carried out in orion2 with no magnetic fields, several protostars form, ranging in mass from ∼1 to 30 M⊙; with magnetic fields, only a single ∼30 M⊙ protostar forms by the end of the simulation. Magnetic fields thus suppress the formation of low-mass Pop III stars, yielding a top-heavy Pop III initial mass function and contributing to the absence of observed Pop III stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

SILCC – VII. Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities

ABSTRACT We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σgas = 10, 30, 50, and 100 $\mathrm{M_\odot \, pc^{-2}}$. The silcc project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density ΣSFR. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $\sigma _\mathrm{z} \propto \Sigma _\mathrm{SFR}^a$, with a ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼2.2 higher than that of the CNM, in agreement with local observations. For $\Sigma _\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\, \mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$ the WIM motions become supersonic.

Astronomy & Astrophysics↗