2025-05-02 | | Total: 9
The calculation of reactive properties is a challenging task in chemical reaction discovery. Machine learning (ML) methods play an important role in accelerating electronic structure predictions of activation energies and reaction enthalpies, and are a crucial ingredient to enable large-scale automated reaction network discovery with >103 reactions. Unfortunately, the predictive accuracy of existing ML models does not yet reach the required accuracy across the space of possible chemical reactions to enable subsequent kinetic simulations that even qualitatively agree with experimental kinetics. Here, we comprehensively assess the underlying reasons for prediction failures within a selection of machine-learned models of reactivity. Models based on difference fingerprints between reactant and product structures lack transferability despite providing good in-distribution predictions. This results in a significant loss of information about the context and mechanism of chemical reactions. We propose a convolutional ML model that uses atom-centered quantum-chemical descriptors and approximate transition state information. Inclusion of the latter improves transferability for out-of-distribution benchmark reactions, making greater use of the limited chemical reaction space spanned by the training data. The model further delivers atom-level contributions to activation energies and reaction enthalpies that provide a useful interpretational tool for rationalizing reactivity.
Unconditional molecular generation is a stepping stone for conditional molecular generation, which is important in \emph{de novo} drug design. Recent unconditional 3D molecular generation methods report saturated benchmarks, suggesting it is time to re-evaluate our benchmarks and compare the latest models. We assess five recent high-performing 3D molecular generation methods (EQGAT-diff, FlowMol, GCDM, GeoLDM, and SemlaFlow), in terms of both standard benchmarks and chemical and physical validity. Overall, the best method, SemlaFlow, has a success rate of 87% in generating valid, unique, and novel molecules without post-processing and 92.4% with post-processing.
We present an implementation of relativistic double-ionization-potential (DIP) equation-of-motion coupled-cluster (EOMCC) with up to 4-hole--2-particle (4h2p) excitations that makes use of the molecular mean-field exact two-component (mmfX2C) framework. We apply mmfX2C-DIP-EOMCC to several neutral atoms and diatomic molecules to obtain the ground and first few excited states of the corresponding di-cation species, and we observe excellent agreement (to within 0.001 eV) between double ionization potentials (IPs) obtained from mmfX2C- and four-component DIP-EOMCC calculations that include 3-hole--1-particle (3h1p) excitations, with either the Dirac--Coulomb or Dirac--Coulomb--Gaunt Hamiltonians. We also compare double IPs for mmfX2C-DIP-EOMCC calculations with the full Dirac--Coulomb--Breit (DCB) Hamiltonian to those from experiment. The mmfX2C-DIP-EOMCC with 3h1p excitations leads to errors in absolute double IPs that are larger than 0.1 eV, whereas the addition of 4h2p excitations reduces these errors dramatically, often by an order of magnitude or more.
We introduce a unified statistical framework for quantifying system-environment coupling by treating the interaction energy VSE as a stochastic variable. Using a reference-particle decomposition, we derive exact, closed-form expressions for the mean and variance of VSE in terms of the single-particle density and up to four-body correlation functions. When VSE is approximately Gaussian, these two moments suffice to compute the free energy shift of the strongly coupled system. To validate our framework, we ran explicit Monte Carlo simulations of the full system-environment configurations across a range of system sizes, generating reference distributions of the interaction energy VSE. We then applied our derived analytical formulas to predict these distributions and found excellent agreement in both the weak- and strong-coupling regimes.
We present a kinematically complete study of dissociative electron attachment (DEA) in linear OCS molecules, focusing on how electrons resonantly attach and trigger dissociation. Near the Franck-Condon regime, DEA is dominated by molecular shape resonances, where transient OCS− states form with high vibrational amplitudes, spectroscopically evident as broad features in DEA cross-sections. As the electron beam energy increases from 5.5 to 6.0 eV, S− population shifts from lower to higher-energy highly dense bending vibrational states, reinforcing our findings on dipole-forbidden vibronic intensity borrowing. Our advanced potential energy curve calculations, employing the Equation-of-motion coupled cluster singles and doubles for electron attachment (EA-EOMCCSD) method, reveal that beyond the shape resonance, non-adiabatic resonant tunneling governs the avoided crossings, dynamically generating three mode-specific vibronic Feshbach resonances before complete dissociation into three distinct kinetic energy bands of S−. Our theoretical results probe most of the experimental observations quantitatively and qualitatively. These insights deepen our fundamental understanding of resonance-mediated dissociation in electron-molecule resonant scattering, with broader implications for quantum mechanics, plasma physics, vibrational revival, astrochemistry, and radiation damage research.
We derive and implement a second-order adjoint method to compute exact gradients and Hessians for a prototypical quantum optimal control problem, that of solving for the minimal energy applied electric field that drives a molecule from a given initial state to a desired target state. For small to moderately sized systems, we demonstrate a vectorized GPU implementation of a second-order adjoint method that computes both Hessians and gradients with wall times only marginally more than those required to compute gradients via commonly used first-order adjoint methods. Pairing our second-order adjoint method with a trust region optimizer (a type of Newton method), we show that it outperforms a first-order method, requiring significantly fewer iterations and wall time to find optimal controls for four molecular systems. Our derivation of the second-order adjoint method allows for arbitrary parameterizations of the controls.
In open systems with strong coupling, the interaction energy between the system and the environment is significant, so thermodynamic quantities cannot be reliably obtained by traditional statistical mechanics methods. The Hamiltonian of mean force H∗β offers an in principle accurate theoretical basis by explicitly accounting for the interaction energy. However, calculating the Hamiltonian of mean force is challenging both theoretically and computationally. We demonstrate that when the condition VarE0(e−βVSE)=0 is met, the dependence of thermodynamic variables can be shifted from {Pβ(xS),H∗β(xS)} to {Pβ(xS),P(VSE)}. This change simplifies thermodynamic measurements. As a central result, we derive a general equality that holds for arbitrary coupling strengths and from which an inequality follows - aligned with Jensen's inequality applied to the Gibbs-Bogoliubov-Feynman bound. This equality, analogous in importance to the Jarzynski equality, offers deeper insight into free energy differences in strongly coupled systems. Finally, by combining our result with said Jarzynski equality, we derive additional relations that further clarify thermodynamic behavior in strongly coupled open systems.
We introduce an efficient method TTN-HEOM for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme to the bexcitonic generalization of the numerically-exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system-bath interaction. These master equations are constructed based on the time-dependent Dirac-Frenkel variational principle which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, Tensor Equations for Non-Markovian Structured Open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: Two fixed-rank methods that require a constant memory footprint during the dynamics, and one adaptive-rank method with variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude-Lorentzian component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable to simulate both dephasing and relaxation dynamics of driven quantum system interacting with structured baths, even those of chemical complexity, with affordable computational cost.
Electronic-structure theory is the foundation of the description of materials including multiscale modeling of their properties and functions. Obviously, without sufficient accuracy at the base, reliable predictions are unlikely at any level that follows. The software package FHI-aims has proven to be a game changer for accurate free-energy calculations because of its scalability, numerical precision, and its efficient handling of density functional theory (DFT) with hybrid functionals and van der Waals interactions. It treats molecules, clusters, and extended systems (solids and liquids) on an equal footing. Besides DFT, FHI-aims also includes quantum-chemistry methods, descriptions for excited states and vibrations, and calculations of various types of transport. Recent advancements address the integration of FHI-aims into an increasing number of workflows and various artificial intelligence (AI) methods. This Roadmap describes the state-of-the-art of FHI-aims and advancements that are currently ongoing or planned.