Source author record

T. J. Millar

T. J. Millar appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2025arXiv

ATOMIUM: Continuum emission and evidence of dust enhancement from binary motion

Low- and intermediate-mass stars on the asymptotic giant branch (AGB) account for a significant portion of the dust and chemical enrichment in their host galaxy. Here we present ALMA observations of the continuum emission at 1.24 mm around a sample of 17 stars from the ATOMIUM survey. From our analysis of the stellar contributions to the continuum flux, we find that the semi-regular variables all have smaller physical radii and fainter monochromatic luminosities than the Mira variables. Comparing these properties with pulsation periods, we find a positive trend between stellar radius and period only for the Mira variables with periods above 300 days and a positive trend between the period and the monochromatic luminosity only for the red supergiants and the most extreme AGB stars with periods above 500 days. We find that the continuum emission at 1.2 mm can be classified into four groups. "Featureless" continuum emission is confined to the (unresolved) regions close to the star for five stars in our sample, relatively uniform extended flux is seen for four stars, tentative elongated features are seen for three stars, and the remaining five stars have unique or unusual morphological features in their continuum maps. These features can be explained by binary companions to 10 out of the 14 AGB stars in our sample. Based on our results we conclude that there are two modes of dust formation: well established pulsation-enhanced dust formation and our newly proposed companion-enhanced dust formation. If the companion is located close to the AGB star, in the wind acceleration region, then additional dust formed in the wake of the companion can increase the mass lost through the dust driven wind. This explains the different dust morphologies seen around our stars and partly accounts for a large scatter in literature mass-loss rates, especially among semiregular stars with small pulsation periods.

preprint2021arXiv

ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars: Motivation, sample, calibration, and initial results

This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. 17 oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (0.025-0.050 arcsec), medium (0.13-0.24 arcsec), and low (about 1 arcsec) angular resolution. The sensitivity per 1.3 km/s channel was 1.5-5 mJy/beam. 13 molecules were designated as primary molecules in the survey: CO, SiO, AlO, AlOH, TiO, TiO2, HCN, SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described; and we highlight one scientific result - the wind kinematics of the ATOMIUM sources. The ATOMIUM sources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion, and the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In 15 sources, some molecular transitions other than 12CO v=0 J=2-1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. Binary interaction with a (sub)stellar companion might (partly) explain the non-monotonic behaviour of the projected velocity field.

preprint2020arXiv

The Role of Ultraviolet Photons in Circumstellar Astrochemistry

Stars with masses between 1 and 8 solar masses (M$_\odot$) lose large amounts of material in the form of gas and dust in the late stages of stellar evolution, during their Asymptotic Giant Branch phase. Such stars supply up to 35% of the dust in the interstellar medium and thus contribute to the material out of which our solar system formed. In addition, the circumstellar envelopes of these stars are sites of complex, organic chemistry with over 80 molecules detected in them. We show that internal ultraviolet photons, either emitted by the star itself or from a close-in, orbiting companion, can significantly alter the chemistry that occurs in the envelopes particularly if the envelope is clumpy in nature. At least for the cases explored here, we find that the presence of a stellar companion, such as a white dwarf star, the high flux of UV photons destroys H$_2$O in the inner regions of carbon-rich AGB stars to levels below those observed and produces species such as C$^+$ deep in the envelope in contrast to the expectations of traditional descriptions of circumstellar chemistry.

preprint2016arXiv

ALMA Observations of a Gap and a Ring in the Protoplanetary Disk around TW Hya

We report the first detection of a gap and a ring in 336 GHz dust continuum emission from the protoplanetary disk around TW Hya, using the Atacama Large Millimeter/Submillimeter Array (ALMA). The gap and ring are located at around 25 and 41 au from the central star, respectively, and are associated with the CO snow line at ~30 au. The gap has a radial width of less than 15 au and a mass deficit of more than 23%, taking into account that the observations are limited to an angular resolution of ~15 au. In addition, the 13CO and C18O J = 3 - 2 lines show a decrement in CO line emission throughout the disk, down to ~10 au, indicating a freeze-out of gas-phase CO onto grain surfaces and possible subsequent surface reactions to form larger molecules. The observed gap could be caused by gravitational interaction between the disk gas and a planet with a mass less than super-Neptune (2M_{Neptune}), or could be the result of the destruction of large dust aggregates due to the sintering of CO ice.

