Stellar Mass Accretion Research Articles

Large Fluctuations within 1 au in Protoplanetary Disks

Protoplanetary disks are often assumed to change slowly and smoothly during planet formation. Here, we investigate the time evolution of isolated disks subject to viscosity and a disk wind. The viscosity is assumed to increase rapidly at around 900 K due to thermal ionization of alkali metals, or thermionic and ion emission from dust, and the onset of magnetorotational instability (MRI). The disks generally undergo large, rapid fluctuations for a wide range of time-averaged mass accretion rates. Fluctuations involve coupled waves in temperature and surface density that move radially in either direction through the inner 1.5 au of the disk. Two types of waves are seen with radial speeds of roughly 50 and 1000 cm s−1, respectively. The pattern of waves repeats with a period of roughly 10,000 yr that depends weakly on the average mass accretion rate. Viscous transport due to MRI is confined to the inner disk. This region is resupplied by mass flux from the outer disk driven by the disk wind. Interior to 1 au, the temperature and surface density can vary by a factor of 2–10 on timescales of years to kiloyears. The stellar mass accretion rate varies by 3 orders of magnitude on a similar timescale. This behavior lasts for at least 1 Myr for initial disks comparable to the minimum-mass solar nebula.

NGC 6302: The Tempestuous Life of a Butterfly

NGC 6302 (The Butterfly Nebula) is an extremely energetic and rapidly expanding bipolar planetary nebula (PN). If the central source is a single star, then its apparent location in an H-R diagram places it among the most massive, hottest, and presumably rapidly evolving of all central stars of PNe. Our proper motion study of NGC 6302, based on Hubble Space Telescope WFC3 images spanning 11 yr, has uncovered at least four different pairs of uniformly expanding internal lobes ejected at various times and orientations over the past two millennia at speeds ranging from 10–600 km s−1. In addition, we find a pair of collimated off-axis flows in constant motion at ∼770 ± 100 km s−1 within which bright [Fe ii] feathers are conspicuous. Combining our results with those previously published, we find that the ensemble of flows has an ionized mass >0.1 M ⊙ and its kinetic energy, between 1046 and 1048 erg, lies at the upper end of gravity-powered PNe ejection processes such as stellar mergers or mass accretion. We assemble our results into a plausible historical timeline of ejections from the nucleus and suggest that the ejections are powered by gravitational infall.

Co-evolution of dust grains and protoplanetary disks

Abstract We propose a new evolutionary process for protoplanetary disks, the co-evolution of dust grains and protoplanetary disks, revealed by dust–gas two-fluid non-ideal magnetohydrodynamics simulations considering the growth of dust grains and associated changes in magnetic resistivity. We found that the dust growth significantly affects disk evolution by changing the coupling between the gas and the magnetic field. Moreover, once the dust grains grow sufficiently large and the adsorption of charged particles on to them becomes negligible, the physical quantities (e.g., density and magnetic field) of the disk are well described by characteristic power laws. In this disk structure, the radial profile of density is steeper and the disk mass is smaller than those of the model ignoring dust growth. We analytically derive these power laws from the basic equations of non-ideal magnetohydrodynamics. The analytical power laws are determined only by observable physical quantities, e.g., central stellar mass and mass accretion rate, and do not include difficult-to-determine parameters, e.g., the viscous parameter α. Therefore, our model is observationally testable and this disk structure is expected to provide a new perspective for future studies on protostar and disk evolution.

Gaia Data Release 3

