We propose a new Brownian bridge construction for our newly devised Coherent States Path Integral Monte Carlo algorithm. We apply it to the numerically exact calculation of the thermodynamic properties of the Helium fluid on a plane at low non zero temperature. We find very good agreement with the conventional plane waves path integral Monte Carlo results.

Fractional anomalous Hall crystals (FAHCs) replicate the topological order of the fractional quantum Hall effect in the continuum without requiring any external magnetic field. They spontaneously break continuous translation symmetry like a Wigner crystal, but are distinguished by each unit cell holding a fixed fractional number of electrons. Until now, these states have been confined to theoretical speculation or engineered models, leaving open the question of whether they can plausibly emerge in actual physical systems. Here, we establish them as energetically competitive candidate states in a realistic material setting. We study rhombohedral pentalayer graphene (R5G) with variational wavefunctions that are exact zero modes of a recently proposed ideal model of R5G. We evaluate their energies using Monte Carlo, after reinstating realistic dispersion and screened Coulomb interactions. We find FAHCs to be energetically competitive with integer anomalous Hall crystals and Fermi liquids, and their stability follows a simple principle. Each crystal maps onto a parent quantum Hall liquid that fixes its interaction energy, while the kinetic energy favors crystal periods that match the finite-momentum minimum of R5G's Mexican-hat dispersion. A weak periodic potential can then selectively lower and pin the commensurate fractional crystals. This picture predicts how the integer and fractional quantum anomalous Hall stability windows evolve with twist angle and displacement field, which we compare to recent experiments. These results support a continuum-and-interactions-first route to fractional anomalous Hall states in rhombohedral graphene.

Topological semimetals host band nodes characterized by quantized invariants that can appear in bulk responses, yet some invariants remain hidden from standard probes. In particular, band nodes can carry secondary topological charges whose transport signatures are still largely unexplored. Here we study three-dimensional nodal-line semimetals in which nodal rings carry both the Berry phase $w_1\pi$ and a $\mathbb{Z}_2$ monopole charge $w_2$. We show that magnetic quantum oscillations, usually treated as a probe of $w_1$, can directly diagnose $w_2$, with the relevant signal selected by the magnetic-field direction. For a field along the ring axis, the inner and outer extremal orbits of the toroidal Fermi surface both encircle the $w_2$-enforced thread and exhibit a topological phase shift $\nu w_2\pi$ in the $\nu$th harmonic, which is accessible through standard phase-resolved quantum-oscillation analysis. By contrast, for a field applied perpendicular to the ring axis, the relevant extremal orbit exhibits the usual $\pi$ phase shift associated with the Berry phase $w_1\pi$, independent of $w_2$. For weak doping, three-dimensional ABC-stacked graphdiyne is predicted to exhibit the proposed oscillations in a field range accessible with present-day high-field facilities.

Certain 3D $U(1)$ spin liquids, such as those arising in dipolar quantum spin ice, have an emergent monopole which is the source of an emergent magnetic field that transforms under symmetries like an electric polarization. As a consequence, these monopoles carry a physical electric charge in their cores and form a plasma at low temperatures. Due to the monopole coupling to emergent gauge fields, the full system behaves as an electrical insulator for DC transport, but can display a sharp plasma resonance analogous to a metal at very low frequencies. This can serve as a clear fingerprint to detect these states in materials. We discuss the optimal conditions to observe this phenomenon in the 3D $U(1)$ spin-liquid candidate Ce$_2$Zr$_2$O$_7$.

Analyzing the dynamics of physical observables in open quantum many-body systems is a fundamental but highly challenging task that has yielded very few exact results. In this work, we identify intrinsic conformal structures that restore exact solvability in $(1+1)$D conformal field theories. For $N$ Majorana fermions with linear mode jumps, the adjoint Lindbladian is triangular on reduced even Majorana monomials, yielding recursive exact Heisenberg evolution. In Wess-Zumino-Witten models admitting conformal Majorana embeddings, this hierarchy gives exact dynamics of affine-current products realized as Majorana bilinears, including regimes where the Kac-Moody current algebra alone does not close. In diagonal rational conformal field theories, Verlinde topological defect lines furnish jump operators whose primary-sector dynamics is exactly diagonal: topological-charge probabilities are conserved, while intersector coherences dephase at rates fixed by the modular $S$ matrix and nonnegative measurement strengths. These examples show that intrinsic conformal structures, such as conformal embeddings and modular data, can organize exactly solvable open conformal dynamics.

Altermagnetic order gives rise to momentum-dependent spin splitting of electronic and magnonic excitations even in the absence of a net magnetization. Here, we investigate the magnetic field dependence of the spin-wave spectrum of altermagnetic $\alpha$-MnTe using inelastic neutron scattering and linear spin-wave theory. An in-plane magnetic field continuously reorients the Néel vector by overcoming the weak crystalline anisotropy, while remaining small compared with the dominant exchange scale. We find that this reorientation leaves the magnon energies and line widths essentially unchanged, but strongly modifies the measured spectral intensity through the transverse-momentum projection. Our results demonstrate a clear separation between the soft orientational degree of freedom of the antiferromagnetic order and the robust exchange-dominated chiral magnon spectrum. This combination establishes $\alpha$-MnTe as a platform for reconfigurable magnon coupling, in which external fields tune how excitations interact with polarized probes without substantially altering their frequency or coherence.

The Casimir effect is typically governed by intrinsic material properties and lacks in situ tunability. We show that, in mesoscopic rings, both the magnitude and sign of the fermion-mediated Casimir interaction can be controlled via the Aharonov-Bohm effect. The resulting interplay between the Aharonov-Bohm phase and the Casimir interaction provides a route to engineer long-range interactions. In particular, this mechanism enables the implementation of non-Clifford SWAP$^{\alpha}$ gates between spatially separated spin qubits, thereby reducing the overhead for universal quantum computation and quantum error correction in spin-qubit architectures.

The Cooper pair momentum in a superconductor is associated with a phase gradient of the superconducting order parameter. In general, this momentum is small compared to the Fermi momentum, which makes it challenging to measure. Josephson junctions, however, enable the creation of large phase gradients and transfer of the Cooper pair momentum to quasiparticles via Andreev reflection. In this work we demonstrate that Andreev bound states propagating along ballistic planar Josephson junctions eject into an adjacent normal region at a phase-controlled angle that scales as $\Theta \sim \sqrt{\Delta/ \mu}$, where $\Delta$ is the superconducting gap and $\mu$ is the chemical potential. This angle parametrically exceeds the conventional Cooper pair momentum scale $\Delta/ \mu$, and thus this phenomenon is sizeable even within the Andreev approximation regime $\Delta / \mu \ll 1$. Our results establish phase-controlled quasiparticle ejection as a kinematic probe of condensate momentum transfer: unlike existing probes that detect the Doppler energy shift, the signal appears as a momentum-space deflection of emitted quasiparticles.

A pair of parallel polarizable planar walls at distance d is considered. The walls are symmetrically charged with a uniform surface charge density, neutralized by mobile point counterions moving between them. The case of repulsive particle images is studied in the strong-coupling (SC) regime. Of interest is the dependence of the effective inter-wall interaction (pressure), mediated by the mobile counterions, as a function of the distance d. It is shown that previous virial SC single-particle theories work well at small d when the dielectric jump is small; for intermediate and large dielectric jumps they are inadequate even in the SC region. Here, we propose a Wigner-type SC theory based on harmonic deviations of particles from their ground-state monolayer or bilayer Wigner structures formed inside the space between the dielectric walls. Our Monte-Carlo simulations are in very good agreement with the Wigner SC predictions, even down to moderate coupling constants ($\Xi > 10$).

We investigate the role of the CaCl2-type (Pnnm) phase in the high-pressure transformation of rutile TiO2, whose experimental signature has remained elusive. While analogies with other rutile-type oxides suggest such an intermediate, TiO2 typically exhibits a direct transformation to higher-coordination phases such as baddeleyite. Using an all-electron density functional framework combined with density functional perturbation theory, we show that rutile TiO2 undergoes a ferroelastic instability characterized by the development of an orthorhombic strain and a double-well energy landscape at 13.5 GPa. This instability is associated with the softening of the C11 - C12 elastic combination and the condensation of a B1g phonon mode, involving coordinated rotations of TiO6 octahedra that lower the symmetry to the Pnnm structure. Despite this clear elastic and dynamical pathway, enthalpy calculations show that the CaCl2-type phase is only weakly stabilized relative to rutile and remains energetically unfavorable compared to competing columbite and baddeleyite phases. Consequently, the Pnnm phase does not emerge as a stable high-pressure polymorph but instead exists as a transient or weakly metastable intermediate. These results demonstrate that the CaCl2-type phase represents the intrinsic ferroelastic response of rutile TiO2, yet is suppressed by thermodynamic competition, providing a consistent and unified explanation for its elusive experimental observation.

We investigate the formation of loop-current ordered phases in half-filled SU($N$) two-leg fermionic ladders. Using a low-energy approach, we uncover the existence of non-perturbative duality symmetries relating four competing orders. Two of these orders correspond to loop-current ordered phases that spontaneously break the time-reversal symmetry and describe charge currents circulating in a staggered pattern either around the plaquettes or along the diagonals of the ladder. These unconventional phases are shown to be dual to conventional (charge and bond) density-wave phases through an exact density-current duality symmetry existing on the lattice. From a perturbative renormalization group approach, we find that these phases for $N>2$ are stabilized in a half-filled SU($N$) two-leg Hubbard ladder with an additional SU($N$) Hund's interaction. The effect of a small doping on these phases is also discussed.

The graphite-to-diamond transition exhibits striking variability under high-pressure, high-temperature (HPHT) conditions, producing diamond, graphitic phases, or metastable, mixed diamond-graphite nanocomposites despite similar synthesis conditions. Existing atomistic models, largely based on idealised single-crystal graphite, do not explain the persistence of partially transformed intermediate states under HPHT conditions. Here, using large-scale molecular dynamics simulations, we show that precursor grain structure governs graphite-to-diamond transformation pathways by decoupling diamond nucleation from cooperative transformation propagation. Grain boundaries first facilitate local sp$^3$ nucleation, after which diamond growth propagates within individual grains but becomes arrested at crystallographically mismatched grain boundaries. As a result, structurally heterogeneous graphite stabilizes kinetically arrested mixed sp$^2$-sp$^3$ states, whereas large or single-crystalline domains favour cooperative bulk transformation into diamond. Our findings identify structural heterogeneity as a missing control parameter alongside pressure and temperature, reframing metastable transformation products as kinetically trapped states arising from precursor microstructure rather than thermodynamic intermediates. Precursor crystallinity therefore emerges as a practical control parameter governing graphite-to-diamond transformation pathways.

Physical Vapor Deposition (PVD) is widely utilized for the production of organic semiconductor devices due to its ability to form thin layers with exceptional properties. Although the layers in the device usually consist of two or more components, there is limited understanding about the fundamental characteristics of such multi-component vapor-deposited glasses. Here, spectroscopic ellipsometry was employed to characterize the densities, thermal stabilities and optical properties of co-vapor deposited NPD and TPD glasses across the entire range of composition. We find that co-deposited NPD and TPD form high density glasses with enhanced thermal stability. The dependences of density and stability upon substrate temperature are correlated, and the birefringence of the co-deposited glasses is determined by the reduced substrate temperature of mixtures. Additionally, we observe that the transformation of a highly stable and dense two component glass into its supercooled liquid initiates from the free surface and propagates into the bulk at constant velocity, like single component PVD glasses. All these features are consistent with the surface equilibration mechanism.

Vapor-deposited glass mixtures of organic semiconductors commonly serve as active layers in organic electronic devices, whose lifetime and performance are strongly influenced by the stability and structure of these mixed glasses. Here, we study the stability and anisotropic molecular packing of six co-deposited organic semiconductor glass mixtures with 50:50 weight ratio, by differential scanning calorimetry and spectroscopic ellipsometry. We find that all six binary systems exhibit high kinetic stability and significantly reduced enthalpy relative to the corresponding liquid-cooled glassy mixtures (ultrastable behavior), even for systems where the glass transition temperatures of the components differ by more than 90 K. Furthermore, we demonstrate that the birefringence of a co-deposited glass mixture, a measure of its anisotropic packing, can be predicted from the birefringence of glasses of the two pure components. These results for stability and structure are expected to be applicable to other co-deposited organic semiconductor glass mixtures, so long as the two components mix well in the glass and individually can form ultrastable glasses. Therefore, our findings are significant for designing novel electronic devices with enhanced device lifetime and increased operational efficiency.

X-ray diffraction (XRD) is fundamental to structural materials characterization, yet transforming a raw powder pattern into a refined crystal structure still demands considerable domain expertise. We present AGAPI-XRD, a hybrid framework integrating DiffractGPT generative structure prediction, database pattern matching against JARVIS-DFT and COD, and automated Rietveld refinement and ALIGNN-FF relaxation through a unified API at this https URL. First, we used the AGAPI-XRD pipeline to evaluate the crystal structure of a variety of minerals in the RRUFF database that were experimentally characterized using powder x-ray diffraction. Next, we benchmarked the lattice parameter prediction fidelity of the AGAPI-XRD pipeline using a subset of the Alexandria PBE-hull dataset and the subset of RRUFF minerals that have known lattice parameters. AGAPI-XRD returns valid lattice parameters for 79.7\% of the RRUFF benchmark minerals and for 94.8--98.1\% of the Alexandria subset, while identifying a candidate structure for 93.8\% of RRUFF minerals. For this benchmark, pattern matching delivers the highest accuracy for known phases, while DiffractGPT extends structure generation to complex materials absent from existing databases. Together, AGAPI-XRD advances accessible, end-to-end automated crystal structure determination from powder XRD data.

Two-dimensional van der Waals heterobilayers provide an attractive platform for the development of next-generation spintronic devices. Here, first-principles calculations are performed to investigate the structural, electronic and magnetic properties of MSe2/WTe2 (M = V, Cr, Mn, Fe, and Co) van der Waals heterobilayers. The pristine WSe2/WTe2 heterobilayer in AA'-configuration is found to be energetically favorable and exhibits type-II band alignment with a band gap of 0.70 eV, and this provides an ideal platform for controlling carrier transport. Substituting W with 3d transition metal atoms, induces long-range magnetic ordering and reconstructs the spin-resolved electronic band structure. The formation of the heterointerface generates pronounced charge redistribution and an intrinsic built-in electric field, leading to interface-induced electronic reconstruction. MnSe2/WTe2 heterobilayer exhibits half-metallicity, whereas FeSe2/WTe2 heterobilayer simultaneously exhibits half-metallicity and spin-resolved type-II band alignment. Interfacial electronic reconstruction further produces a substantial perpendicular magnetic anisotropy, driving MnSe2 from an in-plane easy axis with MAE value of 1.10 meV in the isolated monolayer to a robust out-of-plane easy axis with MAE value of 20.8 meV in the heterobilayer. Among all the structures, CoSe2/WTe2 heterobilayer exhibits maximum Curie temperature (273.87 K). The combined results establish that interface engineering makes MSe2/WTe2 heterobilayers as a promising candidates for next-generation low-dimensional spintronic applications.

The choice of crystalline orientation and offcut angle is non-trivial for low-symmetry $\beta\text{-Ga}_2\text{O}_3$, where anisotropy impacts bulk and thin film synthesis, material properties, and power device fabrication and performance. Scalable (100)-oriented $\beta\text{-Ga}_2\text{O}_3$ wafers are desirable for electronic devices but are not typically used due to 10-30x slower growth rates compared to other orientations. Here we report molecular beam epitaxy (MBE) growth rates equal to the fast growth direction by using (100) $Ga_2O_3$ wafers with large grown-in offcuts. The offcuts (up to 13.4°) are directly grown by Edge-defined Film-fed Growth (EFG) of 2D ribbons with rotated seed crystals, avoiding material loss from crystal boule offcut methods while maintaining high crystalline quality. Chemical-mechanical polishing produces epitaxy-ready substrates, and step flow growth is observed across all offcut angles. We measure an unintentional n-type doping density of $2{\times}10^{15} cm^{-3}$, one of the lowest values reported for MBE-grown films. Planar Schottky barrier diodes on these epilayers without edge termination have an on/off ratio ~10$^5$ and an average breakdown field of 1.56 MV/cm, comparable to or exceeding similar devices fabricated on other orientations. Overall, these results illustrate the importance of both crystal face and offcut angle and validate the use of the scalable (100)-oriented $\beta\text{-Ga}_2\text{O}_3$ wafers.

