Showing 1–50 of 180 results for author: Reichman, D R

Coupling of Electronic Transitions to Ferroelectric Order in a 2D Semiconductor

Authors: Chun-Ying Huang, Daniel G. Chica, Zhi-Hao Cui, Taketo Handa, Morgan Thinel, Nicholas Olsen, Yufeng Liu, Michael E. Ziebel, Guiying He, Yinming Shao, Connor A. Occhialini, Jonathan Pelliciari, Dmitri N. Basov, Matthew Sfeir, Abhay Pasupathy, Valentina Bisogni, David R. Reichman, Xavier Roy, Xiaoyang Zhu

Abstract: A ferroelectric material often exhibits a soft transvers optical (TO) phonon mode which governs it phase transition. Charge coupling to this ferroelectric soft mode may further mediate emergent physical properties, including superconductivity and defect tolerance. However, direct experimental evidence for such coupling is scarce. Here we show that a photo-launched coherent phonon couples strongly… ▽ More A ferroelectric material often exhibits a soft transvers optical (TO) phonon mode which governs it phase transition. Charge coupling to this ferroelectric soft mode may further mediate emergent physical properties, including superconductivity and defect tolerance. However, direct experimental evidence for such coupling is scarce. Here we show that a photo-launched coherent phonon couples strongly to electronic transitions across the bandgap in the van der Waals (vdW) two-dimensional (2D) ferroelectric semiconductor NbOI2. Using terahertz time-domain spectroscopy and first-principles calculations, we identify this mode as the TO phonon responsible for ferroelectric order. This exclusive coupling occurs only with above-gap electronic transition and is absent in the valence band as revealed by resonant inelastic X-ray scattering. Our findings suggest a new role of the soft TO phonon mode in electronic and optical properties of ferroelectric semiconductors. △ Less

Submitted 11 October, 2024; originally announced October 2024.

Comments: 18 pages, 5 figures, 13 pages SI

First-Principles Phase-Field Modeling

Authors: Jaehyeok Jin, David R. Reichman

Abstract: Phase-field methods offer a versatile computational framework for simulating large-scale microstructure evolution. However, the applicability and predictability of phase-field models are inherently limited by their ad hoc nature, and there is currently no bottom-up theory available that enables truly first-principles predictive modeling of large-scale non-equilibrium processes. Here, we present a… ▽ More Phase-field methods offer a versatile computational framework for simulating large-scale microstructure evolution. However, the applicability and predictability of phase-field models are inherently limited by their ad hoc nature, and there is currently no bottom-up theory available that enables truly first-principles predictive modeling of large-scale non-equilibrium processes. Here, we present a bottom-up framework that provides a route to the construction of mesoscopic phase-field models entirely based on atomistic information. By introducing a molecular coarse-grained system as an intermediate step, we demonstrate the approach on the example of ice nucleation dynamics, with a spatiotemporal scale-up of nearly $10^8$ times compared to the microscopic model. Our framework offers a unique approach for incorporating atomistic details into mesoscopic models, systematically bridging the gap between microscopic particle-based simulations and field-theoretic models. △ Less

Submitted 10 October, 2024; originally announced October 2024.

Comments: 6 pages, 4 figures, see ancillary file for Supplemental Material

Beyond CCSD(T) accuracy at lower scaling with auxiliary field quantum Monte Carlo

Authors: Ankit Mahajan, James H. Thorpe, Jo S. Kurian, David R. Reichman, Devin A. Matthews, Sandeep Sharma

Abstract: We introduce a black-box auxiliary field quantum Monte Carlo (AFQMC) approach to perform highly accurate electronic structure calculations using configuration interaction singles and doubles (CISD) trial states. This method consistently provides more accurate energy estimates than coupled cluster singles and doubles with perturbative triples (CCSD(T)), often regarded as the gold standard in quantu… ▽ More We introduce a black-box auxiliary field quantum Monte Carlo (AFQMC) approach to perform highly accurate electronic structure calculations using configuration interaction singles and doubles (CISD) trial states. This method consistently provides more accurate energy estimates than coupled cluster singles and doubles with perturbative triples (CCSD(T)), often regarded as the gold standard in quantum chemistry. This level of precision is achieved at a lower asymptotic computational cost, scaling as $O(N^6)$ compared to the $O(N^7)$ scaling of CCSD(T). We provide numerical evidence supporting these findings through results for challenging main group and transition metal-containing molecules. △ Less

Submitted 3 October, 2024; originally announced October 2024.

Comments: 17 pages, 11 figures

Simulating anharmonic vibrational polaritons beyond the long wavelength approximation

Authors: Dipti Jasrasaria, Arkajit Mandal, David R. Reichman, Timothy C. Berkelbach

Abstract: In this work we investigate anharmonic vibrational polaritons formed due to strong light-matter interactions in an optical cavity between radiation modes and anharmonic vibrations beyond the long-wavelength limit. We introduce a conceptually simple description of light-matter interactions, where spatially localized cavity radiation modes couple to localized vibrations. Within this theoretical fram… ▽ More In this work we investigate anharmonic vibrational polaritons formed due to strong light-matter interactions in an optical cavity between radiation modes and anharmonic vibrations beyond the long-wavelength limit. We introduce a conceptually simple description of light-matter interactions, where spatially localized cavity radiation modes couple to localized vibrations. Within this theoretical framework, we employ self-consistent phonon theory and vibrational dynamical mean-field theory to efficiently simulate momentum-resolved vibrational-polariton spectra, including effects of anharmonicity. Numerical simulations in model systems demonstrate the accuracy and applicability of our approach. △ Less

Submitted 12 September, 2024; originally announced September 2024.

Comments: 10 pages, 6 figures

Correlation Functions From Tensor Network Influence Functionals: The Case of the Spin-Boson Model

Authors: Haimi Nguyen, Nathan Ng, Lachlan P. Lindoy, Gunhee Park, Andrew J. Millis, Garnet Kin-Lic Chan, David R. Reichman

Abstract: We investigate the application of matrix product state (MPS) representations of the influence functionals (IF) for the calculation of real-time equilibrium correlation functions in open quantum systems. Focusing specifically on the unbiased spin-boson model, we explore the use of IF-MPSs for complex time propagation, as well as IF-MPSs for constructing correlation functions in the steady state. We… ▽ More We investigate the application of matrix product state (MPS) representations of the influence functionals (IF) for the calculation of real-time equilibrium correlation functions in open quantum systems. Focusing specifically on the unbiased spin-boson model, we explore the use of IF-MPSs for complex time propagation, as well as IF-MPSs for constructing correlation functions in the steady state. We examine three different IF approaches: one based on the Kadanoff-Baym contour targeting correlation functions at all times, one based on a complex contour targeting the correlation function at a single time, and a steady state formulation which avoids imaginary or complex times, while providing access to correlation functions at all times. We show that within the IF language, the steady state formulation provides a powerful approach to evaluate equilibrium correlation functions. △ Less

Submitted 22 June, 2024; originally announced June 2024.

Benchmarking the Exponential Ansatz for the Holstein model

Authors: Junjie Yang, Zhi-Hao Cui, Ankit Mahajan, Huanchen Zhai, David R. Reichman, Garnet Kin-Lic Chan

Abstract: Polarons are quasiparticles formed as a result of lattice distortions induced by charge carriers. The single-electron Holstein model captures the fundamentals of single polaron physics. We examine the power of the exponential ansatz for the polaron ground-state wavefunction in its coupled cluster, canonical transformation, and (canonically transformed) perturbative variants across the parameter sp… ▽ More Polarons are quasiparticles formed as a result of lattice distortions induced by charge carriers. The single-electron Holstein model captures the fundamentals of single polaron physics. We examine the power of the exponential ansatz for the polaron ground-state wavefunction in its coupled cluster, canonical transformation, and (canonically transformed) perturbative variants across the parameter space of the Holstein model. Our benchmark serves to guide future developments of polaron wavefunctions beyond the single-electron Holstein model. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 10 pages, 6 figures

Structure and dynamics of electron-phonon coupled systems using neural quantum states

Authors: Ankit Mahajan, Paul J. Robinson, Joonho Lee, David R. Reichman

Abstract: In this work, we use neural quantum states (NQS) to describe the high-dimensional wave functions of electron-phonon coupled systems. We demonstrate that NQS can accurately and systematically learn the underlying physics of such problems through a variational Monte Carlo optimization of the energy with minimal incorporation of physical information even in highly challenging cases. We assess the abi… ▽ More In this work, we use neural quantum states (NQS) to describe the high-dimensional wave functions of electron-phonon coupled systems. We demonstrate that NQS can accurately and systematically learn the underlying physics of such problems through a variational Monte Carlo optimization of the energy with minimal incorporation of physical information even in highly challenging cases. We assess the ability of our approach across various lattice model examples featuring different types of couplings. The flexibility of our NQS formulation is demonstrated via application to ab initio models parametrized by density functional perturbation theory consisting of electron or hole bands coupled linearly to dispersive phonons. We compute accurate real-frequency spectral properties of electron-phonon systems via a novel formalism based on NQS. Our work establishes a general framework for exploring diverse ground state and dynamical phenomena arising in electron-phonon systems, including the non-perturbative interplay of correlated electronic and electron-phonon effects in systems ranging from simple lattice models to realistic models of materials parametrized by ab initio calculations. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 17 pages, 13 figures

