Source author record

Serena Viti

Serena Viti 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

49works
7topics
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

49 published item(s)

preprint2026arXiv

High-energy astrochemistry in the molecular interstellar medium

In the past decade, there has been a significant shift in astrochemistry with a renewed focus on the role of non-thermal processes on the molecular interstellar medium, in particular energetic particles (such as cosmic ray particles and fast electrons) and X-ray radiation. This has been brought about in large part due to new observations of interstellar complex organic molecules (iCOMS) in environments that would inhibit their formation, such as cold, dense gas in prestellar cores or in the highly energetic environments in galactic centers. In parallel, there has been a plethora of new laboratory investigations on the role of high-energy radiation and electrons on the chemistry of astrophysical ices, demonstrating the ability of this radiation to induce complex chemistry. In recent years, theoretical models have also begun to include newer cosmic-ray-driven processes in both the gas and ice phases. In this review, we unify aspects of the chemistry driven by X-ray radiation and energetic particles into a ``high-energy astrochemistry'', defining this term and reviewing the underlying chemical processes. We conclude by examining various laboratories where high-energy astrochemistry is at play and identify future issues to be tackled.

preprint2024arXiv

The ALCHEMI atlas: principal component analysis reveals starburst evolution in NGC 253

Molecular lines are powerful diagnostics of the physical and chemical properties of the interstellar medium (ISM). These ISM properties, which affect future star formation, are expected to differ in starburst galaxies from those of more quiescent galaxies. We investigate the ISM properties in the central molecular zone of the nearby starburst galaxy NGC 253 using the ultra-wide millimeter spectral scan survey from the ALMA Large Program ALCHEMI. We present an atlas of velocity-integrated images at a 1".6 resolution of 148 unblended transitions from 44 species, including the first extragalactic detection of HCNH$^+$ and the first interferometric images of C$_3$H$^+$, NO, HCS$^+$. We conduct a principal component analysis (PCA) on these images to extract correlated chemical species and to identify key groups of diagnostic transitions. To the best of our knowledge, our dataset is currently the largest astronomical set of molecular lines to which PCA has been applied. The PCA can categorize transitions coming from different physical components in NGC 253 such as i) young starburst tracers characterized by high-excitation transitions of HC$_3$N and complex organic molecules (COMs) versus tracers of on-going star formation (radio recombination lines) and high-excitation transitions of CCH and CN tracing PDRs, ii) tracers of cloud-collision-induced shocks (low-excitation transitions of CH$_3$OH, HNCO, HOCO$^+$, and OCS) versus shocks from star-formation-induced outflows (high-excitation transitions of SiO), as well as iii) outflows showing emission from HOC$^+$, CCH, H$_3$O$^+$, CO isotopologues, HCN, HCO$^+$, CS, and CN. Our findings show these intensities vary with galactic dynamics, star formation activities, and stellar feedback.

preprint2022arXiv

ALMA imaging of the cold molecular and dusty disk in the type 2 active nucleus of the Circinus galaxy

