Showing papers in "arXiv: Chemical Physics in 2020"

CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations.

Abstract: CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2k to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.

NWChem: Past, Present, and Future

Abstract: Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principledriven methodologies to model complex chemical and materials processes. Over the last few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach and outlook.

Machine Learning Force Fields

Abstract: In recent years, the use of Machine Learning (ML) in computational chemistry has enabled numerous advances previously out of reach due to the computational complexity of traditional electronic-structure methods. One of the most promising applications is the construction of ML-based force fields (FFs), with the aim to narrow the gap between the accuracy of ab initio methods and the efficiency of classical FFs. The key idea is to learn the statistical relation between chemical structure and potential energy without relying on a preconceived notion of fixed chemical bonds or knowledge about the relevant interactions. Such universal ML approximations are in principle only limited by the quality and quantity of the reference data used to train them. This review gives an overview of applications of ML-FFs and the chemical insights that can be obtained from them. The core concepts underlying ML-FFs are described in detail and a step-by-step guide for constructing and testing them from scratch is given. The text concludes with a discussion of the challenges that remain to be overcome by the next generation of ML-FFs.

Photoelectrochemical water splitting in separate oxygen and hydrogen cells

Abstract: Solar water splitting provides a promising path for sustainable hydrogen production and solar energy storage. One of the greatest challenges towards large-scale utilization of this technology is reducing the hydrogen production cost. The conventional electrolyzer architecture, where hydrogen and oxygen are co-produced in the same cell, gives rise to critical challenges in photoelectrochemical (PEC) water splitting cells that directly convert solar energy and water to hydrogen. Here we overcome these challenges by separating the hydrogen and oxygen cells. The ion exchange in our cells is mediated by auxiliary electrodes, and the cells are connected to each other only by metal wires, enabling centralized hydrogen production. We demonstrate hydrogen generation in separate cells with solar-to-hydrogen conversion efficiency of 7.5%, which can readily surpass 10% using standard commercial components. A basic cost comparison shows that our approach is competitive with conventional PEC systems, enabling safe and potentially affordable solar hydrogen production.

The DIRAC code for relativistic molecular calculations.

Abstract: DIRAC is a freely distributed general-purpose program system for 1-, 2- and 4-component relativistic molecular calculations at the level of Hartree--Fock, Kohn--Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, coupled cluster and electron propagator theory. At the self-consistent-field level a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding, and frozen density embedding models. DIRAC was one of the earliest codes for relativistic molecular calculations and remains a reference in its field.

Perspective: Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM)

Abstract: An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system--bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article, presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.

Attosecond-resolved photoionization of chiral molecules

Abstract: Chiral light-matter interactions have been investigated for two centuries, leading to the discovery of many chiroptical processes used for discrimination of enantiomers. Whereas most chiroptical effects result from a response of bound electrons, photoionization can produce much stronger chiral signals that manifest as asymmetries in the angular distribution of the photoelectrons along the light propagation axis. Here we implement a self-referenced attosecond photoelectron interferometry to measure the temporal profile of the forward and backward electron wavepackets emitted upon photoionization of camphor by circularly polarized laser pulses. We found a delay between electrons ejected forward and backward, which depends on the ejection angle and reaches 24 attoseconds. The asymmetric temporal shape of electron wavepackets emitted through an autoionizing state further reveals the chiral character of strongly-correlated electronic dynamics.

Atom-Tunneling in Chemistry

Abstract: Quantum mechanical tunneling of atoms is increasingly found to play an important role in many chemical transformations. Experimentally, atom-tunneling can be indirectly detected by temperature-independent rate constants at low temperature or by enhanced kinetic isotope effects. On the contrary, using computational investigations the influence of tunneling on the reaction rates can directly be monitored. The tunnel effect, for example, changes reaction paths and branching ratios, enables chemical reactions in an astrochemical environment that would be impossible by thermal transition, and influences biochemical processes.

Mathematical Model of a Direct Methanol Fuel Cell.

