Showing 1–50 of 56 results for author: Yorke, H W

Makemake + Sedna: A Continuum Radiation Transport and Photoionization Framework for Astrophysical Newtonian Fluid Dynamics

Authors: Rolf Kuiper, Harold W. Yorke, Andrea Mignone

Abstract: Astrophysical fluid flow studies often encompass a wide range of physical processes to account for the complexity of the system under consideration. In addition to gravity, a proper treatment of thermodynamic processes via continuum radiation transport and/or photoionization is becoming the state of the art. We present a major update of our continuum radiation transport module, MAKEMAKE, and a new… ▽ More Astrophysical fluid flow studies often encompass a wide range of physical processes to account for the complexity of the system under consideration. In addition to gravity, a proper treatment of thermodynamic processes via continuum radiation transport and/or photoionization is becoming the state of the art. We present a major update of our continuum radiation transport module, MAKEMAKE, and a newly developed module for photoionization, SEDNA, coupled to the magnetohydrodynamics code PLUTO. These extensions are currently not publicly available; access can be granted on a case-by-case basis. We explain the theoretical background of the equations solved, elaborate on the numerical layout, and present a comprehensive test suite for radiation-ionization hydrodynamics. The grid-based radiation and ionization modules support static one-dimensional, two-dimensional, and three-dimensional grids in Cartesian, cylindrical, and spherical coordinates. Each module splits the radiation field into two components, one originating directly from a point source - solved using a ray-tracing scheme - and a diffuse component - solved with a three-dimensional flux-limited diffusion (FLD) solver. The FLD solver for the continuum radiation transport makes use of either the equilibrium one-temperature approach or the linearization two-temperature approach. The FLD solver for the photoionization module enables accounting for the temporal evolution of the radiation field from direct recombination of free electrons into hydrogen's ground state as an alternative to the on-the-spot approximation. A brief overview of completed and ongoing scientific studies is given to explicitly illustrate the multipurpose nature of the numerical framework presented. △ Less

Submitted 25 September, 2020; originally announced September 2020.

Journal ref: The Astrophysical Journal Supplement Series, 2020, Volume 250, Issue 1, id.13, 28 pp

Python Radiative Transfer Emission code (PyRaTE): non-LTE spectral lines simulations

Authors: Aris Tritsis, Harold W. Yorke, Konstantinos Tassis

