Source author record

K. Willacy

K. Willacy 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

8works
4topics
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

8 published item(s)

preprint2016arXiv

Chemistry as a diagnostic of prestellar core geometry

We present a new method for assessing the intrinsic 3D shape of prestellar cores from molecular column densities. We have employed hydrodynamic simulations of contracting, isothermal cores considering three intrinsic geometries: spherical, cylindrical/filamentary and disk-like. We have coupled our hydrodynamic simulations with non-equilibrium chemistry. We find that a) when cores are observed very elongated (i.e. for aspect ratios $\le$ 0.15) the intrinsic 3D geometry can be probed by their 2D molecular emission maps, since these exhibit significant qualitative morphological differences between cylindrical and disk-like cores. Specifically, if a disk-like core is observed as a filamentary object in dust emission, then it will be observed as two parallel filaments in $\rm{N_2H^{+}}$; b) for cores with higher aspect ratios (i.e. 0.15 $\sim$ 0.9) we define a metric $Δ$ that quantifies whether a molecular column density profile is centrally peaked, depressed or flat. We have identified one molecule ($\rm{CN}$) for which $Δ$ as a function of the aspect ratio probes the 3D geometry of the core; and c) for cores with almost circular projections (i.e. for aspect ratios $\sim$ 1), we have identified three molecules ($\rm{OH}$, $\rm{CO}$ and $\rm{H_2CO}$) that can be used to probe the intrinsic 3D shape by close inspection of their molecular column density radial profiles. We alter the temperature and the cosmic-ray ionization rate and demonstrate that our method is robust against the choice of parameters.

preprint2015arXiv

The composition of the protosolar disk and the formation conditions for comets

Conditions in the protosolar nebula have left their mark in the composition of cometary volatiles, thought to be some of the most pristine material in the solar system. Cometary compositions represent the end point of processing that began in the parent molecular cloud core and continued through the collapse of that core to form the protosun and the solar nebula, and finally during the evolution of the solar nebula itself as the cometary bodies were accreting. Disentangling the effects of the various epochs on the final composition of a comet is complicated. But comets are not the only source of information about the solar nebula. Protostellar disks around young stars similar to the protosun provide a way of investigating the evolution of disks similar to the solar nebula while they are in the process of evolving to form their own solar systems. In this way we can learn about the physical and chemical conditions under which comets formed, and about the types of dynamical processing that shaped the solar system we see today. This paper summarizes some recent contributions to our understanding of both cometary volatiles and the composition, structure and evolution of protostellar disks.

preprint2014arXiv

A Herschel [C II] Galactic plane survey II: CO-dark H2 in clouds

ABRIDGED: Context: HI and CO large scale surveys of the Milky Way trace the diffuse atomic clouds and the dense shielded regions of molecular hydrogen clouds. However, until recently, we have not had spectrally resolved C+ surveys to characterize the photon dominated interstellar medium, including, the H2 gas without C, the CO-dark H2, in a large sample of clouds. Aims: To use a sparse Galactic plane survey of the 1.9 THz [C II] spectral line from the Herschel Open Time Key Programme, Galactic Observations of Terahertz C+ (GOT C+), to characterize the H2 gas without CO in a statistically significant sample of clouds. Methods: We identify individual clouds in the inner Galaxy by fitting [CII] and CO isotopologue spectra along each line of sight. We combine these with HI spectra, along with excitation models and cloud models of C+, to determine the column densities and fractional mass of CO-dark H2 clouds. Results: We identify 1804 narrow velocity [CII] interstellar cloud components in different categories. About 840 are diffuse molecular clouds with no CO, 510 are transition clouds containing [CII] and 12CO, but no 13CO, and the remainder are dense molecular clouds containing 13CO emission. The CO-dark H2 clouds are concentrated between Galactic radii 3.5 to 7.5 kpc and the column density of the CO-dark H2 layer varies significantly from cloud-to-cloud with an average 9X10^(20) cm-2. These clouds contain a significant fraction of CO-dark H2 mass, varying from ~75% for diffuse molecular clouds to ~20% for dense molecular clouds. Conclusions: We find a significant fraction of the warm molecular ISM gas is invisible in HI and CO, but is detected in [CII]. The fraction of CO-dark H2 is greatest in the diffuse clouds and decreases with increasing total column density, and is lowest in the massive clouds.

preprint2014arXiv

Effect of OH depletion on measurements of the mass-to-flux ratio in molecular cloud cores

The ratio of mass and magnetic flux determines the relative importance of magnetic and gravitational forces in the evolution of molecular clouds and their cores. Its measurement is thus central in discriminating between different theories of core formation and evolution. Here we discuss the effect of chemical depletion on measurements of the mass-to-flux ratio using the same molecule (OH) both for Zeeman measurements of the magnetic field and the determination of the mass of the region. The uncertainties entering through the OH abundance in determining separately the magnetic field and the mass of a region have been recognized in the literature. It has been proposed however that, when comparing two regions of the same cloud, the abundance will in both cases be the same. We show that this assumption is invalid. We demonstrate that when comparing regions with different densities, the effect of OH depletion in measuring changes of the mass-to-flux ratio between different parts of the same cloud can even reverse the direction of the underlying trends (for example, the mass-to-flux ratio may appear to decrease as we move to higher density regions). The systematic errors enter primarily through the inadequate estimation of the mass of the region.

