2026-03-09 | | Total: 16
We present two models with explicit long-range electrostatics in the form of Coulomb interactions. Both models include point charges depending on their local atomic environments, and the second model also conserves a total charge of an atomic system. We combine the proposed long-range models with local Moment Tensor Potential and demonstrate that they reduce the training errors of the MTP models fitted on the same training sets including the CH$_3$COO$^-$+4-methylphenol and CH$_3$COO$^-$+4-methylimidazole organic dimers (non-periodic systems) and the NaCl crystal (periodic system). For the organic dimers, the proposed models also give qualitatively correct predictions of the binding curves. Furthermore, in this study we introduce a method for calculating phonon spectra of isotropic materials only via these long-range models fitted to energies, forces, and stresses. The developed long-range model with point charges dependent on atomic environments and conserving total charge is capable of predicting the correct value of the LO-TO splitting in the $Γ$-point in the isotropic NaCl. For this system, we also predict dielectric constant from dipole moment fluctuations calculated with molecular dynamics simulations conducted with the developed long-range model. The calculated dielectric constant is in good agreement with experiment. Finally, we demonstrate the broader applicability of the introduced approach by computing the phonon spectrum of uniaxial tetragonal PbTiO$_3$. Although the method is formally derived for isotropic materials, we show that it is also perspective for uniaxial materials (e.g., PbTiO$_3$) as the spectrum obtained with our long-range interatomic potential corresponds to the one calculated with density functional theory.
Methane is the simplest hydrocarbon, yet it exhibits an extraordinarily complicated series of crystal phases. Notably, the non-plastic phases have large unit cells with nearly, but not quite cubic symmetry. Furthermore, although non-polar molecules interact very weakly, their reorganisation across phase transitions is very sluggish. Here, we demonstrate that these complex structures can be understood as simple packing of near-spherical supermolecular clusters of methane molecules: the departure from cubic symmetry arising from the non-spherical nature of the molecules. We use molecular dynamics based on density functional theory calculations to simulate the finite-temperature crystal structures of methane, finding that the complex Phase A is based around a 13-molecule regular icosahedron, with 8 additional molecules forming the 21-molecule unit cell. Similarly, Phase B is based on a body-centred cubic bcc packing of 17-molecule Z16 polyhedra, with the remaining 12 molecules per cell in tetrahedral interstices. We demonstrate that the favored intermolecular separation depends sensitively on molecular orientation, leading to hindered rotation and suppressed entropy. The structures are determined by a trade-off between efficient packing and entropy.
Particle in Cell (PIC) simulations have become a vital tool for the investigation of kinetic processes in plasma physics. Many of the systems investigated with PIC simulations contain particles with relativistic velocities. The correct integration and the knowledge of possible sources of errors in relativistic particle trajectories is of importance to accurately judge the validity of the simulation results. Over the past few decades, various new integration schemes for relativistic particle trajectories in PIC simulations have been proposed. These are aimed at improving numerical accuracy in specific scenarios. This article presents a comprehensive comparison of particle pushers with a focus on explicit schemes. An important class of these schemes is found to be generalisable to arbitrary high order. A comparison of the fourth order variants of these schemes with their second order counterpart is also presented.
Semiconducting devices face persistent challenges in operating at high pressure, as the band theory predicts that materials transition to a more metallic state under compression. However, our findings with carbon chains in KFI substrates reveal a conditional deviation from this norm. We not only witness the transition from polyyne (semiconductor) to cumulene (metal) at medium pressure, but we also observe an unexpected re-entrance of the polyyne at high pressures, where the band gap in the polyyne increases with pressure. In addition, the synthesis of long cumulene chains has posed a longstanding challenge in the quest for high-temperature organic superconductivity. We have identified critical conditions for synthesizing extended cumulene chains within zeolite frameworks, highlighting the interplay between unconventional charge density waves and significant torsions. The KFI zeolite facilitates the formation of carbon chains exceeding 5,000 atoms, in stark contrast to around 100 other zeolites that are limited to ~10 atoms. The cumulene@KFI system demonstrates a superconducting transition temperature reaching ~62 K, surpassing the highest reported values for bulk iron-based superconductors. This interplay between carbon structures and superconductivity not only advances our understanding of charge density waves but also heralds a new era in the study of novel applications
Transforming rovibronic Hamiltonians of molecular systems from the $ΛS$ (Hund's case a) basis to the adiabatic $Ω$ representation is widely used to "remove" spin-orbit coupling (SOC) and enable single-state treatments of spectra and dynamics. We show that this simplification is only apparent: the SOC elimination necessarily generates sizeable non-adiabatic couplings (NACs) from the nuclear kinetic energy operator. Neglecting these spin-orbit-induced NACs causes severe errors in rovibronic energies and transition properties. Using an analytically tractable two electronic state model and high-accuracy variational benchmarks, we derive the exact conditions for numerical equivalence between $Ω$ and $ΛS$ formulations and quantify how missing NAC terms and bond-length-dependent spin factors degrade predictions. We implement a complete $Ω$-representation workflow in Duo for diatomics, fully transforming all Hamiltonian terms and enabling side-by-side $Ω$ vs $ΛS$ calculations. For common single-state pipelines (e.g., LEVEL), we provide diagnostics that flag unsafe regimes and practical remedies to restore accuracy. The results deliver actionable guidance for spectroscopy, photophysics, and kinetics: $Ω$-based single-state approximations are reliable only when interacting states are well separated in the Franck-Condon region; otherwise, explicit non-adiabatic terms are required - even for "forbidden" transitions.
This study investigates the use of machine learning based mesh adaptivity, specifically mesh movement methods (UM2N), with depth integrated non-hydrostatic shallow water models. Motivation for this comes from the need for models which balance efficiency and accuracy for use in probabilistic coastal hazard assessment. Implementations are built on the discontinuous Galerkin finite-element (DG-FE) based software, Thetis, which leverages the partial differential equation (PDE) framework Firedrake for automated code generation. Verification on benchmark test cases and validation against laboratory measurements of coastal hazards, focusing on tsunami propagation, run-up, and inundation is performed. In these tests, the UM2N-driven meshes help resolve key non-hydrostatic dynamics and yield numerical solutions in close agreement with reference computations and measured data. Numerical results indicate that the UM2N surrogate based approach significantly accelerates conventional mesh movement techniques and has high robustness over long integration periods and under strongly nonlinear wave conditions.
Accurate quantum chemical calculations are critical for understanding molecular properties, yet their computational cost remains a major challenge. Full Configuration Interaction (FCI) provides exact solutions but is prohibitively expensive for large systems. To address this, quantum computers are expected to be useful, but developing practical quantum computers is still ongoing. Here we introduce an efficient Configuration Interaction (CI) computation algorithm based on classical mechanics, which we call Simulated Bifurcation-based CI (SBCI), because we derive this algorithm from a quantum inspired algorithm for combinatorial optimization called Simulated Bifurcation. Applying it to FCI computations of representative molecular systems and comparing the results with those by a standard method, we demonstrate that SBCI can reduce computation costs such as computation times and/or required memory sizes, while keeping high accuracy comparable to the standard method. Thus, SBCI will be promising for accelerating high-precision electronic structure calculations without compromising reliability.
Many HPC applications that solve differential equations rely on the Runge-Kutta family of methods for time integration. Among these methods, the fourth-order accurate RK4 scheme is especially popular. This time integration scheme requires applications to evaluate four intermediate stages to take one time step. Depending on the complexity of the problem being solved, the evaluation of these intermediate stages can be computationally expensive. In this paper we develop explicit fourth-order accurate Multistep Runge-Kutta (MSRK) methods. The advantage of such methods is that they re-use data from previous time steps, thus requiring fewer intermediate stage evaluations and potentially speeding up applications. We outline a procedure to obtain and tune the method's coefficients by adjusting their stability regions in an attempt to maximize the size that a time step can take. We validate and evaluate our new methods in the context of Numerical Relativity applications using the EinsteinToolkit. We believe, however, that these methods and results should generalize to other applications using explicit Runge-Kutta methods.
Next-generation microreactors enable remote deployment and semi-autonomous operation, but compact, sealed, heterogeneous cores limit conventional safeguard approaches that rely on access and bulk accountancy. Limited inspection access and complex internal geometry reduce sensitivity to localized anomalies such as missing fuel. Here we demonstrate missing-fuel detection in microreactor scale geometries using muon scattering tomography under realistic cosmic-ray conditions. We introduce $μ$TRec, a physics-informed framework that reconstructs event-level curved muon trajectories by combining a Gaussian multiple Coulomb scattering model with Bayesian updating, then maps scattering density through voxel wise M-values for core integrity verification. We evaluate a representative hexagonal core containing 61 fuel flakes with embedded control drums and shutdown rods, using both idealized 5 GeV muons and zenith-angle-dependent 0-60 GeV cosmic-ray spectra. A single missing fuel flake is detected with $3\times 10^{6}$ muons at 50 mm voxel resolution. Incorporating per-muon momentum further increases detectability by up to 149.85% for laser-driven sources and 105.11% for cosmic-ray sources relative to momentum-agnostic reconstruction. The approach remains robust under practical detector limits, with only an 8.88% reduction in detectability for 10 mm spatial resolution and 10% energy resolution. Compared with PoCA, $μ$TRec delivers 326.13% to 392.14% higher detectability at equal muon counts, enabling faster defect identification.
We investigate the impact of tokeniser pretraining on the accuracy and efficiency of physics emulation. Modern high-resolution simulations produce vast volumes of data spanning diverse physical regimes and scales. Training foundation models to learn the dynamics underlying such data enables the modelling of complex multiphysics phenomena, especially in data-limited settings. The emerging class of physics foundation models typically aims to learn two tasks jointly: (i) extracting compact representations of high-resolution spatiotemporal data, and (ii) capturing governing physical dynamics. However, learning both tasks from scratch simultaneously can impede the effectiveness of either process. We demonstrate that pretraining the tokeniser with an autoencoding objective prior to training the dynamics model enhances computational efficiency for downstream tasks. Notably, the magnitude of this benefit depends on domain alignment: pretraining on the same physical system as the downstream task yields the largest improvements, while pretraining on other systems provides moderate gains. In-domain pretraining reduces VRMSE by 64% after 10,500 training steps compared to training from scratch. To our knowledge, this is the first systematic investigation of tokeniser pretraining for physics foundation models. We further introduce flexible spatiotemporal compression operations that extend causal convolutions to support runtime-adjustable compression ratios, enabling efficient adaptation to diverse downstream tasks. Our findings provide practical guidance for training efficient physics emulators and highlight the importance of strategic pretraining data selection.
We present the first application of spectral nudging in a probabilistic ensemble forecasting framework, combining the physics-based ECMWF Integrated Forecasting System ensemble (IFS-ENS) with forecasts from the probabilistic machine-learned AIFS-ENS ensemble. Large scales of virtual temperature and vorticity are relaxed toward the machine-learned forecasts, while mesoscale structures remain governed by the physics-based model. This hybrid ensemble shows substantial improvements in large-scale forecast skill, with gains in predictive skill extended by up to two days in the tropics and by approximately half a day in the extra-tropics relative to IFS-ENS. Despite nudging being applied only to upper-air fields, improvements are also found in several near-surface parameters. Tropical cyclone track forecasts improve significantly, consistent with improved representation of the large-scale steering flow, without degrading storm intensity or ensemble spread. These results demonstrate that spectral nudging can be successfully extended to ensemble prediction and provide an attractive pathway for combining machine-learned and physics-based weather prediction systems.
We investigate the Continuous-Time Koopman Autoencoder (CT-KAE) as a lightweight surrogate model for long-horizon ocean state forecasting in a two-layer quasi-geostrophic (QG) system. By projecting nonlinear dynamics into a latent space governed by a linear ordinary differential equation, the model enforces structured and interpretable temporal evolution while enabling temporally resolution-invariant forecasting via a matrix exponential formulation. Across 2083-day rollouts, CT-KAE exhibits bounded error growth and stable large-scale statistics, in contrast to autoregressive Transformer baselines which exhibit gradual error amplification and energy drift over long rollouts. While fine-scale turbulent structures are partially dissipated, bulk energy spectra, enstrophy evolution, and autocorrelation structure remain consistent over long horizons. The model achieves orders-of-magnitude faster inference compared to the numerical solver, suggesting that continuous-time Koopman surrogates offer a promising backbone for efficient and stable hybrid physical-machine learning climate models.
The acoustic signature of thermal spray processes is known to vary with changes in operating conditions, which also influence particle in-flight velocity and distribution. Building on this idea, the present work first develops an analytical model that links chamber and nozzle parameters to far-field acoustic levels using gas-dynamics relations and simplified acoustic power propagation. The model is then calibrated to reduce systematic error associated with neglected turbulence effects and to improve agreement across operating conditions. In addition, a numerical framework is implemented to complement the analytical model and to resolve supersonic jet flow and in-flight particle transport. The second part of the study uses unsteady compressible simulations with hybrid turbulence modeling such as Unsteady Reynolds-Averaged Navier-Stokes (URANS) and Delayed Detached Eddy Simulation (DDES) to capture the development of the shock-containing jet and the associated near-field pressure fluctuations. Far-field sound is predicted using the Ffowcs Williams-Hawkings acoustic analogy, while a Lagrangian approach tracks particles injected at the nozzle exit to quantify velocity evolution, radial spreading, and downstream flux distributions. The influence of operating conditions (e.g., chamber pressure and temperature) is assessed, and predictions are evaluated against published microphone spectra and particle flux measurements. Overall, the combined analytical and numerical approach captures how changes in nozzle operating conditions affect jet unsteadiness and mixing, leading to measurable shifts in acoustic level and spectral content. These results suggest that aeroacoustic signatures could be used as a non-intrusive pathway to monitor and potentially control thermal spray operating conditions.
Data-driven surrogate models improve the efficiency of simulating continuous dynamical systems, yet their autoregressive rollouts are often limited by instability and spectral blow-up. While global regularization techniques can enforce contractive dynamics, they uniformly damp high-frequency features, introducing a contraction-dissipation dilemma. Furthermore, long-horizon trajectory optimization methods that explicitly correct drift are bottlenecked by memory constraints. In this work, we propose Jacobian-Adaptive Weighting for Stability (JAWS), a probabilistic regularization strategy designed to mitigate these limitations. By framing operator learning as Maximum A Posteriori (MAP) estimation with spatially heteroscedastic uncertainty, JAWS dynamically modulates the regularization strength based on local physical complexity. This allows the model to enforce contraction in smooth regions to suppress noise, while relaxing constraints near singular features to preserve gradients, effectively realizing a behavior similar to numerical shock-capturing schemes. Experiments demonstrate that this spatially-adaptive prior serves as an effective spectral pre-conditioner, which reduces the base operator's burden of handling high-frequency instabilities. This reduction enables memory-efficient, short-horizon trajectory optimization to match or exceed the long-term accuracy of long-horizon baselines. Evaluated on the 1D viscous Burgers' equation, our hybrid approach improves long-term stability, shock fidelity, and out-of-distribution generalization while reducing training computational costs.
Sampling molecular conformations from the Boltzmann distribution is essential for computational chemistry, but iterative diffusion methods are prohibitively slow. Drifting Models offer one-step generation, yet their equilibrium matches the \emph{training} distribution, which may deviate from the true Boltzmann distribution due to sampling bias. We introduce Drifting Models to molecular conformation generation for the first time, establishing a theoretical bridge via the \emph{Drifting Score Identity}: for Gaussian kernels, the drifting field's attraction equals a kernel-weighted average of \emph{any} distribution's score function. Substituting molecular force labels -- which directly encode the Boltzmann score -- yields the \emph{Drifting Force Identity} and decomposes the field into standard drift plus a Boltzmann correction. We further discover a striking phenomenon unique to molecular systems: force incorporation's effectiveness \emph{reverses across representations}. In coordinate space, Force-Interpolated Drifting (FI) dominates by blending physical force directions with data displacements. In distance feature space, Force-Aligned Kernel (FK) achieves superior accuracy by modifying only kernel weights, thereby preserving the manifold of geometrically valid molecules. On MD17 Ethanol, both approaches achieve one-step generation with over 1000x speedup relative to recent score-matching methods with Boltzmann guiding, providing more than million-fold acceleration over traditional molecular dynamics, while ensuring perfect structural validity and distributional accuracy rivaling multi-step methods.
We introduce a data-driven framework for approximating the convex set of $N$-representable two-electron reduced density matrices (2-RDMs). Traditional approaches characterize this set through linear matrix inequalities that define its supporting hyperplanes. Here, we instead learn a vertex-based approximation to its boundary from molecular data and use this information to improve the set defined by low-order positivity constraints, without explicitly constructing higher-order conditions. The resulting semidefinite machine learning approach -- combining an input convex neural network with semidefinite programming -- drives a direct variational calculation of the 2-RDM with enhanced accuracy at computational cost comparable to two-positivity calculations. Applications to the potential energy curves of ${\rm C}_2^{2-}$, ${\rm N}_2$, and ${\rm O}_2^{2+}$ demonstrate these systematic improvements as well as close agreement with complete active space configuration interaction results. Overall, semidefinite machine learning interweaves data-driven boundary information with semidefinite positivity constraints to yield more accurate energies and 2-RDMs without explicit higher-order positivity conditions.