Papers

There are multiple possible origins of transport anisotropies in metals and superconductors. For instance, rotational symmetry can be spontaneously broken in the normal state as a result of electronic nematic order inducing anisotropies in an otherwise $s$-wave superconducting phase. Another possibility is that the dominant source of rotational symmetry breaking is the superconductor itself and its vestiges that may survive in the normal state. We here theoretically analyze the correlations of transport anisotropies in the normal and the corresponding superconducting phase for different scenarios of broken symmetry, either coming solely from the normal state, solely from the superconductor and its vestiges in the metallic regimes, or from both simultaneously. We further include both zero-momentum and finite-momentum pairing; we develop a theory of vestigial order for the latter, characterized by broken rotational and translational symmetry. Our findings reveal that the relative transport anisotropies in the normal and superconducting phases sensitively depend on the scenario, including the form of vestigial order and, in some cases, the parity of the superconducting order parameter. As such, measuring the directional dependence of the critical current and resistivity can provide strong constraints on the origin of rotational symmetry breaking. We demonstrate our findings in minimal models relevant to twisted multilayer graphene, rhombohedral graphene, and twisted transition metal dichalcogenides.

Condensed matter physics continues to seek new frustrated quantum materials that not only deepen our understanding of fundamental physical phenomena but also hold promise for transformative technologies. In this review article, we highlight the unique features of chiral spin topology and review the topological phenomena recently identified in trillium lattice compounds. Based on the unique spin states realized in these systems, we explore the potential for realizing various theoretically proposed chiral quantum phases. We examine representative materials including the magnetic insulating compound K2Ni2(SO4)3 and and the intermetallic EuPtSi discussing both experimental findings and theoretical predictions, while outlining several key questions. Finally, we offer a perspective on promising research directions aimed at uncovering novel emergent behavior in chiral trillium lattice-based materials.

Our work establishes the mathematical equivalence between polymer Self Consistent Field Theory and interfacial SAFT based classical Density Functional Theory by providing an explicit proof that SCFT is the mean field (MF) or saddle point limit of a Legendre dual formulation of iSAFT CDFT. The polymer CDFT ideal chain term is the Legendre conjugate of the single chain generating functional built from standard SCFT propagators, so the Legendre transform converts the iSAFT CDFT Euler Lagrange equations into SCFT fixed point relations. Within this Functional Field Theory, Wertheim TPT1 association free energy can be incorporated into the excess functional, with the mass action variables providing additive SCFT fields yet leaving the propagators unchanged. Linearization about the MF solution produces a reversible mapping for fluctuation corrections, suggesting a universal response that couples the single chain correlation kernel to the iSAFT direct correlation operator, allowing the inclusion of short range structure and association cooperatively beyond mean field. Long range interactions map similarly, showing how functional derivatives can yield SCFT field contributions, while second derivatives yield the long range blocks of the iSAFT correlation operator, providing a unified treatment of polymers. Mathematically, this framework aims to clarify that SCFT and polymer DFT solvers share similar diffusion based propagators, and how iSAFT hard sphere, attraction, and association modules contribute additional field and hessian terms. Specifically, this work provides a Legendre duality proof that places SAFT CDFTs and SCFT on the same MF foundation, and a reversible theoretical framework that transfers iSAFT correlation information into SCFT, enabling modular, quantitative modeling of associating, structured polymer systems.

We propose geometrically nonlinear (finite) continuum models of flexomagnetism based on the Cosserat micropolar and its descendent couple-stress theory. These models introduce the magneto-mechanical interaction by coupling the micro-dislocation tensor of the micropolar model with the magnetisation vector using a Lifshitz invariant. In contrast to conventional formulations that couple strain-gradients to the magnetisation using fourth-order tensors, our approach relies on third-order tensor couplings by virtue of the micro-dislocation being a second-order tensor. Consequently, the models permit centrosymmetric materials with a single new flexomagnetic constant, and more generally allow cubic-symmetric materials with two such constants. We postulate the flexomagnetic action-functionals and derive the corresponding governing equations using both scalar and vectorial magnetic potential formulations, and present numerical results for a nano-beam geometry, confirming the physical plausibility and computational feasibility of the models.

