Search

Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples. [ABSTRACT FROM AUTHOR]

In Paper I, the performances of pre-screening (PS), extended PS (EPS), and cumulant (CU) approximations to the fourth-order density matrix were examined in the context of second-order N-electron valence state perturbation theory (NEVPT2). It has been found that the CU, PS, and even EPS approximations with loose thresholds may introduce intruder states. In the present work, the origin of these "false intruder" states introduced by approximated density matrices is discussed. Canonical NEVPT2 implementations employ a rank reduction trick. By analyzing its residual error, we find that the omission of the rank reduction leads to a more stable multireference perturbation theory for incomplete active space reference wave functions. Such a full rank (FR)-NEVPT2 formulation is equivalent to the conventional NEVPT2 method for the complete active space self-consistent field/complete active space configuration interaction reference wave function. A major drawback of the FR-NEVPT2 formulation is the necessity of the fifth-order density matrix. To avoid the construction of the high-order density matrices, the combination of the FR-NEVPT2 with the CU approximation is studied. However, we find that the CU approximation remains problematic as it still introduces intruder states. The question of how to robustly and efficiently perform internally contracted multireference perturbation theories with approximate densities remains a challenging field of investigation. [ABSTRACT FROM AUTHOR]

This paper reports the release of PathSum, a new software suite of state-of-the-art path integral methods for studying the dynamics of single or extended systems coupled to harmonic environments. The package includes two modules, suitable for system–bath problems and extended systems comprising many coupled system–bath units, and is offered in C++ and Fortran implementations. The system–bath module offers the recently developed small matrix path integral (SMatPI) and the well-established iterative quasi-adiabatic propagator path integral (i-QuAPI) method for iteration of the reduced density matrix of the system. In the SMatPI module, the dynamics within the entanglement interval can be computed using QuAPI, the blip sum, time evolving matrix product operators, or the quantum–classical path integral method. These methods have distinct convergence characteristics and their combination allows a user to access a variety of regimes. The extended system module provides the user with two algorithms of the modular path integral method, applicable to quantum spin chains or excitonic molecular aggregates. An overview of the methods and code structure is provided, along with guidance on method selection and representative examples. [ABSTRACT FROM AUTHOR]

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. This paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios. [ABSTRACT FROM AUTHOR]

In this paper, the Hessian matrix of valence bond (VB) self-consistent field (VBSCF) energy with respect to orbitals are evaluated by applying the nonorthogonal orbital based N-body reduced density matrices, which was presented in Paper I. To this end, an automatic formula/code generator (AFCG) is developed; with which the matrix elements between internally contracted excited configurations of VB wave function and the corresponding codes are generated automatically. Compared to the tedious manual formula deducing and implementing, AFCG is much more convenient and efficient, and enables us to avoid troublesome debugging. With the help of AFCG, the Hessian-based Newton-Raphson algorithm is implemented for the VBSCF orbital optimization. Test calculations indicate that the Newton-Raphson algorithm converges quadratically and has much better convergence behavior than the gradient-based LBFGS algorithms. Furthermore, a combined approach with LBFGS and Newton-Raphson algorithms is applied to reduce the total CPU time of the calculation. [ABSTRACT FROM AUTHOR]

Fewest-switches surface hopping (FSSH) has emerged as one of the leading methods for modeling the quantum dynamics of molecular systems. While its original formulation was limited to adiabatic populations, the growing interest in the application of FSSH to coherent phenomena prompts the question of how one should construct a complete density matrix based on FSSH trajectories. A straightforward solution is to define adiabatic coherences based on wavefunction coefficients. In this paper, we demonstrate that inconsistencies introduced in the density matrix through such treatment may lead to a violation of positivity. We furthermore show that a recently proposed coherent generalization of FSSH results in density matrices that satisfy positivity while yielding improved accuracy throughout much (but not all) of parameter space. [ABSTRACT FROM AUTHOR]

The generalized quantum master equation (GQME) provides a general and formally exact framework for simulating the reduced dynamics of open quantum systems. The recently introduced modified approach to the GQME (M-GQME) corresponds to a specific implementation of the GQME that is geared toward simulating the dynamics of the electronic reduced density matrix in systems governed by an excitonic Hamiltonian. Such a Hamiltonian, which is often used for describing energy and charge transfer dynamics in complex molecular systems, is given in terms of diabatic electronic states that are coupled to each other and correspond to different nuclear Hamiltonians. Within the M-GQME approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic density matrix is fully captured by a memory kernel superoperator, which can be obtained from short-lived (compared to the time scale of energy/charge transfer) projection-free inputs. In this paper, we test the ability of the M-GQME to predict the energy transfer dynamics within a seven-state benchmark model of the Fenna–Matthews–Olson (FMO) complex, with the short-lived projection-free inputs obtained via the Ehrenfest method. The M-GQME with Ehrenfest-based inputs is shown to yield accurate results across a wide parameter range. It is also found to dramatically outperform the direct application of the Ehrenfest method and to provide better-behaved convergence with respect to memory time in comparison to an alternative implementation of the GQME approach previously applied to the same FMO model. [ABSTRACT FROM AUTHOR]

