Search

Simulations of condensed matter systems at the hybrid density functional theory level pose significant computational challenges. The elevated costs arise from the non-local nature of the Hartree–Fock exchange (HFX) in conjunction with the necessity to approach the thermodynamic limit. In this work, we address these issues with the development of a new efficient method for the calculation of HFX in periodic systems, employing k-point sampling. We rely on a local atom-specific resolution-of-the-identity scheme, the use of atom-centered Gaussian type orbitals, and the truncation of the Coulomb interaction to limit computational complexity. Our real-space approach exhibits a scaling that is, at worst, linear with the number of k-points. Issues related to basis set diffuseness are effectively addressed through the auxiliary density matrix method. We report the implementation in the CP2K software package, as well as accuracy and performance benchmarks. This method demonstrates excellent agreement with equivalent Γ-point supercell calculations in terms of relative energies and nuclear gradients. Good strong and weak scaling performances, as well as graphics processing unit (GPU) acceleration, make this implementation a promising candidate for high-performance computing. [ABSTRACT FROM AUTHOR]

The Random-Phase approximation (RPA) provides an appealing framework for semi-local density functional theory. In its Resolution-of-the-Identity (RI) approach, it is a very accurate and more cost-effective method than most other wavefunction-based correlation methods. For widespread applications, efficient implementations of nuclear gradients for structure optimizations and data sampling of machine learning approaches are required. We report a well scaling implementation of RI-RPA nuclear gradients on massively parallel computers. The approach is applied to two polymorphs of the benzene crystal obtaining very good cohesive and relative energies. Different correction and extrapolation schemes are investigated for further improvement of the results and estimations of error bars. [ABSTRACT FROM AUTHOR]

First principles molecular dynamics simulation protocol is established using revised functional of Perdew-Burke-Ernzerhof (revPBE) in conjunction with Grimme's third generation of dispersion (D3) correction to describe the properties of water at ambient conditions. This study also demonstrates the consistency of the structure of water across both isobaric (NpT) and isothermal (NVT) ensembles. Going beyond the standard structural benchmarks for liquid water, we compute properties that are connected to both local structure and mass density fluctuations that are related to concepts of solvation and hydrophobicity. We directly compare our revPBE results to the Becke-Lee-Yang-Parr (BLYP) plus Grimme dispersion corrections (D2) and both the empirical fixed charged model (SPC/E) and many body interaction potential model (MB-pol) to further our understanding of how the computed properties herein depend on the form of the interaction potential.

The development of novel double-hybrid density functionals offers new levels of accuracy and is leading to fresh insights into the fundamental properties of matter. Hartree–Fock exact exchange and correlated wave function methods, such as second-order Møller–Plesset (MP2) and direct random phase approximation (dRPA), are usually required to build such functionals. Their high computational cost is a concern, and their application to large and periodic systems is, therefore, limited. In this work, low-scaling methods for Hartree–Fock exchange (HFX), SOS-MP2, and direct RPA energy gradients are developed and implemented in the CP2K software package. The use of the resolution-of-the-identity approximation with a short range metric and atom-centered basis functions leads to sparsity, allowing for sparse tensor contractions to take place. These operations are efficiently performed with the newly developed Distributed Block-sparse Tensors (DBT) and Distributed Block-sparse Matrices (DBM) libraries, which scale to hundreds of graphics processing unit (GPU) nodes. The resulting methods, resolution-of-the-identity (RI)-HFX, SOS-MP2, and dRPA, were benchmarked on large supercomputers. They exhibit favorable sub-cubic scaling with system size, good strong scaling performance, and GPU acceleration up to a factor of 3. These developments will allow for double-hybrid level calculations of large and periodic condensed phase systems to take place on a more regular basis. [ABSTRACT FROM AUTHOR]

Accurate descriptions of intermolecular interactions are of great importance in simulations of molecular liquids. We present an electronic structure method that combines the accuracy of the Harris functional approach with the computational efficiency of approximately linear-scaling density functional theory (DFT). The non-variational nature of the Harris functional has been addressed by constructing a Lagrangian energy functional, which restores the variational condition by imposing stationarity with respect to the reference density. The associated linear response equations may be treated with linear-scaling efficiency in an atomic orbital based scheme. Key ingredients to describe the structural and dynamical properties of molecular systems are the forces acting on the atoms and the stress tensor. These first-order derivatives of the Harris Lagrangian have been derived and implemented in consistence with the energy correction. The proposed method allows for simulations with accuracies close to the Kohn–Sham DFT reference. Embedded in the CP2K program package, the method is designed to enable ab initio molecular dynamics simulations of molecular solutions for system sizes of several thousand atoms. Available subsystem DFT methods may be used to provide the reference density required for the energy correction at near linear-scaling efficiency. As an example of production applications, we applied the method to molecular dynamics simulations of the binary mixtures cyclohexane-methanol and toluene-methanol, performed within the isobaric-isothermal ensemble, to investigate the hydrogen bonding network in these non-ideal mixtures. [ABSTRACT FROM AUTHOR]

