Scatter In Luminosity Research Articles

Galaxy–Absorber Association in the Epoch of Reionization: Galactic Population Luminosity Distribution for Different Absorbers at 10 ≥ z ≥ 5.5

How do galaxies of different luminosities contribute to the metal absorber populations of varying species and strength? We present our analysis of the predicted metal contributions from galaxies as observed in quasar absorption line spectra during the end of the epoch of reionization (10 ≥ z ≥ 5.5). This was done by implementing on-the-fly particle tracking into the latest Technicolor Dawn simulation and then linking C ii, C iv, Si ii, Si iv, O i, and Mg ii absorbers to host galaxies in postprocessing. We define the host galaxy luminosity distribution (HGLD) as the rest-frame ultraviolet luminosity distribution of galaxies contributing ions to an absorber, weighted by the fractional contribution, and compute its dependence on ion and absorber strength. The HGLD shape is predicted to be indistinguishable from the field luminosity function, indicating that there is no relationship between the absorber strength or ion and the luminosity of the dominant contributing galaxy. Switching from galaxy luminosity to stellar mass, the predicted host galaxy mass distributions (HGMDs) indicate that more-massive galaxies contribute a higher fraction of metal ions to absorbers of each species, with the HGMDs of stronger absorbers extending out to higher masses. We conclude that the fraction of absorbing metal ions contributed by galaxies increases weakly with stellar mass, but the scatter in luminosity at fixed stellar mass obscures this relationship. For the same reason, we predict that observational analyses of the absorber–galaxy relationship will uncover stronger trends with stellar mass than with luminosity.

On the incidence of episodic accretion in Class I YSOs from VVV

ABSTRACT Episodic accretion is one of the competing models to explain the observed luminosity spread in young stellar clusters. These short-lived high accretion events could also have a strong impact on planet formation. Observations of high-amplitude variability in young stellar objects (YSOs) due to large changes in the accretion rate provide direct observational evidence for episodic accretion. However, there are still uncertainties in the frequency of these events and if episodic accretion is universal among YSOs. To determine the frequency of outbursts in Class I YSOs, we built a large and robust sample of objects at this evolutionary stage, and searched for high-amplitude near-infrared (ΔKS > 2 mag) variability in the VIRAC2 database of the Vista Variables in the Via Lactea survey. By complementing with near-IR (2MASS and DENIS) and mid-IR (WISE/Neo-WISE) data, we find that from ∼7000 Class I YSOs, 97 objects can be classified as eruptive variable YSOs. The duration of the outbursts vary from a few months to longer than 9 yr, and cover a similar range of amplitudes. Values of ΔKS > 5 mag, however, are only observed in outbursts with duration longer than 9 yr. When considering different effects of completeness and contamination, we estimate that the incidence of episodic accretion in Class I YSOs is between 2 and 3 per cent. Finally, we determine a recurrence time-scale of long-term outbursts (a.k.a FUors) of $\tau =1.75^{+1.12}_{-0.87}$ kyr. The latter value agrees with previous estimates and is in line with the expectations of higher frequency of FUor outbursts during younger stages of evolution.

Don’t Believe the Hype(r): The Yellow Supergiants of Westerlund 1

Yellow hypergiants (YHGs) are often presumed to represent a transitional post-red supergiant (RSG) phase for stars ∼30–40 M ⊙. Here we present visual-wavelength echelle spectra of six YHG candidates in the Galactic cluster Westerlund 1, and we compare them to known YHGs, IRC +10420 and Hen 3-1979. We find that the six YHG candidates do not exhibit any metal emission lines, nor do they show strong Hα emission, and as such do not meet the criteria necessary to be classified as YHGs. In conjunction with their moderate luminosities of = 4.7–5.4 estimated from optical/infrared photometry, we suggest instead that they are normal yellow supergiants (YSGs) with more modest initial masses around 15–20 M ⊙. This adds additional support to the hypothesis that Wd1 is a multiage cluster with an older age than previously assumed, and is not a ∼5 Myr old cluster caught at a very specific transitional point when single-star evolution might yield Wolf–Rayet stars, luminous blue variables, RSGs, and YHGs in the same cluster. Nevertheless, the population of YSGs in Wd1 is very unusual, with YSGs outnumbering RSGs, but with both spanning a large luminosity range. Here, we discuss evolutionary scenarios that might have led to the high fraction of YSGs. The number of YSGs and their significant luminosity spread cannot be explained by simple population synthesis models with single or binary stars. Even with multiple ages or a large age spread, the high YSG/RSG ratio remains problematic. We suggest instead that the objects may experience a prolonged YSG phase due to evolution in triple systems.

