## Abstract

Antiferromagnets exhibit distinctive characteristics such as ultrafast dynamics and robustness against perturbative fields, thereby attracting considerable interest in fundamental physics and technological applications. Recently, it was revealed that the Néel vector can be switched by a current-induced staggered (Néel) spin-orbit torque in antiferromagnets with the parity-time symmetry, and furthermore, a nonsymmorphic symmetry enables the control of Dirac fermions. However, the real-time dynamics of the magnetic and electronic structures remain largely unexplored. Here, we propose a theory of the ultrafast dynamics in antiferromagnetic Dirac semimetals and show that the Néel vector is rotated in the picosecond timescale by the terahertz-pulse-induced Néel spin-orbit torque and other torques originating from magnetic anisotropies. This reorientation accompanies the modulation of the mass of Dirac fermions and can be observed in real time by the magneto-optical effects. Our results provide a theoretical basis for emerging ultrafast antiferromagnetic spintronics combined with the topological aspects of materials.

## Introduction

The control and detection of antiferromagnetic (AFM) orders have been a challenging problem because of the absence of net magnetizations. Since the energy scale of antiferromagnets lies in the terahertz range, which is beyond the scope of conventional electronics, optical pulses are suitable for controlling AFM orders^{1,2,3} through magnon excitations^{4,5,6,7}, inverse magneto-optical effects^{8,9,10}, magnetic anisotropy^{11,12,13,14}, photoinduced phase transitions^{15,16,17,18}, and so forth^{19,20,21,22,23}. For detecting AFM orders, the magneto-optical effects are frequently utilized for collinear and canted antiferromagnets^{24,25,26,27,28,29,30}. The X-ray magnetic linear dichroism with photoemission electron microscopy has succeeded in imaging the AFM domain structure^{31,32,33,34}. While these optical and X-ray measurements have the potential for time-resolved observations of magnetic dynamics^{1,35,36,37,38}, their application to the AFM orders remains difficult and under development.

Several years ago, another mechanism that modulates AFM orders via electric currents was proposed and demonstrated^{39,40}. In antiferromagnets with a combined space-inversion and time-reversal symmetry, electric currents induce a staggered spin density, and thus, a staggered torque termed the Néel spin-orbit torque (NSOT). The NSOT can efficiently switch the staggered magnetization, i.e., the Néel vector, thereby yielding significant progress in AFM spintronics^{41,42,43,44}. Recently, it was found that a nonsymmorphic crystalline symmetry in addition to the parity-time symmetry preserves the crossing of doubly degenerate bands; this implies that Dirac fermions are controlled by the direction of the Néel vector^{45}. The discovery of this close relationship between the magnetic and electronic structures opened a research field called topological AFM spintronics^{46,47}. Some experiments have shown that not only electric current pulses but optical pulses can control the AFM order^{48,49}; these studies performed transport measurements such as the anisotropic magnetoresistance and the planar Hall effect to observe the magnetic structure^{44}. Meanwhile, the topological electronic structure has received less attention, and therefore, the real-time dynamics of the Néel vector and the Dirac fermions induced by the NSOT are still unclear.

From a theoretical viewpoint, two approaches have been established to study the real-time dynamics of magnets. One is based on microscopic models such as the Hubbard model and Kondo lattice model^{18}, which allows considering electron correlation and quantum effects; however, the cluster size is severely limited and approximations are needed to deal with the symmetry-broken ordered states. The other approach is to use the Landau–Lifshitz–Gilbert (LLG) equation of the classical vectors of magnetic moments, which has been widely adopted in spintronics^{50,51,52,53,54,55}. This equation enables the simulation of large-scale magnetic structures, while the electronic degree of freedom is integrated out. Therefore, a unified framework is required to reveal the relationship between the magnetic and electronic structures.

In this work, we propose a theory based on a minimal two-dimensional model of AFM Dirac semimetals to investigate the real-time magnetic and electronic dynamics. We show that the Néel vector is rotated by optically induced NSOT and other torques originating from magnetic anisotropies. This reorientation accompanies the modulation of the mass of Dirac fermions. Furthermore, the magneto-optical effect is found to be promising for the time-resolved measurement of the Néel vector.

## Results

### Theoretical model

We consider a minimal tight-binding model of AFM Dirac semimetals projected onto the two-dimensional square lattice. The Hamiltonian is divided into three parts as \({{{\mathcal{H}}}}={{{{\mathcal{H}}}}}_{{{{\rm{ele}}}}}+{{{{\mathcal{H}}}}}_{{{{\rm{exc}}}}}+{{{{\mathcal{H}}}}}_{{{{\rm{mag}}}}}\). The first term

consists of the nearest-neighbor (*h*_{1}) and next-nearest-neighbor (*h*_{2}) hoppings, and the spin-orbit coupling (*λ*) of itinerant electrons with momentum **k**. Here, the spin and sublattice degrees of freedom of the itinerant electrons are described by the Pauli matrices **σ** and **τ**, respectively. The localized magnetic moments on each sublattice denoted by **S**_{A} and **S**_{B} couple to the itinerant electrons through the exchange interaction

