Search NASA⌕ Search

SEARCH · Search NASA

Results for “STAR ACCRETION”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Real-time discovery of AT2020xnd: a fast, luminous ultraviolet transient with minimal radioactive ejecta

ABSTRACT The many unusual properties of the enigmatic AT2018cow suggested that at least some subset of the empirical class of fast blue optical transients (FBOTs) represents a genuinely new astrophysical phenomenon. Unfortunately, the intrinsic rarity and fleeting nature of these events have made it difficult to identify additional examples early enough to acquire the observations necessary to constrain theoretical models. We present here the Zwicky Transient Facility discovery of AT2020xnd (ZTF20acigmel, the ‘Camel’) at z = 0.243, the first unambiguous AT2018cow analogue to be found and confirmed in real time. AT2018cow and AT2020xnd share all key observational properties: a fast optical rise, sustained high photospheric temperature, absence of a second peak attributable to ejection of a radioactively heated stellar envelope, extremely luminous radio, millimetre, and X-ray emission, and a dwarf-galaxy host. This supports the argument that AT2018cow-like events represent a distinct phenomenon from slower-evolving radio-quiet supernovae, likely requiring a different progenitor or a different central engine. The sample properties of the four known members of this class to date disfavour tidal disruption models but are consistent with the alternative model of an accretion powered jet following the direct collapse of a massive star to a black hole. Contextual filtering of alert streams combined with rapid photometric verification using multiband imaging provides an efficient way to identify future members of this class, even at high redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Ultra high-energy cosmic rays from relativistic outflows in accretion induced collapse of white dwarfs

When a rapidly-rotating, highly magnetized white dwarf (WD) approaches the Chandrashekhar limit through mass accretion, it can undergo an accretion-induced collapse (AIC) to form a proto-neutron star or protomagnetar. The protomagnetar can drive a magnetically-dominated relativistic outflow, whose low entropy can lead to efficient formation of heavy nuclei. In this work, we propose that such relativistic outflows from AIC of WDs can contribute as sources of ultra high-energy cosmic rays (UHECRs). We model the acceleration of heavy nuclei in these relativistic outflows, and show that AICs can dominantly power the observed UHECRs, if a majority of them host relativistic outflows. Accounting for uncertainties in the acceleration mechanisms and AIC rates, AICs can contribute $\sim$ a few $10^{43} - 10^{45}\ {\rm erg \ Mpc}^{-3} {\rm yr}^{-1}$ in UHECR energy generation rate density, assuming iron-like nuclei.

Mukhopadhyay, Mainak [Fermilab; Chicago U., KICP] ↗
DOE OSTI2021

MOA-2009-BLG-319Lb: A Sub-Saturn Planet inside the Predicted Mass Desert

We present an adaptive optics (AO) analysis of images from the Keck II telescope NIRC2 instrument of the planetary microlensing event MOA-2009-BLG-319. The ∼10 yr baseline between the event and the Keck observations allows the planetary host star to be detected at a separation of 66.5 ± 1.7 mas from the source star, consistent with the light-curve model prediction. The combination of the host star brightness and light-curve parameters yields host star and planet masses of M {sub host} = 0.524 ± 0.048M {sub ⊙} and m {sub p} = 67.3 ± 6.2M {sub ⊕} at a distance of D {sub L} = 7.1 ± 0.7 kpc. The star−planet projected separation is 2.03 ± 0.21 au. The planet-to-star mass ratio of this system, q = (3.857 ± 0.029) × 10{sup −4}, places it in the predicted “planet desert” at 10{sup −4} < q < 4 × 10{sup −4} according to the runaway gas accretion scenario of the core accretion theory. Seven of the 30 planets in the Suzuki et al. sample fall in this mass ratio range, and this is the third with a measured host mass. All three of these host stars have masses of 0.5 ≤ M {sub host}/M {sub ⊙} ≤ 0.7, which implies that this predicted mass ratio gap is filled with planets that have host stars within a factor of two of 1M {sub ⊙}. This suggests that runaway gas accretion does not play a major role in determining giant planet masses for stars somewhat less massive than the Sun. Our analysis has been accomplished with a modified DAOPHOT code that has been designed to measure the brightness and positions of closely blended stars. This will aid in the development of the primary method that the Nancy Grace Roman Space Telescope mission will use to determine the masses of microlens planets and their hosts.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Electron and muon dynamics in neutron stars beyond chemical equilibrium