Water isobaric heat capacity is anomalously large under ambient conditions and exhibits a sharp maximum upon supercooling. Using classical and path-integral molecular dynamics with accurate machine-learning interatomic potentials, we show that nuclear quantum effects primarily act by suppressing high-frequency vibrations, while the anomalous temperature dependence of the isobaric heat capacity originates from structural fluctuations, quantified by the second-solvent-shell intruder order parameter. A simple two-state mapping reveals an effective enthalpy scale of about 4 kJ/mol associated with the interconversion of low- and high-density-like local structures, providing a microscopic link between their population changes and the excess heat capacity from supercooled to ambient conditions.

For a weakly interacting Bose gas trapped by an imperfect one-dimensional artificial crystal, we study the effect of its punctual defects, i.e. vacancies, on the ground state properties of the system. In the framework of the mean field approximation, we numerically solve the corresponding Gross-Pitaevskii equation using the ``Gradient Flow with Discrete Normalization'' method, also known as the imaginary time method. The crystal is artificially produced by applying an external Dirac comb potential to the Bose gas where vacancies are created by randomly removing a predetermined number of deltas. We observe that as the number of randomly removed deltas increases, the ground state energy decreases exponentially from its value for the perfect crystal case until it reaches its value when the Bose gas is free. The ground state energy is reported for different magnitudes of the interaction between bosons and several system sizes which we extrapolate to infinity for the crystal with only one vacancy. Also, we observe the presence of an energy gap between the ground state energies of the perfect system and that with a vacancy, which is more noticeable for values of the particle interaction magnitude $ g \leq 0.1$, when the delta strength $P_0 = 10$. In addition, we report the boson distributions within the crystal, %inside a box with periodic boundary conditions, i.e. the probability density functions which show localization features around vacancies which disappear as $g$ increases. From the ground state energy, the chemical potential is obtained immediately.

Materials characterized by competing interactions often exhibit a large number of nearly degenerate periodic states. Here we show that layered UAgBi$_2$ hosts a cascade of field- and temperature-induced magnetic transitions. Based on specific heat, thermal expansion, and neutron diffraction, we construct a phase diagram that reveals at least seven nearly degenerate magnetic states in UAgBi$_2$. The observed multi-step magnetization process can be understood by square-wave structures with distinct propagation vectors $\mathbf{k}$=(0, 0, \textit{k}) in the presence of strong easy-axis anisotropy that confines the moments along the \textit{c} axis. Our findings are consistent with a magnetic devil's staircase described by the axial next-nearest neighbor Ising (ANNNI) model and place UAgBi$_2$ as a rare realization of the devil's staircase in a 5\textit{f}-electron system.

We revisit the anomalous Hall effect in magnetic conductors, and its generalization to finite frequencies, using a formalism based on microscopic notions of polarization, magnetization, and free charges and currents. The electronic degrees of freedom are treated within second-quantized field theory, where the Hamiltonian features a static and cell-periodic magnetic field that encodes the magnetic order in the crystal and breaks time-reversal symmetry. We study the dynamics of bound and free charge carriers at the microscopic level as they respond to a spatially uniform electric field at finite frequency. The conductivity tensor describing the long-wavelength response is a sum of three terms, including a Kubo term associated with the polarization response, along with the metallic Drude term and the anomalous Hall conductivity that are associated with the longitudinal and transverse parts of the free current response, respectively. We also present numerical calculations of these contributions for the ferromagnetic body-centered cubic phase of iron.

The light metal Al was predicted to have strong orbital Hall and Rashba effects from its bulk and interface, respectively. In this letter, we report experimental evidence that neither the bulk nor the interface of the Al contributed a detectable torque on adjacent Co and FePt layers with significant spin Hall effect, spin-orbit coupling, and resistivity mismatch with Al. These results suggest minimal orbital-spin conversion in Al and negligible orbital transport and spin-vorticity torque in Co/Al, Al/Co, and Al/FePt bilayers. Our findings suggest poor generality and/or effectiveness of torque contributions by the orbital Hall effect, the interfacial orbital Rashba effect, and the spin-vorticity coupling.

We investigate the dynamics of a dark-bright soliton in harmonically trapped two-component Bose-Einstein condensates and reveal an interesting spontaneous spatial symmetry breaking driven by nonlinear interactions. When the interaction parameter crosses a threshold value, we find that the dark-bright soliton's motion demonstrates a transition from symmetric periodic oscillation about the origin to asymmetric oscillations offset from the origin. In particular, at the transition point, the effective soliton mass, determined by the ratio of inertial mass to physical mass, diverges. The underlying mechanism is uncovered by constructing trial wave functions and employing the Lagrangian variational method to obtain an effective potential in the quasiparticle picture, which changes from a single well to a double well. The anomalous ``superheavy soliton'' phenomenon is a direct consequence of the dark-bright soliton's physical mass vanishing at the transition point. We obtain the phase diagram of this spatial symmetry-breaking transition. Possible implications of our finding for quantum metrology are discussed.

Universal scaling behavior appears across a wide range of nonlinear systems despite substantial differences in their governing equations and physical mechanisms. We develop a renormalization-group (RG) framework that identifies two complementary RG mechanisms underlying such universality. First, scale invariance generates RG fixed points corresponding to asymptotic self-similar solutions. Second, repeated RG transformations eliminate non-scale-invariant irrelevant structures, causing broad classes of equations to flow toward the same fixed points and thereby form universality classes. The framework applies to finite-time singularities, long-time intermediate asymptotics, stochastic Edwards--Wilkinson growth, nonlinear diffusion, density-dependent biological diffusion, and fluid-interface dynamics. In each case, it reproduces known scaling behavior and identifies the associated universality class through explicit irrelevance criteria. A central feature of the framework is its predictive character. Once a scale-invariant fixed point is identified, the theory predicts entire families of nonlinear equations sharing the same asymptotic self-similar solution. While the diffusion class is partially supported by existing mathematical RG results, most universality classes identified here have not previously been established and therefore constitute falsifiable predictions. These results provide a unified RG perspective on universality in nonlinear systems and show that universality emerges from the same fundamental RG principles that underlie critical phenomena. In contrast to critical phenomena, where observable behavior is typically governed by unstable fixed points requiring fine tuning, self-similar dynamics are generally selected through dynamically stable RG fixed points.

We study the quantum properties of two- and three-dimensional spin textures -- $k\pi$-skyrmions and hopfion rings -- within the cluster mean-field approximation (CMFA). By combining the CMFA with a symmetrization procedure, we achieve two key advances: the accurate computation of quantum fluctuations in large spin textures and reliable access to metastable states. These challenges are generally insurmountable using standard methods, which are severely limited by the curse of dimensionality and typically restricted to ground-state properties. Exploiting the cylindrical symmetry of the studied magnetic configurations, we construct one-dimensional chain-like clusters that can be efficiently simulated using the density matrix renormalization group method, while inter-cluster interactions are treated at the mean-field level. The resulting spatial profiles of quantum features such as the local variation of the magnetization length in hopfion rings reveal limitations of the classical micromagnetic model and indicate the necessity of its extension. We demonstrate that the recently proposed regularized micromagnetic equation provides a suitable framework for this purpose.

Yielding, plasticity, and necking are central to the mechanical performance of materials, yet a concise unified physical picture of how these nonlinear responses arise remains lacking. We develop a minimal theory for glassy polymers based on a classical volume-dependent relaxation time following the Doolittle equation, and derive the constitutive relation using the Onsager variational principle. Surprisingly, this simple theory explains yielding, plasticity, and neck initiation under constant strain rate loading via a shear unfreezing mechanism: as the sample is stretched, volume-increasing activated molecular mobility drives shear deformation from an initially frozen state to an unfrozen state. The theory yields an analytical expression for the yielding stress as a function of strain rate and temperature. It also predicts a phase diagram for necking initiation in the same parameter space, providing a mechanism beyond the classical Considère criterion. Our results establish a unified framework for nonlinear tensile behavior in glassy materials.

Microwave kinetic inductance detectors (MKIDs) are generally read out with microwave readout tones of high enough amplitude to adequately suppress the noise contribution of the first-stage amplifier. At high readout power, the detector's resonant frequency is altered as a result of the dependence of the kinetic inductance on the internally circulating microwave current. With the tone placed below the resonant frequency, the nonlinear frequency shift results in a positive feedback effect that can significantly enhance the responsivity of the detector, to both optical and microwave power. We report a factor of 10 enhancement in optical response by tuning the readout power and frequency to close to the resonator's bifurcation point. A corresponding decrease in the bandwidth of the resonator is observed under these conditions. We show that the strength of the feedback effect can be easily selected by adjusting the excitation, and provide a map of possible operational states to do so. Operation of MKIDs in this mode could be used to improve sensitivity when non-intrinsic noise sources are significant.

The interplay between band topology and strong correlations is central to modern condensed matter physics, but exact solutions are rare. Here, we present an exactly solvable model on the Kagome lattice by combining a Kane-Mele-type spin-orbit coupling with the Hatsugai-Kohmoto interaction. At 1/3 filling, we uncover a continuous topological quantum phase transition driven by electron interaction. A weakly correlated $\mathbb{Z}_2$ topological insulator gives way to a strongly correlated insulator that, while $\mathbb{Z}_2$-trivial, hosts a nontrivial spin Chern number $C_s=2$. The transition exhibits critical scaling consistent with the universality class of two-dimensional Dirac fermions. At half-filling, the same model yields a non-Fermi-liquid to Mott-insulator transition, demonstrating that correlation-driven topological and Mott transitions can be unified within a single solvable framework. Our results establish the Kagome Hatsugai-Kohmoto model as a valuable benchmark for interacting topological systems.

At continuous electronic phase transitions, Coulomb interactions can modify the relation between length, energy, and temperature, but experimentally disentangling their effects on spatial versus dynamical criticality has remained difficult, since finite-temperature scaling alone measures only the combined exponent $\kappa = 1/(z\gamma)$. Here, we introduce two advances that resolve this limitation. First, by combining temperature scaling with independent current scaling, we separately extract the dynamical exponent $z$ and the localization-length exponent $\gamma$ at the quantum Hall plateau transition -- rather than inferring one from an assumed value of the other. Second, using dual-graphite-gated graphene devices in which the effective Coulomb interaction range is tuned geometrically by the ratio of the magnetic length $l_B$ to the graphite-gate distance $d$, we track this separation across both screened and unscreened interaction regimes within the same device platform. Temperature scaling gives $\kappa \simeq 0.21$ in the screened regime and $\kappa \simeq 0.41$ in the unscreened regime; combining this with current scaling reveals that screening changes $z$ from $\simeq 1$ in the unscreened regime to $\simeq 2$ in the screened regime. In contrast, $\gamma$ remains close to $2.4$ throughout. Our results establish that gate-controlled screening selectively modifies the interaction-dependent dynamical sector of the quantum Hall transition, leaving the localization-length exponent $\gamma$ unchanged within experimental uncertainty. More broadly, this work establishes geometric screening as a versatile tool for controlling interactions and disentangling interaction and disorder effects in correlated two-dimensional systems, including fractional quantum Hall states, moiré materials, and other strongly localized electronic phases.

Continuous time crystals (CTCs) represent a nonequilibrium quantum phase that spontaneously breaks time-translation symmetry without periodic external driving, manifesting as persistent, long-lived oscillations under steady pumping. The lifetime is constrained by the instability of the limit cycle phase, balanced between nonlinear feedback and energy dissipation, which have rarely been studied in experiments before. Here, we report an observation of an ultralong-lived Rydberg-atom CTC in a driven-dissipative many-body atomic system. By harnessing long-range interactions and engineering a dissipative environment that stabilizes the limit-cycle dynamics, we suppress heating and decay effects that typically destroy time-crystalline order. The key factor underlying the ultralong-lived CTC is the closing of the Liouvillian gap and the near-zero real part of the system's Liouvillian eigenspectrum. Through systematic optimization, we achieve an oscillatory lifetime exceeding 16.95 hours-orders of magnitude longer than previous CTC realizations. Our work establishes a robust platform for exploring long-lived autonomous nonequilibrium phases and paves the way for applications in quantum sensing and continuous-time quantum information processing.

We study terahertz field-driven second-order nonlinear electron transport phenomena, including the nonlinear Hall effect (NLHE), in two-dimensional tellurene flakes. The dc current excited by linearly polarized terahertz (THz) radiation in Hall bar samples is investigated in directions both along and perpendicular to the $c$-axis of tellurene. As expected for second-order transport phenomena, the current scales as the square of the in-plane electric field of the radiation $\bf E$, and depends on its orientation. The current results from a combination of three contributions, including the NLHE, the Nonlinear Longitudinal (NLL) and Nonlinear Diagonal (NLD) currents. We established the equivalence between NLH, NLL, and NLD transport currents and Linear photogalvanic effect (LPGE) contributions induced by the absorption of linearly polarized and unpolarized THz radiation. All contributions can be controlled by a gate voltage and have opposite signs for electron and hole conductivity. The magnitude of the current increases drastically when the samples are cooled from room temperature to 4.2 K. It also increases with decreasing radiation frequency. These results are well described by the developed phenomenological and microscopic theories. We show that the THz radiation-induced electric current originates from microscopic mechanisms such as skew scattering, side jump, and the Berry curvature dipole.

The classical dichotomy between crystalline order and amorphous disorder is increasingly challenged by novel states that lack conventional crystalline symmetries while retaining crystal-like properties. Quasicrystals occupy a distinctive position within this expanding framework by possessing long-range order without translational periodicity, thereby permitting arbitrary $N$-fold rotational symmetry. Paradoxically, far from their unique symmetry center, high-symmetry quasicrystals closely resemble disordered patterns, raising the question of how deterministic order can be detected. Here we show that increasing rotational symmetry progressively suppresses local statistical signatures of quasiperiodicity, while preserving its underlying exact long-range order. This order is thus concealed below an emergent crossover length that grows linearly with $N$. Therefore, as $N \rightarrow \infty$, the disorder-like regime expands without bound, defining a symmetry-controlled geometric critical point at which deterministic order and randomness become statistically indistinguishable over any finite observation window. For finite $N$, however, quasiperiodic order becomes detectable beyond this crossover, revealing a second emergent length scale that we identify as the size of a \textit{statistical unit cell} -- finite patches over which statistical properties recur despite the absence of conventional translational periodicity. In one dimension, the statistical-unit-cell size coincides with the crossover length, whereas in two dimensions it grows as $N^2$, remaining smaller than the size of typical approximants and establishing a hierarchy of emergent length scales. Together, the disorder-to-order crossover and statistical unit cells provide a quantitative framework connecting crystals, quasicrystals, and amorphous matter, showing how apparent disorder can emerge from purely deterministic geometry.

The attention matrix of a causal transformer is row-stochastic, iterated over depth, and non-normal by construction. For non-normal operators, eigenvalues control only asymptotic behavior; finite-depth behavior is controlled by resolvent quantities such as pseudospectra and Kreiss constants. We test, under pre-registered criteria, whether this resolvent view predicts anything about trained transformers that eigenvalues miss. Two structural facts organize the analysis: the mask pins the Kreiss constant of every causal stochastic matrix at $\sqrt{n}$, and deflating the mask-forced Perron projector factorizes the depth deviation dynamics exactly into a product of deflated operators. Across GPT-2, Pythia-410m, and Llama-3-8B, learned non-normality proves to be signed. A routing minority carries excess transient reserve that tracks previous-token function and doubles when induction heads engage, while the sink majority is suppressed below matched shuffle nulls, so that attention sinks act as transient dampers. On depth products, eigenvalue predictions of surviving deviations err by seven to eleven orders of magnitude, an error absent in matched nulls. Checkpoint censuses date this organization to a consolidation phase after circuit formation, and a clamping intervention on Llama-3-8B establishes a causal chain from three massive activation dimensions through sink attention to transient damping; LayerNorm models implement the same functions elsewhere. A cross-validated contest concludes that resolvent features are required for depth-transient persistence and routing-head identity, and that no single-operator summary of any kind predicts per-head causal criticality.

Aqueous corrosion of metals involves multiple interconnected processes. Thus, computer simulation of corrosion as a whole needs to be able to describe the individual processes and how they influence each other. Atomistic simulations are designed to obtain detailed information for small regions of space over short times. Thus there are limits to the understanding that can be obtained by atomistic simulations alone. Here is presented a method that uses generalised rate equations to extend the length and time scales that can be accessed. It is shown to reduce to either the phase field or cellular automata methods in certain limits. The generalised kinetic equations can reproduce the behaviour described by both the Nernst-Planck and Butler-Volmer equations, which are frequently used to describe corrosion. In addition, the method can describe local rearrangements of atoms such as chemical reactions. Example results are shown for illustrative 1D and 2D problems, with good agreement being found with other methods.