To address the challenge of performance portability and facilitate the implementation of electronic structure solvers, we developed the basic matrix library (BML) and Parallel, Rapid O(N), and Graph-based Recursive Electronic Structure Solver (PROGRESS) library. The BML implements linear algebra operations necessary for electronic structure kernels using a unified user interface for various matrix formats (dense and sparse) and architectures (CPUs and GPUs). Focusing on density functional theory and tight-binding models, PROGRESS implements several solvers for computing the single-particle density matrix and relies on BML. In this paper, we describe the general strategies used for these implementations on various computer architectures, using OpenMP target functionalities on GPUs, in conjunction with third-party libraries to handle performance critical numerical kernels. We demonstrate the portability of this approach and its performance in benchmark problems. [ABSTRACT FROM AUTHOR]

Time-Resolved Serial Femtosecond Crystallography (TR-SFX) conducted at X-ray Free Electron Lasers (XFELs) has become a powerful tool for capturing macromolecular structural movies of light-initiated processes. As the capabilities of XFELs advance, we anticipate that a new range of coherent control and structural Raman measurements will become achievable. Shorter optical and x-ray pulse durations and increasingly more exotic pulse regimes are becoming available at free electron lasers. Moreover, with high repetition enabled by the superconducting technology of European XFEL (EuXFEL) and Linac Coherent Light Source (LCLS-II) , it will be possible to improve the signal-to-noise ratio of the light-induced differences, allowing for the observation of vibronic motion on the sub-Angstrom level. To predict and assign this coherent motion, which is measurable with a structural technique, new theoretical approaches must be developed. In this paper, we present a theoretical density matrix approach to model the various population and coherent dynamics of a system, which considers molecular system parameters and excitation conditions. We emphasize the use of the Wigner transform of the time-dependent density matrix, which provides a phase space representation that can be directly compared to the experimental positional displacements measured in a TR-SFX experiment. Here, we extend the results from simple models to include more realistic schemes that include large relaxation terms. We explore a variety of pulse schemes using multiple model systems using realistic parameters. An open-source software package is provided to perform the density matrix simulation and Wigner transformations. The open-source software allows us to define any arbitrary level schemes as well as any arbitrary electric field in the interaction Hamiltonian. [ABSTRACT FROM AUTHOR]

Seniority-zero geminal wavefunctions are known to capture bond-breaking correlation. Among this class of wavefunctions, Richardson–Gaudin states stand out as they are eigenvectors of a model Hamiltonian. This provides a clear physical picture, clean expressions for reduced density matrix (RDM) elements, and systematic improvement (with a complete set of eigenvectors). Known expressions for the RDM elements require the computation of rapidities, which are obtained by first solving for the so-called eigenvalue based variables (EBV) and then root-finding a Lagrange interpolation polynomial. In this paper, we obtain expressions for the RDM elements directly in terms of the EBV. The final expressions can be computed at the same cost as the rapidity expressions. Therefore, except, in particular, circumstances, it is entirely unnecessary to compute rapidities at all. The RDM elements require numerically inverting a matrix, and while this is usually undesirable, we demonstrate that it is stable, except when there is degeneracy in the single-particle energies. In such cases, a different construction would be required. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the time evolution of quantum coherence—the off-diagonal elements of the density matrix of a multistate quantum system—from the perspective of the Wigner–Moyal formalism. This approach provides an exact phase space representation of quantum mechanics. We consider the coherent evolution of nuclear wavepackets in a molecule with two electronic states. For harmonic potentials, the problem is analytically soluble for both a fully quantum mechanical description and a semiclassical description. We highlight the serious deficiencies of the semiclassical treatment of coherence for general systems and illustrate how even qualitative accuracy requires higher order terms in the Moyal expansion to be included. The model provides an experimentally relevant example of a molecular Schrödinger's cat state. The alive and dead cats of the exact two-state quantum evolution collapse into a "zombie" cat in the semiclassical limit—an averaged behavior, neither alive nor dead, leading to significant errors. The inclusion of the Moyal correction restores a faithful simultaneously alive and dead representation of the cat that is experimentally observable. [ABSTRACT FROM AUTHOR]

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). In this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H2Z2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces. [ABSTRACT FROM AUTHOR]

In recent years, the performance of different entanglement indicators obtained directly from tomograms has been assessed in continuous-variable and hybrid quantum systems. In this paper, we carry out this task in the case of spin systems. We compute the entanglement indicators from actual experimental data obtained from three liquid-state nuclear magnetic resonance (NMR) experiments and compare them with standard entanglement measures calculated from the corresponding density matrices, both experimentally reconstructed and numerically computed. The gross features of entanglement dynamics and spin squeezing properties are found to be reproduced by these entanglement indicators. However, the extent to which these indicators and spin squeezing track the entanglement during time evolution of the multipartite systems in the NMR experiments is very sensitive to the precise nature and strength of interactions as well as the manner in which the full system is partitioned into subsystems. We also use the IBM quantum computer to implement equivalent circuits that capture the dynamics of the multipartite system in one of the NMR experiments and carry out a similar comparative assessment of the performance of tomographic indicators. This exercise shows that these indicators can estimate the degree of entanglement without necessitating detailed state reconstruction procedures, establishing the advantage of the tomographic approach. [ABSTRACT FROM AUTHOR]

One-particle Green's functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green's function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations. [ABSTRACT FROM AUTHOR]

