2026-03-23 | | Total: 9
Machine learning interatomic potentials (MLIPs) trained on large, chemically diverse datasets are revolutionizing computational chemistry, enabling molecular dynamics simulations of battery electrolytes with near-DFT accuracy over 10,000 times faster than DFT. While previous MLIP training datasets with suitable elemental coverage for electrolytes have been based on inorganic materials, the Open Molecules 2025 (OMol25) dataset provides large-scale molecular DFT MLIP training data with broad elemental coverage and specifically samples tens of millions of electrolyte configurations. Here, we integrate computational modeling with experimental validation to systematically assess the ability of large-scale MLIPs pre-trained on materials data or on OMol25 to accurately resolve nanoscale structural organization and ion-solvation characteristics in Na-ion battery electrolytes across diverse physicochemical conditions and compositional regimes. We find that the OMol25-trained Universal Model of Atoms (UMA-OMol) predicts experimentally measured densities and X-ray structure factors in substantially better agreement compared to state-of-the-art models trained only on inorganic materials data. Using UMA-OMol, we further analyze systematic trends in solvation structure as a function of cation identity, anion chemistry, salt concentration, and solvent topology. We observe that increasing system temperature amplifies the heterogeneity within the solvation environment, perturbing cation-solvent interactions and promoting the formation of contact ion pairs (CIPs). Moreover, subtle variations in the solvent topology of glyme-based electrolytes cause pronounced changes in ion-correlations and solvation structure. The experimental agreement and microscopic insights shown here position OMol25-trained MLIPs as a practical route to predictive, high-throughput electrolyte simulations.
The $Δ$SCF DFT approach defines the system energy as a function of orbital occupancy. Inspired by Landau Fermi liquid theory, we develop an occupancy extrapolation (OE) method that captures excited-state energies via a Taylor expansion of the energy with respect to occupation fluctuation from a reference state. OE retains the physics of $Δ$SCF while offering a physical interpretation of excitation energies as sums of quasiparticle energies and their generalized screened interactions. It yields accurate valence, Rydberg, and charge-transfer excitation energies at $O(N^3)$ cost, avoids separate SCF calculations for each excited state, and enables efficient large-scale excited-state simulations from ground-state calculations.
We present the positron coupled cluster singles and doubles (POS-CCSD) method to calculate positron binding energies in molecules. This framework treats electrons and positrons on an equal footing and includes up to simultaneous double-electron-single-positron excitations. We benchmark the approach by computing binding energies for atomic anions and several polar and non-polar polyatomic systems, comparing the results with independent theoretical studies and, where available, experimental data. The fully converged results for H$^{-}$ are in excellent agreement with quantum Monte Carlo and multi-reference configuration interaction results. Quantitative agreement with experiments is not reached in the present study due to the slow convergence of the binding energy with respect to the size of the orbital bases for the electrons and the positron. However, the POS-CCSD results underscore the critical role of electron correlation in the description of electron-positron systems required for a balanced description of these complex systems. In addition, we examine nuclear relaxation effects following positron attachment in LiH.
Transition-metal complexes (TMCs) are promising photosensitizers for Type~I photodynamic therapy (PDT), where electron-transfer processes can generate reactive oxygen species under hypoxic conditions. Yet identifying candidates with the required ground- and excited-state redox energetics remains challenging across the vast chemical space of TMCs. Here, we develop a data-efficient active learning (AL) framework for the discovery of Type~I active TMC photosensitizers by combining a chemically structured design space of over 2.1 million Ru(II), Os(II), and Ir(III) complexes with targeted DFT calculations and pretrained atomistic representations. With only 300 quantum-chemical evaluations, the approach efficiently enriches candidates within a mechanistically defined optimal redox region. Analysis of the viable complexes reveals chemical design principles linking metal identity, ligand framework, substituent pattern, and physicochemical properties to Type~I photoreactivity, including a pronounced preference for Os(II)-based complexes and electronically asymmetric ligand environments along with combination of electronic donating and accepting substituents. More broadly, the strategy presented herein provides a scalable, mechanism-guided route for the rational design of transition-metal photocatalysts for applications spanning biomedicine, solar energy conversion, and photoredox chemistry.
The analysis of nuclear magnetic resonance parameters, such as the indirect nuclear spin-spin coupling constants, in terms of contributions from localised molecular orbitals is a commonly used approach for gaining a deeper understanding of experimentally observed trends in these parameters. In the vast majority of these studies, contributions from pairs of one occupied and one virtual orbital are calculated and analyzed. Analyses in terms of two pairs of an occupied and a virtual orbital, that would allow for the study of coupling pathways, are much more seldom, as they require calculating the coupling constants as a sum over all excited states. Previous studies have shown that, for the often dominating Fermi-contact contribution to the coupling constants, more or less all excited states have to be calculated when employing a Davidson algorithm, because the most high-lying excited states can also make a significant contribution to the Fermi-contact term. In this study we investigated therefore, whether by employing a Lanczos algorithm one can obtain converged values of the Fermi-contact contribution to the indirect nuclear spin-spin coupling constants already with a significantly smaller percentage of the total number of excited states included in the sum-over-states expression. To this purpose we have extended the recent implementation of a Lanczos algorithm for the RPA/TDHF or TDDFT eigenvalue problem in the Dalton program (L. Zamok et al. J. Chem. Phys. 156, 014102 (2022)). The new procedure was tested on 17 molecules containing first, second and third row atoms. We find that, for most coupling constants, less than 50% of the excited pseudo states are necessary for converging the Fermi-contact term with an error of less than 0.5 Hz. For the few exceptions, typically for molecules with third-row atoms, around 60% were necessary.
Pulsed Dynamic Nuclear Polarization (DNP) is currently receiving substantial interest as a means to enhance the sensitivity of nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI) by orders of magnitude. It has also received much attention as a central ingredient in many modalities of electron spin-involved quantum sensing. Relative to spin engineering associated with NMR, the design of efficient pulsed DNP experiments with a broad experimental scope are challenged by large electron-nuclear spin systems, large electron spin-involved interactions, and instrumental non-idealities and limitations. All of this may challenge traditional NMR-like theoretical and numerical pulse sequence engineering. Exploiting state-of-the-art instrumentation and taking advantage of the high sensitivity of DNP relative to NMR, we here demonstrate the use of combinations of Bayesian machine learning methods and constrained random walk procedures to design pulse sequences \textit{in situ}, by experiments, directly on the spin systems responding to spectrometer instructions. For trityl and nitroxide samples, it is demonstrated that efficient broadband DNP pulse sequences can be designed in situ with experimental protocols benchmarked against in silico analogs.
Electromagnetic (EM) coupling between molecular exciton and plasmon has been studied using in Rayleigh scattering or extinction spectroscopy. However, evaluating EM coupling involving subradiant plasmon is challenging because this resonance does not manifest clearly in far-field spectra. In this study, we developed a method to evaluate such coupling using EM enhancement factors (FR) derived from ultrafast surface-enhanced fluorescence (ultrafast SEF). This SEF, which appears as a broad background in surface-enhanced resonant Raman scattering (SERRS) spectra, were measured using silver nanoparticle dimers containing dye molecules within their nanogaps. Our results show that the spectral peaks of FR for subradiant resonances appear near the dips in Rayleigh scattering spectra. Furthermore, these FR peaks exhibit blue-shifts during the quenching processes of both ultrafast SEF and SERRS. We examined these static and temporal spectral properties using a coupled oscillator model composed of radiant plasmons, subradiant plasmons, and molecular excitons. The static properties were reproduced by increasing the linewidths of the radiant plasmon resonance, while the temporal properties were captured by decreasing the EM coupling energies between the exciton and both plasmon oscillators. These findings indicate that this methodology is a powerful tool for evaluating EM coupling between subradiant plasmons and molecular excitons.
Hydroxyapatite (HA) on a magnesium (Mg) surface is studied using density functional theory, to help understand the effect of HA coating and alloying in the surfaces of Mg-based biodegradable implants. We determine the adsorption energies and structural changes of a single layer of HA on pure Mg(0001) and on sparsely calcium (Ca) or zinc (Zn) doped Mg(0001) and find that both Zn and Ca doping improves the adsorption, except in a few positions of HA relative to the dopant position. All adsorption configurations, whether with pure or doped Mg surfaces, show deformation of the surface and HA layer. For Ca doping, we found that for a certain adsorption configuration, the dopant Ca atom moves out of the Mg surface and into the HA layer, leaving behind a Mg vacancy in the top layer of the Mg surface. Plots of electron density changes show that electrons accumulate around the Ca dopant and the neighboring Mg atoms, while in Zn doping this is less pronounced. Overall, our results demonstrate that the dopant choice and relative position of HA influence the interaction between HA and Mg-surfaces, and affect both adsorption energies and atomic and electronic structures.
Electrostatic interactions fundamentally govern the structure, stability, and dynamics of charged (bio)matter, yet the impact of heterogeneous and anisotropic charge distributions on the behavior of protein solutions remains elusive. Here, we introduce a versatile multiscale framework that directly connects molecular-level electrostatics to collective properties via a colloid-inspired coarse-grained modeling combined with neural network-assisted optimization. Using monoclonal antibodies as model system, our inverse design approach identifies charge patterns capable of reliably reproducing experimental structure factors, osmotic compressibility and collective diffusion coefficients in a wide region of protein concentrations. Close inspection of our data further uncovers how specific physical features and spatial arrangements of localized charge patches significantly influence the solution structure. This transferable strategy provides a predictive pathway to decode and control charge-driven interactions in complex biomolecules and, more generally, in heterogeneously-charged soft matter systems, with immediate relevance to protein formulation and biomaterials engineering.