We aim to shed light on the physical properties and kinematics of the molecular material in the nucleus of one of the closest type 2 active galaxies. To this end, we obtained high angular resolution ALMA observations of the nucleus of the Circinus galaxy. The observations map the emission at 350GHz and 690GHz with spatial resolutions of ~3.8pc and ~2.2pc, respectively. The continuum emission traces cold ($T\lesssim100$K) dust in a circumnuclear disk with spiral arms on scales of 25pc, plus a marginally resolved nuclear emission peak. The latter is not extended in polar direction as claimed based on earlier ALMA observations. A significant amount (of the order of 40%) of the 350GHz emission is not related to dust, but most likely free-free emission instead. We detect CO(3-2) and CO(6-5) as well as HCO$^+$(4-3), HCN(4-3), and CS(4-3). The CO emission is extended, showing a spiral pattern, similar to the extended dust emission. Towards the nucleus, CO is excited to higher transitions and its emission is self-absorbed, leading to an apparent hole in the CO(3-2) but not the CO(6-5) emission. On the other hand, the high gas density tracers HCO$^+$, HCN, and CS show a strong, yet unresolved (($\lesssim4$pc) concentration of the emission at the nucleus, pointing at a very small 'torus'. The kinematics are dominated by rotation and point at a geometrically thin disk down to the resolution limit of our observations. In contrast to several other AGNs, no HCN enhancement is found towards the nucleus. The Circinus nucleus is therefore composed of at least two distinct components: (1) an optically thin, warm outflow of ionised gas containing clouds of dust; and (2) a cold molecular and dusty disk. These findings support the most recent radiative transfer calculations of the obscuring structures in AGNs, which find a similar two-component structure. (Abridged)

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

Energizing Star Formation: The Cosmic Ray Ionization Rate in NGC 253 Derived From ALCHEMI Measurements of H$_3$O$^+$ and SO

The cosmic ray ionization rate (CRIR) is a key parameter in understanding the physical and chemical processes in the interstellar medium. Cosmic rays are a significant source of energy in star formation regions, which impacts the physical and chemical processes which drive the formation of stars. Previous studies of the circum-molecular zone (CMZ) of the starburst galaxy NGC 253 have found evidence for a high CRIR value; $10^3-10^6$ times the average cosmic ray ionization rate within the Milky Way. This is a broad constraint and one goal of this study is to determine this value with much higher precision. We exploit ALMA observations towards the central molecular zone of NGC 253 to measure the CRIR. We first demonstrate that the abundance ratio of H$_3$O$^+$ and SO is strongly sensitive to the CRIR. We then combine chemical and radiative transfer models with nested sampling to infer the gas properties and CRIR of several star-forming regions in NGC 253 due to emission from their transitions. We find that each of the four regions modelled has a CRIR in the range $(1-80)\times10^{-14}$ s$^{-1}$ and that this result adequately fits the abundances of other species that are believed to be sensitive to cosmic rays including C$_2$H, HCO$^+$, HOC$^+$, and CO. From shock and PDR/XDR models, we further find that neither UV/X-ray driven nor shock dominated chemistry are a viable single alternative as none of these processes can adequately fit the abundances of all of these species.

preprint2022arXiv

FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed.

preprint2022arXiv

History-independent tracers: Forgetful molecular probes of the physical conditions of the dense interstellar medium

Molecular line emission is a powerful probe of the physical conditions of astrophysical objects but can be complex to model, and it is often unclear which transitions would be the best targets for observers who wish to constrain a given parameter. We therefore produce a list of molecular species for which the gas history can be ignored, removing a major modelling complexity. We then determine the best of these species to observe when attempting to constrain various physical parameters. To achieve this, we use a large set of chemical models with different chemical histories to determine which species have abundances at 1 MYr that are insensitive to the initial conditions. We then use radiative transfer modelling to produce the intensity of every transition of these molecules. We finally compute the mutual information between the physical parameters and all transitions and transition ratios in order to rank their usefulness in determining the value of a given parameter. We find 48 species that are insensitive to the chemical history of the gas, 23 of which have collisional data available. We produce a ranked list of all the transitions and ratios of these species using their mutual information with various gas properties. We show mutual information is an adequate measure of how well a transition can constrain a physical parameter by recovering known probes and demonstrating that random forest regression models become more accurate predictors when high-scoring features are included. Therefore, this list can be used to select target transitions for observations in order to maximize knowledge about those physical parameters.

preprint2022arXiv

Interstellar Hydrocarbons: Degradation Chemistry in Diffuse Clouds

Observations of diffuse clouds showed that they contain a number of simple hydrocarbons (e.g. CH, C$_2$H, (l- and c-)C$_3$H$_2$, and C$_4$H) in abundances that may be difficult to understand on the basis of conventional gas-phase chemical models. Recent experimental results revealed that the photo-decomposition mechanism of hydrogenated amorphous carbon (HAC) and of solid hexane releases a range of hydrocarbons into the gas, containing up to 6 C-atoms for the case of HAC decomposition. These findings motivated us to introduce a new potential input to interstellar chemistry; the "top-down" or degradation scheme, as opposed to the conventional "build-up" or synthesis scheme. In this work, we demonstrate the feasibility of the top-down approach in diffuse clouds using gas-grain chemical models. In order to examine this scheme, we derived an expression to account for the formation of hydrocarbons when HACs are photo-decomposed after their injection from grain mantles. Then, we calculated the actual formation rate of these species by knowing their injected fraction (from experimental work) and the average rate of mantle carbon injection into the ISM (from observations). Our preliminary results are promising and reveal that the degradation scheme can be considered as an efficient mechanism for the formation of some simple hydrocarbons in diffuse clouds. However, an actual proof of the efficiency of this process and its rate constants would require comprehensive experimental determination.

preprint2022arXiv

The Effects of Cosmic Rays on the Chemistry of Dense Cores

Cosmic rays are crucial for the chemistry of molecular clouds and their evolution. They provide essential ionizations, dissociations, heating and energy to the cold, dense cores. As cosmic rays pierce through the clouds they are attenuated and lose energy, which leads to a dependency on the column density of a system. The detailed effects these particles have on the central regions still needs to be fully understood. Here, we revisit how cosmic rays are treated in the UCLCHEM chemical modeling code by including both ionization rate and H2 dissociation rate dependencies alongside the production of cosmic ray induced excited species and we study in detail the effects of these treatments on the chemistry of pre-stellar cores. We find that these treatments can have significant effects on chemical abundances, up to several orders of magnitude, depending on physical conditions. The ionization dependency is the most significant treatment, influencing chemical abundances through increased presence of ionized species, grain desorptions and enhanced chemical reactions. Comparisons to chemical abundances derived from observations show the new treatments reproduce these observations better than the standard handling. It is clear that more advanced treatments of cosmic rays are essential to chemical models and that including this type of dependency provides more accurate chemical representations.

preprint2022arXiv

UCLCHEMCMC: A MCMC Inference tool for Physical Parameters of Molecular Clouds

We present the publicly available, open source code UCLCHEMCMC, designed to estimate physical parameters of an observed cloud of gas by combining Monte Carlo Markov Chain (MCMC) sampling with chemical and radiative transfer modeling. When given the observed values of different emission lines, UCLCHEMCMC runs a Bayesian parameter inference, using a MCMC algorithm to sample the likelihood and produce an estimate of the posterior probability distribution of the parameters. UCLCHEMCMC takes a full forward modeling approach, generating model observables from the physical parameters via chemical and radiative transfer modeling. While running UCLCHEMCMC, the created chemical models and radiative transfer code results are stored in an SQL database, preventing redundant model calculations in future inferences. This means that the more UCLCHEMCMC is used, the more efficient it becomes. Using UCLCHEM and RADEX, the increase of efficiency is nearly two orders of magnitude, going from 5185.33 \pm 1041.96 s for ten walkers to take one thousand steps when the database is empty, to 68.89 \pm 45.39 s when nearly all models requested are in the database. In order to demonstrate its usefulness we provide an example inference of UCLCHEMCMC to estimate the physical parameters of mock data, and perform two inferences on the well studied prestellar core, L1544, one of which show that it is important to consider the substructures of an object when determining which emission lines to use.

preprint2022arXiv

Using Statistical Emulation and Knowledge of Grain-Surface Diffusion for Bayesian Inference of Reaction Rate Parameters: An Application to a Glycine Network

There exists much uncertainty surrounding interstellar grain-surface chemistry. One of the major reaction mechanisms is grain-surface diffusion for which the the binding energy parameter for each species needs to be known. However, these values vary significantly across the literature which can lead to debate as to whether or not a particular reaction takes place via diffusion. In this work we employ Bayesian inference to use available ice abundances to estimate the reaction rates of the reactions in a chemical network that produces glycine. Using this we estimate the binding energy of a variety of important species in the network, by assuming that the reactions take place via diffusion. We use our understanding of the diffusion mechanism to reduce the dimensionality of the inference problem from 49 to 14, by demonstrating that reactions can be separated into classes. This dimensionality reduction makes the problem computationally feasible. A neural network statistical emulator is used to also help accelerate the Bayesian inference process substantially. The binding energies of most of the diffusive species of interest are found to match some of the disparate literature values, with the exceptions of atomic and diatomic hydrogen. The discrepancies with these two species are related to limitations of the physical and chemical model. However, the use of a dummy reaction of the form H + X -> HX is found to somewhat reduce the discrepancy with the binding energy of atomic hydrogen. Using the inferred binding energies in the full gas-grain version of UCLCHEM results in almost all the molecular abundances being recovered.

preprint2021arXiv

Constraining chemical networks inAstrochemistry

Databases of gas and surface chemical reactions are a key tool for scientists working in a wide range of physical sciences. In Astrochemistry, databases of chemical reactions are used as inputs to chemical models to determine the abundances of the interstellar medium. Gas chemistry and, in particular, grain surface chemistry and its treatment in gas-grain chemical models are however areas of large uncertainty. Many reactions - especially on the dust grains - have not been systematically experimentally studied. Moreover, experimental measurements are often not easily translated to the rate equation approach most commonly used in astrochemical modelling. Reducing the degree of uncertainty intrinsic in these databases is therefore a prime problem, but has so far been approached mainly by ad hoc procedures of essentially trial and error. In this chapter we review the problem of the determination of accurate and complete chemical networks in the wider context of Astrochemistry and explore the possibility of using statistical methods and machine learning (ML) techniques to reduce the uncertainty in chemical networks.

preprint2020arXiv

A chemical study of carbon fractionation in external galaxies

In the interstellar medium carbon exists in the form of two stable isotopes $^{12}$C and $^{13}$C and their ratio is a good indicator of nucleosynthesis in galaxies. However, chemical fractionation can potentially significantly alter this ratio and in fact observations of carbon fractionation within the same galaxy has been found to vary from species to species. In this paper we theoretically investigate the carbon fractionation for selected abundant carbon-bearing species in order to determine the conditions that may lead to a spread of the $^{12}$C/$^{13}$C ratio in external galaxies. We find that carbon fractionation is sensitive to almost all the physical conditions we investigated, it strongly varies with time for all species but CO, and shows pronounced differences across species. Finally we discuss our theoretical results in the context of the few observations of the $^{12}$C/$^{13}$C ratio in both local and higher redshift galaxies.

preprint2019arXiv

The CO(3-2)/CO(1-0) luminosity line ratio in nearby star-forming galaxies and AGN from xCOLD GASS, BASS and SLUGS

We study the r31=L'CO(3-2)/L'CO(1-0) luminosity line ratio in a sample of nearby (z < 0.05) galaxies: 25 star-forming galaxies (SFGs) from the xCOLD GASS survey, 36 hard X-ray selected AGN host galaxies from BASS and 37 infrared luminous galaxies from SLUGS. We find a trend for r31 to increase with star-formation efficiency (SFE). We model r31 using the UCL-PDR code and find that the gas density is the main parameter responsible for variation of r31, while the interstellar radiation field and cosmic ray ionization rate play only a minor role. We interpret these results to indicate a relation between SFE and gas density. We do not find a difference in the r31 value of SFGs and AGN host galaxies, when the galaxies are matched in SSFR (<r31>= 0.52 +/- 0.04 for SFGs and <r31> = 0.53 +/- 0.06 for AGN hosts). According to the results of UCL-PDR models, the X-rays can contribute to the enhancement of the CO line ratio, but only for strong X-ray fluxes and for high gas density (nH > 10$^4$ cm-3). We find a mild tightening of the Kennicutt-Schmidt relation when we use the molecular gas mass surface density traced by CO(3-2) (Pearson correlation coefficient R=0.83), instead of the molecular gas mass surface density traced by CO(1-0) (R=0.78), but the increase in correlation is not statistically significant (p-value=0.06). This suggests that the CO(3-2) line can be reliably used to study the relation between SFR and molecular gas for normal SFGs at high redshift, and to compare it with studies of low-redshift galaxies, as is common practice.

preprint2019arXiv

Tracing shock type with chemical diagnostics: an application to L1157

Aims: The physical structure of a shock wave may take a form unique to its shock type, implying that the chemistry of each shock type is unique as well. We aim to investigate the different chemistries of J-type and C-type shocks in order to identify unique molecular tracers of both shock types. We apply these diagnostics to the protostellar outflow L1157 to establish whether the B2 clump could host shocks exhibiting type-specific behaviour. Of particular interest is the L1157-B2 clump, which has been shown to exhibit bright emission in S-bearing species and HNCO. Methods: We simulate, using a parameterised approach, a planar, steady-state J-type shock wave using UCLCHEM. We compute a grid of models using both C-type and J-type shock models to determine the chemical abundance of shock-tracing species as a function of distance through the shock and apply it to the L1157 outflow. We focus on known shock-tracing molecules such as H2O, HCN, and CH3 OH. Results: We find that a range of molecules including H2O and HCN have unique behaviour specific to a J-type shock, but that such differences in behaviour are only evident at low $v_{s}$ and low $n_{H}$. We find that CH3OH is enhanced by shocks and is a reliable probe of the pre-shock gas density. However, we find no difference between its gas-phase abundance in C-type and J-type shocks. Finally, from our application to L1157, we find that the fractional abundances within the B2 region are consistent with both C-type and J-type shock emission.

preprint2016arXiv

H2S in the L1157-B1 Bow shock

Sulfur-bearing molecules are highly reactive in the gas phase of the ISM. However, the form in which most of the sulfur is locked onto interstellar dust grains is unknown. By taking advantage of the short time-scales of shocks in young molecular outflows, one could track back the main form of sulfur in the ices. In this paper, six transitions of H$_2$S and its isotopologues in the L1157-B1 bowshock have been detected using data from the Herschel-CHESS survey and the IRAM-30m ASAI large program. These detections are used to calculate the properties of H$_2$S gas in L1157-B1 through use of a rotation diagram and to explore the possible carriers of sulfur on the grains. The isotopologue detections allow the first calculation of the H$_2$S deuteration fraction in an outflow from a low mass protostar. The fractional abundance of H$_2$S in the region is found to be 6.0$\times$10$^{-7}$ and the deuteration fraction is 2$\times$10$^{-2}$. In order to investigate the form of sulfur on the grains, a chemical model is run with four different networks, each with different branching ratios for the freeze out of sulfur bearing species into molecules such as OCS and H$_2$S. It is found that the model best fits the data when at least half of each sulfur bearing species hydrogenates when freezing. We therefore conclude that a significant fraction of sulfur in L1157-B1 is likely to be locked in H$_2$S on the grains.

preprint2016arXiv

Radiative transfer meets Bayesian statistics: where does a galaxy's [CII] emission come from?

The [CII] 158$μ$m emission line can arise in all phases of the ISM, therefore being able to disentangle the different contributions is an important yet unresolved problem when undertaking galaxy-wide, integrated [CII] observations. We present a new multi-phase 3D radiative transfer interface that couples Starburst99, a stellar spectrophotometric code, with the photoionisation and astrochemistry codes Mocassin and 3D-PDR. We model entire star forming regions, including the ionised, atomic and molecular phases of the ISM, and apply a Bayesian inference methodology to parametrise how the fraction of the [CII] emission originating from molecular regions, $f_{[CII],mol}$, varies as a function of typical integrated properties of galaxies in the local Universe. The main parameters responsible for the variations of $f_{[CII],mol}$ are specific star formation rate (sSFR), gas phase metallicity, HII region electron number density ($n_e$), and dust mass fraction. For example, $f_{[CII],mol}$ can increase from 60% to 80% when either $n_e$ increases from 10$^{1.5}$ to 10$^{2.5}$cm$^{-3}$, or SSFR decreases from $10^{-9.6}$ to $10^{-10.6}$ yr$^{-1}$. Our model predicts for the Milky Way that $f_{[CII],mol}$$=75.8\pm5.9$%, in agreement with the measured value of 75%. When applying the new prescription to a complete sample of galaxies from the Herschel Reference Survey (HRS), we find that anywhere from 60 to 80% of the total integrated [CII] emission arises from molecular regions.

preprint2016arXiv

The Spatial Distribution of Complex Organic Molecules in the L1544 Pre-stellar Core

The detection of complex organic molecules (COMs) toward cold sources such as pre-stellar cores (with T<10 K), has challenged our understanding of the formation processes of COMs in the interstellar medium. Recent modelling on COM chemistry at low temperatures has provided new insight into these processes predicting that COM formation depends strongly on parameters such as visual extinction and the level of CO freeze out. We report deep observations of COMs toward two positions in the L1544 pre-stellar core: the dense, highly-extinguished continuum peak with Av>=30 mag within the inner 2700 au; and a low-density shell with average Av~7.5-8 mag located at 4000 au from the core's center and bright in CH3OH. Our observations show that CH3O, CH3OCH3 and CH3CHO are more abundant (by factors ~2-10) toward the low-density shell than toward the continuum peak. Other COMs such as CH3OCHO, c-C3H2O, HCCCHO, CH2CHCN and HCCNC show slight enhancements (by factors <=3) but the associated uncertainties are large. This suggests that COMs are actively formed and already present in the low-density shells of pre-stellar cores. The modelling of the chemistry of O-bearing COMs in L1544 indicates that these species are enhanced in this shell because i) CO starts freezing out onto dust grains driving an active surface chemistry; ii) the visual extinction is sufficiently high to prevent the UV photo-dissociation of COMs by the external interstellar radiation field; and iii) the density is still moderate to prevent severe depletion of COMs onto grains.

preprint2015arXiv

A new study of an old sink of sulfur in hot molecular cores: the sulfur residue

Sulfur appears to be depleted by an order of magnitude or more from its elemental abundance in star-forming regions. In the last few years, numerous observations and experiments have been performed in order to to understand the reasons behind this depletion without providing a satisfactory explanation of the sulfur chemistry towards high-mass star-forming cores. Several sulfur-bearing molecules have been observed in these regions, and yet none are abundant enough to make up the gas-phase deficit. Where, then, does this hidden sulfur reside? This paper represents a step forward in our understanding of the interactions among the various S-bearing species. We have incorporated recent experimental and theoretical data into a chemical model of a hot molecular core in order to see whether they give any indication of the identity of the sulfur sink in these dense regions. Despite our model producing reasonable agreement with both solid-phase and gas-phase abundances of many sulfur-bearing species, we find that the sulfur residue detected in recent experiments takes up only ~6 per cent of the available sulfur in our simulations, rather than dominating the sulfur budget.

preprint2015arXiv

Chemical Tracers of Pre-Brown Dwarf Cores Formed Through Turbulent Fragmentation

A gas-grain time dependent chemical code, UCL\_CHEM, has been used to investigate the possibility of using chemical tracers to differentiate between the possible formation mechanisms of brown dwarfs. In this work, we model the formation of a pre-brown dwarf core through turbulent fragmentation by following the depth-dependent chemistry in a molecular cloud through the step change in density associated with an isothermal shock and the subsequent freefall collapse once a bound core is produced. Trends in the fractional abundance of molecules commonly observed in star forming cores are then explored to find a diagnostic for identifying brown dwarf mass cores formed through turbulence. We find that the cores produced by our models would be bright in CO and NH$_3$ but not in HCO$^+$. This differentiates them from models using purely freefall collapse as such models produce cores that would have detectable transitions from all three molecules.

preprint2015arXiv

Effective destruction of CO by cosmic rays: implications for tracing H$_2$ gas in the Universe

We report on the effects of cosmic rays (CRs) on the abundance of CO in $\rm H_2$ clouds under conditions typical for star-forming galaxies in the Universe. We discover that this most important molecule for tracing H$_2$ gas is very effectively destroyed in ISM environments with CR energy densities $\rm U_{CR}\sim(50-10^{3})\times U_{CR,Gal}$, a range expected in numerous star-forming systems throughout the Universe. This density-dependent effect operates volumetrically rather than only on molecular cloud surfaces (i.e. unlike FUV radiation that also destroys CO), and is facilitated by: a) the direct destruction of CO by CRs, and b) a reaction channel activated by CR-produced He$^{+}$. The effect we uncover is strong enough to render Milky-Way type Giant Molecular Clouds (GMCs) very CO-poor (and thus CO-untraceable), even in ISM environments with rather modestly enhanced average CR energy densities of $\rm U_{CR}\sim(10-50)\times\rm U_{CR,Gal}$. We conclude that the CR-induced destruction of CO in molecular clouds, unhindered by dust absorption, is perhaps the single most important factor controlling the CO-visibility of molecular gas in vigorously star-forming galaxies. We anticipate that a second order effect of this CO destruction mechanism will be to make the H$_2$ distribution in the gas-rich disks of such galaxies appear much clumpier in CO $J$=1--0, 2--1 line emission than it actually is. Finally we give an analytical approximation of the CO/H$_2$ abundance ratio as a function of gas density and CR energy density for use in galaxy-size or cosmological hydrodynamical simulations, and propose some key observational tests.

