Papers

Many studies of non-relativistic spin-splitting (NRSS), or altermagnetism, have focused on idealized, perfectly ordered crystals, relying on symmetry-based approaches to identify candidate materials. Here, we theoretically investigate how local short-range ordering (SRO) influences NRSS of energy bands in partially ordered collinear antiferromagnetic iron oxyfluoride (FeOF). Using the cluster expansion method, we identify four nearly degenerate structures (energy difference $\leq 8$ meV per formula unit) that represent distinct snapshots of local plane-to-plane O/F correlations. Our density functional theory (DFT) results show robust NRSS along the $\Gamma$-M direction in all four structures, despite the absence of long-range order. The magnitude and character of the splitting depend sensitively on the specific direction of anion correlations, effects that are not fully captured in high-symmetry average structures. Notably, two configurations ($Pmc2_1$ and $Pm$) exhibit $\Gamma$-point spin splitting absent in ordered FeF$_2$ and a virtual crystal approximation model of FeOF. We further predict distinct magneto-optical Kerr effect (MOKE) signatures, enabling experimental detection of SRO-driven electronic structure changes. These results highlight heteroanionic compounds as a promising design space for NRSS antiferromagnets, with experimentally synthesized FeOF already exhibiting a substantially higher Néel temperature (315\,K) than FeF$_2$ (79\,K).

Spin-transfer torque magnetic random-access memory (STT-MRAM) relies on nanoscale magnetic tunnel junctions (MTJs) as its fundamental building blocks. Next-generation STT-MRAM requires strategies that simultaneously improve switching energy efficiency and device endurance. Here, we present the first study of perpendicular STT-MRAM writing assisted by radio-frequency (RF) spin torque. We show that applying a small-amplitude RF pulse prior to a direct-current (DC) writing pulse enhances the MTJ switching probability, with the efficiency gain increasing at lower RF frequencies. This RF+DC writing scheme enables shorter DC pulses, thereby improving device endurance. Analytical and numerical modeling qualitatively reproduces the experimental trends, while quantitative discrepancies indicate that realistic MTJ properties beyond idealized models play an important role in RF-assisted switching.

Accelerating the discovery of structural materials is essential for applications in hard and refractory alloys, hypersonic platforms, nuclear systems, and other extreme environment technologies. Progress is often constrained by slow synthesis and characterization cycles and the need for extensive mechanical testing across large compositional spaces. Here, we propose a rapid screening strategy based on random material libraries, in which thousands of distinct compositions are embedded within a single specimen, mapped by EDS, and subsequently characterized. Using nanoindentation as a representative case, we show that such libraries enable dense composition property mapping while reducing the number of samples required to explore high dimensional composition spaces compared to traditional synthesis and test workflows. An experimentally calibrated Monte Carlo framework is developed to quantify practical limits, including particle size, EDS noise and resolution, positional accuracy, and nanoindenter motion costs. The simulations identify regimes where random libraries provide orders of magnitude acceleration over classical workflows. Finally, we demonstrate experimental navigation of these libraries using automated indentation. Together, these results establish random libraries as a general route to high throughput characterization in structurally critical material systems.

As an ultrawide bandgap semiconductor, beta-Ga2O3 has been attractive for its strong tolerance to irradiation damage and high n-type conductivity through ion implantation. Homoepitaxial (010) e̱ta}-Ga2O3 films grown by MOCVD were implanted with Ge to study the post-implantation damage and lattice recovery after thermal annealing. Box profiles of 100 or 50 nm at concentration of 5*10^19 or 3*10^19 cm^-3 were formed, with maximum displacement per atom (DPA) of 1.2 or 2.0. Lattice recovery was investigated using X-ray diffraction (XRD) for anneals from 100 C to 1050 C. A gamma-phase related peak was observed for all implant conditions. All samples showed strain relaxation of beta-phase peak at temperature below 500 C, with no significant change for the gamma-phase related peak. For lower damage implants, films recovered fully to epitaxial beta-phase after sequential annealing to 900 C. For the higher damage implant, the gamma-phase associated peak annealed out with increasing temperature, but a new diffraction peak formed at slightly smaller lattice spacing; full recovery of the lattice was not observed until annealing at 1050 C. The newly formed diffraction peak is identified as beta-(20-4), beta-(512), or beta-(71-2), each potentially arising from the conversion of gamma-phase to beta-phase via a common oxygen sub-lattice.

