Researcher profile

Jonathan Tennyson

Jonathan Tennyson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

28 published item(s)

preprint2026arXiv

ExoMol line lists -- LXIII: ExoMol line lists for 12 isotopologues of CO$_2$

Extensive rovibrational line lists are constructed for 12 isotopologues of carbon dioxide: $^{12}$C$^{16}$O$_2$, $^{13}$C$^{16}$O$_2$, $^{12}$C$^{17}$O$_2$, $^{13}$C$^{17}$O$_2$, $^{12}$C$^{18}$O$_2$, $^{13}$C$^{18}$O$_2$, $^{16}$O$^{12}$C$^{17}$O, $^{16}$O$^{12}$C$^{18}$O, $^{16}$O$^{13}$C$^{17}$O, $^{16}$O$^{13}$C$^{18}$O, $^{17}$O$^{12}$C$^{18}$O, and $^{17}$O$^{13}$C$^{18}$O. The variational program TROVE was employed together with an exact kinetic energy operator, accurate empirical potential energy surface (Ames-2) and the ab initio dipole moment surface Ames-2021-40K. Empirical energy levels from the most recent MARVEL analyses, as well as from the HITRAN and CDSD databases, are used to replace calculated values where available. The line lists are further supplemented by assigning AFGL quantum numbers using machine-learning based estimators. The resulting data were employed to generate opacities with four radiative transfer codes, TauREx, ARCiS, NEMESIS, and petitRADTRANS, both for individual isotopologues and for CO$_2$ at terrestrial isotopic natural abundance. All line lists and associated data are available at www.exomol.com.

preprint2025arXiv

ExoAtom: A Database of Atomic Spectra in ExoMol Format

We present the ExoAtom database, www.exomol.com/exoatom, an extension of the ExoMol database to provide atomic line lists in the ExoMol format. ExoAtom is designed for detailed astrophysical, planetary, and laboratory applications. ExoAtom currently includes atomic data for 80 neutral atoms and 74 singly charged ions. These data are extracted from both the NIST and Kurucz databases, with 79/71 atoms/ions sourced from NIST and 38/37 atoms/ions sourced from Kurucz. ExoAtom uses the file types .all, .def, .states, .trans and .pf as fundamental components for structuring atomic data in a consistent hierarchy. The .states file contains quantum numbers, uncertainties, lifetimes, etc. The .trans file specifies Einstein A coefficients and their associated wavenumbers. The .pf file provides partition functions over a wide grid of temperatures. Post-processing of the ExoAtom data is provided by the program PyExoCross. Future development of ExoAtom will include additional ionization stages.

preprint2024arXiv

The impact of spectral line wing cut-off: Recommended standard method with application to MAESTRO opacity database

When computing cross-sections from a line list, the result depends not only on the line strength, but also the line shape, pressure-broadening parameters, and line wing cut-off (i.e., the maximum distance calculated from each line centre). Pressure-broadening can be described using the Lorentz lineshape, but it is known to not represent the true absorption in the far wings. Both theory and experiment have shown that far from the line centre, non-Lorentzian behaviour controls the shape of the wings and the Lorentz lineshape fails to accurately characterize the absorption, leading to an underestimation or overestimation of the opacity continuum depending on the molecular species involved. The line wing cut-off is an often overlooked parameter when calculating absorption cross sections, but can have a significant effect on the appearance of the spectrum since it dictates the extent of the line wing that contributes to the calculation either side of every line centre. Therefore, when used to analyse exoplanet and brown dwarf spectra, an inaccurate choice for the line wing cut-off can result in errors in the opacity continuum, which propagate into the modeled transit spectra, and ultimately impact/bias the interpretation of observational spectra, and the derived composition and thermal structure. Here, we examine the different methods commonly utilized to calculate the wing cut-off and propose a standard practice procedure (i.e., absolute value of 25~cm$^{-1}$ for $P\leqslant$~200~bar and 100~cm$^{-1}$ for $P >$ ~200~bar) to generate molecular opacities which will be used by the open-access {\tt MAESTRO} (Molecules and Atoms in Exoplanet Science: Tools and Resources for Opacities) database. The pressing need for new measurements and theoretical studies of the far-wings is highlighted.