Quantum master equations provide a general framework for describing the dynamics of electronic observables within a complex molecular system. One particular family of such equations is based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. In this paper, we show how different choices of projection operators, as well as whether one starts out with the time-convolution or the time-convolutionless forms of the generalized quantum master equation, give rise to four different types of such off-diagonal quantum master equations (OD-QMEs), namely, time-convolution and time-convolutionless versions of a Pauli-type OD-QME for only the electronic populations and an OD-QME for the full electronic density matrix (including both electronic populations and coherences). The fact that those OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrödinger picture electronic coherences from them. To address this, we also extend a procedure for extracting Schrödinger picture electronic coherences from interaction picture populations recently introduced by Trushechkin in the context of time-convolutionless Pauli-type OD-QME to the other three types of OD-QMEs. The performance of the aforementioned four types of OD-QMEs is explored in the context of the Garg–Onuchic–Ambegaokar benchmark model for charge transfer in the condensed phase across a relatively wide parameter range. The results show that time-convolution OD-QMEs can be significantly more accurate than their time-convolutionless counterparts, particularly in the case of Pauli-type OD-QMEs, and that rather accurate Schrödinger picture coherences can be obtained from interaction picture electronic inputs. [ABSTRACT FROM AUTHOR]

It is important that any dynamics method approaches the correct population distribution at long times. In this paper, we derive a one-body reduced density matrix dynamics for electrons in energetic contact with a bath. We obtain a remarkable equation of motion which shows that in order to reach equilibrium properly, rates of electron transitions depend on the density matrix. Even though the bath drives the electrons towards a Boltzmann distribution, hole blocking factors in our equation of motion cause the electronic populations to relax to a Fermi-Dirac distribution. These factors are an old concept, but we show how they can be derived with a combination of time-dependent perturbation theory and the extended normal ordering of Mukherjee and Kutzelnigg for a general electronic state. The resulting non-equilibrium kinetic equations generalize the usual Redfield theory to many-electron systems, while ensuring that the orbital occupations remain between zero and one. In numerical applications of our equations, we show that relaxation rates of molecules are not constant because of the blocking effect. Other applications to model atomic chains are also presented which highlight the importance of treating both dephasing and relaxation. Finally, we show how the bath localizes the electron density matrix. [ABSTRACT FROM AUTHOR]

Using the formulas and techniques developed in Papers I and II of this series, the recently developed second-order perturbation theory based on a valence bond self-consistent field reference function (VBPT2) has been extended by using the internally contracted correction wave function. This ansatz strongly reduces the size of the interaction space compared to the uncontracted wave function and thus improves the capability of the VBPT2 method dramatically. Test calculations show that internally contracted VBPT2 using only a small number of reference valence bond functions, can give results as accuracy as the VBPT2 method and other more sophisticated methods such as full configuration interaction and multireference configuration interaction. [ABSTRACT FROM AUTHOR]

A full configuration interaction calculation (FCI) ultimately defines the innate molecular orbital description of a molecule. Its density matrix and the natural orbitals obtained from it quantify the difference between having N-dominantly occupied orbitals in a reference determinant for a wavefunction to describe N-correlated electrons and how many of those N-electrons are left to the remaining virtual orbitals. The latter provides a measure of the multi-determinantal character (MDC) required to be in a wavefunction. MDC is further split into a weak correlation part and a part that indicates stronger correlation often called multi-reference character (MRC). If several virtual orbitals have high occupation numbers, then one might argue that these additional orbitals should be allowed to have a larger role in the calculation, as in MR methods, such as MCSCF, MR-CI, or MR-coupled-cluster (MR-CC), to provide adequate approximations toward the FCI. However, there are problems with any of these MR methods that complicate the calculations compared to the uniformity and ease of application of single-reference CC calculations (SR-CC) and their operationally single-reference equation-of-motion (EOM-CC) extensions. As SR-CC theory is used in most of today's "predictive" calculations, an assessment of the accuracy of SR-CC at some truncation of the cluster operator would help to quantify how large an issue MRC actually is in a calculation, and how it might be alleviated while retaining the convenient SR computational character of CC/EOM-CC. This paper defines indices that identify MRC situations and help assess how reliable a given calculation is. [ABSTRACT FROM AUTHOR]

Accurate prediction of water properties in its gas and condensed phases, including the interaction of water with surfaces, is of prime importance for many scientific disciplines. However, accurate simulation of all water properties together within semilocal approximations of the density functional theory possesses great challenges. The Strongly Constrained and Appropriately Normed semilocal density functional, which satisfies 17 known exact constraints and includes the intermediate range van der Waals interaction, performs quite well for different properties of water including the correct energy ordering of isomers. Despite its impressive performance, the energy overestimation for water isomers, ice lattice energies, and volume underestimation for ice are noticeable. However, it is recently shown that [S. Jana et al., J. Chem. Theory Comput. 16(2), 974–987 (2020)] meta-generalized gradient approximations based on the density matrix expansion [i.e., Tao-Mo (TM) and revised TM (revTM)] can achieve quite a good accuracy for the diverse properties of water. In this paper, we assess the performance of the dispersion corrected counterparts of the TM and revTM functionals. It is shown that the dispersion corrected counterparts of both methods are also quite accurate for diverse water properties, especially for the water–solid interactions. Moreover, the extent of accuracy of TM-based functionals is also analyzed from the viewpoint of the density and functional-driven error. Finally, a comparison in the performance of the dispersion corrected functionals is exhibited. It is shown that the "Optimized Power" damping function together with Grimme's D3 correction and revTM functional is in excellent agreement for the water adsorption on carbon nanostructure materials and ice-lattice mismatch problem without deviating accuracy of other water properties compared to its bare functional. [ABSTRACT FROM AUTHOR]

