Floquet Bessel Tuning of Chern Transport in Altermagnet–Topological-Insulator Interfaces

1. Introduction

Altermagnets (AMs) [1,2] provide a route to spin-split electronic structure and Berry-curvature control without net magnetization. This property is particularly attractive for nanoscale spintronic devices, where stray magnetic fields, Joule heating, and optical Floquet excitation are undesirable. In an AM/topological-insulator (AM–TI) heterostructure, the interfacial exchange field generated by a thin altermagnetic layer can act as an anisotropic Dirac mass for the TI surface states. If this mass is controlled mechanically, a nanoscale Hall-bar device can in principle switch between distinct Chern sectors using strain rather than magnetic-field reversal. The relevant experimental observable is the transverse Hall response ${\sigma }_{xy}$, defined by $j_x={\sigma }_{xy}{E}_y$, together with its thermoelectric counterpart ${\alpha }_{xy}$. In altermagnets, the momentum-dependent spin splitting, constrained by crystal symmetry, generates Berry-curvature multipoles that can be reoriented by strain. Materials such as ${\mathrm{RuO}}_2$, MnTe, and $\alpha$${\mathrm{-Fe}}_2{\mathrm{O}}_3$ exhibit strain-tunable anomalous Hall responses at frequencies far below conventional THz optical drives [3,4,5,6]. Mechanical control therefore offers a low-dissipation route to topological transport modulation without optical pumping or global magnetization reversal.

We consider a mechanically driven AM–TI nanoscale heterostructure in which the crystalline phase of two altermagnetic electrodes is periodically modulated. Piezoelectric or surface-acoustic-wave actuation induces a controlled time dependence in the interfacial exchange field acting on the TI surface states. The resulting structure operates as a Chern valve, where the sign structure of the Dirac mass is tuned in momentum space by independently controlling the crystalline phases of the left and right electrodes. A concrete device geometry and driving protocol are shown in Figure 1a.

At low energies, the TI surface states are described by the Dirac Hamiltonian

$$\mathcal{H}\left(\mathbf{k}\right)=v_{\mathrm{F}}(k_x{\sigma }_y-k_y{\sigma }_x)+m(\theta ,\phi ){\sigma }_z,$$

(1)

where $v_{\mathrm{F}}$ is the Fermi velocity, ${\sigma }_i$ are Pauli matrices in spin space, and $\theta =\mathrm{atan}2(k_y,k_x)$ denotes the polar angle in momentum space. The parameter $\phi$ represents the in-plane crystalline phase of the altermagnetic electrode, setting the orientation of the interfacial exchange anisotropy relative to the TI momentum coordinates.

We model the interfacial exchange field as

$$m(\theta ,\phi )=m_0+m_{2}cos(2\theta -2\phi )+m_{4}cos(4\theta -4\phi ),$$

(2)

where $m_2$ and $m_4$ are effective second- and fourth-order angular harmonics generated by interface symmetry reduction and strain-induced anisotropy. These coefficients characterize interfacial symmetry rather than a bulk spin-group classification.

Chiral edge channels appear when the masses induced by the two electrodes satisfy $m_L\left(\theta \right)m_R\left(\theta \right)<0$ for some $\theta$. In this case the occupied Dirac band acquires a nonzero Chern number $\mathcal{C}$ leading to quantized Hall conductance. Independent control of ${\phi }_L$ and ${\phi }_R$ therefore enables topological switching without reversing a global magnetization. Time dependence is introduced via

$$\phi \left(t\right)={\phi }_0+\mathcal{A}cos(\Omega t),$$

(3)

with modulation amplitude $\mathcal{A}$ and drive frequency $\Omega$. In the experimentally relevant regime of slow mechanical driving, the dynamics are governed by an adiabatic Floquet Hamiltonian obtained by cycle averaging. This produces a Bessel-function renormalization of each angular harmonic, $m_n\to m_n{J}_0\left(n\mathcal{A}\right)$, allowing selective suppression of individual components at the zeros of $J_0$.

Unlike static phase control considered previously, the amplitude dependence of the Bessel factors enables controlled suppression of specific altermagnetic harmonics, producing re-entrant transitions between distinct Chern sectors (e.g., $\mathcal{C}=2\to 1\to 0$) within a single device geometry.