Topological phases in narrow-gap materials often arise from a delicate interplay of spin-orbit coupling (SOC) and electron correlation, making their prediction sensitive to the underlying electronic-structure approximation. Electronic topology is commonly characterized using effective one-body Hamiltonians, but this approach becomes delicate near gap closings and does not reveal how a band inversion is realized in the interacting spectrum. Here we introduce a multi-messenger analysis that combines three complementary readouts of a single interacting Green's function: the Fu-Kane invariant of the topological Hamiltonian, the symmetry-resolved spectral function, and the symmetry-resolved orbital occupation. Together, they identify the topological class, resolve band connectivity, and track the orbital character of the inversion. We implement this framework within fully self-consistent relativistic GW, treating SOC and electronic correlation on equal footing while eliminating dependence on a density-functional starting point. Applied to strain-tunable $\alpha$-Sn, the approach distinguishes its competing zero-gap, Dirac-semimetal, and topological-insulator regimes and provides an internally consistent description of their common band-inversion mechanism. More broadly, this framework offers a frequency-resolved route to diagnosing correlated topology across both gapped and near-critical states.

Reactive dopant atoms embedded in inert host metal surfaces define the active sites in single-atom alloys (SAAs), yet SAA synthesis remains challenging. To address this, we elucidate how dopant adatoms deposited on Cu and Ag surfaces become incorporated into the metal and identify periodic trends from early to late transition metals (TMs) using density functional theory. Adatoms diffuse nearly freely across terraces, as diffusion barriers are small, whereas direct incorporation into terraces is unfavourable. In line with conventional wisdom, step edges and kink sites strongly facilitate dopant incorporation, confirming their critical role in alloy formation. Attachment of adatoms to steps and kinks from the lower terrace is favoured. Incorporation then proceeds either from this attached state or when adatoms approach a step edge from above, where reactions often proceed without barrier. Incorporation barriers are generally lower for early and central TMs, increase towards late TMs, and are slightly higher on Cu than on Ag surfaces. Repulsive interactions between Pd adatoms and dopants explain the experimental observation that a dopant-rich brim on the upper terrace of Cu surfaces inhibits incorporation from above. In contrast, attractive interactions, as found for Ru, anchor diffusing adatoms (even on terraces) and promote the formation of adatom islands, yet hinder incorporation next to the dopant and may impede the growth of embedded dopant clusters. By rationalising periodic trends and experimental observations, we show how specific surface sites and adatom–dopant interactions shape dopant incorporation, offering guidance on the surface environments most conducive to SAA synthesis for different dopant elements.

The directed motion of active colloids is governed by spatial variations in surface chemistry and interfacial stress, yet these properties remain extremely difficult to measure directly. We introduce a physics-informed neural network framework that infers the slip distribution driving propulsion from partial observations of the surrounding flow. By combining sparse fluid velocity measurements with the Stokes equations and boundary constraints, the method reconstructs both the near-surface slip and the full velocity and pressure fields. Validation against analytical solutions and Boundary Element Method calculations for canonical active colloid models shows quantitative agreement in both unbounded and confined geometries. Crucially, the framework recovers the surface slip even when no flow data are available near the particle, demonstrating that accessible bulk measurements encode the interfacial stresses responsible for active motion. These results establish physics-informed inference as a powerful tool for characterizing and ultimately controlling interfacially driven transport in colloidal active matter.

This review provides an introduction to the essential physics of soft adhesion, including the thermodynamics of adhesion and wetting, the mechanics of contact with deformable materials, and the material properties that most affect interfacial interactions with soft solid gels and elastomers. Throughout, we emphasize both foundational physics and current experimental and theoretical research in these areas. We conclude with a practical overview of standard experimental test methods for characterizing soft adhesion. The physical understanding developed herein provides the basis for understanding the mechanics of contact with soft materials.

