Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STAR EVOLUTION”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

How galaxies populate halos in extreme density environments: An analysis of the halo occupation distribution in SDSS

Recent works have shown that the properties of galaxy populations in simulated dark matter halos vary with large-scale environments. These results suggest a variation in the halo occupation distribution (HOD) in extreme density environments since the dynamical and astrophysical conditions prevailing in these regions may significantly affect the formation and evolution of their halos and residing galaxies, influencing the mean number of satellite galaxies. To analyse these effects, we identify cosmic voids and future virialised structures (FVS) in the Sloan Digital Sky Server Data Release 12 (SDSS-DR12) and estimate the HOD within these super-structures using group catalogues as dark matter halo proxies. Our goal is to use observational galaxy data to characterise the HOD within cosmic voids and FVS, explore the different properties of these galaxies populations and compare them with the general results outside these super-structures. We use a publicly available observational galaxy catalogue with information on redshifts, positions, magnitudes and other astrophysical features to build a volume complete galaxy sample and identify cosmic voids and FVS. Using a publicly available galaxy group catalogue as a proxy to dark matter halos, we compute the HOD within both types of super-structures for different absolute magnitude thresholds. We also study the dependence on the results on the main void and FVS properties, density and volume. We also analysed the main characteristics of the stellar content of galaxies inside these extreme-density regions such as the mean stellar age and the galaxy light concentration index. In all cases, we compare the results with those derived from the Field sample, defined by objects outside both types of environments. Inside cosmic voids, we find a strong decrease in HOD concerning the Field results. In the most extreme cases, the mean number of satellites fall to ~ 50%. Inside FVS, the HOD shows a significant increase to the Field, with a ~ 40% excess in the mean number of satellites. These results are present for the different galaxy luminosity ranges explored. In both environments, the differences with respect to the Field increases for the extreme values of the density environments. However, we obtain no signs of variations related to intrinsic characteristics of the super-structures, indicating that the effects depend mainly on the density of the large-scale environment. In addition, we find that the cumulative distribution of the mean age of stars of the central galaxy also varies in the different regions, this suggests that the history of the formation of the dark matter halos may be different. Finally, we explore the HOD for the 25% youngest (oldest) galaxies, based on the mean age of their stars. We find that for the low-mass groups the youngest galaxies are only present inside voids and are generally central galaxies. On the other hand, for the high-mass groups the FVS environments show the same increase in the HOD concerning the Field as previously mentioned. We find that cosmic voids lack a significant fraction of galaxies with the oldest stellar population.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Observing the onset of pressure-driven K-shell delocalization

The gravitational pressure in many astrophysical objects exceeds one gigabar (one billion atmospheres), creating extreme conditions where the distance between nuclei approaches the size of the K shell. This close proximity modifies these tightly bound states and, above a certain pressure, drives them into a delocalized state. Both processes substantially affect the equation of state and radiation transport and, therefore, the structure and evolution of these objects. Still, our understanding of this transition is far from satisfactory and experimental data are sparse. Here, in this work, we report on experiments that create and diagnose matter at pressures exceeding three gigabars at the National Ignition Facility where 184 laser beams imploded a beryllium shell. Bright X-ray flashes enable precision radiography and X-ray Thomson scattering that reveal both the macroscopic conditions and the microscopic states. The data show clear signs of quantum-degenerate electrons in states reaching 30 times compression, and a temperature of around two million kelvins. At the most extreme conditions, we observe strongly reduced elastic scattering, which mainly originates from K-shell electrons. We attribute this reduction to the onset of delocalization of the remaining K-shell electron. With this interpretation, the ion charge inferred from the scattering data agrees well with ab initio simulations, but it is significantly higher than widely used analytical models predict.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Orbital Decay of Double White Dwarfs: Beyond Gravitational-wave Radiation Effects