Ab Initio Polaron Wave Functions

Authors: Paul J. Robinson, Joonho Lee, Ankit Mahajan, David R. Reichman

Abstract: In this work we demonstrate that accurate ground state wave functions may be constructed for polarons in a fully ab initio setting across the wide range of couplings associated with both the large and small polaron limits. We present a single general unitary transformation approach which encompasses an ab initio version of the Lee-Low-Pines theory at weak coupling and the coherent state Landau-Pek… ▽ More In this work we demonstrate that accurate ground state wave functions may be constructed for polarons in a fully ab initio setting across the wide range of couplings associated with both the large and small polaron limits. We present a single general unitary transformation approach which encompasses an ab initio version of the Lee-Low-Pines theory at weak coupling and the coherent state Landau-Pekar framework at strong coupling while interpolating between these limits in general cases. We show that perturbation theory around these limits may be performed in a facile manner to assess the accuracy of the approach, as well as provide an independent route to the ab initio properties of polarons. We test these ideas on the case of LiF, where the electron-polaron is expected to be large and relatively weakly coupled, while the hole-polaron is expected to be a strongly coupled small polaron. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 12 pages, 3 figures

Scalable Ab Initio Electronic Structure Methods with Near Chemical Accuracy for Main Group Chemistry

Authors: Yujing Wei, Sibali Debnath, John L. Weber, Ankit Mahajan, David R. Reichman, Richard A. Friesner

Abstract: This study evaluates the precision of widely recognized quantum chemical methodologies, CCSD(T), DLPNO-CCSD(T) and localized ph-AFQMC, for determining the thermochemistry of main group elements. DLPNO-CCSD(T) and localized ph-AFQMC, which offer greater scalability compared to canonical CCSD(T), have emerged over the last decade as pivotal in producing precise benchmark chemical data. Our investiga… ▽ More This study evaluates the precision of widely recognized quantum chemical methodologies, CCSD(T), DLPNO-CCSD(T) and localized ph-AFQMC, for determining the thermochemistry of main group elements. DLPNO-CCSD(T) and localized ph-AFQMC, which offer greater scalability compared to canonical CCSD(T), have emerged over the last decade as pivotal in producing precise benchmark chemical data. Our investigation includes closed-shell, neutral molecules, focusing on their heat of formation and atomization energy sourced from four specific small molecule datasets. Firstly, we selected molecules from the G2 and G3 datasets, noted for their reliable experimental heat of formation data. Additionally, we incorporate molecules from the W4-11 and W4-17 sets, which provide high-level theoretical reference values for atomization energy at 0 K. Our findings reveal that both DLPNO-CCSD(T) and ph-AFQMC methods are capable of achieving a root-mean-square deviation (RMSD) of less than 1 kcal/mol across the combined dataset, aligning with the threshold for chemical accuracy. Moreover, we make efforts to confine the maximum deviations within 2 kcal/mol, a degree of precision that significantly broadens the applicability of these methods in fields such as biology and materials science. △ Less

Submitted 29 July, 2024; v1 submitted 2 May, 2024; originally announced May 2024.

Journal ref: J. Phys. Chem. A 2024, 128, 5796-5807

Investigating the Collective Nature of Cavity Modified Chemical Kinetics under Vibrational Strong Coupling

Authors: Lachlan P. Lindoy, Arkajit Mandal, David R. Reichman

Abstract: In this paper we develop quantum dynamical methods capable of treating the dynamics of chemically reacting systems in an optical cavity in the vibrationally strong-coupling (VSC) limit at finite temperatures and in the presence of a dissipative solvent in both the few and many molecule limits. In the context of two simple models we demonstrate how reactivity in the {\em collective} VSC regime does… ▽ More In this paper we develop quantum dynamical methods capable of treating the dynamics of chemically reacting systems in an optical cavity in the vibrationally strong-coupling (VSC) limit at finite temperatures and in the presence of a dissipative solvent in both the few and many molecule limits. In the context of two simple models we demonstrate how reactivity in the {\em collective} VSC regime does not exhibit altered rate behavior in equilibrium, but may exhibit resonant cavity modification of reactivity when the system is explicitly out of equilibrium. Our results suggest experimental protocols that may be used to modify reactivity in the collective regime and point to features not included in the models studied which demand further scrutiny. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Coarse-Graining in Space versus Time

Authors: Jaehyeok Jin, David R. Reichman, Jeppe C. Dyre, Ulf R. Pedersen

Abstract: Understanding the structure and dynamics of liquids is pivotal for the study of larger spatiotemporal processes, particularly for glass-forming materials at low temperatures. The so-called thermodynamic scaling relation, validated for many molecular systems through experiments, offers an efficient means to explore a vast range of time scales along a one-dimensional phase diagram. Isomorph theory p… ▽ More Understanding the structure and dynamics of liquids is pivotal for the study of larger spatiotemporal processes, particularly for glass-forming materials at low temperatures. The so-called thermodynamic scaling relation, validated for many molecular systems through experiments, offers an efficient means to explore a vast range of time scales along a one-dimensional phase diagram. Isomorph theory provides a theoretical framework for thermodynamic scaling based on strong virial-potential energy correlations, but this approach is most successful for simple point particles. In particular, isomorph theory has resisted extension to complex molecular liquids due to the existence of high-frequency intramolecular interactions. To elucidate the microscopic origin of density scaling for molecular systems, we employ two distinct approaches for coarse-graining in space or in time. The former eliminates fast degrees of freedom by reducing a molecule to a center-of-mass-level description, while the latter involves temporally averaged fluctuations or correlation functions over the characteristic time scale. We show that both approaches yield a consistent density scaling coefficient for ortho-terphenyl, which is moreover in agreement with the experimental value. Building upon these findings, we derive the density scaling relationship exhibiting a single-parameter phase diagram from fully atomistic simulations. Our results unravel the microscopic nature underlying thermodynamic scaling and shed light on the role of coarse-graining for assessing the slow fluctuations in molecular systems, ultimately enabling the extension of systematic bottom-up approaches to larger and more complex molecular liquids that are experimentally challenging to probe. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 24 pages, 18 figures

Tensor network influence functionals in the continuous-time limit: connections to quantum embedding, bath discretization, and higher-order time propagation

Authors: Gunhee Park, Nathan Ng, David R. Reichman, Garnet Kin-Lic Chan

Abstract: We describe two developments of tensor network influence functionals (in particular, influence functional matrix product states (IF-MPS)) for quantum impurity dynamics within the fermionic setting of the Anderson impurity model. The first provides the correct extension of the IF-MPS to continuous time by introducing a related mathematical object, the boundary influence functional MPS. The second c… ▽ More We describe two developments of tensor network influence functionals (in particular, influence functional matrix product states (IF-MPS)) for quantum impurity dynamics within the fermionic setting of the Anderson impurity model. The first provides the correct extension of the IF-MPS to continuous time by introducing a related mathematical object, the boundary influence functional MPS. The second connects the dynamics described by a compressed IF-MPS to that of a quantum embedding method with a time-dependent effective bath undergoing nonunitary dynamics. Using these concepts, we implement higher-order time propagators for the quench dynamics of the Anderson impurity model within the boundary IF-MPS formalism. The calculations illustrate the ability of the current formulation to efficiently remove the time step error in standard discrete-time IF-MPS implementations as well as to interface with state vector propagation techniques. They also show the advantages of IF-MPS dynamics, with its associated highly compact effective bath dynamics, over state vector propagation with a static bath discretization. △ Less

Submitted 2 July, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Journal ref: Physical Review B 110, 045104 (2024)

Dynamical Facilitation Governs the Equilibration Dynamics of Glasses

Authors: Rahul N. Chacko, François P. Landes, Giulio Biroli, Olivier Dauchot, Andrea J. Liu, David R. Reichman

