A stochastic orbital approach to the resolution of identity (RI) approximation for 4 index electron repulsion integrals (ERIs) is presented. The stochastic RI-ERIs are then applied to second order Møller-Plesset perturbation theory (MP2) utilizing a multiple stochastic orbital approach. The introduction of multiple stochastic orbitals results in an O(N_AO^3 ) scaling for both the stochastic RI-ERIs and stochastic RI-MP2, NAO being the number of basis functions. For a range of water clusters we demonstrate that this method exhibits a small prefactor and observed scalings of O(Ne^2.4) for total energies and O(Ne^3.1) for forces (Ne being the number of correlated electrons), outperforming MP2 for clusters with as few as 21 water molecules.

Quasiparticle (QP) excitations are extremely important for understanding and predicting charge transfer and transport in molecules, nanostructures and extended systems. Since density functional theory (DFT) within Kohn-Sham (KS) formulation does not provide reliable QP energies, a many-body perturbation technique within the GW approximation are essential. The steep computational scaling of GW prohibits its use in extended, open boundary, systems with thousands of electrons and more. Recently, a stochastic formulation of GW has been proposed [Phys. Rev. Lett. 113, 076402 (2014)] which scales nearly linearly with the system size, as illustrated for a series of silicon nanocrystals exceeding 3000 electrons. Here, we implement the stochastic GW (sGW) approach to study the ionization potential (IP) of a subset of molecules taken from the "GW 100" benchmark. We show that sGW provides a reliable results in comparison to GW WEST code and to experimental results, numerically establishing its validity. For completeness, we also provide a detailed review of sGW and a summary of the numerical algorithm.

The second-order Matsubara Green’s function method (GF2) is a robust temperature dependent quantum chemistry approach, extending beyond the random-phase approximation. However, till now the scope of GF2 applications was quite limited as they require computer resources which rise steeply with system size. In each step of the self-consistent GF2 calculation there are two parts: the estimation of the self-energy from the previous step’s Green’s function, and updating the Green’s function from the self-energy. The first part formally scales as the fifth power of the system size while the second has a much gentler cubic scaling. Here, we develop a stochastic approach to GF2 (sGF2) which reduces the fifth power scaling of the first step to merely quadratic, leaving the overall sGF2 scaling as cubic. We apply the method to linear hydrogen chains containing up to 1000 electrons, showing that the approach is numerically stable, efficient and accurate. The stochastic errors are very small, of the order of 0.1% or less of the correlation energy for large systems, with only a moderate computational effort. The first iteration of GF2 is an MP2 calculation that is done in linear scaling, hence we obtain an extremely fast stochastic MP2 (sMP2) method as a by-product. While here we consider finite systems with large band gaps where at low temperatures effects are negligible, the sGF2 formalism is temperature dependent and general and can be applied to finite or periodic systems with small gaps at finite temperatures.

We develop a new scheme for determining molecular partial atomic charges (PACs) with external electrostatic potential (ESP) closely mimicking that of the molecule. The PACs are the ‘minimal corrections’ to a reference set of PACs necessaryfor reproducing exactly the tensor components of the Cartesian zero-, first- and second- molecular electrostatic multipoles. We evaluate the quality of ESP reproduction when ‘minimally correcting’(MC) Mulliken, Hirshfeld or iterative-Hirshfeld reference PACs. In all these cases, the MC-PACs significantly improve the ESP while preserving the reference PACs’invariance under the molecular symmetry operations. When iterative-Hirshfeld PACs are used as reference, the MC-PACs yield ESPs of comparable quality to those of the ChElPG charge fitting method.

An ab initio Langevin dynamics approach is developed based on stochastic density functional theory (sDFT) within a new embedded fragment formalism. The forces on the nuclei generated by sDFT contain a random component natural to Langevin dynamics and its standard deviation is used to estimate the friction term on each atom by satisfying the fluctuation–dissipation relation. The overall approach scales linearly with system size even if the density matrix is not local and is thus applicable to ordered as well as disordered extended systems. We implement the approach for a series of silicon nanocrystals (NCs) of varying size with a diameter of up to 3nm corresponding to Ne = 3000 electrons and generate a set of configurations that are distributed canonically at a fixed temperature, ranging from cryogenic to room temperature. We also analyze the structure properties of the NCs and discuss the reconstruction of the surface geometry.

Charge carrier localization in extended atomic systems has been described previously as being driven by disorder, point defects, or distortions of the ionic lattice. Here we show for the first time by means of first-principles computations that charge carriers can spontaneously localize due to a purely electronic effect in otherwise perfectly ordered structures. Optimally tuned range-separated density functional theory and many-body perturbation calculations within the GW approximation reveal that in trans-polyacetylene and polythiophene the hole density localizes on a length scale of several nanometers. This is due to exchange-induced translational symmetry breaking of the charge density. Ionization potentials, optical absorption peaks, excitonic binding energies, and the optimally tuned range parameter itself all become independent of polymer length as it exceeds the critical localization length. Moreover, we find that lattice disorder and the formation of a polaron result from the charge localization in contrast to the traditional view that lattice distortions precede charge localization. Our results can explain experimental findings that polarons in conjugated polymers form instantaneously after exposure to ultrafast light pulses.