In a previous paper [J. R. Mannouch and J. O. Richardson, J. Chem. Phys. 153, 194109 (2020)], we derived a new partially linearized mapping-based classical-trajectory technique called the spin partially linearized density matrix (spin-PLDM) approach. This method describes the dynamics associated with the forward and backward electronic path integrals using a Stratonovich–Weyl approach within the spin-mapping space. While this is the first example of a partially linearized spin-mapping method, fully linearized spin-mapping is already known to be capable of reproducing dynamical observables for a range of nonadiabatic model systems reasonably accurately. Here, we present a thorough comparison of the terms in the underlying expressions for the real-time quantum correlation functions for spin-PLDM and fully linearized spin mapping in order to ascertain the relative accuracy of the two methods. In particular, we show that spin-PLDM contains an additional term within the definition of its real-time correlation function, which diminishes many of the known errors that are ubiquitous for fully linearized approaches. One advantage of partially linearized methods over their fully linearized counterparts is that the results can be systematically improved by re-sampling the mapping variables at intermediate times. We derive such a scheme for spin-PLDM and show that for systems for which the approximation of classical nuclei is valid, numerically exact results can be obtained using only a few "jumps." Additionally, we implement focused initial conditions for the spin-PLDM method, which reduces the number of classical trajectories that are needed in order to reach convergence of dynamical quantities, with seemingly little difference to the accuracy of the result. [ABSTRACT FROM AUTHOR]

We present a new partially linearized mapping-based approach for approximating real-time quantum correlation functions in condensed-phase nonadiabatic systems, called the spin partially linearized density matrix (spin-PLDM) approach. Within a classical trajectory picture, partially linearized methods treat the electronic dynamics along forward and backward paths separately by explicitly evolving two sets of mapping variables. Unlike previously derived partially linearized methods based on the Meyer–Miller–Stock–Thoss mapping, spin-PLDM uses the Stratonovich–Weyl transform to describe the electronic dynamics for each path within the spin-mapping space; this automatically restricts the Cartesian mapping variables to lie on a hypersphere and means that the classical equations of motion can no longer propagate the mapping variables out of the physical subspace. The presence of a rigorously derived zero-point energy parameter also distinguishes spin-PLDM from other partially linearized approaches. These new features appear to give the method superior accuracy for computing dynamical observables of interest when compared with other methods within the same class. The superior accuracy of spin-PLDM is demonstrated in this paper through application of the method to a wide range of spin-boson models as well as to the Fenna–Matthews–Olsen complex. [ABSTRACT FROM AUTHOR]

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo (FCIQMC) algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24 000 central processing unit cores. In this paper, we describe the core functionalities of NECI and its recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing us to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs, and it is licensed under GPL-3.0. [ABSTRACT FROM AUTHOR]

We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange–correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications. [ABSTRACT FROM AUTHOR]

There are three essential problems in computational relativistic chemistry: Electrons moving at relativistic speeds, close lying states, and dynamical correlation. Currently available quantum-chemical methods are capable of solving systems with one or two of these issues. However, there is a significant class of molecules in which all the three effects are present. These are the heavier transition metal compounds, lanthanides, and actinides with open d or f shells. For such systems, sufficiently accurate numerical methods are not available, which hinders the application of theoretical chemistry in this field. In this paper, we combine two numerical methods in order to address this challenging class of molecules. These are the relativistic versions of coupled cluster methods and the density matrix renormalization group (DMRG) method. To the best of our knowledge, this is the first relativistic implementation of the coupled cluster method externally corrected by DMRG. The method brings a significant reduction of computational costs as we demonstrate on the system of TlH, AsH, and SbH. [ABSTRACT FROM AUTHOR]

The time dependent density matrix renormalization group (TD-DMRG) has become one of the cutting edge methods of quantum dynamics for complex systems. In this paper, we comparatively study the accuracy of three time evolution schemes in the TD-DMRG, the global propagation and compression method with the Runge-Kutta algorithm (P&C-RK), the time dependent variational principle based methods with the matrix unfolding algorithm (TDVP-MU), and with the projector-splitting algorithm (TDVP-PS), by performing benchmarks on the exciton dynamics of the Fenna-Matthews-Olson complex. We show that TDVP-MU and TDVP-PS yield the same result when the time step size is converged and they are more accurate than P&C-RK4, while TDVP-PS tolerates a larger time step size than TDVP-MU. We further adopt the graphical processing units to accelerate the heavy tensor contractions in the TD-DMRG, and it is able to speed up the TDVP-MU and TDVP-PS schemes by up to 73 times. [ABSTRACT FROM AUTHOR]

In this paper, we consider the momentum operator of a quantum particle directed along the displacement of two of its neighbors. A modified open-path path integral molecular dynamics is presented to sample the distribution of this directional momentum distribution, where we derive and use a new estimator for this distribution. Variationally enhanced sampling is used to obtain this distribution for an example molecule, malonaldehyde, in the very low temperature regime where deep tunneling happens. We find no secondary feature in the directional momentum distribution and that its absence is due to quantum entanglement through a further study of the reduced density matrix. [ABSTRACT FROM AUTHOR]

This paper presents techniques for Fock matrix construction that are designed for high performance on shared and distributed memory parallel computers when using Gaussian basis sets. Four main techniques are considered. (1) To calculate electron repulsion integrals, we demonstrate batching together the calculation of multiple shell quartets of the same angular momentum class so that the calculation of large sets of primitive integrals can be efficiently vectorized. (2) For multithreaded summation of entries into the Fock matrix, we investigate using a combination of atomic operations and thread-local copies of the Fock matrix. (3) For distributed memory parallel computers, we present a globally accessible matrix class for accessing distributed Fock and density matrices. The new matrix class introduces a batched mode for remote memory access that can reduce the synchronization cost. (4) For density fitting, we exploit both symmetry (of the Coulomb and exchange matrices) and sparsity (of 3-index tensors) and give a performance comparison of density fitting and the conventional direct calculation approach. The techniques are implemented in an open-source software library called GTFock. [ABSTRACT FROM AUTHOR]

