Stochastic orbital techniques offer reduced computational scaling and memory requirements to describe ground and excited states at the cost of introducing controlled statistical errors. Such techniques often rely on two basic operations, stochastic trace estimation and stochastic resolution of identity, both of which lead to statistical errors that scale with the number of stochastic realizations (\$N\_\\textbackslashxi\\$) as \$\textbackslashsqrt\N\_\\textbackslashxi\ˆ\-1\\\$. Reducing the statistical errors without significantly increasing \$N\_\\textbackslashxi\\$ has been challenging and is central to the development of efficient and accurate stochastic algorithms. In this work, we build upon recent progress made to improve stochastic trace estimation based on the ubiquitous Hutchinson's algorithm and propose a two-step approach for the stochastic resolution of identity, in the spirit of the Hutch++ method. Our approach is based on employing a randomized low-rank approximation followed by a residual calculation, resulting in statistical errors that scale much better than \$\textbackslashsqrt\N\_\\textbackslashxi\ˆ\-1\\\$. We implement the approach within the second-order Born approximation for the self-energy in the computation of neutral excitations and discuss three different low-rank approximations for the two-body Coulomb integrals. Tests on a series of hydrogen dimer chains with varying lengths demonstrate that the Hutch++-like approximations are computationally more efficient than both deterministic and purely stochastic (Hutchinson) approaches for low error thresholds and intermediate system sizes. Notably, for arbitrarily large systems, the Hutchinson-like approximation outperforms both deterministic and Hutch++-like methods.

We report high-level calculations of the excited states of [2,2]-paracyclophane (PCP), which was recently investigated experimentally by ultrafast pump-probe experiments on oriented single crystals [Haggag et al., ChemPhotoChem 6 e202200181 (2022)]. PCP, in which the orientation of the two benzene rings and their range of motion are constrained, serves as a model for studying benzene exciplex formation. The character of the excimer state and the state responsible for the brightest transition are similar to those in benzene dimer. The constrained structure of PCP allows one to focus on the most important degree of freedom, the inter-ring distance. The calculations explain the main features of the transient absorption spectral evolution. This brightest transition of the excimer is polarized along the inter-fragment axis. The absorption of light polarized in the plane of the rings reveals the presence of other absorbing states of Rydberg character, with much weaker intensities. We also report new transient absorption data obtained by a broadband 8 fs pump, which time-resolve strong modulations of the excimer absorption. The combination of theory and experiment provides a detailed picture of the evolution of the electronic structure of the PCP excimer in the course of a single molecular vibration.

Electron-transfer reactions are ubiquitous in chemistry and biology. The electrons quantum nature allows its transfer across long distances. In the well-known harpoon mechanism, electron-transfer results in Coulombic attraction between initially neutral reactants that leads to dramatic increase in the reaction rate. Here we present a different mechanism, in which electron-transfer from a neutral reactant to a multiply charged cation results in strong repulsion that encodes the electron-transfer distance in the kinetic energy release. 3D coincidence-imaging allows to identify such “inverse” harpoon products, predicted by non adiabatic molecular dynamics simulations to occur between H2 and HCOH2+ following double-ionization of isolated methanol molecules. Detailed comparison of measured and simulated data indicates that while the relative probability of long-range electron-transfer events is correctly predicted, theory overestimates the electron-transfer distance.

We consider an arbitrary quantum mechanical system, initially in its ground-state, exposed to a time-dependent electromagnetic pulse with a carrier frequency ω₀ and a slowly varying envelope of finite duration. By working out a solution to the time-dependent Schrödinger equation in the high-ω₀ limit, we find that, to the leading order in ω−10, a perfect self-cancellation of the system’s linear response occurs as the pulse switches off. Surprisingly, the system’s observables are, nonetheless, describable in terms of a combination of its linear density response function and nonlinear functions of the electric field. An analysis of a jellium slab and jellium sphere models reveals a very high surface sensitivity of the considered setup, producing a richer excitation spectrum than accessible within the conventional linear response regime. On this basis, we propose a new spectroscopic technique, which we provisionally name the Nonlinear High-Frequency Pulsed Spectroscopy (NLHFPS). Combining the advantages of the extraordinary surface sensitivity, the absence of constraints by the traditional dipole selection rules, and the clarity of theoretical interpretation utilizing the linear response time-dependent density functional theory, NLHFPS has a potential to evolve into a powerful characterization method for nanoscience and nanotechnology

