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

O. Sipilä

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

ResearcherUnclaimed source record
astro-ph.GA

Catalog footprint

What is connected

14works
1topics
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)

preprint2026arXiv

New Estimate for the Cosmic Ray-Induced $\rm H_2$ Photodissociation Rate in the Interstellar Medium

In the interstellar medium, cosmic rays (CRs) generate a field of ultraviolet (UV) photons via the excitation and subsequent radiative decay of $\rm H_2$ molecules. This UV field is a major agent of ionization and dissociation in the inner regions of molecular clouds that are shielded from the effects of the interstellar radiation field. In particular, the dissociation of $\rm H_2$, by far the most abundant molecule in interstellar clouds, leads to the production of atomic hydrogen which then takes part in the production of a multitude of molecules, in particular complex organics on the surfaces of interstellar dust grains. Precise knowledge of the rates of CR-induced dissociation processes is thus crucial for constructing reliable chemical models. For the present paper, we have derived a new value of $k_{\rm diss, CR}(\mbox{$\rm H_2$})=0.831ζ$ for the rate of $\rm H_2$ dissociation, where $ζ$ is the CR ionization rate of $\rm H_2$. This prediction contrasts a previous value from the Leiden database which overestimated the rate due to an inconsistent treatment of the $\rm H_2$ abundances and photodissociation cross sections. By running a series of chemical models, we show that the overestimated dissociation rate has a large effect on the results of chemical simulations, with the abundance of methanol being overestimated by over one order of magnitude. Hence, we strongly recommend the adoption of our new estimate $k_{\rm diss, CR}(\mbox{$\rm H_2$})=0.831ζ$ in all chemical models that include this process. Our newly derived value corresponds to $\rm H_2$ being purely in the para form ($J^{\prime\prime} = 0$). However, in the interiors of molecular clouds the $\rm H_2$ ortho-to-para ratio is low and using the rate for para-$\rm H_2$ is an adequate approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Deuteration of c-C$_3$H$_2$ towards the pre-stellar core L1544

Context: In the centre of pre-stellar cores, the deuterium fractionation is enhanced due to the cold temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to probe the evolution and the kinematics in the earliest stages of star formation. Aims: We analyse emission maps of cyclopropenylidene, c-C$_3$H$_2$, to study the distribution of the deuteration throughout the prototypical pre-stellar core L1544. Methods: We use single-dish observations of c-C$_3$H$_2$, c-H$^{13}$CC$_2$H, c-C$_3$HD, and c-C$_3$D$_2$ towards the pre-stellar core L1544, performed at the IRAM 30m telescope. We derive the column density and deuterium fraction maps, and compare these observations with non-LTE radiative transfer simulations. Results: The highest deuterium fractions are found close to the dust peak at the centre of L1544, where the increased abundance of H$_2$D$^+$ ions drives the deuteration process. The peak values are N(c-C$_3$HD)/N(c-C$_3$H$_2)=0.17\pm0.01$, N(c-C$_3$D$_2$)/N(c-C$_3$H$_2)=0.025\pm0.003$ and N(c-C$_3$D$_2$)/N(c-C$_3$HD$)=0.16\pm0.03$, which is consistent with previous single point observations. The distributions of c-C$_3$HD and c-C$_3$D$_2$ indicate that the deuterated forms of c-C$_3$H$_2$ in fact trace the dust peak and not the c-C$_3$H$_2$ peak. Conclusions: The N(c-C$_3$D$_2$)/N(c-C$_3$HD) map confirms that the process of deuteration is more efficient towards the centre of the core and demonstrates that carbon-chain molecules are still present at high densities. This is likely caused by an increased abundance of He$^+$ ions destroying CO, which increases the amount of carbon atoms in the gas phase.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

H2CS deuteration maps towards the pre-stellar core L1544

