Peer Revenue Publish or profit

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Physics Chemical Physics

physics.chem-ph Watch Chemical Physics (sign in to use)

Wednesday 9 September

Chemical Physics

Efficient Deterministic-Stochastic Representation of the Coulomb Operator in Real Space

We present an efficient mixed deterministic-stochastic approach for constructing the Coulomb operator in real space that preserves accuracy across the full spectral range. The dominant long-range components of the interaction are captured deterministically via a compact, low-rank approximation, constructed using Chebyshev-filtered subspace iteration, while the remaining spectral tail is treated using unbiased probing with a small number of stochastic vectors. The resulting operator is benchmarked within self-consistent field calculations for the extended Hubbard Hamiltonian on periodic, perturbed, and non-periodic three-dimensional lattices. We account for stochastic bias in observables by implementing a jackknife correction. The method systematically improves with deterministic rank and stochastic sample count, while substantially reducing the computational cost of constructing the Coulomb matrix.
Read the full abstract

All authors Nafiz Faiaz, Tucker Allen, Dimitri Bazile, Kajsa Williams, Daniel Neuhauser

Chemical Physics

Resolving the ground state intersection problem in coupled cluster theory

A physically correct description of ground state conical intersections is essential for simulating non-adiabatic dynamics of radiationless decay to the electronic ground state. Coupled cluster theory offers a particularly balanced treatment of static and dynamic correlation in excited states. However, coupled cluster methods and other single-reference approaches fail to describe ground state intersection regions, and several workarounds have been proposed to address this shortcoming. Here we present convex similarity constrained coupled cluster theory (CVX-SCC), a framework that resolves the ground state intersection problem by construction. We have applied this theory to two different coupled cluster methods, and tested the resulting models on ethylene, uracil, PSB3, and HeH$_2$, showing that both variants produce physically correct ground state conical intersections. With future development of nuclear gradients and non-adiabatic coupling vectors, this approach will provide single-reference electronic structures suitable for non-adiabatic dynamics simulations all the way to the ground state.
Read the full abstract

Chemical Physics

Modeling electrolytes in nanopores by Monte Carlo simulations and the Bazant–Storey–Kornyshev model

We study electrolyte-filled cylindrical nanopores through Monte Carlo (MC) simulations and the Bazant–Storey–Kornyshev model, for different ionic valencies, sizes, and bulk concentrations, pore radii and surface charge densities. Our model accounts for finite ion size through a Stern layer and we use de~Souza and Bazant's mechanical equilibrium principle to derive a boundary condition for the outer Helmholtz plane. For 1:1 electrolytes and moderate surface charge densities, we find that both the classical Poisson–Boltzmann–Stern (PB–Stern) and our Bazant–Storey–Kornyshev–Boltzmann–Stern (BSKB–Stern) model fit MC data well—for 2:1 and 3:1 electrolytes, the BSKB–Stern outperforms the PB-Stern model. Conversely, the BSKB–Stern model poorly fits MC data in other scenarios. First, BSKB–Stern does not capture drying and apparent like-charge attraction in 3:1 electrolytes and small surface charge densities. Second, BSKB–Stern only fits MC data for finely-tuned ionic diameters; for other ionic diameters, the ionic charge densities show oscillations or extended near-surface regions caused by ionic packing and strong Coulomb interactions, not captured by the BSKB–Stern model.
Read the full abstract

All authors Nader Nekoubin, David Fertig, Dezső Boda, Mathijs Janssen

Chemical Physics

Analytic Gradients and Nonadiabatic Couplings for Device-Resident DMRG-QD-NEVPT2 Through Conical Intersections on a Consumer GPU

Nonadiabatic dynamics through a conical intersection requires both static and dynamic correlation and, at every geometry, an excited-state gradient and interstate nonadiabatic coupling (NACME). Analytic multireference derivatives at this level have been confined to datacenter-scale hardware. We report device-resident \mbox{DMRG-QD-NEVPT2} with analytic gradients and NACMEs within the $8$\,GB of a consumer GPU: a density-matrix-renormalization-group (DMRG) active-space reference supplies the static correlation, and quasi-degenerate SC-NEVPT2 dresses it through a multi-state effective Hamiltonian that remains valid through a conical intersection, where single-state perturbation theory develops a vanishing denominator and fails. Every gradient and NACME is one reverse-mode transpose of a single contraction graph: the DMRG sweep enters as a gauge-free differentiable node, and all two-electron work is AO-direct through Cholesky factors. The DMRG reference is exact within the active space ($\sim\!10^{-15}$\,Ha). The corrected adiabats close smoothly to a real $S_0/S_1$ conical intersection with a non-vanishing, $1/\Delta E$-divergent NACME—where linear-response TDDFT returns zero identically—and the dressing restores the differential dynamic correlation that ionic $\pi\pi^*$ states require. DMRG-NEVPT2 and quasi-degenerate NEVPT2 are each established; what is new is the combined stack—multi-state and through a conical intersection, with analytic gradients and couplings—realized device-resident on commodity hardware.
Read the full abstract

