Search

In recent years there has been a rapid growth in the development and application of new stochastic methods in electronic structure. These methods are quite diverse, from many-body wave function techniques in real space or determinant space to being used to sum perturbative expansions. This growth has been spurred by the more favorable scaling with the number of electrons and often better parallelization over large numbers of central processing unit (CPU) cores or graphical processing units (GPUs) than for high-end non-stochastic wave function based methods. This special issue of the Journal of Chemical Physics includes 33 papers that describe recent developments and applications in this area. As seen from the articles in the issue, stochastic electronic structure methods are applicable to both molecules and solids and can accurately describe systems with strong electron correlation. This issue was motivated, in part, by the 2019 Telluride Science Research Center workshop on Stochastic Electronic Structure Methods that we organized. Below we briefly describe each of the papers in the special issue, dividing the papers into six subtopics. [ABSTRACT FROM AUTHOR]

To address the challenge of performance portability and facilitate the implementation of electronic structure solvers, we developed the basic matrix library (BML) and Parallel, Rapid O(N), and Graph-based Recursive Electronic Structure Solver (PROGRESS) library. The BML implements linear algebra operations necessary for electronic structure kernels using a unified user interface for various matrix formats (dense and sparse) and architectures (CPUs and GPUs). Focusing on density functional theory and tight-binding models, PROGRESS implements several solvers for computing the single-particle density matrix and relies on BML. In this paper, we describe the general strategies used for these implementations on various computer architectures, using OpenMP target functionalities on GPUs, in conjunction with third-party libraries to handle performance critical numerical kernels. We demonstrate the portability of this approach and its performance in benchmark problems. [ABSTRACT FROM AUTHOR]

Dilute magnetic semiconductors (DMSs) with both charge and spin degrees of freedom have emerged as promising candidates in the spintronic industry. However, the Curie temperature below room temperature and uncertainty about the origin of ferromagnetism hinder the application of DMSs. To address these issues, we explored a better SnO2-based co-doped method (Fe, V + VSn) using ab initio calculations. The calculation results show that the Sn13FeVO32 (Fe, V + VSn) has a high Curie temperature (716 K), good ferromagnetic properties, stronger covalency of bonds, and better optical transparency in the visible light range. In addition, the holes or electrons generated by the complexes in the (Fe, V + VO/VSn) co-doped system cause a spin-polarized double exchange effect in the Fe-3d, V-3d, and O-2p orbitals, which leads to magnetism of the co-doped systems. The static dielectric constant ɛ1(0) of the system increases after doping. Among them, Sn14FeVO31 (Fe, V + VO) has the largest ɛ1(0), indicating that Sn14FeVO31 has the strongest polarization ability and better photocatalytic properties. In Sn14FeVO31, the imaginary part of the dielectric function and the absorption spectrum all have new peaks in the low-energy region, which are caused by the jump of electrons from the guide band of the spin-polarized impurity energy level. This paper proposes a new method for preparing dilute magnetic semiconductors in spin electronic devices with high room temperature ferromagnetic properties and excellent optical properties through the (Fe, V + VO/VSn) co-doped SnO2. [ABSTRACT FROM AUTHOR]

In this paper, a global and full-dimensional potential energy surface at the 2A″ ground state for the Al + O2 → AlO + O reaction was constructed, for the first time, based on extensive electronic structure calculations using the doubly hybrid functional XYG3 and potential energy surface fittings by neural networks. Details of the reaction paths have been analyzed. It was found that both two intermediates, the cyclic-AlO2 and the linear-OAlO, were able to dissociate to the AlO + O products, and the isomerization process between these two intermediates was controlled by conical intersections between two 2A″ states. Ro-vibrational state resolved integral cross sections have also been calculated at collision energies from 1.0 to 10.0 kcal/mol. The results support the harpooning mechanism in this metal-oxidant-involved reaction. [ABSTRACT FROM AUTHOR]

Gibbs energies for reactions involving aqueous ions are challenging to predict due to the large solvation energies of such ions. A stringent test would be the ab initio reproduction of the aqueous-phase chelate effect, an entropic effect in reactions of very small enthalpy changes. This paper examines what is required to achieve such a reproduction for the paradigmatic reaction M(NH3)42+ + 2 en → M(en)22+ + 4 NH3 (en = 1,2-ethylenediamine), for which ΔrxnG* and ΔrxnH* are −2.3 and +1.6 kcal mol−1, respectively, if M = Zn. Explicit solvation via simulation was avoided in order to allow sufficiently accurate electronic structure models; this required the use of continuum solvation models (CSMs), and a great deal of effort was made in attempting to lower the relative errors of ΔsolvG*[M(NH3)42+] vs ΔsolvG*[M(en)22+] from the CSMs available in Gaussian software. CSMs in ADF and JDFTx software were also tested. A uniform 2.2 kcal mol−1 accuracy in ΔrxnG* for all three metal-atom choices M = {Zn, Cd, Hg} was eventually achieved, but not from any of the known CSMs tested, nor from cavity size reoptimization, nor from semicontinuum modeling: post facto solvation energy corrections [one per solute type, NH3, en, M(NH3)42+, M(en)22+] were needed. It is hoped that this study will aid (and encourage) further CSM development for coordination-complex ions. [ABSTRACT FROM AUTHOR]

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). In this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H2Z2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces. [ABSTRACT FROM AUTHOR]

The Pauli kinetic energy functional and its functional derivative, termed Pauli potential, play a crucial role in the successful implementation of orbital-free density functional theory for electronic structure calculations. However, the exact forms of these two quantities are not known. Therefore, perforce, one employs the approximate forms for the Pauli functional or Pauli potential for performing orbital-free density functional calculations. In the present study, we developed a feed-forward neural network-based representation for the Pauli potential using a 1-dimensional (1-D) model system. We expanded density in terms of basis functions, and the coefficients of the expansion were used as input to a feed-forward neural network. Using the neural network-based representation of the Pauli potential, we calculated the ground-state densities of the 1-D model system by solving the Euler equation. We calculated the Pauli kinetic energy using the neural network-based Pauli potential employing the exact relation between the Pauli kinetic energy functional and the potential. The sum of the neural network-based Pauli kinetic energy and the von Weizsäcker kinetic energy resulted in an accurate estimation of the total kinetic energy. The approach presented in this paper can be employed for the calculation of Pauli potential and Pauli kinetic energy, obviating the need for a functional derivative. The present study is an important step in the advancement of application of machine learning-based techniques toward the orbital-free density functional theory-based methods. [ABSTRACT FROM AUTHOR]