Abstract: Convincing evidence of domain growth in the heating of ultrastable glasses suggests that the equilibration dynamics of super-cooled liquids could be driven by a nucleation and growth mechanism. We investigate this possibility by simulating the equilibration dynamics of a model glass during both heating and cooling between poorly and well-annealed states. Though we do observe the growth of domains… ▽ More Convincing evidence of domain growth in the heating of ultrastable glasses suggests that the equilibration dynamics of super-cooled liquids could be driven by a nucleation and growth mechanism. We investigate this possibility by simulating the equilibration dynamics of a model glass during both heating and cooling between poorly and well-annealed states. Though we do observe the growth of domains during heating, we find that domains are absent during cooling. This absence is inconsistent with classical nucleation theory. By comparing the equilibration dynamics of our glass with that of two models with kinetic constraints, we demonstrate that dynamical facilitation generically leads to heating driven by domain growth and cooling without domains. Our results provide strong evidence that dynamical facilitation, not nucleation and interfacial-tension-driven domain growth, is the driving mechanism for the equilibration dynamics of glass-formers. △ Less

Submitted 10 May, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

Comments: 21 pages, 16 figures. Supplemental Movies too large for arXiv. Will provide download link once available. Available on request until then

Journal ref: Physical Review X 14.3 (2024): 031012

Hierarchical Framework for Predicting Entropies in Bottom-Up Coarse-Grained Models

Authors: Jaehyeok Jin, David R. Reichman

Abstract: The thermodynamic entropy of coarse-grained (CG) models stands as one of the most important properties for quantifying the missing information during the CG process and for establishing transferable (or extendible) CG interactions. However, performing additional CG simulations on top of model construction often leads to significant additional computational overhead. In this work, we propose a simp… ▽ More The thermodynamic entropy of coarse-grained (CG) models stands as one of the most important properties for quantifying the missing information during the CG process and for establishing transferable (or extendible) CG interactions. However, performing additional CG simulations on top of model construction often leads to significant additional computational overhead. In this work, we propose a simple hierarchical framework for predicting the thermodynamic entropies of various molecular CG systems. Our approach employs a decomposition of the CG interactions, enabling the estimation of the CG partition function and thermodynamic properties a priori. Starting from the ideal gas description, we leverage classical perturbation theory to systematically incorporate simple yet essential interactions, ranging from the hard sphere model to the generalized van der Waals model. Additionally, we propose an alternative approach based on multiparticle correlation functions, allowing for systematic improvements through higher-order correlations. Numerical applications to molecular liquids validate the high fidelity of our approach, and our computational protocols demonstrate that a reduced model with simple energetics can reasonably estimate the thermodynamic entropy of CG models without performing any CG simulations. Overall, our findings present a systematic framework for estimating not only the entropy but also other thermodynamic properties of CG models, relying solely on information from the reference system. △ Less

Submitted 13 February, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: 36 pages, 6 figures

Journal ref: J. Phys. Chem. B 128 (13), 3182-3199 (2024)

Variational Lang-Firsov approach plus Møller-Plesset perturbation theory with applications to ab initio polariton chemistry

Authors: Zhi-Hao Cui, Arkajit Mandal, David R. Reichman

Abstract: We apply the Lang-Firsov (LF) transformation to electron-boson coupled Hamiltonians and variationally optimize the transformation parameters and molecular orbital coefficients to determine the ground state. Møller-Plesset (MP-$n$, with $n = 2$ and $4$) perturbation theory is then performed on top of the optimized LF mean-field state to improve the description of electron-electron and electron-boso… ▽ More We apply the Lang-Firsov (LF) transformation to electron-boson coupled Hamiltonians and variationally optimize the transformation parameters and molecular orbital coefficients to determine the ground state. Møller-Plesset (MP-$n$, with $n = 2$ and $4$) perturbation theory is then performed on top of the optimized LF mean-field state to improve the description of electron-electron and electron-boson correlations. The method (LF-MP) is applied to several electron-boson coupled systems, including the Hubbard-Holstein model, diatomic molecule dissociation ($\text{H}_2$, $\text{HF}$), and the modification of proton transfer reactions (malonaldehyde and aminopropenal) via the formation of polaritons in an optical cavity. We show that with a correction for the electron-electron correlation, the method gives quantitatively accurate energies comparable to exact diagonalization or coupled-cluster theory. The effect of multiple photon modes, spin polarization, and the comparison to the coherent state MP theory are also discussed. △ Less

Submitted 2 February, 2024; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 16 pages, 8 figures

Journal ref: J. Chem. Theory Comput. (2024)

Perturbative Expansion in Reciprocal Space: Bridging Microscopic and Mesoscopic Descriptions of Molecular Interactions

Authors: Jaehyeok Jin, David R. Reichman

Abstract: Determining the Fourier representation of various molecular interactions is important for constructing density-based field theories from a microscopic point of view, enabling a multiscale bridge between microscopic and mesoscopic descriptions. However, due to the strongly repulsive nature of short-ranged interactions, interparticle interactions cannot be formally defined in Fourier space, which re… ▽ More Determining the Fourier representation of various molecular interactions is important for constructing density-based field theories from a microscopic point of view, enabling a multiscale bridge between microscopic and mesoscopic descriptions. However, due to the strongly repulsive nature of short-ranged interactions, interparticle interactions cannot be formally defined in Fourier space, which renders coarse-grained approaches in $\textit{k}$-space somewhat ambiguous. In this paper, we address this issue by designing a perturbative expansion of pair interactions in reciprocal space. Our perturbation theory, starting from reciprocal space, elucidates the microscopic origins underlying zeroth-order (long-range attractions) and divergent repulsive interactions from higher-order contributions. We propose a systematic framework for constructing a faithful Fourier space representation of molecular interactions, capturing key structural correlations in various systems, including simple model systems and molecular coarse-grained models of liquids. Building upon the Ornstein-Zernike equation, our approach can be combined with appropriate closure relations, and to further improve the closure approximations, we develop a bottom-up parametrization strategy for inferring the bridge function from microscopic statistics. By incorporating the bridge function into the Fourier representation, our findings suggest a systematic, bottom-up approach to performing coarse-graining in reciprocal space, leading to the systematic construction of a bottom-up classical field theory of complex aqueous systems. △ Less

Submitted 21 November, 2023; v1 submitted 12 September, 2023; originally announced September 2023.

Comments: 54 pages, 11 figures

Journal ref: J. Phys. Chem. B 128 (4), 1061-1078 (2024)

Response properties in phaseless auxiliary field quantum Monte Carlo

Authors: Ankit Mahajan, Jo S. Kurian, Joonho Lee, David R. Reichman, Sandeep Sharma

Abstract: We present a method for calculating first-order response properties in phaseless auxiliary field quantum Monte Carlo (AFQMC) through the application of automatic differentiation (AD). Biases and statistical efficiency of the resulting estimators are discussed. Our approach demonstrates that AD enables the calculation of reduced density matrices (RDMs) with the same computational cost scaling as en… ▽ More We present a method for calculating first-order response properties in phaseless auxiliary field quantum Monte Carlo (AFQMC) through the application of automatic differentiation (AD). Biases and statistical efficiency of the resulting estimators are discussed. Our approach demonstrates that AD enables the calculation of reduced density matrices (RDMs) with the same computational cost scaling as energy calculations, accompanied by a cost prefactor of less than four in our numerical calculations. We investigate the role of self-consistency and trial orbital choice in property calculations. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: 10 pages, 6 figures

Competing Generalized Wigner Crystal States in Moiré Heterostructures

Authors: Shu Fay Ung, Joonho Lee, David R. Reichman

Abstract: We present a comprehensive study of generalized Wigner crystals across various filling factors for system sizes up to 162 holes employing Hartree-Fock theory and explicitly correlated wave function approaches. While we find broad agreement with the behavior observed in experiments and classical Monte Carlo simulations, we highlight the fact that the Hartree-Fock energy landscape appears to be rema… ▽ More We present a comprehensive study of generalized Wigner crystals across various filling factors for system sizes up to 162 holes employing Hartree-Fock theory and explicitly correlated wave function approaches. While we find broad agreement with the behavior observed in experiments and classical Monte Carlo simulations, we highlight the fact that the Hartree-Fock energy landscape appears to be remarkably complex, exhibiting many competing states, both ordered and disordered, separated by energies of a fraction of $\sim$ 1 meV/hole. We demonstrate which of the located states are metastable by performing a stability analysis at the Hartree-Fock level. Correlated wave function methods furthermore reveal small correlation energies that are nevertheless large enough to tip the balance of state ordering found within Hartree-Fock theory. △ Less

Submitted 5 December, 2023; v1 submitted 6 August, 2023; originally announced August 2023.

Comments: 17 pages, 12 figures

Journal ref: Phys. Rev. B 108, 245113 (2023)

The Design of New Practical Constraints in Auxiliary-Field Quantum Monte Carlo

Authors: John L. Weber, Hung Vuong, Richard A. Friesner, David R. Reichman

