Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Characterizing the Variation and Covariation of Cloud Microphysical Properties and Implications for Simulation of Subgrid-scale Warm-Rain Processes in Earth System Models (Final DOE-ASR Report)

Warm marine boundary layer (MBL) clouds constitute an important component in the global climate system, and precipitation plays a central role in controlling the water budget, radiative effects, and lifetime of these MBL clouds. Unfortunately, because of the relatively coarse effective grid resolution of the current generation of Earth system models (ESMs), the variety of cloud microphysical processes occurring inside an ESM grid cell are often oversimplified or unconstrained by observations. For example, the warm rain processes (e.g., autoconversion and accretion) are usually parameterized as nonlinear functions of grid-mean cloud properties. Because of the nonlinear nature of these functions, neglecting variability within the ESM grid volume can lead to substantial biases in precipitation production, cloud cover, and surface radiative fluxes. In state-of-the art ESMs, the influence of subgrid-scale variability is represented as an enhancement factor (EF) coefficient to the autoconversion, and accretion rates calculated from the model variables. However, EF is typically taken to be a constant or even used as a knob to tune model cloud properties to match observations, an ad hoc approach that may yield a desired cloud outcome yet introduce compensating errors. In this project, we used the combination of in situ cloud microphysics measurements from the ACE-ENA field campaign and large-eddy simulations (LES) to characterize and understand subgrid-scale variations and co-variations of cloud microphysical properties and use the results to evaluate and improve the representation of subgrid warm-rain processes in ESMs, in particular the EF used to tune the autoconversion and accretion processes. In this final report, we summarize our research activities and main findings in Section 2, provide a list of publications (Section 3) and presentations (Section 4) resulted from our research, and briefly discuss the student activities supported by this project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Periodic Eruptive Variability of the Isolated Pre-main-sequence Star V347 Aurigae

V347 Aurigae is associated with the small dark cloud L1438 and appears to be an isolated pre-main-sequence star located at distance d ≈ 200 pc. Multi-epoch, archival photometry reveals periodic brightness variations with amplitude V ≈ 2.0 mag occurring on timescales of ~160 days that have persisted for decades. Regular-cadence, optical imaging of the source with the Zwicky Transient Facility shows that a small reflection nebula illuminated by V347 Aur also fluctuates in brightness, at times fading completely. Multi-epoch, Keck/HIRES data suggests the presence of two distinct spectral components: a prominent emission-line-dominated spectrum with a heavily veiled continuum correlated with the bright photometric state, and an M-type absorption line spectrum associated with quiescence. All spectra exhibit strong Balmer and He i line emission, consistent with accretion, as well as high velocity emission arising from the forbidden transitions of [O i], [N ii], and [S ii] that are generally associated with collimated jets and disk winds. There is no evidence in existing high-dispersion spectroscopy or high-resolution imaging for binarity of V347 Aur. The repeating outburst events are possibly linked to accretion instabilities induced by an undetected companion or a structure within the circumstellar disk that periodically increases the mass accretion rate. V347 Aur is perhaps analogous to an EXor-type variable, though more regularly recurring.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Role of Early Giant-planet Instability in Terrestrial Planet Formation

The terrestrial planets are believed to have formed by violent collisions of tens of lunar- to Mars-size protoplanets at time t < 200 Myr after the protoplanetary gas disk dispersal (t {sub 0}). The solar system giant planets rapidly formed during the protoplanetary disk stage and, after t {sub 0}, radially migrated by interacting with outer disk planetesimals. An early (t < 100 Myr) dynamical instability is thought to have occurred with Jupiter having gravitational encounters with a planetary-size body, jumping inward by ∼0.2–0.5 au, and landing on its current, mildly eccentric orbit. Here we investigate how the giant-planet instability affected the formation of the terrestrial planets. We study several instability cases that were previously shown to match many solar system constraints. We find that resonances with giant planets help to remove solids available for accretion near ∼1.5 au, thus stalling the growth of Mars. It does not matter, however, whether the giant planets are placed on their current orbits at t {sub 0} or whether they realistically evolve in one of our instability models; the results are practically the same. The tight orbital spacing of Venus and Earth is difficult to reproduce in our simulations, including cases where bodies grow from a narrow annulus at 0.7–1 au, because protoplanets tend to spread radially during accretion. The best results are obtained in the narrow-annulus model when protoplanets emerging from the dispersing gas nebula are assumed to have (at least) the mass of Mars. This suggests efficient accretion of the terrestrial protoplanets during the first ∼10 Myr of the solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2021

