Observation of parity-time symmetry breaking in a single-spin system

See allHide authors and affiliations

Science  31 May 2019:
Vol. 364, Issue 6443, pp. 878-880
DOI: 10.1126/science.aaw8205

Breaking symmetry with single spins

The energetics of quantum systems are typically described by Hermitian Hamiltonians. The exploration of non-Hermitian physics in classical parity-time (PT)–symmetric systems has provided fertile theoretical and experimental ground to develop systems exhibiting exotic behavior. Wu et al. now demonstrate that non-Hermitian physics can be found in a solid-state quantum system. They developed a protocol, termed dilation, which transformed a PT-symmetric Hamiltonian into a Hermitian one. This allowed them to investigate PT-symmetric physics with a single nitrogen-vacancy center in diamond. The results provide a starting point for exploiting and understanding the exotic properties of PT-symmetric Hamiltonians in quantum systems.

Science, this issue p. 878


Steering the evolution of single spin systems is crucial for quantum computing and quantum sensing. The dynamics of quantum systems has been theoretically investigated with parity-time–symmetric Hamiltonians exhibiting exotic properties. Although parity-time symmetry has been explored in classical systems, its observation in a single quantum system remains elusive. We developed a method to dilate a general parity-time–symmetric Hamiltonian into a Hermitian one. The quantum state evolutions ranging from regions of unbroken to broken PT symmetry have been observed with a single nitrogen-vacancy center in diamond. Owing to the universality of the dilation method, our result provides a route for further exploiting and understanding the exotic properties of parity-time symmetric Hamiltonian in quantum systems.

In quantum mechanics, the real energies of a system are guaranteed by a fundamental axiom associated with the Hermiticity of physical observables. However, a class of non-Hermitian Hamiltonians satisfying parity-time (PT) symmetry can still exhibit real eigenenergies (1). An alternative formulation of quantum mechanics can be established when the axiom of Hermiticity is replaced by the condition of PT symmetry (2, 3). The physics associated with PT-symmetric Hamiltonian has aroused considerable experimental interest (4, 5). The optical analog of PT-symmetric quantum mechanics was first proposed (6) and then extended to other systems, such as electronics (79), microwaves (10), mechanics (11), acoustics (1214), and optical systems with atomic media (1517). Experimental study on PT-symmetric classical optical systems has stimulated many applications such as unidirectional light transport (18, 19) and single-mode lasers (20, 21).

Experimentally investigating PT-symmetric physics in quantum systems is challenging, because quantum systems are governed by Hermitian Hamiltonians. A possible approach is to realize a PT-symmetric Hamiltonian in an open quantum system. However, it is difficult in general to realize a controllable PT-symmetric Hamiltonian by controlling the environment (22). Some progress has been made with this approach in the system of light-matter quasiparticles (23, 24). A lossy Hamiltonian has been constructed to simulate the quantum dynamics under PT-symmetric Hamiltonians (25). However, the additional dissipation introduced in this protocol is usually detrimental to quantum features such as coherence and entanglement. Instead of engineering PT-symmetric Hamiltonians, nonunitary evolution operators have been mimicked in optics (26, 27). Alternatively, two theoretical approaches have been developed to dilate a PT-symmetric Hamiltonian into a Hermitian Hamiltonian with a higher-dimensional Hilbert space (28, 29). However, these dilation methods are limited to the cases of unbroken PT-symmetric Hamiltonians. Thus, observing the transition from the unbroken to the broken PT symmetry through the exceptional point in a single quantum system, such as a single spin, remains elusive.

For the PT-symmetric Hamiltonian Hs and the quantum state |ψ(t) satisfying the Schrödinger type equation, iddt|ψ(t)=Hs|ψ(t), we introduce a dilated state |Ψ(t)=|ψ(t)|+η(t)|ψ(t)|+, governed by a dilated Hermitian Hamiltonian Hs,a(t), wherein |=(|0i|1)/2 and |+=i(|0+i|1)/2 are the eigenstates of σy of the introduced ancilla and η(t) is an appropriate linear operator. By measuring state |Ψ(t), the evolution of |ψ(t)can be obtained in the subspace where the ancilla is |.