We demonstrate the emergence of a robust tetrahedral magnetic ground state in monolayer graphene doped to the van Hove singularity (vHS). This noncoplanar, gapped spin configuration - featuring four equally inclined moments - has been previously identified as a candidate instability. Here, not only do we confirm its stability across all finite interactions using fully self-consistent, real-space-resolved calculations, but we also go beyond earlier work by charting the full surrounding phase diagram. In doing so, we unravel a cascade of symmetry-broken magnetic states - pseudo-tetrahedral, planar, collinear, and modulated textures - which we classify using spin structure factors and vector order parameters. These results stem from unrestricted Hartree-Fock simulations on large supercells with dense k-point sampling, enabling us to resolve interaction-driven magnetic and charge inhomogeneities. Our findings connect directly with recent ARPES and doping experiments near the vHS in graphene, and establish the tetrahedral state as the central correlated instability in this regime, offering predictive insight into emergent magnetism in correlated Dirac materials.

Topologically protected Majorana zero modes (MZMs) have attracted intense interest due to their potential application in fault-tolerant quantum computation (TQC). Graphene nanoribbons, with tunable edge terminations and compatibility with planar device architectures, offer a promising alternative to semiconductor nanowires. Here we present a comprehensive theoretical study of finite graphene strips with armchair, zigzag, and nearly square geometries, proximitized by an s-wave superconductor and subject to Rashba spin-orbit coupling, Zeeman fields, and disorder. Using exact diagonalization of the Bogoliubov-de Gennes tight-binding Hamiltonian, we analyze Majorana polarization, low-energy spectra, and real-space wavefunctions to identify the non-trivial topological phases supporting MZMs and distinguish them from from partially separated Andreev bound states (psABS) or the quasi-Majoranas. We systematically chart the robustness of these modes across geometries and disorder regimes, finding that armchair strips with short zigzag edges provide the most stable platform. Our results unify polarization diagnostics with spatial wavefunction analysis and disorder effects, yielding concrete design guidelines for graphene-based topological superconductors.

We present a historical review of anyon and exclusion statistics, introduced in the 1980s and 1990s respectively, and then turn to developments in the recently introduced inclusion statistics. In contrast to exclusion statistics, where particles tend to be more exclusive than usual fermions, inclusion statistics particles tend to be more gregarious than usual bosons and manifest an enhanced propensity to form condensates. Inclusion and exclusion statistics are related through a duality transformation, generalizing the well-known Bose-Fermi duality. We conclude with a review of the Calogero model realization of exclusion statistics and its extension to inclusion statistics.

We uncover a class of nonlinear odd viscoelastic effects in three spatial dimensions. We show that these dissipationless effects arise upon combining geometric distortions in two orthogonal directions, yielding momentum flow in the third direction. We demonstrate that the effect arises from nontrivial geometric tensors in quantum states, and can be scaled up with integer topological invariants. We further show that the effect fingerprints the multiband Hilbert-space geometry of the underlying quantum states, as encoded in the nonmetricity and three-state quantum geometric tensors. Our findings unravel the role of multistate geometry in viscoelastic phenomena, paving a path for experimental observation of uncharted nonlinear odd viscoelastic responses in quantum systems.

To study the impact of active systems on their surroundings, we introduce a model that couples an active nematic fluid to an isotropic substrate fluid via friction. We numerically show that as the active layer develops turbulence, the substrate inherits the chaotic behaviour, exhibiting a novel form of turbulence driven by locally generated stochastic forcing from the active layer. In particular, the short-length-scale flow structures in the active layer are filtered out, so the system behaves as a de facto low-pass filter. We derive analytically the transfer function between the two layers and use it to predict the large-q decay of the substrate energy spectrum, and to investigate how tensorial quantities, such as the strain rate and the active stresses, are transmitted between the active layer and the substrate. Our analysis agrees with recent experiments measuring velocity-velocity correlations in mixtures of active and passive microtubules, and it may have implications for traction force microscopy measurements in cellular layers.

In recent years tantalizing signs for a novel phase have been reported that is chirally symmetric but nevertheless exhibits massive bound states. The necessary condition for such a phase, referred to as Symmetric Mass Generation (SMG), is the cancellation of all (continuous and discrete) 't~Hooft anomalies. In 3+1 dimensions this occurs in systems containing a multiple of 16 massless Weyl fermions. SMG was originally discovered in lower dimensional condensed matter systems. We present results investigating four dimensional field theories with gauge group SU(3). Our findings suggest that SU(3) with $N_f=8$ fundamental fermions exhibits an SMG phase not only on the lattice but also in the infinite cutoff continuum limit. If confirmed, SMG could provide a new UV completion of the standard model and give rise to new scenarios for beyond standard model physics.