With a He– H 3 + interaction potential obtained from advanced electronic structure calculations, we computed the vibration-rotation-tunneling (VRT) states of this complex for total angular momenta J from 0 to 9, both for the vibrational ground state and for the twofold degenerate v2 = 1 excited state of H 3 + . The potential has three equivalent global minima with depth De = 455.3 cm−1 for He in the plane of H 3 + , three equatorial saddle points that separate these minima with barriers of 159.5 cm−1, and two axial saddle points with energies of 243.1 cm−1 above the minima. The dissociation energies calculated for the complexes of He with ortho- H 3 + (o H 3 + ) and para- H 3 + (p H 3 + ) are D0 = 234.5 and 236.3 cm−1, respectively. Wave function plots of the VRT states show that they may be characterized as weakly hindered internal rotor states, delocalized over the three minima in the potential and with considerable amplitude at the barriers. Most of them are dominated by the jk = 10 and 11 rotational ground states of o H 3 + and p H 3 + , with the intermolecular stretching mode excited up to v = 4 inclusive. However, we also found excited internal rotor states: 33 in He–o H 3 + , and 22 and 21 in He–p H 3 + . The VRT levels and wave functions were used to calculate the frequencies and line strengths of all allowed v2 = 0 → 1 rovibrational transitions in the complex. Theoretical spectra generated with these results are compared with the experimental spectra in Paper II [Salomon et al., J. Chem. Phys. 156, 144308 (2022)] and are extremely helpful in assigning these spectra. This comparison shows that the theoretical energy levels and spectra agree very well with the measured ones, which confirms the high accuracy of our ab initio He– H 3 + interaction potential and of the ensuing calculations of the VRT states. [ABSTRACT FROM AUTHOR]

Recent work has demonstrated the promise of using machine-learned surrogates, in particular, Gaussian process (GP) surrogates, in reducing the number of electronic structure calculations (ESCs) needed to perform surrogate model based (SMB) geometry optimization. In this paper, we study geometry meta-optimization with GP surrogates where a SMB optimizer additionally learns from its past "experience" performing geometry optimization. To validate this idea, we start with the simplest setting where a geometry meta-optimizer learns from previous optimizations of the same molecule with different initial-guess geometries. We give empirical evidence that geometry meta-optimization with GP surrogates is effective and requires less tuning compared to SMB optimization with GP surrogates on the ANI-1 dataset of off-equilibrium initial structures of small organic molecules. Unlike SMB optimization where a surrogate should be immediately useful for optimizing a given geometry, a surrogate in geometry meta-optimization has more flexibility because it can distribute its ESC savings across a set of geometries. Indeed, we find that GP surrogates that preserve rotational invariance provide increased marginal ESC savings across geometries. As a more stringent test, we also apply geometry meta-optimization to conformational search on a hand-constructed dataset of hydrocarbons and alcohols. We observe that while SMB optimization and geometry meta-optimization do save on ESCs, they also tend to miss higher energy conformers compared to standard geometry optimization. We believe that further research into characterizing the divergence between GP surrogates and potential energy surfaces is critical not only for advancing geometry meta-optimization but also for exploring the potential of machine-learned surrogates in geometry optimization in general. [ABSTRACT FROM AUTHOR]

Electronic structure calculations have the potential to predict key matter transformations for applications of strategic technological importance, from drug discovery to material science and catalysis. However, a predictive physicochemical characterization of these processes often requires accurate quantum chemical modeling of complex molecular systems with hundreds to thousands of atoms. Due to the computationally demanding nature of electronic structure calculations and the complexity of modern high-performance computing hardware, quantum chemistry software has historically failed to operate at such large molecular scales with accuracy and speed that are useful in practice. In this paper, novel algorithms and software are presented that enable extreme-scale quantum chemistry capabilities with particular emphasis on exascale calculations. This includes the development and application of the multi-Graphics Processing Unit (GPU) library LibCChem 2.0 as part of the General Atomic and Molecular Electronic Structure System package and of the standalone Extreme-scale Electronic Structure System (EXESS), designed from the ground up for scaling on thousands of GPUs to perform high-performance accurate quantum chemistry calculations at unprecedented speed and molecular scales. Among various results, we report that the EXESS implementation enables Hartree–Fock/cc-pVDZ plus RI-MP2/cc-pVDZ/cc-pVDZ-RIFIT calculations on an ionic liquid system with 623 016 electrons and 146 592 atoms in less than 45 min using 27 600 GPUs on the Summit supercomputer with a 94.6% parallel efficiency. [ABSTRACT FROM AUTHOR]

Heavy-atom alkaline earth monofluoride molecules are considered as prospective systems to study spatial parity or spatial parity and time-reversal symmetry violating effects such as the nuclear anapole moment or the electron electric dipole moment. A comprehensive and highly accurate theoretical study of the electronic structure properties and transition energies in such systems can simplify the preparation and interpretation of the experiments. However, almost no attempts to calculate quantum electrodynamics (QED) effects' contribution into characteristics of these neutral heavy-atom molecules have been performed. Recently, we have formulated and implemented such an approach to calculate QED contributions to transition energies of molecules [L. V. Skripnikov, J. Chem. Phys. 154, 201101 (2021)]. In this paper, we perform a benchmark theoretical study of the transition energies in the Ba+ cation and BaF molecule. The deviation of the calculated values from the experimental ones is of the order 10 cm−1 and is more than an order of magnitude better than the "chemical accuracy," 350 cm−1. The achievement of such an agreement has been provided, in particular, by the inclusion of the QED effects. The latter appeared to be not less important than the high-order correlation effects beyond the coupled cluster with single, double, and perturbative triple cluster amplitude level. We compare the role of QED effects for transition energies with heavier molecules—RaF and E120F, where E120 is the superheavy Z = 120 homolog of Ra. [ABSTRACT FROM AUTHOR]

In this paper, we present CTRAMER (Charge-Transfer RAtes from Molecular dynamics, Electronic structure, and Rate theory)—an open-source software package for calculating interfacial charge-transfer (CT) rate constants in organic photovoltaic (OPV) materials based on ab initio calculations and molecular dynamics simulations. The software is based on identifying representative donor/acceptor geometries within interfacial structures obtained from molecular dynamics simulation of donor/acceptor blends and calculating the corresponding Fermi's golden rule CT rate constants within the framework of the linearized-semiclassical approximation. While the methods used are well established, the integration of these state-of-the-art tools originating from different disciplines to study photoinduced CT processes with explicit treatment of the environment, in our opinion, makes this package unique and innovative. The software also provides tools for investigating other observables of interest. After outlining the features and implementation details, the usage and performance of the software are demonstrated with results from an example OPV system. [ABSTRACT FROM AUTHOR]

