Showing papers in "Journal of Chemical Physics in 2014"
TL;DR: This work implements an implicit solvation model that has a firm theoretical foundation into the widely used density-functional code Vienna ab initio Software Package and finds that solvation reduces the surface energies of the nanocrystals and increases the energy barrier of the SN2 reaction.
Abstract: Solid-liquid interfaces are at the heart of many modern-day technologies and provide a challenge to many materials simulation methods. A realistic first-principles computational study of such systems entails the inclusion of solvent effects. In this work, we implement an implicit solvation model that has a firm theoretical foundation into the widely used density-functional code Vienna ab initio Software Package. The implicit solvation model follows the framework of joint density functional theory. We describe the framework, our algorithm and implementation, and benchmarks for small molecular systems. We apply the solvation model to study the surface energies of different facets of semiconducting and metallic nanocrystals and the SN2 reaction pathway. We find that solvation reduces the surface energies of the nanocrystals, especially for the semiconducting ones and increases the energy barrier of the SN2 reaction.
1,378 citations
TL;DR: This Perspective outlines the history, basic concepts, and present status of KS-DFT in chemical physics, and offers suggestions for its future development.
Abstract: Since its formal inception in 1964–1965, Kohn-Sham density-functional theory (KS-DFT) has become the most popular electronic structure method in computational physics and chemistry. Its popularity stems from its beautifully simple conceptual framework and computational elegance. The rise of KS-DFT in chemical physics began in earnest in the mid 1980s, when crucial developments in its exchange-correlation term gave the theory predictive power competitive with well-developed wave-function methods. Today KS-DFT finds itself under increasing pressure to deliver higher and higher accuracy and to adapt to ever more challenging problems. If we are not mindful, however, these pressures may submerge the theory in the wave-function sea. KS-DFT might be lost. I am hopeful the Kohn-Sham philosophical, theoretical, and computational framework can be preserved. This Perspective outlines the history, basic concepts, and present status of KS-DFT in chemical physics, and offers suggestions for its future development.
1,136 citations
TL;DR: A systematic examination of the computational expense and accuracy of Symmetry-Adapted Perturbation Theory for the prediction of non-covalent interaction energies is provided and three standards are recommended as the gold, silver, and bronze standard of SAPT.
Abstract: A systematic examination of the computational expense and accuracy of Symmetry-Adapted Perturbation Theory (SAPT) for the prediction of non-covalent interaction energies is provided with respect to both method [SAPT0, DFT-SAPT, SAPT2, SAPT2+, SAPT2+(3), and SAPT2+3; with and without CCD dispersion for the last three] and basis set [Dunning cc-pVDZ through aug-cc-pV5Z wherever computationally tractable, including truncations of diffuse basis functions]. To improve accuracy for hydrogen-bonded systems, we also include two corrections based on exchange-scaling (sSAPT0) and the supermolecular MP2 interaction energy (δMP2). When considering the best error performance relative to computational effort, we recommend as the gold, silver, and bronze standard of SAPT: SAPT2+(3)δMP2/aug-cc-pVTZ, SAPT2+/aug-cc-pVDZ, and sSAPT0/jun-cc-pVDZ. Their respective mean absolute errors in interaction energy across the S22, HBC6, NBC10, and HSG databases are 0.15 (62.9), 0.30 (4.4), and 0.49 kcal mol−1 (0.03 h for adenine·thymine complex).
554 citations
TL;DR: An extensive structural study on solid H2S at pressure ranges of 10-200 GPa is performed through an unbiased structure prediction method based on particle swarm optimization algorithm and is able to establish stable metallic structures violating an earlier proposal of elemental decomposition into sulfur and hydrogen.
Abstract: Hydrogen sulfide (H2S) is a prototype molecular system and a sister molecule of water (H2O). The phase diagram of solid H2S at high pressures remains largely unexplored arising from the challenges in dealing with the pressure-induced weakening of S–H bond and larger atomic core difference between H and S. Metallization is yet achieved for H2O, but it was observed for H2S above 96 GPa. However, the metallic structure of H2S remains elusive, greatly impeding the understanding of its metallicity and the potential superconductivity. We have performed an extensive structural study on solid H2S at pressure ranges of 10–200 GPa through an unbiased structure prediction method based on particle swarm optimization algorithm. Besides the findings of candidate structures for nonmetallic phases IV and V, we are able to establish stable metallic structures violating an earlier proposal of elemental decomposition into sulfur and hydrogen [R. Rousseau, M. Boero, M. Bernasconi, M. Parrinello, and K. Terakura, Phys. Rev. Lett. 85, 1254 (2000)]. Our study unravels a superconductive potential of metallic H2S with an estimated maximal transition temperature of ∼80 K at 160 GPa, higher than those predicted for most archetypal hydrogen-containing compounds (e.g., SiH4, GeH4, etc.).
549 citations
TL;DR: This work builds upon the previously developed many-body dispersion framework and proposes an effective range-separation of the coupling between the atomic response functions that extends the already broad applicability of the MBD method to non-metallic materials with highly anisotropic responses, such as layered nanostructures.
Abstract: An accurate determination of the electron correlation energy is an essential prerequisite for describing the structure, stability, and function in a wide variety of systems. Therefore, the development of efficient approaches for the calculation of the correlation energy (and hence the dispersion energy as well) is essential and such methods can be coupled with many density-functional approximations, local methods for the electron correlation energy, and even interatomic force fields. In this work, we build upon the previously developed many-body dispersion (MBD) framework, which is intimately linked to the random-phase approximation for the correlation energy. We separate the correlation energy into short-range contributions that are modeled by semi-local functionals and long-range contributions that are calculated by mapping the complex all-electron problem onto a set of atomic response functions coupled in the dipole approximation. We propose an effective range-separation of the coupling between the atomic response functions that extends the already broad applicability of the MBD method to non-metallic materials with highly anisotropic responses, such as layered nanostructures. Application to a variety of high-quality benchmark datasets illustrates the accuracy and applicability of the improved MBD approach, which offers the prospect of first-principles modeling of large structurally complex systems with an accurate description of the long-range correlation energy.
509 citations
TL;DR: New analysis tools for excited states are introduced including state averaged natural transition orbitals, which give a compact description of a number of states simultaneously, and natural difference orbitals (defined as the eigenvectors of the difference density matrix), which reveal details about orbital relaxation effects.
Abstract: A variety of density matrix based methods for the analysis and visualization of electronic excitations are discussed and their implementation within the framework of the algebraic diagrammatic construction of the polarization propagator is reported. Their mathematical expressions are given and an extensive phenomenological discussion is provided to aid the interpretation of the results. Starting from several standard procedures, e.g., population analysis, natural orbital decomposition, and density plotting, we proceed to more advanced concepts of natural transition orbitals and attachment/detachment densities. In addition, special focus is laid on information coded in the transition density matrix and its phenomenological analysis in terms of an electron-hole picture. Taking advantage of both the orbital and real space representations of the density matrices, the physical information in these analysis methods is outlined, and similarities and differences between the approaches are highlighted. Moreover, new analysis tools for excited states are introduced including state averaged natural transition orbitals, which give a compact description of a number of states simultaneously, and natural difference orbitals (defined as the eigenvectors of the difference density matrix), which reveal details about orbital relaxation effects.
364 citations
TL;DR: The electronic and thermoelectric properties of one to four monolayers of MoS2, MoSe2, WS2, and WSe2 are calculated and the transition from few layers to bulk is discussed.
Abstract: The electronic and thermoelectric properties of one to four monolayers of MoS2, MoSe2, WS2, and WSe2 are calculated For few layer thicknesses, the near degeneracies of the conduction band K and Σ valleys and the valence band Γ and K valleys enhance the n-type and p-type thermoelectric performance The interlayer hybridization and energy level splitting determine how the number of modes within kBT of a valley minimum changes with layer thickness In all cases, the maximum ZT coincides with the greatest near-degeneracy within kBT of the band edge that results in the sharpest turn-on of the density of modes The thickness at which this maximum occurs is, in general, not a monolayer The transition from few layers to bulk is discussed Effective masses, energy gaps, power-factors, and ZT values are tabulated for all materials and layer thicknesses
321 citations
TL;DR: The computational effort in finding the minimum energy path with DFT was reduced by a factor ranging from 50% to an order of magnitude by using an IDPP path as the initial path.
Abstract: A method is presented for generating a good initial guess of a transition path between given initial and final states of a system without evaluation of the energy. An objective function surface is constructed using an interpolation of pairwise distances at each discretization point along the path and the nudged elastic band method then used to find an optimal path on this image dependent pair potential (IDPP) surface. This provides an initial path for the more computationally intensive calculations of a minimum energy path on an energy surface obtained, for example, by ab initio or density functional theory. The optimal path on the IDPP surface is significantly closer to a minimum energy path than a linear interpolation of the Cartesian coordinates and, therefore, reduces the number of iterations needed to reach convergence and averts divergence in the electronic structure calculations when atoms are brought too close to each other in the initial path. The method is illustrated with three examples: (1) rotation of a methyl group in an ethane molecule, (2) an exchange of atoms in an island on a crystal surface, and (3) an exchange of two Si-atoms in amorphous silicon. In all three cases, the computational effort in finding the minimum energy path with DFT was reduced by a factor ranging from 50% to an order of magnitude by using an IDPP path as the initial path. The time required for parallel computations was reduced even more because of load imbalance when linear interpolation of Cartesian coordinates was used.
284 citations
TL;DR: According to calculations, magnesium exhibits a tendency towards the growth of smooth surfaces as it exhibits lower diffusion barriers than lithium and sodium, and as an hcp metal it favors higher-coordinated configurations in contrast to the bcc metals Li and Na.
Abstract: Lithium and magnesium exhibit rather different properties as battery anode materials with respect to the phenomenon of dendrite formation which can lead to short-circuits in batteries. Diffusion processes are the key to understanding structure forming processes on surfaces. Therefore, we have determined adsorption energies and barriers for the self-diffusion on Li and Mg using periodic density functional theory calculations and contrasted the results to Na which is also regarded as a promising electrode material in batteries. According to our calculations, magnesium exhibits a tendency towards the growth of smooth surfaces as it exhibits lower diffusion barriers than lithium and sodium, and as an hcp metal it favors higher-coordinated configurations in contrast to the bcc metals Li and Na. These characteristic differences are expected to contribute to the unequal tendencies of these metals with respect to dendrite growth.
278 citations
TL;DR: This work found that the collective inclusion of Exx and vdW as resulting from a large-scale AIMD simulation of (H2O)128 significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, SOO(Q), and corresponding oxygen- oxygengen radial distribution function, gOO(r), that are now in quantitative agreement with the best available experimental data.
Abstract: In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local van der Waals/dispersion (vdW) interactions, via a fully self-consistent density-dependent dispersion correction, and an approximate treatment of nuclear quantum effects, via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx and vdW as resulting from a large-scale AIMD simulation of (H2O)128 significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, SOO(Q), and corresponding oxygen-oxygen radial distribution function, gOO(r), that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of a hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions. This increasingly more accurate description of the underlying hydrogen bond network in liquid water also yields higher-order correlation functions, such as the oxygen-oxygen-oxygen triplet angular distribution, POOO(θ), and therefore the degree of local tetrahedrality, as well as electrostatic properties, such as the effective molecular dipole moment, that are in much better agreement with experiment.
266 citations
TL;DR: The use of Mie potentials to describe the group-group interaction is shown to allow accurate simultaneous descriptions of the fluid-phase behavior and second-order thermodynamic derivative properties of the pure fluids based on a single set of group parameters.
Abstract: A generalization of the recent version of the statistical associating fluid theory for variable range Mie potentials [Lafitte et al., J. Chem. Phys. 139, 154504 (2013)] is formulated within the framework of a group contribution approach (SAFT-γ Mie). Molecules are represented as comprising distinct functional (chemical) groups based on a fused heteronuclear molecular model, where the interactions between segments are described with the Mie (generalized Lennard-Jonesium) potential of variable attractive and repulsive range. A key feature of the new theory is the accurate description of the monomeric group-group interactions by application of a high-temperature perturbation expansion up to third order. The capabilities of the SAFT-γ Mie approach are exemplified by studying the thermodynamic properties of two chemical families, the n-alkanes and the n-alkyl esters, by developing parameters for the methyl, methylene, and carboxylate functional groups (CH3, CH2, and COO). The approach is shown to describe accurately the fluid-phase behavior of the compounds considered with absolute average deviations of 1.20% and 0.42% for the vapor pressure and saturated liquid density, respectively, which represents a clear improvement over other existing SAFT-based group contribution approaches. The use of Mie potentials to describe the group-group interaction is shown to allow accurate simultaneous descriptions of the fluid-phase behavior and second-order thermodynamic derivative properties of the pure fluids based on a single set of group parameters. Furthermore, the application of the perturbation expansion to third order for the description of the reference monomeric fluid improves the predictions of the theory for the fluid-phase behavior of pure components in the near-critical region. The predictive capabilities of the approach stem from its formulation within a group-contribution formalism: predictions of the fluid-phase behavior and thermodynamic derivative properties of compounds not included in the development of group parameters are demonstrated. The performance of the theory is also critically assessed with predictions of the fluid-phase behavior (vapor-liquid and liquid-liquid equilibria) and excess thermodynamic properties of a variety of binary mixtures, including polymer solutions, where very good agreement with the experimental data is seen, without the need for adjustable mixture parameters.
TL;DR: The excited states of a diverse set of molecules are examined using a collection of newly implemented analysis methods to expose the particular power of natural difference orbitals for the description of orbital relaxation effects and the utility of a wide array of additional analysis methods is highlighted.
Abstract: The excited states of a diverse set of molecules are examined using a collection of newly implemented analysis methods. These examples expose the particular power of three of these tools: (i) natural difference orbitals (the eigenvectors of the difference density matrix) for the description of orbital relaxation effects, (ii) analysis of the one-electron transition density matrix in terms of an electron-hole picture to identify charge resonance and excitonic correlation effects, and (iii) state-averaged natural transition orbitals for a compact simultaneous representation of several states. Furthermore, the utility of a wide array of additional analysis methods is highlighted. Five molecules with diverse excited state characteristics are chosen for these tasks: pyridine as a prototypical small heteroaromatic molecule, a model system of six neon atoms to study charge resonance effects, hexatriene in its neutral and radical cation forms to exemplify the cases of double excitations and spin-polarization, respectively, and a model iridium complex as a representative metal organic compound. Using these examples a number of phenomena, which are at first sight unexpected, are highlighted and their physical significance is discussed. Moreover, the generality of the conclusions of this paper is verified by a comparison of single- and multireference ab initio methods.
TL;DR: Numerically how the approximation to quantum dynamics depends on this friction is explored, with a particular emphasis on vibrational spectroscopy: a broad range of frictions approaching optimal damping give similar results, and these results are immune to both the resonance problem of RPMD and the curvature problem of CMD.
Abstract: Two of the most successful methods that are presently available for simulating the quantum dynamics of condensed phase systems are centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD). Despite their conceptual differences, practical implementations of these methods differ in just two respects: the choice of the Parrinello-Rahman mass matrix and whether or not a thermostat is applied to the internal modes of the ring polymer during the dynamics. Here, we explore a method which is halfway between the two approximations: we keep the path integral bead masses equal to the physical particle masses but attach a Langevin thermostat to the internal modes of the ring polymer during the dynamics. We justify this by showing analytically that the inclusion of an internal mode thermostat does not affect any of the established features of RPMD: thermostatted RPMD is equally valid with respect to everything that has actually been proven about the method as RPMD itself. In particular, because of the choice of bead masses, the resulting method is still optimum in the short-time limit, and the transition state approximation to its reaction rate theory remains closely related to the semiclassical instanton approximation in the deep quantum tunneling regime. In effect, there is a continuous family of methods with these properties, parameterised by the strength of the Langevin friction. Here, we explore numerically how the approximation to quantum dynamics depends on this friction, with a particular emphasis on vibrational spectroscopy. We find that a broad range of frictions approaching optimal damping give similar results, and that these results are immune to both the resonance problem of RPMD and the curvature problem of CMD.
TL;DR: This work reports results from molecular dynamics simulations that indicate that NPG maintains its ultrahigh permeability even at low pressures, allowing a permeate water flux of 6.1 × 10−15 l/h bar per pore.
Abstract: Nanoporous graphene (NPG) shows tremendous promise as an ultra-permeable membrane for water desalination thanks to its atomic thickness and precise sieving properties. However, a significant gap exists in the literature between the ideal conditions assumed for NPG desalination and the physical environment inherent to reverse osmosis (RO) systems. In particular, the water permeability of NPG has been calculated previously based on very high pressures (1000–2000 bars). Does NPG maintain its ultrahigh water permeability under real-world RO pressures (<100 bars)? Here, we answer this question by drawing results from molecular dynamics simulations. Our results indicate that NPG maintains its ultrahigh permeability even at low pressures, allowing a permeate water flux of 6.0 l/h-bar per pore, or equivalently 1041 ± 20 l/m2-h-bar assuming a nanopore density of 1.7 × 1013 cm−2.
TL;DR: In this paper, the authors introduced a database (HAB11) of electronic coupling matrix elements (Hab) for electron transfer in 11 π-conjugated organic homo-dimer cations.
Abstract: We introduce a database (HAB11) of electronic coupling matrix elements (Hab) for electron transfer in 11 π-conjugated organic homo-dimer cations. High-level ab inito calculations at the multireference configuration interaction MRCI+Q level of theory, n-electron valence state perturbation theory NEVPT2, and (spin-component scaled) approximate coupled cluster model (SCS)-CC2 are reported for this database to assess the performance of three DFT methods of decreasing computational cost, including constrained density functional theory (CDFT), fragment-orbital DFT (FODFT), and self-consistent charge density functional tight-binding (FODFTB). We find that the CDFT approach in combination with a modified PBE functional containing 50% Hartree-Fock exchange gives best results for absolute Hab values (mean relative unsigned error = 5.3%) and exponential distance decay constants β (4.3%). CDFT in combination with pure PBE overestimates couplings by 38.7% due to a too diffuse excess charge distribution, whereas the ec...
TL;DR: The implementation of an efficient program of the algebraic diagrammatic construction method for the polarisation propagator in third-order perturbation theory (ADC(3)) for the computation of excited states is reported and it is reported that ADC(3) has a much larger range of applicability due to its more favourable scaling of O(N(6)) with system size.
Abstract: The implementation of an efficient program of the algebraic diagrammatic construction method for the polarisation propagator in third-order perturbation theory (ADC(3)) for the computation of excited states is reported. The accuracies of ADC(2) and ADC(3) schemes have been investigated with respect to Thiel's recently established benchmark set for excitation energies and oscillator strengths. The calculation of 141 vertical excited singlet and 71 triplet states of 28 small to medium-sized organic molecules has revealed that ADC(3) exhibits mean error and standard deviation of 0.12 ± 0.28 eV for singlet states and −0.18 ± 0.16 eV for triplet states when the provided theoretical best estimates are used as benchmark. Accordingly, the ADC(2)-s and ADC(2)-x calculations revealed accuracies of 0.22 ± 0.38 eV and −0.70 ± 0.37 eV for singlets and 0.12 ± 0.16 eV and −0.55 ± 0.20 eV for triplets, respectively. For a comparison of CC3 and ADC(3), only non-CC3 benchmark values were considered, which comprise 84 singlet states and 19 triplet states. For these singlet states CC3 exhibits an accuracy of 0.23 ± 0.21 eV and ADC(3) an accuracy of 0.08 ± 0.27 eV, and accordingly for the triplet states of 0.12 ± 0.10 eV and −0.10 ± 0.13 eV, respectively. Hence, based on the quality of the existing benchmark set it is practically not possible to judge whether ADC(3) or CC3 is more accurate, however, ADC(3) has a much larger range of applicability due to its more favourable scaling of O(N6) with system size.
TL;DR: This work extends previous work on singlet exciton fission in isolated dimers to the case of crystalline materials, focusing on pentacene as a canonical and concrete example, demonstrating significant charge-transfer character in the low-lying excited states.
Abstract: We extend our previous work on singlet exciton fission in isolated dimers to the case of crystalline materials, focusing on pentacene as a canonical and concrete example. We discuss the proper interpretation of the character of low-lying excited states of relevance to singlet fission. In particular, we consider a variety of metrics for measuring charge-transfer character, conclusively demonstrating significant charge-transfer character in the low-lying excited states. The impact of this electronic structure on the subsequent singlet fission dynamics is assessed by performing real-time master-equation calculations involving hundreds of quantum states. We make direct comparisons with experimental absorption spectra and singlet fission rates, finding good quantitative agreement in both cases, and we discuss the mechanistic distinctions that exist between small isolated aggregates and bulk systems.
TL;DR: A size-extensive construction of SIC orbitals is proposed which, unlike earlier constructions, makes SIC computationally efficient, and a true spin-density functional.
Abstract: Standard spin-density functionals for the exchange-correlation energy of a many-electron ground state make serious self-interaction errors which can be corrected by the Perdew-Zunger self-interaction correction (SIC). We propose a size-extensive construction of SIC orbitals which, unlike earlier constructions, makes SIC computationally efficient, and a true spin-density functional. The SIC orbitals are constructed from a unitary transformation that is explicitly dependent on the non-interacting one-particle density matrix. When this SIC is applied to the local spin-density approximation, improvements are found for the atomization energies of molecules.
TL;DR: Comparing the performance between the two classes of vdW-inclusive methods for adsorption at equilibrium adsorbate-substrate separations, it is concluded that both methods allow for an accurate treatment of adsorptive energies, equilibrium distances, and binding curves.
Abstract: Exploring the role of van der Waals (vdW) forces on the adsorption of molecules on extended metal surfaces has become possible in recent years thanks to exciting developments in density functional theory (DFT). Among these newly developed vdW-inclusive methods, interatomic vdW approaches that account for the nonlocal screening within the bulk [V. G. Ruiz, W. Liu, E. Zojer, M. Scheffler, and A. Tkatchenko, Phys. Rev. Lett. 108, 146103 (2012)] and improved nonlocal functionals [J. Klimes, D. R. Bowler, and A. Michaelides, J. Phys.: Condens. Matter 22, 022201 (2010)] have emerged as promising candidates to account efficiently and accurately for the lack of long-range vdW forces in most popular DFT exchange-correlation functionals. Here we have used these two approaches to compute benzene adsorption on a range of close-packed (111) surfaces upon which it either physisorbs (Cu, Ag, and Au) or chemisorbs (Rh, Pd, Ir, and Pt). We have thoroughly compared the performance between the two classes of vdW-inclusive methods and when available compared the results obtained with experimental data. By examining the computed adsorption energies, equilibrium distances, and binding curves we conclude that both methods allow for an accurate treatment of adsorption at equilibrium adsorbate-substrate distances. To this end, explicit inclusion of electrodynamic screening in the interatomic vdW scheme and optimized exchange functionals in the case of nonlocal vdW density functionals is mandatory. Nevertheless, some discrepancies are found between these two classes of methods at large adsorbate-substrate separations.
TL;DR: It is demonstrated that only the TS/HI method leads to accurate results for systems where both electrostatic and dispersion interactions are important, as illustrated for Li-intercalated graphite and for molecular adsorption on the surfaces in ionic solids and in the cavities of zeolites.
Abstract: Recently we have demonstrated that the applicability of the Tkatchenko-Scheffler (TS) method for calculating dispersion corrections to density-functional theory can be extended to ionic systems if the Hirshfeld method for estimating effective volumes and charges of atoms in molecules or solids (AIM’s) is replaced by its iterative variant [T. Bucko, S. Lebegue, J. Hafner, and J. Angyan, J. Chem. Theory Comput.9, 4293 (2013)]. The standard Hirshfeld method uses neutral atoms as a reference, whereas in the iterative Hirshfeld (HI) scheme the fractionally charged atomic reference states are determined self-consistently. We show that the HI method predicts more realistic AIM charges and that the TS/HI approach leads to polarizabilities and C 6 dispersion coefficients in ionic or partially ionic systems which are, as expected, larger for anions than for cations (in contrast to the conventional TS method). For crystalline materials, the new algorithm predicts polarizabilities per unit cell in better agreement with the values derived from the Clausius-Mosotti equation. The applicability of the TS/HI method has been tested for a wide variety of molecular and solid-state systems. It is demonstrated that for systems dominated by covalent interactions and/or dispersion forces the TS/HI method leads to the same results as the conventional TS approach. The difference between the TS/HI and TS approaches increases with increasing ionicity. A detailed comparison is presented for isoelectronic series of octet compounds, layered crystals, complex intermetallic compounds, and hydrides, and for crystals built of molecules or containing molecular anions. It is demonstrated that only the TS/HI method leads to accurate results for systems where both electrostatic and dispersion interactions are important, as illustrated for Li-intercalated graphite and for molecular adsorption on the surfaces in ionic solids and in the cavities of zeolites.
TL;DR: Out of the three water models, TIP5P water shows the largest degree of tetrahedral ordering and is often the one that is least perturbed by the presence of the wall, while in the hydrophilic regime there are noticeable differences in structure and dynamics between all stages of wall roughness.
Abstract: We perform molecular dynamics simulations to observe the structure and dynamics of water using different water models (TIP3P, TIP4P, TIP5P) at ambient conditions, constrained by planar walls, which are either modeled by smooth potentials or regular atomic lattices, imitating the honeycomb-structure of graphene. We implement walls of different hydroaffinity, different lattice constant, and different types of interaction with the water molecules. We find that in the hydrophobic regime the smooth wall generally represents a good abstraction of the atomically rough walls, while in the hydrophilic regime there are noticeable differences in structure and dynamics between all stages of wall roughness. For a small lattice constant however the smooth and the atomically rough wall still share a number of structural and dynamical similarities. Out of the three water models, TIP5P water shows the largest degree of tetrahedral ordering and is often the one that is least perturbed by the presence of the wall.
TL;DR: A new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of abinitio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest methods, and explores the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost.
Abstract: We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as “cloning,” in analogy to the “spawning” procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, “trains,” as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.
TL;DR: It is shown how one can pair the orbitals so that the role of the Brueckner orbitals at the CCD level is retained at the p-CCD level, and ways of extending CCD to accurately describe strongly correlated systems are explored.
Abstract: Coupled cluster theory with single and double excitations accurately describes weak electron correlation but is known to fail in cases of strong static correlation. Fascinatingly, however, pair coupled cluster doubles (p-CCD), a simplified version of the theory limited to pair excitations that preserve the seniority of the reference determinant (i.e., the number of unpaired electrons), has mean field computational cost and is an excellent approximation to the full configuration interaction (FCI) of the paired space provided that the orbital basis defining the pairing scheme is adequately optimized. In previous work, we have shown that optimization of the pairing scheme in the seniority zero FCI leads to a very accurate description of static correlation. The same conclusion extends to p-CCD if the orbitals are optimized to make the p-CCD energy stationary. We here demonstrate these results with numerous examples. We also explore the contributions of different seniority sectors to the coupled cluster doubles (CCD) correlation energy using different orbital bases. We consider both Hartree-Fock and Brueckner orbitals, and the role of orbital localization. We show how one can pair the orbitals so that the role of the Brueckner orbitals at the CCD level is retained at the p-CCD level. Moreover, we explore ways of extending CCD to accurately describe strongly correlated systems.
TL;DR: A newly developed publicly available genetic algorithm for global structure optimisation within atomic scale modeling and the implementation is described and benchmarked through a detailed statistical analysis employing averages across many independent runs of the GA.
Abstract: We present a newly developed publicly available genetic algorithm (GA) for global structure optimisation within atomic scale modeling. The GA is focused on optimizations using first principles calculations, but it works equally well with empirical potentials. The implementation is described and benchmarked through a detailed statistical analysis employing averages across many independent runs of the GA. This analysis focuses on the practical use of GA’s with a description of optimal parameters to use. New results for the adsorption of M8 clusters (M = Ru, Rh, Pd, Ag, Pt, Au) on the stoichiometric rutile TiO2(110) surface are presented showing the power of automated structure prediction and highlighting the diversity of metal cluster geometries at the atomic scale.
TL;DR: In this article, surface roughness influences the adhesion between elastic solids and a Tabor number which depends on the length scale or magnification is introduced to give information about the nature of adhesion at different length scales.
Abstract: We discuss how surface roughness influences the adhesion between elastic solids. We introduce a Tabor number which depends on the length scale or magnification, and which gives information about the nature of the adhesion at different length scales. We consider two limiting cases relevant for (a) elastically hard solids with weak (or long ranged) adhesive interaction (DMT-limit) and (b) elastically soft solids with strong (or short ranged) adhesive interaction (JKR-limit). For the former cases we study the nature of the adhesion using different adhesive force laws (F ∼ u−n, n = 1.5–4, where u is the wall-wall separation). In general, adhesion may switch from DMT-like at short length scales to JKR-like at large (macroscopic) length scale. We compare the theory predictions to results of exact numerical simulations and find good agreement between theory and simulation results.
TL;DR: This work focuses on first-principles crystal structure prediction at high pressures, paying particular attention to novel high pressure structures uncovered by efficient structure prediction methods.
Abstract: Crystal structure prediction at high pressures unbiased by any prior known structure information has recently become a topic of considerable interest. We here present a short overview of recently developed structure prediction methods and propose current challenges for crystal structure prediction. We focus on first-principles crystal structure prediction at high pressures, paying particular attention to novel high pressure structures uncovered by efficient structure prediction methods. Finally, a brief perspective on the outstanding issues that remain to be solved and some directions for future structure prediction researches at high pressure are presented and discussed.
TL;DR: The results suggest that it may be possible to observe in computer simulations spontaneous crystallization of TIP4P/2005 at about 60 K below the melting point, and that for all water models γ decreases as the temperature decreases.
Abstract: In this work, we evaluate by means of computer simulations the rate for ice homogeneous nucleation for several water models such as TIP4P, TIP4P/2005,TIP4P/ICE, and mW (following the same procedure as in Sanz et al. [J. Am. Chem. Soc. 135, 15008 (2013)]) in a broad temperature range. We estimate the ice-liquid interfacial free-energy, and conclude that for all water models γ decreases as the temperature decreases. Extrapolating our results to the melting temperature, we obtain a value of the interfacial free-energy between 25 and 32 mN/m in reasonable agreement with the reported experimental values. Moreover, we observe that the values of γ depend on the chosen water model and this is a key factor when numerically evaluating nucleation rates, given that the kinetic prefactor is quite similar for all water models with the exception of the mW (due to the absence of hydrogens). Somewhat surprisingly the estimates of the nucleation rates found in this work for TIP4P/2005 are slightly higher than those of the mW model, even though the former has explicit hydrogens. Our results suggest that it may be possible to observe in computer simulations spontaneous crystallization of TIP4P/2005 at about 60 K below the melting point.
TL;DR: The ability of the CPM electrode to respond to local charge fluctuations in the electrolyte is seen to significantly lower the energy (and barrier) for the approach of Li(+) ions to the electrode surface.
Abstract: A major challenge in the molecular simulation of electric double layer capacitors (EDLCs) is the choice of an appropriate model for the electrode. Typically, in such simulations the electrode surface is modeled using a uniform fixed charge on each of the electrode atoms, which ignores the electrode response to local charge fluctuations in the electrolyte solution. In this work, we evaluate and compare this Fixed Charge Method (FCM) with the more realistic Constant Potential Method (CPM), [S. K. Reed et al., J. Chem. Phys. 126, 084704 (2007)], in which the electrode charges fluctuate in order to maintain constant electric potential in each electrode. For this comparison, we utilize a simplified LiClO4-acetonitrile/graphite EDLC. At low potential difference (ΔΨ ⩽ 2 V), the two methods yield essentially identical results for ion and solvent density profiles; however, significant differences appear at higher ΔΨ. At ΔΨ ⩾ 4 V, the CPM ion density profiles show significant enhancement (over FCM) of “inner-sphere...
TL;DR: In this paper, quantum Monte Carlo calculations were performed to study the cohesive energy of carbon allotropes, including sp 3-bonded diamond, sp 2-bunded graphene, sp-sp 2 hybridized graphynes, and sp-bonding carbyne, and it was shown that γ-graphyne, the most energetically stable graphyne, turns out to be 6.766(6) eV/atom, which is smaller than that of graphene by 0.698(12)
Abstract: We have performed quantum Monte Carlo calculations to study the cohesion energetics of carbon allotropes, including sp 3-bonded diamond, sp 2-bonded graphene, sp–sp 2 hybridized graphynes, and sp-bonded carbyne. The computed cohesive energies of diamond and graphene are found to be in excellent agreement with the corresponding values determined experimentally for diamond and graphite, respectively, when the zero-point energies, along with the interlayer binding in the case of graphite, are included. We have also found that the cohesive energy of graphyne decreases systematically as the ratio of sp-bonded carbon atoms increases. The cohesive energy of γ-graphyne, the most energetically stable graphyne, turns out to be 6.766(6) eV/atom, which is smaller than that of graphene by 0.698(12) eV/atom. Experimental difficulty in synthesizing graphynes could be explained by their significantly smaller cohesive energies. Finally, we conclude that the cohesive energy of a newly proposed graphyne can be accurately estimated with the carbon–carbon bond energies determined from the cohesive energies of graphene and three different graphynes considered here.
TL;DR: Six recently developed exchange functionals for pairing with two different versions of van der Waals density functionals (vdW-DF) are tested for weakly bonded solids, and the overall performance for lattice constants is good.
Abstract: Six recently developed exchange functionals for pairing with two different versions of van der Waals density functionals (vdW-DF) are tested for weakly bonded solids. The test, using 26 layered weakly bonded compounds, benchmarks the lattice constants against experimental data and the interlayer binding energies against reference data from the random-phase approximation (RPA). The investigated functionals tend to give interlayer binding energies higher than the RPA benchmark, and the overall performance for lattice constants is good. The exchange functionals optB86b and cx13 paired with the original vdW-DF and the B86R functional paired with vdW-DF2 are found to give particularly good results for equilibrium geometries.