Abstract: A one dimensional (1-D), isothermal model for a direct methanol fuel cell (DMFC) is presented. This model accounts for the kinetics of the multi-step methanol oxidation reaction at the anode. Diffusion and crossover of methanol are modeled and the mixed potential of the oxygen cathode due to methanol crossover is included. Kinetic and diffusional parameters are estimated by comparing the model to data from a 25 cm2 DMFC. This semi-analytical model can be solved rapidly so that it is suitable for inclusion in real-time system level DMFC simulations.

Delta-Machine Learning for Potential Energy Surfaces: A PIP approach to bring a DFT-based PES to CCSD(T) Level of Theory

Abstract: ``$\Delta$-machine learning" refers to a machine learning approach to bring a property such as a potential energy surface (PES) based on low-level (LL) density functional theory (DFT) energies and gradients to close to a coupled cluster (CC) level of accuracy. Here we present such an approach that uses the permutationally invariant polynomial (PIP) method to fit high-dimensional PESs. The approach is represented by a simple equation, in obvious notation $V_{LL{\rightarrow}CC}=V_{LL}+\Delta{V_{CC-LL}}$, and demonstrated for \ce{CH4}, \ce{H3O+}, and $trans$ and $cis$-$N$-methyl acetamide (NMA), \ce{CH3CONHCH3}. For these molecules, the LL PES, $V_{LL}$, is a PIP fit to DFT/B3LYP/6-31+G(d) energies and gradients, and $\Delta{V_{CC-LL}}$ is a precise PIP fit obtained using a low-order PIP basis set and based on a relatively small number of CCSD(T) energies. For \ce{CH4} these are new calculations adopting an aug-cc-pVDZ basis, for \ce{H3O+} previous CCSD(T)-F12/aug-cc-pVQZ energies are used, while for NMA new CCSD(T)-F12/aug-cc-pVDZ calculations are performed. With as few as 200 CCSD(T) energies, the new PESs are in excellent agreement with benchmark CCSD(T) results for the small molecules, and for 12-atom NMA training is done with 4696 CCSD(T) energies.

OrbNet: Deep Learning for Quantum Chemistry Using Symmetry-Adapted Atomic-Orbital Features

Abstract: We introduce a machine learning method in which energy solutions from the Schrodinger equation are predicted using symmetry adapted atomic orbitals features and a graph neural-network architecture. \textsc{OrbNet} is shown to outperform existing methods in terms of learning efficiency and transferability for the prediction of density functional theory results while employing low-cost features that are obtained from semi-empirical electronic structure calculations. For applications to datasets of drug-like molecules, including QM7b-T, QM9, GDB-13-T, DrugBank, and the conformer benchmark dataset of Folmsbee and Hutchison, \textsc{OrbNet} predicts energies within chemical accuracy of DFT at a computational cost that is thousand-fold or more reduced.

Influence of Ti$_3^+$ Defect-type on Heterogeneous Photocatalytic H$_2$ Evolution Activity of TiO$_2$

Abstract: Reduced titanium dioxide has recently attracted large attention, particularly for its unique co-catalyst-free H$_2$ heterogeneous photocatalytic application. The enhanced photocatalytic activity of the reduced TiO$_2$ was previously ascribed to the introduction of point crystal defects (mainly Ti$_3^+$ centers), which result in the formation of intrinsic co-catalytic centers and enhanced visible light absorption. In this work, we systematically investigate the effect of different defects in the TiO$_{2-x}$ lattice on photocatalytic H$_2$ evolution. To introduce different types of defects, thermal annealing in air (oxidative), Ar (inert), Ar/H$_2$ (reducing), and H$_2$ (reducing) atmospheres were performed on commercially available anatase nanopowder. Then, the powders were characterized by scanning electron microscopy (SEM), X-ray diffraction (XRD) and high-resolution transmission electron microscopy (HR-TEM) to clarify the effect of treatment on material properties. Furthermore, the defect types were characterized by electron paramagnetic resonance (EPR) spectroscopy. We show that thermal annealing in different atmospheres can form different amounts of different defect types in the TiO$_2$ structure. The highest photocatalytic activation is achieved by annealing the anatase powder in a reducing atmosphere for an appropriate temperature/annealing time. By combining the results from H$_2$ generation and EPR analysis we show that the simultaneous presence of two types of defects, i.e. surface exposed Ti$_3^+$ and lattice embedded Ti$_3^+$ centers, in an optimum low concentration, is the determining factor for an optimized photocatalytic H$_2$ evolution rate. In fact, annealing anatase powder under the so-reported optimized conditions in reducing atmosphere leads to the generation of a considerable amount of H$_2$, with rates as high as 338 $\mu$molh-1g-1.

