Harmonic Infrared and Raman Spectra in Molecular Environments Using the Polarizable Embedding Model.

We present a fully analytic approach to calculate infrared (IR) and Raman spectra of molecules embedded in complex molecular environments modeled using the fragment-based polarizable embedding (PE) model. We provide the theory for the calculation of analytic second-order geometric derivatives of molecular energies and first-order geometric derivatives of electric dipole moments and dipole-dipole polarizabilities within the PE model. The derivatives are implemented using a general open-ended response theory framework, thus allowing for an extension to higher-order derivatives.

Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems.

The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized.

On Resolution-of-the-Identity Electron Repulsion Integral Approximations and Variational Stability.

The definiteness of the Mulliken and Dirac electron repulsion integral (ERI) matrices is examined for different classes of resolution-of-the-identity (RI) ERI approximations with particular focus on local fitting techniques. For global RI, robust local RI, and nonrobust local RI we discuss the definiteness of the approximated ERI matrices as well as the resulting bounds of Hartree, exchange, and total energies. Lower bounds of Hartree and exchange energy contributions are crucial as their absence may lead to variational instabilities, causing severe convergence problems or even convergence to a spurious state in self-consistent-field optimizations.

Explicitly correlated second-order Møller-Plesset perturbation theory in a Divide-Expand-Consolidate (DEC) context.