Finally, cyclic modulation of the relative phase $\Delta \phi ={\phi }_R-{\phi }_L$ generates quantized adiabatic pumping characterized by a winding number defined over one driving period [16,23]. The resulting topology arises from cyclic evolution in parameter space rather than quasienergy band inversions, establishing AM/TI interfaces as a mechanically driven Floquet platform for Hall and thermoelectric switching.

2. Theory and Methods

2.1. Theory

2.1.1. Floquet Formalism in Extended Hilbert Space

For a time-periodic Hamiltonian $\mathcal{H}\left(t\right)=\mathcal{H}(t+T)$ with $T=2\pi /\Omega$, the dynamics can be formulated in the extended Floquet (Sambe) space [7,8]. The corresponding Floquet Hamiltonian reads

$${\mathcal{H}}_F=\mathcal{H}\left(t\right)-i\hslash {\partial }_t,$$

(4)

which, in Fourier representation, couples photon sectors separated by integer multiples of $\hslash \Omega$.

In the regime relevant here, the mechanical frequency satisfies $\hslash \Omega \ll \Delta E$, where the minimum instantaneous gap is defined as

$$\Delta E=\underset{\mathbf{k},t}{min}2\sqrt{{\left(v_{F}k\right)}^2+m^2(\theta ,t)}.$$

(5)

The minimum occurs near $k=0$ in proximity to mass sign changes. For representative ${\mathrm{RuO}}_2{\mathrm{/Bi}}_2{\mathrm{Se}}_3$ parameters [3,28], we obtain $\Delta E=0.2$–3 meV. With $\Omega /\left(2\pi \right)=10$ MHz, $\hslash \Omega \approx 0.041\mu$eV, yielding $\hslash \Omega /\Delta E\sim {10}^{-4}$–${10}^{-5}$. In this limit, inter-sector hybridization in Sambe space is perturbatively small and does not induce quasienergy gap closings. The effective Hamiltonian reduces to the cycle average

$${\mathcal{H}}_F^{\left(0\right)}={\displaystyle \frac{1}{T}}{\int }_0^{T}dt\mathcal{H}\left(t\right),$$

(6)

with corrections of order ${(\hslash \Omega /\Delta E)}^2$.

2.1.2. Cycle-Averaged Mass Renormalization

For the phase modulation $\phi \left(t\right)={\phi }_0+\mathcal{A}cos(\Omega t)$, time averaging of each angular harmonic gives

$${\displaystyle \frac{1}{T}}{\int }_0^{T}dtcos\left[n(\theta -{\phi }_0-\mathcal{A}cos\Omega t)\right]=J_0\left(n\mathcal{A}\right)cos\left[n(\theta -{\phi }_0)\right],$$

(7)

where $J_0$ is the zeroth-order Bessel function [16,17]. The effective Dirac mass becomes

$$\overline{m}\left(\theta \right)=m_0+m_2{J}_0\left(2\mathcal{A}\right)cos(2\theta -2{\phi }_0)+m_4{J}_0\left(4\mathcal{A}\right)cos(4\theta -4{\phi }_0).$$

(8)

Under $C_2:\theta \to \theta +\pi$, both harmonics are invariant. Under $C_4:\theta \to \theta +\pi /2$, the fourth-order harmonic remains invariant whereas the second-order term changes sign. The coexistence of $m_2$ and $m_4$ therefore signals reduced interfacial symmetry rather than a bulk spin-group classification.

2.1.3. Adiabatic Pumping

For a fixed momentum angle ${\theta }^{*}$, the complex mass product $z\left(t\right)=m_L({\theta }^{*},t)m_R({\theta }^{*},t)$ traces a closed trajectory in the complex plane over one driving period. The winding number [16,23]

$$\mathcal{W}={\displaystyle \frac{1}{2\pi }}{\int }_0^{T}dt{\displaystyle \frac{d}{dt}}arg\left[z\left(t\right)\right]$$

(9)

counts the number of encirclements of the origin. Here time parametrizes cyclic evolution of $\phi \left(t\right)$ while ${\theta }^{*}$ is fixed; the topological invariant is therefore defined in parameter space rather than momentum space.

2.2. Methods

2.2.1. Model and Parameters

The TI surface states are described by a two-band Dirac Hamiltonian $\mathcal{H}\left(\mathbf{k}\right)=\mathbf{d}\left(\mathbf{k}\right)·\sigma$ with $\mathbf{d}=(-\hslash v_F{k}_y,\hslash v_F{k}_x,m(\theta ,\phi ))$ and $m(\theta ,\phi )=m_0+m_{2}cos(2\theta -2\phi )+m_{4}cos(4\theta -4\phi )$. Representative parameters for ${\mathrm{RuO}}_2{\mathrm{/Bi}}_2{\mathrm{Se}}_3$ heterostructures [3,28] are $v_F=5.0\times {10}^5$ m/s, $m_0=0.2$ meV, $m_2=8.0$ meV, $m_4=3.0$ meV, temperature $T=15$ K, and chemical potential $\mu =1.5$ meV.