Hereafter, we define the uniform and staggered magnetizations as **m** = (**S**_{A} + **S**_{B})/2 and **n** = (**S**_{A} − **S**_{B})/2, respectively. The vector **n** is termed the Néel vector. The exchange Hamiltonian is an extension of the model for **S**_{A} = − **S**_{B} = **n** proposed by Šmejkal et al.^{45}. While both the parity (\({{{\mathcal{P}}}}\)) and time-reversal (\({{{\mathcal{T}}}}\)) symmetries are broken because of the presence of the localized magnetic moments, the parity-time (\({{{\mathcal{PT}}}}\)) symmetry is preserved when **m** = 0. This guarantees the Kramers degeneracy of the energy bands in the whole Brillouin zone. In addition to the \({{{\mathcal{PT}}}}\) symmetry, this model is invariant under a nonsymmorphic glide symmetry operation only when the Néel vector is directed along the 〈100〉 directions, which causes the crossings of the doubly degenerate energy bands, i.e., the gapless Dirac points, to be protected at the Brillouin zone boundary^{45}, as shown in Fig. 1a. The energy gap at the Dirac points depends on the value of *J*_{exc} and the direction of **n**(*t*). The spin-orbit coupling and the AFM ordering of the localized moments lead to the spin-momentum locking of the itinerant electrons. The expectation value of the staggered spin moments of the itinerant electrons in the lower bands is displayed in Fig. 1b, which indicates that the net staggered moment 〈*τ*^{z}*σ*^{x}〉 takes a negative value when **n** = (1, 0, 0) because of the AFM exchange interaction.

To stabilize the initial state in which the Néel vector is aligned with the [100] direction, we introduce magnetic anisotropy energy

where *K*_{z} and *K*_{xy} represent coefficients of the easy-plane and biaxial anisotropies, respectively. Here, other magnetic interactions such as the exchange interaction between the localized moments, which can be implemented in a straightforward manner, are omitted because the AFM state is stabilized by the indirect exchange mediated by the itinerant electrons at half filling.

The time-dependent electric field **F**(*t*) parallel to the two-dimensional plane is incorporated in the model via the Peierls substitution. The time evolutions of the itinerant electrons and the localized magnetic moments are described by the von Neumann and LLG equations, respectively; these equations are solved simultaneously by the fourth-order Runge–Kutta method. The parameter values and units are presented in “Methods”.

### Application of a constant electric field