Abstract: We formulate and characterize a new constraint for Auxiliary Field Quantum Monte Carlo (AFQMC) applicable for general fermionic systems, which allows for the accumulation of phase in the random walk but disallows walkers with a magnitude of phase greater than $π$ with respect to the trial wave function. For short imaginary times, before walkers accumulate sizable phase values, this approach is equ… ▽ More We formulate and characterize a new constraint for Auxiliary Field Quantum Monte Carlo (AFQMC) applicable for general fermionic systems, which allows for the accumulation of phase in the random walk but disallows walkers with a magnitude of phase greater than $π$ with respect to the trial wave function. For short imaginary times, before walkers accumulate sizable phase values, this approach is equivalent to exact free projection, allowing one to observe the accumulation of bias associated with the constraint and thus estimate its magnitude \textit{a priori}. We demonstrate the stability of this constraint over arbitrary imaginary times and system sizes, highlighting the removal of noise due to the fermionic sign problem. Benchmark total energies for a variety of weakly and strongly correlated molecular systems reveal a distinct bias with respect to standard phaseless AFQMC, with a comparative increase in accuracy given a sufficient quality of trial wave function for the set of studied cases. We then take this constraint, referred to as linecut (lc-) AFQMC, and systematically release it (lcR-AFQMC), providing a route to obtain a smooth bridge between constrained AFQMC and the exact free projection results. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: 26 pages, 11 figures, 3 tables

Correlated Anharmonicity and Dynamic Disorder Control Carrier Transport in Halide Perovskites

Authors: Maximilian J. Schilcher, David J. Abramovitch, Matthew Z. Mayers, Liang Z. Tan, David R. Reichman, David A. Egger

Abstract: Halide pervoskites are an important class of semiconducting materials which hold great promise for optoelectronic applications. In this work we investigate the relationship between vibrational anharmonicity and dynamic disorder in this class of solids. Via a multi-scale model parameterized from first-principles calculations, we demonstrate that the non-Gaussian lattice motion in halide perovskites… ▽ More Halide pervoskites are an important class of semiconducting materials which hold great promise for optoelectronic applications. In this work we investigate the relationship between vibrational anharmonicity and dynamic disorder in this class of solids. Via a multi-scale model parameterized from first-principles calculations, we demonstrate that the non-Gaussian lattice motion in halide perovskites is microscopically connected to the dynamic disorder of overlap fluctuations among electronic states. This connection allows us to rationalize the emergent differences in temperature-dependent mobilities of prototypical MAPbI$_3$ and MAPbBr$_3$ compounds across structural phase-transitions, in agreement with experimental findings. Our analysis suggests that the details of vibrational anharmonicity and dynamic disorder can complement known predictors of electronic conductivity and can provide structure-property guidelines for the tuning of carrier transport characteristics in anharmonic semiconductors. △ Less

Submitted 2 August, 2023; v1 submitted 23 May, 2023; originally announced May 2023.

Microscopic Theory of Multimode Polariton Dispersion in Multilayered Materials

Authors: Arkajit Mandal, Ding Xu, Ankit Mahajan, Joonho Lee, Milan E. Delor, David R. Reichman

Abstract: We develop a microscopic theory for the multimode polariton dispersion in materials coupled to cavity radiation modes. Starting from a microscopic light-matter Hamiltonian, we devise a general strategy for obtaining simple matrix models of polariton dispersion curves based on the structure and spatial location of multi-layered 2D materials inside the optical cavity. Our theory exposes the connecti… ▽ More We develop a microscopic theory for the multimode polariton dispersion in materials coupled to cavity radiation modes. Starting from a microscopic light-matter Hamiltonian, we devise a general strategy for obtaining simple matrix models of polariton dispersion curves based on the structure and spatial location of multi-layered 2D materials inside the optical cavity. Our theory exposes the connections between seemingly distinct models that have been employed in the literature and resolves an ambiguity that has arisen concerning the experimental description of the polaritonic band structure. We demonstrate the applicability of our theoretical formalism by fabricating various geometries of multi-layered perovskite materials coupled to cavities and demonstrating that our theoretical predictions agree with the experimental results presented here. △ Less

Submitted 19 March, 2023; originally announced March 2023.

Stochastic Real-Time Second-Order Green's Function Theory for Neutral Excitations in Molecules and Nanostructures

Authors: Leopoldo Mejía, Jia Yin, David R. Reichman, Roi Baer, Chao Yang, Eran Rabani

Abstract: We present a real-time second-order Green's function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of $O(N_{\rm e}^3$), where $N_{\rm e}$ is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals (ERI). To improve the time-propagation a… ▽ More We present a real-time second-order Green's function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of $O(N_{\rm e}^3$), where $N_{\rm e}$ is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals (ERI). To improve the time-propagation and the spectral resolution, we adopt the dynamic mode decomposition (DMD) technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Journal ref: J. Chem. Theory Comput. 2023, 19, 16, 5563-5571

Remembering the work of Phillip L. Geissler: A coda to his scientific trajectory

Authors: Gregory R. Bowman, Stephen J. Cox, Christoph Dellago, Kateri H. DuBay, Joel D. Eaves, Daniel A. Fletcher, Layne B. Frechette, Michael Grünwald, Katherine Klymko, JiYeon Ku, Ahmad K. Omar, Eran Rabani, David R. Reichman, Julia R. Rogers, Andreana M. Rosnik, Grant M. Rotskoff, Anna R. Schneider, Nadine Schwierz, David A. Sivak, Suriyanarayanan Vaikuntanathan, Stephen Whitelam, Asaph Widmer-Cooper

Abstract: Phillip L. Geissler made important contributions to the statistical mechanics of biological polymers, heterogeneous materials, and chemical dynamics in aqueous environments. He devised analytical and computational methods that revealed the underlying organization of complex systems at the frontiers of biology, chemistry, and materials science. In this retrospective, we celebrate his work at these… ▽ More Phillip L. Geissler made important contributions to the statistical mechanics of biological polymers, heterogeneous materials, and chemical dynamics in aqueous environments. He devised analytical and computational methods that revealed the underlying organization of complex systems at the frontiers of biology, chemistry, and materials science. In this retrospective, we celebrate his work at these frontiers. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Journal ref: Ann. Rev. Phys. Chem. 74, 11.1-11.27 (2023)

Quantum simulation of exact electron dynamics can be more efficient than classical mean-field methods

Authors: Ryan Babbush, William J. Huggins, Dominic W. Berry, Shu Fay Ung, Andrew Zhao, David R. Reichman, Hartmut Neven, Andrew D. Baczewski, Joonho Lee

Abstract: Quantum algorithms for simulating electronic ground states are slower than popular classical mean-field algorithms such as Hartree-Fock and density functional theory, but offer higher accuracy. Accordingly, quantum computers have been predominantly regarded as competitors to only the most accurate and costly classical methods for treating electron correlation. However, here we tighten bounds showi… ▽ More Quantum algorithms for simulating electronic ground states are slower than popular classical mean-field algorithms such as Hartree-Fock and density functional theory, but offer higher accuracy. Accordingly, quantum computers have been predominantly regarded as competitors to only the most accurate and costly classical methods for treating electron correlation. However, here we tighten bounds showing that certain first quantized quantum algorithms enable exact time evolution of electronic systems with exponentially less space and polynomially fewer operations in basis set size than conventional real-time time-dependent Hartree-Fock and density functional theory. Although the need to sample observables in the quantum algorithm reduces the speedup, we show that one can estimate all elements of the $k$-particle reduced density matrix with a number of samples scaling only polylogarithmically in basis set size. We also introduce a more efficient quantum algorithm for first quantized mean-field state preparation that is likely cheaper than the cost of time evolution. We conclude that quantum speedup is most pronounced for finite temperature simulations and suggest several practically important electron dynamics problems with potential quantum advantage. △ Less

Submitted 3 January, 2023; originally announced January 2023.

Comments: 31 pages, 2 tables, 1 figure

Journal ref: Nat. Comms 14: 4058 (2023)

Instantons and the quantum bound to chaos

Authors: Vijay Ganesh Sadhasivam, Lars Meuser, David R. Reichman, Stuart C. Althorpe

Abstract: We investigate why the Lyapunov exponents $λ$ of out-of-time-ordered correlators (OTOCs) satisfy a universal bound $λ< 2 πk_B T/{\hbar}$ by probing imaginary-time path-integral space using ring-polymer molecular dynamics (RPMD) for a barrier-crossing model. We find that the RPMD OTOC satisfies the same bound as the quantum OTOC, which is caused by the stability of quantum thermal fluctuations arou… ▽ More We investigate why the Lyapunov exponents $λ$ of out-of-time-ordered correlators (OTOCs) satisfy a universal bound $λ< 2 πk_B T/{\hbar}$ by probing imaginary-time path-integral space using ring-polymer molecular dynamics (RPMD) for a barrier-crossing model. We find that the RPMD OTOC satisfies the same bound as the quantum OTOC, which is caused by the stability of quantum thermal fluctuations around the instanton on the barrier. We expect that similar instantons (or other delocalised structures) will be found in many other systems with exponentially growing OTOCs. △ Less

