Final Stellar Mass Research Articles

From the CMF to the IMF: beyond the core-collapse model

ABSTRACT Observations have indicated that the pre-stellar core mass function (CMF) is similar to the stellar initial mass function (IMF), except for an offset towards larger masses. This has led to the idea that there is a one-to-one relation between cores and stars, such that the whole stellar mass reservoir is contained in a gravitationally bound pre-stellar core, as postulated by the core-collapse model, and assumed in recent theoretical models of the stellar IMF. We test the validity of this assumption by comparing the final mass of stars with the mass of their progenitor cores in a high-resolution star formation simulation that generates a realistic IMF under physical condition characteristic of observed molecular clouds. Using a definition of bound cores similar to previous works we obtain a CMF that converges with increasing numerical resolution. We find that the CMF and the IMF are closely related in a statistical sense only; for any individual star there is only a weak correlation between the progenitor core mass and the final stellar mass. In particular, for high-mass stars only a small fraction of the final stellar mass comes from the progenitor core, and even for low-mass stars the fraction is highly variable, with a median fraction of only about 50 per cent. We conclude that the core-collapse scenario and related models for the origin of the IMF are incomplete. We also show that competitive accretion is not a viable alternative.

The Core Mass Function in the Orion Nebula Cluster Region: What Determines the Final Stellar Masses?

Abstract Applying dendrogram analysis to the CARMA-NRO C18O (J = 1–0) data having an angular resolution of ∼8″, we identified 692 dense cores in the Orion Nebula Cluster region. Using this core sample, we compare the core and initial stellar mass functions in the same area to quantify the step from cores to stars. About 22% of the identified cores are gravitationally bound. The derived core mass function (CMF) for starless cores has a slope similar to Salpeter’s stellar initial mass function (IMF) for the mass range above 1 M ⊙, consistent with previous studies. Our CMF has a peak at a subsolar mass of ∼0.1 M ⊙, which is comparable to the peak mass of the IMF derived in the same area. We also find that the current star formation rate is consistent with the picture in which stars are born only from self-gravitating starless cores. However, the cores must gain additional gas from the surroundings to reproduce the current IMF (e.g., its slope and peak mass), because the core mass cannot be accreted onto the star with 100% efficiency. Thus, the mass accretion from the surroundings may play a crucial role in determining the final stellar masses of stars.

The influence of infall on the properties of protoplanetary discs

Context. The properties of protoplanetary discs determine the conditions for planet formation. In addition, planets can already form during the early stages of infall. Aims. We constrain physical quantities such as the mass, radius, lifetime, and gravitational stability of protoplanetary discs by studying their evolution from formation to dispersal. Methods. We perform a population synthesis of protoplanetary discs with a total of 50 000 simulations using a 1D vertically integrated viscous evolution code, studying a parameter space of final stellar mass from 0.05 to 5 M⊙. Each star-and-disc system is set up shortly after the formation of the protostar and fed by infalling material from the parent molecular cloud core. Initial conditions and infall locations are chosen based on the results from a radiation-hydrodynamic population synthesis of circumstellar discs. We also consider a different infall prescription based on a magnetohydrodynamic (MHD) collapse simulation in order to assess the influence of magnetic fields on disc formation. The duration of the infall phase is chosen to produce a stellar mass distribution in agreement with the observationally determined stellar initial mass function. Results. We find that protoplanetary discs are very massive early in their lives. When averaged over the entire stellar population, the discs have masses of ~0.3 and 0.1 M⊙ for systems based on hydrodynamic or MHD initial conditions, respectively. In systems characterised by a final stellar mass ~1 M⊙, we find disc masses of ~0.7 M⊙ for the “hydro” case and ~0.2 M⊙ for the “MHD” case at the end of the infall phase. Furthermore, the inferred total disc lifetimes are long, ≈5–7 Myr on average. This is despite our choice of a high value of 10−2 for the background viscosity α-parameter. In addition, we find that fragmentation is common in systems that are simulated using hydrodynamic cloud collapse, with more fragments of larger mass formed in more massive systems. In contrast, if disc formation is limited by magnetic fields, fragmentation may be suppressed entirely. Conclusions. Our work draws a picture quite different from the one often assumed in planet formation studies: protoplanetary discs are more massive and live longer. This means that more mass is available for planet formation. Additionally, when fragmentation occurs, it can affect the disc’s evolution by transporting large amounts of mass radially. We suggest that the early phases in the lives of protoplanetary discs should be included in studies of planet formation. Furthermore, the evolution of the central star, including its accretion history, should be taken into account when comparing theoretical predictions of disc lifetimes with observations.