preprint2016arXiv

ALMA-resolved salt emission traces the chemical footprint and inner wind morphology of VY CMa

(abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative.

preprint2016arXiv

Candidate Water Vapor Lines to Locate the H$_{2}$O Snowline through High-Dispersion Spectroscopic Observations I. The Case of a T Tauri Star

Inside the H$_{2}$O snowline of protoplanetary disks, water evaporates from the dust-grain surface into the gas phase, whereas it is frozen out on to the dust in the cold region beyond the snowline. H$_{2}$O ice enhances the solid material in the cold outer part of a disk, which promotes the formation of gas-giant planet cores. We can regard the H$_{2}$O snowline as the surface that divides the regions between rocky and gaseous giant planet formation. Thus observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In this paper, we find candidate water lines to locate the H$_{2}$O snowline through future high-dispersion spectroscopic observations. First, we calculate the chemical composition of the disk and investigate the abundance distributions of H$_{2}$O gas and ice, and the position of the H$_{2}$O snowline. We confirm that the abundance of H$_{2}$O gas is high not only in the hot midplane region inside the H$_{2}$O snowline but also in the hot surface layer of the outer disk. Second, we calculate the H$_{2}$O line profiles and identify those H$_{2}$O lines which are promising for locating the H$_{2}$O snowline: the identified lines are those which have small Einstein $A$ coefficients and high upper state energies. The wavelengths of the candidate H$_{2}$O lines range from mid-infrared to sub-millimeter, and they overlap with the regions accessible to ALMA and future mid-infrared high dispersion spectrographs (e.g., TMT/MICHI, SPICA).

preprint2016arXiv

Chemistry in AGB stars: successes and challenges

Emission and absorption line observations of molecules in late-type stars are a vital component in our understanding of stellar evolution, dust formation and mass loss in these objects. The molecular composition of the gas in the circumstellar envelopes of AGB stars reflects chemical processes in gas whose properties are strong functions of radius with density and temperature varying by more than ten and two orders of magnitude, respectively. In addition, the interstellar UV field plays a critical role in determining not only molecular abundances but also their radial distributions. In this article, I shall briefly review some recent successful approaches to describing chemistry in both the inner and outer envelopes and outline areas of challenge for the future.

preprint2014arXiv

ALMA hints at the presence of two companions in the disk around HD 100546

HD 100546 is a well-studied Herbig Be star-disk system that likely hosts a close-in companion with compelling observational evidence for an embedded protoplanet at 68 AU. We present ALMA observations of the HD 100546 disk which resolve the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390+/-20 AU in radius) whereas the continuum emission is more compact (230+/-20 AU in radius) suggesting radial drift of the mm-sized grains. The CO emission is similar in extent to scattered light images indicating well-mixed gas and um-sized grains in the disk atmosphere. Assuming azimuthal symmetry, a single-component power-law model cannot reproduce the continuum visibilities. The visibilities and images are better reproduced by a double-component model: a compact ring with a width of 21 AU centered at 26 AU and an outer ring with a width of 75+/-3 AU centered at 190+/-3 AU. The influence of a companion and protoplanet on the dust evolution is investigated. The companion at 10 AU facilitates the accumulation of mm-sized grains within a compact ring, ~ 20 - 30 AU, by ~ 10 Myr. The injection of a protoplanet at 1 Myr hastens the ring formation (~ 1.2 Myr) and also triggers the development of an outer ring (~ 100 - 200 AU). These observations provide additional evidence for the presence of a close-in companion and hint at dynamical clearing by a protoplanet in the outer disk.

preprint2014arXiv

Complex organic molecules along the accretion flow in isolated and externally irradiated protoplanetary disks