Multiexciton generation, by which more than a single electron–hole pair is generated on optical excitation, is a promising paradigm for pushing the efficiency of solar cells beyond the Shockley–Queisser limit of 31%. Utilizing this paradigm, however, requires the onset energy of multiexciton generation to be close to twice the band gap energy and the efficiency to increase rapidly above this onset. This challenge remains unattainable even using confined nanocrystals, nanorods or nanowires. Here, we show how both goals can be achieved in a nanorod heterostructure with type-II band offsets. Using pseudopotential atomistic calculation on a model type-II semiconductor heterostructure we predict the optimal conditions for controlling multiexciton generation efficiencies at twice the band gap energy. For a finite band offset, this requires a sharp interface along with a reduction of the exciton cooling and may enable a route for breaking the Shockley–Queisser limit.

Despite its importance, it is still not known by which mechanism Li-doped MgO catalyzes the oxidative coupling of methane to ethane. Nevertheless, it is commonly assumed that the mechanism goes through catalytic H abstraction from methane via a Li+O– surface defect. In this paper we use first-principles density functional theory calculations to show that the reaction is significantly more exothermic when the Li+O– defect is situated on a step edge instead of on the flat surface. We find that the reaction on the step is exothermic by 0.25 eV, whereas it is endothermic by 0.3 eV on the flat surface. The presence of the Li dopant in the step edge is crucial for the exothermicity of the reaction. These findings suggest that surface steps which include lithium defects could be responsible for the catalytic behavior of Li/MgO. Following the binding of hydrogen to the Li+O– defect on the step edge the methyl radical can either depart to the gas phase or bind to an adjacent step-edge oxygen atom, increasing the exothermicity of the overall process to 0.8 eV. Activation energies of 0.2 eV for the first pathway and 0.5–0.8 eV for the second were calculated.

In exact density functional theory, the total ground-state energy is a series of linear segments between integer electron points, a condition known as “piecewise linearity.” Deviation from this condition is indicative of poor predictive capabilities for electronic structure, in particular of ionization energies, fundamental gaps, and charge transfer. In this article, we take a new look at the deviation from linearity (i.e., curvature) in the solid-state limit by considering two different ways of approaching it: a large finite system of increasing size and a crystal represented by an increasingly large reference cell with periodic boundary conditions. We show that the curvature approaches vanishing values in both limits, even for functionals which yield poor predictions of electronic structure, and therefore cannot be used as a diagnostic or constructive tool in solids. We find that the approach towards zero curvature is different in each of the two limits, owing to the presence of a compensating background charge in the periodic case. Based on these findings, we present a new criterion for functional construction and evaluation, derived from the size-dependence of the curvature, along with a practical method for evaluating this criterion. For large finite systems, we further show that the curvature is dominated by the self-interaction of the highest occupied eigenstate. These findings are illustrated by computational studies of various solids, semiconductor nanocrystals, and long alkane chains.

The applicability of Mulliken’s theory for photoinduced as well as intramolecular charge-transfer states is examined for several systems of interest by comparing its predictions to TDDFT excitation energies, obtained using functionals appropriate for charge-transfer (CT) states. The results show that it is possible to estimate the energy of the CT state of a donor–acceptor pair on the basis of information on the separate donor and acceptor moieties, along with structural data, within 0.3 eV of TDDFT values. The novelty and usefulness of the proposed method lies mainly in PET applications where the TDDFT determination of the CT state is challenging.

We develop a stochastic formulation of the optimally tuned range-separated hybrid density functional theory that enables significant reduction of the computational effort and scaling of the nonlocal exchange operator at the price of introducing a controllable statistical error. Our method is based on stochastic representations of the Coulomb convolution integral and of the generalized Kohn–Sham density matrix. The computational cost of the approach is similar to that of usual Kohn–Sham density functional theory, yet it provides a much more accurate description of the quasiparticle energies for the frontier orbitals. This is illustrated for a series of silicon nanocrystals up to sizes exceeding 3000 electrons. Comparison with the stochastic GW many-body perturbation technique indicates excellent agreement for the fundamental band gap energies, good agreement for the band edge quasiparticle excitations, and very low statistical errors in the total energy for large systems. The present approach has a major advantage over one-shot GW by providing a self-consistent Hamiltonian that is central for additional postprocessing, for example, in the stochastic Bethe–Salpeter approach.