The traditional description of the orbital evolution of compact-object binaries, like double white dwarfs (DWDs), assumes that the system is driven only by gravitational-wave (GW) radiation. However, the high magnetic fields with intensities of up to gigagausses measured in WDs alert a potential role of the electromagnetic (EM) emission in the evolution of DWDs. We evaluate the orbital dynamics of DWDs under the effects of GW radiation, tidal synchronization, and EM emission by a unipolar inductor generated by the magnetic primary and the relative motion of the nonmagnetic secondary. We show that the EM emission can affect the orbital dynamics for magnetic fields larger than megagausses. We applied the model to two known DWDs, SDSS J0651+2844 and ZTF J1539+5027, for which the GW radiation alone does not fully account for the measured orbital decay rate. We obtain upper limits to the primary's magnetic field strength, over which the EM emission causes an orbital decay faster than observed. The contribution of tidal locking and the EM emission is comparable, and together they can contribute up to 20% to the measured orbital decay rate. We show that the gravitational waveform for a DWD modeled as purely driven by GWs and including tidal interactions and EM emission can have large relative dephasing detectable in the mHz regime of frequencies relevant for space-based detectors like LISA. Therefore, including physics besides GW radiation in the waveform templates is essential to calibrate the GW detectors using known sources, e.g., ZTF J1539+5027, and to infer binary parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Revision of Faraday rotation measure constraints on the primordial magnetic field using the IllustrisTNG simulation

ABSTRACT Previously derived Faraday rotation constraints on the volume-filling intergalactic magnetic field (IGMF) have used analytical models that made a range of simplifying assumptions about magnetic field evolution in the intergalactic medium and did not consider the effect of baryonic feedback on large-scale structures. In this work, we revise existing Faraday rotation constraints on the IGMF using a numerical model of the intergalactic medium from the IllustrisTNG cosmological simulation that includes a sophisticated model of the baryonic feedback. We use the IllustrisTNG model to calculate the rotation measure and compare the resulting mean and median of the absolute value of the rotation measure with data from the NRAO VLA Sky Survey (NVSS). The numerical model of the intergalactic medium includes a full magnetohydrodynamic model of the compressed primordial magnetic field as well as a model of the regions where the magnetic field is not primordial, but is rather produced by the process of baryonic feedback. Separating these two types of regions, we are able to assess the influence of the primordial magnetic field on the Faraday rotation signal. We find that by correcting for regions of compressed primordial field and accounting for the fact that part of the intergalactic medium is occupied by magnetic fields spread by baryonic feedback processes rather than by the primordial field relaxes the Faraday rotation bound by a factor of ≃3. This results in B0 < 1.8 × 10−9 G for large correlation length IGMFs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Galactic echoes

Gaia has revealed a variety of substructures in the phase space of stars in the Solar neighbourhood, including the vertical ‘Snail’ $(z, v_z)$ in space. Such substructures are often interpreted as the incompletely phase-mixed response of the disc stars to a single perturbation, such as an impulsive encounter with a satellite galaxy. In this paper, we consider the possibility that such structures contain manifestations of phase-space echoes. First established in plasma physics in the 1960s, echoes arise when a collisionless system is perturbed twice: the macroscopic responses to both perturbations mix to small scales in phase space, whereupon they couple non-linearly, producing a third macroscopic ‘echo’ response without the need for a third perturbation. We derive the galactic analogue of the plasma echo theory using angle-action variables and apply it to a one-dimensional model of vertical motion in the Milky Way. We verify the predicted echo behaviour using idealized test particle simulations, both with and without the inclusion of diffusion through orbital scattering off molecular clouds. While we conclude that the Gaia Snail itself is unlikely a (pure) echo effect, the basic physics we uncover is sufficiently generic that we expect phase-space echoes to be common in disc galaxies.

Galaxy: evolution↗
DOE OSTI2022

The High Fraction of Thin Disk Galaxies Continues to Challenge ΛCDM Cosmology

Any viable cosmological framework has to match the observed proportion of early- and late-type galaxies. In this contribution, we focus on the distribution of galaxy morphological types in the standard model of cosmology (Lambda cold dark matter, ΛCDM). Using the latest state-of-the-art cosmological ΛCDM simulations known as Illustris, IllustrisTNG, and EAGLE, we calculate the intrinsic and sky-projected aspect ratio distribution of the stars in subhalos with stellar mass M * > 10 10 M ⊙ at redshift z = 0. There is a significant deficit of intrinsically thin disk galaxies, which however comprise most of the locally observed galaxy population. Consequently, the sky-projected aspect ratio distribution produced by these ΛCDM simulations disagrees with the Galaxy And Mass Assembly (GAMA) survey and Sloan Digital Sky Survey at ≥12.52σ (TNG50-1) and ≥14.82σ (EAGLE50) confidence. The deficit of intrinsically thin galaxies could be due to a much less hierarchical merger-driven build-up of observed galaxies than is given by the ΛCDM framework. It might also arise from the implemented sub-grid models, or from the limited resolution of the above-mentioned hydrodynamical simulations. We estimate that an 8 5 times better mass resolution realization than TNG50-1 would reduce the tension with GAMA to the 5.58σ level. Finally, we show that galaxies with fewer major mergers have a somewhat thinner aspect ratio distribution. Given also the high expected frequency of minor mergers in ΛCDM, the problem may be due to minor mergers. In this case, the angular momentum problem could be alleviated in Milgromian dynamics because of a reduced merger frequency arising from the absence of dynamical friction between extended dark matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The incidence of X-ray selected AGN in nearby galaxies

