Simulating anharmonic vibrational polaritons beyond the long wavelength approximation.

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.

Investigating the collective nature of cavity-modified chemical kinetics under vibrational strong coupling.

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

Correlation functions from tensor network influence functionals: The case of the spin-boson model.

We investigate the application of matrix product state (MPS) representations of the influence functionals (IFs) 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.

Benchmarking the exponential ansatz for the Holstein model.

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.

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

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 past 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 data sets.

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

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.

Variational Lang-Firsov Approach Plus Møller-Plesset Perturbation Theory with Applications to Ab Initio Polariton Chemistry.

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-, with = 2 and 4) perturbation theory is then applied 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 (H, HF), and the modification of proton transfer reactions (malonaldehyde and aminopropenal) via the formation of polaritons in an optical cavity.

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

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 (CG) approaches in -space somewhat ambiguous. In this paper, we address this issue by designing a perturbative expansion of pair interactions in reciprocal space.

Electron and Spin Delocalization in [Co Se (PEt ) ] Superatoms.

Molecular clusters can function as nanoscale atoms/superatoms, assembling into superatomic solids, a new class of solid-state materials with designable properties through modifications on superatoms. To explore possibilities on diversifying building blocks, here we thoroughly studied one representative superatom, Co Se (PEt ) . We probed its structural, electronic, and magnetic properties and revealed its detailed electronic structure as valence electrons delocalize over inorganic [Co Se ] core while ligands function as an insulated shell.

Instantons and the quantum bound to chaos.

The rate at which information scrambles in a quantum system can be quantified using out-of-time-ordered correlators. A remarkable prediction is that the associated Lyapunov exponent [Formula: see text] that quantifies the scrambling rate in chaotic systems obeys a universal bound [Formula: see text]. Previous numerical and analytical studies have indicated that this bound has a quantum-statistical origin.

Response properties in phaseless auxiliary field quantum Monte Carlo.

We present a method for calculating first-order response properties in phaseless auxiliary field quantum Monte Carlo by applying 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 with the same computational cost scaling per sample as energy calculations, accompanied by a cost prefactor of less than four in our numerical calculations.

Expanding the Design Space of Constraints in Auxiliary-Field Quantum Monte Carlo.

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

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

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 (), where 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. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths.

Finding defects in glasses through machine learning.

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. Due to their extremely low density, it is very hard to directly identify them in computer simulations.

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

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.

Ultrafast imaging of polariton propagation and interactions.

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.

Quantum dynamical effects of vibrational strong coupling in chemical reactivity.

Recent experiments suggest that ground state chemical reactivity can be modified when placing molecular systems inside infrared cavities where molecular vibrations are strongly coupled to electromagnetic radiation. This phenomenon lacks a firm theoretical explanation. Here, we employ an exact quantum dynamics approach to investigate a model of cavity-modified chemical reactions in the condensed phase.

Microscopic Theory of Multimode Polariton Dispersion in Multilayered Materials.

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

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond.

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within "chemical accuracy" (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N) computational cost.

Evaluating the evidence for exponential quantum advantage in ground-state quantum chemistry.

Due to intense interest in the potential applications of quantum computing, 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, for generic chemical problems where heuristic quantum state preparation might be assumed to be efficient. The availability of exponential quantum advantage then centers on whether features of the physical problem that enable efficient heuristic quantum state preparation also enable efficient solution by classical heuristics.

Data-Efficient Machine Learning Potentials from Transfer Learning of Periodic Correlated Electronic Structure Methods: Liquid Water at AFQMC, CCSD, and CCSD(T) Accuracy.

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 and provide a data-efficient approach to obtain machine-learned condensed-phase potential energy surfaces using AFQMC, CCSD, and CCSD(T) from a very small number (≤200) of energies by leveraging a transfer learning scheme starting from lower-tier electronic structure methods. We demonstrate the effectiveness of this approach for liquid water by performing both classical and path integral molecular dynamics simulations on these machine-learned potential energy surfaces.

Remembering the Work of Phillip L. Geissler: A Coda to His Scientific Trajectory.

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.

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

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.