Semiempirical electronic structure methods such as Density-Functional Tight-Binding (DFTB) offer a computationally efficient approach to molecular and materials simulations, bridging the gap between first-principles accuracy and classical force field speed while retaining full access to electronic properties. However, DFTB calculations based on self-consistent charge (SCC) schemes can still suffer from slow convergence, particularly for complex molecular and materials systems, making the iterative procedure a significant bottleneck in large-scale simulations and high-throughput workflows. We present a machine learning approach that accelerates DFTB simulations by predicting optimal initial atomic charges. Using element-specific models based on the Smooth Overlap of Atomic Positions descriptor and kernel ridge regression, we train charge models on reference calculations and demonstrate that ML-predicted initial charges consistently and significantly improve SCC convergence across diverse chemical systems including organic molecules, biomolecules, water clusters, transition metal oxides and solid electrolytes.

A semi-empirical kinetic framework is formulated for predicting phase selection in multicomponent concentrated alloys under rapid solidification. The approach is based on the critical cooling rate required to suppress competing crystalline pathways and combines topology-dependent ranking of BCC-, FCC-, and HCP-like crystallization pathways with a correction for glass-forming ability. The formulation includes a topology-dependent viscosity correction with a smoothed BCC multiplier and a continuous correction factor for glass-forming ability based on mixing enthalpy, excess entropy, and atomic-size dispersion. Comparison with experimental and computational data shows that the kinetic criterion captures changes in the lattice type expected from the valence electron concentration criterion, describes kinetic suppression of phase separation, and identifies competitive multiphase crystallization. The model also distinguishes alloys with high and low glass-forming ability. The proposed framework provides a practical approach for preliminary evaluation of kinetic phase competition in rapidly solidified multicomponent melts.

A common view in physical chemistry literature is that the Eyring representation, in which a linear fit of $\ln(k/T)$ versus $1/T$ is attempted, is more fundamental than the Arrhenius representation, $\ln(k)$ versus $1/T$. This perception is typically motivated by its derivation from statistical mechanics and quantum mechanics, and by the interpretation of the intercept in terms of the activation entropy $\Delta S^\ddagger$, whereas the Arrhenius equation and its prefactor are often regarded as purely phenomenological. However, harmonic approximation models yield exact linearity in Arrhenius plots but not in Eyring plots, although for real experimental data both generally appear equally linear within typical experimental accuracy. Furthermore, the impression that the Eyring formulation is inherently quantum mechanical arises from the presence of the Planck constant in the prefactor, whereas this term results from normalization conventions in the partition function. This also highlights an interpretational issue in $\Delta G^\ddagger$, which is based on partition functions of different dimensionality between reactant and transition state. This dimensional mismatch can be reformulated in an alternative representation that improves interpretability and reduces to an Arrhenius-type expression in which the prefactor is directly related to an entropy of activation. In this framework, both activation enthalpy and entropy obtained from an Arrhenius fit are arguably more physically relevant than the corresponding Eyring fit values.

Photoacoustic (PA) signal generation in colloidal suspensions of optically absorbing nanoparticles is dominated by the thermal expansion of water for gold nanoparticles, but remains mostly unexplored for organic nanoparticles. Here, we derive a model where the PA generation efficiency scales with particle size and thermoelastic contrast with water. The model is validated using solid lipid nanoparticles labeled with several BODIPY dyes. This experimental validation paves the way for quantitative PA characterization of nanomaterials and rational design of PA contrast agents.

Polymerization-induced phase separation (PIPS) provides a powerful route to generate structured polymeric materials by coupling chemical conversion with thermodynamic demixing. PIPS in liquid-state systems underlies dispersion polymerization, serving as a cornerstone technique for microparticle production, yet is constrained by solvent compatibility and limited range of morphologies. Here, we establish an elastically mediated PIPS regime that bridges these two limits by conducting controlled polymerization within a deformable elastomeric network. This approach, termed Dispersion Polymerization in an Elastomeric Solvent (DiPolES), serves as a solid-state analogue of dispersion polymerization in which an elastomeric network simultaneously serves as solvent and physical stabilizer. Using photoiniferter-mediated polymerization of methyl methacrylate (MMA) within poly(dimethyl siloxane) (PDMS) elastomeric solvent, DiPolES enables robust fabrication of elastomeric composites containing uniform PMMA microparticles with tunable size (0.85 to 3 {\mu}m) and shape (spheroidal and ellipsoidal). The strategy is generalizable beyond the PDMS/MMA system and is applicable to diverse monomers, such as acrylonitrile and 2-vinyl pyridine, which can be extracted from the elastomeric solvent, enabling high-yield production of microparticles. Real-time imaging and compositional analysis reveal that particle formation proceeds through rapid nucleation at low monomer conversion, followed by growth accompanied by cavitation of the surrounding network. Monomer loading governs the particle size, while solvent elasticity modulates the transition from isolated uniform spheroids to heterogeneous clusters. Interestingly, applying uniaxial strain during DiPolES enables production of ellipsoidal particles without any post-processing.

The conventional multipole framework provides a standard symmetry-based language for one-body electronic degrees of freedom, but fails to distinguish physically distinct sectors in many-body operator space. We show that canonical symmetries compatible with rotations provide this structure through two additional 2 labels associated with the particle-number gauge transformation G{\pi}/2 and the particle-hole transformation CA. These labels separate operators with different body numbers and particle-number changes, leading to a sixteenfold classification that organizes selection rules for nonzero expectation values, induced multipoles, and symmetry-allowed couplings in many-body multipole space.

The finite-rate cooling dynamics of the triangular-lattice $J_1-J_2$ Ising antiferromagnet is studied under local Metropolis updates. Although an antiferromagnetic next-nearest-neighbor coupling selects a stripe phase in equilibrium, the simulations show that this phase is not automatically reached on finite time scales. A kinetic stripe-formation time $n^*(L,J_2/J_1)$ is defined from the probability of obtaining a globally stripe-ordered final state. This time shifts to much slower cooling as the system size increases and to faster cooling as $J_2/J_1$ increases. The size dependence is compatible with a coarsening-controlled process, with an effective growth at least quadratic in L over the simulated range. Real-space morphology and fixed-temperature diagnostics show that failed trajectories are not simply disordered states: they often contain locally stripe-ordered domains separated by residual walls or competing orientations. In the weak-$J_2/J_1$ regime, the system can restore the local nearest-neighbor frustrated constraint while still failing to select a global stripe sector. These results separate three processes that are usually conflated: energetic degeneracy lifting by $J_2/J_1$, local constraint restoration, and global stripe-orientation selection under local dynamics.

The speed of optically driven magnetization dynamics is fundamentally determined by how efficiently angular momentum can be transferred between electronic, spin and lattice degrees of freedom. In rare-earth magnets, this process is typically slow because optical excitation primarily addresses itinerant electrons, whereas the magnetic moment resides in localized 4f states. Here we show that selective optical excitation of localized magnetic states can overcome this limitation. Using femtosecond pump-probe magneto-optical spectroscopy of ferrimagnetic gadolinium iron garnet, we resonantly excite an intra-4f transition of Gd3+ at 4.65 eV and resolve the ensuing dynamics of the antiferromagnetically coupled Gd and Fe sublattices. Direct excitation of the 4f manifold induces an ultrafast demagnetization of the Gd sublattice with a characteristic time of 38 fs, more than two orders of magnitude faster than in elemental gadolinium and even faster than the response of the Fe sublattice in the same material. By contrast, off-resonant excitation strongly suppresses the acceleration of the Gd dynamics while leaving the Fe response largely unchanged. These results demonstrate that the ultrafast magnetic response of rare-earth systems is governed not only by intrinsic material properties but also by the optical excitation pathway. Selective access to localized magnetic states therefore provides a powerful photonic handle for engineering angular-momentum flow and controlling magnetism far from equilibrium.

We develop a statistical-mechanical theory of water in unsaturated soil whose outcome is a continuum field equation for the pore-occupancy g(r,x,t), the fraction of pores of radius r that are water-filled at position x and time t. The theory is built across three scales: microscopic inter-pore transfers set by Hagen-Poiseuille rates and a driving potential (the difference of pore-class chemical potentials, taken in capillary-gravitational form but open to adsorptive, osmotic, or thermal refinement); a mesoscale master equation relaxing the occupancy toward the equilibrium step g_eq=H(r*-r); and, on contracting the averaging volume to a point, the continuum balance d_t g + div F = C[g] - E - T, of which everything else is a limit, a moment, or a boundary resolution. The kinetic equation is an Onsager gradient flow descending the Gibbs free energy, with an H-theorem for the isothermal unforced system and mass conservation as its zeroth moment. A single dimensionless group, the pore-resolved Damkohler number Da(r,x), organizes the behavior and unifies phenomenologies long modelled separately. A Chapman-Enskog reduction identifies Richards' equation as the quasi-static (Da->0) limit, with matric potential and hydraulic conductivity K emerging only there and K vanishing below the percolation threshold; capillary-bundle and critical-path models are its diagonal and spectral limits. Hysteresis is the holonomy of a forcing bundle, a geometric phase rather than per-pore bistability, with a falsifiable loop-area law H ~ I^2. Preferential flow is what the same equation does where Da>1, so the Richards/preferential-flow dichotomy becomes a continuous Da-controlled crossover. Out of the quasi-static limit g(r) is the irreducible state variable. All inputs are geometric properties of the pore network, measurable from micro-CT and calibrated against no macroscopic data.

We measured coherent ultrafast dynamics of exciton complexes in a single strongly-confined InAs quantum dot embedded in a GaAs nanopost. Such a photonic structure combines a wave guiding with a cavity effect and assures an enhanced light-matter coupling. Coherence properties of an exciton-biexciton system hosted by a quantum dot are assessed with four-wave mixing microscopy. Our results show that this broad-band photonic structure is an excellent asset to probe coherent couplings in a small set of solid state quantum systems and to investigate the coherence dynamics within the level structure of their excited states.

We investigate the evolution of spin correlations and the nature of the multi-stage magnetic transition in the quasi-two-dimensional easy-plane van der Waals magnet $\rm CrCl_3$. By combining broadband ferromagnetic resonance (FMR) spectroscopy and DC SQUID magnetometry on mechanically exfoliated micro-flakes with non-local dynamical mean-field theory (DFT+DMFT) calculations, we analyze both long- and short range magnetic order in CrCl$_3$. Experimentally, SQUID and FMR measurements confirm the presence of the crossover to a spin polarized phase with the subsequent transition into an antiferromagnetic ground state upon cooling, but show robust short-range correlations at temperatures far above the magnetic ordering temperatures. Theoretically, we show the existence of highly stable local magnetic moments at room temperature, with a giant room temperature lifetime $\tau$ of 130--300 ps due to a wide Mott bandgap. Below room temperature, a rapid growth of the in-plane correlation length $\xi$ signals the formation of strong short range magnetic order consistent with the experimental observations. We also obtain a temperature-driven crossover of the interlayer exchange interaction, which changes from positive (ferromagnetic) at high temperatures to negative (antiferromagnetic) in the low-temperature ordered phase.

When interaction among atoms or ions is strong enough, they often arrange periodically, forming a crystal. The arrangement patterns of atoms or ions can encode information, a concept that has enabled devices such as ferroelectric memories. It has been found that not only atoms or ions but also electrons in condensed matter can crystallize when Coulomb interaction is strong enough. Typical examples are charge-ordered states in solids, where different valences, or different electron numbers, of an ion spontaneously form a spatial pattern on the lattice. In such electron crystals, information is expected to be encoded into the electron-ordering patterns. Here, we demonstrate electronic manipulation and readout of charge-ordering directions in a paramagnetic semiconductor LuFe$_2$O$_4$. By applying current pulses at room temperature, we observed that the non-reciprocal resistivity of LuFe$_2$O$_4$ is modulated along with a sign reversal, which disappears above the charge-ordering temperature. A numerical calculation incorporating inter-band Berry curvature affected by the charge ordering is consistent with the experimental results. By applying the observed phenomenon, we also demonstrate a non-reciprocal resistance memory operation in the charge-ordered LuFe$_2$O$_4$. This result opens the door to realizing charge-ordering electronics.

The low-energy dynamics of many-body systems is governed by gapless modes whose properties are dictated by symmetry. Their existence follows from Goldstone's theorem, while their effective description at zero temperature is determined by the pattern of symmetry breaking. At finite temperature, an analogous role is played by hydrodynamics, which describes the universal behavior of many-body systems over long times and large distances. These principles are well understood for uniform symmetries, which act homogeneously in spacetime and lead to a correspondence between massless Goldstone modes and broken generators, as well as gapless hydrodynamic modes and conserved charges. However, this simple picture changes in the presence of nonuniform symmetries, whose generators do not commute with spacetime translations. The low-energy implications of these symmetries remain less understood, as they do not introduce additional gapless modes but instead constrain the dynamics of the existing degrees of freedom. In this thesis, we develop a unified framework for many-body systems with nonuniform symmetries and show that their effects can be understood in terms of additional fields that are not independent at low energies and can be eliminated, leading to kinematic constraints that reshape the infrared dynamics. In systems with spontaneous symmetry breaking, this mechanism modifies the effective theory and often softens the dispersion relations of the remaining modes. At finite temperature, it manifests in hydrodynamics as constraints on macroscopic currents. As a result, nonuniform symmetries give rise to qualitatively new physical phenomena, including modified spectra of collective excitations, exemplified by transverse Tkachenko oscillations in quantum vortex crystals, and unconventional transport phenomena, such as anomalously slow diffusion and softened sound modes in dipole-conserving systems.

The interplay of tangential activity and excluded-volume interactions in ring polymers adsorbed to a surface, consistently results in the overall swelling of the ring configurations. This is in strong contrast to the case for three dimensional linear and ring polymers, where activity induces frequently collapsed structures. By means of Brownian Dynamic simulations, we investigate how the scaling properties of such active rings can be universally characterized by an activity-dependent Flory exponent, as a generalization of the equilibrium behavior. At high activity, an effective persistence length characterizes the conformations of active flexible rings.

The superconductivity in RuO2 emerges under strain. RuO2 is also an altermagnet candidate. The nature of superconductivity and its relation with neighboring orders, however, are not understood. To address this problem, we grew epitaxial RuO2 films on TiO2(100) and TiO2(110) single crystal substrates and studied the electronic transport and emergent superconductivity along various crystallographic directions. We show that the superconducting transition strongly depends on the growth orientation and the crystallographic direction of the transport in the RuO2 films. We also observe a strong violation of Pauli paramagnetic limit with in-plane applied magnetic field which we attribute to strong spin-orbit scattering. These results offer opportunities for epitaxially engineered superconductors.

The Rosenzweig--Porter (RP) random matrix ensemble has emerged as a minimal model for the integrability-to-chaos crossover in quantum many-body systems. Its phase diagram features a region with fractal eigenstates, exhibiting intermediate spectral and localization properties between the fully localized and fully delocalized regimes. In this work, we explore several generalizations of the RP model and determine their level statistics at the scale of the Thouless energy $E_T$, which characterizes the crossover. Using tools from free probability theory and the replica method, we compute the full counting statistics in the limit of large system size, and show that it takes a simple, universal scaling form around $E_T$, shared across all variations of the model. We validate our analytical predictions using exact numerical diagonalization of large samples, and large-deviation algorithms that resolve the full counting statistics down to probabilities as low as $10^{-40}$. We also contrast our predictions with measurements on the quantum random energy model, which is the simplest model displaying many-body localization.

I critically overview my research on strongly correlated fermion systems for almost five decades. It concentrated on: (i) the first derivation of what is now called the t--J model, comprising both the limit of Anderson kinetic exchange of spin--spin interaction in the Mott--Hubbard insulator and taking into account real-space pairing, subsequently applied to high-temperature superconductivity; (ii) the concept of spin-dependent heavy mass of quasiparticles in heavy-fermion systems, and (iii) the first nontrivial model of statistical thermodynamics of the Mott--Hubbard transition. Those three features, together with the specific quantum critical phenomena provide, in my view, fundamental components of the theory of strongly correlated fermions established in the 1960s. Some related questions such as introduction of atomicity in the chemical bonding (iv), and specific properties of correlated nanosystems within the rigorous EDABI (v) approach are also briefly elaborated at the end.