General-relativistic instability in hylotropic supermassive stars

Context. The formation of supermassive black holes by direct collapse would imply the existence of supermassive stars (SMSs) and their collapse through the general-relativistic (GR) instability into massive black hole seeds. However, the final mass of SMSs is weakly constrained by existing models, in spite of the importance of this value for the consistency of the direct collapse scenario. Aims. We estimate the final masses of spherical SMSs within the whole parameter space that is relevant to these objects. Methods. We built analytical stellar structures (hylotropes) that mimic existing numerical SMS models, accounting for full stellar evolution with rapid accretion. From these hydrostatic structures, we determine ab initio the conditions for GR instability and compare the results with the predictions for full stellar evolution. Results. We show that hylotropic models predict the onset of GR instability with a high level of precision. The mass of the convective core appears as a decisive quantity. The lower it is, the larger the total mass required for GR instability. The typical conditions for GR instability feature a total mass of ≳105 M⊙ with a core mass of ≳104 M⊙. If the core mass remains below 104 M⊙, total masses in excess of 106 − 107 M⊙ can be reached. Conclusions. Our results confirm that spherical SMSs forming in primordial, atomically cooled haloes collapse at masses below 500 000 M⊙. On the other hand, accretion rates in excess of 1000 M⊙ yr−1, leading to final stellar masses of ≳106 M⊙, are required for massive black hole formation in metal-rich gas. Thus, the different channels of direct collapse imply distinct final masses for the progenitor of the black hole seed.

EDGE: a new approach to suppressing numerical diffusion in adaptive mesh simulations of galaxy formation

ABSTRACT We introduce a new method to mitigate numerical diffusion in adaptive mesh refinement (AMR) simulations of cosmological galaxy formation, and study its impact on a simulated dwarf galaxy as part of the ‘EDGE’ project. The target galaxy has a maximum circular velocity of $21\, \mathrm{km}\, \mathrm{s}^{-1}$ but evolves in a region that is moving at up to $90\, \mathrm{km}\, \mathrm{s}^{-1}$ relative to the hydrodynamic grid. In the absence of any mitigation, diffusion softens the filaments feeding our galaxy. As a result, gas is unphysically held in the circumgalactic medium around the galaxy for $320\, \mathrm{Myr}$, delaying the onset of star formation until cooling and collapse eventually triggers an initial starburst at z = 9. Using genetic modification, we produce ‘velocity-zeroed’ initial conditions in which the grid-relative streaming is strongly suppressed; by design, the change does not significantly modify the large-scale structure or dark matter accretion history. The resulting simulation recovers a more physical, gradual onset of star formation starting at z = 17. While the final stellar masses are nearly consistent ($4.8 \times 10^6\, \mathrm{M}_{\odot }$ and $4.4\times 10^6\, \mathrm{M}_{\odot }$ for unmodified and velocity-zeroed, respectively), the dynamical and morphological structure of the z = 0 dwarf galaxies are markedly different due to the contrasting histories. Our approach to diffusion suppression is suitable for any AMR zoom cosmological galaxy formation simulations, and is especially recommended for those of small galaxies at high redshift.

The highly variable time evolution of star-forming cores identified with dendrograms