Deuteration is a crucial tool to understand the complexity of interstellar chemical processes, especially when they involve the interplay of gas-phase and grain-surface chemistry. In the case of multiple deuteration, comparing observation with the results of chemical modelling is particularly effective to study how molecules are inherited in the different stages within the process of star and planet formation. We aim to study the the D/ H ratio in H2CS across the prototypical pre-stellar core L1544. This study allows us to test current gas-dust chemical models involving sulfur in dense cores. We present here single-dish observations of H2CS, HDCS and D2CS with the IRAM 30m telescope. We analyse their column densities and distributions, and compare these observations with gas-grain chemical models. The deuteration maps of H2CS in L1544 are compared with the deuteration maps of methanol, H2CO, N2H+ and HCO+ towards the same source. Furthermore, the single and double deuteration of H2CS towards the dust peak of L1544 is compared with H2CO and c-C3H2. The difference between the deuteration of these molecules in L1544 is discussed and compared with the prediction of chemical models. The maximum deuterium fractionation for the first deuteration of H2CS is N(HDCS)/N(H2CS)$\sim$30$\%$ and is located towards the north-east at a distance of about 10000 AU from the dust peak. While for c-C3H2 the first and second deuteration have a similar efficiency, for H2CS and H2CO the second deuteration is more efficient, leading to D2CX/HDCX$\sim$100$\%$ (with X= O or S). Our results imply that the large deuteration of H2CO and H2CS observed in protostellar cores as well as in comets is likely inherited from the pre-stellar phase. However, the comparison with state-of-the-art chemical models suggests that the reaction network for the formation of the doubly deuterated H2CS and H2CO it is not complete yet.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effect of grain size distribution and size-dependent grain heating on molecular abundances in starless and pre-stellar cores

We present a new gas-grain chemical model to constrain the effect of grain size distribution on molecular abundances in starless and pre-stellar cores. We introduce grain-size dependence simultaneously for cosmic-ray (CR)-induced desorption efficiency and for grain equilibrium temperatures. We keep explicit track of ice abundances on a set of grain populations. We find that the size-dependent CR desorption efficiency affects ice abundances in a highly non-trivial way that depends on the molecule. Species that originate in the gas phase follow a simple pattern where the ice abundance is highest on the smallest grains (the most abundant in the distribution). Some molecules, such as HCN, are instead concentrated on large grains throughout the time evolution, while others (like $\rm N_2$) are initially concentrated on large grains, but at late times on small grains, due to grain-size-dependent competition between desorption and hydrogenation. Most of the water ice is on small grains at high medium density ($n({\rm H_2}) \gtrsim 10^6 \, \rm cm^{-3}$), where the water ice fraction, with respect to total water ice reservoir, can be as low as $\sim 10^{-3}$ on large (> 0.1 $μ$m) grains. Allowing the grain equilibrium temperature to vary with grain size induces strong variations in relative ice abundances in low-density conditions where the interstellar radiation field and in particular its ultraviolet component are not attenuated. Our study implies consequences not only for the initial formation of ices preceding the starless core stage, but also for the relative ice abundances on the grain populations going into the protostellar stage. In particular, if the smallest grains can lose their mantles due to grain-grain collisions as the core is collapsing, the ice composition in the beginning of the protostellar stage could be very different to that in the pre-collapse phase.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effect of multilayer ice chemistry on gas-phase deuteration in starless cores

Aims. We aim to investigate whether a multilayer ice model can be as successful as a bulk ice model in reproducing the observed abundances of various deuterated gas-phase species toward starless cores. Methods. We calculate abundances for various deuterated species as functions of time adopting fixed physical conditions. We also estimate abundance gradients by adopting a modified Bonnor-Ebert sphere as a core model. In the multilayer ice scenario, we consider desorption from one or several monolayers on the surface. Results. We find that the multilayer model predicts abundances of $\rm DCO^+$ and $\rm N_2D^+$ that are about an order of magnitude lower than observed, caused by the trapping of CO and $\rm N_2$ into the grain mantle. As a result of the mantle trapping, deuteration efficiency in the gas phase increases and we find stronger deuterium fractionation in ammonia than what has been observed. Another distinguishing feature of the multilayer model is that $\rm D_3^+$ becomes the main deuterated ion at high density. The bulk ice model is generally easily reconciled with observations. Conclusions. Our results underline that more theoretical and experimental work is needed to understand the composition and morphology of interstellar ices, and the desorption processes that can act on them. With the current constraints, the bulk ice model appears to be better in reproducing observations than the multilayer ice model. According to our results, the $\rm H_2D^+$ to $\rm N_2D^+$ abundance ratio is higher than 100 in the multilayer model, while only a few $\times$ 10 in the bulk model, and so observations of this ratio could provide information on the ice morphology in starless cores. Observations of the abundance of $\rm D_3^+$ compared to $\rm H_2D^+$ and $\rm D_2H^+$ would provide additional constraints for the models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The NH$_2$D hyperfine structure revealed by astrophysical observations