In trapped Bose-Einstein condensates, interaction quenches which are abrupt changes of the interaction strength typically implemented via Feshbach tuning, are a practical and widely used protocol to address far-from-equilibrium collective modes. Using both numerical Gross Pitaevskii and analytical schemes we study these interaction-quench-induced collective modes in a harmonically trapped two-dimensional Bose–Einstein condensate contrasting the behavior found at low and high energies. In the low-lying regime, we characterize realistic circumstances in which there is a breakdown of the expected scale invariance so that the collective excitations follow hydrodynamic theory instead of the predictions given by SO(2,1) conformal symmetry. In the high energy regime, we focus on important trap effects associated with acoustic oscillations which have been of interest experimentally. This comprehensive analysis of the collective excitations in trapped two-dimensional Bose-Einstein condensates is experimentally accessible. Through their frequencies and damping, this reflects an important built-in spectroscopy of such many-body states.

For long-range percolation on $\mathbb{Z}$ with translation-invariant edge kernel $J$, it is a classical theorem of Aizenman and Newman (1986) that the phase transition is discontinuous when $J(x,y)$ is of order $|x-y|^{-2}$ and that there is no phase transition at all when $J(x,y)=o(|x-y|^{-2})$. We prove a strengthened version of this theorem for the hierarchical lattice, where the relevant threshold is at $|x-y|^{-2d} \log\log |x-y|$ rather than $|x-y|^{-2}$: There is a continuous phase transition for kernels of larger order, a discontinuous phase transition for kernels of exactly this order, and no phase transition at all for kernels of smaller order. As such, $|x-y|^{-2d} \log\log |x-y|$ is essentially the only kernel that produces a discontinuous phase transition. We also prove a hierarchical analogue of the “$M^2\beta=1$” conjecture of Imbrie and Newman (1988), which gives an exact formula for the density of the infinite cluster at the point of discontinuous phase transition and remains open in the Euclidean setting.

This manuscript reports a semiconductor-sensitized thermal cell (STC) that converts ambient heat into electrical power while simultaneously reducing its own temperature under isothermal conditions. Using a printable semiconductor–electrolyte architecture, we fabricate $4\,\mathrm{cm} \times 4\,\mathrm{cm}$ devices that generate up to approximately $0.2\,\mathrm{mW}$ at temperatures of $40$–$55\,^\circ\mathrm{C}$. During continuous discharge, the STC exhibits a transient temperature decrease followed by thermal equilibration with the environment. In contrast, periodic on–off discharge produces sustained cooling of approximately $1\,^\circ\mathrm{C}$ relative to a non-discharging reference. Notably, parallel integration of four STCs yields a nonlinear enhancement of cooling (approximately $5\,^\circ\mathrm{C}$) without a corresponding increase in electrical output. The observed behavior can be understood within a macroscopic energy-balance framework, in which time modulation of electrochemical heat consumption prevents the establishment of thermal steady state. These results demonstrate sustained isothermal cooling induced by heat-to-electricity conversion at practical device scales, and highlight semiconductor-sensitized thermal cells as a platform for coupled energy harvesting and thermal management.