First, we discuss real-time dynamics induced by a constant electric field that is switched on at time *t* = 0 and applied in the [110] direction, i.e., \({{{\bf{F}}}}(t)={F}_{0}\Theta (t)(\cos \pi /4,\sin \pi /4,0)\) with Θ being the step function. Figure 2a shows the real-time profiles of the in-plane *x*-component of the Néel vector and the out-of-plane *z*-component of the uniform magnetization. Figure 2b displays the excited electron density defined by the number density of electrons that occupy single-particle energy levels above the chemical potential *μ* of the initial state. The energy bandgap shown in Fig. 2c is defined by the difference between the single-particle energies of the lowest level above *μ* and the highest level below *μ*; this difference is equal to the direct gap when there is the particle-hole symmetry with *h*_{2} = 0. As time evolves, *n*^{x} abruptly decreases towards zero at time *t* = *t*_{1}, *t*_{2}, or *t*_{3}, which is indicated by arrows in Fig. 2a and is hereafter termed a reorientation time *t*_{r}. Meanwhile, the *z* component of the uniform magnetization appears with values of the order of *m*^{z} ~ 0.03; this means that the magnetic structure is slightly canted as illustrated in Fig. 2e. As the field magnitude increases from *F*_{0} = 0.005 to 0.009, the reorientation time *t*_{r} becomes short. The excited electron density shown in Fig. 2b increases as time evolves. Figure 2b shows some plateaus in which the excited electron density is constant in time and the energy bandgap opens simultaneously (Fig. 2c). Once the excited electron density becomes greater than approximately 0.01 (indicated by the shade in Fig. 2b), the Néel vector starts to rotate. Note that the bandgap opens through the Néel-vector dynamics induced by the constant field, which is in contrast to the gap opening of the Floquet–Bloch bands under the high-frequency driving of electrons^{56}. These phenomena induced by the constant electric field occur in a timescale of the order of *t* ~ 10,000*ħ*/*h*_{1}, which corresponds to 6.6 ps for *h*_{1} = 1 eV. Therefore, the Néel vector and the energy bandgap can be controlled by the external field in such a short timescale.

The mechanism of the reorientation of the Néel vector is understood as follows. The external electric field induces the nonequilibrium staggered spin density of the itinerant electrons through the spin-orbit coupling, and it generates torques in the [001] direction (see Fig. 2d). Then, the two sublattice magnetic moments are slightly canted towards the [001] direction, which means *m*^{z} > 0. The easy-plane anisotropy, whose energy per sublattice moment is estimated to be \({K}_{z}{({m}^{z})}^{2}\), induces another effective field (0, 0, − 2*K*_{z}*m*^{z}) that rotates the localized moments in the anticlockwise direction as depicted in Fig. 2e. Here, we have shown the NSOT-driven reorientation of the Néel vector on the basis of the microscopic model in which the real-time dynamics of the itinerant electrons and the localized moments are explicitly considered.

Figure 2f shows the polar plots of the in-plane angle of the Néel vector *φ* as a function of time multiplied by the field magnitude for different values of *F*_{0}. The Néel vector rotates towards the [010] direction for sufficiently large field magnitudes *F*_{0} ≳ 0.004. For the weak magnitude *F*_{0} ≲ 0.004, the trajectories of *φ*(*F*_{0}*t*) are on a universal curve independent of *F*_{0}; this implies that the timescale of the dynamics is determined by the vector potential **A** ~ **F***t* rather than the electric field **F**. This is because the NSOT is induced by the electric current due to the diamagnetic response. The reorientation is regarded as a deviation from the weak-field universal trajectory. Even before the reorientation, the Néel vector slightly deviates from the *x* axis because of the NSOT; once the excited electron density exceeds approximately 0.01, the NSOT-induced staggered torque shown in Fig. 2d overcomes the other torques originating from the magnetic anisotropies. These behaviors are summarized in Fig. 2g, where the reorientation time is plotted as a function of *F*_{0} for different values of the biaxial anisotropy *K*_{xy} and the Gilbert damping constant *α*. Overall, *t*_{r} is inversely proportional to *F*_{0}, which means that the vector potential ~ *F*_{0}*t*_{r} governs the reorientation dynamics; *t*_{r} becomes shorter as *K*_{xy} or *α* decreases. It is also found that the product *t*_{r}*K*_{xy} is a function of *F*_{0}/*K*_{xy} for *K*_{xy }≥ 0.01 and *α* = 1. The reorientation time shows the step-like feature indicated by the dotted lines, which reflects a discontinuous behavior of *φ*(*F*_{0}*t*) with respect to *F*_{0} as seen in the inset of Fig. 2f. Although the anisotropy energies and the Gilbert damping constant are highly material dependent, this reorientation can occur in the timescale of subpicoseconds with a sufficiently large *F*_{0} (≳ 0.01) irrespective of *K*_{xy} and *α*.

### Irradiation of a monocycle pulse

We show the real-time dynamics induced by a monocycle pulse of which the frequency is in the terahertz range but is off-resonant (see “Methods” for parameters). The electric field is parallel to the [110] direction. Figure 3a shows the in-plane components of **n** and the out-of-plane component of **m** for *F*_{0} = 0.02. The Néel vector rotates towards the [010] direction during pulse irradiation, and it stops the rotation after irradiation. This reorientation accompanies the modulation of the bandgap as shown in Fig. 3b. When the bandgap is opened, the excited electron density does not increase as in the case of the constant field, which means that the pulse is off-resonant.

The dashed lines in Fig. 3a shows *n*^{y} and *m*^{z} for the same field amplitude but different polarization, i.e., \({{{{\bf{F}}}}}_{0}\parallel [\bar{1}10]\), whose signs are opposite to those in the **F**_{0}∥[110] case (solid lines). This polarization dependence is also seen in the case of the constant field (not shown), and it implies that the direction of the Néel vector is determined by the polarization and inverted by the mirror operation on **F**_{0} against a plane perpendicular to **n**(*t* = 0).

The pulse amplitude dependence of the Néel vector dynamics is summarized in Fig. 3c. When the amplitude is weak (*F*_{0} < 0.01), the Néel vector cannot climb over the potential barrier of the magnetic anisotropy. However, the in-plane angle exceeds 45^{∘} for a sufficiently large *F*_{0} (≥ 0.01) for which the excited electron density becomes greater than approximately 0.01, and the Néel vector remains in the [010] direction after irradiation for even larger *F*_{0} (≥ 0.02).

Figure 3d, e show the in-plane angular velocity of the Néel vector and the *z* component of the uniform magnetization, both of which are divided by *F*_{0}. For weak *F*_{0} < 0.01, the time profiles of *d**φ*/*d*(*F*_{0}*t*) and *m*^{z}/*F*_{0} are on each universal curve, whereas the deviations from the curves are found for *F*_{0 }≥ 0.01 and they lead to the reorientation towards the [010] direction. These time profiles are quite similar because *d**φ*/*d**t* is approximated by 2*K*_{z}*m*^{z} as long as ∥**m**∥ ≪ ∥**n**∥.

### Optical conductivity and magneto-optical effect

Finally, we show the optical longitudinal and transverse (Hall) conductivities for different Néel vector angles, and we discuss the resulting magneto-optical Voigt effect, assuming that the system is in the collinear AFM ground states. Figure 4a, b display the real parts of the longitudinal (*σ*_{xx}) and transverse (*σ*_{xy}) conductivities as functions of the in-plane Néel vector angle *φ*. The peaks of \({{\mathrm{Re}}}\,{\sigma }_{xx}(\omega )\) for *ω* > 0 are attributed to the Van Hove singularity; the Drude peak is also seen at *ω* = 0 for **n**∥[100] and [010] when *h*_{2} ≠ 0. Because of the two-fold symmetry around the *z* axis, *σ*_{xx} and *σ*_{xy} are periodic in *φ*. Here, *σ*_{yy} is equivalent to *σ*_{xx} under the rotation *φ* ↦ *φ* + *π*/2, and *σ*_{xy} is equal to *σ*_{yx}. The sign of the transverse conductivity is unchanged by *λ* ↦ − *λ* or *J*_{exc} ↦ − *J*_{exc} (not shown), which implies the conductivity is an even function of *λ* and *J*_{exc}. Since the present model holds the D_{4h} symmetry in the limit of *J*_{exc}**n** → 0, the conductivity linear in the electric field can be expanded with respect to a time-reversal-breaking field **B** ~ *J*_{exc}**n** as \({\sigma }_{ij}({{{\bf{B}}}})={\sigma }_{ij}^{(0)}+{\sigma }_{ijk}^{(1)}{B}_{k}+{\sigma }_{ijkl}^{(2)}{B}_{k}{B}_{l}\), where *σ*^{(1)} vanishes when **B**⊥[001] and