The 1$_{11}$-1$_{01}$ lines of ortho and para--NH$_2$D (o/p-NH$_2$D), respectively at 86 and 110 GHz, are commonly observed to provide constraints on the deuterium fractionation in the interstellar medium. In cold regions, the hyperfine structure due to the nitrogen ($^{14}$N) nucleus is resolved. To date, this splitting is the only one which is taken into account in the NH$_2$D column density estimates. We investigate how the inclusion of the hyperfine splitting caused by the deuterium (D) nucleus affects the analysis of the rotational lines of NH$_2$D. We present 30m IRAM observations of the above mentioned lines, as well as APEX o/p-NH$_2$D observations of the 1$_{01}$-0$_{00}$ lines at 333 GHz. The hyperfine spectra are first analyzed with a line list that only includes the hyperfine splitting due to the $^{14}$N nucleus. We find inconsistencies between the line widths of the 1$_{01}$-0$_{00}$ and 1$_{11}$-1$_{01}$ lines, the latter being larger by a factor of $\sim$1.6$\pm0.3$. Such a large difference is unexpected given the two sets of lines are likely to originate from the same region. We next employ a newly computed line list for the o/p-NH$_2$D transitions, where the hyperfine structure induced by both nitrogen and deuterium nuclei is included. With this new line list, the analysis of the previous spectra leads to linewidths which are compatible. Neglecting the hyperfine structure owing to D leads to overestimate the linewidths of the o/p-NH$_2$D lines at 3 mm. The error for a cold molecular core is about 50\%. This error propagates directly to the column density estimate. It is therefore recommended to take into account the hyperfine splittings caused by both the $^{14}$N and D nuclei in any analysis relying on these lines.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Understanding the $\rm C_3H_2$ cyclic-to-linear ratio in L1544

Aims. We aim to understand the high cyclic-to-linear $\rm C_3H_2$ ratio ($32 \pm 4$) observed toward L1544 by Spezzano et al. (2016). Methods. We combine a gas-grain chemical model with a physical model for L1544 to simulate the column densities of cyclic and linear $\rm C_3H_2$ observed toward L1544. The most important reactions for the formation and destruction of both forms of $\rm C_3H_2$ are identified, and their relative rate coefficients are varied to find the best match to the observations. Results. We find that the ratio of the rate coefficients of $\rm C_3H_3^+ + e^- \longrightarrow C_3H_2 + H$ for cyclic and linear $\rm C_3H_2$ must be $\sim 20$ in order to reproduce the observations, depending on the branching ratios assumed for the $\rm C_3H_3^+ + e^- \longrightarrow C_3H + H_2$ reaction. In current astrochemical networks it is assumed that cyclic and linear $\rm C_3H_2$ are formed in a 1:1 ratio in the aforementioned reactions. Laboratory studies and/or theoretical calculations are needed to confirm the results of our chemical modeling based on observational constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Benchmarking spin-state chemistry in starless core models

Aims. We aim to present simulated chemical abundance profiles for a variety of important species, with special attention given to spin-state chemistry, in order to provide reference results against which present and future models can be compared. Methods. We employ gas-phase and gas-grain models to investigate chemical abundances in physical conditions corresponding to starless cores. To this end, we have developed new chemical reaction sets for both gas-phase and grain-surface chemistry, including the deuterated forms of species with up to six atoms and the spin-state chemistry of light ions and of the species involved in the ammonia and water formation networks. The physical model is kept simple in order to facilitate straightforward benchmarking of other models against the results of this paper. Results. We find that the ortho/para ratios of ammonia and water are similar in both gas-phase and gas-grain models, at late times in particular, implying that the ratios are determined by gas-phase processes. We derive late-time ortho/para ratios of ~0.5 and ~1.6 for ammonia and water, respectively. We find that including or excluding deuterium in the calculations has little effect on the abundances of non-deuterated species and on the ortho/para ratios of ammonia and water, especially in gas-phase models where deuteration is naturally hindered owing to the presence of abundant heavy elements. Although we study a rather narrow temperature range (10-20 K), we find strong temperature dependence in, e.g., deuteration and nitrogen chemistry. For example, the depletion timescale of ammonia is significantly reduced when the temperature is increased from 10 to 20 K; this is because the increase in temperature translates into increased accretion rates, while the very high binding energy of ammonia prevents it from being desorbed at 20 K.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the stability of non-isothermal Bonnor-Ebert spheres. II. The effect of gas temperature on the stability