Searching for Active Low-mass Stars in the CMa Star-forming Region: Multi-band Photometry with T80S

An exotic environment surrounds the young stellar groups associated with the Canis Major (CMa) OB1/R1 region, which probably was formed under feedback from at least three supernova events having occurred a few million years ago. We use astrometric data from the Gaia-DR2 to confirm the membership of the stars in CMa R1, based on proper motion and parallax, which revealed 514 new members and candidates. The mean age of 5 Myr estimated from the color–magnitude diagram characterizes the sources as likely pre-main-sequence candidates. In total, a sample of 694 stars detected with the T80-South telescope was analyzed according to different color–color diagrams, which were compared with theoretical colors from evolutionary models, aiming to reveal the objects that exhibit color excess due to accretion processes. Accretion and magnetic activity were also explored on the basis of empirical flux–flux relation, such as F {sub 660} and F {sub 861} that are related to Hα and Ca ii triplet emission, respectively. A low fraction (∼3 percent) of the sample have Hα excess and other colors expected for stars exhibiting chromospheric activity. The number of Class I and Class II objects, identified by the infrared (Wide-field Infrared Survey Explorer) colors, indicates a disk fraction of ∼6 percent, which is lower than that expected for stellar clusters with a similar age. Such a large sample of objects associated with CMa R1 without evidences of circumstellar accretion can be interpreted as a lack of disk-bearing stars, unusual for young star-forming regions. However, this may be explained as the result of supernova events.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Evidence that the Hot Jupiter WASP-77 A b Formed Beyond Its Parent Protoplanetary Disk’s H2O Ice Line

Idealized protoplanetary disk and giant planet formation models have been interpreted to suggest that a giant planet's atmospheric abundances can be used to infer its formation location in its parent protoplanetary disk. It has recently been reported that the hot Jupiter WASP-77 A b has subsolar atmospheric carbon and oxygen abundances with a solar C/O abundance ratio. Assuming solar carbon and oxygen abundances for its host star WASP-77 A, WASP-77 A b's atmospheric carbon and oxygen abundances possibly indicate that it accreted its envelope interior to its parent protoplanetary disk's H 2 O ice line from carbon-depleted gas with little subsequent planetesimal accretion or core erosion. We show that the photospheric abundances of carbon and oxygen in WASP-77 A are supersolar with a subsolar C/O abundance ratio, implying that WASP-77 A b's atmosphere has significantly substellar carbon and oxygen abundances with a superstellar C/O ratio. Our result possibly indicates that WASP-77 A b's envelope was accreted by the planet beyond its parent protoplanetary disk's H 2 O ice line. While numerous theoretical complications to these idealized models have now been identified, the possibility of nonsolar protoplanetary disk abundance ratios confound even the most sophisticated protoplanetary disk and giant planet formation models. We therefore argue that giant planet atmospheric abundance ratios can only be meaningfully interpreted relative to the possibly nonsolar mean compositions of their parent protoplanetary disks as recorded in the photospheric abundances of their dwarf host stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

NS 1987A in SN 1987A

The possible detection of a compact object in the remnant of SN 1987A presents an unprecedented opportunity to follow its early evolution. The suspected detection stems from an excess of infrared emission from a dust blob near the compact object’s predicted position. The infrared excess could be due to the decay of isotopes like 44 Ti, accretion luminosity from a neutron star or black hole, magnetospheric emission or a wind originating from the spin down of a pulsar, or to thermal emission from an embedded, cooling neutron star. In this work it is shown that the last possibility is the most plausible as the other explanations are disfavored by other observations and/or require fine-tuning of parameters. Not only are there indications that the dust blob overlaps the predicted location of a kicked compact remnant, but its excess luminosity also matches the expected thermal power of a 30 yr old neutron star. Furthermore, models of cooling neutron stars within the minimal cooling paradigm readily fit both NS 1987A and Cas A, the next-youngest known neutron star. If correct, a long heat transport timescale in the crust and a large effective stellar temperature are favored, implying relatively limited crustal n- 1 S 0 superfluidity and an envelope with a thick layer of light elements, respectively. If the locations do not overlap, then pulsar spin down or accretion might be more likely, but the pulsar’s period and magnetic field or the accretion rate must be rather finely tuned. In this case, NS 1987A may have enhanced cooling and/or a heavy-element envelope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Wide Dust Gaps in Protoplanetary Disks Induced by Eccentric Planets: A Mass-eccentricity Degeneracy

