Chemical Processes in Protoplanetary Disks

We have developed a high-resolution combined physical and chemical model of a protoplanetary disk surrounding a typical T Tauri star. Our aims were to use our model to calculate the chemical structure of disks on small scales (submilliarcsecond in the inner disk for objects at the distance of Taurus, ~140 pc) to investigate the various chemical processes thought to be important in disks and to determine potential molecular tracers of each process. Our gas-phase network was extracted from the UMIST Database for Astrochemistry to which we added gas–grain interactions including freezeout and thermal and non-thermal desorption (cosmic-ray-induced desorption, photodesorption, and X-ray desorption), and a grain-surface network. We find that cosmic-ray-induced desorption has the least effect on our disk chemical structure while photodesorption has a significant effect, enhancing the abundances of most gas-phase molecules throughout the disk and affecting the abundances and distribution of HCN, CN, and CS, in particular. In the outer disk, we also see enhancements in the abundances of H2O and CO2. X-ray desorption is a potentially powerful mechanism in disks, acting to homogenize the fractional abundances of gas-phase species across the depth and increasing the column densities of most molecules, although there remain significant uncertainties in the rates adopted for this process. The addition of grain-surface chemistry enhances the fractional abundances of several small complex organic molecules including CH3OH, HCOOCH3, and CH3OCH3 to potentially observable values (i.e., a fractional abundance of greater than 10-11).

"Irish Perspectives on EC Law by Mary Catherine Lucey and Cathrina Keville (eds.), Dublin, Thomson Round Hall, 2003 (ISBN 1-85800-280-X)"

Book Review

A commentary on climate change, stone decay dynamics and the ‘greening’ of natural stone buildings: new perspectives on ‘deep wetting’

Environmental controls on stone decay processes are rapidly changing as a result of changing climate. UKCP09 projections for the 2020s (2010–2039) indicate that over much of the UK seasonality of precipitation will increase. Summer dryness and winter wetness are both set to increase, the latter linked to projected precipitation increases in autumn and spring months. If so, this could increase the time that stone structures remain wet and possibly the depth of moisture penetration, and it appears that building stone in Northern Ireland has already responded through an increased incidence of algal ‘greening’.This paper highlights the need for understanding the effects of climate change through a series of studies of largely sandstone structures. Current and projected climatic trends are therefore considered to have aesthetic, physical and chemical implications that are not currently built into our models of sandstone decay, especially with respect to the role played by deep-seated wetness on sandstone deterioration and decay progression and the feedbacks associated with, for example surface algal growth. In particular,it is proposed that algal biofilms will aid moisture retention and further facilitate moisture and dissolved salt penetration to depth. Thus, whilst the outer surface of stone may continue to experience frequent wetting and drying associated with individual precipitation events, the latter is less likely to be complete, and the interiors of building blocks may only experience wetting/drying in response to seasonal cycling. A possible consequence of deeper salt penetration could be a delay in the onset of surface deterioration,but more rapid and effective retreat once it commences as decay mechanisms ‘tap into a reservoir of deep salt’.

Chemical evolution of protoplanetary disks - the effects of viscous accretion, turbulent mixing and disk winds

We calculate the chemical evolution of protoplanetary disks considering radial viscous accretion, vertical turbulent mixing, and vertical disk winds. We study the effects on the disk chemical structure when different models for the formation of molecular hydrogen on dust grains are adopted. Our gas-phase chemistry is extracted from the UMIST Database for Astrochemistry (Rate06) to which we have added detailed gas-grain interactions. We use our chemical model results to generate synthetic near- and mid-infrared local thermodynamic equilibrium line emission spectra and compare these with recent Spitzer observations. Our results show that if H2 formation on warm grains is taken into consideration, the H2O and OH abundances in the disk surface increase significantly. We find that the radial accretion flow strongly influences the molecular abundances, with those in the cold midplane layers particularly affected. On the other hand, we show that diffusive turbulent mixing affects the disk chemistry in the warm molecular layers, influencing the line emission from the disk and subsequently improving agreement with observations. We find that NH3, CH3OH, C2H2, and sulfur-containing species are greatly enhanced by the inclusion of turbulent mixing. We demonstrate that disk winds potentially affect the disk chemistry and the resulting molecular line emission in a manner similar to that found when mixing is included.

Electron-impact excitation of Ni II: Effective collision strengths for optically allowed fine-structure transitions