Are Brightest Cluster Galaxies Special?

A longstanding question when studying Brightest Cluster Galaxies (BCGs) and their evolution is whether their luminosities are statistical extremes of the luminosity distribution of other cluster member galaxies, or whether their luminosities follow a different,“special” distribution. To explore this question, we use the two Tremaine–Richstone (TR) statistics that sometimes lead to different conclusions about BCG specialness. We run Monte Carlo simulations of galaxy clusters with varying galaxy luminosity functions and BCG properties, and explore the sensitivity of the TR statistics to the adjusted BCG properties. We identify cases where the two statistics lead to different conclusions. The first of the TR statistics, which relies on the spread of BCG luminosities, is more reliable for determining specialness than the second, which utilizes the spread in the luminosity gap between the BCG and second ranked galaxy.

Revisiting the [C ii] 158 μm line-intensity mapping power spectrum from the EoR using non-uniform line-luminosity scatter

Abstract Detecting the line-intensity mapping (LIM) signal from the galaxies of the epoch of reionization is an emerging tool to constrain their role in reionization. Ongoing and upcoming experiments target the signal fluctuations across the sky to reveal statistical and astrophysical properties of these galaxies via signal statistics, e.g. the power spectrum. Here, we revisit the [C ii]$_{158 \, \mu \text{m}}$ LIM power spectrum under non-uniform line–luminosity scatter, which has a halo-mass variation of statistical properties. Line–luminosity scatter from a cosmological hydrodynamic and radiative transfer simulation of galaxies at $z$ = 6 is considered in this study. We test the robustness of different model frameworks that interpret the impact of the line-luminosity scatter on the signal statistics. We use a simple power-law model to fit the scatter and demonstrate that the mean luminosity–halo mass correlation fit cannot preserve the mean intensity of the LIM signal (hence the clustering power spectrum) under non-uniform scatter. In our case, the mean intensity changes by ∼48 per cent compared to the mean correlation fit in contrast to the general case with semi-analytical scatter. However, we find that the prediction for the mean intensity from the most-probable fit can be modelled robustly, considering the generalized and more realistic non-uniform scatter. We also explore the possibility of diminishing luminosity bias under non-uniform scatter, affecting the clustering power spectrum, although this phenomenon might not be statistically significant. Therefore, we should adopt appropriate approaches that can consistently interpret the LIM power spectrum from observations.

Differences between the globular cluster systems of the Virgo and Fornax galaxy clusters

ABSTRACT It is well-known that Globular cluster systems are different among galaxies. Here, we test to which degree these differences remain on the scale of galaxy clusters by comparing the globular clusters (GCs) in optical surveys of the Virgo galaxy cluster (ACSVCS) and the ACS Fornax galaxy cluster (ACSFCS) in Kolmogorov–Smirnoff Tests. Both surveys were obtained with the Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope, and contain thousands of GCs in dozens of galaxies each. Also, well over 100 point sources in the Chandra X-ray Observatory source catalogue were attributed to the GCs in both optical catalogues, and interpreted as low-mass X-ray binaries (LMXBs). Thus, the optical and X-ray data are as uniform as possible. Our main findings are as follows: (1) The spread in luminosities and half-light radii is larger in the ACSVCS than in the ACSFCS. (2) The ratio between the half-light radii for the F475W-passband and the F850LP-passband is on average smaller in the ACSVCS. (3) The distribution of the LMXBs with the luminosity of the GCs is different between both surveys. These findings are significant. The first finding could be a consequence of a wider spread in the distances of the GCs in the ACSVCS, but the others must have internal reasons in the GCs. Thus, the GC systems are also different on a galaxy cluster scale.