Despite advancements in electron optics and spectrometer design over the past twenty years, atomic-resolution valence-electron energy-loss spectroscopy imaging remains challenging due to the delocalization of inelastic electron scattering. In this study, we used an energy-filtered spectrometer equipped with a hybrid-pixel direct electron detector and spherical aberration-corrected scanning transmission electron microscopy to analyze many-electron excitations and interband transitions in a single-crystal calcium niobate, CaNb2O6, with spatial resolution ranging from the nanometers to the atomic scale. In the low-loss region above the bandgap at about 3.8 eV, we observed volume plasmons, around 6 eV and 15 eV energy loss, as well as a mix of strongly correlated plasmons and excitons, known as plexcitons, at approximately 7.3 eV energy loss. Additionally, we employed an on-axis EELS setup for atomic-resolution zero-loss peak (ZLP) imaging and visualized energy- and atom-resolved images of plexcitons and VPs, which showed contrast reversal relative to high-angle annular dark-field images. To investigate elastic contrast preservation, we also analyzed the effect of the collection angle and minimized its influence to produce delocalized VP images. In fact, the ZLP and VEELS images obtained using the weak-beam setup demonstrate that, in both cases, the contrast resembles Z-contrast. Moreover, we found that [NbO6] octahedra directly contributed to the lateral maps of interband transitions in the range from 3.2 eV to 3.5 eV energy loss. These findings demonstrate that Cs-STEM-EELS, which examines atomic-scale contrast associated with low-energy losses, can be a powerful tool for visualizing the structure, bonding, and electronic properties of complex crystalline nanostructures, including individual atomic sites, interstitial sites, and point defects.

Classical spin liquids are correlated paramagnets in which local constraints generate extensive degeneracy and emergent gauge structures, often observable as pinch-point singularities in spin structure factors. Here we introduce pinch-curve spin liquids, in which the pinch singularities form one-dimensional algebraic curves in momentum space. Inversion symmetry protects these curves by reducing the singularity condition to two real algebraic constraints in three dimensions, and the geometry of the pinch locus is algebraically programmable. We identify elementary mechanisms for generating straight and curved pinch loci, construct lattice spin models that realize them, and test the predicted structure factors using Monte Carlo simulations. We further show that pinch curves can host an infrared Gauss-law transition: the leading local constraint and the associated anisotropic scaling of the structure factor change, even though the singular locus remains one-dimensional.

Chaotic many-body spectra are expected to densely fill their energy window. We show that constrained spin chains with chiral symmetry evade this expectation by hosting an exponentially large manifold of symmetry-protected exact zero modes separated from the surrounding spectrum by a sharp gap at zero energy. The gap is generated by chaotic level repulsion, with width set by the number of zero modes times the mean level spacing. We verify this mechanism in an East-West kinetically constrained chain, develop a minimal random-matrix description, and show how the gap can be detected through linear-response spectroscopy.

Majorana zero modes provide a route to fault-tolerant qubits by encoding information non-locally in fermion parity. Their sensitivity to noise is expected to decrease exponentially with increasing separation between the Majoranas, a suppression known as topological protection. Kitaev chains engineered in quantum dot-superconductor arrays provide a tunable platform in which separated Majorana zero modes can emerge at the ends of the chain, even in two-site chains. These minimal-chain modes are known as poor man's Majoranas and retain characteristic Majorana properties, including near-zero energy and equal electron-hole character, but have only limited protection. A key outstanding challenge is to move beyond identifying such modes in electrical transport measurements and achieve coherent qubit control in the time domain. Here, we demonstrate a Majorana parity qubit by realizing coherent coupling between two-site Kitaev chains. Since total fermion parity is conserved, the system separates into global even and odd parity manifolds. We observe coherent parity oscillations in both manifolds with equal oscillation frequencies at the Majorana sweet spot, as predicted for isolated Majorana zero modes. We further show that the oscillation frequency and coherence depend systematically on inter-chain coupling and quantum-dot detunings, in close agreement with our model for short, partially protected chains. Our results establish the first coherent control of a Majorana qubit, encoded in the fermion parity of Majorana zero modes in minimal Kitaev chains.

Understanding the dynamical properties of coupled phase oscillator systems with heterogeneous oscillator frequencies has been a long-standing challenge of complex systems theory. While the seminal work of Ott and Antonsen dramatically improved our theoretical understanding of coupled phase oscillators for a small family of oscillator frequency distributions, we here present a mean-field reduction method for arbitrary frequency distributions. Our method leverages the drastic dimensionality reduction obtained for Lorentzian frequency distributions, and combines it with a data-driven multi-ensemble approach. As such, the method renders the Ott-Antonsen equations directly applicable to empirical distributions of phase oscillator frequencies, often achieving a drastic dimensionality reduction and allowing to study real-world physical and biological systems by means of stability, sensitivity, and bifurcation analyses.

We discuss two different cases of strongly correlated fermions statistics. The first of them is the non-Fermi liquid (NFL) case, $i.e.$, that of fermions with exclusion of doubly occupancy of quasimomentum states $\{\bf k\}$ with opposite spins ($\uparrow\downarrow$). The second is that of statistical spin liquid (SSL) case, in which the fermion spins hop around and mix with the holes (unoccupied) states. For both cases we calculate the system entropy and the corresponding statistical distribution function, analyzed for a two-dimensional square-lattice filling $n\in [0,1]$ and relative temperature $k_B T/|t|$, where $t<0$ is the nearest neighbor hopping integral. We are particularly interested in the situation when the system of itinerant fermions reduce to the Mott-insulator state for the half-filled band ($n\to1$). This limiting situation signals a qualitative difference between the present SSL statistics and that of uncorrelated fermions representing normal Fermi-liquid state.

Salinity gradient energy offers a sustainable route to convert ionic chemical potential differences into usable power. Capacitive mixing enables this conversion without membranes, but suffers from limited ion regulation at electrode interfaces. Here we show that by grafting electrode surfaces with titrating polymers, the performance can be substantially improved. Using Grand Canonical Monte Carlo simulations with exact image-charge Ewald summations, we demonstrate how the coupled effects of ion adsorption and charge regulation in response to an external potential can be harnessed. Grafted electrodes are shown to deliver substantially more energy relative to bare surfaces, driven by charge regulation effects that exploit the pH difference that typically exists between rivers and the ocean. While the effect is in principle maximized at high grafting densities and moderate chain lengths, the performance is fairly robust to variations of these parameters, within reasonable bounds. Complementary classical polymer Density Functional Theory calculations confirm these trends, validating the mechanistic framework. This work also establishes a practical approach to harvest electrical energy during wastewater neutralization, where acidic (or alkaline) effluents serve as complementary reservoirs, and offers a promising strategy to couple environmental remediation with renewable energy recovery.

Long-chain branching profoundly modifies the nonlinear oscillatory response of entangled polymer melts by introducing arm-retraction pathways absent in linear polymers. We present a molecular tube theory that explains the characteristic maximum of the Nonlinearity Index (NLI) observed experimentally in long-chain-branched polymers. The theory extends the recently developed nonlinear tube-orientation description of linear polymers by incorporating branch-point force transmission and delayed arm retraction. The backbone initially develops nonlinear orientation as in the corresponding linear polymer, whereas long-arm retraction subsequently relaxes the stored branch-point tension and progressively erases backbone orientational memory. This competition produces a characteristic NLI maximum followed by a post-peak decay. The theory predicts two distinct nonlinear regimes corresponding to sparse and dense long-chain branching and introduces an architecture parameter governing the height and width of the nonlinear peak. The resulting framework provides a molecular interpretation of nonlinear Fourier rheology and directly links the nonlinear harmonic response to polymer architecture.

Controlling the energy scales of a quantum system is essential for defining robust qubits. In silicon spin qubits, the nearly degenerate conduction-band valleys create a leakage channel from the single-spin computational basis, posing a challenge to scaling and to shuttling-based architectures. Here, we measure the relevant energy scales of single-electron spin qubits in buried silicon quantum wells co-designed for low disorder and high valley splitting. Across a linear array of four quantum dots with an average orbital energy of 2.4(2) meV, we report an average single-electron valley splitting of 0.40(6) meV and an average two-electron singlet-triplet splitting of 0.24(7) meV. In three dots, we observe a strong correlation between valley splitting and orbital energy, with an average linear coefficient of $\approx 0.22$ (meV/meV), demonstrating that electrostatic confinement can increase the valley splitting by several hundred microelectronvolts. In contrast, the remaining dot exhibits the highest valley splitting of 0.76(2) meV and low correlation, suggesting excellent characteristics for spin-qubit operation. Our findings demonstrate that strong confinement can be exploited in buried quantum wells to effectively enhance the valley splitting, thereby establishing a viable path toward the realization of shuttling and sparse-occupation-based architectures in low-disorder heterostructures.

We present the H-cloud formalism, a cloud-accessible and open-source finite-element framework for electromagnetic modelling of applied superconductors. The proposed method expresses the nonlinear electromagnetic response of type-II superconductors in a curl-conforming discretisation based on Nédélec finite elements, where the tangential applied-field boundary condition, nonlinear E-J power law, and fully implicit time-discrete residual are stated explicitly at the variational level, all within a scripted Python finite-element workflow. The weak form is used as the basis for forward simulation and for extension to adjoint differentiation and PDE-constrained optimisation, while keeping the governing assumptions, boundary conditions, and solver structure fully visible to the user. The implementation is realised in Firedrake with UFL and PETSc-backed nonlinear solvers, allowing the identical script to run on local machines and in browser-accessible environments such as Google Colab without reformulating the problem. The method is verified on the canonical magnetisation benchmark of a cylindrical superconductor under Bean-like penetration conditions and then benchmarked against an independently constructed COMSOL model for a practical high temperature superconducting Bi2212 wire. Across matched mesh studies, the open-source workflow reproduces the commercial-reference magnetisation loops to within approximately \(1\%\) , with relative peak errors below 1.5%, while cloud execution preserves the same numerical solution at rather modest additional runtime considering the use of (freely available) reduced hardware resources. The proposed framework provides a rigorous, reproducible, and portable route for superconducting simulation, benchmarking, and future optimisation-led modelling of applied and functional superconductors, shareable and executable into open cloud environments.

Strongly constrained quantum systems, in which local rules forbid most configurations, play a central role in condensed matter and lattice gauge theory. Their thermalization is often thought to be delicate: extensive conservation laws and dynamically frozen states can shatter the Hilbert space into many disconnected sectors. A natural question is whether, once the frozen states are removed, the dynamics within a single sector still thermalizes. We address this in the height-conserving quantum dimer model on the square lattice, whose local plaquette flips conserve an emergent height field. Resolving the winding numbers, the four sublattice heights, and lattice momentum , we isolate the dominant connected Krylov component of each fragmented sector and analyze its spectral spectral statistics, entanglement, and connectivity. The two standard chaos diagnostics then show different behavior:across momentum sectors the level-spacing statistics range from near-Poisoon to Wigner-Dyson, yet in every sector the eigenstate entanglement entropy collapses onto a narrow, dome-shaped curve characteristic of eigenstate thermalization. Only a handful of low-entanglement outliers interrupt this thermal pattern, in selected sectors. Thus, strong kinematic constraints can lead to a situation where spectral correlations and eigenstate thermalization need not follow the same universal signatures -- a manifestation of constrained quantum chaos.

We expand on the ideas developed by C. Jarzynski in Physica A 552, 122077 (2020), where an integral fluctuation theorem was derived with the aim of obtaining thermodynamic inequalities stronger than those implied by the Jarzynski equality. Restricting ourselves to the quantum setting, we derive the corresponding detailed fluctuation theorem and additional detailed and integral fluctuation theorems; we also provide a clear physical interpretation of the stochastic quantities defined in the previous reference. Furthermore, we show that their averages are given by the nonadiabaticity parameter (i.e., the relative entropy between the final state after a finite-time driving protocol and the corresponding adiabatically evolved state) and the excess work (also known as inner friction). We elaborate on the inequalities derived from the fluctuation theorems and discuss their connection to irreversibility and formulations of the Second Law.

Achieving valley splittings well in excess of the thermal energy of electrons and avoiding valley excitations is essential for the consistent initialization, operation and readout of gate-defined Si spin qubits. In this work, we present a device-level optimization strategy for pushing valley splittings to between 1 and 5 meV, well beyond values reported in nearly all previous theoretical studies. Using device-scale simulations that incorporate atomistic alloy disorder through a 1D tight-binding theory, we demonstrate that our proposed approach yields large valley splittings with a tight distribution across disorder realizations, a key requirement for reproducible qubit performance at scale. The approach rests on an unorthodox Si/SiGe heterostructure design combining a narrow quantum well, a small Ge spike, and a pure-Ge cap. We corroborate these predictions with targeted atomistic density functional theory calculations. These results offer a clear path forward for scalable Si/SiGe spin qubit devices and, if realized experimentally, effectively eliminate valley splitting as an existential problem for large scale SiGe-based quantum processors.

Non-local magic has recently emerged as a fundamental resource for characterizing genuinely non-local non-stabilizer correlations. However, its direct calculation is an intractable numerical problem, except for small systems, and its understanding remains limited. We derive a representation of non-local magic in terms of the Walsh--Hadamard autocorrelations of the entanglement spectrum. Our representation makes the underlying harmonic structure explicit and enables a systematic analysis of its behaviors for various scenarios. We prove that non-local magic can be upper-bounded by an entanglement entropy and we derive exact analytical results for broad classes of quantum states, characterizing the scaling of non-local magic for volume-law states, as well as ground states of one-dimensional gapped and critical systems. Our results identify the spectral organization of the entanglement spectrum as the key ingredient governing non-local magic and provide a framework for further systematic analytical investigation.

By integrating the quantum Monte Carlo technique into the similarity renormalization group (SRG), we have developed a stochastic SRG framework (SRGQMC) capable of both free-space two-body and in-medium many-body evolutions. This approach circumvents the combinatorial tensor-space explosion of many-body flow equations by mapping continuous unitary transformations onto an ensemble of signed random walkers. We benchmark the SRGQMC against deterministic free-space SRG evolutions of realistic nucleon-nucleon (NN) interactions, as well as against in-medium SRG (IMSRG) many-body calculations with the Richardson pairing model at two- and three-body levels [IMSRG(2)/(3)]. While a deterministic extension to the four-body level [IMSRG(4)] remains unfeasible due to prohibitive computational costs, we have achieved the first IMSRG(4) calculation by using the stochastic technique, demonstrating a substantial improvement toward the full configuration-interaction limit. This stochastic framework provides a practical pathway to higher-order IMSRG calculations.

Molecule-surface interactions are central to many research and technological areas, spanning from heterogeneous catalysis and polymer science to electrochemistry. Of particular relevance are metallic nanogaps used in molecular electronics and near-field spectroscopy. Due to the buried nature of these double interfaces, few methods exist to monitor side-specific interactions and relative molecular orientation inside the gap. In this work, we introduce plasmon-enhanced nonlinear vibrational spectroscopy as an efficient tool to investigate surface molecular adsorption within metallic nanojunctions. By exploiting simultaneous vibrational sum- and difference-frequency generation in dual-resonant nanocavities, we resolve molecular orientation and preferential binding to one of the two gold surfaces, with few-molecule sensitivity. We also discover that the non-resonant (electronic) second-order nonlinear response is not an intrinsic property of the metal surface, but is instead governed by the molecule-surface interaction. Our findings provide a powerful analytical tool, easily implementable as an add-on to Raman spectroscopy, thanks to commercially available mid-infrared quantum cascade lasers.

Interfacial charge-transfer processes play a critical role in plasmon-enhanced spectroscopy, yet the energetic conditions governing charge-transfer resonance within molecule-metal nanocavities remain poorly understood. Here, plasmonic nanoparticle-on-mirror junctions incorporating systematically engineered biphenylthiol derivatives monolayers were used to investigate how frontier orbital alignment influences chemical enhancement mechanisms. A molecular library spanning a broad range of electronically tuned acceptor states was examined using surface-enhanced Raman scattering (SERS), vibrational sum frequency generation (vSFG) spectroscopy, and density functional theory calculations. By combining different excitation wavelengths with controlled variation of substrate composition and molecular electronic structure, the energetic relationship between plasmon-enhanced charge transfer excitation and molecular orbital alignment was quantitatively evaluated. The results reveal that charge-transfer enhancement of Raman scattering is governed by a well-defined interfacial resonance condition dependent on substrate work function and excitation energy. We further probe a subset of molecular-metal systems by nanocavity-enhanced vSFG and identify the same resonance conditions as in SERS, consistent with expectations. These findings establish an experimentally accessible framework for probing and engineering charge-transfer processes in plasmonic molecular junctions and provide mechanistic insight relevant to molecular plasmonics, charge carrier photophysics, and nanoscale interfacial spectroscopy.