Abstract: We describe PyRaTE, a new, non-local thermodynamic equilibrium (non-LTE) line radiative transfer code developed specifically for post-processing astrochemical simulations. Population densities are estimated using the escape probability method. When computing the escape probability, the optical depth is calculated towards all directions with density, molecular abundance, temperature and velocity va… ▽ More We describe PyRaTE, a new, non-local thermodynamic equilibrium (non-LTE) line radiative transfer code developed specifically for post-processing astrochemical simulations. Population densities are estimated using the escape probability method. When computing the escape probability, the optical depth is calculated towards all directions with density, molecular abundance, temperature and velocity variations all taken into account. A very easy-to-use interface, capable of importing data from simulations outputs performed with all major astrophysical codes, is also developed. The code is written in Python using an `embarrassingly parallel' strategy and can handle all geometries and projection angles. We benchmark the code by comparing our results with those from RADEX (van der Tak et al. 2007) and against analytical solutions and present case studies using hydrochemical simulations. The code is available on GitHub (https://github.com/ArisTr/PyRaTE). △ Less

Submitted 23 May, 2018; originally announced May 2018.

Comments: 10 pages, 13 figures. Accepted for publication in MNRAS. Code link: https://github.com/ArisTr/PyRaTE

Science Impacts of the SPHEREx All-Sky Optical to Near-Infrared Spectral Survey II: Report of a Community Workshop on the Scientific Synergies Between the SPHEREx Survey and Other Astronomy Observatories

Authors: Olivier Doré, Michael W. Werner, Matthew L. N. Ashby, Lindsey E. Bleem, Jamie Bock, Jennifer Burt, Peter Capak, Tzu-Ching Chang, Jonás Chaves-Montero, Christine H. Chen, Francesca Civano, I. Ilsedore Cleeves, Asantha Cooray, Brendan Crill, Ian J. M. Crossfield, Michael Cushing, Sylvain de la Torre, Tiziana DiMatteo, Niv Dvory, Cora Dvorkin, Catherine Espaillat, Simone Ferraro, Douglas Finkbeiner, Jenny Greene, Jackie Hewitt , et al. (38 additional authors not shown)

Abstract: SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$μ$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$σ$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx… ▽ More SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$μ$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$σ$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop. △ Less

Submitted 24 May, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: 50 pages, 24 figures, more details at http://spherex.caltech.edu

Protostellar Outflows and Radiative Feedback from Massive Stars. II. Feedback, Star Formation Efficiency, and Outflow Broadening

Authors: Rolf Kuiper, Neal J. Turner, Harold W. Yorke

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, 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 mas… ▽ More 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, 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 are 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 degree during early protostellar accretion phases, increasing rapidly up to 65 degree at the onset of radiation pressure feedback. △ Less

Submitted 16 September, 2016; originally announced September 2016.

Comments: accepted for publication at ApJ

Supermassive star formation via episodic accretion: protostellar disc instability and radiative feedback efficiency

Authors: Yuya Sakurai, Eduard I. Vorobyov, Takashi Hosokawa, Naoki Yoshida, Kazuyuki Omukai, Harold W. Yorke

Abstract: The formation of SMSs is a potential pathway to seed SMBHs in the early universe. A critical issue for forming SMSs is stellar UV feedback, which may limit the stellar mass growth via accretion. In this paper we study the evolution of an accreting SMS and its UV emissivity under conditions of realistic variable accretion from a self-gravitating circumstellar disc. First we conduct a 2D hydrodynami… ▽ More The formation of SMSs is a potential pathway to seed SMBHs in the early universe. A critical issue for forming SMSs is stellar UV feedback, which may limit the stellar mass growth via accretion. In this paper we study the evolution of an accreting SMS and its UV emissivity under conditions of realistic variable accretion from a self-gravitating circumstellar disc. First we conduct a 2D hydrodynamical simulation to follow the long-term protostellar accretion until the stellar mass exceeds $10^4~M_\odot$. The disc fragments due to gravitational instability, creating a number of small clumps that rapidly migrate inward to fall onto the star. The resulting accretion history is thus highly time-dependent: short episodic accretion bursts are followed by longer, relative quiescent phases. We show that the circumstellar disc for the so-called direct collapse model is more unstable and generates greater variability over shorter timescales than normal Pop III cases. We conduct a post-process stellar evolution calculation using the obtained accretion history. Our results show that, regardless of the strong variability of the accretion rates, the stellar radius monotonically increases with almost constant effective temperature at $T_{\rm eff} \simeq 5000$ K as the stellar mass increases. The resulting UV feedback is too weak to hinder mass accretion due to the low flux of stellar UV photons, thus verifying our implicit assumption of no stellar feedback during the hydrodynamic simulations. The insensitivity of stellar evolution to variable accretion is attributed to the fact that typical timescales of variability, $\lesssim 10^3$ years, are too short to affect the stellar structure. We argue that this evolution will continue until the SMS eventually collapses to produce a massive black hole by the general relativistic instability after the stellar mass reaches $\gtrsim 10^5~M_\odot$. △ Less

Submitted 19 November, 2015; originally announced November 2015.

Comments: 9 pages, 6 figures, submitted to MNRAS

Formation of Massive Primordial Stars: Intermittent UV Feedback with Episodic Mass Accretion

Authors: Takashi Hosokawa, Shingo Hirano, Rolf Kuiper, Harold W. Yorke, Kazuyuki Omukai, Naoki Yoshida

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 Msun $\lesssim M_* \lesssim$ 1000 Msun.… ▽ More 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 Msun $\lesssim M_* \lesssim$ 1000 Msun. The formation of very massive ($\gtrsim$ 250 Msun) stars is possible under weak UV feedback, whereas ordinary massive (a few x 10 Msun) 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. (2014), 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 Msun/yr, 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 HII 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. △ Less

Submitted 25 April, 2016; v1 submitted 5 October, 2015; originally announced October 2015.

Comments: 28 pages, 27 figures, accepted for publication in ApJ. Movies are available at https://www.youtube.com/playlist?list=PLy0BOLTBcHhb-Q6cG8hwmjUEvTzcH4AMF

Formation of primordial supermassive stars by burst accretion

Authors: Y. Sakurai, T. Hosokawa, N. Yoshida, H. W. Yorke

Abstract: A promising formation channel of SMBHs at redshift 6 is the so-called DC model, which posits that a massive seed BH forms through gravitational collapse of a $\sim 10^5~M_\odot$ SMS. We study the evolution of such a SMS growing by rapid mass accretion. In particular, we examine the impact of time-dependent mass accretion of repeating burst and quiescent phases that are expected to occur with a sel… ▽ More A promising formation channel of SMBHs at redshift 6 is the so-called DC model, which posits that a massive seed BH forms through gravitational collapse of a $\sim 10^5~M_\odot$ SMS. We study the evolution of such a SMS growing by rapid mass accretion. In particular, we examine the impact of time-dependent mass accretion of repeating burst and quiescent phases that are expected to occur with a self-gravitating circumstellar disk. We show that the stellar evolution with such episodic accretion differs qualitatively from that expected with a constant accretion rate, even if the mean accretion rate is the same. Unlike the case of constant mass accretion, whereby the star expands roughly following $R_* \simeq 2.6 \times 10^3 R_\odot (M_*/100~M_\odot)^{1/2}$, the protostar can substantially contract during the quiescent phases between accretion bursts. The stellar effective temperature and ionizing photon emissivity increase accordingly as the star contracts, which can cause strong ionizing feedback and halt the mass accretion onto the star. With a fixed duration of the quiescent phase $Δt_{\rm q}$, such contraction occurs in early evolutionary phases, i.e. for $M_* \lesssim 10^3~M_\odot$ with $Δt_{\rm q} \simeq 10^3$ yr. For later epochs and larger masses but the same $Δt_{\rm q}$, contraction is negligible even during quiescent phases. With larger quiescent times $Δt_{\rm q}$, however, the star continues to contract during quiescent phases even for the higher stellar masses. We show that such behavior is well understood by comparing the interval time and the thermal relaxation time for a bloated surface layer. We conclude that the UV radiative feedback becomes effective if the quiescent phase associated by the burst accretion is longer than $\sim 10^3$ yr, which is possible in an accretion disk forming in the direct collapse model. △ Less

Submitted 14 May, 2015; originally announced May 2015.

Comments: 10 pages, 7 figures, submitted to MNRAS

Early structure formation from primordial density fluctuations with a blue-tilted power spectrum

Authors: Shingo Hirano, Nick Zhu, Naoki Yoshida, David Spergel, Harold W. Yorke

Abstract: While observations of large-scale structure and the cosmic microwave background (CMB) provide strong constraints on the amplitude of the primordial power spectrum (PPS) on scales larger than 10~Mpc, the amplitude of the power spectrum on sub-galactic length scales is much more poorly constrained. We study early structure formation in a cosmological model with a blue-tilted PPS. We assume that the… ▽ More While observations of large-scale structure and the cosmic microwave background (CMB) provide strong constraints on the amplitude of the primordial power spectrum (PPS) on scales larger than 10~Mpc, the amplitude of the power spectrum on sub-galactic length scales is much more poorly constrained. We study early structure formation in a cosmological model with a blue-tilted PPS. We assume that the standard scale-invariant PPS is modified at small length scales as $P(k) \sim k^{m_{\rm s}}$ with $m_{\rm s} > 1$. We run a series of cosmological hydrodynamic simulations to examine the dependence of the formation epoch and the characteristic mass of primordial stars on the tilt of the PPS. In models with $m_{\rm s} > 1$, star-forming gas clouds are formed at $z > 100$, when formation of hydrogen molecules is inefficient because the intense CMB radiation destroys chemical intermediates. Without efficient coolant, the gas clouds gravitationally contract while keeping a high temperature. The protostars formed in such "hot" clouds grow very rapidly by accretion to become extremely massive stars that may leave massive black holes with a few hundred solar-masses at $z > 100$. The shape of the PPS critically affects the properties and the formation epoch of the first generation of stars. Future experiments of the CMB polarization and the spectrum distortion may provide important information on the nature of the first stars and their formation epoch, and hence on the shape of the small-scale power spectrum. △ Less

Submitted 8 October, 2015; v1 submitted 20 April, 2015; originally announced April 2015.

Comments: 8 pages, 5 figures, 1 table, accepted for publication in ApJ

Primordial Star Formation under the Influence of Far Ultraviolet Radiation: 1540 Cosmological Halos and the Stellar Mass Distribution

Authors: Shingo Hirano, Takashi Hosokawa, Naoki Yoshida, Kazuyuki Omukai, Harold W. Yorke

Abstract: We perform a large set of cosmological simulations of early structure formation and follow the formation and evolution of 1540 star-forming gas clouds to derive the mass distribution of primordial stars. The star formation in our cosmological simulations is characterized by two distinct populations, the so-called Population III.1 stars and primordial stars formed under the influence of far ultravi… ▽ More We perform a large set of cosmological simulations of early structure formation and follow the formation and evolution of 1540 star-forming gas clouds to derive the mass distribution of primordial stars. The star formation in our cosmological simulations is characterized by two distinct populations, the so-called Population III.1 stars and primordial stars formed under the influence of far ultraviolet (FUV) radiation (Population III.2D stars). In this work, we determine the stellar masses by using the dependences on the physical properties of star-forming cloud and/or the external photodissociating intensity from nearby primordial stars, which are derived from the results of two-dimensional radiation hydrodynamic simulations of protostellar feedback. The characteristic mass of the Pop III stars is found to be a few hundred solar masses at z ~ 25, and it gradually shifts to lower masses with decreasing redshift. At high redshifts z > 20, about half of the star-forming gas clouds are exposed to intense FUV radiation and thus give birth to massive Pop III.2D stars. However, the local FUV radiation by nearby Pop III stars becomes weaker at lower redshifts, when typical Pop III stars have smaller masses and the mean physical separation between the stars becomes large owing to cosmic expansion. Therefore, at z < 20, a large fraction of the primordial gas clouds host Pop III.1 stars. At z =< 15, the Pop III.1 stars are formed in relatively cool gas clouds due to efficient radiative cooling by H_2 and HD molecules; such stars have masses of a few x 10 Msun. Since the stellar evolution and the final fate are determined by the stellar mass, Pop III stars formed at different epochs play different roles in the early universe. △ Less

Submitted 7 February, 2015; v1 submitted 7 January, 2015; originally announced January 2015.

Comments: 22 pages, 19 figures, published in MNRAS

Journal ref: MNRAS, 2015, Volume 448, Issue 1, Pages. 568-587

Protostellar Outflows and Radiative Feedback from Massive Stars

Authors: Rolf Kuiper, Harold W. Yorke, Neal J. Turner

Abstract: We carry out radiation hydrodynamical simulations of the formation of massive stars in the super-Eddington regime including both their radiative feedback and protostellar outflows. The calculations start from a prestellar core of dusty gas and continue until the star stops growing. The accretion ends when the remnants of the core are ejected, mostly by the force of the direct stellar radiation in… ▽ More We carry out radiation hydrodynamical simulations of the formation of massive stars in the super-Eddington regime including both their radiative feedback and protostellar outflows. The calculations start from a prestellar core of dusty gas and continue until the star stops growing. The accretion ends when the remnants of the core are ejected, mostly by the force of the direct stellar radiation in the polar direction and elsewhere by the reradiated thermal infrared radiation. How long the accretion persists depends on whether the protostellar outflows are present. We set the mass outflow rate to 1% of the stellar sink particle's accretion rate. The outflows open a bipolar cavity extending to the core's outer edge, through which the thermal radiation readily escapes. The radiative flux is funneled into the polar directions while the core's collapse proceeds near the equator. The outflow thus extends the "flashlight effect", or anisotropic radiation field, found in previous studies from the few hundred AU scale of the circumstellar disk up to the 0.1 parsec scale of the core. The core's flashlight effect allows core gas to accrete on the disk for longer, in the same way that the disk's flashlight effect allows disk gas to accrete on the star for longer. Thus although the protostellar outflows remove material near the core's poles, causing slower stellar growth over the first few free-fall times, they also enable accretion to go on longer in our calculations. The outflows ultimately lead to stars of somewhat higher mass. △ Less

Submitted 19 December, 2014; originally announced December 2014.

Comments: accepted for publication at ApJ

Effect of OH depletion on measurements of the mass-to-flux ratio in molecular cloud cores

Authors: K. Tassis, K. Willacy, H. W. Yorke, N. J. Turner

Abstract: The ratio of mass and magnetic flux determines the relative importance of magnetic and gravitational forces in the evolution of molecular clouds and their cores. Its measurement is thus central in discriminating between different theories of core formation and evolution. Here we discuss the effect of chemical depletion on measurements of the mass-to-flux ratio using the same molecule (OH) both for… ▽ More The ratio of mass and magnetic flux determines the relative importance of magnetic and gravitational forces in the evolution of molecular clouds and their cores. Its measurement is thus central in discriminating between different theories of core formation and evolution. Here we discuss the effect of chemical depletion on measurements of the mass-to-flux ratio using the same molecule (OH) both for Zeeman measurements of the magnetic field and the determination of the mass of the region. The uncertainties entering through the OH abundance in determining separately the magnetic field and the mass of a region have been recognized in the literature. It has been proposed however that, when comparing two regions of the same cloud, the abundance will in both cases be the same. We show that this assumption is invalid. We demonstrate that when comparing regions with different densities, the effect of OH depletion in measuring changes of the mass-to-flux ratio between different parts of the same cloud can even reverse the direction of the underlying trends (for example, the mass-to-flux ratio may appear to decrease as we move to higher density regions). The systematic errors enter primarily through the inadequate estimation of the mass of the region. △ Less

Submitted 25 August, 2014; originally announced August 2014.

Comments: 5 pages, 3 figures, accepted for publication in MNRAS

Formation of Primordial Supermassive Stars by Rapid Mass Accretion

Authors: Takashi Hosokawa, Harold W. Yorke, Kohei Inayoshi, Kazuyuki Omukai, Naoki Yoshida

Abstract: Supermassive stars (SMSs) forming via very rapid mass accretion (Mdot >~ 0.1 Msun/yr) could be precursors of supermassive black holes observed beyond redshift of about 6. Extending our previous work, we here study the evolution of primordial stars growing under such rapid mass accretion until the stellar mass reaches 10^{4 - 5} Msun. Our stellar evolution calculations show that a star becomes supe… ▽ More Supermassive stars (SMSs) forming via very rapid mass accretion (Mdot >~ 0.1 Msun/yr) could be precursors of supermassive black holes observed beyond redshift of about 6. Extending our previous work, we here study the evolution of primordial stars growing under such rapid mass accretion until the stellar mass reaches 10^{4 - 5} Msun. Our stellar evolution calculations show that a star becomes supermassive while passing through the "supergiant protostar" stage, whereby the star has a very bloated envelope and a contracting inner core. The stellar radius increases monotonically with the stellar mass, until =~ 100 AU for M_* >~ 10^4 Msun, after which the star begins to slowly contract. Because of the large radius the effective temperature is always less than 10^4 K during rapid accretion. The accreting material is thus almost completely transparent to the stellar radiation. Only for M_* >~ 10^5 Msun can stellar UV feedback operate and disturb the mass accretion flow. We also examine the pulsation stability of accreting SMSs, showing that the pulsation-driven mass loss does not prevent stellar mass growth. Observational signatures of bloated SMSs should be detectable with future observational facilities such as the James Webb Space Telescope. Our results predict that an inner core of the accreting SMS should suffer from the general relativistic instability soon after the stellar mass exceeds 10^5 Msun. An extremely massive black hole should form after the collapse of the inner core. △ Less

Submitted 14 October, 2013; v1 submitted 20 August, 2013; originally announced August 2013.

Comments: 14 pages, 13 figures, accepted for publication in ApJ

One Hundred First Stars : Protostellar Evolution and the Final Masses

Authors: Shingo Hirano, Takashi Hosokawa, Naoki Yoshida, Hideyuki Umeda, Kazuyuki Omukai, Gen Chiaki, Harold W. Yorke

Abstract: We perform a large set of radiation hydrodynamics simulations of primordial star formation in a fully cosmological context. Our statistical sample of 100 First Stars show that the first generation of stars have a wide mass distribution M_popIII = 10 ~ 1000 M_sun. We first run cosmological simulations to generate a set of primordial star-forming gas clouds. We then follow protostar formation in eac… ▽ More We perform a large set of radiation hydrodynamics simulations of primordial star formation in a fully cosmological context. Our statistical sample of 100 First Stars show that the first generation of stars have a wide mass distribution M_popIII = 10 ~ 1000 M_sun. We first run cosmological simulations to generate a set of primordial star-forming gas clouds. We then follow protostar formation in each gas cloud and the subsequent protostellar evolution until the gas mass accretion onto the protostar is halted by stellar radiative feedback. The accretion rates differ significantly among the primordial gas clouds which largely determine the final stellar masses. For low accretion rates the growth of a protostar is self-regulated by radiative feedback effects and the final mass is limited to several tens of solar masses. At high accretion rates the protostar's outer envelope continues to expand and the effective surface temperature remains low; such protostars do not exert strong radiative feedback and can grow in excess to one hundred solar masses. The obtained wide mass range suggests that the first stars play a variety of roles in the early universe, by triggering both core-collapse supernovae and pair-instability supernovae as well as by leaving stellar mass black holes. We find certain correlations between the final stellar mass and the physical properties of the star-forming cloud. These correlations can be used to estimate the mass of the first star from the properties of the parent cloud or of the host halo, without following the detailed protostellar evolution. △ Less

Submitted 19 February, 2014; v1 submitted 20 August, 2013; originally announced August 2013.

Comments: 25 pages, 24 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.781:60-81,2014

On the simultaneous evolution of massive protostars and their host cores

Authors: Rolf Kuiper, Harold W. Yorke

Abstract: Studies of the evolution of massive protostars and the evolution of their host molecular cloud cores are commonly treated as separate problems. However, interdependencies between the two can be significant. Here, we study the simultaneous evolution of massive protostars and their host molecular cores using a multi-dimensional radiation hydrodynamics code that incorporates the effects of the therma… ▽ More Studies of the evolution of massive protostars and the evolution of their host molecular cloud cores are commonly treated as separate problems. However, interdependencies between the two can be significant. Here, we study the simultaneous evolution of massive protostars and their host molecular cores using a multi-dimensional radiation hydrodynamics code that incorporates the effects of the thermal pressure and radiative acceleration feedback of the centrally forming protostar. The evolution of the massive protostar is computed simultaneously using the stellar evolution code STELLAR, modified to include the effects of variable accretion. The interdependencies are studied in three different collapse scenarios. For comparison, stellar evolutionary tracks at constant accretion rates and the evolution of the host cores using pre-computed stellar evolutionary tracks are computed. The resulting interdependencies of the protostellar evolution and the evolution of the environment are extremely diverse and depend on the order of events, in particular the time of circumstellar accretion disk formation with respect to the onset of the bloating phase of the star. Feedback mechanisms affect the instantaneous accretion rate and the protostar's radius, temperature and luminosity on timescales equal or smaller than 5 kyr, corresponding to the accretion timescale and Kelvin-Helmholtz contraction timescale, respectively. Nevertheless, it is possible to approximate the overall protostellar evolution in many cases by pre-computed stellar evolutionary tracks assuming appropriate constant average accretion rates. △ Less

Submitted 27 May, 2013; originally announced May 2013.

Comments: accepted for publication at ApJ

On the effects of optically thick gas (disks) around massive stars

Authors: Rolf Kuiper, Harold W. Yorke

Abstract: Numerical simulations have shown that the often cited radiation pressure barrier to accretion onto massive stars can be circumvented, when the radiation field is highly anisotropic in the presence of a circumstellar accretion disk with high optical depth. Here, these studies of the so-called flashlight effect are expanded by including the opacity of the innermost dust-free but potentially opticall… ▽ More Numerical simulations have shown that the often cited radiation pressure barrier to accretion onto massive stars can be circumvented, when the radiation field is highly anisotropic in the presence of a circumstellar accretion disk with high optical depth. Here, these studies of the so-called flashlight effect are expanded by including the opacity of the innermost dust-free but potentially optically thick gas regions around forming massive stars. In addition to frequency-dependent opacities for the dust grains, we use temperature- and density-dependent Planck- and Rosseland mean opacities for the gas. The simulations show that the innermost dust-free parts of the accretion disks are optically thick to the stellar radiation over a substantial fraction of the solid angle above and below the disk's midplane. The temperature in the shielded disk region decreases faster with radius than in a comparison simulation with a lower constant gas opacity, and the dust sublimation front is shifted to smaller radii. The shielding by the dust-free gas in the inner disk thus contributes to an enhanced flashlight effect, which ultimately results in a smaller opening angle of the radiation pressure driven outflow and in a much longer timescale of sustained feeding of the circumstellar disk by the molecular cloud core. We conclude that it is necessary to properly account for the opacity of the inner dust-free disk regions around forming massive stars in order to correctly assess the effectiveness of the flashlight effect, the opening angle of radiation pressure driven outflows, and the lifetime and morphological evolution of the accretion disk. △ Less

Submitted 27 November, 2012; originally announced November 2012.

Comments: accepted at ApJ

Journal ref: ApJ (2013), vol. 763 (2), pp. 104

Protostellar Feedback and Final Mass of the Second-Generation Primordial Stars

Authors: Takashi Hosokawa, Naoki Yoshida, Kazuyuki Omukai, Harold W. Yorke

Abstract: The first stars in the universe ionized the ambient primordial gas through various feedback processes. "Second-generation" primordial stars potentially form from this disturbed gas after its recombination. In this Letter, we study the late formation stage of such second-generation stars, where a large amount of gas accretes onto the protostar and the final stellar mass is determined when the accre… ▽ More The first stars in the universe ionized the ambient primordial gas through various feedback processes. "Second-generation" primordial stars potentially form from this disturbed gas after its recombination. In this Letter, we study the late formation stage of such second-generation stars, where a large amount of gas accretes onto the protostar and the final stellar mass is determined when the accretion terminates. We directly compute the complex interplay between the accretion flow and stellar ultraviolet (UV) radiation, performing radiation-hydrodynamic simulations coupled with stellar evolution calculations. Because of more efficient H2 and HD cooling in the pre-stellar stage, the accretion rates onto the star are ten times lower than in the case of the formation of the first stars. The lower accretion rates and envelope density result in the occurrence of an expanding bipolar HII region at a lower protostellar mass M_* \simeq 10Msun, which blows out the circumstellar material, thereby quenching the mass supply from the envelope to the accretion disk. At the same time the disk loses mass due to photoevaporation by the growing star. In our fiducial case the stellar UV feedback terminates mass accretion onto the star at M_* \simeq 17Msun. Although the derived masses of the second-generation primordial stars are systematically lower than those of the first generation, the difference is within a factor of only a few. Our results suggest a new scenario, whereby the majority of the primordial stars are born as massive stars with tens of solar masses, regardless of their generations. △ Less

Submitted 10 October, 2012; originally announced October 2012.

Comments: 5 pages, 4 figures, to be published in ApJL

Rapidly Accreting Supergiant Protostars: Embryos of Supermassive Black Holes?

Authors: Takashi Hosokawa, Kazuyuki Omukai, Harold W. Yorke

Abstract: Direct collapse of supermassive stars (SMSs) is a possible pathway for generating supermassive black holes in the early universe. It is expected that an SMS could form via very rapid mass accretion with Mdot ~ 0.1 - 1 Msun/yr during the gravitational collapse of an atomic-cooling primordial gas cloud. In this paper we study how stars would evolve under such extreme rapid mass accretion, focusing o… ▽ More Direct collapse of supermassive stars (SMSs) is a possible pathway for generating supermassive black holes in the early universe. It is expected that an SMS could form via very rapid mass accretion with Mdot ~ 0.1 - 1 Msun/yr during the gravitational collapse of an atomic-cooling primordial gas cloud. In this paper we study how stars would evolve under such extreme rapid mass accretion, focusing on the early evolution until the stellar mass reaches 1000 Msun. To this end we numerically calculate the detailed interior structure of accreting stars with primordial element abundances. Our results show that for accretion rates higher than 0.01 Msun/yr, stellar evolution is qualitatively different from that expected at lower rates. While accreting at these high rates the star always has a radius exceeding 100 Rsun, which increases monotonically with the stellar mass. The mass-radius relation for stellar masses exceeding ~ 100 Msun follows the same track with R_* \propto M_*^0.5 in all cases with accretion rates > 0.01 Msun/yr; at a stellar mass of 1000 Msun the radius is about 7000 Rsun (~= 30 AU). With higher accretion rates the onset of hydrogen burning is shifted towards higher stellar masses. In particular, for accretion rates exceeding Mdot > 0.1 Msun/yr, there is no significant hydrogen burning even after 1000 Msun have accreted onto the protostar. Such "supergiant" protostars have effective temperatures as low as Teff ~= 5000 K throughout their evolution and because they hardly emit ionizing photons, they do not create an HII region or significantly heat their immediate surroundings. Thus, radiative feedback is unable to hinder the growth of rapidly accreting stars to masses in excess of 1000 Msun, as long as material is accreted at rates Mdot > 0.01 Msun/yr. △ Less

Submitted 12 March, 2012; originally announced March 2012.

Comments: 11 pages, 10 figures

Non-Equilibrium Chemistry of Dynamically Evolving Prestellar Cores: II. Ionization and Magnetic Field

Authors: Konstantinos Tassis, Karen Willacy, Harold W. Yorke, Neal Turner

Abstract: We study the effect that non-equilibrium chemistry in dynamical models of collapsing molecular cloud cores has on measurements of the magnetic field in these cores, the degree of ionization, and the mean molecular weight of ions. We find that OH and CN, usually used in Zeeman observations of the line-of-sight magnetic field, have an abundance that decreases toward the center of the core much faste… ▽ More We study the effect that non-equilibrium chemistry in dynamical models of collapsing molecular cloud cores has on measurements of the magnetic field in these cores, the degree of ionization, and the mean molecular weight of ions. We find that OH and CN, usually used in Zeeman observations of the line-of-sight magnetic field, have an abundance that decreases toward the center of the core much faster than the density increases. As a result, Zeeman observations tend to sample the outer layers of the core and consistently underestimate the core magnetic field. The degree of ionization follows a complicated dependence on the number density at central densities up to 10^5 cm^{-3} for magnetic models and 10^6 cm^{-3} in non-magnetic models. At higher central densities the scaling approaches a power-law with a slope of -0.6 and a normalization which depends on the cosmic-ray ionization rate ζ and the temperature T as (ζT)^1/2. The mean molecular weight of ions is systematically lower than the usually assumed value of 20 - 30, and, at high densities, approaches a value of 3 due to the asymptotic dominance of the H3+ ion. This significantly lower value implies that ambipolar diffusion operates faster. △ Less

Submitted 22 May, 2012; v1 submitted 17 November, 2011; originally announced November 2011.

Comments: 7 pages, 5 figures, accepted for publication in ApJ

Non-Equilibrium Chemistry of Dynamically Evolving Prestellar Cores: I. Basic Magnetic and Non-Magnetic Models and Parameter Studies

Authors: Konstantinos Tassis, Karen Willacy, Harold W. Yorke, Neal J. Turner

Abstract: We combine dynamical and non-equilibrium chemical modeling of evolving prestellar molecular cloud cores, and explore the evolution of molecular abundances in the contracting core. We model both magnetic cores, with varying degrees of initial magnetic support, and non-magnetic cores, with varying collapse delay times. We explore, through a parameter study, the competing effects of various model par… ▽ More We combine dynamical and non-equilibrium chemical modeling of evolving prestellar molecular cloud cores, and explore the evolution of molecular abundances in the contracting core. We model both magnetic cores, with varying degrees of initial magnetic support, and non-magnetic cores, with varying collapse delay times. We explore, through a parameter study, the competing effects of various model parameters in the evolving molecular abundances, including the elemental C/O ratio, the temperature, and the cosmic-ray ionization rate. We find that different models show their largest quantitative differences at the center of the core, whereas the outer layers, which evolve slower, have abundances which are severely degenerate among different dynamical models. There is a large range of possible abundance values for different models at a fixed evolutionary stage (central density), which demonstrates the large potential of chemical differentiation in prestellar cores. However, degeneracies among different models, compounded with uncertainties induced by other model parameters, make it difficult to discriminate among dynamical models. To address these difficulties, we identify abundance ratios between particular molecules, the measurement of which would have maximal potential for discrimination among the different models examined here. In particular, we find that the ratios between NH3 and CO; NH2 and CO; NH3 and HCO+ are sensitive to the evolutionary timescale, and that the ratio between HCN and OH is sensitive to the C/O ratio. Finally, we demonstrate that measurements of the central deviation (central depletion or enhancement) of abundances of certain molecules are good indicators of the dynamics of the core. △ Less

Submitted 26 April, 2012; v1 submitted 16 November, 2011; originally announced November 2011.

Comments: 20 pages, 15 figures, accepted for publication in ApJ

Protostellar Feedback Halts the Growth of the First Stars in the Universe

Authors: Takashi Hosokawa, Kazuyuki Omukai, Naoki Yoshida, Harold W. Yorke

Abstract: The first stars fundamentally transformed the early universe by emitting the first light and by producing the first heavy elements. These effects were predetermined by the mass distribution of the first stars, which is thought to have been fixed by a complex interplay of gas accretion and protostellar radiation. We performed radiation-hydrodynamics simulations that followed the growth of a primord… ▽ More The first stars fundamentally transformed the early universe by emitting the first light and by producing the first heavy elements. These effects were predetermined by the mass distribution of the first stars, which is thought to have been fixed by a complex interplay of gas accretion and protostellar radiation. We performed radiation-hydrodynamics simulations that followed the growth of a primordial protostar through to the early stages as a star with thermo-nuclear burning. The circumstellar accretion disk was evaporated by ultraviolet radiation from the star when its mass was 43 times that of the Sun. Such massive primordial stars, in contrast to the often postulated extremely massive stars, may help explain the fact that there are no signatures of the pair-instability supernovae in abundance patterns of metal-poor stars in our galaxy. △ Less

Submitted 15 November, 2011; originally announced November 2011.

Comments: to appear in Science (published in online Science Express), combined with SOM, additional images and movies are available at http://www-tap.scphys.kyoto-u.ac.jp/~hosokawa/firststarstop_e.html

A New Recipe for Obtaining Central Volume Densities of Prestellar Cores from Size Measurements

Authors: Konstantinos Tassis, Harold W. Yorke

Abstract: We propose a simple analytical method for estimating the central volume density of prestellar molecular cloud cores from their column density profiles. Prestellar cores feature a flat central part of the column density and volume density profiles of the same size indicating the existence of a uniform density inner region. The size of this region is set by the thermal pressure force which depends o… ▽ More We propose a simple analytical method for estimating the central volume density of prestellar molecular cloud cores from their column density profiles. Prestellar cores feature a flat central part of the column density and volume density profiles of the same size indicating the existence of a uniform density inner region. The size of this region is set by the thermal pressure force which depends only on the central volume density and temperature of the core, and can provide a direct measurement of the central volume density. Thus a simple length measurement can immediately yield a central density estimate independent of any dynamical model for the core and without the need for fitting. Using the radius at which the column density is 90% of the central value as an estimate of the size of the flat inner part of the column density profile yields an estimate of the central volume density within a factor of 2 for well resolved cores. △ Less

Submitted 12 May, 2011; originally announced May 2011.

Comments: 4 pages, 4 figures, accepted for publication in ApJ Letters

Herschel observations of Extra-Ordinary Sources: Methanol as a probe of physical conditions in Orion KL

Authors: S. Wang, E. A. Bergin, N. R. Crockett, P. F. Goldsmith, D. C. Lis, J. C. Pearson, P. Schilke, T. A. Bell, C. Comito, G. A. Blake, E. Caux, C. Ceccarelli, J. Cernicharo, F. Daniel, M. -L. Dubernet, M. Emprechtinger, P. Encrenaz, M. Gerin, T. F. Giesen, J. R. Goicoechea, H. Gupta, E. Herbst, C. Joblin, D. Johnstone, W. D. Langer , et al. (23 additional authors not shown)

Abstract: We have examined methanol emission from Orion KL with of the {\em Herschel}/HIFI instrument, and detected two methanol bands centered at 524 GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band 1a) is dominated by the isolated $Δ$J$=$0, K$=-4\rightarrow$-3, v$_t$$=$0 Q branch, and includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is do… ▽ More We have examined methanol emission from Orion KL with of the {\em Herschel}/HIFI instrument, and detected two methanol bands centered at 524 GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band 1a) is dominated by the isolated $Δ$J$=$0, K$=-4\rightarrow$-3, v$_t$$=$0 Q branch, and includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is dominated by the $Δ$J$=$0, K$=7\rightarrow$6, v$_t$$=$0 Q branch transitions which are mostly blended. We have used the isolated E-type v$_t$$=$0 methanol transitions to explore the physical conditions in the molecular gas. With HIFI's high velocity resolution, the methanol emission contributed by different spatial components along the line of sight toward Orion KL (hot core, low velocity flow, and compact ridge) can be distinguished and studied separately. The isolated transitions detected in these bands cover a broad energy range (upper state energy ranging from 80 K to 900 K), which provides a unique probe of the thermal structure in each spatial component. The observations further show that the compact ridge is externally heated. These observations demonstrate the power of methanol lines as probes of the physical conditions in warm regions in close proximity to young stars. △ Less

Submitted 4 January, 2011; originally announced January 2011.

Comments: 8 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

Detection of OH$^+$ and H$_2$O$^+$ towards Orion~KL

Authors: H. Gupta, P. Rimmer, J. C. Pearson, S. Yu, E. Herbst, N. Harada, E. A. Bergin, D. A. Neufeld, G. J. Melnick, R. Bachiller, W. Baechtold, T. A. Bell, G. A. Blake, E. Caux, C. Ceccarelli, J. Cernicharo, G. Chattopadhyay, C. Comito, S. Cabrit, N. R. Crockett, F. Daniel, E. Falgarone, M. C. Diez-Gonzalez, M. -L. Dubernet, N. Erickson , et al. (54 additional authors not shown)

Abstract: We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H… ▽ More We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H$_2$O$^+$ are observed purely in absorption, showing a narrow component at the source velocity of 9 kms$^{-1}$, and a broad blueshifted absorption similar to that reported recently for HF and {\it para}-H$_{2}^{18}$O, and attributed to the low velocity outflow of Orion KL. We estimate column densities of OH$^+$ and H$_2$O$^+$ for the 9 km s$^{-1}$ component of $9 \pm 3 \times 10^{12}$cm$^{-2}$ and $7 \pm 2 \times 10^{12}$cm$^{-2}$, and those in the outflow of $1.9 \pm 0.7 \times 10^{13}$cm$^{-2}$ and $1.0 \pm 0.3 \times 10^{13}$cm$^{-2}$. Upper limits of $2.4\times 10^{12}$cm$^{-2}$ and $8.7\times 10^{12}$cm$^{-2}$ were derived for the column densities of {\it ortho} and {\it para}-H$_3$O$^+$ from transitions near 985 and 1657GHz. The column densities of the three ions are up to an order of magnitude lower than those obtained from recent observations of W31C and W49N. The comparatively low column densities may be explained by a higher gas density despite the assumption of a very high ionization rate. △ Less

Submitted 8 September, 2010; originally announced September 2010.

Comments: 5 pages, 1 figure. Accepted to and to be published in the Herschel HIFI special issue of A\&A Letters

Report number: CL#10-3059

2-D Magnetohydrodynamic Simulations of Induced Plasma Dynamics in the Near-Core Region of a Galaxy Cluster

Authors: Ioannis G. Mikellides, Konstantinos Tassis, Harold W. Yorke

Abstract: We present results from numerical simulations of the cooling-core cluster A2199 produced by the two-dimensional (2-D) resistive magnetohydrodynamics (MHD) code MACH2. In our simulations we explore the effect of anisotropic thermal conduction on the energy balance of the system. The results from idealized cases in 2-D axisymmetric geometry underscore the importance of the initial plasma density in… ▽ More We present results from numerical simulations of the cooling-core cluster A2199 produced by the two-dimensional (2-D) resistive magnetohydrodynamics (MHD) code MACH2. In our simulations we explore the effect of anisotropic thermal conduction on the energy balance of the system. The results from idealized cases in 2-D axisymmetric geometry underscore the importance of the initial plasma density in ICM simulations, especially the near-core values since the radiation cooling rate is proportional to ${n_e}^2$. Heat conduction is found to be non-effective in preventing catastrophic cooling in this cluster. In addition we performed 2-D planar MHD simulations starting from initial conditions deliberately violating both thermal balance and hydrostatic equilibrium in the ICM, to assess contributions of the convective terms in the energy balance of the system against anisotropic thermal conduction. We find that in this case work done by the pressure on the plasma can dominate the early evolution of the internal energy over anisotropic thermal conduction in the presence of subsonic flows, thereby reducing the impact of the magnetic field. Deviations from hydrostatic equilibrium near the cluster core may be associated with transient activity of a central active galactic nucleus and/or remnant dynamical activity in the ICM and warrant further study in three dimensions. △ Less

Submitted 3 September, 2010; originally announced September 2010.

Comments: 16 pages, 13 figures, accepted for publication in MNRAS

Gas morphology and energetics at the surface of PDRs: new insights with Herschel observations of NGC 7023

Authors: C. Joblin, P. Pilleri, J. Montillaud, A. Fuente, M. Gerin, O. Berné, V. Ossenkopf, J. Le Bourlot, D. Teyssier, J. R. Goicoechea, F. Le Petit, M. Röllig, M. Akyilmaz, A. O. Benz, F. Boulanger, S. Bruderer, C. Dedes, K. France, R. Güsten, A. Harris, T. Klein, C. Kramer, S. D. Lord, P. G. Martin, J. Martin-Pintado , et al. (15 additional authors not shown)

Abstract: We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach co… ▽ More We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species. △ Less

Submitted 7 August, 2010; originally announced August 2010.

Comments: Accepted for publication in Astronomy and Astrophysics Letters (Herschel HIFI special issue), 5 pages, 5 figures

Herschel observations of EXtra-Ordinary Sources (HEXOS): Observations of H2O and its isotopologues towards Orion KL

Authors: G. J. Melnick, V. Tolls, D. A. Neufeld, E. A. Bergin, T. G. Phillips, S. Wang, N. R. Crockett, T. A. Bell, G. A. Blake, S. Cabrit, E. Caux, C. Ceccarelli, J. Cernicharo, C. Comito, F. Daniel, M. -L. Dubernet, M. Emprechtinger, P. Encrenaz, E. Falgarone, M. Gerin, T. F. Giesen, J. R. Goicoechea, P. F. Goldsmith, E. Herbst, C. Joblin , et al. (27 additional authors not shown)

Abstract: We report the detection of more than 48 velocity-resolved ground rotational state transitions of H2(16)O, H2(18)O, and H2(17)O - most for the first time - in both emission and absorption toward Orion KL using Herschel/HIFI. We show that a simple fit, constrained to match the known emission and absorption components along the line of sight, is in excellent agreement with the spectral profiles of al… ▽ More We report the detection of more than 48 velocity-resolved ground rotational state transitions of H2(16)O, H2(18)O, and H2(17)O - most for the first time - in both emission and absorption toward Orion KL using Herschel/HIFI. We show that a simple fit, constrained to match the known emission and absorption components along the line of sight, is in excellent agreement with the spectral profiles of all the water lines. Using the measured H2(18)O line fluxes, which are less affected by line opacity than their H2(16)O counterparts, and an escape probability method, the column densities of H2(18)O associated with each emission component are derived. We infer total water abundances of 7.4E-5, 1.0E-5, and 1.6E-5 for the plateau, hot core, and extended warm gas, respectively. In the case of the plateau, this value is consistent with previous measures of the Orion-KL water abundance as well as those of other molecular outflows. In the case of the hot core and extended warm gas, these values are somewhat higher than water abundances derived for other quiescent clouds, suggesting that these regions are likely experiencing enhanced water-ice sublimation from (and reduced freeze-out onto) grain surfaces due to the warmer dust in these sources. △ Less

Submitted 30 July, 2010; originally announced July 2010.

Comments: 5 pages, 2 figures

Reversal of infall in SgrB2(M) revealed by Herschel/HIFI observations of HCN lines at THz frequencies

Authors: Rainer Rolffs, Peter Schilke, Claudia Comito, E. A. Bergin, F. F. S. van der Tak, D. C. Lis, S. -L. Qin, K. M. Menten, R. Guesten, T. A. Bell, G. A. Blake, E. Caux, C. Ceccarelli, J. Cernicharo, N. R. Crockett, F. Daniel, M. -L. Dubernet, M. Emprechtinger, P. Encrenaz, M. Gerin, T. F. Giesen, J. R. Goicoechea, P. F. Goldsmith, H. Gupta, E. Herbst , et al. (32 additional authors not shown)

Abstract: To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to… ▽ More To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity (HCN to H$^{13}$CN). This is most evident in the HCN 12--11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field changes from infall in the outer part to expansion in the inner part. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HIFI were critically important. △ Less

Submitted 29 July, 2010; originally announced July 2010.

Comments: A&A, HIFI special issue, accepted

Herschel observations of deuterated water towards Sgr B2(M)

Authors: Claudia Comito, Peter Schilke, Rainer Rolffs, D. C. Lis, A. Belloche, E. A. Bergin, T. G. Phillips, T. A. Bell, N. R. Crockett, S. Wang, G. A. Blake, E. Caux, C. Ceccarelli, J. Cernicharo, F. Daniel, M. -L. Dubernet, M. Emprechtinger, P. Encrenaz, M. Gerin, T. F. Giesen, J. R. Goicoechea, P. F. Goldsmith, H. Gupta, E. Herbst, C. Joblin , et al. (31 additional authors not shown)

Abstract: Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before wit… ▽ More Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from a low abundance ($2.5\times 10^{-11}$) in the outer envelope at temperatures below 100~K through a medium abundance ($1.5\times 10^{-9}$) in the inner envelope/outer core, at temperatures between 100 and 200~K, and finally a high abundance ($3.5\times 10^{-9}$) at temperatures above 200~K in the hot core. △ Less

Submitted 29 July, 2010; originally announced July 2010.

Comments: A&A HIFI special issue, accepted

A Sample of [CII] Clouds Tracing Dense Clouds in Weak FUV Fields observed by Herschel

Authors: Jorge L. Pineda, Thangasamy Velusamy, William D. Langer, Paul F. Goldsmith, Di Li., Harold W. Yorke, Jet Propulsion Laboratory

Abstract: The [CII] fine--structure line at 158um is an excellent tracer of the warm diffuse gas in the ISM and the interfaces between molecular clouds and their surrounding atomic and ionized envelopes. Here we present the initial results from Galactic Observations of Terahertz C+ (GOTC+), a Herschel Key Project devoted to study the [CII] fine structure emission in the galactic plane using the HIFI instrum… ▽ More The [CII] fine--structure line at 158um is an excellent tracer of the warm diffuse gas in the ISM and the interfaces between molecular clouds and their surrounding atomic and ionized envelopes. Here we present the initial results from Galactic Observations of Terahertz C+ (GOTC+), a Herschel Key Project devoted to study the [CII] fine structure emission in the galactic plane using the HIFI instrument. We use the [CII] emission together with observations of CO as a probe to understand the effects of newly--formed stars on their interstellar environment and characterize the physical and chemical state of the star-forming gas. We collected data along 16 lines--of--sight passing near star forming regions in the inner Galaxy near longitudes 330 degrees and 20 degrees. We identify fifty-eight [CII] components that are associated with high--column density molecular clouds as traced by 13CO emission. We combine [CII], 12CO, and 13CO observations to derive the physical conditions of the [CII]--emitting regions in our sample of high--column density clouds based on comparison with results from a grid of Photon Dominated Region (PDR) models. From this unbiased sample, our results suggest that most of [CII] emission originates from clouds with H2 volume densities between 10e3.5 and 10e5.5 cm^-3 and weak FUV strength (CHI_0=1-10). We find two regions where our analysis suggests high densities >10e5 cm^-3 and strong FUV fields (CHI=10e4-10e6), likely associated with massive star formation. We suggest that [CII] emission in conjunction with CO isotopes is a good tool to differentiate between regions of massive star formation (high densities/strong FUV fields) and regions that are distant from massive stars (lower densities/weaker FUV fields) along the line--of--sight △ Less

Submitted 28 July, 2010; originally announced July 2010.

Comments: To be published in A&A HIFI Special Edition

Herschel observations of EXtra-Ordinary Sources: The Terahertz spectrum of Orion KL seen at high spectral resolution

Authors: N. R. Crockett, E. A. Bergin, S. Wang, D. C. Lis, T. A. Bell, G. A. Blake, A. Boogert, B. Bumble, S. Cabrit, E. Caux, C. Ceccarelli, J. Cernicharo, C. Comito, F. Daniel, M. -L. Dubernet, M. Emprechtinger, P. Encrenaz, E. Falgarone, M. Gerin, T. F. Giesen, J. R. Goicoechea, P. F. Goldsmith, H. Gupta, R. Gusten, P. Hartogh , et al. (38 additional authors not shown)

Abstract: We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but… ▽ More We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands. △ Less

Submitted 26 July, 2010; originally announced July 2010.

Comments: Accepted for publication in the Herschel HIFI special issue of Astronomy and Astrophysics Letters, 5 pages, 3 figures

[CII] observations of H$_2$ molecular layers in transition clouds

Authors: T. Velusamy, W. D. Langer, J. L. Pineda, P. F. Goldsmith, D. Li., H. W. Yorke

Abstract: We present the first results on the diffuse transition clouds observed in [CII] line emission at 158 microns (1.9 THz) towards Galactic longitudes near 340deg (5 LOSs) and 20deg (11 LOSs) as part of the GOT C+ survey. Out of the total 146 [CII] velocity components detected by profile fitting we identify 53 as diffuse molecular clouds with associated $^{12}$CO emission but without $^{13}$CO emissio… ▽ More We present the first results on the diffuse transition clouds observed in [CII] line emission at 158 microns (1.9 THz) towards Galactic longitudes near 340deg (5 LOSs) and 20deg (11 LOSs) as part of the GOT C+ survey. Out of the total 146 [CII] velocity components detected by profile fitting we identify 53 as diffuse molecular clouds with associated $^{12}$CO emission but without $^{13}$CO emission and characterized by A$_V$ < 5 mag. We estimate the fraction of the [CII] emission in the diffuse HI layer in each cloud and then determine the [CII] emitted from the molecular layers in the cloud. We show that the excess [CII] intensities detected in a few clouds is indicative of a thick H$_2$ layer around the CO core. The wide range of clouds in our sample with thin to thick H$_2$ layers suggests that these are at various evolutionary states characterized by the formation of H$_2$ and CO layers from HI and C$^+$, respectively. In about 30% of the clouds the H$_2$ column densities (''dark gas'') traced by the [CII] is 50% or more than that traced by $^{12}$CO emission. On the average about 25% of the total H$_2$ in these clouds is in an H$_2$ layer which is not traced by CO. We use the HI, [CII], and $^{12}$CO intensities in each cloud along with simple chemical models to obtain constraints on the FUV fields and cosmic ray ionization rates. △ Less

Submitted 20 July, 2010; originally announced July 2010.

Comments: To be published in HIFI A&A special edition

C$^+$ detection of warm dark gas in diffuse clouds

Authors: W. D. Langer, T. Velusamy, J. L. Pineda, P. F. Goldsmith, D. Li, H. W. Yorke

Abstract: We present the first results of the Herschel open time key program, Galactic Observations of Terahertz C$^+$ (GOT C+) survey of the [CII] fine-structure line at 1.9 THz (158 microns) using the HIFI instrument on Herschel. We detected 146 interstellar clouds along sixteen lines-of-sight towards the inner Galaxy. We also acquired HI and CO isotopologue data along each line-of-sight for analysis of t… ▽ More We present the first results of the Herschel open time key program, Galactic Observations of Terahertz C$^+$ (GOT C+) survey of the [CII] fine-structure line at 1.9 THz (158 microns) using the HIFI instrument on Herschel. We detected 146 interstellar clouds along sixteen lines-of-sight towards the inner Galaxy. We also acquired HI and CO isotopologue data along each line-of-sight for analysis of the physical conditions in these clouds. Here we analyze 29 diffuse clouds (A$_{V}$ < 1.3 mag.) in this sample characterized by having [CII] and HI emission, but no detectable CO. We find that [CII] emission is generally stronger than expected for diffuse atomic clouds, and in a number of sources is much stronger than anticipated based on their HI column density. We show that excess [CII] emission in these clouds is best explained by the presence of a significant diffuse warm H$_2$, dark gas, component. This first [CII] 158 micron detection of warm dark gas demonstrates the value of this tracer for mapping this gas throughout the Milky Way and in galaxies. △ Less

Submitted 18 July, 2010; originally announced July 2010.

Comments: To be published in A&A HIFI Special Edition

Herschel observations of EXtra-Ordinary Sources: the present and future of spectral surveys with Herschel/HIFI

Authors: E. A. Bergin, T. G. Phillips, C. Comito, N. R. Crockett, D. C. Lis, P. Schilke, S. Wang, T. A. Bell, G. A. Blake, B. Bumble, E. Caux, S. Cabrit, C. Ceccarelli, J. Cernicharo, F. Daniel, Th. de Graauw, M. -L. Dubernet, M. Emprechtinger, P. Encrenaz, E. Falgarone, M. Gerin, T. F. Giesen, J. R. Goicoechea, P. F. Goldsmith, H. Gupta , et al. (34 additional authors not shown)

Abstract: We present initial results from the Herschel GT key program: Herschel observations of EXtra-Ordinary Sources (HEXOS) and outline the promise and potential of spectral surveys with Herschel/HIFI. The HIFI instrument offers unprecedented sensitivity, as well as continuous spectral coverage across the gaps imposed by the atmosphere, opening up a largely unexplored wavelength regime to high-resolution… ▽ More We present initial results from the Herschel GT key program: Herschel observations of EXtra-Ordinary Sources (HEXOS) and outline the promise and potential of spectral surveys with Herschel/HIFI. The HIFI instrument offers unprecedented sensitivity, as well as continuous spectral coverage across the gaps imposed by the atmosphere, opening up a largely unexplored wavelength regime to high-resolution spectroscopy. We show the spectrum of Orion KL between 480 and 560 GHz and from 1.06 to 1.115 THz. From these data, we confirm that HIFI separately measures the dust continuum and spectrally resolves emission lines in Orion KL. Based on this capability we demonstrate that the line contribution to the broad-band continuum in this molecule-rich source is ~20-40% below 1 THz and declines to a few percent at higher frequencies. We also tentatively identify multiple transitions of HD18O in the spectra. The first detection of this rare isotopologue in the interstellar medium suggests that HDO emission is optically thick in the Orion hot core with HDO/H2O ~ 0.02. We discuss the implications of this detection for the water D/H ratio in hot cores. △ Less

Submitted 13 July, 2010; originally announced July 2010.

Comments: 6 pages, 5 figures, accepted and to be published in Herschel HIFI special issue of A&A Letters

Herschel observations of extra-ordinary sources: Detection of Hydrogen Fluoride in absorption towards Orion~KL

Authors: T. G. Phillips, E. A. Bergin, D. C. Lis, D. A. Neufeld, T. A. Bell, S. Wang, N. R. Crockett, M. Emprechtinger, G. A. Blake, E. Caux, C. Ceccarelli, J. Cernicharo, C. Comito, F. Daniel, M. -L. Dubernet, P. Encrenaz, M. Gerin, T. F. Giesen, J. R. Goicoechea, P. F. Goldsmith, E. Herbst, C. Joblin, D. Johnstone, W. D. Langer, W. D. Latter , et al. (30 additional authors not shown)

Abstract: We report a detection of the fundamental rotational transition of hydrogen fluoride in absorption towards Orion KL using Herschel/HIFI. After the removal of contaminating features associated with common molecules ("weeds"), the HF spectrum shows a P-Cygni profile, with weak redshifted emission and strong blue-shifted absorption, associated with the low-velocity molecular outflow. We derive an esti… ▽ More We report a detection of the fundamental rotational transition of hydrogen fluoride in absorption towards Orion KL using Herschel/HIFI. After the removal of contaminating features associated with common molecules ("weeds"), the HF spectrum shows a P-Cygni profile, with weak redshifted emission and strong blue-shifted absorption, associated with the low-velocity molecular outflow. We derive an estimate of 2.9 x 10^13 cm^-2 for the HF column density responsible for the broad absorption component. Using our best estimate of the H2 column density within the low-velocity molecular outflow, we obtain a lower limit of ~1.6 x 10^-10 for the HF abundance relative to hydrogen nuclei, corresponding to 0.6% of the solar abundance of fluorine. This value is close to that inferred from previous ISO observations of HF J=2--1 absorption towards Sgr B2, but is in sharp contrast to the lower limit of 6 x 10^-9 derived by Neufeld et al. (2010) for cold, foreground clouds on the line of sight towards G10.6-0.4. △ Less

Submitted 13 July, 2010; originally announced July 2010.

Comments: 5 pages, 3 figures, paper to be published in the Herschel special issue of A&A letters

Herschel observations of extra-ordinary sources: Detecting spiral arm clouds by CH absorption lines

Authors: S. -L. Qin, P. Schilke, C. Comito, T. Möller, R. Rolffs, H. S. P. Müller, A. Belloche, K. M. Menten, D. C. Lis, T. G. Phillips, E. A. Bergin, T. A. Bell, N. R. Crockett, G. A. Blake, S. Cabrit, E. Caux, C. Ceccarelli, J. Cernicharo, F. Daniel, M. -L. Dubernet, M. Emprechtinger, P. Encrenaz, E. Falgarone, M. Gerin, T. F. Giesen , et al. (45 additional authors not shown)

Abstract: We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds sho… ▽ More We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds show `spiral arm' cloud characteristics, while the absorption component with the broadest line width and highest column density corresponds to the gas from the Sgr~B2 envelope. The observations show that each `spiral arm' harbors multiple velocity components, indicating that the clouds are not uniform and that they have internal structure. This line-of-sight through almost the entire Galaxy offers unique possibilities to study the basic chemistry of simple molecules in diffuse clouds, as a variety of different cloud classes are sampled simultaneously. We find that the linear relationship between CH and H$_2$ column densities found at lower $A_V$ by UV observations does not continue into the range of higher visual extinction. There, the curve flattens, which probably means that CH is depleted in the denser cores of these clouds. △ Less