Chemical Physics

More is not always better: Dissociative photoionization limits the EUV absorbing photoacid generator pentafluorophenyl triflate in photolithography

Pentafluorophenyl triflate has been explored as a highly absorbing neutral photoacid generator (PAG) candidate for next generation chemically amplified resists used in extreme ultraviolet (EUV) lithography. Although increased fluorination enhances EUV absorption, this study demonstrates that such an approach does not necessarily improve photoacid generation efficiency. Using photoelectron-photoion coincidence (PEPICO) spectroscopy at the 92 eV photon energy of the EUV scanners in combination with quantum chemical calculations, the dissociative photoionization of pentafluorophenyl triflate was systematically investigated. The photoionization mass spectrum reveals extensive fragmentation, with the parent ion contributing only 3.1 % of the total signal and CF$_3^+$ representing the dominant product ion. Computed appearance energies align well with experimental trends and support a sequential fragmentation pathway involving loss of SO$_2$, CF$_3$, and CO. Crucially, none of the major dissociation channels yield precursors capable of forming triflic acid, the strong photoacid required for efficient deprotection reactions in chemically amplified resists. Combined with previous dissociative electron attachment studies indicating similarly unfavorable fragmentation, the results demonstrate that despite its high EUV absorption cross section, pentafluorophenyl triflate is unsuitable as a PAG for EUV lithography. The findings highlight the importance of understanding fundamental photoionization and electron interaction mechanisms to guide the rational design of next generation high performance EUV photoresists.
Read the full abstract

All authors Fabian Holzmeier, Michiel J. van Setten, Ragnar Björnsson, Robert Richter, Filipe Ferreira da Silva, Oddur Ingólfsson

Chemical Physics

Comprehensive Study of L-Menthol and Octanoic Acid as a Hydrophobic Eutectic Solvent

Hydrophobic eutectic solvents (HES) based on natural compounds represent promising green alternatives to conventional solvents. In this work, we investigate the physicochemical, structural, and dynamical properties of an ES formed by L-menthol and octanoic acid using a combined experimental and molecular dynamics simulation approach. Five compositions with molar ratios from 1:3 to 3:1 were studied with molecular dynamics simulation in the temperature range 15 degrees C to 35 degrees C. Experimental measurements of density and viscosity in the temperature range from 5 degrees C to 35 degrees C were complemented with results obtained from MD simulations employing the OPLS force field. Structural analyses based on radial distribution functions and Kirkwood-Buff integrals reveal that the dominant interactions in the mixture are hydrogen bonds between L-menthol and octanoic acid molecules. Dynamic properties, including self-diffusion coefficients and hydrogen-bond lifetimes, indicate that intermolecular hydrogen bonds between the two components are stronger and longer-lived than bonds between identical species. These findings provide molecular-level insight into the structure and transport properties of menthol-based ESs relevant for green solvent applications.
Read the full abstract

All authors Bernarda Lovrinčević, Martina Požar, Romana Cerc Korošec, Klara Švegelj, Peter Ogrin, Marija Bešter-Rogač, Tomaž Urbič

Tuesday 8 September

Chemical Physics · Cross-listed from cond-mat.stat-mech

The Role of Concentration in Determining NaCl Nucleation Mechanism: A Story of Pathways Coexistence

Deciphering the microscopic details of crystal nucleation remains a key open problem in chemical physics. For aqueous sodium chloride (NaCl), the extent to which nucleation involves amorphous ion aggregation, and the supersaturation regime in which two-step-like behavior emerges, remains actively debated. In this work, we examine NaCl nucleation at 12, 13, and 14 mol/kg ($S = 3.2, 3.5, 3.8$) using unbiased path sampling with $\infty$RETIS. Through a combination of reaction coordinate analysis, joint and conditional probability distributions, diffusion tensor analysis, and conditional free energy landscapes, we identify nucleation pathways without a priori mechanistic assumptions. While all concentrations are governed by the same reaction coordinate, indicating that nucleus growth and structural ordering are central to crossing the barrier, the nature of reactive trajectories changes with supersaturation. At lower concentrations, nucleation mainly occurs via tightly coupled increases in cluster size and ordering, compatible with a one-step-like mechanism. As concentration increases, amorphous aggregation becomes increasingly decoupled from crystallization, and pathways involving substantial amorphous growth prior to crystallization become progressively more probable. Importantly, these pathways coexist within a single broad reaction channel rather than proceeding through distinct metastable intermediate states. Together, these results reconcile earlier, contradictory reports on NaCl nucleation mechanisms and support a view of nucleation as an ensemble of competing reactive pathways whose relative probabilities vary continuously with supersaturation. More broadly, our findings illustrate that substantial mechanistic diversity can emerge within a classical nucleation framework, without requiring distinct non-classical descriptions.
Read the full abstract