This work presents a finite element method for simulating dynamic processes that involve the coupled evolution of dislocation motion and crack propagation. The method numerically solves the Concurrent Atomistic-Continuum (CAC) formulation of the conservation of linear momentum. A crystalline material is discretized at the unit-cell level using 6-node prism elements whose geometry allows dislocations and cracks to nucleate and propagate along element facets. Nanoscale simulations of single-crystal Cu, Fe, and Si demonstrate the initiation and propagation of dislocations and cracks, and these results are reproduced by the finite element method in excellent agreement with fully atomistic molecular dynamics simulations. Mesoscale simulations of single-crystal Cu further demonstrate the ability of the method to capture size-dependent brittle and ductile behavior. Under plane-strain conditions the Cu model fractures in a brittle manner, while a fully three-dimensional model exhibits curved and intersecting dislocations that blunt the crack tip and prevent crack propagation, resulting in ductile behavior. The accuracy, efficiency, and applicability of the method are discussed.

Infinite projected entangled-pair states (iPEPS) provide a powerful tensor network ansatz for two-dimensional quantum many-body systems in the thermodynamic limit. In this paper we introduce an approach to accurately compute the energy variance of an iPEPS, enabling systematic extrapolations of the ground-state energy to the exact zero-variance limit. It is based on the contraction of a large cell of tensors using the corner transfer matrix renormalization group (CTRMG) method, to evaluate the correlator between pairs of local Hamiltonian terms. We show that the accuracy of this approach is substantially higher than that of previous methods, and we demonstrate the usefulness of variance extrapolation for the Heisenberg model, for a free fermionic model, and for the Shastry-Sutherland model. Finally, we apply the approach to compute $\langle \langle \mathcal{L}^\dagger \mathcal{L} \rangle \rangle$ for an open quantum system described by the Liouvillian $\mathcal{L}$, in order to assess the quality of the steady-state solution and to locate first-order phase transitions, using the dissipative quantum Ising model as an example.

Out-of-time-order correlators (OTOCs) are central probes of quantum scrambling, and their generalizations have recently become key primitives for both benchmarking quantum advantage and learning the structure of Hamiltonians. Yet their behavior has lacked a unified algorithmic interpretation. We show that higher-order OTOCs naturally fit within the framework of quantum signal processing (QSP): each $\mathrm{OTOC}^{(k)}$ measures the $2k$-th Fourier component of the phase distribution associated with the singular values of a spatially resolved truncated propagator. This explains the contrasting sensitivities of time-ordered correlators (TOCs) and higher-order OTOCs to causal-cone structure and to chaotic, integrable, or localized dynamics. Based on this understanding, we further generalize higher-order OTOCs by polynomial transformation of the singular values of the spatially resolved truncated propagator. The resultant signal allows us to construct frequency-selective filters, which we call OTOC spectroscopy. This extends conventional OTOCs into a mode-resolved tool for probing scrambling and spectral structure of quantum many-body dynamics.

We investigate the phonon magnetochiral effect (PMCE) in chiral Weyl semimetals by deriving the nonreciprocal dispersion relations of both acoustic and non-polar optical phonons in the presence of a magnetic field. Using a semiclassical Boltzmann kinetic framework that incorporates Berry curvature, orbital magnetic moment, and node-dependent electronic structure, we obtain analytic expressions for the magnetic-field-induced corrections to the phonon dynamical matrix. Inequivalent Weyl nodes with distinct Fermi velocities, Fermi energies, and relaxation times generate a dynamical chiral imbalance that alters the phonon dispersion. For acoustic phonons, the formalism yields the magnetic-field-dependent corrections to the longitudinal mode, while for optical phonons we identify an optical analogue of the PMCE that produces a corresponding shift in the optical branch. Together, these results provide a unified theoretical description of how band-geometric properties of Weyl fermions influence both acoustic and optical phonon dispersions in chiral Weyl semimetals.

