T. J.
Millar
*
Astrophysics Research Centre, Queen’s University Belfast, University Road, Belfast BT7 1NN, UK. E-mail: tom.millar@qub.ac.uk
First published on 7th July 2023
Fifty years on from the first detailed chemical kinetic modelling of astronomical sources, I provide some introductory comments on the history of astrochemistry, summarise some personal views on the topics covered in this discussion meeting, and conclude with some thoughts on its future development. I have left out the jokes.
The earliest attempt to describe the gas phase chemical processes that might operate in interstellar clouds was published by Bates and Spitzer2 in an effort to understand the CH and CH+ observations. With the exception of ion–neutral reactions, they identified most of the important chemical processes likely to operate in the diffuse environment in which these species were found. Indeed, finding that their calculated abundances did not reproduce the observations, they even considered the possibility that the thermal desorption of CH4 from grains might be their source. Of course, one should recognise that a major reason for model failures was their lack of knowledge of the physical properties of the diffuse clouds in which these molecules were found.
The first detailed numerical model of chemical kinetics in interstellar clouds was published by Herbst and Klemperer3 and is rightly regarded as a classic in our field. The authors included many of the processes at work and worried about some that still concern us today – the elemental abundance of oxygen, the carbon-to-oxygen (C/O) elemental abundance ratio, the role of grains as a source and sink of interstellar molecules, and the extent to which a global chemical kinetic model can be applied to individual molecular clouds. Although some arguments concerning time scales are somewhat dated, this is, I would argue, the foundational paper for our subject and one which led directly to the plethora of experimental, theoretical, gas–grain modelling approaches to astrochemistry presented at this meeting.
Starting around 1985, spectral line surveys of high-mass star-formation systems began to detect in a large number of objects, a plethora of interstellar complex organic molecules (iCOMs), very loosely defined as containing at least six atoms including one element in addition to carbon and hydrogen. Many of the sources were rather hot and dense by interstellar standards, ∼100–300 K, n(H2) ∼ 106–109 cm−3. These species, typified by their smallest member CH3OH, are often saturated and tend to have larger abundances, often by several orders of magnitude, in these sources compared to those in cold (10 K) clouds. Observations of their D-substituted analogues, however, showed a surprising result – molecular D/H abundance ratios that are much larger, by factors of 103–105, than the cosmic D/H ratio of a few 10−5. Indeed, the largest enhancement I am aware of occurs for triply-deuterated ammonia, ND3, which has an abundance relative to NH3 of about 10−3, some 1011 times larger than the statistical D/H ratio. Such large fractionation can only occur at very low temperatures and the accepted paradigm is that iCOMs and their deuterated versions are formed in molecular ices on cold dust grains before being desorbed into hot gas through the effect of stellar heating or shock waves. Such observations led to a number of models in which H atom collisions with a cold ice surface was assumed to convert atoms such as N, C, O and S to NH3, CH4, H2O and H2S in a water-rich ice mantle. Subsequent heating due to the presence of a newly-formed hot, massive star liberated large abundances of these simple hydrogenated molecules to the gas phase where subsequent gas-phase chemistry produced the iCOMs.4–6 The lack of any detailed description of the actual chemistry on these icy grains and the cry of the modellers that ‘all’ abundance anomalies could be explained by surface chemistry led Malcolm Walmsley to note that grain chemistry ‘was the last refuge of the scoundrel’.
Fig. 1 Number of molecules detected in interstellar and circumstellar environments as a function of the number of atoms. |
In the last five years there has been a substantial increase in the number of detections of complex hydrocarbons mainly driven through dedicated spectral line surveys of the cold, dark Taurus Molecular Cloud 1 (TMC-1). McGuire and colleagues have concentrated on using sophisticated analysis of rotational emission lines as part of their GOTHAM survey to identify a number of cyclic hydrocarbons including benzonitrile,7 1- and 2-cyanocyclopentadiene, C5H5CN,8,9 indene,10 and 2-cyanoindene.11 Cernicharo and co-workers, using the very sensitive NANOCOSMOS receiver, have detected many more hydrocarbons in their QUIJOTE survey of TMC-1 over the past few years, including several 5- and 6-membered ring species and branched chains. Interesting species include 1- and 2-ethynyl cyclopentadiene (c-C5H5CCH), ethynyl benzene (C6H5CCH),12 and ethynyl cyclopropenylidene (c-C3HCCH).13
This discussion meeting was convened at a very appropriate time – and with a very timely agenda – given the recent development of very sensitive facilities such as the JWST (James Webb Space Telescope) with both higher sensitivity and spectral resolution than any other IR space observatory; as well as high spatial resolution instruments such as the ALMA (Atacama Large Millimeter/submillimeter Array). The early results from the JWST are breath-taking with an impact that affects every discipline in astronomy. The ALMA, which has been in operation for a decade, also continues to offer new views of the molecular universe and many of the observational papers at this meeting reflected the current dominant positions of these two facilities.
The meeting was arranged around four main themes: (1) observational astrochemistry; (2) laboratory gas-phase studies; (3) experimental studies of dust and ices; (4) computational astrochemistry. A number of poster papers were also on display during the meeting. My comments on the presentations are not meant to be comprehensive but reflect a number of personal factors including a lack of time to assimilate everything in the papers, presentations, extensive discussions and posters, a subjective bias toward certain topics, and a degree of ignorance in some areas. Nevertheless, I do commend the proceedings to the reader; they contain a wealth of information on both the rude health and the underlying concerns – the known unknowns – of our subject and its importance in answering fundamental questions in astrophysics, including the formation of stars and planetary systems. In particular, Cecilia Ceccarelli’s Spiers Memorial Lecture (https://doi.org/10.1039/d3fd00106g), is a masterly introduction to the field, outlining recent successes, identifying complexities yet to be addressed in a comprehensive manner, and pointing to the future. For my part, the most important message of the meeting is that the foundation of astrochemistry is spectroscopy. Without the decades of advances in laboratory spectroscopy, and especially in astronomical instrumentation, our exploration of chemical processes in space would not be ‘much more complex’ than those considered 50 years ago.
Kamp derived similar results from fitting a combination of slab and 2D thermochemical disk models to her observed spectra. She noted that line lists for 13C12CH2 were incomplete, preventing a full analysis of the spectra and of the (likely large) 13C/12C ratio. Again differences within and between source types were evident. All T Tauri disks contain H2O whereas it is only tentatively detected in the low mass disks. She also stressed the need to go beyond the slab retrieval models which treat only one molecule at a time and assume constant density and temperature. In addition, she noted that the underlying molecular line data and partition functions needed for such analyses are often missing. Her comparison between column densities derived with both slab and thermochemical models show that the latter typically predict values a factor of ten larger than the former. Clearly a more accurate approach is needed but again it was stressed that databases currently lack fundamental information – Einstein A-coefficients, partition functions at high T and P – as well as spectroscopy to provide hot band and isotopologue line lists for hydrocarbons. Finally, it is worth re-iterating that these are the early results from the JWST; we can expect a host of new species and new insights into the chemistry and physics of these planet-forming systems in the coming years. As always with new facilities, it will be a challenge for both laboratory and theoretical approaches to keep abreast of new observational developments.
Three papers presented high angular resolution ALMA observations on spatial scales ranging from the Solar System to 105 times larger. On the largest scale, Huang’s observations (https://doi.org/10.1039/d3fd00007a) used molecular tracers of shock waves to probe the similarities and differences between molecular gas in two galaxies – NGC 253, a starburst galaxy, and NGC 1068, an active galactic nucleus harbouring a supermassive black hole which generates high X-ray and cosmic-ray fluxes. Detailed observations were interpreted by models of shock chemistry finding evidence for the presence of both slow shocks, via observations of HNCO, and fast shocks, via the detection of large abundances of SiO which arise from the sputtering of silicate dust grains. Wilkins and Blake (https://doi.org/10.1039/d3fd00003f) presented a pilot spectral line search in the 1.2 mm window toward 11 giant molecular clouds (GMCs) at a distance of 4–8 kpc from the Galactic Centre, while Bianchi (https://doi.org/10.1039/d3fd00018d) investigated molecular line emission in SVS13-A, proto-binary system, in particular from the streamers connecting the binary system to their larger scale molecular envelope. This source has a rich array of molecular lines generated by iCOMs in the hot gas in the innermost region. Her study of HDO and SO2 found a large, warm (95–160 K), column density of the former, indicating that it has been released from dust grain mantles, most likely by shock processes generated in the infalling gas. Astrochemistry and the astrophysics that can be derived from it are likely to be dominated by the ALMA and the JWST for the next decade or two – an exciting prospect.
The remainder of the session was devoted to spectroscopy and the study of transient species. This included papers by Brünken (https://doi.org/10.1039/d3fd00015j) on the fragmentation products of cyclic ions; Gupta who used a cryogenic ion trap at 4 K to record the rotational spectrum of c-C3H2D (https://doi.org/10.1039/d3fd00068k); and Jacovella (https://doi.org/10.1039/d3fd00016h) who measured the rotational spectrum of the cross-branched molecule norbornadiene (C8H7), and its cyano-derivatives, which have much larger electric dipole moments, as a pre-cursor to an (unsuccessful) astronomical search. Puzzarini (https://doi.org/10.1039/d3fd00052d) presented a strategy for determining the rotational spectra of transient molecules using quantum calculations to predict properties and thereby guide laboratory searches and line identification.
‘Nothing puzzles me more than molecules and ices yet nothing bothers me less as I never think about them’.
Nevertheless it is the job of the astrochemist to both worry and puzzle over the known and unknown unknowns of ice chemistry.
Four major questions on the cold (∼10 K) ice chemistry that is thought to produce iCOMs, continue to exercise the community: what process removes these species from grains at temperatures of ∼6–10 K to the cold gas phase in which they are detected?; do these processes return whole or fragmented species?; what are the binding energies of surface iCOMs since these determine the efficiency of thermal desorption processes?; and, given the fact that experiments and theory are limited to relatively few species, can we generalise results from the few to the many? A number of answers to these questions have been suggested in the literature and some were covered here.
Brown (https://doi.org/10.1039/d3fd00024a) discussed experiments which showed that infrared radiation could desorb CO from CO–water ice mixtures, noting that CO didn’t desorb when irradiated with 4.6 μm photons but did so when the irradiation was at the wavelengths of the H2O vibrational modes at 2.9 and 12 μm. If a single photon process and if the efficiency is high, this could be an important desorption mechanism in dark clouds which are much more transparent to IR radiation than to UV photons, which have long been thought to drive photodesorption. Brown identified two mechanisms: a fast process of indirect resonant desorption in which an excited H2O molecule removes a neighbouring CO molecule, and a slow, photon-induced process in which energy is accumulated in the heat bath of water ice. It should be noted that the water ice layers are thin enough that the IR photons also heat the cold, metal substrate and that this may also contribute to desorption. Although poorly defined, Brown suggested yields of CO molecules removed per incident photon to be on the order of a few 10−5.
Like Brown, Bertin (https://doi.org/10.1039/d3fd00004d) studied the photodesorption of ices, in his case with vacuum UV (VUV) photons and ices composed of either pure formic acid (HCOOH), or methyl formate (HCOOCH3), or with these species embedded in CO or H2O ices.
In addition to determining desorption yields, around 10−5 per incident VUV photon, another key issue is whether molecules are returned intact or as fragments. Bertin found similar results to previous experiments on methanol-containing ice, namely, that species detected in the gas phase are similar to the gas-phase photodesorption products of the molecule. Were this to be a universal rule it would be very convenient for modellers, but it is not – recent experiments on acetonitrile, CH3CN, show that it photodesorbed intact.16
Gudipati (https://doi.org/10.1039/d3fd00048f) discussed an alternative desorption mechanism, one driven by the conversion of water ice from the amorphous to crystalline form. His experiments doped water ice with the volatiles CO, CO2 and O2 at 10 K and traced the outgassing species via mass spectrometry as the ice was heated to about 200 K. He showed that their abundance in the ice reduce by about a factor of 100 due to outgassing when the phase change occurs.
Finally, Sameera (https://doi.org/10.1039/d3fd00033h) discussed the important issues of binding energies and reaction barriers on ices. It has been known for some time, and certainly since the pioneering experimental work of Collings et al.,17 that the binding energies of a molecule in water ice depends on the specific site in which it sits. More recently, Tinacci et al.18 used theoretical calculations of NH3 on a model cluster containing 100 H2O molecules to show that binding energies followed a double peaked profile with that at higher energies well fit by experimental TPD measurements. Rather surprisingly, however, there was a second distribution with much lower binding energies comparable to those of CO.
Sameera concentrated on CH3OH in water ice using a theoretical model to consider the reactions:
CH3OH + OH → CH2OH + H2O | (1) |
CH3OH + OH → CH3O + H2O | (2) |
Many of the key issues raised here, including the effects of binding energies on grain chemistry, are discussed at length in Ceccarelli’s Spiers Memorial Lecture (https://doi.org/10.1039/d3fd00106g).
Garrod and Herbst (https://doi.org/10.1039/d3fd00014a) presented a very detailed time-dependent model of iCOM formation in molecular clouds based on their identification that several detected and candidate interstellar molecules containing the amine group, –NH2, have proton affinities (PA) that are greater than that of NH3, 853.6 kJ mol−1. Because of its relatively large abundance in molecular clouds and its large PA, astrochemical models tend to include proton transfer reactions between ammonia and molecular ions as an important step in forming the neutral molecule. An example is:
H2COH+ + NH3 → NH4+ + H2CO | (3) |
The advantage of this proton transfer reaction is that it preserves the underlying neutral molecule, whereas the alternative loss of molecular ions, dissociative recombination with electrons, can destroy the molecule in a large fraction of cases. In the case of protonated formaldehyde, reaction with electrons has five product channels with that leading to formaldehyde occurring only around 30% of the time.19 For protonated methanol, the situation is much worse with only 5% of recombinations re-cycling the neutral.20 As a result, stable neutral molecules with PAs greater than that of ammonia are destroyed faster than those with PAs less than ammonia. Garrod’s paper shows that this destruction process has a particular affect on the abundances of many amine species, including methylamine, ethylamine, urea and glycine. The latter molecule, yet to be identified in interstellar clouds, is a cornerstone in the argument that the precursors of life on Earth may be widely dispersed in the Galaxy – Garrod’s results suggest that it may be more difficult to detect glycine than we currently envisage. One should note, however, that for reactions in which the difference between PAs is large, most importantly for reactions involving H3+, the excess reaction energy can lead to dissociative products rather than proton transfer, although not in every case.21
Viti (https://doi.org/10.1039/d3fd00008g) and Van de Sande (https://doi.org/10.1039/d3fd00039g) discussed novel ways of approaching the gas–grain modelling of molecular cloud and circumstellar chemistry, respectively. The former presented an alternative approach of using machine learning to model some of the uncertainties in complex chemical kinetic networks while the latter used a sophisticated physical model to discuss the effects of clumpy structures in the expanding envelopes around late-type stars. Viti’s paper used both Bayesian and machine learning approaches to identify both those binding energies with the largest impact on calculated abundances and those species for which determining binding energies are a priority. Van de Sande’s paper includes the words ‘disentangling’ and ‘when adding complexity’ in its title – a bold claim. As a co-author, I recuse myself at this point and leave it as an exercise for the reader as to whether the content lives up to its title!
Bromley (https://doi.org/10.1039/d3fd00055a) addressed the important topic of dust nucleation in the context of arguments – in essence that the dust destruction time scale is shorter than the dust formation time scale in evolved stars and supernovae – that a significant fraction of silicate dust must nucleate in the low density interstellar medium. Bromley concentrated on the creation of nanosilicates, combining quantum calculations on the lowest energy structures of small seed particles with their infrared absorption characteristics. He showed that the JWST could be sensitive enough to detect such particles if they composed 5–10% of the silicate grain population.
Finally, Madhusudhan (https://doi.org/10.1039/d3fd00075c) described the atmospheric physics and chemistry on a Hycean world, that is a sub-Neptunian exoplanet whose surface is covered by water oceans and with a H2-rich atmosphere. Although no such planet has yet been detected, several candidate systems have been identified orbiting M dwarf stars. In addition to evolution time scales of several Gyr, Madhusudhan’s paper showed one of the important differences between exoplanetary chemistry and astrochemistry, namely that the degree of physical complexity included in the former is of a much higher order than in the latter. In particular, the atmospheric models, internal planetary structures and delivery of exogenous materials are all informed by detailed observations and models of planet Earth and other solar-system bodies. Of course major uncertainties remain but he showed that the three essential ingredients for life – an energy source, liquid water and bio-essential elements – are all present in such worlds and that the JWST should be able to detect specific atmospheric molecules that are characteristic of Hycean worlds, an exciting prospect.
Advances in chemical kinetics have also been dramatic over the past 50 years. Many important low-temperature experimental studies are now routinely supplemented by theoretical calculations of potential energy surfaces and reaction rate coefficients and product branching ratios. Such research, built on advanced computational techniques and computing power, is now routine and, although expensive and time-consuming, it can answer important questions in astrochemistry.
In Section 2, I noted that in recent years, observational searches coupled with high resolution spectroscopy have identified more than 100 iCOMs and complex hydrocarbon species in interstellar clouds. Their observed abundances have led to innovative chemical models including both gas-phase and surface reactions. For example, Garrod’s basic model, as described here, has over 22000 reactions among 1385 species of which roughly half are in the gas phase with the rest equally divided between surface ice and bulk ice. Yet for many species, and as we heard in the discussions, abundance calculations rest on ‘thin ice’ (pardon the pun).
Let me use, as an example, some results from my own, as yet unpublished work,22 discovered while updating the UMIST Database for Astrochemistry chemical network, Rate12.23 Propene (CH3CHCH2), is a key intermediary in the synthesis of large hydrocarbons. It was detected in TMC-1 with a relatively large fractional abundance, 4 × 10−9, by Marcelino et al.24 and was thought to form via two sequential radiative association reactions:
C3H3+ + H2 → C3H5+ + hν | (4) |
C3H5+ + H2 → C3H7+ + hν | (5) |
Fig. 2 The fractional abundance of propene relative to H2 as a function of time for a model with properties similar to TMC-1 (T = 10 K, n(H2) = 104 cm−3). |
This fall in its calculated abundance by about 10 orders of magnitude leads to a ‘propene catastrophe’ in the models since it is known experimentally to be a reactive molecule even at very low temperatures. For example, it has fast measured rate coefficients in reactions with CH (15 K),26 CN (23 K),27 C2 (77 K),28 C2H (15 K)29 and C4H (39 K),30 where the value in parentheses refers to the lowest laboratory temperature at which the rate coefficient has been measured. The products of these reactions include many of the hydrocarbons and their cyano-derivatives that have been detected in TMC-1, including c-C5H6, CH3C4H, H2CCCHCCH, CH2CHCHCH2, CH3C6H, CH3CHCHCN, CH2C(CH3)CN and CH2CHCH2CN. Their calculated abundances are now more than 104 times less than observed – and the misery goes on, since these species are also involved in the synthesis of other complex molecules including c-C5H5CN, c-C5H4CCH2, c-C5H5CCH and H2CCCHC3N. Identifying efficient, low-temperature formation pathways to propene is one of the most important tasks for astrochemistry today. Such ‘single-point failures’ also exist for other species in the networks. Indeed it is likely that somewhere in the region of 70–100 detected molecules are not yet included in even the most complex networks because their synthesis – in the gas phase or in ices – is unknown.
To conclude, let me mention four generic areas that we will need to cultivate in the coming years if we are to ensure that the field of astrochemistry continues both its remarkable growth over the past 10–20 years and its relevance to astrophysics far and near. The first is synergy. To me, one of the central messages of this discussion meeting – funded by the Royal Society of Chemistry – is that research in core theoretical and experimental physical chemistry disciplines, including spectroscopy and gas-phase kinetics, plays a central role in astrochemistry. Collaborations between observers and modellers within the astronomy community were evident in the papers describing JWST and ALMA results and noticeably increase their impact.
Synergy across disciplines was also evident in and central to, many of the presentations and posters in this meeting and will increase in importance. For example, it is clear that many of the assumptions that underlie the grain-surface chemical kinetic models are too simplistic. Close co-operation between experimental and theoretical approaches on topics such as binding energies, surface and bulk diffusion, photodesorption yields, reaction products and desorption mechanisms, will be needed to ensure that theoretical models of interstellar chemistry become more predictive rather than reactive in nature. In a similar fashion, the astrochemical community needs to better identify critical reactions and processes for the chemical physics community to study.
The second is strategy. These synergistic links between chemistry and astronomy are often the result of individual actions, that is, where one researcher makes a measurement or calculation and reaches out to someone in the other field. There is an element of chance about these and indeed such interactions, while useful to a certain degree, may not have as much impact as one would like. We therefore need to find ways in which we can support research in a much more strategic manner. In the Netherlands, the Dutch Astrochemistry Network has been funded since 2010 and has supported many cross-disciplinary projects and, just as importantly, people. It not only brings together expertise in chemistry and astronomy but it also allows for the exchange of PhD students between institutions.
Such an approach is problematic in the UK where chemistry and astrophysics are funded by two separate research councils and where, in addition, many of the EU’s cross-disciplinary programmes are ever more difficult to access. In the UK the joint RSC/RAS body, the Astrophysical Chemistry Group, provides a forum – but no funding – through which interested parties can interact.
There are major issues facing many of us in astrochemistry – can the astronomical community help shield experimental groups, often expensive to run, from the growing commodification of modern university operations, where the bottom line is, too often, profit? This is not a one-way support of course. Chemistry in its broadest sense is universal and as we move further into the era of the ALMA and the JWST, to more detailed explorations of exoplanet atmospheres and the planets, moons and minor bodies of our solar system, astronomers will need the chemical disciplines to help interpret and understand observational results.
The third area is succession. In my view, the paper by Herbst and Klemperer,3 published 50 years ago, represents the birth of astrochemistry as a scientific discipline. Many of the scientists involved in the development of the foundations of the subject in its first 20 years, have retired or are close to retirement. They/we need to be active in planning their/our replacements. This means getting involved in departmental/university politics to ensure that the best people are recruited to, and promoted within, our discipline.
The final area is students. Our view of the molecular universe has improved enormously in the past 50 years. The development of high-resolution facilities, such as discussed here, provides us with an amazing panoply of images and spectra, coupled with a recognition of astrochemistry as central to answering many fundamental questions in astronomy, from the formation and evolution of galaxies to the formation of stars and planets. It was very noticeable at this meeting that the field has a wonderful array of talented, smart people who already possess a mature, informed, view of their research. The quality and impact of presentations and posters by PhD students and postdoctoral researchers was impressive as was their engagement in the discussion sessions. We need to ensure that, collectively, we harness those skills and encourage, excite, and enthuse our early career scientists. Where we can, let us support their intellectual growth, their exchange among institutions and across national boundaries, and let us find ways to employ them in full-time positions.
This journal is © The Royal Society of Chemistry 2023 |