We examine the possibility of using a Metropolis algorithm for computing the exchange energy in a large molecular system. Following ideas set forth in a recent publication (Baer, Neuhauser, and Rabani, Phys. Rev. Lett. 111, 106402 (2013)) we focus on obtaining the exchange energy per particle (ExPE, as opposed to the total exchange energy) to a predefined statistical error and on determining the numerical scaling of the calculation achieving this. For this we assume that the occupied molecular orbitals (MOs) are known and given in terms of a standard Gaussian atomic basis set. The Metropolis random walk produces a sequence of pairs of three-dimensional points (x,x'), which are distributed in proportion to $\rho(x,x')^2$, where $\rho(x,x')$ is the density matrix. The exchange energy per particle is then simply the average of the Coulomb repulsion energy U_C(|x–x'|) over these pairs. To reduce the statistical error we separate the exchange energy into a short-range term that can be calculated deterministically in a linear scaling fashion and a long-range term that is treated by the Metropolis method. We demonstrate the method on water clusters and silicon nanocrystals showing the magnitude of the ExPE standard deviation is independent of system size. In the water clusters a longer random walk was necessary to obtain full ergodicity as Metropolis walkers tended to get stuck for a while in localized regions. We developed a diagnostic tool that can alert a user when such a situation occurs. The calculation effort scales linearly with system size if one uses an atom screening procedure that can be made numerically exact. In systems where the MOs can be localized efficiently the ExPE can even be computed with “sublinear scaling” as the MOs themselves can be screened.

Density functional theory with optimally tuned range-separated hybrid (OT-RSH) functionals has been recently suggested [Refaely-Abramson et al. Phys. Rev. Lett. 2012, 109, 226405] as a nonempirical approach to predict the outer-valence electronic structure of molecules with the same accuracy as many-body perturbation theory. Here, we provide a quantitative evaluation of the OT-RSH approach by examining its performance in predicting the outer-valence electron spectra of several prototypical gas-phase molecules, from aromatic rings (benzene, pyridine, and pyrimidine) to more complex organic systems (terpyrimidinethiol and copper phthalocyanine). For a range up to several electronvolts away from the frontier orbital energies, we find that the outer-valence electronic structure obtained from the OT-RSH method agrees very well (typically within  0.1–0.2 eV) with both experimental photoemission and theoretical many-body perturbation theory data in the GW approximation. In particular, we find that with new strategies for an optimal choice of the short-range fraction of Fock exchange, the OT-RSH approach offers a balanced description of localized and delocalized states. We discuss in detail the sole exception found—a high-symmetry orbital, particular to small aromatic rings, which is relatively deep inside the valence state manifold. Overall, the OT-RSH method is an accurate DFT-based method for outer-valence electronic structure prediction for such systems and is of essentially the same level of accuracy as contemporary GW approaches, at a reduced computational cost.

We study the role of the effective mass, band mixing, and phonon emission on multiexciton generation in IV–VI nanocrystals. A four-band k · p effective mass model, which allows for an independent variation of these parameters, is adopted to describe the electronic structure of the nanocrystals. Multiexciton generation efficiencies are calculated using a Green’s function formalism, providing results that are numerically similar to impact excitation. We find that multiexciton generation efficiencies are maximized when the effective mass of the electron and hole are small and similar. Contact with recent experimental results for multiexciton generation in PbS and PbSe is made.

A fast stochastic method for calculating the second order Møller-Plesset (MP2) correction to the correlation energy of large systems of electrons is presented. The approach is based on reducing the exact summation over occupied and unoccupied states to a time-dependent trace formula amenable to stochastic sampling. We demonstrate the abilities of the method to treat systems with thousands of electrons using hydrogen passivated silicon spherical nanocrystals represented on a real space grid, much beyond the capabilities of present day MP2 implementations.

A fast method is developed for calculating the random phase approximation (RPA) correlation energy for density functional theory. The correlation energy is given by a trace over a projected RPA response matrix, and the trace is taken by a stochastic approach using random perturbation vectors. For a fixed statistical error in the total energy per electron, the method scales, at most, quadratically with the system size; however, in practice, due to self-averaging, it requires less statistical sampling as the system grows, and the performance is close to linear scaling. We demonstrate the method by calculating the RPA correlation energy for cadmium selenide and silicon nanocrystals with over 1500 electrons. We find that the RPA correlation energies per electron are largely independent of the nanocrystal size. In addition, we show that a correlated sampling technique enables calculation of the energy difference between two slightly distorted configurations with scaling and a statistical error similar to that of the total energy per electron.