In this paper, we present collision strengths and Maxwellian averaged effective collision strengths for the electron-impact excitation of Ni II. Attention is expressly concentrated on the optically allowed fine-structure transitions between the 3d 9, 3d 84s, and 3d 74s 2 even parity levels and the 3d 84p and 3d 74s 4p odd parity levels. The parallel RMATRXII R-matrix package has been recently extended to allow for the inclusion of relativistic fine-structure effects. This suite of codes has been utilized in conjunction with the parallel PSTGF and PSTGICF programs in order to compute converged total collision strengths for the allowed transitions with which this study is concerned. All 113 LS terms identified with the 3d 9, 3d 84s, 3d 74s 2, 3d 84p, and 3d 74s 4p basis configurations were included in the target wavefunction representation, giving rise to a sophisticated 295 jj-level, 1930 coupled channel scattering complex. Maxwellian averaged effective collision strengths have been computed at 30 individual electron temperatures ranging from 30 to 1,000,000 K. This range comfortably encompasses all temperatures significant to astrophysical and plasma applications. The convergence of the collision strengths is exhaustively investigated and comparisons are made with previous theoretical works, where significant discrepancies exist for the majority of transitions. We conclude that intrinsic in achieving converged collision strengths and thus effective collision strengths for the allowed transitions is the combined inclusion of contributions from the (N + 1) partial waves extending to a total angular momentum value of L = 50 and further contributions from even higher partial waves accomplished by employing a "top-up" procedure.

Electron-impact excitation of CrII: A theoretical calculation of effective collision strengths for optically allowed transitions.

In this paper, we present electron-impact excitation collision strengths and Maxwellian averaged effective collision strengths for the complicated iron-peak ion Cr II. We consider specifically the allowed lines for transitions from the 3d(5) and 3d(4)4s even parity configuration states to the 3d(4)4p odd parity configuration levels. The parallel suite of R-Matrix packages, RMATRX II, which have recently been extended to allow for the inclusion of relativistic effects, were used to compute the collision cross sections. A total of 108 LS pi/280 J pi levels from the basis configurations 3d(5), 3d(4)4s, and 3d(4)4p were included in the wavefunction representation of the target including all doublet, quartet, and sextet terms. Configuration interaction and correlation effects were carefully considered by the inclusion of seven more configurations and a pseudo-corrector (4d) over bar type orbital. The 10 configurations incorporated into the Cr II model thus listed are 3d(5), 3d(4)4s, 3d(4)4p, 3d(3)4s(2), 3d(3)4p(2), 3d(3)4s4p, 3d(4)(4d) over bar, 3d(3)4s (4d) over bar, 3d(3)4p (4d) over bar, and 3d(3)(4d) over bar (2), constituting the largest Cr II target model considered to date in a scattering calculation. The Maxwellian averaged effective collision strengths are computed for a wide range of electron temperatures 2000-100,000 K which are astrophysically significant. Care has been taken to ensure that the partial wave contributions to the collision strengths for these allowed lines have converged with "top-up" from the Burgess-Tully sum rule incorporated. Comparisons are made with the results of Bautista et al. and significant differences are found for some of the optically allowed lines considered.

Collisions Strengths and Effective Collision Strengths for transitions within the ground state configuration of S III

We have carried out a 29-state R-matrix calculation in order to calculate collision strengths and effective collision strengths for the electron impact excitation of S III. The recently developed parallel RMATRX II suite of codes have been used, which perform the calculation in intermediate coupling. Collision strengths have been generated over an electron energy range of 0-12 Ryd, and effective collision strength data have been calculated from these at electron temperatures in the range 1000-100,000 K. Results are here presented for the fine-structure transitions between the ground-state configurations of 3s(2)3p(2) P-3(0,1,2), D-1(2), and S-1(0), and the values given resolve a discrepancy between two previous R-matrix calculations.

Electron-impact excitation of Sc II: collision strengths and effective collision strengths for fine-structure transitions.

Accurate fine-structure atomic data for the Fe-peak elements are essential for interpreting astronomical spectra. There is a severe paucity of data available for Sc II, highlighted by the fact that no collision strengths are readily available for this ion. We present electron-impact excitation collision strengths and Maxwellian averaged effective collision strengths for Sc II. The collision strengths were calculated for all 3916 transitions amongst 89 jj levels (arising from the 3d4s, 3d2, 4s2, 3d4p, 4s4p, 3d5s, 3d4d, 3d5p, 4p2 and 3d4f configurations), resulting in a 944 coupled channel problem. The R-matrix package RMATRXII was utilized, along with the transformation code FINE and the external region code PSTGF, to calculate the collision strengths for a range of incident electron energies in the 0 to 8.3 Rydberg region. Maxwellian averaged effective collision strengths were then produced for 27 temperatures lying within the astrophysically significant range of 30 to 105 K.
The collision strengths and effective collision strengths were produced for two different target models. The purpose was to systematically examine the effect of including open 3p correlation terms into the configuration interaction expansion for the wavefunction. The first model consisted of all 36 CI terms that could be generated with the 3p core closed. The second model incorporated an additional six configurations which allowed for single-electron excitations from within the 3p core. Comparisons are made between the two models and the results of Bautista et al., obtained by private communication. It is concluded that the first model produced the most reliable set of collision and effective collision strengths for use in astrophysical and plasma applications.