-
PDF
- Split View
-
Views
-
Annotate
-
Cite
Cite
Zheng Liu, Ye-Xin Wang, Yu-Hui Fang, Si-Xue Qin, Zhe-Ming Wang, Shang-Da Jiang, Song Gao, Electric field manipulation enhanced by strong spin-orbit coupling: promoting rare-earth ions as qubits, National Science Review, Volume 7, Issue 10, October 2020, Pages 1557–1563, https://doi.org/10.1093/nsr/nwaa148
- Share Icon Share
Abstract
Quantum information processing based on magnetic ions has potential for applications as the ions can be modified in their electronic properties and assembled by a variety of chemical methods. For these systems to achieve individual spin addressability and high energy efficiency, we exploited the electric field as a tool to manipulate the quantum behaviours of the rare-earth ion which has strong spin-orbit coupling. A Ce:YAG single crystal was employed with considerations to the dynamics and the symmetry requirements. The Stark effect of the Ce3+ ion was observed and measured. When demonstrated as a quantum phase gate, the electric field manipulation exhibited high efficiency which allowed up to 57 π/2 operations before decoherence with optimized field direction. It was also utilized to carry out quantum bang-bang control, as a method of dynamic decoupling, and the refined Deutsch-Jozsa algorithm. Our experiments highlighted rare-earth ions as potentially applicable qubits because they offer enhanced spin-electric coupling which enables high-efficiency quantum manipulation.
INTRODUCTION
Quantum computation offers accelerated ways of solving problems such as database searching [1] and prime factor decomposition [2]. Recently, there is an uprising trend to employ magnetic molecules as quantum bits (qubits), because of their advantages including monodispersed size and chemically controllable properties [3]. In these magnetic molecules where the electron spin is utilized, the electronic structure and quantum coherence are tunable by modifying the environment of the spin [4], and the dimension of the Hilbert space can be increased by designing the magnetic coupling between the spin carriers. The quantum phase memory time of the electron spin has been largely increased to nearly a millisecond by molecular design using nuclear spin-free ligands [5], indicating a bright future for quantum information processing based on electron spin qubits.
However, the magnetic field is hardly desirable when controlling electron spins in the quantum information processing context as the addressability of a single spin is difficult to realize because of the poor locality of the magnetic field. The electric field (E field), which can be much more concentrated and efficient, offers a promising alternative. Coupling between the E field and spins is referred to as spin-electric coupling and observed in the Stark effect. Related magnetic resonance research was first published in the 1960s [6,7], and later reviewed by Mims [8].
Electric control of spins has been proposed to realize electric dipole-induced spin resonance in semiconductors [9], and numerous experiments have been performed to control magnetic quantum dots [10], single electron spins [11], diamond defects [12] and the single nuclear spin in a molecular magnet [13]. Coherent control of an electron spin ensemble has been demonstrated in pnictogen-doped silicon semiconductors [14–16], piezo-diluted magnetic semiconductors [17] and molecular magnets [18]. Most recently, A. Morello and co-workers reported nuclear electric resonance using a high-spin 123Sb nuclear doped in silicon, indicating the importance of electric control of the spin in quantum information processing [19].
The electric control of electron spin was not observed for some of the molecular magnets reported as qubits, and was observed to be inefficient in others, especially those with light elements. This is mainly because the spin-electric coupling is weak, and hinders their further application as qubits. The present study aims to solve this problem using rare-earth ions in a crystal. We report largely enhanced spin-electric coupling with which coherent control by the E field is illustrated as a phase gate. Using this electric phase gate, quantum bang-bang decoupling is realized with microwave (mw) pulses, and the refined Deutsch-Jozsa (D-J) algorithm is demonstrated.
Background: the role of spin-orbit coupling
where |$\overline{\overline {\boldsymbol{g}}} = {g_e}\;\overline{\overline {\bf 1}} + 2\lambda \overline{\overline {\boldsymbol \Lambda} } ,\;\overline{\overline {\boldsymbol{D}}} = {\lambda ^2}\overline{\overline {\boldsymbol{\Lambda }}} ,{\rm{and}}\;\overline{\overline {\boldsymbol{\Lambda }}} $| is a 3 × 3 matrix containing the orbital angular momentum eigenstates and their energies determined by the crystal field term. The effect of |${\hat{H}_E}$| is then included via |$\overline{\overline {\boldsymbol{\Lambda }}} $| in both the Zeeman and zero-field splitting terms and scales up with |$\lambda$|. Cases with weak crystal field and strong spin-orbit coupling require more complicated theoretical treatments, but a similar basic principle still holds, indicating that strong spin-orbit coupling is desirable in seeking materials for efficientEfield control of a spin system. A recent study concluded that, with similar structures, spin-orbit coupling is not the primary factor leading to decoherence [20]. This ensures the possibility of employing strong spin-orbit coupling systems as qubits while maintaining long quantum coherence time. We therefore propose that rare-earth ions are desirable for strong spin-electric coupling.
The Stark effect has been observed in several bulk systems including pnictogen-doped silicon semiconductors [16] and piezoelectric materials [17], where the spin-orbit coupling of the spin carriers is weak, but the applied E field dramatically redistributes the electron wave function. In contrast, in the present study, the E field is used to manipulate the quantum coherence properties of the spins localized on the rare-earth ions individually.
Spin-orbit coupling may also help via spin-spin coupling. Antisymmetric exchange, also known as the Dzyaloshinskii-Moriya (D-M) interaction, arises from the electron spin anisotropy as a result of large spin-orbit coupling. This kind of exchange results in spin canting and net electric polarization, which provides strong spin-electric coupling. R. Sessoli and co-workers tuned the magnetic interaction between an organic radical and a Mn2+ ion in a single crystal in a chiral space group [21]. By applying the ac E field to the sample continuously, the electron paramagnetic resonance (EPR) spectra showed a shift of the effective g-factor up to −2.5 × 10−6.
RESULTS AND DISCUSSION
Observation of the spin-electric coupling
The Stark effect can be observed with two methods. A straightforward method is to apply the E field to the sample continuously to modify the energy levels of the spin carrier. The slight influence can be observed with a strong dc E field or, via signal modulation, with an ac E field. Another approach, which is more efficient and common for molecular magnets, is to measure via quantum phase evolution with pulsed EPR [18].
Previous researches using pulsed EPR were performed with powder samples. To evaluate the spin-electric coupling in more detail, a single crystal was used in our experiments. The emergence of the Stark effect in magnetic resonance requires that the spin carriers must not be located at any inversion centre. On the other hand, the number of inequivalent parameters needed to model the effect would be smaller if the spin centres were of higher symmetry. Based on these considerations, a yttrium aluminium garnet (YAG) single crystal doped with Ce3+ was employed. YAG crystallizes in the space group |$Ia\bar{3}d$| with cubic unit cells. Ce3+ was doped into the YAG crystal at the Y3+ positions [22]. There are six magnetically inequivalent sites [23] and each Ce3+ is identically coordinated by eight oxygen atoms with a local symmetry of D2. The three C2 axes serve as principal axes for the g-factor. Ce3+ ions appear in inversion pairs as a result of the central symmetry of the crystal (Fig. 1a). To eliminate the electron spin dipolar interactions, the concentration of Ce3+ ions was reduced to less than 0.1%, and the average distance between Ce3+ ions was more than 3 nm. To enhance the E field strength and acquire a significant signal, the Ce:YAG single crystal was cut to 0.5 mm in thickness, with a 2 × 6 mm2 area to fit the sample into the cavity. The crystal was mounted so that the (111) face was perpendicular to the E field and parallel with the B0 field. This geometry was fixed for all our experiments.