The generalized quantum master equation (GQME) provides a powerful framework for simulating electronically nonadiabatic molecular dynamics. Within this framework, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density matrix is fully captured by a memory kernel superoperator. In this paper, we consider two different procedures for calculating the memory kernel of the GQME from projection-free inputs obtained via the combination of the mapping Hamiltonian (MH) approach and the linearized semiclassical (LSC) approximation. The accuracy and feasibility of the two procedures are demonstrated on the spin-boson model. We find that although simulating the electronic dynamics by direct application of the two LSC-based procedures leads to qualitatively different results that become increasingly less accurate with increasing time, restricting their use to calculating the memory kernel leads to an accurate description of the electronic dynamics. Comparison with a previously proposed procedure for calculating the memory kernel via the Ehrenfest method reveals that MH/LSC methods produce memory kernels that are better behaved at long times and lead to more accurate electronic dynamics. [ABSTRACT FROM AUTHOR]

The Caldeira-Leggett model provides a compact characterization of a thermal environment in terms of a spectral density function, which has led to a variety of numerically exact quantum methods for reduced density matrix propagation. Since spectral densities are often computed from classical molecular dynamics simulations, we investigate in this paper whether quantum effects should be accounted for in the calculations. Therefore, we reformulate the recently developed Fourier method for spectral density calculations from semiclassical simulations which approximately allow for quantum effects. We propose two possible protocols based on either correlation functions or expectation values. These protocols are tested on a generic Calderra-Leggett model for the linearized semiclassical initial-value representation (LSC-IVR), the thawed Gaussian wave packet dynamics (TGWD), and hybrid schemes combining the two with the more accurate Herman-Kluk formula. Surprisingly, spectral densities from the LSC-IVR method, which treats the dynamics completely classically, are found to be extremely accurate, even in the quantum regime, where this method does not give a correct description of the correlation functions and expectation values. In contrast, the TGWD method turns out as too inaccurate for spectral density calculations, and the hybrid schemes perform well only if the system is close to the classical regime. This implies that, if the bath has a Caldeira-Leggett form, spectral densities are insensitive to quantum effects and any effort to approximately account for them rather leads to errors. Hence, in this case, spectral densities can be computed from classical simulations and used in a reduced quantum simulation as well. [ABSTRACT FROM AUTHOR]

The molecular dynamics following the electronic B Π u 3 0 + ⟵ X Σ g + 1 photoexcitation of the iodine molecule embedded in solid krypton are studied quantum mechanically using the Gaussian variant of the multiconfigurational time-dependent Hartree method (G-MCTDH). The accuracy of the Gaussian wave packet approximation is validated against numerically exact MCTDH simulations for a fully anharmonic seven-dimensional model of the I2Kr18 cluster in a crystal Kr cage. The linear absorption spectrum, time-evolving vibrational probability densities, and I2 energy expectation value are accurately reproduced by the numerically efficient G-MCTDH approach. The reduced density matrix of the chromophore is analyzed in the coordinate, Wigner and energy representations, so as to obtain a multifaceted dynamical view of the guest-host interactions. Vibrational coherences extending over the bond distance range 2.7 Å < RI–I < 4.0 Å are found to survive for several vibrational periods, despite extensive dissipation. The present results prepare the ground for the simulation of time-resolved coherent Raman spectroscopy of the I2-krypton system addressed in Paper II. [ABSTRACT FROM AUTHOR]

Quantum chemistry is increasingly performed using large cluster computers consisting of multiple interconnected nodes. For a fixed molecular problem, the efficiency of a calculation usually decreases as more nodes are used, due to the cost of communication between the nodes. This paper empirically investigates the parallel scalability of Hartree-Fock calculations. The construction of the Fock matrix and the density matrix calculation are analyzed separately. For the former, we use a parallelization of Fock matrix construction based on a static partitioning of work followed by a work stealing phase. For the latter, we use density matrix purification from the linear scaling methods literature, but without using sparsity. When using large numbers of nodes for moderately sized problems, density matrix computations are network-bandwidth bound, making purification methods potentially faster than eigendecomposition methods. [ABSTRACT FROM AUTHOR]

Recently it has been shown how to employ the adiabatic connection (AC) formalism to obtain correlation energy for multireference wavefunctions [K. Pernal, Phys. Rev. Lett. 120, 013001 (2018)]. Approximations to the exact AC formulation have been based on assuming that a one-electron reduced density matrix is constant along the AC path and by employing the extended random phase approximation. In this paper, the importance of these approximations is examined by comparing approximate AC integrands with their exact counterparts obtained for the hydrogen molecule in its ground and excited states. Encouraging results obtained for H2 indicate that AC is a viable and promising approach to a correlation energy problem not only for ground but also for excited states of electronic systems. [ABSTRACT FROM AUTHOR]

A simple molecular orbital treatment of local current distributions inside single molecular junctions is developed in this paper. Using the first-order perturbation theory and nonequilibrium Green's function techniques in the framework of Hückel theory, we show that the leading contributions to local current distributions are directly proportional to the off-diagonal elements of transition density matrices. Under the orbital approximation, the major contributions to local currents come from a few dominant molecular orbital pairs which are mixed by the interactions between the molecule and electrodes. A few simple molecular junctions consisting of single- and multi-ring conjugated systems are used to demonstrate that local current distributions inside molecular junctions can be decomposed by partial sums of a few leading contributing transition density matrices. [ABSTRACT FROM AUTHOR]