It is important that any dynamics method approaches the correct population distribution at long times. In this paper, we derive a one-body reduced density matrix dynamics for electrons in energetic contact with a bath. We obtain a remarkable equation of motion which shows that in order to reach equilibrium properly, rates of electron transitions depend on the density matrix. Even though the bath drives the electrons towards a Boltzmann distribution, hole blocking factors in our equation of motion cause the electronic populations to relax to a Fermi-Dirac distribution. These factors are an old concept, but we show how they can be derived with a combination of time-dependent perturbation theory and the extended normal ordering of Mukherjee and Kutzelnigg for a general electronic state. The resulting non-equilibrium kinetic equations generalize the usual Redfield theory to many-electron systems, while ensuring that the orbital occupations remain between zero and one. In numerical applications of our equations, we show that relaxation rates of molecules are not constant because of the blocking effect. Other applications to model atomic chains are also presented which highlight the importance of treating both dephasing and relaxation. Finally, we show how the bath localizes the electron density matrix. [ABSTRACT FROM AUTHOR]

In the full optimized reaction space and valence-complete active space self-consistent field (vCAS) methods, a set of active orbitals is defined as the union of the valence orbitals on the atoms, all possible configurations involving the active orbitals are generated, and the orbitals and configuration coefficients are self-consistently optimized. Such wave functions have tremendous flexibility, which makes these methods incredibly powerful but can also lead to inconsistencies in the description of the electronic structure of molecules. In this paper, the problems that can arise in vCAS calculations are illustrated by calculations on the BH and BF molecules. BH is well described by the full vCAS wave function, which accounts for molecular dissociation and 2s–2p near-degeneracy in the boron atom. The same is not true for the full vCAS wave function for BF. There is mixing of core and active orbitals at short internuclear distances and swapping of core and active orbitals at large internuclear distances. In addition, the virtual 2π orbitals, which were included in the active space to account for the 2s–2p near degeneracy effect, are used instead to describe radial correlation of the electrons in the F2pπ-like pairs. Although the above changes lead to lower vCAS energies, they lead to higher vCAS+1+2 energies as well as irregularities and/or discontinuities in the potential energy curves. All of the above problems can be addressed by using the spin-coupled generalized valence bond-inspired vCAS wave function for BF, which includes only a subset of the atomic valence orbitals in the active space. [ABSTRACT FROM AUTHOR]

The nonequilibrium Fermi's golden rule describes the transition between a photoexcited bright donor electronic state and a dark acceptor electronic state, when the nuclear degrees of freedom start out in a nonequilibrium state. In a previous paper [X. Sun and E. Geva, J. Chem. Theory Comput. 12, 2926 (2016)], we proposed a new expression for the nonequilibrium Fermi's golden rule within the framework of the linearized semiclassical approximation and based on the Condon approximation, according to which the electronic coupling between donor and acceptor is assumed constant. In this paper we propose a more general expression, which is applicable to the case of non-Condon electronic coupling. We test the accuracy of the new non-Condon nonequilibrium Fermi's golden rule linearized semiclassical expression on a model where the donor and acceptor potential energy surfaces are parabolic and identical except for shifts in the equilibrium energy and geometry, and the coupling between them is linear in the nuclear coordinates. Since non-Condon effects may or may not give rise to conical intersections, both possibilities are examined by considering the following: (1) A modified Garg-Onuchic-Ambegaokar model for charge transfer in the condensed phase, where the donor-acceptor coupling is linear in the primary-mode coordinate, and for which non-Condon effects do not give rise to a conical intersection; (2) the linear vibronic coupling model for electronic transitions in gas phase molecules, where non-Condon effects give rise to conical intersections. We also present a comprehensive comparison between the linearized semiclassical expression and a progression of more approximate expressions, in both normal and inverted regions, and over a wide range of initial nonequilibrium states, temperatures, and frictions. [ABSTRACT FROM AUTHOR]

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. This paper explains the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science. [ABSTRACT FROM AUTHOR]

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-principle-driven methodologies to model complex chemical and materials processes. Over the past 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. [ABSTRACT FROM AUTHOR]

The vibrational spectroscopy of adsorbates is becoming an important investigation tool for catalysis and material science. This paper presents a semiclassical molecular dynamics method able to reproduce the vibrational energy levels of systems composed by molecules adsorbed on solid surfaces. Specifically, we extend our divide-and-conquer semiclassical method for power spectra calculations to gas-surface systems and interface it with plane-wave electronic structure codes. The Born–Oppenheimer classical dynamics underlying the semiclassical calculation is full dimensional, and our method includes not only the motion of the adsorbate but also those of the surface and the bulk. The vibrational spectroscopic peaks related to the adsorbate are accounted together with the most coupled phonon modes to obtain spectra amenable to physical interpretations. We apply the method to the adsorption of CO, NO, and H2O on the anatase-TiO2 (101) surface. We compare our semiclassical results with the single-point harmonic estimates and the classical power spectra obtained from the same trajectory employed in the semiclassical calculation. We find that CO and NO anharmonic effects of fundamental vibrations are similarly reproduced by the classical and semiclassical dynamics and that H2O adsorption is fully and properly described in its overtone and combination band relevant components only by the semiclassical approach. [ABSTRACT FROM AUTHOR]

We develop several configuration interaction approaches for characterizing the electronic structure of an adsorbate on a metal surface (at least in model form). When one can separate the adsorbate from the substrate, these methods can achieve a reasonable description of adsorbate on-site electron-electron correlation in the presence of a continuum of states. While the present paper is restricted to the Anderson impurity model, there is hope that these methods can be extended to ab initio Hamiltonians and provide insight into the structure and dynamics of molecule-metal surface interactions. [ABSTRACT FROM AUTHOR]

There are many ways to numerically represent chemical systems in order to compute their electronic structure. Basis functions may be localized in real-space (atomic orbitals), in momentum-space (plane waves), or in both components of phase-space. Such phase-space localized basis functions in the form of wavelets have been used for many years in the electronic structure. In this paper, we turn to a phase-space localized basis set first introduced by Wilson [Generalized Wannier Functions (Cornell University, 1987), preprint]. We provide the first full study of this basis and its numerical implementation. To calculate electronic energies of a variety of small molecules and states, we utilize the sum-of-products form, Gaussian quadratures, and introduce methods for selecting sample points from a grid of phase-space localized Wilson basis. Both full configuration interaction and Hartree-Fock implementations are discussed and implemented numerically. As with many grid based methods, describing both tightly bound and diffuse orbitals is challenging so we have considered augmenting the Wilson basis set as projected Slater-type orbitals. We have also compared the Wilson basis set against the recently introduced wavelet transformed Gaussians (gausslets). Throughout, we give comments on the implementation and use small atoms and molecules to illustrate the convergence properties of the Wilson basis. [ABSTRACT FROM AUTHOR]