A time-dependent formulation for electron-hole excitations in extended finite systems, based on the Bethe-Salpeter equation (BSE), is developed using a stochastic wave function approach. The time-dependent formulation builds on the connection between time-dependent Hartree-Fock (TDHF) theory and the configuration-interaction with single substitution (CIS) method. This results in a time-dependent Schrödinger-like equation for the quasiparticle orbital dynamics based on an effective Hamiltonian containing direct Hartree and screened exchange terms, where screening is described within the random-phase approximation (RPA). To solve for the optical-absorption spectrum, we develop a stochastic formulation in which the quasiparticle orbitals are replaced by stochastic orbitals to evaluate the direct and exchange terms in the Hamiltonian as well as the RPA screening. This leads to an overall quadratic scaling, a significant improvement over the equivalent symplectic eigenvalue representation of the BSE. Application of the time-dependent stochastic BSE (TDsBSE) approach to silicon and CdSe nanocrystals up to size of  3000 electrons is presented and discussed.

We develop a new smoothing or extrapolating method, based on discrete Laguerre functions, for systematically analyzing the stochastic signal of shifted-contour auxiliary-field Monte Carlo. We study the statistical errors and extrapolation errors using full configuration-interaction energies for the doubly stretched water molecule. The only free parameter is the order N of the fit. We show that low N emphasizes stability while higher N enable improved extrapolation, at the cost of increased statistical errors. Typically, one should use low order for signals based on a small number of iterations while higher order is efficacious for signals based on large number of iterations. We provide a heuristic algorithm for determining the order to be used and show its utility.

A stochastic approach to time-dependent density functional theory is developed for computing the absorption cross section and the random phase approximation (RPA) correlation energy. The core idea of the approach involves time-propagation of a small set of stochastic orbitals which are first projected on the occupied space and then propagated in time according to the time-dependent Kohn-Sham equations. The evolving electron density is exactly represented when the number of random orbitals is infinite, but even a small number ( 16) of such orbitals is enough to obtain meaningful results for absorption spectrum and the RPA correlation energy per electron. We implement the approach for silicon nanocrystals using real-space grids and find that the overall scaling of the algorithm is sublinear with computational time and memory.

We develop a formalism to calculate the quasiparticle energy within the GW many-body perturbation correction to the density functional theory. The occupied and virtual orbitals of the Kohn-Sham Hamiltonian are replaced by stochastic orbitals used to evaluate the Green function G, the polarization potential W, and, thereby, the GW self-energy. The stochastic GW (sGW) formalism relies on novel theoretical concepts such as stochastic time-dependent Hartree propagation, stochastic matrix compression, and spatial or temporal stochastic decoupling techniques. Beyond the theoretical interest, the formalism enables linear scaling GW calculations breaking the theoretical scaling limit for GW as well as circumventing the need for energy cutoff approximations. We illustrate the method for silicon nanocrystals of varying sizes with Ne > 3000 electrons.

We develop a method in which the electronic densities of small fragments determined by Kohn-Sham density functional theory (DFT) are embedded using stochastic DFT to form the exact density of the full system. The new method preserves the scaling and the simplicity of the stochastic DFT but cures the slow convergence that occurs when weakly coupled subsystems are treated. It overcomes the spurious charge fluctuations that impair the applications of the original stochastic DFT approach. We demonstrate the new approach on a fullerene dimer and on clusters of water molecules and show that the density of states and the total energy can be accurately described with a relatively small number of stochastic orbitals.

Using a combination of density functional theory and quantum master equations approach, we study the effect of electromagnetic (EM) coupling on the nonequilibrium steady-state behavior of a recently introduced gated molecular junction. This junction was demonstrated in a previous publication to exhibit sharp current switching near a certain critical DC field Ez*, which induces intramolecular charge transfer, and here, we analyze the steady-state population and current when an AC EM field (EMF) is present. The AC EMF at frequency $ømega_0$ produces pronounced population and current features at gate fields Ez = Ez* ± $\hbar ømega_0/ez$ (where $e_z$ is the dipole of the charge-transfer state) and thus allows additional sharp switching capability at lower gate fields. We found that even when EMF is absent, the EM coupling itself changes the overall steady-state population and current distributions because it allows for relaxation via spontaneous emission

Despite its importance, it is still not known by which mechanism Li-doped MgO catalyzes the oxidative coupling of methane to ethane. Nevertheless, it is commonly assumed that the mechanism goes through catalytic H abstraction from methane via a Li+O– surface defect. In this paper we use first-principles density functional theory calculations to show that the reaction is significantly more exothermic when the Li+O– defect is situated on a step edge instead of on the flat surface. We find that the reaction on the step is exothermic by 0.25 eV, whereas it is endothermic by 0.3 eV on the flat surface. The presence of the Li dopant in the step edge is crucial for the exothermicity of the reaction. These findings suggest that surface steps which include lithium defects could be responsible for the catalytic behavior of Li/MgO. Following the binding of hydrogen to the Li+O– defect on the step edge the methyl radical can either depart to the gas phase or bind to an adjacent step-edge oxygen atom, increasing the exothermicity of the overall process to 0.8 eV. Activation energies of 0.2 eV for the first pathway and 0.5–0.8 eV for the second were calculated.