Chemical Physics

Quantum-accurate atomistic modeling of enzyme catalysis using a machine learned potential

Electronic rearrangements associated with bond forming/breaking in catalytic enzymes require quantum mechanical (QM) treatment beyond classical molecular mechanics (MM). Hybrid QM/MM methods enable tractable simulations but require system-specific setup and are sensitive to the QM region choice and treatment of the QM/MM interface. We demonstrate quantum-accurate treatment of all-atom, complete enzymes in explicit solvent comprising up to 54k atoms and 1 microsecond of total simulation time using the machine-learned interatomic potential (MLIP) eSEN-omol. We reproduce experimental barrier trends for Claisen rearrangement in chorismate mutase, resolve critical intermediate states in PETase catalyzed polymer depolymerization, and distinguish mechanistic alternatives for metal-activated phosphoryl transfer in nucleoside diphosphate kinase. We realize 1000x speedups relative to typical QM/MM calculations without system-specific tuning. These results establish MLIPs as a practical route to QM-accurate simulations of enzyme catalysis.
Read the full abstract

All authors Meng Gao, Armin Shayesteh Zadeh, Aniruddha Seal, Siva Dasetty, Siddarth K. Achar, Misko Dzamba, Benjamin K. Miller, Leif D. Jacobson, C. Lawrence Zitnick, Brandon M. Wood, Zachary W. Ulissi, Daniel S. Levine, Andrew L. Ferguson

Chemical Physics

Multipole splats for optimized and inverted effective potentials

In modern density functional theory, effective Hamiltonians are constructed to reproduce densities and forces of electrons at equilibrium. However, their nonlocal potentials leave systematic errors in spectral properties and real-time dynamics. Access to local optimized or inverted effective potentials would remove this limitation, but the numerical fragility in finite orbital bases has long prevented their wide adoption. Here, we introduce multipole splats, a class of trial potentials that carry the correct asymptotic decay required to support the unoccupied spectrum. By connecting the computation of effective potentials to variational and supervised variants of Hamiltonian learning, we recast both problems as stable nonlinear optimization formulated directly in standard orbital basis sets and applicable to any hybrid functional approximation. The resulting solver allows us to resolve spatial profiles of exchange-correlation potential errors during molecular dissociation and accurately reconstruct key excited states without empirical asymptotic corrections. We also show the deviation from the ionization potential theorem for different exchange-correlation approximations on a dataset of molecular systems. Multipole splats extract insights from established approximations and provide capacity for robust dataset generation for downstream processing and learning.
Read the full abstract

Chemical Physics

Balanced and broadly-normed meta-generalized gradient approximation

Conventional semi-local density functional approximations for the exchange-correlation energy systematically underestimate fundamental band gaps. Advanced meta-generalized gradient approximations (meta-GGAs) can partially overcome this limitation through their orbital dependence via the kinetic energy density. In principle, the flexible meta-GGA form should allow for balanced performance across diverse properties. However, atomization energies, reaction barrier heights, and lattice constants, as well as band gaps, exhibit different sensitivities to the large parameter space of meta-GGA functionals. In this work, we explore the space of appropriate norms that guide the r$^2$SCAN-type functional construction through the interpolation function connecting two fundamental paradigms: the one-electron limit and the uniform electron gas limit. We introduce local modifications to a strongly-smoothed interpolation function and investigate their impact on functional performance. To further probe the limits of the meta-GGA framework, we include spin-unpolarized bonded systems among the set of norms. The resulting bonded-norm bn-r$^2$SCAN meta-GGA achieves balance: band-gap predictions approaching the accuracy of the ultra-nonlocal LAK meta-GGA and lattice constants of nearly r$^2$SCAN accuracy.
Read the full abstract

All authors Yanyong Wang, Chandra Shahi, John P. Perdew, Adrienn Ruzsinszky

Chemical Physics · Cross-listed from math-ph

The inverse problem of time-dependent density functional theory on the torus

We prove that, for any given sufficiently regular one-particle density, there exists a unique external potential for which the solution to the many-body Schrödinger equation has this prescribed density. This implies in particular that one can exactly reproduce the time-dependent density of an interacting system using non-interacting electrons, which is at the heart of time-dependent density functional theory. Our result covers the Coulomb interaction and densities arising from extended nuclei.
Read the full abstract

Chemical Physics · Cross-listed from cond-mat.mtrl-sci

Time-Dependent Density Functional Theory with Coulomb Interactions

We provide the first fully rigorous justification of time-dependent density functional theory in the continuum: Given a time-dependent density, we prove that there exists an external potential, unique up to a time-dependent constant, such that the corresponding Schrödinger equation reproduces the given density. This potential can be obtained with an iteration scheme. Our argument covers Coulomb interactions and does not need any Taylor expansion in time. It instead requires analyticity in space for the density and potential, which is compatible with extended nuclei.
Read the full abstract