Submitted 24 January, 2023; v1 submitted 20 December, 2022; originally announced December 2022.

Finding defects in glasses through machine learning

Authors: Simone Ciarella, Dmytro Khomenko, Ludovic Berthier, Felix C. Mocanu, David R. Reichman, Camille Scalliet, Francesco Zamponi

Abstract: Structural defects control the kinetic, thermodynamic and mechanical properties of glasses. For instance, rare quantum tunneling two-level systems (TLS) govern the physics of glasses at very low temperature. Because of their extremely low density, it is very hard to directly identify them in computer simulations. We introduce a machine learning approach to efficiently explore the potential energy… ▽ More Structural defects control the kinetic, thermodynamic and mechanical properties of glasses. For instance, rare quantum tunneling two-level systems (TLS) govern the physics of glasses at very low temperature. Because of their extremely low density, it is very hard to directly identify them in computer simulations. We introduce a machine learning approach to efficiently explore the potential energy landscape of glass models and identify desired classes of defects. We focus in particular on TLS and we design an algorithm that is able to rapidly predict the quantum splitting between any two amorphous configurations produced by classical simulations. This in turn allows us to shift the computational effort towards the collection and identification of a larger number of TLS, rather than the useless characterization of non-tunneling defects which are much more abundant. Finally, we interpret our machine learning model to understand how TLS are identified and characterized, thus giving direct physical insight into their microscopic nature. △ Less

Submitted 16 May, 2023; v1 submitted 11 December, 2022; originally announced December 2022.

Journal ref: Nature Communications 14, 4229 (2023)

Machine learning potentials from transfer learning of periodic correlated electronic structure methods: Application to liquid water with AFQMC, CCSD, and CCSD(T)

Authors: Michael S. Chen, Joonho Lee, Hong-Zhou Ye, Timothy C. Berkelbach, David R. Reichman, Thomas E. Markland

Abstract: Obtaining the atomistic structure and dynamics of disordered condensed phase systems from first principles remains one of the forefront challenges of chemical theory. Here we exploit recent advances in periodic electronic structure to show that, by leveraging transfer learning starting from lower tier electronic structure methods, one can obtain machine learned potential energy surfaces for liquid… ▽ More Obtaining the atomistic structure and dynamics of disordered condensed phase systems from first principles remains one of the forefront challenges of chemical theory. Here we exploit recent advances in periodic electronic structure to show that, by leveraging transfer learning starting from lower tier electronic structure methods, one can obtain machine learned potential energy surfaces for liquid water from the higher tier AFQMC, CCSD, and CCSD(T) approaches using $\le$200 energies. By performing both classical and path integral molecular dynamics simulations on these machine learned potential energy surfaces we uncover the interplay of dynamical electron correlation and nuclear quantum effects across the entire liquid range of water while providing a general strategy for efficiently utilizing periodic correlated electronic structure methods to explore disordered condensed phase systems. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Real time evolution of Anderson impurity models via tensor network influence functionals

Authors: Nathan Ng, Gunhee Park, Andrew J. Millis, Garnet Kin-Lic Chan, David R. Reichman

Abstract: In this work we present and analyze two tensor network-based influence functional approaches for simulating the real-time dynamics of quantum impurity models such as the Anderson model. Via comparison with recent numerically exact simulations, we show that such methods accurately capture the long-time non-equilibrium quench dynamics. The two parameters that must be controlled in these tensor netwo… ▽ More In this work we present and analyze two tensor network-based influence functional approaches for simulating the real-time dynamics of quantum impurity models such as the Anderson model. Via comparison with recent numerically exact simulations, we show that such methods accurately capture the long-time non-equilibrium quench dynamics. The two parameters that must be controlled in these tensor network influence functional approaches are a time discretization (Trotter) error and a bond dimension (tensor network truncation) error. We show that the actual numerical uncertainties are controlled by an intricate interplay of these two approximations which we demonstrate in different regimes. Our work opens the door to using these tensor network influence functional methods as general impurity solvers. △ Less

Submitted 21 November, 2022; v1 submitted 18 November, 2022; originally announced November 2022.

Bipolaronic superconductivity out of a Coulomb gas

Authors: J. Sous, C. Zhang, M. Berciu, D. R. Reichman, B. V. Svistunov, N. V. Prokof'ev, A. J. Millis

Abstract: Employing unbiased sign-problem-free quantum Monte Carlo, we investigate the effects of long-range Coulomb forces on BEC of bipolarons using a model of bond phonon-modulated electron hopping. In absence of long-range repulsion, this model was recently shown to give rise to small-size, light-mass bipolarons that undergo a superfluid transition at high values of the critical transition temperature… ▽ More Employing unbiased sign-problem-free quantum Monte Carlo, we investigate the effects of long-range Coulomb forces on BEC of bipolarons using a model of bond phonon-modulated electron hopping. In absence of long-range repulsion, this model was recently shown to give rise to small-size, light-mass bipolarons that undergo a superfluid transition at high values of the critical transition temperature $T_\mathrm{c}$. We find that $T_\mathrm{c}$ in our model even with the long-range Coulomb repulsion remains much larger than that of Holstein bipolarons, and can be on the order of or greater than the typical upper bounds on phonon-mediated $T_\mathrm{c}$ based on the Migdal-Eliashberg and McMillan approximations. Our work points to a physically simple mechanism for superconductivity in the low-density regime that may be relevant to current experiments on dilute superconductors. △ Less

Submitted 1 December, 2022; v1 submitted 25 October, 2022; originally announced October 2022.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. B 108, L220502 (2023)

The Generalized Green's function Cluster Expansion: A Python package for simulating polarons

Authors: Matthew R. Carbone, Stepan Fomichev, Andrew J. Millis, Mona Berciu, David R. Reichman, John Sous

Abstract: We present an efficient implementation of the Generalized Green's function Cluster Expansion (GGCE), which is a new method for computing the ground-state properties and dynamics of polarons (single electrons coupled to lattice vibrations) in model electron-phonon systems. The GGCE works at arbitrary temperature and is well suited for a variety of electron-phonon couplings, including, but not limit… ▽ More We present an efficient implementation of the Generalized Green's function Cluster Expansion (GGCE), which is a new method for computing the ground-state properties and dynamics of polarons (single electrons coupled to lattice vibrations) in model electron-phonon systems. The GGCE works at arbitrary temperature and is well suited for a variety of electron-phonon couplings, including, but not limited to, site and bond Holstein and Peierls (Su-Schrieffer-Heeger) couplings, and couplings to multiple phonon modes with different energy scales and coupling strengths. Quick calculations can be performed efficiently on a laptop using solvers from NumPy and SciPy, or in parallel at scale using the PETSc sparse linear solver engine. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 3 pages, software can be found open source under the BSD-3-clause license at github.com/x94carbone/GGCE

Journal ref: J. Open Source Softw. 8, 5115 (2023)

Quantum Dynamics of Vibrational Polariton Chemistry

Authors: Lachlan P. Lindoy, Arkajit Mandal, David R. Reichman

Abstract: We employ an exact quantum mechanical simulation technique to investigate a model of cavity-modified chemical reactions in the condensed phase. The model contains the coupling of the reaction coordinate to a generic solvent, cavity coupling to either the reaction coordinate or a non-reactive mode, and the coupling of the cavity to lossy modes. Thus, many of the most important features needed for r… ▽ More We employ an exact quantum mechanical simulation technique to investigate a model of cavity-modified chemical reactions in the condensed phase. The model contains the coupling of the reaction coordinate to a generic solvent, cavity coupling to either the reaction coordinate or a non-reactive mode, and the coupling of the cavity to lossy modes. Thus, many of the most important features needed for realistic modeling of the cavity modification of chemical reactions are included. We find that when a molecule is coupled to an optical cavity it is essential to treat the problem quantum mechanically in order to obtain a quantitative account of alterations to reactivity. We find sizable and sharp changes in the rate constant that are associated with quantum mechanical state splittings and resonances. The features that emerge from our simulations are closer to those observed in experiments than are previous calculations, even for realistically small values of coupling and cavity loss. This work highlights the importance of a fully quantum treatment of vibrational polariton chemistry. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Microscopic observation of two-level systems in a metallic glass model

Authors: Felix C. Mocanu, Ludovic Berthier, Simone Ciarella, Dmytro Khomenko, David R. Reichman, Camille Scalliet, Francesco Zamponi

