Papers

$\alpha$-RuCl$_3$ is a leading material for proximate Kitaev magnetism. We address the origin of the broad, $\Gamma$-point centered excitation continuum observed by inelastic neutron scattering at elevated temperatures in this compound. Using stochastic Landau-Lifshitz dynamics augmented with quantum-equivalent corrections, we reproduce the temperature-dependent dynamical spin structure factor across both the correlated and conventional paramagnetic regimes. A meta-analysis of 38 published exchange parameter sets identifies those most consistent with the full temperature evolution. A Bayesian optimization procedure is used to derive parameters that capture the low-energy star-like momentum dependence and the overall bandwidth of the continuum. Rescaling temperatures by the Curie–Weiss scale produces a collapse of spectral measures, demonstrating that the high-$T$ dynamics are governed by correlated paramagnetism below $\theta_{\mathrm{CW}}$ rather than by the Kitaev crossover to fractionalized excitations. Complementary 24-site exact diagonalization clarifies finite-size systematics at low temperature and the proximity to zigzag/incommensurate ordering. Beyond $\alpha$-RuCl$_3$, our simulation pipeline provides a reproducible, data-driven framework to infer effective spin models in magnets that exhibit broad continua.

The Kramers-Wannier transformation of the 1+1d transverse-field Ising model exchanges the paramagnetic and ferromagnetic phases and, at criticality, manifests as a non-invertible symmetry. Extending such self-duality symmetries beyond gauging of abelian groups in tensor-product Hilbert spaces has, however, remained challenging. In this work, we construct a generalized 1+1d Ising model based on a finite-dimensional semisimple Hopf algebra $H$ that enjoys an anomaly-free non-invertible symmetry $\mathrm{Rep}(H)$. We provide an intuitive diagrammatic formulation of both the Hamiltonian and the symmetry operators using a non-(co)commutative generalization of ZX-calculus built from Hopf-algebraic data. When $H$ is self-dual, we further construct a generalized Kramers-Wannier duality operator that exchanges the paramagnetic and ferromagnetic phases and becomes a non-invertible symmetry at the self-dual point. This enlarged symmetry mixes with lattice translation and, in the infrared, flows to a weakly integral fusion category given by a $\mathbb{Z}_2$ extension of $\mathrm{Rep}(H)$. Specializing to the Kac-Paljutkin algebra $H_8$, the smallest self-dual Hopf algebra beyond abelian group algebras, we numerically study the phase diagram and identify four of the six $\mathrm{Rep}(H_8)$-symmetric gapped phases, separated by Ising critical lines and meeting at a multicritical point. We also realize all six $\mathrm{Rep}(H_8)$-symmetric gapped phases on the lattice via the $H$-comodule algebra formalism, in agreement with the module-category classification of $\mathrm{Rep}(H_8)$. Our results provide a unified Hopf-algebraic framework for non-invertible symmetries, dualities, and the tensor product lattice models that realize them.

We provide explicit constant-depth local unitary circuits that realize general anyon permutations in Kitaev's quantum double models. This construction can be naturally understood through a correspondence between anyon permutation symmetries of two-dimensional topological orders and self-dualities in one-dimensional systems, where local gates implement self-duality transformations on the boundaries of microscopic regions. From this holographic perspective, general anyon permutations in the $D(G)$ quantum double correspond to compositions of three classes of one-dimensional self-dualities, including gauging of certain subgroups of $G$, stacking with $G$ symmetry-protected topological phases, and outer automorphisms of the group $G$. We construct circuits realizing the first class by employing self-dual unitary gauging maps, and present transversal circuits for the latter two classes.

Altermagnetism (AM) and its associated spin-transport phenomena are typically linked to spin-split electronic band structures in bulk materials. However, the crystal surface has a reduced symmetry with respect to the bulk, which can induce AM at the surface of conventional antiferromagnets (AFMs) $\unicode{x2013}$ a local effect which cannot be detected using bulk properties. In this work we define the symmetry conditions necessary for surface AM and show how it can be topologically protected, rendering it a robust effect. We provide a minimal model for one trivial and two topological examples of surface AM. We show that the spin spectral density, accessible by spin- and angle-resolved photoemission spectroscopy, can exhibit a $d$-wave-like altermagnetic character at the surface, even when the full band structure is completely spin degenerate. Our topological model describes the Dirac semimetal CuMnAs, which provides an existing realization of our theory. Our results identify crystal surfaces as a platform to realize robust, topology- and symmetry-driven unconventional magnetism beyond the bulk classification of magnetic materials.