The Ce:YAG crystal structure and the device used in our experiments. (a) Locations of the six magnetically inequivalent Ce3+ ions marked Ce-1 to Ce-6. Each Ce3+ is coordinated by eight oxygen atoms. Their principle axes are shown as coloured arrows (red, green and blue for x, y and z, respectively). Inversion centres are located at (0.5, 0.5, 0.5) and (0.5±0.25, 0.5±0.25, 0.5±0.25). The inversion pairs involving Ce-1 to Ce-6 generated by (0.5, 0.5, 0.5) are also shown. (b) Ce:YAG single crystal of size 2 × 6 × 0.5 mm3 is mounted on the electrode, and connected to 0.2 mm gold wires with a 50 Ω match resistor.

Spin echo oscillation upon the application of an E field. (a) The pulse sequence employed in the experiment. The rectangular E field pulse is applied between the mw π/2 and π pulses of the standard Hahn echo pulse sequence. (b) The transient spin echo oscillates with increasing frequencies upon enhancing the E field strength from 0 to 106 V/m. (c) The quantum phase evolution frequency, derived from the Fourier transformation of (b), shows a linear dependence on the E field strength. The error bars lie within the size of the points. (d) The quantum phase evolution frequency is simulated with varying directions of the E and B0 fields. The most efficient frequency for quantum phase control is optimized to be 0.96 MHz when E//z and B0//(0.44, 0.90, 0) with respect to the local Ce3+ coordinates. θ and φ denote the polar and azimuthal angles. The experiment was conducted at 10 K, with the microwave frequency 9.667 GHz.
With strong spin-orbit coupling, the Ce3+ with 4f1 electron configuration presents a J = 5/2 state. As a result of the large crystal field splitting, only the ground doublets can be observed in X-band EPR measurements. For simplicity, an effective spin 1/2 model is used, where the values of the highly anisotropic g-tensor depend on the crystal field parameters. Hence, the spin-electric coupling, enhanced by large spin-orbit coupling and significantly influencing the crystal field parameters, can be detected through the g-factor shift in our experiments.
In this model, the Hamiltonian of the Ce3+ ion in the external static magnetic field is |$\hat{H} = {\mu _{\boldsymbol B}}\ {\boldsymbol B}^{T}\overline{\overline {\boldsymbol g}} \hat{\boldsymbol S}$|. By cw-EPR measurements and simulations, the principal values of the g-tensor are determined to be gxx = 1.85, gyy = 0.90 and gzz = 2.74, which agree with the previous report [23].
where parameters |${T_{xyz}}$|, |${T_{yxz}}$| and |${T_{zxy}}$| can be determined by detecting the phase evolution with different applied E fields, as detailed in SI. In our work, five different amplitudes from 0 to 106 V/m were applied to the Ce:YAG single crystal. It can be clearly seen that the spin echo phase evolution frequency showed a linear response to the external E field, as shown in Fig. 2c. It is worth noting that the rising and falling edges of the E field pulse are about 40 ns. The induced magnetic field is at the mG level, comparable to the static and dynamic inhomogeneity from the magnet providing the external magnetic field (B0 field) in the spectrometer, and thus is neglected.
It is important to quantify the effect to optimize the manipulation speed of this quantum phase gate. The Stark effect shows strong anisotropy when the E field pulse is applied in different directions in the Ce3+ local coordinate system. By fitting the four experimental curves with a calculated formula, the parameters are determined to be |${T_{xyz}} = \ 3.30( 2 )\ \times {10^{ - 8}}\ {\rm{m}}/{\rm{V}}$|, |${T_{yxz}} = \ 8.76( 6 )\ \times {10^{ - 8}}\ {\rm{m}}/{\rm{V}}$| and |${T_{zxy}} = \ 12.1( 1 )\ \times {10^{ - 8}}\ {\rm{m}}/{\rm{V}}$|. This measurement, as well as all below, was conducted at 10 K.
With this set of parameters, we can optimize the efficiency of the E field used to manipulate the spin quantum phase. As the Stark effect is highly anisotropic, we can find the directions of the E field and the B0 field in the Ce3+ local coordinates that give the largest coupling strength. Through simulation, a coupling constant of 1.6 MHz·m/(T·MV) is acquired when the B0 field is along (0.44, 0.90, 0) and the E field is along (0, 0, 1) in the Ce3+ local coordinates as shown in Fig. 2d. With |$E\ = \ 1\ {\rm{MV}}/{\rm{m}}$| and |${B_0} = \ 0.6\ {\rm{T}}$|, the coupling efficiency can be enhanced to 0.96 MHz with proper orientation.
Different from the effect of crystal-field and hyperfine-coupling terms, the energy shift induced by the E field, originating from the Zeeman term, is proportional to the B0 field, as in Eq. 3. Therefore, the manipulation efficiency of this E field pulse can be enhanced by increasing the B0 field. A strongB0 field provides plenty of other advantages, e.g. resolution enhancement, a high-quality initialized state, and a long phase memory time [24]. Conclusively, in the present research, it is safe to work with a relatively weak E field.
Electric field as a phase gate
This efficient quantum phase evolution offers us the possibility of demonstrating several quantum manipulations by the E field. Because of the higher operating rate, we can use the mw pulse to reverse the phase factor while the electron spin is undergoing the evolution in the external E field. In this demonstration, the electron spin is firstly prepared in a superposition state in the xy plane of the Bloch sphere by a mw pulse. Then the E field pulse could drive phase evolution of this superposition state, as if there were an additional magnetic field along the z axis. During this evolution driven by the E field, the system can be ‘kicked’ by short mw π pulses to another state within the xy plane. Based on the pulse sequence used in the former experiment in Stark effect detection, these microwave π pulses were added symmetrically about the π pulse in the Hahn-echo sequence. During the E field pulse, a series of mw π pulses at arbitrarily chosen moments were applied to reverse the evolution of the spin polarization, as shown in Fig. 3a. This kind of operation is similar to the experiments with fullerene qubits [25]. The spin polarization evolves in the xy plane under the E field, while a mw π pulse drives the spin rapidly around the B1 field along a half-circle, as illustrated in Fig. 3b.