Understanding structure-property relationships in materials is fundamental in condensed matter physics and materials science. Over the past few years, machine learning (ML) has emerged as a powerful tool for advancing this understanding and accelerating materials discovery. Early ML approaches primarily focused on constructing and screening large material spaces to identify promising candidates for various applications. More recently, research efforts have increasingly shifted toward generating crystal structures using end-to-end generative models. This review analyzes the current state of generative modeling for crystal structure prediction and de novo generation. It examines crystal representations, outlines the generative models used to design crystal structures, and evaluates their respective strengths and limitations. Furthermore, the review highlights experimental considerations for evaluating generated structures and provides recommendations for suitable existing software tools. Emerging topics, such as modeling disorder and defects, integration in advanced characterization, incorporating synthetic feasibility constraints, and model explainability are explored. Ultimately, this work aims to inform both experimental scientists looking to adapt suitable ML models to their specific circumstances and ML specialists seeking to understand the unique challenges related to inverse materials design and discovery.

Recent progress in generalized symmetry and topological holography has shown that, in conformal field theory (CFT), topological data from one dimensional higher can play a key role in determining local dynamics. Based on this insight, a fixed-point (FP) tensor complex (TC) for CFT has recently been constructed. In this work, we develop a TC renormalization (TCR) algorithm adapted to this CFT-based structure, forming a renormalization-group (RG) framework with generalized symmetry. We show that the full FP tensor can emerge from the RG flow starting with only the three-point function of the primary fields. Remarkably, even when starting solely from topological data, the RG process can still reconstruct the full FP tensor–a method we call as topological bootstrap. This approach deepens the connection between the topological and dynamical aspects of CFT and suggests pathways toward a fully algebraic description of gapless quantum states, with potential extensions to higher dimensions.

We study large deviations, over a long time window $T \to \infty$, of the dynamical observables $A_n = \int_{0}^{T} x^n(t) dt$, $n=3,4,\dots$, where $x(t)$ is a centered stationary Gaussian process in continuous time. We show that, for short-correlated processes the probability density of $A_n$ exhibits an anomalous scaling $P(A_n,T) \sim \exp[-T^{\mu} f_n(\Delta A_n T^{-\nu})]$ at $T\to \infty$ while keeping $\Delta A_n T^{-\nu}$ constant. Here $\Delta A_n$ is the deviation of $A_n$ from its ensemble average. The anomalous exponents $\mu$ and $\nu$ depend on $n$ and are smaller than $1$, whereas the rate function $f_n(z)$ exhibits a first-order dynamical phase transition (DPT) which resembles condensation transitions observed in many systems. The same type of anomaly and DPT, with the same $\mu$ and $\nu$, was previously uncovered for the Ornstein-Uhlenbeck process - the only stationary Gaussian process which is also Markovian. We also uncover an anomalous behavior and a similar DPT in the long-correlated Gaussian processes. However, the anomalous exponents $\mu$ and $\nu$ are determined in this case not only by $n$ but also by the power-law long-time decay $\sim |t|^{-\alpha}$ of the covariance. The different anomalous scaling behavior is a consequence of a faster-than-linear scaling with $T$ of the variance of $A_n$. Finally, for sufficiently long-ranged correlations, $\alpha<2/n$, the DPT disappears, giving way to a smooth crossover between the regions of typical, Gaussian fluctuations and large deviations. The basic mechanism behind the DPT is the existence of strongly localized optimal paths of the process conditioned on very large $A_n$ and coexistence between the localized and delocalized paths of the conditioned process. Our theoretical predictions are corroborated by replica-exchange Wang-Landau simulations where we could probe probability densities down to $10^{-200}$.

We give a short proof that low-temperature dynamics for the Sherrington-Kirkpatrick model have mixing time exponential in the system size, based on the recently proved existence of gapped spin configurations by (Minzer-Sah-Sawhney 2023, Dandi-Gamarnik-Zdeborová 2023). This result is in contrast with a well established physics prediction which posits a stretched exponential mixing time of order $e^{N^{1/3 \pm o(1)}}$. Our proof clarifies that this prediction cannot apply to mixing from worst case initial conditions, but should presumably be understood to concern dynamics from a suitably random initialization.