A neutron star harbors $\mathcal{O}$(10 56 ) electrons in its core, and almost the same number of muons, with muon decay prohibited by Pauli blocking. However, as macroscopic properties of the star such as its mass, rotational velocity, or magnetic field evolve over time, the equilibrium lepton abundances (dictated by the weak interactions) change as well. Scenarios where this can happen include spin-down, accretion, magnetic field decay, and tidal deformation. Here, we discuss the mechanisms by which a star disrupted in one of these ways re-establishes lepton chemical equilibrium. In most cases, the dominant processes are out-of-equilibrium Urca reactions, the rates of which we compute for the first time. If, however, the equilibrium muon abundance decreases, while the equilibrium electron abundance increases (or decreases less than the equilibrium muon abundance), outward diffusion of muons plays a crucial role as well. This is true in particular for stars older than about 10 4 yrs whose core has cooled to ≲ 20 keV. The muons decay in a region where Pauli blocking is lifted, and we argue that these decays lead to a flux of Script $\mathcal{O}$(10 MeV) neutrinos. Realistically, however, this flux will remain undetectable for the foreseeable future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simulating highly eccentric common envelope jet supernova impostors

We conduct three-dimensional hydrodynamical simulations of eccentric common envelope jet supernova (CEJSN) impostors, i.e. a neutron star that crosses through the envelope of a red supergiant star on a highly eccentric orbit and launches jets as it accretes mass from the envelope. Because of numerical limitations, we apply a simple prescription where we inject the assumed jets’ power into two opposite conical regions inside the envelope. We find the outflow morphology to be very complicated, clumpy, and non-spherical, having a large-scale symmetry only about the equatorial plane. The outflow morphology can substantially differ between simulations that differ by their jets’ power. We estimate by simple means the light curve to be very bumpy, to have a rise time of one to a few months, and to slowly decay in about a year to several years. These eccentric CEJSN impostors will be classified as ‘gap’ objects, i.e. having a luminosity between those of classical novae and typical supernovae (termed also ILOTs for intermediate luminosity optical transients). We strengthen a previous conclusion that CEJSN impostors might account for some peculiar ILOTs, in particular those that might repeat over time-scales of months to years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

GRB 171205A: Hypernova and Newborn Neutron Star

GRB 171205A is a low-luminosity, long-duration gamma-ray burst (GRB) associated with SN 2017iuk, a broad-line type Ic supernova (SN). It is consistent with having been formed in the core collapse of a widely separated binary, which we have called the binary-driven hypernova of type III. The core collapse of the CO star forms a newborn NS (νNS) and the SN explosion. Fallback accretion transfers mass and angular momentum to the νNS, here assumed to be born non-rotating. The accretion energy injected into the expanding stellar layers powers the prompt emission. The multiwavelength power-law afterglow is explained by the synchrotron radiation of electrons in the SN ejecta, powered by energy injected by the spinning νNS. We calculate the amount of mass and angular momentum gained by the νNS, as well as the νNS rotational evolution. The νNS spins up to a period of 47 ms, then releases its rotational energy powering the synchrotron emission of the afterglow. The paucity of the νNS spin explains the low-luminosity characteristic and that the optical emission of the SN from the nickel radioactive decay outshines the optical emission from the synchrotron radiation. From the νNS evolution, we infer that the SN explosion had to occur at most 7.36 h before the GRB trigger. Therefore, for the first time, the analysis of the GRB data leads to the time of occurrence of the CO core collapse leading to the SN explosion and the electromagnetic emission of the GRB event.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Evolution of Inclined Binary Black Holes in the Disks of Active Galactic Nuclei