preprint2022arXiv

A method for the variational calculation of hyperfine-resolved rovibronic spectra of diatomic molecules

An algorithm for the calculation of hyperfine structure and spectra of diatomic molecules based on the variational nuclear motion is presented. Hyperfine coupling terms considered are Fermi-contact, nuclear spin-electron spin dipole-dipole, nuclear spin-orbit, nuclear spin-rotation and nuclear electric quadrupole interactions. Initial hyperfine-unresolved wavefunctions are obtained for given set of potential energy curves and associated couplings by variation solution of the nuclear-motion Schrödinger equation. Fully hyperfine-resolved parity-conserved rovibronic Hamiltonian matrices for a given final angular momentum, $\bm{F}$, are constructed and then diagonalized to give hyperfine-resolved energies and wavefunctions.Electric transition dipole moment curves can then be used to generate a hyperfine-resolved line list by applying rigorous selection rules.The algorithm is implemented in \textsc{Duo}, which is a general program for calculating spectra of diatomic molecules. This approach is tested for NO and MgH,and the results are compared to experiment and shown to be consistent with those given the well-used effective Hamiltonian code PGOPHER.

preprint2022arXiv

Cross-sections for heavy atmospheres: H$_2$O continuum

Most of the exoplanets detected up to now transit in front of their host stars, allowing for the generation of transmission spectra; the study of exoplanet atmospheres relies heavily upon accurate analysis of these spectra. Recent discoveries mean that the study of atmospheric signals from low-mass, temperate worlds are becoming increasingly common. The observed transit depth in these planets is small and more difficult to analyze. Analysis of simulated transmission spectra for two small, temperate planets (GJ 1214 b and K2-18 b) is presented, giving evidence for significant differences in simulated transit depth when the water vapor continuum is accounted for when compared to models omitting it. These models use cross-sections from the CAVIAR lab experiment for the water self-continuum up to 10,000 cm$^{-1}$; these cross-sections exhibit an inverse relationship with temperature, hence lower-temperature atmospheres are the most significantly impacted. Including the water continuum strongly affects transit depths, increasing values by up to 80 ppm, with the differences for both planets being detectable with the future space missions Ariel and JWST. It is imperative that models of exoplanet spectra move toward adaptive cross-sections, increasingly optimized for H$_2$O-rich atmospheres. This necessitates including absorption contribution from the water vapor continuum into atmospheric simulations.

preprint2022arXiv

Cross-sections for heavy atmospheres: H$_2$O self-broadening

The discovery of super-Earth and mini-Neptune exoplanets means that atmospheric signals from low-mass, temperate exoplanets are being increasingly studied. The signal acquired as the planet transits its host star, known as the transit depth, is smaller for these planets and, as such, more difficult to analyze. The launch of the space telescopes James Webb (JWST) & Ariel will give rise to an explosion in the quality and quantity of spectroscopic data available for an unprecedented number of exoplanets in our galaxy. Accurately extracting the information content, thereby permitting atmospheric science, of such data-sets will require robust models and techniques. We present here the analysis of simulated transmission spectra for water-rich atmospheres, giving evidence for non-negligible differences in simulated transit depths when self-broadening of H$_2$O is correctly accounted for, compared with the currently typically accepted standard of using H$_2$ and He-broadened cross-sections. Our case-study analysis is carried out on two super-Earths, focusing on water-based atmospheres, ranging from H$_2$-rich to H$_2$O-rich. The transit depth is considerably affected, increasing values by up to 60 ppm, which is shown to be detectable with JWST and Ariel. The differences are most pronounced for the lighter (i.e. $mu$ $\sim$ 4) atmospheres. Our work illustrates that it is imperative that the field of exoplanet spectroscopy moves toward adapted cross-sections, increasingly optimized for high-mu atmospheres for studies of super-Earths and mini-Neptunes.

preprint2022arXiv

ExoMol photodissociation cross sections I: HCl and HF

