2026-05-08 | | Total: 14
The electronic structure calculations remain a major bottleneck in ab initio nonadiabatic molecular dynamics. We develop an efficient TDDFT-based FSSH implementation in the GPU4PySCF package for medium-sized molecular systems. Our approach combines density fitting, TDDFT with minimal auxiliary basis sets (TDDFT-ris), and an approximate Z-vector solver to reduce the computational cost of TDDFT excited states and derivative coupling calculations. These approximations introduce negligible errors in realistic FSSH workloads while maintaining high computational efficiency. Benchmark results show that, for 73-atom systems with a triple-$ζ$ basis set, individual electronic structure calculations are completed within one minute on a single NVIDIA A100 GPU.
Multiconfiguration pair-density functional theory (MC-PDFT) offers an efficient and accurate framework for computing electronic energies in strongly correlated molecular systems, with the quality of the on-top functional being a key determinant of its predictive accuracy. Here we introduce FunctionalAgent, an agentic system for fully automated functional development. FunctionalAgent orchestrates a team of specialized sub-agents to decompose the development process into dataset construction, active-space generation, MCSCF calculation and descriptor generation, loss-function construction, and functional fitting, optimization, and evaluation, thereby linking all stages into a closed-loop automated workflow. Using FunctionalAgent, we developed MC26, a hybrid meta-GGA on-top functional that achieves improved overall accuracy on the training set compared with other methods evaluated on the same benchmark dataset. We further introduce COF26, a new functional form that, owing to the optimized training process, achieves the best performance on both the training and test sets.
Accurate prediction of adiabatic $0$-$0$ excited-state energies is crucial for modeling molecular photophysical processes. Here, we benchmark computational strategies for evaluating excited-state energies and singlet-triplet gaps obtained using different geometry-optimization strategies, including time-dependent density functional theory (TDDFT), spin-unrestricted Kohn-Sham (UKS) DFT for triplet states (${\rm T}_1$), and state-specific orbital-optimized UKS (ssUKS) DFT for singlet excited states (${\rm S}_1$). Zero-point vibrational energy corrections are evaluated consistently at the optimized geometries and combined with ADC(2) excitation energies for comparison with experimental anion photoelectron spectroscopy data for a representative set of molecules. Among the protocols considered, adiabatic $0$-$0$ energies evaluated at TDDFT-optimized ${\rm S}_1$ and ${\rm T}_1$ geometries show the best agreement with experiment, with a mean absolute error below 0.1 eV. Replacing these geometries with UKS-optimized ${\rm T}_1$ and ssUKS-optimized ${\rm S}_1$ structures yields comparable accuracy. Vertical excitation energies are substantially more sensitive to the choice of geometry than the corresponding ${\rm S}_1$-${\rm T}_1$ gaps, which are comparatively more robust because of partial error cancellation. As a larger case study, we examine rubrene and find that UKS/ssUKS-based geometries remain useful for evaluating singlet-fission energetics. Overall, UKS/ssUKS-based workflows provide an efficient and accurate route to excited-state geometry optimization and to the evaluation of adiabatic $0$-$0$ energies for states with dominant single-determinant character.
Recent reports suggest that chemical reaction rates can change when reactants are placed inside an optical cavity. These effects have been attributed to the hybridization of molecular vibrational modes with cavity modes into polaritons, but the underlying mechanism remains debated. Recently, attempts to reproduce the key experiments have sometimes failed, which poses also ambiguity and impedes the determination of the possible mechanism. Without a reliable theoretical framework, polaritonic chemistry -- which seeks to use optical resonators as catalysts to control reactions -- has reached a pivotal stage. Standardized protocols for reproducible cavity experiments are therefore urgently needed. Here, we identify pitfalls in approaches that monitor reaction progress with UV/Vis spectroscopy. Using the Transfer Matrix Method, we analyze a model pseudo-first-order reaction and assess how transient cavity thickness variations, cavity inhomogeneity, and fitting protocols influence the extracted rate constant. We find that changes in cavity thickness upon reactant introduction can strongly distort apparent kinetics when monitoring at a single wavelength, an artifact that can be mitigated by spectral smoothing. Additionally, we demonstrate that, unlike in many previous studies, the asymptotic extinction should be treated as a fitting parameter rather than fixed to the final experimental value. By identifying these pitfalls, our work lays the foundation for more robust analyses and reliable measurements in polaritonic chemistry.
Machine learning interatomic potentials (MLIPs) offer first-principles accuracy with reduced computational cost, but their transferability across different thermodynamic states remains questionable, particularly for fluid systems where molecules experience local environments far from crystalline equilibrium. Here, we demonstrate that diversifying the density of training configurations, rather than temperature, is the most effective strategy for building thermodynamically transferable MLIPs within a fixed computational budget. We first show that foundation MLIPs trained on solid-state databases accurately describe liquid-like densities but fail at gas-like conditions, while molecular-database-trained models exhibit the opposite behavior. Controlled from-scratch training and distillation experiments confirm that density-diverse datasets resolve both failure modes, whereas temperature-diverse datasets cannot compensate for missing density regimes. Coordination number analysis reveals the physical origin of this behavior: local coordination topology is more susceptible to density than temperature, leading to further structural diversity. These results establish density diversity as a design principle for thermodynamically transferable MLIPs and provide a validation framework for assessing the thermodynamic coverage of both foundation and from-scratch models, enabling reliable atomistic simulation of fluid-phase processes across diverse operating conditions.
Iridium-based photosensitizers have attracted significant attention in photodynamic therapy (PDT) due to their exceptional photophysical properties and chemical stability, as well as tunable phosphorescence emission spectrum and high triplet state production yields. Photosensitizers with large two-photon absorption (TPA) and mitochondrial targeting capabilities are particularly promising for clinical PDT, as they enable deeper tissue penetration and reduced damage to normal cells. In this study, we theoretically studied photophysical, photodynamic properties and photosensitization reaction mechanism of a series of iridium-based photosensitizers with modified C^N and N^N ligands (a2-a6, b1/b1-r and b2/b2-r) by TDDFT/DFT methods. The photophysical properties, including one- and two-photon absorption spectra, frontier molecular orbitals, and singlet and triplet excitation energies, were calculated. Additionally, rate constants for intersystem crossing, fluorescence, and phosphorescence, along with water solubility and lipophilicity metrics (logP), were determined to assess both efficacy and biocompatibility. The results elucidate the modulation roles of the chelated ligands and ancillary ligands in TP-PDT efficiency, indicating that the asymmetric iso-fused-benzene ring modification to the N^N ligand is a robust design strategy for comprehensively enhancing photosensitization performance. Complexes a2, b2 and b1-r show greater promise as candidates for two-photon PDT photosensitizers, owing to their large TPA cross-sections, extended triplet state lifetimes, and balanced water solubility and lipophilicity. Notably, the b1-r complex can undergo both Type I and Type II PDT photosensitization mechanisms, which will help address the issue of drug resistance arising from the hypoxic environment in deep-seated tumors.
We propose a mixed quantum-classical hydrodynamic framework to model short-time inertial effects in the non-adiabatic evolution of a quantum solute coupled to a classical polar solvent. Drawing upon the work of Burghardt and Bagchi [Chem. Phys. 329 (2006), 343], we employ the Hamiltonian approach to incorporate consistent backreaction and preserve quantum decoherence beyond standard Ehrenfest dynamics. The solvent is treated as an ideal polar fluid and the quantum solute state is correlated to both the position and molecular orientation coordinates of the liquid. This approach retains essential solute-solvent correlations while significantly reducing the computational complexity of previous approaches. We further incorporate dissipative terms to capture both inertial effects and polarization relaxation. After establishing the general setting for non-local dielectric continua, the Marcus local approximation is integrated into the model thereby extending traditional solvation theory to account for collective fluid sloshing on fast timescales.
The fluctuations of ions in polar solvents remain poorly understood theoretically due to the complex coupling between ionic motion and solvent polarization. Indeed, while all-atom resolution can be achieved in numerical simulations, analytical approaches require suitable levels of coarse-graining. In this work, we describe ions and solvent molecules as interacting Brownian particles and use stochastic density functional theory to derive a generalized Langevin equation for the ionic charge density, explicitly accounting for solvent-mediated memory effects. In the regime where there is a clear timescale separation between fast solvent and slow ion dynamics, we obtain simple expressions for dynamical charge structure factors, which are validated by BD simulations. For slow solvents, we predict an emerging two-step relaxation in ionic dynamics. These results provide a mesoscopic approach for ion-solvent dynamics and open pathways to study fluctuation-induced phenomena in electrolytes.
Quantum simulation has begun to penetrate the field of quantum chemistry in hopes of efficiently calculating ground state energies and approximating real-time evolution. With modern research highlighting nonadiabatic dynamics, tunably approximating deep circuits representing potential landscapes becomes crucial for simulating real quantum systems. Variationally approximating unitaries allows for shallower circuits and accuracy tunable to hardware fidelity, so long as the observable quantities are preserved. We show the variational compression of Trotter terms preserve reaction rate coefficients via classical emulation of a hybrid quantum-classical optimization method, as well as fast-forwarded adiabatic dynamics on quantum hardware. Compressed circuits can be incorporated with product-formula-based time evolution to approximate dynamics of a particle in two coupled harmonic potentials, allowing tunability when removing high-cost qubit interactions. Approximate rate coefficients are recovered after substituting terms in a nonadiabatic dynamic process, giving proof-of-principle for observable preservation under variational optimization. Attention is paid to minimizing qubit and gate-count resources.
Conserved quantities increasingly underpin the inference of physical models. Recently new conserved quantities have been found in this context, that currently lack an interpretation. Here, we show that irreversible reactions in CRNs and Markov Chains lead to emergent conservation laws and broken cycles. Linearly dependent currents - characterized by the "co-production index" - arise due to irreversible reactions. We derive a law relating conserved quantities, broken cycles, and co-production. This resolves a recent conundrum posed by a machine-discovered candidate for a non-integer conservation law. Our findings introduce heretofore overlooked extensions to a widely used index law for CRNs and Markov Chains that undercounts conservation laws. This furnishes new tools and immediate applications for the inference and analysis of models based on conservation laws.
Long-range electrostatics and polarization remain central obstacles to extending machine learning interatomic potentials (MLIPs) to ionic, polar, and interfacial systems. Here, we introduce a semi-local framework for learning electrostatics from energies and forces using polarizable atomic multipoles. Local equivariant descriptors predict environment-dependent latent monopoles, dipoles, and quadrupoles, while residual non-local charge transfer and polarization are captured by non-self-consistent linear response in induced charges and dipoles. Across four diverse benchmarks and four short-range MLIP architectures, the multipole hierarchy and response terms systematically improve potential energy surface accuracy, with the largest gains in systems where long-range effects are essential. More importantly, the learned latent variables recover physically meaningful electrical responses: accurate Born effective charge tensors, emergent polarizabilities, infrared spectra in close agreement with experiments, and semi-quantitative Raman spectra for bulk water and hybrid MAPbI$_3$ perovskite. This systematically improvable, physically transparent framework enables MLIPs trained on standard energy and force labels to predict polarization-sensitive observables.
We introduce a scalable, translationally invariant variational theory for ab initio polarons that remains applicable across coupling regimes without resorting to supercells. Our approach combines a momentum-projected Toyozawa-type wavefunction with a low-rank factorization of the electron-phonon kernel, enabling near-linear scaling with the number of $\mathbf{k}$-points while capturing both delocalized and self-trapped carriers. Benchmarks for the Fröhlich model, LiF, and anatase and rutile TiO$_2$ yield accurate polaron binding energies, thermodynamic-limit band structures, and transparent real-space measures of polaron extent. For LiF, comparison with first-principles diagrammatic Monte Carlo (DiagMC) reveals close agreement for the weak-coupling electron-polaron ground state and band structure. However, in the hole-polaron of LiF, which is in the strong-coupling regime, we found a significant bias in DiagMC results. These results establish momentum-projected variational wavefunctions as a systematically improvable route to thermodynamic limit studies of polarons in real materials.
The development of accurate exchange-correlation (XC) functionals remains a longstanding challenge in density functional theory (DFT). The vast majority of XC functionals have been hand designed by human researchers combining physical insight, exact constraints, and empirical fitting. Recent advances in large language models enable a systematic, automated alternative to this human-driven design loop. This report presents an agentic search system in which an LLM proposes structured functional-form changes guided by evolutionary history. The system attempts to improve functional performance through an iterative plan-execute-summarize loop, where improvements are measurable by optimizing functional parameters against a standard thermochemistry dataset, then evaluating performance on a held-out subset. The strongest discovered functional, SAFS26-a (Seed Agentic Functional Search 2026), improves upon the gold-standard ωB97M-V baseline by ~9%. These results also surface a cautionary lesson for AI-assisted science: models powerful enough to discover genuine improvements are equally capable of exploiting unphysical shortcuts to game the benchmark; domain expertise translated into explicitly enforced constraints remains essential to keeping results scientifically grounded.
CDFCI is a shared-memory parallel numerical program for computing low-lying eigenpairs of large-scale, non-relativistic fermionic Hamiltonians. The software is designed to handle a broad class of many-body quantum models, including both ab initio electronic structure Hamiltonians and lattice-based Hamiltonians arising in condensed matter physics. CDFCI combines an efficient coordinate-descent-based selected configuration interaction algorithm with dedicated parallelization strategies, achieving high performance on modern multi-core architectures. Benchmark results on representative quantum chemistry and condensed matter test cases demonstrate that CDFCI attains state-of-the-art accuracy with competitive performance compared to established selected configuration interaction (such as CIPSI or SHCI) and DMRG implementations. The software is open-source, extensively documented, and provides a Python interface for seamless integration with PySCF and other many-body simulation workflows.