Chemical Physics

Orbital choice in constructing model Hamiltonians

Model Hamiltonians (e.g., the Hubbard or Heisenberg models) provide a simple yet physically meaningful description of complex electronic phenomena through a small number of parameters, such as electron hopping integrals or magnetic exchange couplings. Their construction typically exploits the local nature of electron correlation effects and therefore relies on localized orbitals. Although many localization schemes are available, the resulting localized orbitals depend not only on the chosen localization functional but also on the orbital space to which the localization procedure is applied. These choices are often treated as technical details and are rarely discussed explicitly. Here, we demonstrate that different, chemically reasonable choices of localized orbitals can lead to substantially different model Hamiltonian parameters, even when they produce nearly identical electronic energies. Using Density Matrix Downfolding, we systematically investigate how the localization functional and, more importantly, the choice of model orbital space affect effective Hamiltonians derived from ab initio calculations for representative $\pi$-conjugated systems and transition metal complexes.
Read the full abstract

Chemical Physics · Cross-listed from physics.comp-ph

Structure-Informed Bayesian Inference of Anomalous Transport and Hidden Molecular Trapping in Amorphous Media

Molecular diffusion in fluctuating amorphous and macromolecular media governs key transport processes across soft-matter physics, energy storage, and biological membranes. Extracting localized trapping states from single-particle tracking trajectories remains a fundamental challenge; because thermal structural breathing continuously reconfigures pore boundaries, conventional geometric algorithms suffer from severe systematic biases, erroneously merging distinct localized states during cyclic molecular returns. Here, we address this deadlock by shifting the paradigm from local geometric recurrence to a structure-informed Bayesian regularization. Leveraging discrete Morse theory, we extract the time-invariant topological skeleton of the fluctuating host matrix to construct robust, gas-specific physical priors that account for individual molecular dimensions. Trajectory steps are sequentially partitioned via a two-stage probabilistic refinement that dynamically adapts to the transport landscape. Benchmarked against a rigorous environment where synthetic particles explore the actual interconnected matrix graph, our approach eliminates systemic biases, restricting macroscopic trapping parameter deviations to just a few percent under optimal linear $O(N)$ computational scaling. Applied to hydrogen and methane transport within a type-I kerogen matrix, serving as a prototype for highly tortuous, flexible macromolecular networks, the method successfully decodes the hidden microscopic mechanisms of confined diffusion. To ensure immediate broad impact, the documented open-source code and data are made publicly available, offering an accessible strategy readily adaptable to a broad spectrum of tracking phenomena, from ion transport in battery polymers to protein trafficking within cellular environments.
Read the full abstract

All authors Andrey Ananev, Maria Potapova, Nikolay Kondratyuk, Timur Vostroknutov, Aleksey Khlyupin

Chemical Physics

Examination of Cluster Based Methods to Fermionic Systems

We present the cluster mean field(cMF) and its extensions into perturbation theory(cPT) and coupled cluster(cCC) methods for fermionic systems. We discuss how the commutation rules for fermionic systems modifies the method as well as other computational details. We then perform calculations on the Hubbard model and a few select conjugated molecules to demonstrate the efficacy of the method. We will also use these examples to compare cluster methods to CAS-SCF.
Read the full abstract

Chemical Physics · Cross-listed from cond-mat.str-el

Examining Convergence of Cluster Perturbation Theory and Cluster Coupled Cluster for J1-J2 Heisenberg Spin Lattices

In previous work, it was shown that using the cluster mean field(cMF) wavefunction as the zeroth order wavefunction for perturbation theory and coupled cluster produces a reasonable approximation for the J 1-J 2 Heisenberg model in the thermodynamic limit(TDL). However, this work was limited to cPT4 and cCCSD. The central claim of cMF is that it takes strongly correlated systems and makes them weakly correlated so that PT weakly correlated methods can be used. To confirm that this claim is true, it is important to show that the PT and CC series converge. In this paper, we extend the previous calculations to cPT7 and cCCSDTQ5 where we show that the energy is converging. We also discuss how to perform cluster based calculations efficiently in the thermodynamic limit, which is crucial for these more expensive calculations.
Read the full abstract

Chemical Physics

Competing Ring-Opening and Hofmann Elimination Pathways in Aqueous TEMPO Catholytes: A First-Principles Study

