2026-03-23 | | Total: 26
Atomic-scale disorder can create hidden optical anisotropy even in crystals that are structurally cubic on average. Here, we show that CuInSnS$_4$ single crystals host locally symmetry-broken environments arising from intrinsic In/Sn cation disorder, which affect vibrational and excitonic properties in markedly different ways. Combining polarization- and temperature-dependent Raman spectroscopy, infrared near-field microscopy, steady-state and time-resolved photoluminescence, and first-principles calculations, we find that phonons remain largely symmetry-averaged and locally homogeneous on the nanoscale. In contrast, photoluminescence reveals a lower-energy band-tail emission with pronounced polarization anisotropy following a well-defined angular symmetry, highlighting the strong sensitivity of excitonic states to local symmetry breaking. This phonon-exciton decoupling reveals that intrinsic disorder can localize excitons while preserving vibrational coherence and dielectric homogeneity, thereby opening new opportunities for polarization-sensitive light sources, anisotropic photodetectors, and exciton-based optical functionalities even in nominally cubic multinary semiconductors.
Strain engineering provides a powerful route for tuning the electronic properties of two-dimensional (2D) materials, but exploring the full multidimensional strain space with density functional theory (DFT) is computationally prohibitive due to the nonlinear coupling between normal and shear components. In this work, we introduce a Transformer-based, multi-target surrogate model framework that achieves DFT-level bandgap prediction accuracy, reaching a mean absolute error of 0.0103 eV while retaining full interpretability through attention-weight analysis. The learned self-attention map consistently identifies shear strain as the interaction center that influences both bandgap and phonon stability, an insight not readily captured by classical feature-importance metrics. This work establishes attention-based architectures as physically interpretable surrogate models for multi-property prediction, offering a generalizable strategy for accelerating deep elastic strain engineering in materials informatics.
In solid-state batteries, ceramic solid electrolytes are penetrated by dendrites when plating above a critical current density $J_\mathrm{crit}$. A dendrite will propagate by metal deposition at a pre-existing dendrite tip if the mechanical energy required to crack the ceramic open is less than the electrical energy (Joule heating) wasted by forcing the current to detour around the dendrite to the flat electrode surface. Based on this principle of minimal power dissipation, a dependence of $J_\mathrm{crit}\propto c_\mathrm{max}^{3/2}$ is derived. $c_\mathrm{max}$ is the length of the longest preexisting, sufficiently thin interfacial defect. Consequentially, scattering of $J_\mathrm{crit}$ between samples must follow a Weibull-distribution, similar to the tensile strength of ceramic components.
Thermomigration is the driving force for hydrogen transport due to a temperature gradient. It can compete with hydrogen transport induced by stress gradients. While stress-driven hydrogen migration is well established, thermomigration remains comparatively underexplored, largely due to limited mechanistic understanding and a scarcity of experimental data. In this work, we develop a thermodynamically consistent framework for hydrogen transport, incorporating a mechanistic model for thermomigration. This is implemented within a finite element framework using an effective chemical potential. Using case studies of iron and nickel heat exchangers and zirconium alloy nuclear fuel cladding, we quantify the competing and synergistic effects of thermomigration and stress-driven transport. We show that thermomigration often dominates hydrogen redistribution in heat-carrying components, even in the presence of significant thermal incompatibility stresses. However, stress-driven transport is shown to become decisive near sharp stress concentrators. A graphical method is introduced to rapidly identify the dominant transport mechanism without requiring fully coupled simulations. The results provide practical guidance for assessing hydrogen redistribution and embrittlement risk in heat-carrying structural components.
Atomic electron tomography (AET) determines the three-dimensional (3D) coordinates and chemical identities of individual atoms from a series of scanning transmission electron microscopy images taken at different tilt angles. However, under the low dose conditions required to mitigate beam damage, the reduced signal-to-noise ratio forces a trade off among accuracy, robustness, and throughput, which ultimately limits the broader application of AET. Here, we introduce a physics aware, two stage neural networks (PANN) that incorporates physical constraints throughout its workflow to achieve accurate AET under low-dose imaging. First, a global local 3D ResUNet refines the initial reconstruction and corrects geometric distortions in the volume. Second, the local density around each identified atom is encoded using 3D Zernike moments. These feature descriptors, along with the atomic coordinates are then processed by a graph attention Transformer to classify the elemental species. We benchmark the PANN workflow using a dataset of 42,588 reconstructed volumes, covering diverse noise models, materials morphologies, and dose settings. Under low dose conditions, PANN significantly improves performance, reducing the atomic coordinates error and leading to an increase in the atomic recovery rate. The framework's performance on experimental lose dose AET data across nanoparticles of varying morphology and composition demonstrate robust generalization. We anticipate this approach will extend the applicability of AET, particularly in investigating materials sensitive to electron dose or chemical state, including halide perovskites, zeolite, and quantum dot.
Differential scanning calorimetry (DSC) is a standard tool for studying precipitation and phase transformations in aluminum alloys, yet its relation to mechanical performance has so far remained mostly indirect. Here, we demonstrate that DSC curves themselves act as fingerprints that directly encode mechanical properties. Four representative 6xxx series alloys (Al-Mg-Si) were subjected to different natural and artificial aging regimens, followed by DSC heat-flow measurements and tensile testing. Machine learning models trained on the thermograms predicted yield strength, ultimate tensile strength, and uniform elongation in five-fold grouped cross-validation, with the best model (Lasso) achieving R^2 values of 0.93, 0.86, and 0.87 and mean absolute errors of 14.3 MPa, 11.1 MPa, and 1.5 percent, respectively. Leave-one-alloy-out evaluation with sparse calibration using anchor samples further demonstrated generalization across alloy chemistries. While direct prediction on unseen alloy data degraded performance substantially, inclusion of as few as one to two anchor conditions from the target alloy recovered predictive accuracy, approaching that of the standard cross-validation. Feature importance analysis revealed that the 230 to 270 C region, associated with precipitation of the primary hardening phase beta'', contributed most strongly to predictive accuracy, providing direct mechanistic validation of the model. These findings establish DSC as a diagnostic tool that can serve as a rapid proxy for mechanical property evaluation, enabling accelerated alloy screening, process optimization, and integration of thermal analysis into data-driven manufacturing.
Hydroxyapatite (HA) on a magnesium (Mg) surface is studied using density functional theory, to help understand the effect of HA coating and alloying in the surfaces of Mg-based biodegradable implants. We determine the adsorption energies and structural changes of a single layer of HA on pure Mg(0001) and on sparsely calcium (Ca) or zinc (Zn) doped Mg(0001) and find that both Zn and Ca doping improves the adsorption, except in a few positions of HA relative to the dopant position. All adsorption configurations, whether with pure or doped Mg surfaces, show deformation of the surface and HA layer. For Ca doping, we found that for a certain adsorption configuration, the dopant Ca atom moves out of the Mg surface and into the HA layer, leaving behind a Mg vacancy in the top layer of the Mg surface. Plots of electron density changes show that electrons accumulate around the Ca dopant and the neighboring Mg atoms, while in Zn doping this is less pronounced. Overall, our results demonstrate that the dopant choice and relative position of HA influence the interaction between HA and Mg-surfaces, and affect both adsorption energies and atomic and electronic structures.
The interaction of oxygen with cobalt and cobalt-based alloys has been a very important topic in the field of spintronics as it leads to enhanced orbital anisotropy and interfacial Dzyaloshinskii-Moriya interaction (DMI), which are crucial in the context of applications such as magnetic tunnel junctions (MTJs) based data storage and domain wall (DW) motion. To understand the complex and interesting relationship between oxygen and ferromagnetic (FM)/heavy metal (HM) interfaces, we studied controlled oxygen ion implantation in a cobalt layer located in a Pt/Co 1.2 /Pt (nm) multilayer with a specific structure. At high implantation fluence, the perpendicular anisotropy was lost, as verified by in-plane hysteresis measurements. Under low magnetic field conditions, the DW dynamics of Co/Pt multilayers were analyzed, highlighting key parameters such as DW velocity, roughness amplitude, and roughness exponent. After O+-ion implantation, the DW velocity increased by more than 50 times, rising from 5 um/s to 300 um/s compared with the as-deposited multilayer. The fundamental cause of this improvement is the structural and magnetic changes brought by the implantation, which successfully lower the energy barriers preventing DW movements. The results show how oxygen implantation can be used to precisely tailor the ferromagnetic interfaces, leading to promised improvements in the functionality of next-generation spintronic devices.
Stable and remarkable valley polarization effect is the key to utilizing valley degree of freedom in valleytronic devices. According to first-principles calculations and symmetry analysis, we reveal that valley polarization effect in monolayer V2Se2O altermagnet is correlated with the net magnetic moment between magnetic V atoms under uniaxial strain, thereby proposing two strategies for achieving giant valley polarization effect. Firstly, substituting one V atom in V2Se2O with Cr to construct a ferrimagnetic monolayer VCrSe2O enhances the net magnetic moment between magnetic atoms, thereby realizing a giant valley polarization effect. Applying uniaxial strain along either the a-axis or b-axis significantly increases the value of valley polarization, which exhibits a nearly linear relationship with the net magnetic moments between the magnetic atoms. Secondly, constructing a van der Waals heterostructure composed of V2Se2O and α-SnO monolayers breaks mirror symmetry, thereby inducing a net magnetic moment, which in turn causes a remarkable valley polarization effect. Compressing the interlayer distance of the heterostructure can increase the net magnetic moment between V atoms, then enhancing the value of valley polarization to nearly 400 meV. This work reveals that valley polarization in monolayer altermagnet is correlated with the net magnetic moment between magnetic atoms. Finally, we propose two strategies to achieve giant valley polarization based on monolayer altermagnets, providing theoretical guidance for the potential applications of ferrimagnetic monolayers and altermagnet-based heterostructures in valleytronics.
In this paper, we examine the properties of spin-polarized interfaces consisting of single-molecule magnet bis(cyclopentadienyl)cobalt(II) (cobaltocene) and two-dimensional magnetic materials, semiconducting CrI$_3$ and metallic Fe$_3$GeTe$_2$, using first-principles density functional theory based calculations. Our calculated adsorption energies indicate the stability of these hetero-interfaces with the observation of hybridization of electronic states across the interface. Magnetic exchange interaction parameters have been obtained from both total energy differences and the Liechtenstein-Katsnelson-Antropov-Gubanov (LKAG) formalism in the basis of maximally localized Wannier functions (MLWFs). Analysis of these parameters shows a strong directional anisotropy in the magnetic substrate-molecule interaction in agreement with the nature of orbital hybridization. Additionally, possible exchange mechanisms are proposed based on orbital-resolved exchange and hopping parameters. We also show that the molecular adsorption may enhance the intralayer exchange interactions, with some exchange parameters reaching up to a 3-fold increase in magnitude compared to the freestanding case. Finally, we observe a 100 % spin polarization at the Fermi level in the cobaltocene/CrI$_3$ interface, which makes it particularly promising for spin-transport applications.
The orthorhombic crystal structure of the MgSnN$_2$ compound with Pna2$_1$ symmetry has been investigated as a low-cost, non-toxic material for photovoltaic (PV) applications using density functional theory (DFT) and spectroscopic limited maximum efficiency (SLME) calculations. A detailed analysis of the electronic and optical properties was performed using the mBJ semilocal exchange functional. The bandgap of MgSnN$_2$ is found to be 2.45 eV. SLME photovoltaic analysis suggests that a thin film of MgSnN$_2$ with a thickness of 2 $μ$m can reach an efficiency of 13.17% at room temperature. This efficiency was further improved through the simulation of a multi-junction device, where the tandem configuration increased the efficiency from 12.80% (single-junction) to 22.42%. Furthermore, introducing cation disorder can further reduce the bandgap, enhancing its suitability for solar cell applications.
Recent advances in nanoscale fabrication enable atomic-scale manipulation of two-dimensional (2D) materials by introducing engineered pores and perforations. This provides new opportunities to tailor functional properties of 2D materials for applications such as selective ion transport, desalination membranes, and molecular filtration. Despite this progress, the auxetic mechanical behavior of perforated 2D materials has received little attention. In this work, large-scale reactive molecular dynamics (MD) simulations, validated against experimental measurements and first-principles calculations, are employed to investigate the mechanical response of perforated monolayer titanium-based MXene metamaterials. Architectures containing rectangular perforations with straight ligaments and sinusoidally curved ligaments are systematically examined under uniaxial tension and compression over a range of geometric parameters and temperatures, from the onset of deformation to fracture. The results demonstrate that MXene metamaterials exhibit a tunable negative Poisson's ratio (NPR), which can be controlled through the perforation geometry and surface termination. Atomistic stress analysis reveals alternating in-plane shear stresses at the junctions that induce rotational deformation of the ligaments. This rotating-junction mechanism is coupled with out-of-plane deflections arising from the low bending rigidity of atomically thin materials, producing complex three-dimensional deformations. Comparison with graphene metamaterials indicates that the perforation geometry governs qualitative auxetic trends, whereas intrinsic material properties determine quantitative responses. These findings identify MXenes as a versatile candidate for the design of tunable 2D mechanical metamaterials and provide atomistic insight into the interplay between geometry, bending rigidity, and auxetic deformation mechanisms.
We report emergent magnetic behavior in heterostructures composed of (111)-oriented La$_{0.7}$Sr$_{0.3}$MnO$_3$ (LSMO) and (00$l$)-oriented Bi$_2$Te$_3$ (BT), controlled by interfacial reconstructions. When BT is deposited directly onto LSMO, an intermediate interfacial layer forms between the two materials. Polarized Neutron Reflectometry modeling reveals that this reconstructed region stabilizes a secondary magnetically ordered phase that is coupled to the underlying ferromagnetic LSMO layer. As a consequence, the heterostructures exhibit unconventional self-crossing magnetic hysteresis loops at room temperature, characterized by a reversal of the net magnetization at low applied fields. In contrast, the introduction of a tellurium seed layer results in a sharper LSMO-BT interface, while preserving the anomalous hysteresis behavior and enhancing the saturation magnetization. Element-specific X-ray absorption spectroscopy suggests that the emergent magnetic phase originates from the chemical reconstruction of manganese species. These results demonstrate that interface engineering in magnetic oxide-topological insulator heterostructures provides a pathway to control emergent magnetic coupling and emergent magnetic states in oxide-topological insulator heterostructures.
Nonlinear phononics provides a powerful ultrafast route to control lattice excitations, enabling access to hidden quantum orders, phononic computing, and quantum transduction. However, dynamic control of anharmonic phonon interactions remains limited, as these interactions are typically fixed by the equilibrium crystal lattice and lack external tunability. Emergent ferrons in ferroelectrics, which are collective oscillations of the spontaneous electric polarization, may offer a promising platform to overcome this limitation by combining intrinsic phononic nonlinearity with direct electrical control of the ferroelectric order parameter. Here we report electrically controllable nonlinear ferron upconversion in the van der Waals ferroelectric NbOI2. We show that resonant THz excitation of a 3.1 THz ferron drives coherent upconversion to a 7.0 THz optical phonon. Using two-dimensional THz spectroscopy, we directly resolve off-diagonal coupling features and establish the nonlinear upconversion pathway. Supported by first-principles calculations and analytical modeling, we identify the microscopic origin as a cubic anharmonic lattice coupling. Importantly, in situ electric-field switching enables nonvolatile control of both the ferron dynamics and the associated upconversion process. The phase reversal and hysteretic behavior across the coercive fields establish that the ferron-mediated nonlinear phononic interaction is strongly dependent on the underlying ferroelectric order parameter. These results introduce ferron upconversion as a new and universal regime of nonlinear phononics in ferroelectrics and establish an electrically programmable platform for coherent lattice control, paving the way for ferronic information processing and quantum phononic transduction.
The ST1 diamond color center was experimentally demonstrated to involve a substitutional oxygen atom (O$_C$) and carbon vacancy (V$_C$), has a spin singlet ground-state, and a metastable electron spin ancilla: a triplet. ST1's structure was left unsolved for more than a decade. With embedded multiconfigurational quantum mechanical theory, we investigate O$_C$-V$_C$-derived diamond defects, specifically both 0 and +2-charged coupled O$_C$V$_C$, and O$_C$ surrounded by V$_C$s along the [110] axis (V$_C$O$_C$V$_C$). We found both O$_C$V$_{C}^{2+}$ (C$_{3v}$) and V$_C$O$_C$V$_{C}^{2+}$ (C$_{2v}$) to have a spin-singlet ground state (1$^1$A$_1$) and metastable spin triplets. We demonstrate ST1 to be V$_C$O$_C$V$_{C}^{2+}$. The calculated vertical excitation energies of V$_C$O$_C$V$_{C}^{2+}$'s first (1$^1$B$_2$) and second (2$^1$A$_1$) bright spin-singlet excited states closely match ST1's experimental zero phonon line (2.2-2.3 eV). O$_C$V$_{C}^{2+}$ ($^1$E) absorbs much higher (2.8 eV). The two O lone pairs favor V$_C$O$_C$V$_C$ over O$_C$V$_C$, in a similar manner as the single N lone pair favors formation of N$_C$V$_C$ centers.
Machine-learning electronic Hamiltonians achieve orders-of-magnitude speedups over density-functional theory, yet current models omit long-range Coulomb interactions that govern physics in polar crystals and heterostructures. We derive closed-form long-range Hamiltonian matrix elements in a nonorthogonal atomic-orbital basis through variational decomposition of the electrostatic energy, deriving a variationally consistent mapping from the electron density matrix to effective atomic charges. We implement this framework in HamGNN-LR, a dual-channel architecture combining E(3)-equivariant message passing with reciprocal-space Ewald summation. Benchmarks demonstrate that physics-based long-range corrections are essential: purely data-driven attention mechanisms fail to capture macroscopic electrostatic potentials. Benchmarks on polar ZnO slabs, CdSe/ZnS heterostructures, and GaN/AlN superlattices show two- to threefold error reductions and robust transferability to systems far beyond training sizes, eliminating the characteristic staircase artifacts that plague short-range models in the presence of built-in electric fields.
In the Comment by Bordag et al. (arXiv:2506.10792), concerns are raised regarding the validity of the results presented in Phys. Rev. B 111, 115428 (2025) (arXiv:2403.02279), where the theoretical descriptions of the electric conductivity of graphene obtained from the Kubo formula and from quantum field theory via the polarization tensor are compared. In this Reply, we show that these concerns arise from misinterpretations of Phys. Rev. B 111, 115428 (2025), in which the results are either inaccurately represented or applied outside the domain of validity of the model. We address the comments concerning the derivation of the Luttinger formula for the electric conductivity from the Kubo formula and clarify why the results of Phys. Rev. B 111, 115428 (2025) cannot be arbitrarily extended to make claims on the gauge invariance. We further demonstrate that our findings are fully consistent with the established and widely accepted literature cited in the Comment. We confirm that the model for electric conductivity discussed in Phys. Rev. B 111, 115428 (2025) correctly predicts a vanishing electric current in the absence of an external electric field, as physically required, and in contrast with the model advocated by the Authors of the Comment. We also show that the electric permittivity does not exhibit a double pole in $ω$, contrary to the claim made in the Comment. Finally, we emphasize that the inclusion of losses is a standard and well-established approach in the study of transport properties of materials, including graphene, and we take the opportunity to correct a few minor typographical errors in Phys. Rev. B 111, 115428 (2025). We show and maintain that all results derived in Phys. Rev. B 111, 115428 (2025) are fully valid and correct.
We propose a framework to explore the internal charge distribution of mesoscopic quasiparticles by inelastic x-ray scattering, while also accounting for the conventional scattering from electrons. Specifically, we investigate a new contribution of intrinsic and optically pumped excitons (bound electron-hole pairs) to the x-ray scattering spectrum of transition metal dichalcogenides (TMDCs). The optical excitation leads to the creation of Wannier exciton populations, adding new quasi-elastic processes beyond the conventional electronic features to the x-ray scattering spectra. Differential spectra (with and without optical pumping) can be used to isolate and identify the internal charge distribution of the optically pumped excitons in the scattering response, potentially offering insights into many-body interactions and quasi-particle dynamics in 2D systems.
We present an in-depth analysis of the orbital magnetism by means of the so-called modern theory based on the Berry phase across distinct classes of materials-d transition metals, sp metals, and transition metal dichalcogenides-highlighting the microscopic nature of band structure characteristics. We adopt a gauge-covariant formulation of the modern theory proposed in [Lopez et al. Phys. Rev. B 85, 014435 (2012)], which enables the calculation of orbital magnetism in a controlled manner in any chosen gauge of Wannier functions and gives the total contribution as a gauge-invariant measurable. This captures consistently the contributions due to the anomalous position, velocity, and orbital angular momentum of Wannier basis, as well as the contributions due to Hamiltonian such that their sum is gauge-invariant. For d transition metals, we find that the atom-centered approximation captures the majority of the total contribution given by modern theory, which we attribute to localized nature of d electrons. However, 5d metals tend to exhibit larger deviation between the two methods than 3d metals do, as 5d electrons are more delocalized than 3d electrons. On the other hand, sp metals exhibit a strong deviation between the two methods, where large kinetic energy of sp electrons is important. Finally, in 1H-MoS2, we find that the valley orbital moment far exceeds the atomic limit of d electrons due to coherent hybridization between valence and conduction bands in direct band gaps. Our work elucidates the interplay of the chemical nature of electronic orbitals and the effect of band structures in a consistent manner and highlights the role of Berry phase in orbital magnetism. The results suggest a promising direction of orbitronics beyond controlling atomic orbitals, in which the orbital magnetism can be greatly enhanced by exploiting Berry phase.
We propose a magnon-driven anomalous Hall effect in altermagnets, arising from the coupling between coherently excited chiral magnons and chiral electronic motion. Using density-matrix perturbation theory and symmetry analysis, we show that the resulting Hall conductivity is solely determined by the chiralithy of the Néel-order precession, in sharp contrast to the anomalous Hall effect from the equilibrium Néel order. It then has distinct symmetry requirements from the latter and can exist even when the latter is forbidden by symmetry. The magnon-driven anomalous Hall effect is exemplified in a minimal lattice model with the same symmetry of the altermagnet CrSb, which hosts no static anomalous Hall effect. Our results reveal a direct interplay between chiral magnons and chiral electronic motion, paving the way of probing magnon chirality and to control electronic chirality through magnons.
We investigate how anisotropic strain reorganizes commensurate moiré superlattices and electronic structure in twisted bilayer graphene across a finite range of twist angles. Motivated by experiments demonstrating robust magic angle phenomenology under angular disorder and heterostrain, we construct exact commensurate moiré supercells generated by general anisotropic strain applied to the top graphene layer. At fixed effective moiré deformation, with particular focus on the $\pm 6.008^\circ$ twist angle, we uncover two distinct and symmetry inequivalent commensurate geometries: tilted two dimensional moiré superlattices and quasi one dimensional stripe like patterns. Anisotropic strain qualitatively reshapes the low energy band structure by reducing the number of Dirac points in the moiré Brillouin zone, leading to sharply different electronic and magnetic field responses in these regimes. Strikingly, tilted two dimensional structures near pristine angles preserve bandwidths and AA region localization comparable to the unstrained case, providing a natural explanation for the persistence of magic angle physics over a finite distortion window, whereas quasi one dimensional moiré patterns exhibit dimensional reduction and immediate Hofstadter butterfly splitting at infinitesimal magnetic fields. Our results identify anisotropic strain as a unifying geometric mechanism controlling commensurate moiré physics beyond the pristine twist angle limit.
Converting mid-infrared (MIR) radiation to visible or near-infrared wavelengths is essential for imaging and sensing, yet achieving sensitive, low-power, and scalable detection remains challenging. Lanthanide nanocrystals provide an alternative through ratiometric luminescence but are typically constrained by Boltzmann statistics, which tie population distributions to lattice temperature and limit signal contrast. Here we show that MIR irradiation rebalances dissipative relaxation pathways, driving lanthanide emitters into a non-Boltzmann steady state that enables non-thermal control of population distributions. This allows emission behaviors inaccessible under thermal equilibrium. We exploit this regime to achieve linear MIR detection with respect to MIR power across 6.8 to 8.6 micrometers. The ratiometric response is intrinsically independent of the pump power, enabling operation at an ultralow excitation power of 10 uW, several orders of magnitude lower than conventional approaches. Using standard silicon photodetectors, we then demonstrate room-temperature MIR imaging with detection limits approaching 4 nW um-2. Our results establish lanthanide nanoparticles as an efficient platform for MIR conversion and sensing in nanophotonic systems.
Spin-orbit coupling (SOC) gives rise to complex magnetic states such as spin liquids, skyrmion crystals, and topological spin-wave excitations. We consider exchange interactions in multi-orbital Mott insulators where SOC is strong on ligand ions. SOC on the ligands enables electron hopping accompanied by spin flips and fluctuations in the orbital state of the ligand hole. These processes generate anisotropic exchange interactions and greatly increase the number of possible exchange paths. The number grows further with the inclusion of hopping between ligands, which mediates interactions between more distant spins. We propose an effective method to calculate exchange interactions at arbitrary separations between spins. Applying it to monolayer CrI$_3$, we obtain anisotropic interactions between nearest-neighbor and next-nearest-neighbor Cr spins, as well as single-ion anisotropy induced by long-range hopping. In this material, magnetic anisotropy stabilizes long-range ferromagnetic order and opens a magnon gap at the Dirac points, which defines a nontrivial magnon band topology. Using Hubbard model parameters from first-principles calculations, the resulting spectrum agrees well with the spin-wave dispersion observed experimentally in bulk CrI$_3$, except that the calculated Dirac gap is much smaller.
Energy renormalizations of resident carriers and excitons are studied theoretically, and compared with recent experiments of electrostatically-doped WSe$_2$ monolayers. The calculated energy renormalization of resident carriers, subjected to strong out-of-plane magnetic field, reveals the importance of dynamical screening in transition metal dichalcogenides. The energy renormalization of tightly bound excitons is analyzed through the exchange interaction between the electron (or hole) component of the exciton and resident carriers that share the same spin and valley quantum numbers. Our theory explains the weak energy shift of excitonic resonances despite the strong energy renormalization of resident carriers. We identify the dependence of the energy renormalization on the envelope function of a tightly-bound exciton, showing that unlike free electron-hole pairs, this energy renormalization is not the added renormalizations of a resident electron and resident hole.
Memory is a foundational concept across disciplines, from neurobiology and electronics to artificial intelligence and quantum gravity. In materials, memory effects typically arise from ferroic orders, such as ferroelectricity and ferromagnetism, where information is stored in charge or spin degrees of freedom. Here, we report a surprising discovery of a nonvolatile superlattice memory effect in monolayer TaIrTe4, a dual quantum spin Hall insulator, where information is encoded through sharply contrasting lattice periodicities. In particular, in a pristine monolayer, we observe the spontaneous emergence of a long-period superlattice that can be programmed ON and OFF in a nonvolatile manner by electrostatic tuning of low-energy electronic states. This switching toggles the system between two structural configurations with unit cell areas differing by nearly two orders of magnitude. Mechanistically, our results reveal two independent and distinct instabilities, one in the lattice and the other in the QSH electrons, which are coupled, leading to electrostatic control of lattice configurations with nonvolatile memory. This finding is enabled by combining linear and nonlinear transport measurements, Raman spectroscopy, and scanning tunneling microscopy, which probe complementary aspects of the underlying orders. Remarkably, this nonvolatile memory effect stabilizes a spontaneous superlattice with a periodicity on the few-nanometer scale that remains robust across a wide doping range, persists over days, and survives above 70 K. Combined with the QSH topology, this stability offers a promising route to nonvolatile memory control of topological flat bands and their filling enabled quantum states. Our preliminary data indeed show the emergence of new insulating states at fractional superlattice fillings, which can be clearly switched ON and OFF together with the superlattice.