2025-07-15 | | Total: 16
Transition state (TS) structures define the critical geometries and energy barriers underlying chemical reactivity, yet their fleeting nature renders them experimentally elusive and drives the reliance on costly, high-throughput density functional theory (DFT) calculations. Here, we introduce TS-GEN, a conditional flow-matching generative model that maps samples from a simple Gaussian prior directly to transition-state saddle-point geometries in a single, deterministic pass. By embedding both reactant and product conformations as conditioning information, TS-GEN learns to transport latent noise to true TS structures via an optimal-transport path, effectively replacing the iterative optimization common in nudged-elastic band or string-method algorithms. TS-GEN delivers unprecedented accuracy, achieving a root-mean-square deviation of 0.004 ˚A (vs. 0.103 ˚A for prior state-of-the-art) and a mean barrier-height error of 1.019 kcal/mol (vs. 2.864 kcal/mol), while requiring only 0.06 s GPU time per inference. Over 87% of generated TSs meet chemical-accuracy criteria (<1.58 kcal/mol error), substantially outpacing existing methods. TS-GEN also exhibits strong transferability to out-of-distribution reactions from a larger database. By uniting sub-angstrom precision, sub-second speed, and broad applicability, TS-GEN will be highly useful for high-throughput exploration of complex reaction networks, paving the way to the exploration of novel chemical reaction mechanisms.
Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na2SO4(aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in 0.1-2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium-sulfate ion pairing at the low concentration limit (0.1 M) which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. We demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns long simulations.
We present a detailed derivation and discussion of cavity Born-Oppenheimer coupled cluster (CBO-CC) theory and address cavity-modified electron correlation in the vibrational strong coupling regime. Methodologically, we combine the recently proposed cavity reaction potential (CRP) approach with the Lagrangian formulation of CC theory and derive a self-consistent CRP-CC method at the singles and doubles excitations level (CRP-CCSD). The CRP-CC approach is formally similar to implicit solvation CC models and provides access to the CBO-CC electronic ground state energy minimized in cavity coordinate space on a CC level of theory. A hierarchy of linearisation schemes (lCRP-CCSD) similar to canonical CC theory systematically lifts the self-consistent nature of the CRP-CCSD approach and mitigates numerical cost by approximating electron correlation effects in energy minimization. We provide a thorough comparison of CRP-CCSD, lCRP-CCSD and CRP-Hartee-Fock methods for a cavity-modified Menshutkin reaction, pyridine+CH3Br, and cavity-induced collective electronic effects in microsolvation energies of selected methanol-water clusters. We find lCRP-CCSD methods to provide excellent results compared to the self-consistent CRP-CCSD approach in the few-molecule limit. We furthermore observe significant differences between mean-field and correlated results in both reactive and collective scenarios. Our work emphasizes the non-trivial character of electron correlation under vibrational strong coupling and provides a starting point for further developments in ab initio vibro-polaritonic chemistry beyond the mean-field approximation.
Wearable sensors are critical components in smart sports systems for real-time monitoring of athletic performance, physiological conditions, and psychological states. In this work, a sodium alginate/gelatin-based triboelectric nanogenerator (SG-TENG) was developed for mechanical energy harvesting and real-time sensing in track and field applications. The SA/gelatin composite film exhibits high transparency, good flexibility, and uniform morphology, enabling stable triboelectric output. The SG-TENG delivers high electrical performance, including an open-circuit voltage (Voc) of 156.6 V, short-circuit current (Isc) of 46.9 uA, and transferred charge (Qsc) of 139.6 nC, with a maximum power of 13.5 mW under optimal loading. Its output characteristics are closely related to mechanical parameters such as frequency, force, displacement, and contact area. Furthermore, the device demonstrates energy storage capability by charging capacitors under various dynamic conditions. Integrated into a running shoe, the SG-TENG enables self-powered gait monitoring and can distinguish between walking, running, and jumping. It also shows potential for inferring psychological and physiological states from gait features, offering a promising solution for battery-free, multifunctional sensing in sports performance evaluation and health monitoring.
We present a domain-based local pair natural orbital Møller--Plesset second order perturbation theory (DLPNO-MP2) for periodic systems, working within an LCAO formalism within the Tubromole program package. This approach, Megacell-DLPNO-MP2, embeds a supercell correlation treatment within a megacell and does not involve periodic image summation for the Coulomb integrals. Working in a basis of well-localised Wannier functions, periodicity is instead imposed through rigorous translational symmetry of Hamiltonian integrals and wavefunction parameters. The accuracy of the method is validated through comparison with a complementary periodic DLPNO-MP2 method that employs Born--von Kármán boundary conditions, described in paper I of this series. The PNO approximations are shown to be equivalent in the two approaches and entirely consistent with molecular DLPNO-MP2 calculations. The Megacell-DLPNO-MP2 method displays sub-linear scaling with respect to supercell size at the asymptotic limit and example calculations are presented with up to 15000 basis functions in the correlation treatment.
We present domain-based local pair natural orbital Møller--Plesset second order perturbation theory (DLPNO-MP2) with Born--von Kármán boundary (BvK) conditions. The approach is based on well-localised Wannier functions in a LCAO formalism and extends the molecular DLPNO-MP2 implementation Tubromole program package to periodic systems. The PNOs are formed through a PAO-OSV-PNO cascade, using BvK projected atomic orbitals and orbital specific virtuals as intermediaries in an analogous manner to the molecular scheme. Our chargeless and surface-dipole corrected local density fitting approach is shown to be numerically stable and to ensure convergent lattice summations over the periodic images for the two- and three-index Coulomb integrals. Through careful benchmarking, we show that the DLPNO approximations in the BvK-DLPNO-MP2 methods are entirely consistent with those of molecular DLPNO-MP2 calculations, and with an alternative periodic approach Megacell-DLPNO-MP2, reported in Paper II of this series. Smooth convergence to the canonical correlation energy with tightening PNO threshold is observed. Reference MP2 correlation energies are provided for a set of 2D and 3D periodic systems using a triple-zeta basis and supercell sizes up to 11×11 and 7×7×7.
The gas-phase hydrogen exchange reaction (HER) is the most fundamental chemical process for benchmarking quantum reaction dynamics. In this Letter, we focus on controlling HER by means of strong light-matter coupling inside a resonant cavity, an approach often called polariton chemistry. In particular, we focus on the isotopic variation of HER involving collisions between molecular hydrogen H2 and deuterium atom D, i.e., H2+D→HD+H. We find that the asymmetry introduced by the different isotopes, despite being small, enables strong cavity-induced modifications of reaction rates. Outside of the cavity the reaction is as usual D+H2→DH+H. However, inside the cavity another type of reactions take place where D+H2→DH+H+Ephoton, where Ephoton=ℏωcav. Our results show that HER is an ideal platform to make a significant step toward closing the gap between theory and experiment in polariton chemistry.
Currently, the pre-treatment for membrane filtration systems is receiving more attention that can be remove large particles and debris or enhance permeation flux. In this study, commercial nylon 6 fabric is modified via hydrolysis amidation reaction in alkaline media with NaOH and Dimethylene diamine. The increase of hydrophilicity amino groups on the surface of treated fabric which increases the hydrophilicity of the fabric. This opens up new possibilities for the treatment process in commercial membrane filtration systems. Besides, amino enriched nylon 6 fabric can capture carboxylic functionalized carbon-based nanomaterials forwards to the synthesis of thin film nanocomposite membranes based on carboxylated multi-walled carbon nanotubes.
Nuclear Magnetic Resonance (NMR) spectroscopy is a central characterization method for molecular structure elucidation, yet interpreting NMR spectra to deduce molecular structures remains challenging due to the complexity of spectral data and the vastness of the chemical space. In this work, we introduce DiffNMR, a novel end-to-end framework that leverages a conditional discrete diffusion model for de novo molecular structure elucidation from NMR spectra. DiffNMR refines molecular graphs iteratively through a diffusion-based generative process, ensuring global consistency and mitigating error accumulation inherent in autoregressive methods. The framework integrates a two-stage pretraining strategy that aligns spectral and molecular representations via diffusion autoencoder (Diff-AE) and contrastive learning, the incorporation of retrieval initialization and similarity filtering during inference, and a specialized NMR encoder with radial basis function (RBF) encoding for chemical shifts, preserving continuity and chemical correlation. Experimental results demonstrate that DiffNMR achieves competitive performance for NMR-based structure elucidation, offering an efficient and robust solution for automated molecular analysis.
We present assembly-theory, a Rust package for computing assembly indices of covalently bonded molecular structures. This is a key complexity measure of assembly theory, a recent theoretical framework quantifying selection across diverse systems, most importantly chemistry. assembly-theory is designed for researchers and practitioners alike, providing (i) extensible, high-performance implementations of assembly index calculation algorithms, (ii) comprehensive benchmarks against which current and future algorithmic improvements can be tested, and (iii) Python bindings and RDKit-compatible data loaders to support integration with existing computational pipelines.
We present a machine learning model for predicting the electric dipole moment of diatomic molecules using only the atomic properties of the constituent atoms. Our model can screen the entire periodic table and identify the molecules with the largest dipole moment for applications in cold molecular sciences, or find the molecule with the largest dipole moment that contains a given atom. Similarly, our model identifies useful trends that explain the dipole moment of molecules, improving our intuition in chemical physics. Finally, we condense our model into an analytical expression to predict the dipole moment in terms of atomic properties.
An optical phonon at 2.5 meV is observed in a thin film of the chiral metal halide (R-MBA)2PbI4, but is absent in the racemic counterpart, as revealed by femtosecond transient absorption spectroscopy. This experimental result indicates the chiral origin of the 2.5 meV mode and supports recent theoretical predictions of chirality transfer from the organic to the inorganic layers, with implications for the spin polarization properties of hybrid metal halides and perovskites.
Fast and accurate generation of molecular conformers is desired for downstream computational chemistry and drug discovery tasks. Currently, training and sampling state-of-the-art diffusion or flow-based models for conformer generation require significant computational resources. In this work, we build upon flow-matching and propose two mechanisms for accelerating training and inference of generative models for 3D molecular conformer generation. For fast training, we introduce the SO(3)-Averaged Flow training objective, which leads to faster convergence to better generation quality compared to conditional optimal transport flow or Kabsch-aligned flow. We demonstrate that models trained using SO(3)-Averaged Flow can reach state-of-the-art conformer generation quality. For fast inference, we show that the reflow and distillation methods of flow-based models enable few-steps or even one-step molecular conformer generation with high quality. The training techniques proposed in this work show a path towards highly efficient molecular conformer generation with flow-based models.
The excitation energy of 229mTh is extremely low at 8.4eV; thus, this isotope exhibits changes in its decay modes depending on the chemical state, specifically the outermost electronic states. However, the reported half-lives of γ-ray transition are not consistent among the previous experiments. In this study, we investigate the chemical states of 229mTh by DFT calculations. Based on these results, we evaluate the relationship between the experimental half-life for each sample and the electronic state of Th. The calculation results indicate that ion trap method, CaF2 model and MgF2 one decay only via the γ-ray transition, whereas LiSrAlF6 decay via the γ-ray transition and has a possibility of decay via internal conversion, and electron bridge.
Proton transfer reactions are fundamental to many chemical and biological systems, where quantum effects such as tunneling, delocalization, and zero-point motion play key kinetic control roles. However, classical methods capable of accurately capturing these phenomena scale prohibitively with system size. Here, we develop and demonstrate quantum computing algorithms based on the Nuclear-Electronic Orbital framework, treating the transferring proton quantum mechanically. We assess the potential of current quantum devices for simulating proton transfer kinetics with high accuracy. We first construct a deep initial ansätze within a truncated orbital space by employing the frozen natural orbital approximation. Then, to balance circuit depth against state fidelity, we implement an adaptive form of approximate quantum compiling. Using resulting circuits at varying compression levels transpiled for the ibm_fez device, we compute barrier heights and delocalised proton densities along the proton transfer pathway using a realistic hardware noise model. We find that, although current quantum hardware introduces significant noise relative to the demanding energy tolerances involved, our approach allows substantial circuit simplification while maintaining energy barrier estimates within 13% of the reference value. Despite present hardware limitations, these results offer a practical means of approximating key circuit segments in near-term devices and early fault-tolerant quantum computing systems.
Protein-ligand binding affinity prediction is essential for drug discovery and toxicity assessment. While machine learning (ML) promises fast and accurate predictions, its progress is constrained by the availability of reliable data. In contrast, physics-based methods such as absolute binding free energy perturbation (AB-FEP) deliver high accuracy but are computationally prohibitive for high-throughput applications. To bridge this gap, we introduce ToxBench, the first large-scale AB-FEP dataset designed for ML development and focused on a single pharmaceutically critical target, Human Estrogen Receptor Alpha (ERα). ToxBench contains 8,770 ERα-ligand complex structures with binding free energies computed via AB-FEP with a subset validated against experimental affinities at 1.75 kcal/mol RMSE, along with non-overlapping ligand splits to assess model generalizability. Using ToxBench, we further benchmark state-of-the-art ML methods, and notably, our proposed DualBind model, which employs a dual-loss framework to effectively learn the binding energy function. The benchmark results demonstrate the superior performance of DualBind and the potential of ML to approximate AB-FEP at a fraction of the computational cost.