Submitted 12 July, 2010; originally announced July 2010.

Comments: Accepted for publication in A&A, HIFI Special Issue

Herschel/HIFI measurements of the ortho/para ratio in water towards Sagittarius B2(M) and W31C

Authors: D. C. Lis, T. G. Phillips, P. F. Goldsmith, D. A. Neufeld, E. Herbst, C. Comito, P. Schilke, H. S. P. Müller, E. A. Bergin, M. Gerin, T. A. Bell, M. Emprechtinger, J. H. Black, G. A. Blake, F. Boulanger, S. Cabrit, E. Caux, C. Ceccarelli, J. Cernicharo, A. Coutens, N. R. Crockett, F. Daniel, E. Dartois, M. De Luca, M. -L. Dubernet , et al. (52 additional authors not shown)

Abstract: We present Herschel/HIFI observations of the fundamental rotational transitions of ortho- and para-H$_2^{16}$O and H$_2^{18}$O in absorption towards Sagittarius~B2(M) and W31C. The ortho/para ratio in water in the foreground clouds on the line of sight towards these bright continuum sources is generally consistent with the statistical high-temperature ratio of 3, within the observational uncertain… ▽ More We present Herschel/HIFI observations of the fundamental rotational transitions of ortho- and para-H$_2^{16}$O and H$_2^{18}$O in absorption towards Sagittarius~B2(M) and W31C. The ortho/para ratio in water in the foreground clouds on the line of sight towards these bright continuum sources is generally consistent with the statistical high-temperature ratio of 3, within the observational uncertainties. However, somewhat unexpectedly, we derive a low ortho/para ratio of $2.35 \pm 0.35$, corresponding to a spin temperature of $\sim$27~K, towards Sagittarius~B2(M) at velocities of the expanding molecular ring. Water molecules in this region appear to have formed with, or relaxed to, an ortho/para ratio close to the value corresponding to the local temperature of the gas and dust. △ Less

