Publications by Author: Baer, Roi

Submitted
Haggag, O. S. ; Baer, R. ; Ruhman, S. ; Krylov, A. Revisiting the benzene excimer using [2, 2] paracyclophane model system: Experiment and theory. ChemRxiv 2024 Submitted. Publisher's VersionAbstract

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.

Stochastic density functional theory combined with Langevin dynamics for warm dense matter
Hadad, E. R. ; Roy, A. ; Rabani, E. ; Redmer, R. ; Baer, R. Stochastic density functional theory combined with Langevin dynamics for warm dense matter. arXiv.2401.11336 Submitted. Publisher's VersionAbstract

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.

2024
Weak second-order quantum state diffusion unraveling of the Lindblad master equation
Baer, R. ; Adhikari, S. Weak second-order quantum state diffusion unraveling of the Lindblad master equation. Journal of Chemical Physics 2024, 160, 064107. Publisher's VersionAbstract

Abstract Simulating mixed-state evolution in open quantum systems is crucial for various chemical physics, quantum optics, and computer science applications. These simulations typically follow the Lindblad master equation dynamics. An alternative approach known as quantum state diffusion unraveling is based on the trajectories of pure states generated by random wave functions, which evolve according to a nonlinear Itô-Schrödinger equation (ISE). This study introduces weak first- and second-order solvers for the ISE based on directly applying the Itô-Taylor expansion with exact derivatives in the interaction picture. We tested the method on free and driven Morse oscillators coupled to a thermal environment and found that both orders allowed practical estimation with a few dozen iterations. The variance was relatively small compared to the linear unraveling and did not grow with time. The second-order solver delivers much higher accuracy and stability with bigger time steps than the first-order scheme, with a small additional workload. However, the second-order algorithm has quadratic complexity with the number of Lindblad operators as opposed to the linear complexity of the first-order algorithm.

AdhikariWeak2024.pdf
2022
An “inverse” harpoon mechanism
Gope, K. ; Livshits, E. ; Bittner, D. M. ; Baer, R. ; Strasser, D. An “inverse” harpoon mechanism. Science Advances 2022, 8 eabq8084. Publisher's VersionAbstract

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.

Sequential and concerted C-C and C-O bond dissociation in the Coulomb explosion of 2-propanol
Bittner, D. ; Gope, K. ; Livshits, E. ; Baer, R. ; Strasser, D. Sequential and concerted C-C and C-O bond dissociation in the Coulomb explosion of 2-propanol. Journal of Chemical Physics 2022, 157, 074309. Publisher's Version
Linear Weak Scalability of Density Functional Theory Calculations without Imposing Electron Localization
Fabian, M. D. ; Shpiro, B. ; Baer, R. Linear Weak Scalability of Density Functional Theory Calculations without Imposing Electron Localization. J. Chem. Theory Comput. 2022, acs.jctc.1c00829. Publisher's VersionAbstract

Linear scaling density functional theory (DFT) approaches to the electronic structure of materials are often based on the tendency of electrons to localize in large atomic and molecular systems. However, in many cases of actual interest, such as semiconductor nanocrystals, system sizes can reach a substantial extension before significant electron localization sets in, causing a considerable deviation from linear scaling. Herein, we address this class of systems by developing a massively parallel DFT approach which does not rely on electron localization and is formally quadratic scaling yet enables highly efficient linear wall-time complexity in the weak scalability regime. The method extends from the stochastic DFT approach described in Fabian et al. (WIRES: Comp. Mol. Sci. 2019, e1412) but is entirely deterministic. It uses standard quantum chemical atomcentered Gaussian basis sets to represent the electronic wave functions combined with Cartesian real-space grids for some operators and enables a fast solver for the Poisson equation. Our main conclusion is that when a processor-abundant high-performance computing (HPC) infrastructure is available, this type of approach has the potential to allow the study of large systems in regimes where quantum confinement or electron delocalization prevents linear scaling.