Aims. We investigate the stability of non-isothermal Bonnor-Ebert spheres with a model that includes a self-consistent calculation of the gas temperature. This way we can discard the assumption of equality between the dust and gas temperatures, and study the stability as the gas temperature changes with chemical evolution of the gas. Methods. We use a gas-grain chemical model including a time-dependent treatment of depletion onto grain surfaces, which strongly influences the gas temperature as the main coolant, CO, depletes from the gas. Dust and gas temperatures are solved with radiative transfer. For comparison with previous work, we assume that the cores are deeply embedded in a larger external structure, corresponding to visual extinction $A_{\rm V}^{\rm ext}=10$ mag. Results. We find that the critical non-dimensional radius $ξ_1$ derived here is similar to our previous work where we assumed $T_{\rm dust}=T_{\rm gas}$; the $ξ_1$ values lie below the isothermal critical value $ξ_0\sim6.45$, but the difference is less than 10%. Chemical evolution does not affect notably the stability condition of low-mass cores (<0.75 $M_\odot$). For higher masses the decrease of cooling owing to CO depletion causes substantial temporal changes in the temperature and density profiles of the cores. In the mass range 1-2 $M_\odot$ , $ξ_1$ decreases with chemical evolution, whereas above 3 $M_\odot$ , $ξ_1$ instead increases. We also find that decreasing $A_{\rm V}^{\rm ext}$ increases the gas temperature especially when the gas is chemically old, causing $ξ_1$ to increase with respect to models with higher $A_{\rm V}^{\rm ext}$. The derived $ξ_1$ values are close to $ξ_0$. The density contrast between the core center and edge varies between 8 to 16 depending on core mass and the chemical age of the gas, compared to the constant value $\sim$ 14.1 for the isothermal BES.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spin-state chemistry of deuterated ammonia

Aims. We aim to develop a chemical model that contains a consistent description of spin-state chemistry in reactions involving chemical species with multiple deuterons. We apply the model to the specific case of deuterated ammonia, to derive values for the various spin-state ratios. Methods. We apply symmetry rules in the complete scrambling assumption to calculate branching ratio tables for reactions between chemical species that include multiple protons and/or deuterons. Reaction sets for both gas-phase and grain-surface chemistry are generated using an automated routine that forms all possible spin-state variants of any given reaction with up to six H/D atoms. Single-point and modified Bonnor-Ebert models are used to study the density and temperature dependence of ammonia and its isotopologs, and the associated spin-state ratios. Results. We find that the spin-state ratios of the ammonia isotopologs are, at late times, very different from their statistical values. The ratios are rather insensitive to variations in the density, but present strong temperature dependence. We derive high peak values ($\sim$ 0.1) for the deuterium fraction in ammonia, in agreement with previous (gas-phase) models. The deuterium fractionation is strongest at high density, corresponding to a high degree of depletion, and also presents temperature dependence. We find that in the temperature range 5 to 20 K, the deuterium fractionation peaks at $\sim$ 15 K while most of the ortho/para (and meta/para for $\rm ND_3$) ratios present a minimum at 10 K (ortho/para $\rm NH_2D$ has instead a maximum at this temperature). Conclusions. Owing to the density and temperature dependence found in the abundances and spin-state ratios of ammonia and its isotopologs, it is evident that observations of ammonia and its deuterated forms can provide important constraints on the physical structure of molecular clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

HD depletion in starless cores

Aims: We aim to investigate the abundances of light deuterium-bearing species such as HD, H2D+ and D2H+ in a gas-grain chemical model including an extensive description of deuterium and spin state chemistry, in physical conditions appropriate to the very centers of starless cores. Methods: We combine a gas-grain chemical model with radiative transfer calculations to simulate density and temperature structure in starless cores. The chemical model includes deuterated forms of species with up to 4 atoms and the spin states of the light species H2, H2+ and H3+ and their deuterated forms. Results: We find that HD eventually depletes from the gas phase because deuterium is efficiently incorporated to grain-surface HDO, resulting in inefficient HD production on grains. HD depletion has consequences not only on the abundances of e.g. H2D+ and D2H+, whose production depends on the abundance of HD, but also on the spin state abundance ratios of the various light species, when compared with the complete depletion model where heavy elements do not influence the chemistry. Conclusions: While the eventual HD depletion leads to the disappearance of light deuterium-bearing species from the gas phase in a relatively short timescale at high density, we find that at late stages of core evolution the abundances of H2D+ and D2H+ increase toward the core edge and the disributions become extended. The HD depletion timescale increases if less oxygen is initially present in the gas phase, owing to chemical interaction between the gas and the dust predecing the starless core phase. Our results are greatly affected if H2 is allowed to tunnel on grain surfaces, and therefore more experimental data not only on tunneling but also on the O + H2 surface reaction in particular is needed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radial molecular abundances and gas cooling in starless cores