The accretion disks that fuel active galactic nuclei (AGNs) may house numerous stars and compact objects, formed in situ or captured from nearby star clusters. Embedded neutron stars and black holes may form binaries and eventually merge, emitting gravitational waves detectable by LIGO/VIRGO. AGN disks are a particularly promising environment for the production of high-mass gravitational-wave events involving black holes in the pair-instability mass gap, and may facilitate electromagnetic counterparts to black hole binary mergers. However, many orders of magnitude separate the typical length scales of binary formation and those on which gravitational waves can drive binary inspirals, making binary mergers inside the disk uncertain. Previous hydrodynamical simulations of binaries have either been restricted to two dimensions entirely, or focused on binaries aligned with the midplane of the disk. Herein we present the first three-dimensional, high-resolution, local-shearing-box, inviscid hydrodynamical simulations of disk-embedded binaries over a range of orbital inclinations. We find that retrograde binaries can shrink up to 4 times as quickly as prograde binaries, and that all binaries not perfectly aligned (or anti-aligned) with the AGN disk are driven into alignment. An important consequence of this is that initially retrograde binaries will traverse the inclinations where von Zeipel–Lidov–Kozai oscillations can drive binary eccentricities to large values, potentially facilitating mergers. We also find that interactions with the AGN disk may excite eccentricities in retrograde binaries and cause the orbits of embedded binaries to precess.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Gravitational-wave signature of core-collapse supernovae

Here, we calculate the gravitational-wave (GW) signatures of detailed 3D core-collapse supernova sim ulations spanning a range of massive stars. Most of the simulations are carried out to times late enough to capture more than 95% of the total GW emission. We find that the f/g-mode and f-mode of proto-neutron star oscillations carry away most of the GW power. The f-mode frequency inex orably rises as the proto-neutron star (PNS) core shrinks. We demonstrate that the GW emission is excited mostly by accretion plumes onto the PNS that energize modal oscillations and also high frequency (“haze”) emission correlated with the phase of violent accretion. The duration of the major phase of emission varies with exploding progenitor and there is a strong correlation between the total GW energy radiated and the compactness of the progenitor. Moreover, the total GW emissions vary by as much as three orders of magnitude from star to star. For black-hole formation, the GW signal tapers off slowly and does not manifest the haze seen for the exploding models. For such failed models, we also witness the emergence of a spiral shock motion that modulates the GW emission at a frequency near ~100 Hertz that slowly increases as the stalled shock sinks. We find significant angular anisotropy of both the high- and low-frequency (memory) GW emissions, though the latter have very little power.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Design and Scaling of an Omega-EP Experiment to Study Cold Streams Feeding Early Galaxies

Galaxies form in the centers of dark matter halos, and grow by accreting gas from the halos. The gas supply to the galaxy is the bottleneck for star formation and ultimately sets its stellar mass and star formation rate. As gas falls from the cosmic web into the halos, the highly supersonic infalling gas (Mach ~20 for Milky Way-type halos) eventually forms a shock. If cooling is inefficient, this shock quickly stabilizes as a quasi-static slowly expanding accretion shock at the outer edge of the halo. If, however, cooling is important, a recent theory suggests that filaments of gas from the cosmic web will freefall unshocked through the halos, and shock only at the outskirts of the central galaxies. Furthermore, these "cold streams" allow for a steady supply of gas into the galaxies and for more efficient star formation. The cold dense filament flowing into a hot less dense environment is potentially Kelvin–Helmholtz unstable. This instability may hinder the ability of the stream to deliver gas deeply enough into the halo. A design of a well-scaled laser experiment on Omega-EP meant for studying this phenomena is presented in the current work, with relevant theory informed by 2D hydrodynamic simulations of the experiment. We establish the hydrodynamic scaling analysis between the cosmological system and its experimental analog, presenting the experiment as an adiabatic upper limit to informing the role of mixing due to the Kelvin–Helmholtz instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Stellar Rotation of T Tauri Stars in the Orion Star-forming Complex