The Diversity of Environments around Luminous Quasars at Redshift z ∼ 6

Abstract Significant clustering around the rarest luminous quasars is a feature predicted by dark matter theory combined with number density matching arguments. However, this expectation is not reflected by observations of quasars residing in a diverse range of environments. Here we assess the tension in the diverse clustering of visible i-band dropout galaxies around luminous z ∼ 6 quasars. Our approach uses a simple empirical method to derive the median luminosity–to–halo mass relation, L c (M h ), for both quasars and galaxies under the assumption of lognormal luminosity scatter, Σ Q and Σ G . We show that higher Σ Q reduces the average halo mass hosting a quasar of a given luminosity, thus introducing at least a partial reversion to the mean in the number count distribution of nearby Lyman-break galaxies. We generate a large sample of mock Hubble Space Telescope fields of view centered across rare z ∼ 6 quasars by resampling pencil beams traced through the dark matter component of the BlueTides cosmological simulation. We find that diverse quasar environments are expected for Σ Q > 0.4, consistent with numerous observations and theoretical studies. However, we note that the average number of galaxies around the central quasar is primarily driven by galaxy evolutionary processes in neighboring halos, as embodied by our parameter Σ G , instead of a difference in the large-scale structure around the central quasar host, embodied by Σ Q . We conclude that models with Σ G > 0.3 are consistent with current observational constraints on high-z quasars, and that such a value is comparable to the scatter estimated from hydrodynamical simulations of galaxy formation.

Jet–Accretion System in the Nearby mJy Radio Galaxies

Abstract It is generally thought that FRII radio galaxies host thin optically thick disks, while FRIs are powered by advection-dominated accretion flows. Sources with an efficient engine are optically classified as high-excitation radio galaxies (HERGs) and those with an inefficient motor as low-excitation radio galaxies (LERGs). Recently, the study of radio galaxies down to mJy fluxes has cast serious doubts on the LERG-FRI and HERG-FRII correspondence, revealing that many LERGs show FRII radio morphologies. The FR catalogs recently compiled by Capetti et al. and Baldi et al. have allowed us to explore this issue in the local (z ≤ 0.15) mJy universe. Our statistical study shows that the majority of nearby mJy objects are in a late stage of their life. FRII-LERGs appear more similar to the old FRI-LERGs than to the young FRII-HERGs. FRII-LERGs may be aged HERGs that, having exhausted their cold fuel, have changed their accretion regime or are a separate LERG class particularly efficient in launching jets. Exploiting the empirical relations that convert L [O III] and L 1.4 GHz into accretion power and jet kinetic power, respectively, we observed that LERGs with similar masses and accretion rates seem to expel jets of different powers. We speculate that intrinsic differences related to the black hole properties (spin and magnetic field at its horizon) can determine the observed spread in jet luminosity. In this view, FRII-LERGs should have the fastest spinning black holes and/or the most intense magnetic fluxes. On the contrary, compact LERGs (i.e., FR0s) should host extremely slow black holes and/or weak magnetic fields.

An Improved Hertzsprung–Russell Diagram for the Orion Trapezium Cluster