Context. The Gaia Radial Velocity Spectrometer (RVS) provides the unique opportunity of a spectroscopic analysis of millions of stars at medium resolution (λ/Δλ ∼ 11 500) in the near-infrared (845−872 nm). This wavelength range includes the Ca II infrared triplet (IRT) at 850.03, 854.44, and 866.45 nm, which is a good indicator of magnetic activity in the chromosphere of late–type stars. Aims. Here we present the method devised for inferring the Gaia stellar activity index from the analysis of the Ca II IRT in the RVS spectrum, together with its scientific validation. Methods. The Gaia stellar activity index is derived from the Ca II IRT excess equivalent width with respect to a reference spectrum, taking the projected rotational velocity (vsini) into account. We performed scientific validation of the Gaia stellar activity index by deriving a R′IRT index, which is largely independent of the photospheric parameters, and considering the correlation with the R′HK index for a sample of stars. A sample of well-studied pre-main-sequence (PMS) stars is considered to identify the regime in which the Gaia stellar activity index may be affected by mass accretion. The position of these stars in the colour–magnitude diagram and the correlation with the amplitude of the photometric rotational modulation is also scrutinised. Results.Gaia DR3 contains a stellar activity index derived from the Ca II IRT for some 2 × 106 stars in the Galaxy. This represents a ‘gold mine’ for studies on stellar magnetic activity and mass accretion in the solar vicinity. Three regimes of the chromospheric stellar activity are identified, confirming suggestions made by previous authors based on much smaller R′HK datasets. The highest stellar activity regime is associated with PMS stars and RS CVn systems, in which activity is enhanced by tidal interaction. Some evidence of a bimodal distribution in main sequence (MS) stars with Teff ≳ 5000 K is also found, which defines the two other regimes, without a clear gap in between. Stars with 3500 K ≲ Teff ≲ 5000 K are found to be either very active PMS stars or active MS stars with a unimodal distribution in chromospheric activity. A dramatic change in the activity distribution is found for Teff ≲ 3500 K, with a dominance of low activity stars close to the transition between partially- and fully convective stars and a rise in activity down into the fully convective regime.

Predicting Stellar Mass Accretion: An Optimized Echo State Network Approach in Time Series Modeling

Modeling the dynamics of the formation and evolution of protostellar disks as well as the history of stellar mass accretion typically involve the numerical solution of complex systems of coupled differential equations. The resulting mass accretion history of protostars is known to be highly episodic due to recurrent instabilities and also exhibits short timescale flickering. By leveraging the strong predictive abilities of neural networks, we extract some of the critical temporal dynamics experienced during the mass accretion including periods of instability. Particularly, we utilize a novel form of the Echo-State Neural Network (ESN), which has been shown to efficiently deal with data having inherent nonlinearity. We introduce the use of Optimized-ESN (Opt-ESN) to make model-independent time series forecasting of mass accretion rate in the evolution of protostellar disks. We apply the network to multiple hydrodynamic simulations with different initial conditions and exhibiting a variety of temporal dynamics to demonstrate the predictability of the Opt-ESN model. The model is trained on simulation data of $\sim 1-2$ Myr, and achieves predictions with a low normalized mean square error ($\sim 10^{-5}$ to $10^{-3}$) for forecasts ranging between 100 and 3800 yr. This result shows the promise of the application of machine learning based models to time-domain astronomy.

A Direct Measurement of Galaxy Major and Minor Merger Rates and Stellar Mass Accretion Histories at Z < 3 Using Galaxy Pairs in the REFINE Survey

We measure the role of major and minor mergers in forming the stellar masses of galaxies over redshifts 0 < z < 3 using a combination of ∼3.25 deg2 of the deepest ground-based near-infrared imaging taken to date (Ultra Deep Survey, Ultra-VISTA, and VIDEO) as part of the collated REFINE survey. We measure the pair fraction and merger fractions for galaxy mergers of different mass ratios, and quantify the merger rate with newly measured timescales derived from the Illustris simulation as a function of redshift and merger mass ratio. For a M * > 1011 M ⊙ selection, we find that over 0 < z < 3 major mergers with mass ratios greater than 1:4 occur times on average, while minor mergers down to ratios of 1:10 occur on average times per galaxy. We also quantify the role of major and minor mergers in galaxy formation, whereby the increase in mass due to major mergers is while minor mergers account for an increase of using a M * > 1011 M ⊙ selection. We thus find that major mergers add more stellar mass to galaxies than minor mergers over this epoch. Overall, mergers will more than double the mass of massive galaxies over this epoch when selecting by stellar mass. We however find a lower increase in stellar mass when selecting by a constant number density. Finally, we compare our results to simulations, finding that minor mergers are overpredicted in Illustris and in semi-analytical models, suggesting a mismatch between observations and theory in this fundamental aspect of galaxy assembly.

A semi-analytical model for the temporal evolution of the episodic disc-to-star accretion rate during star formation