Abstract: The low-temperature quasi-universal behavior of amorphous solids has been attributed to the existence of spatially-localized tunneling defects found in the low-energy regions of the potential energy landscape. Computational models of glasses can be studied to elucidate the microscopic nature of these defects. Recent simulation work has demonstrated the means of generating stable glassy configurati… ▽ More The low-temperature quasi-universal behavior of amorphous solids has been attributed to the existence of spatially-localized tunneling defects found in the low-energy regions of the potential energy landscape. Computational models of glasses can be studied to elucidate the microscopic nature of these defects. Recent simulation work has demonstrated the means of generating stable glassy configurations for models that mimic metallic glasses using the swap Monte Carlo algorithm. Building on these studies, we present an extensive exploration of the glassy metabasins of the potential energy landscape of a variant of the most widely used model of metallic glasses. We carefully identify tunneling defects and reveal their depletion with increased glass stability. The density of tunneling defects near the experimental glass transition temperature appears to be in good agreement with experimental measurements. △ Less

Submitted 4 January, 2023; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 17 pages, 12 figures

Journal ref: J. Chem. Phys. 158, 014501 (2023)

Dynamical Heterogeneity in Glass-Forming Liquids

Authors: Giulio Biroli, Kunimasa Miyazaki, David R. Reichman

Abstract: We review the phenomena of dynamical heterogeneity in glass-forming systems and its description within replica and mean-field theories of the glass transition. We review the phenomena of dynamical heterogeneity in glass-forming systems and its description within replica and mean-field theories of the glass transition. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: To appear as a chapter in the book "Spin Glass Theory and Far Beyond - Replica Symmetry Breaking after 40 years"

Modern computational studies of the glass transition

Authors: Ludovic Berthier, David R. Reichman

Abstract: The physics of the glass transition and amorphous materials continues to attract the attention of a wide research community after decades of effort. Supercooled liquids and glasses have been studied numerically since the advent of molecular dynamics and Monte Carlo simulations in the last century. Computer studies have greatly enhanced both experimental discoveries and theoretical developments and… ▽ More The physics of the glass transition and amorphous materials continues to attract the attention of a wide research community after decades of effort. Supercooled liquids and glasses have been studied numerically since the advent of molecular dynamics and Monte Carlo simulations in the last century. Computer studies have greatly enhanced both experimental discoveries and theoretical developments and constitute an active and continually expanding research field. Our goal in this review is to provide a modern perspective on this area. We describe the need to go beyond canonical methods to attack a problem that is notoriously difficult in terms of time scales, length scales, and physical observables. We first summarise recent algorithmic developments to achieve enhanced sampling and faster equilibration using replica exchange methods, cluster and swap Monte Carlo algorithms, and other techniques. We then review some major recent advances afforded by these novel tools regarding the statistical mechanical description of the liquid-to-glass transition as well as the mechanical, vibrational and thermal properties of the glassy solid. We finally describe some important challenges for future research. △ Less

Submitted 4 January, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 17 pages, 5 figures, to be published in Nature Reviews Physics

Journal ref: Nat. Rev. Phys. 5, 102 (2023)

Is there evidence for exponential quantum advantage in quantum chemistry?

Authors: Seunghoon Lee, Joonho Lee, Huanchen Zhai, Yu Tong, Alexander M. Dalzell, Ashutosh Kumar, Phillip Helms, Johnnie Gray, Zhi-Hao Cui, Wenyuan Liu, Michael Kastoryano, Ryan Babbush, John Preskill, David R. Reichman, Earl T. Campbell, Edward F. Valeev, Lin Lin, Garnet Kin-Lic Chan

Abstract: The idea to use quantum mechanical devices to simulate other quantum systems is commonly ascribed to Feynman. Since the original suggestion, concrete proposals have appeared for simulating molecular and materials chemistry through quantum computation, as a potential ``killer application''. Indications of potential exponential quantum advantage in artificial tasks have increased interest in this ap… ▽ More The idea to use quantum mechanical devices to simulate other quantum systems is commonly ascribed to Feynman. Since the original suggestion, concrete proposals have appeared for simulating molecular and materials chemistry through quantum computation, as a potential ``killer application''. Indications of potential exponential quantum advantage in artificial tasks have increased interest in this application, thus, it is critical to understand the basis for potential exponential quantum advantage in quantum chemistry. Here we gather the evidence for this case in the most common task in quantum chemistry, namely, ground-state energy estimation. We conclude that evidence for such an exponential advantage across chemical space has yet to be found. While quantum computers may still prove useful for quantum chemistry, it may be prudent to assume exponential speedups are not generically available for this problem. △ Less

Submitted 14 November, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Journal ref: Nat Commun 14, 1952 (2023)

Twenty Years of Auxiliary-Field Quantum Monte Carlo in Quantum Chemistry: An Overview and Assessment on Main Group Chemistry and Bond-Breaking

Authors: Joonho Lee, Hung Q. Pham, David R. Reichman

Abstract: In this work, we present an overview of the phaseless auxiliary-field quantum Monte Carlo (ph- AFQMC) approach from a computational quantum chemistry perspective, and present a numerical assessment of its performance on main group chemistry and bond-breaking problems with a total of 1004 relative energies. While our benchmark study is somewhat limited, we make recommendations for the use of ph-AFQ… ▽ More In this work, we present an overview of the phaseless auxiliary-field quantum Monte Carlo (ph- AFQMC) approach from a computational quantum chemistry perspective, and present a numerical assessment of its performance on main group chemistry and bond-breaking problems with a total of 1004 relative energies. While our benchmark study is somewhat limited, we make recommendations for the use of ph-AFQMC for general main-group chemistry applications. For systems where single determinant wave functions are qualitatively accurate, we expect the accuracy of ph-AFQMC in conjunction with a single determinant trial wave function to be between that of coupled-cluster with singles and doubles (CCSD) and CCSD with perturbative triples (CCSD(T)). For these applications, ph-AFQMC should be a method of choice when canonical CCSD(T) is too expensive to run. For systems where multi-reference (MR) wave functions are needed for qualitative accuracy, ph-AFQMC is far more accurate than MR perturbation theory methods and competitive with MR configuration interaction (MRCI) methods. Due to the computational efficiency of ph-AFQMC compared to MRCI, we recommended ph-AFQMC as a method of choice for handling dynamic correlation in MR problems. We conclude with a discussion of important directions for future development of the ph-AFQMC approach. △ Less

Submitted 27 September, 2022; v1 submitted 2 August, 2022; originally announced August 2022.

Response to "Exponential challenges in unbiasing quantum Monte Carlo algorithms with quantum computers"

Authors: Joonho Lee, David R. Reichman, Ryan Babbush, Nicholas C. Rubin, Fionn D. Malone, Bryan O'Gorman, William J. Huggins

Abstract: A recent preprint by Mazzola and Carleo numerically investigates exponential challenges that can arise for the QC-QMC algorithm introduced in our work, "Unbiasing fermionic quantum Monte Carlo with a quantum computer." As discussed in our original paper, we agree with this general concern. However, here we provide further details and numerics to emphasize that the prospects for practical quantum a… ▽ More A recent preprint by Mazzola and Carleo numerically investigates exponential challenges that can arise for the QC-QMC algorithm introduced in our work, "Unbiasing fermionic quantum Monte Carlo with a quantum computer." As discussed in our original paper, we agree with this general concern. However, here we provide further details and numerics to emphasize that the prospects for practical quantum advantage in QC-QMC remain open. The exponential challenges in QC-QMC are dependent on (1) the choice of QMC methods, (2) the underlying system, and (3) the form of trial and walker wavefunctions. While one can find difficult examples with a specific method, a specific system, and a specific walker/trial form, for some combinations of these choices, the approach is potentially more scalable than other near-term quantum algorithms. Future research should aim to identify examples for which QC-QMC enables practical quantum advantage. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Two-Dimensional Spectroscopy of Two-Dimensional Materials

Authors: Lachlan P Lindoy, Yao-Wen Chang, David R Reichman

Abstract: In this work we provide an exact and efficient numerical approach to simulate multi-time correlation functions in the Mahan-Nozières-De Dominicis model, which crudely mimics the spectral properties of doped two-dimensional semiconductors such as monolayer transition metal dichalcogenides. We apply this approach to study the coherent two-dimensional electronic spectra of the model. We show that sev… ▽ More In this work we provide an exact and efficient numerical approach to simulate multi-time correlation functions in the Mahan-Nozières-De Dominicis model, which crudely mimics the spectral properties of doped two-dimensional semiconductors such as monolayer transition metal dichalcogenides. We apply this approach to study the coherent two-dimensional electronic spectra of the model. We show that several experimentally observed phenomena, such as peak asymmetry and coherent oscillations in the waiting-time dependence of the trion-exciton cross peaks of the two-dimensional rephasing spectrum, emerge naturally in our approach. Additional features are also present which find no correspondence with experimentally expected behavior. We trace these features to the infinite hole mass property of the model. We use this understanding to construct an efficient approach which filters out configurations associated with the lack of exciton recoil, enabling the connection to previous work and providing a route to the construction of realistic two-dimensional spectra over a broad doping range in two-dimensional semiconductors. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: 8 Pages, 6 Figures