We present an ab initio calculation of the band-to-band internal-conversion rate of the $\hbar\omega_{\rm nuc} \approx 8.35$ eV isomeric transition in $^{229}$ThO$_2$. Because the nuclear transition energy exceeds the electronic band gap of ThO$_2$, the isomer can decay nonradiatively by resonantly promoting a valence electron into the conduction band. We formulate this process as a Brillouin-zone sum over vertical interband transitions weighted by local Th-centered hyperfine matrix elements, which are evaluated directly from all-electron full-potential linearized augmented-plane-wave Bloch spinors. A finite nuclear magnetization model is included to regularize the short-range hyperfine interaction and to account for the Bohr-Weisskopf effect. After applying scissor shifts to span the experimentally reported ThO$_2$ band gaps, we find calculated internal-conversion lifetimes in the range of $1-16~\mu{\rm s}$. The lifetime increases strongly as the band gap approaches $\omega_{\rm nuc}$ because the resonant interband phase space at the nuclear transition energy is reduced. For the larger reported ThO$_2$ gaps, the calculated lifetime is comparable to the measured conversion-electron Mössbauer lifetime [Nature 648, 300 (2025)]. Our analysis implies that choosing solid-state hosts with band-gap values slightly lower than $\omega_{\rm nuc}$ can optimize solid-state nuclear clock performance with internal-conversion electron readout.

We derive a low-energy $\mathbf{k}\cdot\mathbf{p}$ effective Hamiltonian for monolayer osmium carbide (OsC) and ruthenium carbide (RuC) in a planar hexagonal configuration. First-principles calculations indicate that both monolayers are dynamically stable and exhibit features of a two-dimensional quantum spin Hall (QSH) phase, characterized by a nontrivial $\mathbb{Z}_2$ topological invariant. Using symmetry analysis at the $\Gamma$ point, we construct a multiband $\mathbf{k}\cdot\mathbf{p}$ Hamiltonian including spin-orbit coupling and reduce it to a four-band low-energy model through Löwdin partitioning. The effective Hamiltonian has a block-diagonal form, with two blocks related by time-reversal symmetry, analogous to the Bernevig--Hughes--Zhang (BHZ) model. In contrast to the standard BHZ form, the symmetry-allowed off-diagonal coupling contains quadratic momentum-dependent terms, which modify the low-energy dispersion near the $\Gamma$ point. The fitted parameters reproduce the ab initio band structures in the low-energy region, yielding a compact model for analyzing the electronic and topological properties of monolayer OsC and RuC.

Large language model (LLM) agents are increasingly expected to play a central role in AI-driven scientific discovery. Equipped with broad knowledge, flexible reasoning, and tool use, they have the potential to autonomously explore and solve scientific problems by repeatedly proposing hypotheses, testing them, and revising their beliefs in the light of the evidence. In current agents, however, these hypotheses, tests, and belief updates are buried in unstructured logs, and no mechanism lets the agent or the human researcher audit that process. Here we propose the Hypothesis Evolution Protocol (HEP), an agent harness that provides hypothesis generation, evaluation, and evolution as explicit, auditable operations. On materials-science research tasks, a HEP-equipped agent operates the hypothesis--test--evidence--belief cycle that planning-style agents lack, generalizes across research questions, and exploits the protocol more fully as the base LLM becomes more capable. These results mark a step toward auditable AI scientists, whose scientific reasoning can be inspected, verified, and built upon.

We present a novel real-time framework for the decay of metastable states in quantum field theories using Wigner functions. The framework introduces a nonperturbative nucleation rate formula that captures both quantum tunneling and over-the-barrier nucleation, alongside steps to evaluate it perturbatively. We apply it to a simple thermal example with direct relevance to current analog experiments. Our derived one-loop nucleation rate fundamentally differs from the widely cited high-temperature result by Linde: The prefactor contains quantum effects and also asymptotes to a differing form at high temperatures, where the quantum effects become negligible. Rather, the result is a generalization of Affleck's rate formula to quantum field theories, asymptoting to the effective field theory approach, and Langer's rate, at high temperatures. The example also reveals that the high-temperature side of the ``quantum-to-classical'' transition of thermal vacuum decay is still inherently quantum mechanical, even though the bounce background possesses the classical, $\mathrm{O}(d)\times S^1$, symmetry.

We employ large-scale density matrix renormalization group calculations to study the quantum phases of dipolar molecules confined in bent (zigzag) endofullerene chains, as a function of the chain angle $\gamma$. For LiF, ferroelectric order persists across the full range $60^\circ < \gamma 180^\circ$, with the critical effective dipole moment increasing as the chain bends and parallel alignment becomes less favorable. Near the equilateral configuration ($\gamma = 60^\circ$), geometric frustration drives a transition to an antiferroelectric Néel-ordered phase in which neighboring dipoles anti-align along the chain axis. We show that capturing this reorientation requires including dipolar couplings beyond the nearest-neighbor approximation, since next-nearest-neighbor interactions become equally strong at $\gamma = 60^\circ$. For confined water, o-D$_2$O reproduces both ordered phases, whereas p-H$_2$O -- owing to its large rotational constants -- develops no order at any chain angle despite the enhanced coordination of the bent geometry. Because a zigzag chain is the narrowest stripe of a two-dimensional lattice, these results suggest that engineered endofullerene layers could host a rich variety of dipole-ordered quantum phases beyond the ferroelectric ordering observed in previous work.

Cooperative multi-agent reinforcement learning is well suited to problems with large parameter spaces and exploitable local structure, such as the tuning of electrostatically-defined quantum-dot arrays. However, if parameter cross-talk is strong, a non-stationary environment from the perspective of any individual agent can destabilize learning - the same effect that plagues manual tuning of such systems. We propose using a factored representation of the action space, learned online, to decouple agents and minimize their interference. Our framework, QADAPT, uses this factorization to efficiently learn shared policies based on local measurements and rewards. With this modular strategy, we achieve zero-shot generalization to unseen quantum device sizes and maintain an approximately constant number of convergence steps to reach target regimes. This work provides a scalable route toward the rapid calibration of large-scale quantum processors.

We investigate the forces acting on a proton vortex in superconducting neutron-star matter composed of neutrons, protons, and electrons, accounting for Fermi-liquid interactions in the neutron-proton subsystem. While the force arising from electron scattering by the vortex magnetic field is well known, we demonstrate that normal neutrons also exert a force on proton vortices due to Fermi-liquid coupling with superconducting protons. Using a kinetic approach based on Landau Fermi-liquid theory, we show that neutron quasiparticles scatter off a proton vortex through the spatially varying condensate momentum, giving rise to a longitudinal force proportional to the relative neutron-vortex velocity. In contrast to the electron contribution, this force has no transverse component and vanishes in the absence of neutron-proton Fermi-liquid interaction. Within a simple model, we analyze the main properties of this force and provide physically motivated estimates of its magnitude.

Hybrid devices integrating quantum dots with Josephson junctions are gaining interest because they combine spin-based quantum computing with circuit quantum electrodynamics (circuit QED) methods. In particular, Andreev spin qubits have shown significant experimental progress including strong two-qubit coupling, and are predicted to exhibit all-to-all connectivity. Here we propose superconducting singlet-triplet (SST) qubits that rely on parallel-aligned double quantum dots in Josephson junctions. While Andreev spin qubits require spin-orbit interaction to unlock the spin degree-of-freedom, SST qubits do not require spin-orbit interaction, making the advantages of hybrid devices available to a wider range of materials. Similar to Andreev spin qubits, the qubit states couple to the superconducting phase across the junction, which allows for control and readout using circuit QED, and supports all-to-all connectivity. Only $N$ flux lines are required to perform any single- and two-qubit gate among $N$ qubits, and thus the overhead of control lines is small. Finally, linear protection from charge or flux noise makes these qubits interesting candidates for a future quantum processor.

We consider quasi-one-dimensional (Q1D) continuous waves (CWs) in the two-dimensional (2D) optical system with the cubic-quintic nonlinearity and a Q1D potential trough. In the case of a smooth trough profile, we confirm the known modulational instability (MI) of Q1D CWs with the transverse structure corresponding to the 1D ground state (GS) in the potential trough, and demonstrate the MI of CWs with the dipole-mode (DM) transverse structure, corresponding to the lowest 1D excited state in the potential trough. The CWs of both GS and DM types remain nearly stable close to the edges of their existence regions. Stable stationary states in the form of periodic chains of 2D solitons, trapped in the potential trough, are produced in a numerical form. The dynamics of the soliton chains excited by a localized kick is studied too. For the potential trough with the singular delta-functional profile, we find two species of exact analytical solutions for CWs, one of which is completely stable.

Financial and economic history is strewn with bubbles and crashes, booms and busts, crises and upheavals of all sorts. Understanding the origin of these events is arguably one of the most important problems in economic theory: are economies intrinsically unstable, and can one ``stabilize unstable economies''? In this review I argue, from a physicist's vantage point, that the concept of equilibrium -- so central to mainstream economic thinking -- is likely to be the exception rather than the rule in large, complex, interacting systems. Drawing on a series of stylized ``toy'' models, I show how excess volatility, endogenous crises and crashes, inflation swells and persistent inequalities can all emerge naturally from genuinely out-of-equilibrium dynamics, without invoking large exogenous shocks. Three generic mechanisms recur throughout: trapping in a multiplicity of history-dependent equilibria; the impossibility of dynamically reaching equilibrium, leading to oscillations and chaos; and the spontaneous evolution towards fragile, marginally stable states -- the self-organized criticality paradigm. I stress that these are phenomenological scenarios rather than calibrated theories: there is, at this stage, no ``smoking gun''. But the burden of proof, I contend, should be on the equilibrium camp.

Tunneling is essential in the initialization, measurement, and control of quantum dot qubits. In silicon, such tunneling connects not only the qubit states but also valley minima in the conduction band on opposite sides of the Brillouin zone, with large consequences for the quantum dot behavior. Here we present a full characterization of the intravalley and intervalley tunnel couplings, including their complex phases -- the valley phases. These phases are shown to control measurable parameters, including the ratios of the gaps at anticrossings between quantum states of a double quantum dot. The valley phases themselves evolve as a function of the quantum dot gate voltages and depend on the underlying atomic structure of the quantum well. Knowledge of the valley phases completes the picture and fills a key gap in our understanding of sample-wide variations of valley couplings and the physical parameters that depend on them, including spin-orbit coupling, valley-orbit mixing, and Landé $g$-factors.

Floquet engineering, which involves controlling systems through time-periodic driving, is a method for coherently manipulating quantum materials and realizing dynamical states with novel functionalities. Most research in solid-state systems has focused on the use of AC-\textit{electric} fields as the controlling drive. In this study, we investigate the effects of AC-\textit{magnetic} fields on two-dimensional (2D) Dirac electrons and report the emergence of new states and new transport phenomena. In a magnetic field that temporarily changes its direction, the 2D Dirac electrons form a new localized state with a flat band dispersion, dubbed as a $\pi$-Landau level. Its wave function is a superposition of the clockwise and counterclockwise cyclotron orbits with time-periodic amplitudes, resulting in a novel closed trajectory shaped like a figure eight. Then, what would be the counterpart of the Hall effect in AC-magnetic fields? We find that a DC-current in the transverse direction, \textit{i.e.} a homodyne Hall current, is generated when an additional AC-electric field is applied. In the case of Dirac electrons, several electronic states contribute to this phenomenon including the $\pi$-Landau level. However, when the chemical potential $\mu$ is near the Dirac point, the dominant contribution comes from the low-energy electrons and we numerically find the homodyne Hall current to behave as $I_y=-\frac{e}{h}\mu$ per valley and spin. We explain this phenomenon through the high-frequency effective Floquet Hamiltonian which resembles the chiral Landau level Hamiltonian of three-dimensional Weyl Hamiltonian exhibiting chiral anomaly. We discuss the experimental feasibility and conclude that it is possible to realize this new exotic state using techniques such as THz metamaterial enhancement of magnetic fields.

A point-like defect in a uniform current-carrying conductor induces a dipole in the electrochemical potential, which counteracts the original transport field. If the mean free path of the carriers is much smaller than the size of the defect, the dipole results from the purely diffusive motion of the carriers around the defect. In the opposite limit, ballistic carriers scatter from the defect -- for this situation, Rolf Landauer postulated the emergence of residual resistivity dipoles that are independent of the defect size and thus impose a fundamental limit on the resistance of the parent conductor. Here, we study resistivity dipoles around holes of different sizes in two-dimensional Bi films on Si(111). Using scanning tunneling potentiometry to image the dipoles, we find a crossover from linear to constant scaling behavior of their amplitudes with defect size, manifesting the transition from diffusive to Landauer dipoles. The extracted parameters of the transition allow us to estimate the Fermi wave vector and the carrier mean free path in our Bi films.

Aluminum-based Josephson junctions are currently the main sources of nonlinearity for control and manipulation of superconducting qubits. A constriction-based junction provides an alternative source of nonlinearity that promises new types of protected qubits and the possibility of high-temperature and high-frequency operation through the use of superconductors with larger energy gaps. Junctions made from such superconductors have been challenging to incorporate into superconducting qubits because of difficulty controlling junction parameters and have had extremely low lifetimes, which limited their utility. Here we demonstrate that junctions made using titanium nitride (TiN) are a promising and controllable qubit platform. We use TiN junctions to build superconducting qubits based on quantum phase slips through the junction. We operate the qubit at zero flux where the qubit frequency (~17 GHz) is mainly determined by the inductance of the qubit. We perform readout and coherent control of the superconducting qubit, and measure qubit lifetimes >60 $\mu$s. Finally, we demonstrate operation of the qubit at temperatures exceeding 300 mK. Our results add the TiN-based junction as a tool for superconducting quantum information processing and opens avenues for new classes of superconducting qubits.

Theoretical studies of X-ray photoemission from simple metals have traditionally focused on the frequency domain, aiming to reproduce experimental spectra. Here, we investigate the same problem in the time domain in search of physical insight and methodological advances. Our results reveal prominent aspects of the problem that are inconspicuous in the frequency domain. The calculated $\mathcal{F}(t)$ exhibits a weakly damped harmonic oscillation that modulates the Doniach-Sunjic power-law decay of the photoemission rate, a behavior arising from the coherent interference between two classes of particle-hole excitations. From a methodological perspective, we advance the time-dependent numerical renormalization-group (NRG) approach by exploring the eNRG method, a real-space variant that is more flexible than Wilson's construction. Giving special attention to strong core-hole potentials, we compare the photocurrents obtained from two complementary time-dependent eNRG algorithms with (i) an analytical expression for $\mathcal{F}(t)$ that becomes highly accurate at moderately long times and (ii) results from numerical diagonalization of a tight-binding Hamiltonian, which covers the time interval in which our analytical expression is less precise. Anticipating extensions to correlated-impurity models, we identify the sources of deviation and discuss the virtues and drawbacks of the two algorithms.

The sequential response of frustrated materials-ranging from crumpled sheets and amorphous media to metamaterials-reveals their memory effects and emergent computational potential. Despite their spatial extension, most studies rely on a single global stimulus, such as compression, effectively reducing the problem to scalar driving. Here, we introduce vectorial driving of frustrated materials by applying multiple spatially localized stimuli to explore path-dependent, sequential responses. We uncover a wealth of phenomena absent in scalar driving, including non-Abelian responses, mixed-mode behavior, and chiral loop transients. We show that fold singularities connect three states -- ancestor, descendant, and sibling. This recurring pattern serves as the elementary building block of all sequential paths. We then introduce three levels of description of sequential, path-dependent responses. At the most fundamental level, path-dependent transition graphs (pt-graphs) and strain maps capture the response under arbitrary vectorial driving and connect pathways to the underlying singularities. They provide a complete description analogous to transition graphs (t-graphs) for scalar driving. However, as pt-graphs and strain maps become unwieldy for high-dimensional driving, we introduce b-graphs -- graphs whose nodes and transitions encode the systems response to binarized vectorial driving. These present a less complete but much simpler second-level description by restricting attention to binary input and their induced transitions. Finally, we introduce graph-based motifs that enable a systematic analysis of b-graphs. As statistical measures of pathway complexity, these motifs can be obtained in systems of any size or complexity. Our work paves the way for strategies to explore, harness, and understand complex materials and memory, while advancing embodied intelligence and in-materia computing.

