Introduction

Thermoelectric (TE) properties, which determine heat-to-electricity energy conversion of materials, lie at the center toward TE applications. Traditionally, semimetals are usually regarded as bad TE materials, because they naturally possess two types of carriers, which counterbalance each other’s contribution to the induced TE voltage, leading to a reduced thermopower compare to semiconductors which usually has unipolar transports1. Interestingly, recent TE transport measurements on Weyl semimetals have shown that both types of carriers can contribute to the transverse TE voltage constructively through Nernst effect when a magnetic field is applied perpendicular to the applied temperature gradient, because electrons and holes are forced to deflect in opposite directions by the magnetic force2,3,4,5,6,7. Taking advantage of Nernst effect, four times higher thermopower factor compared to conventional Seebeck response is achieved on NbP4, which is comparable to that of the state-of-the-art thermoelectric materials (see Supplementary Table 1 for a summary of thermoelectric coefficients of typical materials).

Different from steady-state TE response, which is usually slow because it involves the heat transport dynamics of lattice, the photothermoelectric (PTE) response upon photoexcitation can be ultrafast as it is dominated by the transient carrier temperature, which is induced by the highly nonequilibrium carriers excited by ultrafast laser pulses8,9,10,11,12. The instantaneous transient states after the photoexcitation can support transient ambipolar transport and provide ambipolar TE effect. As such response is usually on the picosecond timescale in semimetals9,12,13,14, such transient current is difficult to be measured by conventional transport measurement. Instead, the picosecond transient current emits electronmagnetic wave in terahertz (THz) frequency range, thus can be characterized by measuring its THz emission15,16,17.

The ultrafast PTE response is an important aspect of photoexcited carrier dynamics, which is crucial for high field/speed electronic device applications. With ultrahigh electron mobility and ultrafast carrier dynamics, topological semimetals are suitable for such applications11,13,18,19. In these devices, the carriers are accelerated by the electric field which is similar to the photoexcited state after photoexcitation9,14,20,21,22,23,24. It is equally important for optoelectronic-related applications, such as light-harvesting and detection, where light interaction with materials is involved and the simultaneous transient optical process involves the TE response of the materials8,10,11,25,26,27.

In this work, we reveal the transient ambipolar TE current response of typical Dirac semimetal Cd3As2 after ultrafast photoexcitation under magnetic field. The transient thermal gradient is created by the thickness gradient of the as-grown Cd3As2 sample after ultrafast photoexcitation. As a contactless approach, the THz emission from the transient TE current is studied at room temperature. Our results clearly distinguish the transverse Seebeck response and the Nernst response from Cd3As2 after photoexcitation. According to the THz emission amplitudes, the Nernst effect is one order of magnitude larger than the Seebeck effect with a relatively weak magnetic field B ~ 0.4 T, and the large enhancement is attributed to the constructive ambipolar transport nature under magnetic field. Furthermore, the efficiency of THz emission from Cd3As2 is comparable to that of typical semiconductor THz sources. This work may open opportunities for high-speed TE devices, taking advantage of their PTE effects.

Results

THz emissions from a Cd3As2 film with a thickness gradient

The Cd3As2 samples measured in this work are [112]-oriented thin films grown by molecular beam epitaxy28,29,30. As shown in our previous study28, the Fermi level of Cd3As2 is about 250 meV above the Dirac node, the electron doping density is on the order of 1013 cm−2, and the mobility is about 5 × 103 cm2 V−1 s−1 at room temperature. The Cd3As2 film has a radial thickness gradient due to the nature of the growth method. The central area is thicker, and the thickness gradually decreases from the center to the outer part on a circular substrate. The gradient is characterized by relative transmittance scanning with a 633-nm He-Ne laser. Figure 1a shows the thickness variation of a Cd3As2 sample deducted from the transmittance measurement. The maximum thickness of the Cd3As2 film is estimated to be around 50 nm. The relative thickness gradient is about 10% per millimeter.