2.2.2. Momentum-Space Discretization

Berry-curvature and transport integrals are evaluated on a polar $(k,\theta )$ grid with $N_k=300$ radial points up to $k_{max}={10}^9{\mathrm{m}}^{-1}$ and $N_{\theta }=601$ angular points. All integrals use composite trapezoidal quadrature [31]. Doubling $N_k$ and $N_{\theta }$ changes ${\sigma }_{xy}^{\mathrm{tot}}$ and ${\alpha }_{xy}^{\mathrm{tot}}$ by less than $5\times {10}^{-4}$ for representative parameter sets.

2.2.3. Berry Curvature

For the two-band model the Berry curvature of the lower band is

$${\Omega }_z\left(\mathbf{k}\right)=-{\displaystyle \frac{1}{2}}{\displaystyle \frac{\mathbf{d}·({\partial }_{k_x}\mathbf{d}\times {\partial }_{k_y}\mathbf{d})}{{\left|\mathbf{d}\right|}^3}}.$$

(10)

The Chern number is calculated as

$$\mathcal{C}={\displaystyle \frac{1}{2\pi }}\int d^{2}k{\Omega }_z\left(\mathbf{k}\right),$$

(11)

Numerically, Berry curvature and Chern-sector diagnostics were obtained using two complementary procedures. For the transport calculations we used the continuum two-band expression above, with analytic derivatives of $\mathbf{d}\left(\mathbf{k}\right)$ and a small norm regularization near gap closings. For the Floquet Chern phase diagrams we also used a gauge-invariant Fukui–Hatsugai–Suzuki link-variable discretization of the occupied-band spinors, which avoids gauge ambiguities in the numerical evaluation of $\mathcal{C}$. Since the mass depends on $\theta$, derivatives ${\partial }_{k_x}m$ and ${\partial }_{k_y}m$ are evaluated via ${\partial }_{\theta }m$ together with ${\partial }_{k_x}\theta =-k_y/k^2$ and ${\partial }_{k_y}\theta =k_x/k^2$. To avoid numerical instabilities near gap closings, $\left|\mathbf{d}\right|$ is regularized as $\left|\mathbf{d}\right|\ge max({\epsilon }_{\mathrm{reg}}max|\mathbf{d}|,{10}^{-30})$ with ${\epsilon }_{\mathrm{reg}}=5\times {10}^{-6}$.

2.2.4. Hall and Thermoelectric Response

The anomalous Hall conductivity follows from Berry-curvature integration weighted by the Fermi–Dirac occupation [32,33],

$${\sigma }_{xy}(\theta ;\mu ,T)=-{\displaystyle \frac{e^2}{\hslash }}{\int }_0^{k_{max}}{\displaystyle \frac{kdk}{{\left(2\pi \right)}^2}}{\Omega }_z(k,\theta )f(E_{-}(k,\theta )-\mu ),$$

(12)

with $f(E-\mu )=1/[1+exp((E-\mu )/k_{B}T)]$ and $E_{-}(k,\theta )=-\sqrt{{(\hslash v_{F}k)}^2+m^2(\theta ,\phi )}$. Angular averaging yields ${\sigma }_{xy}^{\mathrm{tot}}$.

The thermoelectric Hall coefficient is obtained from the Mott relation,

$${\alpha }_{xy}=-{\displaystyle \frac{{\pi }^2{k}_B^{2}T}{3e}}{\displaystyle \frac{\partial {\sigma }_{xy}}{\partial \mu }},$$

(13)

where ${\partial }_{\mu }{\sigma }_{xy}$ is evaluated by central finite differences with $\Delta \mu =0.5$ meV. The plotted dimensionless Hall response is ${\tilde{\sigma }}_{xy}={\sigma }_{xy}/(e^2/h)$.

2.2.5. Chern-Valve Criterion