Stochastic Vector Techniques in Ground-State Electronic Structure
Baer, R. ; Neuhauser, D. ; Rabani, E. Stochastic Vector Techniques in Ground-State Electronic Structure. Annu. Rev. Phys. Chem. 2022, 73, annurev–physchem–090519–045916. Publisher's VersionAbstract

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 secondorder 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.

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets
Shpiro, B. ; Fabian, M. D. ; Rabani, E. ; Baer, R. Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets. J. Chem. Theory Comput. 2022, 18, 1458–1466. Publisher's VersionAbstract

We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atomcentered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trpzip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (∼1.2 × 10−3Eh/a0) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.

High frequency limit of spectroscopy
Nazarov, V. U. ; Baer, R. High frequency limit of spectroscopy. Journal of Chemical Physics 2022, 157, 084112.Abstract

We consider an arbitrary quantum mechanical system, initially in its ground-state, exposed to a time-dependent electromagnetic pulse with a carrier frequency ω0 and a slowly varying envelope of finite duration. By working out a solution to the time-dependent Schrödinger equation in the high-ω0 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

2021
Tempering stochastic density functional theory
Nguyen, M. ; Li, W. ; Li, B. (Y. ); Baer, R. ; Rabani, E. ; Neuhauser, D. Tempering stochastic density functional theory. J. Chem. Phys. 2021, 5.0063266. Publisher's VersionAbstract

We introduce a tempering approach with stochastic density functional theory (sDFT), labeled t-sDFT, which reduces the statistical errors in the estimates of observable expectation values. This is achieved by rewriting the electronic density as a sum of a "warm" component complemented by "colder" correction(s). Since the "warm" component is larger in magnitude but faster to evaluate, we use many more stochastic orbitals for its evaluation than for the smaller-sized colder correction(s). This results in a significant reduction of the statistical fluctuations and the bias compared to sDFT for the same computational effort. We the method's performance on large hydrogen-passivated silicon nanocrystals (NCs), finding a reduction in the systematic error in the energy by more than an order of magnitude, while the systematic errors in the forces are also quenched. Similarly, the statistical fluctuations are reduced by factors of around 4-5 for the total energy and around 1.5-2 for the forces on the atoms. Since the embedding in t-sDFT is fully stochastic, it is possible to combine t-sDFT with other variants of sDFT such as energy-window sDFT and embedded-fragmented sDFT.

nguyen2021temperin.pdf
Two pathways and an isotope effect in H3+ formation following double ionization of methanol
Gope, K. ; Livshits, E. ; Bittner, D. M. ; Baer, R. ; Strasser, D. Two pathways and an isotope effect in H3+ formation following double ionization of methanol. Natural Sciences 2021, ntls.10022. Publisher's VersionAbstract

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.

ntls.10022.pdf
Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction
Chen, M. ; Baer, R. ; Neuhauser, D. ; Rabani, E. Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction. J. Chem. Phys. 2021, 154, 204108. Publisher's VersionAbstract

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn–Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could√be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1/ Nχ according to the central limit theorem, where Nχ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

chen2021stochastic.pdf
2020
Time-resolving the ultrafast H2 roaming chemistry and H3+ formation using extreme-ultraviolet pulses
Livshits, E. ; Luzon, I. ; Gope, K. ; Baer, R. ; Strasser, D. Time-resolving the ultrafast H2 roaming chemistry and H3+ formation using extreme-ultraviolet pulses. Communications Chemistry 2020, 3 49. Publisher's Version livshits2020time.pdf
Absence of Triplets in Single-Photon Double Ionization of Methanol
Gope, K. ; Livshits, E. ; Bittner, D. M. ; Baer, R. ; Strasser, D. Absence of Triplets in Single-Photon Double Ionization of Methanol. The Journal of Physical Chemistry Letters 2020, 11, 8108–8113. Publisher's VersionAbstract

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.

