Source author record

Laura K. McKemmish

Laura K. McKemmish appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

12 published item(s)

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

Assessing Spectra and Thermal Inversions due to TiO in Hot Jupiter Atmospheres

Recent detections of thermal inversions in the dayside atmospheres of some hot Jupiters are motivating new avenues to understand the interplay between their temperature structures and other atmospheric conditions. In particular, TiO has long been proposed to cause thermal inversions in hot Jupiters, depending on other factors such as stellar irradiation, C/O, and vertical mixing. TiO also has spectral features in the optical and near-infrared that have been detected. However, interpretations of TiO signatures rely on the accuracy of TiO opacity used in the models. The recently reported Toto TiO line list provides a new opportunity to investigate these dependencies, which is the goal of the present work. First, we investigate how the Toto line list affects observable transmission and emission spectra of hot Jupiters at low and high resolution. The improvement in the Toto line list compared to a previous line list results in observable differences in the model spectra, particularly in the optical at high resolution. Secondly, we explore the interplay between temperature structure, irradiation and composition with TiO as the primary source of optical opacity, using 1D self-consistent atmospheric models. Among other trends, we find that the propensity for thermal inversions due to TiO peaks at C/O$\sim$0.9, consistent with recent studies. Using these models, we further assess metrics to quantify thermal inversions due to TiO, compared to frequently-used Spitzer photometry, over a range in C/O, irradiation, metallicity, gravity and stellar type.

preprint2020arXiv

Computational insight into diatomic molecules as probes to measure the variation of the proton-to-electron mass ratio

Astrophysical molecular spectroscopy is an important means of searching for new physics through probing the variation of the proton-to-electron mass ratio, $μ$. New molecular probes could provide tighter constraints on the variation of $μ$ and better direction for theories of new physics. Here we summarise our previous paper \citep{19SyMoCu.CN} for astronomers, highlighting the importance of accurate estimates of peak molecular abundance and temperature as well as spectral resolution and sensitivity of telescopes in different regions of the electromagnetic spectrum. Whilst none of the 11 astrophysical diatomic molecules we investigated showed enhanced sensitive rovibronic transitions at observable intensities for astrophysical environments, we have gained a better understanding of the factors that contribute to high sensitivities. From our results, CN, CP, SiN and SiC have shown the most promise of all astrophysical diatomic molecules for further investigation, with further work currently being done on CN.

preprint2020arXiv

Mixed ramp-Gaussian basis sets for core-dependent properties: STO-RG and STO-R2G for Li-Ne

The traditional Gaussian basis sets used in modern quantum chemistry lack an electron-nuclear cusp, and hence struggle to accurately describe core electron properties. A recently introduced novel type of basis set, mixed ramp-Gaussians, introduce a new primitive function called a ramp function which addresses this issue. This paper introduces three new mixed ramp-Gaussian basis sets - STO-R, STO-RG and STO-R2G, made from a linear combination of ramp and Gaussian primitive functions - which are derived from the single-core-zeta Slater basis sets for the elements Li to Ne. This derivation is done in an analogous fashion to the famous STO-$n$G basis sets. The STO-RG basis functions are found to outperform the STO-3G basis functions and STO-R2G outperforms STO-6G, both in terms of wavefunction fit and other key quantities such as the one-electron energy and the electron-nuclear cusp. The second part of this paper performs preparatory investigations into how standard all-Gaussian basis sets can be converted to ramp-Gaussian basis sets through modifying the core basis functions. Using a test case of the 6-31G basis set for carbon, we determined that the second Gaussian primitive is less important when fitting a ramp-Gaussian core basis function directly to an all-Gaussian core basis function than when fitting to a Slater basis function. Further, we identified the basis sets that are single-core-zeta and thus should be most straightforward to convert to mixed ramp-Gaussian basis sets in the future.

preprint2020arXiv

Non-detection of TiO and VO in the atmosphere of WASP-121b using high-resolution spectroscopy