ABSTRACT We investigate the time evolution of dense cores identified in molecular cloud simulations using dendrograms, which are a common tool to identify hierarchical structure in simulations and observations of star formation. We develop an algorithm to link dendrogram structures through time using the three-dimensional density field from magnetohydrodynamical simulations, thus creating histories for all dense cores in the domain. We find that the population-wide distributions of core properties are relatively invariant in time, and quantities like the core mass function match with observations. Despite this consistency, an individual core may undergo large (>40 per cent), stochastic variations due to the redefinition of the dendrogram structure between time-steps. This variation occurs independent of environment and stellar content. We identify a population of short-lived (<200 kyr) overdensities masquerading as dense cores that may comprise $\sim\!20$ per cent of any time snapshot. Finally, we note the importance of considering the full history of cores when interpreting the origin of the initial mass function; we find that, especially for systems containing multiple stars, the core mass defined by a dendrogram leaf in a snapshot is typically less than the final system stellar mass. This work reinforces that there is no time-stable density contour that defines a star-forming core. The dendrogram itself can induce significant structure variation between time-steps due to small changes in the density field. Thus, one must use caution when comparing dendrograms of regions with different ages or environment properties because differences in dendrogram structure may not come solely from the physical evolution of dense cores.

EDGE: The Origin of Scatter in Ultra-faint Dwarf Stellar Masses and Surface Brightnesses

Abstract We demonstrate how the least luminous galaxies in the universe, ultra-faint dwarf galaxies, are sensitive to their dynamical mass at the time of cosmic reionization. We select a low-mass ( ) dark matter halo from a cosmological volume, and perform zoom hydrodynamical simulations with multiple alternative histories using “genetically modified” initial conditions. Earlier-forming ultra-faints have higher stellar mass today, due to a longer period of star formation before their quenching by reionization. Our histories all converge to the same final dynamical mass, demonstrating the existence of extended scatter (≥1 dex) in stellar masses at fixed halo mass due to the diversity of possible histories. One of our variants builds less than 2% of its final dynamical mass before reionization, rapidly quenching in situ star formation. The bulk of its final stellar mass is later grown by dry mergers, depositing stars in the galaxy’s outskirts and hence expanding its effective radius. This mechanism constitutes a new formation scenario for highly diffuse ( , ), metal-poor ( ), ultra-faint ( ) dwarf galaxies within the reach of next-generation low surface brightness surveys.

Disruption of giant molecular clouds and formation of bound star clusters under the influence of momentum stellar feedback

Abstract Energetic feedback from star clusters plays a pivotal role in shaping the dynamical evolution of giant molecular clouds (GMCs). To study the effects of stellar feedback on the star formation efficiency of the clouds and the dynamical response of embedded star clusters, we perform a suite of isolated GMC simulations with star formation and momentum feedback subgrid models using the moving-mesh hydrodynamics code Arepo. The properties of our simulated GMCs span a wide range of initial mass, radius, and velocity configurations. We find that the ratio of the final stellar mass to the total cloud mass, ϵint, scales strongly with the initial cloud surface density and momentum feedback strength. This correlation is explained by an analytic model that considers force balancing between gravity and momentum feedback. For all simulated GMCs, the stellar density profiles are systematically steeper than that of the gas at the epochs of the peaks of star formation, suggesting a centrally concentrated stellar distribution. We also find that star clusters are always in a sub-virial state with a virial parameter ∼0.6 prior to gas expulsion. Both the sub-virial dynamical state and steeper stellar density profiles prevent clusters from dispersal during the gas removal phase of their evolution. The final cluster bound fraction is a continuously increasing function of ϵint. GMCs with star formation efficiency smaller than 0.5 are still able to form clusters with large bound fractions.

Radiation pressure in super star cluster formation

The physics of star formation at its extreme, in the nuclei of the densest and the most massive star clusters in the universe - potential massive black hole nurseries - has for decades eluded scrutiny. Spectroscopy of these systems has been scarce, whereas theoretical arguments suggest that radiation pressure on dust grains somehow inhibits star formation. Here, we harness an accelerated Monte Carlo radiation transport scheme to report a radiation hydrodynamical simulation of super star cluster formation in turbulent clouds. We find that radiation pressure reduces the global star formation efficiency by 30-35%, and the star formation rate by 15-50%, both relative to a radiation-free control run. Overall, radiation pressure does not terminate the gas supply for star formation and the final stellar mass of the most massive cluster is $\sim 1.3\times10^6\,M_\odot$. The limited impact as compared to in idealized theoretical models is attributed to a radiation-matter anti-correlation in the supersonically turbulent, gravitationally collapsing medium. In isolated regions outside massive clusters, where the gas distribution is less disturbed, radiation pressure is more effective in limiting star formation. The resulting stellar density at the cluster core is $\ge 10^{8}\,M_{\odot}\,\textrm{pc}^{-3}$, with stellar velocity dispersion $\gtrsim 70\,\text{km}\,\text{s}^{-1}$. We conclude that the super star cluster nucleus is propitious to the formation of very massive stars via dynamical core collapse and stellar merging. We speculate that the very massive star may avoid the claimed catastrophic mass loss by continuing to accrete dense gas condensing from a gravitationally-confined ionized phase.