Fig. 1: THz emissions from Cd3As2 epitaxial film with a thickness gradient.
figure 1

a The position of a cut piece (white dash line region) in an as-grown wafer, and the mapping of relative thickness of the Cd3As2 film determined by two-dimensional scanning of laser transmittance. The blue circle marks the excitation laser spot and the arrow illustrates the direction of thickness gradient on the spot. b Schematic diagram of experimental setup, where θ denotes the azimuthal angle of the sample with respect to y-axis. Inset: the crystal structure of Cd3As2 along (112) plane. c Schematic diagram of THz emissions from transient current generated by 800-nm optical excitation. EN and E// here denote the instantaneous electrical field of THz emission from Nernst current (JN) under magnetic field and from the in-plane current (J//), respectively. J⊥ denotes the out-of-plane current. d Typical THz waveform generated by J//, J⊥ and JN in (c). The two THz waveforms from J⊥ and JN are obtained by sample tilting and by applying an in-plane magnetic field respectively.

Full size image

Because the Cd3As2 has strong light absorption at the 800-nm excitation wavelength (~50% for 50-nm thick film31), the light intensity decays quickly when penetrating into Cd3As2, which builds up transient temperature gradient of absorbed light energy and thus the electron temperature gradient along the direction perpendicular to the surface. On the other hand, as a result of the radial thickness gradient of Cd3As2 film, the ultrafast pulse excitation can also produce a temperature gradient along the in-plane direction. This effect has been previously experimentally verified in Bi2Te3-based TE thin films16. The temperature gradient, along both out-of-plane and in-plane directions, can generated an out-of-plane (J⊥) and in-plane (J//) TE current through Seebeck effect, respectively. In our previous studies11,28, we have already demonstrated a PTE current in picosecond timescale generated in Cd3As2 when excited by ultrafast pulses. The picosecond currents emit electromagnetic waves in THz frequency, which can be detected in the electro-optical (E-O) sampling geometry.

The measurement of THz emission from ultrafast laser excitation of Cd3As2 is schematically depicted in Fig. 1b, which is based on a typical E-O sampling using [110] orientated ZnTe crystal. In a default experimental configuration at room temperature, an 800-nm 150-femtosecond laser beam excites Cd3As2 film under normal incidence with a magnetic field (when applicable) applied along the x-axis, and the E-O sampling is configured to detect y-polarized THz signal (with ZnTe[110] along x-axis). More experimental details of transient THz emission setup can be found in the method section and the detail of THz sampling geometry can be found in Supplementary Note 2.

As shown in Fig. 1c, the THz field emitted by the in-plane TE current (J//) is determined by the following relation: E//(t) ∝ dJ///dt, which propagates along the light path (z-axis). The THz field emitted by the out-of-plane TE current J⊥ cannot be detected in the E-O sampling geometry. However, when the sample is tilted, J⊥ would have projection perpendicular to z-axis, which contributes to detectable THz emission (see Supplementary Note 3). Alternatively, when an x-direction magnetic field is applied, a y-direction current JN is generated from Nernst effect as determined later on. Both approaches can emit detectable THz in the aforementioned E-O sampling geometry.

Figure 1d shows typical THz emission signals detected in E-O sampling with and without magnetic field. According to the magnitudes of THz emission signals (S(J//) and S(J⊥)) emitted from J// and J⊥, respectively, we can find that J⊥ is stronger than J// when no magnetic field is applied. Furthermore, when an x-direction magnetic field of B = 0.4 T is applied, the THz emission signal (S(JN, 0.4 T)) for normal incidence is an order of magnitude larger than that for 45°-incidence excitation at 0 T (S(J⊥, 0 T)), indicating a Nernst effect that is much larger than the Seebeck effect as result of ambipolar transport and possibly anomalous Nernst effect, which will be discussed in the discussion section.

THz emission from photo-Seebeck effect

First, we study the characteristics of THz emission from Cd3As2 without applying magnetic field. Figure 2a shows the peak amplitude of THz waveform (Smax(t), the waveforms are shown in Supplementary Fig. 4) increases with the excitation intensity (I). It can be fitted by a model of photo-Seebeck current (see Supplementary Note 5), which approximately follows a power law: ETHz ∝ I0.66. This intensity dependence is different from that of an optical rectification effect or any other second-order effects17,32,33,34,35,36,37,38,39,40,41, which typically exhibits linear excitation intensity dependence. Here we can rule out the effect from absorption saturation of the excitation beam, because a power-dependent transmittance measurement of the 800-nm excitation beam is flat within the power range as shown in Supplementary Fig. 4a. Furthermore, the THz emission is independent of the polarization of the excitation beam as shown in Fig. 2b. This is consistent with the isotropic optical properties on Cd3As2 (112) plane for near-infrared photon42,43,44. However, if the excitation polarization is fixed and the sample is rotated around the surface normal axis (varying the azimuthal angle θ, as shown in the inset of Fig. 2c), the amplitude of the THz waveform changes with θ as shown in Fig. 2c. The peak signal of THz emission exhibits a cosinoidal dependence on θ as shown in Fig. 2d. According to the cosinoidal fitting shown in Fig. 2d, we find the THz polarization direction is along the direction of thickness gradient of Cd3As2 film with accuracy of 6.2°. This azimuthal angle dependence of THz emission is consistent with that the THz emission is generated via the in-plane Seebeck current J// induced by thickness gradient16 as discussed above.

Fig. 2: THz emission from Cd3As2 without applying magnetic field.
figure 2

a The laser intensity dependence of the peak field of emitted THz pulse. The red line is the simulation of power dependence according to the TE response presented in Supplementary Note 5. b Typical THz waveforms excited by laser with different polarized direction. c Typical THz waveforms at different azimuthal angle θ, the thickness gradient is roughly set along y-axis for θ = 0°, with unintentional deviation of θ0. d The azimuthal angle dependence of the peak field of emitted THz pulse, and the red line is the cosinoidal fitting cos (θ+θ0), the deviation (θ0) is 6.2° according to the fitting. Insets in panels (b, c): the schematic illustrations of the corresponding measurements.

Full size image

THz emission by photo-Nernst effect

In the next, we study the characteristics of the THz emission measurement results when a 0.4-T in-plane magnetic field is applied. After applying the magnetic field along x-direction, the THz emission is about an order of magnitude larger. Figure 3a shows the excitation intensity dependence of THz emission amplitude in the presence of a magnetic field. Different from that at 0 T which follows the power law of ETHz ∝ I0.66, the excitation intensity dependence under magnetic field follows a different power law of ETHz ∝ I0.5. The THz emission is also independent of the polarization of the excitation beam as shown in Fig. 3b. Both power and polarization dependence support the response is due to photothermal effect instead of nonlinear optical effect, which is similar to the case without magnetic field. Figure 3c shows the dependence of the THz waveforms at different azimuthal angle θ. The THz response at 0.4 T magnetic field is composed of a θ-dependent cosinoidal oscillation component and a θ-independent component. The θ-independent component is one order of magnitude larger than the θ-dependent component. The θ-dependent component is from J//, the oscillation amplitude of the θ-dependent component is about the same for those taken with and without magnetic field as shown in Fig. 3d. The θ-independent component is induced by the magnetic field due to photo-Nernst effect.

Fig. 3: THz emission from Cd3As2 with a magnetic field of ~0.4 T at x-direction.
figure 3

a The laser intensity dependence of the peak field of emitted THz pulse. The red line is a power-law fitting. b Typical THz waveforms excited by laser with different polarization. c Typical THz waveforms at different azimuthal angle θ. d The azimuthal angle dependences of the peak field of THz pulses with and without magnetic field, and the lines are the cosinoidal fitting. The data are normalized by the average amplitude at time zero (Savg(0), 0.4 T). Insets in panels (b, c): the schematic illustrations of the corresponding measurements.

Full size image

In the following, we systematically studied the magnetic field dependence of THz emission. First, we minimized the amplitude of nonmagnetic signal by adjusting θ, therefore only the magnetic field-dependent THz component can be detected. We find the THz waveform can flip over and remain the same amplitude as the direction of the magnetic field is reversed as shown in Fig. 4a. However, if θ is set so that there is a strong nonmagnetic signal, the directly detected THz signals are obviously asymmetric under reverse magnetic fields, unless the nonmagnetic component of the signal is subtracted (Supplementary Fig. 5). This further confirms that the THz emission under magnetic field contains a pure magnetic component (S(t, B) − S(t, 0)) and a nonmagnetic component (S(t, 0)) that are independent of each other, and the direction of magnetic component is affected by the direction of the magnetic field. Moreover, when the magnetic field is 40°-away from the x-axis in xz-plane, the THz amplitude approximately equals to the THz amplitude under in-plane magnetic field multiplying by cos40°, as show in Fig. 4b. It implies that only the in-plane projection of the magnetic field influences the detected THz emission.

Fig. 4: Magnetic field dependences of THz emission.
figure 4

a The THz emissions at reverse magnetic fields of 0.4 T. The dash line is the reversion of waveform at B = ‒0.4 T. b The THz emissions when a 0.1 T-magnetic field is applied at x-axis and 40° away from x-axis in xz-plane, respectively. The dash line is the THz waveform under x-directional magnetic field multiplying by cos40°. c The x- and y-polarization components of the THz emissions with magnetic field B‖x and B‖y. d The amplitude of net magnetic THz signal at time zero as a function of the intensity of magnetic field. Insets in panels (a, b): the schematic illustrations of the corresponding measurements.

Full size image

The polarization of THz emission can be determined by measuring the orthogonal polarization THz components Sx and Sy along x- and y-axis, respectively (see the detailed method in Supplementary Note 3). As shown in Fig. 4c, for B‖x, Sy dominates the signal and Sx is close to zero, so the THz emission is y-polarized; for B‖y, the magnetic THz becomes x-polarized. We can summarize that the polarization of THz emission is perpendicular to the in-plane magnetic field. The residual signal of Sx at B‖x and the THz amplitudes difference between B‖x and B‖y should be due to experimental errors, eg. ZnTe [001] or B not exactly along x/y-axis, and the inhomogeneity of ZnTe crystal. Practically, these errors are difficult to be removed completely. Finally, as shown in Fig. 4d, the pure magnetic THz component is linearly proportional to the in-plane magnetic field (THz waveforms presented in Supplementary Fig. 6).

Discussion

The experimental evidences described above are consistent with the interpretations that the observed THz emission is dominated by transient PTE current. Due to the thickness gradient of the sample and strong absorption to the excitation light, the electron temperature gradient is generated after photoexcitation along both in-plane and out-of-plane directions which are responsible for detectable THz emission without and with magnetic field respectively. According to the above experimental results, the THz emission polarized along y-axis, which is detectable in the E-O sampling configuration, can be given by the following empirical equation:

$${{{{\bf{E}}}}}_{{{{\rm{THz}}}}}\propto {I}^{0.66}\hat{{{\bf{y}}}}+b{I}^{0.5}(\hat{{{\bf{z}}}}\times {{{\bf{B}}}})$$
(1)

where b is a constant that corresponds to the ratio of Seebeck and Nernst contributions. The first term represents the contribution of THz emission generated from a transverse Seebeck effect due to thickness gradient that is independent of the magnetic field (as shown in Fig. 5a); the second term represents the contribution from the Nernst effect that is dependent on the magnetic field. These two terms exhibit different power-law dependences according to Figs. 2a and 3a, and the second term has a linear dependence on B according to Fig. 4d. The direction of the TE current is determined by the thermal gradient of the sample, which does not rely on the polarization of excitation light (as observed in Figs. 2b and 3b).

Fig. 5: Mechanisms of PTE THz emissions and the comparison with semiconductor THz sources.
figure 5

Schematics of (a) transverse Seebeck effect and (b) Nernst effect, where JS is the Seebeck current, JS// is the in-plane component, JS⊥ and JN denote the out-of-plane component and the Nernst current respectively, which are composed by both the electron current and hole current. Comparison of the optimized THz emission from Cd3As2 nano-films with the (c) ZnTe (110) and (d) InSb (111) wafers with thicknesses of 0.5 mm, respectively. Insets in panels (c, d): the corresponding setups.

Full size image

The sublinear excitation intensity dependence and the independence of the excitation polarization help to rule out other nonlinear optical effects which are usually considered as the mechanisms of THz emission. These effects include optical current injection, shift current, optical rectification, depletion field effect, photo-Dember effect and so on15,17,33,34,35,39,41. Although the second-order nonlinear tensor is zero due to the inversion symmetry of bulk Cd3As2, nonzero third-order tensor can exist around an interface or by the assistance of a perpendicular DC electric field. However, those second-order nonlinear optical responses should have linear power dependence instead of sublinear power dependence. The deduction of sublinear power dependence of the TE response, either the transvers Seebeck effect or the Nernst effect, is presented in Supplementary Note 5, by assuming that the intensity of transient currents is dominated by the laser-induced TE effects. The polarization independence is also consistent with the polarization-independent absorption of Cd3As2.

In the next, we analyze the TE response more quantitatively. The transient thermal electric current along y-axis (Jy (B)), which emits THz emission that is detectable when setting [001] of ZnTe crystal along x-axis, is composed of two components:

$${J}_{y}(B)={J}_{S//}(B)+{J}_{N}(B)={\alpha }_{{yy}}(B)\frac{d{T}_{e}}{{dy}}+{\alpha }_{{yz}}(B)\frac{d{T}_{e}}{{dz}}$$
(2)

where JS// and JN are the transverse Seebeck and the Nernst currents respectively, Te is the instantaneous electron temperature, and αyy, αyz are the TE conductivities of Seebeck and Nernst effect. When a magnetic field B is applied along x direction, the TE conductivities can be written as: \({{\alpha }}_{{yy}}={{\alpha }}_{0}\frac{1}{1+{({\mu }B)}^{2}}\) and \({{\alpha }}_{{yz}}={{\alpha }}_{0}\frac{{\mu }B}{1+{({\mu }B)}^{2}}\), as derived from the linearized Boltzmann equation45, where α0 is the zero-field TE conductivity and μ is the carrier mobility. Taking μ ~5000 cm2 V-1 S-1 according to our previous study28, and B ~0.4 T in the measurement, (μB)2 ~ 0.04 ≪ 1. Therefore,

$${J}_{y}(B)\approx {\alpha }_{0}\frac{d{T}_{e}}{{dy}}+\mu B{\alpha }_{0}\frac{d{T}_{e}}{{dz}}$$
(3)

which is fully consistent with empirical Eq. 1 summarized from the experimental observation. According to Eq. 3, the second term (Nernst current JN) should be smaller than the longitudinal Seebeck current JS⊥ = \({{\alpha }}_{0}\frac{d{T}_{e}}{{dz}}\), this is inconsistent with the experimental observation that JN is over an order of magnitude larger than JS⊥ as shown in Fig. 1d and Supplementary Fig. 2d.

This large Nernst signal is due to the ambipolar transport nature during the transient process as illustrated in Fig. 5b. To describe the Nernst response, Eq. 3 would become \({J}_{y}(B)\approx {\alpha }_{e0}\frac{d{T}_{e}}{{dy}}-{\alpha }_{h0}\frac{d{T}_{h}}{{dy}}+{\mu }_{e}B{\alpha }_{e0}\frac{d{T}_{e}}{{dz}}+{\mu }_{h}B{\alpha }_{h0}\frac{d{T}_{h}}{{dz}}\), where e and h denote the electron and hole respectively. In the right side of this equation, the first two terms are contributions from electron and hole to JS// respectively, and the last two terms are contributions to JN. These two types of carriers have the same longitudinal thermal diffusion direction, counterbalancing each other’s contribution to the Seebeck current. However, distinct from the Seebeck effect, the electrons and holes are deflected in opposite directions under magnetic field, and thus, their contributions to the Nernst current are added together as shown in Fig. 5b. The lifetime scale of photoexcited electrons and holes are consistent with recent experimental observation that the photoexcited electrons and holes have an approximate lifetime of 3 ps before recombination in a Tr-ARPES measurement46. The picosecond photoexcited carrier lifetime is sufficient to support the transient ambipolar transport on picosecond timescale and the subsequently enhanced THz emission. Usually, it is not easy to observe the ambipolar transport in low-temperature electrical transport measurement, because one type of carriers predominated in heavily doped semimetals. Previously, the Nernst signal enhancement had been observed in compensated semimetals such as bismuth and graphite hosting almost equal concentrations of electrons and holes under transport measurement without photoexcitation45,47,48. The transient PTE response and suitable photoexcited carrier lifetime of Cd3As2 provide an ideal experimental platform to observe the enhancement on Nernst effect as demonstrated in this work.

Additionally, the anomalous Nernst effect may also contribute to the large enhancement of the PTE signal under magnetic field49,50. The anomalous Nernst effect arises from the nonzero Berry curvature associated with the Weyl nodes without the need to apply a real magnetic field. However, it should be absent in Dirac semimetal without applying magnetic field due to restoration of time-reversal symmetry. When a magnetic field in applied, the magnetic field can split a Dirac node into a pair of Weyl nodes with opposite chirality, and then a nonzero Berry curvature is generated. It has been found that the Berry curvature can induce an anomalous Nernst effect to add on the conventional one6,49,50,51. At low magnetic field and low temperature, it is predicted theoretically that the total Nernst effect is characterized by an almost steplike profile at the low magnetic field and the response is dominated by the anomalous Nernst effect49,50. Although in this work, we cannot clearly sort out the contribution from conventional and anomalous Nernst effects, the large enhancement factor implies very possible contribution from anomalous Nernst response, which calls for further studies.

At last, we compared the THz emission from Cd3As2 by photo-Nernst effect at the weak magnetic field with that from the typical semiconductor THz sources. Figure 5c compares the THz intensity from the response with ZnTe, which is the most widely used THz crystal through optical rectification52. The intensity of THz emission from Cd3As2 at 0.4 T can reach ~70% of that from 0.5 mm-thick ZnTe optimized under the same excitation condition. Figure 5d further compare the THz emission from Cd3As2 with that from a 0.5 mm-thick InSb (111), which is one of the best narrow bandgap semiconductors for THz source with the dominated mechanism of photo-Dember effect34,35,53. The intensity of THz emission from Cd3As2 at 0.14 T can reach ~30% of that from InSb under a reflective detection geometry with 2-μm excitation, confirming the thermoelectric THz emission can also be triggered by longer wavelength due to the zero-bandgap nature of Cd3As2. The comparison measurements indicate that, based on the photoexcited transient Nernst effect, the THz emission efficiency from Cd3As2 with the thickness of tens nanometer is already comparable with typical semiconductor THz sources with submillimeter thickness. When integrated with ferromagnetic materials with careful design54,55,56,57,58, it would provide a highly efficient compact THz source solution based on photo-Nernst effect.

In summary, we have studied THz emissions generated from transient PTE current in Cd3As2 by femtosecond laser excitation under zero and weak magnetic field at room temperature. The THz emissions are generated from a thickness-gradient-assisted transverse Seebeck and Nernst effect, respectively. The magnitude of photo-Nernst current is an order of magnitude larger than Seebeck current at low magnetic field indicating ambipolar transport nature after the photoexcitation. Because of the Nernst effect, THz emission efficiency is comparable to ZnTe and InSb crystals. The THz emission of Cd3As2 can be triggered by ultrafast laser with wavelength over broad-spectrum range due to its semimetallic nature. In addition, the THz emission can be versatility controlled through optical, electric, magnetic and thermal approaches. Compared with its two-dimensional Dirac semimetal counterpart, graphene, which has shown strong and fast PTE responses10,25,59, Cd3As2 has larger steady-state TE coefficients and smaller thermal conductivity49,60,61,62,63,64, which are both more favorable for TE conversion than graphene and thus promises better PTE response at transient state. We expect the transient PTE dynamics revealed in this work provides indispensable device physics for high speed/field electronic and optoelectronic device applications based on three-dimensional topological semimetals.

Methods

Sample preparation

Cd3As2 thin films with about 50-nm thickness were grown in a PerkinElmer (Waltham, MA) 425B molecular beam epitaxy system. Fresh cleaved 2-inch mica was used as the substrate. The data presented in the main text are obtained from two samples with substrate thickness of 30 µm and 70 µm, respectively. The absorbed molecules on the surface of the substrate were removed by annealing at 300 °C for 30 min. Then a CdTe buffer layer with a thickness of approximately 10 nm was deposited on the substrate to assist Cd3As2 nucleation. After that, the Cd3As2 layer was grown on the buffer layer by evaporating the bulk material (99.9999%, American Elements Inc., Los Angeles, CA) at 170 °C. The growth was in situ monitored by the reflection high-energy electron diffraction system. The Cd3As2 layer has high quality and excellent electric properties, which have been systematically studied in previous works28,29,30. The relative thickness mapping of Cd3As2 film is determined by two-dimensional transmittance scanning with a 633-nm cw laser. The laser spot size is about 1 mm, and the scanning step is 0.5 mm.

Transient THz emission setup

In a default experimental configuration, 150-fs laser pulses with a central wavelength of 800 nm from a 250-kHz Ti-sapphire amplifier (RegA) system65 were used to excite the samples and to probe the emitted THz pluses through a 0.5 mm-thick ZnTe (110) crystal using a standard E-O sampling technique66. The excitation beam propagates along the z-direction. It passes through a 3-mm-diameter aperture to shape a circular excitation beam, and the full width at half maximum is about 1.6 mm. The excitation power is measured just in front of the sample and chopped with a frequency about 1 kHz. A half waveplate is mounted into the setup to adjust the polarization direction of the excitation beam before incident onto a sample. The incident angle could be adjusted by tilting the sample in xz-plane, and it is defined as 0° for normal incidence on the Cd3As2 side. The THz pulse was emitted along the z-direction, and then it was collected and focused onto the ZnTe crystal using a pair of parabolic mirrors, which is collimated with a y-polarized sampling beam. Because the E-O sampling is THz polarization dependent67, in general measurements, the ZnTe [001] axis is set along x-direction to detect y-polarized THz component. To acquire the x-component, ZnTe [001] axis could be set along y-directions. In magneto measurements, a pair of permanent magnets was applied along x-direction. The intensity of the magnetic field could be controlled by change the spacing between the two magnets. The maximum magnetic field can reach 0.4 T in this setup. In out-of-plane magnetic field measurement, the magnetic field is applied with an approximately 40°-angle with respect to the x-axis in xz-plane, and the intensity is lowered to 0.1 T to get a larger space between the permanent magnets to permit the unhindered passing through of the THz wave. In comparative measurement with InAs (111), a reflective setup was used and the excitation beam with a wavelength of 2 μm was from an optical parametric amplification. The light is incidence on the 45°-tilted sample along x-direction. the magnetic field can only reach 0.14 T and the direction of the magnetic field has to set along y-direction as limited by the reflection geometry. All measurements were carried out in atmosphere condition at room temperature.