ABSTRACT We present the identification and analysis of an unbiased sample of active galactic nuclei (AGN) that lie within the local galaxy population. Using the MPA-JHU catalogue (based on SDSS DR8) and 3XMM DR7 we define a parent sample of 25 949 local galaxies (z ≤ 0.33). After confirming that there was strictly no AGN light contaminating stellar mass and star-formation rate calculations, we identified 917 galaxies with central, excess X-ray emission likely originating from an AGN. We analysed their optical emission lines using the BPT diagnostic and confirmed that such techniques are more effective at reliably identifying sources as AGN in higher mass galaxies: rising from 30 per cent agreement in the lowest mass bin to 93 per cent in the highest. We then calculated the growth rates of the black holes powering these AGN in terms of their specific accretion rates (∝LX/M*). Our sample exhibits a wide range of accretion rates, with the majority accreting at rates $\le 0.5\ \mathrm{ per \, cent}$ of their Eddington luminosity. Finally, we used our sample to calculate the incidence of AGN as a function of stellar mass and redshift. After correcting for the varying sensitivity of 3XMM, we split the galaxy sample by stellar mass and redshift and investigated the AGN fraction as a function of X-ray luminosity and specific black hole accretion rate. From this we found the fraction of galaxies hosting AGN above a fixed specific accretion rate limit of 10−3.5 is constant (at $\approx 1\ \mathrm{ per \, cent}$) over stellar masses of 8 < log M*/M⊙ < 12 and increases (from $\approx 1\ \mathrm{ per \, cent}$ to 10 per cent) with redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Oxidation of the Interiors of Carbide Exoplanets

Astrophysical measurements have shown that some stars have sufficiently high carbon-to-oxygen ratios such that the planets they host would be mainly composed of carbides instead of silicates. We studied the behavior of silicon carbide in the presence of water under the high pressure–temperature conditions relevant to planetary interiors in the laser-heated diamond-anvil cell. When reacting with water, silicon carbide converts to silica (stishovite) and diamond at pressures up to 50 GPa and temperatures up to 2500 K: $\mathrm{SiC}+2{{\rm{H}}}_{2}{\rm{O}}\to {\mathrm{SiO}}_{2}+{\rm{C}}+2{{\rm{H}}}_{2}$. Therefore, if water can be incorporated into carbide planets during their formation or through later delivery, they could be oxidized and have mineralogy dominated by silicates and diamond in their interiors. The reaction could produce CH 4 at shallower depths and H 2 at greater depths that could be degassed from the interior, causing the atmospheres of the converted carbon planets to be rich in reducing gases. Excess water after the reaction can be stored in dense silica polymorphs in the interiors of the converted carbon planets. Such conversion of mineralogy to diamond and silicates would decrease the density of carbon-rich planet, making the converted planets distinct from silicate planets in mass–radius relations for the 2–8 Earth mass range.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Host galaxies of ultrastrong Mg ii absorbers at z ∼ 0.5