(Abridged) The birth environment of the Sun will have influenced the conditions in the pre-solar nebula, including the attainable chemical complexity, important for prebiotic chemistry. The formation and distribution of complex organic molecules (COMs) in a disk around a T Tauri star is investigated for two scenarios: (i) an isolated disk, and (ii) a disk irradiated externally by a nearby massive star. The chemistry is calculated along the accretion flow from the outer disk inwards using a comprehensive network. Two simulations are performed, one beginning with complex ices and one with simple ices only. For the isolated disk, COMs are transported without major alteration into the inner disk where they thermally desorb into the gas reaching an abundance representative of the initial assumed ice abundance. For simple ices, COMs efficiently form on grain surfaces under the conditions in the outer disk. Gas-phase COMs are released into the molecular layer via photodesorption. For the irradiated disk, complex ices are also transported inwards; however, they undergo thermal processing caused by the warmer conditions in the irradiated disk which tends to reduce their abundance along the accretion flow. For simple ices, grain-surface chemistry cannot synthesise COMs in the outer disk because the necessary grain-surface radicals, which tend to be particularly volatile, are not sufficiently abundant on the grain surfaces. Gas-phase COMs are formed in the inner region of the irradiated disk via gas-phase chemistry induced by the desorption of strongly bound molecules such as methanol; hence, the abundances are not representative of the initial molecular abundances injected into the outer disk. These results suggest that the composition of comets formed in isolated disks may differ from those formed in externally irradiated disks with the latter composed of more simple ices.

preprint2014arXiv

Complex organic molecules in protoplanetary disks

(Abridged) Protoplanetary disks are vital objects in star and planet formation, possessing all the material which may form a planetary system orbiting the new star. We investigate the synthesis of complex organic molecules (COMs) in disks to constrain the achievable chemical complexity and predict species and transitions which may be observable with ALMA. We have coupled a 2D model of a protoplanetary disk around a T Tauri star with a gas-grain chemical network including COMs. We compare compare synthesised line intensities and calculated column densities with observations and determine those COMs which may be observable in future. COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances 1e-6 - 1e-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, 1e-12 - 1e-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 several T Tauri star-disk systems. 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. 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.

preprint2014arXiv

Photodissociation and chemistry of N$_2$ in the circumstellar envelope of carbon-rich AGB stars

The envelopes of AGB stars are irradiated externally by ultraviolet photons; hence, the chemistry is sensitive to the photodissociation of N$_2$ and CO, which are major reservoirs of nitrogen and carbon, respectively. The photodissociation of N$_2$ has recently been quantified by laboratory and theoretical studies. Improvements have also been made for CO photodissociation. For the first time, we use accurate N$_2$ and CO photodissociation rates and shielding functions in a model of the circumstellar envelope of the carbon-rich AGB star, IRC +10216. We use a state-of-the-art chemical model of an AGB envelope, the latest CO and N$_2$ photodissociation data, and a new method for implementing molecular shielding functions in full spherical geometry with isotropic incident radiation. We compare computed column densities and radial distributions of molecules with observations. The transition of N$_2$ $\to$ N (also, CO $\to$ C $\to$ C$^+$) is shifted towards the outer envelope relative to previous models. This leads to different column densities and radial distributions of N-bearing species, especially those species whose formation/destruction processes largely depend on the availability of atomic or molecular nitrogen, for example, C$_n$N ($n$=1, 3, 5), C$_n$N$^-$ ($n$=1, 3, 5), HC$_n$N ($n$=1, 3, 5, 7, 9), H$_2$CN and CH$_2$CN. The chemistry of many species is directly or indirectly affected by the photodissociation of N$_2$ and CO, especially in the outer shell of AGB stars where photodissociation is important. Thus, it is important to include N$_2$ and CO shielding in astrochemical models of AGB envelopes and other irradiated environments. In general, while differences remain between our model of IRC +10216 and the observed molecular column densities, better agreement is found between the calculated and observed radii of peak abundance.

preprint2014arXiv

The irradiation of water ice by C$^+$ ions in the cosmic environment

