2026-04-16 | | Total: 100
Non-Hermitian band descriptions capture how loss, gain, and environmental coupling reshape quantum matter, yet most experimental tests rely on wave-based or dynamical probes. Here we establish a new equilibrium route to exceptional physics in Dirac materials: in the weakly non-Hermitian regime, the thermodynamic density of states and the quantum capacitance exhibit a universal equilibrium approach to the exceptional point. In our minimal non-reciprocal graphene model, the hopping imbalance reduces the Dirac velocity as $v_F=v\sqrt{1-β^2}$, implying that the low-energy density of states, the thermodynamic density of states, and the quantum capacitance all scale as $(1-β^2)^{-1}$ as $|β|\to 1^-$. Consequently, at charge neutrality the quantum capacitance remains linear in temperature but with a diverging prefactor, while the inverse response softens linearly on approaching the exceptional point. In a magnetic field, this manifests as a collapse of the Landau-level spacing and a corresponding crowding of thermally active levels. Complementarily, the biorthogonal Bloch states exhibit a Petermann factor $K=(1-β^2)^{-1}$, which isolates the irreducibly non-Hermitian effect of eigenvector non-orthogonality. These results identify quantum capacitance as an experimentally accessible bulk equilibrium probe of effective non-Hermiticity in Dirac materials.
Investigating the interplay among topology and electron-electron interactions is an intriguing research quest which has recently gathered steam across the community of condensed-matter physics. In the present work, we study the fate of a three-dimensional Berry-dipole semimetal, lying at the topological quantum critical point separating a Hopf insulator from a trivial insulator, in the presence of long-range Coulomb interactions. Utilizing large-$N_f$ analysis at three spatial dimensions and an $ε$-expansion within the renormalization-group scheme, we uncover the emergence of a spatially \textit{anisotropic} non-Fermi liquid with enhanced Berry-dipole moment. We further derive the corresponding scaling relations of certain physical observables as functions of the probed energy and temperature scale, and we provide a simple observational criterion for distinguishing the onset of the topological anisotropic non-Fermi liquid from a Berry-dipole semimetal.
We construct upper bounds on entanglement entropies of many-body quantum states that have fixed energy expectation values with respect to geometrically local Hamiltonians. Our focus is on entanglement entropies of subsystems that make up approximately half of the full system. The upper bound on the von Neumann entanglement entropy is half the sum of the thermal entropies of two fictitious systems at the same temperature as one another, with an additional area-law contribution in some systems. The effective temperature is chosen such that the sum of the thermal energies of the two fictitious systems matches the constraint on the energy of the state in the original problem; at subextensive energies, this temperature decreases with increasing system size. Our upper bounds on Rényi entanglement entropies take an analogous form. As a first application we show that ground-state Schmidt ranks in frustration-free (FF) systems are upper bounded by the ground-state degeneracies of Hamiltonians acting on subsystems. Ground-state von Neumann and Rényi entanglement entropies therefore follow an area law when the zero-temperature thermal entropies of subsystems scale with surface areas, rather than with subsystem volumes. This result holds independently of the spectral gap. For physical models of quantum matter, which have well-defined specific heat capacities (and are not necessarily FF), our bounds provide a way to convert this thermodynamic data into constraints on pure-state entanglement at both subextensive and extensive energies. We also show that our upper bounds on half-system entanglement entropies are optimal, up to subleading corrections, in wide varieties of systems. Our results relate physical thermodynamic properties to the structure of many-body Hilbert space at low energies.
Granular flows are ubiquitous in nature and industrial applications, yet a complete continuum theory remains a long-standing challenge. The leading empirical approach, μ(I) rheology, lacks microscopic foundations and becomes multivalued in dense, slowly sheared flows where nonlocal corrections are required. Exploiting state-of-the-art high-speed X-ray tomography to investigate microscopic dynamics of dense granular flows in a Couette geometry, we establish a new, universal constitutive law spanning quasi-static to inertial regimes based on structural relaxation, resolving the fundamental difficulty in the original μ(I) framework. By further establishing a non-equilibrium statistical framework for granular flows, we demonstrate an intrinsic analogy between driven granular matter and hard-sphere liquids owing to their identical Carnahan-Starling equation of state, naturally explaining our rheological approach and the emergence of glassy behaviors. Our framework unifies granular rheology with the broader physics of disordered systems and provides a complete, microscopically-based theoretical framework for dense granular flow.
Materials discovery is fundamental to advance next-generation technologies as well as for sustainable and circular economy. Beyond computational screening, generative models are efficient at finding materials with desired properties, via multi-modal learning using multiscale data. This perspective examines the landscape of generative design for inorganic materials and discusses the integration of multi-modal learning with high-throughput experimental validation. We contextualize these challenges through the lens of a generative design framework as a unified approach to address the data-driven inverse design of functional materials. The central idea of the framework is constructed around a foundation AI model for inorganic materials interlinked deeply with various property databases and high-throughput experiments via a machine learning driven closed loop, which enables the framework to solve key challenges in functional materials. We argue that domain-specific implementations of such integrated workflows represent a promising pathway toward the unresolved challenge of data-driven inverse design for atom-engineered inorganic functional materials.
Long-term chemical durability of glass, crucial for immobilizing nuclear waste, is governed by glass properties such as composition, surface geometry, as well as external factors like thermodynamic conditions and surrounding medium. Despite decades of research, there are no models that account for these intrinsic and extrinsic factors to predict the dissolution rates of glass compositions. To address this challenge, we evaluate the role of natural language embeddings capturing the synthesis and testing conditions in enhancing the predictability of glass dissolution. Evaluating the approach on hand-curated ~700 datapoints extracted from the literature, we reveal that the machine learning (ML) model including natural language embeddings (NLP-ML) outperforms classical ML model in predicting glass dissolution rate. Furthermore, we developed a generalizable ML model by transforming the compositional features to structural descriptors of glass alongside NLP-derived features, enabling extrapolation capability to glass compositions with completely new elements absent in the training data. Evaluating this model on a completely new dataset of glass compositions 34 chemical components in contrast to the training dataset that had only 28 components, we demonstrate that the model indeed exhibits generalizability to glass compositions that are out-of-distribution. Altogether, this integrated approach offers a pathway towards high-fidelity glass dissolution prediction and accelerate the discovery of novel glass compositions with tailored durability for sustainable nuclear waste management.
We investigate the thermal responses of a harmonic oscillator chain coupled at its boundaries to heat baths held at different temperatures. This setup sustains a steady energy flux, continuously dissipating heat into both reservoirs. By introducing slow variations in the bath temperatures, we quantify the resulting excess heat currents and thereby obtain the nonequilibrium heat capacity matrix at fixed but arbitrary temperature differences. We demonstrate the existence of a well-defined thermodynamic limit for long chains. The specific heat associated with energy exchanges with a single bath depends on the difference in friction coefficients governing the system-bath couplings. That thermokinetic effect is typical for nonequilibrium response. When the couplings with the thermal baths acquire temperature dependence, the specific heat correspondingly inherits a nontrivial temperature dependence, in sharp contrast with equilibrium. Our results provide the first explicit determination of specific heat(s) in a locally interacting, spatially extended driven system. Beyond its exact solvability, the model may offer a natural nonequilibrium extension of the Dulong-Petit law, capturing the high-temperature behavior of driven molecules.
The supercurrent diode effect (SDE), characterized by unequal critical currents in opposite directions, has been observed with or without magnetic fields, yet mechanisms enabling zero-field SDE without explicit symmetry breaking remain underexplored. Here we investigate a Josephson junction with strong electron-electron interaction modeled by a Hubbard $U$ term and an odd number of electrons. We find that strong correlations induce spontaneous breaking of time-reversal and mirror symmetries, forming a $\varphi$-junction with degenerate energy minima at $\pm\varphi$, resulting in zero-field Josephson diode effect (JDE) without magnetic order. Spin-orbit coupling breaks SU(2) symmetry but does not determine diode polarity, contrasting with magneto-chiral mechanisms. We further show that applying a tiny Zeeman field enables controllable JDE with sizable efficiency due to the enhancement by the strong magnetic correlation, and the JDE strength peaks when the field induces a level-crossing transition. These findings establish strong electron correlation as a distinct mechanism for nonreciprocal superconducting transport, broadening the understanding of SDE origins.
The interplay between charge and spin dynamics is at the heart of strongly correlated materials. Inspired by recent quantum simulation experiments, we develop a conserving diagrammatic method to describe the Fermi-Hubbard model for strong repulsion and small hole doping away from the half-filled anti-ferromagnetic ground state. We show that doping leads to four hole pockets in the Brillouin zone formed by magnetic polarons, which become increasingly damped with hole concentration. Likewise, the magnon spectrum of the anti-ferromagnet softens and dampens with doping due to hole-induced magnetic frustration. This gives rise to a suppression of the anti-ferromagnetic correlations in agreement with recent experiments. We then calculate the response of the system to a lattice modulation and recover the qualitative difference between in-phase and out-of-phase modulations seen in experiments, which was interpreted as signs of pseudogap physics. Our results indicate that the complex competition between spin and charge degrees of freedom and the emergence of the pseudogap phase may be usefully analyzed for small dopings, where systematic theories can be developed.
The current-driven resistive transition is central to superconducting single-photon detectors, transition-edge sensors, and fluxonic devices. Depending on sample uniformity, dimensions, and heat removal, it can be driven by phase-slip events, flux-flow instabilities (FFI), or normal-domain formation. Here, we investigate the influence of two types of microstructural defects on vortex dynamics in MgB$_2$ films: columnar growth in textured films and buffer-layer roughness in single-crystal films. The current-voltage ($I$-$V$) curves measured at $T \approx 0.25 T_\mathrm{c}$ for both films exhibit multiple steps. Time-dependent Ginzburg-Landau simulations reproduce the major features of the experimental $I$-$V$ curves and suggest that the resistive transitions for both films are mediated by the formation and growth of normal domains rather than FFI. The single-crystal film with buffer-layer roughness exhibits superconductivity breakdown at higher currents and pinning activation energies approximately twice those of the textured film, along with more pronounced multi-step features in the $I$-$V$ curves. These features are attributed to the combination of stronger pinning induced by lateral variations of the superconducting order parameter along the MgO buffer layer and its lower thermal boundary resistance. Our results show that both the film microstructure and the film-buffer interface are critical for the resistive transition, offering insights for superconducting devices requiring controlled dissipation at high transport currents.
The novel kagome family RTi3Bi4 (R: rare-earth) offers a unique platform for exploring distinctive physical phenomena such as anisotropy, spin density wave, and anomalous Hall effect. In particular, the magnetic frustration and behavior of magnetic anisotropy in antiferromagnetic (AFM) kagome materials are of great interest for the fundamental studies and hold promise for next-generation device applications. Here, we report a tunable bifurcation of magnetic anisotropic and bi-oriented AFM order observed in the quasi-1D kagome antiferromagnet GdTi3Bi4. The magnetic domain evolutions during two plateau transition processes are directly visualized, unveiling a pronounced in-plane anisotropy along the a-axis. Temperature-dependent characterization reveals a bifurcation transition of anisotropy at approximately 2 K, where the a-axis anisotropy splits into two special orientations, revealing a hidden bi-oriented in-plane AFM order deviating from the high-symmetry direction by 7 degree. More intriguingly, the characteristics of the bifurcated anisotropy are clearly illustrated through vector magnetic field modulation, revealing three distinct in-plane domain phases in the transverse magnetic field phase diagram. Our results not only provide valuable insights into the tunable bifurcation of magnetic anisotropic in GdTi3Bi4, but also pave a novel pathway for AFM spintronics development.
Charge density waves (CDWs) are a widespread collective electronic order in quantum materials, furnishing key insights into symmetry breaking and competing phases. However, their dynamic control with external fields remains a pivotal challenge. Here, we report deterministic and hysteretic switching of unidirectional CDW orientation via in-plane magnetic field rotation in magnetic kagome metal GdTi3Bi4. Atomically resolved spectroscopy shows two types of 3a0*1a0 CDW domains, Q1 and Q2 oriented 60 degree apart along two distinct crystallographic directions and separated by atomically sharp domain walls. Rotating the magnetic field drives reversible transitions between these CDW configurations, exhibiting a robust C2-symmetric phase diagram with pronounced hysteresis. This hysteretic switching is mediated by a field-dependent reorientation of underlying antiferromagnetic spins, revealing a tunable energy landscape with stable and metastable states and modulates the electronic charge order via spin-lattice coupling. Our findings not only demonstrate the switching of CDW configurations by in-plane magnetic field but also reveal the mechanism of coupling between CDW and magnetic fields, offering new insights into CDW manipulation and versatile platform for developing a spin-mediated multistate spin-charge coupling memory and programmable quantum devices.
Twist-engineering of topological phases in two-dimensional materials offers a powerful route to modulate electronic structure beyond conventional strain or chemical control. In particular, group 15 (pnictogens) monolayers such as bismuthene provide an ideal platform due to their strong intrinsic spin-orbit coupling (SOC) and robust topological character. Here, we investigate a previously unexplored heterostructure consisting of a $β$-bismuthene monolayer rotated by 30$^\circ$ on a planar bismuthene layer stabilized on a SiC(0001) substrate. Using first-principles calculations, we demonstrate that this specific rotational alignment induces a unique interlayer orbital hybridization which, combined with the strong SOC and the naturally broken inversion symmetry, gives rise to a pronounced Rashba spin-splitting, absent in the isolated monolayers. The topological nature of the system is confirmed through the calculation of the Z2 topological invariant and Spin Hall Conductivity (SHC), revealing a robust Quantum Spin Hall (QSH) phase with an enhanced topological response compared to the individual layers. Furthermore, we explore the chemical tunability of this system via Sb substitution, showing that the gradual reduction of SOC systematically narrows the band gap while preserving the non-trivial topology. Our results establish large-angle twisted group 15 heterostructures as a versatile platform for engineering spin-orbit-driven phenomena and advancing topological spintronics.
Many-body effects in condensed matter yield novel quantum states when the electronic density of states is enhanced. A vivid example is flat bands, which suppress kinetic energy and let interactions dominate, when they are filled with an integer number of electrons in moire systems. Yet flat bands and commensurate fillings are not the only conditions for correlated phenomena. Situations may occur where the band structure develops locally enhanced density of states, leading to strong correlations even at non-integer fillings, although such cases often yield pseudogaps that make detection elusive. Here we demonstrate that small-angle twisted monolayer-bilayer graphene combines moire-induced global flat band and additional local band flattening. Their coexistence allows direct comparison of correlated effects. The global route stabilizes commensurate states, while the local mechanism produces nearly flat bands, lifting degeneracy and generating symmetry breaking at non-integer fillings, yet without opening a global gap. Because there is no global gapped signature, the system remains metallic, but the effect reveals itself in anomalous Hall responses, signaling time-reversal symmetry breaking and valley polarization. Our results demonstrate dual-flatness as a guiding principle, extending moire physics beyond commensurate fillings and identifying topological transport as a probe of gapless correlated metals.
The coexistence of topological states with different dimensionalities in a single crystalline system offers a unique platform to study the interplay of distinct fermionic excitations. Here, integrating first-principles calculations with symmetry analysis, we propose the three-dimensional boron allotrope $P6_3$-$\text{B}_{30}$ as an ideal, structurally stable candidate for exploring multidimensional topological physics. Benefiting from the practically negligible spin-orbit coupling of the light-element framework, $P6_3$-$\text{B}_{30}$ operates as a pristine spinless topological semimetal. We show that the combined time-reversal and twofold screw symmetry ($\mathcal{T}S_{2z}$) enforces a robust two-dimensional nodal surface on the $k_z = π$ plane via a Kramers-like degeneracy. Concurrently, the system hosts a diverse set of zero-dimensional Weyl fermions -- including an unconventional double-Weyl point ($\mathcal{C} = -2$), conventional Type-I WPs ($\mathcal{C} = -1$), and completely tilted Type-II WPs ($\mathcal{C} = +1$) -- emerging at the high-symmetry points $Γ$ and K, as well as along the H-K path, protected by $C_6$ and $C_3$ crystalline rotational symmetries. Crucially, the substantial momentum-space separation between the nodal surface and Weyl points allows for their unambiguous independent resolution. Calculations of the (100) surface states reveal distinct, nontrivial Fermi arcs connecting Weyl nodes of opposite chirality. This work establishes $P6_3$-$\text{B}_{30}$ as a compelling material platform for investigating the physics of multidimensional hybrid topological fermions and their interplay.
Kitaev chains realized in quantum dots coupled via superconducting segments provide a controllable platform for engineering Majorana zero modes (MZMs). In these systems, subgap states in the hybrid region mediate the effective coupling between quantum dots and determine the emergence of sweet-spots where MZMs are strongly localized. However, existing minimal treatments often oversimplify the mesoscopic hybrid region. We perform a full microscopic treatment of this hybrid segment, capturing the quasiparticle continuum and spin-split Andreev bound states (ABSs), and show that it fundamentally alters the minimal picture. We derive analytical expressions for the renormalized couplings and sweet-spot conditions, establishing a direct link between microscopic chain parameters and Majorana optimization and identifying experimentally relevant regimes for improved device performance. Critically, we find that parity-crossings of the ABS, marking the onset of an odd-parity spin-polarized regime in the segment, identify the optimal operating windows where MZMs are simultaneously well localized with a large gap to excited states.
We study a Su-Schrieffer-Heeger chain coupled to a single mode photonic cavity. Considering an off-resonant regime we use the high-frequency expansion in order to obtain an effective fermionic Hamiltonian with cavity-mediated interactions. We characterize the effects of the cavity on topology in a finite size chain by studying three different markers adapted for interacting systems: correlation functions between edges in a chain with open boundary conditions, and a winding number based on the single-particle Green's function and bulk electric polarization via the many-body formula by Resta for a chain with periodic boundary conditions. There is excellent agreement between the winding number and polarization approaches to compute the phase diagram, with the presence of the edge states being confirmed through the calculations of the two-point correlation function. Our approach provides an alternative perspective on cavity-modified topological phases through a study of an effective interacting electronic Hamiltonian and complements methods that treat the full light-matter Hamiltonian directly.
Topologically nontrivial magnetic textures such as skyrmions offer promising opportunities for spintronic applications. In recent years, it has been shown that the magnetic properties of layered materials can be affected by depositing chiral molecules on the surface, while the influence of chiral overlayers on skyrmion properties such as their stability and interactions remains largely unexplored. To address this challenge, we employ wide-field nitrogen-vacancy (NV) magnetometry to directly image skyrmions in chiral-molecule-functionalized magnetic thin films, enabling quantitative mapping of magnetic stray fields over extended areas under ambient conditions. Using pixel-resolved optically detected magnetic resonance (ODMR) combined with controlled magnetic fields, we reproducibly nucleate and probe skyrmion states in CoFeB ferromagnetic samples, enabling quantitative investigation of their properties. We find evidence for enantioselective and magnetic-field-polarity-dependent modifications of skyrmion diameter, spacing, and shape, pointing to a possibility of molecular control of topological spin textures via magneto-chiral coupling.
Nonlinear modal interactions in resonant systems govern a wide range of phenomena, with broad relevance across modern physics and engineering. Yet, experimentally determining the strength of nonlinear coupling in multimode resonators remains highly challenging. Here, we introduce a multi-tone spectroscopy method for identifying nonlinear coupling coefficients directly from experimental data. Our approach employs dual-tone excitation near selected resonances which, in combination with additional probing tones at higher-order modes, generates sideband responses associated with specific modal couplings. These spectral signatures are analyzed using an inverse reconstruction procedure to quantitatively determine the corresponding nonlinear coupling strengths in the frequency domain. Using this method, we determine ten pairwise nonlinear coupling parameters across five modes of highly tensioned nanostrings, enabling the reconstruction of fully experimental, device-specific nonlinear reduced-order models. Our experimentally derived models show excellent agreement with values obtained numerically using finite element based nonlinear reduced-order models. Our method is generic and can be used for the characterization of diverse modal and intermodal couplings in mechanical and hybrid resonant systems.
Quantum geometry, a quantum mechanical quantity comprised of Berry curvature and quantum metric, describes the geometric structure of the electronic bands in solids. The correlation between nontrivial quantum geometry and quantum materials leads to new findings in condensed matter systems. Here we demonstrate that altermagnets, with spontaneously broken time-reversal (T)- half-lattice-translation and parity-time symmetry, host both T-odd and T-even quantum geometric quantities that simultaneously manifest themselves despite the vanishing net magnetization. Consequently, giant room-temperature third-order electrical transport responses with sizable quantum geometric contributions are observed in (101)-oriented RuO2 thin films, an altermagnetic candidate; in particular, the third-order Hall effect is intimately correlated with altermagnetic order and can serve as a promising tool for detecting the Neel vector. Our work not only supports the existence of altermagnetism in 8-nm-thick RuO2 thin films, but also shows altermagnets as a versatile platform for exploring quantum geometry and constructing quantum electronic and spintronic devices.
We study the role of heterogeneous volatility in a networked wealth dynamics model and its impact on the wealth condensation transition. Extending the Bouchaud--Mézard framework, we introduce binary volatility in networks and investigate how its configuration affects the effective power-law tail exponent of the wealth distribution. Using a stochastic block model, we control the mixing between volatility groups and show that the effective exponent is governed not only by the global parameter $Λ=2J/β^2$ but also by the volatility configuration in the network. We find that local interactions between nodes with different volatility induce a neutralization of group-wise exponents, which lowers the aggregate tail exponent and can drive a condensation transition across $γ_{\rm c}=2$. Our results identify volatility mixing as another control mechanism for wealth condensation and highlight the importance of noise heterogeneity in nonequilibrium systems on networks.
The new family of superconducting nickelates centered around La$_{3}$Ni$_{2}$O$_{7}$ possesses attractive features, such as the high transition temperature and the presence of an antiferromagnetic ground state at ambient pressure, suggesting an unconventional pairing mechanism. In the nonsuperconducting state, the possibility of different density-wave orders with opposite pressure dependencies is discussed, whose relationships and microscopic origins are largely unknown. However, sample-quality issues, such as impurity-phase formation or oxygen vacancies, impede the progress in the field. Here, we employ high-pressure synthesis and hydrostatic high-pressure transport techniques to investigate bilayer nickelates with controlled band width and filling, and perform a systematic study on their impact on the superconductivity and other characteristic properties. While increasing the tilting of the NiO$_6$ octahedra shifts the superconducting phase to higher pressure, simultaneous hole doping reverts this trend. We also observe up to three distinct anomalies in the nonsuperconducting state which are possibly related to density-wave formation.
Structural disorder in the magnetically frustrated pyrochlore system leads to intriguing magnetic states. We present the thermodynamic behavior and short range magnetic correlations in Tb$_{2}$Zr$_{2}$O$_{7}$ and Tb$_{2}$Zr$_{1.5}$Ti$_{0.5}$O$_{7}$ compounds. The parent compound Tb$_{2}$Zr$_{2}$O$_{7}$ has defect fluorite structure, which evolves toward the pyrochlore phase on Ti doping at Zr site. There is no long range magnetic order down to 0.4 K, and a magnetic field dependent spin freezing evolves below 1.25 K and 1.05 K for the parent and doped compounds, respectively. The ac susceptibility measurements indicate slow spin relaxation process below 20 K in these compounds. Inelastic neutron scattering reveals broad diffuse scattering, indicative of short range correlations at low temperature, owing to local structural distortions and persistent spin fluctuations. These results suggest a correlated, disorder influenced magnetic state in Tb$_{2}$Zr$_{2}$O$_{7}$, Tb$_{2}$Zr$_{1.5}$Ti$_{0.5}$O$_{7}$ compounds.
We derive the asymptotic first passage time (FPT) distribution for space-dependent variable-order time-fractional diffusion, where the fractional exponent $α(x)$ varies with position. For any sufficiently smooth $α(x)$ on a finite domain with absorbing and reflecting boundaries, we show that the survival probability decays as $Ψ(t)\sim C\,t^{-α_*}/(\ln t)^ν$, where $α_*$ is the minimum value of the fractional exponent and $ν$ is determined by the location and shape of the minimum. For a constant fractional exponent $ν=0$ and this provides a theoretical prediction that can identify spatially heterogeneous anomalous transport in experiments. We validate the theory against exact Laplace-space solutions and Monte Carlo simulations for linear and nonlinear profiles of $α(x)$.
We developed a machine learning interatomic potential (MLIP) for Ge-rich GeSbTe alloys of interest for applications in phase change memories embedded in microcontrollers. The MLIP was generated by fitting with a neural network method a large database of energies and forces computed within density functional theory of elemental, binary, stoichiometric and non-stoichiometric ternary alloys in the Ge-Sb-Te phase diagram. The MLIP is demonstrated to be highly transferable to large regions of the phase diagram around the compositions included in the dataset. The MLIP is then exploited to simulate the crystallization with phase separation of three Ge-rich alloys on the Ge-Sb$_2$Te$_3$ and Ge- Ge$_2$Sb$_2$Te$_5$ tie-lines that correspond to the set process of the memory cell. The transformation on the ns time scale and at 600 K, comparable to the operation conditions of the memory, yields crystalline cubic GeTe slightly Sb-doped and amorphous GeSb and Ge. These metastable phases differ from the thermodynamically stable products and form due to kinetics effects on the short time span of the set operation in phase change memories.