QUESTDB: a database of highly-accurate excitation energies for the electronic structure community

Abstract: We describe our efforts of the past few years to create a large set of more than 500 highly-accurate vertical excitation energies of various natures ($\pi \to \pi^*$, $n \to \pi^*$, double excitation, Rydberg, singlet, doublet, triplet, etc) in small- and medium-sized molecules. These values have been obtained using an incremental strategy which consists in combining high-order coupled cluster and selected configuration interaction calculations using increasingly large diffuse basis sets in order to reach high accuracy. One of the key aspect of the so-called QUEST database of vertical excitations is that it does not rely on any experimental values, avoiding potential biases inherently linked to experiments and facilitating theoretical cross comparisons. Following this composite protocol, we have been able to produce theoretical best estimate (TBEs) with the aug-cc-pVTZ basis set for each of these transitions, as well as basis set corrected TBEs (i.e., near the complete basis set limit) for some of them. The TBEs/aug-cc-pVTZ have been employed to benchmark a large number of (lower-order) wave function methods such as CIS(D), ADC(2), CC2, STEOM-CCSD, CCSD, CCSDR(3), CCSDT-3, ADC(3), CC3, NEVPT2, and others (including spin-scaled variants). In order to gather the huge amount of data produced during the QUEST project, we have created a website [this https URL] where one can easily test and compare the accuracy of a given method with respect to various variables such as the molecule size or its family, the nature of the excited states, the type of basis set, etc. We hope that the present review will provide a useful summary of our effort so far and foster new developments around excited-state methods.

Intermolecular interactions in optical cavities: an ab initio QED study

Abstract: Intermolecular bonds are weak compared to covalent bonds, but they are strong enough to influence the properties of large molecular systems. In this work, we investigate how strong light-matter coupling inside an optical cavity can modify these intermolecular forces. We perform a detailed comparison between currently available ab initio electron-photon methodologies. The electromagnetic field inside the cavity can modulate the ground state properties of weakly bound complexes. Controlling the field polarization, the interactions can be stabilized or destabilized, and electron densities, dipole moments, and polarizabilities can be altered. We demonstrate that electron-photon correlation is fundamental to describe intermolecular interactions in strong light-matter coupling. This work proposes optical cavities as a novel tool to manipulate and control ground state properties, solvent effects, and intermolecular interactions for molecules and materials.

Niobium disulphide (NbS$_2$)-based (heterogeneous) electrocatalysts for an efficient hydrogen evolution reaction

Abstract: The design of efficient and cost-effective catalysts for the hydrogen evolution reaction (HER) is the key for molecular hydrogen (H2) production from electrochemical water splitting. Transition metal dichalcogenides (MX2), most notably group-6 MX2 (e.g., MoS2 and WS2), are appealing catalysts for the HER alternative to the best, but highly expensive, Pt-group elements. However, their HER activity is typically restricted to their edge sites rather than their basal plane. Furthermore, their semiconducting properties hinder an efficient electron transfer to the catalytic sites, which impedes a high rate of H2 production. Herein, we exploit liquid-phase exfoliation-produced metallic (1H, 2H and 3R) NbS2 nanoflakes, belonging to the class of metallic layered group-5 MX2, to overcome the abovementioned limitations. Both chemical treatment with hygroscopic Li salt and electrochemical in operando self-nanostructuring are exploited to improve the NbS2 nanoflake HER activity. The combination of NbS2 with other MX2, in our case MoSe2, also provides heterogeneous catalysts accelerating the HER kinetics of the individual counterparts. The designed NbS2-based catalysts exhibit an overpotential at a cathodic current of 10 mA cm-2 (n10) as low as 0.10 and 0.22 V vs. RHE in 0.5 M H2SO4 and 1 M KOH, respectively. In 0.5 M H2SO4, the HER activity of the NbS2-based catalysts is also superior to those of the Pt/C benchmark at current densities higher than 80 mA cm-2. Our work provides general guidelines for a scalable and cost-effective exploitation of NbS2, as well as the entire MX2 portfolio, for attaining a viable H2 production through electrochemical routes.