This paper presents the nonadiabatic molecular dynamics simulation combined with the spin-flip time-dependent density functional theory (SF-TDDFT). In contrast to the conventional single-reference electronic structure methods, which have difficulty in describing the S0/S1 conical intersections, the SF-TDDFT can yield the correct topology of crossing points. Thus, one expects that the method can take naturally into account the S1 → S0 nonadiabatic transitions. We adopt Tully's fewest switch surface hopping algorithm by introducing the analytic SF-TDDFT nonadiabatic coupling vector. We apply the proposed method to the photoisomerization reactions of E-azomethane, methanimine, and ethene molecules and reproduce the results of previous studies based on the multireference methods. The proposed approach overcomes the ad hoc treatment of S1 → S0 transition at the single-reference calculation level and affords both the dynamics on the S1 state and the recovery of the S0 state with modest computational costs. [ABSTRACT FROM AUTHOR]

The bulk properties such as lattice constants, bulk moduli, and cohesive energies of alkali, alkaline-earth, and transition metals are studied within the framework of the recently developed meta-GGA (meta-Generalized Gradient Approximation) level semilocal exchange-correlation functionals. To establish the applicability, broadness, and accuracy of meta-GGA functionals, we also put the results obtained using PBE (Perdew-Burke-Ernzerhof) [J. P. Perdew et al., Phys. Rev. Lett. 77, 3865 (1996)] and PBE reparameterized for solid [J. P. Perdew et al., Phys. Rev. Lett. 100, 136406 (2008)] GGA functionals. The interesting feature of the present paper is that it measures the accuracy of the recently developed TM (Tao-Mo), TMTPSS [TM exchange with Tao-Perdew-Staroverov-Scuseria (TPSS)] [J. Tao and Y. Mo, Phys. Rev. Lett. 117, 073001 (2016)] correlation, and strongly constrained and appropriately normed [J. Sun et al., Phys. Rev. Lett. 115, 036402 (2015)] functionals to calculate the aforementioned properties. Not only that, we also include other (popular) meta-GGA functionals in order to have a closer look at the performance of the meta-GGA functionals too. The present systematic investigation shows that the TM functional is accurate in describing the lattice constants while for cohesive energies and bulk moduli, the PBE and modified TPSS perform better compared to others. [ABSTRACT FROM AUTHOR]

Density functional theory (DFT) is the most successful and widely used approach for computing the electronic structure of matter. However, for tasks involving large sets of candidate molecules, running DFT separately for every possible compound of interest is forbiddingly expensive. In this paper, we propose a neural network based machine learning algorithm which, assuming a sufficiently large training sample of actual DFT results, can instead learn to predict certain properties of molecules purely from their molecular graphs. Our algorithm is based on the recently proposed covariant compositional networks framework and involves tensor reduction operations that are covariant with respect to permutations of the atoms. This new approach avoids some of the representational limitations of other neural networks that are popular in learning from molecular graphs and yields promising results in numerical experiments on the Harvard Clean Energy Project and QM9 molecular datasets. [ABSTRACT FROM AUTHOR]

Nanoplates such as graphene and MoS2 are promising materials due to their excellent electronic and mechanical properties. The preparation of such nanoplates is, however, still challenging due to the large free energy barrier that multilayer nanoplates need to overcome during exfoliation. In the case of a Bernal-stacked bilayer graphene, the binding energy between two graphene layers is about 17.8 meV per atom such that harsh chemical and/or mechanical treatment is usually necessary. In this paper, we perform extensive molecular dynamics simulations for a generic model of nanoplates and illustrate that when the shear is applied to the nanoplate bilayer solution, the nanoplate bilayer may exfoliate readily. In our simulations, the free energy barrier that two nanoplate layers need to overcome reaches up to 21.8kBT, where kB and T denote the Boltzmann constant and temperature, respectively. This implies that without external stimuli, the nanoplate bilayer would hardly exfoliate. Upon the application of shear, however, the transition between different Bernal stacked conformations occurs, which provides multiple intermediate states for exfoliation and facilitates the shear exfoliation. We also find that if one were to increase the affinity between the solvent and nanoplates slightly, the free energy barrier would be decreased significantly. [ABSTRACT FROM AUTHOR]

Photodetachment spectroscopy of H2B-7 and its deuterated isotopomer probing the energetically low-lying electronic states of the respective neutral cluster is theoretically investigated in this paper. The theoretical methodology is based on detailed quantum chemistry calculations of electronic state energies, construction of a vibronic coupling model in the diabatic electronic basis, and nuclear dynamics calculations from first principles using time-dependent and time-independent quantum mechanical methods. The theoretical model consists of five coupled electronic states and fifteen vibrational modes. Several reduced dimensional calculations are performed to identify the relevant vibrational modes contributing to the vibronic structure of electronic bands and the impact of non-adiabatic coupling on them. The low-energy part of the spectrum of both H2B7 and its deuterated analogue is assigned by examining the vibronic wavefunctions and the results are compared with the experimental findings. The nonadiabatic decay dynamics of the electronic excited states of the neutral clusters is examined at length. [ABSTRACT FROM AUTHOR]

Provided the initial state, the Runge-Gross theorem establishes that the time-dependent (TD) external potential of a system of non-relativistic electrons determines uniquely their TD electronic density, and vice versa (up to a constant in the potential). This theorem requires the TD external potential and density to be Taylor-expandable around the initial time of the propagation. This paper presents an extension without this restriction. Given the initial state of the system and evolution of the density due to some TD scalar potential, we show that a perturbative (not necessarily weak) TD potential that induces a non-zero divergence of the external force-density, inside a small spatial subset and immediately after the initial propagation time, will cause a change in the density within that subset, implying that the TD potential uniquely determines the TD density. In this proof, we assume unitary evolution of wavefunctions and first-order differentiability (which does not imply analyticity) in time of the internal and external force-densities, electronic density, current density, and their spatial derivatives over the small spatial subset and short time interval. [ABSTRACT FROM AUTHOR]