For electrode phases ${\phi }_L$ and ${\phi }_R={\phi }_L+\Delta \phi$, sign-change sectors are defined by $m_L\left(\theta \right)m_R\left(\theta \right)<0$. Each connected sector corresponds to a domain wall in the effective Dirac mass and therefore to one chiral edge channel (Jackiw–Rebbi mechanism [29,30]). Connected components on the periodic $\theta$ grid determine the channel number $N_{\mathrm{ch}}$. To suppress numerical noise, a minimum angular width $w_{min}=6^{\circ }$ and a small threshold $m_L{m}_R<-{\eta }^2$ with $\eta ={10}^{-6}$ eV are imposed.

2.2.6. Phase Modulation and Diagrams

Periodic modulation $\phi \left(t\right)={\phi }_0+\mathcal{A}cos(\Omega t)$ is implemented at the level of the cycle-averaged Hamiltonian, $m_n\to m_n{J}_0\left(n\mathcal{A}\right)$. Phase diagrams are obtained by scanning $\Delta \phi \in [-\pi /2,\pi /2]$ with $N_{\Delta \phi }=241$ points and $\mu$ over $[-6,6]$ meV or $[-15,15]$ meV with $N_{\mu }=151$ points. Contours of ${\alpha }_{xy}^{\mathrm{tot}}=0$ are overlaid on $N_{\mathrm{ch}}$ maps.

3. Results

Figure 1b shows $\overline{m}\left(\theta \right)$ (scaled by the magneto-elastic coupling $J_1\sim 0.3$ meV [28]) for $\mathcal{A}=0$ (black), $0.6$ (red), and $1.2$ (blue), at phase offsets $\Delta \phi =0$ (solid) and $\pi /6$ (dashed). For $\Delta \phi =0$ the mass exhibits a symmetric double-well angular structure. A finite phase offset shifts the extrema but preserves the overall nodal pattern. Increasing $\mathcal{A}$ suppresses the $C_2$ and $C_4$ harmonics through the Bessel renormalization, progressively flattening the angular mass profile and reducing the number of sign-change sectors.

Figure 1c displays the Chern number $\mathcal{C}(\Delta \phi )$ for the same amplitudes. In the undriven case ($\mathcal{A}=0$), a broad plateau with $\left|\mathcal{C}\right|=2$ extends over $\Delta \phi \in [-\pi /3,\pi /3]$, with sharp transitions at the phase boundaries. At $\mathcal{A}=0.6$, close to the first zero of $J_0\left(4\mathcal{A}\right)$, the $\left|\mathcal{C}\right|=2$ region narrows significantly. For $\mathcal{A}=1.2$, near the first zero of $J_0\left(2\mathcal{A}\right)$, the response becomes predominantly trivial ($\mathcal{C}=0$), consistent with strong suppression of the leading harmonic. The sequence illustrates re-entrant topological switching controlled by the drive amplitude at fixed phase offset. The global structure is summarized in Figure 1d, which maps $\mathcal{C}(\mathcal{A},\Delta \phi )$ over $\mathcal{A}\in [0,1.4]$ and $\Delta \phi \in [-\pi /2,\pi /2]$. Topological sectors with $\mathcal{C}=\pm 2$ dominate at small $\mathcal{A}$, while intermediate amplitudes stabilize $\mathcal{C}=\pm 1$ regions. As $\mathcal{A}$ approaches the Bessel zeros, these sectors contract and eventually give way to an extended trivial phase. The resulting “eye”-shaped structure reflects the coexistence of symmetry-inequivalent $C_2$ and $C_4$ harmonics: tuning $\mathcal{A}$ renormalizes each component differently and enables transitions such as $\mathcal{C}=2\to 1\to 0$ within the same device.

Panel (e) shows the winding number $\mathcal{W}\left(\mathcal{A}\right)$ for two representative driving protocols. For co-rotating phases ($\Delta \phi =0$), $\mathcal{W}=0$ throughout the explored amplitude range. For mixed-phase driving ($\Delta \phi =\pi /3$), $\mathcal{W}$ develops a quantized plateau at $\mathcal{W}=2$ over a finite interval of $\mathcal{A}$, and returns to zero outside this window. The right subpanels provide complementary diagnostics: inside the plateau, $arg\left[m_L\left(t\right)m_R\left(t\right)\right]$ exhibits a net $4\pi$ phase accumulation per cycle and the complex trajectory $z\left(t\right)=m_L\left(t\right)m_R\left(t\right)$ forms a closed loop encircling the origin twice; outside the plateau no net encirclement occurs.

