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.

Engineered metallic nanoparticles, which are found in numerous applications, are usually stabilized by organic ligands influencing their interfacial properties. We found that the ligands affect tremendously the electrochemical peak oxidation potentials of the nanoparticles. In this work, identical gold nanoparticles were ligand-exchanged and carefully analyzed to enable a precise and highly reproducible comparison. The peak potential difference between gold nanoparticles stabilized by various ligands, such as 2- and 4-mercaptobenzoic acid, can be as high as 71 mV, which is substantial in energetic terms. A detailed study supported by density functional theory (DFT) calculations aimed to determine the source of this interesting effect. The DFT simulations of the ligand adsorption modes on Au surfaces were used to calculate the redox potentials through the thermodynamic cycle method. The DFT results of the peak potential shift were in good agreement with the experimental results for a few ligands, but showed some discrepancy, which was attributed to kinetic effects. The kinetic rate constant of the oxidation of Au nanoparticles stabilized by 4-mercaptobenzoic acid was found to be twice as large as that of the Au nanoparticles stabilized by citrate, as calculated from Laviron’s theory and the Tafel equation. Finally, these findings could be applied to some novel applications such as determining the distribution of nanoparticle population in a dispersion as well as monitoring the ligand exchange between nanoparticles.

This study overviews and extends a recently developed stochastic finite-temperature Kohn-Sham density functional theory to study warm dense matter using Langevin dynamics, specifically under periodic boundary conditions. The method's algorithmic complexity exhibits nearly linear scaling with system size and is inversely proportional to the temperature. Additionally, a novel linear-scaling stochastic approach is introduced to assess the Kubo-Greenwood conductivity, demonstrating exceptional stability for DC conductivity. Utilizing the developed tools, we investigate the equation of state, radial distribution, and electronic conductivity of Hydrogen at a temperature of 30,000K. As for the radial distribution functions, we reveal a transition of Hydrogen from gas-like to liquid-like behavior as its density exceeds 4 g/cm³. As for the electronic conductivity as a function of the density, we identified a remarkable isosbestic point at frequencies around 7eV, which may be an additional signature of a gas-liquid transition in Hydrogen at 30,000K.