Complex organic molecules in protoplanetary disks

Context. Protoplanetary disks are vital objects in star and planet formation, possessing all the material, gas and dust, which may form a planetary system orbiting the new star. Small, simple molecules have traditionally been detected in protoplanetary disks; however, in the ALMA era, we expect the molecular inventory of protoplanetary disks to significantly increase.

Aims. We investigate the synthesis of complex organic molecules (COMs) in protoplanetary disks to put constraints on the achievable chemical complexity and to predict species and transitions which may be observable with ALMA.

Methods. We have coupled a 2D steady-state physical model of a protoplanetary disk around a typical T Tauri star with a large gas-grain chemical network including COMs. We compare the resulting column densities with those derived from observations and perform ray-tracing calculations to predict line spectra. We compare the synthesised line intensities with current observations and determine those COMs which may be observable in nearby objects. We also compare the predicted grain-surface abundances with those derived from cometary comae observations.

Results. We find COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances similar to 10(-6)-10(-4) that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, similar to 10(-12)-10(-7). Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards T Tauri star-disk systems. There is poor agreement with HC3(N) lines observed towards LkCa 15 and GO Tau and we discuss possible explanations for these discrepancies. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging, even with ALMA "Full Science" capabilities. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

Cumulenic and heterocumulenic anions : potential interstellar species?

A recent theoretical investigation by Terzieva & Herbst of linear carbon chains, C-n where n greater than or equal to 6, in the interstellar medium has shown that these species can undergo efficient radiative association to form the corresponding anions. An experimental study by Barckholtz, Snow & Bierbaum of these anions has demonstrated that they do not react efficiently with molecular hydrogen, leading to the possibility of detectable abundances of cumulene-type anions in dense interstellar and circumstellar environments. Here we present a series of electronic structure calculations which examine possible anionic candidates for detection in these media, namely the anion analogues of the previously identified interstellar cumulenes CnH and Cn-1CH2 and heterocumulenes CnO (where n = 2-10). The extraordinary electron affinities calculated for these molecules suggest that efficient radiative electron attachment could occur, and the large dipole moments of these simple (generally) linear molecules point to the possibility of detection by radio astronomy.

Modelling the warm H(2) infrared emission of the Helix nebula cometary knots

Molecular hydrogen emission is commonly observed in planetary nebulae. Images taken in infrared H(2) emission lines show that at least part of the molecular emission is produced inside the ionized region. In the best studied case, the Helix nebula, the H(2) emission is produced inside cometary knots (CKs), comet-shaped structures believed to be clumps of dense neutral gas embedded within the ionized gas. Most of the H(2) emission of the CKs seems to be produced in a thin layer between the ionized diffuse gas and the neutral material of the knot, in a mini-photodissociation region (mini-PDR). However, PDR models published so far cannot fully explain all the characteristics of the H(2) emission of the CKs. In this work, we use the photoionization code AANGABA to study the H(2) emission of the CKs, particularly that produced in the interface H(+)/H(0) of the knot, where a significant fraction of the H(2) 1-0 S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such a region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H(2) infrared lines. We find that the temperature derived from H(2) observations, even of a single knot, will depend very strongly on the observed transitions, with much higher temperatures derived from excited levels. We also proposed that the separation between the H alpha and [N II] peak emission observed in the images of CKs may be an effect of the distance of the knot from the star, since for knots farther from the central star the [N II] line is produced closer to the border of the CK than H alpha.

Silicon monosulphide radiative association

Rate coefficients for radiative association of silicon and sulphur atoms to form silicon monosulphide (SiS) molecule are estimated. The radiative association is due mainly to approach in the E(1)Sigma(+) and A(1)Pi states of SiS. For temperatures ranging from similar to 1000 to similar to 14 000 K, the rate coefficients are found to vary from 8.43 x 10(-17) to 2.69 x 10(-16) cm(3) s(-1). Our calculated rate coefficient is higher than the values used in modelling the chemistry of Type Ia supernovae.

The formation of CN and CN+ by direct radiative association