Due to their improved accuracy, double-hybrid density functionals emerged as an important method for molecular electronic-structure calculations. The high computational costs of double-hybrid calculations in the condensed phase and the lack of efficient gradient implementations thereof inhibit a wide applicability for periodic systems. We present an implementation of forces and stress tensors for double-hybrid density functionals within the Gaussian and plane-waves electronic structure framework. The auxiliary density matrix method is used to reduce the overhead of the Hartree–Fock kernel providing an efficient and accurate methodology to tackle condensed phase systems. First applications to water systems of different densities and molecular crystals show the efficiency of the implementation and pave the way for advanced studies. Finally, we present large benchmark systems to discuss the performance of our implementation on modern large-scale computers. [ABSTRACT FROM AUTHOR]

Linear-response time-dependent density functional theory (LR-TDDFT) for core level spectroscopy using standard local functionals suffers from self-interaction error and a lack of orbital relaxation upon creation of the core hole. As a result, LR-TDDFT calculated x-ray absorption near edge structure spectra needed to be shifted along the energy axis to match experimental data. We propose a correction scheme based on many-body perturbation theory to calculate the shift from first-principles. The ionization potential of the core donor state is first computed and then substituted for the corresponding Kohn–Sham orbital energy, thus emulating Koopmans's condition. Both self-interaction error and orbital relaxation are taken into account. The method exploits the localized nature of core states for efficiency and integrates seamlessly in our previous implementation of core level LR-TDDFT, yielding corrected spectra in a single calculation. We benchmark the correction scheme on molecules at the K- and L-edges as well as for core binding energies and report accuracies comparable to higher order methods. We also demonstrate applicability in large and extended systems and discuss efficient approximations. [ABSTRACT FROM AUTHOR]

Water adsorption and dissociation processes on pristine low-index TiO2 interfaces are important but poorly understood outside the well-studied anatase (101) and rutile (110). To understand these, we construct three sets of machine learning potentials that are simultaneously applicable to various TiO2 surfaces, based on three density-functional-theory approximations. Here we show the water dissociation free energies on seven pristine TiO2 surfaces, and predict that anatase (100), anatase (110), rutile (001), and rutile (011) favor water dissociation, anatase (101) and rutile (100) have mostly molecular adsorption, while the simulations of rutile (110) sensitively depend on the slab thickness and molecular adsorption is preferred with thick slabs. Moreover, using an automated algorithm, we reveal that these surfaces follow different types of atomistic mechanisms for proton transfer and water dissociation: one-step, two-step, or both. These mechanisms can be rationalized based on the arrangements of water molecules on the different surfaces. Our finding thus demonstrates that the different pristine TiO2 surfaces react with water in distinct ways, and cannot be represented using just the low-energy anatase (101) and rutile (110) surfaces. Electronic structure methods are vital, yet they are often too computationally expensive. Here, the authors develop machine learned density matrices to fully represent electronic structures in a computationally cheap and accurate way. [ABSTRACT FROM AUTHOR]

We study self-interaction effects in solvated and strongly-correlated cationic molecular clusters, with a focus on the solvated hydroxyl radical. To address the self-interaction issue, we apply the DC-r2SCAN method, with the auxiliary density matrix approach. Validating our method through simulations of bulk liquid water, we demonstrate that DC-r2SCAN maintains the structural accuracy of r2SCAN while effectively addressing spin density localization issues. Extending our analysis to solvated cationic molecular clusters, we find that the hemibonded motif in the [CH3S∴CH3SH]+ cluster is disrupted in the DC-r2SCAN simulation, in contrast to r2SCAN that preserves the (three-electron-two-center)-bonded motif. Similarly, for the [SH∴SH2]+ cluster, r2SCAN restores the hemibonded motif through spin leakage, while DC-r2SCAN predicts a weaker hemibond formation influenced by solvent–solute interactions. Our findings demonstrate the potential of DC-r2SCAN combined with the auxiliary density matrix method to improve electronic structure calculations, providing insights into the properties of solvated cationic molecular clusters. This work contributes to the advancement of self-interaction corrected electronic structure theory and offers a computational framework for modeling condensed phase systems with intricate correlation effects. [ABSTRACT FROM AUTHOR]