Motivated by the recent discovery of superconductivity in the bilayer nickelate La$_3$Ni$_2$O$_7$ (LNO) under pressure, we study a mixed-dimensional (mixD) bilayer $t_\parallel$-$J_\parallel$-$J_\perp$ model, which has been proposed as an effective low-energy description of LNO. Using neural quantum states (NQS), and in particular Gutzwiller-projected Hidden Fermion Pfaffian State, we access the ground-state properties on large lattices up to $8\times 8\times 2$ sites. We show that this model exhibits superconductivity across a wide range of dopings and couplings, and analyze the pairing behavior in detail. We identify a crossover from tightly bound, Bose-Einstein-condensed interlayer pairs at strong interlayer exchange to more spatially extended Bardeen-Cooper-Schrieffer-like pairs as the interlayer exchange is decreased. Furthermore, upon tuning the intralayer exchange, we observe a sharp transition from interlayer $s$-wave pairing to intralayer $d$-wave pairing, consistent with a first-order change in the pairing symmetry. We verify that our simulations are accurate by comparing with matrix product state simulations on coupled ladders. Our results represent the first simulation of a fermionic multi-orbital system with NQS, and provide the first evidence for superconductivity in two-dimensonal bilayers using high-precision numerics. These findings provide insight into superconductivity in bilayer nickelates and cold atom quantum simulation platforms.

In enzymatic recycling or biodegradation of semi-crystalline plastic waste, crystalline spherulites embedded into an amorphous matrix hinder and slow down depolymerisation. When the enzymatic depolymerisation temperature exceeds the glass transition temperature, these spherulites tend to grow. The depolymerisation process is thus controlled by a competition between erosion of the amorphous matrix from the particle surface and the growth of recalcitrant spherulites within the particle bulk and at its surface. We present a geometric model that captures this competition, together with an algorithm to solve the equations numerically. Our algorithm introduces a new extension of Voronoi/Delaunay tessellation in space. We extract the parameters for the model from experimental data on the enzymatic depolymerization by hydrolase LCC-ICCG of PET bottle flakes and textile waste, in order to make a prediction of the observed degradation yield as a function of time. Both the final yield and the degradation kinetics are correctly predicted. Most importantly, the model clarifies how and to which extent nucleating agents, impurities, additives, and/or rapid crystal growth present in the waste can undermine pretreatment efforts aiming to initiate depolymerisation from a material with a low initial crystallinity.

In this paper we apply the canonical spin wave theory developed for itinerant Kondo lattice magnets in the strong coupling regime to Kondo ferromagnets, and address two general questions pertaining to their magnetic excitations. First, we compute corrections to the strong coupling (i.e., double-exchange) spin wave dispersion of itinerant ferromagnets. We show that the spin wave dispersion beyond the strong coupling limit can be mapped to the spin wave dispersion of a Heisenberg ferromagnet with farther neighbor exchange couplings, and discuss how this affects instabilities towards antiferromagnetism. Second, we examine the effect of including electronic spin-orbit coupling in the spin wave theory of Kondo ferromagnets. Including spin-orbit coupling is natural and straightforward in the formulation of the canonical spin wave expansion. Our key result is to demonstrate that the linear spin wave Hamiltonian of the itinerant Kondo ferromagnet can be mapped to the spin wave Hamiltonian of a Heisenberg ferromagnet with easy-axis Ising anisotropy and antisymmetric Dzyaloshinskii-Moriya exchange interaction. We show that in the case of the Kane-Mele honeycomb lattice Kondo ferromagnet this leads to topological magnon bands, and discuss the implications of this result for itinerant ferromagnets more broadly.