Rate coefficients for direct radiative association of carbon and nitrogen atoms to form CN, and of carbon ions and nitrogen atoms to form CN+ ions, are calculated for temperatures in the range of 300 to 14,700 K. For the CN molecule, the rate coefficients can be represented by the standard expression, k(CN)(T) = 7.87 x 10(-19)(T/300)(0.056) exp (-96.0/T) cm(3) s(-1) for temperatures between 300 and 2700 K and k(CN)(T) = 1.37 x 10(-18)(T/300)-0.128 exp (-520.1/T) cm(-3) s(-1) at T > 2700 K. For the CN+ ion, the corresponding expression is k(CN+)(T) = 1.08 x 10(-18)(T/300)(0.071) exp (-57.5/T) cm(-3) s(-1) for the temperature range studied. Calculated rate coefficients k(CN) are about 2 orders of magnitude lower than the canonical value used in the modeling of the chemistry of various astrophysical environments.

SPH simulations of clumps formation by dissipative collisions of molecular clouds - II. Magnetic case

We performed computer simulations of interstellar cloud-cloud collisions using the three-dimensional smoothed particle magnetohydrodynamics method. In order to study the role of the magnetic field on the process of collision-triggered fragmentation, we focused our attention on head-on supersonic collisions between two identical spherical molecular-clouds. Two extreme configurations of the magnetic field were adopted: parallel and perpendicular to the initial clouds motion. The initial magnetic field strength was approximately 12.0 muG. In the parallel case, much more of the collision debris were retained in the shocking region than in the non-magnetic case where gas escaped freely throughout the symmetry plane. Differently from the non-magnetic case, eddy-like vortices were formed. The regions of highest vorticity and the the regions of highest density are offset. We found clumps formation only in the parallel case, however, they were larger, hotter and less dense than in the analogous non-magnetic case. In the perpendicular case, the compressed field works as a magnetic wall, preventing a stronger compression of the colliding clouds. This last effect inhibits direct contact of the two clouds. In both cases, we found that the field lines show a chaotic aspect in large scales. Also, the field magnitude is considerably amplified in the shock layer. However, the field distribution is almost coherent in the higher density regions.

SPH simulations of clumps formation by dissipative collision of molecular clouds - I. Non magnetic case

Computer experiments of interstellar cloud collisions were performed with a new smoothed-particle-hydrodynamics (SPH) code. The SPH quantities were calculated by using spatially adaptive smoothing lengths and the SPH fluid equations of motion were solved by means of a hierarchical multiple time-scale leapfrog. Such a combination of methods allows the code to deal with a large range of hydrodynamic quantities. A careful treatment of gas cooling by H, H(2), CO and H II, as well as a heating mechanism by cosmic rays and by H(2) production on grains surface, were also included in the code. The gas model reproduces approximately the typical environment of dark molecular clouds. The experiments were performed by impinging two dynamically identical spherical clouds onto each other with a relative velocity of 10 km s(-1) but with a different impact parameter for each case. Each object has an initial density profile obeying an r(-1)-law with a cutoff radius of 10 pc and with an initial temperature of 20 K. As a main result, cloud-cloud collision triggers fragmentation but in expense of a large amount of energy dissipated, which occurred in the head-on case only. Off-center collision did not allow remnants to fragment along the considered time (similar to 6 Myr). However, it dissipated a considerable amount of orbital energy. Structures as small as 0.1 pc, with densities of similar to 10(4) cm(-3), were observed in the more energetic collision.

Formation of SO, SO+, and S-2 by radiative association

Rate coefficients for radiative association of SO, SO+, and S-2 are estimated. For temperatures ranging from 300 to 14,000 K, the direct radiative association rate coefficients are found to vary with temperature from 1.73 x 10(-19) to 7.29 x 10(-19) cm(3) s(-1) and from 1.49 x 10(-21) to 3.70 x 10(-19) cm(3) s(-1) for S-2 and SO, respectively. The rate coefficients for formation through the inverse predissociation for S-2 are found to vary from 3.59 x 10(-18) to 1.44 x 10(-20) cm(3) s(-1). For SO+, the direct rate coefficient varies rapidly with temperature from 3.62 x 10(-27) cm(3) s(-1) at 2000 K to 2.34 x 10(-20) cm(3) s(-1) at 14,000 K. The direct radiative association rate coefficients increase with the increase in temperature, but the inverse predissociation rate coefficients decrease with the increase in temperature.