ABSTRACT We develop a semi-analytical formalism for the determination of the evolution of the stellar mass accretion rate for specified density and velocity profiles that emerge from the runaway collapse of a prestellar cloud core. In the early phase, when the infall of matter from the surrounding envelope is substantial, the star accumulates mass primarily because of envelope-induced gravitational instability in a protostellar disc. In this phase, we model the envelope mass accretion rate from the isothermal free-fall collapse of a molecular cloud core. The disc gains mass from the envelope, and transports matter to the star via a disc accretion mechanism that includes episodic gravitational instability and mass accretion bursts according to the Toomre Q-criterion. In a later phase, mass is accreted on to the star due to gravitational torques within the spiral structures in the disc, in a manner that analytical theory suggests has a mass accretion rate ∝t−6/5. Our model provides a self-consistent evolution of the mass accretion rate by joining the spherical envelope accretion (dominant at the earlier stage) with the disc accretion (important at the later stage), and accounts for the presence of episodic accretion bursts at appropriate times. We show using a simple example that the burst mode can provide a good match to the observed distribution of bolometric luminosities. Our framework reproduces key elements of detailed numerical simulations of disc accretion and can aid in developing intuition about the basic physics as well as to compare theory with observations.

CO isotopolog line fluxes of viscously evolving disks

Context. Protoplanetary disks are thought to evolve viscously, where the disk mass – the reservoir available for planet formation – decreases over time as material is accreted onto the central star over a viscous timescale. Observations have shown a correlation between disk mass and the stellar mass accretion rate, as expected from viscous theory. However, this happens only when using the dust mass as a proxy of the disk mass; the gas mass inferred from CO isotopolog line fluxes, which should be a more direct measurement, shows no correlation with the stellar mass accretion rate. Aims. We investigate how 13CO and C18O J = 3−2 line fluxes, commonly used as gas mass tracers, change over time in a viscously evolving disk and use them together with gas disk sizes to provide diagnostics of viscous evolution. In addition, we aim to determine if the chemical conversion of CO through grain-surface chemistry combined with viscous evolution can explain the CO isotopolog observations of disks in Lupus. Methods. We ran a series of thermochemical DALI models of viscously evolving disks, where the initial disk mass is derived from observed stellar mass accretion rates. Results. While the disk mass, Mdisk, decreases over time, the 13CO and C18O J = 3−2 line fluxes instead increase over time due to their optically thick emitting regions growing in size as the disk expands viscously. The C18O 3–2 emission is optically thin throughout the disk for only for a subset of our models (M*≤ 0.2 M⊙ and αvisc ≥ 10−3, corresponding to Mdisk(t = 1 Myr) ≤ 10−3 M⊙). For these disks the integrated C18O flux decreases with time, similar to the disk mass. Observed 13CO and C18O 3–2 fluxes of the most massive disks (Mdisk ≳ 5 × 10−3 M⊙) in Lupus can be reproduced to within a factor of ~2 with viscously evolving disks in which CO is converted into other species through grain-surface chemistry with a moderate cosmic-ray ionization rate of ζcr ~ 10−17 s−1. The C18O 3–2 fluxes for the bulk of the disks in Lupus (with Mdisk ≲ 5 × 10−3 M⊙) can be reproduced to within a factor of ~2 by increasing ζcr to ~ 5 × 10−17−10−16 s−1, although explaining the stacked upper limits requires a lower average abundance than our models can produce. In addition, increasing ζcr cannot explain the observed 13CO fluxes for lower mass disks, which are more than an order of magnitude fainter than what is predicted. In our models the optically thick 13CO emission originates from a layer higher up in the disk (z∕r ~ 0.25−0.4) where photodissociation stops the conversion of CO into other species. Reconciling the 13CO fluxes of viscously evolving disks with the observations requires either efficient vertical mixing or low mass disks (Mdust ≲ 3 × 10−5 M⊙) being much thinner and/or smaller than their more massive counterparts. Conclusions. The 13CO model flux predominantly traces the disk size, but the C18O model flux traces the disk mass of our viscously evolving disk models if chemical conversion of CO is included. The discrepancy between the CO isotopolog line fluxes of viscously evolving disk models and the observations suggests that CO is efficiently vertically mixed or that low mass disks are smaller and/or colder than previously assumed.

Magnetic Braking of Accreting T Tauri Stars: Effects of Mass Accretion Rate, Rotation, and Dipolar Field Strength