Full configuration interaction (FCI) restricted to a pairing space yields size-extensive correlation energies but its cost scales exponentially with molecular size. Restricting the variational two-electron reduced-density-matrix (2-RDM) method to represent the same pairing space yields an accurate lower bound to the pair FCI energy at a mean-field-like computational scaling of O(r3) where r is the number of orbitals. In this paper, we show that localized molecular orbitals can be employed to generate an efficient, approximately size-extensive pair 2-RDM method. The use of localized orbitals eliminates the substantial cost of optimizing iteratively the orbitals defining the pairing space without compromising accuracy. In contrast to the localized orbitals, the use of canonical Hartree-Fock molecular orbitals is shown to be both inaccurate and non-size-extensive. The pair 2-RDM has the flexibility to describe the spectra of one-electron RDM occupation numbers from all quantum states that are invariant to time-reversal symmetry. Applications are made to hydrogen chains and their dissociation, n-acene from naphthalene through octacene, and cadmium telluride 2-, 3-, and 4-unit polymers. For the hydrogen chains, the pair 2-RDM method recovers the majority of the energy obtained from similar calculations that iteratively optimize the orbitals. The localized-orbital pair 2-RDM method with its mean-field-like computational scaling and its ability to describe multireference correlation has important applications to a range of strongly correlated phenomena in chemistry and physics. [ABSTRACT FROM AUTHOR]

Noise-assisted energy transfer can be explained geometrically in terms of the set of one-electron reduced density matrices (1-RDMs) [R. Chakraborty and D. A. Mazziotti, Phys. Rev.A 91, 010101(R) (2015)]. In this paper, we examine the geometric picture of quantum noise for the seven-chromophore Fenna-Matthews-Olson (FMO) complex. Noise expands the feasible set of orbital occupation trajectories to the target state through the violation of the pure-state N-representability conditions on the 1-RDM, known as the generalized Pauli constraints. While the generalized Pauli constraints are not explicitly known for seven-electron systems, we are able to treat a seven-exciton model of the FMO complex through the use of generalized Pauli constraints for p qubits which are known for arbitrary p. In the model, we find that while dephasing noise alone produces a trajectory of ensemble states that neither expands the set of 1-RDMs nor reaches the reaction center, the inclusion of both dephasing and dissipation expands the set of 1-RDMs and exhibits an efficient energy transfer to the reaction center. The degree to which the noise expands the set of 1-RDMs, violating the generalized Pauli constraints, is quantified by the distance of the 1-RDM outside its pure set to the distance of the 1-RDM inside its ensemble set. The geometric picture of energy transfer has applications to general quantum systems in chemistry and physics. [ABSTRACT FROM AUTHOR]

The density-matrix renormalization group (DMRG) method, which can deal with a large active space composed of tens of orbitals, is nowadays widely used as an efficient addition to traditional complete active space (CAS)-based approaches. In this paper, we present the DMRG algorithm with a multi-level (ML) control of the active space based on chemical intuition-based hierarchical orbital ordering, which is called as ML-DMRG with its self-consistent field (SCF) variant MLDMRG- SCF. Ground and excited state calculations of H2O, N2, indole, and Cr2 with comparisons to DMRG references using fixed number of kept states (M) illustrate that ML-type DMRG calculations can obtain noticeable efficiency gains. It is also shown that the orbital re-ordering based on hierarchical multiple active subspaces may be beneficial for reducing computational time for not only ML-DMRG calculations but also DMRG ones with fixed M values. [ABSTRACT FROM AUTHOR]

In this contribution, we report two different methodologies for characterizing the electronic structure reorganization occurring when a chromophore undergoes an electronic transition. For the first method, we start by setting the theoretical background necessary to the reinterpretation through simple tensor analysis of (i) the transition density matrix and (ii) the natural transition orbitals in the scope of reduced density matrix theory. This novel interpretation is made more clear thanks to a short compendium of the one-particle reduced density matrix theory in a Fock space. The formalism is further applied to two different classes of excited states calculation methods, both requiring a single-determinant reference, that express an excited state as a hole-particle mono-excited configurations expansion, to which particle-hole correlation is coupled (time-dependent Hartree-Fock/timedependent density functional theory) or not (configuration interaction single/Tamm-Dancoff approximation). For the second methodology presented in this paper, we introduce a novel and complementary concept related to electronic transitions with the canonical transition density matrix and the canonical transition orbitals. Their expression actually reflects the electronic cloud polarisation in the orbital space with a decomposition based on the actual contribution of one-particle excitations from occupied canonical orbitals to virtual ones. This approach validates our novel interpretation of the transition density matrix elements in terms of the Euclidean norm of elementary transition vectors in a linear tensor space. A proper use of these new concepts leads to the conclusion that despite the different principles underlying their construction, they provide two equivalent excited states topological analyses. This connexion is evidenced through simple illustrations of (in)organic dyes electronic transitions analysis. [ABSTRACT FROM AUTHOR]