Abstract We present a study of the Trapezium cluster in Orion. We analyze flux-calibrated Very Large Telescope/Multi-Unit Spectroscopic Explorer spectra of 361 stars to simultaneously measure the spectral types, reddening, and the optical veiling due to accretion. We find that the extinction law from Cardelli et al. with a total-to-selective extinction value of R V = 5.5 is more suitable for this cluster. For 68% of the sample the new spectral types are consistent with literature spectral types within two subclasses but, as expected, we derive systematically later types than the literature by one to two subclasses for the sources with significant accretion levels. Here we present an improved Hertzsprung–Russell (H-R) diagram of the Trapezium cluster, in which the contamination by optical veiling on spectral types and stellar luminosities has been properly removed. A comparison of the locations of the stars in the H-R diagram with the non-magnetic and magnetic pre-main-sequence evolutionary tracks indicates an age of 1–2 Myr. The magnetic pre-main-sequence evolutionary tracks can better explain the luminosities of the low-mass stars. In the H-R diagram, the cluster exhibits a large luminosity spread (σ(Log L ⋆/L ⊙) ∼ 0.3). By collecting a sample of 14 clusters/groups with different ages, we find that the luminosity spread tends to be constant (σ(Log L ⋆/L ⊙) ∼ 0.2–0.25) after 2 Myr, which suggests that age spread is not the main cause of the luminosity spread. There are ∼0.1 dex larger luminosity spreads for the younger clusters, e.g., the Trapezium cluster, than the older clusters, which can be explained by the starspots, accretion history, and circumstellar disk orientations.

Stochastic Processes as the Origin of the Double Power-law Shape of the Quasar Luminosity Function

Abstract The quasar luminosity function (QLF) offers insight into the early coevolution of black holes and galaxies. It has been characterized observationally up to redshift z ∼ 6 with clear evidence of a double power-law shape, in contrast to the Schechter-like form of the underlying dark-matter halo mass function. We investigate a physical origin for the difference in these distributions by considering the impact of stochasticity induced by the processes that determine the quasar luminosity for a given host halo and redshift. We employ a conditional luminosity function and construct the relation between median quasar magnitude versus halo mass with log-normal in luminosity scatter Σ, and duty-cycle ϵ DC, and focus on high redshift z ≳ 4. We show that, in order to reproduce the observed QLF, the Σ = 0 abundance matching requires all of the brightest quasars to be hosted in the rarest most massive dark-matter halos (with an increasing in halo mass). Conversely, for Σ > 0 the brightest quasars can be overluminous outliers hosted in relatively common dark-matter halos. In this case, the median quasar magnitude versus halo mass relation, M UV,c, flattens at the high-end, as expected in self-regulated growth due to feedback. We sample the parameter space of Σ and ϵ DC and show that M UV,c flattens above for . Models with ϵ DC ∼ 1 instead require a high mass threshold close to . We investigate the impact of ϵ DC and Σ on measurements of clustering and find there is no luminosity dependence on clustering for Σ > 0.3, consistent with recent observations from Subaru HSC.

On the Gaia DR2 distances for Galactic luminous blue variables

ABSTRACT We examine parallaxes and distances for Galactic luminous blue variables (LBVs) in the Gaia second data release (DR2). The sample includes 11 LBVs and 14 LBV candidates. For about half of the sample, DR2 distances are either similar to commonly adopted literature values, or the DR2 values have large uncertainties. For the rest, reliable DR2 distances differ significantly from values in the literature, and in most cases the Gaia DR2 distance is smaller. Two key results are that the S Doradus instability strip may not be as clearly defined as previously thought, and that there exists a population of LBVs at relatively low luminosities. LBVs seem to occupy a wide swath from the end of the main sequence at the blue edge to ∼8000 K at the red side, with a spread in luminosity reaching as low as log(L/L⊙) ≈ 4.5. The lower-luminosity group corresponds to effective single-star initial masses of 10–25 M⊙, and includes objects that have been considered as confirmed LBVs. We discuss implications for LBVs including (1) their instability and origin in binary evolution, (2) connections to some supernova (SN) impostors such as the class of SN 2008S-like objects, and (3) LBVs that may be progenitors of SNe with dense circumstellar material across a wide initial mass range. Although some of the Gaia DR2 distances for LBVs have large uncertainty, this represents the most direct and consistent set of Galactic LBV distance estimates available in the literature.

Testing asteroseismology with Gaia DR2: hierarchical models of the Red Clump