The Hamiltonian, Hs,a(t), can be designed flexibly according to realistic physical systems, as Hs,a(t) is not uniquely determined [see supplementary text SM1 (30)]. For example, Hs,a(t) can be chosen asHs,a(t)=Λ(t)I+Γ(t)σz(1)where Λ(t)={Hs(t)+[iddtη(t)+η(t)Hs(t)]η(t)}M1(t), and Γ(t)=i[Hs(t)η(t)η(t)Hs(t)iddtη(t)]M1(t). The time-dependent operator M(t)=η(t)η(t)+I, where σx, σy, and σz are Pauli operators and I is the identity matrix. We note that this derivation of Hs,a(t) holds for an arbitrary Hamiltonian Hs, i.e., our method can be applied to dilate a general PT-symmetric Hamiltonian (SM1).

We investigate the PT-symmetric Hamiltonian of the formHs=[ir11ir](2)where r is a real number. The eigenvalues of Hs are E=±1r2. In the region |r|<1, the eigenvalues E are real, and the system is in an unbroken-symmetry region. Especially, the Hamiltonian Hs is Hermitian with r=0. When |r|>1, the imaginary part of E appears, and the system is in a broken-symmetry region. The point |r|=1 is known as the exceptional point. By taking η(0)=η0I, with η0 being a real number, the corresponding Hs,a can be written in terms of Pauli operators as followsHs,a(t)=A1(t)σxI+A2(t)Iσz+A3(t)σyσz+A4(t)σzσz(3)where A1(t), A2(t), A3(t), and A4(t) are real-valued functions determined by Hs (SM1 and fig. S3).

A single nitrogen-vacancy (NV) center in diamond (Fig. 1A) is used to demonstrate our proposal. The Hamiltonian of the NV center is HNV=2π(DSz2+ωeSz+QIz2+ωnIz+ASzIz), where Sz (Iz) is the spin operator of the electron (nuclear) spin, D=2.87 GHz is the electronic zero-field splitting, Q=4.95 MHz is the nuclear quadrupolar interaction, and A=2.16 MHz is the hyperfine interaction. A magnetic field is applied along the NV symmetry axis ( [1 1 1] crystal axis), yielding the electron (nuclear) Zeeman frequencies ωe (ωn). A subspace of the total system is utilized to form a two-qubit system, which is spanned by the four energy levels |mS,mI=|0,1,|0,0,|1,0, and |1,1 relabeled with |0e|1n, |0e|0n, |1e|0n, and |1e|1n (Fig. 1B).

Fig. 1 Constructing a PT-symmetric Hamiltonian in a NV center.

(A) Schematic atomic structure of the NV center. (B) Energy-level structure of the NV center in the ground multiplet states. Two microwave pulses (blue arrows) and two radiofrequency pulses (orange arrows) are applied. (C) Quantum circuit of the experiment. The electron (nuclear) spin is taken as the system (ancilla) qubit. X (Y) denotes the nuclear spin rotation around the x (y) axis. The initial state is prepared by rotations Y(θ) and X(π/2). The two-qubit system evolves under the dilation Hamiltonian Hs,a. The populations of the four energy levels are measured after the rotation X(π/2).

To realize the Hs,a(t) in Eq. 3, we apply on the electron spin of the two-spin system two selective microwave pulses with time-dependent amplitudes, frequencies, and phases. The pulses can be described by the following Hamiltonian Hc=2π2Ω1(t)cos[0tω1(τ)dτ+ϕ1(t)]Sx|1nn1| +2π2Ω2(t)cos[0tω2(τ)dτ+ϕ2(t)]Sx|0nn0|(4) where Ω1,2(t), ω1,2(t) and ϕ1,2(t) correspond to the amplitudes, angular frequencies, and phases of the two pulses. The evolution of the two-qubit system follows Hs,a in Eq. 3 (in the interaction picture) when the frequencies of the two microwave pulses ω1(t)=ωMW1+2A4(t) and ω2(t)=ωMW22A4(t), the amplitudes Ω1(t)=Ω2(t)=A12(t)+A32(t)/π and the phases ϕ1(t)=ϕ2(t)=arctan2(A3(t),A1(t)) (see SM2 and fig. S4 for the derivation and numerical results).

The experiment was performed on an optically detected magnetic resonance setup. The static magnetic field was set to 506 G, and the NV center was polarized into the state |0e|1n via laser pulses (31). Then the initial state was prepared to |Ψ=|0e|n+η0|0e|+n by the rotations Y(θ) followed by X(π/2) applied on the nuclear spin (Fig. 1C), where θ=2arctanη0. Two selective MW pulses applied to the electron spin were generated by an arbitrary waveform generator. Finally, the nuclear spin rotation X(π/2) transformed the state |Ψ(t)=|ψ(t)e|n+η(t)|ψ(t)e|+n into |Φ(t)=|ψ(t)e|1n+η(t)|ψ(t)e|0n for measurement. Populations of the electron spin states (denoted as P0 and P1) when nuclear spin is |1n were obtained by a proper renormalization of the experimentally measured data (SM3).