Photon initiated chemistry, \textit{i.e.} the interaction of light with chemical species, is a key factor in the evolution of the atmosphere of exoplanets. For planets orbiting stars in UV-rich environments, photodissociation induced by high energy photons dominates the atmosphere composition and dynamics. The rate of photodissociation can be highly dependent on atmospheric temperature, as increased temperature leads to increased population of vibrational excited states and the consequent lowering of the photodissociation threshold. This paper inaugurates a new series of papers presenting computed temperature-dependent photodissociation cross sections with rates generated for different stellar fields. Cross sections calculations are performed by solving the time-independent Schrödinger equation for each electronic state involved in the process. Here photodissociation cross sections for hydrogen chloride and hydrogen fluoride are computed for a grid of 34 temperatures between 0 and 10~000 K. Use of different radiation fields shows that for the Sun and cooler stars the photodissociation rate can increase exponentially for molecular temperatures above 1000 K; conversely the photodissociation rates in UV rich fields instead are almost insensitive to the temperature of the molecule. Furthermore, these rates show extreme sensitivity to the radiation model used for cool stars, suggesting that further work on these may be required. The provision of an ExoMol database of cross sections is discussed.

preprint2022arXiv

Original Research By Young Twinkle Students (ORBYTS): Ephemeris Refinement of Transiting Exoplanets III

We report photometric follow-up observations of thirteen exoplanets (HATS-1 b, HATS2 b, HATS-3 b, HAT-P-18 b, HAT-P-27 b, HAT-P-30 b, HAT-P-55 b, KELT-4A b, WASP-25 b, WASP-42 b, WASP-57 b, WASP-61 b and WASP-123 b), as part of the Original Research By Young Twinkle Students (ORBYTS) programme. All these planets are potentially viable targets for atmospheric characterisation and our data, which were taken using the LCOGT network of ground-based telescopes, will be combined with observations from other users of ExoClock to ensure that the transit times of these planets continue to be well-known, far into the future.

preprint2022arXiv

Sub-percent accuracy for the intensity of a near infrared water line at 10670 cm$^{-1}$: experiment and analysis

Laser measurements of the intensity of (201) 3$_{22}$ -- (000) 2$_{21}$ near-infrared water absorption line at 10670.1 \cm\ are made using three different Herriott cells. These measurements determine the line intensity with an standard deviation below of 0.3~\% by consideration of the new geometrically derived formula for the opt ical path length without approximations. This determination together with the current accepted value leads to an overall uncertainty of 0.7~\% of the experimentally assessed line intensity which is compared with previous {\it ab initio} predictions. It is found that steady improvements in the both the dipole moment surface (DMS) and potential energy surface (PES) used in the theoretical studies leads to systematic better agreement with the observation, with the most recent prediction agreeing closely with the experiment.

preprint2022arXiv

Twinkle -- a small satellite spectroscopy mission for the next phase of exoplanet science

With a focus on off-the-shelf components, Twinkle is the first in a series of cost competitive small satellites managed and financed by Blue Skies Space Ltd. The satellite is based on a high-heritage Airbus platform that will carry a 0.45 m telescope and a spectrometer which will provide simultaneous wavelength coverage from 0.5-4.5 $\rm{μm}$. The spacecraft prime is Airbus Stevenage while the telescope is being developed by Airbus Toulouse and the spectrometer by ABB Canada. Scheduled to begin scientific operations in 2025, Twinkle will sit in a thermally-stable, sun-synchronous, low-Earth orbit. The mission has a designed operation lifetime of at least seven years and, during the first three years of operation, will conduct two large-scale survey programmes: one focused on Solar System objects and the other dedicated to extrasolar targets. Here we present an overview of the architecture of the mission, refinements in the design approach, and some of the key science themes of the extrasolar survey.

preprint2020arXiv

A Science Gateway for Atomic and Molecular Physics

We describe the creation of a new Atomic and Molecular Physics science gateway (AMPGateway). The gateway is designed to bring together a subset of the AMP community to work collectively to make their codes available and easier to use by the partners as well as others. By necessity, a project such as this requires the developers to work on issues of portability, documentation, ease of input, as well as making sure the codes can run on a variety of architectures. Here we outline our efforts to build this AMP gateway and future directions.

preprint2020arXiv

An update to the MARVEL dataset and ExoMol line list for 12C2