Figure 1f presents the anomalous Hall conductance ${\tilde{\sigma }}_{xy}(\Delta \phi )$ and the thermoelectric Hall coefficient ${\tilde{\alpha }}_{xy}(\Delta \phi )$. The suppression of ${\tilde{\sigma }}_{xy}$ at large $\mathcal{A}$ follows the collapse of the corresponding Chern sectors near Bessel zeros. The thermoelectric response, obtained from the Mott relation ${\alpha }_{xy}\propto {\partial }_{\mu }{\sigma }_{xy}$, is enhanced when the chemical potential is gate-tracked to remain near energies where $|{\partial }_{\mu }{\sigma }_{xy}|$ is maximal. This optimization sharpens the response without altering the underlying topological structure.

The parameter range used in Figure 1(b)–(f) is compatible with strained AM/TI nanoscale heterostructures. A realistic implementation could consist of a few-quintuple-layer ${\mathrm{Bi}}_2{\mathrm{Se}}_3$ or ${\mathrm{Bi}}_2{\mathrm{Te}}_3$ channel (∼5–10 QL) contacted by thin ${\mathrm{RuO}}_2$ or MnTe altermagnetic electrodes with thicknesses in the ∼5–30 nm range.[19,28] The layer thicknesses play distinct physical roles. The TI should be thick enough to suppress hybridization between the top and bottom surface states, while remaining sufficiently thin for electrostatic gating and interfacial transport to dominate over bulk leakage channels; this motivates the few-quintuple-layer regime. If the TI becomes too thin, hybridization opens an additional finite-size gap competing with the altermagnet-induced Dirac mass, whereas excessively thick or bulk-conducting films progressively shunt the surface Hall response through trivial bulk transport. The altermagnetic layer, in turn, must be thick enough to stabilize the compensated ordered phase and maintain a robust exchange coupling at the interface, while remaining thin enough for strain to be transferred coherently from the piezoelectric or surface-acoustic-wave actuator. In this sense, the effective mass $m(\theta ,\phi )$ should be understood primarily as an interfacial exchange mass controlled by the first few atomic layers adjacent to the TI rather than as a bulk volume-averaged magnetic parameter.

Independent strain control of $\Delta \phi$ has been demonstrated in ${\mathrm{RuO}}_2$-based and MnTe-based platforms,[5,6,26] while integrated piezoelectric or surface-acoustic-wave (SAW) actuators can generate the angular excursions required to reach $\mathcal{A}\sim 0.1$–1.0.[24,25] Such actuation can be implemented on chip using thin-film AlN, ScAlN, or ${\mathrm{LiNbO}}_3$ platforms compatible with nanoscale Hall-bar geometries and MHz operation.[36,37,38,39] Operating at MHz frequencies places the system deep in the adiabatic regime ($\hslash \Omega /\Delta E\sim {10}^{-4}$–${10}^{-5}$), which is particularly advantageous for thin-film nanoscale devices because it minimizes heating, reduces photon-assisted interband processes, and avoids the strong dissipation typically associated with THz Floquet platforms.[34,35]

Experimentally, the proposed mechanism could be identified through the strain-amplitude dependence of the anomalous Hall plateau structure, the collapse of topological sectors near Bessel zeros, and the gate dependence of the thermoelectric Hall response. Since the effect originates from interfacial Berry-curvature engineering, the signal is also expected to scale with TI-film thickness and to weaken as bulk transport progressively dominates.

4. Conclusions

These results establish AM/TI interfaces as a realistic platform for mechanically programmable Hall and thermoelectric switching at MHz frequencies, with $(\mathcal{A},\Delta \varphi )$ control windows accessible through piezoelectric strain and electrostatic gating in nanoscale thin-film heterostructures. The present bilayer geometry enables programmable topological control through independent tuning of the two altermagnetic phases $({\phi }_L,{\phi }_R)$ combined with Bessel-function amplitude renormalization. Because the mechanism relies on interfacial Berry-curvature engineering rather than global magnetization reversal or optical pumping, it is naturally suited to low-dissipation nano-spintronic architectures operating in the adiabatic regime.

Several extensions of the present framework remain open. Multilayer AM–TI heterostructures with repeated interfaces could enable spatially selective Berry-curvature engineering and coupled chiral transport channels. Likewise, multifrequency driving protocols may provide a route toward synthetic topological band structures controlled entirely by mechanical phase modulation.[16] Finally, combining the present mechanism with proximitized superconducting interfaces could offer a platform for studying strain-controlled topological superconducting boundaries and related emergent edge states.[40]