Abstract
Noise is ubiquitous in real quantum systems, leading to non-Hermitian quantum dynamics, and may affect the fundamental states of matter. Here we report in an experiment a quantum simulation of the two-dimensional non-Hermitian quantum anomalous Hall (QAH) model using the nuclear magnetic resonance processor. Unlike the usual experiments using auxiliary qubits, we develop a stochastic average approach based on the stochastic Schrödinger equation to realize the non-Hermitian dissipative quantum dynamics, which has advantages in saving the quantum simulation sources and simplifying the implementation of quantum gates. We demonstrate the stability of dynamical topology against weak noise and observe two types of dynamical topological transitions driven by strong noise. Moreover, a region where the emergent topology is always robust regardless of the noise strength is observed. Our work shows a feasible quantum simulation approach for dissipative quantum dynamics with stochastic Schrödinger equation and opens a route to investigate non-Hermitian dynamical topological physics.
Similar content being viewed by others
Introduction
As a fundamental notion beyond the celebrated Landau-Ginzburg-Wilson framework1, the topological quantum matter has stimulated extensive studies in recent years, with tremendous progress having been achieved in searching for various types of topological states2,3,4,5,6,7. A most important feature of topological matter is the bulk-surface correspondence2,3,4, which relates the bulk topology to boundary states and provides the foundation of most experimental characterizations and observations of topological quantum phases, such as via transport measurements8,9,10 and angle-resolved photoemission spectroscopy11,12,13.
Despite the fact that topological phases are defined at the ground state at equilibrium, quantum quenches in recent studies provide a nonequilibrium way to investigate topological physics14,15,16,17,18,19,20,21,22,23,24,25,26,27,28. Particularly, as a momentum-space counterpart of the bulk-boundary correspondence, the dynamical bulk-surface correspondence was proposed29,30,31,32,33,34, which relates the bulk topology of an equilibrium phase to a nontrivial dynamical topological phase emerging on certain momentum subspaces called band-inversion surfaces (BISs) when quenching the system across topological transitions. This dynamical topology enables a broadly applicable way to characterize and detect topological phases by quantum dynamics and has triggered many experimental studies in quantum simulations, such as in ultracold atoms35,36, nitrogen-vacancy defects in diamond37,38,39, nuclear magnetic resonance (NMR)40, and superconducting circuits41.
The quench-induced dynamical topological phase has been mainly studied in Hermitian systems, while the system is generally non-Hermitian when coupled to the environment42. Recently, the interplay between non-Hermiticity and topology has attracted considerable attention43,44, with rich phenomena being uncovered, such as the exotic topological phases driven by exceptional points45,46,47,48, the anomalous bulk-boundary correspondence46,49,50,51, and the non-Hermitian skin effect52. Experimental observations of the non-Hermitian topological physics have been reported in classical systems with gain and loss, like the photonic systems53,54, the active mechanical metamaterial55, as well as topolectrical circuits56, and in quantum simulators, like the nitrogen-vacancy center57,58, where the non-Hermitian effects are engineered by coupling to auxiliary qubits.
As an important source of dissipation and non-Hermiticity, the dynamical noise is ubiquitous and inevitable in the real quantum simulations, especially for the quantum quench dynamics, and can be described by the stochastic Schrödinger equation59,60. Without the necessity of applying auxiliary qubits, the quantum simulation using the stochastic Schrödinger equation may enable a direct and more efficient way to explore non-Hermitian dynamical phases, hence facilitating the discovery of non-Hermitian topological physics with minimal quantum simulation sources. In particular, the controllable noise can provide a fundamental scheme to explore non-Hermitian dissipative quantum dynamics, and the noise effects on the quench-induced dynamical topological phase give rise to rich nonequilibrium topological physics61. However, the experimental study is currently lacking.
In this article, we report the experimental observation of quench-induced non-Hermitian dynamical topological states by simulating a noising two-dimensional (2D) quantum anomalous Hall (QAH) model on an NMR quantum simulator. Unlike previous experiments using auxiliary qubits57,58, we achieve with advantages the non-Hermitian quench dynamics via simulating the stochastic Schrödinger equation and by averaged measurements over different noise configurations59,60,61. We observe the dynamical topology emerging in the non-Hermitian dissipative quench dynamics on BISs by measuring the time-averaged spin textures in momentum space, and identify two types of dynamical topological transitions classified by distinct dynamical exceptional points by varying the noise strength. Moreover, the existence of a sweet spot region with the emergent topology being robust under arbitrarily strong noise is experimentally verified. Our experiment demonstrates a feasible technique in simulating dynamical topological physics with minimal sources.
Results
Non-Hermitian QAH model
We consider the non-Hermitian 2D QAH model with the magnetic dynamical white noise described by the Hamiltonian
where \({{{{\mathcal{H}}}}}_{{{{\rm{QAH}}}}}({{{\bf{k}}}})={{{\bf{h}}}}({{{\bf{k}}}})\cdot {{{\boldsymbol{\sigma }}}}\) describes the QAH phase62,63, with Bloch vector \({{{\bf{h}}}}=({\xi }_{\rm{so}}\sin {k}_{x},{\xi }_{\rm{so}}\sin {k}_{y},{m}_{z}-{\xi }_{0}\cos {k}_{x}-{\xi }_{0}\cos {k}_{y})\). Here ξ0 (or ξso) simulates the spin-conserved (spin-flipped) hopping coefficient, and mz is the magnetic field. The white noise wi(k, t) of strength \(\sqrt{{w}_{i}}\) couples to the Pauli matrix σi and satisfies \({\langle \langle {w}_{i}({{{\bf{k}}}},t)\rangle \rangle }_{{{{\rm{noise}}}}}=0\) and \({\langle \langle {w}_{i}({{{\bf{k}}}},t){w}_{j}({{{\bf{k}}}},t^{\prime} )\rangle \rangle }_{{{{\rm{noise}}}}}={w}_{i}{\delta }_{ij}\delta (t-t^{\prime} )\), where 〈〈⋅〉〉noise is the stochastic average over different noise configurations. Without noise, the Hamiltonian \({{{{\mathcal{H}}}}}_{{{{\rm{QAH}}}}}\) hosts nontrivial QAH phase for 0 < ∣mz∣ < ∣ξ0∣ with Chern number C1 = sgn(mz), and the phase is trivial for ∣mz∣ > ∣ξ0∣ or mz = 062. The random noise can change the topology of the QAH model, and plays a vital role on the quantum dynamics induced in the present system. We start with the simple situation with a single noise configuration. In this case the quantum dynamics governed by the stochastic Schrödinger equation \({{{\rm{i}}}}{\partial }_{t}\left|\psi ({{{\bf{k}}}},t)\right\rangle ={{{\mathcal{H}}}}({{{\bf{k}}}},t)\left|\psi ({{{\bf{k}}}},t)\right\rangle\) describes a random unitary evolution, which can be further converted into the so-called Itô form59,60 in simulation (see “Methods” for details)
Here \({{{{\mathcal{H}}}}}_{{{{\rm{eff}}}}}={{{{\mathcal{H}}}}}_{{{{\rm{QAH}}}}}-({{{\rm{i}}}}/2){\sum }_{i}{w}_{i}\) is the effective non-Hermitian Hamiltonian, such that the increment of a Wiener process \({W}_{i}({{{\bf{k}}}},t)\equiv (1/\sqrt{{w}_{i}})\int\nolimits_{0}^{t}{{{\rm{d}}}}s\,{w}_{i}({{{\bf{k}}}},s)\) is independent from the wavefunction function \(\left|\psi (t)\right\rangle\), for which we have the Itô rules dtdWi(t) = 0 and dWi(t)dWj(t) = δijdt, and the corresponding expectation value is zero. The formal solution of the above equation reads \(\left|\psi (t)\right\rangle =U(t)\left|\psi (0)\right\rangle\) with
where \({{{\mathcal{T}}}}\) denotes the time ordering. Note that while the equation (3) describes a random unitary evolution in the regime with single noise configuration, after the noise configuration averaging the non-Hermitian dissipative quantum dynamics emerges and is captured by the master equation
where \(\rho ({{{\bf{k}}}},t)\equiv {\langle \langle \left|\psi ({{{\bf{k}}}},t)\right\rangle \left\langle \psi ({{{\bf{k}}}},t)\right|\rangle \rangle }_{{{{\rm{noise}}}}}\) is the stochastic averaged density matrix; see “Methods” for details. The configuration averaging is a key point for the present quantum simulation of non-Hermitian dynamical topological phases.
Quantum simulation approach
We next develop the quantum simulation approach by introducing the discrete Stochastic Schrödinger equation for the non-Hermitian dissipative quantum dynamics, since the continuous evolution cannot be directly emulated with digital quantum simulators. Specifically, we discretize the continuous time as tn = nτ with small-time step τ, where the integer n ranges from zero to the total number of time steps M. The increment of the Wiener process can be simulated by random numbers \({{\Delta }}{W}_{i}({t}_{n})={N}_{i}({t}_{n})\sqrt{\tau }\) for each noise configuration, and we obtain the discretized stochastic Schrödinger equation
with \(\tilde{{{{\mathcal{H}}}}}({{{\bf{k}}}},{t}_{n})={{{{\mathcal{H}}}}}_{{{{\rm{eff}}}}}({{{\bf{k}}}})+{\sum }_{i}\sqrt{{w}_{i}}{\sigma }_{i}{N}_{i}({{{\bf{k}}}},{t}_{n})/\sqrt{\tau }\). Here Ni(tn) is sampled from the standard normal distribution to match the expectation and variance of dWi, and the wavefunction is normalized in each time step. The corresponding unitary evolution operator from time tn to tn+1 reads
leading to the discrete equation of motion
in the linear order of τ after stochastic average, which describes the desired non-Hermitian quantum dynamics. We shall analyse the quality of this discretization versus time step τ in the experiment. The stochastic average of a physical operator \(\hat{{{{\bf{O}}}}}\) at time tn can now be obtained by
This formalism can be directly simulated in experiment.
The above presents the essential idea for simulating the non-Hermitian systems based on the stochastic Schrödinger equation. This method is fundamentally different from that applied in the previous experiments57,58 using auxiliary qubits, where the non-Hermiticity is obtained from a Hermitian Hamiltonian in the extended Hilbert space by tracing the auxiliary degrees of freedom and careful designs of the quantum circuit with complex unitary operations are required64,65. In contrast, our temporal average approach based on the stochastic Schrödinger equation saves the resources of qubits and avoids the implementation of complex gates, which benefits the experimental platforms in various scenarios. Moreover, this quantum simulation approach can be directly extended to exploring higher dimensional non-Hermitian topological phases and phase transitions.
Non-Hermitian dynamical topological phases
Before presenting the experiment, in this section, we briefly introduce the non-Hermitian dynamical topological phases emerging in the quench dynamics described by Eq. (4) and to be studied in this work.
The system is initially prepared at the fully polarized ground state \({\rho }_{0}=\left|\downarrow \right\rangle \left\langle \downarrow \right|\) of a deep trivial Hamiltonian with ∣mz∣ ≫ ∣ξ0∣. After quenching mz to a nontrivial value at time tn = 0, the system starts to evolve under the post-quench Hamiltonian \({{{\mathcal{H}}}}({{{\bf{k}}}},{t}_{n})\); see Fig. 1a. Without noise, the spin polarization \(\langle {{{\boldsymbol{\sigma }}}}({{{\bf{k}}}},{t}_{n})\rangle \equiv \rm{Tr}[{{{\boldsymbol{\sigma }}}}\mathop{\prod }\nolimits_{i = 0}^{n-1}U({t}_{n-i},{t}_{n-i-1}){\rho }_{0}\mathop{\prod }\nolimits_{i = 0}^{n-1}{U}^{{\dagger} }({t}_{i+1},{t}_{i})]\) precesses with respect to the Hamiltonian vector h; see Fig. 1a. The post-quench QAH phase can be determined by the dynamical topology emerging on BISs29, identified as the momentum subspaces with hz = 0, where the initial state is perpendicular to the SO field hso ≡ (hx, hy), leading to vanishing time-averaged spin polarizations.
In the presence of non-Hermiticity, the precession axis for each noise configuration is distorted, leading to the deformation for the BISs and dissipative effect. To characterize the noise effect, the spin polarization needs to be stochastically averaged as
over different noise configurations (Fig. 1b). Compared to the spin polarization 〈σ(k, tn)〉 without noise, the stochastic averaged s(k, tn) follows the non-Hermitian dynamics and exhibits dephasing and amplitude decaying effects. We compensate the amplitude decay by rescaling s(k, tn), leading to the rescaled spin polarization \(\tilde{{{{\bf{s}}}}}({{{\bf{k}}}},{t}_{n})\equiv {{{{\bf{s}}}}}_{0}({{{\bf{k}}}})+{{{{\bf{s}}}}}_{+}({{{\bf{k}}}}){{{{\rm{e}}}}}^{-{{{\rm{i}}}}\omega ({{{\bf{k}}}}){t}_{n}}+{{{{\bf{s}}}}}_{-}({{{\bf{k}}}}){{{{\rm{e}}}}}^{+{{{\rm{i}}}}\omega ({{{\bf{k}}}}){t}_{n}}\), where the coefficients s0,± and oscillation frequency ω are extracted from the experimental data by fitting; see “Methods”. Similar to the noiseless case, the time average
vanishes on the deformed BISs (dubbed as dBISs)61, with the number of steps M being large enough to minimize the error. The non-Hermitian dynamical topological phase is captured by the dynamical invariant \({{{\mathcal{W}}}}\equiv \frac{1}{2\pi }{\oint }_{{{{\rm{dBIS}}}}}{{{\bf{g}}}}({{{\bf{k}}}}){{{\rm{d}}}}{{{\bf{g}}}}({{{\bf{k}}}})\), which describes the winding of dynamical field \({{{\bf{g}}}}({{{\bf{k}}}})=(1/{{{{\mathcal{N}}}}}_{{{{\bf{k}}}}}){\partial }_{{k}_{\perp }}(\overline{{\tilde{s}}_{x}({{{\bf{k}}}})},\overline{{\tilde{s}}_{y}({{{\bf{k}}}})})\) on the dBISs. Here k⊥ is perpendicular to the dBISs and \({{{{\mathcal{N}}}}}_{{{{\bf{k}}}}}\) is a normalization factor. Under the dynamical noise, the non-Hermitian dynamical topological phases and phase transitions may be induced, as studied in the experiment presented below.
Experimental setup
The demonstration is performed on the NMR quantum simulator. The sample is the 13C-labeled chloroform dissolved in acetone-d6, with 13C and 1H nuclei denoted as two qubits. The 2D QAH model is simulated by the qubit 13C, while the other qubit 1H enhances the signal by Overhauser effect (see Fig. 1c and “Methods”). In the double-rotating frame, the total Hamiltonian of this sample is
where J = 215 Hz is the coupling strength, Bi is the amplitude of the control pulse, and ϕi is the phase. We firstly initialize the system into the fully polarized state \(\left|\downarrow \right\rangle\) using the nuclear Overhauser effect66. Then we quench mz to the nontrivial region with ∣mz∣ < 2ξ0 and allow the system to evolve under the effective Hamiltonian \(\tilde{{{{\mathcal{H}}}}}\), in which the non-Hermitian constant term i∑iwi can be ignored. The evolution is realized by the Trotter approximation combined with control pulse optimizations as follows.
We study the non-Hermitian dissipative quantum dynamics from time t = 0–30 ms. For each noise configuration, numerical results show that the discrete evolution approximates the continuous evolution of the stochastic Schrödinger equation quite well, when the total number of time steps is greater than 100; see Fig. 2. In the experiment, we discretize the time into 300 segments, such that the Hamiltonian in each interval is approximately time-independent. As the interval τ is sufficiently small, the evolution in the n-th step can be realized using the first-order Trotter decomposition:
with \({\eta }_{x,y}={\xi }_{\rm{so}}\sin {k}_{x,y}+\sqrt{{w}_{x,y}}{N}_{x,y}({t}_{n})/\sqrt{\tau }\) and \({\eta }_{z}={m}_{z}-{\xi }_{0}\cos {k}_{x}-{\xi }_{0}\cos {k}_{y}+\sqrt{{w}_{z}}{N}_{z}({t}_{n})/\sqrt{\tau }\). Here ξ0 is set to 1 kHz. Each term on the right-hand side represents a single-qubit rotation with a rotating angle 2ηiτ along axis σi, which can be experimentally realized by tuning the amplitude and phase of the control pulse in Eq. (11) (z-rotation can be indirectly realized via x- and y-rotations), with further pulse optimization techniques to reduce control errors; see Fig. 1c.
We measure the spin polarization 〈σ(k, t)〉 for single noise configuration at every 20τ interval. After averaging over all noise configurations, we obtain the stochastic averaged spin polarization s(k, t), from which the rescaled spin polarization \(\tilde{{{{\bf{s}}}}}({{{\bf{k}}}},t)\) can be constructed by fitting. We repeat the above procedures for the whole momentum space to obtain the time-averaged spin textures \(\overline{\tilde{{{{\bf{s}}}}}({{{\bf{k}}}})}\).
Experimental results
We start from the weak noise regime, where the noise strength is chosen as wx = 0.05ξ0, wy = 0, and wz = 0.01ξ0 with ξso = 0.2ξ0. The system is quenched to the topological phase with mz = 1.2ξ0. In Fig. 3a, we plot the spin polarization 〈σ(t)〉 at the momentum k = (1.286, −0.257) for four different noise configurations. For each noise configuration, no notable decay exists in the spin polarization, manifesting the unitary evolution. However, after averaged over all noise configurations, the system clearly exhibits the non-Hermitian dissipative quantum dynamics; see Fig. 3b.
Figure 3c shows the measured time-averaged spin textures \(\overline{{\tilde{s}}_{i}({{{\bf{k}}}})}\) with fixed ky = −0.257 and kx ∈ [−1.8, 1.8], obtained by rescaling the stochastic averaged spin polarization s(k, t). The momenta with vanishing values represent dBIS points. To obtain the 2D time-averaged spin texture, we discretize the whole momentum space kx, ky ∈ [−1.8, 1.8] into a 15 × 15 lattice and repeat the above measurements. The results are shown in Fig. 3d, from which the dBIS momenta can be identified. Although the corresponding shape is slightly deformed from the ideal BIS with hz = 0 in the absence of noise (see “Methods”), it is obvious that under weak noise, the dynamical field g(k) can be defined everywhere on dBIS and characterizes the nontrivial non-Hermitian dynamical topological phase (Fig. 3e). Indeed, this emergent dynamical topology is robust against the weak noise and is protected by the finite minimal oscillation frequency on the dBISs, serving as a bulk gap for the dynamical topological phase. The experimental minimum oscillation frequency on dBISs is given by \({\omega }_{\min }=0.4175\) kHz, close to the theoretical value 0.4063 kHz (Fig. 3b). Further, this non-Hermitian dynamical topological phase may break down under strong noise, with two types of dynamical transition being observed below.
We now increase the noise strength to a strong regime with wx = 0.1ξ0, wy = 0.05ξ0, and wz = 0.45ξ0. The averaged spin polarization is measured in the same way as in the weak noise regime. However, the quench dynamics are essentially different, where the spin polarization s(t) at certain momenta, for instance kx = −1.286 and ky = −0.257, displays pure decay without oscillation; see Fig. 4a, b. For these momenta, the dynamical field g vanishes. In Fig. 4c, we show the corresponding spin textures. From the result for \(\overline{{\tilde{s}}_{z}}\), we find that singularities emerges on the dBISs and interrupt their continuity. Thus the dBIS breaks down, while the deformation of the shape of dBIS is small, and the non-Hermitian dynamical topological phase transition occurs. In Fig. 4d, we increase the noise strength to wx = 1.6ξ0, wy = 0, wz = 0.8ξ0 and set a strong SO coupling coefficient with ξso = 2ξ0. A qualitatively different dynamical transition is uncovered, where the dBISs are dramatically deformed by the noise and are connected to the topological charge at k = 0. Due to this singularity, the dynamical topology also breaks down. The above two qualitatively different phenomena are referred to as type-I and type-II dynamical transitions, respectively, which we examine below in more detail.
We notice that the equilibrium topological phase transition usually corresponds to the close of energy gap. In the nonequilibrium regime, the analogous quantity is the oscillation frequency. Here we observe the corresponding momentum distribution in Fig. 5a. One can see that the oscillation frequency is in general nonzero but may vanish on certain dBISs momenta when these two types of dynamical transition occur, i.e., \({\omega }_{\min }({{{{\bf{k}}}}}_{c})\to 0\). Indeed, the momenta (kc) with just vanishing oscillation frequency are exceptional points of the Liouvillian superoperator, on which the eigenvectors \({{{{\bf{s}}}}}_{\pm }^{L(R)}\) coalesce61. Thus the dynamical transitions are driven by exceptional points with vanishing oscillation frequency on dBISs. To further distinguish these two types of dynamical transition and the corresponding exceptional points, we treat the Liouvillian superoperator as a three-level system; see “Methods”. The coefficient s+ of rescaled dynamical spin polarization \(\tilde{{{{\bf{s}}}}}({{{\bf{k}}}},t)\) contains the information of corresponding eigenvectors \({{{{\bf{s}}}}}_{\pm }^{L(R)}\). Like the spin-1 system, we measure the Liouvillian polarization \(\langle {L}_{\alpha }\rangle \equiv {{{{\bf{s}}}}}_{+}^{{\dagger} }{L}_{\alpha }{{{{\bf{s}}}}}_{+}\) to characterize the Liouvillian superoperator. Here the operator Lα is defined as
and Lz = i[Ly, Lx], which satisfies [Lα, Lβ] = iϵαβγLγ. The measured momentum distribution of these quantities in the experiment is shown in Fig. 5b, c, from which an important feature of exceptional points is observed that the component 〈Lx〉 ≈ 0 and 〈Ly〉 ≈ 0 vanish on these points while 〈Lz〉 is in general nonzero (e.g., see Fig. 5c). Therefore, the exceptional points are actually the singularities in the two-component vector field (〈Lx〉, 〈Ly〉).
With this observation and to characterize the exceptional points, we consider the Liouvillian polarization on a small loop \({{{\mathcal{S}}}}\) enclosing the exceptional points, as shown in Fig. 6a, b. Although the component 〈Lz〉 is nonzero on this loop, the trajectory projected on the 〈Lx〉-〈Ly〉 plane indeed defines a winding number61
which distinguishes the two types of dynamical transitions. We observe that for type-I transition, the winding number NE = 0 is trivial, while the winding NE = 1 is nontrivial for the type-II dynamical transition. Consequently, these distinct exceptional points on dBISs shows the fundamental difference between the type-I and type-II dynamical transitions. Moreover, regardless of the shape and size of the loop \({{{\mathcal{S}}}}\), the winding number NE only depends on the topological properties of the enclosed exceptional points as long as the loop does not cross any other singular points; see Fig. 6a, b. Here we note that the topological charges are always singularities of the field (〈Lx〉, 〈Ly〉) and have nontrivial winding number61 (Fig. 6c, d). The loop \({{{\mathcal{S}}}}\) should be introduced without enclosing any non-exceptional charge momentum in characterizing the dynamical transitions and corresponding exceptional points. This also tells that the type-II dynamical transition is similar to the equilibrium topological phase transition, in which the topological charges serve as singular points and the transition occurs when they pass through the BISs29,30. On the other hand, the type-I transition is a peculiar feature of the quench-induced non-Hermitian dynamical topological phase transition.
Although the non-Hermitian dynamical topological phase may typically be destroyed in the strong noise regime, a quite interesting feature of the present system is the existence of a sweet spot region satisfying61
in which regime the dynamical topology is always robust at any finite noise strength, as characterized by the taper-type region in Fig. 7. In particular, for the central line with wx = wy = wz, we experimentally increase the noise strength wi in each direction from 0.5ξ0 to a very large value wi ≃ 10ξ0 (points O1,2,3) and measure the corresponding time-averaged dynamical spin textures. We observe that although the noise strength is much large compared with all other energy scales, the dBIS in \(\overline{{\tilde{s}}_{z}}\) remains stable, without suffering singularities. Inside the taper-type region the dynamical topology is well-defined on the dBIS, in sharp contrast to outside points (P1,2). The experimental confirmation of this sweet spot region may offer guidance in designing noise-tolerant topological devices.
Discussion
We have experimentally reported the quantum simulation of non-Hermitian quantum dynamics for a 2D QAH model coupled to dynamical noise based on a stochastic average approach of the stochastic Schrödinger equation, and simulated non-Hermitian dynamical topological phases and phase transitions. Our method does not require the ancillary qubits and careful designs of complex unitary gates, hence saving the simulation sources and avoiding the implementation of complex gates in the experiment. The dynamical topological physics driven by dynamical noise has been observed, including the stability of non-Hermitian dynamical topological states protected by the minimal oscillation frequency of quench dynamics under weak noise and two basic types of dynamical topological transitions driven by strong noise and classified by distinct exceptional points. Moreover, a sweet spot region is observed, where the non-Hermitian dynamical topological phase survives at arbitrarily strong noise.
Our experiment has shown an advantageous quantum simulation approach to explore the non-Hermitian dynamical topological physics, in which only a minimal number of qubits are used. This approach is directly applicable to high dimensions by taking into account more, but still minimal number of qubits, in which the rich phenomena are expected, and also to other digital quantum simulators.
Methods
Stratonovich stochastic Schrödinger equation
We consider the non-Hermitian 2D QAH model (1) with the magnetic dynamical white noise wi(k, t). Since the dynamical white noise is in some sense infinite, the dynamical equation \({\partial }_{t}\left|\psi ({{{\bf{k}}}},t)\right\rangle =-{{{\rm{i}}}}{{{\mathcal{H}}}}({{{\bf{k}}}},t)\left|\psi ({{{\bf{k}}}},t)\right\rangle\) cannot be considered as an ordinary differential equation. Instead, it should be regarded as an integral equation
where \({W}_{i}({{{\bf{k}}}},t)=(1/\sqrt{{w}_{i}})\int\nolimits_{0}^{t}{{{\rm{d}}}}s\,{w}_{i}({{{\bf{k}}}},s)\) is a Wiener process. For brevity, the symbols of integration are usually dropped, leading to the stochastic Schrödinger equation
In general, there are two definitions of stochastic integration, i.e., the Stratonovich form
and the Itô form
The basic difference is that the integrand f(t) and the increment dW(t) are independent of each other in the Itô form, namely 〈〈f(t)dW(t)〉〉noise = f(t)〈〈dW(t)〉〉noise = 0, while they are not independent in the Stratonovich form. The Schrödinger equation (17) must be interpreted as a Stratonovich stochastic differential equation59,60, such that the quantum mechanical probability is preserved, i.e., d〈ψ(t)∣ψ(t)〉 = 0.
Converting into the Itô form
Since the wavefunction \(\left|\psi (t)\right\rangle\) and the increment dWi(t) are not independent in the Stratonovich form, it is usually convenient to convert the Stratonovich stochastic Schrödinger equation (17) into the Itô form, which takes the form
Due to \((1/2)(\left|\psi (t+{{{\rm{d}}}}t)\right\rangle +\left|\psi (t)\right\rangle )=[{{{\bf{1}}}}-({{{\rm{i}}}}/2)({{{{\mathcal{H}}}}}_{{{{\rm{eff}}}}}{{{\rm{d}}}}t+{\sum }_{i}{\alpha }_{i}{{{\rm{d}}}}{W}_{i}(t))]\left|\psi (t)\right\rangle\), we have the following relation between the Stratonovich integral and the Itô integral
where we have used the Itô rules dtdWi(t) = 0 and dWi(t)dWj(t) = δijdt for the increment of a Wiener process. Substituting this into the Itô stochastic Schrödinger equation (20), we obtain
Compared with the original Stratonovich stochastic Schrödinger equation (17), it is easy to find
In the main text, we have shown that the formal solution of the Itô stochastic Schrödinger equation (20) is given by a unitary evolution U(t) [see Eq. (3)]. To prove that U(t) is indeed the solution of the Itô equation, we shall note that
where the terms other than dt and dWidWi = dt vanish according to the Itô rules. Thus we recover the Itô stochastic Schrödinger equation, i.e., \({{{\rm{d}}}}U(t)=U(t+{{{\rm{d}}}}t)-U(t)=-{{{\rm{i}}}}[{{{{\mathcal{H}}}}}_{{{{\rm{eff}}}}}{{{\rm{d}}}}t+{\sum }_{i}\sqrt{{w}_{i}}{\sigma }_{i}{{{\rm{d}}}}{W}_{i}(t)]U(t)\).
Non-Hermitian dissipative quantum dynamics
We now consider the equation of motion for the stochastic density operator \(\varrho (t)=\left|\psi (t)\right\rangle \left\langle \psi (t)\right|\), namely
Since the increments dWi(t) are independent of ϱ(t) in the Itô form, after average over different noise configurations the last term vanishes and we arrive at the Lindblad master equation (4) for the stochastic averaged density matrix ρ(t) ≡ 〈〈ϱ(t)〉〉noise, which describes the non-Hermitian dissipative quantum dynamics.
Stochastically averaged spin dynamics
In this section, we show the stochastically averaged spin dynamics. According to the master equation (4), the stochastically averaged spin polarization s(k, t) is governed by the equation of motion
with the Liouvillian superoperator
The solution to this dissipative quantum dynamics can be written as
with the coefficients \({{{{\bf{s}}}}}_{\alpha }({{{\bf{k}}}})=[{{{{\bf{s}}}}}_{\alpha }^{L}({{{\bf{k}}}})\cdot {{{\bf{s}}}}({{{\bf{k}}}},0)]{{{{\bf{s}}}}}_{\alpha }^{R}\) for α = 0, ±. Here \({{{{\bf{s}}}}}_{\alpha }^{L(R)}\) satisfying \({{{{\bf{s}}}}}_{\alpha }^{L}({{{\bf{k}}}})\cdot {{{{\bf{s}}}}}_{\beta }^{R}({{{\bf{k}}}})={\delta }_{\alpha \beta }\) are the left (right) eigenvectors of the Liouvillian superoperator \({{{{\mathcal{L}}}}}^{T}{{{{\bf{s}}}}}_{\alpha }^{L}=-{\lambda }_{\alpha }{{{{\bf{s}}}}}_{\alpha }^{L}\), \({{{\mathcal{L}}}}{{{{\bf{s}}}}}_{\alpha }^{R}=-{\lambda }_{\alpha }{{{{\bf{s}}}}}_{\alpha }^{R}\) with eigenvalues λ0 and λ± = λ1 ± iω, respectively. The oscillation frequency is denoted as ω.
In experiments, the coefficients sα, decay rates λ0,1, and oscillation frequency can be extracted by fitting the experimental data. By ignoring λ0,1, we obtain the rescaled spin polarization \(\tilde{{{{\bf{s}}}}}({{{\bf{k}}}},t)\).
NMR sample
The experiment is performed on the nuclear magnetic resonance processor (NMR). The sample we used is the 13C-labeled chloroform dissolved in rmacetone—d6. The 13C spin is used as the working qubit, which is controlled by radio-frequency (RF) fields. The 1H is decoupled throughout the experiment by Overhauser effect which can enhance the signal strength of 13C.
Overhauser effect
Applying a weak RF field at the Larmor frequency of one nuclear spin for a sufficient duration may enhance the longitudinal magnetization of the others, this is the steady-state nuclear Overhauser effect (NOE). In modern NMR, the steady-state NOE is mainly exploited in heteronuclear spin systems, where the enhancement of magnetization is useful and dramatic.
For an ensemble of heteronuclear systems made up with a nuclei I with gyromagnetic ratio γI and a nuclei S with gyromagnetic ratio γS, with ∣γI∣ > ∣γS∣, the thermal equilibrium state of the heteronuclear system is
where βI/βS = γI/γS, \(\frac{1}{4}{\hat{{{{\rm{I}}}}}}_{z}=\hat{{\sigma }_{z}}\otimes \hat{1}\), \(\frac{1}{4}{\hat{{{{\rm{S}}}}}}_{z}=\hat{1}\otimes \hat{{\sigma }_{z}}\). Assume that a continuous RF field is applied at the I-spin Larmor frequency, inducing transitions across two pairs of energy levels. After sufficient time, the RF field equalizes the populations across the irradiated transitions. At that time, the populations settle into steady-state values, which do not change any more, as long as the RF field is left on. The steady-state spin density operator is
By comparing with thermal equilibrium Eq. (29), the S-spin magnetization is enhanced by factor ϵNOE. For our experiment I = 1H and S = 13C.
Noise configurations
For the stochastic average, it is clear that the more noise configurations are considered, the more reasonable result we obtain, as shown in Fig. 8. On the other hand, the large number of noise configurations takes a lot of time. We have performed numerical simulations, and found that the average of 5000 noise configurations can precisely approximate the non-Hermitian dissipative quantum dynamics; see Fig. 8d. However, in NMR experiments, as the relaxation time is in the magnitude of seconds, a complete implementation of all 5000 noise configurations requires an extremely long-running time that we cannot afford.
An alternative method to solve the issue is to reduce the number of noise configurations by numerical simulation prior to the implementation of experiments. We test different number of noise configurations, and plot their average dynamics in comparison with the ideal dynamics of the non-Hermitian Hamiltonian; see Fig. 8a-d. The simulated results show that with the increase of the number of noise configurations, the stochastic averaged spin polarization 〈〈〈σ(k, t)〉〉〉noise would eventually approach to the spin polarization s(k, t) solved by the Lindblad master equation61. The opposite is that with the decrease of the number of noise configurations, the performance of the approximation becomes more fluctuating (Fig. 8e). But the 〈〈〈σ(k, t)〉〉〉noise always fluctuates above and below the theoretical spin polarization s(k, t). After a sufficient number of averaging, the stochastic averaged spin polarization that in the opposite side of theoretical value will be offset by each other. We randomly generated 5000 noise configurations N(tn) that satisfy the normal distribution and separate these noise configurations into two subgroups in which the noise has an opposite effect on 〈〈〈σ(k, t)〉〉〉noise. Then we use numerical simulations to select two noise configurations from these two subgroups respectively such that the 〈〈〈σ(k, t)〉〉〉noise obtained from these four noise configurations precisely approximate the one obtained from the 5000 configurations (Fig. 8f). From experimental results and the corresponding fidelities, it can be concluded that the experiment is in excellent accordance with the simulations. And the theory and experiment results of each group of noise are in good agreement (Fig. 8g) So, it is somehow reasonable to utilize four noise configurations to replace a full description of the non-Hermitian dynamics under 5000 noise configurations. We would like to emphasize that the above numerical simulations to reduce the number of noise configurations does not affect the applicability of the method. In many other quantum systems such as the superconducting circuits or nitrogen-vacancy centers in diamond, the implementation of experiments takes a much shorter time, so they can realize the stochastic average with a larger number of noise configurations.
Experimental results vs. theoretical results
In this section, we show the agreement of our experimental results with the theoretical calculations. In Fig. 9, we compare the experimental spin textures with theoretical ones. Although the resolution of experimental data is lower than that of numerical calculations, the experimental results and the theoretical simulations reach the same conclusion. In Fig. 10, we show the numerical calculations for exceptional points and the corresponding winding numbers, which are consistent with our experimental results (Fig. 7).
Data availability
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
References
Landau, L. & Lifshitz, E. Statistical Physics, Course Theoretical Physics, vol. 5 (Pergamon Press 1999).
Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Rev. Mod. Phys. 82, 3045 (2010).
Qi, X.-L. & Zhang, S.-C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 1057 (2011).
Chiu, C.-K., Teo, J. C., Schnyder, A. P. & Ryu, S. Classification of topological quantum matter with symmetries. Rev. Mod. Phys. 88, 035005 (2016).
Ando, Y. & Fu, L. Topological crystalline insulators and topological superconductors: from concepts to materials. Annu. Rev. Condens. Matter Phys. 6, 361–381 (2015).
Sato, M. & Ando, Y. Topological superconductors: a review. Rep. Prog. Phys. 80, 076501 (2017).
Xie, B. et al. Higher-order band topology. Nat. Rev. Phys. 3, 520–532 (2021).
Konig, M. et al. Quantum spin hall insulator state in hgte quantum wells. Science 318, 766–770 (2007).
Chang, C.-Z. et al. Experimental observation of the quantum anomalous hall effect in a magnetic topological insulator. Science 340, 167–170 (2013).
He, Q. L. et al. Chiral majorana fermion modes in a quantum anomalous hall insulator–superconductor structure. Science 357, 294–299 (2017).
Hsieh, D. et al. A topological dirac insulator in a quantum spin hall phase. Nature 452, 970–974 (2008).
Chen, Y. et al. Experimental realization of a three-dimensional topological insulator, bi2te3. Science 325, 178–181 (2009).
Xia, Y. et al. Observation of a large-gap topological-insulator class with a single dirac cone on the surface. Nat. Phys. 5, 398–402 (2009).
Vajna, S. & Dóra, B. Topological classification of dynamical phase transitions. Phys. Rev. B 91, 155127 (2015).
Caio, M., Cooper, N. R. & Bhaseen, M. Quantum quenches in chern insulators. Phys. Rev. Lett. 115, 236403 (2015).
Budich, J. C. & Heyl, M. Dynamical topological order parameters far from equilibrium. Phys. Rev. B 93, 085416 (2016).
Wilson, J. H., Song, J. C. & Refael, G. Remnant geometric hall response in a quantum quench. Phys. Rev. Lett. 117, 235302 (2016).
Gong, Z. & Ueda, M. Topological entanglement-spectrum crossing in quench dynamics. Phys. Rev. Lett. 121, 250601 (2018).
Fläschner, N. et al. Observation of dynamical vortices after quenches in a system with topology. Nat. Phys. 14, 265–268 (2018).
Song, B. et al. Observation of symmetry-protected topological band with ultracold fermions. Sci. Adv. 4, eaao4748 (2018).
Qiu, X., Deng, T.-S., Hu, Y., Xue, P. & Yi, W. Fixed points and dynamic topological phenomena in a parity-time-symmetric quantum quench. iScience 20, 392–401 (2019).
Huang, B. & Liu, W. V. Floquet higher-order topological insulators with anomalous dynamical polarization. Phys. Rev. Lett. 124, 216601 (2020).
Ünal, F. N., Bouhon, A. & Slager, R.-J. Topological euler class as a dynamical observable in optical lattices. Phys. Rev. Lett. 125, 053601 (2020).
Hu, H. & Zhao, E. Topological invariants for quantum quench dynamics from unitary evolution. Phys. Rev. Lett. 124, 160402 (2020).
Hu, H., Yang, C. & Zhao, E. Quench dynamics of hopf insulators. Phys. Rev. B 101, 155131 (2020).
Lu, Y.-H., Wang, B.-Z. & Liu, X.-J. Ideal weyl semimetal with 3d spin-orbit coupled ultracold quantum gas. Sci. Bull. 65, 2080–2085 (2020).
Wang, Z.-Y. et al. Realization of an ideal weyl semimetal band in a quantum gas with 3d spin-orbit coupling. Science 372, 271–276 (2021).
Mizoguchi, T., Kuno, Y. & Hatsugai, Y. Detecting bulk topology of quadrupolar phase from quench dynamics. Phys. Rev. Lett. 126, 016802 (2021).
Zhang, L., Zhang, L., Niu, S. & Liu, X.-J. Dynamical classification of topological quantum phases. Sci. Bull. 63, 1385–1391 (2018).
Zhang, L., Zhang, L. & Liu, X.-J. Dynamical detection of topological charges. Phys. Rev. A 99, 053606 (2019).
Zhang, L., Zhang, L. & Liu, X.-J. Unified theory to characterize floquet topological phases by quench dynamics. Phys. Rev. Lett. 125, 183001 (2020).
Yu, X.-L. et al. Quantum dynamical characterization and simulation of topological phases with high-order band inversion surfaces. PRX Quantum 2, 020320 (2021).
Li, L., Zhu, W. & Gong, J. Direct dynamical characterization of higher-order topological phases with nested band inversion surfaces. Sci. Bull. 66, 1502–1510 (2021).
Zhang, L., Jia, W. & Liu, X.-J. Universal topological quench dynamics for \({{\mathbb{Z}}}_{2}\) topological phases. Sci. Bull. https://www.sciencedirect.com/science/article/abs/pii/S2095927322001554 (2022).
Sun, W. et al. Uncover topology by quantum quench dynamics. Phys. Rev. Lett. 121, 250403 (2018).
Yi, C.-R. et al. Observing topological charges and dynamical bulk-surface correspondence with ultracold atoms. Phys. Rev. Lett. 123, 190603 (2019).
Wang, Y. et al. Experimental observation of dynamical bulk-surface correspondence in momentum space for topological phases. Phys. Rev. A 100, 052328 (2019).
Ji, W. et al. Quantum simulation for three-dimensional chiral topological insulator. Phys. Rev. Lett. 125, 020504 (2020).
Chen, B. et al. Digital quantum simulation of floquet topological phases with a solid-state quantum simulator. Photon. Res. 9, 81–87 (2021).
Xin, T. et al. Quantum phases of three-dimensional chiral topological insulators on a spin quantum simulator. Phys. Rev. Lett. 125, 090502 (2020).
Niu, J. et al. Simulation of higher-order topological phases and related topological phase transitions in a superconducting qubit. Sci. Bull. 66, 1168–1175 (2021).
Breuer, H.-P. et al. The Theory of Open Quantum Systems (Oxford University Press on Demand, 2002).
Ashida, Y., Gong, Z. & Ueda, M. Non-hermitian physics. Adv. Phys. 69, 249–435 (2020).
Bergholtz, E. J., Budich, J. C. & Kunst, F. K. Exceptional topology of non-hermitian systems. Rev. Mod. Phys. 93, 015005 (2021).
Heiss, W. The physics of exceptional points. J. Phys. A 45, 444016 (2012).
Lee, T. E. Anomalous edge state in a non-hermitian lattice. Phys. Rev. Lett. 116, 133903 (2016).
Leykam, D., Bliokh, K. Y., Huang, C., Chong, Y. D. & Nori, F. Edge modes, degeneracies, and topological numbers in non-hermitian systems. Phys. Rev. Lett. 118, 040401 (2017).
Xu, Y., Wang, S.-T. & Duan, L.-M. Weyl exceptional rings in a three-dimensional dissipative cold atomic gas. Phys. Rev. Lett. 118, 045701 (2017).
Xiong, Y. Why does bulk boundary correspondence fail in some non-hermitian topological models. J. Phys. Commun. 2, 035043 (2018).
Kunst, F. K., Edvardsson, E., Budich, J. C. & Bergholtz, E. J. Biorthogonal bulk-boundary correspondence in non-hermitian systems. Phys. Rev. Lett. 121, 026808 (2018).
Yao, S., Song, F. & Wang, Z. Non-hermitian chern bands. Phys. Rev. Lett. 121, 136802 (2018).
Yao, S. & Wang, Z. Edge states and topological invariants of non-hermitian systems. Phys. Rev. Lett. 121, 086803 (2018).
Xiao, L. et al. Non-hermitian bulk–boundary correspondence in quantum dynamics. Nat. Phys. 16, 761–766 (2020).
Weidemann, S. et al. Topological funneling of light. Science 368, 311–314 (2020).
Ghatak, A., Brandenbourger, M., Van Wezel, J. & Coulais, C. Observation of non-hermitian topology and its bulk–edge correspondence in an active mechanical metamaterial. Proc. Natl Acad. Sci. 117, 29561–29568 (2020).
Helbig, T. et al. Generalized bulk–boundary correspondence in non-hermitian topolectrical circuits. Nat. Phys. 16, 747–750 (2020).
Wu, Y. et al. Observation of parity-time symmetry breaking in a single-spin system. Science 364, 878–880 (2019).
Zhang, W. et al. Observation of non-hermitian topology with nonunitary dynamics of solid-state spins. Phys. Rev. Lett. 127, 090501 (2021).
Gardiner, C. & Zoller, P. Quantum Noise (Springer-Verlag, 2004).
Gardiner, C. & Zoller, P. The Quantum World of Ultra-Cold Atoms and Light Book I: Foundations of Quantum Optics, vol. 2 (World Scientific Publishing Company, 2014).
Zhang, L., Zhang, L. & Liu, X.-J. Quench-induced dynamical topology under dynamical noise. Phys. Rev. Res. 3, 013229 (2021).
Liu, X.-J., Law, K. T. & Ng, T. K. Realization of 2d spin-orbit interaction and exotic topological orders in cold atoms. Phys. Rev. Lett. 112, 086401 (2014).
Wu, Z. et al. Realization of two-dimensional spin-orbit coupling for bose-einstein condensates. Science 354, 83–88 (2016).
Günther, U. & Samsonov, B. F. Naimark-dilated p t-symmetric brachistochrone. Phys. Rev. Lett. 101, 230404 (2008).
Kawabata, K., Ashida, Y. & Ueda, M. Information retrieval and criticality in parity-time-symmetric systems. Phys. Rev. Lett. 119, 190401 (2017).
Schirmer, R. E. & Noggle, J. H. The Nuclear Overhauser Effect; Chemical Applications (Academic Press, 1971).
Acknowledgements
This work is supported by the National Key Research and Development Program of China (2019YFA0308100, 2021YFA1400900), the National Natural Science Foundation of China (12075110, 11825401, 11975117, 11905099, 11875159, and U1801661), the Guangdong Basic and Applied Basic Research Foundation (2019A1515011383, 2021B1515020070), the Guangdong International Collaboration Program (2020A0505100001), the Science, Technology and Innovation Commission of Shenzhen Municipality (ZDSYS20170303165926217, KQTD20190929173815000, JCYJ20200109140803865, JCYJ20170412152620376, RCYX20200714114522109, and JCYJ20180302174036418), the Open Project of Shenzhen Institute of Quantum Science and Engineering (Grant No.SIQSE202003), the Pengcheng Scholars, the Guangdong Innovative and Entrepreneurial Research Team Program (2019ZT08C044), and the Guangdong Provincial Key Laboratory (2019B121203002). L.Z. also acknowledges support from Agencia Estatal de Investigación (the R&D project CEX2019-000910-S, funded by MCIN/AEI/10.13039/501100011033, Plan National FIDEUA PID2019-106901GB-I00, FPI), Fundació Privada Cellex, Fundació Mir-Puig, and Generalitat de Catalunya (AGAUR Grant No. 2017 SGR 1341, CERCA program).
Author information
Authors and Affiliations
Contributions
D.L. and X.L. supervised the experiments. L.Z. and X.L. elaborated the theoretical framework. Z.L. and X.L. wrote the computer code and accomplished the NMR experiments. All authors analyzed the data, discussed the results, and wrote the manuscript. Z.L., L.Z., and X.L. who made equal contributions to this work are considered “co-first authors".
Corresponding authors
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
Lin, Z., Zhang, L., Long, X. et al. Experimental quantum simulation of non-Hermitian dynamical topological states using stochastic Schrödinger equation. npj Quantum Inf 8, 77 (2022). https://doi.org/10.1038/s41534-022-00587-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41534-022-00587-3
This article is cited by
-
Non-Hermitian non-equipartition theory for trapped particles
Nature Communications (2024)
-
Noisy intermediate-scale quantum computers
Frontiers of Physics (2023)