The tidal perturbation of embedded protoplanets on their natal disks has been widely attributed to be the cause of gap-ring structures in submillimeter images of protoplanetary disks around T Tauri stars. Numerical simulations of this process have been used to propose scaling of characteristic dust-gap width/gap-ring distance with respect to planet mass. Applying such scaling to analyze observed gap samples yields a continuous mass distribution for a rich population of hypothetical planets in the range of several Earth to Jupiter masses. In contrast, the conventional core-accretion scenario of planet formation predicts a bimodal mass function due to (1) the onset of runaway gas accretion above ~20 Earth masses and (2) suppression of accretion induced by gap opening. Here, we examine the dust disk response to the tidal perturbation of eccentric planets as a possible resolution of this paradox. Based on simulated gas and dust distributions, we show the gap-ring separation of Neptune-mass planets with small eccentricities might become comparable to that induced by Saturn-mass planets on circular orbits. Furthermore, this degeneracy may obliterate the discrepancy between the theoretical bimodal mass distribution and the observed continuous gap width distribution. Despite damping due to planet–disk interaction, modest eccentricity may be sustained either in the outer regions of relatively thick disks or through resonant excitation among multiple super Earths. Moreover, the ring-like dust distribution induced by planets with small eccentricities is axisymmetric even in low viscosity environments, consistent with the paucity of vortices in Atacama Large Millimeter/submillimeter Array images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Exploring Changes in Quasar Spectral Energy Distributions across C IV Parameter Space

We examine the UV/X-ray properties of 1378 quasars in order to link empirical correlations to theoretical models of the physical mechanisms dominating quasars as a function of mass and accretion rate. The clarity of these correlations is improved when (1) using C IV broad emission line equivalent width (EQW) and blueshift (relative to systemic) values calculated from high signal-to-noise ratio reconstructions of optical/UV spectra and (2) removing quasars expected to be absorbed based on their UV/X-ray spectral slopes. In addition to using the traditional C IV parameter space measures of C IV EQW and blueshift, we define a "C IV ∥ distance" along a best-fit polynomial curve that incorporates information from both C IV parameters. We find that the C IV ∥ distance is linearly correlated with both the optical-to-X-ray slope, α ox , and broad-line He II EQW, which are known spectral energy distribution indicators, but does not require X-ray or high spectral resolution UV observations to compute. The C IV ∥ distance may be a better indicator of the mass-weighted accretion rate, parameterized by L/L Edd , than the C iv EQW or blueshift alone, as those relationships are known to break down at the extrema. Conversely, there is only a weak correlation with the X-ray energy index (Γ), an alternate L/L Edd indicator. We find no X-ray or optical trends in the direction perpendicular to the C IV distance that could be used to reveal differences in accretion disk, wind, or corona structure that could be widening the C IV EQW–blueshift distribution. A different parameter (such as metallicity) not traced by these data must come into play.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Angular-momentum Transport in Proto-neutron Stars and the Fate of Neutron Star Merger Remnants

Both the core collapse of rotating massive stars, and the coalescence of neutron star (NS) binaries result in the formation of a hot, differentially rotating NS remnant. The timescales over which differential rotation is removed by internal angular-momentum transport processes (viscosity) have key implications for the remnant’s long-term stability and the NS equation of state (EOS). Guided by a nonrotating model of a cooling proto-NS, we estimate the dominant sources of viscosity using an externally imposed angular-velocity profile Ω(r). Although the magneto-rotational instability provides the dominant source of effective viscosity at large radii, convection and/or the Tayler–Spruit dynamo dominate in the core of merger remnants where dΩ/dr ≥ 0. Furthermore, the viscous timescale in the remnant core is sufficiently short that solid-body rotation will be enforced faster than matter is accreted from rotationally supported outer layers. Guided by these results, we develop a toy model for how the merger remnant core grows in mass and angular momentum due to accretion. We find that merger remnants with sufficiently massive and slowly rotating initial cores may collapse to black holes via envelope accretion, even when the total remnant mass is less than the usually considered threshold ≈1.2 M TOV for forming a stable solid-body rotating NS remnant (where M TOV is the maximum nonrotating NS mass supported by the EOS). This qualitatively new picture of the post-merger remnant evolution and stability criterion has important implications for the expected electromagnetic counterparts from binary NS mergers and for multimessenger constraints on the NS EOS

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Contracting and Expanding Binary Black Holes in 3D Low-mass AGN Disks: The Importance of Separation