Submitted 8 July, 2010; originally announced July 2010.

Comments: Accepted for publication in A&A, HIFI Special Issue; 5 pages, 3 figures

Herschel/HIFI discovery of interstellar chloronium (H$_2$Cl$^+$)

Authors: D. C. Lis, J. C. Pearson, D. A. Neufeld, P. Schilke, H. S. P. Müller, H. Gupta, T. A. Bell, C. Comito, T. G. Phillips, E. A. Bergin, C. Ceccarelli, P. F. Goldsmith, G. A. Blake, A. Bacmann, A. Baudry, M. Benedettini, A. Benz, J. Black, A. Boogert, S. Bottinelli, S. Cabrit, P. Caselli, A. Castets, E. Caux, J. Cernicharo , et al. (80 additional authors not shown)

Abstract: We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S)… ▽ More We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models. △ Less

Submitted 8 July, 2010; originally announced July 2010.

Comments: Accepted for publication in A&A, HIFI Special Issue; 4 pages, 2 figures, 1 table

Journal ref: Astron. Astrophys. 521 (2010) Art. No. L9

Herschel observations of ortho- and para-oxidaniumyl (H2O+) in spiral arm clouds toward Sgr B2(M)

Authors: P. Schilke, C. Comito, H. S. P. Mueller, E. A. Bergin, E. Herbst, D. C. Lis, D. A. Neufeld, T. G. Phillips, T. A. Bell, G. A. Blake, S. Cabrit, E. Caux, C. Ceccarelli, J. Cernicharo, N. R. Crockett, F. Daniel, M. -L. Dubernet, M. Emprechtinger, P. Encrenaz, M. ~Gerin, T. F. Giesen, J. R. Goicoechea, P. F. Goldsmith, H. Gupta, C. Joblin , et al. (33 additional authors not shown)