We investigate finite-temperature magnetic and topological phase transitions in two-dimensional honeycomb ferromagnets using an extended self-consistent renormalized spin-wave theory (SRSWT) that incorporates higher-order corrections from the Holstein–Primakoff expansion. Focusing on the combined effects of single-ion anisotropy, Zeeman field, next-nearest-neighbor (NNN) exchange, and Dzyaloshinskii–Moriya interaction, we analyze how these parameters influence the magnetization curves and magnon spectra. This work serves two main goals. First, we critically examine the limitations of SRSWT, showing that in the absence of external or interaction tuning, the theory tends to overestimate magnon self-energy corrections, often predicting first-order magnetic transitions with multivalued magnetization and metastable solution branches (i.e., self-consistent but thermodynamically unstable states). Second, we demonstrate that topological transitions – signaled by magnon gap closings at the Dirac points – can be tuned to occur below the magnetic transition temperature and within the thermodynamically stable regime. In particular, we identify two practical tuning strategies: applying an external Zeeman field of appropriate sign depending on the anisotropy strength, and introducing a small antiferromagnetic NNN exchange coupling. These findings not only clarify the predictive scope and limitations of SRSWT but also provide experimentally relevant guidance for realizing thermally driven topological transitions in two-dimensional honeycomb magnetic insulators.

We develop a perturbative input-output formalism to compute the reflectance and transmittance spectra of cavity magnonics systems subject to a Floquet modulation. The method exploits the strong hierarchy between the magnetic-dipole couplings transverse (drive field) and parallel (modulation field) to the static bias field, which naturally introduces the small parameter $\epsilon = (2Ns)^{-1/2}$ associated with the total spin $Ns$ of the ferromagnet. By organizing the cavity and magnon fields in a systematic expansion in $\epsilon$, we obtain compact analytic expressions for the spectra up to second order. Using these results, we reproduce the characteristic sideband structure observed in recent Floquet cavity electromagnonics experiments. Furthermore, accounting for the Zeeman interaction between the modulation field and the fully polarized ground state - a contribution typically neglected in previous treatments - we predict an additional magnon detuning of approximately $0.8\,\mathrm{GHz}$, independent of both modulation frequency and sample size and determined solely by the spatial volume occupied by the modulation field. This identifies a measurable and previously overlooked shift relevant for the interpretation and design of cavity magnonics experiments.

A self-contained treatment of the Bogoliubov-Valatin transformation for homogeneous fermionic Hamiltonians is presented. The aim is to provide a quick reference that may also serve as supplementary material for a graduate-level course, and that can be understood with quantum mechanics knowledge up to the level of the second quantization's rules. The objective of the transformation is to cast a quadratic Hamiltonian into a diagonal form that resembles the Hamiltonian of a system of non-interacting particles. To obtain this, the first step consists in putting its coefficient matrix into its canonical form; the transformation can always be performed on fermionic Hamiltonians, only some care must be taken when this form is singular. Having explained how to cast a general matrix into its standard form, a complete description of the transformation is provided; a novel procedure is proposed here for the singular matrix case.

In Ref. [1], arXiv:2511.03384, Markos and Hlubina argue that "contrary to the expectations of Hirsch" [2] the conventional theory of superconductivity correctly describes the dynamics of the Meissner effect. Here I point out the flaws in their arguments that render them invalid, and propose an experiment to shed further light on these issues.

We provide an overview of the Gross-Pitaevskii-Poisson equation (GPPE) that is used to model self-gravitating superfluid systems, which include gravitationally collapsed boson and axion stars and dark-matter haloes. We outline how this framework can be used to develop minimal models for neutron stars and for pulsars and their glitches. We account not only for vortices in the neutron superfluid inside these stars, but also for the flux tubes in the proton-superconductor subsystem, using a coupled model with the neutron superfluid, proton superconductor, the Maxwell equations for the vector potential ${\bf A}$, and the Poisson equation for self-gravity.

We introduce a symmetry-adapted framework for simulating quantum dynamics based on Pauli propagation. When a quantum circuit possesses a symmetry, many Pauli strings evolve redundantly under actions of the symmetry group. We exploit this by merging Pauli strings related through symmetry transformations. This procedure, formalized as the symmetry-merging Pauli propagation algorithm, propagates only a minimal set of orbit representatives. Analytically, we show that symmetry merging reduces space complexity by a factor set by orbit sizes, with explicit gains for translation and permutation symmetries. Numerical benchmarks of all-to-all Heisenberg dynamics confirm improved stability, particularly under truncation and noise. Our results establish a group-theoretic framework for enhancing Pauli propagation, supported by open-source code demonstrating its practical relevance for classical quantum-dynamics simulations.