ABSTRACT From a sample of 109 candidate ultrastrong Mg ii (USMg ii; having rest equivalent width of Mg ii, W2796 > 3.0 Å) systems at z = 0.4–0.6, we confirm 27 and identify host galaxies of 20 systems based on associated nebular line emission from our SALT observations or from Sloan Digital Sky Survey (SDSS) fiber spectra. The measured impact parameter, [O ii] luminosity, star formation rate, B-band luminosity, and stellar mass are in the ranges 7.3 ≤ D[kpc] ≤ 79, $0.2\le L_{[\mathrm{ O}\,~\small {\rm II}]}[ 10^{41}~\mathrm{ erg} \mathrm{ s}^{-1}]\le 4.5$, 2.59 ≤ SFR[M⊙yr−1] ≤ 33.51, $0.15L_B^{*}\le L_B\le 1.63L_B^{*}$, and 10.21 ≤ log[M*/ M ⊙] ≤ 11.62, respectively. The impact parameters found are larger than that predicted by the W2796 versus D relationship of the general population of Mg ii absorbers. At a given D, USMg ii host galaxies are more luminous and massive compared to typical Mg ii absorbers. However, the measured SFRs are slightly lower than that of main-sequence galaxies with same M⋆ at z ∼ 0.5. We report a correlation between $L_{[\mathrm{ O}\,\small {\rm II}]}$ and W2796 for the full population of Mg ii absorbers, driven mainly by the host galaxies of weak Mg ii absorbers that tend to have low $L_{[\mathrm{ O}\,\small {\rm II}]}$ and large impact parameters. We find at least ∼33 per cent of the USMg ii host galaxies (with a limiting magnitude of mr < 23.6) are isolated and the large W2796 in these cases may originate from gas flows (infall/outflow) in single haloes of massive but not starburst galaxies. We also find galaxy interactions could be responsible for large velocity widths in at least ∼17 per cent cases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

3D Hydrodynamics of Pre-supernova Outbursts in Convective Red Supergiant Envelopes

Eruptive mass loss likely produces the energetic outbursts observed from some massive stars before they become core-collapse supernovae (SNe). The resulting dense circumstellar medium may also cause the subsequent SNe to be observed as Type IIn events. The leading hypothesis of the cause of these outbursts is the response of the envelope of the red supergiant (RSG) progenitor to energy deposition in the months to years prior to collapse. Early theoretical studies of this phenomenon were limited to 1D, leaving the 3D convective RSG structure unaddressed. Using FLASH's hydrodynamic capabilities, we explore the 3D outcomes by constructing convective RSG envelope models and depositing energies less than the envelope binding energies on timescales shorter than the envelope dynamical time deep within them. We confirm the 1D prediction of an outward-moving acoustic pulse steepening into a shock, unbinding the outermost parts of the envelope. However, we find that the initial 2–4 km s -1 convective motions seed the intrinsic convective instability associated with the high-entropy material deep in the envelope, enabling gas from deep within the envelope to escape and increasing the amount of ejected mass compared to an initially "quiescent" envelope. The 3D models reveal a rich density structure, with column densities varying by ≈10× along different lines of sight. Our work highlights that the 3D convective nature of RSG envelopes impacts our ability to reliably predict the outburst dynamics, the amount, and the spatial distribution of the ejected mass associated with deep energy deposition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Stellar mass as a galaxy cluster mass proxy: application to the Dark Energy Survey redMaPPer clusters

Abstract We introduce a galaxy cluster mass observable, μ⋆, based on the stellar masses of cluster members, and we present results for the Dark Energy Survey (DES) Year 1 (Y1) observations. Stellar masses are computed using a Bayesian model averaging method, and are validated for DES data using simulations and COSMOS data. We show that μ⋆ works as a promising mass proxy by comparing our predictions to X-ray measurements. We measure the X-ray temperature–μ⋆ relation for a total of 129 clusters matched between the wide-field DES Y1 redMaPPer catalogue and Chandra and XMM archival observations, spanning the redshift range 0.1 < $z$ < 0.7. For a scaling relation that is linear in logarithmic space, we find a slope of α = 0.488 ± 0.043 and a scatter in the X-ray temperature at fixed μ⋆ of $\sigma _{{\rm ln} T_\mathrm{ X}|\mu _\star }= 0.266^{+0.019}_{-0.020}$ for the joint sample. By using the halo mass scaling relations of the X-ray temperature from the Weighing the Giants program, we further derive the μ⋆-conditioned scatter in mass, finding $\sigma _{{\rm ln} M|\mu _\star }= 0.26^{+ 0.15}_{- 0.10}$. These results are competitive with well-established cluster mass proxies used for cosmological analyses, showing that μ⋆ can be used as a reliable and physically motivated mass proxy to derive cosmological constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Oxygen isotope ratios in zircon and garnet: A record of assimilation and fractional crystallization in the Dinkey Dome peraluminous granite, Sierra Nevada, California