In this paper, we employ an evolutionary algorithm along with the full-potential density functional theory (DFT) computations to perform a comprehensive search for the stable structures of stoichiometric (WS2)n nano-clusters (n = 1 - 9), within three different exchange-correlation functionals. Our results suggest that n = 5 and 8 are possible candidates for the low temperature magic sizes of WS2 nano-clusters while at temperatures above 500 Kelvin, n = 7 exhibits a comparable relative stability with n = 8. The electronic properties and energy gap of the lowest energy isomers were computed within several schemes, including semilocal Perdew-Burke-Ernzerhof and Becke- Lee-Yang-Parr functionals, hybrid B3LYP functional, many body based DFT+GW approach, ΔSCF method, and time dependent DFT calculations. Vibrational spectra of the lowest lying isomers, computed by the force constant method, are used to address IR spectra and thermal free energy of the clusters. Time dependent density functional calculation in a real time domain is applied to determine the full absorption spectra and optical gap of the lowest energy isomers of the WS2 nano-clusters. [ABSTRACT FROM AUTHOR]

Accurate and computationally efficient exchange-correlation functionals are critical to the successful application of linear-scaling density functional theory (DFT). Local and semi-local functionals of the density are naturally compatible with linear-scaling approaches, having a general form which assumes the locality of electronic interactions and which can be efficiently evaluated by numerical quadrature. Presently, the most sophisticated and flexible semi-local functionals are members of the meta-generalized-gradient approximation (meta-GGA) family, and depend upon the kinetic energy density, τ, in addition to the charge density and its gradient. In order to extend the theoretical and computational advantages of τ-dependent meta-GGA functionals to large-scale DFT calculations on thousands of atoms, we have implemented support for τ-dependent meta-GGA functionals in the ONETEP program. In this paper we lay out the theoretical innovations necessary to implement τ-dependent meta-GGA functionals within ONETEP's linear-scaling formalism. We present expressions for the gradient of the τ-dependent exchange-correlation energy, necessary for direct energy minimization. We also derive the forms of the τ-dependent exchange-correlation potential and kinetic energy density in terms of the strictly localized, self-consistently optimized orbitals used by ONETEP. To validate the numerical accuracy of our self-consistent meta-GGA implementation, we performed calculations using the B97M-V and PKZB meta-GGAs on a variety of small molecules. Using only a minimal basis set of self-consistently optimized local orbitals, we obtain energies in excellent agreement with large basis set calculations performed using other codes. Finally, to establish the linear-scaling computational cost and applicability of our approach to large-scale calculations, we present the outcome of self-consistent meta-GGA calculations on amyloid fibrils of increasing size, up to tens of thousands of atoms. [ABSTRACT FROM AUTHOR]

We investigate the sensitivity of X-ray absorption spectra, simulated using a general method, to properties of molecular excited states. Recently, Averbukh and co-workers [M. Ruberti et al., J. Chem. Phys. 140, 184107 (2014)] introduced an efficient and accurate L² method for the calculation of excited state valence photoionization cross-sections based on the application of Stieltjes imaging to the Lanczos pseudo-spectrum of the algebraic diagrammatic construction (ADC) representation of the electronic Hamiltonian. In this paper, we report an extension of this method to the calculation of excited state core photoionization cross-sections. We demonstrate that, at the ADC(2)x level of theory, ground state X-ray absorption spectra may be accurately reproduced, validating the method. Significantly, the calculated X-ray absorption spectra of the excited states are found to be sensitive to both geometric distortions (structural dynamics) and the electronic character (electronic dynamics) of the initial state, suggesting that core excitation spectroscopies will be useful probes of excited state non-adiabatic dynamics. We anticipate that the method presented here can be combined with ab initio molecular dynamics calculations to simulate the time-resolved X-ray spectroscopy of excited state molecular wavepacket dynamics. [ABSTRACT FROM AUTHOR]

Previous work has shown how a symmetrical quasi-classical (SQC) windowing procedure can be used to quantize the initial and final electronic degrees of freedom in the Meyer-Miller (MM) classical vibronic (i.e, nuclear + electronic) Hamiltonian, and that the approach provides a very good description of electronically non-adiabatic processes within a standard classical molecular dynamics framework for a number of benchmark problems. This paper explores application of the SQC/MM approach to the case of very weak non-adiabatic coupling between the electronic states, showing (as anticipated) how the standard SQC/MM approach used to date fails in this limit, and then devises a new SQC windowing scheme to deal with it. Application of this new SQC model to a variety of realistic benchmark systems shows that the new model not only treats the weak coupling case extremely well, but it is also seen to describe the "normal" regime (of electronic transition probabilities ≳ 0.1) even more accurately than the previous "standard" model. [ABSTRACT FROM AUTHOR]

Double-hybrid density functional approximations (DH-DFAs) provide an accurate description of the electronic structure of molecules by semiempirically mixing density functional and wavefunction theory. In this paper, we investigate the properties of the potential used in such approximations. By using the optimized effective potential approach, the consistent Kohn-Sham (KS) potential for a double-hybrid functional (including the second-order perturbational contribution) can be generated. This potential is shown to provide an improved description of orbital energies as vertical ionization potentials (IPs), relative to the perturbation-free KS potential typically used. Based on this observation, we suggest that DH-DFAs should be constructed in such a way that the potential provides accurate orbital energies. As a proof of principle, the B2-PLYP functional is reparameterized to obtain the IP-optimized B2IP-PLYP functional, using a small set of vertical IPs and atomization energies as reference data. This functional is shown to outperform B2-PLYP in a wide range of benchmarks and is en par with the related B2GP-PLYP. In particular, it is shown to be the most reliable choice in electronically difficult and multireference cases. [ABSTRACT FROM AUTHOR]