Aqueous redox-flow batteries based on TEMPO derivatives are promising for large-scale energy storage, but their practical use is limited by the chemical instability of the oxidized N -oxoammonium state. In this work, we investigate the degradation of five TEMPO derivatives using ab initio molecular dynamics combined with enhanced sampling. Two proposed degradation mechanisms, ring opening and Hofmann elimination, are examined and their corresponding activation free energies are compared. For all derivatives considered, ring opening exhibits a lower activation free energy than Hofmann elimination, identifying it as the kinetically preferred degradation pathway. The magnitude of the ring-opening barrier, however, varies significantly between molecules, showing that different functionalizations strongly influence its stability toward degradation. The predicted preference for ring opening is consistent with available experimental studies, which have identified or inferred ring-opening degradation for several TEMPO-based catholytes. These results provide an atomistic picture of degradation pathways that are difficult to resolve experimentally and highlight the importance of molecular structure in controlling the kinetic stability of TEMPO derivatives in aqueous electrolytes.
Read the full abstract

All authors Line Mouaffac, Maiwenn Souetre, Guillaume Jeanmairet, Mathieu Salanne

Chemical Physics

TSBench: A physics-grounded benchmark for evaluating LLM understanding of chemical reaction mechanisms

Understanding a chemical reaction requires mapping a symbolic reactant-product description onto the three-dimensional pathway through which atoms rearrange, yet chemistry benchmarks for large language models (LLMs) largely probe factual knowledge and text-based reasoning. Here we introduce TSBench, a benchmark in which an LLM agent uses structure-editing tools to construct three-dimensional transition-state (TS) guesses verified by an automated quantum-chemical pipeline, yielding a physics-grounded pass/fail verdict. Across 546 evaluations of seven frontier LLMs on 78 elementary reactions, the aggregate success rate rose from 50.4% to 66.8% under diagnosis-driven revision; the best models approached 90% on the simplest reactions, yet performance dropped sharply with mechanistic complexity. Most failed attempts produced locally plausible saddle points whose reaction paths led to the wrong reactant-product pair, revealing local geometric intuition without a robust grasp of the global reaction coordinate. TSBench establishes a mechanism-level yardstick for LLM agents in mechanism-sensitive tasks such as synthesis planning and autonomous experimentation.
Read the full abstract

Chemical Physics · Cross-listed from quant-ph

Explicit block encodings of rate matrices for simulating polymerization kinetics on quantum computers

Predicting how molecular weight distribution and monomer sequence evolve during polymerization is central to polymer science, yet classical approaches face a trade-off between molecular resolution and computational cost: for copolymers, the number of distinguishable species grows exponentially with chain length. Quantum computing offers a potential alternative, provided the non-unitary rate matrices governing the kinetics can be embedded into unitary quantum circuits, a task known as block encoding. Here we construct explicit block-encoding circuits for two kinetic models of living polymerization: Model A, single-monomer polymerization, whose lower-bidiagonal rate matrix is encoded via a sparse-oracle construction and a two-term linear combination of unitaries (LCU) decomposition; and Model B, two-monomer copolymerization, where a bijective labeling of polymer species by an integer index (the m-index) yields a structured sparse matrix encoded via either a five-term LCU or a sparse-oracle construction. Numerical simulations with the sparse-oracle encodings reproduce the classical time evolution for reactivity ratios drawn from reported olefin copolymerization systems spanning near-random ($r_1 r_2 \simeq 1$) and blocky ($r_1 r_2 > 1$) microstructures, and the LCU encodings are verified by explicit reconstruction of the encoded matrix block. Resource estimation shows that both implementations require only $O(\log N)$ qubits in the matrix dimension $N$ (an exponential memory saving over the classical state space), with gate counts growing gradually, reaching $10^4$ to $10^5$ gates at $10^3$ system qubits. These results establish a concrete quantum circuit foundation for simulating polymerization kinetics on fault-tolerant quantum hardware, and a first step toward exploiting exponential state-space compression for high-dimensional polymer reaction networks.
Read the full abstract

All authors Yuhei Ikeda, Hokuto Iwakiri, Soichiro Nishio, Kentaro Matsumoto

Chemical Physics · Cross-listed from cond-mat.mes-hall

From London to Morse via Binnig, Quate, and Gerber

In their landmark paper introducing the atomic force microscope [Phys. Rev. Lett. 56, 930 (1986)], Binnig, Quate, and Gerber presciently anticipated that the technique would ultimately be capable of probing interactions running the gamut from weak van der Waals interactions to strong covalent bonding. They also highlighted that the tip-sample forces central to AFM are present, and often highly influential, in scanning tunnelling microscopy; indeed, this realisation directly inspired the invention of the force microscope. In this perspective for the Forty Years of AFM special issue, we review selected aspects of two decades of work from our group at the University of Nottingham that span the force range highlighted by BQG and are united by a common, central theme: the probe as active participant rather than passive observer. Our selection of results also covers length- and correlation-scales from the microscopic right down to the single chemical bond limit, tracking a spectrum of interactions from van der Waals/Hamaker forces, through hydrogen bonding, to covalent bonds and, finally, atom-by-atom assembly of metal clusters via vertical tip-sample transfer. Echoing BQG's own observations on the prevalence of probe-sample forces in STM, we also discuss recent evidence that tip-induced heterogeneity underpins first-passage dynamics in molecular diffusion and highlight the challenges in acquiring non-invasive measurements of diffusion barriers for adsorbed molecules that are readily perturbed by the probe. We close with a perspective on machine learning's growing role in automating tip-driven atomic and molecular manipulation.
Read the full abstract