An integral scheme for the efficient evaluation of two-center integrals over contracted solid harmonic Gaussian functions is presented. Integral expressions are derived for local operators that depend on the position vector of one of the two Gaussian centers. These expressions are then used to derive the formula for three-index overlap integrals where two of the three Gaussians are located at the same center. The efficient evaluation of the latter is essential for local resolution-of-the-identity techniques that employ an overlap metric.We compare the performance of our integral scheme to the widely used Cartesian Gaussian-based method of Obara and Saika (OS). Non-local interaction potentials such as standard Coulomb, modified Coulomb, and Gaussian-type operators, which occur in range-separated hybrid functionals, are also included in the performance tests. The speed-up with respect to the OS scheme is up to three orders of magnitude for both integrals and their derivatives. In particular, our method is increasingly efficient for large angular momenta and highly contracted basis sets. [ABSTRACT FROM AUTHOR]

We have studied polarized Au(100) and Au(111) electrodes immersed in electrolyte solution by implementing finite‐field methods in density functional theory‐based molecular dynamics simulations. This allows us to directly compute the Helmholtz capacitance of electric double layer by including both electronic and ionic degrees of freedom, and the results turn out to be in excellent agreement with experiments. It is found that the electronic response of Au electrode makes a crucial contribution to the high Helmholtz capacitance and the instantaneous adsorption of Cl can lead to a charge inversion on the anodic polarized Au(100) surface. These findings point out ways to improve popular semi‐classical models for simulating electrified solid–liquid interfaces and to identify the nature of surface charges therein which are difficult to access in experiments. [ABSTRACT FROM AUTHOR]