This paper addresses the definition of correlation energy within 4-component relativistic atomic and molecular calculations. In the nonrelativistic domain the correlation energy is defined as the difference between the exact eigenvalue of the electronic Hamiltonian and the Hartree-Fock energy. In practice, what is reported is the basis set correlation energy, where the "exact" value is provided by a full Configuration Interaction (CI) calculation with some specified one-particle basis. The extension of this definition to the relativistic domain is not straightforward since the corresponding electronic Hamiltonian, the Dirac-Coulomb Hamiltonian, has no bound solutions. Present-day relativistic calculations are carried out within the no-pair approximation, where the Dirac-Coulomb Hamiltonian is embedded by projectors eliminating the troublesome negative-energy solutions. Hartree-Fock calculations are carried out with the implicit use of such projectors and only positive-energy orbitals are retained at the correlated level, meaning that the Hartree-Fock projectors are frozen at the correlated level. We argue that the projection operators should be optimized also at the correlated level and that this is possible by full Multiconfigurational Self-Consistent Field (MCSCF) calculations, that is, MCSCF calculations using a no-pair full CI expansion, but including orbital relaxation from the negative-energy orbitals. We show by variational perturbation theory that the MCSCF correlation energy is a pure MP2-like correlation expression, whereas the corresponding CI correlation energy contains an additional relaxation term. We explore numerically our theoretical analysis by carrying out variational and perturbative calculations on the two-electron rare gas atoms with specially tailored basis sets. In particular, we show that the correlation energy obtained by the suggested MCSCF procedure is smaller than the no-pair full CI correlation energy, in accordance with the underlying minmax principle and our theoretical analysis. We also show that the relativistic correlation energy, obtained from no-pair full MCSCF calculations, scales at worst as X-2 with respect to the cardinal number X of our correlation-consistent basis sets optimized for the two-electron atoms. This is better than the X-1 scaling suggested by previous studies, but worse than the X-3 scaling observed in the nonrelativistic domain. The well-known 1/Z- expansion in nonrelativistic atomic theory follows from coordinate scaling. We point out that coordinate scaling for consistency should be accompanied by velocity scaling. In the nonrelativistic domain this comes about automatically, whereas in the relativistic domain an explicit scaling of the speed of light is required. This in turn explains why the relativistic correlation energy to the lowest order is not independent of nuclear charge, in contrast to nonrelativistic theory. [ABSTRACT FROM AUTHOR]

Quantum molecular dynamics requires an accurate representation of the molecular potential energy surface from a minimal number of electronic structure calculations, particularly for nonadiabatic dynamics where excited states are required. In this paper, we employ pseudospectral sampling of time-dependent Gaussian basis functions for the simulation of non-adiabatic dynamics. Unlike other methods, the pseudospectral Gaussian molecular dynamics tests the Schrödinger equation with N Dirac delta functions located at the centers of the Gaussian functions reducing the scaling of potential energy evaluations from O(N²) to O(N). By projecting the Gaussian basis onto discrete points in space, the method is capable of efficiently and quantitatively describing the nonadiabatic population transfer and intra-surface quantum coherence. We investigate three model systems: the photodissociation of three coupled Morse oscillators, the bound state dynamics of two coupled Morse oscillators, and a two-dimensional model for collinear triatomic vibrational dynamics. In all cases, the pseudospectral Gaussian method is in quantitative agreement with numerically exact calculations. The results are promising for nonadiabatic molecular dynamics in molecular systems where strongly correlated ground or excited states require expensive electronic structure calculations. [ABSTRACT FROM AUTHOR]

This Special Topic section on Advanced Electronic Structure Methods for Solids and Surfaces contains a collection of research papers that showcase recent advances in the high accuracy prediction of materials and surface properties. It provides a timely snapshot of a growing field that is of broad importance to chemistry, physics, and materials science. [ABSTRACT FROM AUTHOR]

Ab initio electronic structure calculations of two-dimensional layered structures are typically performed using codes that were developed for three-dimensional structures, which are periodic in all three directions. The introduction of a periodicity in the third direction (perpendicular to the layer) is completely artificial and may lead in some cases to spurious results and to difficulties in treating the action of external fields. In this paper we develop a new approach, which is "native" to quasi-2D materials, making use of basis function that are periodic in the plane, but atomic-like in the perpendicular direction. We show how some of the basic tools of ab initio electronic structure theory -- density functional theory, GW approximation and Bethe-Salpeter equation -- are implemented in the new basis. We argue that the new approach will be preferable to the conventional one in treating the peculiarities of layered materials, including the long range of the unscreened Coulomb interaction in insulators, and the effects of strain, corrugations, and external fields. [ABSTRACT FROM AUTHOR]

Electronic structure calculations of atomistic systems based on density functional theory involve solving the Poisson equation. In this paper, we present a plane-wave based algorithm for solving the generalized Poisson equation subject to periodic or homogeneous Neumann conditions on the boundaries of the simulation cell and Dirichlet type conditions imposed at arbitrary subdomains. In this way, source, drain, and gate voltages can be imposed across atomistic models of electronic devices. Dirichlet conditions are enforced as constraints in a variational framework giving rise to a saddle point problem. The resulting system of equations is then solved using a stationary iterative method in which the generalized Poisson operator is preconditioned with the standard Laplace operator. The solver can make use of any sufficiently smooth function modelling the dielectric constant, including density dependent dielectric continuum models. For all the boundary conditions, consistent derivatives are available and molecular dynamics simulations can be performed. The convergence behaviour of the scheme is investigated and its capabilities are demonstrated. [ABSTRACT FROM AUTHOR]

In this paper, the size-dependent changes in energetic, vibrational, and electronic properties of C-O gas molecule interacting with surface Pd atom of a variety of AuPd nanoalloy structures are investigated by means of first principles calculations. The variation in C-O adsorption energies, C-O vibration frequencies (?C-O), and Pd d-bond centers (ed) on a series of non-supported Aun-1-Pd1 nanoparticles (with n varying from 13 to 147) and on two semi-finite surfaces are inspected with cluster size. We demonstrate for the first time that, with small AuPd bimetallic three-dimensional clusters as TOh38, one can reach cluster size convergence even for such a sensitive observable as the adsorption energy on a metal surface. Indeed, the results show that the adsorbate-induced perturbation is extremely local and it only concerns the isolated Pd interacting with the reactive gas molecule. Except for 13 atom clusters, in which molecular behaviour is predominant, no finite-size effects are observed for surface Pd atom substituted in AuPd free nanoclusters above 38 atoms. [ABSTRACT FROM AUTHOR]

Functionalized materials consisting of inorganic substrates with organic adsorbates play an increasing role in emerging technologies like molecular electronics or hybrid photovoltaics. For such applications, the adsorption geometry of the molecules under operating conditions, e.g., ambient temperature, is crucial because it influences the electronic properties of the interface, which in turn determine the device performance. So far detailed experimental characterization of adsorbates at room temperature has mainly been done using a combination of complementary methods like photoelectron spectroscopy together with scanning tunneling microscopy. However, this approach is limited to ensembles of adsorbates. In this paper, we show that the characterization of individual molecules at room temperature, comprising the determination of the adsorption configuration and the electrostatic interaction with the surface, can be achieved experimentally by atomic force microscopy (AFM) and Kelvin probe force microscopy (KPFM).We demonstrate this by identifying two different adsorption configurations of isolated copper(II) meso-tetra (4-carboxyphenyl) porphyrin (Cu-TCPP) on rutile TiO2 (110) in ultra-high vacuum. The local contact potential difference measured by KPFM indicates an interfacial dipole due to electron transfer from the Cu-TCPP to the TiO2. The experimental results are verified by state-of-the-art first principles calculations. We note that the improvement of the AFM resolution, achieved in this work, is crucial for such accurate calculations. Therefore, high resolution AFM at room temperature is promising for significantly promoting the understanding of molecular adsorption. [ABSTRACT FROM AUTHOR]