Abstract: H2O+ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H2O+ in most velocity intervals is 4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the for… ▽ More H2O+ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H2O+ in most velocity intervals is 4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the formation temperature and current physical temperature of the gas hosting H2O+ is discussed, but no firm conclusion is reached. In the velocity interval 0 to 60 km/s, an ortho/para ratio of below unity is found, but if this is due to an artifact of contamination by other species or real is not clear. △ Less

Submitted 5 July, 2010; originally announced July 2010.

Comments: A&A, HIFI special issue, accepted

Journal ref: Astron. Astrophys. 521 (2010) Art. No. L11

Do Lognormal Column-Density Distributions in Molecular Clouds Imply Supersonic Turbulence?

Authors: K. Tassis, D. A. Christie, A. Urban, J. L. Pineda, T. Ch. Mouschovias, H. W. Yorke, H. Martel

Abstract: Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supers… ▽ More Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supersonic turbulence. The first star-forming cloud is simulated using smoothed particle hydrodynamics (SPH); in this case gravity, opposed only by thermal-pressure forces, drives the evolution. The second cloud is magnetically subcritical with subsonic turbulence, simulated using nonideal MHD; in this case the evolution is due to gravitationally-driven ambipolar diffusion. The third cloud is isothermal, self-gravitating, and has a smooth density distribution analytically approximated with a uniform inner region and an r^-2 profile at larger radii. We show that in all three cases the column-density distributions are lognormal. Power-law tails develop only at late times (or, in the case of the smooth analytic profile, for strongly centrally concentrated configurations), when gravity dominates all opposing forces. It therefore follows that lognormal column-density distributions are generic features of diverse model clouds, and should not be interpreted as being a consequence of supersonic turbulence. △ Less