Aims: We aim to simulate radial profiles of molecular abundances and the gas temperature in cold and heavily shielded starless cores by combining chemical and radiative transfer models. Methods: A determination of the dust temperature in a modified Bonnor-Ebert sphere is used to calculate initial radial molecular abundance profiles. The abundances of selected cooling molecules corresponding to two different core ages are then extracted to determine the gas temperature at two time steps. The calculation is repeated in an iterative process yielding molecular abundances consistent with the gas temperature. Line emission profiles for selected substances are calculated using simulated abundance profiles. Results: The gas temperature is a function of time; the gas heats up as the core gets older because the cooling molecules are depleted onto grain surfaces. The contributions of the various cooling molecules to the total cooling power change with time. Radial chemical abundance profiles are non-trivial: different species present varying degrees of depletion and in some cases inward-increasing abundances profiles, even at t > 10^5 years. Line emission simulations indicate that cores of different ages can present significantly different line emission profiles, depending on the tracer species considered. Conclusions: Chemical abundances and the associated line cooling power change as a function of time. Most chemical species are depleted onto grain surfaces at densities exceeding ~10^5 cm^-3. Notable exceptions are NH_3 and N2H^+; the latter is largely undepleted even at n_H~10^6 cm-3. On the other hand, chemical abundances are not significantly developed in regions of low gas density even at t~10^5 years, revealed by inward-increasing abundance gradients. The gas temperature can be significantly different from the dust temperature; this may have implications on core stability.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the stability of non-isothermal Bonnor-Ebert spheres

Aims: We aim to derive a stability condition for non-isothermal Bonnor-Ebert spheres and compare the physical properties of critical non-isothermal and isothermal gas spheres. These configurations can serve as models for prestellar cores before gravitational collapse. Methods: A stability condition for non-isothermal spheres is derived by constructing an expression for the derivative of boundary pressure with respect to core volume. The temperature distribution is determined by means of radiative transfer calculations. Based on the stability analysis, we derive the physical parameters of critical cores for the mass range 0.1 - 5.0 M_sun. In addition, the properties of roughly Jupiter-mass cores are briefly examined. Results: At the low-mass end the critical non-isothermal sphere has lower central density and a slightly larger physical radius than the corresponding isothermal sphere (i.e. one with the same mass and average temperature). The temperature decrease towards the core centre becomes steeper towards smaller masses as the central density becomes higher. The slope depends on the adopted dust model. We find that the critical dimensionless radius increases above the isothermal value xi_0 = 6.45 for very low-mass cores (< 0.2 M_sun). However, in the mass-range studied here the changes are within 5% from the isothermal value. Conclusions: The density structures of non-isothermal and isothermal Bonnor-Ebert spheres for a given mass are fairly similar. However, the present models predict clear differences in the average temperatures for the same physical radius. Especially for low-mass cores, the temperature gradient probably has implications on the chemistry and the observed line emission. We also find that hydrostatic Jupiter-mass cores with radii less than 100 AU would have very high boundary pressures compared with typical pressures in the interstellar space.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modelling line emission of deuterated H_3^+ from prestellar cores

Context: The depletion of heavy elements in cold cores of interstellar molecular clouds can lead to a situation where deuterated forms of H_3^+ are the most useful spectroscopic probes of the physical conditions. Aims: The aim is to predict the observability of the rotational lines of H_2D^+ and D_2H^+ from prestellar cores. Methods: Recently derived rate coefficients for the H_3^+ + H_2 isotopic system were applied to the "complete depletion" reaction scheme to calculate abundance profiles in hydrostatic core models. The ground-state lines of H_2D^+(o) (372 GHz) and D_2H^+(p) (692 GHz) arising from these cores were simulated. The excitation of the rotational levels of these molecules was approximated by using the state-to-state coefficients for collisions with H_2. We also predicted line profiles from cores with a power-law density distribution advocated in some previous studies. Results: The new rate coefficients introduce some changes to the complete depletion model, but do not alter the general tendencies. One of the modifications with respect to the previous results is the increase of the D_3^+ abundance at the cost of other isotopologues. Furthermore, the present model predicts a lower H_2D^+ (o/p) ratio, and a slightly higher D_2H^+ (p/o) ratio in very cold, dense cores, as compared with previous modelling results. These nuclear spin ratios affect the detectability of the submm lines of H_2D^+(o) and D_2H^+(p). The previously detected H_2D^+ and D_2H^+ lines towards the core I16293E, and the H_2D^+ line observed towards Oph D can be reproduced using the present excitation model and the physical models suggested in the original papers.

0 citations0 reviewsNo code linkedOpen access