Fundamental gap renormalization due to electronic polarization is a basic phenomenon in molecular crystals. Despite its ubiquity and importance, all conventional approaches within density-functional theory completely fail to capture it, even qualitatively. Here, we present a new screened range-separated hybrid functional, which, through judicious introduction of the scalar dielectric constant, quantitatively captures polarization-induced gap renormalization, as demonstrated on the prototypical organic molecular crystals of benzene, pentacene, and C60. This functional is predictive, as it contains system-specific adjustable parameters that are determined from first principles, rather than from empirical considerations.

A recently introduced molecular junction, for which the gate acts as an on/off switch for intrajunction electron transfer between localized donor and acceptor sites is studied. We demonstrate that a Landauer + density functional (DFT) approach is fundamentally flawed for describing the electronic conductance in this system. By comparing the Landauer + DFT conductance to that predicted by the Redfield quantum master equations, we point out several effects that cannot be explained by the former approach. The molecular junction is unique in the small number of conductance channels and their sharp response to the gate.

We develop an alternative formulation in the energy-domain to calculate the second order Møller–Plesset (MP2) perturbation energies. The approach is based on repeatedly choosing four random energies using a nonseparable guiding function, filtering four random orbitals at these energies, and averaging the resulting Coulomb matrix elements to obtain a statistical estimate of the MP2 correlation energy. In contrast to our time-domain formulation, the present approach is useful for both quantum chemistry and real-space/plane wave basis sets. The scaling of the MP2 calculation is roughly linear with system size, providing a useful tool to study dispersion energies in large systems. This is demonstrated on a structure of 64 fullerenes within the SZ basis as well as on silicon nanocrystals using real-space grids.

We study how shape affects multiexciton generation rates in a semiconducting nanocrystal by considering CdSe nanorods with varying diameters and aspect ratios. The calculations employ an atomistic semiempirical pseudopotential model combined with an efficacious stochastic approach applied to systems containing up to 20 000 atoms. The effect of nanorod diameter and aspect ratio on multiexciton generation rates is analyzed in terms of the scaling of the density of trion states and the scaling of the Coulomb couplings. Both show distinct scaling from spherical nanocrystals leading to a surprising result where the multiexciton generation rates are roughly independent of the nanorod length.

We formulate the Kohn-Sham density functional theory (KS-DFT) as a statistical theory in which the electron density is determined from an average of correlated stochastic densities in a trace formula. The key idea is that it is sufficient to converge the total energy per electron to within a predefined statistical error in order to obtain reliable estimates of the electronic band structure, the forces on nuclei, the density and its moments, etc. The fluctuations in the total energy per electron are guaranteed to decay to zero as the system size increases. This facilitates “self-averaging” which leads to the first ever report of sublinear scaling KS-DFT electronic structure. The approach sidesteps calculation of the density matrix and thus, is insensitive to its evasive sparseness, as demonstrated here for silicon nanocrystals. The formalism is not only appealing in terms of its promise to far push the limits of application of KS-DFT, but also represents a cognitive change in the way we think of electronic structure calculations as this stochastic theory seamlessly converges to the thermodynamic limit.

Perdew et al. discovered two different properties of exact Kohn–Sham density functional theory (DFT): (i) The exact total energy versus particle number is a series of linear segments between integer electron points. (ii) Across an integer number of electrons, the exchange-correlation potential “jumps” by a constant, known as the derivative discontinuity (DD). Here we show analytically that in both the original and the generalized Kohn–Sham formulation of DFT the two properties are two sides of the same coin. The absence of a DD dictates deviation from piecewise linearity, but the latter, appearing as curvature, can be used to correct for the former, thereby restoring the physical meaning of orbital energies. A simple correction scheme for any semilocal and hybrid functional, even Hartree–Fock theory, is shown to be effective on a set of small molecules, suggesting a practical correction for the infamous DFT gap problem. We show that optimally tuned range-separated hybrid functionals can inherently minimize both DD and curvature, thus requiring no correction, and that this can be used as a sound theoretical basis for novel tuning strategies.