The forces acting on the atoms as well as the stress tensor are crucial ingredients for calculating the structural and dynamical properties of systems in the condensed phase. Here, these derivatives of the total energy are evaluated for the second-order Møller-Plesset perturbation energy (MP2) in the framework of the resolution of identity Gaussian and plane waves method, in a way that is fully consistent with how the total energy is computed. This consistency is non-trivial, given the different ways employed to compute Coulomb, exchange, and canonical four center integrals, and allows, for example, for energy conserving dynamics in various ensembles. Based on this formalism, a massively parallel algorithm has been developed for finite and extended system. The designed parallel algorithm displays, with respect to the system size, cubic, quartic, and quintic requirements, respectively, for the memory, communication, and computation. All these requirements are reduced with an increasing number of processes, and the measured performance shows excellent parallel scalability and efficiency up to thousands of nodes. Additionally, the computationally more demanding quintic scaling steps can be accelerated by employing graphics processing units (GPU's) showing, for large systems, a gain of almost a factor two compared to the standard central processing unit-only case. In this way, the evaluation of the derivatives of the RI-MP2 energy can be performed within a few minutes for systems containing hundreds of atoms and thousands of basis functions. With good time to solution, the implementation thus opens the possibility to perform molecular dynamics (MD) simulations in various ensembles (microcanonical ensemble and isobaric-isothermal ensemble) at the MP2 level of theory. Geometry optimization, full cell relaxation, and energy conserving MD simulations have been performed for a variety of molecular crystals including NH3, CO2, formic acid, and benzene. [ABSTRACT FROM AUTHOR]

Water is a ubiquitous liquid that displays a wide range of anomalous properties and has a delicate structure that challenges experiment and simulation alike. The various intermolecular interactions that play an important role, such as repulsion, polarization, hydrogen bonding, and van der Waals interactions, are often difficult to reproduce faithfully in atomistic models. Here, electronic structure theories including all these interactions at equal footing, which requires the inclusion of non-local electron correlation, are used to describe structure and dynamics of bulk liquid water. Isobaricisothermal (NpT) ensemble simulations based on the Random Phase Approximation (RPA) yield excellent density (0.994 g/ml) and fair radial distribution functions, while various other density functional approximations produce scattered results (0.8-1.2 g/ml). Molecular dynamics simulation in the microcanonical (NVE) ensemble based on Møller-Plesset perturbation theory (MP2) yields dynamical properties in the condensed phase, namely, the infrared spectrum and diffusion constant. At the MP2 and RPA levels of theory, ice is correctly predicted to float on water, resolving one of the anomalies as resulting from a delicate balance between van der Waals and hydrogen bonding interactions. For several properties, obtaining quantitative agreement with experiment requires correction for nuclear quantum effects (NQEs), highlighting their importance, for structure, dynamics, and electronic properties. A computed NQE shift of 0.6 eV for the band gap and absorption spectrum illustrates the latter. Giving access to both structure and dynamics of condensed phase systems, non-local electron correlation will increasingly be used to study systems where weak interactions are of paramount importance. [ABSTRACT FROM AUTHOR]

We describe the calculation of Raman spectra for periodic systems via ab initio molecular dynamics (AIMD) utilizing the Gaussian and plane wave method in the program package CP2K. The electric -dipole-electric-dipole polarizability tensor has been implemented for an arbitrary shape of the simulation cell. In addition, a computationally efficient approach for its decomposition into local contributions is presented. As an example for the application of computational Raman spectroscopy to liquids, the Raman spectra of S-methyloxirane in the liquid phase have been calculated together with Raman spectra obtained from static calculations employing the double-harmonic approximation. The comparison to experimental data illustrates that a very good agreement between experiment and simulated spectra can be obtained employing AIMD, which takes into account anharmonicities and dynamical effects at ambient conditions. [ABSTRACT FROM AUTHOR]

We discuss our implementation of linear-response time-dependent density functional theory (LR-TDDFT) for core level near-edge absorption spectroscopy. The method is based on established LR-TDDFT approaches to X-ray absorption spectroscopy (XAS) with additional accurate approximations for increased efficiency. We validate our implementation by reproducing benchmark results at the K-edge and showing that spin–orbit coupling effects at the L2,3-edge are well described. We also demonstrate that the method is suitable for extended systems in periodic boundary conditions and measure a favorable sub-cubic scaling of the calculation cost with system size. We finally show that GPUs can be efficiently exploited and report speedups of up to a factor 2. [ABSTRACT FROM AUTHOR]

Recently, the nonlocal van der Waals (vdW) density functionals [M. Dion, H. Rydberg, E. Schröder, D. C. Langreth, and B. I. Lundqvist, Phys. Rev. Lett. 92, 246401 (2004)] have attracted considerable attention due to their good performance for systems where weak interactions are important. Since the physics of dispersion is included in these functionals, they are usually more accurate and show less erratic behavior than the semilocal and hybrid methods. In this work, several variants of the vdW functionals have been tested on rare-gas dimers (from He2 to Kr2) and solids (Ne, Ar, and Kr) and their accuracy compared to standard semilocal approximations, supplemented or not by an atom-pairwise dispersion correction [S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, J. Chem. Phys. 132, 154104 (2010)]. An analysis of the results in terms of energy decomposition is also provided. [ABSTRACT FROM AUTHOR]

The ability to reproduce the experimental structure of water around the sodium and potassium ions is a key test of the quality of interaction potentials due to the central importance of these ions in a wide range of important phenomena. Here, we simulate the Na+ and K+ ions in bulk water using three density functional theory functionals: (1) the generalized gradient approximation (GGA) based dispersion corrected revised Perdew, Burke, and Ernzerhof functional (revPBE-D3) (2) the recently developed strongly constrained and appropriately normed (SCAN) functional (3) the random phase approximation (RPA) functional for potassium. We compare with experimental X-ray diffraction (XRD) and X-ray absorption fine structure (EXAFS) measurements to demonstrate that SCAN accurately reproduces key structural details of the hydration structure around the sodium and potassium cations, whereas revPBE-D3 fails to do so. However, we show that SCAN provides a worse description of pure water in comparison with revPBE-D3. RPA also shows an improvement for K+, but slow convergence prevents rigorous comparison. Finally, we analyse cluster energetics to show SCAN and RPA have smaller fluctuations of the mean error of ion–water cluster binding energies compared with revPBE-D3. [ABSTRACT FROM AUTHOR]

Structural information is extracted from the all-particle (non-Born-Oppenheimer) wave function by calculating radial and angular densities derived from n-particle densities. As a result, one- and two-dimensional motifs of classical molecular structure can be recognized in quantum mechanics. Numerical examples are presented for three- (H-, Ps-, H2+), four- (Ps2, H2), and five-particle (H2D+) systems. [ABSTRACT FROM AUTHOR]

We introduce a method for the all-electron calculation of the NMR chemical shifts and the EPR g tensor using the Gaussian and augmented-plane-wave method. The presented approach is based on the generalized density functional perturbation theory. The method is validated by comparison with other theoretical methods for a selection of small molecules. We also present two exemplary applications that involve the calculation of the chemical shifts of a hydrated adenine and the g tensor for the E1′ center in α-quartz using a quantum mechanical/molecular mechanical approach. [ABSTRACT FROM AUTHOR]

In an article [A. Aguado and P. A. Madden, J. Chem. Phys. 119, 7471 (2003)] published in this journal, Ewald summation expressions were derived for the energy, interatomic forces, pressure tensor, electrostatic field, and electrostatic field gradients in simulation system composed of molecules with charges, induced dipoles, and quadrupoles. In this letter we propose an alternative formulation of the reciprocal space terms generalized to higher multipoles, providing, at the same time, a few important corrections for previously published derivations. The present expressions, more compact than the ones proposed in the original work, provide a straightforward approach to implement an Ewald summation scheme for multipole interactions in codes where a standard Ewald summation is already available. A major result of the present derivation is an increase in the computational efficiency compared to the previous implementation of the several different electrostatic terms. [ABSTRACT FROM AUTHOR]

Ab initio molecular dynamics simulations with hybrid density functionals have so far found little application due to their computational cost. In this work, an implementation of the Hartree–Fock exchange is presented that is specifically targeted at ab initio molecular dynamics simulations of medium sized systems. We demonstrate that our implementation, which is available as part of the CP2K/Quickstep program, is robust and efficient. Several prescreening techniques lead to a linear scaling cost for integral evaluation and storage. Integral compression techniques allow for in-core calculations on systems containing several thousand basis functions. The massively parallel implementation respects integral symmetry and scales up to hundreds of CPUs using a dynamic load balancing scheme. A time-reversible multiple time step scheme, exploiting the difference in computational efficiency between hybrid and local functionals, brings further time savings. With extensive simulations of liquid water, we demonstrate the ability to perform, for several tens of picoseconds, ab initio molecular dynamics based on hybrid functionals of systems in the condensed phase containing a few thousand Gaussian basis functions. [ABSTRACT FROM AUTHOR]

The direct energy functional minimization problem in electronic structure theory, where the single-particle orbitals are optimized under the constraint of orthogonality, is explored. We present an orbital transformation based on an efficient expansion of the inverse factorization of the overlap matrix that keeps orbitals orthonormal. The orbital transformation maps the orthogonality constrained energy functional to an approximate unconstrained functional, which is correct to some order in a neighborhood of an orthogonal but approximate solution. A conjugate gradient scheme can then be used to find the ground state orbitals from the minimization of a sequence of transformed unconstrained electronic energy functionals. The technique provides an efficient, robust, and numerically stable approach to direct total energy minimization in first principles electronic structure theory based on tight-binding, Hartree–Fock, or density functional theory. For sparse problems, where both the orbitals and the effective single-particle Hamiltonians have sparse matrix representations, the effort scales linearly with the number of basis functions N in each iteration. For problems where only the overlap and Hamiltonian matrices are sparse the computational cost scales as O(M2N), where M is the number of occupied orbitals. We report a single point density functional energy calculation of a DNA decamer hydrated with 4003 water molecules under periodic boundary conditions. The DNA fragment containing a cis-syn thymine dimer is composed of 634 atoms and the whole system contains a total of 12 661 atoms and 103 333 spherical Gaussian basis functions. [ABSTRACT FROM AUTHOR]

A smooth ℓ1-norm based function to obtain a sparse representation of the orbital coefficients is introduced. This sparseness function is further parametrized with respect to unitary transformations among the occupied orbitals. Thus the function can be straightforwardly included in an optimization scheme or used on the fly during self-consistent field iterations to induce or maintain the sparsity of the orbital coefficients. As practical examples, we induce sparsity in the orbital coefficients of liquid water and bulk silicon. We also report the sparsity of the orbital coefficients of 1024 water molecules along a short Born–Oppenheimer molecular dynamics trajectory. It is observed that, after a stabilization period, the sparsity of the orbitals can be kept stable along the dynamics with small additional computational effort. [ABSTRACT FROM AUTHOR]

We present a library of Gaussian basis sets that has been specifically optimized to perform accurate molecular calculations based on density functional theory. It targets a wide range of chemical environments, including the gas phase, interfaces, and the condensed phase. These generally contracted basis sets, which include diffuse primitives, are obtained minimizing a linear combination of the total energy and the condition number of the overlap matrix for a set of molecules with respect to the exponents and contraction coefficients of the full basis. Typically, for a given accuracy in the total energy, significantly fewer basis functions are needed in this scheme than in the usual split valence scheme, leading to a speedup for systems where the computational cost is dominated by diagonalization. More importantly, binding energies of hydrogen bonded complexes are of similar quality as the ones obtained with augmented basis sets, i.e., have a small (down to 0.2 kcal/mol) basis set superposition error, and the monomers have dipoles within 0.1 D of the basis set limit. However, contrary to typical augmented basis sets, there are no near linear dependencies in the basis, so that the overlap matrix is always well conditioned, also, in the condensed phase. The basis can therefore be used in first principles molecular dynamics simulations and is well suited for linear scaling calculations. [ABSTRACT FROM AUTHOR]

A combined linear-response–frozen electron-density model has been implemented in a molecular-dynamics scheme derived from an extended Lagrangian formalism. This approach is based on a partition of the electronic charge distribution into a frozen region described by Kim-Gordon theory [J. Chem. Phys. 56, 3122 (1972); J. Chem. Phys. 60, 1842 (1974)] and a response contribution determined by the instantaneous ionic configuration of the system. The method is free from empirical pair potentials and the parametrization protocol involves only calculations on properly chosen subsystems. We apply this method to a series of alkali halides in different physical phases and are able to reproduce experimental structural and thermodynamic properties with an accuracy comparable to Kohn-Sham density-functional calculations. [ABSTRACT FROM AUTHOR]

The performance of density functional theory methods for the modeling of condensed aqueous systems is hard to predict and validation by ab initio molecular simulation of liquid water is absolutely necessary. In order to assess the reliability of these tests, the effect of temperature on the structure and dynamics of liquid water has been characterized with 16 simulations of 20 ps in the temperature range of 280–380 K. We find a pronounced influence of temperature on the pair correlation functions and on the diffusion constant including nonergodic behavior on the time scale of the simulation in the lower temperature range (which includes ambient temperature). These observations were taken into account in a consistent comparison of a series of density functionals (BLYP, PBE, TPSS, OLYP, HCTH120, HCTH407). All simulations were carried out using an ab initio molecular dynamics approach in which wave functions are represented using Gaussians and the density is expanded in an auxiliary basis of plane waves. Whereas the first three functionals show similar behavior, it is found that the latter three functionals yield more diffusive dynamics and less structure.© 2005 American Institute of Physics. [ABSTRACT FROM AUTHOR]

A recently introduced formulation of time dependent linear response density functional theory within the plane-wave pseudopotential framework [J. Hutter, J. Chem. Phys. 118, 3928 (2003)] is applied to the study of solvent shift and intensity enhancement effects of the [sup 1]A[sub 2] n→π[sup *] electronic transition in acetone, treating solute and solvent at the same level of theory. We propose a suitable formalism for computing transition intensities based on the modern theory of polarization, which is applicable to condensed-phase and finite systems alike. The gain in intensity brought about by thermal fluctuations is studied in molecular acetone at room temperature, and in gas-phase (CH[sub 3])[sub 2]CO·(H[sub 2]O)[sub 2] at 25 K. The latter system is characterized by the appearance of relatively intense features in the low-energy region of the spectrum, attributable to spurious solvent→solute charge-transfer excitations created by deficiencies in the DFT methodology. The n→π[sup *] transition can be partially isolated from the charge-transfer bands, yielding a blueshift of 0.17 eV with respect to gas-phase acetone. This analysis is then carried over to a solution of acetone in water, where further complications are encountered in the from of a solute→solvent charge transfer excitations overlapping with the n→π[sup *] band. The optically active occupied states are found to be largely localized on either solute or solvent, and using this feature we were again able to isolate the physical n→π[sup *] band and compute the solvatochromic shift. The result of 0.19 eV is in good agreement with experiment, as is the general increase in the mean oscillator strength of the transition. The unphysical charge transfers are interpreted in terms of degeneracies in the spectrum of orbital energies of the aqueous acetone solution.© 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

An efficient method for optimizing single-determinant wave functions of medium and large systems is presented. It is based on a minimization of the energy functional using a new set of variables to perform orbital transformations. With this method convergence of the wave function is guaranteed. Preconditioners with different computational cost and efficiency have been constructed. Depending on the preconditioner, the method needs a number of iterations that is very similar to the established diagonalization-DIIS approach, in cases where the latter converges well. Diagonalization of the Kohn-Sham matrix can be avoided and the sparsity of the overlap and Kohn-Sham matrix can be exploited. If sparsity is taken into account, the method scales as O(MN[SUP2]), where M is the total number of basis functions and N is the number of occupied orbitals. The relative performance of the method is optimal for large systems that are described with high quality basis sets, and for which the density matrices are not yet sparse. We present a benchmark calculation on a DNA crystal containing 2312 base pairs, solvent and counter ions (2388 atoms), using a TZV(2d,2p) basis (38 688 basis functions) and conclude that the electronic structure of systems of this size can now be studied routinely. [ABSTRACT FROM AUTHOR]

An efficient formulation of time-dependent linear response density functional theory for the use within the plane wave basis set framework is presented. The method avoids the transformation of the Kohn–Sham matrix into the canonical basis and references virtual orbitals only through a projection operator. Using a Lagrangian formulation nuclear derivatives of excited state energies within the Tamm–Dancoff approximation are derived. The algorithms were implemented into a pseudo potential/plane wave code and applied to the calculation of adiabatic excitation energies, optimized geometries and vibrational frequencies of three low lying states of formaldehyde. An overall good agreement with other time-dependent density functional calculations, multireference configuration interaction calculations and experimental data was found. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Three frequently used gradient-corrected density functionals (B, BP, and BLYP) are applied in an ab initio molecular dynamics simulation of liquid water in order to evaluate their performance for the description of condensed aqueous systems. A comparison of structural characteristics (radial distribution functions) and dynamical properties (vibrational spectra, orientational relaxation, and self-diffusion) leads to the conclusion that hydrogen bonding is too weak in the usual local density approximation corrected for exchange only according to Becke (B), whereas adding the gradient correction for correlation according to Perdew (BP) yields effective hydrogen bonds in the liquid that are too strong. The combination of B with the semilocal correlation functional according to Lee, Yang, and Parr (BLYP) yields the best agreement with experiment. The computational method, which is the basis for the determination of (adiabatic) electronic structure in the ab initio molecular dynamics simulation, has been validated by an extensive series of test calculations for the water dimer, which will also be presented here. © 1996 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The ground- and first excited state of s-tetrazine arising from a π*←n excitation (1Ag,1B3u) have been studied using the complete active space (CASSCF) and the second order multiconfiguration perturbation theory (CASPT2) ab initio methods. The focus of this study is on the effect of the electronic excitation on the molecular structure and on those electronic properties which are important to model the solvatochromatic behavior of the molecule in polymer matrices as used in permanent hole burning experiments. Since the accurate computation of excited state molecular properties represents a major challenge for today’s numerical quantum chemistry, some technical aspects are also considered. The present study shows that the change in geometry upon electronic excitation is small. This is in partial contradiction with the experimental studies which however disagree among themselves [see K. K. Innes, I. G. Ross, and W. R. Moomaw, J. Mol. Spectrosc. 132, 492 (1988), and R. E. Smalley, L. Wharton, and D. H. Levi, ibid. 66, 375 (1977)]. This study also confirms that the first excited state equilibrium structure is of D2h symmetry. In an earlier theoretical study it was found that the D2h symmetry structure may represent a saddle point rather than a minimum on the excited state potential surface [see A. C. Scheiner and H. F. Schaefer III, J. Chem. Phys. 87, 3539 (1987)]. In the first excited state, we observe an increase of the mean polarizability of s-tetrazine along with an enhanced anisotropy. The change in the polarizability is almost exclusively in the ‘‘in-plane’’ components of the tensor; the polarizability in the vertical direction is nearly unchanged. This observation questions recent experimental results reported for this molecule [see S. Heitz, D. Weidnauer, and A. Hese, J. Chem. Phys. 95, 7952 (1991)]. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

New methods for integrating the Car–Parrinello equations with ultrasoft pseudopotentials are introduced. In particular, the difficulties associated with the generalized orthonormality constraint condition 〈[lowercase_phi_synonym]i|S({RI})|[lowercase_phi_synonym]j=δij are addressed. It is shown that the equations of motion can be integrated using the velocity Verlet/RATTLE scheme, and a new method, the constrained nonorthogonal orbital method, that eliminates the need to enforce this constraint explicitly is introduced. In this new scheme, the generalized orthonormality constraint is satisfied implicitly, thus allowing the freedom to choose simpler constraint conditions. We show that usual N3 scaling associated with the calculation of the Lagrange multipliers in the constraint force can be reduced to an N2 calculation by the use of a simple set of norm or length constraints on the electronic orbitals without sacrificing accuracy. The constrained nonorthogonal orbital method is shown to be considerably simpler to implement and more efficient than the standard approach to the ultrasoft pseudopotential problem. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In this paper we present new formulas related to the parametrization of molecular orbitals by exponential transformations. We show that it is not necessary to use series expansions of the exponentials and that the first derivative of the energy with respect to independent parameters can be calculated at arbitrary positions in an efficient manner. These formulas can be implemented such that computational costs are of order N2M and storage requirements of order NM, where N is the number of occupied orbitals and M the size of the basis set. We discuss possible applications to Hartree–Fock calculations and ab initio molecular dynamics. [ABSTRACT FROM AUTHOR]

The question of the lowest energy structure of the C6 molecule has been addressed using high level of theory ab initio methods such as the coupled cluster and multiconfiguration perturbation theory approaches. The results show that the linear cumulenic chain and the monocyclic six-membered ring structure with D3h symmetry are very close in energy, indicating that both structures may be observed experimentally. The hexagonal six-membered ring structure of D6h symmetry was confirmed as a saddle point geometry. These results are in partial contradiction with two earlier ab initio studies that either find the distorted hexagon structure [Raghavachari, Whiteside, and Pople, J. Chem. Phys. 85, 6623 (1986)] or the linear cumulenic chain [V. Parasuk and J. Almlöf, J. Chem. Phys. 91, 1137 (1989)] to be the lowest energy structure. [ABSTRACT FROM AUTHOR]

Understanding lithium solvation and transport in ionic liquids is important due to their possible application in electrochemical devices. Using first-principles simulations aided by a metadynamics approach we study the free-energy landscape for lithium ions at infinite dilution in ethylammonium nitrate, a protic ionic liquid. We analyze the local structure of the liquid around the lithium cation and obtain a quantitative picture in agreement with experimental findings. Our simulations show that the lowest two free energy minima correspond to conformations with the lithium ion being solvated either by three or four nitrate ions with a transition barrier between them of 0.2 eV. Other less probable conformations having different solvation pattern are also investigated. [ABSTRACT FROM AUTHOR]

A novel strategy for post-synthesis amine borane functionalization of MOFs under gas-solid phase transformation, utilizing gaseous diborane, is reported. The covalently confined amine borane derivative decorated on the framework backbone is stable when preserved at low temperature, but spontaneously liberates soft chemical hydrogen at room temperature, leading to the development of an unusual borenium type species (-NH=BH2+) ion-paired with a hydroborate anion. Furthermore, the unsaturated amino borane (-NH=BH2) and the μ-iminodiborane (-μ-NHB2H5) were detected as final products. A combination of DFT based molecular dynamics simulations and solid state NMR spectroscopy, utilizing isotopically enriched materials, were undertaken to unequivocally elucidate the mechanistic pathways for H2 liberation. [ABSTRACT FROM AUTHOR]

Solar-light driven water splitting is a promising way for the sustainable production of molecular hydrogen, the latter representing an efficient carrier for energy storage and conversion into common liquid fuels. In search of novel catalysts for high-performance water splitting devices, Co-pyrphyrin (CoPy) has been recently synthesized and successfully used as a homogeneous water reduction catalyst. We investigate the adsorption of this molecule on the rutile TiO2(110) surface as a possible first step towards the design of a heterogeneous water reduction system. We find that the adsorption of the molecule is stabilized by the interaction of the cyano groups with the under-coordinated Ti centers present at the surface. This interaction induces the rehybridization of the molecular orbitals localized on the cyano groups and the realignment of the lowest unoccupied molecular states. Moreover, the highest occupied molecular orbital of CoPy@rutile(110) is localized on CoPy and the energy gap turns out to be significantly smaller than the gap of pristine rutile(110). This implies that direct or indirect injection of electrons from CoPy to the rutile(110) surface is in principle possible upon the absorption of light in the visible range. On the other hand, the electronic properties of the Co(ii) center are not modified by the adsorption, which suggests that CoPy and its derivatives may be used in water electrolysis for hydrogen production also in the adsorbed state. [ABSTRACT FROM AUTHOR]