The spectrum of dicarbon (C2) is important in astrophysics and for spectroscopic studies of plasmas and flames. The C2 spectrum is characterized by many band systems with new ones still being actively identified; astronomical observations involve eight of these bands. Recently, Furtenbacher et al. (2016, Astrophys. J. Suppl., 224, 44) presented a set of 5699 empirical energy levels for 12C2, distributed among 11 electronic states and 98 vibronic bands, derived from 42 experimental studies and obtained using the MARVEL (Measured Active Rotational-Vibrational Energy Levels) procedure. Here, we add data from 13 new sources and update data from 5 sources. Many of these data sources characterize high-lying electronic states, including the newly detected 3Pig state. Older studies have been included following improvements in the MARVEL procedure which allow their uncertainties to be estimated. These older works in particular determine levels in the C1Pig state, the upper state of the insufficiently characterized Deslandres-d'Azambuja (C1Pig-A1Piu) band. The new compilation considers a total of 31323 transitions and derives 7047 empirical(MARVEL) energy levels spanning 20 electronic and 142 vibronic states. These new empirical energy levels are used here to update the 8states C2 ExoMol line list. This updated line list is highly suitable for high-resolution cross-correlation studies in astronomical spectroscopy of, for example, exoplanets, as 99.4% of the transitions with intensities over 10^(-18) cm/molecule at 1000 K have frequencies determined by empirical energy levels.

preprint2020arXiv

Calculated line lists for H$_{\text{2}} \mkern0.1mu ^{\text{16}}$O and H$_{\text{2}} \mkern0.1mu ^{\text{18}}$O with extensive comparisons to theoretical and experimental sources including the HITRAN2016 database