Committee neural network potentials control generalization errors and enable active learning

Abstract: It is well known in the field of machine learning that committee models improve accuracy, provide generalization error estimates, and enable active learning strategies. In this work, we adapt these concepts to interatomic potentials based on artificial neural networks. Instead of a single model, multiple models that share the same atomic environment descriptors yield an average that outperforms its individual members as well as a measure of the generalization error in the form of the committee disagreement. We not only use this disagreement to identify the most relevant configurations to build up the model's training set in an active learning procedure, but also monitor and bias it during simulations to control the generalization error. This facilitates the adaptive development of committee neural network potentials and their training sets, while keeping the number of ab initio calculations to a minimum. To illustrate the benefits of this methodology, we apply it to the development of a committee model for water in the condensed phase. Starting from a single reference ab initio simulation, we use active learning to expand into new state points and to describe the quantum nature of the nuclei. The final model, trained on 814 reference calculations, yields excellent results under a range of conditions, from liquid water at ambient and elevated temperatures and pressures to different phases of ice, and the air-water interface - all including nuclear quantum effects. This approach to committee models will enable the systematic development of robust machine learning models for a broad range of systems.

On the role of gradients for machine learning of molecular energies and forces

Abstract: The accuracy of any machine learning potential can only be as good as the data used in the fitting process. The most efficient model therefore selects the training data that will yield the highest accuracy compared to the cost of obtaining the training data. We investigate the convergence of prediction errors of quantum machine learning models for organic molecules trained on energy and force labels, two common data types in molecular simulations. When training and predicting on different geometries corresponding to the same single molecule, we find that the inclusion of atomic forces in the training data increases the accuracy of the predicted energies and forces 7-fold, compared to models trained on energy only. Surprisingly, for models trained on sets of organic molecules of varying size and composition in non-equilibrium conformations, inclusion of forces in the training does not improve the predicted energies of unseen molecules in new conformations. Predicted forces, however, also improve about 7-fold. For the systems studied, we find that force labels and energy labels contribute equally per label to the convergence of the prediction errors. Choosing to include derivatives such as atomic forces in the training set or not should thus depend on, not only on the computational cost of acquiring the force labels for training, but also on the application domain, the property of interest, and the desirable size of the machine learning model. Based on our observations we describe key considerations for the creation of datasets for potential energy surfaces of molecules which maximize the efficiency of the resulting machine learning models.

Polaritonic normal modes in Transition State Theory

Abstract: A series of experiments demonstrate that strong light-matter coupling between vibrational excitations in isotropic solutions of molecules and resonant infrared optical microcavity modes leads to modified thermally-activated kinetics. However, Feist and coworkers [\emph{Phys. Rev. X.}, \textbf{9}, 021057(2019)] have recently demonstrated that, within transition state theory, the effects of strong light-matter coupling with reactive modes are electrostatic and essentially independent of light-matter resonance or even of the formation of vibrational polaritons. To analyze this puzzling theoretical result in further detail, we revisit it under a new light, invoking a normal mode analysis of the transition state and reactant configurations for an ensemble of an arbitrary number of molecules in a cavity, obtaining simple analytical expressions that produce similar conclusions as Feist. While these effects become relevant in optical microcavities if the molecular dipoles are anisotropically aligned, or in cavities with extreme confinement of the photon modes, they become negligible for isotropic solutions in microcavities. It is concluded that further studies are necessary to track the origin of the experimentally observed kinetics.

Recent developments in MAS DNP-NMR of materials