Submitted 14 June, 2010; originally announced June 2010.

Comments: 6 pages, 6 figures, accepted for publication in MNRAS

Clouds, filaments and protostars: the Herschel Hi-GAL Milky Way

Authors: S. Molinari, B. Swinyard, J. Bally, M. Barlow, J. P. Bernard, P. Martin, T. Moore, A. Noriega-Crespo, R. Plume, L. Testi, A. Zavagno, A. Abergel, B. Ali, L. Anderson, P. André, J. P. Baluteau, C. Battersby, M. T. Beltrán, M. Benedettini, N. Billot, J. Blommaert, S. Bontemps, F. Boulanger, J. Brand, C. Brunt , et al. (99 additional authors not shown)

Abstract: We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2° x 2° tiles approximately centered at l=30° and l=59°. The two regions are extremely rich in intense and highly structure… ▽ More We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2° x 2° tiles approximately centered at l=30° and l=59°. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59° field; a A_V value between 5 and 10 is found for the l=30° field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows. △ Less

Submitted 18 May, 2010; originally announced May 2010.

Comments: A&A, accepted

Evolution of Massive Protostars via Disk Accretion

Authors: Takashi Hosokawa, Harold W. Yorke, Kazuyuki Omukai

Abstract: Mass accretion onto (proto-)stars at high accretion rates > 10^-4 M_sun/yr is expected in massive star formation. We study the evolution of massive protostars at such high rates by numerically solving the stellar structure equations. In this paper we examine the evolution via disk accretion. We consider a limiting case of "cold" disk accretion, whereby most of the stellar photosphere can radiate f… ▽ More Mass accretion onto (proto-)stars at high accretion rates > 10^-4 M_sun/yr is expected in massive star formation. We study the evolution of massive protostars at such high rates by numerically solving the stellar structure equations. In this paper we examine the evolution via disk accretion. We consider a limiting case of "cold" disk accretion, whereby most of the stellar photosphere can radiate freely with negligible backwarming from the accretion flow, and the accreting material settles onto the star with the same specific entropy as the photosphere. We compare our results to the calculated evolution via spherically symmetric accretion, the opposite limit, whereby the material accreting onto the star contains the entropy produced in the accretion shock front. We examine how different accretion geometries affect the evolution of massive protostars. For cold disk accretion at 10^-3 M_sun/yr the radius of a protostar is initially small, about a few R_sun. After several solar masses have accreted, the protostar begins to bloat up and for M \simeq 10 M_sun the stellar radius attains its maximum of 30 - 400 R_sun. The large radius about 100 R_sun is also a feature of spherically symmetric accretion at the same accreted mass and accretion rate. Hence, expansion to a large radius is a robust feature of accreting massive protostars. At later times the protostar eventually begins to contract and reaches the Zero-Age Main-Sequence (ZAMS) for M \simeq 30 M_sun, independent of the accretion geometry. For accretion rates exceeding several 10^-3 M_sun/yr the protostar never contracts to the ZAMS. The very large radius of several 100s R_sun results in a low effective temperature and low UV luminosity of the protostar. Such bloated protostars could well explain the existence of bright high-mass protostellar objects, which lack detectable HII regions. △ Less