The E field phase gate demonstration of quantum bang-bang control and the dynamic decoupling realized by mw pulses. (a) (I) The phase evolution within the xy plane of the Bloch sphere driven by theE field phase gate. (II) to (V), Four trains of mw pulses to kick the spin echo. The arrows indicate the moment at which the mw pulses are applied. The quantum phase evolution reverses when applying the microwave pulses. (b) The Bloch sphere demonstration of this bang-bang control. The bold purple curve represents the spin phase evolution induced by the E field phase gate, and the thin green curves indicate the mw ‘kicking’ operation. The straight red arrow is the B1 field direction. (c) Successive mw pulses, which were continuously implemented during the plateaus, can lock and release the phase evolution at arbitrary positions. The experiment was conducted at 10 K, with the microwave frequency 9.667 GHz.
In the above demonstration, the mw pulses operate much faster than the E field pulses. In the original proposal of bang-bang control, rapid spin flips are employed by short mw pulses to prevent the unwanted phase evolution [26]. As in Fig. 3c, when the system evolves in the xy plane driven by the E field and a series of mw π pulses is applied, one can see that the spin evolution is locked, and the spin is decoupled from the E field during the mw pulses. This is an example of the quantum Zeno effect, in which the quantum phase evolution is supressed by repeated measurements or controllable interactions with the environment [27]. The spin evolution is released at the end of each mw pulse, which can be locked again by the next. It is necessary to mention that, unlike the projective quantum Zeno effect experiment which locks the system in the eigenstate [28], this experiment can lock the system in any state.
Phase gate application: D-J algorithm
This efficient quantum phase gate could also help to demonstrate the refined D-J algorithm. The D-J problem is to determine whether a given oracle function |$f\!\!:{\rm{}}{\{ {0,1} \}^n} \to \{ {0,1} \}$|, which takes an n-digit binary input and produces 0 or 1 as the output, is a balanced function (0 for half of the inputs and 1 for the other half) or a constant function (uniformly 0 or 1 for all inputs). A classical computer requires |${2^{n - 1}} + 1$| times of evaluation to solve this problem for the n-bit input situation, whereas quantum computing can solve it in a single try. This has been achieved with 19F nuclear spins [29], but not yet with electron spins. A refined version of the D-J algorithm proposed by D. Collins and co-workers uses one of the input qubits as the output, allowing the n-bit problem to be solved using n qubits [30]. Two types of quantum phase gates are required in our demonstration of the refined D-J algorithm with |$n\ = \ 1$|. The first, the Hadamard gate regarding the basis set |$|\ {s_z} = \ \pm 1/2 \rangle$|, which can be performed by a π/2 mw pulse, is used to generate a uniform superposition of the two possible inputs and to convert the output to a readable eigenstate in the end. The second, the f-controlled control gates are used to encode the functions for calculation. Here, the algorithm can be demonstrated in this two-level system with the E field quantum phase gate with a field strength of |$0.8\ \times {10^6}$| V/m.
There are only four possible oracle functions for the single input situation in the D-J problem: |${f_1} (x) = 0$|, |${f_2}\! (x) = 1$|,|${f_3}\! (x) = 1 - x$| and |${f_4}\! (x) = x$|. These four functions can be encoded as 4 f-controlled gates realized by:

The refined D-J algorithm (n = 1) using the E field pulses as the f-controlled gates. (a) The pulse sequence of the D-J algorithm demonstration. The initial π/2 pulse is applied to prepare the superposition state of |$\frac{1}{{\sqrt 2 }}( {| {0 \rangle + } |1 \rangle} )$|. The green blocks indicate the four f-controlled gates, as shown in the insets of (b–e). The following π/2 pulse is to convert the evolved superposition state into eigenstates |${|0 \rangle}$| (for constant functions, |${f_1}$| and |${f_2}$|) or |${|1 \rangle}$| (for balanced functions, |${f_3}$| and |${f_4}$|). The final Hahn echo pulse sequence is to read out the eigenstates |${|0 \rangle}$| or |${|1 \rangle}$|. The spin echo results are shown in (b–e), among which (b) and (c) with positive phases correspond to eigenstate |${|0 \rangle}$|, whereas (d) and (e) with negative phases indicate eigenstate |${|1 \rangle}$|. The experiment was conducted at 10 K, with the microwave frequency 9.667 GHz.
CONCLUSION
We have demonstrated that an E field pulse applied to the Ce:YAG can be a highly efficient quantum phase gate. The phase gate is utilized to perform quantum bang-bang control and the refined Deutsch-Jozsa algorithm. The example of Ce3+ with 4f1 configuration showed that rare-earth ions, with the aid of a proper symmetry, can enhance the spin-electric coupling via their strong spin-orbit interaction. The optimized manipulation condition with the X-band frequency and 106 V/m field strength allows over 50 π/2 phase gate operations within the phase memory time of up to 15 μs at 10 K. The application of E field quantum phase gates allows for achieving effective manipulations on the rare-earth qubit.
We would like to highlight the importance of spin-orbit coupling in the research on the E field coupling with the spin centre. The manipulation in our experiments, resulting from the strong spin-orbit coupling of the rare-earth ion, is much more efficient compared to those with other separated electron spin centres. The quantum coherence time, of course, is still an important factor to determine whether any application is possible, but faster manipulations allow us to be less picky about this. The E field in our experiments cannot exceed 106 V/m because of the sample size, setting a limit to the maximal number of manipulations before decoherence. The 1 MHz coupling constant might be enhanced further in other setups such as single spin break-junction devices, where the E field can easily reach 108 V/m with a few volts. Such an E field is focused enough to control the electron spin individually.
Shortly after an earlier version of our manuscript was submitted to arXiv [32], a work by J.-J. Liu and co-workers was published investigating the E field effect on the quantum coherence involving a clock transition in the HoW10 molecular nanomagnet, which provides a platform with even longer phase memory time for operation [33]. With other advantages of rare-earth ions such as ultralong optical coherence times, optical readout [34,35], and the possibility of device fabrication by chemical modifications and assembly, the high-efficiency quantum phase gate achieved with enhanced spin-electric coupling in this research indicates that electron spin qubits based on rare-earth ions make promising candidates for quantum information processing towards application.
FUNDING
This work was supported by the National Basic Research Program of China (2018YFA0306003 and 2017YFA0204903), the National Natural Science Foundation of China (21822301, 11805024 and 11847301) and the Beijing Academy of Quantum Information Sciences (Y18G23).
AUTHOR CONTRIBUTIONS
Z.L. performed the measurements, assisted by Y.-X.W. and Y.-H.F. Z.L. and S.-D.J. processed the data and performed theoretical analysis. S.-X.Q. provided the methods to simulate the spin evolutions. Z.-M.W. performed the single crystal face index. The project was conceived by S.G. and S.-D.J. Z.L. and S.-D.J. designed the experiments and wrote the manuscript, assisted by Y.-X.W. All the authors revised the manuscript.
Conflict of interest statement. None declared.
REFERENCES