We present a large-scale study of stellar rotation for T Tauri stars in the Orion star-forming complex. We use the projected rotational velocity ($v\sin (i)$) estimations reported by the APOGEE-2 collaboration as well as individual masses and ages derived from the position of the stars in the HR diagram, considering Gaia-EDR3 parallaxes and photometry plus diverse evolutionary models. We find an empirical trend for $v\sin (i)$ decreasing with age for low-mass stars (0.4M ⊙ < M $\star$ < 1.2M ⊙ ). Our results support the existence of a mechanism linking $v\sin (i)$ to the presence of accreting protoplanetary disks, responsible for regulating stellar rotation on timescales of about 6 Myr, which is the timescale in which most of the T Tauri stars lose their inner disk. Our results provide important constraints to models of rotation in the early phases of evolution of young stars and their disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Spatially resolved gas and stellar kinematics in compact starburst galaxies

The kinematics of galaxies provide valuable insights into their physics and assembly history. Kinematics are governed not only by the gravitational potential, but also by merger events and stellar feedback processes such as stellar winds and supernova explosions. We aim to identify what governs the kinematics in a sample of SDSS-selected nearby starburst galaxies, by obtaining spatially resolved measurements of the gas and stellar kinematics. We obtained near-infrared integral-field K-band spectroscopy with VLT/SINFONI for 15 compact starburst galaxies. We derived the integrated as well as spatially resolved stellar and gas kinematics. The stellar kinematics were derived from the CO absorption bands, and Paα and Brγ emission lines were used for the gas kinematics. Based on the integrated spectra, we find that the majority of galaxies have gas and stellar velocity dispersion that are comparable. A spatially resolved comparison shows that the six galaxies that deviate show evidence for a bulge or stellar feedback. Two galaxies are identified as mergers based on their double-peaked emission lines. In our sample, we find a negative correlation between the ratio of the rotational velocity over the velocity dispersion (v rot /σ) and the star formation rate surface density. We propose a scenario where the global kinematics of the galaxies are determined by gravitational instabilities that affect both the stars and gas. This process could be driven by mergers or accretion events. Effects of stellar feedback on the ionised gas are more localised and detected only in the spatially resolved analysis. The mass derived from the velocity dispersion provides a reliable mass even if the galaxy cannot be spatially resolved. The technique used in this paper is applicable to galaxies at low and high redshift with the next generation of infrared-focussed telescopes (JWST and ELT).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Collapsar Gamma-Ray Bursts Grind Their Black Hole Spins to a Halt

Abstract The spin of a newly formed black hole (BH) at the center of a massive star evolves from its natal value due to two competing processes: accretion of gas angular momentum that increases the spin and extraction of BH angular momentum by outflows that decreases the spin. Ultimately, the final, equilibrium spin is set by a balance between both processes. In order for the BH to launch relativistic jets and power aγ-ray burst (GRB), the BH magnetic field needs to be dynamically important. Thus, we consider the case of a magnetically arrested disk (MAD) driving the spin evolution of the BH. By applying the semianalytic MAD BH spin evolution model of Lowell et al. to collapsars, we show that if the BH accretes ∼20% of its initial mass, its dimensionless spin inevitably reaches small values,a≲ 0.2. For such spins, and for mass accretion rates inferred from collapsar simulations, we show that our semianalytic model reproduces the energetics of typical GRB jets,L jet ∼ 10 50 erg s −1 . We show that our semianalytic model reproduces the nearly constant power of typical GRB jets. If the MAD onset is delayed, this allows powerful jets at the high end of the GRB luminosity distribution,L jet ∼ 10 52 erg s −1 , but the final spin remains low,a≲ 0.3. These results are consistent with the low spins inferred from gravitational wave detections of binary BH mergers. In a companion paper by Gottlieb et al., we use GRB observations to constrain the natal BH spin to bea≃ 0.2.

Astronomy & Astrophysics↗
DOE OSTI2023

A resolved study of the inner regions of nearby galaxies with an excess of young massive stars: missing link in the AGN–starburst connection?