A new approach for obtaining interacting quantum atoms-defined components of binding energy of intermolecular interactions, which bypasses the use of standard six-dimensional integrals and two-particle reduced density matrix (2-RDM) reconstruction, is proposed. To examine this approach, three datasets calculated within the density functional theory framework using the def2-TZVP basis have been explored. The first two, containing 53 weakly bound bimolecular associates and 13 molecular clusters taken from the crystal, were used in protocol refinement, and the third one containing other 20 bimolecular and three cluster systems served as a validation reference. In addition, to verify the performance of the proposed approach on an exact 2-RDM, calculations within the coupled cluster formalism were performed for part of the first set systems using the cc-pVTZ basis set. The process of optimization of the proposed parametric model is considered, and the role of various energy contributions in the formation of non-covalent interactions is discussed with regard to the obtained trends. [ABSTRACT FROM AUTHOR]

Accurate multireference electronic structure calculations are important for constructing potential energy surfaces. Still, even in the case of low-scaling methods, their routine use is limited by the steep growth of the computational and storage costs as the active space grows. This is primarily due to the occurrence of three- and higher-body density matrices or, equivalently, their cumulants. This work examines the effect of various cumulant truncation schemes on the accuracy of the driven similarity renormalization group second-order multireference perturbation theory. We test four different levels of three-body reduced density cumulant truncations that set different classes of cumulant elements to zero. Our test cases include the singlet–triplet gap of CH2, the potential energy curves of the X Σ g + 1 and A Σ u + 3 states of N2, and the singlet–triplet splittings of oligoacenes. Our results show that both relative and absolute errors introduced by these cumulant truncations can be as small as 0.5 kcal mol−1 or less. At the same time, the amount of memory required is reduced from O ( N A 6 ) to O ( N A 5 ) , where NA is the number of active orbitals. No additional regularization is needed to prevent the intruder state problem in the cumulant-truncated second-order driven similarity renormalization group multireference perturbation theory methods. [ABSTRACT FROM AUTHOR]

We develop an accurate and numerically efficient non-adiabatic path-integral approach to simulate the non-linear spectroscopy of exciton–polariton systems. This approach is based on the partial linearized density matrix approach to model the exciton dynamics with explicit propagation of the phonon bath environment, combined with a stochastic Lindblad dynamics approach to model the cavity loss dynamics. Through simulating both linear and polariton two-dimensional electronic spectra, we systematically investigate how light–matter coupling strength and cavity loss rate influence the optical response signal. Our results confirm the polaron decoupling effect, which is the reduced exciton–phonon coupling among polariton states due to the strong light–matter interactions. We further demonstrate that the polariton coherence time can be significantly prolonged compared to the electronic coherence outside the cavity. [ABSTRACT FROM AUTHOR]

The equivalence in one-electron quantum baths between the practical implementation of density matrix embedding theory (DMET) and the more recent Householder-transformed density matrix functional embedding theory has been shown previously in the standard but special case where the reference full-size (one-electron reduced) density matrix, from which the bath is constructed, is idempotent [S. Yalouz et al., J. Chem. Phys. 157, 214112 (2022)]. We prove mathematically that the equivalence remains valid when the density matrix is not idempotent anymore, thus allowing for the construction of correlated (one-electron) quantum baths. A density-matrix functional exactification of DMET is derived within the present unified quantum embedding formalism. Numerical examples reveal that the embedding cluster can be quite sensitive to the level of density-matrix functional approximation used for computing the reference density matrix. [ABSTRACT FROM AUTHOR]

Electron occupations that arise from pure quantum states are restricted by a stringent set of conditions that are said to generalize the Pauli exclusion principle. These generalized Pauli constraints (GPCs) define the boundary of the set of one-electron reduced density matrices (1-RDMs) that are derivable from at least one N-electron wavefunction. In this paper, we investigate the sparsity of the Slater-determinant representation of the wavefunction that is a necessary, albeit not sufficient, condition for its 1-RDM to lie on the boundary of the set of pure N-representable 1-RDMs or in other words saturate one of the GPCs. The sparse wavefunction, we show, is exact not only for 3 electrons in 6 orbitals but also for 3 electrons in 8 orbitals. For larger numbers of electrons and/or orbitals in the lowest spin state, the exact wavefunction does not generally saturate one of the GPCs, and hence, the sparse representation is typically an approximation. Because the sparsity of the wavefunction is a necessary but not sufficient condition for saturation of one of the GPCs, optimization of the sparse wavefunction Ansatz to minimize the ground-state energy does not necessarily produce a wavefunction whose 1-RDM exactly saturates one of the GPCs. While the sparse Ansatz can be employed with arbitrary orbitals or optimized orbitals, in this paper, we explore the Ansatz with the natural orbitals from full configuration interaction, which yields an upper bound to the ground-state energy that equals the exact energy for a given basis set if the full-configuration-interaction wavefunction saturates the Ansatz’s GPC. With calculations on the boron isoelectronic sequence, the dinitrogen cation N 2 + , hydrogen chains, and cyclic conjugated π systems, we examine the quality of the sparse wavefunction Ansatz from the amount of correlation energy recovered. [ABSTRACT FROM AUTHOR]

We present a machine learning approach to calculating electronic specific heat capacities for a variety of benchmark molecular systems. Our models are based on data from density matrix quantum Monte Carlo, which is a stochastic method that can calculate the electronic energy at finite temperature. As these energies typically have noise, numerical derivatives of the energy can be challenging to find reliably. In order to circumvent this problem, we use Gaussian process regression to model the energy and use analytical derivatives to produce the specific heat capacity. From there, we also calculate the entropy by numerical integration. We compare our results to cubic splines and finite differences in a variety of molecules in which Hamiltonians can be diagonalized exactly with full configuration interaction. We finally apply this method to look at larger molecules where exact diagonalization is not possible and make comparisons with more approximate ways to calculate the specific heat capacity and entropy. [ABSTRACT FROM AUTHOR]