We present a first principles molecular dynamics (FPMD) study of the interaction of low energy, positively charged, carbon (C+) projectiles with amorphous solid water clusters at 30 K. Reactions involving the carbon ion at an initial energy of 11 eV and 1.7 eV with 30-molecule clusters have been investigated. Simulations indicate that the neutral isoformyl radical, COH, and carbon monoxide, CO, are the dominant products of these reactions. All these reactions are accompanied by the transfer of a proton from the reacting water molecule to the ice, where it forms a hydronium ion. We find that COH is formed either via a direct, "knock-out", mechanism following the impact of the C+ projectile upon a water molecule or by creation of a COH_2^+ intermediate. The direct mechanism is more prominent at higher energies. CO is generally produced following the dissociation of COH. More frequent production of the formyl radical, HCO, is observed here than in gas phase calculations. A less commonly occurring product is the dihydroxymethyl, CH(OH)_2, radical. Although a minor result, its existence gives an indication of the increasing chemical complexity which is possible in such heterogeneous environments.

preprint2013arXiv

Molecular line emission from a protoplanetary disk irradiated externally by a nearby massive star

Star formation often occurs within or nearby stellar clusters. Irradiation by nearby massive stars can photoevaporate protoplanetary disks around young stars (so-called proplyds) which raises questions regarding the ability of planet formation to take place in these environments. We investigate the two-dimensional physical and chemical structure of a protoplanetary disk surrounding a low-mass (T Tauri) star which is irradiated by a nearby massive O-type star to determine the survivability and observability of molecules in proplyds. Compared with an isolated star-disk system, the gas temperature ranges from a factor of a few (in the disk midplane) to around two orders of magnitude (in the disk surface) higher in the irradiated disk. Although the UV flux in the outer disk, in particular, is several orders of magnitude higher, the surface density of the disk is sufficient for effective shielding of the disk midplane so that the disk remains predominantly molecular in nature. We also find that non-volatile molecules, such as HCN and H2O, are able to freeze out onto dust grains in the disk midplane so that the formation of icy planetesimals, e.g., comets, may also be possible in proplyds. We have calculated the molecular line emission from the disk assuming LTE and determined that multiple transitions of atomic carbon, CO (and isotopologues, 13CO and C18O), HCO+, CN, and HCN may be observable with ALMA, allowing characterization of the gas column density, temperature, and optical depth in proplyds at the distance of Orion (~400 pc).

preprint2013arXiv

Photodetachment as destruction mechanism for CN- and C3N- anions in circumstellar envelopes

Absolute photodetachment cross sections of two anions of astrophysical importance CN- and C3N- were measured to be (1.18 +- (0.03)_stat (0.17)_sys) * 10^-17 cm^2 and (1.43 +- (0.14)_stat (0.37)_sys) * 10^-17 cm^2 respectively at the ultraviolet wavelength of 266 nm (4.66 eV). These relatively large values of the cross sections imply that photodetachment can play a major role in the destruction mechanisms of these anions particularly in photon-dominated regions. We have therefore carried out model calculations using the newly measured cross sections to investigate the abundance of these molecular anions in the cirumstellar envelope of the carbon-rich star IRC+10216. The model predicts the relative importance of the various mechanisms of formation and destruction of these species in different regions of the envelope. UV photodetachment was found to be the major destruction mechanism for both CN- and C3N- anions in those regions of the envelope, where they occur in peak abundance. It was also found that photodetachment plays a crucial role in the degradation of these anions throughout the circumstellar envelope.

preprint2012arXiv

Chemical Processes in Protoplanetary Disks II. On the importance of photochemistry and X-ray ionization