The range-separated Baer-Neuhauser-Livshits functional with optimized range-separation parameter \gamma was employed to predict ionization energies of alkanes and oligothiophenes. For all systems negative orbital energies of neutral species are consistent with explicitly calculated states of cations. For \sigma-systems excellent agreement with experiment is obtained while for conjugated π-systems IPs are underestimated.

We present a method for obtaining outer-valence quasiparticle excitation energies from a density-functional-theory-based calculation, with an accuracy that is comparable to that of many-body perturbation theory within the GW approximation. The approach uses a range-separated hybrid density functional, with an asymptotically exact and short-range fractional Fock exchange. The functional contains two parameters, the range separation and the short-range Fock fraction. Both are determined nonempirically, per system, on the basis of the satisfaction of exact physical constraints for the ionization potential and frontier-orbital many-electron self-interaction, respectively. The accuracy of the method is demonstrated on four important benchmark organic molecules: perylene, pentacene, 3,4,9,10-perylene-tetracarboxylic-dianydride (PTCDA), and 1,4,5,8-naphthalene-tetracarboxylic-dianhydride (NTCDA). We envision that for the outer-valence excitation spectra of finite systems the approach could provide an inexpensive alternative to GW, opening the door to the study of presently out of reach large-scale systems.

Excitation gaps are of considerable significance in electronic structure theory. Two different gaps are of particular interest. The fundamental gap is defined by charged excitations, as the difference between the first ionization potential and the first electron affinity. The optical gap is defined by a neutral excitation, as the difference between the energies of the lowest dipole-allowed excited state and the ground state. Within many-body perturbation theory, the fundamental gap is the difference between the corresponding lowest quasi-hole and quasi-electron excitation energies, and the optical gap is addressed by including the interaction between a quasi-electron and a quasi-hole. A long-standing challenge has been the attainment of a similar description within density functional theory (DFT), with much debate on whether this is an achievable goal even in principle. Recently, we have constructed and applied a new approach to this problem. Anchored in the rigorous theoretical framework of the generalized Kohn–Sham equation, our method is based on a range-split hybrid functional that uses exact long-range exchange. Its main novel feature is that the range-splitting parameter is not a universal constant but rather is determined from first principles, per system, based on satisfaction of the ionization potential theorem. For finite-sized objects, this DFT approach mimics successfully, to the best of our knowledge for the first time, the quasi-particle picture of many-body theory. Specifically, it allows for the extraction of both the fundamental and the optical gap from one underlying functional, based on the HOMO–LUMO gap of a ground-state DFT calculation and the lowest excitation energy of a linear-response time-dependent DFT calculation, respectively. In particular, it produces the correct optical gap for the difficult case of charge-transfer and charge-transfer-like scenarios, where conventional functionals are known to fail. In this perspective, we overview the formal and practical challenges associated with gap calculations, explain our new approach and how it overcomes previous difficulties, and survey its application to a variety of systems.

In conventional spectroscopy, transitions between electronic levels are governed by the electric dipole selection rule because electric quadrupole, magnetic dipole, and coupled electric dipole-magnetic dipole transitions are forbidden in a far field. We demonstrated that by using nanostructured electromagnetic fields, the selection rules of absorption spectroscopy could be fundamentally manipulated. We also show that forbidden transitions between discrete quantum levels in a semiconductor nanorod structure are allowed within the near-field of a noble metal nanoparticle. Atomistic simulations analyzed by an effective mass model reveal the breakdown of the dipolar selection rules where quadrupole and octupole transitions are allowed. Our demonstration could be generalized to the use of nanostructured near-fields for enhancing light-matter interactions that are typically weak or forbidden.