High-Temperature Superconductors (HTS) such as YBa2Cu3O7-delta (YBCO) are essential for next-generation Tokamak fusion reactors, where Rare-Earth Barium Copper Oxides (REBCO) form the functional layers in HTS magnets. Because YBCO's superconductivity depends strongly on oxygen stoichiometry and defect structure, atomistic simulations can provide crucial insight into radiation-damage mechanisms and pathways to maintain material performance. In this work, we develop and benchmark four Machine-Learned Interatomic Potentials (MLPs) for YBCO: the Atomic Cluster Expansion (ACE), the Message-Passing Atomic Cluster Expansion (MACE), the Gaussian Approximation Potential (GAP), and the Tabulated Gaussian Approximation Potential (tabGAP), trained on an extensive Density Functional Theory (DFT) database explicitly designed to include irradiation-damaged-like configurations. The resulting models achieve DFT-level accuracy across a wide range of atomic environments, faithfully capturing the interatomic forces relevant to radiation damage processes. Among the tested models, MACE delivers the highest accuracy, although at greater computational cost, while ACE and tabGAP provide an excellent balance between efficiency and fidelity. These machine-learned potentials establish a robust foundation for large-scale molecular dynamics simulations of radiation-induced defect evolution in complex superconducting materials

In living systems, DNA undergoes continuous and rhythmic mechanical remodeling through condensation, looping, and disentangling to regulate gene expression, segregate chromosomes, and guide morphogenesis. Here, we demonstrate a purely mechanical route to rhythmic DNA reorganization in a minimal active composite of microtubules, kinesin motors, and DNA. We embed a DNA polymer in an active turbulent microtubule-kinesin fluid, creating a self-morphing material. The active flows stretch and entangle the DNA, forming a self-organized viscoelastic network that resists active stresses and affects flow over large length scales. This mechanical feedback loop progressively amplifies velocity correlations and drives a nonequilibrium phase transition tuned by DNA contour length: from disordered flow to synchronized, millimeter-scale oscillations with vortices. We rationalize the phase transition with an active-gel model that predicts a growing length scale and an oscillatory instability emerging from the interplay between activity, orientational order, and self-generated viscoelasticity, rather than chemical signaling. The dependence of the oscillation frequency on system size and activity quantitatively agrees with experiment. Thus, flow-driven DNA remodeling provides a minimal physical route to autonomous, system-spanning oscillations in three dimensions and suggests design principles for programmable soft matter that coordinates, actuates, and reshapes itself.

The Kibble-Zurek mechanism (KZM) successfully predicts the density of topological defects deposited by the phase transitions, but it is not clear why. Its key conjecture is that, near the critical point of the second-order phase transition, critical slowing down will result in a period when the system is too sluggish to follow the potential that is changing faster than its reaction time. The correlation length at the freeze-out instant $\hat t$ when the order parameter catches up with the post-transition broken symmetry configuration is then decisive, determining when the mosaic of broken symmetry domains locks in topological defects. To understand why the KZM works so well we analyze Landau-Ginzburg model and show why temporal evolution of the order parameter plays such a key role. The analytical solutions we obtain suggest novel, hitherto unexplored, experimentally accessible observables that can shed light on symmetry breaking dynamics while testing the conjecture on which the KZM is based.

Charge Kondo circuits consist of metallic islands connected by single-mode quantum point contacts (QPCs). The island's charging energy makes these circuits tunable quantum simulators of various strongly interacting models. Here we propose a circuit that realizes the Kondo effect with effective Luttinger-liquid interactions, and show that it undergoes a localization transition in which the QPC transmission is fully suppressed below a critical value. Experimental signatures include a diverging charge susceptibility and an entropy step. Our findings open a path toward realizing localization transitions in more exotic settings.

Photon-pressure circuits are the circuit implementation of the cavity optomechanical Hamiltonian and discussed for qubit readout, low-frequency quantum photonics and dark matter axion detection. Due to the enormous design flexibility of superconducting circuits, photon-pressure systems provide fascinating possibilities to explore unusual parameter regimes of the optomechanical Hamiltonian. Here, we report the realization of a photon-pressure platform, in which a GHz circuit interacts with a MHz circuit via a magnetic-flux-tunable combination of dispersive and dissipative photon-pressure. In addition, both coupling rates are considerably enhanced by nonlinearities of the GHz-mode, which leads to the multi-photon coupling rates scaling stronger with the pump photon number $n_\mathrm{c}$ than the usual $\sqrt{n_\mathrm{c}}$ dependence. We demonstrate that interference of the two interaction paths leads to a Fano-like response in photon-pressure induced transparency, and that the dynamical backaction is considerably modified compared to the dispersive case, including a parametric instability caused by a red-detuned pump tone.