Two-dimensional PtSe2 exhibits outstanding intrinsic properties such as high carrier mobility, tunable bandgap, broadband absorption and air stability, making it ideal for (opto)electronic applications. In particular, thick PtSe2 is semimetallic and well suited for ultrafast optoelectronics in the infrared domain. However, achieving PtSe2 films of high crystalline quality with controlled properties on low-cost and insulating substrates remains challenging. Here, highly crystalline semimetallic PtSe2 films are grown by molecular beam epitaxy on sapphire substrates. It is shown how an optimized post-growth annealing remarkably improves the out-of-plane crystallinity and leads to record sheet conductances, up to 1.6 mS. In-depth structural analyses reveal the strong influence of the domain arrangement within the films on their electrical properties. Films that are mainly composed of vertically single crystalline domains exhibit high sheet conductance (1.1 - 1.6 mS), whereas films that contain superimposed twisted domains present low sheet conductance (0.5 - 0.6 mS). Moreover, it is demonstrated that the A1g Raman peak width, in addition to the commonly used Eg peak width, are both effective metrics for evaluating the quality of PtSe2: films with narrower Eg and A1g peaks exhibit higher in-plane and out-of-plane crystalline quality, respectively, as well as higher sheet conductance. Finally, coplanar waveguides integrating a semimetallic PtSe2 channel are fabricated on a 2-inch sapphire substrate to demonstrate optoelectronic devices operating at the 1.55 {\mu}m telecom wavelength. This includes photodetectors with a record 60 GHz bandwidth and the first PtSe2-based optoelectronic mixer with a bandwidth above 30 GHz.

Two-dimensional coherent spectroscopy (2DCS) provides insights into the nonlinear response of correlated lattice systems. We simulate multipulse excitations in the Hubbard model using nonequilibrium dynamical mean-field theory to extract the 2DCS signal of weak-coupling antiferromagnets with and without local potential disorder. By comparing calculations with static and dynamic Hartree terms, and analyzing the waiting-time dependence of the signal, we identify the contribution of the collective amplitude mode to the spectroscopic features and the relevant underlying processes. With broadband pulses, the rephasing and nonrephasing peaks at the gap energy are found to be of predominant amplitude mode character. Using narrow-band pulses, we also demonstrate a strong enhancement of these amplitude mode-related signals at a pulse frequency of half the gap size.

Locally noncentrosymmetric materials are attracting significant attention due to the unique phenomena associated with sublattice degrees of freedom. The recently discovered heavy-fermion superconductor CeRh$_2$As$_2$ has emerged as a compelling example of this class, garnering widespread interest for its remarkable temperature-magnetic-field phase diagram, which features a field-induced first-order superconductor-to-superconductor phase transition with nontrivial dependence on the field direction and high critical fields, as well as antiferromagnetic and potentially higher multipole orders. To investigate the complex interplay of the ordered phases in CeRh$_2$As$_2$, we develop a theoretical framework based on symmetry analysis applied to a Bogoliubov--de Gennes Hamiltonian and Landau methods. This approach allows us to propose probable symmetries of the superconducting states and elucidate their close relationship with magnetism. Among other results, we find that the near degeneracy of two pairing symmetries is naturally explained if and only if intralayer spin-orbit coupling is large compared to interlayer hopping. Intriguingly, we find that the first-order transition can be interpreted as a transition between coexistence phases of the same superconducting order parameters, albeit with distinct admixtures. This line may end in a critical end point below the superconducting critical temperature. Our approach accurately reproduces current experimental phase diagrams for varying temperature as well as out-of-plane and in-plane magnetic field, both if the transition to a magnetic phase occurs below the superconducting critical temperature and if it occurs above. Furthermore, we calculate the magnetic susceptibility and the specific heat and compare these quantities to recent experimental results.

In these notes, we present a rigorous and self-contained introduction to the fundamental concepts and methods of quantum many-body theory. The text is designed to provide a solid theoretical foundation for the study of interacting quantum systems, combining clarity with mathematical precision. Core topics are developed systematically, with detailed derivations and comprehensive proofs that aim to make the material accessible to graduate students and beginning PhD students. Special attention is given to formal consistency and pedagogical structure, so as to guide the reader through both the conceptual and technical aspects of the subject. This work is intended as a reliable starting point for further exploration and research in modern quantum many-body physics.

We perform the variational Monte Carlo calculation for recently proposed chiral superconducting states driven by strong Coulomb interactions. We compare the resulting energetics of these electronic phases for the electron dispersion relation $E_k = c_2 k^2+c_4 k^4$. Motivated by the recent discovery of chiral superconductivity in rhombohedral graphene systems, we apply our analysis to relevant parameter regimes. We demonstrate that topological chiral superconducting phases (including a spin-unpolarized state) can be energetically favored over the spin-valley polarized Fermi liquid above the density of Wigner crystal phase. Our results show that the preference for chiral superconductivity is strongest when $c_2$ lies between zero and a negative value corresponding to a Fermi sea on the verge of forming a hole pocket around $k=0$. This finding suggests that superconductivity can arise from pure repulsive Coulomb interactions in systems with an almost flat band bottom, without relying on the pairing instability of a Fermi surface. This mechanism opens a new pathway to superconductivity beyond the conventional BCS mechanism.

Generative inverse design is reshaping the discovery of functional crystalline materials. Yet current generative models face challenges in simultaneously achieving stability, novelty, and precise controllability in a single trained model. We address these challenges with a key physical insight: the diversity of crystals is governed by their crystallographic information (CI), namely composition, space group and lattice, whereas only a few stable atomic configurations remain once the CI is fixed. Built on this insight, we introduce SCGEN (stable and controllable crystal structure generation), a physics-grounded generative model with two components. A variational autoencoder samples diverse, physically plausible CI, and a symmetry- and Wyckoff-position-constrained optimizer locates stable atomic positions via universal machine-learning potentials. Benchmarked on roughly two million structures, SCGEN reaches state-of-the-art stability while preserving comparable novelty, and it satisfies any specified composition, space group, lattice or joint constraint with 100% success and no task-specific retraining. Applied to photocatalytic water splitting, property-guided optimization with SCGEN generates 200,000 candidate structures and identify top 22 stable, active and synthesizable photocatalysts. By decoupling CI generation from coordinate optimization, SCGEN establishes a physics-grounded inverse-design paradigm that yields synthesis-ready crystals on demand, rather than structures requiring post hoc repair, relaxation, or retraining.

Recently, several spin chain models have been discovered that admit solutions in terms of "free fermions in disguise." A graph-theoretical treatment of such models was also established, giving sufficient conditions for free fermionic solvability. These conditions involve a particular property of the so-called frustration graph of the Hamiltonian, namely that it must be claw-free. Additionally, one set of sufficient conditions also requires the absence of so-called even holes. In this paper, we present a model with disguised free fermions where the frustration graph contains both claws and even holes. Special relations between coupling constants ensure that the free fermionic property still holds. Notably, the central elements associated with the even holes can be removed by fixing the gauge, revealing our model to be an integrable deformation within the original algebra of free fermions in disguise. The transfer matrix of this model can be factorized in a special case, thereby proving the conjectured free fermionic nature of a special quantum circuit published recently by two of the present authors.

In a two-dimensional kagome monolayer, a nontrivial intrinsic Berry curvature may arise in the $d$-electron system from the interaction with a non-collinear spin order induced by an underlying antiferromagnetic exchange. This opens the route for a quantum anomalous Hall effect in the multi-orbital system, even without an external magnetic field, explicit spin-orbit coupling or relativistic effects. For spin orders with an out-of-plane component, the scalar spin chirality is finite, and the integration of the Berry curvature over the Brillouin zone may yield integer Hall conductivities in units of $e^2/h$. For a Fermi level within a nontrivial gap, the canted configuration offers, at least in principle, the possibility of a maximal Chern number, $C=\pm 5$. Candidate materials are considered in this paper. In existing materials, the electron hopping is generally highly anisotropic, leading to a quantum anomalous Hall effect with smaller Chern numbers. A topological phase transition between Hall plateaus of opposite $C$ can be driven by flipping the out-of-plane component of the spin order, alluding to the potential of this system to applications in quantum information.

We study the activity of 'living' droplets, which confine tens of swimming meso-organisms in 3D using a superhydrophobic substrate. With few swimmers, the droplet oscillates at its inherent resonant frequency. We observe deviations from this classical regime as the level of confinement or crowding increases and develop scaling law models to successfully describe our results. We report a difference in swimming kinematics in crowded 3D environments compared to quasi-2D. Our work reveals mechanisms for bio-inspired droplet actuation with implications for mesoscale robotics, fluidics, and sensing.

Valley, as a ubiquitous degree of freedom in lattices, has found wide applications in both electronic and classical-wave devices in recent years. However, achieving valley-polarized states, a prerequisite for valley-based operations, still remains challenging. Here, we propose and experimentally demonstrate a universal non-Hermitian mechanism for valley filtering using only uniform background dissipation, which creates a propagation length contrast between valleys through their intrinsic group velocity differences. We implement this concept in an acoustic crystal, observing switchable and robust valley polarization of sound through large-scale field mapping. Remarkably, our approach is solely based on uniform loss, without the need for any special lattice structures, tailored excitations, or external fields. We further provide designs of our non-Hermitian valley filter on photonic and electronic platforms. Our results offer a simple and effective solution to valley-polarized state generation and may advance the development of novel valley-based devices in both classical and quantum regimes.

We investigate the effects of Van Hove singularities in the electronic band structure of two-dimensional semiconductors on the energy dispersion of excitons and positive trions. In particular, we study valence band energy dispersions which possess (i) a typical logarithmic Van Hove singularity, (ii) a line high-order Van Hove singularity (HOVHS) from a Mexican-hat dispersion or (iii) a point HOVHS such as a monkey saddle. We find that the density of states (DOS) of excitons and trions containing such singularities is dramatically enhanced and shows, in general, how the HOVHS in the valence band can strongly affect and be mirrored in the DOS of excitons and trions. This leads to new states that govern the optical properties of the system. In addition, we study a set of materials, InSe, GaSe and $\alpha$-SnAs, from a class of materials in which the topmost valence band has an inverted Mexican-hat shape. The most favourable exciton occurs when the singularity is at the $\Gamma$-point, as in the example of monolayer $\alpha$-SnAs, which hosts a HOVHS. Our work thus provides a pathway to engineer specific bound states in two-dimensional materials that host such singularities, thereby opening new avenues for potential applications.

Understanding the fundamental properties that dictate photoexcited polarons in materials is critical to tuning their properties. Theoretical models of polarons have only recently been extended to the excited state. Experimental measurements of polaron formation and transport have been widely undertaken across a range of materials, from photocatalysts and superconductors to soft conducting polymers. Here, we map thermalized excited state experimental measurements of quantities such as polaron strength onto phase diagrams of the Holstein, Hubbard-Holstein, and t-J-Holstein models. This work demonstrates that tuning electron-phonon coupling strength, electron localization, and spin exchange can be leveraged to suppress or control polarons in transition metal oxides. We find that the t-J-Holstein model best describes the measured iron oxides and could be generally applied to a wide range of systems that exhibit polaron formation in the excited state. This work combines experimental data with ground state models to provide a qualitative parameter space for informing photoexcited polaron design, under which excited state polaronic behavior can be classified within ground-state calculable models.

We develop theoretical methods to predict the effects of paramagnetic impurities on the microwave response of conventional spin-singlet superconductors. Our focus is on superconducting devices and resonators with low concentrations of impurities and exchange interactions with conduction electrons. We connect the sub-gap quasiparticle spectrum generated by pair-breaking to the frequency and temperature dependence of the conductivity for superconductors operating at microwave frequencies. We report theoretical results for superconducting device performance -- dissipation, quality factor and frequency shift anomalies -- based on self-consistent calculations of the current response and penetraion of the electromagnetic field at the vacuum-superconducting interface. Key results include the prediction of a non-monotonic anomaly in the low-frequency superfluid fraction and penetration depth at very low temperatures related to the sub-gap quasiparticle spectrum. Dissipation of microwave power is predicted from intra- and inter- impurity band transitions at GHz frequencies at low temperatures, including a physical mechanism responsible for residual resistance. We predict anomalies in the resonant frequency, $f(T)$, and quality factor, $Q(T)$, of high-Q SRF cavities operating in the GHz range at low-temperatures that are sensitive to non-magnetic and paramagnetic impurity disorder.

We study polaritonic bound states in the continuum (BIC) created in GaN waveguides. The existence of symmetry-protected BICs is confirmed by the suppression of light emission and the observation of a polarization vortex in momentum space. Upon increasing the pumping, polariton population accumulates at the BIC and we observe polariton lasing from the blueshifted BIC states. The assessment of the polariton BIC emission energy and of its real and momentum space wavefunctions as a function of pumping power, i.e. of polariton density, indicates the formation of a bright soliton above the lasing threshold. Soliton formation at the BIC is induced by the combination of negative mass BIC and of repulsive polariton-polariton and polariton-reservoir interactions.

We performed de Haas-van Alphen (dHvA) oscillation measurements and band-structure calculations for YNi$_2$B$_2$C. Our improved band structure successfully explained the origins of the large dHvA frequencies $\beta$ and $\zeta$, which were inexplicable in previous works. By comparing experimental effective masses with band masses, we determined the electron-phonon coupling for each orbit. The results showed a clear Fermi-surface-sheet dependence of the electron-phonon coupling strength, especially highlighting that the coupling for the band-28 sheet is very weak, almost absent for the orbit with $B \parallel c$. This finding is consistent with previous observations of dHvA oscillations from this orbit in the mixed state down to very low fields. Amidst growing interest in high-temperature superconductivity driven by electron-phonon coupling in hydrides under high pressure, this study provides foundational data pivotal to precisely understanding electron-phonon coupling.

Significant efforts have been devoted to studying the properties of quantum droplets, an ultra low-temperature phase of bosonic quantum matter that emerges as a consequence of the Lee-Huang-Yang fluctuating correction. However, the temporal dynamics of droplet formation for heteronuclear bosonic mixtures is only partially understood. Here, we numerically analyze the droplet formation process for homonuclear and heteronuclear binary bosonic mixtures in one dimension, using a tight-binding model and real-time evolution with a novel, highly robust integration algorithm. We proceed with a systematic scan of interaction intensities, mass ratios, and initial conditions that allows us to characterize quantitative criteria for droplet formation and equilibrium prospects. Noticeably, most droplets readily form across the entire parameter space, although only a small fraction achieves a stable equilibrium configuration within the simulation horizon. We attribute this equilibrium deficiency to damping from a breathing mode, which we extract directly from the width oscillations at late times. Within the balanced subset of results, the Lee-Huang-Yang contribution dominates the binding energy, while the density profile at late times is better described by a soliton-like shape rather than the flat-topped profiles characteristic of larger droplets. Heteronuclear droplets achieve deeper bindings and exhibit larger breathing oscillations in comparison to their homonuclear counterparts, which is consistent with the role of mass-imbalanced kinetic terms. Our analysis distinguishes physical results from artifacts of incomplete relaxation and provides a methodological framework for interpreting real-time quantum droplet simulations in the constrained dissipation regime characteristic of one-dimensional geometries.

High-throughput powder X-ray diffraction (XRD) simulations are a key prerequisite for generating large datasets used in the development of machine-learning models for XRD-based materials analysis. However, the widely used pymatgen powder XRD calculator, implemented entirely in Python, can be computationally inefficient for large-scale workloads, limiting throughput. We present XRD-Rust, a Rust-accelerated implementation of the pymatgen powder XRD calculator that maintains compatibility with existing Python-based workflows. The method retains pymatgen for crystal structure handling while reimplementing the computationally intensive parts of the XRD calculation in Rust, with optional further acceleration via SIMD vectorization and multi-threaded execution across reflections using the Rayon library. Performance benchmarking on two large crystallographic datasets, the Materials Cloud Three-Dimensional Structure Database (MC3D, 33 142 structures) and the Crystallography Open Database (COD, 515 181 structures), demonstrates substantial speedups. For MC3D, XRD-Rust achieves a median serial SIMD speedup of 10.8x (median absolute deviation, MAD, 1.8x), increasing to 15.1x (MAD 3.8x) with 8-thread parallel execution, with a maximum runtime reduction from 40.5 s to 0.9 s. For COD, the median serial SIMD acceleration reaches 10.7x (MAD 4.2x), while 8-thread parallel execution yields 19.5x (MAD 10.0x), with a maximum reduction from 1437 min to 1 min. SIMD vectorization alone provides additional performance gains ranging from a few percent to several tens of percent, depending on the workload and execution mode. Overall, these benchmarks demonstrate that XRD-Rust significantly accelerates powder XRD simulations compared to the original pymatgen implementation, enabling efficient high-throughput dataset generation and improving performance in interactive diffraction analysis applications.