is a symmetric tensor. The off-diagonal components of *σ*(**B**) can be written as \({\sigma }_{xy}({{{\bf{B}}}})={\sigma }_{yx}({{{\bf{B}}}})\propto \sin 2\varphi\) for \({{{\bf{n}}}}=(\cos \varphi ,\sin \varphi ,0)\). This is consistent with the numerical result of *σ*_{xy}, although the present values of *λ* (=0.8) and *J*_{exc} (=0.6) are beyond the perturbative regime. From these properties of *σ*_{xx} and *σ*_{xy}, we can identify the direction in which the Néel vector rotates by measuring the longitudinal optical conductivity and magneto-optical signal. Figure 4c, d shows the magneto-optical rotation angle and ellipticity calculated from *σ*_{xx} and *σ*_{xy} (see “Methods”). The rotation angle *θ* exhibits large values up to 1^{∘}. In real materials, the values of *λ* and *J*_{exc} may be smaller than the values adopted here, which causes a reduction in the rotation angle down to a few millidegrees. Even so, the magneto-optical rotation can be measured using modern optical apparatus.

The transverse conductivity *σ*_{xy} appears when *λ* ≠ 0 and *J*_{exc} ≠ 0 even though the magnetic structure is collinear AFM. Recently, the transverse response in collinear antiferromagnets is attracting considerable attention, and it has been discussed for the case in which the conductivity tensor is antisymmetric^{57,58}; our calculations demonstrate the transverse optical response attributed to the symmetric tensor in collinear antiferromagnets.

## Discussion

We investigated the real-time dynamics of the magnetic and electronic structures induced by the constant and pulse electric fields, and we revealed the microscopic mechanism of the NSOT-driven reorientation of the Néel vector. Our method can be applied to any tight-binding model coupled to localized moments. The experimental verification of the present results is necessary for a better understanding of the NSOT-induced dynamics and for further progress in ultrafast AFM spintronics combined with the topological aspects of the materials. Candidate materials include orthorhombic CuMnAs^{45,59,60,61} and MnPd_{2}^{62}, which have been studied in the field of AFM spintronics and are expected to host Dirac fermions. Note that the reorientation of the Néel vector demonstrated in this study is not directly related to what has been observed in experiments because of the discrepancy in the timescale. To make more quantitative discussions on, e.g., the timescale of the reorientation and the required magnitude of the field, we should consider the crystal structure and dimensionality of materials, obtain realistic estimates of model parameters as well as the Gilbert damping constant from the ab initio calculations and experiments, and take other dissipation processes caused by the electron–electron interaction into account; these are left for future work.

Recent intense terahertz pulse sources, whose oscillation period is 1 ps and peak amplitude is more than 1 MV cm^{−1}^{63}, are promising to drive the reorientation of the Néel vector in the picosecond timescale. In our calculations, a typical value of *F*_{0} is of order 0.01, which corresponds to 0.26 MV cm^{−1} (see “Methods”). There are many earlier studies on the electrical modulation and detection of the Néel vector^{44}; their time resolution was limited to ≳1 ns. However, it may be possible to observe the Néel vector dynamics in real time via optical measurements because the magnetic structure is closely related to the electronic structure in the present model. One approach is the time- and angle-resolved photoemission spectroscopy; this can directly access the energy band structure, which reflects the direction of the Néel vector as the positions and energy gaps of the Dirac points or nodal lines. Another approach is to use magneto-optical effects. This system exhibits the magneto-optical rotation with angles up to 1^{∘}. The time-resolved magneto-optical measurements can be carried out even in the subpicosecond timescale, which provides a sufficient time resolution to observe the real-time dynamics discussed in this study.

## Methods

### Real-time dynamics

We calculate the real-time dynamics of the itinerant electrons and the localized magnetic moments as follows^{64,65}. The quantum state of the electrons is exactly described by a one-body density matrix \({\hat{\rho }}_{{{{\bf{k}}}}}\) because the electronic Hamiltonian \({{{{\mathcal{H}}}}}_{{{{\rm{ele}}}}}\) contains no many-body interaction. Once the localized magnetic moments are given, the Hamiltonians in Eqs. (1) and (2) can be diagonalized for each **k** as \(({{{{\mathcal{H}}}}}_{{{{\rm{ele}}}}}+{{{{\mathcal{H}}}}}_{{{{\rm{exc}}}}})\left|{{{\bf{k}}}}\nu \right\rangle ={\varepsilon }_{{{{\bf{k}}}}\nu }\left|{{{\bf{k}}}}\nu \right\rangle\), where \(\left|{{{\bf{k}}}}\nu \right\rangle\) is obtained by the unitary transform as \(\left|{{{\bf{k}}}}\nu \right\rangle ={\sum }_{\gamma = {{{\rm{A}}}},{{{\rm{B}}}}}{\sum }_{s = \uparrow ,\downarrow }\left|{{{\bf{k}}}}\gamma s\right\rangle {U}_{\nu ,\gamma s}^{* }\). In the initial state in which the external field is absent, the density matrix is given by \({\tilde{\rho }}_{{{{\bf{k}}}}}={\sum }_{\nu }\left|{{{\bf{k}}}}\nu \right\rangle {n}_{{{{\bf{k}}}}\nu }\left\langle {{{\bf{k}}}}\nu \right|\), where *n*_{kν} = Θ(*μ* − *ϵ*_{kν}) is the step function with *μ* being the chemical potential chosen such that the number of electrons is *N*_{ele}. The time evolution of the density matrix is governed by the von Neumann equation