gope2020absence.pdf
Dopant levels in large nanocrystals using stochastic optimally tuned range-separated hybrid density functional theory
Lee, A. J. ; Chen, M. ; Li, W. ; Neuhauser, D. ; Baer, R. ; Rabani, E. Dopant levels in large nanocrystals using stochastic optimally tuned range-separated hybrid density functional theory. Physical Review B 2020, 102, 035112. Publisher's Version lee2020dopant.pdf
Linear-Response Time-Dependent Density Functional Theory with Stochastic Range-Separated Hybrids
Zhang, X. ; Lu, G. ; Baer, R. ; Rabani, E. ; Neuhauser, D. Linear-Response Time-Dependent Density Functional Theory with Stochastic Range-Separated Hybrids. Journal of Chemical Theory and Computation 2020, 16, 1064–1072. Publisher's VersionAbstract

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.

zhang2020linear.pdf
Range-separated stochastic resolution of identity: Formulation and application to second-order Green’s function theory
Dou, W. ; Chen, M. ; Takeshita, T. Y. ; Baer, R. ; Neuhauser, D. ; Rabani, E. Range-separated stochastic resolution of identity: Formulation and application to second-order Green’s function theory. The Journal of Chemical Physics 2020, 153, 074113. Publisher's VersionAbstract

We develop a range-separated stochastic resolution of identity (RS-SRI) approach for the four-index electron repulsion integrals, where the larger terms (above a predefined threshold) are treated using a deterministic RI and the remaining terms are treated using a SRI. The approach is implemented within a second-order Green’s function formalism with an improved O(N3) scaling with the size of the basis set, N. Moreover, the RS approach greatly reduces the statistical error compared to the full stochastic version [T. Y. Takeshita et al., J. Chem. Phys. 151, 044114 (2019)], resulting in computational speedups of ground and excited state energies of nearly two orders of magnitude, as demonstrated for hydrogen dimer chains and water clusters.

dou2020range.pdf
Efficient Langevin dynamics for "noisy" forces
Arnon, E. ; Rabani, E. ; Neuhauser, D. ; Baer, R. Efficient Langevin dynamics for "noisy" forces. J. Chem. Phys. 2020, 152, 161103. Publisher's VersionAbstract

Efficient Boltzmann-sampling using first-principles methods is challenging for extended systems due to the steep scaling of electronic structure methods with the system size. Stochastic approaches provide a gentler system-size dependency at the cost of introducing "noisy" forces, which serve to limit the efficiency of the sampling. In the first-order Langevin dynamics (FOLD), efficient sampling is achievable by combining a well-chosen preconditioning matrix S with a time-step-bias-mitigating propagator (Mazzola et al., Phys. Rev. Lett., 118, 015703 (2017)). However, when forces are noisy, S is set equal to the force-covariance matrix, a procedure which severely limits the efficiency and the stability of the sampling. Here, we develop a new, general, optimal, and stable sampling approach for FOLD under noisy forces. We apply it for silicon nanocrystals treated with stochastic density functional theory and show efficiency improvements by an order-of-magnitude.

arnon2020efficient.pdf
2019
Energy window stochastic density functional theory
Chen, M. ; Baer, R. ; Neuhauser, D. ; Rabani, E. Energy window stochastic density functional theory. The Journal of Chemical Physics 2019, 151, 114116. Publisher's Version chen2019energy.pdf
Stochastic Resolution of Identity for Real-Time Second-Order Green’s Function: Ionization Potential and Quasi-Particle Spectrum
Dou, W. ; Takeshita, T. Y. ; Chen, M. ; Baer, R. ; Neuhauser, D. ; Rabani, E. Stochastic Resolution of Identity for Real-Time Second-Order Green’s Function: Ionization Potential and Quasi-Particle Spectrum. Journal of Chemical Theory and Computation 2019. Publisher's VersionAbstract

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.

dou2019stochastic.pdf

Pages