We investigate the impact of photochemistry and X-ray ionization on the molecular composition of, and ionization fraction in, a protoplanetary disk surrounding a typical T Tauri star. We use a sophisticated physical model, which includes a robust treatment of the radiative transfer of UV and X-ray radiation, and calculate the time-dependent chemical structure using a comprehensive chemical network. In previous work, we approximated the photochemistry and X-ray ionization, here, we recalculate the photoreaction rates using the explicit UV wavelength spectrum and wavelength-dependent reaction cross sections. We recalculate the X-ray ionization rate using our explicit elemental composition and X-ray energy spectrum. We find photochemistry has a larger influence on the molecular composition than X-ray ionization. Observable molecules sensitive to the photorates include OH, HCO+, N2H+, H2O, CO2 and CH3OH. The only molecule significantly affected by the X-ray ionization is N2H+ indicating it is safe to adopt existing approximations of the X-ray ionization rate in typical T Tauri star-disk systems. The recalculation of the photorates increases the abundances of neutral molecules in the outer disk, highlighting the importance of taking into account the shape of the UV spectrum in protoplanetary disks. A recalculation of the photoreaction rates also affects the gas-phase chemistry due to the adjustment of the H/H2 and C+/C ratios. The disk ionization fraction is not significantly affected by the methods adopted to calculate the photochemistry and X-ray ionization. We determine there is a probable 'dead zone' where accretion is suppressed, present in a layer, Z/R <~ 0.1 - 0.2, in the disk midplane, within R \approx 200 AU.

preprint2012arXiv

The UMIST database for astrochemistry 2012

We present the fifth release of the UMIST Database for Astrochemistry (UDfA). The new reaction network contains 6173 gas-phase reactions, involving 467 species, 47 of which are new to this release. We have updated rate coefficients across all reaction types. We have included 1171 new anion reactions and updated and reviewed all photorates. In addition to the usual reaction network, we also now include, for download, state-specific deuterated rate coefficients, deuterium exchange reactions and a list of surface binding energies for many neutral species. Where possible, we have referenced the original source of all new and existing data. We have tested the main reaction network using a dark cloud model and a carbon-rich circumstellar envelope model. We present and briefly discuss the results of these models.

preprint2011arXiv

Carbon chemistry in Galactic Bulge Planetary Nebulae

Galactic Bulge Planetary Nebulae show evidence of mixed chemistry with emission from both silicate dust and PAHs. This mixed chemistry is unlikely to be related to carbon dredge up, as third dredge-up is not expected to occur in the low mass Bulge stars. We show that the phenomenon is widespread, and is seen in 30 nebulae out of 40 of our sample, selected on the basis of their infrared flux. HST images and UVES spectra show that the mixed chemistry is not related to the presence of emission-line stars, as it is in the Galactic disk population. We also rule out interaction with the ISM as origin of the PAHs. Instead, a strong correlation is found with morphology, and the presence of a dense torus. A chemical model is presented which shows that hydrocarbon chains can form within oxygen-rich gas through gas-phase chemical reactions. The model predicts two layers, one at $A_V\sim 1.5$ where small hydrocarbons form from reactions with C$^+$, and one at $A_V\sim 4$, where larger chains (and by implication, PAHs) form from reactions with neutral, atomic carbon. These reactions take place in a mini-PDR. We conclude that the mixed chemistry phenomenon occurring in the Galactic Bulge Planetary Nebulae is best explained through hydrocarbon chemistry in an UV-irradiated, dense torus.

preprint2011arXiv

Discovery of interstellar anions in Cepheus star-forming region

We report the detection of microwave emission lines from the hydrocarbon anion C6H- and its parent neutral C6H in the star-forming region L1251A (in Cepheus), and the pre-stellar core L1512 (in Auriga). The carbon-chain-bearing species C4H, HC3N, HC5N, HC7N and C3S are also detected in large abundances. The observations of L1251A constitute the first detections of anions and long-chain polyynes and cyanopolyynes (with more than 5 carbon atoms) in the Cepheus Flare star-forming region, and the first detection of anions in the vicinity of a protostar outside of the Taurus molecular cloud complex, highlighting a wider importance for anions in the chemistry of star formation. Rotational excitation temperatures have been derived from the HC3N hyperfine structure lines, and are found to be 6.2 K for L1251A and 8.7 K for L1512. The anion-to-neutral ratios are 3.6% and 4.1%, respectively, which are within the range of values previously observed in the interstellar medium, and suggest a relative uniformity in the processes governing anion abundances in different dense interstellar clouds. This research contributes towards the growing body of evidence that carbon chain anions are relatively abundant in interstellar clouds throughout the Galaxy, but especially in the regions of relatively high density and high depletion surrounding pre-stellar cores and young, embedded protostars.