The rotational evolution of accreting pre-main-sequence stars is influenced by its magnetic interaction with its surrounding circumstellar disk. Using the PLUTO code, we perform 2.5D magnetohydrodynamic, axisymmetric, time-dependent simulations of star-disk interaction---with an initial dipolar magnetic field structure, and a viscous and resistive accretion disk---in order to model the three mechanisms that contribute to the net stellar torque: accretion flow, stellar wind, and magnetospheric ejections (periodic inflation and reconnection events). We investigate how changes in the stellar magnetic field strength, rotation rate, and mass accretion rate (changing the initial disk density) affect the net stellar torque. All simulations are in a net spin-up regime. We fit semi-analytic functions for the three stellar torque contributions, allowing for the prediction of the net stellar torque for our parameter regime, and the possibility of investigating spin-evolution using 1D stellar evolution codes. The presence of an accretion disk appears to increase the efficiency of stellar torques compared to isolated stars, for cases with outflow rates much smaller than accretion rates, because the star-disk interaction opens more of the stellar magnetic flux compared to that from isolated stars. In our parameter regime, a stellar wind with a mass loss rate of $\approx 1 \%$ of the mass accretion rate is capable of extracting $\lesssim 50 \%$ of the accreting angular momentum. These simulations suggest that achieving spin-equilibrium in a representative T Tauri case within our parameter regime, e.g., BP Tau, would require a wind mass loss rate of $\approx 25\%$ of the mass accretion rate.

Magnetic Braking of Accreting T Tauri Stars: Effects of Mass Accretion Rate, Rotation, and Dipolar Field Strength

Abstract The rotational evolution of an accreting pre-main-sequence star is influenced by its magnetic interaction with its surrounding circumstellar disk. Using the PLUTO code, we perform 2.5D magnetohydrodynamic, axisymmetric, time-dependent simulations of star–disk interaction—with an initial dipolar magnetic field structure, and a viscous and resistive accretion disk—in order to model the three mechanisms that contribute to the net stellar torque: accretion flow, stellar wind, and magnetospheric ejections (periodic inflation and reconnection events). We investigate how changes in the stellar magnetic field strength, rotation rate, and mass accretion rate (changing the initial disk density) affect the net stellar torque. All simulations are in a net spin-up regime. We fit semi-analytic functions for the three stellar torque contributions, allowing for the prediction of the net stellar torque for our parameter regime, as well as the possibility of investigating spin evolution using 1D stellar evolution codes. The presence of an accretion disk appears to increase the efficiency of stellar torques compared to isolated stars, for cases with outflow rates much smaller than accretion rates, because the star–disk interaction opens more of the stellar magnetic flux compared to that from isolated stars. In our parameter regime, a stellar wind with a mass-loss rate of ≈1% of the mass accretion rate is capable of extracting ≲50% of the accreting angular momentum. These simulations suggest that achieving spin equilibrium in a representative T Tauri case within our parameter regime, e.g., BP Tau, would require a wind mass-loss rate of ≈25% of the mass accretion rate.

Migration jumps of planets in transition discs

Context. Transition discs form a special class of protoplanetary discs that are characterised by a deficiency of disc material close to the star. In a subgroup, inner holes in these discs can stretch out to a few tens of au while there is still mass accretion onto the central star observed at the same time. Aims. We analyse the proposition that this type of wide transition disc is generated by the interaction of the disc with a system of embedded planets. Methods. We performed two-dimensional hydrodynamics simulations of a flat disc. Different equations of state were used including locally isothermal models and more realistic cases that consider viscous heating, radiative cooling, and stellar heating. Two massive planets (with masses of between three and nine Jupiter masses) were embedded in the disc and their dynamical evolution due to disc–planet interaction was followed for over 100 000 yr. The simulations account for mass accretion onto the star and planets. We included models with parameters reminiscent of the system PDS 70. To assess the observability of features in our models we performed synthetic ALMA observations. Results. For systems with a more massive inner planet, there are phases where both planets migrate outward engaged in a 2:1 mean motion resonance via the Masset-Snellgrove mechanism. In sufficiently massive discs, the resulting formation of a vortex and the interaction with it can trigger rapid outward migration of the outer planet where its distance can increase by tens of au in a few thousand years. After another few thousand years, the outer planet rapidly migrates back inwards into resonance with the inner planet. We call this emerging composite phenomenon a migration jump. Outward migration and the migration jumps are accompanied by a high mass accretion rate onto the star. The synthetic images reveal numerous substructures depending on the type of dynamical behaviour. Conclusions. Our results suggest that the outward migration of two embedded planets is a prime candidate for the explanation of the observed high stellar mass accretion rate in wide transition discs. The models for PDS 70 indicate it is not currently undergoing a migration jump but might very well be in a phase of outward migration.

