2026-07-30 | | Total: 11
Oscillator Ising machines (OIMs) and dynamical Ising machines (DIMs) encode binary spins in phase states stabilized by second-harmonic injection (SHI). In a coupled network, the competition between SHI and the instantaneous local network field reshapes each oscillator's conditional energy landscape. We show that this competition drives a transition between monostable and bistable regimes through an astroid spinodal boundary. Near this boundary, the barrier scales as $ΔE_i\proptoμ_i^{3/2}$ at generic smooth points and as $ΔE_i\proptoμ_i^{2}$ at the longitudinal cusp. OIMs and DIMs obey the same spinodal geometry, with their conditional landscapes related by a reversal of the transverse field. Finally, the first-harmonic conditional landscape is mathematically equivalent, up to an additive constant, to the Stoner--Wohlfarth energy of a uniaxial magnetic particle.
We present SPECTRE, a new equilibrium code based on the Multi-Region relaxed MHD model for robustly calculating 3D equilibria with general magnetic topology, allowing for flux surfaces, magnetic islands, and chaos. The code builds on a previous MRxMHD solver, SPEC, but performs significantly better thanks to a new formulation of force, the use of a stable trust-region-based least squares minimization scheme, as well as several additional features. SPECTRE is verified through application to configurations with known equilibrium solutions, in vacuum and with finite beta, both in the fixed boundary and the free boundary mode. Notably, these include vacuum equilibria of a quasi-axisymmetric (QA) device in fixed-boundary mode, and of W7-X in the free-boundary mode, along with a classical stellarator finite-beta free-boundary case. Finally, the solver is applied to a modern optimized finite-beta quasi-isodynamic (QI) configuration, where we demonstrate calculation of a strongly-shaped equilibrium with a core island which is in agreement with a HINT calculation.
Precipitation-strengthened Cu-Ni-Al alloys are of interest for technological applications because coherent, $\mathrm{L}1_2$-ordered $γ'$-$\mathrm{Ni}_3\mathrm{Al}$ precipitates can confer high mechanical strength while allowing the material to retain many of the good transport properties characteristic of elemental Cu. In this work, we study the thermodynamics and phase stability of the pseudobinary $\textrm{Cu}_x (\textrm{Ni}_{3/4} \textrm{Al}_{1/4})_{1-x}$ system, $0 \leq x \leq 1$. We use a computational modelling framework combining first-principles electronic structure calculations with a concentration wave analysis from which atom-atom effective pair interactions are extracted for use in atomistic Monte Carlo simulations. Our modelling reveals three distinct, composition-dependent regimes of phase behaviour, in qualitative agreement with the experimentally determined phase diagram. At low Cu content, Cu is soluble in the $\mathrm{L}1_2$-ordered $\mathrm{Ni}_3\mathrm{Al}$ phase, with a single identifiable phase transition corresponding to chemical ordering between Ni and Al. At intermediate compositions, this high-temperature ordering is followed at lower temperatures by phase separation of Cu and $\mathrm{L}1_2$-ordered $\mathrm{Ni}_3\mathrm{Al}$. Finally, at high Cu content, $\mathrm{L}1_2$-ordered $\mathrm{Ni}_3\mathrm{Al}$ precipitates directly from the solid solution, with no clearly identifiable secondary transition. We relate these phase transformations to features of the underlying electronic structures of the considered alloys. Overall, this work demonstrates a computationally efficient workflow capturing both chemical ordering and coherent precipitation in multicomponent substitutional alloys, with relevance to the study of phenomena such as precipitation strengthening.
Spin-dependent transport in ferromagnets underpins the development of high-density spintronic memories. Spin-dependent transport in strong spin-orbit-coupled ferromagnets exhibits a significant anisotropy. Both the overall spin polarization during charge transport and the magnetic spin Hall conductivity are found to exhibit pronounced anisotropy when the magnetization is tilted away from the crystallographic easy axis or when the electric field is rotated relative to the crystal axes. These anisotropic responses originate primarily from spin-orbit coupling, which is identified as the key driver of the large anisotropy observed in ferromagnet. Furthermore, strain tunability of the magnetic spin Hall anisotropy is demonstrated, with tensile strain progressively enhancing the oscillatory amplitude of the spin Hall conductivity. These findings establish strong spin-orbit-coupled ferromagnets as a platform for anisotropic spin-current generation and field-free spintronic devices that exploit intrinsic material anisotropy for improved performance and energy efficiency.
Intermittency in turbulence manifests as intense vortices and sharp peaks of dissipation. Causing the breakdown of Kolmogorov's simple self-similar theory, it leads to anomalous scaling, multifractality and so far remains beyond the scope of a complete theoretical description. How intermittent flow structures influence these different measurements is not known quantitatively. With a simple filtering procedure-thresholding vorticity and inverting the Biot-Savart law to generate filtered velocity fields-we show the effects of intermittent flow structures can be disentangled. As extreme vorticity contributions to the velocity field are filtered out, the energy spectrum scaling persists, while the bottleneck is flattened, and structure function scalings tend towards their Kolmogorov values. The approach is more rapid for transverse exponents, revealing the selective importance of intensely swirling flow regions. Similarly, the extent of multifractality reduces as intermittency is filtered, shrinking the range of roughness singularity exponents. The residual fields are curiously more multifractal, but their structure begins to break away from an underlying turbulence skeleton. The effects on vortex stretching and strain self-amplification are quantified. Our work shows that a Biot-Savart approach can selectively remove the effects of intermittency from turbulence, and hence from its scalings.
Data-driven subgrid-scale closures for large eddy simulation are of significant interest in many engineering and geoscience applications. In this context, several important questions remain about the role of rotational equivariance as an inductive bias for learned tensorial mappings. We investigate whether equivariance improves accuracy, parameter efficiency, and generalization for subgrid-scale modelling at realistic filter ratios. For turbulent channel flow, we compare data-augmented non-equivariant architectures to those with equivariance as an inductive bias. We compare both pointwise and nonlocal versions of these two model classes. All models are evaluated at matched parameter counts across spatiotemporal, anisotropy, and Reynolds number generalization. We show that non-augmented models learn a small degree of equivariance directly from turbulence data, especially when that data is more isotropic. The equivariant nonlocal architecture attains the highest correlation coefficient on every generalization test at approximately half the parameter count of its non-equivariant counterpart, while the pointwise architectures do not improve on the analytical Clark baseline. Additionally, the equivariant model is more data-efficient than a non-equivariant model. The benefit of equivariance grows with the receptive field of the model, indicating that equivariance and nonlocality are both useful for the subgrid-scale closure task at realistic dataset size, parameter counts, and filter size.
Magnetic material parameters such as the exchange stiffness and magnetic anisotropy govern the behavior and functionality of magnetic systems, yet their local inference from magnetization data remains challenging, particularly in strongly fluctuating regimes with polycrystalline or multiphase microstructure, where conventional texture-based methods become unreliable. We introduce a magnetization-only framework for inferring material parameters from thermally driven magnetization dynamics. Using micromagnetic simulations, we extract statistical quantities such as temporal mean and latent entropy from the magnetization dynamics, fit models to these descriptors, and invert the models to infer material parameters. We show that this framework enables material-parameter inference as well as grain-boundary detection in a heterogeneous sample. Among the descriptors considered, latent entropy yields more accurate parameter estimates than the temporal mean. Our results establish latent entropy as an efficient descriptor for inferring magnetic material parameters from dynamical magnetization data and point toward its use for experimental parameter extraction at high temperatures and, more broadly, under strongly fluctuating conditions.
This study considers extracting causally important vortical structures from the extreme vortex gust-airfoil interaction at a chord-based Reynolds number of $5000$. This extraction is achieved by decomposing a given vortical flow snapshot into its informative and residual components based on the contribution to an arbitrary future target variable with convolutional information-theoretic learning. For the current vortex-airfoil interactions that exhibit transient and multiscale flow characteristics, we first examine the important vortical structures with respect to a future lift coefficient. While the vortex cores are primarily highlighted before vortex impingement, the emerging shear layers are additionally captured after the massive separation, which is evident from a comparison to an instantaneous force-element analysis. We further take the scale-dependent energy transfer as a future variable of interest to examine its impact on the extracted informative structures compared to the lift-associated structures. They are distinct from the lift-based structures in the early stage of the gust encounter yet become similar after impingement, revealing an analogy between informative structures across different transient aerodynamic mechanisms. The present data-driven approach selectively extracts the specific important flow structures responsible for the physics of interest, which can support studying a range of transient aerodynamic flows from the causal, data-driven perspective.
Giant piezoelectricity in ferroelectrics is commonly associated with phase-boundary instabilities, among which the polymorphic phase boundary (PPB) is a prominent example conventionally attributed to the coexistence of ferroelectric phases. Here, using large-scale molecular dynamics simulations of the lead-free (K,Na)NbO3-(Bi,Na)ZrO3 solid solutions, we show that the PPB hosts a percolating multifractal polar domain which governs the dielectric and piezoelectric responses. By quantifying the global fractal dimension and multifractal spectrum width of this polar network, we identify fractal connectivity and multiscale heterogeneity as microstructural order parameters for the PPB. The maximum reversible piezoelectric response occurs when the fractal-domain volume fraction approaches the three-dimensional percolation threshold, suggesting that near-critical polar connectivity enables giant reversible electromechanical coupling. In this mechanism, the fractal backbone preserves polar memory and provides the restoring force required for reversibility, while the surrounding nonfractal regions supply the polar compliance needed for large polarization rotation and strain. These results establish percolating multifractal polar domains as a microscopic mechanism for PPB-enhanced piezoelectricity and suggest fractal connectivity as a design parameter for high-performance piezoelectrics.
The imposition of boundary velocities in finite element models of port-Hamiltonian elastodynamics typically relies on Lagrange multipliers, yielding Differential-Algebraic Equations (DAEs). Alternatively, weak imposition methods that maintain an Ordinary Differential Equation (ODE) structure often exhibit poor accuracy at Dirichlet boundaries. To address these limitations, this paper introduces an additive kinematic decomposition at the continuous level, splitting the displacement and velocity fields into a relative dynamic component that vanishes on the boundary and a prescribed lifting function extending into the interior domain. This decomposition induces a distributed port that maps the effects of the boundary actuation inside the domain. By incorporating this mapping into suitable virtual power principles, we derive lifted port-Hamiltonian system (PHS) models that, upon finite element discretization, reduce to ODE systems in which Dirichlet boundary velocities are strongly imposed. The framework is applied to derive 2-field and 4-field formulations suited to distinct PHS geometric representations. Furthermore, we show that under specific shape functions, standard FEM schemes are recovered, demonstrating that the lifting framework in the discrete models is equivalent to the classic algebraic matrix partitioning in computational mechanics practice. The energy-balance properties and computational performance of the proposed methodology are verified through numerical simulations.
The interaction of nanosecond laser pulses with metallic materials involves multiple complex physical processes. It is challenging to construct a self-consistent model capable of uniformly describing all stages. This work establishes a multi-physics coupling model for pure iron, encompassing laser energy deposition, solid-liquid phase transition, gas-liquid interfacial kinetic transport, plasma expansion and ionization, and spectral radiation. The numerical solution adopts a partition method, utilizing an implicit compact difference scheme for the target and a Mac-Cormack explicit scheme for the plasma. The simulations elucidate the emergence of plasma shielding and its inhibitory effect on the evaporation process, thereby confirming that 81.6% of the early-stage ablation products are transported through a supersonic expansion mode. The model successfully captures the complete evolution of the plasma plume from a high-temperature, highly ionized state to a low-temperature, neutral atomic state. Based on this, spectral calculations demonstrate the dynamic evolution of radiative characteristics from an early stage featuring a strong continuum background dominated by ion lines to a later stage where the continuum attenuates, atomic lines become prominent, and self-absorption appears. The emergence of self-absorption proves the ability of the model to effectively capture the optical thickness effects arising from spatial inhomogeneity within the plasma. Through systematic comparison between experimentally measured spectra and calculated results from the PrismSPECT and NIST LIBS spectral programs, the model presented here achieves the highest comprehensive scores in quantitative evaluations of multiple channels. This validates the necessity and superiority of the full-chain self-consistent modeling approach, especially in describing plasma inhomogeneity and radiation transport.