preprint2010arXiv

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 (sub-milli-arcsecond 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 freeze out 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 homogenise 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 >~ 1.0E-11).

preprint2010arXiv

Complex Molecule Formation in Grain Mantles

Context: Complex molecules such as ethanol and dimethyl ether have been observed in a number of hot molecular cores and hot corinos. Attempts to model the molecular formation process using gas phase only models have so far been unsuccessful. Aims : To demonstrate that grain surface processing is a viable mechanism for complex molecule formation in these environments. Methods: A variable environment parameter computer model has been constructed which includes both gas and surface chemistry. This is used to investigate a variety of cloud collapse scenarios. Results: Comparison between model results and observation shows that by combining grain surface processing with gas phase chemistry complex molecules can be produced in observed abundances in a number of core and corino scenarios. Differences in abundances are due to the initial atomic and molecular composition of the core/corino and varying collapse timescales. Conclusions: Grain surface processing, combined with variation of physical conditions, can be regarded as a viable method for the formation of complex molecules in the environment found in the vicinity of a hot core/corino and produce abundances comparable to those observed.

preprint2009arXiv

The abundance of HNCO and its use as a diagnostic of environment

We aim to investigate the chemistry and gas phase abundance of HNCO and the variation of the HNCO/CS abundance ratio as a diagnostic of the physics and chemistry in regions of massive star formation. A numerical-chemical model has been developed which self-consistently follows the chemical evolution of a hot core. The model comprises of two distinct stages. An initial collapse phase is immediately followed by an increase in temperature which represents the switch on of a central massive star and the subsequent evolution of the chemistry in a hot, dense gas cloud (the hot core). During the collapse phase, gas species are allowed to accrete on to grain surfaces where they can participate in further reactions. During the hot core phase surface species thermally desorb back in to the ambient gas and further chemical evolution takes place. For comparison, the chemical network was also used to model a simple dark cloud and photodissociation regions. Our investigation reveals that HNCO is inefficiently formed when only gas-phase formation pathways are considered in the chemical network with reaction rates consistent with existing laboratory data. Using currently measured gas phase reaction rates, obtaining the observed HNCO abundances requires its formation on grain surfaces. However our model shows that the gas phase HNCO in hot cores is not a simple direct product of the evaporation of grain mantles. We also show that the HNCO/CS abundance ratio varies as a function of time in hot cores and can match the range of values observed. This ratio is not unambiguously related to the ambient UV field as been suggested - our results are inconsistent with the hypothesis of Martin et al (2008). In addition, our results show that this ratio is extremely sensitive to the initial sulphur abundance.

preprint2007arXiv

Molecular Hydrogen Emission from Protoplanetary Disks II. Effects of X-ray Irradiation and Dust Evolution

Detailed models for the density and temperature profiles of gas and dust in protoplanetary disks are constructed by taking into account X-ray and ultraviolet (UV) irradiation from a central T Tauri star, as well as dust size growth and settling toward the disk midplane. The spatial and size distributions of dust grains in the disks are numerically computed by solving the coagulation equation for settling dust particles. The level populations and line emission of molecular hydrogen are calculated using the derived physical structure of the disks. X-ray irradiation is the dominant heating source of the gas in the inner disk region and in the surface layer, while the far UV heating dominates otherwise. If the central star has strong X-ray and weak UV radiation, the H2 level populations are controlled by X-ray pumping, and the X-ray induced transition lines could be observable. If the UV irradiation is strong, the level populations are controlled by thermal collisions or UV pumping, depending on the properties of the dust grains in the disks. As the dust particles evolve in the disks, the gas temperature at the disk surface drops because the grain photoelectric heating becomes less efficient, while the UV radiation fields become stronger due to the decrease of grain opacity. This makes the H2 level populations change from local thermodynamic equilibrium (LTE) to non-LTE distributions, which results in changes to the line ratios of H2 emission. Our results suggest that dust evolution in protoplanetary disks could be observable through the H2 line ratios. The emission lines are strong from disks irradiated by strong UV and X-rays and possessing small dust grains; such disks will be good targets in which to observe H2 emission.