Approximations to galaxy star formation rate histories: properties and uses of two examples

Galaxies evolve via a complex interaction of numerous different physical processes, scales and components. In spite of this, overall trends often appear. Simplified models for galaxy histories can be used to search for and capture such emergent trends, and thus to interpret and compare results of galaxy formation models to each other and to nature. Here, two approximations are applied to galaxy integrated star formation rate histories, drawn from a semi-analytic model grafted onto a dark matter simulation. Both a lognormal functional form and principal component analysis (PCA) approximate the integrated star formation rate histories fairly well. Machine learning, based upon simplified galaxy halo histories, is somewhat successful at recovering both fits. The fits to the histories give fixed time star formation rates which have notable scatter from their true fixed time rates at final time, especially for quiescent and "green valley" galaxies, and more so for the PCA fit. For classifying galaxies into subfamilies sharing similar integrated histories, both approximations are better than using final stellar mass or specific star formation rate. Several subsamples from the simulation illustrate how these simple parameterizations can provide points of contact for comparisons between different galaxy formation samples, or more generally, models. As a side result, the halo masses of simulated galaxies with early peak star formation rate (according to the lognormal fit) lie on one of two tracks. The small fraction of galaxies with a lower halo mass at peak star formation rate appear to stall in their halo growth, even though they are central in their host halos.

First results from the IllustrisTNG simulations: the galaxy colour bimodality

We introduce the first two simulations of the IllustrisTNG project, a next generation of cosmological magnetohydrodynamical simulations, focusing on the optical colors of galaxies. We explore TNG100, a rerun of the original Illustris box, and TNG300, which includes 2x2500^3 resolution elements in a volume twenty times larger. Here we present first results on the galaxy color bimodality at low redshift. Accounting for the attenuation of stellar light by dust, we compare the simulated (g-r) colors of 10^9 < M*/Msun < 10^12.5 galaxies to the observed distribution from the Sloan Digital Sky Survey (SDSS). We find a striking improvement with respect to the original Illustris simulation, as well as excellent quantitative agreement in comparison to the observations, with a sharp transition in median color from blue to red at a characteristic M* ~ 10^10.5 Msun. Investigating the build-up of the color-mass plane and the formation of the red sequence, we demonstrate that the primary driver of galaxy color transition in the TNG model is supermassive blackhole feedback in its low-accretion state. Across the entire population we measure a median color transition timescale dt_green of ~1.6 Gyr, a value which drops for increasingly massive galaxies. We find signatures of the physical process of quenching: at fixed stellar mass, redder galaxies have lower SFRs, gas fractions, and gas metallicities; their stellar populations are also older and their large-scale interstellar magnetic fields weaker than in bluer galaxies. Finally, we measure the amount of stellar mass growth on the red sequence. Galaxies with M* > 10^11 Msun which redden at z<1 accumulate on average ~25% of their final z=0 mass post-reddening; at the same time, ~18% of such massive galaxies acquire half or more of their final stellar mass while on the red sequence.

Some observational tests of a minimal galaxy formation model