with the initial condition \({\hat{\rho }}_{{{{\bf{k}}}}}(t=0)={U}_{{{{\bf{k}}}}}^{-1}{\tilde{\rho }}_{{{{\bf{k}}}}}{U}_{{{{\bf{k}}}}}\). The expectation value is defined by \(\langle \cdot \rangle ={N}^{-1}{\sum }_{{{{\bf{k}}}}}{{\mathrm{Tr}}}\,[\ \cdot \ {\hat{\rho }}_{{{{\bf{k}}}}}]\), where *N* represents the number of **k** points in the Brillouin zone. The localized moments on each sublattice denoted by **S**_{A} and **S**_{B} are treated as classical unit vectors, and therefore, the Hamiltonian \({{{\mathcal{H}}}}(t)\) depends on the classical spins and the external field. The classical spins follow the LLG equation

for *γ* = A, B, where *α* denotes the Gilbert damping constant and \({{{{\bf{b}}}}}_{\gamma }\equiv -\langle \partial {{{\mathcal{H}}}}/\partial {{{{\bf{S}}}}}_{\gamma }\rangle\) represents the effective field given by

Here, **σ**_{A/B} = (1 ± *τ*^{z})**σ**/2 denotes the sublattice spin density. The time-dependent electric field **F**(*t*) is incorporated in the model via the Peierls substitution \({{{\bf{A}}}}(t)=-{\int}^{t}{{{\bf{F}}}}(t^{\prime} ) dt^{\prime}\), where \({{{\bf{A}}}}(t)=-{\int}^{t}{{{\bf{F}}}}(t^{\prime} ) dt^{\prime}\), *e*, and *a* denote the vector potential, elementary charge, and lattice constant, respectively. Throughout this paper, *h*_{1}, *e*, *a*, and *ħ* are set to unity. The physical quantities of time, electric-field magnitude, and conductivity are expressed in units of *ħ*/*h*_{1} = 0.66 fs, *h*_{1}/(*e**a*) = 26.3 MV cm^{−1}, and *e*^{2}/(*ħ**a*) = 6.4 × 10^{3} S cm^{−1}, respectively, for *h*_{1} = 1 eV and *a* = 0.38 nm. In most calculations, we use *h*_{2} = 0.08, *λ* = 0.8, *J*_{exc} = 0.6, *K*_{z} = 0.1, *K*_{xy} = 0.01, and *α* = 1.

The easy-plane anisotropy coefficient *K*_{z} determines the in-plane angular velocity of the Néel vector *d**φ*/*dt*. For the above parameters, the difference in the total energies of the ground states is \(\langle {{{\mathcal{H}}}}\rangle {| }_{{{{\bf{n}}}} = (0,0,1)}-\langle {{{\mathcal{H}}}}\rangle {| }_{{{{\bf{n}}}} = (1,0,0)}=0.13{h}_{1}\). This may be considerably larger than the anisotropy energy of real materials (≲1 meV)^{60,61}. However, we confirmed that the reorientation of the Néel vector occurs within 10 ps even for smaller values: *K*_{z} = *K*_{xy} = 0.001, *λ* = 0.08, *J*_{exc} = 0.06, and *F*_{0} = 0.001.

Equations (5) and (6) are simultaneously solved using the fourth-order Runge–Kutta method with a time step *δ**t* = 0.01. We perform the numerical calculation with *N* = 512 × 512-point mesh in the Brillouin zone, and we fix the number of electrons to *N*_{ele} = 2*N* (half-filled). The external field **F**(*t*) = − ∂_{t}**A**(*t*) is given by **A**(*t*) = − **F**_{0}*t*Θ(*t*) for the constant electric field and \({{{\bf{A}}}}(t)=-{{{{\bf{F}}}}}_{0}\cdot (t-{t}_{{{{\rm{c}}}}})\exp [-{(t-{t}_{{{{\rm{c}}}}})}^{2}/(2{t}_{{{{\rm{w}}}}}^{2})]\) for the monocycle pulse. We adopt *t*_{c} = 2000 and *t*_{w} = 300.

### Optical conductivity and magneto-optical effect

We calculate the optical conductivity from

where \({{{{\bf{J}}}}}_{{{{\bf{k}}}}\mu \nu }=\left\langle {{{\bf{k}}}}\mu \right|{{{{\bf{J}}}}}_{{{{\bf{k}}}}}\left|{{{\bf{k}}}}\nu \right\rangle\) denotes the matrix element of the electric current operator \({{{{\bf{J}}}}}_{{{{\bf{k}}}}}=-\delta {{{\mathcal{H}}}}({{{\bf{k}}}}-{{{\bf{A}}}})/\delta {{{\bf{A}}}}{| }_{{{{\bf{A}}}} = 0}\) and Γ represents a broadening factor. The dielectric tensor is given by