An ab initio variational grand-canonical electronic structure mean-field method, based on the Gibbs–Peierls–Bogoliubov minimum principle for the Gibbs free energy, is applied to the di-lithium (Li+Li) system at temperatures around T \approx 10,000 K and electronic chemical potential of μ \approx -0.1Eh. The method is an extension of the Hartree–Fock approach to finite temperatures. We first study the Li2 molecule at a frozen inter-nuclear distance of R = 3 \AA as a function of temperature. The mean-field electronic structure changes smoothly as temperature increases, up to 104 K, where a sharp spontaneous spin-polarization emerges as the variational mean-field solution. Further increase in the temperature extinguishes this polarization. We analyze the mean-field behavior using a correlated single-site Hubbard model and show it arises from an attempt of the mean-field to mimic the polarization of the spin–spin correlation function of the exact solution. Next, we keep constant the temperature at 104 K and examine the electronic structure as a function of inter-nuclear distance R. At R = 3.7 \AA, a crossing between two free energy states occurs: One state is “spin-unpolarized” (becomes lower in energy when R \ge 3.7 \AA), while the other is “spin polarized”. This crossing causes near-discontinuous jumps in calculated properties of the system and is associated with using the noninteracting electron character of our mean-field approach. Such problems will likely plague FT-DFT calculations as well. We use second-order perturbation theory (PT2) to study effects of electron correlation on the potential of mean force between the two colliding Li atoms. We find that PT2 correlation free energy at  104 K is larger than at 0 K and tends to restore the spin-polarized state as the lowest free energy solution.

We present a molecular junction composed of a donor (polyacetylene strands) and an acceptor (malononitrile) connected together via a benzene ring and coupled weakly to source and drain electrodes on each side, for which a gate electrode induces intramolecular charge transfer, switching reversibly the character of conductance. Using a new brand of density functional theory, for which orbital energies are similar to the quasiparticle energies, we show that the junction displays a single, gate-tunable differential conductance channel in a wide energy range. The gate field must align parallel to the displacement vector between donors and acceptor to affect their potential difference; for strong enough fields, spontaneous intramolecular electron transfer occurs. This event radically affects conductance, reversing the charge of carriers, enabling a spin-polarized current channel. We discuss the physical principles controlling the operation of the junction and find interplay of quantum interference, charging, Coulomb blockade, and electron-hole binding energy effects. We expect that this switching behavior is a generic property for similar donor-acceptor systems of sufficient stability.

In recent generalized Kohn-Sham (GKS) schemes for density functional theory (DFT) Hartree-Fock type exchange is important. In plane waves and grid approaches the high cost of exchange energy calculations makes these GKS considerably more expensive than Kohn-Sham DFT calculations. We develop a stochastic approach for speeding up the calculation of exchange for large systems. We show that stochastic error per particle does not grow and can even decrease with system size (at a given number of iterations). We discuss several alternative approaches and explain how these ideas can be included in the GKS framework.

A stochastic method is developed to calculate the multiexciton generation (MEG) rates in semiconductor nanocrystals (NCs). The numerical effort scales near-linearly with system size allowing the study of MEG rates up to diameters and exciton energies previously unattainable using atomistic calculations. Illustrations are given for CdSe NCs of sizes and energies relevant to current experimental setups, where direct methods require treatment of over 1011 states. The approach is not limited to the study of MEG and can be applied to calculate other correlated electronic processes.

Forster Resonance Energy Transfer (FRET) between fluorescent proteins (FPs) is widely used to construct fluorescent sensor proteins, to study intracellular protein-protein interactions and to monitor conformational changes in multidomain proteins. Although FRET depends strongly on the orientation of the transition dipole moments (TDMs) of the donor and acceptor fluorophores, this orientation dependence is currently not taken into account in FRET sensor design. Similarly, studies that use FRET to derive structural constrains typically assume a kappa(2) of 2/3 or use the TDM of green fluorescent protein, as this is the only FP for which the TDM has been determined experimentally. Here we used time-dependent density functional theory (TD-DFT) methods to calculate the TDM for a comprehensive list of commonly used fluorescent proteins. The method was validated against higher levels of calculation. Validation with model compounds and the experimentally determined TDM of GFP shows that the TDM is mostly determined by the structure of the pi-conjugated fluorophore and is insensitive to non-conjugated side chains or the protein surrounding. Our calculations not only provide TDM for most of the currently used FPs, but also suggest an empirical rule that can be used to obtain the TDMs for newly developed fluorescent proteins in the future.