Dark matter simulations can serve as a basis for creating galaxy histories via the galaxy-dark matter connection. Here, one such model by Becker (2015) is implemented with several variations on three different dark matter simulations. Stellar mass and star formation rates are assigned to all simulation subhalos at all times, using subhalo mass gain to determine stellar mass gain. The observational properties of the resulting galaxy distributions are compared to each other and observations for a range of redshifts from 0-2. Although many of the galaxy distributions seem reasonable, there are noticeable differences as simulations, subhalo mass gain definitions, or subhalo mass definitions are altered, suggesting that the model should change as these properties are varied. Agreement with observations may improve by including redshift dependence in the added-by-hand random contribution to star formation rate. There appears to be an excess of faint quiescent galaxies as well (perhaps due in part to differing definitions of quiescence). The ensemble of galaxy formation histories for these models tend to have more scatter around their average histories (for a fixed final stellar mass) than the two more predictive and elaborate semi-analytic models of Guo et al (2013) and Henriques et al (2015), and require more basis fluctuations (using PCA) to capture 90 percent of the scatter around their average histories. The codes to plot model predictions (in some cases alongside observational data) are publicly available to test other mock catalogues at https://github.com/jdcphysics/validation . Information on how to use these codes is in the appendix.

PROTOSTELLAR OUTFLOWS AND RADIATIVE FEEDBACK FROM MASSIVE STARS. II. FEEDBACK, STAR-FORMATION EFFICIENCY, AND OUTFLOW BROADENING

ABSTRACT We perform two-dimensional axially symmetric radiation hydrodynamic simulations to assess the impact of outflows and radiative force feedback from massive protostars by varying when the protostellar outflow starts, and to determine the ratio of ejection to accretion rates and the strength of the wide-angle disk wind component. The star-formation efficiency, i.e., the ratio of final stellar mass to initial core mass, is dominated by radiative forces and the ratio of outflow to accretion rates. Increasing this ratio has three effects. First, the protostar grows slower with a lower luminosity at any given time, lowering radiative feedback. Second, bipolar cavities cleared by the outflow become larger, further diminishing radiative feedback on disk and core scales. Third, the higher momentum outflow sweeps up more material from the collapsing envelope, decreasing the protostar's potential mass reservoir via entrainment. The star-formation efficiency varies with the ratio of ejection to accretion rates from 50% in the case of very weak outflows to as low as 20% for very strong outflows. At latitudes between the low-density bipolar cavity and the high-density accretion disk, wide-angle disk winds remove some of the gas, which otherwise would be part of the accretion flow onto the disk; varying the strength of these wide-angle disk winds, however, alters the final star-formation efficiency by only ±6%. For all cases, the opening angle of the bipolar outflow cavity remains below 20° during early protostellar accretion phases, increasing rapidly up to 65° at the onset of radiation pressure feedback.

Variations of the stellar initial mass function in semi-analytical models: implications for the mass assembly and the chemical enrichment of galaxies in the gaea model

In this work, we investigate the implications of the Integrated Galaxy-wide stellar Initial Mass Function (IGIMF) approach in the framework of the semi-analytic model GAEA (GAlaxy Evolution and Assembly), which features a detailed treatment of chemical enrichment and stellar feedback. The IGIMF provides an analytic description of the dependence of the stellar IMF shape on the rate of star formation in galaxies. We find that our model with a universal IMF predicts a rather flat [$\alpha$/Fe]-stellar mass relation. The model assuming the IGIMF, instead, is able to reproduce the observed increase of $\alpha$-enhancement with stellar mass, in agreement with previous studies. This is mainly due to the fact that massive galaxies are characterized by larger star formation rates at high-redshift, leading to stronger $\alpha$-enhancement with respect to low-mass galaxies. At the same time, the IGIMF hypothesis does not affect significantly the trend for shorter star formation timescales for more massive galaxies. We argue that in the IGIMF scenario the [$\alpha$/Fe] ratios are good tracers of the highest star formation events. The final stellar masses and mass-to-light-ratio of our model massive galaxies are larger than those estimated from the synthetic photometry assuming a universal IMF, providing a self-consistent interpretation of similar recent results, based on dynamical analysis of local early type galaxies.

FORMATION OF MASSIVE PRIMORDIAL STARS: INTERMITTENT UV FEEDBACK WITH EPISODIC MASS ACCRETION