Resonant Cavity Modification of Ground State Chemical Kinetics

Authors: Lachlan P. Lindoy, Arkajit Mandal, David R. Reichman

Abstract: Recent experiments have suggested that ground state chemical kinetics can be suppressed or enhanced by coupling the vibrational degrees of freedom of a molecular system with a radiation mode inside an optical cavity. Experiments show that the chemical rate is strongly modified when the photon frequency is close to characteristic vibrational frequencies. The origin of this remarkable effect remains… ▽ More Recent experiments have suggested that ground state chemical kinetics can be suppressed or enhanced by coupling the vibrational degrees of freedom of a molecular system with a radiation mode inside an optical cavity. Experiments show that the chemical rate is strongly modified when the photon frequency is close to characteristic vibrational frequencies. The origin of this remarkable effect remains unknown. In this work, we develop an analytical rate theory for cavity-modified ground state chemical kinetics based on the Pollak-Grabert-Hänggi rate theory. Unlike previous work, our theory covers the complete range of solvent friction values, from the energy-diffusion limited to the spatial-diffusion limited regimes. We show that the chemical reaction rate can either be enhanced or suppressed depending on the bath friction; when bath friction is weak chemical kinetics is enhanced as opposed to the case of strong bath friction, where chemical kinetics is suppressed. Further, we show that the photon frequency at which maximum modification of chemical rate is achieved is close to the reactant well, and hence resonant rate modification occurs. In the strong friction limit the {\it resonant} photon frequency is instead close to the barrier frequency, as obtained using the Grote-Hynes rate theory. Finally, we observe that the rate changes (as a function of photon frequency) are much sharper and more sizable in the weak friction limit than in the strong friction limit, and become increasingly sharp with decreasing well frequency. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Ultrafast imaging of polariton propagation and interactions

Authors: Ding Xu, Arkajit Mandal, James M. Baxter, Shan-Wen Cheng, Inki Lee, Haowen Su, Song Liu, David R. Reichman, Milan Delor

Abstract: Semiconductor excitations can hybridize with cavity photons to form exciton-polaritons (EPs) with remarkable properties, including light-like energy flow combined with matter-like interactions. To fully harness these properties, EPs must retain ballistic, coherent transport despite matter-mediated interactions with lattice phonons. Here we develop a nonlinear momentum-resolved optical approach tha… ▽ More Semiconductor excitations can hybridize with cavity photons to form exciton-polaritons (EPs) with remarkable properties, including light-like energy flow combined with matter-like interactions. To fully harness these properties, EPs must retain ballistic, coherent transport despite matter-mediated interactions with lattice phonons. Here we develop a nonlinear momentum-resolved optical approach that directly images EPs in real space on femtosecond scales in a range of polaritonic architectures. We focus our analysis on EP propagation in layered halide perovskite microcavities. We reveal that EP-phonon interactions lead to a large renormalization of EP velocities at high excitonic fractions at room temperature. Despite these strong EP-phonon interactions, ballistic transport is maintained for up to half-exciton EPs, in agreement with quantum simulations of dynamic disorder shielding through light-matter hybridization. Above 50% excitonic character, rapid decoherence leads to diffusive transport. Our work provides a general framework to precisely balance EP coherence, velocity, and nonlinear interactions. △ Less

Submitted 15 June, 2023; v1 submitted 2 May, 2022; originally announced May 2022.

Bipolaronic high-temperature superconductivity

Authors: C. Zhang, J. Sous, D. R. Reichman, M. Berciu, A. J. Millis, N. V. Prokof'ev, B. V. Svistunov

Abstract: Electron-lattice interactions play a prominent role in quantum materials, making a deeper understanding of direct routes to phonon-mediated high-transition-temperature ($T_{\mathrm{c}}$) superconductivity desirable. However, it has been known for decades that weak electron-phonon coupling gives rise to low values of $T_{\mathrm{c}}$, while strong electron-phonon coupling leads to lattice instabili… ▽ More Electron-lattice interactions play a prominent role in quantum materials, making a deeper understanding of direct routes to phonon-mediated high-transition-temperature ($T_{\mathrm{c}}$) superconductivity desirable. However, it has been known for decades that weak electron-phonon coupling gives rise to low values of $T_{\mathrm{c}}$, while strong electron-phonon coupling leads to lattice instability or formation of bipolarons, generally assumed to be detrimental to superconductivity. Thus, the route to high-$T_{\mathrm{c}}$ materials from phonon-mediated mechanisms has heretofore appeared to be limited to raising the phonon frequency as in the hydrogen sulfides. Here we present a simple model for phonon-mediated high-$T_{\mathrm{c}}$ superconductivity based on superfluidity of light bipolarons. In contrast to the widely studied Holstein model where lattice distortions modulate the electron's potential energy, we investigate the situation where lattice distortions modulate the electron hopping. This physics gives rise to small-size, yet light bipolarons, which we study using an exact sign-problem-free quantum Monte Carlo approach, demonstrating a new route to phonon-mediated high-$T_\mathrm{c}$ superconductivity. We find that $T_\mathrm{c}$ in our model generically and significantly exceeds typical upper bounds based on Migdal-Eliashberg theory or superfluidity of Holstein bipolarons. The key ingredient in this bipolaronic mechanism that gives rise to high $T_\mathrm{c}$ is the combination of light mass and small size of bipolarons. Our work establishes principles towards the design of high-$T_{\mathrm{c}}$ superconductors via functional material engineering. △ Less

Submitted 1 February, 2023; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: 6 pages main text + 12 pages appendices, 5 figures main text + 11 figures appendices

Journal ref: Phys. Rev. X 13, 011010 (2023)

Cumulant methods for electron-phonon problems. II. The self-consistent cumulant expansion

Authors: Paul J. Robinson, Ian S. Dunn, David R. Reichman

Abstract: In this work we present a self-consistent cumulant expansion (SC-CE) and investigate its accuracy for the one-dimensional Holstein model with and without phonon dispersion. We show that for finite lattices sizes, the numerical integration of the SC-CE equations becomes unstable at long times. This defect is partially ameliorated when studying systems in the thermodynamic limit, enabling the demons… ▽ More In this work we present a self-consistent cumulant expansion (SC-CE) and investigate its accuracy for the one-dimensional Holstein model with and without phonon dispersion. We show that for finite lattices sizes, the numerical integration of the SC-CE equations becomes unstable at long times. This defect is partially ameliorated when studying systems in the thermodynamic limit, enabling the demonstration that the SC-CE corrects many deficits of the standard perturbative CE in the (non-dispersive) Holstein model. The natural phonon damping that arises in the more realistic dispersed Holstein model renders the SC-CE stable, allowing for a complete assessment of the method. Here we find that self-consistency dramatically corrects many of the failures found in the perturbative CE, but also introduces some unphysical features. Finally, we comment on the potential use of SC-CE as a tool for calculating Green's functions in generic many-body problems. △ Less

Submitted 8 June, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 11 pages, 4 figures

Journal ref: Phys. Rev. B 105, 224305 (2022)

Cumulant methods for electron-phonon problems. I. Perturbative expansions

Authors: Paul J. Robinson, Ian S. Dunn, David R. Reichman

Abstract: In this work we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CE for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as w… ▽ More In this work we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CE for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as well as the approach to the thermodynamic limit, drawing distinctions and connections between the behavior of systems in and away from the thermodynamic limit that enable a greater understanding of the ability of the CE to capture real-frequency information across the full range of wave vectors. We find that for zero electronic momentum, the spectral function is well described by the second-order CE at low and high temperatures. However, for non-zero electronic momenta, the CE is only accurate at high temperature. We analyze the fourth-order cumulant, and find that while it improves the description of the short-time dynamics encoded in the one-particle Green's function, it can introduce divergences in the time domain as well as unphysical negative spectral weight in the spectral function. When well-behaved, the fourth-order CE does provide notable accurate corrections to the second-order CE. Finally, we use our results to comment on the use of the CE as a tool for calculating transport behavior in the realistic ab initio modeling of materials. △ Less

Submitted 8 June, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 16 pages, 7 figures

Journal ref: Phys. Rev. B 105, 224304 (2022)

A Localized-Orbital Energy Evaluation for Auxiliary-Field Quantum Monte Carlo

Authors: John L. Weber, Hung Vuong, Pierre A. Devlaminck, James Shee, Joonho Lee, David R. Reichman, Richard A. Friesner

