No Arabic abstract
We present a radiative transfer model, which is applicable to the study of submillimetre spectral line observations of protostellar envelopes. The model uses an exact, non-LTE, spherically symmetric radiative transfer `Stenholm method, which numerically solves the radiative transfer problem by the process of `Lambda-iteration. We also present submillimetre spectral line data of the Class 0 protostars NGC1333-IRAS2 and Serpens SMM4. We examine the physical constraints which can be used to limit the number and range of parameters used in protostellar envelope models, and identify the turbulent velocity and tracer molecule abundance as the principle sources of uncertainty in the radiative transfer modelling. We explore the trends in the appearance of the predicted line profiles as key parameters in the models are varied. We find that the separation of the two peaks of a typical infall profile is dependent not on the evolutionary status of the collapsing protostar, but on the turbulent velocity dispersion in the envelope. We also find that the line shapes can be significantly altered by rotation. Fits are found for the observed line profiles of IRAS2 and SMM4 using plausible infall model parameters. The density and velocity profiles in our best fit models are inconsistent with a singular isothermal sphere model. We find better agreement with a form of collapse which assumes non-static initial conditions. We also find some evidence that the infall velocities are retarded from free-fall towards the centre of the cloud, probably by rotation, and that the envelope of SMM4 is rotationally flattened.
This paper investigates small-scale structures of dense gas and dust around the low-mass protostellar binary NGC1333-IRAS2 using millimeter-wavelength aperture-synthesis observations from the OVRO and BIMA interferometers. The detected 3 mm continuum emission is consistent with models of the envelope around IRAS2A, based on previously reported submillimeter continuum images, down to the 3, or 500 AU, resolution of the interferometer data. Our data constrain the contribution of an unresolved point source to 22 mJy. Within the accuracy of the parameters describing the envelope model, the point source flux has an uncertainty by up to 25%. We interpret this point source as a cold disk of mass gtrsim 0.3 M_odot. The same envelope model also reproduces aperture-synthesis line observations of the optically thin isotopic species C34S and H13CO+. The more optically thick main isotope lines show a variety of components in the protostellar environment: N2H+ is closely correlated with dust concentrations as seen at submillimeter wavelengths and is particularly strong toward the starless core IRAS2C. We hypothesize that N2H+ is destroyed through reactions with CO that is released from icy grains near the protostellar sources IRAS2A and B. CS, HCO+, and HCN have complex line shapes apparently affected by both outflow and infall. In addition to the east-west jet from IRAS2A, a north-south velocity gradient near this source indicates a second, perpendicular outflow. This suggests the presence of a binary companion within 0.3 (65 AU) from IRAS2A as driving source of this outflow. Alternative explanations of the velocity gradient such as rotation in a circumstellar envelope or a single, wide-angle outflow are less likely.
We report on spectro-imaging observations employing Spitzer IRS and Herschel PACS, aiming to constrain the physical conditions around SMM3 and SMM4 in Serpens. The combined power of both instruments provides an almost complete wavelength coverage between 5 and 200 micron at an angular resolution of 10. We detect line emission from all major molecular (H2, CO, H2O and OH) and many atomic ([OI], [CII], [FeII], [SiII] and [SI]) coolants. Line emission tends to peak at distances of 10 - 20 from the protostellar sources, at positions of known outflow shocks. The only exception is [CII] which likely traces a PDR excited from the neighboring source SMM6. Excitation analysis indicates that H2 and CO originate from gas at two distinct rotational temperatures of 300 K and 1000 K, while H2O and OH emission corresponds to rotational temperatures of 100 - 200 K. The morphological and physical association between CO and H2 suggests a common excitation mechanism which allows direct comparisons between the two molecules. The CO/H2 abundance ratio varies from 10^-5 in the warm gas up to 10^-4 in the hotter regions. The occurrence of J-shocks is suggested by the strong atomic/ionic (except for [CII]) emission as well as a number of line ratio diagnostics. Both C- and J-shocks can account for the observed molecular emission, however J-shocks are strongly advocated by the atomic emission and provide simpler and more homogeneous solutions for CO and H2. C-shocks describe better the emission from H2O and OH. The variations in the CO/H2 abundance ratio for gas at different temperatures can be interpreted by their reformation rates in dissociative J-type shocks, or the simultaneous influence of both C and J shocks.
The evolution of star-forming core analogues undergoing inside-out collapse is studied with a multi-point chemodynamical model which self-consistently computes the abundance distribution of chemical species in the core. For several collapse periods the output chemistry of infall tracer species such as HCO+, CS, and N2H+, is then coupled to an accelerated Lambda-iteration radiative transfer code, which predicts the emerging molecular line profiles using two different input gas/dust temperature distributions. We investigate the sensitivity of the predicted spectral line profiles and line asymmetry ratios to the core temperature distribution, the time-dependent model chemistry, as well as to ad hoc abundance distributions. The line asymmetry is found to be strongly dependent on the adopted chemical abundance distribution. In general, models with a warm central region show higher values of blue asymmetry in optically thick HCO+ and CS lines than models with a starless core temperature profile. We find that in the formal context of Shu-type inside-out infall, and in the absence of rotation or outflows, the relative blue asymmetry of certain HCO+ and CS transitions is a function of time and, subject to the foregoing caveats, can act as a collapse chronometer. The sensitivity of simulated HCO+ line profiles to linear radial variations, subsonic or supersonic, of the internal turbulence field is investigated in the separate case of static cores.
The role of accretion disks in the formation of low-mass stars has been well assessed by means of high angular resolution observations at various wavelengths. These findings confirm the prediction that conservation of angular momentum during the collapse leading to the formation of a star is bound to produce flattening and rotation of the collapsing core. What about high-mass stars? At present, several authors have reported on detections of disks around high-mass YSOs. Notwithstanding these important results, the presence of disks rotating about high-mass stars is not sufficient by itself to prove unambiguously the accretion model: what is needed is iron-clad evidence of infall. Such evidence is very difficult to find, as the free-fall velocity becomes significant only very close to the accreting star, i.e., over a region of a few 0.01 pc ($sim$2000 au), which is very difficult to access and disentangle from the surrounding quiescent or rotating material. In this chapter we discuss how to characterize the infall of material in a sample of 36 high-mass accretion disk candidates covering a broad range of luminosities, from 10$^3$ $L_odot$ to 10$^6$ $L_odot$, compiled by Beltran & de Wit (2016) with the next generation Very Large Array (ngVLA).
We aim to study dust properties of massive star forming regions in the outer Galaxy, in a direction opposite to the Galactic center. We present observations of six outer Galaxy point sources IRAS 01045+6505, 01420+6401, 05271+3059, 05345+3556, 20222+3541 and 20406+4555, taken with the Submillimeter Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT) at 450 and 850 micron. Single temperature greybody models are fitted to the Spectral Energy Distribution of the detected sub-mm cores to derive dust temperature, dust emissivity index and optical depth at 250 micron. The observed radial intensity profiles of the sub-mm cores were fitted with power laws to derive the indices describing the density distribution. At a resolution of 15 all six IRAS point sources show multiple emission peaks. Only four out of fourteen detected sub-mm cores show associated mid-infrared emission. For the sub-mm cores we derive dust temperatures of 32+-5 K and dust emissivity indices between 0.9 and 2.5. The density profiles of the sub-mm cores can be fitted by a single power law distribution with indices -1.5+-0.3, with most cores showing an index of -1.5. This is consistent with most observations of massive star forming regions and supports predictions of models of star formation which consider non-thermal support against gravitational collapse.