preprint2014arXiv

A high resolution study of complex organic molecules in hot cores

We present the results of a line identification analysis using data from the IRAM Plateau de Bure Inferferometer, focusing on six massive star-forming hot cores: G31.41+0.31, G29.96-0.02, G19.61-0.23, G10.62-0.38, G24.78+0.08A1 and G24.78+0.08A2. We identify several transitions of vibrationally excited methyl formate (HCOOCH$_3$) for the first time in these objects as well as transitions of other complex molecules, including ethyl cyanide (C$_2$H$_5$CN), and isocyanic acid (HNCO). We also postulate a detection of one transition of glycolaldehyde (CH$_2$(OH)CHO) in two new hot cores. We find G29.96-0.02, G19.61-0.23, G24.78+0.08A1 and 24.78+0.08A2 to be chemically very similar. G31.41+0.31, however, is chemically different: it manifests a larger chemical inventory and has significantly larger column densities. We suggest that it may represent a different evolutionary stage to the other hot cores in the sample, or it may surround a star with a higher mass. We derive column densities for methyl formate in G31.41+0.31, using the rotation diagram method, of $\times$10$^{17}$ cm$^{-2}$ and a T$_{rot}$ of $\sim$170 K. For G29.96-0.02, G24.78+0.08A1 and G24.78+0.08A2, glycolaldehyde, methyl formate and methyl cyanide all seem to trace the same material and peak at roughly the same position towards the dust emission peak. For G31.41+0.31, however, glycolaldehyde shows a different distribution to methyl formate and methyl cyanide and seems to trace the densest, most compact inner part of hot cores.

