Source author record

Eric Keto

Eric Keto appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

Clustered Formation of Massive Stars within an Ionized Rotating Disk

We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich envelope. The MCMC fits to the data give total stellar masses $M_\star$ from 120 to $200~M_\odot$, with much smaller ionized-gas masses $M_\mathrm{ion-gas} = 0.2$ to $0.25~M_\odot$. The stellar mass is distributed within the gravitational radius $R_g\approx 1000$ to 1500 au, where the ionized gas is bound. The viewing inclination angle from the face-on orientation is $i = 49$ to $56~°$. Radial motions at radii $r > R_g$ converge to $v_{r,0} \approx 8.7$ km/s, or about the speed of sound of ionized gas, indicating that this gas is marginally unbound at most. From additional constraints on the ionizing-photon rate and far-IR luminosity of the region, we conclude that the stellar cluster consists of a few massive stars with $M_\mathrm{star} = 32$ to $60~M_\odot$, or one star in this range of masses accompanied by a population of lower-mass stars. Any active accretion of ionized gas onto the massive (proto)stars is residual. The inferred cluster density is very large, comparable to that reported at similar scales in the Galactic Center. Stellar interactions are likely to occur within the next Myr.

preprint2022arXiv

The Central 1000 au of a Pre-stellar Core Revealed with ALMA. II. Almost Complete Freeze-out

Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particular CO, are heavily depleted in pre-stellar cores due to the efficient molecular freeze-out onto the surface of cold dust grains. However, N-bearing species such as NH$_3$ and, in particular, its deuterated isotopologues, appear to maintain high abundances where CO molecules are mainly in solid phase. Thanks to ALMA, we present here the first clear observational evidence of NH$_2$D freeze-out toward the L1544 pre-stellar core, suggestive of the presence of a"complete-depletion zone" within a $\simeq$1800 au radius, in agreement with astrochemical pre-stellar core model predictions. Our state-of-the-art chemical model coupled with a non-LTE radiative transfer code demonstrates that NH$_2$D becomes mainly incorporated in icy mantles in the central 2000 au and starts freezing-out already at $\simeq$7000 au. Radiative transfer effects within the pre-stellar core cause the NH$_2$D(1$_{11}$-1$_{01}$) emission to appear centrally concentrated, with a flattened distribution within the central $\simeq$3000 au, unlike the 1.3 mm dust continuum emission which shows a clear peak within the central $\simeq$1800 au. This prevented NH$_2$D freeze-out to be detected in previous observations, where the central 1000 au cannot be spatially resolved.

preprint2021arXiv

An observational correlation between magnetic field, angular momentum and fragmentation in the envelopes of Class 0 protostars?

To assess the potential role of magnetic fields in regulating the envelope rotation and the fragmentation of Class 0 protostars, we carried out observations of the dust polarized emission at 0.87 mm with the SMA, in the envelopes of a large sample of 20 Class 0 protostars. We estimate the mean magnetic field orientation over the central 1000 au envelope scales and compared it to that of the protostellar outflow in order to study the relation between their misalignment and the kinematics of the circumstellar gas. We discover a strong relationship between the misalignment of the magnetic field orientation with the outflow and the amount of angular momentum observed at similar scales in the protostellar envelope, revealing a potential link between the kinetic and the magnetic energy at envelope scales. The relation could be driven by favored B misalignments in more dynamical envelopes or a dependence of the envelope dynamics with the large-scale B initial configuration. Comparing the trend with the presence of fragmentation, we observe that single sources are mostly associated with conditions of low angular momentum in the inner envelope and good alignment of the magnetic field with protostellar outflows, at intermediate scales. Our results suggest that the properties of the magnetic field in protostellar envelopes bear a tight relationship with the rotating-infalling gas directly involved in the star and disk formation: we find that it may not only influence the fragmentation of protostellar cores into multiple stellar systems, but also set the conditions establishing the pristine properties of planet-forming disks.

preprint2020arXiv

A Turbulent-Entropic Instability and the Fragmentation of Star-Forming Clouds

The kinetic energy of supersonic turbulence within interstellar clouds is subject to cooling by dissipation in shocks and subsequent line radiation. The clouds are therefore susceptible to a condensation process controlled by the specific entropy. In a form analogous to the thermodynamic entropy, the entropy for supersonic turbulence is proportional to the log of the product of the mean turbulent velocity and the size scale. We derive a dispersion relation for the growth of entropic instabilities in a spherical self-gravitating cloud and find that there is a critical maximum dissipation time scale, about equal to the crossing time, that allows for fragmentation and subsequent star formation. However, the time scale for the loss of turbulent energy may be shorter or longer, for example with rapid thermal cooling or the injection of mechanical energy. Differences in the time scale for energy loss in different star-forming regions may result in differences in the outcome, for example, in the initial mass function.

