2026-07-28 | | Total: 17
We provide recommended values for the second virial coefficient, $B(T)$, and its uncertainty, for molecular nitrogen and oxygen. The temperature range covered is $20-3000$ K for nitrogen and $20-2000$ K for oxygen. The recommendations are based on tuning previously published state-of-the-art ab initio pair potentials so that the $B(T)$ calculated from the potentials match selected high-accuracy experimental data; for nitrogen the tuning utilizes values of $B$ derived from literature density data with greatly reduced uncertainty by analyzing the data with the aid of ab initio calculated higher virial coefficients. Quantum effects on $B$ are fully included with the path-integral Monte Carlo method. The resulting $B(T)$ have uncertainties similar to those of the best experimental data, but cover a much wider temperature range.
Exact simulation of high-resolution NMR spectra requires block diagonalization of the spin Hamiltonian, whose dimension grows exponentially with the number of spins $N$; symmetry is the principal tool for taming this growth, yet which permutation groups can occur as the full symmetry group of a scalar-coupled spin system has lacked an exhaustive treatment. Formulating the spin system as an undirected edge-weighted complete graph, we prove an exact realizability criterion: a subgroup of $S_N$ is realizable if and only if it coincides with its symmetrized (undirected) Wielandt 2-closure. In particular, purely rotational symmetry of a single spin ring is impossible, yet chiral spin systems do exist as multi-orbit twisted stacks, and we determine the minimal spin count $μ^{*}(C_n)$ for every cyclic group, including the counter-intuitive realizations $C_8$ and $C_9$ at $N = 12$. A sequential symmetrization algorithm, completed by an orbit-partition decomposition, yields a provably exhaustive enumeration of all realizable symmetry types up to $N = 14$: the apparently new sequence $a(N) = 1, 1, 3, 8, 11, 27, 36, 90, 131, 282, 394, 948, 1316, 2866$ with the tower law $a(N) = a(N-1) + f(N)$ - a catalogue of 6112 entries in all, organized by canonical identifiers and a structural grammar extending the Pople nomenclature. Finally, we present a hierarchical methodology for exact block diagonalization without physical approximations: factorization by the conserved total spin projection, Schur-Weyl contraction of magnetically equivalent composites, orbit-weight deduplication of the spin configurations, and isotypic projection over the representations of the factor group, with a uniform treatment of non-abelian groups and complex characters.
Fanpy is a Python library for developing new wavefunction methods. It enables users to quickly convert mathematical expressions into working code through a modular design based on the Flexible Ansatz for N-electron Configuration Interaction (FANCI) theory. This architecture facilitates a straightforward extension of the codebase. Here we present version 2.0 of the Fanpy package. This release includes several new wavefunction implementations, including coupled- cluster-inspired geminal approaches. A new analysis module enables a more detailed inspection of computational results and lays the groundwork for future features. In addition, the PySCF interface has been redesigned, and an interface to the PyCI package has been introduced to offload computationally expensive components. Finally, we introduce an improved software development environment, including automated testing and issue tracking.
One of the distinguishing aspects of CP2K is its seamless integration of diverse structural and transition-state optimization techniques with advanced sampling approaches including Monte Carlo, molecular dynamics, and metadynamics, enabling the efficient exploration of complex potential- and free-energy landscapes, including rare events. These capabilities are combined with a broad hierarchy of energy and force evaluation methods, ranging from classical and machine-learned interaction potentials and mixed quantum-classical multiscale and semiempirical schemes, to highly accurate quantum-mechanical electronic-structure approaches. At the heart of the latter lies the Gaussian and plane-wave framework, along with its augmented all-electron generalization, which have been described in detail in our previous code review [T. D. Kühne et al., J. Chem. Phys. 152, 194103 (2020)]. Building on this foundation, the present work revisits the methods within CP2K that turn electronic structure into dynamics, transport, and spectroscopic response. Particular emphasis is placed on the coupling between static response calculations and nuclear motion: spectra may be evaluated at optimized structures, averaged over thermally sampled configurations, obtained from time-correlation functions along ab-initio or path integral molecular trajectories, or followed in real time together with electronic and nuclear dynamics. The same modular structure also enables equilibrium and biased transport simulations, from Kubo-type linear response to open-boundary approaches under external potentials, highlighting CP2K's unique capability to unify quantum chemistry with quantum and statistical mechanics within a versatile, holistic simulation environment.
In this work, we present a novel perspective on the coupling force employed to compensate the interface artifacts prevalent in adaptive resolution simulations (AdResS) of open many-particle systems. We show that a substantial part of this "thermodynamic force" can be framed in terms of the theoretical model of the Liouville-type hierarchy for open systems. The correspondence is made explicit for the case of a simple atomistic fluid, for which a one-dimensional integral expression is derived. This enables the analysis for dependencies of the thermodynamic force on important simulation parameters, which is taken to inspire both simplifications for the numerical calculation of the thermodynamic force and new criteria for its validation that are adequate to the interfacial nature of the problem. The theoretical claims are then verified in a simulation study of the atomistic supercritical Lennard-Jones fluid at different thermodynamic states.
Heterogeneous density functional theory (DFT) calculations, particularly plane-wave implementations, introduce systematic formation energy errors ranging from tens to hundreds of meV/atom, depending on the selection of exchange-correlation functionals, kinetic energy cutoffs, pseudopotentials, and dispersion corrections. As demonstrated by the MatPES dataset, identical structures can exhibit an average energy discrepancy of 107 meV/atom between PBE and r2SCAN calculations. Such method-dependent discrepancies hinder the integration of multi-source DFT data, greatly limiting the scale and quality of datasets for training robust materials AI models. Here, we resolve this fundamental data silo barrier via graph-based transfer learning. Leveraging 380,190 structurally paired PBE-r2SCAN entries from the MatPES database, we train a structure-aware graph neural network to predict cross-functional energy residuals and align inconsistent DFT energy scales. By adopting GPTFF model architecture, the model converts conventional PBE energies to r2SCAN-level accuracy with a mean absolute error of 14.3 meV/atom, compared with 18.2 meV/atom achieved by CHGNet. This versatile approach effectively upgrades massive legacy PBE datasets to high-precision r2SCAN standards. It enables reliable predictions of phase stability, battery voltage profiles, and reaction thermodynamics, while allowing the integration of multi-source DFT data to advance the development of high-performance materials foundation models.
Cold polar molecules offer fascinating prospects for ultracold chemistry and quantum physics, including new platforms for quantum simulation or computation. However, their inherent collision properties remain largely unknown. It has proven extremely hard to experimentally probe collisions between two dipolar molecules at sufficiently low energies and high precision, as it appears fundamentally impossible to merge two beams of molecules with significant dipole moments. Here we report measurements of state-to-state cross sections for collisions between ammonia isotopologues at energies between 0.3 and 100 cm$^{-1}$ using a novel beam merger. We experimentally observed a local maximum in the cross sections that indicates the effective dipole moments can switch off at low collision energies. Scattering calculations reproduced this maximum in good agreement and explained the observed scaling with the parity splitting energies in the molecular energy level structures. Measurements of the correlated energy transfer in both collision partners yielded direct evidence of the suppression of the dipole-dipole interaction at low energies. Our results demonstrate how collisions between an important class of polar molecules evolve from the high temperature limit towards the ultracold regime in a counterintuitive way, have major consequences for the feasibility of future experiments and the interpretation of previous work, and offer distinctive opportunities to control cold molecular collisions with external fields.
We characterise the interplay between cation-translational and anion-rotational dynamics on the model system $\rm α\text{-}Li_{2}SO_{4}$ through virtual isotopic substitution molecular dynamics (VISMD) simulations. We independently control in turn the lithium translational- and sufate rotational- diffusivity by varying the masses of lithium and oxygen, respectively, showing that both diffusivities exhibit power-law relationships with the masses of each individual component. Changing either diffusivity produces an effect on the other; i.e. there exists a bidirectional coupling between anion and cation diffusive dynamics. From our variable-temperature/variable-mass dataset we demonstrate that the apparent activation energy for lithium diffusion depends strongly on the isotope mass, decreasing with accelerated sulfate dynamics; an observation which provides clear evidence of mode-coupling between cation and anion dynamics.
We scale a CUDA-Q-native pipeline coupling a generative quantum eigensolver (GQE) to quantum-selected configuration interaction (QSCI) across active spaces of 14 to 44 qubits, applied to the extreme-ultraviolet (EUV) photoresist chemistry of monoalkyltin oxo-hydroxides. A GPT-2 policy emits UCCSD operator sequences; sampled bitstrings become determinants, diagonalised classically, and a cross-circuit generalised-eigenvalue refinement makes every reported GQE+QSCI energy a variational upper bound. Every rung from 14 to 40 qubits carries an exact CASCI or FCI reference, up to 166 million determinants for SnO at 32 qubits. The pipeline is chemically accurate, below 1.6 mHa, through 30 qubits on methyltin trihydroxide and through 32 on SnO, on the best seed at the top rungs. Circuit depth rather than training length is the scaling lever; the refined subspace grows near-linearly with the operator count while staying a vanishing fraction of the determinant space, 0.017% at the 32-qubit SnO rung. It also runs on the 54-qubit IQM Emerald processor, at the shallow depths its routed two-qubit gates allow, reaching +0.330 mHa for SnO at 14 qubits from a CCSD-amplitude-ordered pool prefix and +3.92 mHa for the industrial n-butyltin ligand at 22 qubits from depth-truncated trained circuits under per-circuit readout self-calibration, 81% of the active-space correlation; classical configuration recovery on those counts tightens the 22-qubit result to +0.18 to 0.21 mHa. For the methyl resist, ionisation collapses the classical UCCSD(T) Sn-C bond dissociation energy from 72.6 to 21.2 kcal/mol, the switch that flips solubility on exposure. Against that, the 40-qubit result is support-limited at 22.8 mHa, the full trained ansatz on hardware awaits better fidelities, and classical subspace expansion reaches the 32 to 40-qubit spaces with no quantum sampler, so that boundary is mapped, not beaten.
We investigate the Uhlenbeck-Ford (UF) model as a reference system for free-energy calculations in two-dimensional (2D) fluids. The 2D virial coefficients are computed exactly up to tenth order and combined with molecular simulation data to construct highly accurate numerical representations of the equation of state and the excess Helmholtz free energy.We then determine the phase diagram of the model in order to establish the thermodynamic stability limits of the fluid phase and thereby identify the range of applicability of the UF model as a fluid reference system. In the course of this analysis, we identify the solid, hexatic, and fluid phases, and show that the fluid remains the only thermodynamically stable phase, independent of density, for scaling parameters up to $p\lesssim 70$. Finally, we demonstrate the practical applicability of the 2D UF model as a reference system through thermodynamic integration calculations of the free energy of a two-dimensional Lennard-Jones fluid.
We present a method that enables efficient simulations of coarse-grained electrolyte solutions inside a narrow slit pore with an anisotropic dielectric permittivity tensor. The electrostatic equations for polarizable surfaces are solved using a 2D periodic Green's function method combined with a slab-corrected anisotropic 3D Ewald summation. We apply this approach in Monte Carlo simulations to study 1:1 electrolytes confined between both polarizable dielectric and metallic surfaces. Our results show that dielectric anisotropy aggressively reshapes the double-layer structure. While a coordinate stretching transformation demonstrates that individual ion-image interactions depend strictly on the bulk-like parallel permittivity, the suppression of the perpendicular permittivity dramatically amplifies direct in-plane ion-ion correlations. Under strong anisotropy, these lateral correlations dominate the thermodynamics completely, rendering the structural profiles of mutually opposing dielectric and metallic boundaries practically identical by forcing the smaller cations directly into the contact plane of the larger anions.
Attention enables context modeling via query-key scoring with softmax normalization. Driven by industrial long-context demands, mainstream research has converged toward sparsity and efficiency--yet softmax's independence assumption persists. For scientific tasks unburdened by long-token constraints, however, richer structured coupling may often be essential, making tailored attention both viable and more appropriate. To this end, we propose Variational-Ising-Attention (VIA), which augments softmax normalization with an interacting Ising model; attention patterns emerge from learnable pairwise couplings via variational mean-field inference, redefining attention from a ranking over isolated items to a collective state over interacting entities. We instantiate VIA on retrosynthesis reaction center prediction, a task inherently governed by cooperative bond-breaking constraints. Comprehensive experiments across model variants, coupled with mechanistic analyses, demonstrate that VIA consistently and substantially outperforms standard softmax attention. More broadly, our findings suggest that for scientific problems, the optimal solution is not general-purpose efficiency, but appropriately tailored attention aligned with intrinsic domain structure. This work provides a theoretically grounded and empirically validated instantiation of this paradigm.
In this work, we introduce a multiconfigurational mixed quantum-classical many-body approach for simulating the finite-temperature correlated multi-exciton dynamics in the presence of phonon-induced static and dynamic disorder. In this mixed quantum-classical approach, the excitonic subsystem is described using a multiconfigurational wavefunction that extends beyond the mean-field limit, while the phonons are evolved quasi-classically. Using this approach, we simulate a multi-excitonic dissipative system and show how the interplay between phonon-induced dynamic disorder and exciton-exciton many-body interactions determines excitation-dependent excitonic transport and spatial correlations. Our results show that while the mean-field approach produces semi-quantitatively accurate diffusive dynamics, it does not capture the spatial correlations as expected. We find that a mean-field path approximation, where we generate pre-computed trajectories using our mean-field mixed quantum-classical approach and then perform multiconfigurational dynamics, can reproduce the spatial correlations to a good accuracy, positioning this approach as an efficient method for capturing spatial correlations in complex systems.
Understanding the effects of doping on the crystal structure and optical properties of semiconductor materials is crucial for advancing next-generation semiconductor and photonic technologies. Although various studies have focused on doped hematite ($α$-Fe$_2$O$_3$), the relationship between dynamical stability and optical properties remains insufficiently explored. This study presents a comprehensive first-principles investigation that simultaneously evaluates the phonon dispersion characteristics and frequency-dependent optical response of B-doped, Y-doped, and (B, Y)-co-doped $α$-Fe$_2$O$_3$, providing deeper insights into the underlying mechanisms. We examined the finite-temperature vibrational properties, dielectric function, and optical characteristics to comprehend the lattice dynamics and light-matter interactions under electromagnetic radiation. Vibrational thermodynamics reveal that pristine and Y-doped hematite maintain dynamic stability, while B-doped hematite exhibits imaginary phonon modes indicating lattice instability due to distortions in the Fe--O framework. Notably, Y co-doping with B helps suppress these soft modes, restoring structural stability through lattice relaxation and improved interatomic forces. B doping enhances low-energy absorption by introducing additional states in the valence band, while Y doping alters orbital hybridization, leading to a broader dispersion. In the optical regime, doped hematite displays dominant interband transitions below $2$ eV and strong absorption between 1.80 eV and 4 eV. The (B, Y) co-doping combines the low-energy benefits with an improved optical response profile. In summary, doping significantly enhances lattice vibrations, light-matter interactions, and optical responses, providing an effective strategy for tailoring hematite for diverse applications in photoactive, optoelectronic, and photonic technologies.
Multi-ionic materials pose a distinct representational challenge in machine learning-driven materials design. Different from single-molecule or composition-based materials, their properties arise from how charged building blocks aggregate into specific assemblies. Here, we show how pretrained machine-learned interatomic potentials (MLIPs) can bypass full crystal-structure prediction and support pre-synthesis screening from stoichiometric ionic clusters using multi-ionic integrated explosives (MIXs) as a synthesis-facing example. This strategy combines a stoichiometric ionic-cluster representation, which represents each candidate material by a non-periodic, stoichiometry-preserved formula-unit cluster, with multi-task fine-tuning (MT-FT), which adapts a pretrained atomistic backbone while retaining the energy--force objective as physical regularization for the sparse detonation-velocity labels. With the pretrained backbone regularized by MT-FT, this surrogate provides a cross-validated screen across only 25 structurally curated perovskite-type energetic materials (PEMs) with experimentally derived Kamlet--Jacobs (K--J) detonation velocities. Representation probes show that the learned descriptors implicitly retain site-aware ionic organization, density information, and coarse packing compatibility, implying why non-periodic clusters can remain predictive before full crystal structures are known. The surrogate extends known PEMs chemistry to three newly synthesized ABX$_4$ materials with both unseen ABX$_4$ stoichiometry and an unseen ethylenediammonium B-site cation, yielding three-point concordance with K--J reference velocities and a mean absolute error (MAE) of 92~m$\cdot$s$^{-1}$ without retraining. Together, these results establish stoichiometry-preserved cluster learning as a synthesis-facing screening strategy for data-scarce multi-ionic materials.
We show how finite-difference time-domain (FDTD) simulations can be extended to model ultrafast nonlinear microscopy, enabling the prediction of spatially-resolved pump--probe signals in arbitrary electromagnetic environments. Focusing on polariton transport in strongly coupled light--matter systems, we develop a perturbative framework to study the ultrafast propagation of hybrid light--matter excitations in nanophotonic structures. We first apply the framework to a standard distributed Bragg reflector (DBR) cavity, reproducing established results for polariton transport from a multimode Tavis--Cummings model. We then consider the full modal landscape of the same cavity, including guided modes below the light line that are typically neglected in single-mode-family descriptions. Exploiting these modes, we design compact mode converters that transfer radiative cavity polaritons into photon-like guided polaritons and back, utilizing the guided modes for low-loss propagation. Despite molecular dephasing, this enables transport of molecular excitation over a millimeter, an order of magnitude beyond current transport experiments. We further compute the pump--probe differential transmission signal, providing an experimental signature of the mechanism. Our results show that the modal landscape of a photonic cavity can be engineered to bypass limitations commonly assumed to be intrinsic to the transport of molecular polaritons.
Quantum-centric workflows are a promising route to improving the accuracy of property predictions in computational chemistry and materials science. By integrating quantum sampling algorithms with classical solvers, electronic structure calculations have recently demonstrated their potential even on noisy intermediate-scale quantum devices. In principle, the method of Vibrational Configuration Interaction (VCI) is suitable for integration with quantum sampling algorithms as well. However, demonstrations of computational workflows for quantum-centric, vibrational property predictions are still lacking. Here, we introduce a methodology for performing anharmonic vibrational structure calculations that can be deployed in a hybrid, quantum-classical mode. Starting from a quartic force field, the approach combines a Vibrational Self-Consistent Field (VSCF) with VCI in either Full, Selected (S-VCI), or Symmetry-Adapted (SA-VCI) mode. In S-VCI, an Epstein-Nesbet perturbative screening significantly reduces the configuration space while retaining high predictive accuracy. A state-list input enables the integration of externally generated vibrational configurations as a seed space. As a proof-of-concept, we demonstrate a hybrid, quantum-classical computational workflow, in which a quantum sampling algorithm provides the seed. Our vibrational wave function analysis package ViBra, equipped with a graphical interface, is available at https://github.com/raphafe96/ViBra.