The deformation and rupture of a lipid vesicle due to the forced normal approach of an inclusion are essential for optimizing the design of magnetic giant unilamellar vesicles [magGUVs, Malik et al., Nanoscale 17, 13720 (2025)], with implications for active colloid-membrane interactions and cellular-scale chemical delivery. Here, we investigate vesicles propelled by a force-driven rigid inclusion and reveal a robust elastohydrodynamic mechanism: the inclusion outpaces the vesicle, sustaining a thinning film that drains symmetrically and self-similarly, largely independent of initial shape. For soft membranes and small inclusions, coupling drives a monotonic tension increase that can exceed the lysis tension. Evaluating the maximal tension over a delivery distance, we map an operating window in vesicle reduced area and size relative to the inclusion.

Despite decades of research, the universal nature of fluctuations in disordered quantum systems remains poorly understood. Here, we present extensive numerical evidence that fluctuations in two-dimensional (2D) Anderson localization belongs to the (1+1)-dimensional Kardar-Parisi-Zhang (KPZ) universality class. In turn, by adopting the KPZ framework, we gain fresh insight into the structure and phenomenology of Anderson localization itself. We analyze both localized eigenstates and time-evolved wave packets, demonstrating that the fluctuation of their logarithmic density follows the KPZ scaling. Moreover, we reveal that the internal structure of these eigenstates exhibits glassy features characteristic of the directed polymer problem, including the emergence of dominant paths together with pinning and avalanche behavior. Localization is not isotropic but organized along preferential branches of weaker confinement, corresponding to these dominant paths. For localized wave packets, we further demonstrate that their spatial profiles obey a stretched-exponential form consistent with the KPZ scaling, while remaining fully compatible with the single-parameter scaling (SPS) hypothesis, a cornerstone of Anderson localization theory. Altogether, our results establish a unified KPZ framework for describing fluctuations and microscopic organization in 2D Anderson localization, revealing the glassy nature of localized states and providing new understanding into the universal structure of disordered quantum systems.

