Coexistence of 1,3-butadiene conformers in ionization energies and Dyson orbitals

The minimum-energy structures on the torsional potential-energy surface of 1,3-butadiene have been studied quantum mechanically using a range of models including ab initio Hartree-Fock and second-order M?ller-Plesset theories, outer valence Green's function, and density-functional theory with a hybrid functional and statistical average orbital potential model in order to understand the binding-energy (ionization energy) spectra and orbital cross sections observed by experiments. The unique full geometry optimization process locates the s-trans-1,3-butadiene as the global minimum structure and the s-gauche-1,3-butadiene as the local minimum structure. The latter possesses the dihedral angle of the central carbon bond of 32.81 degrees in agreement with the range of 30 degrees-41 degrees obtained by other theoretical models. Ionization energies in the outer valence space of the conformer pair have been obtained using Hartree-Fock, outer valence Green's function, and density-functional (statistical average orbital potentials) models, respectively. The Hartree-Fock results indicate that electron correlation (and orbital relaxation) effects become more significant towards the inner shell. The spectroscopic pole strengths calculated in the Green's function model are in the range of 0.85-0.91, suggesting that the independent particle picture is a good approximation in the present study. The binding energies from the density-functional (statisticaly averaged orbital potential) model are in good agreement with photoelectron spectroscopy, and the simulated Dyson orbitals in momentum space approximated by the density-functional orbitals using plane-wave impulse approximation agree well with those from experimental electron momentum spectroscopy. The coexistence of the conformer pair under the experimental conditions is supported by the approximated experimental binding-energy spectra due to the split conformer orbital energies, as well as the orbital momentum distributions of the mixed conformer pair observed in the orbital cross sections of electron momentum spectroscopy.     

The optimal P3M algorithm for computing electrostatic energies in periodic systems

We optimize Hockney and Eastwood's particle-particle particle-mesh algorithm to achieve maximal accuracy in the electrostatic energies (instead of forces) in three-dimensional periodic charged systems. To this end we construct an optimal influence function that minimizes the root-mean-square (rms) errors of the energies. As a by-product we derive a new real-space cutoff correction term, give a transparent derivation of the systematic errors in terms of Madelung energies, and provide an accurate analytical estimate for the rms error of the energies. This error estimate is a useful indicator of the accuracy of the computed energies and allows an easy and precise determination of the optimal values of the various parameters in the algorithm (Ewald splitting parameter, mesh size, and charge assignment order).     

A calculation of exciton energies in periodic systems with helical symmetry: Application to a hydrogen fluoride chain

Ab initio Hartree–Fock crystal orbital results were used as input for the calculation of exciton energies in the Tamm–Dancoff and random-phase approximation for polymers with helical symmetry. The calculations were applied to a hydrogen fluoride chain. We show that the Tamm–Dancoff method is a good approximation to the random-phase theory. Furthermore, the influence of the band–band interaction in the exciton calculation has been investigated.     

Perturbative treatment of the electron-correlation contribution to the diagonal Born-Oppenheimer correction

A perturbative scheme for the treatment of electron-correlation effects on the diagonal Born-Oppenheimer correction (DBOC) is suggested. Utilizing the usual Moller-Plesset partitioning of the Hamiltonian formulas for first and second orders (termed as MP1 and MP2) are obtained by expanding the wave function in the corresponding coupled-cluster expressions for the DBOC[J. Gauss et al., J. Chem. Phys. 125, 144111 (2006)]. The obtained expressions are recast in terms of one- and two-particle density matrices in order to take advantage of existing analytic second-derivative implementations for many-body methods. Test calculations show that both MP1 and MP2 recover large fractions (on average 90% and 95%, respectively) of the coupled-cluster singles and doubles (CCSD) electron-correlation corrections to the DBOC and thus render the suggested MP treatments cost-effective (though still accurate) alternatives to high-level coupled cluster (CC) treatments. The applicability of the MP1 and MP2 schemes for treating DBOC is demonstrated in calculations for the atomization energies of benzene, naphthalene, anthracene, and tetracene. The corresponding corrections are surprisingly large (about 0.6 kJmol for benzene, 1.1 kJmol for naphthalene, 1.5 kJmol for anthracene, and 1.8 kJmol for tetracene) with the electron-correlation corrections reducing the corresponding Hartree-Fock self-consistent field values by 25%-30%.     

Tautomeric forms of azolide anions: vertical electron detachment energies and Dyson orbitals

