Stochastic Resolution of Identity for Real-Time Second-Order Green’s Function: Ionization Potential and Quasi-Particle Spectrum

Stochastic Resolution of Identity to CC2 for Large Systems: Excited State Properties

We apply a stochastic resolution of identity approximation (sRI) to the CC2 method for the excitation energy calculations. A set of stochastic orbitals are employed to decouple the crucial 4-index electron repulsion integrals and optimize the contraction steps in CC2 response theory. The CC2 response for excitations builds upon sRI-CC2 ground-state calculations, which scales as O(N3), where N is a measure for the system size. Overall, the current algorithm for excited states also allows a sharp scaling reduction from original O(N5) to O(N3). We test the sRI-CC2 for different molecular systems and basis sets, and we show that our sRI-CC2 method can accurately reproduce the results of the deterministic CC2 approach. Our sRI-CC2 exhibits an experimental scaling of O(N2.59) for a series of olefin chains, allowing us to calculate systems with nearly thousands of electrons.

Stochastic Real-Time Second-Order Green's Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green's function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(Ne3), where Ne is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

Interpolative Separable Density Fitting for Accelerating Two-Electron Integrals: A Theoretical Perspective

Low-rank approximations have long been considered an efficient way to accelerate electronic structure calculations associated with the evaluation of electron repulsion integrals (ERIs). As an accurate and efficient algorithm for compressing the ERI tensor, the interpolative separable density fitting (ISDF) decomposition has recently attracted great attention in this context. In this perspective, we introduce the ISDF decomposition from the theoretical aspects and technique details. The ISDF decomposition can construct a fully separable low-rank approximation (tensor hypercontraction factorization) of ERIs in real space with a cubic cost, offering great flexibility for accelerating high-scaling electronic structure calculations. We review the typical applications of ISDF in hybrid functionals, time-dependent density functional theory, and GW approximation. Finally, we discuss the promising directions for future development of ISDF.

Time-Dependent Second-Order Green's Function Theory for Neutral Excitations

We develop a time-dependent second-order Green's function theory (GF2) for calculating neutral excited states in molecules. The equation of motion for the lesser Green's function (GF) is derived within the adiabatic approximation to the Kadanoff-Baym (KB) equation, using the second-order Born approximation for the self-energy. In the linear response regime, we recast the time-dependent KB equation into a Bethe-Salpeter-like equation (GF2-BSE), with a kernel approximated by the second-order Coulomb self-energy. We then apply our GF2-BSE to a set of molecules and atoms and find that GF2-BSE is superior to configuration interaction with singles (CIS) and/or time-dependent Hartree-Fock (TDHF), particularly for charge-transfer excitations, and is comparable to CIS with perturbative doubles (CIS(D)) in most cases.

Simulations of nonradiative processes in semiconductor nanocrystals

The description of carrier dynamics in spatially confined semiconductor nanocrystals (NCs), which have enhanced electron-hole and exciton-phonon interactions, is a great challenge for modern computational science. These NCs typically contain thousands of atoms and tens of thousands of valence electrons with discrete spectra at low excitation energies, similar to atoms and molecules, that converge to the continuum bulk limit at higher energies. Computational methods developed for molecules are limited to very small nanoclusters, and methods for bulk systems with periodic boundary conditions are not suitable due to the lack of translational symmetry in NCs. This perspective focuses on our recent efforts in developing a unified atomistic model based on the semiempirical pseudopotential approach, which is parameterized by first-principle calculations and validated against experimental measurements, to describe two of the main nonradiative relaxation processes of quantum confined excitons: exciton cooling and Auger recombination. We focus on the description of both electron-hole and exciton-phonon interactions in our approach and discuss the role of size, shape, and interfacing on the electronic properties and dynamics for II-VI and III-V semiconductor NCs.

Reduced Scaling of Optimal Regional Orbital Localization via Sequential Exhaustion of the Single-Particle Space

Wannier functions have become a powerful tool in the electronic structure calculations of extended systems. The generalized Pipek-Mezey Wannier functions exhibit appealing characteristics (e.g., reaching an optimal localization and the separation of the σ–π orbitals) compared with other schemes. However, when applied to giant nanoscale systems, the orbital localization suffers from a large computational cost overhead if one is interested in localized states in a small fragment of the system. Herein, we present a swift, efficient, and robust approach for obtaining regionally localized orbitals of a subsystem within the generalized Pipek-Mezey scheme. The proposed algorithm introduces a reduced work space and sequentially exhausts the entire orbital space until the convergence of the localization functional. It tackles systems with ∼10000 electrons within 0.5 h with no loss in localization quality compared to the traditional approach. Regionally localized orbitals with a higher extent of localization are obtained via judiciously extending the subsystem’s size. Exemplifying on large bulk and a 4 nm wide slab of diamond with an NV^– center, we demonstrate the methodology and discuss how the choice of the localization region affects the excitation energy of the defect. Furthermore, we show how the sequential algorithm is easily extended to stochastic methodologies that do not provide individual single-particle eigenstates. It is thus a promising tool to obtain regionally localized states for solving the electronic structure problems of a subsystem embedded in giant condensed systems.

Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-CDF (cumulative distribution function) method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with demonstrated 92% strong scalability going from 64 to 4096 processors for a fixed job size.

