Grain processes in massive star formation

Observational evidence suggests that stars greater than 100 M(solar) exist in the Galaxy and Large Magellanic Cloud (LMC), however classical star formation theory predicts stellar mass limits of only approx. 60 M(solar). A protostellar accretion flow consists of inflowing gas and dust. Grains are destroyed as they are near the central protostar creating a dust shell or cocoon. Radiation pressure acting on the grain can halt the inflow of material thereby limiting the amount of mass accumulated by the protostar. We first consider rather general constraints on the initial grain to gas ratio and mass accretion rates that permit inflow. We further constrain these results by constructing a numerical model. Radiative deceleration of grains and grain destruction processes are explicitly accounted for in an iterative solution of the radiation-hydrodynamic equations. Findings seem to suggest that star formation by spherical accretion requires rather extreme preconditioning of the grain and gas environment

Physical conditions in photodissociation regions: Application to galactic nuclei

Infrared and sub-millimeter observations are used in a simple procedure to determine average physical properties of the neutral interstellar medium in Galactic photodissociation regions as well as in ensembles of clouds which exist in the nuclei of luminous infrared galaxies. The relevant observations include the Infrared Astronomy Satellite (IRAS) infrared continuum measurements, infrared spectroscopy of the fine-structure lines of SiII 35 microns, OI 63 microns, and CII 158 microns, and the 2.6 mm CO (J=1-0) rotational transition. The diagnostic capabilities of the OI 145 microns line is also addressed. Researchers attribute these emission lines as well as the continuum to the atomic/molecular photodissociation region on the surfaces of molecular clouds which are illuminated by strong ultraviolet fields. They use the theoretical photodissociation region models of Tielens and Hollenbach (1985, Ap. J., 291, 722) to construct simple diagrams which utilize line ratios and line to continuum ratios to determine the average gas density n, the average incident far-ultraviolet flux G sub o, and the temperature of the atomic gas T

The Chemistry of Interstellar OH+, H2O+, and H3O+: Inferring the Cosmic Ray Ionization Rates from Observations of Molecular Ions

We model the production of OH+, H2O+, and H3O+ in interstellar clouds, using a steady state photodissociation region code that treats the freeze-out of gas species, grain surface chemistry, and desorption of ices from grains. The code includes PAHs, which have important effects on the chemistry. All three ions generally have two peaks in abundance as a function of depth into the cloud, one at A_V<~1 and one at A_V~3-8, the exact values depending on the ratio of incident ultraviolet flux to gas density. For relatively low values of the incident far ultraviolet flux on the cloud ({\chi}<~ 1000; {\chi}= 1= local interstellar value), the columns of OH+ and H2O+ scale roughly as the cosmic ray primary ionization rate {\zeta}(crp) divided by the hydrogen nucleus density n. The H3O+ column is dominated by the second peak, and we show that if PAHs are present, N(H3O+) ~ 4x10^{13} cm^{-2} independent of {\zeta}(crp) or n. If there are no PAHs or very small grains at the second peak, N(H3O+) can attain such columns only if low ionization potential metals are heavily depleted. We also model diffuse and translucent clouds in the interstellar medium, and show how observations of N(OH+)/N(H) and N(OH+)/N(H2O+) can be used to estimate {\zeta}(crp)/n, {\chi}/n and A_V in them. We compare our models to Herschel observations of these two ions, and estimate {\zeta}(crp) ~ 4-6 x 10^-16 (n/100 cm^-3) s^-1 and \chi/n = 0.03 cm^3 for diffuse foreground clouds towards W49N

Herschel PACS and SPIRE spectroscopy of the Photodissociation Regions associated with S 106 and IRAS 23133+6050

Photodissociation regions (PDRs) contain a large fraction of all of the interstellar matter in galaxies. Classical examples include the boundaries between ionized regions and molecular clouds in regions of massive star formation, marking the point where all of the photons energetic enough to ionize hydrogen have been absorbed. In this paper we determine the physical properties of the PDRs associated with the star forming regions IRAS 23133+6050 and S 106 and present them in the context of other Galactic PDRs associated with massive star forming regions. We employ Herschel PACS and SPIRE spectroscopic observations to construct a full 55-650 {\mu}m spectrum of each object from which we measure the PDR cooling lines, other fine- structure lines, CO lines and the total far-infrared flux. These measurements are then compared to standard PDR models. Subsequently detailed numerical PDR models are compared to these predictions, yielding additional insights into the dominant thermal processes in the PDRs and their structures. We find that the PDRs of each object are very similar, and can be characterized by a two-phase PDR model with a very dense, highly UV irradiated phase (n $\sim$ 10^6 cm^(-3), G$_0$ $\sim$ 10^5) interspersed within a lower density, weaker radiation field phase (n $\sim$ 10^4 cm^(-3), G$_0$ $\sim$ 10^4). We employed two different numerical models to investigate the data, firstly we used RADEX models to fit the peak of the $^{12}$CO ladder, which in conjunction with the properties derived yielded a temperature of around 300 K. Subsequent numerical modeling with a full PDR model revealed that the dense phase has a filling factor of around 0.6 in both objects. The shape of the $^{12}$CO ladder was consistent with these components with heating dominated by grain photoelectric heating. An extra excitation component for the highest J lines (J > 20) is required for S 106.Comment: 20 pages, 10 figures, A&A Accepte