Asteroseismology provides fundamental stellar parameters independent of distance, but subject to systematics under calibration. Gaia DR2 has provided parallaxes for a billion stars, which are offset by a parallax zero-point. Red Clump (RC) stars have a narrow spread in luminosity, thus functioning as standard candles to calibrate these systematics. This work measures how the magnitude and spread of the RC in the Kepler field are affected by changes to temperature and scaling relations for seismology, and changes to the parallax zero-point for Gaia. We use a sample of 5576 RC stars classified through asteroseismology. We apply hierarchical Bayesian latent variable models, finding the population level properties of the RC with seismology, and use those as priors on Gaia parallaxes to find the parallax zero-point offset. We then find the position of the RC using published values for the zero-point. We find a seismic temperature insensitive spread of the RC of ~0.03 mag in the 2MASS K band and a larger and slightly temperature-dependent spread of ~0.13 mag in the Gaia G band. This intrinsic dispersion in the K band provides a distance precision of ~1% for RC stars. Using Gaia data alone, we find a mean zero-point of -41 $\pm$ 10 $\mu$as. This offset yields RC absolute magnitudes of -1.634 $\pm$ 0.018 in K and 0.546 $\pm$ 0.016 in G. Obtaining these same values through seismology would require a global temperature shift of ~-70 K, which is compatible with known systematics in spectroscopy.

The Effect of Spots on the Luminosity Spread of the Pleiades

Abstract Cool spots on the surface of magnetically active stars modulate their observed brightnesses and temperatures, thereby affecting the stellar locus on the H-R diagram. Recent high-precision space-based photometric surveys reveal the rotational modulation from spots on stars in young clusters, including K2 monitoring of the 125 Myr old Pleiades cluster. However, light curves reveal only the asymmetries in the visible spot distributions rather than the total sizes of spots on stellar surfaces, which causes a discrepancy between the spot coverage measured by photometric and spectroscopic observations. In this paper, we simulate photometric variability introduced by randomly distributed starspots on a 125 Myr old coeval cluster. Our simulation results show that randomly distributed small spots on the stellar surface would explain the discrepancy that the photometric observations only reveal 10%–40% of the spot coverage measured by spectra. The colors and luminosities of photospheres are modeled for a range of photospheric temperature, spot coverage, and spot temperature. The colors and luminosities of a simulated population are then compared to the luminosity spread of Pleiades members, excluding the 25% of candidates that are identified as non-members with Gaia DR2 astrometry. The observed luminosities of Pleiades members have a standard deviation of 0.05 dex, which could be entirely explained by spots with a star-to-star standard deviation of spot coverage of 10%, but with an average coverage area that is not well constrained.

Type Ia supernova Hubble diagram with near-infrared and optical observations

Context.Type Ia Supernovae (SNe Ia) have been used as standardizable candles in the optical wavelengths to measure distances with an accuracy of ~7% out to redshiftz~ 1.5. There is evidence that in the near-infrared (NIR) wavelengths SNe Ia are even better standard candles, however, NIR observations are much more time-consuming.Aims.We aim to test whether the NIR peak magnitudes could be accurately estimated with only a single observation obtained close to maximum light, provided that the time ofBband maximum, theB−Vcolor at maximum and the optical stretch parameter are known.Methods.We present multi-epochUBV RIand single-epochJandHphotometric observations of 16 SNe Ia in the redshift rangez= 0.037 − 0.183, doubling the leverage of the current SN Ia NIR Hubble diagram and the number of SNe beyond redshift 0.04. This sample was analyzed together with 102 NIR and 458 optical light curves (LCs) of normal SNe Ia from the literature.Results.The analysis of 45 NIR LCs with well-sampled first maximum shows that a single template accurately describes the LCs if its time axis is stretched with the optical stretch parameter. This allows us to estimate the peak NIR magnitudes of SNe with only few observations obtained within ten days fromB-band maximum. The NIR Hubble residuals show weak correlation with ΔM15and the color excessE(B−V), and for the first time we report a potential dependence on theJmax−Hmaxcolor. With these corrections, the intrinsic NIR luminosity scatter of SNe Ia is estimated to be ~0.10 mag, which is smaller than what can be derived for a similarly heterogeneous sample at optical wavelengths. Analysis of both NIR and optical data shows that the dust extinction in the host galaxies corresponds to a lowRV≃ 1.8–1.9.Conclusions.We conclude that SNe Ia are at least as good standard candles in the NIR as in the optical and are potentially less affected by systematic uncertainties. We extended the NIR SN Ia Hubble diagram to its nonlinear part atz~ 0.2 and confirmed that it is feasible to accomplish this result with very modest sampling of the NIR LCs, if complemented by well-sampled optical LCs. With future facilities it will be possible to extend the NIR Hubble diagram beyond redshiftz≃ 1, and our results suggest that the most efficient way to achieve this would be to obtain a single observation close to the NIR maximum.