preprint2014arXiv

An Alternative Accurate Tracer of Molecular Clouds: The "$X_{\rm CI}$-Factor"

We explore the utility of CI as an alternative high-fidelity gas mass tracer for Galactic molecular clouds. We evaluate the X$_{\rm CI}$-factor for the 609 $μ$m carbon line, the analog of the CO X-factor, which is the ratio of the H$_2$ column density to the integrated $^{12}$CO(1-0) line intensity. We use 3D-PDR to post-process hydrodynamic simulations of turbulent, star-forming clouds. We compare the emission of CI and CO for model clouds irradiated by 1 and 10 times the average background and demonstrate that CI is a comparable or superior tracer of the molecular gas distribution for column densities up to $6 \times 10^{23}$ cm$^{-2}$. Our results hold for both reduced and full chemical networks. For our fiducial Galactic cloud we derive an average $X_{\rm CO}$ of $3.0\times 10^{20}$ cm$^{-2}$K$^{-1}$km$^{-1}$s and $X_{\rm CI}$ of $1.1\times 10^{21}$ cm$^{-2}$K$^{-1}$km$^{-1}$s.

preprint2014arXiv

Detectability of Glycine in Solar-type System Precursors

Glycine (NH2CH2COOH) is the simplest amino acid relevant for life. Its detection in the interstellar medium is key to understand the formation mechanisms of pre-biotic molecules and their subsequent delivery onto planetary systems. Glycine has extensively been searched for toward hot molecular cores, although these studies did not yield any firm detection. In contrast to hot cores, low-mass star forming regions, and in particular their earliest stages represented by cold pre-stellar cores, may be better suited for the detection of glycine as well as more relevant for the study of pre-biotic chemistry in young Solar System analogs. We present 1D spherically symmetric radiative transfer calculations of the glycine emission expected to arise from the low-mass pre-stellar core L1544. Water vapour has recently been reported toward this core, indicating that a small fraction of the grain mantles in L1544 (~0.5%) has been injected into the gas phase. Assuming that glycine is photo-desorbed together with water in L1544, and considering a solid abundance of glycine on ices of ~1E-4 with respect to water, our calculations reveal that several glycine lines between 67 GHz and 80 GHz have peak intensities larger than 10 mK. These results show for the first time that glycine could reach detectable levels in cold objects such as L1544. This opens up the possibility to detect glycine, and other pre-biotic species, at the coldest and earliest stages in the formation of Solar-type systems with near-future instrumentation such as the Band 2 receivers of ALMA.

preprint2014arXiv

Deuterium chemistry of dense gas in the vicinity of low-mass and massive star forming regions