All authors Sofia Alonso Perez, Matthew O. Blunt, Frederick Carlisle, Neil R. Champness, Janette L. Dunn, Matthew Edmondson, Rowan Evers, Connor Fields, Subhashis Gangopadhyay, James Hayton, Samuel P. Jarvis, Filipe Junqueira, Lev Kantorovich, Brian Kiraly, Natalio Krasnogor, Ioannis Lekkas, Morten Møller, Philip Moriarty, Chris Pakes, Emmanuelle Pauliac-Vaujour, Oliver Phillips, Adrian Radocea, Philipp Rahe, Mohammad Abdur Rashid, Hongqian Sang, Alex Saywell, Nikhil Seeja Sivakumar, Peter Sharp, Andrew Stannard, Julian Stirling, Adam Sweetman, Simon Taylor, Richard A.J. Woolley

Chemical Physics · Cross-listed from cs.LG

MLIP Detective: Active Failure Mode Discovery Beyond Benchmark Scores for Machine-Learning Interatomic Potentials

Universal machine-learning interatomic potentials (u-MLIPs) aim to generalize across diverse configurations. Benchmarks enable reproducible evaluation but may not expose failures outside their predefined scope. Here, we show that physics-informed search can complement benchmark-based evaluation by uncovering hidden failure modes. We introduce MLIP Detective, an agentic framework for active failure mode discovery. Starting from benchmark evidence, MLIP Detective generates falsifiable, physics-informed failure hypotheses, screens them with inexpensive simulations, and escalates only the most suspicious cases to human experts together with proposed verification protocols. Without issue-specific prompting, MLIP Detective identified and characterized a systematic anomaly in MACE-MPA-0: the model predicted some relaxed adsorbate-surface systems involving O- or F-containing adsorbates to be higher in energy than their corresponding separated fragments. Using cross-model comparisons, MLIP Detective further inferred a likely training-data origin for the anomaly, consistent with recent reports.
Read the full abstract

All authors Ryuhei Okuno, Nontawat Charoenphakdee, Kaoru Hisama, Yuta Tsuboi

Chemical Physics

Fixed-Dimensional Latent Flow for Generating Variable-Size 3D Molecules

Molecular size is coupled to composition, structure, and function, yet most 3D molecular generators require a predefined atom count. We introduce Equivariant-Free Transformer-Autoencoded Latent Flow Matching, a two-stage framework that samples a fixed-dimensional latent vector using flow matching and uses an autoregressive Transformer to determine molecular size, atom types, coordinates, and chemical attributes. Canonical atom ordering and rigid-pose alignment enable Transformers without equivariant layers, while decoded attributes guide bond reconstruction. On PCQM4Mv2, unconditional generation yields 87.9% unique, novel molecules passing sanitization and PoseBusters checks, exceeding baselines with lower end-to-end training and sampling time and higher end-to-end throughput. Across ten target HOMO-LUMO gaps, internal ranking retains 30% of screened candidates and increases the density functional theory-verified hit rate within 0.1 eV from 25.0% to 52.4%, while largely preserving novelty and diversity. These results demonstrate fixed-dimensional latent generation with autoregressive decoding as a practical approach to molecular design without prespecifying size.
Read the full abstract

All authors Weichi Yao, Cameron Gruich, Bryan R. Goldsmith, Yixin Wang

Chemical Physics · Cross-listed from cond-mat.mtrl-sci

Band Structure Modulation of ZrO2 Nanoparticles for Control of CO Adsorption Properties: A Combined Density Functional Theory - Density Functional Tight Binding Study

We present a combined density functional theory (DFT) and density functional tight binding (DFTB) study of zirconia (ZrO2) nanoparticles of experimentally relevant sizes of several nanometers and their interactions with the CO molecule. A hybrid DFTB - Force Field (DFTB-FF) framework is developed, whereby band structure calculations rely on an existing Slater-Koster framework, while the accuracy of structural optimization and adsorption properties is controlled by the introduction of classical long-range interatomic potentials into DFTB instead of the traditional repulsive potentials. Additionally, coordination-dependent Zr-C potentials are introduced to account for the distinct local chemical environments of bulk-like facet sites and under-coordinated tip and edge sites, thereby improving the description of CO adsorption. This hybrid DFTB-FF approach substantially improves the robustness of geometry optimization and provides a practical strategy for extending the applicability of DFTB to complex oxide nanostructures. The calculations reveal that termination stoichiometry can be used to engineer intrinsic, p-type, or n-type electronic structures and thereby tune the adsorption activity of zirconia nanoparticles. While stoichiometric nanoparticles do not activate the C-O bond, low-coordinated sites in Zr-rich (n-type) nanoparticles exhibit chemisorption accompanied by charge donation into a CO antibonding LUMO-derived orbital, resulting in C-O bond activation. These results demonstrate that stoichiometry-controlled electronic structure and under-coordinated surface sites introduced by nanostructuring play a key role in governing the adsorption strength and reactivity of zirconia nanoparticles.
Read the full abstract

