2025-01-13 | | Total: 67
We address the longstanding challenge in quantum statistical mechanics of reconciling unitary dynamics with irreversible relaxation. In classical chaos, the unitary evolution operator develops Ruelle-Pollicott (RP) resonances inside the unit circle in the continuum limit, leading to mixing. In contrast, the quantum theory of many-body RP resonances and their link to irreversibility remain underdeveloped. We relate the spectral form factor to the sum of autocorrelation functions and, in generic many-body systems without conservation laws, argue that all quantum RP resonances converge inside the unit disk. While we conjecture this picture to be general, we analytically prove the emergence of irreversibility in the random phase model (RPM), a paradigmatic Floquet quantum circuit model, in the limit of large local Hilbert space dimension. To this end, we couple it to local environments and compute the exact time evolution of autocorrelation functions, the dissipative form factor, and out-of-time-order correlation functions. Although valid for any dissipation strength, we then focus on weak dissipation to clarify the origin of irreversibility in unitary systems. When the dissipationless limit is taken after the thermodynamic limit, the unitary quantum map develops an infinite tower of RP resonances -- chaotic systems display so-called anomalous relaxation. We identify the exact RP resonances in the RPM and prove that the same RP resonances are obtained from operator truncation. We also show that the OTOC in the RPM can undergo a two-stage relaxation and that during the second stage, the approach to the stationary value is again controlled by the leading RP resonance. Finally, we demonstrate how conservation laws, many-body localization, and nonlocal interactions merge the leading RP resonance into the unit circle, thereby suppressing anomalous relaxation.
We review the recent development of machine-learning (ML) force-field frameworks for Landau-Lifshitz-Gilbert (LLG) dynamics simulations of itinerant electron magnets, focusing on the general theory and implementations of symmetry-invariant representations of spin configurations. The crucial properties that such magnetic descriptors must satisfy are differentiability with respect to spin rotations and invariance to both lattice point-group symmetry and internal spin rotation symmetry. We propose an efficient implementation based on the concept of reference irreducible representations, modified from the group-theoretical power-spectrum and bispectrum methods. The ML framework is demonstrated using the s-d models, which are widely applied in spintronics research. We show that LLG simulations based on local fields predicted by the trained ML models successfully reproduce representative non-collinear spin structures, including 120$^\circ$, tetrahedral, and skyrmion crystal orders of the triangular-lattice s-d models. Large-scale thermal quench simulations enabled by ML models further reveal intriguing freezing dynamics and glassy stripe states consisting of skyrmions and bi-merons. Our work highlights the utility of ML force-field approach to dynamical modeling of complex spin orders in itinerant electron magnets.
This paper provides an overview of the research on the metastable behavior of the Ising model. We analyze the transition times from the set of metastable states to the set of the stable states by identifying the critical configurations that the system crosses with high probability during this transition and by computing the energy barrier that the system must overcome to reach the stable state starting from the metastable one. We describe the dynamical phase transition of the Ising model evolving under Glauber dynamics across various contexts, including different lattices, dimensions and anisotropic variants. The analysis is extended to related models, such as long-range Ising model, Blume-Capel and Potts models, as well as to dynamics like Kawasaki dynamics, providing insights into metastability across different systems.
Graphene spintronics offers a promising route to achieve low power 2D electronics for next generation classical and quantum computation. As device length scales are reduced to the limit of the electron mean free path, the transport mechanism crosses over to the ballistic regime. However, ballistic transport has yet to be shown in a graphene spintronic device, a necessary step towards realising ballistic spintronics. Here, we report ballistic injection of spin polarised carriers via one-dimensional contacts between magnetic nanowires and a high mobility graphene channel. The nanowire-graphene interface defines an effective constriction that confines charge carriers over a length scale smaller than that of their mean free path. This is evidenced by the observation of quantised conductance through the contacts with no applied magnetic field and a transition into the quantum Hall regime with increasing field strength. These effects occur in the absence of any constriction in the graphene itself and occur across several devices with transmission probability in the range T = 0.08 - 0.30.
We are interested in rationalizing the phenomenon of like-charge attraction between charged bodies, such as a pair of colloids, in the strong coupling regime. The two colloids are modelled as uniformly charged parallel plates, neutralized by mobile counterions. In an earlier work [Palaia et al., J. Phys. Chem. B 126, 3143 (2022)], we developed an effective-field theory for symmetric plates, stemming from the ground-state description that holds at infinite couplings. Here, we generalize the approach to the asymmetric case, where the plates bear charges of the same sign, but of different values. In the symmetric situation, the mobile ions, which are localized in the vicinity of the two plates, share equally between both of them. Here, the sharing is non-trivial, depending both on the coupling parameter and the distance between the plates. We thus introduce a counterion occupation parameter, that is determined variationally to ensure minimum of the free energy. The resulting analytical results for the pressure as a function of the plate-plate distance $d$ agree well with our Monte Carlo data, in a large interval of strong and intermediate coupling constants $\Xi$. We show in particular that within this description, there exists a range of large distances at which the attractive pressure features a $1/d^2$ behavior.
We study the stability of vortices in a binary system of Bose-Einstein condensates, with their wave functions modeled by a set of coupled, time-dependent Gross-Pitaevskii equations. Beginning with an effective two-dimensional system, we identify miscible and immiscible regimes characterized by the inter- and intra-atomic interactions and the initial configuration of the system. We then consider a binary system of Bose-Einstein condensates placed in a rotating harmonic trap and study the single vortex state in this system. We derive an approximate form for the energy of a single vortex in the binary system and the critical angular velocity for the global stability of a vortex at the center of the trap. We also compute the metastability onset angular velocity for the local stability of a vortex at the center of the trap. Numerical solutions to the Gross-Pitaevskii equations support these expressions. These rotational results inform us of a novel subphase within the miscible regime of the binary condensate system. We thus demonstrate the non-trivial aspects of vortex stability in interacting binary Bose-Einstein condensates as a result of their non-linear interactions.
Titanium and its alloys are technologically important materials that display a rich phase behaviour. In order to enable large-scale, realistic modelling of Ti and its alloys on the atomistic scale, Machine Learning Interatomic Potentials (MLIPs) are crucial, but rely on databases of atomic configurations. We report databases of such configurations that represent the {\alpha}, \b{eta}, {\omega} and liquid phases of Ti and the Ti-6Al-4V alloy, where we provide total energy, force and stress values evaluated by Density Functional Theory (DFT) using the PBE exchange-correlation functional. We have also leveraged and extended a data reduction strategy, via non-diagonal supercells, for the vibrational properties of Ti and sampling of atomic species within bulk crystalline data for Ti-6Al-4V. These configurations may be used to fit MLIP models that can accurately model the phase behaviour of Ti and Ti-6Al-4V across a broad range of thermodynamic conditions. To validate models, we assembled a set of benchmark protocols, which can be used to rapidly develop and evaluate MLIP models. We demonstrated the utility of our databases and validation tools by fitting models based on the Gaussian Approximation Potential (GAP) and Atomic Cluster Expansion (ACE) frameworks.
Maximum-entropy methods, rooted in the inverse Ising/Potts problem from statistical mechanics, have become indispensable tools for modeling pairwise interactions in disciplines such as bioinformatics, ecology, and neuroscience. Despite their remarkable success, these methods often overlook high-order interactions that may be crucial in complex systems. Conversely, while modern machine learning approaches can capture such interactions, existing interpretable frameworks are computationally expensive, making it impractical to assess the relevance of high-order interactions in real-world scenarios. Restricted Boltzmann Machines (RBMs) offer a computationally efficient alternative by encoding statistical correlations via hidden nodes in a bipartite neural network. Here, we present a method that maps RBMs exactly onto generalized Potts models with interactions of arbitrary high order. This approach leverages large-$N$ approximations, facilitated by the simple architecture of the RBM, to enable the efficient extraction of effective many-body couplings with minimal computational cost. This mapping also enables the development of a general formal framework for the extraction of effective higher-order interactions in arbitrarily complex probabilistic models. Additionally, we introduce a robust formalism for gauge fixing within the generalized Potts model. We validate our method by accurately recovering two- and three-body interactions from synthetic datasets. Additionally, applying our framework to protein sequence data demonstrates its effectiveness in reconstructing protein contact maps, achieving performance comparable to state-of-the-art inverse Potts models. These results position RBMs as a powerful and efficient tool for investigating high-order interactions in complex systems.
Active fluids exhibit chaotic flows at low Reynolds number known as active turbulence. Whereas the statistical properties of the chaotic flows are increasingly well understood, the nature of the transition from laminar to turbulent flows as activity increases remains unclear. Here, through simulations of a minimal model of unbounded active nematics, we find that the transition to active turbulence is discontinuous. We show that the transition features a jump in the mean-squared velocity, as well as bistability and hysteresis between laminar and chaotic flows. From distributions of finite-time Lyapunov exponents, we identify the transition at a value $A^*\approx 4900$ of the dimensionless activity number. Below the transition to chaos, we find subcritical bifurcations that feature bistability of different laminar patterns. These bifurcations give rise to oscillations and to chaotic transients, which become very long close to the transition to turbulence. Overall, our findings contrast with the continuous transition to turbulence in channel confinement, where turbulent puffs emerge within a laminar background. We propose that, without confinement, the long-range hydrodynamic interactions of Stokes flow suppress the spatial coexistence of different flow states, and thus render the transition discontinuous.
Recent numerical works have revealed the instability of many-body localized (MBL) phase in disordered quantum many-body systems with finite system sizes and over finite timescales. This instability arises from Griffith regions that occur at the thermodynamic limit, which rapidly thermalize and affect the surrounding typical MBL regions, introducing an avalanche mechanism into the system. Here, we consider the $\mathbb{Z}_2$-preserving interacting Ising Majorana chain model, which exhibits a more complex phase diagram, where an ergodic phase emerges between two MBL phases with different long-range order properties. We calculate the dynamic characteristics of the model when coupled to an infinite bath under perturbation, and through scaling behavior of the slowest thermalization rate, we find how critical disorder strengths in finite-size systems are affected by the avalanche mechanism. We also employe the embedded inclusion model and use the time evolution of mutual information between each spin and the artificial Griffith region to probe the diffusion of the thermal bubble. We observe that in finite-sized systems, the critical disorder strength gradually drifts away from the central. Our work demonstrate that both MBL paramagnetic phase and MBL spin-glass phase are unstable at finite sizes.
We report the bias voltage controlled inversions of tunneling magnetoresistance (TMR) in magnetic tunnel junctions composed of Fe3GaTe2 electrodes and hBN tunneling barrier, observed at room temperature. The polarity reversal of TMR occurs consistently at around 0.625 V across multiple devices and temperatures, highlighting the robustness of the effect. To understand this behavior, we developed a theoretical model incorporating spin-resolved density of states (DOS) at high energy levels. By adjusting the DOS weighting at different k points to account for misalignment between the crystal structure of electrodes in experimental devices, we improved agreement between experimental and theoretical inversion voltages. Our results provide valuable insight into the voltage-controlled spin injection and detection in two-dimensional magnetic tunnel junctions, with implications for the development of energy-efficient spintronic devices.
We study non-Abelian fractional quantum Hall state in double layer systems at total filling factor $1/2$. Recent progresses in two-dimensional van der Waals materials made it possible to explore the regime with very small interlayer distance. Numerical calculations suggests interlayer phase coherence can develop between the layers such that the electrons may redistribute between them without changing the Hall response. It corresponds to spontaneous breaking of the U(1) symmetry associated with the particle number difference in the layers. This state manifests itself as superfluid in counterflow measurement and has characteristic Hall response when current is passed through one layer and voltages in both layers are measured. As the interlayer distance increases, a phase transition into the Halperin 331 state occurs. We also discuss similar physics for bosonic systems with specially designed interactions.
We conduct a systematic analysis of cavity effects on the decay dynamics of an open magnonic system. The Purcell effect on the magnon oscillator decay is thoroughly examined for both driven and non-driven scenarios. Analytical conditions are determined to distinguish between strong and weak coupling regimes, corresponding to oscillatory and pure decay behaviors respectively. Additionally, our theory also predicts the decay of the photon mode within the cavity-magnonic open system, demonstrating excellent agreement with existing experimental data. Our findings and methodologies can provide valuable insights for advancing research in cavity magnonic quantum control, quantum information processing, and the development of magnonic quantum devices.
Here we investigate epitaxial Hf0.5Zr0.5O2 ferroelectric thin films as potential candidates to be used as non-volatile electric-field-modulated thermal memories. The electric-field dependence of the thermal conductivity of metal/Hf0.5Zr0.5O2/YSZ devices is found to be hysteretic, resembling the polarization vs electric field hysteresis loops, being maximum (minimum) at large applied positive (negative) voltages from the top metallic electrode. This dynamic thermal response is compatible with the coupling between the ferroelectric polarization and the oxygen ion migration, in which the oxygen vacancies are the main phonon scattering sources and the polarization acts as an electrically active ion migration barrier that creates the hysteresis. This new mechanism enables two non-volatile thermal states: high (ON) and low (OFF) thermal conductivity, with an ON/OFF ratio of 1.6. Both the ON and OFF states exhibit high stability over time, though the switching speed is limited by ion mobility in the YSZ electrode.
Transition metal dichalcogenides are a family of quasi-two-dimensional materials that display a high technological potential due to their wide range of electronic ground states, e.g., from superconducting to semiconducting, depending on the chemical composition, crystal structure, or electrostatic doping. Here, we unveil that by tuning a single parameter, the hydrostatic pressure P, a cascade of electronic phase transitions can be induced in the few-layer transition metal dichalcogenide 1T'-WS2, including superconducting, topological, and anomalous Hall effect phases. Specifically, as P increases, we observe a dual phase transition: the suppression of superconductivity with the concomitant emergence of an anomalous Hall effect at P=1.15 GPa. Remarkably, upon further increasing the pressure above 1.6 GPa, we uncover a reentrant superconducting state that emerges out of a state still exhibiting an anomalous Hall effect. This superconducting state shows a marked increase in superconducting anisotropy with respect to the phase observed at ambient pressure, suggesting a different superconducting state with a distinct pairing symmetry. Via first-principles calculations, we demonstrate that the system concomitantly transitions into a strong topological phase with markedly different band orbital characters and Fermi surfaces contributing to the superconductivity. These findings position 1T'-WS2 as a unique, tunable superconductor, wherein superconductivity, anomalous transport, and band features can be tuned through the application of moderate pressures.
We develop a continuum theory to describe the collective dynamics of deformable epithelial cells, using two tensor order parameters to distinguish the force-generating active filaments in the cells from their shape. The theory demonstrates how active flows create nematic domains of anisotropic cells, drive active turbulence, and create topological defects. We show that the filament flow-aligning parameter, $\lambda_{\tens{Q}}$, a rheological quantity that determines the response of the filaments to velocity gradients in the active flows, plays a significant, to date unappreciated, role in determining the pattern of extensional and compressional active flows. In a contractile cell layer, local flows are expected to align elongated cells perpendicular to the active filaments. However, with increasing $\lambda_{\tens{Q}}$, long-range correlations in the active turbulent flow field lead to extended regions where this alignment is parallel, as recently observed in experiments on confluent MDCK cell layers. Further, the two order parameter formalism allows us to distinguish defects in the filament director field, which contribute to the active driving, and those in the shape director field, measured in experiments, which are advected by the active flows. By considering the relative orientations of shape and filaments we are able to explain the surprising observation of defects moving towards their heads in contractile cell layers.
The quantum dynamics of a compound sample consisting from a semiconductor double quantum dot (DQD) system non-linearly coupled with a leaking single-mode micro-resonator is theoretically investigated. The focus is on the resonance condition when the transition frequency of the double quantum dot equals to the doubled resonator frequency, respectively, and the resulting interplay among the involved phonon or photon channels. As a result, the steady-state quantum dynamics of this complex non-linear system exhibits a variety of possible effects that have been demonstrated here. Particularly, we have found the relationship among the electrical current through the double quantum dot and the microwave field inside the resonator that is nonlinearly coupled to it, with a corresponding emphasizing on their critical behaviors. Additionally, the quantum correlations of the photon flux generated into the resonator mode vary from super-Poissonian to Poissonian photon statistics, leading to single-qubit lasing phenomena at microwave frequencies.
Glucose oxidation is a key reaction in biosensing and energy conversion, with anodized copper electrodes showing great potential as catalysts for its enhancement. This study explores the relationship between anodization parameters and glucose oxidation on thin copper films, providing valuable insights for optimization. By adjusting anodization parameters, specific species such as Cu2O, Cu(OH)2, and CuO can be selectively formed, each exhibiting different catalytic activities. Applying polarization in 1 M KOH at 0 V (vs Ag|AgCl) generates a highly active CuO surface layer, leading to significant performance improvements over bare copper electrodes. Specifically, sensitivity increases by 55% in the 0.1 mM-0.5 mM range and 73% in the 0.75 mM-2 mM range. Additionally, the electrodes demonstrate an exceptionally low detection limit of just 0.004 mM. These results are attributed to scalable, reproducible manufacturing on cleanroom-compatible platforms. This study not only sheds light on the effects of thin copper film anodization on glucose oxidation but also provides a practical approach for enhancing the efficiency of copper-based systems in biosensing and energy conversion applications.
The Kekulé-Kitaev model, an extension of the Kitaev model, exhibits quantum spin liquid (QSL) properties, which has an exact solution through Kitaev parton construction. In this study, we calculate the dynamical spin structure factor as a linear response and the third-order magnetic susceptibility as a nonlinear response using two-dimensional coherent spectroscopy for Kekulé-Kitaev model. Our results reveal that the few-matter fermion excitations approximation provides reliable results for both linear and nonlinear responses. Notably, while the Kekulé-Kitaev model shows linear and nonlinear responses qualitatively similar to the Kitaev model, it displays distinct behavior under a weak uniform/staggered magnetic field. Specifically, the Kekulé-Kitaev model does not present a non-Abelian phase under a uniform magnetic field, while such phase appears in the presence of a staggered magnetic field. Interestingly, within the non-Abelian phase, signals originating from two non-adjacent fluxes in the nonlinear response are stronger than signals from other flux excitations. Furthermore, We demonstrate that the ground state of the Kekulé-Kitaev model at isotropic coupling is mapped, through unitary spin rotations, to an excited state of the Kitaev model with a uniform flux configuration.
Recent advancements in exchange-correlation functionals within density functional theory highlight the need for rigorous validation across diverse types of materials properties. In this study, we assess the performance of the newly developed meta-GGA r$^2$SCAN and its predecessor, SCAN, in predicting the Néel transition temperature of antiferromagnetic materials. Our analysis includes 48 magnetic materials, spanning both simple and complex systems. Using DFT, we compute the energies of various magnetic configurations and extract exchange interaction parameters through a least-squares fitting approach. These parameters are then used in classical Monte Carlo simulations to estimate the transition temperatures. Our results demonstrate that both SCAN and r$^2$SCAN greatly outperform standard GGA and GGA+$U$ methods, yielding predictions that closely align with experimental values. The Pearson correlation coefficients for SCAN and r$^2$SCAN are 0.97 and 0.98, respectively, when compared to experimental transition temperatures. Additionally, we calculate the energy differences between antiferromagnetic and ferromagnetic configurations to assess the performance of the hybrid HSE06 functional. We found that the HSE06 functional underestimates transition temperatures compared to the meta-GGA functionals and experimental values.
Single-time and two-time correlators are computed exactly in the $1D$ Glauber-Ising model after a quench to zero temperature and on a periodic chain of finite length $N$, using a simple analytical continuation technique. Besides the general confirmation of finite-size scaling in non-equilibrium dynamics, this allows to test the scaling behaviour of the plateau height $C_{\infty}^{(2)}$ to which the two-time auto-correlator converges, when deep into the finite-size regime.
Materials and systems that exhibit persistent spin texture provide a platform for creating robust spin states that can be used in quantum computing, memory storage, and other advanced technological applications. In this paper we show that persistent spin-texture in single-layer phosphorene electrons close to the $\Gamma$ point, subjected to the finite perpendicular electric field, can be achieved by appropriately tuning the extrinsic spin-orbit coupling strength using the tensile strain of about $1.2\%$ in the zigzag direction. This is confirmed by detailed numerical investigations of the effects of strain on the intrinsic and extrinsic spin-orbit coupling, and by the effective spin-orbit Hamiltonian of phosphorene electrons and holes around the $\Gamma$ point, assuming the presence of the perpendicular electric field. Furthermore, the calculated spin relaxation rates due to the Dyakonov-Perel mechanism indicate a giant anisotropy of the in-plane spin, up to $10^5$, which is directly related to the discovered persistent spin texture of phosphorene electrons close to the $\Gamma$ point. %%%%%%%% We also show, that strain can reverse the anisotropy of spin mixing parameter $b^2$ connected to the Elliott-Yafet spin relaxation mechanism which dominates spin relaxation in phosphorene. We find the conditions under which Elliott-Yafet spin lifetime anisotropy can be largely enhanced due to synergy of spin mixing and g-factor anisotropy. %%%%%%%%%% Our results suggest that spin texture in phosphorene can be modulated by strain, enabling its potential usage in the field of spintronics.
We revisit the structure of the phase diagram of the two-component mean-field Bose mixture at finite temperatures, considering both the cases of attractive and repulsive interspecies interactions. In particular, we analyze the evolution of the phase diagram upon driving the system towards collapse and point out its distinctive features in this limit. We provide analytical insights into the global structure of the phase diagram and the properties of the phase transitions between the normal phase and the phases involving Bose-Einstein condensates. \emph{Inter alia} we analytically demonstrate that for sufficiently weak interspecies interactions $a_{12}$ the system generically exhibits a line of quadruple points but has no triple nor tricritical points in the phase diagram spanned by the chemical potentials $\mu_1$, $\mu_2$ and temperature $T$. In contrast, for sufficiently large, positive values of $a_{12}$, the system displays both triple and tricritical points but no quadruple points. As pointed out in recent studies, in addition to the phase transitions involving condensation, the mixture may be driven through a liquid-gas type transition, and we clarify the conditions for its occurrence. We finally discuss the impact of interaction- and mass-imbalance on the phase diagram of the mixture.
Does a machine learning model actually gain an understanding of the material space? We answer this question in the affirmative on the example of the OptiMate model, a graph attention network trained to predict the optical properties of semiconductors and insulators. By applying the UMAP dimensionality reduction technique to its latent embeddings, we demonstrate that the model captures a nuanced and interpretable representation of the materials space, reflecting chemical and physical principles, without any user-induced bias. This enables clustering of almost 10,000 materials based on optical properties and chemical similarities. Beyond this understanding, we demonstrate how the learned material space can be used to identify more sustainable alternatives to critical materials in energy-related technologies, such as photovoltaics. These findings demonstrate the dual utility of machine learning models in materials science: Accurately predicting material properties while providing insights into the underlying materials space. The approach demonstrates the broader potential of leveraging learned materials spaces for the discovery and design of materials for diverse applications, and is easily applicable to any state-of-the-art machine learning model.
To realize practical quantum computers, a large number of quantum bits (qubits) will be required. Semiconductor spin qubits offer advantages such as high scalability and compatibility with existing semiconductor technologies. However, as the number of qubits increases, manual qubit tuning becomes infeasible, motivating automated tuning approaches. In this study, we use U-Net, a neural network method for object detection, to identify charge transition lines in experimental charge stability diagrams. The extracted charge transition lines are analyzed using the Hough transform to determine their positions and angles. Based on this analysis, we obtain the transformation matrix to virtual gates. Furthermore, we identify the single-electron regime by clustering the Hough transform outputs. We also show the single-electron regime within the virtual gate space. These sequential processes are performed automatically. This approach will advance automated control technologies for large-scale quantum devices.