with *ϵ*_{0} = 8.854 pF m^{−1} being the electric constant. Since the optical conductivity tensor is symmetric (*σ*_{xy} = *σ*_{yx}) and anisotropic (*σ*_{xx} ≠ *σ*_{yy}) in the present model, we cannot use the well-known formula of the complex Kerr rotation angle \({\Phi }_{{{{\rm{K}}}}}={\epsilon }_{xy}/[(1-{\epsilon }_{xx})\sqrt{{\epsilon }_{xx}}]\). According to Maxwell’s equations in material media, the electric-field amplitude **F**_{0} satisfies

where \(\hat{{{{\bf{N}}}}}\) represents a complex refractive index vector. From Eq. (11) we have the relation \(({\hat{N}}^{2}-{\epsilon }_{xx})({\hat{N}}^{2}-{\epsilon }_{yy})={\epsilon }_{xy}^{2}\), assuming that the propagation vector of incident light is perpendicular to the two-dimensional *x**y*-plane. The symmetric property of the conductivity and dielectric tensors leads to the eigenmodes of Eq. (11) being linearly polarized, contrary to the antisymmetric case in which the eigenmodes are circularly polarized. The boundary condition is that the tangential components of the electric field **F** and magnetic field \(\hat{{{{\bf{N}}}}}\times {{{\bf{F}}}}/{c}_{0}\) (*c*_{0} is the speed of light) are continuous at the surface, from which we obtain two eigenmodes and complex refractive indices. The incident linearly polarized light is decomposed into a linear combination of the eigenmodes. The amplitude vector of the reflected light is proportional to (*r*_{pp}, *r*_{sp}, 0), where *r*_{pp} and *r*_{sp} are the amplitude reflection coefficients of the incident p-polarized light to the reflected p- and s-polarized light. The magneto-optical rotation angle *θ* and ellipticity *η* are determined by the fit of the electric-field trajectory to an ellipse rotated around the *z* axis by angle *θ*, and the relation \(\eta ={{\mathrm{sgn}}}\,[\arg ({r}_{{{{\rm{sp}}}}}/{r}_{{{{\rm{pp}}}}})]\times | {r}_{{{{\rm{sp}}}}}| /| {r}_{{{{\rm{pp}}}}}|\), respectively.

## Data availability

The data that support the findings of this study are available upon reasonable request to the corresponding author.

## Code availability

The computer codes used to generate the results of this study are available upon reasonable request to the corresponding author.

## References

- 1.
Němec, P., Fiebig, M., Kampfrath, T. & Kimel, A. V. Antiferromagnetic opto-spintronics.

*Nat. Phys.***14**, 229–241 (2018). - 2.
Song, C. et al. How to manipulate magnetic states of antiferromagnets.

*Nanotechnology***29**, 112001 (2018). - 3.
Kampfrath, T., Tanaka, K. & Nelson, K. A. Resonant and nonresonant control over matter and light by intense terahertz transients.

*Nat. Photonics***7**, 680–690 (2013). - 4.
Kampfrath, T. et al. Coherent terahertz control of antiferromagnetic spin waves.

*Nat. Photonics***5**, 31–34 (2011). - 5.
Nishitani, J., Nagashima, T. & Hangyo, M. Terahertz radiation from antiferromagnetic MnO excited by optical laser pulses.

*Appl. Phys. Lett.***103**, 081907 (2013). - 6.
Baierl, S. et al. Terahertz-driven nonlinear spin response of antiferromagnetic nickel oxide.

*Phys. Rev. Lett.***117**, 197201 (2016). - 7.
Li, J. et al. Spin current from sub-terahertz-generated antiferromagnetic magnons.

*Nature***578**, 70–74 (2020). - 8.
Kimel, A. V. et al. Ultrafast non-thermal control of magnetization by instantaneous photomagnetic pulses.

*Nature***435**, 655–657 (2005). - 9.
Kimel, A. V. et al. Inertia-driven spin switching in antiferromagnets.

*Nat. Phys.***5**, 727–731 (2009). - 10.
Satoh, T., Iida, R., Higuchi, T., Fiebig, M. & Shimura, T. Writing and reading of an arbitrary optical polarization state in an antiferromagnet.

*Nat. Photonics***9**, 25–29 (2015). - 11.
Duong, N. P., Satoh, T. & Fiebig, M. Ultrafast Manipulation of Antiferromagnetism of NiO.

*Phys. Rev. Lett.***93**, 117402 (2004). - 12.
Kimel, A. V., Kirilyuk, A., Tsvetkov, A., Pisarev, R. V. & Rasing, T. Laser-induced ultrafast spin reorientation in the antiferromagnet TmFeO

_{3}.*Nature***429**, 850–853 (2004). - 13.
Baierl, S. et al. Nonlinear spin control by terahertz-driven anisotropy fields.

*Nat. Photonics***10**, 715–718 (2016). - 14.
Schlauderer, S. et al. Temporal and spectral fingerprints of ultrafast all-coherent spin switching.

*Nature***569**, 383–387 (2019). - 15.
Ju, G. et al. Ultrafast generation of ferromagnetic order via a laser-induced phase transformation in FeRh thin films.