ABSTRACT We have selected galaxies with very high levels of H α emission (EQW(H α) >700 Å.) in their central regions from the final data release of the MaNGA survey. Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range 9.5 < log M*/(M⊙) < 11.5. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populations of young massive stars, finding that selection in the 2D space of extinction-corrected H α EQW and [S iii]/[S ii] line ratio produces the best results. By comparing stacked spectra covering these regions with stacked spectra covering normal starburst regions with 100 Å<EQW(H α) <200 Å, we obtain the following main results: (1) Clear signatures of excess Wolf–Rayet stars are found in half of the H α excess regions, (2) galaxy regions containing excess Wolf–Rayet stars are more often associated with the presence of high-ionization emission lines characteristic of accreting black holes. Excess [Ne iii] is detected in 4 out of 8 of the WR regions and there are tentative [Fe x] detections in 2 galaxies. (3) Regions of the galaxy with excess Wolf–Rayet stars are located where the interstellar medium has larger ionized gas turbulent velocities and higher neutral gas overdensities. We make a first attempt to constrain changes in the high-mass end of the stellar initial mass function (IMF) using the HR-pyPopStar evolutionary population synthesis models that include high-wavelength-resolution theoretical atmosphere libraries for Wolf–Rayet stars.

Kauffmann, Guinevere↗
DOE OSTI2022

Galactic ArchaeoLogIcaL ExcavatiOns (GALILEO): I. An updated census of APOGEE N-rich giants across the Milky Way

We use the 17th data release of the second phase of the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) to provide a homogenous census of N-rich red giant stars across the Milky Way (MW). We report a total of 149 newly identified N-rich field giants toward the bulge, metal-poor disk, and halo of our Galaxy. They exhibit significant enrichment in their nitrogen abundance ratios ([N/Fe] ≳ +0.5), along with simultaneous depletions in their [C/Fe] abundance ratios ([C/Fe] < +0.15), and they cover a wide range of metallicities (-1.8 < [Fe/H] < -0.7). The final sample of candidate N-rich red giant stars with globular-cluster-like (GC-like) abundance patterns from the APOGEE survey includes a grand total of ~412 unique objects. These strongly N-enhanced stars are speculated to have been stripped from GCs based on their chemical similarities with these systems. Even though we have not found any strong evidence for binary companions or signatures of pulsating variability yet, we cannot rule out the possibility that some of these objects were members of binary systems in the past and/or are currently part of a variable system. In particular, the fact that we identify such stars among the field stars in our Galaxy provides strong evidence that the nucleosynthetic process(es) producing the anomalous [N/Fe] abundance ratios occurs over a wide range of metallicities. This may provide evidence either for or against the uniqueness of the progenitor stars to GCs and/or the existence of chemical anomalies associated with likely tidally shredded clusters in massive dwarf galaxies such as “Kraken/Koala”, Gaia-Enceladus-Sausage, among others, before or during their accretion by the MW. In conclusion, a dynamical analysis reveals that the newly identified N-rich stars exhibit a wide range of dynamical characteristics throughout the MW, indicating that they were produced in a variety of Galactic environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Birth of a Be star: an APOGEE search for Be stars forming through binary mass transfer

ABSTRACT Motivated by recent suggestions that many Be stars form through binary mass transfer, we searched the APOGEE survey for Be stars with bloated, stripped companions. From a well-defined parent sample of 297 Be stars, we identified one mass-transfer binary, HD 15124. The object consists of a main-sequence Be star ($M_{\rm Be}=5.3\pm 0.6\, {\rm M}_{\odot }$) with a low-mass ($M_{\rm donor}=0.92\pm 0.22\, {\rm M}_{\odot }$), subgiant companion on a 5.47-d orbit. The emission lines originate in an accretion disc caused by ongoing mass transfer, not from a decretion disc as in classical Be stars. Both stars have surface abundances bearing imprint of CNO processing in the donor’s core: the surface helium fraction is YHe ≈ 0.6, and the nitrogen-to-carbon ratio is 1000 times the solar value. The system’s properties are well-matched by binary evolution models in which mass transfer begins while a $3-5\, {\rm M}_{\odot }$ donor leaves the main sequence, with the originally less massive component becoming the Be star. These models predict that the system will soon become a detached Be + stripped star binary like HR 6819 and LB-1, with the stripped donor eventually contracting to become a core helium-burning sdO/B star. Discovery of one object in this short-lived (∼1 Myr) evolutionary phase implies the existence of many more that have already passed through it and are now Be + sdO/B binaries. We infer that $(10-60)\, {{\ \rm per\ cent}}$ of Be stars have stripped companions, most of which are $\sim 100\, \times$ fainter than the Be stars in the optical. Together with the dearth of main-sequence companions to Be stars and recent discovery of numerous Be + sdO/B binaries in the UV, our results imply that binarity plays an important role in the formation of Be stars.

