Molecular Line Emission from Accretion Disks Around YSOs

In this work, we model the expected molecular emission from protoplanetary disks, modifying different physical parameters, such as dust grain size, mass accretion rate, viscosity, and disk radius, to obtain observational signatures in these sources. Having in mind possible future observations, we study correlations between physical parameters and observational characteristics. Our aim is to determine the kind of observations that will allow us to extract information about the physical parameters of disks. We also present prospects for molecular line observations of protoplanetary disks, using millimeter and submillimeter interferometers (e.g., SMA or ALMA), based on our results.     

Scientific role of ACA for low-mass star-formation study

We discuss the scientific role of the Atacama Compact Array (ACA), the Japanese contribution to the ALMA project, for low-mass star-formation study. Our recent observations of several low-mass protostellar envelopes in the submillimeter CS (J=7–6) and HCN (J=4–3) lines with the SMA and ASTE have revealed that these submillimeter emissions are more extended than ∼2000 AU and show different velocity structures from those traced by millimeter lines. These results suggest the importance of taking short-spacing informations the ACA can offer. Our comprehensive imaging simulations of these protostellar envelopes, as well as prestellar cores and debris disks, unprecedentedly demonstrate the scientific importance of ACA.     

A massive, dusty toroid with large grains in the pre-planetary nebula IRAS22036+5306

Using the Submillimeter Array (SMA), we have obtained high angular-resolution (∼1″) interferometric maps of the submillimeter (0.88 mm) continuum and CO J=3–2 line from IRAS 22036+5306 (I 22036), a bipolar pre-planetary nebula (PPN) with knotty jets discovered in our HST SNAPshot survey of young PPNe. In addition, we have obtained supporting lower-resolution (∼10″) 2.6 mm continuum and CO, ¹³CO J=1–0 observations with the Owens Valley Radio Observatory (OVRO) interferometer. We find an unresolved source of submillimeter (and millimeter-wave) continuum emission in I 22036, implying a very substantial mass (0.02–0.04M _⊙) of large (i.e., radius ≳1 mm), cold (≲50 K) dust grains associated with I 22036’s toroidal waist. The CO J=3–2 observations show the presence of a very fast (∼220 km s⁻¹), highly collimated, massive (0.03M _⊙) bipolar outflow with a very large scalar momentum (about 10³⁹ g cm s⁻¹), and the characteristic spatio-kinematic structure of bow-shocks at the tips of this outflow. The fast outflow in I 22036, as in most PPNe, cannot be driven by radiation pressure. The large mass of the torus suggests that it has most likely resulted from common-envelope evolution in a binary, however it remains to be seen whether or not the time-scales required for the growth of grains to millimeter sizes in the torus are commensurate with such a formation scenario. The presence of the torus should facilitate the formation of the accretion disk needed to launch the jet. We also find that the ¹³C/¹²C ratio in I 22036 is very high (0.16), close to the maximum value achieved in equilibrium CNO-nucleosynthesis (0.33). The combination of the high circumstellar mass (i.e., in the torus and an extended dust shell inferred from ISO far-infrared spectra) and the high ¹³C/¹²C ratio in I 22036 provides strong support for this object having evolved from a massive (≳4M _⊙) progenitor in which hot-bottom-burning has occurred.     

Tidal torques and galactic warps

%We study the evolution of galactic disks subject to tidal torques motivated by cosmological N-body simulations using analytic and numerical techniques. We find that self-gravitating disks subject to these torques resemble observed warped galaxies. The warps develop at a local surface density of 70 M _⊙ pc^-2 and move out through the disk at a rate that depends on the surface density of the disk. This revised version was published online in August 2006 with corrections to the Cover Date.     

Disk galaxy evolution: from the Milky Way to high-redshift disks