Several decouplings of the electron propagator, including the relatively new P3+ approximation for the self-energy, have been used to calculate vertical electron detachment energies of tautomeric forms of closed-shell, pentagonal, aromatic anions in which ring carbons without bonds to hydrogens appear. This study extends previous work in which the most stable forms of anionic, five-member rings with one to five nitrogens were considered. Whereas the lowest electron detachment energies sometimes are assigned by Koopmans's theorem results to pi orbital vacancies, electron propagator calculations always obtain sigma orbital vacancies for the ground states of the doublet radicals. Higher electron detachment energies that correspond to excited doublets with pi vacancies also are presented. The predicted transition energies are in good agreement with low-intensity peaks in recent anion photoelectron spectra that have been assigned to less stable, tautomeric forms of these anions.     

Rearrangement correction to higher polarizabilities of two-electron systems

In this paper we include the rearrangement correction (discussed in the preceding paper) in a coupled Hartree–Fock (CHF) calculation of atomic hyperpolarizabilities and other related properties. We have studied the effect of these corrections on properties like electric dipole hyperpolarizabilities, uniform electric field quadrupole polarizabilities and shielding factors in two-electron ions and have noticed significant changes in the computed values over the CHF results.     

Analytic calculation of the diagonal Born-Oppenheimer correction within configuration-interaction and coupled-cluster theory

Schemes for the analytic calculation of the diagonal Born-Oppenheimer correction (DBOC) are formulated and implemented for use with general single-reference configuration-interaction and coupled-cluster wave function models. Calculations are reported to demonstrate the convergence of the DBOC with respect to electron-correlation treatment and basis set as well as to investigate the size-consistency error in configuration-interaction calculations of the DBOC. The importance of electron-correlation contributions to the DBOC is illustrated in the computation of the corresponding corrections for the reaction energy and activation barrier of the F + H2 --> FH + H reaction as well as of the atomization energy for trans-butadiene.     

Calculation of quasiparticle energy of molecular systems by the GW method

The quasiparticle energy of the H₂ molecule is calculated by using the GW method, in which the self‐energy operator fully depends on the frequency. The initial Green function G₀ is constructed from the wave function obtained by the Hartree–Fock approximation (HFA) and local density approximation (LDA) in the framework of the density functional theory (DFT). From the results obtained we have shown that the wave function from the DFT–LDA is more effective than that from the HFA for G₀. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem 84: 348–353, 2001     

Accurate lattice energies of organic molecular crystals from periodic turbomole calculations

Accurate lattice energies of organic crystals are important i.e. for the pharmaceutical industry. Periodic DFT calculations with atom‐centered Gaussian basis functions with the Turbomole program are used to calculate lattice energies for several non‐covalently bound organic molecular crystals. The accuracy and convergence of results with basis set size and k‐space sampling from periodic calculations is evaluated for the two reference molecules benzoic acid and naphthalene. For the X23 benchmark set of small molecular crystals accurate lattice energies are obtained using the PBE‐D3 functional. In particular for hydrogen‐bonded systems, a sufficiently large basis set is required. The calculated lattice energy differences between enantiopure and racemic crystal forms for a prototype set of chiral molecules are in good agreement with experimental results and allow the rationalization and computer‐aided design of chiral separation processes. © 2018 Wiley Periodicals, Inc.     

Hyperviral extension for energies of parameter-dependent systems

The energy of a linearly perturbed enclosed system is calculated by way of the hypervirial methodology in the semiclassical limit. Solutions of this equation are compared with the exact as well as with WKB approximate results for the harmonic oscillator model.     

On localized orbitals in infinite,periodic systems

It is shown that in two- or three-dimensional periodic systems with an odd number of electrons per unit cell no unitary transformation of the canonical Hartree - Fock orbitals can produce equivalent, localized orbitals for all electrons. However, in one-dimensional systems such orbitals can be found. The physical implications are discussed.     

Application of semiempirical long-range dispersion corrections to periodic systems in density functional theory