The standard interstellar ratio of deuterium to hydrogen (D/H) atoms is $\sim 1.5 \times 10^{-5}$. However, the deuterium fractionation is in fact found to be enhanced, to different degrees, in cold, dark cores, hot cores around massive star forming regions, lukewarm cores, and warm cores ({\it hereafter}, hot corinos) around low-mass star forming regions. In this paper, we investigate the overall differences in the deuterium chemistry between hot cores and hot corinos. We have modelled the chemistry of dense gas around low-mass and massive star forming regions using a gas-grain chemical model. We investigate the influence of varying the core density, the depletion efficiency of gaseous species on to dust grains, the collapse mode and the final mass of the protostar on the chemical evolution of star forming regions. We find that the deuterium chemistry is, in general, most sensitive to variations of the depletion efficiency on to grain surfaces, in agreement with observations. In addition, the results showed that the chemistry is more sensitive to changes in the final density of the collapsing core in hot cores than in hot corinos. Finally, we find that ratios of deuterated sulphur bearing species in dense gas around hot cores and corinos may be good evolutionary indicators in a similar way as their non deuterated counterparts.

preprint2014arXiv

Herschel HIFI observations of O$_2$ toward Orion: special conditions for shock enhanced emission

We report observations of molecular oxygen (O$_2$) rotational transitions at 487 GHz, 774 GHz, and 1121 GHz toward Orion Peak A. The O2 lines at 487 GHz and 774 GHz are detected at velocities of 10-12 km/s with line widths 3 km/s; however, the transition at 1121 GHz is not detected. The observed line characteristics, combined with the results of earlier observations, suggest that the region responsible for the O$_2$ emission is 9" (6e16 cm) in size, and is located close to the H2 Peak 1position (where vibrationally-excited H$_2$ emission peaks), and not at Peak A, 23" away. The peak O2 column density is 1.1e18/cm2. The line velocity is close to that of 621 GHz water maser emission found in this portion of the Orion Molecular Cloud, and having a shock with velocity vector lying nearly in the plane of the sky is consistent with producing maximum maser gain along the line-of-sight. The enhanced O$_2$ abundance compared to that generally found in dense interstellar clouds can be explained by passage of a low-velocity C-shock through a clump with preshock density 2e4/cm3, if a reasonable flux of UV radiation is present. The postshock O$_2$ can explain the emission from the source if its line of sight dimension is ~10 times larger than its size on the plane of the sky. The special geometry and conditions required may explain why O$_2$ emission has not been detected in the cores of other massive star-forming molecular clouds.

preprint2014arXiv

Understanding the Formation and Evolution of Interstellar Ices: A Bayesian Approach

Understanding the physical conditions of dark molecular clouds and star forming regions is an inverse problem subject to complicated chemistry that varies non-linearly with time and the physical environment. In this paper we apply a Bayesian approach based on a Markov Chain Monte Carlo (MCMC) method for solving the non-linear inverse problems encountered in astrochemical modelling. We use observations for ice and gas species in dark molecular clouds and a time dependent, gas grain chemical model to infer the values of the physical and chemical parameters that characterize quiescent regions of molecular clouds. We show evidence that in high dimensional problems, MCMC algorithms provide a more efficient and complete solution than more classical strategies. The results of our MCMC method enable us to derive statistical estimates and uncertainties for the physical parameters of interest as a result of the Bayesian treatment.

preprint2013arXiv

An improved chemical scheme for the reactions of atomic oxygen and simple unsaturated hydrocarbons - implications for star-forming regions

Recent laboratory experiments have demonstrated that, even though contribution from other reaction channels cannot be neglected, unsaturated hydrocarbons easily break their multiple C-C bonds to form CO after their interactions with atomic oxygen. Here we present an upgraded chemical modelling including a revision of the reactions between oxygen atoms and small unsaturated hydrocarbons for different astrochemical environments. A first conclusion is that towards hot cores/corinos atomic oxygen easily degrades unsaturated hydrocarbons directly to CO or to its precursor species (such as HCCO or HCO) and destroys the double or triple bond of alkenes and alkynes. Therefore, environments rich in atomic oxygen at a relatively high temperature are not expected to be rich in large unsaturated hydrocarbons or polycyclic aromatic hydrocarbons. On the contrary, in O-poor and C-rich objects, hydrocarbon growth can occur to a large extent. In addition, new radical species, namely ketyl and vinoxy radicals, generated from other reaction channels can influence the abundances of hydrocarbons towards hot cores. We, therefore, suggest they should be included in the available databases. Hydrocarbon column densities are calculated in the 10(13)-10(15) cm(-2) range, in good agreement with their observed values, despite the small number of data currently published in the literature.

preprint2013arXiv

Modeling the Atomic-to-Molecular Transition and Chemical Distributions of Turbulent Star-Forming Clouds

We use 3D-PDR, a three-dimensional astrochemistry code for modeling photodissociation regions (PDRs), to post-process hydrodynamic simulations of turbulent, star-forming clouds. We focus on the transition from atomic to molecular gas, with specific attention to the formation and distribution of H, C+, C, H2 and CO. First, we demonstrate that the details of the cloud chemistry and our conclusions are insensitive to the simulation spatial resolution, to the resolution at the cloud edge, and to the ray angular resolution. We then investigate the effect of geometry and simulation parameters on chemical abundances and find weak dependence on cloud morphology as dictated by gravity and turbulent Mach number. For a uniform external radiation field, we find similar distributions to those derived using a one-dimensional PDR code. However, we demonstrate that a three-dimensional treatment is necessary for a spatially varying external field, and we caution against using one-dimensional treatments for non-symmetric problems. We compare our results with the work of Glover et al. (2010), who self-consistently followed the time evolution of molecule formation in hydrodynamic simulations using a reduced chemical network. In general, we find good agreement with this in situ approach for C and CO abundances. However, the temperature and H2 abundances are discrepant in the boundary regions (Av < 5), which is due to the different number of rays used by the two approaches.

preprint2013arXiv

Molecules in the Circumnuclear Disk of the Galactic Center

Within a few parsecs around the central Black Hole Sgr A*, chemistry in the dense molecular cloud material of the circumnuclear disk (CND) can be affected by many energetic phenomena such as high UV-flux from the massive central star cluster, X-rays from Sgr A*, shock waves, and an enhanced cosmic-ray flux. Recently, spectroscopic surveys with the IRAM 30 meter and the APEX 12 meter telescopes of substantial parts of the 80--500 GHz frequency range were made toward selected positions in and near the CND. These datasets contain lines from the molecules HCN, HCO$^+$, HNC, CS, SO, SiO, CN, H$_2$CO, HC$_3$N, N$_2$H$^+$, H$_3$O$^+$ and others. We conduct Large Velocity Gradient analyses to obtain column densities and total hydrogen densities, $n$, for each species in molecular clouds located in the southwest lobe of CND. The data for the above mentioned molecules indicate 10$^5\,$cm$^{-3} \lesssim n <10^6\,$cm$^{-3}$, which shows that the CND is tidally unstable. The derived chemical composition is compared with a chemical model calculated using the UCL\_CHEM code that includes gas and grain reactions, and the effects of shock waves. Models are run for varying shock velocities, cosmic-ray ionization rates, and number densities. The resulting chemical composition is fitted best to an extremely high value of cosmic-ray ionization rate $ζ\sim 10^{-14}\,$s$^{-1}$, 3 orders of magnitude higher than the value in regular Galactic molecular clouds, if the pre-shock density is $n=10^5\,$cm$^{-3}$.

preprint2012arXiv

A lambda=3 mm molecular line survey of NGC1068. Chemical signatures of an AGN environment

