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

Jonathan Holdship

Jonathan Holdship contributes to research discovery and scholarly infrastructure.

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

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
3topics
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

6 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access