Abstract: Solid-state NMR spectroscopy is a powerful technique for the characterization of the atomic-level structure and dynamics of materials. Nevertheless, the use of this technique is often limited by its lack of sensitivity, which can prevent the observation of surfaces, defects or insensitive isotopes. Dynamic Nuclear Polarization (DNP) has been shown to improve by one to three orders of magnitude the sensitivity of NMR experiments on materials under Magic-Angle Spinning (MAS), at static magnetic field $B_0 \ge 5$ T, conditions allowing for the acquisition of high-resolution spectra. The field of DNP-NMR spectroscopy of materials has undergone a rapid development in the last ten years, spurred notably by the availability of commercial DNP-NMR systems. We provide here an in-depth overview of MAS DNP-NMR studies of materials at high $B_0$ field. After a historical perspective of DNP of materials, we describe the DNP transfers under MAS, the transport of polarization by spin diffusion and the various contributions to the overall sensitivity of DNP-NMR experiments. We discuss the design of tailored polarizing agents and the sample preparation in the case of materials. We present the DNP-NMR hardware and the influence of key experimental parameters, such as microwave power, magnetic field, temperature and MAS frequency. We give an overview of the isotopes that have been detected by this technique, and the NMR methods that have been combined with DNP. Finally, we show how MAS DNP-NMR has been applied to gain new insights into the structure of organic, hybrid and inorganic materials with applications in fields, such as health, energy, catalysis,optoelectronics etc.

Deep Learning for UV Absorption Spectra with SchNarc: First Steps Towards Transferability in Chemical Compound Space

Abstract: Machine learning (ML) has shown to advance the research field of quantum chemistry in almost any possible direction and has recently also entered the excited states to investigate the multifaceted photochemistry of molecules. In this paper, we pursue two goals: i) We show how ML can be used to model permanent dipole moments for excited states and transition dipole moments by adapting the charge model of [Chem. Sci., 2017, 8, 6924-6935], which was originally proposed for the permanent dipole moment vector of the electronic ground state. ii) We investigate the transferability of our excited-state ML models in chemical space, i.e., whether an ML model can predict properties of molecules that it has never been trained on and whether it can learn the different excited states of two molecules simultaneously. To this aim, we employ and extend our previously reported SchNarc approach for excited-state ML. We calculate UV absorption spectra from excited-state energies and transition dipole moments as well as electrostatic potentials from latent charges inferred by the ML model of the permanent dipole moment vectors. We train our ML models on CH$_2$NH$_2^+$ and C$_2$H$_4$, while predictions are carried out for these molecules and additionally for CHNH$_2$, CH$_2$NH, and C$_2$H$_5^+$. The results indicate that transferability is possible for the excited states.

QM7-X: A comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules

Abstract: We introduce QM7-X, a comprehensive dataset of 42 physicochemical properties for $\approx$ 4.2 M equilibrium and non-equilibrium structures of small organic molecules with up to seven non-hydrogen (C, N, O, S, Cl) atoms. To span this fundamentally important region of chemical compound space (CCS), QM7-X includes an exhaustive sampling of (meta-)stable equilibrium structures - comprised of constitutional/structural isomers and stereoisomers, e.g., enantiomers and diastereomers (including cis-/trans- and conformational isomers) - as well as 100 non-equilibrium structural variations thereof to reach a total of $\approx$ 4.2 M molecular structures. Computed at the tightly converged quantum-mechanical PBE0+MBD level of theory, QM7-X contains global (molecular) and local (atom-in-a-molecule) properties ranging from ground state quantities (such as atomization energies and dipole moments) to response quantities (such as polarizability tensors and dispersion coefficients). By providing a systematic, extensive, and tightly-converged dataset of quantum-mechanically computed physicochemical properties, we expect that QM7-X will play a critical role in the development of next-generation machine-learning based models for exploring greater swaths of CCS and performing in silico design of molecules with targeted properties.

Perspective: Essentials of Relativistic Quantum Chemistry

Abstract: Relativistic quantum chemistry has evolved into a fertile and large field and is now becoming an integrated part of mainstream chemistry. Yet, given the much-involved physics and mathematics (as compared with nonrelativistic quantum chemistry), it is still necessary to clean up the essentials underlying the relativistic electronic structure theories and methodologies (such that uninitiated readers can pick up quickly the right ideas and tools for further development or application) and meanwhile pinpoint future directions of the field. To this end, the three aspects of electronic structure calculations, i.e., relativity, correlation, and QED, will be highlighted.

