2026-10-20 | | Total: 19
High-accuracy computation of nonsmooth dynamics is a challenging problem in systems with multiple frictional contacts. This study proposes a general physics-informed neural network framework, termed the mixed leading variable physics-informed neural network (MLV-PINN). By using the governing equations and contact-friction laws as physical constraints, this framework provides common solution approaches for computing nonsmooth dynamics not only in discrete systems but also in continuum systems. Specifically, the network represents accelerations, normal contact forces and friction margins as independent outputs and determines them simultaneously, without first reducing the dynamics to discrete velocity-impulse complementary subproblems. A composite physics-informed loss is constructed from the equations of motion in their original form, contact-friction complementary conditions and a friction-margin consistency relation. To improve stick-slip state identification during PINN optimisation, a mixed leading variable combining relative tangential velocity and acceleration is introduced to reformulate the tangential complementary conditions. Four benchmark simulations cover one-dimensional (1-D) nonsmooth dynamics involving single-contact stick-slip, 1-D nonsmooth dynamics involving multi-contact stick-slip, two-dimensional (2-D) nonsmooth dynamics with coupled normal contact-separation and tangential stick-slip, and 2-D nonsmooth dynamics in a continuum system involving moving contact on an elastic beam. Benchmark comparisons show substantial reductions in displacement and velocity errors relative to the corresponding conventional methods at large time steps. These results, together with an ablation study and tests involving friction-velocity slopes, long-time integration and loss weights, illustrate the computational accuracy and support the broad applicability of the proposed framework.
For deeper understanding of quantum non-local correlations a new path integral real-time quantum dynamics approach has been developed and applied to the Monte Carlo simulation of the decoherence of pair entangled states in many-body systems in the framework of the ``replica trick''. Calculations have been carried out for the Heisenberg-Wigner path integral representation of the entangled state correlation function (ESCF) in the modified unified electron gas model as example. The Wiener-Khinchin transformation of ESCF has been applied to random paths in path integral allowing for the Monte Carlo simulation of ESCF. Presented results include the spin-resolved radial distribution functions, ESCFs and estimations of the decoherence time of the entangled electron pairs due to thermalization and interparticle interaction.
Bone scaffolds must reconcile a macroscale mechanical response with microscale geometric feasibility, yet conventional approaches to their design restrict attainable constitutive properties. This work presents a concurrent two-scale topology optimization framework; a conditional rectified flow generates periodic lattice microstructures, which are then refined in latent space by exact finite-element homogenization sensitivities. The resulting unit-cell geometries and their spatial assignments are jointly optimized under mechanical and feature-size constraints. In a distal rabbit femur benchmark, all eight flow-started refinements produce printable cells with maximum relative errors under 1\% across six prescribed engineering moduli. Constructing the implant from a single tiled cell reproduces intact whole-femur compliance within 0.1\% and regional strain energy within 0.35\% under four loading conditions. In human femoral reconstruction, multiscale concurrent optimization produces a cortical mid-shaft scaffold with porosity 63\% while satisfying the compliance constraint within 2\% under walking and stair-climbing loads. It achieves an 84\% porosity on a proximal cancellous scaffold with strain-energy deviations below 0.25\% under the same loads. Analysis of the learned representations reveals distinct architectural families in latent space. Further, the attention heads of the rectified-flow transformer specialize for spatial shifts -- learning the analogue of a finite-difference stencil -- and coordinate symmetries, with ablations showing their influence on generated topologies. These results establish a route to printable, mechanically constrained scaffold designs. The framework is general enough for extending to many other multiscale structural design applications including topology-microstructure co-design in metal additive manufacturing.
We develop a Unified Interface Neural Network (UINN) framework for elliptic interface problems with discontinuous coefficients. A single shared network provides a unified representation of the solution in both subdomains, while a projection--correction construction embeds the solution and flux jump conditions directly into the neural trial space. This design aligns with the classical numerical principle of incorporating interfacial physics into the approximation basis to construct interface-unfitted numerical methods. The jump conditions are therefore satisfied exactly for any network parameters. Compared with methods based on separate phase networks and numerical formulations with interface penalty terms, UINN provides a unified representation with exact enforcement of interfacial jump conditions and interface-penalty-free optimization. Those features streamline the numerical formulations and can enhance training robustness. Numerical experiments on smooth interfaces at modestly large and high coefficient contrasts demonstrate UINN's solution accuracy, near-machine-precision satisfaction of the interface jump conditions, and enhanced robustness.
We present and validate an Implicit Monte Carlo (IMC) radiation-transport solver in the open-source, mesh-agnostic storm library, coupled to the moving-mesh hydrodynamics code rich in three-dimensional Voronoi tessellations. The solver is multigroup, or continuous in frequency, throughout; it transforms photon packets exactly between the comoving and laboratory frames, books every exchange with the material as a laboratory four-momentum so that the IMC path conserves energy to roundoff, formulates the Fleck factor in the fluid's proper time, and accelerates optically thick cells with multigroup Discrete Diffusion Monte Carlo (DDMC) behind an emissivity-preserving interface. We tested the solver on a broad suite of benchmark problems, and it performed well on all of them. The solver is built on the MPI-parallel Monte Carlo framework of storm, which scales to 13,440 cores.
We propose generative marginal learning (GML), a framework for equilibrium sampling and local observable estimation in molecular simulations. By training a conditional diffusion model on large-system simulations, GML learns the marginal distribution of a retained region, enabling configuration generation without explicitly simulating the surrounding environment or imposing periodic boundary conditions on the retained region. We derive a marginal effective potential and the associated conditional mean-force identity, which provide the basis for both denoising and force-alignment training. As an application, we estimate the spatially averaged Irving--Kirkwood stress near a vacancy defect. Numerical results demonstrate that GML substantially improves pair structure and local stress accuracy compared with small periodic molecular dynamics simulations, mitigating finite-size effects in calculating local molecular observables.
Depletion analyses are often driven by Monte Carlo neutron transport methods that introduce uncertainty into the transmutation matrix used to evolve a material composition. The nonlinearity of the matrix exponential applied to this transmutation matrix causes the uncertainty to introduce a systematic bias in the computed isotope vectors. This work derives a second-order correction for this bias that can be computed with quantities obtainable from typical Monte Carlo transport routines. An augmented Bateman system is formulated to evaluate these quantities and an adapted Chebyshev rational approximation method uses transmutation matrix factorizations for efficient evaluation of the second-order correction. Then, the correction is demonstrated using a ten-nuclide verification problem, an ITER port-plug activation problem, and a 9-pin boiling water reactor lattice containing self-shielded gadolinium-bearing fuel.
Density Functional Theory (DFT) is the established workhorse of computational materials science, yet it faces fundamental challenges with strongly correlated electrons, such as in transition metal complexes and bond-breaking processes. While quantum computing offers a theoretically exact path to resolve these problems, current hardware limitations prohibit fully quantum simulations of realistic material systems. To bridge this gap, we present SIESTA-QCOMP, a modular suite embedding quantum methodologies within DFT. Harnessing SIESTA's localised basis set, we implement a projection-based embedding scheme that treats the bulk environment classically while mapping the strongly correlated "impurity" subspace to a quantum processor. The architecture uses a socket-based driver to orchestrate workflows between the Fortran-based SIESTA engine and Python-based frameworks like Qiskit, ensuring high interoperability. We detail the hybrid workflow, currently utilising a Variational Quantum Eigensolver (VQE), and discuss the road map toward methods such as Sample-based Quantum Diagonalisation (SQD). We illustrate the framework on an hypothetical iron porphyrin calculation and argue that such hybrid architectures are not merely interim solutions but essential foundations for integrating quantum utility into materials design, as envisioned by our broader QCOMP4DFT initiative.
Low-order Carleman embeddings replace nonlinear dynamics by a larger linear system, but truncation errors can accumulate even when the underlying nonlinear model is accurately approximated by a quadratic equation. We construct a residual Carleman lift (RC2) of lattice Boltzmann dynamics in which an independently evolved spectral background supplies the coefficients and forcing of a second-order residual lift. In a double-shear-layer test, conventional fixed-rest Carleman reaches $87.3\%$ final error, whereas prepared RC2 achieves $0.136\%$ through vortex roll-up, without restarting the lift or using fine-reference feedback. The results demonstrate control of the finite-time accuracy of a fixed-order linear embedding through nonlinear preparation.
Implicit discretizations of transient heat conduction repeatedly generate sparse linear systems, making diffusion problems a natural test bed for quantum linear solvers. This study evaluates how accurately a variational quantum linear solver (VQLS) recovers the temperature field produced by one backward-Euler step of one-dimensional conduction on IBM quantum hardware. A pre-specified confirmatory experiment compares N=4 interior temperatures encoded on two qubits with N=8 temperatures encoded on three qubits using ten paired hardware blocks; a separately defined N=16, four-qubit extension is evaluated in five independently scheduled exploratory blocks. Local VQLS optimization reduces relative temperature-field error to approximately 10^{-8} for all three cases, establishing an essentially exact ideal reference at this scale. On hardware, mean relative L2 temperature-field error is 8.07% for N=4 and 20.14% for N=8, a paired increase of 12.06 percentage points (95% CI, 10.71-13.42; exact paired sign-flip p=0.001953). The exploratory N=16 mean is 24.75+/-1.84%. A pre-specified N=16 shot-count follow-up then evaluates the same locked circuit at 2,048, 8,192, and 32,768 shots in six independently scheduled blocks. Mean field error remains 23.64%, 24.57%, and 24.37%, respectively; the 32,768-minus-2,048-shot contrast is +0.74 percentage points (95% CI, -0.65 to 2.12; p=0.281). In contrast, ideal multinomial sampling predicts a reduction from 2.68% to 0.75%. Higher shot counts reduce block-to-block variability but do not recover the mean cold-side temperature field. The results show that full-field readout of a strongly nonuniform thermal state is limited by more than ordinary finite-shot sampling under the tested hardware conditions, motivating observable-focused quantum formulations for heat flux, stored energy, and selected temperatures.
Foundation machine-learned interatomic potentials (MLIPs) promise broad transferability, yet their out-of-domain reliability is hard to establish. Here, we use pure and O/H-containing disordered carbon spanning multiple phases, competing hybridisations and porous frameworks to stress-test six popular foundation MLIPs. We first screen all six on density functional theory (DFT) energy and force errors and binding curves, then subject four to melt-quench-anneal molecular dynamics and evaluate the resulting structural observables. Binding curves and O/H-containing trajectories reveal unphysical extrapolation and dynamical failures that low DFT energy and force errors fail to reveal. Of the models tested, MACE-MH-1 is the most consistently robust, with low DFT errors and stable MD. Despite differing architectures, OMAT-based models share structural biases, highlighting the importance of training coverage. All four models generate structures of nanoporous carbon that reproduce the principal experimental and specialised carbon potential signatures. Incorporating O/H leaves short-range order largely unchanged but shifts pore-size distributions towards larger pores. These results emphasise the need for application-specific validation and extend molecular dynamics modelling of nanoporous carbon towards experimentally realistic structures.
Discrete adjoint methods give the sensitivities that aerodynamic design and data assimilation require, but they remain hard to obtain for flow solvers that run on graphics processors. This paper presents CuAD, a source-to-source automatic differentiation tool that reads CUDA kernels and writes their tangent and adjoint versions as ordinary CUDA source. Its reverse mode is organised around the GPU thread and exploits the structure of flow-solver kernels: pure accumulations need no tape, the remaining intermediate values are recomputed rather than stored, and adjoint updates are atomic only where several threads reach the same memory. We evaluate CuAD on four kernels from simple flow models: a viscous Burgers residual, a two-dimensional advection-diffusion residual, a quasi-one-dimensional Euler residual and a Biot-Savart summation over point vortices. All generated derivatives agree with finite differences and pass adjoint dot-product tests to round-off, while the adjoint kernels cost two to six times the primal kernels. Compared with the compiler-based tool Clad under the same compiler, CuAD gives identical results at identical cost on the stencil kernels and is 16 to 470 times faster on the summation kernel because it stores no tape. A steady Burgers study shows that in single precision the adjoint gradient of a finite-difference Jacobian reaches errors of 45 to 75 per cent once the condition number grows to $2.5\times10^{5}$, whereas the CuAD Jacobian keeps the error below $5\times10^{-5}$. Four applications complete the evaluation: adjoint fields of the Burgers shock, the sensitivity of a scalar sensor in a cellular vortex, the inverse design of a transonic nozzle whose area is recovered to $2\times10^{-8}$, and an unsteady adjoint of the roll-up of a vortex sheet over 400 time steps.
A limited understanding of how evaporation couples to crystallization pathways hampers control over the formation of solution-processed perovskite photoactive layers. Identifying these mechanisms-and learning how to select and tune them-is essential for achieving the desired crystal morphology and, ultimately, improving solar-cell performance and stability. To investigate crystallization pathways in evaporating thin films, we develop a simple and general zero-dimensional population-balance model that couples evaporation, chemical reactions, and crystallization. We use the model to explore a range of parameters for a solvent-solute system that follows a two-step crystallization pathway: solvent evaporation first induces formation of a crystalline precursor, which subsequently recrystallizes into the final phase. The model predicts the evolution of the film composition during drying and provides quantitative LaMer curves. Based on these results, we introduce a vocabulary of "supply limitations" and "rate limitations" that enables the interpretation of crystallization kinetics in terms of the governing process conditions.
Solutions of second-order elliptic PDE problems on nonconvex or cracked domains with mixed Dirichlet-Neumann conditions typically lie outside $H^2(Ω)$. Near re-entrant corners, crack tips and points where the boundary condition changes type, the gradient grows like $r^{λ-1}$ with $λ<1$. We study residual-based adaptive refinement in hybrid variational-collocation PINNs. A multiscale Fourier-feature network is trained on a composite loss. The loss combines a variational residual tested against Legendre bubble functions, a strong-form collocation term and penalties for the Dirichlet and Neumann data. The variational term needs only first derivatives of the network, whereas the collocation term requires second derivatives and provides a residual that can be evaluated at any point. After each training stage, candidate points selected by this residual are added to the collocation set. The method is tested on five two dimensional benchmarks, namely a Mode~III crack, a cross-shaped domain with four re-entrant corners with unequal amplitude, a reaction-diffusion problem with four localized peaks, a monotone semilinear problem with a boundary layer and the Motz problem. Compared with uniform point insertion under the same collocation and iteration budget, residual-based refinement reduced the relative $L^2$ error by up to a factor of $32$, the relative $L^\infty$ error by up to a factor of $38$ and the relative $H^1$ semi-norm error by up to a factor of $15$.
We present NEP-kappa 2.0, a unified framework that integrates density-functional theory, material-specific machine-learning interatomic potentials, and pretrained force models for advanced phonon calculations. A common workflow enables harmonic lattice dynamics, quasi-harmonic and stochastic self-consistent phonon treatments, three- and four-phonon transport, Wigner population--coherence transport, and time-dependent phonon dynamics. Applications to 3C-SiC, BAs, and Cs$_2$AgInCl$_6$ demonstrate consistent force-model comparisons, finite-temperature renormalization, higher-order scattering, and intermode coherence. Beyond steady-state transport, mode-resolved simulations of selectively excited 3C-SiC reveal that half of the injected optical-phonon energy reaches acoustic modes within approximately 3.75 ps, primarily through one-optical--two-acoustic scattering. Acoustic branches exhibit response times differing by a factor of five, with competing scattering channels producing a transient energy overshoot. Channel-resolved entropy production and persistent nonthermal occupations distinguish interbranch energy redistribution from intrabranch thermalization. The framework connects equilibrium phonon transport with microscopic nonequilibrium energy dynamics.
Tungsten-tantalum alloys are promising plasma-facing materials for fusion reactors, but simulations of their coupled defect, mechanical, and microstructural evolution require an interatomic potential that remains accurate across diverse atomic environments. We develop a transferable W-Ta potential within the GRaph Atomic Cluster Expansion (GRACE) framework. The primary training dataset was generated automatically by maximizing information entropy in an atomically relevant feature space, avoiding manual enumeration of the structures and processes required for specific downstream applications. Extensive validation demonstrates a consistent description of unary and binary phase stability, elasticity, phonons, high-pressure response, surfaces, point defects, dislocation energetics, transformation pathways, alloy thermodynamics, and melting behavior. As a demonstration, the potential is applied to Ta segregation and shear-coupled migration of a $\Sigma5(120)$ W GB. Simulations reveal strongly site-selective segregation: Ta preferentially occupies a small subset of crystallographically distinct atomic rows. Under shear, segregated Ta impedes disconnection-mediated grain boundary migration. The de-pinning stress initially increases with Ta concentration but plateaus as the dominant pinning environment becomes saturated, while increasing temperature weakens both segregation and pinning. These results establish a broadly validated W-Ta potential and demonstrate how automated entropy-maximized data generation can produce transferable atomistic models capable of resolving mechanisms not explicitly targeted during training.
Terahertz scanning near-field optical microscopy (THz-SNOM) enables spatially and spectrally resolved investigation of low-energy collective and topological states in quantum materials and nano-structured devices by exploiting the strongly localized near field generated by a polarized atomic force microscopy (AFM) tip. Because the THz near-field response is highly sensitive to tip geometry, understanding how geometric parameters control field confinement and decay is essential for selecting tip designs that provide the spatial resolution required for specific experiments. In this work, we develop a boundary element method (BEM) for modeling the electrostatic near-field of a conducting AFM tip in response to an external background field. The tip is represented by a simplified closed, blunt-cone geometry with its surface discretized into panels using increased mesh density near the apex. A modified Green's function is used to calculate interactions between surface panels, while an augmented linear system simultaneously enforces the equipotential boundary condition on the conducting surface and overall charge neutrality to determine the induced surface-charge distribution. We validate the BEM by comparing the resulting near-field decay with that obtained from a three-dimensional quasi-static COMSOL Multiphysics simulation. The BEM accurately reproduces the COMSOL field-decay behavior while requiring orders of magnitude fewer mesh elements and substantially reducing computational cost. These results demonstrate that the BEM provides an accurate and computationally efficient framework for systematically investigating the dependence of THz near-field confinement on AFM tip geometry. This framework can be extended to more complex and experimentally realistic tip geometries, providing a route toward rational tip selection and design for application-specific spatial-resolution requirements.
QRecoupling.jl is a Julia package for numerical, exact and symbolic evaluation of classical and quantum recoupling coefficients and generic finite sums of quantum factorials, with one interface across the classical limit, principal unitary levels and real or complex deformation. A compact factorial rule and its cyclotomic exponents expose factor cancellation, polynomial structure and level support. The package covers quantum Clebsch-Gordan coefficients, 3j and 6j symbols and user-defined sums. Its adaptive numerical evaluation is combined with exact decisions for sums that cancel to zero. Classical 3j and 6j symbols are evaluated at a speed comparable to that of WIGXJPF, which is based on prime factorization tables, at the smallest spins, and about 6 and 10 times faster at spin cutoff 500. Under real, root-of-unity and complex deformation, the largest relative error measured against independent high-precision references is $3.5 \times 10^{-14}$ for spins up to 400. Fusion matrices with 401 channels are computed about 100 times faster than by loops over their entries.
Physics-informed neural networks (PINNs) are trained to minimize the residual of a prescribed equation, so a low training loss does not indicate whether a coefficient in that equation was specified incorrectly. This study evaluates what can be inferred about the encoded parameter from the trained network alone. When the residual is linear in the parameter, written as $A - νB$, the least-squares estimate $\hatν = A^{\top}B/B^{\top}B$ is obtained from a single evaluation of the frozen network, without retraining or reference data. Its accuracy depends on how closely the network satisfies its own equation, quantified by a normalized residual $ρ$. Recovery error is associated with $ρ$ (log-log correlation 0.72 across 36 convection-diffusion networks), and empirical exponents of 0.63 to 1.20 across three equations decompose into one half plus a measurable angular term. For Burgers' equation, a single-model test based on a block-bootstrap standard error flagged 3 of 300 correctly specified networks (observed false-alarm rate 1%), all 2% misspecifications at a practical budget, and two thirds of 1% misspecifications after convergence. For the self-similar blowup rate of the Córdoba-Córdoba-Fontelos equation, where the rate also enters the architecture, the estimator recovers about half of the misspecification, whereas $ρ$ increases approximately quadratically and distinguishes a 0.2% error. When the incorrect rate enters both the equation and the architecture, errors of up to 2% were not detected at the precision reached. At practical budgets, the raw loss can be lower for the misspecified model, a scaling effect that $ρ$ does not share.