The state evolution under Hs is explored by monitoring P0 in the time range of 0 to 8 μs, with the Hermitian case r=0 (Fig. 2A), PT unbroken case r=0.6 (Fig. 2B), exceptional point r=1.0 (Fig. 2C) and PT broken case r=1.4 (Fig. 2D). The measured dependence of the state evolution under Hs with r (Fig. 2F), shows good agreement with the theoretical predictions (Fig. 2E). The transition from the unbroken to the broken PT symmetry through the exceptional point has been clearly demonstrated in Fig. 2. When |r|<1, Hs has real-valued eigenvalues and determines an oscillatory dynamics (Fig. 2B). When |r|>1, the eigenvalues become imaginary, leading to the breakdown of the oscillation (Fig. 2D). Physically, the dynamics of P0 in the region of unbroken PT symmetry is due to the nonvanished oscillatory amplitudes of microwave pulses. While in the region of broken PT symmetry, the amplitudes of microwave pulses decay with time, and the system approaches to a steady state (fig. S4).

Fig. 2 State evolution under Hs.

(A to D) Experimental dynamics of renormalized population P0 when r=0 (A), r=0.6 (B), r=1.0 (C), and r=1.4 (D). All data were collected with 0.5 million averages. Blue dots with error bars are experimental results, and red lines are the theoretical predictions. (E and F) Plot of theoretical results (E) and experimental results (F) for various r values. The color bar stands for P0.

The PT phase transition is also characterized by eigenvalues of the PT-symmetric Hamiltonians derived from experimental data (Fig. 3). The eigenvalues of the Hs are expressed as E±=±1rexp2, where the parameter rexp is obtained by curve fitting the experimental time evolution of the population P0 to theoretical predictions under PT-symmetric Hamiltonian Hs (SM4). When r<1, the eigenvalues E± remain real as r approaches 1 from 0. At the exceptional point r=1, the eigenvalues E± coalesce to 0. When r>1, the system is in the PT-symmetry broken region, and the real parts of the two eigenvalues are zero and the imaginary parts appear. The experimental results show excellent agreement with the corresponding theoretical predictions.

Fig. 3 Experimental observation of the breaking of the PT symmetry.

(A and B) Real and imaginary part of the eigenvalue, E+, of Hs as the function of parameter r. The blue dots with error bars are experimental results, and red lines are theoretical predictions. The interval of 0<r<1 indicates the PT unbroken region and the interval of r>1 indicates the PT broken region. The exceptional point occurs at r=1.

Our work makes NV center a desirable platform for investigating important non-Hermitian physics, such as new topological invariants (3235), quantum thermodynamics (36), and information criticality (29) in this scenario. This platform can be used to investigate various models of decoherence and dissipation in open quantum systems (SM5 and fig. S6). Furthermore, the feature of the exceptional point can be used to improve the sensitivity of quantum sensing (37, 38). This is expected to enhance the performance of NV-based quantum sensor with a variety of applications.

Supplementary Materials

Materials and Methods

Supplementary Text

Figs. S1 to S6

Table S1

References (3941)

References and Notes

  1. Materials, methods, and supplementary text are available as supplementary materials.
Acknowledgments: We thank J. Gong for the helpful discussion. Funding: This work was supported by the National Key R&D Program of China (grant nos. 2018YFA0306600 and 2016YFB0501603), the NNSFC (grant no. 11761131011), the Chinese Academy of Sciences (grant nos. GJJSTD20170001, QYZDY-SSW-SLH004, and QYZDB-SSW-SLH005), and Anhui Initiative in Quantum Information Technologies (grant no. AHY050000). X.R. thanks the Youth Innovation Promotion Association of Chinese Academy of Sciences for their support. Author contributions: J.D. and X. R. proposed the idea and supervised the experiments. X.R., J.G., and Y.W. designed the experiments. Y.W. and W.L. performed the experiments. X.Y. prepared the sample. J.G. and X.S. elaborated the theoretical framework. Y.W. and W.L. carried out the calculations. All authors analyzed the data, discussed the results, and wrote the manuscript. Competing interests: The authors declare no competing interests. Data and materials availability: All data are available in the manuscript or the supplementary materials.
View Abstract

Stay Connected to Science

Navigate This Article