All authors Kexin Chen, William Dawson, Aulia Sukma Hutama, Takahito Nakajima, Keisuke Kameda, Manabu Ihara, Sergei Manzhos

Chemical Physics · Cross-listed from quant-ph

Classical Active-Space Hybrid Quantum Subspace Expansion (CASH-QSE): Quantum Corrections without Remeasuring the Classically Calculable Energy

The wave function prepared by a variational quantum eigensolver (VQE) can contain a substantial classically tractable component, yet its energy contribution is sampled on the quantum device. Representing the wave function as a combination of a classical reference and separately prepared quantum states allows the reference energy to be evaluated without quantum sampling. However, near-linear dependence among these components can amplify measurement errors and destabilize the resulting energy estimate. We introduce Classical Active-Space Hybrid Quantum Subspace Expansion (CASH-QSE), which retains a complete-active-space self-consistent field (CASSCF) reference classically and uses occupation structure to construct quantum states exactly orthogonal to the reference and to one another. The energy follows from an ordinary Hermitian eigenvalue problem without overlap measurements, while the same occupation constraints simplify the measured operators. Using full configuration interaction to guide component selection, we test CASH-QSE along H$_2$O and N$_2$ bond-stretching coordinates spanning weakly to strongly correlated regimes. The benchmarks use STO–3G for both molecules and a restricted cc-pVDZ orbital space for H$_2$O. CASH-QSE reaches chemical accuracy while limiting the largest complete measurement circuits to a few hundred all-to-all logical controlled-NOT gates. In favorable cases, it also reduces the idealized final-energy sampling cost by several orders of magnitude relative to VQE with an adaptive derivative-assembled pseudo-Trotter ansatz (ADAPT-VQE).
Read the full abstract

Chemical Physics

Construction of downfolded Hamiltonians from projective transcorrelation

Projective transcorrelation recovers the short-range electron correlation by a similarity transformation with a geminal function $f(r_{12})$, at the cost of an effective Hamiltonian containing a 3-body term. We replace the term by an effective operator of rank at most two, obtained from the two-body cumulant (2C) approximation to the three-particle reduced density matrix. The 2C 3-body energy is written without cumulants, and the effective one- and two-body interactions are derived as its partial derivatives with respect to the reduced density matrices. The truncation is assessed on the atomization and reaction energies of the HEAT set with CCSD(T), and on the CAS-pTC model, whose downfolded Hamiltonian is held on a qubit register with the number of Pauli strings reduced from the sixth power of the orbital count to the fourth.
Read the full abstract

Monday 7 September

Chemical Physics · Cross-listed from cond-mat.soft

Substrate-induced crystallisation in polymers: Exceptional impact of prefreezing versus heterogeneous nucleation on kinetics and morphology

Crystallisation of liquids often occurs at solid surfaces, via heterogeneous nucleation or prefreezing. Understanding this process is crucial for controlling crystal morphology in thin films and developing efficient nucleating agents - especially for polymers, typically processed in a molten state. However, the efficiency of nucleation and prefreezing with respect to crystallisation kinetics and semicrystalline polymer morphology has not been well studied and remains unclear. Here, we address this challenge by conducting a systematic study of substrate-induced crystallisation in dewetted droplets of the model polymer poly($\varepsilon$-caprolactone) (PCL) on various crystalline and amorphous substrates. Using polarized light optical microscopy, droplet crystallisation temperature was monitored and varied across the substrates, revealing three scenarios: close-to-homogeneous nucleation on glass and silicon; heterogeneous nucleation on ordered and disordered carbon-rich substrates; and prefreezing on graphene. The temperature-dependent nucleation rate was analysed and found to obey classical nucleation theory for all systems except PCL on graphene, where it diverges as all the droplets crystallise simultaneously. The highest in- and out-of-plane PCL crystal orientation is found on graphene, where the epitaxially grown lamellae extend vertically over several hundred nanometres. These results demonstrate that prefreezing has a superior influence on crystallisation kinetics and crystal morphology, making prefreezing substrates perfect nucleating agents.
Read the full abstract

Chemical Physics

Hemoglycin as a Candidate for Permanent Memory