The second-order optical susceptibility of semiconductors $\chi_{ijk}^{(2)}(-2\omega;\omega,\omega)$ finds application in metrology, spectroscopy, telecommunications, material characterization, and quantum information. Pioneering calculations of $\chi_{ijk}^{(2)}(-2\omega;\omega,\omega)$ utilized non-orthogonal Gaussian orbitals centered at atoms. That formulation transitioned into plane-wave-based algorithms as time went by. As of late, nevertheless, multiple tools for calculating optical susceptibilities have recast the problem using Wannier (i.e., localized) orbitals, making a comeback onto frameworks based on localized basis sets. Here, we present an approach for calculating $\chi_{ijk}^{(2)}(-2\omega;\omega,\omega)$ reliant on numerical pseudoatomic orbitals (PAOs) within perturbation theory in the velocity gauge. Its salient feature is a calculation of `Slater-Koster-like' two-center integrals of the momentum operator in between PAOs identified by symmetry. The approach was successfully tested on paradigmatic cubic silicon carbide (3C-SiC) and gallium arsenide, for which linear responses are contributed as well.

A nonrelativistic proof of the spin-statistics theorem is given in terms of the field operators satisfying commutation and anticommutation relations, which are introduced here in the coordinate space as a means to build the permutation symmetry into the brackets of identical particles. An eigenvalue problem of a $\pi$-rotation for a product of two annihilation operators is combined with an analysis on its rotational property to prove the connection that the field operators for integral-spin and half-integral-spin particles obey the commutation and anticommutation relations, respectively.

We show that codimension-two defects in Fermi liquids deform the renormalization group flow via a marginally relevant coupling. The mechanism for generating the flow is distinct from the case of the Kondo problem (codimension-three defects) in that the effective particle-hole asymmetry that leads to the log running is due to the spatial anisotropy generated by the defect. The mechanism for the log generation has a simple geometric explanation which shows that hole fluctuations are suppressed as the incoming momentum is taken to be along the direction of the defect. The RG flow time is shown to scale with the length of the defect. We also show that the dislon, the Goldstone mode localized to the defect, couples in a non-derivative fashion to the bulk fermions and becomes relevant above the dislons' Debye frequency which depends upon the defect tension.

NiO is a promising p-type material for photovoltaics and power electronics, but its temperature limits remain unclear. Using in situ high-temperature X-ray diffraction (HT-XRD) from 30 to 1100 C, we track the structural evolution of NiO thin films in air. The film crystallizes from an amorphous phase to cubic NiO between 300 and 400 C, where the emergence and growth of the (111) diffraction peak correlate with an increase in electrical resistivity. Further increases in temperature lead to improved crystallinity and higher resistivity. At 1100 C, the formation of Ni2O3 is observed, resulting in a highly resistive film. This study establishes a clear correlation between phase evolution, crystallinity, and resistive behavior in NiO thin films.

Randomly self-assembled nanowire networks (NWNs) are dynamical systems in which junctions between two nanowires can be modelled as memristive units viewed as adaptive resistors with memory. Various memristive models have been proposed to capture the complex mechanics of these junctions. Here, we showcase a novel computational framework named Memristive Nanowire Network Simulator (MemNNetSim) to simulate and analyze random memristive NWNs in a unified approach. Implemented using the Python programming language, MemNNetSim allows for the analysis of static and dynamic scenarios of NWNs under arbitrary memristive models. This provides a versatile foundation to build upon in further work, such as reservoir dynamics with NWNs, which has seen increased interest due to the interconnected architecture of NWNs. In this work, we introduce the package, demonstrate its utility in simulating NWNs, and then test advanced scenarios in which it can aid in the exploratory analysis of these systems, particularly in learning how to use NWNs as a physical reservoir in reservoir computing applications.

Advancing the chemical synthesis of crystals is important for both fundamental research and practical applications of quantum materials. While established bulk-phase and thin-film growth methods have enabled enormous progress, synthesizing single crystals suitable for quantum electronic discoveries remains challenging for many emerging materials. Here, we introduce van der Waals (vdW) stacks as nanochemical reactors for single-crystal synthesis and demonstrate their broad applicability in growing both elemental and compound crystals at the micrometer scale. By encapsulating atomically thin reactants that are stacked compactly with inert vdW layers such as hexagonal boron nitride (hBN), we achieve nanoconfined synthesis with the resulting crystals remaining encapsulated. As proof of concept, we synthesized isolated single crystals of elemental tellurium and distinct types of Pd-Te compounds. Structural characterization, including atomic-resolution scanning transmission electron microscopy, confirms the high crystalline quality of the products. We confirm the intrinsic semiconducting gap of tellurium and observe that non-stoichiometric PdTe1-x with a significantly reduced Te content (x ~ 0.18, a regime not previously achieved) retains uniform crystallinity and exhibits superconductivity below a critical temperature of 3.8 K. This nanochemical synthesis is broadly generalizable, chip-integrable, well-suited to a wide range of processing conditions, and compatible with nanofabrication routines for constructing devices. The concept of vdW nanoreactors offers a powerful and versatile pathway to expand the accessible landscape of quantum materials.

Detecting broken time-reversibility at micro- and nanoscale is often difficult when experiments offer limited state resolution. We introduce a lumping method that builds an effective semi-Markov model able to reproduce exactly the full entropy-production statistics of the microscopic dynamics. The mean entropy production stays accurate even when hidden current-carrying cycles are merged, though higher-order information can be unavoidably lost. In these cases, we capture violations of fluctuation theorems consistent with experiments, opening a path to novel inference strategies out of equilibrium.

We show that certain 't Hooft anomalies not detected by standard closed four-dimensional probes can become visible when a quantum field theory is placed on asymptotically locally Euclidean (ALE) spaces. As a concrete example, we use the Eguchi-Hanson (EH) space, whose defining features are its nontrivial second cohomology generated by the self-intersecting two-sphere and its asymptotic boundary $\mathbb{RP}^3$, which carries torsion and thus furnishes additional cohomological data absent on conventional backgrounds. For a theory with symmetry $G_1\times G_2$, we turn on background flux for $G_1$ and probe potential anomalies by performing a global $G_2$ transformation; the resulting anomaly is captured by a five-dimensional mapping torus. The anomaly receives contributions from the four-dimensional characteristic classes on EH space as well as from the $\eta$-invariant associated with the $\mathbb{RP}^3$ boundary. The anomaly detected in this way imposes additional constraints on asymptotically free gauge theories with fixed trivial asymptotic boundary conditions. In particular, infrared composite spectra that match anomalies on standard closed manifolds may nevertheless fail to reproduce the EH anomaly, and thus cannot by themselves furnish a complete symmetry-preserving infrared realization.

Canonical forms are central to the analytical understanding of tensor network states, underpinning key results such as the complete classification of one-dimensional symmetry-protected topological phases within the matrix product state (MPS) framework. Yet, the established theory applies only to uniform MPS with periodic boundary conditions, leaving many physically relevant states beyond its reach. Here we introduce a generalized canonical form for uniform MPS with a boundary matrix, thus extending the analytical MPS framework to a more general setting of wider physical significance. This canonical form reveals that any such MPS can be represented as a block-invertible matrix product operator acting on a structured class of algebraic regular language states that capture its essential long-range and scale-invariant features. Our construction builds on new algebraic results of independent interest that characterize the span and algebra generated by non-semisimple sets of matrices, including a generalized quantum Wielandt's inequality that gives an explicit upper bound on the blocking length at which the fixed-length span stabilizes to an algebra. Together, these results establish a unified theoretical foundation for uniform MPS with boundaries, bridging the gap between periodic and arbitrary-boundary settings, and providing the basis for extending key analytical and classification results of matrix product states to a much broader class of states and operators.

Tensor network methods, most prominently matrix product states (MPS), have become fundamental tools in modern quantum many-body physics. While MPS and extensions like the multiscale entanglement renormalization ansatz (MERA) and tree tensor networks (TTN) efficiently capture area-law entanglement and its logarithmic violations, they inherently struggle to represent highly entangled wavefunctions. Specifically, reaching the volume-law regime typically demands exponential resources within these conventional frameworks. Motivated by this challenge, we propose holographic isometric tensor network states (holographic isoTNS) that simulate quantum lattice models in $D$ spatial dimensions via $(D+1)$-dimensional networks of tensors. The additional dimension substantially enlarges the representational manifold, while isometric constraints on each tensor ensure efficient contractibility. Using one-dimensional systems as testbeds, we analyze the properties of holographic isoTNS. First, we show that randomly initialized holographic isoTNS typically display volume-law entanglement at fixed bond dimension. Second, through analytic constructions and variational optimization, we demonstrate that holographic isoTNS can faithfully describe a broad class of highly entangled yet low-complexity states. In particular, the ansatz can represent arbitrary fermionic Gaussian states, Clifford states, extensions of rainbow states, and certain short-time-evolved states under local evolution. Third, to exploit this expressivity in broader contexts, we implement a time-evolving block decimation (TEBD) algorithm on holographic isoTNS. While the method remains efficient and scalable, error accumulation over TEBD sweeps suggests the need for further algorithmic improvement. Overall, holographic isoTNS broaden the scope of tensor-network methods, opening new avenues to study physics in the volume-law regime.

The traditional multichannel Kondo effect takes place when several gapless metallic electronic channels interact with a localized spin-$S$ impurity, with the number of channels $n$ exceeding the size of the impurity spin, $n>2S$, leading to the emergence of non-Fermi liquid impurity behavior at low temperatures. Here, we show that the effect can be realized even when the electronic degrees of freedom are strongly correlated and gapped. The system under consideration consists of a single spin-$\frac{1}{2}$ impurity coupled isotropically to $n$ spin singlet superconducting channels realized by one-dimensional leads with quasi-long-range superconducting order. The competition between the Kondo and superconducting fluctuations induces multiple distinct ground states and boundary phases depending on the relative strengths of the bulk and boundary interactions. Using the Bethe Ansatz technique, we identify four regimes: an overscreened Kondo phase, a zero-mode phase, a Yu-Shiba-Rusinov (YSR) phase, and a local-moment phase with an unscreened impurity, each with its own experimental characteristic. We describe the renormalization-group flow, the excitation spectrum, and the full impurity thermodynamics in each phase. Remarkably, even in the presence of a bulk mass gap, the boundary critical behavior in the Kondo phase is governed by the same exponents as in the gapless theory with the low-energy impurity sector flowing to the $SU(2)_n$ Wess-Zumino-Witten (WZW) fixed point, and the impurity entropy monotonically decreasing as a function of temperature. In both the overscreened Kondo and zero-mode phases, the residual impurity entropy is $S_{\mathrm{imp}}(T \to 0) = \ln[2\cos(\pi/(n+2))]$. In the YSR and unscreened phases on the other hand the impurity entropy exhibits non-monotonic temperature dependence and is effectively free at low temperatures with $S_{\mathrm{imp}}(T \to 0) = \ln 2$.

We perform a systematic classification of (2+1)d Gross–Neveu–Yukawa-like models built out of one or more 4-component Dirac fermions and $M$ scalar fields, which preserve an O($M$) symmetry rotating the scalars. We then identify the perturbative fixed points of these models in the $4-\epsilon$ expansion. Our classification highlights several targets for the conformal bootstrap and reveals a new fixed point with $M=3$, which we call the "orthogonal Heisenberg" CFT.

We present an attention-based foundation model architecture for learning and predicting quantum states across Hamiltonian parameters, system sizes, and physical systems. Using only basis configurations and physical parameters as inputs, our trained neural network is able to produce highly accurate ground state wavefunctions. For example, we build the phase diagram for the 2D square-lattice $t-V$ model with $N$ particles, from only 18 parameters $(V/t,N)$. Thus, our architecture provides a basis for building a universal foundation model for quantum matter.

We find an unexpected phenomenon of coherently synchronized oscillations in a mirror-symmetric many-body localized system. A synchronization transition of the spin oscillations is found by changing the spin-spin interactions. To understand this phenomenon, an effective Ising model based on local integrals of motion is proposed. We find that the synchronization transition can be understood as a paramagnetic-to-ferromagnetic Ising transition. Based on the Ising model, we theoretically estimate the synchronized frequencies and the synchronization transition points, which agree well with numerical results.

Agentic AI systems increasingly connect large language models (LLMs) to external scientific tools, yet whether and when tool access improves prediction accuracy remains uncharacterized. We present AGAPI (AtomGPT.org API), an open access platform integrating eight open-source LLMs with 18 REST endpoints (28 agent tools, 50 web apps) spanning materials databases, force fields, tight-binding band structures, X-ray diffraction, and protein structure. A three-evaluation residual decomposition on JARVIS-Leaderboard electronic-structure test sets separates agent pipeline fidelity from inherited density functional theory (DFT) functional bias. For bulk modulus and bandgap the agent reproduces JARVIS-DFT entries to numerical precision, so the experimental-reference degradation is functional bias, not agentic malfunction. On memorization-resistant test sets (57 defective supercells, 60 hypothetical compositions), tool-augmented mean absolute error (MAE) is below 0.005 eV versus 1.25 to 1.86 eV tool-free, confirming tools are indispensable where parametric knowledge is unavailable. We further demonstrate autonomous multi-step workflows including 10-operation defect-engineering pipelines. AGAPI is available at https://github.com/atomgptlab/agapi.