The charge transfer integral is a key parameter required by various theoretical models to describe charge transport properties, e.g., in organic semiconductors. The accuracy of this important property depends on several factors, which include the level of electronic structure theory and internal simplifications of the applied formalism. The goal of this paper is to identify the performance of various approximate approaches of the latter category, while using the high level equation-of-motion coupled cluster theory for the electronic structure. The calculations have been performed on the ethylene dimer as one of the simplest model systems. By studying different spatial perturbations, it was shown that while both energy split in dimer and fragment charge difference methods are equivalent with the exact formulation for symmetrical displacements, they are less efficient when describing transfer integral along the asymmetric alteration coordinate. Since the "exact" scheme was found computationally expensive, we examine the possibility to obtain the asymmetric fluctuation of the transfer integral by a Taylor expansion along the coordinate space. By exploring the efficiency of this novel approach, we show that the Taylor expansion scheme represents an attractive alternative to the "exact" calculations due to a substantial reduction of computational costs, when a considerably large region of the potential energy surface is of interest. Moreover, we show that the Taylor expansion scheme, irrespective of the dimer symmetry, is very accurate for the entire range of geometry fluctuations that cover the space the molecule accesses at room temperature. [ABSTRACT FROM AUTHOR]

Tensor factorization of the 2-electron integral matrix is a well-known technique for reducing the computational scaling of ab initio electronic structure methods toward that of Hartree-Fock and density functional theories. The simplest factorization that maintains the positive semidefinite character of the 2-electron integral matrix is the Cholesky factorization. In this paper, we introduce a family of positive semidefinite factorizations that generalize the Cholesky factorization. Using an implementation of the factorization within the parametric 2-RDM method [D. A. Mazziotti, Phys. Rev. Lett. 101, 253002 (2008)], we study several inorganic molecules, alkane chains, and potential energy curves and find that this generalized factorization retains the accuracy and size extensivity of the Cholesky factorization, even in the presence of multi-reference correlation. The generalized family of positive semidefinite factorizations has potential applications to low-scaling ab initio electronic structure methods that treat electron correlation with a computational cost approaching that of the Hartree-Fock method or density functional theory. [ABSTRACT FROM AUTHOR]

The aim of this paper is to shed light on the topology and properties of the nodes (i.e., the zeros of the wave function) in electronic systems. Using the "electrons on a sphere" model, we study the nodes of two-, three-, and four-electron systems in various ferromagnetic configurations (sp, p², sd, pd, p³, sp², and sp³). In some particular cases (sp, p², sd, pd, and p³), we rigorously prove that the non-interacting wave function has the same nodes as the exact (yet unknown) wave function. The number of atomic and molecular systems for which the exact nodes are known analytically is very limited and we show here that this peculiar feature can be attributed to interdimensional degeneracies. Although we have not been able to prove it rigorously, we conjecture that the nodes of the non-interacting wave function for the sp3 configuration are exact. [ABSTRACT FROM AUTHOR]

Changes in the electronic structure of pigments in protein environments and of polar molecules in solution inevitably induce a re-adaption of molecular nuclear structure. Both changes of electronic and vibrational energies can be probed with visible or infrared lasers, such as two-dimensional electronic spectroscopy or vibrational spectroscopy. The extent to which the two changes are correlated remains elusive. The recent demonstration of two-dimensional electronic-vibrational (2DEV) spectroscopy potentially enables a direct measurement of this correlation experimentally. However, it has hitherto been unclear how to characterize the correlation from the spectra. In this paper, we present a theoretical formalism to demonstrate the slope of the nodal line between the excited state absorption and ground state bleach peaks in the spectra as a characterization of the correlation between electronic and vibrational transition energies. We also show the dynamics of the nodal line slope is correlated to the vibrational spectral dynamics. Additionally, we demonstrate the fundamental 2DEV spectral line-shape of a monomer with newly developed response functions. [ABSTRACT FROM AUTHOR]

The electronic structure of eight zinc-centered porphyrin macrocyclic molecules are investigated using density functional theory for ground-state properties, time-dependent density functional theory (TDDFT) for excited states, and Franck-Condon (FC) analysis for further characterization of the UV-vis spectrum. Symmetry breaking was utilized to find the lowest energy of the excited states for many states in the spectra. To confirm the theoretical modeling, the spectroscopic result from zinc phthalocyanine (ZnPc) is used to compare to the TDDFT and FC result. After confirmation of the modeling, five more planar molecules are investigated: zinc tetrabenzoporphyrin (ZnTBP), zinc tetrabenzomonoazaporphyrin (ZnTBMAP), zinc tetrabenzocisdiazaporphyrin (ZnTBcisDAP), zinc tetrabenzotransdiazaporphyrin (ZnTBtransDAP), and zinc tetrabenzotriazaporphyrin (ZnTBTrAP). The two latter molecules are then compared to their phenylated sister molecules: zinc monophenyltetrabenzotriazaporphyrin (ZnMPTBTrAP) and zinc diphenyltetrabenzotransdiazaporphyrin (ZnDPTBtransDAP). The spectroscopic results from the synthesis of ZnMPTBTrAP and ZnDPTBtransDAP are then compared to their theoretical models and non-phenylated pairs. While the Franck-Condon results were not as illuminating for every B-band, the Q-band results were successful in all eight molecules, with a considerable amount of spectral analysis in the range of interest between 300 and 750 nm. The π-π* transitions are evident in the results for all of the Q bands, while satellite vibrations are also visible in the spectra. In particular, this investigation finds that, while ZnPc has a D4h symmetry at ground state, a C4v symmetry is predicted in the excited-state Q band region. The theoretical results for ZnPc found an excitation energy at the Q-band 0-0 transition of 1.88 eV in vacuum, which is in remarkable agreement with published gas-phase spectroscopy, as well as our own results of ZnPc in solution with Tetrahydrofuran that are provided in this paper. [ABSTRACT FROM AUTHOR]

In this paper, we provide detailed insights into the electronic structure of the gas phase biphenylene molecule through core and valence spectroscopy. By comparing results of X-ray Photoelectron Spectroscopy (XPS) measurements with ASCF core-hole calculations in the framework of Density Functional Theory (DFT), we could decompose the characteristic contributions to the total spectra and assign them to non-equivalent carbon atoms. As a difference with similar molecules like biphenyl and naphthalene, an influence of the localized orbitals on the relative XPS shifts was found. The valence spectrum probed by photoelectron spectroscopy at a photon energy of 50 eV in conjunction with hybrid DFT calculations revealed the effects of the localization on the electronic states. Using the transition potential approach to simulate the X-ray absorption spectroscopy measurements, similar contributions from the non-equivalent carbon atoms were determined from the total spectrum, for which the slightly shifted individual components can explain the observed asymmetric features. [ABSTRACT FROM AUTHOR]

The use of broken-symmetry electronic structure methods is required in order to obtain correct behavior of electronically strained open-shell systems, such as transition states, biradicals, and transition metals. This approach often has issues with spin contamination, which can lead to significant errors in predicted energies, geometries, and properties. Approximate projection schemes are able to correct for spin contamination and can often yield improved results. To fully make use of these methods and to carry out exploration of the potential energy surface, it is desirable to develop an efficient second energy derivative theory. In this paper, we formulate the analytical second derivatives for the Yamaguchi approximate projection scheme, building on recent work that has yielded an efficient implementation of the analytical first derivatives. [ABSTRACT FROM AUTHOR]

In this paper, we systematically investigate the electronic structure for the ²Σ+ ground state of the polar alkali-metal-alkaline-earth-metal molecules BaAlk (Alk = Li, Na, K, Rb, and Cs). Potential energy curves and permanent dipole moments (PDMs) are determined using power quantum chemistry complete active space self-consistent field and multi-reference configuration interaction methods. Basic spectroscopic constants are derived from ro-vibrational bound state calculation. From the calculations, it is shown that BaK, BaRb, and BaCs molecules have moderate values of PDM at the equilibrium bond distance (BaK:1.62 D, BaRb:3.32 D, and BaCs:4.02 D). Besides, the equilibrium bond length (4.93 Å and 5.19 Å) and dissociation energy (0.1825 eV and 0.1817 eV) for the BaRb and BaCs are also obtained. [ABSTRACT FROM AUTHOR]

In Paper I of this series [F. A. Bischoff, "Regularizing the molecular potential in electronic structure calculations. I. SCF methods," J. Chem. Phys. 141,184105 (2014)] a regularized molecular Hamilton operator for electronic structure calculations was derived and its properties in SCF calculations were studied. The regularization was achieved using a correlation factor that models the electron-nuclear cusp. In the present study we extend the regularization to correlated methods, in particular the exact solution of the two-electron problem, as well as second-order many body perturbation theory. The nuclear and electronic correlation factors lead to computations with a smaller memory footprint because the singularities are removed from the working equations, which allows coarser grid resolution while maintaining the precision. Numerical examples are given. [ABSTRACT FROM AUTHOR]

Methanol synthesis from CO and H2 over ZnO, which requires high temperatures and high pressures giving rise to a complex interplay of physical and chemical processes over this heterogeneous catalyst surface, is investigated using ab initio simulations. The redox properties of the surrounding gas phase are known to directly impact on the catalyst properties and thus, set the overall catalytic reactivity of this easily reducible oxide material. In Paper III of our series [J. Kiss, J. Frenzel, N. N. Nair, B. Meyer, and D. Marx, J. Chem. Phys. 134, 064710 (2011)] we have qualitatively shown that for the partially hydroxylated and defective ZnO(0001) surface there exists an intricate network of surface chemical reactions. In the present study, we employ advanced molecular dynamics techniques to resolve in detail this reaction network in terms of elementary steps on the defective surface, which is in stepwise equilibrium with the gas phase. The two individual reduction steps were investigated by ab initio metadynamics sampling of free energy landscapes in three-dimensional reaction subspaces. By also sampling adsorption and desorption processes and thus molecular species that are in the gas phase but close to the surface, our approach successfully generated several alternative pathways of methanol synthesis. The obtained results suggest an Eley-Rideal mechanism for both reduction steps, thus involving "near-surface" molecules from the gas phase, to give methanol preferentially over a strongly reduced catalyst surface, while important side reactions are of Langmuir-Hinshelwood type. Catalyst re-reduction by H2 stemming from the gas phase is a crucial process after each reduction step in order to maintain the catalyst's activity toward methanol formation and to close the catalytic cycle in some reaction channels. Furthermore, the role of oxygen vacancies, side reactions, and spectator species is investigated and mechanistic details are discussed based on extensive electronic structure analysis. [ABSTRACT FROM AUTHOR]

The photodesorption of CO from TiO2(110) by femtosecond pulses is investigated with the Surrogate Hamiltonian approach. The aim of the study is to resolve the relaxation mechanism and forecast the lifetime of the exited state based on a microscopic description of the excitation and relaxation processes. The parameters characterizing the system are obtained from ab initio and Density Functional Theory-calculations with one parameter estimated from physical considerations and convergence studies. Two electronic states are considered and the relaxation is assumed to be due to the interaction of the excited adsorbate with electron hole pairs in the surface. Desorption probabilities and velocity distributions of the desorbing molecules are calculated and an exited state lifetime is predicted. Throughout this paper atomic units, i.e., ℏ = me = e = a0 = 1, have been used unless otherwise stated. [ABSTRACT FROM AUTHOR]

Relativistic, infinite order exact two-component, density functional theory electronic structure calculations were performed for MCl4 and MOCl2 of group-4 elements Ti, Zr, Hf, and element 104, Rf, with the aim to predict their behaviour in gas-phase chromatography experiments. RfCl4 and RfOCl2 were shown to be less stable than their lighter homologs in the group, tetrachlorides and oxychlorides of Zr and Hf, respectively. The oxychlorides turned out to be stable as a bent structure, though the stabilization energy with respect to the flat one (C2v) is very small. The trend in the formation of the tetrachlorides from the oxychlorides in group 4 is shown to be Zr < Hf < Rf, while the one in the formation of the oxychlorides from the chlorides is opposite. All the calculated properties are used to estimate adsorption energy of these species on various surfaces in order to interpret results of gas-phase chromatography experiments, as is shown in Paper II. [ABSTRACT FROM AUTHOR]