2026-03-24 | | Total: 159
We study the phase diagram of interacting spinless fermions on the honeycomb bilayer at charge neutrality using large-scale quantum Monte Carlo simulations. In the noninteracting limit, the low-energy spectrum features quadratically dispersing bands that touch at the corners of the hexagonal Brillouin zone. Weak to intermediate interactions induce a splitting of each of the quadratic band touching points into four Dirac points, located along high-symmetry directions of the reciprocal lattice. Strong interactions lead to the formation of a layer-polarized charge density wave, which spontaneously breaks the $\mathbb Z_2$ layer inversion symmetry and opens an insulating gap in the spectrum. We show that the semimetal-to-insulator quantum phase transition as a function of interaction is continuous and characterized by emergent relativistic symmetry. Our results for the values of the correlation-length exponent $ν$, the order-parameter anomalous dimension $η_φ$, and the fermion anomalous dimension $η_ψ$ agree with those of the theoretically predicted 2+1D Gross-Neveu-Ising universality class with eight two-component Dirac fermions within less than 5\%\ deviation. We also determine the crossover scale as a function of interaction strength between the nonrelativistic semimetal state at high temperatures, characterized by dynamical critical exponent $z = 2$, and the Dirac semimetal state at intermediates temperatures, characterized by $z=1$. Further reducing the temperature below the crossover scale at a fixed value of the interaction strength above the quantum critical point results in a classical ordering transition in the 2D Ising universality class.
Sodium-ion batteries require anodes that combine high capacity, low operating voltage, fast Na-ion transport, and mechanical stability, which conventional anodes struggle to deliver. Here, we use the SpookyNet machine-learning force field (MLFF) together with all-electron density-functional theory calculations to characterize Na storage in aminobenzene-functionalized Janus graphene (Na$_x$AB) at room-temperature. Simulations across state of charge reveal a three-stage storage mechanism-site-specific adsorption at aminobenzene groups and Na$_n$@AB$_m$ structure formation, followed by interlayer gallery filling-contrasting the multi-stage pore-, graphite-interlayer-, and defect-controlled behavior in hard carbon. This leads to an OCV profile with an extended low-voltage plateau of 0.15 V vs. Na/Na$^{+}$, an estimated gravimetric capacity of $\sim$400 mAh g$^{-1}$, negligible volume change, and Na diffusivities of $\sim10^{-6}$ cm$^{2}$ s$^{-1}$, two to three orders of magnitude higher than in hard carbon. Our results establish Janus aminobenzene-graphene as a promising, structurally defined high-capacity Na-ion anode and illustrate the power of MLFF-based simulations for characterizing electrode materials.
Global entanglement in quantum many-body systems is inherently nonlocal, raising the question of whether it can be inferred from local observations. We investigate this problem in monitored quantum circuits, where projective measurements generate classical records distributed across spacetime. Using graph neural networks (GNNs), we represent individual quantum trajectories as directed spacetime graphs and reconstruct the half-chain entanglement entropy from local measurement data alone. Because information propagates through the network via local message passing, the architecture directly controls the spacetime region over which correlations can be aggregated. By systematically varying this accessible scale -- through network depth and hierarchical spacetime coarse-graining -- we probe how much measurement information is required to reconstruct global entanglement. We find that prediction accuracy improves as the accessible spacetime region grows and that results from different architectures collapse when expressed in terms of an effective spacetime scale combining depth and coarse-graining. These results demonstrate that the information required to reconstruct global entanglement is organized in spacetime scales and show that graph-based learning architectures provide a controlled operational framework for probing how global quantum correlations emerge from local measurement data.
We predict a supercurrent-driven Néel spin-orbit torque in a superconductor/$d$-wave altermagnet heterostructure, associated with the emergence of spin-triplet correlations. The effect can be understood as a consequence of the supercurrent-induced spin polarization, owing to the interplay between spin-orbit coupling and momentum-dependent spin splitting, as found, for example, in altermagnets. Remarkably, the supercurrent can be tuned by the Néel-vector direction, and the supercurrent-induced torque can both propel magnetic domain walls and reverse the Néel-vector orientation within a domain wall. These findings establish superconductor/altermagnet heterostructures as a versatile platform for the dissipationless control of the Néel vector, with potential applications in racetrack memory, dissipationless superconducting electronics, and unconventional computing.
We establish an exact Wilson-loop conjugation identity, $W(-δ)=W(δ)^*$, for the many-body overlap Wilson-loop $W(δ)$ accumulated along a $U(1)$ flux-threading (twist) cycle parametrized by $θ\in[0,2π]$, where $δ$ denotes the bond-dimerization parameter. A minimal sufficient condition is the existence of a composite antiunitary mapping acting on the flux-threaded Hamiltonian family that implements $(δ,θ)\mapsto(-δ,-θ)$. As a concrete demonstration, we construct such a mapping microscopically in a dimerized staggered Hubbard ring at half filling. We then verify the conjugation identity using the density-matrix renormalization group (DMRG) for gapped, nondegenerate ground states along the twist cycle. Importantly, the identity persists in depinned gapped regimes where the Berry-phase $γ\equiv-\arg W$ is not symmetry-quantized; as a corollary, $γ(-δ)=-γ(δ)$ (mod $2π$). More generally, the same conjugation relation applies to any lattice model whose flux-threaded Hamiltonian family is closed under an orientation reversal of the bond pattern (a suitable permutation of link-hopping parameters) combined with a reversal of the flux orientation.
Ultrathin films of widely different materials exhibit a dramatic enhancement of projectile penetration resistance under high--velocity impact. Despite extensive simulations and experiments, a unifying physical explanation has remained elusive. Here we show that the thickness dependence of the specific penetration energy obeys a universal law, $E_p^*(h)=E_{p,\infty}^*+B h^{-3}$, independent of chemical composition and degree of disorder. The inverse--cube scaling is traced back to a finite--size correction to the effective shear modulus arising from the suppression of long--wavelength nonaffine deformation modes in confined solids. The scaling quantitatively describes impact data for multilayer graphene, graphene oxide, and polymer thin films, revealing a common elastic origin for nanoscale impact resistance.
We report on the fabrication and characterization of patterned high-mobility two-dimensional electron gases (2DEG) formed on SrTiO$_3$ (STO) substrate surfaces by hydrogen plasma exposure. The resulting devices consistently showed high electron mobilities up to 7400 cm$^2$/V$\cdot$s. A large range of electron density was systematically explored by controlled aging of the sample between cooldowns, including the expected range for the STO 2DEG superconducting dome. No superconducting transition was observed down to the base temperature of approximately 10 mK. This suggests suppression of superconductivity in high mobility quasi-two-dimensional SrTiO$_3$ electron gas, likely linked to vertical confinement and electronic orbital rearrangement. We systematically explored electrostatic gate modulation in this 2DEG system and its scaling with electron density and side gate geometry. In contrast with our initial expectation, we observed an improvement of achievable total modulation for larger side gate to channel separation. At low electron density, stochastic channel pinch-off events were observed, creating quasi-ballistic constrictions with irregular conductance quantization. This epitaxy-free and high mobility oxide material platform offers a promising new route towards patterning quantum devices.
The thermodynamic critical field $B_{\rm c}$ provides direct access to the superconducting condensation energy, yet its pressure dependence has been studied much less extensively than that of the transition temperature. Here, muon-spin-rotation/relaxation measurements of the thermodynamic critical field $B_{\rm c}$ of elemental Pb under hydrostatic pressure up to $\simeq2.3$ GPa are reported. From the magnetic-field distribution in the intermediate state, $B_{\rm c}(T)$ is determined and $B_{\rm c}(0)$ is extracted at different pressures. In combination with previously reported high-pressure data for $B_{\rm c}$ and $T_{\rm c}$, it is shown that the pressure dependence of $B_{\rm c}(0)$ follows that of the superconducting gap $Δ(0)$ more closely than that of the transition temperature $T_{\rm c}$. At higher pressures, the logarithmic pressure derivatives of $B_{\rm c}(0)$ and $T_{\rm c}$ are found to converge, indicating that the coupling strengths ratio $α=Δ(0)/k_{\rm B}T_{\rm c}$ becomes nearly pressure independent. This behavior is interpreted as thermodynamic evidence for a pressure-driven crossover from strong- to weak-coupling superconductivity in Pb.
Previous experimental reports of long-range interactions in ionic liquids (ILs) stand in contradiction with theoretical predictions and numerical simulations. To provide insights into the literature discrepancies regarding the experimental ranges of electrostatic screening, claimed with orders of magnitude larger, the interactions between pairs of mica and borosilicate surfaces confining ILs are investigated by two complementary advanced Surface Force Apparatuses. Regardless of differences in confinement geometries (crossed-cylinders, sphere-flat), radii of curvature (cm-mm), and measurement techniques (stepwise vs continuous approach), two ever present force regimes are evidenced. At small surface separations, oscillatory forces reflect IL structuration and layering, while outside this gap, the interaction is monotonic repulsive. In both regimes the spatial extent and force magnitude depend critically on motion conditions, as demonstrated by achieving velocities as low as 9 pm/s with equilibration times up to 90 s. At large separations, fast surface displacements generate long-range interactions (over tens of ion size) creating the illusion of anomalous underscreening, whereas increasingly slow ones shrink both magnitude and range of the repulsion with decay-lengths converging ultimately to a screening length consistent with Poisson-Boltzmann theory with finite ion sizes. The transition from apparent long-range to short-range screening unfolds over nearly two orders of magnitude in time, revealing slow relaxation dynamics reminiscent of aging phenomena. These findings definitely resolve a decade-old controversy on force measurements and reveal rich out-of-equilibrium dynamics. The hydrodynamic contribution to the net force is admittedly crucial to be reduced especially when relaxations span decades in time, but approaching thermodynamic equilibrium during measurements proves essential.
Quantum anomalies arise when symmetries of a classical theory cannot be preserved upon quantization, leading to unconventional topological responses. A prominent example is the parity anomaly of a single two-dimensional Dirac fermion, which enforces a half-quantized Hall response. Anomaly inflow mechanism allows this effect to be observed at the surfaces of three-dimensional topological insulators, however, its realization in a genuinely two-dimensional system has remained elusive. Here we report the observation of a parity-anomalous Hall response at the critical point of a quantum Hall topological phase transition in a synthetic two-dimensional system of ultracold dysprosium atoms. By coupling a continuous spatial dimension to a finite synthetic dimension encoded in atomic spin states, we engineer tunable Chern bands with C = 0 and 1. At the transition, the bulk gap closes at a single Dirac point, where we observe a robust half-quantized Hall drift despite strong non-adiabatic excitations. We show that this response originates from the global structure of the band topology, is protected by an emergent parity symmetry at criticality, and disappears when parity is explicitly broken. Our work establishes synthetic quantum systems as a powerful platform to probe quantum anomalies and their interplay with topology and non-equilibrium dynamics.
Recently, a class of spin chains known as ``free fermions in disguise'' (FFD) has been discovered, which possess hidden free-fermion spectra even though they are not solvable via the standard Jordan-Wigner transformation. In this work, we extend this FFD framework to open quantum systems governed by the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) equation. We establish a general class of exactly solvable open quantum systems within the FFD framework: if the Liouvillian frustration graph is claw-free and has a simplicial clique, the Liouvillian possesses a hidden free-fermion spectrum. In particular, the (even-hole, claw)-free condition automatically guarantees this, enabling exact computation of the Liouvillian gap and an infinite-temperature autocorrelation function. Our results provide the first realization of the FFD mechanism in open quantum systems.
Among two-dimensional magnetic materials, CrSBr has attracted considerable attention owing to its coexistence of ferromagnetic and antiferromagnetic ordering, which depends sensitively on crystallographic orientation. An additional distinguishing feature of CrSBr is its highly anisotropic Fermi surface in momentum space. In this work, we present a comprehensive investigation of magnetoresistance by systematically orienting the bias current and the applied magnetic field along all three crystallographic axes. We demonstrate that the magnetoresistance serves as a direct probe of electronic anisotropy, exhibiting pronounced variations when the current is applied along different crystallographic directions under a magnetic field perpendicular to the sample plane. For in-plane magnetic fields, we observe conventional anisotropic magnetoresistance accompanied by hysteresis, indicative of ferromagnetic behavior. Overall, our study provides a complete picture of electronic transport in CrSBr as a function of bias current and magnetic field orientation with respect to crystallographic directions, thereby opening pathways for future experiments requiring high sensitivity of electrical resistance to magnetic field gradients.
Charge transport in superconducting junctions at finite voltages is governed by Andreev reflections, including multiple Andreev reflections, which are processes that enable multiple charge transfer, a hallmark that shot noise can directly quantify. Since the effective charge extracted from shot noise measurements varies with the transparency of the junction, systematic control of transparency is essential but experimentally challenging. Here, we present shot noise scanning tunneling microscopy measurements enabled by a newly developed amplifier, allowing access to different transparency regimes. We perform shot noise measurements on Pb(111) with tunable transparency at 2.2 K and observe that the shot noise evolves from a single electron tunneling regime to multiple charge transfer regime as transparency increases. Our results are quantitatively consistent with theoretical simulations of Andreev reflections and multiple Andreev reflections for a single-channel system. These results establish junction transparency as the key parameter governing the evolution of charge transport and demonstrate that noise-STM is a powerful platform for investigating microscopic charge transport mechanisms with controlled junction transparency at the atomic scale.
Manganese(II) mobilised by mining activity poses a persistent water-quality challenge, yet the mechanisms by which low-cost sorbents, such as biochar, sequester Mn(II) remain poorly resolved. This study identifies the specific chemical drivers of Mn(II) sequestration by combining fixed-bed column and batch experiments with atomistic molecular dynamics simulations. Oilseed rape straw biochars, produced at 350\textdegree C, 550\textdegree C, and 700\textdegree C, removed 20-50% of dissolved Mn from acidic influent (pH 4, 5 ppm). High-temperature biochar achieved the greatest removal ($\sim$50%) and rapidly increased effluent pH to 9, triggering alkaline precipitation. Conversely, lower-temperature biochars removed 20-30% of Mn while maintaining a near-neutral pH (7-7.5). Enhanced \ce{K+} release in these systems indicates significant cation exchange and non-precipitative pathways. Molecular simulations confirmed that while neutral surfaces show weak Mn(II) association, deprotonated sites drive strong adsorption through inner-sphere complexation ($\sim$50% removal) and outer-sphere association ($\sim$10%). These results establish a mechanistic framework to distinguish between precipitation-led and surface-complexation-led removal. By providing specific chemical criteria for Mn-targeted sequestration, this work enables the rational design of engineered biochars for sustainable water remediation.
For ethylene carbonate on the (100) surface of lithium, we calculate the adsorption energy in two binding motifs as well as the barrier height for a ring-opening decomposition reaction. We validate a scheme for producing results in the thermodynamic limit by correcting results obtained on finite lithium clusters containing only 40-100 atoms, which enables the use of hybrid density functionals, the random-phase approximation, and correlated wavefunction theories such as coupled-cluster theory and auxiliary-field quantum Monte Carlo. We find that the high-level theories agree to within 2-5 kcal/mol and can therefore serve as benchmarks for more affordable methods. Using our reference data, we demonstrate that generalized gradient approximation functionals, such as PBE, are not sufficiently accurate for reaction barrier heights, and we identify $ω$B97X-V as an especially promising functional for the interfacial chemistry of electrolyte solvents at lithium metal anodes.
Nonlinear response functions, formulated as multipoint correlation functions or Volterra kernels, encode the dynamical and spectroscopic properties of physical systems and underpin a wide range of nonlinear transport and optical phenomena. However, their evaluation rapidly becomes prohibitive at high orders because of combinatorial (often factorial) scaling or severe numerical errors. Here, we establish a systematic and efficient framework to compute nonlinear response functions directly from real-time dynamics, without explicitly constructing multipoint correlators or relying on numerically unstable finite-difference methods for order-resolved extraction. Our approach is based on the Gateaux derivative with respect to the external field in function space, which yields a closed hierarchy of tangent equations of motion (TEOM). Propagating the TEOM alongside the original dynamics isolates each perturbative order with high accuracy, providing a term-by-term decomposition of physical contributions. The computational cost scales exponentially with response order in the fully general setting and reduces to polynomial complexity when all perturbation directions are identical; both regimes avoid the factorial scaling of explicit multipoint-correlator evaluations. We demonstrate the power of TEOM by computing frequency-resolved fifth-order response functions for a solid-state electron model and by obtaining nonlinear response functions up to the 49th order with controlled accuracy in a classical Duffing oscillator. We further show that our time-evolution formulation allows optical conductivities to be evaluated directly while remaining numerically stable even near zero frequency. TEOM can be incorporated seamlessly into existing real-time evolution methods, yielding a general framework for computing nonlinear response functions in quantum and classical dynamical systems.
Graphene quantum point contacts (QPCs) in the quantum Hall regime host competing transport mechanisms including chiral edge propagation, valley degeneracy, and gate-induced mode mixing. Their interplay is not visible in conductance alone. Shot noise directly probes the statistics of transmission eigenvalues, revealing microscopic mode partitioning that conductance cannot access. We develop a hybrid framework combining tight-binding simulations of gate-defined graphene QPCs with random matrix theory (RMT) to predict shot noise and Fano factor signatures across different quantum Hall regimes, validated against experimental conductance maps of hBN-encapsulated graphene Hall bars. Three distinct regimes are identified: adiabatic propagation, sharp mode filtering, and multi-mode mixing driven by localized states beneath the split gate. For higher Landau levels ($N_L > 0$), complete mode mixing produces the universal chaotic-cavity limit $F \simeq 1/4$. Strikingly, the zeroth Landau level ($N_L = 0$) converges to $F = 1/3$. This distinct value originates in the sublattice polarization of the $N_L = 0$ edge state: coupling to mixed-sublattice localized states beneath the gate is suppressed, confining transport to an effective single channel ($N = 1$). Complete mixing within this single channel yields a flat transmission eigenvalue distribution and hence exactly $F = 1/3$ from single-channel RMT, numerically coincident with but mechanistically distinct from pseudo-diffusive zero-field graphene transport. The $F = 1/3$ versus $F = 1/4$ crossover is a Landau-level-resolved noise signature absent in conductance, providing a direct discriminator between single-channel and multi-channel chaotic transport in graphene QPCs.
Interlayer coupling plays a critical role in van der Waals materials by governing lattice stability and emergent quantum phases, yet its impact on few-layer hexagonal CoTe$_2$ remains unclear. Here, using first-principles calculations, we systematically investigate monolayer and bilayer CoTe$_2$ with an emphasis on their electronic structures, lattice dynamics, and electron-phonon coupling, and elucidate the underlying mechanisms driven by interlayer interactions. Our results show that monolayer CoTe$_2$ exhibits pronounced dynamical instability at low temperatures, whereas interlayer coupling stabilizes the bilayer crystal structure and gives rise to phonon-mediated superconductivity with a predicted critical temperature of about $4.7$~K. The stabilization and superconductivity in bilayer CoTe$_2$ are primarily attributed to interlayer-coupling-induced Te-$p_z$ charge redistribution and the associated modification of the Fermi surface and electron-phonon coupling. Finally, we discuss how spin-orbit coupling in bilayer CoTe$_2$ weakens the EPC and suppresses superconductivity. Our work clarifies how interlayer coupling can jointly tune structural stability and superconductivity in few-layer CoTe$_2$, providing insights for engineering quantum phases in layered transition-metal dichalcogenides.
Traditional percolation theory assumes static microscopic rules, limiting its ability to describe real-world complex systems where macroscopic order actively regulates local interactions. Here, we introduce feedback percolation, an unified framework that dynamically couples the microscopic activation probability to the macroscopic size of the giant component. We show that this simple feedback mechanism produces a rich variety of behaviors both analytically and numerically. Depending on the feedback functions, the system exhibits explosive discontinuous jumps, hybrid transitions, limit-cycle oscillations, and routes to chaos, absent in classical percolation. Our findings establish that macroscopic feedback provides a unifying physical mechanism for phenomena ranging from self-regulating oscillations to systemic infrastructure collapse.
The two-dimensional Heisenberg spin-glass model is investigated by means of a semiclassical expansion around classical states. At leading order, we obtain an effective quadratic spin-wave Hamiltonian and study the localization properties of its spectrum and eigenfunctions. We find that the nature of the spin-wave excitations, whether they are hydrodynamic or localized modes, depends crucially on the relevance/irrelevance - in the renormalization group sense - of the correlations induced by the underlying classical order in the spin-wave Hamiltonian matrix elements: low-energy excitations around magnetically ordered states are delocalized, whereas those around spin-glass ordered states are localized, albeit weakly. Remarkably, in the magnetically ordered case, spin-wave delocalization is robust with respect to the presence of disorder, even in two spatial dimensions. We interpret this phenomenology by relating the spontaneous breaking of spin-rotation symmetry in the original Heisenberg model to the symmetry and universality class of the resulting quadratic spin-wave Hamiltonian. We conjecture that the hydrodynamic picture can be recovered through the inclusion of interactions among the spin-wave excitations at higher order in the semiclassical expansion, favoring the onset of ergodic behavior.
We comment on the recent work by Yamaguchi and Barré [Phys. Rev. E 107, 054203 (2023)], which uses linear stability analysis of the Vlasov equation to characterize phase transitions in a generalized Hamiltonian Mean Field (gHMF) model. By performing extensive molecular dynamics simulations with $N=10^8$ particles, we demonstrate that the bifurcation analysis of the initial stationary distribution is insufficient to predict either the location or the nature of the phase transition to a quasi-stationary state (qSS). Specifically, we show that for bimodal momentum distributions, the instability threshold identified by the authors does not correspond to a ferromagnetic transition; instead, the system remains in a paramagnetic state characterized by magnetization oscillations with a zero time-average. We find that the true paramagnetic-ferromagnetic transition is discontinuous (first-order) and occurs at significantly larger coupling strengths, characterized by a clear coexistence of states. These results indicate that linear bifurcation and symmetry-breaking phase transitions are distinct phenomena in long-range interacting systems, and that the former lacks the predictive power to describe the long-time fate of the system.
We study experimentally the nonlinear propagation of short pulses of forward volume spin waves in nanometer-thick YIG films. We show that nonlinearity of the spin system can efficiently counteract dispersion broadening of the pulses, leading to the formation of envelope solitons. We demonstrate that in microscopic YIG systems, microwave powers of the order of one milliwatt are sufficient to reach the soliton formation threshold. At powers slightly above this threshold, we achieve transmission of 3-ns spin-wave pulses over distances of up to 50 micrometers without increase in their temporal width. Our results demonstrate a promising way towards high-rate transmission of information in microscopic spin-wave circuits unaffected by detrimental dispersion effects.
Landau theory usually treats free-energy coefficients as intrinsic parameters fixed by thermodynamic variables. We show that externally written microscale fields can survive coarse graining and enter the free-energy functional as spatially prescribed coefficient fields. This defines a nonintrinsic sector of Landau theory. The key condition is a hierarchy of correlation, writing, and frustration lengths. We identify ion-patterned FeRh as a plausible realization.
The Drinfeld center fusion category $\mathcal{Z}(\mathrm{Vec}_G)$ famously models anyons in certain lattice models. Here we demonstrate how its fusion rules may also describe topological order in fractional topological insulator materials, in the vicinity of point defects in the Brillouin zone. Concretely, we prove that $\mathcal{Z}(\mathrm{Vec}_G)$ reflects, locally over a punctured disk in the Brillouin zone, the monodromy (topological order) of gapped quantum states over the parameter space of Bloch Hamiltonians whose classifying space has fundamental group $G$.
Out-of-plane and in-plane electric polarization, which rarely coexist in a two-dimensional (2D) ferroelectric material, offer different advantages in ferroelectricity-based devices. Here, we report the coexistence of in-plane and out-of-plane electric dipoles, along with various phase transitions, in 2D van der Waals layered TlGaS2 single crystal. Quantum paraelectricity was observed along both in-plane and out-of-plane directions of the TlGaS2 crystal. Detailed investigation of the quantum paraelectric soft-mode behavior reveals a close correlation between the electric dipoles and the off-center displacement of Tl1+ ions with 6s2 lone pairs in TlGaS2. Anomalies near temperatures of about 120 K and 60-75 K in dielectric and/or infrared spectra indicate the existence of local or weak long-range structural transitions in TlGaS2. Our results provide important experimental evidence for elucidating the phase transitions and coexistence of in-plane and out-of-plane electric dipoles in 2D layered TlGaS2.