Abstract: Phaseless Auxiliary-Field Quantum Monte Carlo (ph-AFQMC) has recently emerged as a promising method for the production of benchmark-level simulations of medium to large-sized molecules, due to its accuracy and favorable polynomial scaling with system size. Unfortunately the memory footprint of standard energy evaluation algorithms are non-trivial, which can significantly impact timings on graphica… ▽ More Phaseless Auxiliary-Field Quantum Monte Carlo (ph-AFQMC) has recently emerged as a promising method for the production of benchmark-level simulations of medium to large-sized molecules, due to its accuracy and favorable polynomial scaling with system size. Unfortunately the memory footprint of standard energy evaluation algorithms are non-trivial, which can significantly impact timings on graphical processing units (GPUs) where memory is limited. Previous attempts to reduce scaling by taking advantage of the low rank structure of the Coulombic integrals have been successful, but are significantly limited by high prefactors, rendering the utility limited to very large systems. Here, we present a complementary, cubic scaling route to reduce memory and computational scaling based on the low rank of the Coulombic interactions between localized orbitals, focusing on the application to phaseless AFQMC. We show that the error due to this approximation, which we term Localized Orbital AFQMC (LO-AFQMC), is systematic and controllable via a single variable, and is computationally favorable even for small systems. We present results demonstrating a robust retention of accuracy versus both experiment and full ph-AFQMC for a variety of test cases chosen for their potential difficulty for localized orbital based methods, including the singlet-triplet gaps of polyacenes benzene through pentacene, the heats of formation for a set of platonic hydrocarbon cages, and the total energy of ferrocene (Fe(Cp)$_2$). Finally, we reproduce our previous result of the gas phase ionization energy of Ni(Cp)$_2$, agreeing with full ph-AFQMC to within statistical error while using less than a fifteenth of the computer time. △ Less

Submitted 4 April, 2022; v1 submitted 1 February, 2022; originally announced February 2022.

Comments: 36 pages, 8 figures; Supplemental Info, 11 pages, 1 figure

Time dependent second order Green's function theory for neutral excitations

Authors: Wenjie Dou, Joonho Lee, David R. Reichman, Roi Baer, Eran Rabani

Abstract: We develop a time dependent second order Green's function theory (GF2) for calculating neutral excited states in molecules. The equation of motion for the lesser Green's function (GF) is derived within the adiabatic approximation to the Kadanoff-Baym (KB) equation using the second order Born approximation for the self-energy. In the linear response regime, we recast the time dependent KB equation… ▽ More We develop a time dependent second order Green's function theory (GF2) for calculating neutral excited states in molecules. The equation of motion for the lesser Green's function (GF) is derived within the adiabatic approximation to the Kadanoff-Baym (KB) equation using the second order Born approximation for the self-energy. In the linear response regime, we recast the time dependent KB equation into a Bethe-Salpeter-like equation (GF2-BSE), with a kernel approximated by the second order Coulomb self-energy. We then apply our GF2-BSE to a set of molecules and atoms and find that GF2-BSE is superior to configuration interaction with singles (CIS) and/or time dependent Hartree-Fock (TDHF), particularly for charge transfer excitations, and is comparable to CIS with perturbative doubles (CIS(D)) in most cases. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Time Evolution of ML-MCTDH Wavefunctions II: Application of the Projector Splitting Integrator

Authors: Lachlan P. Lindoy, Benedikt Kloss, David R. Reichman

Abstract: The multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) approach suffers from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of a regularization parameter. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MC… ▽ More The multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) approach suffers from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of a regularization parameter. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to $10^6$ bath modes, and find that for these problems the PSI requires roughly 3-4 orders of magnitude fewer Hamiltonian evaluations and 2-3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH, and 2-3/1-2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to $\sim 1.3\times 10^9$ parameters to obtain accurate dynamics. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: Main paper: 17 pages, 7 figures. Supplemental information (for paper I and paper II of the series): 22 pages, 2 figures

Time Evolution of ML-MCTDH Wavefunctions I: Gauge Conditions, Basis Functions, and Singularities

Authors: Lachlan P. Lindoy, Benedikt Kloss, David R. Reichman

Abstract: We derive a family of equations of motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than that used for standard ML-MCTDH. These alternative conditions re… ▽ More We derive a family of equations of motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than that used for standard ML-MCTDH. These alternative conditions result in an expansion of the wavefunction in terms of a set of potentially arbitrary orthonormal functions, rather than in terms of a set of non-orthonormal and potentially linearly dependent functions, as is the case for standard ML-MCTDH. We show that the EOMs used in the projector splitting integrator (PSI) and the invariant EOMs approaches are two special cases of this family obtained from different choices for the dynamic gauge condition, with the invariant EOMs making use of a choice that introduces potentially unbounded operators into the EOMs. As a consequence, all arguments for the existence of parallelizable integration schemes for the invariant EOMs can also be applied to the PSI EOMs. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 10 + $ε$ pages, 2 figures

Strongly correlated ladders in K-doped $p$-terphenyl crystals

Authors: John Sous, Natalia A. Gadjieva, Colin Nuckolls, David R. Reichman, Andrew J. Millis

Abstract: Potassium-doped terphenyl has recently attracted attention as a potential host for high-transition-temperature superconductivity. Here, we elucidate the many-body electronic structure of recently synthesized potassium-doped terphenyl crystals. We show that this system may be understood as a set of weakly coupled one-dimensional ladders. Depending on the strength of the inter-ladder coupling the sy… ▽ More Potassium-doped terphenyl has recently attracted attention as a potential host for high-transition-temperature superconductivity. Here, we elucidate the many-body electronic structure of recently synthesized potassium-doped terphenyl crystals. We show that this system may be understood as a set of weakly coupled one-dimensional ladders. Depending on the strength of the inter-ladder coupling the system may exhibit spin-gapped valence-bond solid or antiferromagnetic phases, both of which upon hole doping may give rise to superconductivity. This terphenyl-based ladder material serves as a new platform for investigating the fate of ladder phases in presence of three-dimensional coupling as well as for novel superconductivity. △ Less

Submitted 20 August, 2021; originally announced August 2021.

Comments: 6 pages main text + 3 pages supplementary information, 4 figures main text

Journal ref: Nano Lett. 21, 9573 (2021)

Unbiasing Fermionic Quantum Monte Carlo with a Quantum Computer

Authors: William J. Huggins, Bryan A. O'Gorman, Nicholas C. Rubin, David R. Reichman, Ryan Babbush, Joonho Lee

Abstract: Many-electron problems pose some of the greatest challenges in computational science, with important applications across many fields of modern science. Fermionic quantum Monte Carlo (QMC) methods are among the most powerful approaches to these problems. However, they can be severely biased when controlling the fermionic sign problem using constraints, as is necessary for scalability. Here we propo… ▽ More Many-electron problems pose some of the greatest challenges in computational science, with important applications across many fields of modern science. Fermionic quantum Monte Carlo (QMC) methods are among the most powerful approaches to these problems. However, they can be severely biased when controlling the fermionic sign problem using constraints, as is necessary for scalability. Here we propose an approach that combines constrained QMC with quantum computing tools to reduce such biases. We experimentally implement our scheme using up to 16 qubits in order to unbias constrained QMC calculations performed on chemical systems with as many as 120 orbitals. These experiments represent the largest chemistry simulations performed on quantum computers (more than doubling the size of prior electron correlation calculations), while obtaining accuracy competitive with state-of-the-art classical methods. Our results demonstrate a new paradigm of hybrid quantum-classical algorithm, surpassing the popular variational quantum eigensolver in terms of potential towards the first practical quantum advantage in ground state many-electron calculations. △ Less

Submitted 5 July, 2021; v1 submitted 30 June, 2021; originally announced June 2021.

Bond-Peierls polaron: Moderate mass enhancement and current-carrying ground state

Authors: Matthew R. Carbone, Andrew J. Millis, David R. Reichman, John Sous

Abstract: We study polarons in the one-dimensional Bond-Peierls electron-phonon model, in which phonons on bonds of a lattice modulate the hopping of electrons between lattice sites, and contrast the results to those known for the Breathing-Mode Peierls problem. By inspecting the atomic limit, we show that polaronic dressing and mass enhancement of Bond-Peierls polarons depend on the momentum dependence of… ▽ More We study polarons in the one-dimensional Bond-Peierls electron-phonon model, in which phonons on bonds of a lattice modulate the hopping of electrons between lattice sites, and contrast the results to those known for the Breathing-Mode Peierls problem. By inspecting the atomic limit, we show that polaronic dressing and mass enhancement of Bond-Peierls polarons depend on the momentum dependence of the phonons. For dispersionless phonons, Bond-Peierls polarons are perfectly localized in the atomic limit, unlike their Breathing-mode counterparts, because the carrier creates a string of phonon excitations that can only be annihilated via processes that retrace the carrier to its original site. However, inclusion of phonon dispersion leads to a non-divergent polaron mass even in the atomic limit and depending on the form of the phonon dispersion may lead to a transition to a non-zero-momentum ground state akin to that found in the Breathing-mode Peierls model. △ Less

Submitted 22 November, 2021; v1 submitted 26 May, 2021; originally announced May 2021.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. B 104, L140307 (2021)