BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Karen Willacy

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.EPastro-phastro-ph.SRphysics.chem-ph

Catalog footprint

What is connected

14works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

14 published item(s)

preprint2022arXiv

Comets in context: Comparing comet compositions with protosolar nebula models

Comets provide a valuable window into the chemical and physical conditions at the time of their formation in the young solar system. We seek insights into where and when these objects formed by comparing the range of abundances observed for nine molecules and their average values across a sample of 29 comets to the predicted midplane ice abundances from models of the protosolar nebula. Our fiducial model, where ices are inherited from the interstellar medium, can account for the observed mixing ratio ranges of each molecule considered, but no single location or time reproduces the abundances of all molecules simultaneously. This suggests that each comet consists of material processed under a range of conditions. In contrast, a model where the initial composition of disk material is `reset', wiping out any previous chemical history, cannot account for the complete range of abundances observed in comets. Using toy models that combine material processed under different thermal conditions we find that a combination of warm (CO-poor) and cold (CO-rich) material is required to account for both the average properties of the Jupiter-family and Oort cloud comets, and the individual comets we consider. This could occur by the transport (either radial or vertical) of ice-coated dust grains in the early solar system. Comparison of the models to the average Jupiter-family and Oort cloud comet compositions suggest the two families formed in overlapping regions of the disk, in agreement with the findings of A'Hearn et al. (2012) and with the predictions of the Nice model (Gomes et al. 2005, Tsiganis et al. 2005).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Linking atmospheric chemistry of the hot Jupiter HD 209458b to its formation location through infrared transmission and emission spectra

The elemental ratios of carbon, nitrogen, and oxygen in the atmospheres of hot Jupiters may hold clues to their formation locations in the protostellar disc. In this work, we adopt gas phase chemical abundances of C, N and O from several locations in a disc chemical kinetics model as sources for the envelope of the hot Jupiter HD 209458b and evolve the planet's atmospheric composition using a 1D chemical kinetics model, treating both vertical mixing and photochemistry. We consider two atmospheric pressure-temperature profiles, one with and one without a thermal inversion. From each of the resulting 32 atmospheric composition profiles, we find that the molecules CH4, NH3, HCN, and C2H2 are more prominent in the atmospheres computed using a realistic non-inverted P-T profile in comparison to a prior equilibrium chemistry based work which used an analytical P-T profile. We also compute the synthetic transmission and emission spectra for these atmospheres and find that many spectral features vary with the location in the disc where the planet's envelope was accreted. By comparing with the species detected using the latest high-resolution ground-based observations, our model suggests HD 209458b could have accreted most of its gas between the CO2 and CH4 icelines with a super solar C/O ratio from its protostellar disc, which in turn directly inherited its chemical abundances from the protostellar cloud. Finally, we simulate observing the planet with the James Webb Space Telescope (JWST) and show that differences in spectral signatures of key species can be recognized. Our study demonstrates the enormous importance of JWST in providing new insights into hot Jupiter's formation environments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Vertical distribution of cyclopropenylidene and propadiene in the atmosphere of Titan

Titan's atmosphere is a natural laboratory for exploring the photochemical synthesis of organic molecules. Significant recent advances in the study of the atmosphere of Titan include: (a) detection of C$_3$ molecules: C$_3$H$_6$, CH$_2$CCH$_2$, c-C$_3$H$_2$, and (b) retrieval of C$_6$H$_6$, which is formed primarily via C$_3$ chemistry, from Cassini-UVIS data. The detection of $c$-C$_3$H$_2$ is of particular significance since ring molecules are of great astrobiological importance. Using the Caltech/JPL KINETICS code, along with the best available photochemical rate coefficients and parameterized vertical transport, we are able to account for the recent observations. It is significant that ion chemistry, reminiscent of that in the interstellar medium, plays a major role in the production of c-C$_3$H$_2$ above 1000 km.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Non-ideal MHD simulations of subcritical prestellar cores with non-equilibrium chemistry

Non-ideal magnetohydrodynamic (MHD) effects are thought to be gravity's closest ally in overcoming the support of magnetic fields and in forming stars. Here, we modify the publicly available version of the adaptive mesh refinement code FLASH (Fryxell et al. 2000; Dubey et al. 2008) to include a detailed treatment of non-ideal MHD and study such effects in collapsing prestellar cores. We implement two very extended non-equilibrium chemical networks, the largest of which is comprised of $\sim$ 300 species and includes a detailed description of deuterium chemistry. The ambipolar-diffusion, Ohmic and Hall resistivities are then self-consistently calculated from the abundances of charged species. We present a series of 2-dimensional axisymmetric simulations where we vary the chemical model, cosmic-ray ionization rate, and grain distribution. We benchmark our implementation against ideal MHD simulations and previously-published results. We show that, at high densities ($n_{\rm{H_2}}>~10^6~\rm{cm^{-3}}$), the ion that carries most of the perpendicular and parallel conductivities is not $\rm{H_3^+}$ as was previously thought, but is instead $\rm{D_3^+}$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Physical and Chemical Structure of the Disk and Envelope of the Class 0/I Protostar L1527

Sub-millimeter spectral line and continuum emission from the protoplanetary disks and envelopes of protostars are powerful probes of their structure, chemistry, and dynamics. Here we present a benchmark study of our modeling code, RadChemT, that for the first time uses a chemical model to reproduce ALMA C$^{18}$O (2-1) and CARMA $^{12}$CO (1-0) and N$_{2}$H$^{+}$ (1-0) observations of L1527, that allow us to distinguish the disk, the infalling envelope and outflow of this Class 0/I protostar. RadChemT combines dynamics, radiative transfer, gas chemistry and gas-grain reactions to generate models which can be directly compared with observations for individual protostars. Rather than individually fit abundances to a large number of free parameters, we aim to best match the spectral line maps by (i) adopting a physical model based on density structure and luminosity derived primarily from previous work that fit SED and 2D imaging data, updating it to include a narrow jet detected in CARMA and ALMA data near ($\leq 75$au) the protostar, and then (ii) computing the resulting astrochemical abundances for 292 chemical species. Our model reproduces the C$^{18}$O and N$_{2}$H$^{+}$ line strengths within a factor of 3.0; this is encouraging considering the pronounced abundance variation (factor $> 10^3$) between the outflow shell and CO snowline region near the midplane. Further, our modeling confirms suggestions regarding the anti-correlation between N$_{2}$H$^{+}$ and the CO snowline between 400 au to 2,000 au from the central star. Our modeling tools represent a new and powerful capability with which to exploit the richness of spectral line imaging provided by modern submillimeter interferometers.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

An Ammonia Spectral Map of the L1495-B218 Filaments in the Taurus Molecular Cloud: II CCS & HC$_7$N Chemistry and Three Modes of Star Formation in the Filaments

We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both evolved and young regions and weak emission in two evolved regions. HC$_7$N emission is observed only in L1495A-N and L1521D. We find that CCS and HC$_7$N intensity peaks do not coincide with NH$_3$ or dust continuum intensity peaks. We also find that the fractional abundance of CCS does not show a clear correlation with the dynamical evolutionary stage of dense cores. Our findings and chemical modeling indicate that the fractional abundances of CCS and HC$_7$N are sensitive to the initial gas-phase C/O ratio, and they are good tracers of young condensed gas only when the initial C/O is close to solar value. Kinematic analysis using multiple lines including NH$_3$, HC$_7$N, CCS, CO, HCN, \& HCO$^+$ suggests that there may be three different star formation modes in the L1495-B218 filaments. At the hub of the filaments, L1495A/B7N has formed a stellar cluster with large-scale inward flows (fast mode), while L1521D, a core embedded in a filament, is slowly contracting due to its self-gravity (slow mode). There is also one isolated core that appears to be marginally stable and may undergo quasi-static evolution (isolated mode).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Probing Planet Forming Zones with Rare CO Isotopologues

The gas near the midplanes of planet-forming protostellar disks remains largely unprobed by observations due to the high optical depth of commonly observed molecules such as CO and H$_2$O. However, rotational emission lines from rare molecules may have optical depths near unity in the vertical direction, so that the lines are strong enough to be detected, yet remain transparent enough to trace the disk midplane. Here we present a chemical model of an evolving T-Tauri disk and predict the optical depths of rotational transitions of $^{12}$C$^{16}$O, $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O. The MRI-active disk is primarily heated by the central star due to the formation of the dead zone. CO does not freeze out in our modeled region within $70$AU around a sunlike star. However, the abundance of CO decreases because of the formation of complex organic molecules (COM), producing an effect that can be misinterpreted as the "snow line". These results are robust to variations in our assumptions about the evolution of the gas to dust ratio. The optical depths of low-order rotational lines of C$^{17}$O are around unity, making it possible to see into the disk midplane using C$^{17}$O. Combining observations with modeled C$^{17}$O$/$H$_2$ ratios, like those we provide, can yield estimates of protoplanetary disks' gas masses.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A search for co-evolving ion and neutral gas species in prestellar molecular cloud cores

Comparison of linewidths of spectral line profiles of ions and neutral molecules have been recently used to estimate the strength of the magnetic field in turbulent star-forming regions. However, the ion (HCO+) and neutral (HCN) species used in such studies may not be necessarily co-evolving at every scale and density and may thus not trace the same regions. Here, we use coupled chemical/dynamical models of evolving prestellar molecular cloud cores including non-equilibrium chemistry, with and without magnetic fields, to study the spatial distribution of HCO+ and HCN, which have been used in observations of spectral linewidth differences to date. In addition, we seek new ion-neutral pairs that are good candidates for such observations because they have similar evolution and are approximately co-spatial in our models. We identify three such good candidate pairs: HCO+/NO, HCO+/CO, and NO+/NO.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-Equilibrium Chemistry of Dynamically Evolving Prestellar Cores: I. Basic Magnetic and Non-Magnetic Models and Parameter Studies

We combine dynamical and non-equilibrium chemical modeling of evolving prestellar molecular cloud cores, and explore the evolution of molecular abundances in the contracting core. We model both magnetic cores, with varying degrees of initial magnetic support, and non-magnetic cores, with varying collapse delay times. We explore, through a parameter study, the competing effects of various model parameters in the evolving molecular abundances, including the elemental C/O ratio, the temperature, and the cosmic-ray ionization rate. We find that different models show their largest quantitative differences at the center of the core, whereas the outer layers, which evolve slower, have abundances which are severely degenerate among different dynamical models. There is a large range of possible abundance values for different models at a fixed evolutionary stage (central density), which demonstrates the large potential of chemical differentiation in prestellar cores. However, degeneracies among different models, compounded with uncertainties induced by other model parameters, make it difficult to discriminate among dynamical models. To address these difficulties, we identify abundance ratios between particular molecules, the measurement of which would have maximal potential for discrimination among the different models examined here. In particular, we find that the ratios between NH3 and CO; NH2 and CO; NH3 and HCO+ are sensitive to the evolutionary timescale, and that the ratio between HCN and OH is sensitive to the C/O ratio. Finally, we demonstrate that measurements of the central deviation (central depletion or enhancement) of abundances of certain molecules are good indicators of the dynamics of the core.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-Equilibrium Chemistry of Dynamically Evolving Prestellar Cores: II. Ionization and Magnetic Field

We study the effect that non-equilibrium chemistry in dynamical models of collapsing molecular cloud cores has on measurements of the magnetic field in these cores, the degree of ionization, and the mean molecular weight of ions. We find that OH and CN, usually used in Zeeman observations of the line-of-sight magnetic field, have an abundance that decreases toward the center of the core much faster than the density increases. As a result, Zeeman observations tend to sample the outer layers of the core and consistently underestimate the core magnetic field. The degree of ionization follows a complicated dependence on the number density at central densities up to 10^5 cm^{-3} for magnetic models and 10^6 cm^{-3} in non-magnetic models. At higher central densities the scaling approaches a power-law with a slope of -0.6 and a normalization which depends on the cosmic-ray ionization rate ζ and the temperature T as (ζT)^1/2. The mean molecular weight of ions is systematically lower than the usually assumed value of 20 - 30, and, at high densities, approaches a value of 3 due to the asymptotic dominance of the H3+ ion. This significantly lower value implies that ambipolar diffusion operates faster.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Spitzer mapping of molecular hydrogen pure rotational lines in NGC 1333: A detailed study of feedback in star formation

We present mid-infrared spectral maps of the NGC 1333 star forming region, obtained with the the Infrared Spectrometer on board the Spitzer Space Telescope. Eight pure H2 rotational lines, from S (0) to S (7), are detected and mapped. The H2 emission appears to be associated with the warm gas shocked by the multiple outflows present in the region. A comparison between the observed intensities and the predictions of detailed shock models indicates that the emission arises in both slow (12 - 24 km/s) and fast (36 - 53 km/s) C-type shocks with an initial ortho-to-para ratio of ~ 1. The present H2 ortho-to-para ratio exhibits a large degree of spatial variations. In the post-shocked gas, it is usually about 2, i.e. close to the equilibrium value (~ 3). However, around at least two outflows, we observe a region with a much lower (~ 0.5) ortho-to-para ratio. This region probably corresponds to gas which has been heated-up recently by the passage of a shock front, but whose ortho-to-para has not reached equilibrium yet. This, together with the low initial ortho-to-para ratio needed to reproduce the observed emission, provide strong evidence that H2 is mostly in para form in cold molecular clouds. The H2 lines are found to contribute to 25 - 50% of the total outflow luminosity, and thus can be used to ascertain the importance of star formation feedback on the natal cloud. From these lines, we determine the outflow mass loss rate and, indirectly, the stellar infall rate, the outflow momentum and the kinetic energy injected into the cloud over the embedded phase. The latter is found to exceed the binding energy of individual cores, suggesting that outflows could be the main mechanism for core disruption.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Carbon isotope fractionation in protoplanetary disks

We investigate the gas-phase and grain-surface chemistry in the inner 30 AU of a typical protoplanetary disk using a new model which calculates the gas temperature by solving the gas heating and cooling balance and which has an improved treatment of the UV radiation field. We discuss inner-disk chemistry in general, obtaining excellent agreement with recent observations which have probed the material in the inner regions of protoplanetary disks. We also apply our model to study the isotopic fractionation of carbon. Results show that the fractionation ratio, 12C/13C, of the system varies with radius and height in the disk. Different behaviour is seen in the fractionation of different species. We compare our results with 12C/13C ratios in the Solar System comets, and find a stark contrast, indicative of reprocessing.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Benzene formation in the inner regions of protostellar disks

Benzene (c-C6H6) formation in the inner 3 AU of a protostellar disk can be efficient, resulting in high abundances of benzene in the midplane region. The formation mechanism is different to that found in interstellar clouds and in protoplanetary nebulae, and proceeds mainly through the reaction between allene (C3H4) and its ion. This has implications for PAH formation, in that some fraction of PAHs seen in the solar system could be native rather than inherited from the interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Carbon molecules in the inner wind of the oxygen-rich Mira IK Tauri

The gas-phase non-equilibrium chemistry of the inner wind of the oxygen-rich Mira variable IK Tau (NML Tau) is investigated using a physio-chemical model describing the periodically shocked gas close to the stellar photosphere. We predict the formation in large amounts of a few carbon species in the inner envelope, in particular CO$_2$, which has been recently detected with ISO in the spectra of several oxygen-rich semi-regular and Mira stars. The theoretical abundances are also in excellent agreement with values derived from millimeter and sub-millimeter observations, pointing to the fact that some carbon species in oxygen stars do form from shock chemistry in the inner layers and travel the envelope as ``parent'' species.

0 citations0 reviewsNo code linkedOpen access