Abstract The 119 Ma Dinkey Dome pluton in the central Sierra Nevada Batholith is a peraluminous granite and contains magmatic garnet and zircon that are complexly zoned with respect to oxygen isotope ratios. Intracrystalline SIMS analysis tests the relative importance of magmatic differentiation processes vs. partial melting of metasedimentary rocks. Whereas δ18O values of bulk zircon concentrates are uniform across the entire pluton (7.7‰ VSMOW), zircon crystals are zoned in δ18O by up to 1.8‰, and when compared to late garnet, show evidence of changing magma chemistry during multiple interactions of the magma with wall rock during crustal transit. The evolution from an early high-δ18O magma [δ18O(WR) = 9.8‰] toward lower values is shown by high-δ18O zircon cores (7.8‰) and lower δ 18O rims (6.8‰). Garnets from the northwest side of the pluton show a final increase in δ18O with rims reaching 8.1‰. In situ REE measurements show zircon is magmatic and grew before garnets. Additionally, δ18O in garnets from the western side of the pluton are consistently higher (avg = 7.3‰) relative to the west (avg = 5.9‰). These δ18O variations in zircon and garnet record different stages of assimilation and fractional crystallization whereby an initially high-δ18O magma partially melted low-δ18O wallrock and was subsequently contaminated near the current level of emplacement by higher δ18O melts. Collectively, the comparison of δ18O zoning in garnet and zircon shows how a peraluminous pluton can be constructed from multiple batches of variably contaminated melts, especially in early stages of arc magmatism where magmas encounter significant heterogeneity of wall-rock assemblages. Collectively, peraluminous magmas in the Sierran arc are limited to small <100 km2 plutons that are intimately associated with metasedimentary wall rocks and often surrounded by later and larger metaluminous tonalite and granodiorite plutons. The general associations suggest that early-stage arc magmas sample crustal heterogeneities in small melt batches, but that with progressive invigoration of the arc, such compositions are more effectively blended with mantle melts in source regions. Thus, peraluminous magmas provide important details of the nascent Sierran arc and pre-batholithic crustal structure.

Geochemistry & Geophysics↗
DOE OSTI2026

Mathematical Modeling of the Potential and Time Dependence of Ir Dissolution from Hydrous Ir Oxide Oxygen Evolution Catalysts

One of the main degradation mechanisms of hydrous iridium oxide acidic oxygen evolution reaction (OER) catalysts is dissolution and loss into the acidic membrane. While degradation models have been proposed, there is a gap in understanding the potential and time dependence of the iridium dissolution reaction and its mechanistic underpinnings. In this work, Ir dissolution rates measured as a function of time and potential via time-resolved inductively-coupled plasma mass spectrometry (ICP-MS) in aqueous acidic electrolyte are used to establish a mathematical model for Ir dissolution. The mathematical model is generated using proposed formation and dissolution reactions for Ir species. Through comparison with the ICP-MS data and existing information on the potential-dependent Ir phase, we find that the potential and time-dependence of dissolution can be modeled as dissolution of an oxide phase, here represented as IrO 2 , with potential dependent kinetics and formation of a passivating species, a process with a rate-limiting step that is not potential dependent. This understanding of the potential dependence of dissolution and passivation kinetics using aqueous electrolyte half-cell measurements can be used to predict the degradation of Ir oxide in operating energy conversion devices relying on the OER, such as proton-exchange membrane water electrolyzers.