Ewald summation is used to apply semiempirical long-range dispersion corrections (Grimme, J Comput Chem 2006, 27, 1787; 2004, 25, 1463) to periodic systems in density functional theory. Using the parameters determined before for molecules and the Perdew-Burke-Ernzerhof functional, structure parameters and binding energies for solid methane, graphite, and vanadium pentoxide are determined in close agreement with observed values. For methane, a lattice constant a of 580 pm and a sublimation energy of 11 kJ mol(-1) are calculated. For the layered solids graphite and vanadia, the interlayer distances are 320 pm and 450 pm, respectively, whereas the graphite interlayer energy is -5.5 kJ mol(-1) per carbon atom and layer. Only when adding the semiempirical dispersion corrections, realistic values are obtained for the energies of adsorption of C(4) alkenes in microporous silica (-66 to -73 kJ mol(-1)) and the adsorption and chemisorption (alkoxide formation) of isobutene on acidic sites in the micropores of zeolite ferrierite (-78 to -94 kJ mol(-1)). As expected, errors due to missing self-interaction correction as in the energy for the proton transfer from the acidic site to the alkene forming a carbenium ion are not affected by the dispersion term. The adsorption and reaction energies are compared with the results from M?ller-Plesset second-order perturbation theory with basis set extrapolation.     

Negativity of ground-state interaction energies of atom-excess-electron systems

Any atom—excess-electron system confined to an arbitrary finite region of space that is free of external fields and which can be characterized by a non-relativistic spin-free hamiltonian has, in the limit of infinite nuclear mass of the atom, an exact ground-state expectation value of the atom—electron interaction hamiltonian that is negative-definite.     

Screened potential in π-electron systems and its application to the estimation of excitation energies

A diagrammatic technique was developed for the estimation of the screened potential of -electron systems. The screened potential was expanded in terms of the polarization propagators which were constructed from either the singlet, , or triplet vertex part, . These vertex parts correspond to the singlet or triplet excitations, respectively, in the Random Phase Approximation (RPA) containing exchange diagrams. The excitation energies were calculated by using the screened potential in the framework of RPA with exchange. The excitation energies of several conjugated molecules with or without a hetero atom are in agreement with the experimental data.     

Estimation, computation, and experimental correction of molecular zero-point vibrational energies

For accurate thermochemical tests of electronic structure theory, accurate true anharmonic zero-point vibrational energies ZPVE(true) are needed. We discuss several possibilities to extract this information for molecules from density functional or wave function calculations and/or available experimental data: (1) Empirical universal scaling of density-functional-calculated harmonic ZPVE(harm)s, where we find that polyatomics require smaller scaling factors than diatomics. (2) Direct density-functional calculation by anharmonic second-order perturbation theory PT2. (3) Weighted averages of harmonic ZPVE(harm) and fundamental ZPVE(fund) (from fundamental vibrational transition frequencies), with weights (3/4, 1/4) for diatomics and (5/8,3/8) for polyatomics. (4) Experimental correction of the PT2 harmonic contribution, i.e., the estimate ZPVE(true)PT2 + (ZPVE(fund)expt - ZPVE(fund)PT2) for ZPVE(true). The (5/8,3/8) average of method 3 and the additive correction of method 4 have been proposed here. For our database of experimental ZPVE(true), consisting of 27 diatomics and 8 polyatomics, we find that methods 1 and 2, applied to the popular B3LYP and the nonempirical PBE and TPSS functionals and their one-parameter hybrids, yield polyatomic errors on the order of 0.1 kcal/mol. Larger errors are expected for molecules larger than those in our database. Method 3 yields errors on the order of 0.02 kcal/mol, but requires very accurate (e.g., experimental, coupled cluster, or best-performing density functional) input harmonic ZPVE(harm). Method 4 is the best-founded one that meets the requirements of high accuracy and practicality, requiring as experimental input only the highly accurate and widely available ZPVE(fund)expt and producing errors on the order of 0.05 kcal/mol that are relatively independent of functional and basis set. As a part of our study, we also test the ability of the density functionals to predict accurate equilibrium bond lengths and angles for a data set of 21 mostly polyatomic molecules (since all calculated ZPVEs are evaluated at the correspondingly calculated molecular geometries).     

Implementation of empirical dispersion corrections to density functional theory for periodic systems