We investigate the dynamics of a dense raft of millimeter-sized granular particles at a vertically vibrated air-water interface, which displays a rich set of patterns and particle dynamics as we vary the vibration amplitude, frequency, and particle packing fraction. While the classical parametric instability with standing waves still occurs over a certain range of parameters, the measured wave dispersion relations indicate an increasing role of the raft's emergent elasticity at higher packing fractions, where the effective surface tension decreases and the out-of-plane bending modulus increases. At higher vibration frequencies and lower amplitudes, we identify a regime without standing waves. Instead, individual particles exhibit thermal-like motion, with transport crossing over from diffusive to sub-diffusive as the packing fraction increases. The particle dynamics also display spatial and temporal heterogeneity, as in supercooled liquids. Starting from this regime, when the vibration amplitude is further increased, a large cavity eventually forms inside the raft, whose size and shape depend on the vibration frequency and the injected vibration energy. The cavitation results in the coexistence of free-surface water waves inside the cavity and thermal-like particle motion in the surrounding raft.

Patterned two-dimensional (2D) magnetic nanostructures constitute geometry-engineered spin systems in which exchange, anisotropy, dipolar coupling, and finite-size effects operate on comparable energy scales. Spatial modulation of continuous magnetic films produces confinement-driven critical behavior, compensation phenomena, metastable switching pathways, and topologically nontrivial textures such as vortices and skyrmions. Computational modeling plays a central role in resolving this complexity, enabling quantitative construction of thermodynamic phase diagrams and analysis of geometry-dependent stability regimes. This review synthesizes theoretical and numerical frameworks for patterned 2D magnetism, including classical spin models, stochastic spin dynamics, rare-event methods, and multiscale parameterization informed by first-principles calculations. Representative systems-nanodot and antidot arrays, artificial spin-ice lattices, exchange-modulated heterostructures, and patterned van der Waals magnets - illustrate how geometry functions as an effective thermodynamic control parameter. Emerging directions in nonequilibrium modeling, multiphysics coupling, and scalable data-centric workflows are discussed in the context of predictive phase mapping. Patterned 2D magnetism thus exemplifies geometry as an effective thermodynamic variable in materials engineering and computational statistical physics, with phase stability and controlled spin textures at the core of next-generation spintronic architectures.

Moiré network formation in graphene bilayers breaks stacking symmetry and generates topological domain walls (TDWs) that host one-dimensional boundary states. Here we show that these TDWs undergo an additional symmetry breaking at the level of network geometry, leading to the spontaneous emergence of chiral configurations through lattice relaxation. Using atomistic structural relaxation, we establish a strain-flexibility phase diagram with three equilibrium morphologies -- straight, mono-chiral, and dual-chiral -- arising from TDW energy minimization within the geometric constraints of the moiré network. Tight-binding calculations show that straight and chiral networks support distinct low-energy electronic regimes. Straight networks retain a pronounced topological crossing point (TCP)-centered local peak pair, nearly balanced edge states on the two sides of a TDW, and uniformly gapped AB/BA domains with smooth low-energy interference patterns. Chiral networks suppress the TCP-centered peak pair, redistribute low-energy spectral weight toward the walls, drive curvature-selective edge asymmetry, and break the smooth domain interference pattern into step-like features. Chirality, therefore, selects between a TCP-centered regime with concentrated low-energy response and a wall-centered regime with directional boundary localization, opening distinct routes toward low-energy states concentrated at TCP regions or directional one-dimensional channels along the walls.

Thermomechanical Controlled Processing (TMCP) is widely used to control the microstructure and properties of linepipe or high strength low alloy steels (HSLA). These steels are often joined by welding and used in demanding environments such as the Arctic. In these materials, the thermal path the steel experiences is critical for understanding microstructural evolution during processing. A key step is the solid-state phase transformation during cooling from the high-temperature austenite to the room-temperature microstructure which significantly influences the final mechanical properties. We used 3D electron backscatter diffraction (EBSD) to explore the relationship between the austenite phase and the room temperature microstructure. A significant result for the present work is the collection, and analysis, of data from a large volume (150 x 150 x 100 um3, with a (200 nm)3 voxel size) which enables analysis of a complete prior austenite grain. This grain is twinned, allowing us to analyse the variants at the twin boundary in this grain, which offers new insights into the mechanisms of the transformation to the low temperature phase by highlighting the significant of the twin boundaries on the variants present. This suggests opportunities to engineer novel microstructures by controlling the high-temperature grain boundary character.

We present a computational approach for exciton calculations in two-dimensional (2D) materials within the Bethe-Salpeter equation (BSE) framework, employing an atomistic description with point-like orbitals. Unlike widespread efficient calculations that rely on classical or effective interaction models, such as the Rytova-Keldysh model, our method incorporates quantum screened interactions. By explicitly computing the 2D dielectric function at the random-phase approximation level, we capture screening effects beyond such approximations with an accuracy akin to first-principles methods. Consequently, we can realistically estimate excitonic binding energies with a bearable computational cost. A detailed account of the various convergence parameters sheds light on a possible cause of the large dispersion of binding energies reported in the literature using first-principles GW/BSE implementations. This work thus provides an alternative pathway towards efficient and faithful dielectric screening and exciton computations in low-dimensional materials.

We investigate optimal band structures in band-converged systems to achieve high zT using numerical calculations based on a virtual spectral conductivity model. We consider a two parabolic band system, in which multiple band parameters can be independently controlled. Despite its simplicity, this model provides theoretical validation of empirical trends observed in thermoelectric materials. Our results provide a physically transparent set of design principles for band-structure engineering, offering quantitative design guidelines for the development of a wide range of thermoelectric materials. The main conclusions are as follows: (i) When a band does not cross the chemical potential and |{\mu}-E_edge |>5k_B T, the contribution of the band to zT is negligibly small; (ii) To suppress the bipolar effect, a band gap E_g satisfying E_g>5k_B T_op, where T_op is the operating temperature, is required; (iii) In band-converged systems, the energy separation between the band edge {\Delta}E should satisfy {\Delta}E~0 to maximize zT when interband scattering is insignificant; (iv) Achieving high spectral conductivity {\Sigma} (high band degeneracy N, density of states effective mass m_DOS^*, and relaxation time {\tau}) near the band edge is essential for achieving high zT.

We study the nonlinear chaotic dynamics in a system of linear oscillators coupled by social network links with an additional stratification of oscillator energies, or frequencies, and supplementary nonlinear interactions. It is argued that this system can be viewed as a model of social stratification in a society with nonlinear interacting agents with energies playing a role of wealth states of society. The Hamiltonian evolution is characterized by two integrals of motion being energy and probability norm. Above a certain chaos border the chaotic dynamics leads to dynamical thermalization with the Rayleigh-Jeans (RJ) distribution over states with given energy or wealth. At low energies, this distribution has RJ condensation of norm at low energy modes. We point out a similarity of this condensation with the wealth inequality in the world countries where about a half of population owns only a couple of percent of the total wealth. In the presence of energy pumping and absorption, the system reveals features of the Kolmogorov-Zakharov turbulence of nonlinear waves.

Motivated by growing evidence for the significance of the third-nearest-neighbor Heisenberg ($J_3$) interaction in candidate Kitaev materials, we investigate the dynamical properties of the Kitaev spin liquid (KSL) under a $J_3$ perturbation, focusing on its spin dynamical structure factor (DSF) and Raman scattering. Within a self-consistent parton mean-field plus random-phase approximation framework, we find that $J_3$ induces coherent, paramagnon-like collective modes that coexist with a high-energy Majorana continuum in the spin DSF. The softening of these modes with increasing $|J_3|$ signals a quantum phase transition to magnetic order. Remarkably, magnetic ordering sets in at a common critical $J_3$ for both ferromagnetic ($K<0$) and antiferromagnetic ($K>0$) Kitaev models, with the resulting ordered states forming exact dual pairs under a four-sublattice duality transformation that maps $(K,J_3) \rightarrow (-K,J_3)$. An external magnetic field further softens the preexisting paramagnon modes, thereby enhancing magnetic order. Perturbative Raman calculations show that while the Kitaev-like Raman vertex probes only itinerant matter Majorana fermions, the response from the $J_3$-like vertex features both matter Majoranas and visons. Four-vison excitations produce a sharp peak accompanied by a two-fermion continuum, whereas two-vison excitations yield a continuum closely resembling the single-matter-fermion density of states. These results provide a unified perspective on the dynamical signatures of $J_3$-perturbed KSL and are helpful for interpreting experimental spectra in candidate Kitaev materials with sizable $J_3$ interactions.

Non-linear oscillators serve important functions in many biological systems, including within the inner ear and neuronal networks. The sustainment of oscillations in noisy environments requires continuous energy dissipation, quantified by the steady-state entropy production rate (EPR). We study an idealized, analytically tractable model of a stochastic circular limit cycle and examine how mutual coupling in pairs and populations alters dissipation. For a single oscillator, the EPR depends on three key factors: intrinsic frequency, tangential velocity fluctuations, and mean tangential velocity. The dynamics are characterized by a dimensionless effective temperature given by the ratio of intrinsic relaxation and diffusion timescales. For radial (amplitude), phase (Kuramoto-like), and Cartesian couplings, we derive analytical expressions for the EPR and confirm them numerically. Varying the effective temperature and system size strongly influences how the EPR depends on coupling strength and, in some cases, results in qualitatively distinct behaviors. Moreover, the coupling types affect the tangential velocity distributions differently. Notably, in all cases studied, Cartesian coupling reduces the EPR relative to the uncoupled system, irrespective of effective temperature and system size. The analysis of idealized non-linear oscillators reveals that different classes of coupling interactions and competing timescales present in the oscillators have distinct effects on energy dissipation.

Chiral active Brownian particles (CABPs) are self-propelled agents with intrinsic rotational dynamics, giving rise to circular trajectories commonly observed in biological and synthetic microswimmers. Understanding how CABPs explore confined environments and locate targets is crucial for characterizing transport, search efficiency, and reaction processes in physical and biological systems. We study the escape dynamics of CABPs from one- and two-dimensional confined domains. In one dimension, we consider intervals with either two absorbing boundaries or a reflecting boundary on one side and an absorbing boundary on the other, and derive closed-form asymptotic solutions in the high-chirality regime, revealing the quantitative scaling of the mean first passage time (MFPT) as a function of particle rotation speed (chirality). In two dimensions, we analyze escape from a disk containing one absorbing arc or two symmetric absorbing arcs. By numerically solving the governing partial differential equations, we compute the MFPT for CABPs to escape the domains as a function of the particle's initial orientation, self-propulsion speed, angular velocity, and domain geometry. Our results show that, depending on the parameters and geometry, the MFPT can exhibit non-monotonic behavior as a function of chirality. A minimal escape time exists at an intermediate value of chirality, where the rotational time scale balances the active swimming time scale, redirecting a particle towards the exit which would otherwise be blocked due to unfavorable initial orientation. Our work offers a comprehensive characterization of CABP escape dynamics in canonical confinements and identifies chirality as a key control parameter for transport and search in confined physical and biological systems.

Chirality, defined by the absence of mirror and inversion symmetries, has attracted considerable attention owing to its unique physical phenomena, including cross-correlated responses such as current-induced magnetization (CIM) and chiral phonons. Recently, it has been established that chirality is characterized by electric toroidal (ET) multipoles: the ET monopole $G_0$ in cubic systems and the ET quadrupole $G_u$ in noncubic systems. In this paper, we investigate achiral-to-chiral (AtC) structural phase transitions driven by atomic displacements and construct $G_{0,u}$ as explicit functions of the displacement order parameter $\eta$ based on a group-theoretical approach. We show that the leading-order dependence of $G_{0,u}(\eta)$ is determined by the symmetry of the parent structure and the character of the displacive mode, providing a symmetry-based classification of AtC transitions beyond a binary distinction between achiral and chiral phases. We also demonstrate that $G_{0,u}(\eta)$ is directly reflected in observable quantities such as CIM and chiral phonon splitting (CPS), both of which scale consistently with $G_{0,u}(\eta)$. We further clarify the microscopic mechanism by which AtC transitions give rise to chiral phonons and CPS through the coupling between $G_{0,u}(\eta)$ and phonon degrees of freedom.

We show that negative Gaussian curvature regularizes the infrared sector of crystalline membranes. In a covariant formulation of embedded elasticity, the Green strain contains a symmetry-required linear coupling between the normal phonon field and extrinsic curvature. Integrating out the in-plane phonons converts this coupling into a finite quadratic contribution to the inverse flexural response. Anomalous roughening is thereby replaced by curvature-controlled saturation, and the mechanism survives on minimal hyperbolic patches. Hyperbolic geometry preempts anomalous elasticity before the flat infrared regime is reached, implying the absence of a crumpling phase at the harmonic level. The same Gaussian-order coupling admits sound propagation in the infrared.

In conventional Friedel oscillations, the real-space charge density oscillations induced by an impurity are characterized by an oscillation period set by the Fermi momentum. In this work, we show that the conventional theory is incomplete when the Bloch wavefunctions carry nontrivial quantum geometry. We demonstrate that in metals with an isolated (nearly) flat band at the Fermi energy, quantum geometry induces a distinct type of oscillations, which we call the \emph{quantum geometric Friedel oscillations} (QGFOs). The period of the QGFOs is set by the momentum space separation of the quantum metric hot spots of the flat band. The conventional and quantum metric-induced oscillations coexist at low temperatures. At higher temperatures, the conventional Friedel oscillations away from the impurity site are set by the thermal length such that the oscillations can be easily washed out by temperature effects. Remarkably, the QGFOs decay length is set by the quantum metric length which is defined by the integration of the quantum metric of the flat band. As a result, the QGFOs can persist even at temperatures much larger than the bandwidth of the flat band. Moreover, the decay length is independent of temperature for a wide range of temperatures which is a manifestation of the quantum metric protection. In conclusion, we show that the quantum metric induces novel Friedel oscillations. Our work suggests that the measurement of the QGFOs is a powerful way to detect the quantum metric length (which is associated with the integral of the quantum metric) and the quantum metric hot spot separations (which are associated with the distribution of the quantum metric in the momentum space).

Accumulative polarization switching driven by sequential sub-coercive electric-field pulses offers a promising route toward low-power ferroelectric memories and neuromorphic devices. However, the kinetic regimes governing this nonequilibrium process remain poorly understood. Here, we employ a phase-field model based on the time-dependent Landau-Ginzburg formalism to investigate pulse-driven accumulative switching in ferroelectric HZO. By systematically varying the initial domain configuration, pulse amplitude, pulse-on time, and pulse-off time, we establish a quantitative link between microscopic domain-wall dynamics and macroscopic polarization accumulation. We show that the effective switched-domain radius follows distinct scaling regimes characterized by the local kinetic exponent. Initially, a local exponent greater than 1 indicates superlinear domain growth driven by enhanced irreversible domain-wall propagation under successive pulses. As switching progresses, a local exponent close to unity marks steady self-similar growth, whereas a local exponent less than 1 signifies decelerating dynamics caused by geometric confinement, depletion of switchable polarization, and relaxation-induced back switching. The transition between these regimes is governed by the competition between field-driven excitation during the pulse-on interval and spontaneous relaxation during the pulse-off interval. The initial domain geometry further influences this transition. Increasing the pulse amplitude or pulse-on duration extends the superlinear regime, whereas longer pulse-off times promote relaxation and suppress accumulation. These findings establish a unified scaling framework for pulse-driven accumulative switching, providing quantitative insight into nonequilibrium ferroelectric domain evolution and design guidelines for HZO-based memory and neuromorphic devices.

