2025-06-19 | | Total: 82
We calculate the spin correlations that result from the motion of a single dopant in the hard-core Fermi-Hubbard model, as the geometry evolves from a square to a triangular lattice. In particular, we consider the square lattice with an additional hopping along one diagonal, whose strength is continuously varied. We use a high-temperature expansion which expresses the partition function as a sum over closed paths taken by the dopant. We sample thousands of diagrams in the space of closed paths using the quantum Monte Carlo approach of Raghavan and Elser [1,2], which is free of finite-size effects and allows us to simulate temperatures as low as T∼0.3|t|, even in cases where there is a sign problem. For the case of a hole dopant, we find a crossover from kinetic ferromagnetism to kinetic antiferromagnetism as the geometry is tuned from square to triangular, which can be observed in current quantum gas microscopes.
The discovery of local-moment magnetism in van der Waals (vdW) semiconductors down to the single-layer limit has led to a paradigm shift in the understanding of two-dimensional (2D) magnets and unleashed their potential for applications in microelectronic and optoelectronic devices. The incorporation of strong electronic and magnetic correlations in 2D vdW metals remains a sought-after platform not only to enable control of emergent quantum phases, such as superconductivity, but also to achieve more theoretically tractable microscopic models of complex materials. To date, however, there is limited success in the discovery of such metallic vdW platforms, and f-electron monolayers remain out of reach. Here we demonstrate that the actinide \beta-UTe_3 can be exfoliated to the monolayer limit. A sizable electronic specific heat coefficient provides the hallmark of strong correlations. Remarkably, \beta-UTe_3 remains ferromagnetic in the half-unit-cell limit with an enhanced ordering temperature of 35 K, a factor of two larger than its bulk counterpart. Our work establishes \beta-UTe_3 as a novel materials platform for investigating and modeling correlated behavior in the monolayer limit and opens numerous avenues for quantum control with, e.g., strain engineering.
We propose a Cooper pair analogue of Dzyaloshinskii-Moriya's anisotropic coupling for d-vectors in an inhomogeneous unconventional superconductor, which can arise when itinerant electrons couple to a local exchange field of frustrated spin moments. Using perturbative Green's function methods, we examine an s-d model on a geometrically frustrated lattice with noncollinear spin texture that breaks time-reversal and parity symmetries, leading to spin-triplet pairing correlations. The effective Josephson tunneling in the presence of such an exchange field leads to an anisotropic Josephson coupling between d-vectors, which prefers a spatially inhomogeneous pairing order. Furthermore, we discuss anomalous vortices arising from the spatially inhomogeneous d-vector in the absence of a magnetic field for nonunitary triplet pairing order. We also propose a Josephson diode effect in which the diode efficiency is proportional to the spin chirality of the local exchange field. Our results can help understand the experiments in superconducting proximitized Mn_3Ge and in 4Hb-TaS_2.
When dense granular gases are continuously excited under microgravity conditions, spatial inhomogeneities of the particle number density can emerge. A significant share of particles may collect in strongly overpopulated regions, called clusters. This dynamical clustering, or gas-cluster transition, is caused by a complex interplay and balance between the energy influx and dissipation in particle collisions. Particle number density, container geometry, and excitation strength influence this transition. We perform Discrete Element Method (DEM) simulations for ensembles of frictional spheres in a cuboid container and apply the Kolmogorov Smirnov test and a caging criterion to the local packing fraction profiles to detect clusters. Machine learning can be used to study the gas-cluster transition, and can be a promising alternative to identify the state of the system for a given set of system parameters without time-consuming complex DEM simulations. We test various machine learning models and identify the best models to predict dynamical clustering of frictional spheres in a specific experimental geometry.
Accurate free-energy estimation is essential in molecular simulation, yet the periodic boundary conditions (PBC) commonly used in computer simulations have rarely been explicitly exploited. Equilibrium methods such as umbrella sampling, metadynamics, and adaptive biasing force require extensive sampling, while non-equilibrium pulling with Jarzynski's equality suffers from poor convergence due to exponential averaging. Here, we introduce a physics-informed, score-based diffusion framework: by mapping PBC simulations onto a Brownian particle in a periodic potential, we derive the Fokker-Planck steady-state score that directly encodes free-energy gradients. A neural network is trained on non-equilibrium trajectories to learn this score, providing a principled scheme to efficiently reconstruct the potential of mean force (PMF). On benchmark periodic potentials and small-molecule membrane permeation, our method is up to one order of magnitude more efficient than umbrella sampling.
Disorder, though naturally present in experimental samples and strongly influencing a wide range of material phenomena, remains underexplored in first-principles studies due to the computational cost of sampling the large supercell and configurational space. The recent development of machine learning techniques, particularly graph neural networks (GNNs), has enabled the efficient and accurate predictions of complex material properties, offering promising tools for studying disordered systems. In this work, we introduce a computational framework that integrates GNNs with Monte Carlo simulations for efficient calculations of thermodynamic properties and ensemble-averaged functional properties of disordered materials. Using the surface-termination-disordered MXene monolayer \ch{Ti3C2T}_{2-x} as a representative system, we investigate the effect of surface termination disorder involving \ch{-F}, \ch{-O}, and termination vacancies on the electrical and optical conductivity spectra. We find that surface termination disorder affects the temperature dependence of electrical conductivity, inducing a peak close to the order-disorder phase transition temperature that reflects the competition between scattering and electron filling effects of the surface termination groups across the phase transition. In contrast, optical conductivity remains robust to local disorder across a wide temperature range and is governed primarily by the global chemical composition of surface terminations. These results demonstrate the utility of our machine-learning-assisted framework for statistically modeling disorder effects and ensemble properties in complex materials, opening new avenues for future studies of disorder-driven phenomena in systems such as high-entropy alloys and disordered magnetic compounds.
Quantifying space exploration is a central question in random walk theory, with direct applications ranging from animal foraging, diffusion-limited reactions, and intracellular transport to stock markets. In particular, the explored domain by one or many simple memoryless (or Markovian) random walks has received considerable attention . However, the physical systems mentioned above typically involve significant memory effects, and so far, no general framework exists to analyze such systems. We introduce the concept of a \emph{flip}, defined most naturally in one dimension, where the visited territory is [x_{\rm min}, x_{\rm max}]: a flip occurs when, after discovering a new site at x_{\rm max}, the walker next discovers x_{\rm min} - 1 instead of x_{\rm max} + 1 (and vice-versa). While it reduces to the classical splitting probability in Markovian systems, we show that the flip probability serves as a key observable for quantifying the impact of memory effects on space exploration. Here, we demonstrate that the flip probability follows a strikingly simple and universal law: it decays inversely with the number of sites visited, as 1/n, independently of the underlying stochastic process. We confirm this behavior through simulations across paradigmatic non-Markovian models and observe it in real-world systems, without relying on model assumptions, including biological tracer motion, DNA sequences, and financial market dynamics. Finally, we reveal the physical mechanism behind this universality and show how it extends to higher-dimensional and fractal domains. Our determination of universal flip statistics lay the groundwork for understanding how memory effects govern random explorations.
In recent decades, it has been emphasized that the evolving structure of networks may be shaped by interaction principles that yield sparse graphs with a vertex degree distribution exhibiting an algebraic tail, and other structural traits that are not featured in traditional random graphs. In this respect, through a mean-field approach, this review tackles the statistical physics of graph models based on the interaction principle of duplication-divergence. Additional sophistications extending the duplication-divergence model are also reviewed as well as generalizations of other known models. Possible research gaps and related prior results are then discussed.
We demonstrate that the machine learning of density functionals allows one to determine simultaneously the equilibrium chemical potential across simulation datasets of inhomogeneous classical fluids. Minimization of an implicit loss function based on an Euler-Lagrange equation yields both the universal one-body direct correlation functional, which is represented locally by a neural network, as well as the system-specific unknown chemical potential values. The method can serve as an efficient alternative to conventional computational techniques of measuring the chemical potential. It also facilitates using canonical training data from Brownian dynamics, molecular dynamics, or Monte Carlo simulations as a basis for constructing neural density functionals, which are fit for accurate multiscale prediction of soft matter systems in equilibrium.
To explore and quantitatively map the cation-size mismatch solubility limits in high-entropy oxides (HEOs), we report on Ca^{2+} substitution in prototypical MgCoNiCuZnO, because, while isovalent, Ca^{2+} is 38% larger than its partners' average ionic radii. Using the thermodynamics-grounded bond-length distribution descriptor, we identify Ca^{2+}-Cu^{2+} interactions as the primary prospective lattice destabilizer. Bulk synthesis confirms only 4% Ca solubility with Cu at 950^oC, modestly rising to 5% after Cu removal at 1150^oC. We then employ far-from-equilibrium pulsed-laser deposition to investigate metastable solubility: epitaxial films incorporate 10% Ca with Cu and a full 20% Ca without, doubling and quadrupling the respective bulk limits. This Ca uptake additionally enables deterministic lattice-parameter control via Ca concentration. Overall, our results demonstrate both the extended solubility possible in HEO systems, particularly when accessing metastable states through quenching from high-energy plasma, and that the specific constellation of solid-solvent cations can be rationally engineered to minimize bond-length distributions when largely misfit cations are added, thus expanding the accessible compositional space.
Platinum diselenide (PtSe2) is an emerging two-dimensional (2D) transition metal dichalcogenide known for its excellent electrical and optical properties, along with remarkable air stability. For PtSe2-based electronic devices, understanding high-field breakdown and heat dissipation is crucial for designing high-performance and energy-efficient systems operating under extreme conditions. In this work, we investigate the breakdown mechanisms of semimetallic PtSe2 at both low and room temperatures. Heat dissipation is quantified via interfacial thermal conductivity (ITC) of PtSe2/SiO2 and PtSe2/h-BN interfaces using Raman thermometry. Our findings indicate that at room temperature, device breakdown is predominantly governed by self-heating effects. Conversely, at low temperatures, the breakdown is mainly driven by carrier multiplication under high electric fields, as further confirmed by Hall measurements.
We investigate the emergence of cohesive flocking in open, boundless space using a multi-agent reinforcement learning framework. Agents integrate positional and orientational information from their closest topological neighbours and learn to balance alignment and attractive interactions by optimizing a local cost function that penalizes both excessive separation and close-range crowding. The resulting Vicsek-like dynamics is robust to algorithmic implementation details and yields cohesive collective motion with high polar order. The optimal policy is dominated by strong aligning interactions when agents are sufficiently close to their neighbours, and a flexible combination of alignment and attraction at larger separations. We further characterize the internal structure and dynamics of the resulting groups using liquid-state metrics and neighbour exchange rates, finding qualitative agreement with empirical observations in starling flocks. These results suggest that flocking may emerge in groups of moving agents as an adaptive response to the biological imperatives of staying together while avoiding collisions.
The Langevin equation accounts for unresolved bath degrees of freedom driving the system toward the bath temperature. Because of this, numerical solutions of the Langevin equation have a long history. Here, we recapitulate, combine, and extend existing Langevin-equation based thermostats, scrutinize their properties and demonstrate their superiority over global kinetic-energy controls. Our work includes compact, asymptotic-analysis based derivations of stochastic thermostats, including the highly accurate Grønbech-Jensen scheme. Proposed extensions include a precise, colored and a lean, momentum-conserving thermostat.
We use the quantum metric to understand the properties of quasicrystals, represented by the one-dimensional (1D) Fibonacci chain. We show that the quantum metric can relate the localization properties of the eigenstates to the scale-invariance of both the chain and its energy spectrum. In particular, the quantum metric incorporates information about distances between the local symmetry centers of each eigenstate, making it much more sensitive to the localization properties of quasicrystals than other measures of localization, such as the inverse participation ratio. We further find that a full description of localization requires us to introduce a new phasonic component to the quantum metric, and a mixed phason-position Chern number. Finally, we show that the sum of both position and phasonic components of the quantum metric is lower-bounded by the gap labels associated with each energy gap of the Fibonacci chain. This establishes a direct link between the spatial localization and fractal energy spectrum of quasicrystals. Taken together, the quantum metric provides a unifying, yet accessible, understanding of quasicrystals, rooted in their scale-invariance and with intriguing consequences also for many-body physics.
In this work, we have performed a comprehensive study of dielectric materials NbOX2 (X=Cl, Br, and I) within the framework of density functional theory, incorporating both conventional and hybrid functionals.Our studies focus on the structure, electronic, elastic, and piezoelectric properties. Piezoelectricity is an innovative avenue to extract energy by manipulating the material's structures via mechanical stress. The use of non-lead-based material for piezoelectricity added an advantage of a greener approach. Among the investigated materials, bulk NbOI2 exhibits the highest piezoelectric response of 6.32 C/m^2, which is around 31% higher than NbOCl2, NbOBr2 and even lead zirconate titanate.
The mysterious quantum-critical Planckian bad metal phase with perfect T-linear resistivity persisting beyond the quasi-particle limit and universal T-linear scattering rate has been observed in various high-Tc cuprate superconductors. Here, we develop a realistic theoretical approach to this phase in an analytically solvable large-N multi-channel Kondo lattice model, derived from a heavy-fermion formulated conventionaL t-J model, known for qualitatively describing cuprates. This phase is originated from critical charge Kondo fluctuations where disordered local bosonic charge fluctuations couple to spinon and heavy conduction-electron Fermi surfaces near a charge-Kondo-breakdown local quantum critical point associated with pseudogap-to-Fermi liquid transition. Our results show excellent agreement with experiments and offer broad implications for other unconventional superconductors.
We report on the interplay between atomic scale inhomogeneity and competing magnetic phases and its effect on the anomalous Hall effect in the layered antiferromagnet MnBi_4Te_7, a natural superlattice hosting coexisting ferromagnetic and antiferromagnetic phases. Using a combination of scanning tunneling microscopy (STM), DC and AC magnetization, and magneto-transport measurements, we reveal that intrinsic Mn Bi antisite defects induce strong interlayer exchange coupling, giving rise to a robust exchange bias observed in both magnetic and Hall responses. The exchange bias undergoes a transition from asymmetric to symmetric behavior between 2 K and 6 K, indicating a temperature driven dynamical reconfiguration of interfacial spin structures. The training effect analysis revealed a stronger contribution of frozen spins at 2 K compared to 6 K, with relaxation amplitude shift from -264 Oe to 306 Oe. This sign reversal indicates a field-induced change in interfacial coupling. The temperature dependence of longitudinal resistivity and magnetization reveals complementary behavior, indicating the coexistence of two distinct spin states near the magnetic transition temperature. The phase fraction based resistivity model captures the distinct scattering mechanisms that govern electronic transport across different magnetic regimes. Our findings offer a direct link between microscopic disorder, interfacial magnetism and macroscopic topological phenomena in magnetic topological insulators.
Motivated by recent experiments, we investigate the spin-1/2 XXZ model on the triangular lattice with strong Ising anisotropy, combining large-scale numerical simulations and analytical methods to uncover unconventional spin dynamics. First, we compute the dynamical spin structure factor using density matrix renormalization group (DMRG) simulations and find excellent agreement with inelastic neutron scattering data on the layered compound \text{K}_2\text{Co}(\text{SeO}_3)_2. The low-energy spectrum reveals a roton-like minimum at the M point, absent in linear spin-wave theory, accompanied by peak intensity and a broad continuum above it. Near the \Gamma point, we observe an approximately linear dispersion with vanishing spectral weight. Second, we compare two analytical frameworks that reproduce the observed features. The first is a hard-core boson approach, which includes: (i) an effective staggered boson model (ESBM) at zero magnetic field, (ii) perturbation theory applied to the one-third magnetization plateau, and (iii) a self-consistent mean-field Schwinger boson theory (SBT). The second framework is based on a variational supersolid quantum dimer model (QDM) ansatz, combined with a single-mode approximation. The SBT captures the broad continuum, the M-point minimum, and linear dispersion at \Gamma, whereas the QDM reproduces the roton minimum and linear dispersion at finite momentum near \Gamma. Remarkably, both the QDM wavefunction and the DMRG ground state exhibit nearly identical structure factors with pronounced transverse photon-like excitations. Together, our comprehensive theoretical and numerical analysis elucidates the microscopic origin of supersolid excitations in the XXZ triangular lattice model and their proximity to a spin liquid phase observed experimentally.
The nature of stacking faults - whether intrinsic or extrinsic - plays a pivotal role in defect-mediated processes in crystalline materials. Yet, current electron microscopy techniques for their reliable analysis remain limited to either conventional fringe-contrast imaging of inclined faults or atomic-resolution imaging of edge-on configurations. Here, we overcome this dichotomy by introducing a high-resolution scanning transmission electron microscopy (HRSTEM) method that enables full structural discrimination of inclined stacking faults in fcc and L1_2 crystals. This approach eliminates a long-standing geometric constraint on high-resolution analysis, providing comprehensive access to stacking faults on all glide planes along the widely used [001] and [110] zone axes. We demonstrate the robustness of the method in a CoNi-based superalloy, achieving clear discrimination of fault types even in overlapping configurations and foil thicknesses exceeding 100 nm. Simulations reveal that fault-induced de-channeling is key to contrast formation and is strongly governed by the fault's depth within the sample. Leveraging this effect, we further establish a route to artificially generate ultrathin TEM lamellae - bounded by the stacking fault itself - thereby enhancing contrast for atomic-scale studies of long-range ordering, compositional fluctuations, and nanoclustering.
We propose a novel type of superconducting spin qubit, dubbed solitonic Andreev spin qubit (SASQ), that combines features of Andreev spin qubits and geometric spin qubits. The two SASQ states are the degenerate spin orientations of an Andreev bound state trapped in a circular Josephson junction with a Corbino disk geometry, created on a 2DEG. The junction is subjected to a weak magnetic flux that induces a fluxoid mismatch between the inner disk and outer ring superconductors. The fluxoid mismatch produces a cancellation of the induced pairing that traps unconventional spin-degenerate Andreev bound states analogous to Jackiw-Rebbi solitons. They are localized at a position around the junction that can be controlled by changing the superconducting phase difference. Moving a soliton with a phase bias induces a holonomic rotation of its spin, by virtue of the spin-orbit coupling in the 2DEG. The holonomic qubit trajectories can densely cover the full Bloch sphere as the soliton revolves around the junction. Effects of non-holonomic (dynamic) qubit drift are also analyzed.
In spin-orbit-coupled systems the charge and spin transport are generally coupled to each other, namely a charge current will induce a spin current and vice versa. In the presence of time-reversal symmetry \mathcal{T}, the cross-coupling transport coefficients describing how one process affects the other are constrained by the famous Onsager reciprocal relations. In this paper, we generalize the Onsager reciprocal relations of charge and spin transport to systems that break the time-reversal symmetry but preserve a combined symmetry of \mathcal{T} and some other symmetry operation \mathcal{O}. We show that the symmetry or antisymmetry of the cross-coupling transport coefficients remains in place provided that the operator \mathcal{O} meets certain conditions. Among many candidate systems where our generalized Onsager relations apply, we focus on a conceptually simple and experimentally realized model in cold atomic systems for explicit demonstration and use these relations to predict highly non-trivial transport phenomena that can be readily verified experimentally.
The equilibration dynamics of ultrastable glasses subjected to heating protocols has attracted recent experimental and theoretical interest. With simulations of the mW water model, we investigate the devitrification and melting dynamics of both conventional quenched (QG) and vapor deposited (DG) amorphous ices under controlled heating ramps. By developing an algorithm to reconstruct hydrogen-bond networks, we show that bond ring statistics correlates with the structural stability of the glasses and allows tracking crystalline and liquid clusters during devitrification and melting. We find that QG melts in the bulk, whereas melting in DG preferentially begins near the free surface. During devitrification, the DG shows an excess of 5-membered rings near the free surface, which is consistent with its tendency to nucleate the crystal phase in this region. Additionally, the DG shows an Avrami exponent exceeding the standard 1+d behavior, while both glasses display the same sub-3d growth of liquid clusters across heating rates, indicating that the DG enhanced exponent stems from its higher kinetic stability.
Advances in photovoltaic technology are a viable route to contribute to cleaner and more sustainable energy solutions, placing perovskite-based materials among the best candidates for solar energy conversion. However, some challenges must be addressed to enhance their performance and stability. Herein, we report an investigation of the AsNCa3 antiperovskite system for its potential in photovoltaic devices. We consider eight distinct crystalline phases, their structural parameters, dynamical stability, and electronic and optical properties. Furthermore, we consider each structural phase's contributions to solar harvesting efficiency by calculating the power conversion efficiency (PCE) using the spectroscopiclimited maximum efficiency (SLME) formalism, which in this case reaches a maximum of 31.2%. All dynamically stable phases exhibit a band gap around 1.3 eV, which lies within the optimal range for single-junction solar cells and yields PCE values comparable to the theoretical maximum PCE for silicon. These results place AsNCa3 antiperovskites as promising candidates for high-efficiency photovoltaic applications. Notably, the PCE is only slightly changed by structural phase modification, suggesting that phase transitions induced by environmental conditions during device operation might not compromise device performance.
In this paper we introduce and study the Maximum-Average Subtensor (p-MAS) problem, in which one wants to find a subtensor of size k of a given random tensor of size N, both of order p, with maximum sum of entries. We are motivated by recent work on the matrix case of the problem in which several equilibrium and non-equilibrium properties have been characterized analytically in the asymptotic regime 1 \ll k \ll N, and a puzzling phenomenon was observed involving the coexistence of a clustered equilibrium phase and an efficient algorithm which produces submatrices in this phase. Here we extend previous results on equilibrium and algorithmic properties for the matrix case to the tensor case. We show that the tensor case has a similar equilibrium phase diagram as the matrix case, and an overall similar phenomenology for the considered algorithms. Additionally, we consider out-of-equilibrium landscape properties using Overlap Gap Properties and Franz-Parisi analysis, and discuss the implications or lack-thereof for average-case algorithmic hardness.
Successive phase transitions in the quasi-kagome compound URhSn were investigated by resonant X-ray scattering (RXS) at the uranium M_4 edge. In the high-temperature phase between 16 K and 54 K, an additional RXS signal was detected superposed onto fundamental reflections in both \pi-\sigma' and \pi-\pi' polarization channels. Upon cooling below 16 K, reported as a ferromagnetic phase along c, substantial enhancements were observed again in the both polarization channels at the 300 reflection, demonstrating a simultaneous emergence of in-plane spin alongside the c-axis ferromagnetic components. The observed behavior can be interpreted by an antiferro-quadrupole (AFQ) order of O_{yz} or O_{zx} characterized by a propagation vector q = 0 in the intermediate phase, which then coexists with a ferromagnetic component below 16 K. The resulting ground state structure breaks the mirror symmetry perpendicular to the kagome plane, identifying the formation of a unique AFQ order with either chirality or polarity in URhSn.