Thermal inversions have long been predicted to exist in the atmospheres of ultra-hot Jupiters. However, detection of two species thought to be responsible -- TiO and VO -- remain elusive. We present a search for TiO and VO in the atmosphere of the ultra-hot Jupiter WASP-121b ($T_\textrm{eq} \gtrsim 2400$ K), an exoplanet already known to show water features in its dayside spectrum characteristic of a temperature inversion as well as tentative evidence for VO at low-resolution. We observed its transmission spectrum with UVES/VLT and used the cross-correlation method -- a powerful tool for the unambiguous identification of the presence of atomic and molecular species -- in an effort to detect whether TiO or VO were responsible for the observed temperature inversion. No evidence for the presence of TiO or VO was found at the terminator of WASP-121b. By injecting signals into our data at varying abundance levels, we set rough detection limits of $[\text{VO}] \lesssim -7.9$ and $[\text{TiO}] \lesssim -9.3$. However, these detection limits are largely degenerate with scattering properties and the position of the cloud deck. Our results may suggest that neither TiO or VO are the main drivers of the thermal inversion in WASP-121b, but until a more accurate line list is developed for VO, we cannot conclusively rule out its presence. Future work will search for finding other strong optically-absorbing species that may be responsible for the excess absorption in the red-optical.

preprint2020arXiv

On the Computation of Dipole Moments: A Recommendation on the Choice of the Basis Set and the Level of Theory

Together with experimental data, theoretically predicted dipole moments represent a valuable tool for different branches in the chemical and physical sciences. With the diversity of levels of theory and basis sets available, a reliable combination must be carefully chosen in order to achieve accurate predictions. In a recent publication (arXiv:1709.05075 [physics.chem-ph]), Hait and Head-Gordon took a first step in this regard by providing recommendations on the best density functionals suitable for these purposes. However, no extensive study has been performed to provide recommendations on the basis set choice. Here, we shed some light into this matter by evaluating the performance of 38 general-purpose basis sets of single up to triple zeta-quality, when coupled with nine different levels of theory, in the computation of dipole moments. The calculations were performed on a data set with 114 small molecules containing second- and third-row elements. We based our analysis in regularised root mean square errors (regularised RMSE) as this procedure ensures relative errors for ionic species and absolute errors for species with small dipole moment values. Our results indicate that the best compromise between accuracy and computational efficiency is achieved by performing the computations with an augmented double zeta-quality basis set (i.e. aug-pc-1, aug-pcseg-1, aug-cc-pVDZ) together with a hybrid functional (e.g. wB97X-V, SOGGA11-X). Augmented triple-zeta basis sets could enhance the accuracy of the computations, but the computational cost of introducing such a basis set is substantial compared with the small improvement provided. These findings also highlight the crucial role that augmentation of the basis set with diffuse functions on both hydrogen and non-hydrogen atoms plays in the computation of dipole moments.

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.

preprint2016arXiv

Ab initio calculations to support accurate modelling of the rovibronic spectroscopy calculations of vanadium monoxide (VO)

Accurate knowledge of the rovibronic near-infrared and visible spectra of vanadium monoxide (VO) is very important for studies of cool stellar and hot planetary atmospheres. Here, the required ab initio dipole moment and spin-orbit coupling curves for VO are produced. This data forms the basis of a new VO line list considering 13 different electronic states and containing over 277 million transitions. Open shell transition, metal diatomics are challenging species to model through ab initio quantum mechanics due to the large number of low-lying electronic states, significant spin-orbit coupling and strong static and dynamic electron correlation. Multi-reference configuration interaction methodologies using orbitals from a complete active space self-consistent-field (CASSCF) calculation are the standard technique for these systems. We use different state-specific or minimal-state CASSCF orbitals for each electronic state to maximise the calculation accuracy. The off-diagonal dipole moment controls the intensity of electronic transitions. We test finite-field off-diagonal dipole moments, but found that (1) the accuracy of the excitation energies were not sufficient to allow accurate dipole moments to be evaluated and (2) computer time requirements for perpendicular transitions were prohibitive. The best off-diagonal dipole moments are calculated using wavefunctions with different CASSCF orbitals.

preprint2016arXiv

ExoMol line lists XVIII. The high temperature spectrum of VO