A recently developed empirical dispersion correction (Grimme et al., J. Chem. Phys. 2010, 132, 154104) to standard density functional theory (DFT‐D3) is implemented in the plane‐wave program package VASP. The DFT‐D3 implementation is compared with an implementation of the earlier DFT‐D2 version (Grimme, J. Comput. Chem. 2004, 25, 1463; Grimme, J. Comput. Chem. 2006, 27, 1787). Summation of empirical pair potential terms is performed over all atom pairs in the reference cell and over atoms in shells of neighboring cells until convergence of the dispersion energy is obtained. For DFT‐D3, the definition of coordination numbers has to be modified with respect to the molecular version to ensure convergence. The effect of three‐center terms as implemented in the original molecular DFT‐D3 version is investigated. The empirical parameters are taken from the original DFT‐D3 version where they had been optimized for a reference set of small molecules. As the coordination numbers of atoms in bulk and surfaces are much larger than in the reference compounds, this effect has to be discussed. The results of test calculations for bulk properties of metals, metal oxides, benzene, and graphite indicate that the original parameters are also suitable for solid‐state systems. In particular, the interlayer distance in bulk graphite and lattice constants of molecular crystals is considerably improved over standard functionals. With the molecular standard parameters (Grimme et al., J. Chem. Phys. 2010, 132, 154104; Grimme, J. Comput. Chem. 2006, 27, 1787) a slight overbinding is observed for ionic oxides where dispersion should not contribute to the bond. For simple adsorbate systems, such as Xe atoms and benzene on Ag(111), the DFT‐D implementations reproduce experimental results with a similar accuracy as more sophisticated approaches based on perturbation theory (Rohlfing and Bredow, Phys. Rev. Lett. 2008, 101, 266106). © 2012 Wiley Periodicals, Inc.     

On almost periodic processes in impulsive fractional-order competitive systems

In this paper an impulsive model for the dynamics of competitive Lotka–Volterra systems using the Caputo fractional-order derivative is developed. The existence and uniqueness of almost periodic solutions are investigated. Applying the fractional Lyapunov method, we give sufficient conditions for global perfect uniform-asymptotic stability of the almost periodic solution and sufficient conditions to save these qualities at the uncertain case. Since the competitive relationship in ecosystem models is relevant in various contexts, including many population, neural nets or chemical kinetics problems, our results can be applied in the investigation of almost periodic processes in a wide range of competitive systems of diverse interest.     

On solving the master equation in spatially periodic systems

We present a new method for solving the master equation for a system evolving on a spatially periodic network of states. The network contains 2(ν) images of a "unit cell" of n states, arranged along one direction with periodic boundary conditions at the ends. We analyze the structure of the symmetrized (2(ν)n) × (2(ν)n) rate constant matrix for this system and derive a recursive scheme for determining its eigenvalues and eigenvectors, and therefore analytically expressing the time-dependent probabilities of all states in the network, based on diagonalizations of n × n matrices formed by consideration of a single unit cell. We apply our new method to the problem of low-temperature, low-occupancy diffusion of xenon in the zeolite silicalite-1 using the states, interstate transitions, and transition state theory-based rate constants previously derived by June et al. [J. Phys. Chem. 95, 8866 (1991)]. The new method yields a diffusion tensor for this system which differs by less than 3% from the values derived previously via kinetic Monte Carlo (KMC) simulations and confirmed by new KMC simulations conducted in the present work. The computational requirements of the new method are compared against those of KMC, numerical solution of the master equation by the Euler method, and direct molecular dynamics. In the problem of diffusion of xenon in silicalite-1, the new method is shown to be faster than these alternative methods by factors of about 3.177 × 10(4), 4.237 × 10(3), and 1.75 × 10(7), respectively. The computational savings and ease of setting up calculations afforded by the new method of master equation solution by recursive reduction of dimensionality in diagonalizing the rate constant matrix make it attractive as a means of predicting long-time dynamical phenomena in spatially periodic systems from atomic-level information.     

Stress tensor in model polymer systems with periodic boundaries

The calculation of the stress tensor from molecular simulations of atomistic model polymer systems employing periodic boundary conditions is discussed. Starting from the dynamical equations governing the motion of sites, correct double summation forms of the atomic and the molecular virial equations are derived, which are valid for flexible, infinitely stiff and rigid chain models even in the presence of interactions between different images of the same parent macromolecule. A new expression for the true instantaneous stress (flux of momentum through the faces of the simulation box) is derived and shown to exhibit large fluctuations when applied in molecular dynamics simulations. A new equation for the thermodynamic stress, cast exclusively in terms of intermolecular forces on interaction sites, is also derived. Application to Monte Carlo simulations shows that the molecular virial expression exhibits the smallest fluctuations among all stress expressions discussed, and thus allows computation of the thermodynamic stress with least uncertainty. A scheme is developed for the calculation of surface tension from intermolecular forces only.     

Three-dimensional Ewald method with correction term for a system periodic in one direction

A three-dimensional Ewald summation formula with a shape-dependent correction term for Coulomb interactions in systems with one-dimensional periodicity is derived. Test molecular dynamics simulations of acetone molecules in cylindrical silica pores show that the formula is efficient only when size of the system in a plane perpendicular to the periodicity direction is small in comparison with the periodicity length.