preprint2020arXiv

CMZoom: Survey Overview and First Data Release

We present an overview of the CMZoom survey and its first data release. CMZoom is the first blind, high-resolution survey of the Central Molecular Zone (CMZ; the inner 500 pc of the Milky Way) at wavelengths sensitive to the pre-cursors of high-mass stars. CMZoom is a 500-hour Large Program on the Submillimeter Array (SMA) that mapped at 1.3 mm all of the gas and dust in the CMZ above a molecular hydrogen column density of 10^23 cm^-2 at a resolution of ~3" (0.1 pc). In this paper, we focus on the 1.3 mm dust continuum and its data release, but also describe CMZoom spectral line data which will be released in a forthcoming publication. While CMZoom detected many regions with rich and complex substructure, its key result is an overall deficit in compact substructures on 0.1 - 2 pc scales (the compact dense gas fraction: CDGF). In comparison with clouds in the Galactic disk, the CDGF in the CMZ is substantially lower, despite having much higher average column densities. CMZ clouds with high CDGFs are well-known sites of active star formation. The inability of most gas in the CMZ to form compact substructures is likely responsible for the dearth of star formation in the CMZ, surprising considering its high density. The factors responsible for the low CDGF are not yet understood but are plausibly due to the extreme environment of the CMZ, having far-reaching ramifications for our understanding of the star formation process across the cosmos.

preprint2020arXiv

Stability and Solution of the Time-Dependent Bondi-Parker Flow

Bondi (1952) and Parker (1958} derived a steady-state solution for Bernouilli's equation in spherical symmetry around a point mass for two cases, respectively, an inward accretion flow and an outward wind. Left unanswered were the stability of the steady-state solution, the solution itself of time-dependent flows, whether the time-dependent flows would evolve to the steady-state, and under what conditions a transonic flow would develop. In a Hamiltonian description, we find that the steady state solution is equivalent to the Lagrangian implying that time-dependent flows evolve to the steady state. We find that the second variation is definite in sign for isothermal and adiabatic flows, implying at least linear stability. We solve the partial differential equation for the time-dependent flow as an initial-value problem and find that a transonic flow develops under a wide range of realistic initial conditions. We present some examples of time-dependent solutions.

preprint2016arXiv

A Hot and Massive Accretion Disk around the High-Mass Protostar IRAS 20126+4104