The trihydrogen ion has a central role in creating complex molecules in the interstellar medium. Therefore, its formation and destruction mechanisms in high photon energy environments involving organic molecules are drawing significant experimental and theoretical attention. Here, we employ a combination of time-resolved ultrafast extreme-ultraviolet pump and near-infrared probe spectroscopy applied to the deuterated CH3OD methanol molecule. Similar to other double-ionization studies, the isotopic labeling reveals two competing pathways for forming trihydrogen: A) H+3 + COH+ and B) H+3 + HCO+. We validate our high-level ab initio nonadiabatic molecular dynamic simulations by showing that it closely reproduces the essential features of the measured kinetic energy release distribution and branching ratios of the two pathways of the deuterated system. The success of ab initio simulation in describing single photon double-ionization allows for an unprecedented peek into the formation pathways for the undeuterated species, beyond present experimental reach. For this case, we find that the kinetic energy release of pathway B shifts to lower energies by more than 0.6 eV due to a dynamical isotope effect. We also determine the mechanism for trihydrogen formation from excited states of the dication and elucidate the isotope effect’s role in the observed dynamics.

Despite the abundance of data concerning single-photon double ionization of methanol, the spin state of the emitted electron pair has never been determined. Here we present the first evidence that identifies the emitted electron pair spin as overwhelmingly singlet when the dication forms in low-energy configurations. The experimental data show that while the yield of the CH2O+ + H3+ Coulomb explosion channel is abundant, the metastable methanol dication is largely absent. According to high-level ab initio simulations, these facts indicate that photoionization promptly forms singlet dication states, where they quickly decompose through various channels, with significant H3+ yields on the low-lying states. In contrast, if we assume that the initial dication is formed in one of the low-lying triplet states, the ab initio simulations exhibit a metastable dication, contradicting the experimental findings. Comparing the average simulated branching ratios with the experimental data suggests a \textgreater3 order of magnitude enhancement of the singlet:triplet ratio compared with their respective 1:3 multiplicities.

Generalized Kohn−Sham density functional theory is a popular computational tool for the ground state of extended systems, particularly within range-separated hybrid (RSH) functionals that capture the long-range electronic interaction. Unfortunately, the heavy computational cost of the nonlocal exchange operator in RSH-DFT usually confines the approach to systems with at most a few hundred electrons. A significant reduction in the computational cost is achieved by representing the density matrix with stochastic orbitals and a stochastic decomposition of the Coulomb convolution (J. Phys. Chem. A 2016, 120, 3071). Here, we extend the stochastic RSH approach to excited states within the framework of linear-response generalized Kohn−Sham time-dependent density functional theory (GKS-TDDFT) based on the plane-wave basis. As a validation of the stochastic GKS-TDDFT method, the excitation energies of small molecules N2 and CO are calculated and compared to the deterministic results. The computational efficiency of the stochastic method is demonstrated with a two-dimensional MoS2 sheet (∼1500 electrons), whose excitation energy, exciton charge density, and (excited state) geometric relaxation are determined in the absence and presence of a point defect.

We develop a stochastic resolution of identity approach to the real-time second-order Green’s function (real-time sRI-GF2) theory, extending our recent work for imaginary-time Matsubara Green’s function [Takeshita et al. J. Chem. Phys. 2019, 151, 044114]. The approach provides a framework to obtain the quasi-particle spectra across a wide range of frequencies and predicts ionization potentials and electron affinities. To assess the accuracy of the real-time sRI-GF2, we study a series of molecules and compare our results to experiments as well as to a many-body perturbation approach based on the GW approximation, where we find that the real-time sRI-GF2 is as accurate as self-consistent GW. The stochastic formulation reduces the formal computatinal scaling from O(Ne5) down to O(Ne3) where Ne is the number of electrons. This is illustrated for a chain of hydrogen dimers, where we observe a slightly lower than cubic scaling for systems containing up to Ne ≈ 1000 electrons.