The distinctive symmetry properties of pure altermagnets make them natural candidates for piezomagnetism. Previous work motivated by the piezomagnetic properties of altermagnets has primarily focused on the spin magnetization response to applied strain. In this paper we study orbital piezomagnetic effects–the orbital magnetization response to applied strain–in minimal lattice models of pure insulating altermagnets in two dimensions. We obtain general microscopic expressions for the orbital magnetization in the presence of strain, as well as the orbital piezomagnetic polarizability, i.e., the defining response characteristic of pure altermagnets. We apply these expressions to three specific tetragonal lattice models, two corresponding to $d$-wave altermagnets and one describing a $g$-wave altermagnet. Whereas the $d$-wave altermagnets are associated with a linear piezomagnetic polarizability, the $g$-wave altermagnet exhibits a nonlinear piezomagnetic effect. Our analysis reveals how the polarizabilities are related to and determined by the Berry curvature of the occupied bands. Connections to materials of current interest are discussed.

Vibrational spectroscopy in the electron microscope can reveal phonon excitations with nanometer spatial resolution, yet routine prediction remains out of reach due to fragmented workflows requiring specialized expertise. Here we introduce PySlice, the first publicly available implementation of the Time Autocorrelation of Auxiliary Wavefunction (TACAW) method, providing an automated framework that produces momentum- and energy-resolved vibrational electron energy-loss spectra directly from atomic structures. By integrating universal machine learning interatomic potentials with TACAW, PySlice eliminates the bottleneck of per-system potential development. Users input atomic structures and obtain phonon dispersions, spectral diffraction patterns, and spectrum images through a unified workflow spanning molecular dynamics, GPU-accelerated electron scattering, and frequency-domain analysis. We outline the formulation behind the code, demonstrate its application to canonical systems in materials science, and discuss its use for advanced analysis and materials exploration. The modular Python architecture additionally supports conventional electron microscopy simulations, providing a general-purpose platform for imaging and diffraction calculations. PySlice makes vibrational spectroscopy prediction routine rather than specialized, enabling computational screening for experimental design, systematic exploration of phonon physics across materials families, and high-throughput generation of simulated data for training of future machine learning models.