*Phys. Rev. Lett.***93**, 197403 (2004). - 16.
Ehrke, H. et al. Photoinduced melting of antiferromagnetic order in La

_{0.5}Sr_{1.5}MnO_{4}measured using ultrafast resonant soft X-ray diffraction.*Phys. Rev. Lett.***106**, 217401 (2011). - 17.
Li, T. et al. Femtosecond switching of magnetism via strongly correlated spin–charge quantum excitations.

*Nature***496**, 69–73 (2013). - 18.
Ishihara, S. Photoinduced ultrafast phenomena in correlated electron magnets.

*J. Phys. Soc. Jpn.***88**, 072001 (2019). - 19.
Takayoshi, S., Aoki, H. & Oka, T. Magnetization and phase transition induced by circularly polarized laser in quantum magnets.

*Phys. Rev. B***90**, 085150 (2014). - 20.
Mentink, J. H., Balzer, K. & Eckstein, M. Ultrafast and reversible control of the exchange interaction in Mott insulators.

*Nat. Commun.***6**, 6708 (2015). - 21.
Ono, A. & Ishihara, S. Double-exchange interaction in optically induced nonequilibrium state: a conversion from ferromagnetic to antiferromagnetic structure.

*Phys. Rev. Lett.***119**, 207202 (2017). - 22.
Werner, P., Eckstein, M., Müller, M. & Refael, G. Light-induced evaporative cooling of holes in the Hubbard model.

*Nat. Commun.***10**, 5556 (2019). - 23.
Takasan, K. & Sato, M. Control of magnetic and topological orders with a DC electric field.

*Phys. Rev. B***100**, 060408 (2019). - 24.
Cheong, S.-w, Fiebig, M., Wu, W., Chapon, L. & Kiryukhin, V. Seeing is believing: visualization of antiferromagnetic domains.

*npj Quantum Mater.***5**, 3 (2020). - 25.
Saidl, V. et al. Optical determination of the Néel vector in a CuMnAs thin-film antiferromagnet.

*Nat. Photonics***11**, 91–96 (2017). - 26.
Tzschaschel, C., Satoh, T. & Fiebig, M. Tracking the ultrafast motion of an antiferromagnetic order parameter.

*Nat. Commun.***10**, 3995 (2019). - 27.
Higo, T. et al. Large magneto-optical Kerr effect and imaging of magnetic octupole domains in an antiferromagnetic metal.

*Nat. Photonics***12**, 73–78 (2018). - 28.
Zheng, Z. et al. Magneto-optical probe of ultrafast spin dynamics in antiferromagnetic CoO thin films.

*Phys. Rev. B***98**, 134409 (2018). - 29.
Xu, J. et al. Imaging antiferromagnetic domains in nickel oxide thin films by optical birefringence effect.

*Phys. Rev. B***100**, 134413 (2019). - 30.
Schreiber, F. et al. Concurrent magneto-optical imaging and magneto-transport readout of electrical switching of insulating antiferromagnetic thin films.

*Appl. Phys. Lett.***117**, 082401 (2020). - 31.
Stöhr, J. et al. Images of the antiferromagnetic structure of a NiO(100) surface by means of X-ray magnetic linear dichroism spectromicroscopy.

*Phys. Rev. Lett.***83**, 1862–1865 (1999). - 32.
Scholl, A. et al. Observation of antiferromagnetic domains in epitaxial thin films.

*Science***287**, 1014–1016 (2000). - 33.
Ohldag, H. et al. Spin reorientation at the antiferromagnetic NiO(001) surface in response to an adjacent ferromagnet.

*Phys. Rev. Lett.***86**, 2878–2881 (2001). - 34.
Grzybowski, M. J. et al. Imaging current-induced switching of antiferromagnetic domains in CuMnAs.

*Phys. Rev. Lett.***118**, 057701 (2017). - 35.
Wen, H., Cherukara, M. J. & Holt, M. V. Time-resolved X-ray microscopy for materials science.

*Annu. Rev. Mater. Res.***49**, 389–415 (2019). - 36.
Choe, S.-B. et al. Vortex core-driven magnetization dynamics.

*Science***304**, 420–422 (2004). - 37.
Litzius, K. et al. Skyrmion Hall effect revealed by direct time-resolved X-ray microscopy.

*Nat. Phys.***13**, 170–175 (2017). - 38.
Dean, M. P. M. et al. Ultrafast energy- and momentum-resolved dynamics of magnetic correlations in the photo-doped Mott insulator Sr

_{2}IrO_{4}.*Nat. Mater.***15**, 601–605 (2016). - 39.
Železný, J. et al. Relativistic Néel-order fields induced by electrical current in antiferromagnets.

*Phys. Rev. Lett.***113**, 157201 (2014). - 40.
Wadley, P. et al. Electrical switching of an antiferromagnet.

*Science***351**, 587–590 (2016). - 41.
Jungwirth, T., Marti, X., Wadley, P. & Wunderlich, J. Antiferromagnetic spintronics.

*Nat. Nanotechnol.***11**, 231–241 (2016). - 42.
Baltz, V. et al. Antiferromagnetic spintronics.

*Rev. Mod. Phys.***90**, 015005 (2018). - 43.
Manchon, A. et al. Current-induced spin-orbit torques in ferromagnetic and antiferromagnetic systems.