Submitted 17 May, 2010; originally announced May 2010.

Comments: 20 pages, 16 figures

Detection of interstellar oxidaniumyl: abundant H2O+ towards the star-forming regions DR21, Sgr B2, and NGC6334

Authors: V. Ossenkopf, H. S. P. Müller, D. C. Lis, P. Schilke, T. A. Bell, S. Bruderer, E. Bergin, C. Ceccarelli, C. Comito, J. Stutzki, A. Bacman, A. Baudry, A. O. Benz, M. Benedettini, O. Berne, G. Blake, A. Boogert, S. Bottinelli, F. Boulanger, S. Cabrit, P. Caselli, E. Caux, J. Cernicharo, C. Codella, A. Coutens , et al. (77 additional authors not shown)

Abstract: We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum bac… ▽ More We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158μ\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2. △ Less

Submitted 14 May, 2010; originally announced May 2010.

Comments: Accepted for publication in A&A

Journal ref: Astron. Astrophys. 518 (2010) Art. No. L111

Toward Understanding Massive Star Formation

Authors: Hans Zinnecker, Harold W. Yorke

Abstract: Although fundamental for astrophysics, the processes that produce massive stars are not well understood. Large distances, high extinction, and short timescales of critical evolutionary phases make observations of these processes challenging. Lacking good observational guidance, theoretical models have remained controversial. This review offers a basic description of the collapse of a massive mol… ▽ More Although fundamental for astrophysics, the processes that produce massive stars are not well understood. Large distances, high extinction, and short timescales of critical evolutionary phases make observations of these processes challenging. Lacking good observational guidance, theoretical models have remained controversial. This review offers a basic description of the collapse of a massive molecular core and a critical discussion of the three competing concepts of massive star formation: - monolithic collapse in isolated cores - competitive accretion in a protocluster environment - stellar collisions and mergers in very dense systems We also review the observed outflows, multiplicity, and clustering properties of massive stars, the upper initial mass function and the upper mass limit. We conclude that high-mass star formation is not merely a scaled-up version of low-mass star formation with higher accretion rates, but partly a mechanism of its own, primarily owing to the role of stellar mass and radiation pressure in controlling the dynamics. △ Less

Submitted 9 July, 2007; originally announced July 2007.

Comments: 139 pages, 18 figures, 5 tables, glossary

Journal ref: Ann.Rev.Astron.Astrophys.45:481-563,2007

ExoPTF Science Uniquely Enabled by Far-IR Interferometry: Probing the Formation of Planetary Systems, and Finding and Characterizing Exoplanets

Authors: David Leisawitz, Tom Armstrong, Chad Bender, Dominic Benford, Daniella Calzetti, John Carpenter, William C. Danchi, Michel Fich, Dale Fixsen, Daniel Y. Gezari, Matt Griffin, Martin Harwit, Alan J. Kogut, William D. Langer, Charles Lawrence, Dan Lester, Lee G. Mundy, Joan Najita, David Neufeld, Goran Pilbratt, Stephen Rinehart, Aki Roberge, Eugene Serabyn, Sachindev Shenoy, Hiroshi Shibai , et al. (6 additional authors not shown)