An accurate line list, VOMYT, of spectroscopic transitions is presented for hot VO. The 13 lowest electronic states are considered. Curves and couplings are based on initial {\it ab initio} electronic structure calculations and then tuned using available experimental data. Dipole moment curves, used to obtain transition intensities, are computed using high levels of theory (e.g. MRCI/aug-cc-pVQZ using state-specific or minimal-state CAS for dipole moments). This line list contains over 277 million transitions between almost 640,000 energy levels. It covers the wavelengths longer than 0.29 $μ$m and includes all transitions from energy levels within the lowest nine electronic states which have energies less than 20,000 \cm{} to upper states within the lowest 13 electronic states which have energies below 50,000 \cm{}. The line lists give significantly increased absorption at infrared wavelengths compared to currently available VO line lists. The full line lists is made available in electronic form via the CDS database and at www.exomol.com.

preprint2016arXiv

Low-temperature chemistry using the R-matrix method

Techniques for producing cold and ultracold molecules are enabling the study of chemical reactions and scattering at the quantum scattering limit, with only a few partial waves contributing to the incident channel, leading to the observation and even full control of state-to-state collisions in this regime. A new R-matrix formalism is presented for tackling problems involving low- and ultra-low energy collisions. This general formalism is particularly appropriate for slow collisions occurring on potential energy surfaces with deep wells. The many resonance states make such systems hard to treat theoretically but offer the best prospects for novel physics: resonances are already being widely used to control diatomic systems and should provide the route to steering ultracold reactions. Our R-matrix-based formalism builds on the progress made in variational calculations of molecular spectra by using these methods to provide wavefunctions for the whole system at short internuclear distances, (a regime known as the inner region). These wavefunctions are used to construct collision energy-dependent R-matrices which can then be propagated to give cross sections at each collision energy. The method is formulated for ultracold collision systems with differing numbers of atoms.

preprint2016arXiv

The {\it ab initio} calculation of spectra of open shell diatomic molecules

The spectra (rotational, rotation-vibrational or electronic) of diatomic molecules due to transitions involving only closed-shell ($^1Σ$) electronic states follow very regular, simple patterns and their theoretical analysis is usually straightforward. On the other hand, open-shell electronic states lead to more complicated spectral patterns and, moreover, often appear as a manifold of closely lying electronic states, leading to perturbations with even larger complexity. This is especially true when at least one of the atoms is a transition metal. Traditionally these complex cases have been analysed using approaches based on perturbation theory, with semi-empirical parameters determined by fitting to spectral data. Recently the needs of two rather diverse scientific areas have driven the demand for improved theoretical models of open-shell diatomic systems based on an \emph{ab initio} approach, these areas are ultracold chemistry and the astrophysics of "cool" stars, brown dwarfs and most recently extrasolar planets. However, the complex electronic structure of these molecules combined with the accuracy requirements of high-resolution spectroscopy render such an approach particularly challenging. This review describes recent progress in developing methods for directly solving the effective Schrödinger equation for open-shell diatomic molecules, with a focus on molecules containing a transtion metal. It considers four aspects of the problem: 1. The electronic structure problem, 2. Non-perturbative treatments of the curve couplings, 3. The solution of the nuclear motion Schrödinger equation, 4. The generation of accurate electric dipole transition intensities. Examples of applications are used to illustrate these issues.

preprint2011arXiv

Quantum entanglement between electronic and vibrational degrees of freedom in molecules

We consider the quantum entanglement of the electronic and vibrational degrees of freedom in molecules with a tendency towards double welled potentials using model coupled harmonic diabatic potential-energy surfaces. The von Neumann entropy of the reduced density matrix is used to quantify the electron-vibration entanglement for the lowest two vibronic wavefunctions in such a bipartite system. Significant entanglement is found only in the region in which the ground vibronic state contains a density profile that is bimodal (i.e., contains two separate local minima). However, in this region two distinct types of entanglement are found: (1) entanglement that arises purely from the degeneracy of energy levels in the two potential wells and which is destroyed by slight asymmetry, and (2) entanglement that involves strongly interacting states in each well that is relatively insensitive to asymmetry. These two distinct regions are termed fragile degeneracy-induced entanglement and persistent entanglement, respectively. Six classic molecular systems describable by two diabatic states are considered: ammonia, benzene, semibullvalene, pyridine excited triplet states, the Creutz-Taube ion, and the radical cation of the "special pair" of chlorophylls involved in photosynthesis. These chemically diverse systems are all treated using the same general formalism and the nature of the entanglement that they embody is elucidated.