Kariuki, Nancy N. [Argonne National Laboratory (AN↗
DOE OSTI2020

Low-energy cross-section measurement of the B 10 ( α , n ) N 13 reaction and its impact on neutron production in first-generation stars

Nucleosynthesis in the first generation of massive stars offers a unique setting to explore the creation of the first heavier nuclei in an environment free of impurities from earlier stellar generations. In later generations of massive stars, hydrogen burning occurs predominantly through the CNO cycles, but without the carbon, nitrogen, and oxygen to catalyze the reaction sequence, first stars would have to rely on the inefficient pp chains for their energy production. Observations of second and third generation stars show pronounced abundances of carbon and oxygen isotopes, which suggests a rapid conversion of the primordial abundances to heavier elements. While the triple-alpha-process primarily facilitates this conversion, there are alternative reaction sequences, such as 2 H(α,γ) 6 Li(α,γ) 10 B(α,n) 13 N, that may play a significant role. To study such alternate reaction pathways for production of carbon and heavier nuclei, a number of new measurements are needed. In this work, new measurements are reported for the 10 B(α,n) 13 N reaction, extending the cross section down to 575 keV incident α-particle energy. The measurements were made using a state-of-the-art deuterated liquid scintillator and a spectrum unfolding technique. An R-matrix analysis was performed in order to facilitate a comparison of the underlying nuclear structure with the reaction measurements. An unexpected upturn is observed in the low-energy S factor that indicates the presence of a new low-energy resonance. Finally, a revised reaction rate is determined that takes into account the present data as well as other previous measurements from the literature that were previously neglected.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution

ABSTRACT We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\, \mathrm{M_\odot }$ red giant branch primary and a $1\, \mathrm{M_\odot }$ companion modelled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, ∼15 per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of ${\sim}2\, {\rm yr}$, after which the envelope would be unbound.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Formation and evolution of protostellar accretion discs – I. Angular-momentum budget, gravitational self-regulation, and numerical convergence

Here, we investigate the formation and early evolution of a protostellar disc from a magnetized prestellar core using non-ideal magnetohydrodynamic (MHD) simulations including ambipolar diffusion and Ohmic dissipation. The dynamical contraction of the prestellar core ultimately leads to the formation of a first hydrostatic core, after ambipolar diffusion decouples the magnetic field from the predominantly neutral gas. The hydrostatic core accumulates angular momentum from the infalling material, evolving into a rotationally supported torus; this ‘first hydrostatic torus’ then forms an accreting protostar and a rotationally supported disc. The disc spreads out by gravitational instability, reaching ~30 au in diameter at ~3 kyr after protostar formation. The total mass and angular momentum of the protostar–disc system are determined mainly by accretion of gas from an infalling pseudo-disc, which has low specific angular momentum because of magnetic braking; their removal from the protostar–disc system by outflow and disc magnetic braking are negligible, in part because the magnetic field is poorly coupled there. The redistribution of angular momentum within the protostar–disc system is facilitated mainly by gravitational instability; this allows formation of relatively large discs even when the specific angular momentum of infalling material is low. We argue that such discs should remain marginally unstable as they grow (with Toomre Q ~ 1–2), an idea that is broadly consistent with recent observational estimates for Class 0/I discs. We discuss the numerical convergence of our results, and show that properly treating the inner boundary condition is crucial for achieving convergence at an acceptable computational cost.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Emergence of the temperature–density relation in the low-density intergalactic medium

We examine the evolution of the phase diagram of the low-density intergalactic medium during the Epoch of Reionization in simulation boxes with varying reionization histories from the Cosmic Reionization on Computers project. The probability density function (PDF) of gas temperature at fixed density exhibits two clear modes: a warm and a cold temperature mode, corresponding to the gas inside and outside of ionized bubbles. We find that the transition between the two modes is ‘universal’ in the sense that its timing is accurately parametrized by the value of the volume-weighted neutral fraction for any reionization history. This ‘universality’ is more complex than just a reflection of the fact that ionized gas is warm and neutral gas is cold: it holds for the transition at a fixed value of gas density, and gas at different densities transitions from the cold to the warm mode at different values of the neutral fraction, reflecting a non-trivial relationship between the ionization history and the evolving gas density PDF. Furthermore, the ‘emergence’ of the tight temperature–density relation in the warm mode is also approximately ‘universally’ controlled by the volume-weighted neutral fraction for any reionization history. In particular, the ‘emergence’ of the temperature–density relation (as quantified by the rapid decrease in its width) occurs when the neutral fraction is 10 –4 ≲ X HI ≲ 10 –3 for any reionization history. Our results indicate that the neutral fraction is a primary quantity controlling the various properties of the temperature–density relation, regardless of reionization history.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Gravitational-wave Constraints on the Equatorial Ellipticity of Millisecond Pulsars

In this work, we present a search for continuous gravitational waves from five radio pulsars, comprising three recycled pulsars (PSR J0437-4715, PSR J0711-6830, and PSR J0737-3039A) and two young pulsars: the Crab pulsar (J0534+2200) and the Vela pulsar (J0835-4510). We use data from the third observing run of Advanced LIGO and Virgo combined with data from their first and second observing runs. For the first time, we are able to match (for PSR J0437-4715) or surpass (for PSR J0711-6830) the indirect limits on gravitational-wave emission from recycled pulsars inferred from their observed spin-downs, and constrain their equatorial ellipticities to be less than 10 -8 . For each of the five pulsars, we perform targeted searches that assume a tight coupling between the gravitational-wave and electromagnetic signal phase evolution. We also present constraints on PSR J0711-6830, the Crab pulsar, and the Vela pulsar from a search that relaxes this assumption, allowing the gravitational-wave signal to vary from the electromagnetic expectation within a narrow band of frequencies and frequency derivatives.

79 ASTRONOMY AND ASTROPHYSICS↗