The microscopic origin of supercurrent nonreciprocity in van der Waals Josephson junctions remains under active debate, particularly regarding the role of twist-angle engineering in layered superconductors. Here, we investigate superconducting transport in twisted NbSe2/NbSe2 vertical Josephson junctions fabricated by dry transfer with controlled crystallographic alignment and chemically clean interfaces. High-resolution transmission electron microscopy is employed to directly determine the twist angle and assess interface quality. While the Josephson coupling strength exhibits a pronounced dependence on twist angle, with characteristic voltages maximized near crystallographically equivalent orientations and suppressed at intermediate angles, the supercurrent diode efficiency remains negligibly small and shows no systematic twist-angle dependence. In contrast, enhanced diode-like responses emerge only in weakly coupled junctions exhibiting interfacial disorder and reduced transparency. Deliberate interface degradation further amplifies the apparent nonreciprocity, yielding diode efficiencies approaching 30% together with an irregular magnetic-field-strength dependence. These results establish a clear decoupling between Josephson coupling and supercurrent nonreciprocity in twisted NbSe2/NbSe2 junctions. Our findings identify interface disorder, rather than twist-angle-controlled momentum matching, as the dominant origin of the observed diode response and provide a critical benchmark for interpreting nonreciprocal superconducting transport in van der Waals Josephson devices.

Traditional rheological tools are often limited in characterizing soft materials under ultra-high strain-rate loading conditions (> 1000 s^-1) due to constraints in spatiotemporal resolution, loading rate, and invasiveness. Recently, inertial microcavitation rheometry (IMR), which utilizes laser-induced inertial cavitation (LIC) to dynamically deform surrounding materials, has emerged as a powerful experimental mechanics technique for probing nonlinear viscoelastic properties under extreme loading conditions. However, conventional IMR relies on computationally expensive iterative inverse fitting procedures, limiting its scalability and real-time applicability. Here, we introduce a new AI-enhanced experimental mechanics framework, called Bubble Dynamics Transformer (BDT), that integrates physics-based cavitation simulations with Transformer neural network architectures to achieve rapid inverse characterization of soft material viscoelasticity from experimentally measured bubble dynamics. The proposed framework directly predicts viscoelastic material parameters from time-resolved bubble radius evolution curves without iterative optimization. The BDT is trained using synthetic datasets generated from physics-based Keller--Miksis cavitation simulations and validated using experimental laser-induced cavitation data obtained from hydrogels and viscous polymer solutions. The proposed AI-driven framework demonstrates excellent agreement with our previous IMR while substantially accelerating constitutive parameter inference. Experimental demonstrations further reveal the capability of the framework to characterize rate-dependent material behavior across a wide range of soft materials, from viscous liquids to various viscoelastic hydrogels, at ultra-high strain rates.

Entanglement entropy is a fundamental measure of quantum correlations and a key resource underpinning advances in quantum information and many-body physics. We uncover a universal relationship between bipartite entanglement entropy and particle number after the barrier in a one-dimensional Fermi-Hubbard system with an external asymmetric potential. Decomposing the von Neumann entropy into number entropy $S_n$ and configurational entropy $S_c$, we show that in the barrier-dominated tunneling regime both components are individually well-defined functions of the post-barrier particle density $n_A$, even though $S_c$ encodes off-diagonal coherences that are not directly accessible from density measurements alone. Using Kolmogorov-Arnold Networks - a novel machine learning architecture - we learn the relationship for entropy and its components across a broad range of interaction strengths and barrier heights with high predictive accuracy. Furthermore, we propose a simple analytical binary-entropy-like expression that quantitatively captures the observed correlation for fixed parameters. Our findings open new avenues for characterizing quantum correlations in transport phenomena and provide a powerful framework for estimating the full von Neumann entropy - including its configurational component - from a single transport observable.

We use a 2-categorical version of (de-)equivariantization to classify (3+1)d topological orders with a finite $G$-symmetry. In particular, we argue that (3+1)d fermionic topological order with $G$-symmetry correspond to $\mathbf{2SVect}$-enriched $G$-crossed braided fusion 2-categories. We then show that the categorical data necessary to define these theories agrees with that arising from a fermionic generalization of the Wang-Wen-Witten construction of bosonic topological theories with $G$-symmetry saturating an anomaly. More generally, we also explain how 2-categorical (de-) equivariantization yields a classification of all braided fusion 2-categories.

The extracellular polymeric substances (EPS) matrix of microbial biofilms exhibits a complex structural heterogeneity that profoundly influences mass transport and metabolic activity. Conventional biofilm models typically assume a homogeneous matrix, thereby neglecting the localized transport resistance introduced by the bacterial capsule, a distinct, low-diffusivity polysaccharide layer surrounding individual cells. In this theoretical study, we develop a multiscale "cell-capsule" continuum model that represents the capsule as a concentric shell enveloping each microbial cell core within the bulk EPS. Utilizing a one-dimensional reaction-diffusion framework coupled with a geometric characterization of capsule spacing and thickness, we quantify how microscale architecture modulates oxygen transport in developing biofilms. Model simulations demonstrate that incorporating a discrete capsular phase introduces a pronounced "resistance-in-series" effect, reducing local oxygen availability by up to 70% compared to conventional homogeneous models. Furthermore, our analysis indicates that capsule thickness and matrix compaction jointly control the effective diffusivity and oxygen effectiveness factor within the biofilm. These results provide critical mechanistic insights into how microscale organization governs macroscale biofilm function, offering a new framework for integrating structural heterogeneity into multiscale biofilm simulations.

The statistics of correlations are central quantities characterizing the collective dynamics of recurrent neural networks. We derive exact expressions for the statistics of correlations of nonlinear recurrent networks in the limit of a large number N of neurons, including systematic 1/N corrections, in the regime of Gaussian quenched disorder. Our approach uses a path-integral representation of the network stochastic dynamics, which reduces the description to a few collective variables and enables efficient computation. This generalizes previous results on linear networks to include a wide family of nonlinear activation functions, which enter as interaction terms in the path integral. These interactions can resolve the instability of the linear theory and yield a strictly positive participation dimension. We present explicit results for power-law activations, revealing scaling behavior controlled by the network coupling. In addition, we introduce a class of activation functions based on Pade approximants and provide analytic predictions for their correlation statistics. Numerical simulations confirm our theoretical results with excellent agreement. We also compare with previous works that have studied the complementary case with annealed disorder, and based on this we propose a new self-consistent equation for the more general case of colored noise.

Measurement-induced entanglement (MIE) captures how local measurements generate long-range quantum correlations and drive dynamical phase transitions in many-body systems. Yet estimating MIE experimentally remains challenging: direct evaluation requires extensive post-selection over measurement outcomes, raising the question of whether MIE is accessible with only polynomial resources. We address this challenge by reframing MIE detection as a data-driven learning problem that assumes no prior knowledge of state preparation. Using measurement records alone, we train a neural network in a self-supervised manner to predict the uncertainty metric for MIE--the gap between upper and lower bounds of the average post-measurement bipartite entanglement. Applied to random circuits with one-dimensional all-to-all connectivity, our method reveals a learnability transition with increasing circuit depth: below a threshold the MIE can be effectively learned with resources that grow only polynomially with system size, whereas above it the required resources grow exponentially. This computational phase transition coincides with the breakdown of efficient classical simulation of the underlying quantum state. We further observe signatures of this transition on current noisy quantum devices. These results highlight the power of data-driven approaches for learning MIE and delineate the practical limits of its classical learnability.

Silver iodide (AgI) thin films offer a compelling platform for studying nonlinear optical phenomena due to their intrinsic noncentrosymmetric lattice and direct band gap. Here, we investigate the nonlinear optical properties of AgI thin films grown by physical vapor deposition that selectively produce zincblende (\zbAgI) and wurtzite (\wzAgI) phases. Using a combination of polarization-resolved second harmonic generation (SHG) and two-photon photoluminescence (2PPL) spectroscopy, we identify clear phase- and morphology-dependent anisotropic nonlinear responses. Triangular \zbAgI $(111)$ flakes exhibit sixfold SHG symmetry and isotropic 2PPL emission, while rod-shaped \wzAgI $(101)$ samples display twofold-symmetric patterns in both SHG and 2PPL, which are explained by theories that integrate explicit polarization analysis using second- and third-order nonlinear susceptibilities. We estimate $\chi^{(2)}_\mathrm{eff}$ of 0.45 pm/V (\zbAgI) and 0.16 pm/V (\wzAgI), respectively, after correcting for multiple internal reflections and generation of SHG within the AgI film and optical interference effect in the dielectric layers. These results identify AgI as a useful single-composition halide platform for studying phase-dependent nonlinearity and establish a realistic methodology for evaluating nonlinear susceptibilities in layered materials or thin films supported by a substrate.

The motion of atoms in programmable optical tweezer arrays offers many new opportunities for neutral atom quantum science. These include inter- and intra-site atom motion for resource-efficient implementations of fermionic and bosonic modes, respectively, as well as tweezer transport for efficient compilation of arbitrary circuits. However, the exploitation of atomic motion for all three purposes and others is limited by the inertia of the atoms. We present a comprehensive architectural blueprint for the use of fermionic metastable helium-3 ($^3$He$^*$) atoms -- the lightest trappable atomic species -- in programmable optical tweezer arrays. This includes a concrete analysis of atomic structure considerations as well as Rydberg-mediated interactions. We show that inter-tweezer hopping of $^3$He$^*$ atoms can be $\gtrsim3\times$ faster than previous demonstrations with lithium-6. We also demonstrate a new toolbox for encoding and manipulating qubits directly in the tweezer trap potential, uniquely enabled by the light mass of $^3$He$^*$. Finally, we provide several examples of new opportunities for fermionic quantum simulation and computation that leverage the transport and inter-tweezer hopping of $^3$He$^*$ atom arrays. These tools present new methods to improve the resource efficiency of neutral atom quantum science that may also enable quantum simulations of lattice gauge theories and quantum chemistry outside the Born-Oppenheimer approximation

In quantum metrology, precision is typically characterized by an ensemble-averaged quantity, the quantum Fisher information (QFI), which averages over the fluctuations of individual measurement records. Here we introduce the conditional quantum Fisher information (CQFI), a trajectory-level version of the QFI that generalizes the classical stochastic Fisher information to the quantum domain. Defined through the symmetric logarithmic derivative and conditioned on a measurement outcome, the CQFI is a random variable whose average recovers the QFI. Using it, we derive a trajectory-level quantum speed limit, illustrated by the quantum-jump unraveling of a driven thermal qubit. Moreover, the CQFI decomposes into incoherent (population) and coherent (basis-rotation) contributions, together with an interference cross-term. This cross-term vanishes on average but can take negative values along single trajectories, providing a local witness of destructive interference between classical and quantum information channels.

Using the modified Villain lattice Hamiltonian formulation of the 1+1d compact boson theory, we construct exactly solvable abelian chiral lattice gauge theories in two spacetime dimensions. As a concrete example, we derive an explicit quadratic lattice Hamiltonian for the "34-50" chiral gauge theory. We further show that $N$ copies of the modified Villain theory realize the $O(N,N;\mathbb{Z})$ T-duality transformations, which we then use to solve and analyze these lattice gauge theories.

We calculate the correction exponents in the chiral Heisenberg model in the $1/N$ expansion. These exponents are related to the slopes of $\beta$ functions at the phase transition point. We present the results at order $1/N^2$ and check that they agree with the results of the $\epsilon$ expansion near $d = 4$. We find that one of the correction exponents diverges as $d \to 3$. We argue that the appearance of the pole is a rather general phenomenon and is associated with operator mixing involving the system of four-fermion operators. After analyzing the operator mixing structure, we propose a resummation procedure which modifies the exponents already at leading order. We also perform calculations directly in the three-dimensional model and find complete agreement with the resummed exponents.

Motivated by ongoing interest in the universal behaviour of the Hubbard model of spinning electrons on honeycomb and $\pi$-flux lattices at the semi-metal -- Mott insulator phase transition, we formulate the 3D chiral Heisenberg model, a theory of relativistic fermions in three spacetime dimensions, as a lattice field theory using domain wall fermions. The contact interaction term preserves an SU(2) global symmetry. We perform numerical simulations using the Rational Hybrid Monte Carlo algorithm on system sizes $L^3\times L_s$ with $L\in\{8,\ldots,24\}$ and domain wall separation $L_s\in\{8,16,24\}$. We locate the phase transition corresponding to spontaneous SU(2)$\to$U(1) breaking, yielding critical exponent estimates $\nu^{-1}=0.63(3)$, $\eta_\Phi=1.42(8)$. These values are considerably removed from estimates obtained from simulations performed in (2+1)D, ie. with the time and spatial directions treated differently, but align more closely with analytic estimates obtained using 3D covariant field theory. We also present first results for the fermion correlator, ultimately needed for the determination of the exponent $\eta_\Psi$, highlighting the need to rotate the fermion source to a common reference direction in isospace in order to obtain a signal.

Dijkgraaf-Witten theories have a wide range of applications in topological phases of matter and the study of generalized global symmetries. We develop a method to construct BF-type Lagrangians for Dijkgraaf-Witten theories with non-abelian gauge group by gauging $H^{(0)}$ symmetries from a BF-Lagrangian of an abelian Dijkgraaf-Witten theory. When $H$ nontrivially permutes the operators of the original theory, the Lagrangian of the $H$-gauged theory is constructed with cohomologies with local coefficients. We analyze the structure of the Lagrangians and their gauge transformations with homotopy theory. We also construct the operator spectrum and verify the Lagrangians by matching elementary linking invariants.

Causality -- that a response cannot precede its cause -- is among the most universal constraints in physics. Yet when a unitary microscopic theory is reduced to an open-system memory-kernel equation, causality is not inherited for free: an upper-half-plane pole of the kernel forces exponential growth of the reduced propagator and is incompatible with any completely positive trace-preserving (CPTP) reduced dynamics. We examine three objects with different causal structure under projection -- the Nakajima-Zwanzig (NZ) memory kernel $\tilde{\mathcal{K}}(z)$, the reduced-state Laplace transform $\tilde{\sigma}(z)$, and the force-fit effective kernel $\tilde{\mathcal{K}}_\text{eff}(z)$ -- using the Kramers-Kronig (KK) relations as diagnostic. Under a real-axis spectral hypothesis on the projected generator, $\tilde{\mathcal{K}}(z)$ lies in a vector-valued Hardy space and obeys (subtracted) KK relations, giving a CPTP-consistency criterion, a passivity-analyticity statement, and a Carleman diagnostic. We prove $\tilde{\sigma}(z)$ analytic in the upper half-plane for any initial state -- unitarity bounds $\|\sigma(t)\|_{\mathrm{op}} \leq 1$, so acausality cannot be blamed on the state alone. Yet the force-fit kernel can develop upper-half-plane poles at simple zeros of $\tilde{\sigma}(z)$: passive baths sit in a robust regime where these zeros stay real, while near-resonant systems enter a fragile regime in which coherence-channel zeros migrate into the upper half-plane, an intrinsic symmetry property already present for factorized states. We verify the full operator-valued KK relation on the extracted $4\times4$ NZ memory kernel of the Jaynes-Cummings model, the relative $L_2$ residual decreasing under refinement ($3.8\% \to 0.95\%$), consistent with exact matrix-valued KK in the continuum limit.

Photocatalytic oxidation (PCO) is a promising strategy for indoor air purification and outdoor pollutant abatement, potentially offering treatment for climate- and health-relevant pollutants such as methane (CH$_4$), nitrogen oxides (NO$_\text{x}$) and volatile organic compounds (VOCs). In this work, we present experiments evaluating the PCO of CH$_4$ (2 to 10 ppm) under varying UV-C light intensities (4 to 59 W/m$^2$), using titanium dioxide (TiO$_2$) as the photocatalyst. At 2 ppm CH$_4$, TiO$_2$ achieves a maximum conversion efficiency of 24.4% and a maximum apparent quantum yield of $0.013$% over the tested UV-C light intensities, demonstrating activity at environmentally relevant concentrations. We develop a model to interpret the experimental results and assess the potential of PCO for ventilation applications. The model is validated against our CH$_4$ data and literature results for formaldehyde (HCHO) and NO$_\text{x}$. While laboratory-scale configurations achieve high conversions (e.g., 24.4% for CH$_4$), ventilation-scale performance is predicted to be limited by thin concentration boundary layers and short residence times, with conversion efficiencies dropping to around $0.017$%. Finally, we estimate the climate impact of CH$_4$ removal in terms of CO$_2$e emission rates, demonstrating that TiO$_2$-based PCO in ventilation applications can yield a net climate benefit (i.e., a net-negative CO$_2$e emissions rate) when the modelled CO$_2$e removal rate exceeds the emissions from catalyst material production and UV operation, particularly when pre-existing UV-C irradiation is leveraged.