Abstract
Topology in quench dynamics gives rise to intriguing dynamic topological phenomena, which are intimately connected to the topology of static Hamiltonians yet challenging to probe experimentally. Here we theoretically characterize and experimentally detect momentum-time skyrmions in parity-time \(({\cal{P}}{\cal{T}})\)-symmetric non-unitary quench dynamics in single-photon discrete-time quantum walks. The emergent skyrmion structures are protected by dynamic Chern numbers defined for the emergent two-dimensional momentum-time submanifolds, and are revealed through our experimental scheme enabling the construction of time-dependent non-Hermitian density matrices via direct measurements in position space. Our work experimentally reveals the interplay of \({\cal{P}}{\cal{T}}\) symmetry and quench dynamics in inducing emergent topological structures, and highlights the application of discrete-time quantum walks for the study of dynamic topological phenomena.
Similar content being viewed by others
Introduction
Topological phases feature a wealth of fascinating properties governed by the geometry of their ground-state wave functions at equilibrium1,2, but topological phenomena also manifest as non-equilibrium quantum dynamics in driven-dissipative3 and Floquet systems4,5,6,7, as well as in quench processes8,9,10,11,12,13,14,15. The experimental detection of these dynamic topological phenomena is challenging, since it requires full control and access of the time-evolved state. In recent experiments with ultracold atoms and superconducting qubits, topological objects such as vortices, links, rings, and skyrmions have been identified in the unitary quench dynamics of topological systems via time-resolved and momentum-resolved tomography16,17,18,19. Here we experimentally establish discrete-time photonic quantum walks (QWs) as another promising arena for engineering and detecting dynamic topological phenomena. Compared to other synthetic systems, the relative ease of introducing loss in photonic systems further enables us to experimentally investigate novel dynamic topological phenomena in the non-unitary regime, where parityâtime \(({\cal{P}}{\cal{T}})\) symmetry plays an important role.
In discrete-time photonic QWs20,21,22,23,24,25, single photons, starting from their initial states, are subject to repeated unitary operations26. While QW dynamics support Floquet topological phases (FTPs)22,23,24,25,27,28, discrete-time QWs can also be viewed as a stroboscopic simulation of quench dynamics between FTPs, during which dynamic topological phenomena should occur. As a first endeavor in this direction, two recent experiments report simulation of dynamic quantum phase transitions using QWs29,30, where non-analyticities in the time evolution of physical observables are associated with changes in dynamic topological order parameters31. However, compared to other synthetic systems exhibiting rich dynamic topological structures in higher dimensions16,17,18,19, dynamic topological structures reported in QWs have been limited to one dimension and characterized by topological invariants defined on one-dimensional manifolds.
Here we theoretically characterize and experimentally detect dynamic skyrmion structures, a two-dimensional topological object, in \({\cal{P}}{\cal{T}}\)-symmetric one-dimensional QWs of single photons. Originally proposed in high-energy physics32 and later experimentally observed in magnetic and optical configurations33,34,35, skyrmions are a type of topologically stable defects featuring a three-component vector field in two dimensions. In QW dynamics, dynamic skyrmions manifest themselves in the momentumâtime spin texture of the time-evolved density matrix, and are protected by quantized dynamic Chern numbers in emergent momentumâtime submanifolds12,13,14. To detect dynamic skyrmions, we devise an experimental scheme where time-dependent momentumâspace density matrices of spatially non-localized states are constructed based on a combination of interference-based measurements and projective measurements in position space. Such a practice allows for direct measurements of the density matrix at each time step, which significantly reduces the systematic error introduced by the least-square algorithm in tomographic measurements. We confirm the emergence of dynamic skyrmion structures when QW dynamics correspond to quenches between distinct FTPs in the \({\cal{P}}{\cal{T}}\)-symmetry-unbroken regime, where the dynamics is coherent despite being non-unitary. Effective coherent dynamics is manifested as temporal oscillatory behavior inherent in off-diagonal density-matrix elements. Such oscillatory phenomena reflect the systemâs ability to fully retrieve information temporarily lost to the environment by \({\cal{P}}{\cal{T}}\) dynamics in the unbroken-symmetry regime36. By contrast, when the system is quenched into the \({\cal{P}}{\cal{T}}\)-symmetry-broken regime, skyrmions are absent in the momentumâtime space, as the dynamics become incoherent. Our work unveils the fascinating relation between emergent topology and \({\cal{P}}{\cal{T}}\)-symmetric non-unitary dynamics, and opens up exploration of higher-dimensional dynamic topological structures using QWs.
Results
Quench dynamics in \({\cal{P}}{\cal{T}}\)-symmetric QWs
We experimentally implement \({\cal{P}}{\cal{T}}\)-symmetric non-unitary QWs on a one-dimensional lattice L \(\left( {L \in {\Bbb Z}} \right)\) with single photons in the cascaded interferometric network illustrated in Fig. 1. The corresponding Floquet operator is
where R(θ) rotates coin states (encoded in the horizontal and vertical polarizations of single photons |Hã and |Vã) by θ about the y-axis, and S moves the photon to neighboring spatial modes depending on its polarization (see âMethodsâ). The loss operator \(M = 1_{\mathrm{w}} \otimes \left( {| + \rangle \langle + | + \sqrt {1 - p} | - \rangle \langle - |} \right)\) enforces a partial measurement in the basis \(| \pm \rangle = \left( {|H\rangle | \pm V\rangle } \right)/\sqrt 2\) at each time step with a success probability pâââ[0, 1]. Here \(1_{\mathrm{w}} = \mathop {\sum}\nolimits_x |x\rangle \langle x|\) with |xã (xâââL) denoting the spatial mode. Note that the non-unitary QW driven by U reduces to a unitary one when pâ=â0.
QWs governed by U stroboscopically simulate non-unitary time evolutions driven by the effective Hamiltonian Heff, with \(U = {\mathrm{e}}^{ - iH_{{\mathrm{eff}}}}\). We define the quasienergy ε and eigenstate |Ïã as U|Ïãâ=âγâ1âeâiε|Ïã, where \(\gamma = (1 - p)^{ - \frac{1}{4}}\). U possesses passive \({\cal{P}}{\cal{T}}\) symmetry with \({\cal{P}}{\cal{T}}\gamma U\left( {{\cal{P}}{\cal{T}}} \right)^{ - 1} = \gamma ^{ - 1}U^{ - 1}\), where \({\cal{P}}{\cal{T}} = \mathop {\sum}\nolimits_x {| - x\rangle } \langle x| \otimes \sigma _3{\cal{K}}\), Ï3â=â|HããH|âââ|VããV|, and \({\cal{K}}\) is the complex conjugation. It follows that ε is entirely real in the \({\cal{P}}{\cal{T}}\)-symmetry-unbroken regime, and can take imaginary values in regimes when \({\cal{P}}{\cal{T}}\) symmetry is spontaneously broken29,37,38,39. U also features topological properties, characterized by winding numbers defined through the global Berry phase40,41,42. We show the topological phase diagram of the system in Fig. 2a, where distinct FTPs are marked by their corresponding winding numbers. The boundaries between \({\cal{P}}{\cal{T}}\)-symmetry-unbroken and \({\cal{P}}{\cal{T}}\)-symmetry-broken regimes are also shown in red-dashed lines, with \({\cal{P}}{\cal{T}}\)-symmetry-broken regimes surrounding topological phase boundaries.
To simulate quench dynamics, we initialize the walker photon in the eigenstate |Ïiã of a Floquet operator \(U^{\mathrm{i}} = {\mathrm{e}}^{ - iH_{{\mathrm{eff}}}^{\mathrm{i}}}\), characterized by coin parameters \((\theta _1^{\mathrm{i}},\theta _2^{\mathrm{i}})\). The walker at the tth time step is given by \(|\psi (t)\rangle = {\mathrm{e}}^{ - iH_{{\mathrm{eff}}}t}|\psi ^{\mathrm{i}}\rangle\), such that the resulting QW can be identified as a sudden quench between \(H_{{\mathrm{eff}}}^{\mathrm{i}}\) and Heff. Adopting notations in typical quench dynamics, we denote U and Heff as Uf and \(H_{{\mathrm{eff}}}^{\mathrm{f}}\) in the following, characterized by coin parameters \((\theta _1^{\mathrm{f}},\theta _2^{\mathrm{f}})\).
Fixed points and emergent skyrmions
Due to the lattice translational symmetry of Ui,f, dynamics in different quasi-momentum k-sectors are decoupled. We denote pre-quench and post-quench Floquet operators in each k-sector as \(U_k^{\mathrm{{{i}}}}\) and \(U_k^{\mathrm{f}}\), respectively, whose eigenstates are \(|\psi _{k, \pm }^{{\mathrm{i}},{\mathrm{f}}}\rangle\). Quasienergies of \(U_k^{{\mathrm{i}},{\mathrm{f}}}\) are denoted as \(\varepsilon _{k, \pm }^{{\mathrm{i}},{\mathrm{f}}}\), with \(\varepsilon _{k, \pm }^{{\mathrm{i}},{\mathrm{f}}} = \pm E_k^{{\mathrm{i}},{\mathrm{f}}}\). Focusing on the case where Ui is in the \({\cal{P}}{\cal{T}}\)-symmetry-unbroken regime, we write the initial state as \(|\psi _{k, - }^{\mathrm{i}}\rangle\) in each k-sector.
By invoking the biorthogonal basis43, the non-unitary time evolution of the system is captured by a non-Hermitian density matrix, which can be written as14
where n(k, t)â=â(n1, n2, n3), Ïâ=â(Ï1, Ï2, Ï3), \({\mathbf{\tau }}_i = \mathop {\sum}\nolimits_{\mu ,\nu = \pm } {|\psi _{k,\mu }^{\mathrm{f}}\rangle } \sigma _i^{\mu \nu }\langle \chi _{k,\nu }^{\mathrm{f}}|\) (iâ=â0, 1, 2, 3), and \(\langle \chi _{k,\mu }^{\mathrm{f}}|\) \(\left( {|\psi _{k,\mu }^{\mathrm{f}}\rangle } \right)\) is the left (right) eigenvector of \(U_k^{\mathrm{f}}\). Here Ï0 is a 2âÃâ2 identity matrix, and Ïi (iâ=â1, 2, 3) is the corresponding standard Pauli matrix.
A key advantage of adopting Eq. (2) is that n(k, t) becomes a real unit vector, which enables us to visualize the non-unitary dynamics on a Bloch sphere. As illustrated in Fig. 2b, when \(E_k^{\mathrm{f}}\) is real, n(k, t) rotates around poles of the Bloch sphere with a period \(t_0 = \pi /E_k^{\mathrm{f}}\). Thus, momenta corresponding to poles of the Bloch sphere are identified as two different kinds of fixed points, where the density matrices do not evolve in time. In contrast, when \(E_k^{\mathrm{f}}\) is imaginary, there are no fixed points in the dynamics, as n(k, t) asymptotically approaches the north pole in the long-time limit (see Fig. 2b).
When Ui and Uf belong with distinct FTPs in the \({\cal{P}}{\cal{T}}\)-symmetry-unbroken regime, fixed points of different kinds necessarily emerge in pairs14,29. Each momentum submanifold between a pair of distinct fixed points can be combined with the S1 topology of the periodic time evolution to form an emergent S2 momentumâtime manifold, which can be mapped to the S2 Bloch sphere of n(k, t). The Chern number characterizing such an S2âââS2 mapping is finite and gives rise to intriguing skyrmion structures in the emergent momentumâtime manifolds.
To probe fixed points and momentumâtime skyrmions, we perform projective measurements and interference-based measurements to construct the Hermitian density matrix Ïâ²(k, t)â=â|Ïk(t)ããÏk(t)|. This is achieved by writing \(\rho {\prime}(k,t) = \frac{1}{2}\mathop {\sum}\limits_{j = 0}^3 {{\mathrm{ }}\mathop {\sum}\limits_{x_1,x_2} {{\mathrm{e}}^{ - ik(x_1 - x_2)}} {\mathrm{ }}} \langle \psi _{x_2}(t)|\sigma _j|\psi _{x_1}(t)\rangle \sigma _j\), where |Ïx(t)ã is the coin state on site x at the tth time step. We experimentally measure \(\langle \psi _{x_2}(t)|\sigma _j|\psi _{x_1}(t)\rangle\) (jâ=â0, 1, 2, 3) for each pair of positions x1 and x2 directly. We then calculate the non-Hermitian density matrix Ï(k, t) from Ïâ²(k, t) and determine n(k, t) through n(k, t)â=âTr[Ï(k, t)Ï].
Whereas the Hermitian density matrix Ïâ²(k, t) is experimentally accessible, visualization of non-unitary dynamics on the Bloch sphere starting from Ïâ²(k, t) requires normalization. It is also difficult to characterize the connection between fixed points and dynamic skyrmions using Ïâ²(k, t). By contrast, such a connection is revealed elegantly with the biorthogonal construction of Ï(k, t), as we detail in âMethodsâ and Supplementary Note 4. This is why we characterize skyrmions using the spin texture n(k, t) associated with Ï(k, t).
Dynamics in the \({\cal{P}}{\cal{T}}\)-symmetry-unbroken regime
We first study fixed points and momentumâtime skyrmions in the \({\cal{P}}{\cal{T}}\)-symmetry-unbroken regime. For comparison, we also experimentally characterize these quantities in unitary dynamics. We initialize the walker on a localized lattice site |xâ=â0ã and in the coin state \(\left| {\psi _ - ^{\mathrm{i}}} \right\rangle _c\). Here, |xã denotes the spatial mode. Importantly, \(\left| {\psi _{k, - }^{\mathrm{i}}} \right\rangle = \left| {\psi _ - ^{\mathrm{i}}} \right\rangle _{\mathrm{c}}\) is an eigenstate of \(U_k^{\mathrm{i}}\) for all k, with the corresponding \((\theta _1^{\mathrm{i}},\theta _2^{\mathrm{i}})\) on black dashed lines in Fig. 2a. Without loss of generality, we choose \((\theta _1^{\mathrm{i}} = \pi /4,\theta _2^{\mathrm{i}} = - \pi /2)\) for both the unitary and non-unitary cases.
For the first case of study, we implement unitary QWs with \(\left| {\psi _ - ^{\mathrm{i}}} \right\rangle _{\mathrm{c}} = (\left|H\right\rangle + i\left|V\right\rangle)/\sqrt 2\) and \((\theta _1^{\mathrm{f}} = - \pi /2,\theta _2^{\mathrm{f}} = \pi /3)\), which simulate quench processes between FTPs with νiâ=â0 and νfâ=ââ2. We have chosen \((\theta _1^{\mathrm{f}},\theta _2^{\mathrm{f}})\) on purple dashed lines, where qusienergy bands are flat. Oscillatory dynamics of n(k, t) in different k-sectors thus feature the same period, as illustrated in Fig. 3a. We identify fixed points of unitary dynamics at high-symmetry points of the Brillioun zone {âÏ, âÏ/2, 0, Ï/2}, where n(k, t) become independent of time.
For the second case of study, we implement non-unitary QWs with pâ=â0.36, \(\left| {\psi _ - ^{\mathrm{i}}} \right\rangle _{\mathrm{c}} = 0.7606\left| H \right\rangle + 0.6492i\left| V \right\rangle\), and \(\left[ {\theta _1^{\mathrm{f}} = - \pi /2,\theta _2^{\mathrm{f}} = {\mathrm{arcsin}}\left( {\frac{1}{\alpha }{\mathrm{cos}}\frac{\pi }{6}} \right)} \right]\) (here \(\alpha = \frac{\gamma }{2}(1 + \sqrt {1 - p} )\)). The post-quench FTP is in the \({\cal{P}}{\cal{T}}\)-symmetry unbroken regime with νfâ=ââ2. As shown in Fig. 3b, dynamics of n(k, t) is still oscillatory, but fixed points are shifted away from the high-symmetry points, consistent with theoretical predictions.
In Fig. 4, we plot n(k, t) in the momentumâtime space. The oscillatory behavior in n(k, t) is then manifested as momentumâtime skyrmions, which are protected by dynamic Chern numbers defined on the corresponding momentumâtime submanifold. By contrast, when the system is quenched between FTPs with the same winding number, skyrmion-lattice structures are no longer present, as shown in the Supplementary Fig. 4. Comparing Fig. 4a, b, we see that dynamic skyrmion structures in the unitary and the \({\cal{P}}{\cal{T}}\)-symmetric non-unitary quench processes are qualitatively similar; albeit, in the non-unitary case, skyrmion structures are slightly deformed due to the shift of fixed points. However, as we show later, skyrmion structures are generically absent when the system is quenched into the \({\cal{P}}{\cal{T}}\)-symmetry-broken regime.
Our configuration is also capable of simulating quench dynamics between topologically non-trivial FTPs. To demonstrate this, we implement non-unitary QWs with \(\left| {\psi _ - ^{\mathrm{i}}} \right\rangle = ( - 0.1111 - 0.6983i)\left| 0 \right\rangle \otimes \left| + \right\rangle - \frac{1}{{\sqrt 2 }}\left| { - 2} \right\rangle \otimes \left| - \right\rangle\) and \(\left( {\theta _1^{\mathrm{f}} = \pi /2,\theta _2^{\mathrm{f}} = - {\mathrm{arcsin}}\left( {\frac{1}{\alpha }{\mathrm{{cos}}}\frac{\pi }{6}} \right)} \right)\). The initial state is the eigenstate of Ui with \((\theta _1^{\mathrm{i}} = - \pi /2,\theta _2^{\mathrm{i}} = - \pi /4)\), which is prepared by performing gate operations prior to QW dynamics (see âMethodsâ). The pre-quench and post-quench FTPs are in the \({\cal{P}}{\cal{T}}\)-symmetry unbroken regime with νiâ=ââ2 and νfâ=â2, respectively. In Fig. 5, we show the corresponding oscillations of n(k, t), as well as the momentumâtime skyrmions. Consistent with theoretical predictions (see âMethodsâ), a total of eight fixed points exist in the dynamics, which divide the first Brillioun zone into eight submanifolds with skyrmions emerging in momentumâtime space in between each adjacent pair of fixed points.
Dynamics in the \({\cal{P}}{\cal{T}}\)-symmetry-broken regime
We now turn to the case where Uf belong with the \({\cal{P}}{\cal{T}}\)-symmetry-broken regime. We first initialize the walker on a localized lattice site in the coin state \(\left( {\left| H \right\rangle + \left| V \right\rangle } \right)/\sqrt 2\) and evolve it under Uf characterized by \(\left( {\theta _1^{\mathrm{f}} = - \pi /2,\theta _2^{\mathrm{f}} = \frac{1}{2}\left( {\pi - {\mathrm{arccos}}\frac{1}{\alpha }} \right)} \right)\), which is in the \({\cal{P}}{\cal{T}}\)-symmetry-broken regime with a completely imaginary quasienergy spectrum. As shown in Fig. 6a, there is no periodical evolution in n(k, t) anymore. Instead, different components of n(k, t) slowly approach a steady state with nâ=â(0, 0, 1) in the long-time limit. This is more clearly seen in momentumâtime space shown in Fig. 6b, where skyrmion structures are absent and vectors in all k-sectors tend to point out of the plane in the long-time limit. We note that dynamics of n(k, t) here is insensitive to the choice of initial state, as the system always relaxes to the steady state at long times.
We then start from the same initial state as in Fig. 3b and evolve it under a different Uf characterized by \((\theta _1^{\mathrm{f}} = - 0.39\pi ,\theta _2^{\mathrm{f}} = 0.3864\pi )\), which is in the \({\cal{P}}{\cal{T}}\)-symmetry-broken regime with a partially real (and partially imaginary) quasienergy spectrum. As shown in Fig. 6c, d, whereas the evolution of n(k, t) is still periodic for k-sectors with real quasienergies, the existence of steady-state approaching k-sectors with imaginary quasienergies deform spin textures in those sectors and destroy the overall dynamic skyrmion structures in the momentumâtime space.
Discussion
By simulating quench dynamics of topological systems using photonic QWs, we have revealed emergent momentumâtime skyrmions, protected by dynamic Chern numbers defined on the momentumâtime submanifolds. These emergent topological phenomena are underpinned by fixed points of dynamics, which can exist for both unitary and non-unitary quench processes if the system is quenched across a topological phase boundary. We have further confirmed the decisive role of \({\cal{P}}{\cal{T}}\)-symmetry on the existence of fixed points and skyrmions in non-unitary dynamics.
Emergent momentumâtime skyrmions and dynamic Chern numbers reported here are intimately connected with topological entanglement-spectrum crossing13, as well as the recently observed dynamic quantum phase transitions and dynamic topological order parameters19,29,30. In particular, whereas dynamic topological order parameters and dynamic Chern numbers originate from completely different topological constructions, both are hinged upon the presence of fixed points in quench dynamics14,19. With the highly flexible control of photonic QW protocols, it would be interesting to investigate dynamic topological phenomena in higher dimensions or associated with other topological classifications in the future44. Our work thus paves the way for a systematic experimental study of dynamic topological phenomena in both unitary and non-unitary dynamics.
Methods
Experimental implementation of U
We implement the coin operator \(R(\theta ) = 1_{\mathrm{w}} \otimes {\mathrm{e}}^{ - i\theta \sigma _2}\), the shift operator \(S = \mathop {\sum}\nolimits_x \left( {\left| {x - 1} \right\rangle \left\langle x \right| \otimes \left| H \right\rangle \left\langle H \right| + \left| {x + 1} \right\rangle \left\langle x \right| \otimes \left| V \right\rangle \left\langle V \right|} \right)\), and the loss operator \(M = 1_{\mathrm{w}} \otimes \left( {\left| + \right\rangle \left\langle + \right| + \sqrt {1 - p} \left| - \right\rangle \left\langle - \right|} \right)\), following the approach outlined in ref.â29. Here, \(\left| \pm \right\rangle = (|H\rangle \pm |V\rangle )/\sqrt 2\), Ï2â=âi(â|HããV|+|VããH|) is the standard Pauli operator under the polarization basis, x (xâââL) denotes the spatial mode, \(1_{\mathrm{w}} = \mathop {\sum}\nolimits_x {|x\rangle } \langle x|\), and the loss parameter pâ=â0.36 for non-unitary QWs in our experiment.
Preparation of quasi-local initial state
To simulate quench dynamics between topologically non-trivial FTPs in Fig. 5, the QW dynamics starts from a quasi-local initial state \(\left| {\psi _ - ^{\mathrm{i}}} \right\rangle = ( - 0.1111 - 0.6983i)\left| 0 \right\rangle \otimes \left| + \right\rangle - \frac{1}{{\sqrt 2 }}\left| { - 2} \right\rangle \otimes \left| - \right\rangle\), which is a Ui eigenstate with νiâ=ââ2. As illustrated in Fig. 1, the state is experimentally prepared in the following way. First, starting from |xâ=â0ã, the polarization of single photons is prepared in \(\left( {\frac{1}{{\sqrt 2 }}, - 0.1111 - 0.6983i} \right)^{\mathrm{T}}\) by tuning the setting angles of the wave plates right after the first PBS. Then we use a BD to split the photons into the spatial modes |xâ=â0ã and |xâ=ââ2ã depending on their polarizations. Finally, a HWP at a fixed angle rotates polarizations of single photons in both spatial modes.
\({{\cal{P}}}{{\cal{T}}}\)-symmetric non-unitary QW
The non-unitary Floquet operator \(U\) in Eq. (1) has passive \({\cal{P}}{\cal{T}}\) symmetry, from which we can define \(\tilde U = \gamma U\) with \(\gamma = (1 - p)^{ - \frac{1}{4}}\). \(\tilde U\) has active \({\cal{P}}{\cal{T}}\) symmetry, with the symmetry operator \({\cal{P}}{\cal{T}} = \mathop {\sum}\nolimits_x {| - x\rangle } \langle x| \otimes \sigma _3{\cal{K}}\), where \({\cal{K}}\) is the complex conjugation. As homogeneous QWs have lattice translational symmetry, we write \(\tilde U\) in momentum space
where \(\alpha = \frac{\gamma }{2}(1 + \sqrt {1 - p} ),\beta = \frac{\gamma }{2}(1 - \sqrt {1 - p} )\).
Eigenvalues of \(\tilde U_k\) are \(\lambda _{k, \pm } = d_0 \mp i\sqrt {1 - d_0^2}\), and the corresponding quasienergy εk,±â=âi ln(λk,±). When \(d_0^2 < 1\) for all k, the quasienergy is real, and the system is in the \({\cal{P}}{\cal{T}}\)-symmetry-unbroken regime. Whereas if \(d_0^2 \ge 1\) for some k, the \({\cal{P}}{\cal{T}}\) symmetry is spontaneously broken and the quasienergy is imaginary in the corresponding momentum range.
Winding numbers of non-unitary QWs
Non-unitary QWs governed by U possess topological properties, which are characterized by winding numbers defined through the global Berry phase νâ=âÏB/2Ï. Here ÏBâ=âÏZ+â+âÏZâ, with generalized Zak phases
The integral above is over the first Brillioun zone and ãÏk,μ| and |Ïk,μã (μâ=â±) are, respectively, the left and right eigenstates of Uk, defined through \(U_k^\dagger |\chi _{k,\mu }\rangle = \lambda _\mu ^ \ast |\chi _{k,\mu }\rangle\) and Uk|Ïk,μãâ=âλμ|Ïk,μã, respectively.
Besides \({\cal{P}}{\cal{T}}\) symmetry, the system also possesses particleâhole symmetry \({\cal{K}}U{\cal{K}}^{ - 1} = U\). According to refs.â27,45, the complete topological characterization of the system should be Z â Z, where two topological invariants exist. Following ref.â46, we define another winding number νⲠby choosing a different time frame with Floquet operator \(U\prime = M^{\frac{1}{2}}R\left( {\frac{{\theta _2}}{2}} \right)SR\left( {\theta _1} \right)SR\left( {\frac{{\theta _2}}{2}} \right)M^{\frac{1}{2}}\), where \(M^{\frac{1}{2}} = 1_{\mathrm{w}} \otimes \left( {\left| + \right\rangle \left\langle + \right| + (1 - p)^{\frac{1}{4}}\left| - \right\rangle \left\langle - \right|} \right)\). Note that νⲠcan be calculated from the global Berry phase associated with Uâ². Combinations of ν and νⲠthen give the bulk topological invariants (ν0, νÏ)â=â[(νâ+âνâ²)/2, (νâνâ²)/2]46, which dictate the correct number of edge states with quasienergies Reεâ=â0 and Reεâ=âÏ, respectively, through the bulk-boundary correspondence. However, for the quench dynamics of homogeneous QWs with no boundaries, we find that it is sufficient to discuss the interplay between existence of momentumâtime skyrmions, the pre-quench and post-quench winding numbers, and \({\cal{P}}{\cal{T}}\) symmetry in a single time frame. Indeed, as we demonstrate experimentally in Supplementary Note 2, whereas spin textures and locations of fixed points are quantitatively different under different time frames, going to a different time frame does not qualitatively change the conclusions in the main text, as long as the appropriate winding number is used.
Fixed points and dynamical Chern numbers in QW dynamics
Non-unitary time evolution of the system is captured by the non-Hermitian density matrix
where the time-evolved state \(|\psi _k(t)\rangle = \mathop {\sum}\nolimits_{\mu = \pm } {c_\mu } {\mathrm{e}}^{ - i\varepsilon _{k,\mu }^{\mathrm{f}}t}|\psi _{k,\mu }^{\mathrm{f}}\rangle\), the associated state \(\langle \chi _k(t)|: = \mathop {\sum}\nolimits_\mu {c_\mu ^ \ast } {\mathrm{e}}^{i\varepsilon _{k,{\kern 1pt} u}^{{\mathrm{f}} \ast }t}\langle \chi _{k,\mu }^{\mathrm{f}}|\), \(c_\mu = \langle \chi _{k,\mu }^{\mathrm{f}}|\psi _{k, - }^{\mathrm{i}}\rangle\), and \(\langle \chi _{k,\mu }^{\mathrm{f}}|\) \(\left( {|\psi _{k,\mu }^{\mathrm{f}}\rangle } \right)\) is the left (right) eigenvector of \(U_k^{\mathrm{f}}\), with the biorthonormal conditions \(\langle \chi _{k,\mu }^{\mathrm{f}}|\psi _{k,\nu }^{\mathrm{f}}\rangle = \delta _{\mu \nu }\) and \(\mathop {\sum}\nolimits_\mu |\psi _{k,\mu }^{\mathrm{f}}\rangle \langle \chi _{k,\mu }^{\mathrm{f}}| = 1\). It follows that n(k, t)â=âTr[Ï(k, t) Ï], with {Ïi} satisfying the standard \({\frak{s}}{\frak{u}}(2)\) commutation relations.
Following the convention of the main text, we denote the final Flouqet operator in each quasimomentum k-sector as \(U_k^{\mathrm{f}}\) and the corresponding quasienergy as \(\pm E_k^{\mathrm{f}}\). When \(E_k^{\mathrm{f}}\) is real, we have
By contrast, when \(E_k^{\mathrm{f}}\) is imaginary, and assuming \({\mathrm{Im}}(E_k^{\mathrm{f}}) > 0\), we have
From these expressions, it is straightforward to visualize dynamics of n(k, t) on a Bloch sphere as illustrated in Fig. 2b and discussed in the main text. In particular, when Uf is in the \({\cal{P}}{\cal{T}}\)-symmetry-unbroken regime, fixed points occur at momenta with cââ=â0 or c+â=â0, which we identify as two different types of fixed points. The total number of these fixed points with c+â=â0 or cââ=â0 which are topologically protected is |νfâνi| each14,29. As a concrete example, for the quench illustrated in Fig. 5 with νiâ=ââ2 and νfâ=â2, we observe eight topologically protected fixed points, with half of them cââ=â0 and the other half c+â=â0. It follows that dynamic skyrmions exist in momentumâtime space in between each adjacent pair of fixed points of different kinds.
Dynamic Chern number
When Uf is in the \({\cal{P}}{\cal{T}}\)-symmetry-unbroken regime, periodic evolution of the density matrix in each k-sector gives rise to a temporal S1 topology. In the presence of fixed points, each momentum submanifold between two adjacent fixed points can be combined with the S1 topology in time to form a momentumâtime submanifold S2, which can be mapped to the Bloch sphere associated with the vector n(k, t). These S2âââS2 mappings define a series of dynamic Chern numbers
where km and kn denote two neighboring fixed points, and \(t_0 = \pi /E_k^{\mathrm{f}}\). For quenches between Hamiltonians with different winding numbers, dynamic Chern numbers are quantized, with values dependent on the nature of fixed points at km and kn: Cmnâ=â1 when c+(km)â=â0 and câ(kn)â=â0; Cmnâ=ââ1 when câ(km)â=â0 and c+(kn)â=â0. When the two fixed points are of the same kind, Cmnâ=â0.
According to its definition in Eq. (8), a finite Chern number in a momentumâtime submanifold corresponds to the existence of momentumâtime skyrmions in the same submanifold.
Constructing density matrix from direct measurements
The non-Hermitian density matrix Ï(k, t) is related to the Hermitian one Ïâ²(k, t):â=â|Ïk(t)ããÏk(t)| through
where we have used the biorthonormal conditions \(\langle \chi _{k,\mu }^{\mathrm{f}}|\psi _{k,\nu }^{\mathrm{f}}\rangle = \delta _{\mu \nu }\) and \(\mathop {\sum}\limits_{\mu = \pm } {|\psi _{k,\mu }^{\mathrm{f}}\rangle } \langle \chi _{k,\mu }^{\mathrm{f}}| = 1\).
As we have discussed in the main text, we first determine Ïâ²(k, t) by measuring \(\left\langle {\psi _{x_2}(t)} \right|\sigma _j\left| {\psi _{x_1}(t)} \right\rangle\) (jâ=â0, 1, 2, 3) for each pair of positions x1 and x2. In the case of x1â=âx2, we perform projective measurements on the polarizations of photons at each position. In the case of x1ââ âx2, we employ interference-based measurements to construct the matrix element \(\langle \psi _{x_2}(t)|\sigma _j|\psi _{x_1}(t)\rangle\) from experimental data. We then construct Ï(k, t) and n(k, t) from Ïâ²(k, t) using Eq. (9).
We note that dynamic skyrmions are also visible in spin textures defined through nâ²(k, t)â=âTr[Ïâ²(k, t)] after enforcing normalization on nâ²(k, t). Here the Pauli vectorâÏ =â(Ï1, Ï2, Ï3). As we demonstrate explicitly in Supplementary Fig. 5, whereas n(k, t) and nâ²(k, t) are distinct from each other, dynamic topological structures inherent in these two spin textures are equivalent and emerge under the same conditions.
Data availability
Experimental data, any related experimental background information not mentioned in the text and other findings of this study are available from the corresponding author upon reasonable request.
References
Hasan, M. Z. & Kan, C. L. Colloqium: topological insulators. Rev. Mod. Phys. 82, 3045â3067 (2010).
Qi, X. L. & Zhang, S. C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 1057â1110 (2011).
Diehl, S., Rico, E., Baranov, M. A. & Zoller, P. Topology by dissipation in atomic quantum wires. Nat. Phys. 7, 971 (2011).
Rudner, M. S., Lindner, N. H., Berg, E. & Levin, M. Anomalous edge states and the bulk-edge correspondence for periodically driven two-dimensional systems. Phys. Rev. X 3, 031005 (2013).
DâAlessio, L. & Rigol, M. Dynamical preparation of Floquet Chern insulators. Nat. Commun. 6, 8336 (2015).
Titum, P., Berg, E., Rudner, M. S., Refael, C. & Lindner, N. H. Anomalous FloquetâAnderson insulator as a nonadiabatic quantized charge pump. Phys. Rev. X 6, 021013 (2016).
Potter, A. C., Morimoto, T. & Vishwanath, A. Classification of interacting topological Floquet phases in one dimension. Phys. Rev. X 6, 041001 (2016).
Caio, M. D., Cooper, N. R. & Bhaseen, M. J. Quantum quenches in Chern insulators. Phys. Rev. Lett. 115, 236403 (2015).
Hu, Y., Zoller, P. & Budich, J. C. Dynamical buildup of a quantized Hall response from nontopological states. Phys. Rev. Lett. 117, 126803 (2016).
Wilson, J. H., Song, J. C. W. & Refael, G. Remnant geometric Hall response in a quantum quench. Phys. Rev. Lett. 117, 235302 (2016).
Wang, C., Zhang, P., Chen, X., Yu, J. & Zhai, H. Scheme to measure the topological number of a Chern insulator from quench dynamics. Phys. Rev. Lett. 118, 185701 (2017).
Yang, C., Li, L. & Chen, S. Dynamical topological invariant after a quantum quench. Phys. Rev. B 97, 060304(R) (2018).
Gong, Z. & Ueda, M. Entanglement-spectrum crossing and momentumâtime Skyrmions in quench dynamics. Phys. Rev. Lett. 121, 250601 (2018).
Qiu, X., Deng, T. -S., Hu, Y., Xue, P. & Yi, W. Fixed points and emergent topological phenomena in a parityâtime-symmetric quantum quench. Preprint at http://arxiv.org/abs/1806.10268 (2018).
Zhang, L., Zhang, L., Niu, S. & Liu, X.-J. Dynamical classification of topological quantum phases. Sci. Bull. 63, 1385â1391 (2018).
Fläschner, N. et al. Observation of dynamical vortices after quenches in a system with topology. Nat. Phys. 14, 265 (2018).
Tarnowski, M. et al. Measuring topology from dynamics by obtaining the Chern number from a linking number. Nat. Commun. 10, 1728 (2019).
Sun, W. et al. Uncover topology by quantum quench dynamics. Phys. Rev. Lett. 121, 250403 (2018).
Guo, X.-Y. et al. Observation of dynamical quantum phase transition by a superconducting qubit simulation. Phys. Rev. Appl. 11, 044080 (2019).
Kitagawa, T., Rudner, M. S., Berg, E. & Demler, E. Exploring topological phases with quantum walks. Phys. Rev. A 82, 033429 (2010).
Cardano, F. et al. Statistical moments of quantum-walk dynamics reveal topological quantum transitions. Nat. Commun. 7, 11439 (2016).
Kitagawa, T. et al. Observation of topologically protected bound states in photonic quantum walks. Nat. Commun. 3, 882 (2012).
Cardano, F. et al. Detection of Zak phases and topological invariants in a chiral quantum walk of twisted photons. Nat. Commun. 8, 15516 (2017).
Barkhofen, S. et al. Measuring topological invariants in disordered discrete-time quantum walks. Phys. Rev. A 96, 033846 (2017).
Wang, X. P. et al. Detecting topological invariants and revealing topological phase transitions in discrete-time photonic quantum walks. Phys. Rev. A 98, 013835 (2018).
Peruzzo, A. et al. Quantum walks of correlated photons. Science 329, 1500â1503 (2010).
Xiao, L. et al. Observation of topological edge states in parityâtime-symmetric quantum walks. Nat. Phys. 13, 1117 (2017).
Zhan, X. et al. Detecting topological invariants in nonunitary discrete-time quantum walks. Phys. Rev. Lett. 119, 130501 (2017).
Wang, K. et al. Simulating dynamic quantum phase transitions in photonic quantum walks. Phys. Rev. Lett. 122, 020501 (2019).
Xu, X. Y. et al. Measuring a dynamical topological order parameter in quantum walks. Preprint at https://arxiv.org/abs/1808.03930 (2018).
Heyl, M. Dynamical quantum phase transitions: a review. Rep. Prog. Phys. 81, 054001 (2018).
Skyrme, T. H. R. A unified field theory of mesons and baryons. Nucl. Phys. 31, 556 (1962).
Mühlbauer, S. et al. Skyrmion lattice in a chiral magnet. Science 323, 915 (2009).
Yu, X. Z. et al. Real-space observation of a two-dimensional skyrmion crystal. Nature 465, 901 (2010).
Tsessses, S. et al. Optical skyrmion lattice in evanescent electromagnetic fields. Science 361, 6406 (2018).
Kawabata, K., Ashida, Y. & Ueda, M. Information retrieval and criticality in parityâtime-symmetric systems. Phys. Rev. Lett. 119, 190401 (2017).
Bender, C. M. & Boettcher, S. Real spectra in non-Hermitian Hamiltonians having PT symmetry. Phys. Rev. Lett. 80, 5243â5246 (1998).
Bender, C. M., Brody, D. C. & Jones, H. F. Complex extension of quantum mechanics. Phys. Rev. Lett. 89, 270401 (2002).
Bender, C. M. Making sense of non-Hermitian Hamiltonians. Rep. Prog. Phys. 70, 947â1018 (2007).
Garrison, J. & Wright, E. Complex geometrical phases for dissipative systems. Phys. Lett. A 128, 177 (1988).
Liang, S.-D. & Huang, G.-Y. Topological invariance and global Berry phase in non-Hermitian systems. Phys. Rev. A. 87, 012118 (2013).
Lieu, S. Topological phases in the non-Hermitian SuâSchiefferâHeeger model. Phys. Rev. B 97, 045106 (2018).
Brody, D. C. Biorthogonal quantum mechanics. J. Phys. A 47, 035305 (2014).
Chang, P.-Y. Topology and entanglement in quench dynamics. Phys. Rev. B 97, 224304 (2018).
Kim, D., Mochizuki, K., Kawakami, N. & Obuse, H. Floquet topological phases driven by PT symmetric nonunitary time evolution. Preprint at http://arxiv.org/abs/1609.09650 (2016).
Asbóth, J. K. & Obuse, H. Bulk-boundary correspondence for chiral symmetric quantum walks. Phys. Rev. B 88, 121406(R) (2013).
Acknowledgements
This work has been supported by the Natural Science Foundation of China (Grant Nos. 11674056 and 11522545) and the Natural Science Foundation of Jiangsu Province (Grant No. BK20160024). W.Y. acknowledges support from the National Key Research and Development Program of China (Grant Nos. 2016YFA0301700 and 2017YFA0304100).
Author information
Authors and Affiliations
Contributions
K.K.W. performed the experiments with contributions from L.X., X.Z. and Z.H.B. W.Y. developed the theoretical aspects and performed the theoretical analysis with contribution from X.Z.Q. and wrote part of the paper. P.X. designed the experiments, analyzed the results, and wrote part of the paper. B.C.S. revised the paper. K.K.W. and X.Z.Q. contributed equally to this work.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Journal peer review information: Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisherâs note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
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
Wang, K., Qiu, X., Xiao, L. et al. Observation of emergent momentumâtime skyrmions in parityâtime-symmetric non-unitary quench dynamics. Nat Commun 10, 2293 (2019). https://doi.org/10.1038/s41467-019-10252-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-019-10252-7