ABSTRACT We present coupled stellar evolution (SE) and 3D radiation-hydrodynamic (RHD) simulations of the evolution of primordial protostars, their immediate environment, and the dynamic accretion history under the influence of stellar ionizing and dissociating UV feedback. Our coupled SE RHD calculations result in a wide diversity of final stellar masses covering 10 ≲ M * ≲ 103 . The formation of very massive (≳250 ) stars is possible under weak UV feedback, whereas ordinary massive (a few ×10 ) stars form when UV feedback can efficiently halt the accretion. This may explain the peculiar abundance pattern of a Galactic metal-poor star recently reported by Aoki et al., possibly the observational signature of very massive precursor primordial stars. Weak UV feedback occurs in cases of variable accretion, in particular when repeated short accretion bursts temporarily exceed 0.01 , causing the protostar to inflate. In the bloated state, the protostar has low surface temperature and UV feedback is suppressed until the star eventually contracts, on a thermal adjustment timescale, to create an H ii region. If the delay time between successive accretion bursts is sufficiently short, the protostar remains bloated for extended periods, initiating at most only short periods of UV feedback. Disk fragmentation does not necessarily reduce the final stellar mass. Quite the contrary, we find that disk fragmentation enhances episodic accretion as many fragments migrate inward and are accreted onto the star, thus allowing continued stellar mass growth under conditions of intermittent UV feedback. This trend becomes more prominent as we improve the resolution of our simulations. We argue that simulations with significantly higher resolution than reported previously are needed to derive accurate gas mass accretion rates onto primordial protostars.

NEUTRON STAR MERGERS AS THE ORIGIN OF r -PROCESS ELEMENTS IN THE GALACTIC HALO BASED ON THE SUB-HALO CLUSTERING SCENARIO

Binary mergers (NSMs) of double neutron star (and black hole-neutron star) systems are suggested to be major sites of r-process elements in the Galaxy by recent hydrodynamical and nucleosynthesis studies. It has been pointed out, however, that the estimated long lifetimes of neutron star binaries are in conflict with the presence of r-process-enhanced halo stars at metallicities as low as [Fe/H] ~ -3. To resolve this problem, we examine the role of NSMs in the early Galactic chemical evolution on the assumption that the Galactic halo was formed from merging sub-halos. We present simple models for the chemical evolution of sub-halos with total final stellar masses between 10^4 M_solar and 2 x 10^8 M_solar. Typical lifetimes of compact binaries are assumed to be 100 Myr (for 95% of their population) and 1 Myr (for 5%), according to recent binary population synthesis studies. The resulting metallcities of sub-halos and their ensemble are consistent with the observed mass-metallicity relation of dwarf galaxies in the Local Group, and the metallicity distribution of the Galactic halo, respectively. We find that the r-process abundance ratios [r/Fe] start increasing at [Fe/H] <= -3 if the star formation efficiencies are smaller for less massive sub-halos. In addition, the sub-solar [r/Fe] values (observed as [Ba/Fe] ~ -1.5 for [Fe/H] < -3) are explained by the contribution from the short-lived (~1 Myr) binaries. Our results indicate that NSMs may have a substantial contribution to the r-process element abundances throughout the Galactic history.

Characterizing simulated galaxy stellar mass histories

Cosmological galaxy formation simulations can now produce rich and diverse ensembles of galaxy histories. These simulated galaxy histories, taken all together, provide an answer to the question "how do galaxies form?" for the models used to construct them. We characterize such galaxy history ensembles both to understand their properties and to identify points of comparison for histories within a given galaxy formation model or between different galaxy formation models and simulations. We focus primarily on stellar mass histories of galaxies with the same final stellar mass, for six final stellar mass values and for three different simulated galaxy formation models (a semi-analytic model built upon the dark matter Millennium simulation and two models from the hydrodynamical OverWhelmingly Large Simulations project). Using principal component analysis (PCA) to classify scatter around the average stellar mass history, we find that one fluctuation dominates for all sets of histories we consider, although its shape and contribution can vary between samples. We correlate the PCA characterization with several $z=0$ galaxy properties (to connect with survey observables) and also compare it to some other galaxy history properties. We then explore separating galaxy stellar mass histories into classes, using the largest PCA contribution, k-means clustering, and simple Gaussian mixture models. For three component models, these different methods often gave similar results. These history classification methods provide a succinct and often quick way to characterize changes in the full ensemble of histories of a simulated population as physical assumptions are varied, to compare histories of different simulated populations to each other, and to assess the relation of simulated histories to fixed time observations.

