2026-04-15 | | Total: 15
We introduce a new equation-of-motion coupled-cluster method based on a pair coupled-cluster doubles (pCCD) reference, termed frozen-pair EOM-CCSD (EOM-fpCCSD). This approach combines the computational efficiency of the pCCD ansatz with a dynamical correlation correction, enabling a reliable description of electronically excited states within the EOM framework. The method has been implemented in the open-source PyBEST software package. Its performance is systematically benchmarked against standard EOM-CCSD and its pair-tailored variant (EOM-ptCCSD), using both canonical Hartree-Fock and pCCD natural orbitals. For charge-transfer (CT) excitations taken from the QUEST database, EOM-fpCCSD yields excitation energies very close to those of EOM-CCSD, outperforming EOM-ptCCSD, as well as to the theoretical best estimates (TBEs). Working within the localized pCCD natural orbital basis allows us to determine the directed CT character, which quantifies the directed charge flow from one molecular domain to another. Numerical results show that EOM-fpCCSD, EOM-CCSD, and EOM-ptCCSD provide nearly identical descriptions of the directed CT character, despite changes in excitation energies. The true advantage of EOM-fpCCSD becomes evident for the challenging QUEST subset of doubly excited states. While EOM-ptCCSD performs similarly to standard EOM-CCSD, EOM-fpCCSD significantly outperforms both methods for these problematic states compared to TBEs. In addition to improving the accuracy of excitation energies, EOM-fpCCSD also converges for several states that standard EOM-CCSD and EOM-ptCCSD fail to converge. These results demonstrate that EOM-fpCCSD offers a promising and computationally efficient route toward a more accurate description of complex electronic excitations.
We present an efficient implementation of the low-cost linear-response coupled-cluster singles and doubles (LR-CCSD) method for computing static and frequency-dependent polarizabilities in systems with significant relativistic and electron-correlation effects. The approach employs X2C-based Hamiltonians (X2CAMF and X2CMP) and incorporates Cholesky decomposition to reduce memory requirements. In the current implementation, costly three- and four-external index integrals are generated on the fly, eliminating the need for their storage. Benchmark results indicate that the X2CMP Hamiltonian provides more consistent performance than X2CAMF, particularly for large and highly augmented basis sets. The proposed FNS++CD-X2CMP-LR-CCSD method shows excellent agreement with four-component reference values across a wide range of systems. Additionally, different strategies for constructing the FNS++ basis were assessed, and an averaged density approach was found to offer a favorable balance between accuracy and computational cost. On average, about 73% of the virtual spinor space is removed, demonstrating the efficiency and consistency of the FNS++ density-based truncation approach. The present implementation enables accurate and scalable relativistic response calculations for large molecular systems, as demonstrated by the calculation of the static polarizability of the Uranium Hexafluoride complex with a triple-zeta basis set more than 1400 basis functions.
There has been a veritable explosion of methods and software to perform machine-learned regression on datasets of electronic energies and forces to develop high-dimensional machine learned potential energy surfaces (ML-PESs). A major, but not deeply-studied aspect is how well different ML-PESs represent the same dataset on which they are trained, beyond the standard fitting precision metrics. Here, this is examined in detail using several ''stress tests'', for two widely applied machine-learned potential approaches. One is based on permutationally invariant polynomial (PIP) linear least square regression and the other is the message-passing neural network PhysNet approach. These potentials and dipole moment surfaces are used in VSCF/VCI calculations of vibrational energies and wavefunctions. The energies from the two PESs are directly compared as are the IR spectra. In addition, tunneling splittings for the hydrogen transfer between two equivalent structures are reported from using three methods: ring polymer instanton theory, diffusion Monte Carlo simulations, and the $Q_{im}$ path method. These calculations require the evaluation of on the order of one billion energies that are widely dispersed in the 15-dimensional configurational space. The two PESs yield results for these quantities in excellent agreement with each other.
Gas hydrates are potential candidates in future energy sources while simultaneously providing structures with extensive applications in carbon capture and storage, gas transport, and important separation processes. Prior research in the field considers the dynamics of the water molecule backbone in particular. We investigated the pressure-enthalpy landscape and mechanical stability envelope of sI methane and carbon dioxide hydrates simulated using DFT. We investigated the effect of the revPBE + DFT-D2 and the SCAN + rVV10 and their treatment of the exchange correlation interactions. We examined the zero pressure material properties, finding that revPBE comparatively underbinds the interactions, causing more flexible structures with large equilibrium volumes. Under pressure, the carbon dioxide molecule was found to align itself parallel to the hexagonal faces of the large cage despite the functional used. Additionally, the property differences are caused by the ability of the carbon dioxide molecule to rotate and disperse the changes in the energy landscape in ways that methane molecules cannot. This computational methodology describes the elastic stability of gas hydrate, marginal stability, and critical differences across important molecular interactions, confirming experimentally observed restrictions in guest molecule rotations and novel pressure behaviors under hydrostatic loads
This Perspective is part of a Special Topic that explored the maturity of nonadiabatic molecular dynamics for predicting photochemical processes. In 2023, a prediction challenge was issued to the community of computational photochemists to simulate the photochemistry of cyclobutanone, photoexcited at 200 nm, and the resulting time-resolved MeV-UED signal. The challenge attracted 15 theoretical predictions from more than 70 researchers, employing a wide range of strategies for electronic structure and nonadiabatic molecular dynamics to predict the time-resolved MeV-UED signal before the experiment had been conducted at SLAC (Stanford, USA). The MeV-UED instrument at Shanghai Jiao Tong University was also used to provide a second independent time-resolved MeV-UED signal for the photochemistry of cyclobutanone. This Perspective discusses the various approaches and strategies used by the participants to predict the photochemistry of cyclobutanone. This work also summarizes the strengths and weaknesses of various methods used for photoexcitation, electronic structure, nonadiabatic dynamics, and calculation of observables, as agreed by the participants during a CECAM workshop dedicated to the results of the challenge and organized in Lausanne in April 2025. This Perspective also collects all the predicted time-resolved MeV-UED signals into a single figure, together with the experimental signal. This challenge (i) demonstrated the qualitative predictive power of nonadiabatic molecular dynamics and (ii) underscore the impact of electronic-structure theory on the outcome of the excited-state dynamics and the need for its careful benchmarking. This effort allowed the community to share practical strategies to perform nonadiabatic dynamics (discussed in the present Perspective) and constitutes a 'calibration' exercise for computational photochemistry.
Transferable excited-state dynamics offer a route to efficient screening of photophysical behavior across molecular systems, but conventional nonadiabatic simulations remain prohibitively expensive. Here we introduce X-MACE, a transferable machine-learning potential for excited-state dynamics that predicts multiple potential energy surfaces, forces and oscillator strengths, and combine it with curvature-driven surface hopping to enable data-efficient screening of photochemical pathways. We apply this framework to fluorescent chromophores as an example application, using green fluorescent protein chromophore variants to demonstrate how subtle structural modifications reshape excited-state relaxation, lifetimes and photoisomerization yields. Fine-tuning a single pretrained model with fewer than 100 reference geometries per derivative yields accurate dynamics across a chemically diverse set of analogues. The screening reveals two governing design principles: steric crowding on the phenolate ring lowers the torsional barrier and accelerates access to twisted conical intersections, whereas conjugation extension stabilizes planar excited-state configurations, suppresses non-radiative decay and prolongs fluorescence. More broadly, this workflow provides a general framework for scalable excited-state screening and interpretable design of photophysical properties.
A path-integral hybrid Monte Carlo approach with enveloping bridging potentials (PIHMC-EBP) is proposed for calculating numerically exact tunneling splittings in molecular systems. The central idea is to construct an approximately barrierless bridging potential that smoothly connects symmetry-related regions of ring-polymer phase space, enabling direct sampling of the free-energy profile from which the relevant splittings are obtained. Two tailored nonlocal updates are designed to enhance the sampling of slow collective motions. Compared with path-integral molecular dynamics using thermodynamic integration, PIHMC-EBP requires neither quadrature nor time-step convergence checks, thereby substantially reducing the manual effort required to analyze the results. Applications to malonaldehyde (and its deuterated isotopologue) and the HCl dimer using state-of-the-art potential energy surfaces provide the most precise tunneling splittings reported to date for both systems, while simultaneously reducing the overall computational cost by several times and three orders of magnitude, respectively. Finally, application to the water dimer yields the first numerically exact path-integral calculations of the ground-state tunneling splittings on three different potential energy surfaces, all obtained simultaneously by reweighting a single set of trajectories.
We extend our previous symmetrized path-integral molecular dynamics approach to calculate tunneling splittings of molecules in rotationally excited states. In this new formalism, the system is rigorously projected onto selected rotational manifolds and states of a chosen symmetry through an Eckart spring, which connects the two end beads of the ring polymer via a permutation--inversion--rotation operation. This method is numerically exact within statistical uncertainty once convergence with respect to all simulation parameters has been achieved. Importantly, it enables the simultaneous extraction of tunneling splittings for multiple total angular-momentum quantum numbers $J$ from a single set of simulations, without additional computational cost relative to the original approach. After validating the formalism by computing the rotational levels of water (beyond the rigid-rotor approximation), we apply it to ammonia and obtain rotationally resolved tunneling splittings in excellent agreement with exact variational benchmarks. Except for small errors due to the underlying potential energy surface, the results capture the experimentally observed trend that the tunneling splitting decreases with $J$.
The theory of hyper-Rayleigh scattering optical activity (HRS-OA) spectroscopy has previously been described within the length formulation of the pure electric-dipole and mixed (electric-dipole/magnetic-dipole and electric-dipole/electric-quadrupole) first hyperpolarizabilities required for the description of this process. In this work, we provide an alternative formulation of these pure and mixed hyperpolarizabilities. This new formulation made use of the velocity form of the electric-dipole and electric-quadrupole moment operators that enter in the quadratic response functions. A one-to-one correspondence is found for the gauge-origin shifts obtained in the two formulations. These relations ensure the origin-independence of the theory also for the velocity formulation. Furthermore, even though the basis set dependence of the velocity formulation is more significant compared to the length one, the former is origin-independent by design. This property makes it particularly suitable for calculations of HRS-OA invariants using approximated (variational or not) wavefunctions.
The interest to foster plasmonic applications at energies in the ultra-violet, has escalated research initiatives in clusters of unconventional plasmonic materials like aluminum and indium,for which the surface-plasmon resonance appears above the ionization potential. Naturally, the quantum mechanical description calls for the incorporation of the ionization process, thereby making the ab initio calculations challenging. We present a robust approach within the time-evolution formalism of the time-dependent density-functional theory to calculate surface plasmon resonance in the continuum of metal clusters. Using the much studied Al$_{13}^-$ as a system of reference, we show that accurate description of the continuum and of the ionization of the cluster allow to capture a broad surface-plasmon in the UV. Application of this approach in aluminum clusters has given the size-dependent evolution from discrete spectral features in Al$_{6}$ to the surface-plasmon in larger clusters in the deep ultra-violet.
Machine learning is revolutionizing chemistry. Beyond the value of predictive models accelerating virtual screening, generative AI aims at enabling inverse design, reversing the compound-to-property prediction paradigm into property-to-compound generation. Chemists now have access to a rich AI toolbox for organic chemistry, including drug discovery. However, the application of these methods to inorganic compounds remains limited by the challenges posed by their intrinsic nature. This Review analyzes how these challenges have been addressed, considering widely diverse systems ranging from molecules to crystals, including transition metal complexes and microporous materials. The analysis focuses on how generative AI methods have evolved towards data-representation-model pipelines that address the full complexity of inorganic compounds, including their chemical composition, geometry, symmetry, and electronic structure. Future directions, like benchmark standardization and the development of synthesizability metrics, are also discussed.
The functionality of catalysts, enzymes, and supramolecular assemblies emerges not from individual molecules alone, but from the subtle interplay between multiple components arranged in complex systems. Designing such systems is a grand challenge, the combinatorial explosion of possible chemical compositions and spatial arrangements makes brute-force exploration infeasible, while many current generative approaches remain limited to isolated molecules. In this work, we introduce a hierarchical generative optimization framework that overcomes this barrier by coupling a genetic algorithm for configurational search with a generative model for molecular design. This closed-loop approach enables simultaneous refinement of geometry and composition, efficiently steering discovery toward systems with targeted functionality. As a proof of concept, we design catalytic environments for the Claisen rearrangement of p-tolyl ether by optimizing surrounding components around a fixed reference transition-state geometry. Despite this constraint during the search phase, post-hoc validation via Climbing-Image Nudged Elastic Band calculations confirm a 30% reduction in activation barrier. Beyond this example, our framework provides a general strategy for data-driven discovery of functional multi-component systems, opening the door to automated design of catalysts, enzyme active sites, and advanced materials. Scientific contribution. The study presents a closed loop generative framework that enables joint optimization of molecular components and their spatial organization in multi-component systems. The method moves generative molecular design beyond single molecules toward larger and more complex systems.
Out-of-Time-Order Correlators (OTOCs) quantify quantum information scrambling, but their connection to localized phase-space structures, such as chemical transition states, requires formal development. We derive a leading-order semiclassical expansion for the local microcanonical OTOC in systems with an index-1 saddle point, expressing the scrambling rate as a coherent sum over unstable periodic orbits on the Normally Hyperbolic Invariant Manifold (NHIM). Valid in the semiclassical limit and the intermediate-time regime before the Ehrenfest time, our derivation utilizes the Normal Form theory of the transition state, which transforms the Hamiltonian near the saddle into an integrable (though generally non-separable) form dependent on conserved actions. We outline the derivation of the microcanonical trace, the semiclassical propagator for integrable systems, the factorization of the stability matrix, and the Schur complement reduction of the stationary phase approximation. Our result extends periodic-orbit trace methods to scrambling observables, yielding a local instability exponent Λ(J) governing the leading semiclassical growth window. As a special case, when the observation time coincides with the intrinsic periods of the contributing orbits, the trace sum reduces to an effective 1.5Λ scaling, resulting from the competition between local hyperbolic growth and wavepacket dilution. This simplified form is conditional; the full expansion retains a coherent sum over orbit periods. Finally, we discuss how the dependence of the instability on transverse actions establishes a theoretical mechanism for mode-selective control of scrambling, and outline a numerical evaluation strategy to test these predictions.
We present XANE(3), a physics-based E(3)-equivariant graph neural network for predicting X-ray absorption near-edge structure (XANES) spectra directly from atomic structures. The model combines tensor-product message passing with spherical harmonic edge features, absorber-query attention pooling, custom equivariant layer normalization, adaptive gated residual connections, and a spectral readout based on a multi-scale Gaussian basis with an optional sigmoidal background term. To improve line-shape fidelity, training is performed with a composite objective that includes pointwise spectral reconstruction together with first- and second-derivative matching terms. We evaluate the model on a dataset of 5,941 FDMNES simulations of iron oxide surface facets and obtain a spectrum mean squared error of $1.0 \times 10^{-3}$ on the test set. The model accurately reproduces the main edge structure, relative peak intensities, pre-edge features, and post-edge oscillations. Ablation studies show that the derivative-aware objective, custom equivariant normalization, absorber-conditioned attention pooling, adaptive gated residual mixing, and global background term each improve performance. Interestingly, a capacity-matched scalar-only variant achieves comparable pointwise reconstruction error but reduced derivative-level fidelity, indicating that explicit tensorial channels are not strictly required for low intensity error on this dataset, although they remain beneficial for capturing finer spectral structure. These results establish XANE(3) as an accurate and efficient surrogate for XANES simulation and offer a promising route toward accelerated spectral prediction, ML-assisted spectroscopy, and data-driven materials discovery.
Quantum mechanics owes much of its extraordinary success to a Hilbertian program of mathematical formalization. Yet, the formalism remains poorly aligned with the practical limitations of computations in finite dimensions and under finite accuracy. In this perspective, we argue that this mismatch points to the need for a new mathematical program: a rigorous constructive theory for effective descriptions to identify essential degrees of freedom. We propose an observation-centered point of view in which signals are treated as the primary objects of analysis, while wave functions and Hamiltonians are reconstructed as auxiliary structures to rationalize the observed data. Our starting point is a signal-based spectral equation that reformulates frequency analysis as an operator problem. We connect this point of view to results on prolate Fourier theory, spectral analysis with finite observation time, and short-time quantum simulation. We highlight a sharp accuracy transition relating necessary observation time to the effective spectral density of a signal for achieving accurate resolution. The resulting framework integrates approximation as a fundamental necessity more directly into the foundations of quantum mechanics and points toward a broader program for the effective description of complex quantum systems, such as those found in the molecular sciences.