ISO Detection of CO⁺ toward the protostar IRAS 16293-2422

In this letter we report the detection of eight high-N rotational transitions of CO+ towards a low mass protostar, IRAS 16293-2422. The source was observed with the Long Wavelength Spectrometer on board the Infrared Space Observatory. This is the first time that CO+ has been detected in a low luminosity source and the first time that high-N lines have been detected in any source. The detection of these lines was not predicted by models and consequently, their interpretation is a challenge. We discuss the possibility that the observed CO+ emission originates in the dense inner regions illuminated by the UV field created in the accretion shock (formed by infalling material), and conclude that this is an improbable explanation. We have also considered the possibility that a strong, dissociative J-shock at ~ 500 AU from the star is the origin of the CO+ emission. This model predicts CO+ column densities in rough agreement with the observations if the magnetic field is ~ 1 mG and the shock velocity is 100 km s-1

Detection of CO⁺ with ISO towards the protostar IRAS16293-242

We observed the low luminosity (and low mass) protostar IRAS16293-2422 with the Long Wavelength Spectrometer on board the Infrared Space Observatory. The observed line spectrum is very reach and shows transitions of several molecules and atoms. Here we report the detection of eight high-N rotational transitions of CO+. This is the first time that CO+ has been detected in a low luminosity source and the first time that high-N lines have been detected in any source. The detection of these lines was not predicted by models and consequently, their interpretation is a challenge. We discuss the possibility that the observed CO+ emission originates in the dense inner regions illuminated by the UV field created in the accretion shock (formed by infalling material), and conclude that this is an improbable explanation. We have also considered the possibility that a strong, dissociative J-shock at ~500 AU from the star is the origin of the CO+ emission. This model predicts CO+ column densities in rough agreement with the observations if the magnetic field is ~1 mG and the shock velocity is 100 km s-1

ISO-LWS observations of IRAS16293-2422

We obtained LWS grating spectra toward IRAS 16293-2422 and the surrounding region, which covers the entire extent of the molecular outflow. The LWS spectra show that the region is relatively uncontaminated by PhotoDissociationRegion (PDR)-like emission, showing only a weak diffuse CII emission. The on-source spectrum revealed the presence of the OI(63μm) line and several lines from CO, H2O and OH molecules. In this work we derive the macroscopic quantities associated with the UV-illuminated emitting gas which surrounds IRAS16293-2422 and compare it with previous studies. We show that the molecular lines originate in a hot (~1600 K), dense (~ 3·104cm-3) and extended (~ 8·1016cm) region, that we interprete as the shock of the wind impacting obliquely with the walls of the cavity created by the wind itself. The OI(63μm) line observed by the Kuiper Airborne Observatory (KAO: Ceccarelli et al. 1997a) at ~ 1.2·1017cm west from the central source is hence interpreted as the head of the shock where the wind strikes the ambient gas. Finally we speculate that the OI(63μm) line emission seen on-source originates in the collapsing envelope that surrounds the central object(s

Two-dimensional AMR simulations of colliding flows

Colliding flows are a commonly used scenario for the formation of molecular clouds in numerical simulations. Due to the thermal instability of the warm neutral medium, turbulence is produced by cooling. We carry out a two-dimensional numerical study of such colliding flows in order to test whether statistical properties inferred from adaptive mesh refinement (AMR) simulations are robust with respect to the applied refinement criteria. We compare probability density functions of various quantities as well as the clump statistics and fractal dimension of the density fields in AMR simulations to a static-grid simulation. The static grid with 2048^2 cells matches the resolution of the most refined subgrids in the AMR simulations. The density statistics is reproduced fairly well by AMR. Refinement criteria based on the cooling time or the turbulence intensity appear to be superior to the standard technique of refinement by overdensity. Nevertheless, substantial differences in the flow structure become apparent. In general, it is difficult to separate numerical effects from genuine physical processes in AMR simulations.Comment: 6 pages, 6 figures, submitted to A&