Stochastic Vector Techniques in Ground-State Electronic Structure

We review a suite of stochastic vector computational approaches for studying the electronic structure of extended condensed matter systems. These techniques help reduce algorithmic complexity, facilitate efficient parallelization, simplify computational tasks, accelerate calculations, and diminish memory requirements. While their scope is vast, we limit our study to ground-state and finite temperature density functional theory (DFT) and second-order perturbation theory. More advanced topics, such as quasiparticle (charge) and optical (neutral) excitations and higher-order processes, are covered elsewhere. We start by explaining how to use stochastic vectors in computations, characterizing the associated statistical errors. Next, we show how to estimate the electron density in DFT and discuss highly effective techniques to reduce statistical errors. Finally, we review the use of stochastic vector techniques for calculating correlation energies within the second-order Møller-Plesset perturbation theory and its finite temperature variational form. Example calculation results are presented and used to demonstrate the efficacy of the methods. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Monte Carlo MP2-F12 for Noncovalent Interactions: The C60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green's functions in real space and imaginary time. These integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol^-1 for the C₆₀ dimer at the center distance of 9.748 Å.

Stochastic evaluation of fourth-order many-body perturbation energies

A scalable, stochastic algorithm evaluating the fourth-order many-body perturbation (MP4) correction to energy is proposed. Three hundred Goldstone diagrams representing the MP4 correction are computer generated and then converted into algebraic formulas expressed in terms of Green's functions in real space and imaginary time. They are evaluated by the direct (i.e., non-Markov, non-Metropolis) Monte Carlo (MC) integration accelerated by the redundant-walker and control-variate algorithms. The resulting MC-MP4 method is efficiently parallelized and is shown to display O(n^5.3) size-dependence of cost, which is nearly two ranks lower than the O(n⁷) dependence of the deterministic MP4 algorithm. It evaluates the MP4/aug-cc-pVDZ energy for benzene, naphthalene, phenanthrene, and corannulene with the statistical uncertainty of 10 mE_h (1.1% of the total basis-set correlation energy), 38 mE_h (2.6%), 110 mE_h (5.5%), and 280 mE_h (9.0%), respectively, after about 10⁹ MC steps.

Spin relaxation in radical pairs from the stochastic Schrödinger equation

We show that the stochastic Schrödinger equation (SSE) provides an ideal way to simulate the quantum mechanical spin dynamics of radical pairs. Electron spin relaxation effects arising from fluctuations in the spin Hamiltonian are straightforward to include in this approach, and their treatment can be combined with a highly efficient stochastic evaluation of the trace over nuclear spin states that is required to compute experimental observables. These features are illustrated in example applications to a flavin-tryptophan radical pair of interest in avian magnetoreception and to a problem involving spin-selective radical pair recombination along a molecular wire. In the first of these examples, the SSE is shown to be both more efficient and more widely applicable than a recent stochastic implementation of the Lindblad equation, which only provides a valid treatment of relaxation in the extreme-narrowing limit. In the second, the exact SSE results are used to assess the accuracy of a recently proposed combination of Nakajima-Zwanzig theory for the spin relaxation and Schulten-Wolynes theory for the spin dynamics, which is applicable to radical pairs with many more nuclear spins. We also analyze the efficiency of trace sampling in some detail, highlighting the particular advantages of sampling with SU(N) coherent states.

Auger Recombination Lifetime Scaling for Type I and Quasi-Type II Core/Shell Quantum Dots

Having already achieved near-unity quantum yields, with promising properties for light-emitting diode, lasing, and charge separation applications, colloidal core/shell quantum dots have great technological potential. The shell thickness and band alignment of the shell and core materials are known to influence the efficiency of these devices. In many such applications, improving the efficiency requires a deep understanding of multiexcitonic states. Herein, we elucidate the shell thickness and band alignment dependencies of the biexciton Auger recombination lifetime for quasi-type II CdSe/CdS and type I CdSe/ZnS core/shell quantum dots. We find that the biexciton Auger recombination lifetime increases with the total nanocrystal volume for quasi-type II CdSe/CdS core/shell quantum dots and is independent of the shell thickness for type I CdSe/ZnS core/shell quantum dots. To perform these calculations and compute Auger recombination lifetimes, we developed a low-scaling approach based on the stochastic resolution of identity. The numerical approach provided a framework for studying the scaling of the biexciton Auger recombination lifetimes in terms of the shell thickness dependencies of the exciton radii, Coulomb couplings, and density of final states in quasi-type II CdSe/CdS and type I CdSe/ZnS core/shell quantum dots.