A theoretical study of the A(2)Sigma(+)-X-2 Pi system of the SiP molecule

The A (2)Sigma(+) and X(2)Pi electronic states of the SiP species have been investigated theoretically at a very high level of correlation treatment (CASSCF/MRSDCI). Very accurate potential energy curves are presented for both states, as well as the associated spectroscopic constants as derived from the vib-rotational energy levels determined by means of the numerical solution of the radial Schrodinger equation. Electronic transition moment function, oscillator strengths, Einstein coefficients for spontaneous emission, and Franck-Condon factors for the A(2)Sigma(+)-X(2)Pi system have been calculated. Dipole moment functions and radiative lifetimes for both states have also been determined. Spin-orbit coupling constants are also reported. The radiative lifetimes for the A(2)Sigma(+) state, taking into account the spin-orbit diagonal correction to the X(2)Pi state, decrease from a value of 138 ms at v' = 0 to 0.48 ms at v' = 8, and, for the X(2)Pi state, from 2.32 s at v = 1 to 0.59 s at v = 5. Vibrational and rotational transitions are expected to be relatively strong.

Radiative association of C and P, and Si and P atoms

The rate coefficients for the formation of carbon monophosphide (CP) and silicon monophosphide (SiP) by radiative association are estimated for temperatures ranging from 300 to 14 100 K. In this temperature range, the radiative association rate coefficients are found to vary from 1.14 x 10(-18) to 1.62 x 10(-18) cm(3) s(-1) and from 3.73 x 10(-20) to 7.03 x 10(-20) cm(3) s(-1) for CP and SiP, respectively. In both cases, rate coefficients increase slowly with the increase in temperature.

Formation of SiC by radiative association

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Radiative association of Ti and O atoms

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

A theoretical study of the A2 ∑+-X2∏ system of the SiP molecule

The A2∑+ and Z2∏ electronic states of the SiP species have been investigated theoretically at a very high level of correlation treatment (CASSCF/MRSDCI). Very accurate potential energy curves are presented for both states, as well as the associated spectroscopic constants as derived from the vib-rotational energy levels determined by means of the numerical solution of the radial Schrödinger equation. Electronic transition moment function, oscillator strengths, Einstein coefficients for spontaneous emission, and Franck-Condon factors for the A2∑+-X2∏ system have been calculated. Dipole moment functions and radiative lifetimes for both states have also been determined. Spin-orbit coupling constants are also reported. The radiative lifetimes for the A2∑+ state, taking into account the spin-orbit diagonal correction to the X2∏ state, decrease from a value of 138 ms at v′ = 0 to 0.48 ms at v′ = 8, and, for the X2∏ state, from 2.32 s at v″ = 1 to 0.59 s at v″ = 5. Vibrational and rotational transitions are expected to be relatively strong.

The formation of AlF by radiative association

The formation of the aluminium monofluoride molecule AlF by radiative association of the Al and F atoms is estimated. The radiative association of Al(P-2) and F(P-2) atoms is found to be dominated by the approach along the A(1) potential energy curve accompanied by spontaneous emission into the X-1 Sigma(+) ground state of the AlF. For temperatures ranging from 300 to 14 000 K, the rate coefficients are found to vary from 1.35 x 10(-17) to 9.31 x 10(-16) cm(3) s(-1), respectively. These values indicate that only a small amount of AlF molecules can be formed by radiative association in the inner envelope of carbon-rich stars and other hostile environments.

Radiative association of Ti and O atoms

Rate coefficients for the radiative association of titanium and oxygen atoms to form the titanium monoxide (TiO) molecule are estimated. The radiative association of Ti(F-3) and O(P-3) atoms is dominated by an approach along the C-3 Delta potential energy curve, accompanied by spontaneous emission into the X-3 Delta ground state of TiO. For temperatures ranging from 300-14 000 K, the total rate coefficients are found to vary from 4.76 x 10(-17) to 9.96 x 10(-17) cm(3) s(-1), respectively.