Explaining the luminosity spread in young clusters: proto and pre-main sequence stellar evolution in a molecular cloud environment

Hertzsprung-Russell diagrams of star forming regions show a large luminosity spread. This is incompatible with well-defined isochrones based on classic non-accreting protostellar evo- lution models. Protostars do not evolve in isolation of their environment, but grow through accretion of gas. In addition, while an age can be defined for a star forming region, the ages of individual stars in the region will vary. We show how the combined effect of a protostellar age spread, a consequence of sustained star formation in the molecular cloud, and time-varying protostellar accretion for individual protostars can explain the observed luminosity spread. We use a global MHD simulation including a sub-scale sink particle model of a star forming region to follow the accretion process of each star. The accretion profiles are used to compute stellar evolution models for each star, incorporating a model of how the accretion energy is distributed to the disk, radiated away at the accretion shock, or incorporated into the outer layers of the protostar. Using a modelled cluster age of 5 Myr we naturally reproduce the lumi- nosity spread and find good agreement with observations of the Collinder 69 cluster, and the Orion Nebular Cluster. It is shown how stars in binary and multiple systems can be externally forced creating recurrent episodic accretion events. We find that in a realistic global molecular cloud model massive stars build up mass over relatively long time-scales. This leads to an important conceptual change compared to the classic picture of non-accreting stellar evolution segmented in to low-mass Hayashi tracks and high-mass Henyey tracks.

The Clustering of Luminous Red Galaxies at z ∼ 0.7 from EBOSS and BOSS Data

Abstract We present the first scientific results from the luminous red galaxy (LRG) sample of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) combined with the high-redshift galaxies of the previous BOSS sample. We measure the small- and intermediate-scale clustering from a sample of more than 97,000 galaxies in the redshift range . We interpret these measurements in the framework of the Halo Occupation Distribution. The bias of this sample of LRGs is 2.30 ± 0.03, with a satellite fraction of 13% ± 3% and a mean halo mass of . These results are consistent with expectations, demonstrating that these LRGs will be reliable tracers of large-scale structure at . The galaxy bias implies a scatter of luminosity at fixed halo mass, , of 0.19 dex. Using the clustering of massive galaxies from BOSS CMASS, BOSS LOWZ, and SDSS, we find that is consistent with observations over the full redshift range that these samples cover. The addition of eBOSS to previous surveys allows the investigation of the evolution of massive galaxies over the past ∼7 Gyr.

Age Spreads and the Temperature Dependence of Age Estimates in Upper Sco

Abstract Past estimates for the age of the Upper Sco Association are typically 11–13 Myr for intermediate-mass stars and 4–5 Myr for low-mass stars. In this study, we simulate populations of young stars to investigate whether this apparent dependence of estimated age on spectral type may be explained by the star formation history of the association. Solar and intermediate mass stars begin their pre-main sequence evolution on the Hayashi track, with fully convective interiors and cool photospheres. Intermediate-mass stars quickly heat up and transition onto the radiative Henyey track. As a consequence, for clusters in which star formation occurs on a timescale similar to that of the transition from a convective to a radiative interior, discrepancies in ages will arise when ages are calculated as a function of temperature instead of mass. Simple simulations of a cluster with constant star formation over several Myr may explain about half of the difference in inferred ages versus photospheric temperature; speculative constructions that consist of a constant star formation followed by a large supernova-driven burst could fully explain the differences, including those between F and G stars where evolutionary tracks may be more accurate. The age spreads of low-mass stars predicted from these prescriptions for star formation are consistent with the observed luminosity spread of Upper Sco. The conclusion that a lengthy star formation history will yield a temperature dependence in ages is expected from the basic physics of pre-main sequence evolution, and is qualitatively robust to the large uncertainties in pre-main sequence evolutionary models.

