BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

J. Harju

J. Harju contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

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
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
preprint2015arXiv

The physical state of selected cold clumps

Context. he study of prestellar cores is essential to understanding the initial stages of star formation. With $Herschel$ more cold clumps have been detected than ever before. For this study we have selected 21 cold clumps from 20 $Herschel$ fields observed as a follow-up on original $Planck$ detections. We have observed these clumps in $^{13}$CO (1-0), C$^{18}$O (1-0), and N$_2$H$^+$ (1-0) lines. Aims. Our aim is to find out if these cold clumps are prestellar. We have examined to what extent independent analysis of the dust and the molecular lines lead to similar conclusions about the masses of these objects. Methods. We calculate the clump masses and densities from the dust continuum and molecular line observations and compare these to each other and to the virial and Bonnor-Ebert masses calculated for each clump. Finally we examine two of the fields with radiative transfer models to estimate CO abundances. Results. When excitation temperatures could be estimated, the column densities derived from molecular line observations were comparable to those from dust continuum data. The median column density estimates are 4.2$\times 10^{21}$cm$^{-2}$ and 5.5$\times 10^{21}$cm$^{-2}$ for the line and dust emission data, respectively. The calculated abundances, column densities, volume densities, and masses all have large uncertainties and one must be careful when drawing conclusions. Abundance of $^{13}$CO was found in modeling the two clumps in the field G131.65$+$9.75 to be close to the usual value of 10$^{-6}$. The abundance ratio of $^{13}$CO and C$^{18}$O was $\sim$10. Molecular abundances could only be estimated with modeling, relying on dust column density data. Conclusions. The results indicate that most cold clumps, even those with dust color temperatures close to 11 K, are not necessarily prestellar.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Chemical and Physical Conditions in Molecular Cloud Core DC 000.4-19.5 (SL42) in Corona Australis

Chemical reactions in starless molecular clouds are heavily dependent on interactions between gas phase material and solid phase dust and ices. We have observed the abundance and distribution of molecular gases in the cold, starless core DC 000.4-19.5 (SL42) in Corona Australis using data from the Swedish ESO Submillimeter Telescope. We present column density maps determined from measurements of C18O(J=2-1,1-0) and N2H+(J=1-0) emission features. Herschel data of the same region allow a direct comparison to the dust component of the cloud core and provide evidence for gas phase depletion of CO at the highest extinctions. The dust color emperature in the core calculated from Herschel maps ranges from roughly 10.7 to 14.0 K. This range agrees with the previous determinations from Infrared Space Observatory and Planck observations. The column density profile of the core can be fitted with a Plummer-like density distribution approaching n(r) ~ r^{-2} at large distances. The core structure deviates clearly from a critical Bonnor-Ebert sphere. Instead, the core appears to be gravitationally bound and to lack thermal and turbulent support against the pressure of the surrounding low-density material: it may therefore be in the process of slow contraction. We test two chemical models and find that a steady-state depletion model agrees with the observed C18O column density profile and the observed N(C18O) versus AV relationship.

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

Reliability of NH3 as the temperature probe of cold cloud cores

The temperature is a central parameter affecting the chemical and physical properties of dense cores of interstellar clouds and their evolution to star formation. The chemistry and the dust properties are temperature dependent and the interpretation of observation requires the knowledge of the temperature and its variations. Measurement of the gas kinetic temperature is possible with molecular line spectroscopy, the ammonia molecule, NH3, being the most commonly used tracer. We want to determine the accuracy of the temperature estimates derived from ammonia spectra. The normal interpretation of NH3 observations assumes that all the hyperfine line components are tracing the same gas volume. In the case of temperature gradients they may be sensitive to different layers and cause errors in the optical depth and gas temperature estimates. We examine a series of spherical cloud models, 1.0 and 0.5 M_Sun Bonnor-Ebert spheres, with different radial temperature profiles. We calculate synthetic NH3 spectra and compare the derived column densities and temperatures to the true values. For high signal-to-noise observations, the estimated gas kinetic temperatures are within ~0.3 K of the real mass averaged temperature and the column densities are correct to within ~10%. When the S/N ratio of the (2,2) spectrum decreases below 10, the temperature errors are of the order of 1K but without a significant bias. When the density of the models is increased by a factor of a few, the results begin to show significant bias because of the saturation of the (1,1) main group. The ammonia spectra are found to be a reliable tracer of the mass averaged gas temperature. Because the radial temperature profiles of the cores are not well constrained, the central temperature could still differ from this value. If the cores are optically very thick, there are no guarantees of the accuracy.

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 &#34;complete depletion&#34; 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
preprint2003arXiv

Centimetre continuum emission from young stellar objects in Cederblad 110

The low-mass star formation region associated with the reflection nebula Cederblad 110 in the Chamaeleon I cloud was mapped with the Australian Telescope Compact Array (ATCA) at 6 and 3.5cm. Altogether 11 sources were detected, three of which are previously known low mass young stellar objects associated with the nebula: the illuminating star IRS2 (Class III, Einstein X-ray source CHX7), the brightest far-infrared source IRS4 (Class I), and the weak X-ray source CHX10a (Class III). The other young stellar objects in the region, including the Class 0 protostar candidate Cha-MMS1, were not detected. The radio spectral index of IRS4 (alpha = 1.7 +/- 0.3) is consistent with optically thick free-free emission arising from a dense ionized region, probably a jet-induced shock occurring in the circumstellar material. As the only Class I protostar with a &#39;thermal jet&#39; IRS4 is the strongest candidate for the central source of the molecular outflow found previously in the region. The emission from IRS2 has a flat spectrum (alpha = 0.05 +/- 0.05) but shows no sign of polarization, and therefore its origin is likely to be optically thin free-free emission either from ionized wind or a collimated jet. The strongest source detected in this survey is a new compact object with a steep negative spectral index (-1.1) and a weak linear polarization (about 2 %), which probably represents a background radio galaxy.

0 citations0 reviewsNo code linkedOpen access