We develop a detailed model of the Milky Way (a `prototypical' disk galaxy) and extend it to other disks with the help of some simple scaling relations, obtained in the framework of Cold Dark Matter models. This phenomenological (`hybrid') approach to the study of disk galaxy evolution allows us to reproduce successfully a large number of observed properties of disk galaxies in the local Universe and up to redshift z ∼ 1. The important conclusion is that, on average, massive disks have formed the bulk of their stars earlier than their lower mass counterparts: the `star formation hierarchy' has been apparently opposite to the `dark matter assembly' hierarchy. It is not yet clear whether `feedback' (as used in semi-analytical models of galaxy evolution) can explain that discrepancy. This revised version was published online in August 2006 with corrections to the Cover Date.     

Evidence that Synchrotron Emission from Nonthermal Electrons Produces the Increasing Submillimeter Spectral Component in Solar Flares

We investigate the origin of the increasing spectra observed at submillimeter wavelengths detected in the flare on 2 November 2003 starting at 17:17 UT. This flare, classified as an X8.3 and 2B event, was simultaneously detected by RHESSI and the Solar Submillimeter Telescope (SST) at 212 and 405 GHz. Comparison of the time profiles at various wavelengths shows that the submillimeter emission resembles that of the high-energy X rays observed by RHESSI whereas the microwaves observed by the Owens Valley Solar Array (OVSA) resemble that of ∼50 keV X rays. Moreover, the centroid position of the submillimeter radiation is seen to originate within the same flaring loops of the ultraviolet and X-ray sources. Nevertheless, the submillimeter spectra are distinct from the usual microwave spectra, appearing to be a distinct spectral component with peak frequency in the THz range. Three possibilities to explain this increasing radio spectra are discussed: (1) gyrosynchrotron radiation from accelerated electrons, (2) bremsstrahlung from thermal electrons, and (3) gyrosynchrotron emission from the positrons produced by pion or radioactive decay after nuclear interactions. The latter possibility is ruled out on the grounds that to explain the submillimeter observations requires 3000 to 2×10⁵ more positrons than what is inferred from X-ray and γ-ray observations. It is possible to model the emission as thermal; however, such sources would produce too much flux in the ultraviolet and soft X-ray wavelengths. Nevertheless we are able to explain both spectral components at microwave and submillimeter wavelengths by gyrosynchrotron emission from the same population of accelerated electrons that emit hard X rays and γ rays. We find that the same 5×10³⁵ electrons inferred from RHESSI observations are responsible for the compact submillimeter source (0.5 arcsec in radius) in a region of 4500 G low in the atmosphere, and for the traditional microwave spectral component by a more extended source (50 arcsec) in a 480 G magnetic field located higher up in the loops. The extreme values in magnetic field and source size required to account for the submillimeter emission can be relaxed if anisotropy and transport of the electrons are taken into account.     

Submillimeter Studies of Prestellar Cores and Protostars: Probing the Initial Conditions for Protostellar Collapse

Improving our understanding of the earliest stages of star formation is crucial to gain insight into the origin of stellar masses, multiple systems, and protoplanetary disks. We discuss recent advances made in this area, thanks to submillimeter mapping observations with large single-dish telescopes and interferometers. Although ambipolar diffusion appears to be too slow cores, there is nevertheless good evidence that the gravitational collapse of isolated protostellar cores is strongly magnetically controlled. We also argue that the beginning of protostellar collapse is much more violent in cluster-forming clouds than in regions of distributed star formation.     

Disks Around Post-Main Sequence Binaries

We propose that at least two stars on or near the AGB have long-lived orbiting disks: HD 44179, the central star in the Red Rectangle, and BM Gem, a carbon-rich star with an oxygen-rich circumstellar envelope. The CO emission from both of these disks has a spike with a width near ∼2 km s⁻¹, indicating disk radii of ∼10¹⁶ cm. The dust in such disks is therefore quite cold (near T ∼ 50 K for the Red Rectangle) and may emit primarily at submillimeter wavelengths. The disks around stars where there is also substantial mass loss may not be easily observable; there could be many as yet undiscovered disks around AGB stars This revised version was published online in September 2006 with corrections to the Cover Date.     

A new evolutionary scenario of intermediate-mass star-formation revealed by multi-wavelength observations of OMC-2/3

We have performed millimeter- and submilli- meter-wave survey observations using the Nobeyama millimeter array (NMA) and the Atacama Submillimeter Telescope Experiment (ASTE) in one of the nearest intermediate-mass (IM) star-forming regions: Orion Molecular Cloud-2/3 (OMC-2/3). Using the high-resolution capabilities offered by the NMA (∼several arcsec), we observed dust continuum and H¹³CO⁺(1–0) emission in 12 pre- and proto-stellar candidates identified previously in single-dish millimeter observations. We unveiled the evolutionary changes with variations of the morphology and velocity structure of the dense envelopes traced by the H¹³CO⁺(1–0) emission. Furthermore, using the high-sensitivity capabilities offered by the ASTE, we searched for large-scale molecular outflows associated with these pre- and proto-stellar candidates observed with the NMA. As a result of the CO(3–2) observations, we detected six molecular outflows associated with the dense gas envelopes traced by H¹³CO⁺(1–0) and 3.3 mm continuum emission. The estimated CO outflow momentum increases with the evolutionary sequence from early to late type of the protostellar cores. We also found that the 24 μm flux increases as the dense gas evolutionary sequence. We propose that the enhancement of the 24 μm flux is caused by the growth of the cavity (i.e. the CO outflow destroys the envelope) as the evolutionary sequence. Our results show that the dissipation of the dense gas envelope plays an essential role in the evolution of the IM protostars. The extremely high-sensitivity and high-angular resolution offered by ALMA will reveal unprecedented details of the inner ∼50 AU of these protostars, which will provide us a break through in the classic scenario of IM star/disk formation.     

Galaxy evolution from a microJansky radio survey

A strong evolution of galaxies is observed for 0<z<1, as evidenced by an increase of almost an order of magnitude in the galaxy star-formation rate density. However, it is known that dust obscuration has affected our understanding of galaxy evolution over this significant fraction of the age of the Universe. In order to study galaxy evolution free from dust induced biases, an ultra deep radio survey – the Phoenix Deep Survey – was initiated. With a detection limit of 60μJy, this homogeneous survey, complemented with multiwavelength (photometric and spectroscopic) observations, is being used to build a consistent picture of galaxy evolution. The ultra-deep radio source counts are presented, and interpreted using luminosity function evolutionary models. The discovery of extremely dusty galaxies from this survey, along with the clustering properties of the sub-mJy radio population, are also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Jet-driven disk accretion in low luminosity AGN?

We explore an accretion model for low luminosity AGN (LLAGN) that attributes the low radiative output to a low mass accretion rate, , rather than a low radiative efficiency. In this model, electrons are assumed to drain energy from the ions as a result of collisionless plasma microinstabilities. Consequently, the accreting gas collapses to form a geometrically thin disk at small radii and is able to cool before reaching the black hole. The accretion disk is not a standard disk, however, because the radial disk structure is modified by a magnetic torque which drives a jet and which is primarily responsible for angular momentum transport. We also include relativistic effects. We apply this model to the well known LLAGN M87 and calculate the combined disk-jet steady-state broadband spectrum. A comparison between predicted and observed spectra indicates that M87 may be a maximally spinning black hole accreting at a rate of ∼10⁻³ M _⊙ yr⁻¹. This is about 6 orders of magnitude below the Eddington rate for the same radiative efficiency. Furthermore, the total jet power inferred by our model is in remarkably good agreement with the value independently deduced from observations of the M87 jet on kiloparsec scales.     

Jet-driven disk accretion in low luminosity AGN?*

We explore an accretion model for low luminosity AGN (LLAGN) that attributes the low radiative output to a low mass accretion rate, , rather than a low radiative efficiency. In this model, electrons are assumed to drain energy from the ions as a result of collisionless plasma microinstabilities. Consequently, the accreting gas collapses to form a geometrically thin disk at small radii and is able to cool before reaching the black hole. The accretion disk is not a standard disk, however, because the radial disk structure is modified by a magnetic torque which drives a jet and which is primarily responsible for angular momentum transport. We also include relativistic effects. We apply this model to the well known LLAGN M87 and calculate the combined disk-jet steady-state broadband spectrum. A comparison between predicted and observed spectra indicates that M87 may be a maximally spinning black hole accreting at a rate of ∼10⁻³  M _⊙ yr⁻¹. This is about 6 orders of magnitude below the Eddington rate for the same radiative efficiency. Furthermore, the total jet power inferred by our model is in remarkably good agreement with the value independently deduced from observations of the M87 jet on kiloparsec scales. ^* This paper has previously been published in Astrophysics and Space Science, vol. 310:3–4.     

Unveiling the chemistry of hot protostellar cores with ALMA

High angular resolution mm-wave observations of the Orion-KL region, made with the IRAM Plateau de Bure interferometer (PdBI), reveal the presence of several cores of size 10³ AU, which have distinct spectral signatures. Complex molecules such as ethanol, vinyl cyanide and dimethyl ether show different distributions and their relative abundance varies from core to core by orders of magnitude. The molecular column densities derived in the cores also differ widely from the beam-averaged column densities observed with large single-dish telescopes. Obviously, the predictions of hot core chemistry models must be checked against high resolution observations. ALMA, which allies sensitivity and high angular resolution, will be a key instrument for this type of studies. The PdBI observations were part of a search for interstellar glycine, also carried out with the IRAM 30-m telescope and the Green Bank Telescope. We derive a 3σ upper limit on the column density of glycine of 1×10¹⁵ cm⁻² per 2″×3″ beam in the Orion Hot Core and Compact Ridge. Based on observations made with the IRAM PdB Interferometer, the IRAM 30-m telescope and the NRAO Green-Bank telescope. IRAM is supported by CNRS, MPG and IGN.     

Accretion of helium and metal-rich gas onto neutron stars and black holes at high luminosities

Ultraluminous X-ray sources fed by Wolf-Rayet star winds and X-ray bursters in ultracompact binaries with He or C white dwarfs have accretion disks whose properties may differ significantly from those of pure Hα-accretion disks. Therefore, we have included the dependence on charge number Z and mean molecular weights μ_e/I in the Shakura and Sunyaev (1973) scaling relations for the key parameters of the disk. Furthermore, we also consider the pseudo-Newtonian potential of Paczyńsky and Wiita (1980). These scaling relations might become useful, e.g., when estimating the illumination efficiency of the outer parts of the disk. We also address the changes in the structure of the boundary (spreading) layer on the surfaces of neutron stars that occurs in the case of H-depleted accretion disks. Published in Russian in Pis’ma v Astronomicheskiĭ Zhurnal, 2006, Vol. 32, No. 4, pp. 288–294. The article was translated by the authors.     

CDM, Feedback and the Hubble Sequence

We have performed TreeSPH simulations of galaxy formation in a standard ΛCDM cosmology, including effects of star formation, energetic stellar feedback processes and a meta-galactic UV field, and obtain a mix of disk, lenticular and elliptical galaxies. The disk galaxies are deficient in angular momentum by only about a factor of two compared to observed disk galaxies. The stellar disks have approximately exponential surface density profiles, and those of the bulges range from exponential to r ^1/4, as observed. The bulge-to-disk ratios of the disk galaxies are consistent with observations and likewise are their integrated B-V colours, which have been calculated using stellar population synthesis techniques. Furthermore, we can match the observed I-band Tully-Fisher (TF) relation, provided that the mass-to-light ratio of disk galaxies is (M/L _I) ≃ 0.6–0.7. The ellipticals and lenticulars have approximately r ^1/4 stellar surface density profiles, are dominated by non-disklike kinematics and flattened due to non-isotropic stellar velocity distributions, again consistent with observations. This revised version was published online in September 2006 with corrections to the Cover Date.     

Magnetic fields in protoplanetary disks

Magnetic fields likely play a key role in the dynamics and evolution of protoplanetary disks. They have the potential to efficiently transport angular momentum by MHD turbulence or via the magnetocentrifugal acceleration of outflows from the disk surface. Magnetically-driven mixing has implications for disk chemistry and evolution of the grain population, and the effective viscous response of the disk determines whether planets migrate inwards or outwards. However, the weak ionisation of protoplanetary disks means that magnetic fields may not be able to effectively couple to the matter. I examine the magnetic diffusivity in a minimum solar nebula model and present calculations of the ionisation equilibrium and magnetic diffusivity as a function of height from the disk midplane at radii of 1 and 5 AU. Dust grains tend to suppress magnetic coupling by soaking up electrons and ions from the gas phase and reducing the conductivity of the gas by many orders of magnitude. However, once grains have grown to a few microns in size their effect starts to wane and magnetic fields can begin to couple to the gas even at the disk midplane. Because ions are generally decoupled from the magnetic field by neutral collisions while electrons are not, the Hall effect tends to dominate the diffusion of the magnetic field when it is able to partially couple to the gas, except at the disk surfaces where the low density of neutrals permits the ions to remain attached to the field lines. For a standard population of 0.1 μm grains the active surface layers have a combined column Σ_active≈2 g cm⁻² at 1 AU; by the time grains have aggregated to 3 μm, Σ_active≈80 g cm⁻². Ionisation in the active layers is dominated by stellar X-rays. In the absence of grains, X-rays maintain magnetic coupling to 10% of the disk material at 1 AU (i.e. Σ_active≈150 g cm⁻²). At 5 AU the Σ_active≈Σ_total once grains have aggregated to 1 μm in size.     

High angular resolution millimeter observations of circumstellar disks

In this lecture, we review the properties of protoplanetary disks as derived from high angular resolution observations at millimeter wavelengths. We discuss how the combination of several different high angular resolution techniques allow us to probe different regions of the disk around young stellar objects and to derive the properties of the dust when combined with sophisticated disk models. The picture that emerges is that the dust in circumstellar disks surrounding pre-main sequence stars is in many cases significantly evolved compared to the dust in molecular clouds and the interstellar medium. It is however still difficult to derive a consistent picture and timeline for dust evolution in disks as the observations are still limited to small samples of objects.We also review the evidence for and properties of disks around high-mass young stellar objects and the implications on their formation mechanisms. The study of massive YSOs is complicated by their short lifetimes and larger average distances. In most cases high angular resolution data at millimeter wavelengths are the only method to probe the structure of disks in these objects.We provide a summary of the characteristics of available high angular resolution millimeter and submillimeter observatories. We also describe the characteristics of the ALMA observatory being constructed in the Chilean Andes. ALMA is going to be the world leading observatory at millimeter wavelengths in the coming decades, the project is now in its main construction phase with early science activities envisaged for 2010 and full science operations for 2012.     

A Study on the Characteristics of the Structure of Vega's Debris Disk

The clumpy structure in the Vega's debris disk was reported at millimeter wavelengths previously, and attributed to the concentration of dust grains trapped in resonances with a potential high-eccentricity planet. However, current imaging at multi-wavelengths with higher sensitivity indicates that the Vega's debris disk has a smooth structure. But a planet orbiting Vega could not be neglected, and the present-day observations may place a severe constraint on the orbital parameters for the potential planet. Herein, we utilize the modi- fied MERCURY codes to numerically simulate the Vega system, which consists of a debris disk and a planet. In our simulations, the initial inner and outer boundaries of the debris disk are assumed to be 80 AU and 120 AU, respectively. The dust grains in the disk have the sizes from 10 μm to 100 μm, and the nearly coplanar orbits. From the outcomes, we show that the evolution of debris disk is consistent with recent observations, if there is no planet orbiting Vega. However, if Vega owns a planet with a high eccentricity (e.g., e = 0.6), the planet's semi- major axis cannot be larger than 60 AU, otherwise, an aggregation phenomenon will occur in the debris disk due to the existence of the postulated planet. In addition, the 2:1 mean motion resonances may play a significant role in forming the structure of debris disk.     

Evolutionary processes in protoplanetary accretion disks: the propagation of axisymmetric shock waves

In this paper we investigate both the global and the local hydrodynamics of axisymmetric accretion disks around young stellar objects under the simultaneous action of viscosity, self-gravity and pressure forces. For simplicity, we take for the global model a polytropic equation of state, make the infinitely thin disk approximation and characterize the surface density and temperature profiles in the disk as power laws in the radial distance r from the protostar. We solve the problem of the general density profile of a Keplerian disk showing that self-gravity could not be an important factor for the fast formation of the rocky cores of giant gaseous planets in our solar system. Under the hypothesis that the unperturbed rotation curve of the disk is nearly Keplerian throughout the radial extent, we can estimate with our polytropic model a lower limit for the resulting masses M_d(r) of stable disks up to 100 AU. These masses are in the range of the so-called minimum mass solar nebular (_d/M_s ≈ 0.01–0.02).By adopting a simplified viscosity model, where the height-integrated turbulent dynamical viscosity ν is a function of the surface density σ like η ∝ σ^Γ, we derive in the local shearing sheet model linearized evolution equations for small density perturbations describing both a diffusion process and the propagation of acoustic density waves. We solve a special initial value problem and calculate the appropriate Green's function. The analytical solutions so obtained describe in the case Γ < 0 the successive formation of quasi-stationary ring-shaped density structures in a disk with a definite mode of maximum instability, whereas in the case Γ > Γ_c the density wave equation describes the propagation of an “overstable” ring-shaped acoustic density wavelet to the outer ranges of the accretion disk. Whereas the group velocity of the wave packet is subsonic, the phase velocities of individual wave crests in the wave packet are supersonic. The mode of maximum instability, the growth rate and the number of growing waves in the wavelet are controlled by Γ and α. Our present knowledge concerning turbulent viscosity in protoplanetary disks is not sufficient to decide whether or not the case Γ > Γ_c is realized.The suggested structuring processes in the linear theory should initiate in the non-linear regime the formation of narrow ring-shaped density shock waves moving through the protoplanetary disk. These non-linear waves could produce extremely spatially and temporally heterogeneous temperature regions in the disk. We speculate that ring-shaped density waves, excited by inner boundary conditions and which have dominated the disk's evolution at early times, are responsible both for the fast growth of dust to planetesimals and at least for the rapid accretion of the rocky cores of giant gaseous planets in the protoplanetary accretion disk (shock wave trigger hypothesis). We derive provisional scaling rules for planetary systems regarding the spacing of orbits as a function of the mass ratio of the protoplanetary disk to the protostar. However, further analytical work and linear as well as nonlinear numerical simulations of density waves excited by inner boundary conditions are needed to consolidate the results and speculations of our linear wave mechanics in the future.     

Polarization of radiation from a strongly magnetized accretion disk: The asymptotic spectral distribution

We calculate the polarization of the radiation from an optically thick accretion disk with a vertical averaged magnetic field. The polarization arises from the scattering of light by free electrons in a magnetized disk plasma. The Faraday rotation of the polarization plane during the propagation of a photon in a medium with a magnetic field is considered as the main effect. We discuss various models of optically thick accretion disks with a vertical averaged magnetic field. Our main goal is to derive simple asymptotic formulas for the polarization of radiation in the case where the Faraday rotation angle Ψ ≫ 1 at the Thomson optical depth τ = 1. The results of our calculations allow the magnetic field strength in the region of the marginally stable orbit near a black hole to be estimated from polarimetric observations, including X-ray observations expected in the future. Since the polarization spectrum of the radiation strongly depends on the accretion disk model, a realistic physical model of the accretion disk can be determined from data on the polarization of its radiation.