We present the DEC-RIMP2-F12 method where we have augmented the Divide Expand-Consolidate resolution-of-the-identity second-order Møller-Plesset perturbation theory method (DEC-RIMP2) [P. Baudin et al., J.

Comparison of Three Efficient Approximate Exact-Exchange Algorithms: The Chain-of-Spheres Algorithm, Pair-Atomic Resolution-of-the-Identity Method, and Auxiliary Density Matrix Method.

We compare the performance of three approximate methods for speeding up evaluation of the exchange contribution in Hartree-Fock and hybrid Kohn-Sham calculations: the chain-of-spheres algorithm (COSX; Neese , F. Chem. Phys.

Efficient linear-scaling second-order Møller-Plesset perturbation theory: The divide-expand-consolidate RI-MP2 model.

The Resolution of the Identity second-order Møller-Plesset perturbation theory (RI-MP2) method is implemented within the linear-scaling Divide-Expand-Consolidate (DEC) framework. In a DEC calculation, the full molecular correlated calculation is replaced by a set of independent fragment calculations each using a subset of the total orbital space. The number of independent fragment calculations scales linearly with the system size, rendering the method linear-scaling and massively parallel.

Wavelet formulation of the polarizable continuum model. II. Use of piecewise bilinear boundary elements.

The simplicity of dielectric continuum models has made them a standard tool in almost any Quantum Chemistry (QC) package. Despite being intuitive from a physical point of view, the actual electrostatic problem at the cavity boundary is challenging: the underlying boundary integral equations depend on singular, long-range operators. The parametrization of the cavity boundary should be molecular-shaped, smooth and differentiable.

Ring planarity problem of 2-oxazoline revisited using microwave spectroscopy and quantum chemical calculations.

In a previous infrared, Raman, and microwave spectroscopic work,1 it was claimed that 2-oxazoline has a planar ring equilibrium conformation, and the ring-puckering potential function V(z) = 22.2(z(4) + 1.31z(2)) cm(-1), where z is a dimensionless reduced coordinate, was derived.

Geometry of the magic number H(+)(H2O)21 water cluster by proxy.

Abundance mass spectra, obtained upon carefully electrospraying solutions of tert-butanol (TB) in water into a mass spectrometer, display a systematic series of peaks due to mixed H(+)(TB)m(H2O)n clusters. Clusters with m + n = 21 have higher abundance (magic number peaks) than their neighbours when m ≤ 9, while for m > 9 they have lower abundance. This indicates that the mixed TB-H2O clusters retain a core hydrogen bonded network analogous to that in the famous all-water H(+)(H2O)21 cluster up to the limit m = 9.

Charge-constrained auxiliary-density-matrix methods for the Hartree-Fock exchange contribution.

Three new variants of the auxiliary-density-matrix method (ADMM) of Guidon, Hutter, and VandeVondele [J. Chem. Theory Comput.

Attractive electron-electron interactions within robust local fitting approximations.

An analysis of Dunlap's robust fitting approach reveals that the resulting two-electron integral matrix is not manifestly positive semidefinite when local fitting domains or non-Coulomb fitting metrics are used. We present a highly local approximate method for evaluating four-center two-electron integrals based on the resolution-of-the-identity (RI) approximation and apply it to the construction of the Coulomb and exchange contributions to the Fock matrix. In this pair-atomic resolution-of-the-identity (PARI) approach, atomic-orbital (AO) products are expanded in auxiliary functions centered on the two atoms associated with each product.

Insights into the dynamics of evaporation and proton migration in protonated water clusters from large-scale Born-Oppenheimer direct dynamics.

Large-scale on-the-fly Born-Oppenheimer molecular dynamics simulations using recent advances in linear scaling electronic structure theory and trajectory integration techniques have been performed for protonated water clusters around the magic number (H(2)O)(n)H(+) , for n = 20 and 21. Besides demonstrating the feasibility and efficiency of the computational approach, the calculations reveal interesting dynamical details. Elimination of water molecules is found to be fast for both cluster sizes but rather insensitive to the initial geometry.

MP2 energy and density for large molecular systems with internal error control using the Divide-Expand-Consolidate scheme.

Divide-Expand-Consolidate (DEC) is a local correlation method where the inherent locality of the electron correlation problem is used to express the correlated calculation on a large molecular system in terms of small independent fragment calculations employing small subsets of local HF orbitals. A crucial feature of the DEC scheme is that the sizes of the local orbital spaces are determined in a black box manner during the calculation. In this way it is ensured that the correlation energy has been determined to a predefined precision compared to a conventional calculation.

Molecular gradient for second-order Møller-Plesset perturbation theory using the divide-expand-consolidate (DEC) scheme.

We demonstrate that the divide-expand-consolidate (DEC) scheme--which has previously been used to determine the second-order Møller-Plesset (MP2) correlation energy--can be applied to evaluate the MP2 molecular gradient in a linear-scaling and embarrassingly parallel manner using a set of local Hartree-Fock orbitals. All manipulations of four-index quantities (describing electron correlation effects) are carried out using small local orbital fragment spaces, whereas two-index quantities are treated for the full molecular system. The sizes of the orbital fragment spaces are determined in a black-box manner to ensure that the error in the DEC-MP2 correlation energy compared to a standard MP2 calculation is proportional to a single input threshold denoted the fragment optimization threshold (FOT).

Analytical GIAO and hybrid-basis integral derivatives: application to geometry optimization of molecules in strong magnetic fields.

Analytical integral evaluation is a central task of modern quantum chemistry. Here we present a general method for evaluating differentiated integrals over standard Gaussian and mixed Gaussian/plane-wave hybrid orbitals. The main idea is to have a representation of basis sets that is flexible enough to enable differentiated integrals to be reinterpreted as standard integrals over modified basis functions.

An efficient density-functional-theory force evaluation for large molecular systems.

An efficient, linear-scaling implementation of Kohn-Sham density-functional theory for the calculation of molecular forces for systems containing hundreds of atoms is presented. The density-fitted Coulomb force contribution is calculated in linear time by combining atomic integral screening with the continuous fast multipole method. For higher efficiency and greater simplicity, the near-field Coulomb force contribution is calculated by expanding the solid-harmonic Gaussian basis functions in Hermite rather than Cartesian Gaussians.

Variational and robust density fitting of four-center two-electron integrals in local metrics.

Density fitting is an important method for speeding up quantum-chemical calculations. Linear-scaling developments in Hartree-Fock and density-functional theories have highlighted the need for linear-scaling density-fitting schemes. In this paper, we present a robust variational density-fitting scheme that allows for solving the fitting equations in local metrics instead of the traditional Coulomb metric, as required for linear scaling.

A ground-state-directed optimization scheme for the Kohn-Sham energy.

Kohn-Sham density-functional calculations are used in many branches of science to obtain information about the electronic structure of molecular systems and materials. Unfortunately, the traditional method for optimizing the Kohn-Sham energy suffers from fundamental problems that may lead to divergence or, even worse, convergence to an energy saddle point rather than to the ground-state minimum--in particular, for the larger and more complicated electronic systems that are often studied by Kohn-Sham theory nowadays. We here present a novel method for Kohn-Sham energy minimization that does not suffer from the flaws of the conventional approach, combining reliability and efficiency with linear complexity.

A unified scheme for the calculation of differentiated and undifferentiated molecular integrals over solid-harmonic Gaussians.

Utilizing the fact that solid-harmonic combinations of Cartesian and Hermite Gaussian atomic orbitals are identical, a new scheme for the evaluation of molecular integrals over solid-harmonic atomic orbitals is presented, where the integration is carried out over Hermite rather than Cartesian atomic orbitals. Since Hermite Gaussians are defined as derivatives of spherical Gaussians, the corresponding molecular integrals become the derivatives of integrals over spherical Gaussians, whose transformation to the solid-harmonic basis is performed in the same manner as for integrals over Cartesian Gaussians, using the same expansion coefficients. The presented solid-harmonic Hermite scheme simplifies the evaluation of derivative molecular integrals, since differentiation by nuclear coordinates merely increments the Hermite quantum numbers, thereby providing a unified scheme for undifferentiated and differentiated four-center molecular integrals.

Linear-scaling implementation of molecular response theory in self-consistent field electronic-structure theory.

A linear-scaling implementation of Hartree-Fock and Kohn-Sham self-consistent field theories for the calculation of frequency-dependent molecular response properties and excitation energies is presented, based on a nonredundant exponential parametrization of the one-electron density matrix in the atomic-orbital basis, avoiding the use of canonical orbitals. The response equations are solved iteratively, by an atomic-orbital subspace method equivalent to that of molecular-orbital theory. Important features of the subspace method are the use of paired trial vectors (to preserve the algebraic structure of the response equations), a nondiagonal preconditioner (for rapid convergence), and the generation of good initial guesses (for robust solution).

Linear-scaling implementation of molecular electronic self-consistent field theory.

A linear-scaling implementation of Hartree-Fock and Kohn-Sham self-consistent field (SCF) theories is presented and illustrated with applications to molecules consisting of more than 1000 atoms. The diagonalization bottleneck of traditional SCF methods is avoided by carrying out a minimization of the Roothaan-Hall (RH) energy function and solving the Newton equations using the preconditioned conjugate-gradient (PCG) method. For rapid PCG convergence, the Lowdin orthogonal atomic orbital basis is used.