New line lists are presented for the two most abundant water isotopologues; H$_{2}$$^{16}$O and H$_{2}$$^{18}$O. The H$_{2}$$^{16}$O line list extends to 25710 cm$^{-1}$ with intensity stabilities provided via ratios of calculated intensities obtained from two different semi-empirical potential energy surfaces. The line list for H$_{2}$$^{18}$O extends to 20000 cm$^{-1}$. The minimum intensity considered for all is $10^{-30}$ cm molecule$^{-1}$ at 296~K, assuming 100\% abundance for each isotopologue. Fluctuation of calculated intensities caused by changes in the underlying potential energy are found to be significant, particularly for weak transitions. Direct comparisons are made against eighteen different sources of line intensities, both experimental and theoretical, many of which are used within the HITRAN2016 database. With some exceptions, there is excellent agreement between our line lists and the experimental intensities in HITRAN2016. In the infrared region, many H$_{2}$$^{16}$O bands which exhibit intensity differences of 5-10\% between to the most recent 'POKAZATEL' line list (Polyansky \textit{et al.}, [Mon. Not. Roy. Astron. Soc. \textbf{480}, 2597 (2018)] and observation, are now generally predicted to within 1\%. For H$_{2}$$^{18}$O, there are systematic differences in the strongest intensities calculated in this work versus those obtained from semi-empirical calculations. In the visible, computed cross sections show smaller residuals between our work and both HITRAN2016 and HITEMP2010 than POKAZATEL. While our line list accurately reproduces HITEMP2010 cross sections in the observed region, residuals produced from this comparison do however highlight the need to update line positions in the visible spectrum of HITEMP2010. These line lists will be used to update many transition intensities and line positions in the HITRAN2016 database.

preprint2020arXiv

Empirical Line Lists in the ExoMol Database

The ExoMol database aims to provide comprehensive molecular line lists for exoplanetary and other hot atmospheres. The data are expanded by inclusion of empirically derived line lists taken from the literature for a series of diatomic molecules, namely CH, NH, OH, AlCl, AlF, OH$^+$, CaF, MgF, KF, NaF, LiCl, LiF, MgH, TiH, CrH, FeH, C$_2$, CP, CN, CaH, and triplet N$_2$. Generally, these line lists are constructed from measured spectra using a combination of effective rotational Hamiltonian models for the line positions and ab initio (transition) dipole moments to provide intensities. This work results in the inclusion of 22 new molecules (36 new isotopologues) in the ExoMol database.

preprint2020arXiv

Empirical normal intensity distribution for overtone vibrational spectra of triatomic molecules

Theoretical calculations are contributing a significantly higher proportion of data to contemporary spectroscopic databases, which have traditionally relied on experimental observations and semi-empirical models. It is now a common procedure to extend calculated line lists to include ro-vibrational transitions between all bound states of the ground electronic state up to the dissociation limit. Advanced ab initio methods are utilized to calculate the potential energy and dipole moment surfaces (PESs and DMSs), and semi-empirical PESs are then obtained by combining ab initio and experimental data. The objective is to reach high accuracy in the calculated transition intensities for all parts of spectrum, i.e. to increase the predictive power of the model. We show that in order to perform this task, one needs, in addition to the standard improvements of the PES and DMS in the spectroscopically accessible regions, to extend the ab initio calculations of the PES towards the united-atom limit along the stretching coordinates. The argument is based on the correlation between the intensities of high-overtone transitions and the repulsive potential wall that has previously been theoretically established for diatomic molecules and is empirically extended here to linear and nonlinear triatomic molecules. We generate partial line lists for water and ozone, and together with an already available line list for carbon dioxide, we derive the normal intensity distribution, which is a direct consequence of this correlation. The normal distribution is not an instrument to compute highly accurate intensities, rather it is a means to analyze the intensities computed by the traditional methods.

preprint2020arXiv

ExoMol line lists -- XL. Ro-vibrational molecular line list for the hydronium ion (H$_3$O$^+$)

A new line list for hydronium (H$_3$$^{16}$O$^+$) is computed. The line list is based on a new ab initio dipole moment surface (CCSD(T)/aug-cc-pVQZ) and a new empirical potential energy surface (PES). The empirical PES of H$_3$O$^+$ was obtained by refining an ab initio surface through a global fit to the experimentally determined ro-vibrational energies collected from the literature covering the ground, $ν_1^{\pm}$, $ν_2^{\pm}$, $2ν_2^{\pm}$, $ν_3^{\pm}$ and $ν_4^{\pm}$ vibrational states. The line list covers the wavenumber range up to 10 000 cm$^{-1}$ (wavelengths $>1 μ$m) and should be complete for temperatures up to $T=1500$ K. This is the first comprehensive line list for H$_3$O$^+$ with extensive wavenumber coverage and accurate transitional probabilities. Prospects of detection of hydronium in spectra of solar system giant planets as well as exoplanets are discussed. The eXeL line list is publicly available from the ExoMol and CDS databases.

preprint2020arXiv

ExoMol molecular line lists -- XXXVII: spectra of acetylene

A new ro-vibrational line list for the ground electronic state of the main isotopologue of acetylene, $^{12}$C$_2$H$_2$, is computed as part of the ExoMol project. The aCeTY line list covers the transition wavenumbers up to 10,000 cm$^{-1}$ ($ λ>1$ $μ$m), with lower and upper energy levels up to 12,000 cm$^{-1}$ and 22,000 cm$^{-1}$ considered, respectively. The calculations are performed up to a maximum value for the vibrational angular momentum, $K_{\rm max}=L_{\rm max}$ = 16, and maximum rotational angular momentum, $J$ = 99. Higher values of $J$ were not within the specified wavenumber window. The aCeTY line list is considered to be complete up to 2200 K, making it suitable for use in characterising high-temperature exoplanet or cool stellar atmospheres. Einstein-A coefficients, which can directly be used to calculate intensities at a particular temperature, are computed for 4.3 billion (4,347,381,911) transitions between 5 million (5,160,803) energy levels. We make comparisons against other available data for $^{12}$C$_2$H$_2$, and demonstrate this to be the most complete line list available. The line list is available in electronic form from the online CDS and ExoMol databases.

preprint2020arXiv

Ions in the Thermosphere of Exoplanets: Observable Constraints Revealed by Innovative Laboratory Experiments

With the upcoming launch of space telescopes dedicated to the study of exoplanets, the \textit{Atmospheric Remote-Sensing Infrared Exoplanet Large-survey} (ARIEL) and the \textit{James Webb Space Telescope} (JWST), a new era is opening in the exoplanetary atmospheric explorations. However, especially in relatively cold planets around later-type stars, photochemical hazes and clouds may mask the composition of the lower part of the atmosphere, making it difficult to detect any chemical species in the troposphere or to understand whether there is a surface or not. This issue is particularly exacerbated if the goal is to study the habitability of said exoplanets and to search for biosignatures.\par This work combines innovative laboratory experiments, chemical modeling and simulated observations at ARIEL and JWST resolutions. We focus on the signatures of molecular ions that can be found in upper atmospheres above cloud decks. Our results suggest that H$_3^+$ along with H$_3$O$^+$ could be detected in the observational spectra of sub-Neptunes based on realistic mixing ratio assumption. This new parametric set may help to distinguish super-Earths with a thin atmosphere from H$_2$-dominated sub-Neptunes, to address the critical question whether a low-gravity planet around a low-mass active star is able to retain its volatile components. These ions may also constitute potential tracers to certain molecules of interest such as H$_2$O or O$_2$ to probe the habitability of exoplanets. Their detection will be an enthralling challenge for the future JWST and ARIEL telescopes.

preprint2020arXiv

Molecular Cross Sections for High Resolution Spectroscopy of Super Earths, Warm Neptunes and Hot Jupiters

High resolution spectroscopy (HRS) has been used to detect a number of species in the atmospheres of hot Jupiters. Key to such detections is accurately and precisely modelled spectra for cross-correlation against the R$\gtrsim$20,000 observations. There is a need for the latest generation of opacities which form the basis for high signal-to-noise detections using such spectra. In this study we present and make publicly available cross sections for six molecular species, H$_2$O, CO, HCN, CH$_4$, NH$_3$ and CO$_2$ using the latest line lists most suitable for low- and high-resolution spectroscopy. We focus on the infrared (0.95-5~$μ$m) and between 500-1500~K where these species have strong spectral signatures. We generate these cross sections on a grid of pressures and temperatures typical for the photospheres of super Earth, warm Neptunes and hot Jupiters using the latest H$_2$ and He pressure broadening. We highlight the most prominent infrared spectral features by modelling three representative exoplanets, GJ~1214~b, GJ~3470~b and HD~189733~b, which encompass a wide range in temperature, mass and radii. In addition, we verify the line lists for H$_2$O, CO and HCN with previous high resolution observations of hot Jupiters. However, we are unable to detect CH$_4$ with our new cross sections from HRS observations of HD~102195~b. These high accuracy opacities are critical for atmospheric detections with HRS and will be continually updated as new data becomes available.

preprint2020arXiv

Original Research By Young Twinkle Students (ORBYTS): Ephemeris Refinement of Transiting Exoplanets

We report follow-up observations of transiting exoplanets that have either large uncertainties (>10 minutes) in their transit times or have not been observed for over three years. A fully robotic ground-based telescope network, observations from citizen astronomers and data from TESS have been used to study eight planets, refining their ephemeris and orbital data. Such follow-up observations are key for ensuring accurate transit times for upcoming ground and space-based telescopes which may seek to characterise the atmospheres of these planets. We find deviations from the expected transit time for all planets, with transits occurring outside the 1 sigma uncertainties for seven planets. Using the newly acquired observations, we subsequently refine their periods and reduce the current predicted ephemeris uncertainties to 0.28 - 4.01 minutes. A significant portion of this work has been completed by students at two high schools in London as part of the Original Research By Young Twinkle Students (ORBYTS) programme.

preprint2020arXiv

Original Research By Young Twinkle Students (ORBYTS): Ephemeris Refinement of Transiting Exoplanets II

We report follow-up observations of four transiting exoplanets, TRES-2b, HAT-P-22b, HAT-P-36b and XO-2b, as part of the Original Research By Young Twinkle Students (ORBYTS) programme. These observations were taken using the Las Cumbres Observatory Global Telescope Network's (LCOGT) robotic 0.4 m telescopes and were analysed using the HOlomon Photometric Software (HOPS). Such observations are key for ensuring accurate transit times for upcoming telescopes, such as the James Webb Space Telescope (JWST), Twinkle and Ariel, which may seek to characterise the atmospheres of these planets. The data have been uploaded to ExoClock and a significant portion of this work has been completed by secondary school students in London.

preprint2020arXiv

Spectroscopy of YO from first principles

We report an ab initio study on the spectroscopy of the open-shell diatomic molecule yttrium oxide, YO. The study considers the six lowest doublet states, $X\,{}^{2}Σ^{+}$, $A\,{}^{2}Π$, $C\,{}^{2}Π$, $A'\,{}^{2}Δ$, $B\,{}^{2}Σ^{+}$, $D\,{}^{2}Σ^{+}$ and a few higher-lying quartet states using high levels of electronic structure theory and accurate nuclear motion calculations. The coupled cluster singles, doubles, and perturbative triples, CCSD(T), and multireference configuration interaction (MRCI) methods are employed in conjunction with a relativistic pseudopotential on the yttrium atom and a series of correlation-consistent basis sets ranging in size from triple-$ζ$ to quintuple-$ζ$ quality. Core-valence correlation effects are taken into account and complete basis set limit extrapolation is performed for CCSD(T). Spin-orbit coupling is included through the use of both MRCI state-interaction with spin-orbit (SI-SO) approach and four-component relativistic equation-of-motion CCSD calculations. Using the ab initio data for bond lengths ranging from 1.0 to 2.5 A, we compute 6 potential energy, 12 spin-orbit, 8 electronic angular momentum, 6 electric dipole moment and 12 transition dipole moment (4 parallel and 8 perpendicular) curves which provide a complete description of the spectroscopy of the system of six lowest doublet states. The Duo nuclear motion program is used to solve the coupled nuclear motion Schrödinger equation for these six electronic states. The spectra of $^{89}$Y$^{16}$O simulated for different temperatures are compared with several available high resolution experimental studies; good agreement is found once minor adjustments are made to the electronic excitation energies.

preprint2020arXiv

The 2020 release of the ExoMol database: molecular line lists for exoplanet and other hot atmospheres

The ExoMol database (www.exomol.com) provides molecular data for spectroscopic studies of hot atmospheres. While the data is intended for studies of exoplanets and other astronomical bodies, the dataset is widely applicable. The basic form of the database is extensive line lists; these are supplemented with partition functions, state lifetimes, cooling functions, Landé g-factors, temperature-dependent cross sections, opacities, pressure broadening parameters, $k$-coefficients and dipoles. This paper presents the latest release of the database which has been expanded to consider 80 molecules and 190 isotopologues totaling over 700 billion transitions. While the spectroscopic data is concentrated at infrared and visible wavelengths, ultraviolet transitions are being increasingly considered in response to requests from observers. The core of the database comes from the ExoMol project which primarily uses theoretical methods, albeit usually fine-tuned to reproduce laboratory spectra, to generate very extensive line lists for studies of hot bodies. The data has recently been supplemented by line lists deriving from direct laboratory observations, albeit usually with the use of ab initio transition intensities. A major push in the new release is towards accurate characterisation of transition frequencies for use in high resolution studies of exoplanets and other bodies.

preprint2020arXiv

The ExoMolOP Database: Cross-sections and k-tables for Molecules of Interest in High-Temperature Exoplanet Atmospheres

A publicly available database of opacities for molecules of astrophysical interest, ExoMolOP, has been compiled for over 80 species, based on the latest line list data from the ExoMol, HITEMP and MoLLIST databases. These data are generally suitable for characterising high temperature exoplanet or cool stellar/substellar atmospheres, and have been computed at a variety of pressures and temperatures, with a few molecules included at room-temperature only from the HITRAN database. The data are formatted in different ways for four different exoplanet atmosphere retrieval codes; ARCiS, TauREx, NEMESIS and petitRADTRANS, and include both cross-sections (at R~=~$\fracλ{Δλ}$~=~15,000) and k-tables (at R~=~$\fracλ{Δλ}$~=~1000) for the 0.3~-~50$μ$m wavelength region. Opacity files can be downloaded and used directly for these codes. Atomic data for alkali metals Na and K are also included, using data from the NIST database and the latest line shapes for the resonance lines. Broadening parameters have been taken from the literature where available, or from those for a known molecule with similar molecular properties where no broadening data are available. The data are available from www.exomol.com.

preprint2020arXiv

The high-temperature rotation-vibration spectrum and rotational clustering of silylene (SiH$_2$)

A rotation-vibration line list for the electronic ground state ($\tilde{X}^{1}A_{1}$) of SiH$_2$ is presented. The line list, named CATS, is suitable for temperatures up to 2000 K and covers the wavenumber range 0 - 10,000 cm$^{-1}$(wavelengths $>1.0$ $μ$m) for states with rotational excitation up to $J=52$. Over 310 million transitions between 593 804 energy levels have been computed variationally with a new empirically refined potential energy surface, determined by refining to 75 empirical term values with $J\leq 5$ and a newly computed high-level ab initio dipole moment surface. This is the first, comprehensive high-temperature line list to be reported for SiH$_2$ and it is expected to aid the study of silylene in plasma physics, industrial processes and possible astronomical detection. Furthermore, we investigate the phenomenon of rotational energy level clustering in the spectrum of SiH$_2$. The CATS line list is available from the ExoMol database (www.exomol.com) and the CDS database.

preprint2020arXiv

Use of the complete basis set limit for computing highly accurate ab initio dipole moments

Calculating dipole moments with high-order basis sets is generally only possible for the light molecules, such as water. A simple, yet highly effective strategy of obtaining high-order dipoles with small, computationally less expensive basis sets is described. Using the finite field method for computing dipoles, energies calculated with small basis sets can be extrapolated to produce dipoles that are comparable to those obtained in high order calculations. The method reduces computational resources by approximately 50% (allowing the calculation of reliable dipole moments for larger molecules) and simultaneously improves the agreement with experimentally measured infrared transition intensities. For atmospherically important molecules which are typically too large to consider the use of large basis sets, this procedure will provide the necessary means of improving calculated spectral intensities by several percent.

preprint2019arXiv

Analysis of first overtone bands of isotopologues of CO and SiO in stellar spectra

The first overtone ($Δv$ = 2) bands of the monosubstituted isotopologues of CO at 2.3 $μ$m in the spectrum of Arcturus (K2 III), and of the monosubstituted isotopologues of SiO at 4 $μ$m in the spectrum of the red giant HD196610 (M6 III) are modelled. To investigate problems involving the computation of the first overtone bands of isotopologues of CO and SiO in spectra of late-type stars and to determine isotopic abundances. We use fits of theoretical synthetic spectra to the observed stellar molecular bands of CO and SiO to determine abundances of isotopes of C, O and Si. Fits of synthetic spectra of the \CO first overtone bands at 2.3 $μ$m computed with three available line lists (Goorvitch, HITEMP2010 and Li et al.) to the observed spectrum of Arcturus provide the same carbon abundance [C]=$-0.6$ and isotopic ratio of carbon $^{12}$C/$^{13}$C=10 $\pm$ 2. However, the quality of fits to the observed spectrum differ for three line lists used. Furthermore, the derived oxygen isotopic ratio $^{16}$O/$^{18}$O = 2000 $\pm$ 500 is larger than that known in the solar system where $^{16}$O/$^{18}$O = 500. The silicon isotopic ratio in the atmosphere of the red giant HD196610 is revised. Using the ExoMol SiO line list with appropriate statistical weights for the SiO isotopologues the `non-solar' ratio $^{28}$Si:$^{29}$Si:$^{30}$Si = 0.86$\pm$0.03:0.12$\pm$0.02:0.02$\pm$0.01 is obtained. We found that a) the computed isotopic carbon and silicon ratios determined by the fits to the observed spectrum depend on the adopted abundance of C and Si, respectively; b) Correct treatment of the nuclear spin degeneracies parameter is of crucial importance for the use of nowadays HITRAN/ExoMol line lists in the astrophysical computations.

preprint2019arXiv

UKRmol+: a suite for modelling of electronic processes in molecules interacting with electrons, positrons and photons using the R-matrix method

UKRmol+ is a new implementation of the UK R-matrix electron-molecule scattering code. Key features of the implementation are the use of quantum chemistry codes such as Molpro to provide target molecular orbitals; the optional use of mixed Gaussian -- B-spline basis functions to represent the continuum and improved configuration and Hamiltonian generation. The code is described, and examples covering electron collisions from a range of targets, positron collisions and photionisation are presented. The codes are freely available as a tarball from Zenodo.