LIGO/Virgo has detected several binary black hole (BBH) merger events that may have originated in the accretion disks of active galactic nuclei (AGNs). These events require individual black hole masses that fall within the pair instability supernova mass gap, and therefore these black holes may have been grown from hierarchical mergers. AGN disks are a prime environment for hierarchical mergers and thus a potential location for the progenitors of BBH gravitational wave events. Understanding how a BBH embedded in an AGN disk interacts with the surrounding environment is thus crucial for determining whether this interaction can lead to its merger. However, there are few high-fidelity simulations of this process, and almost all are two-dimensional. We present the results from 3D, high-resolution, local shearing-box simulations of an embedded BBH interacting with an AGN disk. In these first simulations of their kind, we focus on determining the mass accretion rate and the orbital evolution rate at different BBH separations. We find that circular, equal-mass BBHs with separations greater than 10% of their Hill radius contract while accreting at a super-Eddington rate. At smaller separations, however, our 3D simulations find that BBHs expand their orbits. This result suggests that it may be difficult for an AGN disk to push a BBH to merger, but we discuss several mechanisms, including MHD turbulence and radiative and mechanical feedback, which could alleviate this difficulty.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Influence of Disk Composition on the Evolution of Stars in the Disks of Active Galactic Nuclei

Abstract Disks of gas accreting onto supermassive black holes, powering active galactic nuclei (AGN), can capture stars from nuclear star clusters or form stars in situ via gravitational instability. The density and thermal conditions of these disks can result in rapid accretion onto embedded stars, dramatically altering their evolution in comparison to stars in the interstellar medium. Theoretical models predict that, when subjected to sufficiently rapid accretion, fresh gas replenishes hydrogen in the cores of these stars as quickly as it is burned into helium, reaching a quasi-steady state. Such massive, long-lived (“immortal”) stars may be capable of dramatically enriching AGN disks with helium, and would increase the helium abundance in AGN broad-line regions relative to that in the corresponding narrow-line regions and hosts. We investigate how the helium abundance of AGN disks alters the evolution of stars embedded therein. We find, in agreement with analytical arguments, that stars at a given mass are more luminous at higher helium mass fractions, and so undergo more radiation-driven mass loss. We further find that embedded stars tend to be less massive in disks with higher helium mass fractions, and that immortal stars are less common in such disks. Thus, disk composition can alter the rates of electromagnetic and gravitational wave transients as well as further chemical enrichment by embedded stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Metal-weak Milky Way Stellar Disk Hidden in the Gaia–Sausage–Enceladus Debris: The APOGEE DR17 View

We have for the first time identified the early stellar disk in the Milky Way by using a combination of elemental abundances and kinematics. Using data from APOGEE DR17 and Gaia we select stars in the Mg–Mn–Al–Fe plane with elemental abundances indicative of an accreted origin and find stars with both halo-like and disk-like kinematics. The stars with halo-like kinematics lie along a lower sequence in [Mg/Fe], while the stars with disk-like kinematics lie along a higher sequence. Combined with astroseismic observations, we determine that the stars with halo-like kinematics are old, 9–11 Gyr, and that the more evolved stellar disk is about 1–2 Gyr younger. We show that the in situ fraction of stars on deeply bound orbits is not small, in fact the inner Galaxy likely harbors a genuine in situ population together with an accreted one. In addition, we show that the selection of the Gaia–Sausage–Enceladus in the E n -L z plane is not very robust. In fact, radically different selection criteria give almost identical elemental abundance signatures for the accreted stars.

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 OSTI2024

The Santa Barbara Binary-disk Code Comparison