r2SCAN-D4: Dispersion corrected meta-generalized gradient approximation for general chemical applications

Abstract: We combine a regularized variant of the strongly constrained and appropriately normed semilocal density functional [J. Sun, A. Ruzsinszky, and J. P. Perdew, Phys. Rev. Lett. 115, 036402 (2015)] with the latest generation semi-classical London dispersion correction. The resulting density functional approximation r2SCAN-D4 has the speed of generalized gradient approximations while approaching the accuracy of hybrid functionals for general chemical applications. We demonstrate its numerical robustness in real-life settings and benchmark molecular geometries, general main group and organo-metallic thermochemistry, as well as non-covalent interactions in supramolecular complexes and molecular crystals. Main group and transition metal bond lengths have errors of just 0.8%, which is competitive with hybrid functionals for main group molecules and outperforms them for transition metal complexes. The weighted mean absolute deviation (WTMAD2) on the large GMTKN55 database of chemical properties is exceptionally small at 7.5 kcal/mol. This also holds for metal organic reactions with an MAD of 3.3 kcal/mol. The versatile applicability to organic and metal-organic systems transfers to condensed systems, where lattice energies of molecular crystals are within chemical accuracy (errors <1 kcal/mol).

Quantum Chemistry Simulations of Dominant Products in Lithium-Sulfur Batteries

Abstract: Quantum chemistry simulations of four industrially relevant molecules are reported. Dissociation curves and dipole moments are reported for lithium hydride (LiH), hydrogen sulfide (H2S), lithium hydrogen sulfide (LiSH) and lithium sulfide (Li2S). Herein, we demonstrate the ability to calculate dipole moments using up to 21 qubits on a quantum simulator for a lithium sulfur salt molecule, and demonstrate the ability to calculate the dipole moment of the LiH molecule on the IBM Q Valencia device using four qubits. This is the first example to the best of our knowledge of dipole moment calculations being performed on quantum hardware.

On the Origin of Ground-State Vacuum-Field Catalysis: Equilibrium Consideration

Abstract: Recent experiments suggest that vibrational strong coupling (VSC) may significantly modify ground-state chemical reactions and their rates even without external pumping. The intrinsic mechanism of this "vacuum-field catalysis" remains largely unclear. Generally, modifications of thermal reactions in the ground electronic states can be caused by equilibrium or non-equilibrium effects. The former are associated with modifications of the reactant equilibrium distribution as expressed by the transition state theory of chemical reaction rates, while the latter stem from the dynamics of reaching and leaving transition configurations. Here, we examine the VSC effect in a cavity environment on chemical rates as calculated by transition state theory. Our approach is to examine the effect of coupling to cavity mode(s) on the potential of mean force (PMF) associated with the reaction coordinate. Within the context of classical nuclei and classical photons, we find that while the PMF can be affected by the cavity environment, this effect is negligible for the usual cavities used to examine VSC situations.

Generalization of the hierarchical equations of motion theory for efficient calculations with arbitrary correlation functions

Abstract: The hierarchical equations of motion (HEOM) theory is one of the standard methods to rigorously describe open quantum dynamics coupled to harmonic environments. Such a model is used to capture non-Markovian and non-perturbative effects of environments appearing in ultra-fast phenomena. In the regular framework of the HEOM theory, the environment correlation functions are restricted into linear combinations of exponential functions. In this article, we present a new formulation of the HEOM theory including treatments of non-exponential correlation functions, which enables us to describe general environmental effects more efficiently and stably than the original theory and other generalizations. The library and its Python binding we developed to perform simulations based on our approach, named LibHEOM and PyHEOM respectively, are provided as supplementary material.

Accurate Determination of the Charge Transfer Efficiency of Photoanodes for Solar Water Splitting.