We aimed to study the molecular composition of the interstellar medium (ISM) surrounding an Active Galactic Nucleus (AGN), by making an inventory of molecular species and their abundances, as well as to establish a chemical differentiation between starburst galaxies and AGN. We used the IRAM-30 m telescope to observe the central 1.5-2 kpc region of NGC1068, covering the frequencies between 86.2 GHz and 115.6 GHz. Using Boltzmann diagrams, we calculated the column densities of the detected molecules. We used a chemical model to reproduce the abundances found in the AGN, to determine the origin of each detected species, and to test the influence of UV fields, cosmic rays, and shocks on the ISM. We identified 24 different molecular species and isotopologues, among which HC3N, SO, N2H+, CH3CN, NS, 13CN, and HN13C are detected for the first time in NGC1068. We obtained the upper limits to the isotopic ratios 12C/13C=49, 16O/18O=177 and 32S/34S=5. Our chemical models suggest that the chemistry in the nucleus of NGC1068 is strongly influenced by cosmic rays, although high values of both cosmic rays and far ultraviolet (FUV) radiation fields also explain well the observations. The gas in the nucleus of NGC1068 has a different chemical composition as compared to starburst galaxies. The distinct physical processes dominating galaxy nuclei (e.g. C-shocks, UV fields, X-rays, cosmic rays) leave clear imprints in the chemistry of the gas, which allow to characterise the nucleus activity by its molecular abundances.

preprint2012arXiv

Chemical complexity in NGC1068

We aimed to study the chemistry of the circumnuclear molecular gas of NGC1068, and to compare it with those of the starburst galaxies M82 and NGC253. Using the IRAM-30m telescope, we observed the inner 2 kpc of NGC1068 between 86.2 GHz and 115.6 GHz. We identified 35 spectral features, corresponding to 24 different molecular species. Among them, HC3N, SO, N2H+, CH3CN, NS, 13CN, and HN13C are detected for the first time in NGC1068. Assuming local thermodynamic equilibrium (LTE), we calculated the column densities of the detected molecules, as well as the upper limits to the column densities of some undetected species. The comparison among the chemistries of NGC1068, M82, and NGC253, suggests that, apart from X-rays, shocks also determine the chemistry of NGC1068. We propose the column density ratio between CH3CCH and HC3N as a prime indicator of the imprints of starburst and AGN environments in the circumnuclear interstellar medium. This ratio is, at least, 64 times larger in M82 than in NGC1068, and, at least, 4 times larger in NGC253 than in NGC1068. Finally, we used the UCL_CHEM and UCL_PDR chemical codes to constrain the origin of the species, as well as to test the influence of UV radiation fields and cosmic rays on the observed abundances.

preprint2012arXiv

On the formation of glycolaldehyde in dense molecular cores

Glycolaldehyde is a simple monosaccharide sugar linked to prebiotic chemistry. Recently it was detected in a molecular core in the star-forming region G31.41+0.31 at a reasonably high abundance. We investigate the formation of glycolaldehyde at 10K to determine whether it can form efficiently under typical dense core conditions. Using an astrochemical model, we test five different reaction mechanisms that have been proposed in the astrophysical literature, finding that a gas-phase formation route is unlikely. Of the grain-surface formation routes, only two are efficient enough at very low temperatures to produce sufficient glycolaldehyde to match the observational estimates, with the mechanism culminating in CH3OH + HCO being favoured. However, when we consider the feasibility of these mechanisms from a reaction chemistry perspective, the second grain-surface route looks more promising, H3CO + HCO.

preprint2012arXiv

Small-scale structure in the interstellar medium: time-varying interstellar absorption towards κ Velorum

Ultra-high spectral resolution observations of time-varying interstellar absorption towards κ Vel are reported, using the Ultra-High Resolution Facility on the Anglo-Australian Telescope. Detections of interstellar Ca I, Ca II, K I, Na I and CH are obtained, whilst an upper limit on the column density is reported for C_2. The results show continued increases in column densities of K I and Ca I since observations ~ 4 yr earlier, as the transverse motion of the star carried it ~ 10 AU perpendicular to the line of sight. Line profile models are fitted to the spectra and two main narrow components (A & B) are identified for all species except CH. The column density N(K I) is found to have increased by 82 +10-9 % between 1994 and 2006, whilst N(Ca I) is found to have increased by 32 +- 5 % over the shorter period of 2002-2006. The line widths are used to constrain the kinetic temperature to T_k,A < 671 +18-17 K and T_k,B < 114 +15-14 K. Electron densities are determined from the Ca I / Ca II ratio, which in turn place lower limits on the total number density of n_A > 7 * 10^3 cm^-3 and n_B > 2 * 10^4 cm^-3. Calcium depletions are estimated from the Ca I / K I ratio. Comparison with the chemical models of Bell et al. (2005) confirms the high number density, with n = 5 * 10^4 cm^-3 for the best-fitting model. The first measurements of diffuse interstellar bands (DIBs) towards this star are made at two epochs, but only an upper limit of < 40 % is placed on their variation over ~ 9 years. The DIBs are unusually weak for the measured E(B-V) and appear to exhibit similar behaviour to that seen in Orion. The ratio of equivalent widths of the λ5780 to λ5797 DIBs is amongst the highest known, which may indicate that the carrier of λ5797 is more sensitive to UV radiation than to local density.

preprint2012arXiv

The influence of cosmic rays in the circumnuclear molecular gas of NGC1068

We surveyed the circumnuclear disk of the Seyfert galaxy NGC1068 between the frequencies 86.2 GHz and 115.6 GHz, and identified 17 different molecules. Using a time and depth dependent chemical model we reproduced the observational results, and show that the column densities of most of the species are better reproduced if the molecular gas is heavily pervaded by a high cosmic ray ionization rate of about 1000 times that of the Milky Way. We discuss how molecules in the NGC1068 nucleus may be influenced by this external radiation, as well as by UV radiation fields.

preprint2011arXiv

A study of methyl formate in astrochemical environments

Several complex organic molecules are routinely detected in high abundances towards hot cores and hot corinos. For many of them, their paths of formation in space are uncertain, as gas phase reactions alone seem to be insufficient. In this paper, we investigate a possible solid-phase route of formation for methyl formate (HCOOCH3). We use a chemical model updated with recent results from an experiment where simulated grain surfaces were irradiated with 200 keV protons at 16 K, to simulate the effects of cosmic ray irradiation on grain surfaces. We find that this model may be sufficient to reproduce the observed methyl formate in dark clouds, but not that found in hot cores and corinos.

preprint2011arXiv

Chemistry of dense clumps near moving Herbig-Haro objects

Localised regions of enhanced emission from HCO+, NH3 and other species near Herbig-Haro objects (HHOs) have been interpreted as arising in a photochemistry stimulated by the HHO radiation on high density quiescent clumps in molecular clouds. Static models of this process have been successful in accounting for the variety of molecular species arising ahead of the jet; however recent observations show that the enhanced molecular emission is widespread along the jet as well as ahead. Hence, a realistic model must take into account the movement of the radiation field past the clump. It was previously unclear as to whether the short interaction time between the clump and the HHO in a moving source model would allow molecules such as HCO+ to reach high enough levels, and to survive for long enough to be observed. In this work we model a moving radiation source that approaches and passes a clump. The chemical picture is qualitatively unchanged by the addition of the moving source, strengthening the idea that enhancements are due to evaporation of molecules from dust grains. In addition, in the case of several molecules, the enhanced emission regions are longer-lived. Some photochemically-induced species, including methanol, are expected to maintain high abundances for ~10,000 years.