*Rev. Mod. Phys.***91**, 035004 (2019). - 44.
Amin, O. J., Edmonds, K. W. & Wadley, P. Electrical control of antiferromagnets for the next generation of computing technology.

*Appl. Phys. Lett.***117**, 010501 (2020). - 45.
Šmejkal, L., Železný, J., Sinova, J. & Jungwirth, T. Electric control of dirac quasiparticles by spin-orbit torque in an antiferromagnet.

*Phys. Rev. Lett.***118**, 106402 (2017). - 46.
Šmejkal, L., Mokrousov, Y., Yan, B. & MacDonald, A. H. Topological antiferromagnetic spintronics.

*Nat. Phys.***14**, 242–251 (2018). - 47.
Tsai, H. et al. Electrical manipulation of a topological antiferromagnetic state.

*Nature***580**, 608–613 (2020). - 48.
Olejník, K. et al. Terahertz electrical writing speed in an antiferromagnetic memory.

*Sci. Adv.***4**, eaar3566 (2018). - 49.
Kašpar, Z. et al. Quenching of an antiferromagnet into high resistivity states using electrical or ultrashort optical pulses.

*Nat. Electron.***4**, 30–37 (2021). - 50.
Maekawa, S., Valenzuela, S. O., Saitoh, E. & Kimura, T. (eds.)

*Spin Current*2nd edn.(Oxford University Press, Oxford, 2017). - 51.
Blachowicz, T. & Ehrmann, A.

*Spintronics: Theory, Modelling, Devices*. (De Gruyter, Berlin, 2019). - 52.
Abert, C. Micromagnetics and spintronics: models and numerical methods.

*Eur. Phys. J. B***92**, 120 (2019). - 53.
Hals, K. M. D., Tserkovnyak, Y. & Brataas, A. Phenomenology of current-induced dynamics in antiferromagnets.

*Phys. Rev. Lett.***106**, 107206 (2011). - 54.
Roy, P. E., Otxoa, R. M. & Wunderlich, J. Robust picosecond writing of a layered antiferromagnet by staggered spin-orbit fields.

*Phys. Rev. B***94**, 014439 (2016). - 55.
Lopez-Dominguez, V., Almasi, H. & Amiri, P. K. Picosecond electric-field-induced switching of antiferromagnets.

*Phys. Rev. Appl.***11**, 024019 (2019). - 56.
Oka, T. & Aoki, H. Photovoltaic Hall effect in graphene.

*Phys. Rev. B***79**, 081406 (2009). - 57.
Šmejkal, L., González-Hernández, R., Jungwirth, T. & Sinova, J. Crystal time-reversal symmetry breaking and spontaneous Hall effect in collinear antiferromagnets.

*Sci. Adv.***6**, eaaz8809 (2020). - 58.
Naka, M. et al. Anomalous Hall effect in

*κ*-type organic antiferromagnets.*Phys. Rev. B***102**, 075112 (2020). - 59.
Tang, P., Zhou, Q., Xu, G. & Zhang, S.-C. Dirac fermions in an antiferromagnetic semimetal.

*Nat. Phys.***12**, 1100–1104 (2016). - 60.
Kim, Y., Kang, K., Schleife, A. & Gilbert, M. J. Voltage-induced switching of an antiferromagnetically ordered topological Dirac semimetal.

*Phys. Rev. B***97**, 134415 (2018). - 61.
Xu, S.-G. et al. Dirac fermions in the antiferromagnetic spintronics material CuMnAs.

*Phys. Rev. B***102**, 125123 (2020). - 62.
Shao, D.-f, Gurung, G., Zhang, S.-h & Tsymbal, E. Y. Dirac nodal line metal for topological antiferromagnetic spintronics.

*Phys. Rev. Lett.***122**, 077203 (2019). - 63.
Fülöp, J. A., Tzortzakis, S. & Kampfrath, T. Laser-driven strong-field terahertz sources.

*Adv. Opt. Mater.***8**, 1900681 (2020). - 64.
Chern, G.-W., Barros, K., Wang, Z., Suwa, H. & Batista, C. D. Semiclassical dynamics of spin density waves.

*Phys. Rev. B***97**, 035120 (2018). - 65.
Koshibae, W., Furukawa, N. & Nagaosa, N. Real-time quantum dynamics of interacting electrons: self-organized nanoscale structure in a spin-electron coupled system.

*Phys. Rev. Lett.***103**, 266402 (2009).

## Acknowledgements

One of the authors, Sumio Ishihara, passed away on 7 November 2020 during the preparation of the manuscript. We thank H. Matsueda and Y. Masaki for their valuable comments and critical reading of the manuscript. This work was supported by JSPS KAKENHI Grant Nos. JP19K23419, JP20K14394, JP17H02916, JP18H05208, and JP20H00121. The numerical calculations were performed using the facilities of the Supercomputer Center, the Institute for Solid State Physics, the University of Tokyo.

## Author information

### Affiliations

### Contributions

A.O. conceived the idea, performed the calculations, and wrote the manuscript. A.O. and S.I. discussed and interpreted the results.

### Corresponding author

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

**Publisher’s note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Ono, A., Ishihara, S. Ultrafast reorientation of the Néel vector in antiferromagnetic Dirac semimetals.
*npj Comput Mater* **7, **171 (2021). https://doi.org/10.1038/s41524-021-00641-2

Received:

Accepted:

Published: