Skip to main content

    Peter Reinhardt

    We try to decompose the interaction energy of the hydrogen bonded dimer systems HOH⋯NH3, H2O⋯HNH2, and H2O⋯HNH3+ into several well-defined quantities, using the fragment-based density-functional calculations, Symmetry-adapted perturbation... more
    We try to decompose the interaction energy of the hydrogen bonded dimer systems HOH⋯NH3, H2O⋯HNH2, and H2O⋯HNH3+ into several well-defined quantities, using the fragment-based density-functional calculations, Symmetry-adapted perturbation theory, and as well the Heitler–London approach to interaction energies. The theoretical considerations and the numerical results show that despite the different nature of the decomposition schemes common terms beyond the electrostatic
    ABSTRACT The periodic Hartree-Fock approach is used in all-electron calculations on bulk rutile (TiO2) and various slabs modeling the (001), (100), and (110) surface to study possible geometric and electronic surface relaxations,... more
    ABSTRACT The periodic Hartree-Fock approach is used in all-electron calculations on bulk rutile (TiO2) and various slabs modeling the (001), (100), and (110) surface to study possible geometric and electronic surface relaxations, convergence of properties with slab thickness, and the relative stability of the different surfaces. Results are discussed by means of a Mulliken population analysis, band structures, and the density of states. We find displacements of surface atoms ranging from 5 to nearly 40 pm; the (001) surface appears to be significantly unstable with respect to the other two [(100) and (110)], in accord with experimental findings. Specific features of the density of states projected onto different Ti sites and the relative stability of the different surfaces agree with a possible dissociation behavior of the rutile structure into TiO2 chains, which are calculated to be stable.
    Research Interests:
    ABSTRACT The dressed diagonal approximation to the self-consistent – size consistent CI, corrected for off-diagonal Fock matrix elements in localized orbitals is developed and applied to the ammoniac dimer system. A quite correct... more
    ABSTRACT The dressed diagonal approximation to the self-consistent – size consistent CI, corrected for off-diagonal Fock matrix elements in localized orbitals is developed and applied to the ammoniac dimer system. A quite correct correlation energy can be obtained for this system, with a significantly reduced dependence of the results on the choice of the localization procedure. When calculating an interaction energy, the choice of the monomer orbitals and the application of the Boys–Bernardi counterpoise procedure shows in this case an unusual behavior: the correlation energy does not increase with the size of the atomic basis sets. Nevertheless a reasonable potential curve can be obtained.
    ABSTRACT The geometrical structure of the polar, Al-terminated surface of corundum has been investigated by different theoretical methods ranging from empirical pair potentials to periodic Hartree–Fock or density functional theory... more
    ABSTRACT The geometrical structure of the polar, Al-terminated surface of corundum has been investigated by different theoretical methods ranging from empirical pair potentials to periodic Hartree–Fock or density functional theory calculations. The different methods agree in predicting a structure for the bulk, which is close to the experimental one, and an inward relaxation of more than 50% in agreement with experiment. However, a rather strong disagreement exists concerning the quantitative amount of relaxation of the surface layers. Accurate density functional theory calculations for slab models considering up to 18 layers with 30 atoms in the unit cell strongly suggest that the degree of relaxation of the outermost Al atomic layer is much larger than the values obtained by other methods. Reasons for the different description are given and implications for the use of pair potentials discussed.
    Human adult germline stem cells (haGSCs) were established from human testicular biopsies and were claimed to be pluripotent. Recently, the gene expression profile of haGSCs demonstrated that these cells presented with a fibroblast rather... more
    Human adult germline stem cells (haGSCs) were established from human testicular biopsies and were claimed to be pluripotent. Recently, the gene expression profile of haGSCs demonstrated that these cells presented with a fibroblast rather than a pluripotent identity. Nevertheless, haGSCs were reported to generate teratomas. In this report, we address this discrepancy. Instead of using haGSCs, which are no longer available for the stem cell community, we used a human testicular fibroblastic cell (hTFC) line that presents with a gene expression profile highly similar to that of haGSCs. Indeed, as shown by microarray analysis, the similarity between hTFCs and haGSCs is comparable to human embryonic stem cell (hESC) lines derived by different laboratories. We argue that the almost identical gene expression profile of hTFCs and haGSCs should result in a very similar if not identical differentiation potential. Strikingly, hTFCs were not able to generate teratomas after injection into nude ...
    A method is proposed for the solution of the self-consistent field equations that can lead to localized occupied and virtual molecular orbitals, avoiding the need for solving for the canonical molecular orbitals. The method starts with... more
    A method is proposed for the solution of the self-consistent field equations that can lead to localized occupied and virtual molecular orbitals, avoiding the need for solving for the canonical molecular orbitals. The method starts with strongly localized “guess molecular orbitals”, it is nonperturbative and proceeds through the diagonalization of single configuration interaction matrices which may be rendered size-consistent through
    ABSTRACT Starting with self-consistent fields (SCFs), localized orbitals should facilitate the calculation of the correlation energy in extended, and in particular periodic, systems. This idea is exploited on model ring systems (H4n+2).... more
    ABSTRACT Starting with self-consistent fields (SCFs), localized orbitals should facilitate the calculation of the correlation energy in extended, and in particular periodic, systems. This idea is exploited on model ring systems (H4n+2). It is shown that for insulating materials [(H2)2n+1, presenting a large gap in the band structure], most of the energy lowering brought by the orders larger than 2 in the canonical many-body perturbation expansion are due to the local-hole–local-particle interaction and that the localized Epstein–Nesbet second-order energies are close to the best correlation-energy estimates. The situation is completely different for small-gap (metalliclike) systems, such as cyclic H4n+2, where the localized second-order approach misses a large fraction of the correlation energy, involving the propagation of the holes and of the particles and implying specific higher-order diagrams. © 1998 American Institute of Physics.
    ABSTRACT In this work we test two ab initio methodologies which allow the decomposition of the total intermolecular interaction energy into physically meaningful contributions, namely, the symmetry adapted perturbation theory (SAPT) and... more
    ABSTRACT In this work we test two ab initio methodologies which allow the decomposition of the total intermolecular interaction energy into physically meaningful contributions, namely, the symmetry adapted perturbation theory (SAPT) and the use of localized orbitals within a Møller–Plesset perturbation scheme. The accuracy of the two different methods is compared to supermolecular results, within MP2 and coupled-cluster theory within single and double excitations, with perturbative estimates of the amplitudes of triple excitations [CCSD(T)]. Some relations between the different approaches are conjectured from theoretical considerations, and are confirmed by numerical results. The corresponding calculations have been performed for three model dimers: two NH3⋯H2O dimers, with NH3 acting once as a proton acceptor and once as a proton donor, and the NH4+⋯H2O considered as a prototype of the ion–molecule interaction. We may conclude that third-order terms in SAPT help significantly to reproduce the Hartree–Fock induction and the relaxed, total dispersion in the LMP2 decomposition. © 2003 American Institute of Physics.
    Starting from localized bond or lone-pair Hartree-Fock molecular orbitals, one may define contracted doubly excited functions for each pair of bond molecular orbitals. These functions are obtained from local single- and... more
    Starting from localized bond or lone-pair Hartree-Fock molecular orbitals, one may define contracted doubly excited functions for each pair of bond molecular orbitals. These functions are obtained from local single- and double-configuration interaction (CISD) of moderate size. Then one may build a contracted CISD matrix for the whole molecule, spanned by the Hartree-Fock determinant and these contracted doubly excited functions, the number of which is indeed moderate, as scaling at most as the square of the number of bonds. The calculation of the off-diagonal elements of this matrix is straightforward. Its diagonalization provides an upper bound to the lowest CISD eigenvalue. The well-known size-consistency error may be overcome through self-consistent dressings such as coupled-electron pair approximations, and cutoff criteria will lead to linear scaling. Numerical tests on a series of covalent and ionic systems show that the results are very close to that of coupled-cluster calculations. Possible improvements of this already efficient algorithm are suggested.
    ABSTRACT Localized molecular Hartree–Fock orbitals have been determined by means of an iterative procedure consisting of orthogonalization and configuration interaction employing single excitations. For ring systems the rotational... more
    ABSTRACT Localized molecular Hartree–Fock orbitals have been determined by means of an iterative procedure consisting of orthogonalization and configuration interaction employing single excitations. For ring systems the rotational symmetry has been included explicitly to obtain Wannier-like orbitals suited for a posteriori correlation calculations using only the most important contributions within a limited region around one reference ring segment. Applications involving different estimates of the correlation energy include as model systems (H2)2n+1, the ionic LinHn, and a weakly bound beryllium ring as well as the strongly covalent molecule (CH2)36 forming a closed ring. In all cases, the localized and canonical MP2 results are close, and the localized Epstein–Nesbet second-order gives a good estimate of more expensive MP4 or CEPA-0 values. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 70: 167–180, 1998
    This article recalls basic properties of Slater type Orbital (STO) basis sets, in particular, with respect to satisfying Kato's electron-nucleus cusp condition. It is shown that a suitable starting point is the set of... more
    This article recalls basic properties of Slater type Orbital (STO) basis sets, in particular, with respect to satisfying Kato's electron-nucleus cusp condition. It is shown that a suitable starting point is the set of hydrogen-like orbitals or, more flexibly, their generalization as Coulomb Sturmians. A case study is presented on the Hartree-Fock density obtained for a water molecule, N2, and
    ABSTRACT In the present paper it is proposed—without a formal physical demonstration—that the probability to select the lifetime of an activated complex possibly yielding ozone and formed by reactions involving distinguishable isotopes is... more
    ABSTRACT In the present paper it is proposed—without a formal physical demonstration—that the probability to select the lifetime of an activated complex possibly yielding ozone and formed by reactions involving distinguishable isotopes is not equal to the probability to select the lifetime of the same complex but formed by reaction involving indistinguishable isotopes.The origin of this effect is supposed to be linked to the quantum mechanical principle according to which it is not permitted, in a scattering process, to separate the probabilities describing the recoil from the incident particle, if the two particles are indistinguishable. We ascribe the mass independent isotopic fractionation factor to the lifetime ratio of the complexes formed by reactions involving distinguishable and indistinguishable isotopes.In order to illustrate the consequences of this principle, all the accompanying scattering calculations of O+O2→O3*→O+O2 were performed in a thermal gas with oxygen isotopes having the same mass (16 amu). It is numerically shown that the origin of the experimentally observed mass-independent isotope fractionation in ozone is well accounted for within this single assumption.The effect is a priori not restricted to ozone or to oxygen isotopes and may have played a role in several other isotopic anomalies found in meteorites.
    All electron, first principles, periodic calculations have been carried out to investigate the relaxation of the Al-terminated corundum surface. Different theoretical formalisms, Hartree–Fock (HF) and density functional theory using... more
    All electron, first principles, periodic calculations have been carried out to investigate the relaxation of the Al-terminated corundum surface. Different theoretical formalisms, Hartree–Fock (HF) and density functional theory using either the LDA or the hybrid B3LYP exchange-correlation functional, have been employed. In the past, many ab initio studies of surface structure have been performed using the experimental, rather than the
    The hydration free energy, structure, and dynamics of the zinc divalent cation are studied using a polarizable force field in molecular dynamics simulations. Parameters for the Zn2+ are derived from gas-phase ab initio calculation of the... more
    The hydration free energy, structure, and dynamics of the zinc divalent cation are studied using a polarizable force field in molecular dynamics simulations. Parameters for the Zn2+ are derived from gas-phase ab initio calculation of the Zn2+-water dimer. The Thole-based dipole polarization is adjusted on the basis of the constrained space orbital variations (CSOV) calculation, while the symmetry adapted perturbation theory (SAPT) approach is also discussed. The vdW parameters of Zn2+ have been obtained by comparing the AMOEBA Zn2+ water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn2+ solvation in bulk water are subsequently performed with the polarizable force field. The calculated first-shell water coordination number, water residence time, and free energy of hydration are consistent with experimental and previous theoretical values. The study is supplemented with extensive reduced variati...
    Research Interests:
    As for generating localized Hartree−Fock orbitals, we propose a potentially linear-scaling singles-CI scheme to construct fragment-localized density functional theory (DFT) orbitals for molecular systems as water clusters. Due to the use... more
    As for generating localized Hartree−Fock orbitals, we propose a potentially linear-scaling singles-CI scheme to construct fragment-localized density functional theory (DFT) orbitals for molecular systems as water clusters. Due to the use of a deformation step instead of a localization step, the influence of the environment on each separate molecule can be studied in detail. The generated orbital set for the whole molecular system is strictly equivalent to a set of canonical orbitals and is a subsequent energy decomposition of intermolecular interactions into electrostatic, exchange repulsion, and orbital interaction, well beyond dimer systems. Beyond this, the correspondence of the individual orbitals to the initial monomer orbitals permits to assess how an interaction deforms an electron density. We show this for dipole moments, which may be decomposed into monomer contributions, polarization, and charge-transfer contribution. Applications to a water and an ammonia dimer and chains...
    Research Interests:
    In this contribution, we propose accurate intermolecular benchmarks for the 10 standardized water dimers initially proposed by van Duijneveldt et al., following a study by Smith et al. (J Chem Phys, 1990, 92, 1240). Using the popular... more
    In this contribution, we propose accurate intermolecular benchmarks for the 10 standardized water dimers initially proposed by van Duijneveldt et al., following a study by Smith et al. (J Chem Phys, 1990, 92, 1240). Using the popular triple-zeta aug-cc-pVTZ basis set, symmetry-adapted perturbation theory (SAPT) computations have been performed and compared with supermolecular post-Hartree-Fock techniques up to most recent explicitly correlated methods. Effects of approximated inclusion of third-order corrections to induction and exchange-induction energies are discussed. As SAPT results are close to the best available ab initio ones, detailed analysis of the available individual SAPT contributions to the interaction energy confirms that electron correlation acts not only on its long-range dispersion part but also on its separated physical components. This permits to assess the validity of the usual HF+Dispersion model as approximation to add pure correlation effects to orbital-based...
    Research Interests:
    Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is... more
    Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is computed in the spirit of the sum of interacting fragment ab initio (SIBFA) force field by separately evaluating each one of its components: electrostatic (Coulomb), exchange repulsion, polarization, and charge transfer intermolecular interaction energies, in order to reproduce reference constrained space orbital variation (CSOV) energy decomposition calculations at the B3LYP/aug-cc-pVTZ level. The use of an auxiliary basis set restricted to spherical Gaussian functions facilitates the rotation of the fitted densities of rigid fragments and enables a fast and accurate density fitting evaluation of Coulomb and exchange-repulsion energy, the latter using the overlap model introduced by Wheatley and Price [Mol. Phys. 69, 50718 (1990)]. The SIBFA energy scheme for polarization and charge transfer has been implemented using the electric fields and electrostatic potentials generated by the fitted densities. GEM-0 has been tested on ten stationary points of the water dimer potential energy surface and on three water clusters (n=16,20,64). The results show very good agreement with density functional theory calculations, reproducing the individual CSOV energy contributions for a given interaction as well as the B3LYP total interaction energies with errors below kBT at room temperature. Preliminary results for Coulomb and exchange-repulsion energies of metal cation complexes and coupled cluster singles doubles electron densities are discussed
    For the first time the coexistence of a sigma- and a pi-complex in the C(4)H(4)O:HCl system has been observed, in the same supersonic expansion of a molecular jet seeded with argon (or helium) or in a flow-cooled cell at 240 K. This is an... more
    For the first time the coexistence of a sigma- and a pi-complex in the C(4)H(4)O:HCl system has been observed, in the same supersonic expansion of a molecular jet seeded with argon (or helium) or in a flow-cooled cell at 240 K. This is an exception to the third of the Legon-Miller rules which claims the sigma-structure to be the only one to exist. On the grounds of energetic considerations and band contour simulations, two observed bands at 2787.7 and 2795.5 cm(-1) of the nu(s) HCl stretching frequency are assigned to the two complexes, recorded as Fourier transform infrared spectra with a resolution between 0.2 and 0.5 cm(-1). Complementary calculations show that the use of the standard second-order Moller-Plesset perturbation theory may be erroneous for such a complex, due of the overestimation of the dispersion contribution with respect to the electrostatic term. It is finally established that only a balanced version of the second-order Moller-Plesset perturbation method, spin-component scaled-MP2, or a higher level of theory like a coupled-cluster approach, can provide a reliable energetic analysis for this complex.
    ABSTRACT We present a Monte Carlo scheme to simulate equilibrium structures and thermal desorption (TD) spectra of monolayer one- or two-component metal adsorbate films supported on a close-packed metal substrate surface. The scheme is... more
    ABSTRACT We present a Monte Carlo scheme to simulate equilibrium structures and thermal desorption (TD) spectra of monolayer one- or two-component metal adsorbate films supported on a close-packed metal substrate surface. The scheme is used to determine the miscibility properties of an example system, Ag-Cu/Ru(0001), by tuning the model parameters such that experimental TD spectra of this system are well reproduced. The corresponding set of parameters is shown to be unique. We discuss the role of lateral interactions and of algorithmic details. Finally, our results are compared with recently published data for the surface of a Ag-Cu bulk alloy very dilute in Ag [1].
    ABSTRACT We present electronic and structural parameters for bulk rutile (TiO2) determined by means of different theoretical methods, namely, the periodic Hartree-Fock approach, subsequent post-Hartree-Fock density functional correlation... more
    ABSTRACT We present electronic and structural parameters for bulk rutile (TiO2) determined by means of different theoretical methods, namely, the periodic Hartree-Fock approach, subsequent post-Hartree-Fock density functional correlation correction, and a linear-combination-of-atomic-orbitals approach based entirely on density functionals. Differences between the use of effective core potentials, a frozen-core treatment, and all-electron calculations are investigated. © 1996 John Wiley & Sons, Inc.
    Selecting excitations in localized orbitals to calculate long-range correlation contributions to range-separated density-functional theory can reduce the overall computational effort significantly. Beyond simple selection schemes of... more
    Selecting excitations in localized orbitals to calculate long-range correlation contributions to range-separated density-functional theory can reduce the overall computational effort significantly. Beyond simple selection schemes of excited determinants, the dispersion-only approximation, which avoids counterpoise-corrected monomer caculations, is shown to be particularly interesting in this context, which we apply to the random-phase approximation. The approach has been tested on dimers of formamide, water, methane and benzene.
    We consider several spin-unrestricted random-phase approximation (RPA) variants for calculating correlation energies, with and without range separation, and test them on datasets of atomization energies and reaction barrier heights. We... more
    We consider several spin-unrestricted random-phase approximation (RPA) variants for calculating correlation energies, with and without range separation, and test them on datasets of atomization energies and reaction barrier heights. We show that range separation greatly improves the accuracy of all RPA variants for these properties. Moreover, we show that a RPA variant with exchange, hereafter referred to as RPAx-SO2, first proposed by Szabo and Ostlund [J. Chem. Phys. 67, 4351 (1977)] in a spin-restricted closed-shell formalism, and extended here to a spin-unrestricted formalism, provides on average the most accurate range-separated RPA variant for atomization energies and reaction barrier heights. Since this range-separated RPAx-SO2 method had already been shown to be among the most accurate range-separated RPA variants for weak intermolecular interactions [J. Toulouse et al., J. Chem. Phys. 135, 084119 (2011)], this works confirms range-separated RPAx-SO2 as a promising method for general chemical applications.