El-Badry, Kareem (ORCID:0000000268711752)↗
DOE OSTI2022

The role of the turbulence driving mode for the initial mass function

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating, and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ∼1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Diffstar: a fully parametric physical model for galaxy assembly history

We present Diffstar, a smooth parametric model for the in situ star formation history (SFH) of galaxies. The Diffstar model is distinct from traditional SFH models because it is parametrized directly in terms of basic features of galaxy formation physics. Diffstar includes ingredients for: the halo mass assembly history; the accretion of gas into the dark matter halo; the fraction of gas that is eventually transformed into stars, ϵ ms ; the time-scale over which this transformation occurs, τ cons ; and the possibility that some galaxies will experience a quenching event at time t q , and may subsequently experience rejuvenated star formation. We show that our model is sufficiently flexible to describe the average stellar mass histories of galaxies in both the IllustrisTNG (TNG) and UniverseMachine (UM) simulations with an accuracy of ~0.1 dex across most of cosmic time. We use Diffstar to compare TNG to UM in common physical terms, finding that: (i) star formation in UM is less efficient and burstier relative to TNG; (ii) UM galaxies have longer gas consumption time-scales, relative to TNG; (iii) rejuvenated star formation is ubiquitous in UM, whereas quenched TNG galaxies rarely experience sustained rejuvenation; and (iv) in both simulations, the distributions of ϵ ms , τ cons , and t q share a common characteristic dependence upon halo mass, and present significant correlations with halo assembly history. Finally, we conclude with a discussion of how Diffstar can be used in future applications to fit the SEDs of individual observed galaxies, as well as in forward-modelling applications that populate cosmological simulations with synthetic galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Formation and evolution of protostellar accretion discs – II. From 3D simulation to a simple semi-analytic model of Class 0/I discs

Here we use a 3D radiative non-ideal magnetohydrodynamic simulation to investigate the formation and evolution of a young protostellar disc from a magnetized pre-stellar core. The simulation covers the first ${\sim }10\, {\rm kyr}$ after protostar formation and shows a massive, weakly magnetized disc with radius that initially grows and then saturates at ${\sim }30\, {\rm au}$. The disc is gravitationally unstable with prominent large-amplitude spiral arms. We use our simulation results and a series of physical arguments to construct a predictive and quantitative physical picture of Class 0/I protostellar disc evolution from several aspects, including (i) the angular-momentum redistribution in the disc, self-regulated by gravitational instability to make most of the disc marginally unstable; (ii) the thermal profile of the disc, well-approximated by a balance between radiative cooling and accretion heating; and (iii) the magnetic-field strength and magnetic-braking rate inside the disc, regulated by non-ideal magnetic diffusion. Using these physical insights, we build a simple 1D semi-analytic model of disc evolution. We show that this 1D model, when coupled to a computationally inexpensive simulation for the evolution of the surrounding pseudo-disc, can be used reliably to predict disc evolution in the Class 0/I phase. The predicted long-term evolution of disc size, which saturates at ${\sim }30\, {\rm au}$ and eventually shrinks, is consistent with a recent observational survey of Class 0/I discs. Such hierarchical modelling of disc evolution circumvents the computational difficulty of tracing disc evolution through Class 0/I phase with direct, numerically converged simulations.

79 ASTRONOMY AND ASTROPHYSICS↗