Variable Accretion onto Protoplanet Host Star PDS 70

Abstract The PDS 70 system has been subject to many studies in the past year following the discovery of two accreting planets in the gap of its circumstellar disk. Nevertheless, the mass accretion rate onto the star is still not well known. Here, we determined the stellar mass accretion rate and its variability based on Transiting Exoplanet Survey Satellite and High-Accuracy Radial velocity Planetary Searcher (HARPS) observations. The stellar light curve shows a strong signal with a 3.03 ± 0.06 days period, which we attribute to stellar rotation. Our analysis of the HARPS spectra shows a rotational velocity of , indicating that the inclination of the rotation axis is 50° ± 8°. This implies that the rotation axes of the star and its circumstellar disk are parallel within the measurement error. We apply magnetospheric accretion models to fit the profiles of the Hα line and derive mass accretion rates onto the star in the range of , varying over the rotation phase. The measured accretion rates are in agreement with those estimated from near-UV fluxes using accretion shock models. The derived accretion rates are higher than expected from the disk mass and planets’ properties for the low values of the viscous parameter α suggested by recent studies, potentially pointing to an additional mass reservoir in the inner disk to feed the accretion, such as a dead zone. We find that the He I λ10830 line shows a blueshifted absorption feature, indicative of a wind. The mass-loss rate estimated from the line depth is consistent with an accretion-driven inner disk MHD wind.

Resolving the FU Orionis System with ALMA: Interacting Twin Disks?

Abstract FU Orionis objects are low-mass pre-main sequence stars characterized by dramatic outbursts several magnitudes in brightness. These outbursts are linked to episodic accretion events in which stars gain a significant portion of their mass. The physical processes behind these accretion events are not yet well understood. The archetypal FU Ori system, FU Orionis, is composed of two young stars with detected gas and dust emission. The continuum emitting regions have not been resolved until now. Here, we present 1.3 mm observations of the FU Ori binary system using the Atacama Large Millimeter/submillimeter Array. The disks are resolved at 40 mas resolution. Radiative transfer modeling shows that the emission from FU Ori north (primary) is consistent with a dust disk with a characteristic radius of ∼11 au. The ratio between the major and minor axes shows that the inclination of the disk is ∼37°. FU Ori south is consistent with a dust disk of similar inclination and size. Assuming the binary orbit shares the same inclination angle as the disks, the deprojected distance between the north and south components is 0.″6, i.e., ∼250 au. Maps of 12CO emission show a complex kinematic environment with signature disk rotation at the location of the northern component, and also (to a lesser extent) for FU Ori south. The revised disk geometry allows us to update FU Ori accretion models, yielding a stellar mass and mass accretion rate of FU Ori north of 0.6 M ⊙ and 3.8 × 10−5 M ⊙ yr−1, respectively.

The accretion rates and mechanisms of Herbig Ae/Be stars

ABSTRACT This work presents a spectroscopic study of 163 Herbig Ae/Be stars. Amongst these, we present new data for 30 objects. Stellar parameters such as temperature, reddening, mass, luminosity, and age are homogeneously determined. Mass accretion rates are determined from $\rm H\alpha$ emission line measurements. Our data is complemented with the X-Shooter sample from previous studies and we update results using Gaia DR2 parallaxes giving a total of 78 objects with homogeneously determined stellar parameters and mass accretion rates. In addition, mass accretion rates of an additional 85 HAeBes are determined. We confirm previous findings that the mass accretion rate increases as a function of stellar mass, and the existence of a different slope for lower and higher mass stars, respectively. The mass where the slope changes is determined tobe $3.98^{+1.37}_{-0.94}\, \rm M_{\odot }$. We discuss this break in the context of different modes of disc accretion for low- and high-mass stars. Because of their similarities with T Tauri stars, we identify the accretion mechanism for the late-type Herbig stars with the Magnetospheric Accretion. The possibilities for the earlier-type stars are still open, we suggest the Boundary Layer accretion model may be a viable alternative. Finally, we investigated themass accretion–age relationship. Even using the superior Gaia based data, it proved hard to select a large enough sub-sample to remove the mass dependence in this relationship. Yet, it would appear that the mass accretion does decline with age as expected from basic theoretical considerations.