The spin Hall magnetoresistance (SMR) is widely used to study the interplay between charge and spin currents in bilayers of a magnetic insulator and a normal metal. However, not much is known about the spatial variation of the SMR across the surface of one and the same sample. In this work, we investigate the statistical distribution of the SMR in hundreds of nominally identical Hall bar structures patterned into prototypical yttrium iron garnet (YIG)/Pt heterostructures. We find a Gaussian-distributed SMR with a narrow standard deviation of approximately 10$\,$% of the mean value in each YIG/Pt bilayer studied. However, the variation of the mean SMR between different YIG/Pt samples can be as large as ~30$\,$%, despite nominally identical fabrication conditions. This demonstrates that spatial variations of the SMR amplitude must not be neglected, in particular when comparing different heterostructures. On a microscopic level, local variations of the interface quality captured by the spin mixing conductance are the most likely origin for the observed SMR amplitude variations.

The combination of polystyrene particles and multilayers with chiral interfacial interactions offers an interesting platform to study the interplay between topology, chirality, and curvature. Here, we investigate Pt/Co/Ta multilayers sputtered onto self-assembled polystyrene particles with a diameter of 919 nm, which show a magnetic reversal distinctly different from that of planar films. Using non-invasive scanning transmission X-ray microscopy and in-vacuum magnetic force microscopy (MFM) with low-moment tips, we observe maze domains in the planar films, whereas the curvilinear films host three-dimensional spiraling stripe states. Repeated scanning with a conventional higher-moment MFM tip ruptures these spirals into metastable skyrmions at the top of the particles. Micromagnetic simulations of spherical caps reveal that increasing curvature renders a single spiral stripe energetically favorable over an isolated skyrmion, and that this stabilization is predominantly of magnetostatic origin. Our results identify the spiral states as intrinsic to the curvilinear geometry and the skyrmions as metastable states nucleated by local magnetic fields, underlining the value of non-invasive imaging to distinguish intrinsic curvature effects from probe-induced magnetic textures.

We explore the magnetic properties of Nd$_2$Zr$_2$O$_7$ and Pr$_2$Zr$_2$O$_7$ single crystals subjected to pulsed magnetic fields up to 60 T using magnetization and magnetocaloric-effect (MCE) measurements, with initial temperatures ranging from 2 to 31K. The MCE data exhibit pronounced and unconventional hysteresis loops, in which the sample temperature increases during both the up-sweep and down-sweep of the field. In Nd$_2$Zr$_2$O$_7$, the MCE further displays a striking plateau as a function of time, followed by a rapid temperature rise that begins at the maximum applied field, across pulses with differing peak-field strengths. Our magnetization measurements reveal an inferred temperature of the magnetic subsystem that differs significantly from the directly measured sample temperature and exhibits opposite hysteresis: the temperature is higher on the up-sweep than the down-sweep, unlike the direct measurements. These observations indicate a breakdown of thermal equilibrium between magnetic and lattice degrees of freedom on the timescale of the pulse ($\sim 10^{-1}$s). We interpret the results using a phenomenological model involving two thermally coupled subsystems - the magnetic ions and phonons, and a thermal reservoir, which accounts well for the behavior of Pr$_2$Zr$_2$O$_7$. However, it fails to reproduce the plateau seen in Nd$_2$Zr$_2$O$_7$. Agreement with Nd$_2$Zr$_2$O$_7$ data is improved substantially if we allow the thermal coupling between the magnetic and the lattice subsystems to depend on the product $\frac{HdH}{dt}$. Our results reveal anomalously slow heat transfer between magnetic and lattice subsystems and point toward a novel mechanism for dynamically controlling the heat flow in Nd$_2$Zr$_2$O$_7$ via the rate of magnetic field variation.