Abstract: By providing sensitive sub-arcsecond images and integral field spectroscopy in the 25 - 400 micron wavelength range, a far-IR interferometer will revolutionize our understanding of planetary system formation, reveal otherwise-undetectable planets through the disk perturbations they induce, and spectroscopically probe the atmospheres of extrasolar giant planets in orbits typical of most of the pl… ▽ More By providing sensitive sub-arcsecond images and integral field spectroscopy in the 25 - 400 micron wavelength range, a far-IR interferometer will revolutionize our understanding of planetary system formation, reveal otherwise-undetectable planets through the disk perturbations they induce, and spectroscopically probe the atmospheres of extrasolar giant planets in orbits typical of most of the planets in our solar system. The technical challenges associated with interferometry in the far-IR are greatly relaxed relative to those encountered at shorter wavelengths or when starlight nulling is required. A structurally connected far-IR interferometer with a maximum baseline length of 36 m can resolve the interesting spatial structures in nascent and developed exoplanetary systems and measure exozodiacal emission at a sensitivity level critical to TPF-I mission planning. The Space Infrared Interferometric Telescope was recommended in the Community Plan for Far-IR/Submillimeter Space Astronomy, studied as a Probe-class mission, and estimated to cost 800M dollars. The scientific communities in Europe, Japan, and Canada have also demonstrated a keen interest in far-IR interferometry through mission planning workshops and technology research, suggesting the possibility of an international collaborative effort. △ Less

Submitted 5 July, 2007; originally announced July 2007.

Comments: 7 pages, 4 figures, submitted to the Exoplanet Task Force (AAAC), 30 March 2007

Photon Bubbles in the Circumstellar Envelopes of Young Massive Stars

Authors: N. J. Turner, E. Quataert, H. W. Yorke

Abstract: We show that the optically-thick dusty envelopes surrounding young high-mass stars are subject to the photon bubble instability. The infrared radiation passing through the envelope amplifies magnetosonic disturbances, with growth rates in our local numerical radiation MHD calculations that are consistent with a linear analysis. Modes with wavelengths comparable to the gas pressure scale height g… ▽ More We show that the optically-thick dusty envelopes surrounding young high-mass stars are subject to the photon bubble instability. The infrared radiation passing through the envelope amplifies magnetosonic disturbances, with growth rates in our local numerical radiation MHD calculations that are consistent with a linear analysis. Modes with wavelengths comparable to the gas pressure scale height grow by more than two orders of magnitude in a thousand years, reaching non-linear amplitudes within the envelope lifetime. If the magnetic pressure in the envelope exceeds the gas pressure, the instability develops into trains of propagating shocks. Radiation escapes readily through the low-density material between the shocks, enabling accretion to continue despite the Eddington limit imposed by the dust opacity. The supersonic motions arising from the photon bubble instability can help explain the large velocity dispersions of hot molecular cores, while conditions in the shocked gas are suitable for maser emission. We conclude that the photon bubble instability may play a key role in the formation of massive stars. △ Less

Submitted 28 January, 2007; originally announced January 2007.

Comments: Accepted by the Astrophysical Journal; 18 pages, 4 figures

Journal ref: Astrophys.J.662:1052-1058,2007

Massive stars and the energy balance of the interstellar medium. II. The 35 solar mass star and a solution to the "missing wind problem"

Authors: Tim Freyer, Gerhard Hensler, Harold W. Yorke

Abstract: We continue our numerical analysis of the morphological and energetic influence of massive stars on their ambient interstellar medium for a 35 solar mass star that evolves from the main sequence through red supergiant and Wolf-Rayet phases, until it ultimately explodes as a supernova. We find that structure formation in the circumstellar gas during the early main-sequence evolution occurs as in… ▽ More We continue our numerical analysis of the morphological and energetic influence of massive stars on their ambient interstellar medium for a 35 solar mass star that evolves from the main sequence through red supergiant and Wolf-Rayet phases, until it ultimately explodes as a supernova. We find that structure formation in the circumstellar gas during the early main-sequence evolution occurs as in the 60 solar mass case but is much less pronounced because of the lower mechanical wind luminosity of the star. Since on the other hand the shell-like structure of the HII region is largely preserved, effects that rely on this symmetry become more important. At the end of the stellar lifetime 1% of the energy released as Lyman continuum radiation and stellar wind has been transferred to the circumstellar gas. From this fraction 10% is kinetic energy of bulk motion, 36% is thermal energy, and the remaining 54% is ionization energy of hydrogen. The sweeping up of the slow red supergiant wind by the fast Wolf-Rayet wind produces remarkable morphological structures and emission signatures, which are compared with existing observations of the Wolf-Rayet bubble S308. Our model reproduces the correct order of magnitude of observed X-ray luminosity, the temperature of the emitting plasma as well as the limb brightening of the intensity profile. This is remarkable, because current analytical and numerical models of Wolf-Rayet bubbles fail to consistently explain these features. A key result is that almost the entire X-ray emission in this stage comes from the shell of red supergiant wind swept up by the shocked Wolf-Rayet wind rather than from the shocked Wolf-Rayet wind itself as hitherto assumed and modeled. This offers a possible solution to what is called the ``missing wind problem'' of Wolf-Rayet bubbles. △ Less

Submitted 5 December, 2005; originally announced December 2005.

Comments: 52 pages, 20 figures, 2 tables, accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.638:262-280,2006

Turbulent Mixing in the Outer Solar Nebula

Authors: N. J. Turner, K. Willacy, G. Bryden, H. W. Yorke

Abstract: The effects of turbulence on the mixing of gases and dust in the outer Solar nebula are examined using 3-D MHD calculations in the shearing-box approximation with vertical stratification. The turbulence is driven by the magneto-rotational instability. The magnetic and hydrodynamic stresses in the turbulence correspond to an accretion time at the midplane about equal to the lifetimes of T Tauri d… ▽ More The effects of turbulence on the mixing of gases and dust in the outer Solar nebula are examined using 3-D MHD calculations in the shearing-box approximation with vertical stratification. The turbulence is driven by the magneto-rotational instability. The magnetic and hydrodynamic stresses in the turbulence correspond to an accretion time at the midplane about equal to the lifetimes of T Tauri disks, while accretion in the surface layers is thirty times faster. The mixing resulting from the turbulence is also fastest in the surface layers. The mixing rate is similar to the rate of radial exchange of orbital angular momentum, so that the Schmidt number is near unity. The vertical spreading of a trace species is well-matched by solutions of a damped wave equation when the flow is horizontally-averaged. The damped wave description can be used to inexpensively treat mixing in 1-D chemical models. However, even in calculations reaching a statistical steady state, the concentration at any given time varies substantially over horizontal planes, due to fluctuations in the rate and direction of the transport. In addition to mixing species that are formed under widely varying conditions, the turbulence intermittently forces the nebula away from local chemical equilibrium. The different transport rates in the surface layers and interior may affect estimates of the grain evolution and molecular abundances during the formation of the Solar system. △ Less

Submitted 10 November, 2005; originally announced November 2005.

Comments: To appear in the Astrophysical Journal; 20 pages, 9 figures

Journal ref: Astrophys.J.639:1218-1226,2006

Photon Bubbles in Young Massive Stars

Authors: N. J. Turner, H. W. Yorke, A. Socrates, O. M. Blaes

Abstract: Spectroscopic studies indicate that gas in the photospheres of young O stars moves at speeds up to the sound speed. We show, using two-dimensional radiation MHD calculations and results from a local linear analysis, that the motions may be due to photon bubble instability if young O stars have magnetic fields. Spectroscopic studies indicate that gas in the photospheres of young O stars moves at speeds up to the sound speed. We show, using two-dimensional radiation MHD calculations and results from a local linear analysis, that the motions may be due to photon bubble instability if young O stars have magnetic fields. △ Less

Submitted 31 March, 2004; originally announced April 2004.

Comments: 4 pages, 3 figures; to appear in "Gravitational Collapse: from Massive Stars to Planets", eds G. Garcia-Segura, G. Tenorio-Tagle, H. Yorke and J. Franco, Rev. Mex. A. A. S. C

Journal ref: Rev.Mex.Astron.Astrofis. 22 (2004) 54-57

Does the Lunar Surface Still Offer Value As a Site for Astronomical Observatories?

Authors: Daniel F. Lester, Harold W. Yorke, John C. Mather

Abstract: Current thinking about the Moon as a destination has revitalized interest in lunar astronomical observatories. Once seen by a large scientific community as a highly enabling site, the dramatic improvement in capabilities for free-space observatories prompts reevaluation of this interest. Whereas the lunar surface offers huge performance advantages for astronomy over terrestrial sites, free-space… ▽ More Current thinking about the Moon as a destination has revitalized interest in lunar astronomical observatories. Once seen by a large scientific community as a highly enabling site, the dramatic improvement in capabilities for free-space observatories prompts reevaluation of this interest. Whereas the lunar surface offers huge performance advantages for astronomy over terrestrial sites, free-space locales such as Earth orbit or Lagrange points offer performance that is superior to what could be achieved on the Moon. While astronomy from the Moon may be cost effective once infrastructure is there, it is in many respects no longer clearly enabling compared to free-space. △ Less

Submitted 5 February, 2004; v1 submitted 14 January, 2004; originally announced January 2004.

Comments: 14 pages, 1 figure; in press Space Policy

Formation of Massive Stars via Accretion

Authors: Harold W. Yorke

Abstract: The collapse of massive molecular clumps can produce high mass stars, but the evolution is not simply a scaled-up version of low mass star formation. Outflows and radiative effects strongly hinder the formation of massive stars via accretion. A necessary condition for accretion growth of a hydrostatic object up to high masses M > 20 M_sun (rather than coalescence of optically thick objects) is t… ▽ More The collapse of massive molecular clumps can produce high mass stars, but the evolution is not simply a scaled-up version of low mass star formation. Outflows and radiative effects strongly hinder the formation of massive stars via accretion. A necessary condition for accretion growth of a hydrostatic object up to high masses M > 20 M_sun (rather than coalescence of optically thick objects) is the formation of and accretion through a circumstellar disk. Once the central object has accreted approximately 10 M_sun it has already evolved to core hydrogen-burning; the resultant main sequence star continues to accrete material as it begins to photoevaporate its circumstellar disk (and any nearby disks) on a timescale of 100,000 years, similar to the accretion timescale. Until the disk(s) is (are) completely photoevaporated, this configuration is observable as an ultra-compact HII region (UCHII). The final mass of the central star (and any nearby neighboring systems) is determined by the interplay between radiation acceleration, UV photoevaporation, stellar winds and outflows, and the accretion through the disk. Several aspects of this evolutionary sequence have been simulated numerically, resulting in a "proof of concept". This scenario places strong constraints on the accretion rate necessary to produce high mass stars and offers an opportunity to test the accretion hypothesis. △ Less

Submitted 22 December, 2003; originally announced December 2003.

Comments: 12 pages, 11 figures, IAU Symposium 221, invited major review