The Herschel Orion Protostar Survey: Luminosity and Envelope Evolution

Abstract The Herschel Orion Protostar Survey obtained well-sampled 1.2–870 μm spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of the bolometric temperature is roughly flat; 29% of the protostars are in Class 0. The median luminosity decreases by a factor of four with increasing bolometric temperature; consequently, the Class 0 protostars are systematically brighter than the Class I protostars, with a median luminosity of 2.3 as opposed to 0.87 . At a given bolometric temperature, the scatter in luminosities is three orders of magnitude. Using fits to the SEDs, we analyze how the luminosities corrected for inclination and foreground reddening relate to the mass in the inner 2500 au of the best-fit model envelopes. The histogram of the envelope mass is roughly flat, while the median-corrected luminosity peaks at 15 for young envelopes and falls to 1.7 for late-stage protostars with remnant envelopes. The spread in luminosity at each envelope mass is three orders of magnitude. Envelope masses that decline exponentially with time explain the flat mass histogram and the decrease in luminosity, while the formation of a range of stellar masses explains the dispersion in luminosity.

Revisiting the pre-main-sequence evolution of stars

Recent theoretical work has shown that the pre-main-sequence (PMS) evolution of stars is much more complex than previously envisioned. Instead of the traditional steady, one-dimensional solution, accretion may be episodic and not necessarily symmetrical, thereby affecting the energy deposited inside the star and its interior structure. Given this new framework, we want to understand what controls the evolution of accreting stars. We use the MESA stellar evolution code with various sets of conditions. In particular, we account for the (unknown) efficiency of accretion in burying gravitational energy into the protostar through a parameter, $\xi$, and we vary the amount of deuterium present. We confirm the findings of previous works that the evolution changes significantly with the amount of energy that is lost during accretion. We find that deuterium burning also regulates the PMS evolution. In the low-entropy accretion scenario, the evolutionary tracks in the H-R diagram are significantly different from the classical tracks and are sensitive to the deuterium content. A comparison of theoretical evolutionary tracks and observations allows us to exclude some cold accretion models ($\xi\sim 0$) with low deuterium abundances. We confirm that the luminosity spread seen in clusters can be explained by models with a somewhat inefficient injection of accretion heat. The resulting evolutionary tracks then become sensitive to the accretion heat efficiency, initial core entropy, and deuterium content. In this context, we predict that clusters with a higher D/H ratio should have less scatter in luminosity than clusters with a smaller D/H. Future work on this issue should include radiation-hydrodynamic simulations to determine the efficiency of accretion heating and further observations to investigate the deuterium content in star-forming regions. (abbrev.)

Self-consistent evolution of accreting low-mass stars and brown dwarfs

We present self-consistent calculations coupling numerical hydrodynamics simulations of collapsing pre-stellar cores and stellar evolution models of accreting objects. We analyse the main impact of consistent accretion history on the evolution and lithium depletion of young low-mass stars and brown dwarfs. These consistent models confirm the generation of a luminosity spread in Herzsprung-Russell diagrams at ages ~1−10 Myr. They also confirm that early accretion can produce objects with abnormal Li depletion, as found in a previous study that was based on arbitrary accretion rates. The results strengthen that objects with anomalously high level of Li depletion in young clusters should be extremely rare. We also find that early phases of burst accretion can produce coeval models of similar mass with a range of different Li surface abundances, and in particular with Li-excess compared to the predictions of non-accreting counterparts. This result is due to a subtle competition between the effect of burst accretion and its impact on the central stellar temperature, the growth of the stellar radiative core and the accretion of fresh Li from the accretion disk. Only consistent models could reveal such a subtle combination of effects. This new result could explain the recent, puzzling observations of Li-excess of fast rotators in the young cluster NGC 2264. Present self-consistent accreting models are available in electronic form.