The Hemoglycin polymer of glycine derived from space is terminated by iron and silicon to form a three-dimensional regular lattice with a high density of silicon tetrahedral vertices. Divalent cation bridges bypass vertices to facilitate electron conduction axially from layer to layer of the lattice. In principle, the exact distribution of Mg or Ca ions around each vertex can encode a very large amount of information. In space where the polymer forms in a molecular cloud the current derived from photon input will fade but in-fall hemoglycin in life forms on habitable planets like Earth, could be contained behind lipid electronic barriers of the nervous system in neurons or glia. Energy minimization of the permutations for magnesium or calcium occupancy showed deeper binding for calcium-dominated configurations. A low energy state with two independent magnesium to calcium substitutions appears to be key to the maintenance of long-term chemical modification. This state could be read via a pulse train that induced the matching spatial distribution of calcium to render adjacent insulating lattice regions conductive and provide a completed conduction path, signaling recognition. Writing of information may involve spatial modulation of the Na/Ca exchanger within the lipid membrane.
Read the full abstract

Chemical Physics

Simulation-Supervised Foundation Models for Retention Time Prediction in High-Performance Liquid Chromatography beyond Experimental Data Coverage

Accurate prediction of high-performance liquid chromatography (HPLC) retention times (RTs) across diverse molecules and chromatographic methods remains challenging because experimental training data cover only a limited region of chemical and method spaces. Here, we develop FUSE-RT (Foundation model Unifying Simulation and Experimental supervision for Retention Time), a multitask foundation model that integrates RT data from 179 chromatographic methods and adapts to unseen molecules and methods using limited target-domain data. To extend transferability beyond experimental coverage, we introduce simulation-to-real (Sim2Real) transfer learning, in which molecular representations learned from large-scale computational data are transferred to experimental RT prediction. Specifically, we use PolyOmics, comprising 39 properties for approximately 21,400 molecules generated by molecular dynamics and density-functional theory calculations, as auxiliary supervision. We evaluate generalization under molecular, method, and joint molecular–method distribution shifts. Simulation-derived supervision substantially improves transfer beyond the experimental molecular domain, particularly under pronounced coverage gaps and few-shot adaptation. Moreover, RT-prediction error decreases systematically with increasing simulation-data size, following a significant power-law relationship. These results establish Sim2Real transfer as a scalable strategy for extending RT prediction beyond the finite coverage of experimental chromatographic data.
Read the full abstract

All authors Stephen Wu, Yufeng Han, Yasuhiro Mito, Yoshiyuki Watabe, Yoshihiro Hayashi, Hikaru Takaya, Takuya Kubo, Ryo Yoshida

Chemical Physics

Low-rank approximation of Moment Tensor Potential enables reducing training set size without loss of accuracy

In this study, we implement a low-rank approximation of Moment Tensor Potential (MTP) based on the tensor train (TT) decomposition. The implemented tensor-factorized MTP (TFMTP) and the original MTP model are actively trained via a MaxVol-based algorithm during molecular dynamics simulations of a four-component molten salt mixture, LiF-NaF-KF (FLiNaK), and geometry optimizations of a five-component equiatomic MoNbTaWV random alloy. We demonstrate that under a 1.5-fold compression, TFMTP requires two times fewer configurations for fitting than the original MTP model, while maintaining an indistinguishable level of accuracy. These actively trained MTP and TFMTP models are further used to evaluate the density and viscosity of FLiNaK at temperatures ranging from 600 to 1200 K, as well as the elastic constants and bulk modulus of the MoNbTaWV alloy at zero temperature. For both atomic systems, the differences in physical properties predicted by MTP and TFMTP are negligible.
Read the full abstract

All authors Anna Bondarenko, Nikita Rybin, Maxim Rakhuba, Ivan S. Novikov

Chemical Physics

Strong Evidence for Formation of Hydroxyl Anion via 2-Particle-1-Hole Feshbach Resonances

This study investigates the formation of the hydroxyl anion (OH$^-$) via dissociative electron attachment in 2-propanol as a model system for studying both organic and inorganic molecules. Using high-level CAP-EOM-EA-CCSD calculations and advanced ToF mass spectrometry, we demonstrate that OH$^-$ formation at electron energies between 7 and 11 eV is dominated by two-particle-one-hole (2p-1h) Feshbach resonances. The potential-energy curves reveal a dense manifold of anionic states coupled through numerous avoided crossings, facilitating nonadiabatic population transfer during C-OH bond dissociation. Survival-probability analysis identifies a subset of six long-lived resonances that persist long enough to drive fragmentation, with states 25 and 28 acting as primary drivers by funneling the attached electron into the localized $\sigma^*(\mathrm{C{-}OH})$ antibonding orbital. These theoretical predictions are confirmed by experimental observations of a prominent OH$^-$ yield peaking at 8.6 eV, supporting a site-specific fragmentation mechanism that generalizes to the broader class of molecules.
Read the full abstract

All authors Siddique Ali, Meenakshi Rana, Soumya Ghosh, Narayan Kundu, Aryya Ghosh, Dhananjay Nandi