Mapping the Conditions for Hydrodynamic Instability on Steady-State Accretion Models of Protoplanetary Disks

Abstract Hydrodynamic instabilities in disks around young stars depend on the thermodynamic stratification of the disk and on the local rate of thermal relaxation. Here, we map the spatial extent of unstable regions for the Vertical Shear Instability (VSI), the Convective Overstability (COS), and the amplification of vortices via the Subcritical Baroclinic Instability (SBI). We use steady-state accretion disk models, including stellar irradiation, accretion heating, and radiative transfer. We determine the local radial and vertical stratification and thermal relaxation rate in the disk, which depends on the stellar mass, disk mass, and mass accretion rate. We find that passive regions of disks—that is, the midplane temperature dominated by irradiation—are COS unstable about one pressure scale height above the midplane and VSI unstable at radii &gt;10 au. Vortex amplification via SBI should operate in most parts of active and passive disks. For active parts of disks (midplane temperature determined by accretion power), COS can become active down to the midplane. The same is true for the VSI because of the vertically adiabatic stratification of an internally heated disk. If hydrodynamic instabilities or other nonideal MHD processes are able to create α-stresses (&gt;10−5) and released accretion energy leads to internal heating of the disk, hydrodynamic instabilities are likely to operate in significant parts of the planet-forming zones in disks around young stars, driving gas accretion and flow structure formation. Thus, hydrodynamic instabilities are viable candidates to explain the rings and vortices observed with the Atacama Large Millimeter/submillimeter Array and Very Large Telescope.

First hydrodynamics simulations of radiation forces and photoionization feedback in massive star formation

Aims. We present the first simulations of the formation and feedback of massive stars which account for radiation forces as well as photoionization feedback (along with protostellar outflows). In two different accretion scenarios modeled, we determine the relative strength of these feedback components and derive the size of the reservoir from which the forming stars gained their masses. Methods. We performed direct hydrodynamics simulations of the gravitational collapse of high-density mass reservoirs toward the formation of massive stars including self-gravity, stellar evolution, protostellar outflows, continuum radiation transport, photoionization, and the potential impact of ram pressure from large-scale gravitational infall. For direct comparison, we executed these simulations with and without the individual feedback components. Results. Protostellar outflows alone limit the stellar mass growth only in an accretion scenario with a finite mass reservoir; when including accretion and ram pressure from large scales (&gt;0.1 pc), protostellar outflows do not limit stellar mass growth at all. Photoionization and HII regions dominate the feedback ladder only at later times, after the star has already contracted down to the zero-age main sequence, and only on large scales. Specifically, photoionization yields a broadening of the bipolar outflow cavities and a reduction of the gravitational infall momentum by about 50%, but does not limit the stellar mass accretion. On the other hand, we find radiation forces restrain the gravitational infall toward the circumstellar disk, impact the gravito-centrifugal equilibrium at the outer edge of the disk, and eventually shut down stellar accretion completely. The most massive star formed in the simulations accreted 95 M⊙ before disk destruction; this mass was drawn-in from an accretion reservoir of ≈240 M⊙ and ≈0.24 pc in radius. Conclusions. In the regime of very massive stars, the final mass of these stars is controlled by their own radiation force feedback.

A single population of red globular clusters around the massive compact galaxy NGC 1277.