We present new spectral line observations of the CH3CN molecule in the accretion disk around the massive protostar IRAS 20126+4104 with the Submillimeter Array that for the first time measure the disk density, temperature, and rotational velocity with sufficient resolution (0.37", equivalent to ~600 AU) to assess the gravitational stability of the disk through the Toomre-Q parameter. Our observations resolve the central 2000 AU region that shows steeper velocity gradients with increasing upper state energy, indicating an increase in the rotational velocity of the hotter gas nearer the star. Such spin-up motions are characteristics of an accretion flow in a rotationally supported disk. We compare the observed data with synthetic image cubes produced by three-dimensional radiative transfer models describing a thin flared disk in Keplerian motion enveloped within the centrifugal radius of an angular-momentum-conserving accretion flow. Given a luminosity of 1.3x10^4 Lsun, the optimized model gives a disk mass of 1.5 Msun and a radius of 858 AU rotating about a 12.0 Msun protostar with a disk mass accretion rate of 3.9x10^{-5} Msun/yr. Our study finds that, in contrast to some theoretical expectations, the disk is hot and stable to fragmentation with Q > 2.8 at all radii which permits a smooth accretion flow. These results put forward the first constraints on gravitational instabilities in massive protostellar disks, which are closely connected to the formation of companion stars and planetary systems by fragmentation.

preprint2015arXiv

A Keplerian-like disk around the forming O-type star AFGL 4176

We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3" resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm (source mm1) has a deconvolved size of 870+/-110 AU x 330+/-300 AU and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperature of 190 K. The first-moment maps, pixel-to-pixel line modeling, assuming local thermodynamic equilibrium (LTE), and position-velocity diagrams of the CH3CN J=13-12 K-line emission all show a velocity gradient along the major axis of the source, coupled with an increase in velocity at small radii, consistent with Keplerian-like rotation. The LTE line modeling shows that where CH3CN J=13-12 is excited, the temperatures in the disk range from ~70 to at least 300 K and that the H2 column density peaks at 2.8x10^24 cm^-2. In addition, we present Atacama Pathfinder Experiment (APEX) 12CO observations which show a large-scale outflow from AFGL 4176 perpendicular to the major axis of mm1, supporting the disk interpretation. Finally, we present a radiative transfer model of a Keplerian disk surrounding an O7 star, with a disk mass and radius of 12 M_sun and 2000 AU, that reproduces the line and continuum data, further supporting our conclusion that our observations have uncovered a Keplerian disk around an O-type star.

preprint2014arXiv

Chemistry and Radiative Transfer of Water in Cold, Dense Clouds

The Herschel Space Observatory's recent detections of water vapor in the cold, dense cloud L1544 allow a direct comparison between observations and chemical models for oxygen species in conditions just before star formation. We explain a chemical model for gas phase water, simplified for the limited number of reactions or processes that are active in extreme cold ($<$ 15 K). In this model, water is removed from the gas phase by freezing onto grains and by photodissociation. Water is formed as ice on the surface of dust grains from O and OH and released into the gas phase by photodesorption. The reactions are fast enough with respect to the slow dynamical evolution of L1544 that the gas phase water is in equilibrium for the local conditions thoughout the cloud. We explain the paradoxical radiative transfer of the H$_2$O ($1_{10}-1_{01}$) line. Despite discouragingly high optical depth caused by the large Einstein A coefficient, the subcritical excitation in the cold, rarefied H$_2$ causes the line brightness to scale linearly with column density. Thus the water line can provide information on the chemical and dynamical processes in the darkest region in the center of a cold, dense cloud. The inverse P-Cygni profile of the observed water line generally indicates a contracting cloud. This profile is reproduced with a dynamical model of slow contraction from unstable quasi-static hydrodynamic equilibrium (an unstable Bonnor-Ebert sphere).

preprint2014arXiv

Core and filament formation in magnetized, self-gravitating isothermal layers

We examine the role of the gravitational instability in an isothermal, self-gravitating layer threaded by magnetic fields on the formation of filaments and dense cores. Using numerical simulation we follow the non-linear evolution of a perturbed equilibrium layer. The linear evolution of such a layer is described in the analytic work of Nagai et al (1998). We find that filaments and dense cores form simultaneously. Depending on the initial magnetic field, the resulting filaments form either a spiderweb-like network (for weak magnetic fields) or a network of parallel filaments aligned perpendicular to the magnetic field lines (for strong magnetic fields). Although the filaments are radially collapsing, the density profile of their central region (up to the thermal scale height) can be approximated by a hydrodynamical equilibrium density structure. Thus, the magnetic field does not play a significant role in setting the density distribution of the filaments. The density distribution outside of the central region deviates from the equilibrium. The radial column density distribution is then flatter than the expected power law of $r^{-4}$ and similar to filament profiles observed with Herschel. Our results does not explain the near constant filament width of $\sim 0.1$pc. However, our model does not include turbulent motions. It is expected that accretion-driven amplification of these turbulent motions provides additional support within the filaments against gravitational collapse. Finally, we interpret the filamentary network of the massive star forming complex G14.225-0.506 in terms of the gravitational instability model and find that the properties of the complex are consistent with being formed out of an unstable layer threaded by a strong, parallel magnetic field.

preprint2014arXiv

From Filaments to Oscillating Starless Cores

Long wavelength sonic oscillations are observed or inferred in many of the small, dark molecular clouds, the starless cores, that are the precursors to protostars. The oscillations provide significant internal energy and the time scale for their dissipation may control the rate of star formation in the starless cores. Despite their importance, their origin is unknown. We explore one hypothesis that the oscillations develop as a starless core forms from a filament. We model this process with a numerical hydrodynamic simulation and generate synthetic molecular line observations with a radiative transfer simulation. Comparison with actual observations shows general agreement suggesting the proposed mechanism is viable.

preprint2014arXiv

Magnetic Fields and Massive Star Formation

Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $μ$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in these objects ranges from ordered hour-glass configurations to more chaotic distributions. By comparing the SMA data with the single dish data at parsec scales, we found that magnetic fields at dense core scales are either aligned within $40^\circ$ of or perpendicular to the parsec-scale magnetic fields. This finding indicates that magnetic fields play an important role during the collapse and fragmentation of massive molecular clumps and the formation of dense cores. We further compare magnetic fields in dense cores with the major axis of molecular outflows. Despite a limited number of outflows, we found that the outflow axis appears to be randomly oriented with respect to the magnetic field in the core. This result suggests that at the scale of accretion disks ($\lsim 10^3$ AU), angular momentum and dynamic interactions possibly due to close binary or multiple systems dominate over magnetic fields. With this unprecedentedly large sample massive clumps, we argue on a statistical basis that magnetic fields play an important role during the formation of dense cores at spatial scale of 0.01 - 0.1 pc in the context of massive star and cluster star formation.

preprint2014arXiv

The dynamics of collapsing cores and star formation

Low-mass stars are generally understood to form by the gravitational collapse of the dense molecular clouds known as starless cores. Continuum observations have not been able to distinguish among the several different hypotheses that describe the collapse because the predicted density distributions are the almost the same, as they are for all spherical self-gravitating clouds. However, the predicted contraction velocities are different enough that the models can be discriminated by comparing the velocities at large and small radii. This can be done by observing at least two different molecular line transitions that are excited at different densities. For example, the spectral lines of the H2O (110 - 101) and C18O (1-0) have critical densities for collisional de-excitation that differ by 5 orders of magnitude. We compare observations of these lines from the contracting starless core L1544 against the spectra predicted for several different hypothetical models of contraction including the Larson-Penston flow, the inside-out collapse of the singular isothermal sphere, the quasi-equilibrium contraction of an unstable Bonnor-Ebert sphere, and the non-equilibrium collapse of an over-dense Bonnor-Ebert sphere. Only the model of the unstable quasi-equilibrium Bonnor-Ebert sphere is able to produce the observed shapes of both spectral lines. This model allows us to interpret other observations of molecular lines in L1544 to find that the inward velocities seen in observations of CS(2-1) and N2H+ are located within the starless core itself, in particular in the region where the density profile follows an inverse square law. If these conclusions were to hold in the analysis of other starless cores, this would imply that the formation of hydrostatic clouds within the turbulent interstellar medium is not only possible but not exceptional and may be an evolutionary phase in low-mass star formation.

preprint2013arXiv

The Dynamics of Ultracompact HII Regions

Many ultracompact HII regions exhibit a cometary morphology in radio continuum emission. In such regions, a young massive star is probably ablating, through its ultraviolet radiation, the molecular cloud clump that spawned it. On one side of the star, the radiation drives an ionization front that stalls in dense molecular gas. On the other side, ionized gas streams outward into the more rarefied environment. This wind is underpressured with respect to the neutral gas. The difference in pressure draws in more cloud material, feeding the wind until the densest molecular gas is dissipated. Recent, time-dependent simulations of massive stars turning on within molecular gas show the system evolving in a direction similar to that just described. Here, we explore a semi-analytic model in which the wind is axisymmetric and has already achieved a steady state. Adoption of this simplified picture allows us to study the dependence of both the wind and its bounding ionization front on the stellar luminosity, the peak molecular density, and the displacement of the star from the center of the clump. For typical parameter values, the wind accelerates transonically to a speed of about 15 km / s, and transports mass outward at a rate of 10^-4 solar masses per year. Stellar radiation pressure acts to steepen the density gradient of the wind.

preprint2012arXiv

First detection of water vapor in a pre-stellar core

Water is a crucial molecule in molecular astrophysics as it controls much of the gas/grain chemistry, including the formation and evolution of more complex organic molecules in ices. Pre-stellar cores provide the original reservoir of material from which future planetary systems are built, but few observational constraints exist on the formation of water and its partitioning between gas and ice in the densest cores. Thanks to the high sensitivity of the Herschel Space Observatory, we report on the first detection of water vapor at high spectral resolution toward a dense cloud on the verge of star formation, the pre-stellar core L1544. The line shows an inverse P-Cygni profile, characteristic of gravitational contraction. To reproduce the observations, water vapor has to be present in the cold and dense central few thousand AU of L1544, where species heavier than Helium are expected to freeze-out onto dust grains, and the ortho:para H2 ratio has to be around 1:1 or larger. The observed amount of water vapor within the core (about 1.5x10^{-6} Msun) can be maintained by Far-UV photons locally produced by the impact of galactic cosmic rays with H2 molecules. Such FUV photons irradiate the icy mantles, liberating water wapor in the core center. Our Herschel data, combined with radiative transfer and chemical/dynamical models, shed light on the interplay between gas and solids in dense interstellar clouds and provide the first measurement of the water vapor abundance profile across the parent cloud of a future solar-type star and its potential planetary system.

preprint2012arXiv

Hierarchical configurations for cross-correlation interferometers with many elements

Array configurations built on a hierarchy of simple elements have excellent properties for cross-correlation imaging interferometers including a smooth distribution of measured Fourier components, high angular resolution, low side lobes, and compact array size. Compared to arrays with a Gaussian distribution of antenna separations, hierarchical arrays (H-arrays) produce beams with higher angular resolution and a tighter concentration of the total power (encircled energy) within a smaller area around the main beam. An attractive feature of H-arrays is their simplicity. The relationships between the Fourier coverage and the array configuration are easy enough to understand that they can be adjusted to achieve different design goals without the need for numerical optimization. H-arrays will be useful for future multi-element interferometers.

preprint2012arXiv

Two Populations of Molecular Clouds in the Antennae Galaxies

Super star clusters --- extremely massive clusters found predominately in starburst environments --- are essential building blocks in the formation of galaxies and thought to dominate star formation in the high-redshift universe. However, the transformation from molecular gas into these ultra-compact star clusters is not well understood. To study this process, we used the Submillimeter Array and the Plateau de Bure Interferometer to obtain high angular resolution (~1.5" or 160 pc) images of the Antennae overlap region in CO(2--1) to search for the molecular progenitors of the super star clusters. We resolve the molecular gas distribution into a large number of clouds, extending the differential cloud mass function down to a 5σcompleteness limit of 3.8x10^5 M_sun. We identify a distinct break in the mass function around log M_mol/M_sun ~ 6.5, which separates the molecular clouds into two distinct populations. The smaller, less massive clouds reside in more quiescent areas in the region, while the larger, more massive clouds cluster around regions of intense star formation. A broken power-law fit to the mass function yields slopes of α= -1.39+/-0.10 and α= -1.44+/-0.14 for the low- and high-mass cloud population, well-matched to the mass function found for super star clusters in the Antennae galaxies. We find large velocity gradients and velocity dispersions at the locations of intense star formation, suggestive of compressive shocks. It is likely that these environmental factors contribute to the formation of the observed massive molecular clouds and super star clusters in the Antennae galaxies.

preprint2011arXiv

Does External Pressure Explain Recent Results for Molecular Clouds?

The recent paper by Heyer et al (2009) indicates that observations of size, linewidth and column density of interstellar clouds do not agree with simple virial equilibrium (VE) as a balance between gravitational and kinetic energies in the sense that the clouds either have too much kinetic energy or too little mass to be bound. This may be explained by violation of VE as suggested by Dobbs et al 2011, by observational underestimation of the masses as suggested by Heyer et al 2009, or by an external pressure acting as an additional confining force as suggested earlier by Heyer et al 2004. The data of Heyer et al. 2009 cannot be explained with a single value for the external pressure, but if different clouds in the sample have different external pressures in the range of Pe/k = E4 to E7 cm-3 K, then most of the clouds could be in pressure virial equilibrium (PVE). In this paper we discuss two consequences of the external pressure. First, we show that the observational data are consistent with the hypothesis (Chieźe 1987) that most clouds are at a critical mass for dynamical stability determined solely by the pressure. Above this mass a cloud is unstable to gravitational collapse or fragmentation. Second, we show that the external pressure modifies the well-known size-linewidth relationship first proposed by Larson (1981) so that the proportionality is no longer constant but depends on the external pressure.

preprint2011arXiv

The standard model of low-mass star formation applied to massive stars: multi-wavelength modelling of IRAS 20126+4104

In order to investigate whether massive stars form similarly to their low-mass counterparts, we have used the standard envelope plus disc geometry successfully applied to low-mass protostars to model the near-IR to sub-millimetre SED and several mid-IR images of the embedded massive star IRAS20126+4104. We have used a Monte Carlo radiative transfer dust code to model the continuum absorption, emission and scattering through two azimuthally symmetric dust geometries, the first consisting of a rotationally flattened envelope with outflow cavities, and the second which also includes a flared accretion disc. Our results show that the envelope plus disc model reproduces the observed SED and images more accurately than the model without a disc, although the latter model more closely reproduces the morphology of the mid-IR emission within a radius of 1.1" or ~1800au. We have put forward several possible causes of this discontinuity, including inner truncation of the disc due to stellar irradiation, or precession of the outflow cavity. Our best fitting envelope plus disc model has a disc radius of 9200 au. We find that it is unlikely that the outer regions of such a disc would be in hydrostatic or centrifugal equilibrium, however we calculate that the temperatures within the disc would keep it stable to fragmentation.

preprint2011arXiv

Time Variability in Simulated Ultracompact and Hypercompact HII Regions

Ultracompact and hypercompact HII regions appear when a star with a mass larger than about 15 solar masses starts to ionize its own environment. Recent observations of time variability in these objects are one of the pieces of evidence that suggest that at least some of them harbor stars that are still accreting from an infalling neutral accretion flow that becomes ionized in its innermost part. We present an analysis of the properties of the HII regions formed in the 3D radiation-hydrodynamic simulations presented by Peters et al. as a function of time. Flickering of the HII regions is a natural outcome of this model. The radio-continuum fluxes of the simulated HII regions, as well as their flux and size variations are in agreement with the available observations. From the simulations, we estimate that a small but non-negligible fraction (~ 10 %) of observed HII regions should have detectable flux variations (larger than 10 %) on timescales of ~ 10 years, with positive variations being more likely to happen than negative variations. A novel result of these simulations is that negative flux changes do happen, in contrast to the simple expectation of ever growing HII regions. We also explore the temporal correlations between properties that are directly observed (flux and size) and other quantities like density and ionization rates.

preprint2010arXiv

An Interferometric Spectral-line Survey of IRC+10216 in the 345 GHz Band

We report a spectral-line survey of the extreme carbon star IRC+10216 carried out between 293.9 and 354.8 GHz with the Submillimeter Array. A total of 442 lines were detected, more than 200 for the first time; 149 are unassigned. Maps at an angular resolution of ~3" were obtained for each line. A substantial new population of narrow lines with an expansion velocity of ~4 km/s (i.e. ~ 30% of the terminal velocity) was detected. Most of these are attributed to rotational transitions within vibrationally excited states, emitted from energy levels above the v=0, J=0 ground state with excitation energy of 1000-3000 K. Emission from these lines appears to be centered on the star with an angular extent of <1". We use multiple transitions detected in several molecules to derive physical conditions in this inner envelope of IRC+10216.

preprint2010arXiv

Circumstellar Structure around Evolved Stars in the Cygnus-X Star Formation Region

We present observations of newly discovered 24 micron circumstellar structures detected with the Multiband Imaging Photometer for Spitzer (MIPS) around three evolved stars in the Cygnus-X star forming region. One of the objects, BD+43 3710, has a bipolar nebula, possibly due to an outflow or a torus of material. A second, HBHA 4202-22, a Wolf-Rayet candidate, shows a circular shell of 24 micron emission suggestive of either a limb-brightened shell or disk seen face-on. No diffuse emission was detected around either of these two objects in the Spitzer 3.6-8 micron Infrared Array Camera (IRAC) bands. The third object is the luminous blue variable candidate G79.29+0.46. We resolved the previously known inner ring in all four IRAC bands. The 24 micron emission from the inner ring extends ~1.2 arcmin beyond the shorter wavelength emission, well beyond what can be attributed to the difference in resolutions between MIPS and IRAC. Additionally, we have discovered an outer ring of 24 micron emission, possibly due to an earlier episode of mass loss. For the two shell stars, we present the results of radiative transfer models, constraining the stellar and dust shell parameters. The shells are composed of amorphous carbon grains, plus polycyclic aromatic hydrocarbons in the case of G79.29+0.46. Both G79.29+0.46 and HBHA 4202-22 lie behind the main Cygnus-X cloud. Although G79.29+0.46 may simply be on the far side of the cloud, HBHA 4202-22 is unrelated to the Cygnus-X star formation region.

preprint2010arXiv

Dynamics of Molecular Clouds

We further develop the model of molecular cloud fragmentation introduced in Field, Blackman and Keto (2007; FBK). We show that external pressure acting on fragments establishes a scale-dependent critical mass. Fragments with masses less than the critical value are confined largely by pressure, while those with masses greater than or equal to the critical value collapse under self gravitation. Both types of fragments are commonly observed. Without specifying the source of the external pressure, and without assuming any other scaling relations, we predict the power - law index in the relation between the rms velocity of supersonic motions and the size of fragments . We then investigate the possibility that the external pressure is due to the kinetic energy of H atoms released by photodissociation of hydrogen molecules in the fragment. This can account approximately for the observed values of external pressure and two additional observations: the value of the scaling coefficient in the power law mentioned above, and the observation of outflowing atomic hydrogen around molecular clouds. A further prediction is HI at fragment edges with column densities of order 1E20 per sq. cm and velocities of a few km/s that should be detectable with high resolution 21 cm observations. Finally, we predict the magnitude of the coefficient of dissipation in the observed supersonic flows.

preprint2010arXiv

From the Convergence of Filaments to Disk-Outflow Accretion: Massive-Star Formation in W33A

Interferometric observations of the W33A massive star-formation region, performed with the Submillimeter Array (SMA) and the Very Large Array (VLA) at resolutions from 5 arcsec (0.1 pc) to 0.5 arcsec (0.01 pc) are presented. Our three main findings are: (1) parsec-scale, filamentary structures of cold molecular gas are detected. Two filaments at different velocities intersect in the zone where the star formation is occurring. This is consistent with triggering of the star-formation activity by the convergence of such filaments, as predicted by numerical simulations of star formation initiated by converging flows. (2) The two dusty cores (MM1 and MM2) at the intersection of the filaments are found to be at different evolutionary stages, and each of them is resolved into multiple condensations. MM1 and MM2 have markedly different temperatures, continuum spectral indices, molecular-line spectra, and masses of both stars and gas. (3) The dynamics of the "hot-core" MM1 indicates the presence of a rotating disk in its center (MM1-Main) around a faint free-free source. The stellar mass is estimated to be approximately 10 Msun. A massive molecular outflow is observed along the rotation axis of the disk.

preprint2010arXiv

H II regions: Witnesses to massive star formation

We describe the first three-dimensional simulation of the gravitational collapse of a massive, rotating molecular cloud that includes heating by both non-ionizing and ionizing radiation. We find that as the first protostars gain sufficient mass to ionize the accretion flow, their H II regions are initially gravitationally trapped, but soon begin to rapidly fluctuate between trapped and extended states, in agreement with observations. Over time, the same ultracompact H II region can expand anisotropically, contract again, and take on any of the observed morphological classes. In their extended phases, expanding H II regions drive bipolar neutral outflows characteristic of high-mass star formation. The total lifetime of H II regions is given by the global accretion timescale, rather than their short internal sound-crossing time. The pressure of the hot, ionized gas does not terminate accretion. Instead the final stellar mass is set by fragmentation-induced starvation. Local gravitational instabilities in the accretion flow lead to the build-up of a small cluster of stars, all with relatively high masses due to heating from accretion radiation. These companions subsequently compete with the initial high-mass star for the same common gas reservoir and limit its mass growth. Our findings show that the most significant differences between the formation of low-mass and high-mass stars are all explained as the result of rapid accretion within a dense, gravitationally unstable, ionized flow.

preprint2010arXiv

Modeling molecular hyperfine line emission

In this paper we discuss two approximate methods previously suggested for modeling hyperfine spectral line emission for molecules whose collisional transitions rates between hyperfine levels are unknown. Hyperfine structure is seen in the rotational spectra of many commonly observed molecules such as HCN, HNC, NH3, N2H+, and C17O. The intensities of these spectral lines can be modeled by numerical techniques such as Lambda-iteration that alternately solve the equations of statistical equilibrium and the equation of radiative transfer. However, these calculations require knowledge of both the radiative and collisional rates for all transitions. For most commonly observed radio frequency spectral lines, only the net collisional rates between rotational levels are known. For such cases, two approximate methods have been suggested. The first method, hyperfine statistical equilibrium (HSE), distributes the hyperfine level populations according to their statistical weight, but allows the population of the rotational states to depart from local thermodynamic equilibrium (LTE). The second method, the proportional method approximates the collision rates between the hyperfine levels as fractions of the net rotational rate apportioned according to the statistical degeneracy of the final hyperfine levels. The second method is able to model non-LTE hyperfine emission. We compare simulations of N2H+ hyperfine lines made with approximate and more exact rates and find that satisfactory results are obtained.

preprint2010arXiv

The Decrease of Specific Angular Momentum and the Hot Toroid Formation: The Massive Clump G10.6-0.4

This is the first paper of our series of high resolution (1") studies of the massive star forming region G10.6--0.4. We present the emission line observations of the hot core type tracers (O$^{13}$CS, OCS, SO$_{2}$) with $\sim$0$"$.5 resolution. By comparing the results to the high--resolution NH$_{3}$ absorption line observation, we confirm for the first time the rotationally flattened hot toroid in the central $<$0.1 pc region, which has a rotational axis perpendicular to its geometrical major axis. In addition, we present the observations of NH$_{3}$, $^{13}$CS, and CH$_{3}$CN with $\sim$1$"$ resolution, and follow the dynamics of the molecular accretion flow from the 0.3 pc radius to the inner 0.03 pc radius. With reference to the rotational axis of the hot toroid, we measure the rotational velocity from the molecular emission in the region. The results are consistent with an envelope with a rapid decrease of the specific angular momentum from the outer to the inner region. These new results improve the current understanding of the molecular accretion flow in an ultracompact (UC) H\textsc{ii} region created by the embedded O-type cluster.

preprint2010arXiv

The Evolution of Cloud Cores and the Formation of Stars

For a number of starless cores, self-absorbed molecular line and column density observations have implied the presence of large-amplitude oscillations. We examine the consequences of these oscillations on the evolution of the cores and the interpretation of their observations. We find that the pulsation energy helps support the cores and that the dissipation of this energy can lead toward instability and star formation. In this picture, the core lifetimes are limited by the pulsation decay timescales, dominated by non-linear mode-mode coupling, and on the order of ~few x 10^5--10^6 yr. Notably, this is similar to what is required to explain the relatively low rate of conversion of cores into stars. For cores with large-amplitude oscillations, dust continuum observations may appear asymmetric or irregular. As a consequence, some of the cores that would be classified as supercritical may be dynamically stable when oscillations are taken into account. Thus, our investigation motivates a simple hydrodynamic picture, capable of reproducing many of the features of the progenitors of stars without the inclusion of additional physical processes, such as large-scale magnetic fields.

preprint2010arXiv

The standard model of star formation applied to massive stars: accretion disks and envelopes in molecular lines

We address the question of whether the formation of high-mass stars is similar to or differs from that of solar-mass stars through new molecular line observations and modeling of the accretion flow around the massive protostar IRAS20126+4104. We combine new observations of NH3(1,1) and (2,2) made at the Very Large Array, new observations of CHCN(13-12) made at the Submillimeter Array, previous VLA observations of NH(3,3), NH(4,4), and previous Plateau de Bure observations of C34S(2-1), C34S(5-4), and CHCN(12-11) to obtain a data set of molecular lines covering 15 to 419 K in excitation energy. We compare these observations against simulated molecular line spectra predicted from a model for high-mass star formation based on a scaled-up version of the standard disk-envelope paradigm developed for accretion flows around low-mass stars. We find that in accord with the standard paradigm, the observations require both a warm, dense, rapidly-rotating disk and a cold, diffuse infalling envelope. This study suggests that accretion processes around 10 M stars are similar to those of solar mass stars.

preprint2009arXiv

Dynamics and Depletion in Thermally Supercritical Starless Cores

In previous studies we identified two classes of starless cores, thermally subcritical and supercritical, distinguished by different dynamical behavior and internal structure. Here we study the evolution of the dynamically-unstable, thermally-supercritical cores by means of a numerical hydrodynamic simulation that includes radiative equilibrium and simple molecular chemistry. We use our non-LTE radiative transfer code MOLLIE to predict observable CO and N2H+ line spectra, including the non-LTE hyperfine ratios of N2H+, during the contraction. These are compared against observations of the starless core L1544.

preprint2009arXiv

Formation of an O-Star Cluster by Hierarchical Accretion in G20.08-0.14 N

Spectral line and continuum observations of the ionized and molecular gas in G20.08-0.14 N explore the dynamics of accretion over a range of spatial scales in this massive star-forming region. Very Large Array observations of NH_3 at 4'' angular resolution show a large-scale (0.5 pc) molecular accretion flow around and into a star cluster with three small, bright HII regions. Higher resolution (0.4'') observations with the Submillimeter Array in hot core molecules (CH_3CN, OCS, and SO_2) and the VLA in NH_3, show that the two brightest and smallest HII regions are themselves surrounded by smaller scale (0.05 pc) accretion flows. The axes of rotation of the large- and small-scale flows are aligned, and the timescale for the contraction of the cloud is short enough, 0.1 Myr, for the large-scale accretion flow to deliver significant mass to the smaller scales within the star formation timescale. The flow structure appears to be continuous and hierarchical from larger to smaller scales. Millimeter radio recombination line (RRL) observations at 0.4" angular resolution indicate rotation and outflow of the ionized gas within the brightest HII region (A). The broad recombination lines and a continuum spectral energy distribution (SED) that rises continuously from cm to mm wavelengths, are both characteristic of the class of HII regions known as "broad recombination line objects". The SED indicates a density gradient inside this HII region, and the RRLs suggest supersonic flows. These observations are consistent with photoevaporation of the inner part of the rotationally flattened molecular accretion flow. We also report the serendipitous detection of a new NH_3 (3,3) maser.

preprint2007arXiv

The formation of massive stars: accretion, disks and the development of hypercompact HII regions

The hypothesis that massive stars form by accretion can be investigated by simple analytical calculations that describe the effect that the formation of a massive star has on its own accretion flow. Within a simple accretion model that includes angular momentum, that of gas flow on ballistic trajectories around a star, the increasing ionization of a massive star growing by accretion produces a three-stage evolutionary sequence. The ionization first forms a small quasi-spherical HII region gravitationally trapped within the accretion flow. At this stage the flow of ionized gas is entirely inward. As the ionization increases, the HII region transitions to a bipolar morphology in which the inflow is replaced by outflow within a narrow range of angle with about the bipolar axis. At higher rates of ionization, the opening angle of the outflow region progressively increases. Eventually, in the third stage, the accretion is confined to a thin region about an equatorial disk. Throughout this early evolution, the HII region is of hypercompact to ultracompact size depending on the mass of the enclosed star or stars. These small HII regions whose dynamics are dominated by stellar gravitation and accretion are different than compact and larger HII regions whose dynamics are dominated by the thermal pressure of the ionized gas.