We have performed numerical calculations of a binary interacting with a gas disk, using 11 different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, 0.2(a 2 Ω B /v) is the kinematic viscosity. The largest discrepancies between codes resulted from the dimensionality of the setup (3D vs. 2D disks). We find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cosmological Simulations of Stellar Halos with Gaia Sausage–Enceladus Analogs: Two Sausages, One Bun?

Observations of the Milky Way’s stellar halo find that it is predominantly comprised of a radially biased population of stars, dubbed the Gaia Sausage–Enceladus, or GSE. These stars are thought to be debris from dwarf galaxy accretion early in the Milky Way’s history. Though typically considered to be from a single merger, it is possible that the GSE debris has multiple sources. To investigate this possibility, we use the IllustrisTNG50 simulation to identify stellar accretion histories in 98 Milky Way analogs—the largest sample for which such an identification has been performed—and find GSE-like debris in 32, with two-merger GSEs accounting for a third of these cases. Distinguishing single-merger GSEs from two-merger GSEs is difficult in common kinematic spaces, but differences are more evident through chemical abundances and star formation histories. This is because single-merger GSEs are typically accreted more recently than the galaxies in two-merger GSEs: the median infall times (with 16th and 84th percentiles) are $5.{9}_{-2.0}^{+3.3}$ and $10.{7}_{-3.7}^{+1.2}$ Gyr ago for single- and two-merger scenarios, respectively. The systematic shifts in abundances and ages that occur as a result suggest that efforts in modeling these aspects of the stellar halo prove ever-important in understanding its assembly.

galactic archaeology↗
DOE OSTI2025

Low-metallicity Nova Explosions: A Site for Weak rp -process Nucleosynthesis

Classical novae are common cataclysmic events involving a binary system of a white dwarf and a main-sequence or red giant companion star. In metal-poor environments, these explosions produce ejecta differently from their solar counterparts due to the accretion of subsolar-metallicity material onto the white dwarf. In particular, it has been suggested that the nucleosynthesis flow in such low-metallicity nova explosions extends up to the Cu–Zn region, much beyond the expected endpoint, around Ca, predicted for solar-metallicity classical novae. This behavior resembles a weak rp-process, and such nuclear activity has never been observed in accreting white dwarf binaries with typical accretion flows. In this work, we study the characteristics of the weak rp-process for four nova models with metallicities Z = 2 × 10 −9 , 10 −7 , 2 × 10 −6 , and 2 × 10 −5 , and explore the impact of the nuclear physics uncertainties via a Monte Carlo sensitivity study. We identify nuclear reactions whose uncertainties affect the production of intermediate-mass nuclei under these conditions. These reactions and relevant nuclear quantities are targets for measurements at stable or radioactive beam facilities to reduce their rate uncertainties.

Classical novae↗
DOE OSTI2025

The R -process Alliance: Enrichment of r -process Elements in a Simulated Milky Way–like Galaxy

We study the formation of stars with varying amounts of heavy elements synthesized by the rapid neutron-capture process (r-process) based on our detailed cosmological zoom-in simulation of a Milky Way–like galaxy with an N-body/smoothed particle hydrodynamics code, asura. Most stars with no overabundance in r-process elements, as well as the strongly r-process-enhanced (RPE) r-II stars ([Eu/Fe] > +0.7), are formed in dwarf galaxies accreted by the Milky Way within the 6 Gyr after the Big Bang. In contrast, over half of the moderately enhanced r-I stars (+0.3 < [Eu/Fe] ≤ +0.7) are formed in the main in situ disk after 6 Gyr. Our results suggest that the fraction of r-I and r-II stars formed in disrupted dwarf galaxies is larger the higher their [Eu/Fe] is. Accordingly, the most strongly enhanced r-III stars ([Eu/Fe] > +2.0) are formed in accreted components. These results suggest that non-r-process-enhanced stars and r-II stars are mainly formed in low-mass dwarf galaxies that hosted either none or a single neutron star merger, while the r-I stars tend to form in the well-mixed in situ disk. We compare our findings with high-resolution spectroscopic observations of RPE metal-poor stars in the halo and dwarf galaxies, including those collected by the R-Process Alliance. We conclude that observed [Eu/Fe] and [Eu/Mg] ratios can be employed in chemical tagging of the Milky Way’s accretion history.

Chemical enrichment↗