preprint2011arXiv

Extreme Cosmic-Ray-Dominated-Regions: a new paradigm for high star formation density events in the Universe

We examine in detail the recent proposal that extreme Cosmic-Ray-Dominated-Regions (CRDRs) characterize the ISM of galaxies during events of high-density star formation, fundamentally altering its initial conditions (Papadopoulos 2010). Solving the coupled chemical and thermal state equations for dense UV-shielded gas reveals that the large cosmic ray energy densities in such systems (U_{CR} (few)x(10^3-10^4) U_{CR,Gal}) will indeed raise the minimum temperature of this phase (where the initial conditions of star formation are set) from ~10K (as in the Milky Way) to (50-100)K. Moreover in such extreme CRDRs the gas temperature remains fully decoupled from that of the dust, with T_{kin} >> T_{dust}, even at high densities (n(H_2)~10^5--10^6 cm^{-3}), quite unlike CRDRs in the Milky Way where T_k T_{dust} when n(H_2) >= 10^5 cm^{-3}. These dramatically different star formation initial conditions will: a) boost the Jeans mass of UV-shielded gas regions by factors of ~10--100 with respect to those in quiescent or less extreme star forming systems, and b) "erase" the so-called inflection point of the effective equation of state (EOS) of molecular gas. Both these effects occur across the entire density range of typical molecular clouds, and may represent {\it a new paradigm for all high-density star formation in the Universe}, with cosmic rays as the key driving mechanism, operating efficiently even in the high dust extinction environments of extreme starbursts...

preprint2011arXiv

Herschel Measurements of Molecular Oxygen in Orion

We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2.

preprint2011arXiv

Multi-line spectral imaging of dense cores in the Lupus molecular cloud

The molecular clouds Lupus 1, 3 and 4 were mapped with the Mopra telescope at 3 and 12 mm. Emission lines from high density molecular tracers were detected, i.e. NH$_3$ (1,1), NH$_3$ (2,2), N$_2$H$^+$ (1-0), HC$_3$N (3-2), HC$_3$N (10-9), CS (2-1), CH$_3$OH (2$_0-1_0$)A$^+$ and CH$_3$OH (2$_{-1}-1_{-1}$)E. Velocity gradients of more than 1 km s$^{-1}$ are present in Lupus 1 and 3 and multiple gas components are present in these clouds along some lines of sight. Lupus 1 is the cloud richest in high density cores, 8 cores were detected in it, 5 cores were detected in Lupus 3 and only 2 in Lupus 4. The intensity of the three species HC$_3$N, NH$_3$ and N$_2$H$^+$ changes significantly in the various cores: cores that are brighter in HC$_3$N are fainter or undetected in NH$_3$ and N$_2$H$^+$ and vice versa. We found that the column density ratios HC$_3$N/N$_2$H$^+$ and HC$_3$N/NH$_3$ change by one order of magnitude between the cores, indicating that also the chemical abundance of these species is different. The time dependent chemical code that we used to model our cores shows that the HC$_3$N/N$_2$H$^+$ and HC$_3$N/NH$_3$ ratios decrease with time therefore the observed column density of these species can be used as an indicator of the chemical evolution of dense cores. On this base we classified 5 out of 8 cores in Lupus 1 and 1 out of 5 cores in Lupus 3 as very young protostars or prestellar cores. Comparing the millimetre cores population with the population of the more evolved young stellar objects identified in the Spitzer surveys, we conclude that in Lupus 3 the bulk of the star formation activity has already passed and only a moderate number of stars are still forming. On the contrary, in Lupus 1 star formation is on-going and several dense cores are still in the pre--/proto--stellar phase. Lupus 4 is at an intermediate stage, with a smaller number of individual objects.

preprint2010arXiv

Chemistry in Evaporating Ices: Unexplored Territory

We suggest that three-body chemistry may occur in warm high density gas evaporating in transient co\textendash desorption events on interstellar ices. Using a highly idealised computational model we explore the chemical conversion from simple species of the ice to more complex species containing several heavy atoms, as a function of density and of adopted three body rate coefficients. We predict that there is a wide range of densities and rate coefficients in which a significant chemical conversion may occur. We discuss the implications of this idea for the astrochemistry of hot cores.

preprint2010arXiv

Formation Pumping of Molecular Hydrogen in Dark Clouds

Many theoretical and laboratory studies predict H2 to be formed in highly excited ro-vibrational states. The consequent relaxation of excited levels via a cascade of infrared transitions might be observable in emission from suitable interstellar regions. In this work, we model H2 formation pumping in standard dense clouds, taking into account the H/H2 transition zone, through an accurate description of chemistry and radiative transfer. The model includes recent laboratory data on H2 formation, as well as the effects of the interstellar UV field, predicting the populations of gas-phase H2 molecules and their IR emission spectra. Calculations suggest that some vibrationally excited states of H2 might be detectable towards lines of sight where significant destruction of H2 occurs, such as X-ray sources, and provide a possible explanation as to why observational attempts resulted in no detections reported to date.

preprint2010arXiv

N2H+ depletion in the massive protostellar cluster AFGL 5142

We aim at investigating with high angular resolution the NH3/N2H+ ratio toward the high-mass star-forming region AFGL 5142 in order to study whether this ratio behaves similarly to the low-mass case, for which the ratio decreases from starless cores to cores associated with YSOs. CARMA was used to observe the 3.2 mm continuum and N2H+(1-0) emission. We used NH3(1,1) and (2,2), HCO+(1-0) and H13CO+(1-0) data from the literature and we performed a time-dependent chemical modeling of the region. The 3.2 mm continuum emission reveals a dust condensation of ~23 Msun associated with the massive YSOs, deeply embedded in the strongest NH3 core (hereafter central core). The N2H+ emission reveals two main cores, the western and eastern core, located to the west and to the east of the mm condensation, and surrounded by a more extended and complex structure of ~0.5 pc. Toward the central core the N2H+ emission drops significantly, indicating a clear chemical differentiation in the region. We found low values of the NH3/N2H+ ratio ~50-100 toward the western/eastern cores, and high values up to 1000 in the central core. The chemical model indicates that density, and in particular temperature, are key parameters in determining the NH3/N2H+ ratio. The high density and temperature reached in the central core allow molecules like CO to evaporate from grain mantles. The CO desorption causes a significant destruction of N2H+, favoring the formation of HCO+. This result is supported by our observations, which show that N2H+ and HCO+ are anticorrelated in the central core. The observed values of the NH3/N2H+ ratio in the central core can be reproduced by our model for times t~4.5-5.3x10^5 yr (central) and t~10^4-3x10^6 yr (western/eastern). The NH3/N2H+ ratio in AFGL 5142 does not follow the same trend as in regions of low-mass star formation mainly due to the high temperature reached in hot cores.

preprint2010arXiv

The magnetic precursor of L1448-mm: Excitation differences between ion and neutral fluids