preprint2011arXiv

The Influence of Deuteration and Turbulent Diffusion on the Observed D/H Ratio

The influence of turbulent mixing on the chemistry of the interstellar medium has so far received little attention. Previous studies of this effect have suggested that it might play an important role in mixing the various phases of the interstellar medium. In this paper we examine the potential effects of turbulent diffusion on the deuterium chemistry within molecular clouds. We find that such mixing acts to reduce the efficiency of deuteration in these clouds by increasing the ionization fraction and reducing freeze-out of heavy molecules. This leads to lower abundances for many deuterated species. We also examine the influence of turbulent mixing on the transition from atomic hydrogen to H2 and from atomic deuterium to HD near the cloud edge. We find that including turbulent diffusion in our models serves to push these transitions deeper into the cloud and helps maintain a higher atomic fraction throughout the cloud envelope. Based on these findings, we propose a new process to account for the significant scatter in the observed atomic D/H ratio for galactic sightlines extending beyond the Local Bubble. Although several mechanisms have been put forward to explain this scatter, they are unable to fully account for the range in D/H values. We suggest a scenario in which turbulent mixing of atomic and molecular gas at the edges of molecular clouds causes the observed atomic D/H ratio to vary by a factor of ~2.

preprint2010arXiv

Circumstellar molecular composition of the oxygen-rich AGB star IK~Tau: II. In-depth non-LTE chemical abundance analysis

Aims: Little information exists on the circumstellar molecular abundance stratifications of many molecules. The aim is to study the circumstellar chemical abundance pattern of 11 molecules and isotopologs ($^{12}$CO, $^{13}$CO, SiS, $^{28}$SiO, $^{29}$SiO, $^{30}$SiO, HCN, CN, CS, SO, SO$_2$) in the oxygen-rich evolved star IK~Tau. Methods: We have performed an in-depth analysis of a large number of molecular emission lines excited in the circumstellar envelope around IK~Tau. The analysis is done based on a non-local thermodynamic equilibrium (non-LTE) radiative transfer analysis, which calculates the temperature and velocity structure in a self-consistent way. The chemical abundance pattern is coupled to theoretical outer wind model predictions including photodestruction and cosmic ray ionization. Not only the integrated line intensities, but also the line shapes, are used as diagnostic tool to study the envelope structure. Results: The deduced wind acceleration is much slower than predicted from classical theories. SiO and SiS are depleted in the envelope, possibly due to the adsorption onto dust grains. For HCN and CS a clear difference with respect to inner wind non-equilibrium predictions is found, either indicating uncertainties in the inner wind theoretical modeling or the possibility that HCN and CS (or the radical CN) participate in the dust formation. The low signal-to-noise profiles of SO and CN prohibit an accurate abundance determination; the modeling of high-excitation SO$_2$ lines is cumbersome, possibly related to line misidentifications or problems with the collisional rates. The SiO isotopic ratios ($^{29}$SiO/$^{28}$SiO and $^{30}$SiO/$^{28}$SiO) point toward an enhancement in $^{28}$SiO compared to results of classical stellar evolution codes. Predictions for H$_2$O lines in the spectral range of the Herschel/HIFI mission are performed. [abbreviated]

preprint2009arXiv

Deuterium chemistry in protoplanetary disks II The inner 30 AU

We present the results of models of the chemistry, including deuterium, in the inner regions of protostellar disks. We find good agreement with recent gas phase observations of several (non--deuterated) species. We also compare our results with observations of comets and find that in the absence of other processing e.g. in the accretion shock at the surface of the disk, or by mixing in the disk, the calculated D/H ratios in ices are higher than measured and reflect the D/H ratio set in the molecular cloud phase. Our models give quite different abundances and molecular distributions to other inner disk models because of the differences in physical conditions in the model disk. This emphasizes how changes in the assumptions about the density and temperature distribution can radically affect the results of chemical models.

preprint2007arXiv

The chemistry of multiply deuterated molecules in protoplanetary disks. I. The outer disk

We present new models of the deuterium chemistry in protoplanetary disks, including, for the first time, multiply deuterated species. We use these models to explore whether observations in combination with models can give us clues as to which desorption processes occur in disks. We find, in common with other authors, that photodesorption can allow strongly bound molecules such as HDO to exist in the gas phase in a layer above the midplane. Models including this process give the best agreement with the observations. In the midplane, cosmic ray heating can desorb weakly bound molecules such as CO and N$_2$. We find the observations suggest that N$_2$ is gaseous in this region, but that CO must be retained on the grains to account for the observed DCO$^+$/HCO$^+$. This could be achieved by CO having a higher binding energy than N$_2$ (as may be the case when these molecules are accreted onto water ice) or by a smaller cosmic ray desorption rate for CO than assumed here, as suggested by recent theoretical work. For gaseous molecules the calculated deuteration can be greatly changed by chemical processing in the disk from the input molecular cloud values. On the grains singly deuterated species tend to retain the D/H ratio set in the molecular cloud, whereas multiply deuterated species are more affected by the disk chemistry. Consequently the D/H ratios observed in comets may be partly set in the parent cloud and partly in the disk, depending on the molecule.