The percolation phase transition in complex network systems attracts much attention and has numerous applications in various research fields. Finite size effects smooth the transition and make it difficult to predict the critical point of appearance or disappearance of the giant connected component. For this end, we introduce the susceptibility of arbitrary random undirected and directed networks and show that a strong increase of the susceptibility is the early warning signal of approaching the transition point. Our method is based on the introduction of `observers', which are randomly chosen nodes monitoring the local connectivity of a network. To demonstrate efficiency of the method, we derive explicit equations determining the susceptibility and study its critical behavior near continuous and mixed-order phase transitions in uncorrelated random undirected and directed networks, networks with dependency links, and $k$-cores of networks. The universality of the critical behavior is supported by the phenomenological Landau theory of phase transitions.

We introduce a special kind of nonreciprocal electromagnetic response which gives rise to backscattering suppression in the bulk, a long-sought feature in topological photonics, as well as nonreciprocal lensing - an effect when the same structure focuses light incident from one direction and defocuses light propagating in the opposite way. We predict this response in spin spirals and in specially designed metamaterials, validating the key predictions.

In black-box combinatorial optimization, objective evaluations are often expensive, so high quality solutions must be found under a limited budget. Factorization machine with quantum annealing (FMQA) builds a quadratic surrogate model from evaluated samples and optimizes it on an Ising machine. However, FMQA requires binary decision variables, and for nonbinary structures such as integer permutations, the choice of binary encoding strongly affects search efficiency. If the encoding fails to reflect the original neighborhood structure, small Hamming moves may not correspond to meaningful modifications in the original solution space, and constrained problems can yield many infeasible candidates that waste evaluations. Recent work combines FMQA with a binary autoencoder (bAE) that learns a compact binary latent code from feasible solutions, yet the mechanism behind its performance gains is unclear. Using a small traveling salesman problem as an interpretable testbed, we show that the bAE reconstructs feasible tours accurately and, compared with manually designed encodings at similar compression, better aligns tour distances with latent Hamming distances, yields smoother neighborhoods under small bit flips, and produces fewer local optima. These geometric properties explain why bAE+FMQA improves the approximation ratio faster while maintaining feasibility throughout optimization, and they provide guidance for designing latent representations for black-box optimization.

We present a multiscale simulation framework to compute the current vs. voltage (I-V ) characteristics of metal/oxide/metal structures building the core of conductive bridging random access memory (CBRAM) cells and to shed light on their resistance switching properties. The approach relies on a finite element model whose input material parameters are extracted either from ab initio or from machine-learned empirical calculations. The applied techniques range from molecular dynamics and nudged elastic band to electronic and thermal quantum transport. Such an approach drastically reduces the number of fitting parameters needed and makes the resulting modeling environment more accurate than traditional ones. The developed computational framework is then applied to the investigation of an Ag/a-SiO2/Pt CBRAM, reproducing experimental data very well. Moreover, the relevance of Joule heating is assessed by considering various cell geometries. It is found that self-heating manifests itself in devices with thin conductive filaments with few-nanometer diameters and at current concentrations in the 10s-microampere range. With the proposed methodology it is now possible to explore the potential of not-yet fabricated memory cells and to reliably optimize their design.

Among the recognized sources of decoherence in superconducting qubits, the spatial inhomogeneity of the superconducting state and the possible presence of magnetic-flux vortices remain comparatively underexplored. Niobium is commonly used as a structural material in transmon qubits that host Josephson junctions, and excess dissipation anywhere in the transmon can become a bottleneck that limits overall quantum performance. The metal/substrate interfacial layer may simultaneously host pair-breaking loss channels (e.g., two-level systems, TLS) and control thermal transport, thereby affecting dissipation and temperature stability. Here, we use quantitative magneto-optical imaging of the magnetic-flux distribution to characterize the homogeneity of the superconducting state and the critical current density, $j_{c}$, in niobium films fabricated under different sputtering conditions. The imaging reveals distinct flux-penetration regimes, ranging from a nearly ideal Bean critical state to strongly nonuniform thermo-magnetic dendritic avalanches. By fitting the measured magnetic-induction profiles, we extract $j_{c}$ and correlate it with film physical properties and with measured qubit internal quality factors. Our results indicate that the Nb/Si interlayer can be a significant contributor to decoherence and should be considered an important factor that must be optimized.

Simulating low-temperature properties of three-dimensional frustrated quantum magnets is challenging due to the sign problem and the system sizes required to mitigate substantial finite-size effects. However, there are many experimental examples of three-dimensional crystals that could host exotic low-temperature states of matter, such as quantum spin liquids. We calculate the ground-state phase diagrams of frustrated quantum spin models on the body-centered cubic lattice using neural quantum states. First, we study the antiferromagnetic $J_1-J_2$ model where we find a direct first-order phase transition between Néel and collinear long-range-ordered phases at $(J_2/J_1)_c = 0.705$, consistent with previous studies. Then, in a tetragonally-distorted variant, proposed as a minimal model of NaCa$_2$Cu$_2$(VO$_4$)$_3$, we find no evidence of a quantum paramagnetic ground state, with a first-order phase transition between Néel and chain phases at $(J_{2ab}/J_1)_c = 1.0375$. Therefore, the ground state of the tetragonally-distorted model does not reproduce the low-temperature magnetic properties of NaCa$_2$Cu$_2$(VO$_4$)$_3$, and the inclusion of other effects is necessary to rationalize experimental observations.

Thermodynamics is commonly presented as a theory of macroscopic systems in stable equilibrium, built upon assumptions of extensivity and scaling with system size. In this paper, we present a universal formulation of the elementary foundations of thermodynamics, in which entropy and energy are defined and employed beyond equilibrium and without assuming extensivity. The formulation applies to all systems – large and small, with many or few particles – and to all states, whether equilibrium or nonequilibrium, by relying on carefully stated operational definitions and existence principles rather than macroscopic idealizations. Key thermodynamic concepts, including adiabatic availability and available energy, are developed and illustrated using the energy-entropy diagram representation of nonequilibrium states, which provides geometric insight into irreversibility and the limits of work extraction for systems of any size. A substantial part of the paper is devoted to the analysis of entropy transfer in non-work interactions, leading to precise definitions of heat interactions and heat-and-diffusion interactions of central importance in mesoscopic continuum theories of nonequilibrium behavior in simple and complex solids and fluids. As a direct consequence of this analysis, Clausius inequalities and the Clausius statement of the second law are derived in forms explicitly extended to nonequilibrium processes. The resulting framework presents thermodynamics as a universal theory whose concepts apply uniformly to all systems, large and small, and provides a coherent foundation for both teaching and modern applications.

Chemical substitution in crystalline quantum materials is a powerful way to explore the consequences of strong spin-orbit coupling on their structural and electronic properties. In this work, we present an investigation of thin films of the La$\textit{Pn}_2$ ($\textit{Pn}$~=~Sb, Bi) class of layered square-net intermetallics. We report the growth of LaBi$_2$ with a pristine layer-by-layer growth mode, classifying it as a good metal displaying superconductivity at $\sim$0.55~K. Compared to LaSb$_2$, we attribute the enhanced metallic behavior and improved growth dynamics of LaBi$_2$ to significant relativistic corrections to its electronic band structure and the resulting impact on both surface energy and intrinsic phonon scattering.

Moiré superlattices in two-dimensional (2D) materials exhibit rich quantum phenomena, but ab initio modelling of these systems remains computationally prohibitive. Existing machine learning methods for accelerating density-functional theory (DFT) can target the prediction of different quantities and often rely on the locality assumption. Here we train a Gaussian process regression SALTED model exclusively on the electron densities of small displaced bilayer structures and then extrapolate electron density prediction to the large supercells required to describe small twist angles between these bilayers. We show the necessity of long-range descriptors to yield reliable band structures and electrostatic properties of large twisted bilayer structures, when these are derived from predicted densities. We demonstrate that the choice of descriptor determines the distribution of residual density errors, which in turn affects the downstream electronic properties. We apply our models to twisted bilayer graphene, hexagonal boron nitride, and transition metal dichalcogenides, focusing on the model's capacity to predict complex phenomena, including flat band formation, bandwidth narrowing, domain-wall electric fields, and spin-orbit coupling effects. Beyond moiré materials, this approach provides a general methodology for electronic structure prediction in large-scale systems with substantial long-range phenomena related to non-local geometric information.

Chirality in active and passive fluids gives rise to odd transport properties, most notably the emergence of robust edge currents that defy standard dissipative dynamics. While these phenomena are well-described by continuum hydrodynamics, a microscopic framework connecting them to their topological origins has remained elusive. Here, we present a lattice model for an isotropic chiral random walk that bridges the gap between classical stochastic diffusion and unitary quantum evolution. By equipping the walker with an internal degree of freedom and a tunable chirality parameter, $p$, we interpolate between a standard diffusive random walk and a deterministic, topologically non-trivial quantum walk. We show that the topological protection characteristic of the unitary limit ($p=1$) remarkably persists into the dissipative regime ($p<1$). This correspondence allows us to theoretically ground the robustness of edge flows in classical chiral systems using the bulk-boundary correspondence of Floquet topological insulators. Our results provide a discrete microscopic description for odd diffusion, offering a powerful toolkit to predict transport in confined geometries and disordered chiral media.

Spin qubits are typically operated in the lowest orbital of a quantum dot to minimize interference from nearby states. In valence-band hole systems, strong spin-orbit coupling links spin and orbital degrees of freedom, strongly influencing the hole $g$-factor, a key parameter for qubit control. We investigate the out-of-plane $g$-factor in Ge quantum dots using excitation (single-particle) and addition (many-body) spectra. Excitation spectra allow us to distinguish the pure Zeeman $g$-factor from orbital contributions to the magnetic field splitting of states despite the strong spin-orbit coupling. This distinction clarifies discrepancies between $g$-factors extracted with the two methods, for different orbital states and different hole numbers. Furthermore, we find gate-tunability of $g$-factors at the level of 15%, highlighting its relevance for all-electric qubit manipulation.

We consider a clean Josephson-junction chain coupled by one of its extremities to a thermal bath through a resistance. Considering the Langevin dynamics in the classical regime, in the case of Josephson energy much smaller than charging energy, we find that heat propagates logarithmically slowly through the system, rather than diffusively, as highlighted by the logarithmic increase in time of a thermalization length we define and by the logarithmically slow increase in time of the energy. This behavior – typical of quantum Anderson or many-body localized systems – is observed here also in a clean classical glassy Hamiltonian system. We argue that this phenomenon might imply strong robustness to the effect of ergodic inclusions for the nonergodic behavior in the charge-quantized regime.

Adiabatic chiral state conversion (CSC) is one of the many counterintuitive effects associated with non-Hermitian physics. In quantum systems, numerous works have demonstrated this phenomenon under both non-Hermitian Hamiltonian and Lindblad evolution. However, despite considerable progress, the physical mechanism behind it has been a subject of debate. In this work, we present a unified framework that explains CSC in any non-Hermitian system, encompassing non-Hermitian Hamiltonian, Lindblad, and hybrid settings. Our framework relies on perturbative, non-adiabatic corrections to adiabatic evolution and consistently predicts CSC with only the lowest-order corrections. We demonstrate its efficacy with models of single and coupled dissipative qubits, obtaining analytical solutions for the conversion fidelity. Our analysis further reveals the role of non-perturbative dynamics, which can be present even in apparently slow trajectories. We show that this property can be utilised to considerably enhance state conversion. Finally, we demonstrate that CSC can be observed in a model without the presence of exceptional points.

During the development of an organism, cells must coordinate and organize to generate the correct shape, structure, and spatial patterns of tissues and organs, a process known as morphogenesis. The morphogenesis of embryonic tissues is supported by multiple processes that induce the precise physical deformations required for tissues to ultimately form organs with complex geometries. Among the most active players shaping the morphogenetic path are fine-tuned changes in cell adhesion. We review here recent advances showing that changes on cell adhesion, a local, pair-wise property defined at the cell-cell contact level has important global consequences for embryonic tissue topology, being determinant in defining both the geometric and material properties of early embryo tissues.

Conduction band-edge spin-orbit splitting (SOS) in monolayer transition metal dichalcogenides determines a competition between bright and dark excitons and sets conditions for spintronics applications of these semiconductors. Here, we report the SOS measurement for electrons in monolayer MoS$_2$, found from the threshold density, $n_*$, for the upper spin-orbit-split band population, which exceeds by an order of magnitude the values expected from the conventional density functional theory (DFT). Theoretically, half of the observed value can be attributed to the exchange enhancement of SOS in a finite-density electron gas, but explaining the rest requires refining the DFT approach. As the conduction band SOS in MoS$_2$ is set by a delicate balance between the contribution of sulphur $p_x$ and $p_y$ orbitals and $d_{z^2}-d_{xz}$ and $d_{z^2}-d_{yz}$ mixing in molybdenum, we use a DFT+U+V framework for fine-tuning the orbital composition of the relevant band-edge states. An optimised choice of Hubbard U/V parameters produces close agreement with the experimentally observed conduction band SOS in MoS$_2$, simultaneously resulting in the valence-band SOS and the quasi-particle band gap which are closer to their values established in the earlier-published experiments.

${\rm{EuTe}}_{4}$ is a polar charge density wave (CDW) material, with giant thermal hysteresis and non-volatile state switching under electric and optical fields, attracting great attention in recent years. However, the in-depth understanding of these anomalous phenomena remains elusive. Herein, via first-principles calculations, we reveal that the polar CDW state in ${\rm{EuTe}}_{4}$ hosts a novel layer-selective nature, wherein multiple energetically close CDW configurations coexist and exhibit low interconversion energy barriers. Monte Carlo simulations indicate that the giant thermal hysteresis in ${\rm{EuTe}}_{4}$ originates from a phase transition mainly driven by the change of configurational entropy, around which the material hosts a metastable CDW state characterized by diverse local polar configurations breaking the out-of-plane translational symmetry. The configurational composition of this metastable CDW state can be effectively controlled by electric and optical fields, thereby enabling non-volatile state switching. Our theoretical findings align well with recent experimental observations in ${\rm{EuTe}}_{4}$ and pave the way for exploring the emerging phenomena and applications of polar CDW in multilayered systems.

We compute exciton binding energies and fine-structure splittings in CdSe nanoplatelets with two zincblende geometries and one wurtzite geometry, finding that the wurtzite structure exhibits the largest bright-bright splitting due to its intrinsic in-plane anisotropy, while the zincblende structures show smaller but finite splittings arising from atomistic symmetry breaking at edges and corners. These results are obtained using a theoretical framework that we developed, which combines DFT single-particle states with screened configuration interaction, a quasi-2D dielectric screening model, and an efficient Coulomb-cutoff scheme that eliminates periodic-image interactions and enables accurate Coulomb and exchange integrals at low computational cost. This methodology provides a transferable and practical route for studying excitons in CdSe nanoplatelets and other quasi-two-dimensional nanomaterials.

Schelling's model of segregation demonstrates that even in the absence of social or governmental interventions, individuals with mild in-group preferences can self-organize into strongly segregated neighborhoods. Many variants of this celebrated model have been proposed by assuming agents tend to increase their satisfaction. Complementary to this traditional, utility-based approach, we model residential moves using satisfaction-independent reaction rates in a spatially extended chemical reaction network. The resulting model exhibits an emergent phenomenon: despite symmetric in-group preferences, the system undergoes a tipping transition at a critical preference level, beyond which one agent type dominates. We characterize this asymmetric phase transition in details using mean-field analysis, numerical simulations and finite size scaling methods. We find that while the transition shares key features with the Ising universality class, such as $\mathbb{Z}_2$ symmetry breaking and similar exponent ratios, the full set of critical exponents does not match any known universality class.

Topological invariants govern many important physical properties in condensed matter systems. In this work, we obtain the complete set of topological invariants for a family of one-dimensional quasicrystals. The first and best-studied member of the family is the Fibonacci chain, while the successive ones are known in the literature as silver, bronze... and collectively as the metallic mean chains. By considering rational approximants, and by making use of the relationship between these chains and two dimensional Quantum Hall problems, we write down a gap labeling scheme for finite systems, and extend it to the quasiperiodic limit. We show, by numerical computations on open chains, that the proposed scheme correctly yields the winding numbers of edge states in each of the gaps, in all of the quasicrystals. We observe that their relationship to the 2D model leads, in the strict 1D chains family, to spectra forming a simplified Hofstadter “butterfly" diagram, with the analogues of Landau levels appearing in the asymptotic limit. These topological invariants for noninteracting 1D chains can be expected to remain pertinent for superconducting systems where Majorana or Aharanov bound states occur, as shown in recent studies of interacting or Josephson coupled Fibonacci chains.

Three-dimensional electron diffraction (3D ED) has emerged as a powerful method for solving the structures of sub-micron-sized particles down to nanoparticles. However, it faces technical challenges when applied to beam-sensitive samples or agglomerated nanoparticles. This study presents a novel approach that combines 4D-STEM tomography with object tracking and segmentation algorithms to overcome these limitations and achieve single-crystalline 3D ED datasets from nanopowder samples. The method and data quality are assessed on brookite TiO2 nanorods and beam-sensitive CsPbBr3 nanoparticles. To finely sample the reciprocal-space, the data acquisition was automated to acquire hundreds of 4D-STEM scans at fine tilt steps using a slightly convergent beam (0.6 to 1 mrad). Compared to conventional 3D ED methodologies, the proposed method provides enhanced signal-to-noise ratio, low illumination time for reducing beam damage, and the ability to analyze multiple particles from a single tomographic dataset. The procedure is optimized to be feasible using commercially available desktops and detectors. This extends the method applicability to systems and samples that were very challenging for conventional 3D ED methods, by eliminating several technical challenges for the data acquisition.

We present the thermopower of EuCd2P2, a material which exhibits a large resistivity peak with significant magnetic field dependence in the temperature range of 10-25 K. In the same region we observe a highly unusual behavior of the thermopower with two sign changes and giant extrema. The overall variation of the thermopower exceeds 4 000 muV/K and takes place in an extremely narrow temperature region of less than 5 K. The anomaly is suppressed completely in a small magnetic field of 0.5 T. We discuss this observation using a simple drift-diffusion picture and taking into account that the temperature gradient inducing the thermopower voltage is accompanied by a gradient of the electrical resistivity. Our simple estimation yields the correct magnitude, shape, and field dependence of the thermopower anomaly observed in EuCd2P2. These results open a new route to giant thermopower values via gradients of electronic properties.