Abstract: The oxygen evolution reaction (OER) at the surface of semiconductor photoanodes involves photo-generated holes that oxidize water. A certain fraction of the holes that reach the surface recombine with electrons from the conduction band, giving rise to the surface recombination loss. The charge transfer efficiency, xt, defined as the ratio between the flux of holes that contribute to the water oxidation reaction and the total flux of holes that reach the surface, is an important parameter that helps to distinguish between bulk and surface recombination losses. However, accurate determination of xt by conventional voltammetry measurements is complicated because only the total current is measured and it is difficult to discern between different contributions to the current. Chopped light measurement and hole scavenger measurement techniques are widely employed to determine xt, but they often lead to errors. Intensity modulated photocurrent spectroscopy (IMPS) is better suited for accurate determination of xt because it provides direct information on both the total photocurrent and the surface recombination current. Careful analysis of IMPS measurements at different light intensities is required to account for nonlinear effects. We compare the xt values obtained by these methods using heteroepitaxial hematite photoanodes. A wide spread of xt values is obtained by different analysis methods and different light sources and light intensities. Statistical analysis of the results show good correlation between different methods for measurements carried out with the same light source, light intensity and potential. However, there is a considerable spread in the results obtained by different methods. For accurate determination of xt, we recommend IMPS measurements with a bias light intensity such that the irradiance is as close as possible to the standard solar spectrum.

Spaced TiO$_2$ Nanotube Arrays Allow for High Performance Hierarchical Supercapacitor Structure.

Abstract: In this work we describe the synthesis and electrochemical properties of nitridated hierarchical TiO$_2$ nanotubes as an electrode for supercapacitors. The hierarchical TiO$_2$ nanostructures are formed by a controlled layer-by-layer TiO2 nanoparticle decoration on self-organized spaced TiO$_2$ nanotubes. These structures are then annealed in NH$_3$ atmosphere at elevated temperature to convert the material to a nitride structure -- this drastically enhances their electron-transport properties. The areal capacitance of hierarchical structures can be tuned by changing the number of decorated TiO$_2$ nanoparticle layers. The capacitance enhancement of the hierarchical structures reaches a maximum when the surface area through nanoparticle deposition is highest and the conductivity via nitridation is optimized.

Metrics for Benchmarking and Uncertainty Quantification: Quality, Applicability, and a Path to Best Practices for Machine Learning in Chemistry

Abstract: This review aims to draw attention to two issues of concern when we set out to make machine learning work in the chemical and materials domain, i.e., statistical loss function metrics for the validation and benchmarking of data-derived models, and the uncertainty quantification of predictions made by them. They are often overlooked or underappreciated topics as chemists typically only have limited training in statistics. Aside from helping to assess the quality, reliability, and applicability of a given model, these metrics are also key to comparing the performance of different models and thus for developing guidelines and best practices for the successful application of machine learning in chemistry.

Interactions between Large Molecules: Puzzle for Reference Quantum-Mechanical Methods

Abstract: Quantum-mechanical methods are widely used for understanding molecular interactions throughout biology, chemistry, and materials science. Quantum diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] are two state-of-the-art and trusted wavefunction methods that have been categorically shown to yield accurate interaction energies for small organic molecules. These methods provide valuable reference information for widely-used semi-empirical and machine learning potentials, especially where experimental information is scarce. However, agreement for systems beyond small molecules is a crucial remaining milestone for cementing the benchmark accuracy of these methods. Approaching such well-converged predictive power in larger molecules has motivated major developments in CCSD(T) as well as DMC algorithms in the past years, resulting in orders of magnitude time-to-solution reductions. Here, we show that CCSD(T) and DMC interaction energies are not in consistent agreement for a set of polarizable supramolecules. Whilst agreement is found for some of the complexes, in a few key systems disagreements of up to 8 kcal/mol remained. This leads to differences of up to 6 orders of magnitude in the corresponding binding association constant at room temperature for systems which are well within the accustomed domain of applicability for both methods. These findings thus indicate that more caution is required when aiming at reproducible non-covalent interactions between extended molecules. Our data contradicts the expectation that the most comprehensive and robust wavefunction methods predict identical non-covalent interactions and indicate an unsolved challenge for benchmark approaches.