The recently reported Grassmann interpolation (G-Int) method [J. A. Tan and K. U. Lao, J. Chem. Phys. 158, 051101 (2023)] has been extended to spin-unrestricted open-shell systems. In contrast to closed-shell systems, where G-Int has to be performed only once since the α and β density matrices are the same, spin-unrestricted open-shell systems require G-Int to be performed twice—one for the α spin and another for the β spin density matrix. In this work, we tested the performance of G-Int to the carbon monoxide radical cation CO●+ and nickelocene complex, which have the doublet and triple ground states, respectively. We found that the Frobenius norm errors associated with the interpolations for the α and β spin density matrices are comparable for a given molecular geometry. These G-Int density matrices, when used as an initial guess for a self-consistent field (SCF) calculation, outperform the conventional SCF guess schemes, such as the superposition of atomic densities, purified superposition of atomic densities, core Hamiltonian, and generalized Wolfsberg–Helmholtz approximation. Depending on the desired accuracy, these G-Int density matrices can be used to directly evaluate the SCF energy without performing SCF iterations. In addition, the spin-unrestricted G-Int density matrices have been used for the first time to directly calculate the atomic charges using the Mulliken and ChElPG population analysis. [ABSTRACT FROM AUTHOR]

We investigate the adiabatic approximation to the exact-exchange kernel for calculating correlation energies within the adiabatic-connection fluctuation–dissipation framework of time-dependent density functional theory. A numerical study is performed on a set of systems having bonds of different character (H2 and N2 molecules, H-chain, H2-dimer, solid-Ar, and the H2O-dimer). We find that the adiabatic kernel can be sufficient in strongly bound covalent systems, yielding similar bond lengths and binding energies. However, for non-covalent systems, the adiabatic kernel introduces significant errors around equilibrium geometry, systematically overestimating the interaction energy. The origin of this behavior is investigated by studying a model dimer composed of one-dimensional, closed-shell atoms, interacting via soft-Coulomb potentials. The kernel is shown to exhibit a strong frequency dependence at small to intermediate atomic separation that affects both the low-energy spectrum and the exchange-correlation hole obtained from the corresponding diagonal of the two-particle density matrix. [ABSTRACT FROM AUTHOR]

A way to improve the accuracy of the spectral properties in density functional theory (DFT) is to impose constraints on the effective, Kohn–Sham (KS), local potential [J. Chem. Phys. 136, 224109 (2012)]. As illustrated, a convenient variational quantity in that approach is the "screening" or "electron repulsion" density, ρrep, corresponding to the local, KS Hartree, exchange and correlation potential through Poisson's equation. Two constraints, applied to this minimization, largely remove self-interaction errors from the effective potential: (i) ρrep integrates to N − 1, where N is the number of electrons, and (ii) ρrep ≥ 0 everywhere. In this work, we introduce an effective "screening" amplitude, f, as the variational quantity, with the screening density being ρrep = f2. In this way, the positivity condition for ρrep is automatically satisfied, and the minimization problem becomes more efficient and robust. We apply this technique to molecular calculations, employing several approximations in DFT and in reduced density matrix functional theory. We find that the proposed development is an accurate, yet robust, variant of the constrained effective potential method. [ABSTRACT FROM AUTHOR]

Since its first definition, back in 1990, the electron localization function (ELF) has settled as one of the most commonly employed techniques to characterize the nature of the chemical bond in real space. Although most of the work using the ELF has focused on the study of ground-state chemical reactivity, a growing interest has blossomed to apply these techniques to the nearly unexplored realm of excited states and photochemistry. Since accurate excited electronic states usually require to account appropriately for electron correlation, the standard single-determinant ELF formulation cannot be blindly applied to them, and it is necessary to turn to correlated ELF descriptions based on the two-particle density matrix (2-PDM). The latter requires costly wavefunction approaches, unaffordable for most of the systems of current photochemical interest. Here, we compare the exact, 2-PDM-based ELF results with those of approximate 2-PDM reconstructions taken from reduced density matrix functional theory. Our approach is put to the test in a wide variety of representative scenarios, such as those provided by the lowest-lying excited electronic states of simple diatomic and polyatomic molecules. Altogether, our results suggest that even approximate 2-PDMs are able to accurately reproduce, on a general basis, the topological and statistical features of the ELF scalar field, paving the way toward the application of cost-effective methodologies, such as time-dependent-Hartree–Fock or time-dependent density functional theory, in the accurate description of the chemical bonding in excited states of photochemical relevance. [ABSTRACT FROM AUTHOR]

Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1–6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1–2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters. [ABSTRACT FROM AUTHOR]

The leading terms in the large-R asymptotics of the functional of the one-electron reduced density matrix for the ground-state energy of the H2 molecule with the internuclear separation R are derived thanks to the solution of the phase dilemma at the R → ∞ limit. At this limit, the respective natural orbitals (NOs) are given by symmetric and antisymmetric combinations of "half-space" orbitals with the corresponding natural amplitudes having the same amplitudes but opposite signs. Minimization of the resulting explicit functional yields the large-R asymptotics for the occupation numbers of the weakly occupied NOs and the C6 dispersion coefficient. The highly accurate approximates for the radial components of the p-type "half-space" orbitals and the corresponding occupation numbers (that decay like R−6), which are available for the first time thanks to the development of the present formalism, have some unexpected properties. [ABSTRACT FROM AUTHOR]

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. We present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost. [ABSTRACT FROM AUTHOR]