High-redshift quasars host galaxies: is there a stellar mass crisis?

We investigate the evolutionary properties of a sample of quasars at 5<z<6.4 using the semi-analytical hierarchical model GAMETE/QSOdust. We find that the observed properties of these quasars are well reproduced by a common formation scenario in which stars form according to a standard IMF, via quiescent star formation and efficient merger-driven bursts, while the central BH grows via gas accretion and BH-BH mergers. Eventually, a strong AGN driven wind starts to clear up the ISM of dust and gas, damping the star formation and un-obscuring the line of sight toward the QSO. In this scenario, all the QSOs hosts have final stellar masses in the range $(4-6)\times 10^{11} M_{sun}$, a factor 3-30 larger than the upper limits allowed by the observations. We discuss alternative scenarios to alleviate this apparent tension: the most likely explanation resides in the large uncertainties that still affect dynamical mass measurements in these high-z galaxies. In addition, during the transition between the starburst-dominated and the active QSO phase, we predict that about 40% of the progenitor galaxies can be classified as Sub Millimeter Galaxies, although their number rapidly decreases with redshift.

PURELY DRY MERGERS DO NOT EXPLAIN THE OBSERVED EVOLUTION OF MASSIVE EARLY-TYPE GALAXIES SINCEz∼ 1

Several studies have suggested that the observed size evolution of massive early-type galaxies (ETGs) can be explained as a combination of dry mergers and progenitor bias, at least since z\sim1. In this paper we carry out a new test of the dry-merger scenario based on recent lensing measurements of the evolution of the mass density profile of ETGs. We construct a theoretical model for the joint evolution of the size and mass density profile slope \gamma' driven by dry mergers occurring at rates given by cosmological simulations. Such dry-merger model predicts a strong decrease of \gamma' with cosmic time, inconsistent with the almost constant \gamma' inferred from observations in the redshift range 0<z<1. We then show with a simple toy model that a modest amount of cold gas in the mergers -- consistent with the upper limits on recent star formation in ETGs -- is sufficient to reconcile the model with measurements of \gamma'. By fitting for the amount of gas accreted during mergers, we find that models with dissipation are consistent with observations of the evolution in both size and density slope, if \sim4% of the total final stellar mass arises from the gas accreted since z\sim1. Purely dry merger models are ruled out at > 99% CL. We thus suggest a scenario where the outer regions of massive ETGs grow by accretion of stars and dark matter, while small amounts of dissipation and nuclear star formation conspire to keep the mass density profile constant and approximately isothermal.

ON THE NATURE OF THE DEEPLY EMBEDDED PROTOSTAR OMC-2 FIR 4

We use mid-infrared to submillimeter data from the Spitzer, Herschel, and APEX telescopes to study the bright sub-mm source OMC-2 FIR 4. We find a point source at 8, 24, and 70 $\mu$m, and a compact, but extended source at 160, 350, and 870 $\mu$m. The peak of the emission from 8 to 70 $\mu$m, attributed to the protostar associated with FIR 4, is displaced relative to the peak of the extended emission; the latter represents the large molecular core the protostar is embedded within. We determine that the protostar has a bolometric luminosity of 37 Lsun, although including more extended emission surrounding the point source raises this value to 86 Lsun. Radiative transfer models of the protostellar system fit the observed SED well and yield a total luminosity of most likely less than 100 Lsun. Our models suggest that the bolometric luminosity of the protostar could be just 12-14 Lsun, while the luminosity of the colder (~ 20 K) extended core could be around 100 Lsun, with a mass of about 27 Msun. Our derived luminosities for the protostar OMC-2 FIR 4 are in direct contradiction with previous claims of a total luminosity of 1000 Lsun (Crimier et al 2009). Furthermore, we find evidence from far-infrared molecular spectra (Kama et al. 2013, Manoj et al. 2013) and 3.6 cm emission (Reipurth et al 1999) that FIR 4 drives an outflow. The final stellar mass the protostar will ultimately achieve is uncertain due to its association with the large reservoir of mass found in the cold core.