Massive galaxies are thought to form in two phases: an initial collapse of gas and giant burst of central star formation, followed by the later accretion of material that builds up their stellar and dark-matter haloes. The systems of globular clusters within such galaxies are believed to form in a similar manner. The initial central burst forms metal-rich (spectrally red) clusters, whereas more metal-poor (spectrally blue) clusters are brought in by the later accretion of less-massive satellites. This formation process is thought to result in the multimodal optical colour distributions that are seen in the globular cluster systems of massive galaxies. Here we report optical observations of the massive relic-galaxy candidate NGC 1277-a nearby, un-evolved example of a high-redshift 'red nugget' galaxy. We find that the optical colour distribution of the cluster system of NGC 1277 is unimodal and entirely red. This finding is in strong contrast to other galaxies of similar and larger stellar mass, the cluster systems of which always exhibit (and are generally dominated by) blue clusters. We argue that the colour distribution of the cluster system of NGC 1277 indicates that the galaxy has undergone little (if any) mass accretion after its initial collapse, and use simulations of possible merger histories to show that the stellar mass due to accretion is probably at most ten per cent of the total stellar mass of the galaxy. These results confirm that NGC 1277 is a genuine relic galaxy and demonstrate that blue clusters constitute an accreted population in present-day massive galaxies.

The effects of protostellar jet feedback on turbulent collapse

We present results of hydrodynamic simulations of massive star forming regions with and without protostellar jets. We show that jets change the normalization of the stellar mass accretion rate, but do not strongly affect the dynamics of star formation. In particular, $M_*(t) \propto f^2 (t-t_*)^2$ where $f = 1 - f_{\rm jet}$ is the fraction of mass accreted onto the protostar, $f_{\rm jet}$ is the fraction ejected by the jet, and $(t-t_*)^2$ is the time elapsed since the formation of the first star. The star formation efficiency is nonlinear in time. We find that jets have only a small effect (of order 25\%) on the accretion rate onto the protostellar disk (the "raw" accretion rate). We show that the small scale structure -- the radial density, velocity, and mass accretion profiles are very similar in the jet and no-jet cases. Finally, we show that the inclusion of jets does drive turbulence but only on small (parsec) scales.

The Class 0 Protostar BHR71: Herschel Observations and Dust Continuum Models

Abstract We use Herschel spectrophotometry of BHR71, an embedded Class 0 protostar, to provide new constraints on its physical properties. We detect 645 (non-unique) spectral lines among all spatial pixels. At least 61 different spectral lines originate from the central region. A CO rotational diagram analysis shows four excitation temperature components, 43, 197, 397, and 1057 K. Low-J CO lines trace the outflow while the high-J CO lines are centered on the infrared source. The low-excitation emission lines of trace the large-scale outflow, while the high-excitation emission lines trace a small-scale distribution around the equatorial plane. We model the envelope structure using the dust radiative transfer code, hyperion, incorporating rotational collapse, an outer static envelope, outflow cavity, and disk. The evolution of a rotating collapsing envelope can be constrained by the far-infrared/millimeter spectral energy distribution along with the azimuthally averaged radial intensity profile, and the structure of the outflow cavity plays a critical role at shorter wavelengths. Emission at 20–40 μm requires a cavity with a constant-density inner region and a power-law density outer region. The best-fit model has an envelope mass of 19 inside a radius of 0.315 pc and a central luminosity of 18.8 . The time since collapse began is 24,630–44,000 years, most likely around 36,000 years. The corresponding mass infall rate in the envelope (1.2 × 10−5 ) is comparable to the stellar mass accretion rate, while the mass-loss rate estimated from the CO outflow is 20% of the stellar mass accretion rate. We find no evidence for episodic accretion.

THE LINK BETWEEN RARE-EARTH PEAK FORMATION AND THE ASTROPHYSICAL SITE OF THE R PROCESS

ABSTRACT The primary astrophysical source of the rare-earth elements is the rapid neutron capture process (r process). The rare-earth peak that is seen in the solar r-process residuals has been proposed to originate as a pile-up of nuclei during the end of the r process. We introduce a new method utilizing Monte Carlo studies of nuclear masses in the rare-earth region, that includes self-consistently adjusting β-decay rates and neutron capture rates, to find the mass surfaces necessary for the formation of the rare-earth peak. We demonstrate our method with two types of astrophysical scenario, one corresponding to conditions typical of hot winds from core-collapse supernovae and stellar-mass accretion disks, and one corresponding to conditions typical of the ejection of the material from the tidal tails of neutron star mergers. In each type of astrophysical condition, this method successfully locates a region of enhanced stability in the mass surface that is responsible for the rare-earth peak. For each scenario, we find that the change in the mass surface has qualitatively different features, thus future measurements can shed light on the type of environment in which the r process occurred.