Shock modelling predicts an electron density enhancement within the magnetic precursor of C-shocks. Previous observations of SiO, H13CO+, HN13C and H13CN toward the young L1448-mm outflow showed an over-excitation of the ion fluid that was attributed to an electron density enhancement in the precursor. We re-visit this interpretation and test if it still holds when we consider different source morphologies and kinetic temperatures for the observed molecules, and also give some insight on the spatial extent of the electron density enhancement around L1448-mm. We estimate the opacities of H13CO+ and HN13C by observing the J=3\to2 lines of rarer isotopologues to confirm that the emission is optically thin. To model the excitation of the molecules, we use the large velocity gradient (LVG) approximation with updated collisional coefficients to i) re- analyse the observations toward the positions where the over-excitation of H13CO+ has previously been observed [i.e. toward L1448- mm at offsets (0,0) and (0,-10)], and ii) to investigate if the electron density enhancement is still required for the cases of extended and compact emission, and for kinetic temperatures of up to 400 K. We also report several lines of SiO, HN13C and H13CO+ toward new positions around this outflow, to investigate the spatial extent of the over-excitation of the ions in L1448-mm. From the isotopologue observations, we find that the emission of H13CO+ and HN13C from the precursor is optically thin if this emission is extended. Using the new collisional coefficients, an electron density enhancement is still needed to explain the excitation of H13CO+ for extended emission and for gas temperatures of\le 400 K toward L1448-mm (0,-10), and possibly also toward L1448-mm (0,0). For compact emission the data cannot be fitted. We do not find any evidence for the over-excitation of the ion fluid toward the newly observed positions around L1448-mm. The observed line emission of SiO, H13CO+ and HN13C toward L1448-mm (0,0) and (0,-10) is consistent with an electron density enhancement in the precursor component, if this emission is spatially extended. This is also true for the case of high gas temperatures (\le400 K) toward the (0,-10) offset. The electron density enhancement seems to be restricted to the southern, redshifted lobe of the L1448-mm outflow. Interferometric images of the line emission of these molecules are needed to confirm the spatial extent of the over-excitation of the ions and thus, of the electron density enhancement in the magnetic precursor of L1448-mm.

preprint2010arXiv

Warm Cores around Regions of Low-Mass Star Formation

Warm cores (or hot corinos) around low-mass protostellar objects show a rich chemistry with strong spatial variations. This chemistry is generally attributed to the sublimation of icy mantles on dust grains initiated by the warming effect of the stellar radiation. We have used a model of the chemistry in warm cores in which the sublimation process is based on extensive laboratory data; these data indicate that sublimation from mixed ices occurs in several well-defined temperature bands. We have determined the position of these bands for the slow warming by a solar-mass star. The resulting chemistry is dominated by the sublimation process and by subsequent gas-phase reactions; strong spatial and temporal variations in certain molecular species are found to occur, and our results are, in general, consistent with observational results for the well-studied source IRAS 16293-2422. The model used is similar to one that describes the chemistry of hot cores. We infer that the chemistry of both hot cores and warm cores may be described by the same model (suitably adjusted for different physical parameters).

preprint2009arXiv

Destruction of Molecular Gas Reservoirs in Early-Type Galaxies by Active Galactic Nucleus Feedback

Residual star formation at late times in early-type galaxies and their progenitors must be suppressed in order to explain the population of red, passively evolving systems we see today. Likewise, residual or newly accreted reservoirs of molecular gas that are fuelling star formation must be destroyed. This suppression of star formation in early-type galaxies is now commonly attributed to AGN feedback wherein the reservoir of gas is heated and expelled during a phase of accretion onto the central supermassive black hole. However, direct observational evidence for a link between the destruction of this molecular gas and an AGN phase has been missing so far. We present new mm-wavelength observations from the IRAM 30m telescope of a sample of low redshift SDSS early-type galaxies currently undergoing this process of quenching of late-time star formation. Our observations show that the disappearance of the molecular gas coincides within less than 100 Myr with the onset of accretion onto the black hole and is too rapid to be due to star formation alone. Since our sample galaxies are not associated to powerful quasar activity or radio jets, we conclude that low-luminosity AGN episodes are sufficient to suppress residual star formation in early-type galaxies. This `suppression mode' of AGN feedback is very different from the `truncation mode' linked to powerful quasar activity during early phases of galaxy formation.

preprint2009arXiv

The molecular emission of the irradiated star forming core ahead of HH 80N

We present a Berkeley-Illinois-Maryland Association (BIMA) Array molecular survey of the star forming core ahead of HH 80N, the optically obscured northern counterpart of the Herbig-Haro objects HH 80/81. Continuum emission at 1.4 mm and 8 micron is detected at the center of the core, which confirms the presence of an embedded very young stellar object in the core. All detected molecular species arise in a ring-like structure, which is most clearly traced by CS(2-1) emission. This molecular ring suggests that strong molecular depletion occurs in the inner part of the core (at a radius of ~0.1 pc and densities higher than ~5 x 10^4 cm^{-3}). Despite of the overall morphology and kinematic similarity between the different species, there is significant molecular differentiation along the ring-like structure. The analysis of the chemistry along the core shows that part of this differentiation may be caused by the UV irradiation of the nearby HH 80N object, that illuminates the part of the core facing HH 80N, which results in an abundance enhancement of some of the detected species.

preprint2009arXiv

Trends in Molecular Emission from Different Extragalactic Stellar Initial Mass Functions

Banerji et al. (2009) suggested that top-heavy stellar Initial Mass Functions (IMFs) in galaxies may arise when the interstellar physical conditions inhibit low-mass star formation, and they determined the physical conditions under which this suppression may or may not occur. In this work, we explore the sensitivity of the chemistry of interstellar gas under a wide range of conditions. We use these results to predict the relative velocity-integrated antenna temperatures of the CO rotational spectrum for several models of high redshift active galaxies which may produce both top-heavy and unbiased IMFs. We find that while active galaxies with solar metallicity (and top-heavy IMFs) produce higher antenna temperatures than those with sub-solar metallicity (and unbiased IMFs) the actual rotational distribution is similar. The high-J to peak CO ratio however may be used to roughly infer the metallicity of a galaxy provided we know whether it is active or quiescent. The metallicity strongly influences the shape of the IMF. High order CO transitions are also found to provide a good diagnostic for high far-UV intensity and low metallicity counterparts of Milky Way type systems both of which show some evidence for having top-heavy IMFs. We also compute the relative abundances of molecules known to be effective tracers of high density gas in these galaxy models. We find that the molecules CO and CS may be used to distinguish between solar and sub-solar metallicity in active galaxies at high redshift whereas HCN, HNC and CN are found to be relatively insensitive to the IMF shape at the large visual magnitudes typically associated with extragalactic sources.

preprint2008arXiv

Timescales for Low-Mass Star Formation in Extragalactic Environments: Implications for the Stellar IMF

We investigate the physical and chemical conditions necessary for low-mass star formation in extragalactic environments by calculating various characteristic timescales associated with star formation for a range of initial conditions. The balance of these timescales indicates whether low-mass star formation is enhanced or inhibited under certain physical conditions. In this study, we consider timescales for free-fall, cooling, freeze-out, desorption, chemistry and ambipolar diffusion and their variations with changes in the gas density, metallicity, cosmic ray ionisation rate and FUV radiation field strength. We find that extragalactic systems with high FUV radiation field strengths and high cosmic ray fluxes considered at a range of metallicities, are likely to have enhanced low-mass star formation unless the magnetic pressure is sufficient to halt collapse. Our results indicate that this is only likely to be the case for high-redshift galaxies approaching solar metallicities. Unless this is true for all high-redshift sources, this study finds little evidence for a high-mass biased IMF at high redshifts.