Abstract The Kubo-Greenwood (KG) formula is often used in conjunction with Kohn-Sham (KS) density functional theory (DFT) to compute the optical conductivity, particularly for warm dense mater. For applying the KG formula, all KS eigenstates and eigenvalues up to an energy cutoff are required and thus the approach becomes expensive, especially for high temperatures and large systems, scaling cubically with both system size and temperature. Here, we develop an approach to calculate the KS conductivity within the stochastic DFT (sDFT) framework, which requires knowledge only of the KS Hamiltonian but not its eigenstates and values. We show that the computational effort associated with the method scales linearly with system size and reduces in proportion to the temperature unlike the cubic increase with traditional deterministic approaches. In addition, we find that the method allows an accurate description of the entire spectrum, including the high-frequency range, unlike the deterministic method which is compelled to introduce a high-frequency cut-off due to memory and computational time constraints. We apply the method to helium-hydrogen mixtures in the warm dense matter regime at temperatures of \textbackslashsim60\textbackslashtext\kK\ and find that the system displays two conductivity phases, where a transition from non-metal to metal occurs when hydrogen atoms constitute \textbackslashsim0.3 of the total atoms in the system.

The conjecture that, as in bulk semiconductors, hot multiexcitons in nanocrystals cool rapidly to the lowest available energy levels is tested here by recording the effects of a single cold “spectator” exciton on the relaxation dynamics of a subsequently deposited hot counterpart. Results in CdSe/CdS nanodots show that a preexisting cold “spectator exciton” allows only half of the photoexcited electrons to relax directly to the band-edge. The rest are blocked in an excited quantum state due to conflicts in spin orientation. The latter fully relax in this sample only after ∼25 ps as the blocked electrons spins flip, prolonging the temporal window of opportunity for harvesting the retained energy more than 100 fold! Common to all quantum-confined nanocrystals, this process will delay cooling and impact the spectroscopic signatures of hot multiexcitons in all envisioned generation scenarios. How the spin-flipping rate scales with particle size and temperature remains to be determined.

A multifaceted agreement between ab initio theoretical predictions and experimental measurements, including branching ratios, channel-specific kinetic energy release, and three-body momentum correlation spectra, leads to the identification of new mechanisms in Coulomb-explosion (CE) induced two- and three-body breakup processes in methanol. These identified mechanisms include direct nonadiabatic Coulomb explosion responsible for CO bond-breaking, a long-range “ inverse harpooning” dominating the production of H2+ + HCOH+, a transient proton migration leading to surprising energy partitioning in three-body fragmentation and other complex dynamics forming products such as H2O+ and H3+. These mechanisms provide general concepts that should be useful for analyzing future time-resolved Coulomb explosion imaging of methanol as well as other molecular systems. These advances are enabled by a combination of recently developed experimental and computational techniques, using weak ultrafast EUV pulses to initiate the CE and a high-level quantum chemistry approach to follow the resulting field-free nonadiabatic molecular dynamics.

Absorption cross-section spectra for gold nanoparticles were calculated using fully quantum Stochastic Density Functional Theory and a classical Finite-Difference Time Domain Maxwell solver. Spectral shifts were monitored as a function of size (1.3–) and shape (octahedron, cubeoctahedron and truncated cube). Even though the classical approach is forced to fit the quantum time-dependent density functional theory at , at smaller sizes there is a significant deviation as the classical theory is unable to account for peak splitting and spectral blueshifts even after quantum spectral corrections. We attribute the failure of classical methods at predicting these features to quantum effects and low density of states in small nanoparticles. Classically, plasmon resonances are modelled as collective conduction electron excitations, but at small nanoparticle size these excitations transition to few or even individual conductive electron excitations, as indicated by our results.

Single-photon Coulomb explosion of methanol is instigated using the broad bandwidth pulse achieved through high-order harmonics generation. Using 3D coincidence fragment imaging of one molecule at a time, the kinetic energy release (KER) and angular distributions of the products are measured in different Coulomb explosion (CE) channels. Two-body CE channels breaking either the C–O or the C–H bonds are described as well as a proton migration channel forming H2O+, which is shown to exhibit higher KER. The results are compared to intense-field Coulomb explosion measurements in the literature. The interpretation of broad bandwidth single-photon CE data is discussed and supported by ab initio calculations of the predominant C–O bond breaking CE channel. We discuss the importance of these findings for achieving time resolved imaging of ultrafast dynamics.

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.

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

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.