Abstract
Superconducting quantum circuits are a natural platform for quantum simulations of a wide variety of important lattice models describing topological phenomena, spanning condensed matter and high-energy physics. One such model is the bosonic analog of the well-known fermionic Kitaev chain, a 1D tight-binding model with both nearest-neighbor hopping and pairing terms. Despite being fully Hermitian, the bosonic Kitaev chain exhibits a number of striking features associated with non-Hermitian systems, including chiral transport and a dramatic sensitivity to boundary conditions known as the non-Hermitian skin effect. Here, using a multimode superconducting parametric cavity, we implement the bosonic Kitaev chain in synthetic dimensions. The lattice sites are mapped to frequency modes of the cavity, and the in situ tunable complex hopping and pairing terms are created by parametric pumping at the mode-difference and mode-sum frequencies, respectively. We experimentally demonstrate important precursors of nontrivial topology and the non-Hermitian skin effect in the bosonic Kitaev chain, including chiral transport, quadrature wavefunction localization, and sensitivity to boundary conditions. Our experiment is an important first step towards exploring genuine many-body non-Hermitian quantum dynamics.
Similar content being viewed by others
Introduction
While the development of universal fault-tolerant quantum computers is still ongoing, analog quantum simulation (AQS) has emerged as a promising approach to study classically intractable quantum systems1,2,3. In AQS, the simulations take place on an artificial quantum system built to have the same Hamiltonian as the system of interest. One appealing aspect of AQS is that quantum degrees of freedom can be represented natively. For instance, while bosonic degrees of freedom with infinite Hilbert spaces can be naturally simulated using oscillator modes on an AQS platform, boson-to-qubit mapping on a small-scale qubit-based computer typically requires a very large overhead in the number of physical qubits and required gates to implement bosonic operators4,5,6,7.
An important class of AQS is lattice models with topological properties, which are intrinsically plagued by the infamous “sign problem” and thus are unsuitable for quantum Monte Carlo methods8,9,10. Topological systems have been a focus of intense research, both theoretically and experimentally, for some time. More recently the study of classical and quantum topological physics in non-Hermitian systems has attracted significant interest11,12. These systems exhibit intriguing and distinct properties owing to their nonorthogonal eigenstates and singularities in the complex eigenvalue spectrum of non-Hermitian matrices13. Compelling phenomena include oscillations between eigenstates14,15, unidirectional invisibility16, high-performance lasers17, and enhanced sensitivity for potential sensing applications18. Remarkably, non-Hermiticity fundamentally alters concepts such as symmetry and energy gaps inherited from Hermitian physics, giving rise to an enriched variety of topological phases with no Hermitian counterpart13.
There are various approaches to non-Hermitian dynamics19,20,21,22. In the quantum regime, non-Hermiticity appears naturally in open quantum systems described by a Lindblad master equation: conditioned on the absence of a quantum jump, state evolution is described by an effective non-Hermitian Hamiltonian23,24. Experimentally resolving potentially interesting effects from such non-Hermitian Hamiltonians, while possible25,26,27, is challenging. By definition, one must post-select on measurement records without a jump, and such trajectories become exponentially rare at long times. Consequently, one must perform many runs of the experiment to acquire enough data for adequate statistics28. An alternate route to obtaining non-Hermitian quantum dynamics is through the use of Hermitian bosonic Hamiltonians. With unitary squeezing and antisqueezing terms, the equations of motion become effectively non-Hermitian despite the Hermiticity of the underlying Hamiltonians29,30,31. These Hamiltonians present an interesting avenue for probing coherent, genuinely quantum non-Hermitian effects, without the use of dissipation32.
Previously, we have demonstrated the feasibility of an AQS platform based on a multimode superconducting parametric cavity by simulating a plaquette of the bosonic Creutz ladder33. The cavity modes share a superconducting quantum interference device (SQUID) which acts as a common boundary condition. Parametric modulation of the boundary condition induces complex “hopping” couplings that allow us to create a programmable graph of connected (coupled) modes, realizing a lattice in synthetic dimensions. By controlling the phases of the complex hopping terms, we showed that our platform can implement interesting features including static gauge fields and topological phenomena. We note similar work has been done in the context of classical non-Hermitian optics34,35 and chiral photon transport36,37.
In this work, we expand our programmable AQS toolbox by introducing pairing terms between modes in the target Hamiltonian in addition to the hopping terms. They allow us to construct a topologically nontrivial Hamiltonian in synthetic dimensions, the bosonic Kitaev chain (BKC) introduced in ref. 38. Working with a 3-site chain, we experimentally demonstrate phase-dependent chiral transport, quadrature wavefunction localization, and a strong sensitivity to boundary conditions. These observations serve as precursors to nontrivial topology and the much sought-after non-Hermitian skin effect (NHSE)32,39,40,41. Our work enlarges the set of topologically nontrivial models that can be simulated under genuinely quantum conditions, further highlighting the potential of our AQS platform.
Results
We now review the features of the BKC38. In analogy with the celebrated fermionic Kitaev chain42, a spinless p-wave topological superconductor in one dimension, the BKC is a 1D bosonic tight-binding model with both nearest-neighbor hopping and pairing terms38. As a consequence of the pairing terms, the system has effectively non-Hermitian equations of motion and supports phase-dependent chiral transport. Furthermore, it has a topologically nontrivial phase where each band of the bulk spectrum of the dynamical matrix has nonzero winding in the complex energy plane. This is accompanied by the remarkable property of the NHSE43,44,45,46: the entire spectrum depends sensitively on boundary conditions, and the wavefunctions under open boundary conditions are localized at the ends of the chain. (See Supplementary Note 3 for a detailed discussion of the symmetry of the model, the topological invariant, and the topological protection against symmetry-preserving disorder).
The BKC is described by the 1D tight-binding Hamiltonian
illustrated in Fig. 1a. Here \({\hat{a}}_{j}\) is a bosonic annihilation operator on site j, and \(t{e}^{i{\varphi }_{t}}\) and iΔ are respectively the complex hopping and pairing strengths between adjacent sites, with magnitudes t and Δ. Without loss of generality, we choose the pairing strength to be purely imaginary.
To demonstrate that the model supports phase-dependent chiral transport, we consider the Heisenberg equations of motion for the Hermitian position and momentum quadratures \({\hat{x}}_{j}=({\hat{a}}_{j}+{\hat{a}}_{j}^{{{{\dagger}}} })/\sqrt{2}\) and \({\hat{p}}_{j}=-i({\hat{a}}_{j}-{\hat{a}}_{j}^{{{{\dagger}}} })/\sqrt{2}\):
In the special case φt = ± 90°, as long as Δ ≠ 0, the equations of motion correspond exactly to two decoupled non-Hermitian Hatano–Nelson chains with asymmetric left-right hopping for the x and p quadratures47. This is an example of effective non-Hermitian dynamics in Hermitian systems31,48. The decoupling of the x and p quadratures, along with the left-right coupling asymmetry in each chain, gives rise to phase-dependent chiral propagation, with quadratures representing the two chiral species. The chirality becomes perfect in the Δ → t limit: for φt = 90° (−90°), the x (p) quadrature propagates to the right and the p (x) quadrature to the left. From Eqs. (2) to (3), we see that away from φt = ± 90° a finite \(\cos {\varphi }_{t}\) term mixes the left and right-moving quadratures, leading to an overall reduction in the chiral nature of the transport38.
Let us consider the spectrum of the dynamical matrix, i.e., the coefficient matrix of Eqs. (2)–(3), taking into account a uniform onsite single-photon loss rate κ. Eqs. (2)–(3) are now understood as the Heisenberg-Langevin equations of the expectation values \(\langle {\hat{x}}_{j}\rangle\) and \(\langle {\hat{p}}_{j}\rangle\), and acquire the damping terms \(-\kappa \langle {\hat{x}}_{j}\rangle /2\) and \(-\kappa \langle {\hat{p}}_{j}\rangle /2\).
We can perform a spatially uniform squeezing transformation to identify a topological phase transition at \(t| \cos {\varphi }_{t}|=\Delta\). For \(t| \cos {\varphi }_{t}| \, > \, \Delta\) the transformed Hamiltonian for the squeezed fields \({\hat{\beta }}_{j}\) is equivalent to a model with only hopping. In the more interesting case of \(t| \cos {\varphi }_{t}| \, < \, \Delta\), the transformed Hamiltonian assumes the form of Eq. (1) with renormalized parameters \({t}^{{\prime} }=t\sin {\varphi }_{t}\), \({\Delta }^{{\prime} }=\sqrt{{\Delta }^{2}-{t}^{2}{\cos }^{2}{\varphi }_{t}}\) and \({\varphi }_{t}^{{\prime} }=9{0}^{\circ }\).
Exploring the second case further, we consider both open and periodic boundary conditions. The complex energy spectrum for the open chain reads
where kn = nπ/(N + 1), n = 1, 2, …, N and N is the number of sites. The open chain spectrum Eq. (4) is independent of the coupling phase φt, and can be obtained through a second, position-dependent squeezing transformation38. This position dependence is reflected in the quadrature wavefunctions which do depend on φt:
where \({\widehat{d}}_{n}\) is the annihilation operator for the eigenstate corresponding to momentum kn, \({\hat{x}}_{j}^{{\prime} }=({\hat{\beta }}_{j}+{\hat{\beta }}_{j}^{{{{\dagger}}} })/\sqrt{2}\), and \({\hat{p}}_{j}^{{\prime} }=-i({\hat{\beta }}_{j}-{\hat{\beta }}_{j}^{{{{\dagger}}} })/\sqrt{2}\). It is clear from the exponential factors that different quadrature components of a given wavefunction are localized at opposite edges of the chain in the topological phase.
The spectrum for periodic boundary conditions is found by Fourier transforming the dynamical matrix,
where kn = 2nπ/N, n = 1, 2, …, N. In the topological phase \(t| \cos {\varphi }_{t}| < \Delta\), Eq. (6) forms an ellipse in the complex energy plane, yielding a nonzero winding number around the point E = −iκ/2 for each band. Comparing with Eq. (4), we see that the spectrum of the system depends drastically on its boundary conditions, a feature which holds for arbitrary system size (see Fig. 1d). Together with the (phase-dependent) localization in the open-chain wavefunctions in Eq. (5), it thus serves as a paradigmatic example of the NHSE43,44.
Equations (4) and (6) imply that, provided
the system is dynamically stable under open boundary conditions (i.e., all energy eigenvalues have nonpositive imaginary parts), but unstable under periodic boundary conditions. Such instability in the presence of loss stems from the effective non-Hermiticity induced by pairing, and provides direct evidence for the nontrivial topology of the system.
We now describe the experimental realization of the BKC on our AQS platform (see Fig. 1b and ref. 33). We program a chain of three sites in the synthetic frequency dimension. For the open chain, the sites are connected by two links where each link is created by two coherent pumps: a pump at the modes’ frequency difference \({\omega }_{j,j+1}^{t}=| {\omega }_{j}-{\omega }_{j+1}|\) to activate the hopping and a pump at the sum frequency \({\omega }_{j,j+1}^{\Delta }={\omega }_{j}+{\omega }_{j+1}\) to activate the pairing. The magnitudes and phases of these pump tones in turn determine the magnitudes and phases of the complex hopping and pairing terms. To impose periodic boundary conditions on the 3-site chain, we create an additional link with two more pump tones that connect the open ends, forming a closed chain (Fig. 1c).
We measure the spectra of both chains using a vector network analyzer (VNA), determining the eigenmode frequencies directly from the reflection coefficients. We further characterize the open chain using phase-sensitive transport measurements. Sending in a tone set at a constant magnitude but with a phase that ramps at a constant rate from −180° to 180°, we probe at various site frequencies to measure signal transport in synthetic dimensions. The phase-sensitivity of the transport converts the phase sweep of the input signal into magnitude variations in the output signals. Additional details on device fabrication and the measurement setup are given in Supplementary Notes 1 and 2.
Twisted-tubes picture
Before discussing the details of our transport measurements, we present a twisted-tubes picture (Fig. 2) explaining the role of individual hopping and pairing phases. Experimentally, while all pump phases are separately tunable, they are difficult to calibrate absolutely at the sample, which is deep in the cryostat. The twisted-tubes picture is an intuitive way to understand how adjusting the pump phases affects the dynamics of the system.
We begin from a single link a–b with hopping \({t}_{ab}{e}^{i{\varphi }_{ab}^{t}}\) and pairing \({\Delta }_{ab}{e}^{i{\varphi }_{ab}^{\Delta }}\). For future convenience, we define the sum and difference phases, \({\varphi }_{ab}^{\pm }=({\varphi }_{ab}^{t}\pm {\varphi }_{ab}^{\Delta })/2\). The transport properties are solved using a phase-dependent input-output theory (see Methods). The transported signal strength from a to b depends on the input phase. We refer to the input quadrature at a that maximizes transport, together with the corresponding output quadrature at b, as the favored quadratures in transport from a to b. On the other hand, the orthogonal quadratures dominate the transported signal in the opposite direction. This is represented in Fig. 2 as a pair of interleaved directional tubes with (squeezed) elliptic cross sections, where red (blue) tubes transport to the right (left). The favored quadratures in a given transport direction are along the major axis of the corresponding tube. We can twist the tubes by varying \({\varphi }_{ab}^{-}\) and \({\varphi }_{ab}^{+}\), which individually rotate the a and b ends, respectively (Fig. 2a).
In a 3-site chain with two connected links, the relative orientation of the two tubes at site b is quantified by the phase \(\Theta={\varphi }_{ab}^{+}+{\varphi }_{bc}^{-}\) (Fig. 2b), which is invariant under local gauge transformations (see Methods). When the two tubes are maximally misaligned (\(\Theta={0}^{\circ }\,{{{{{{\mathrm{mod}}}}}}}\,\,18{0}^{\circ }\), Fig. 2c), transport along the chain from a to c or from c to a is independent of the input phase. On the other hand, when the two tubes are exactly aligned (\(\Theta=9{0}^{\circ }\,{{{{{{\mathrm{mod}}}}}}}\,\,18{0}^{\circ }\), Fig. 2d), we achieve the maximum transport magnitude for the favored input phase and maximum suppression of the orthogonal phase. The values \(\Theta=9{0}^{\circ }\,{{{{{{\mathrm{mod}}}}}}}\,\,18{0}^{\circ }\) and \(\Theta={0}^{\circ }\,{{{{{{\mathrm{mod}}}}}}}\,\,18{0}^{\circ }\) thus represent two limiting cases of the 3-site chain, which we refer to as the “chiral” chain and the “trivial” chain, respectively. In the translationally invariant BKC in Eq. (1), Θ becomes exactly φt and the chiral chain is topologically nontrivial.
We calibrate the 3-site chain by injecting a signal into the central site (b) and observing the transport to the ends while twisting the b end of the a–b tube by varying \({\varphi }_{ab}^{+}\), effectively changing Θ (see Fig. 3). The phase is calibrated by sweeping it and comparing the complete response to theory, similar to Fig. 4. The calibrated phases are then used for subsequent measurements. Furthermore, we calibrate the magnitudes to be reflection and transport coefficients by normalizing them to an input of unit magnitude during the fitting process (see Methods). Figure 3 shows the measured phase-dependent transport for two extremal cases \({\varphi }_{ab}^{+}={0}^{\circ }\) and 90°. The approximately sinusoidal shape of the magnitudes of the transport coefficients is a manifestation of the transport sensitivity to the input phase, with the maxima (minima) corresponding to the favored (suppressed) phases. At \({\varphi }_{ab}^{+}=9{0}^{\circ }\) (red), the transport magnitudes to a and c are in phase, indicating that the same input phases are favored or suppressed. This corresponds to a trivial chain as the ends of the tubes are completely misaligned at b.
More interestingly, at \({\varphi }_{ab}^{+}={0}^{\circ }\) (blue), the twisted tubes are aligned and we realize a chiral chain. The experimental signature of this is that the transport magnitudes to the ends become completely out of phase. That is, when a certain input phase is enhanced in transport to a, its transport to c is highly suppressed, and vice versa. This corresponds to the chiral chain with the ends of the two tubes aligned at b.
We investigate the chain chirality by performing phase-sensitive transport measurements along the length of the chain while varying \({\varphi }_{ab}^{+}\). The results, in Fig. 4, reveal the predicted chiral transport properties. The chiral regime emerges gradually as we slowly vary \({\varphi }_{ab}^{+}\) away from the trivial case \({\varphi }_{ab}^{+}=9{0}^{\circ }\), reaching maximum chirality at \({\varphi }_{ab}^{+}={0}^{\circ }\).
Figure 5 illustrates the line cuts corresponding to the trivial and chiral chains. In the trivial chain (red), the transported signals are nearly reciprocal and show only minimal changes in response to the varying input phase.
The constant transport magnitude can serve as a baseline for the chiral transport. Conversely, in the chiral chain (blue curves), the transport magnitudes between the ends exhibit strong sensitivity to the input phase. In Fig. 5a, for instance, the transport from the site a to site c is enhanced for ϕ = 0°, compared to the baseline, while the transport of the orthogonal input phases at ϕ = −90° and 90° is highly suppressed. The ratios of enhanced and suppressed magnitudes, compared to the baseline, are expected to be \(\frac{t+\Delta }{t-\Delta }\) and \(\frac{t-\Delta }{t+\Delta }\), respectively, which give us a rough estimate of Δ ≈ 0.4t.
An additional signature of chirality, illustrated in Fig. 5a, is the flattening of the transport phase. Specifically, referring to the right panel of (a) and the left panel of (b), we see that in the trivial case (red), the output phase is approximately linear, the same as the input phase. Conversely in the chiral case (blue), the output phase has a stairstep shape, approximately locking to the phase of the preferred quadrature before jumping by ±180° (which is the same quadrature with opposite amplitude). We see that the observed behavior agrees well with the theoretical predictions.
We further examine the non-Hermitian topology by extracting the x and p quadrature wavefunctions (see Methods). The trivial and chiral chains exhibit a striking difference in the spatial support of the x and p wavefunctions (Fig. 5c, d). In the trivial case, the quadratures are delocalized with nearly equal weights on both ends. In stark contrast, we observe the characteristics of the NHSE in the chiral case: the x and p wavefunctions are strongly localized at the right and left ends, respectively. This demonstrates how we can control non-Hermitian topological effects by changing the gauge-invariant phase of our chain.
We examine the sensitivity of the chiral chain to boundary conditions by connecting the chain ends (Fig. 1d). These measurements are done on a separate device from the above measurements (see Methods). We measure the reflection coefficient around the frequency of site b while varying both hopping and pairing phases of the a–c link (Fig. 6a–c). The dominant pattern is the spectra braiding as a function of \({\varphi }_{ac}^{t}\), which is determined by the loop phase33. However, for certain phase conditions, we observe discontinuities in the central branch of the spectrum (Fig. 6b–d), indicating that the chain is approaching dynamical instability (see Methods).
It is remarkable that the instability is determined solely by the link phases (the pump magnitudes are constant), revealing the transition to the chiral regime of the closed chain. In the twisted-tubes picture, we start with misaligned tubes (Fig. 6e), then rotate both a and c of the a–c tube until it aligns with both the a–b and b–c tubes, forming two directional loops (Fig. 6f). Sufficient alignment of the link phases is necessary to realize nontrivial winding, and consequently, instabilities under periodic boundary conditions (Fig. 1d).
We can give an intuitive picture of the instability in terms of circulating gain. In the absence of onsite loss, when the tubes are aligned, an initial excitation traverses the loop, being amplified indefinitely, resulting in dynamical instability. If the tubes are misaligned, the circulating signal will be amplified through one link but then deamplified through the next, allowing for a stable steady state. Local loss, as in our chain, simply shifts the eigenvalues down in the complex plane (see Fig. 1), increasing the instability threshold for Δ. (The threshold is \(\Delta=\frac{\kappa }{2}\)). We have confirmed that for higher values of Δ, the system becomes unstable, leading to coupled parametric oscillations of the modes.
Discussion
We have demonstrated on our AQS platform that nontrivial non-Hermitian topological systems can be realized using parametric down-conversion. Unlike dissipation-induced non-Hermiticity, this approach allows us to preserve some of the Hamiltonian symmetries, such as time-reversal symmetry. Therefore, this platform can be used to explore the rich topological phases and symmetries of non-Hermitian systems. Furthermore, since the dynamics are the results of a coherent process, as opposed to dissipation, in a Hermitian Hamiltonian31, our AQS platform can implement genuine quantum dynamics with effective non-Hermiticity. An interesting future direction for our platform is to explore the interplay between interactions, topology and non-Hermiticity49,50,51,52,53. In fact, the BKC exhibits nonlinear dynamics in the above-threshold regime, where we have observed coupled parametric oscillations in the system. This would serve as a first example of non-Hermitian non-linear quantum dynamics.
Furthermore, the 3-site chain can also be used in interesting applications. We can utilize the chiral features as a phase-dependent quantum amplifier38,54,55. Alternatively, vacuum squeezing can be used to realize entangled multimode states, which is a complex resource that can play a central role in continuous-variable quantum information processing, e.g., Gaussian boson sampling56. In addition, non-Hermitian systems have a wide range of applications from quantum sensing57 to entanglement creation and control58.
We discussed future directions to improve our AQS platform in earlier work33. Briefly, we can improve the hardware efficiency by increasing the length of the parametric cavity. This would increase the number of modes that can be utilized as lattice sites in AQS. In addition, a single parametric cavity can be used as a sublattice in a network of coupled cavities. The flexibility of controlling local coupling phases allows us to simulate topological models requiring nontrivial phase conditions, which gives our platform advantages over competing platforms.
Methods
Input-output theory
Here we describe in detail the input-output theory used to quantitatively study the transport in our system. Additional information on the fitting procedure is provided in Supplementary Note 4.
We denote the cavity mode j as \({\hat{a}}_{j}\) and its bare frequency as \({\omega }_{j}^{(0)}\). To couple modes j and \({j}^{{\prime} }\), we apply a beam-splitter pump with frequency \({\omega }_{j{j}^{{\prime} }}^{t}\approx | {\omega }_{j}^{(0)}-{\omega }_{{j}^{{\prime} }}^{(0)}|\) and a down-conversion pump with frequency \({\omega }_{j{j}^{{\prime} }}^{\Delta }\, \approx \, {\omega }_{j}^{(0)}+{\omega }_{{j}^{{\prime} }}^{(0)}\); in the rotating-wave approximation, they generate the hopping and pairing terms in the Hamiltonian, respectively. More concretely, we choose a set of frequencies ωj, such that \({\omega }_{j}\, \approx \, {\omega }_{j}^{(0)}\) for all modes, and \({\omega }_{j{j}^{{\prime} }}^{t}=| {\omega }_{j}-{\omega }_{{j}^{{\prime} }}|\), \({\omega }_{j{j}^{{\prime} }}^{\Delta }={\omega }_{j}+{\omega }_{{j}^{{\prime} }}\) for all pumps; since there are only as many free variables ωj as the number of modes, generally the pump frequencies are not independent of each other.
We should distinguish between signal and idler frequencies in the presence of down-conversion pumps. Given the probe detuning Ω, an input signal at the frequency ωj + Ω is coupled to other signal frequencies \({\omega }_{{j}^{{\prime} }}+\Omega\) by beam-splitter pumps, and to other idler frequencies \({\omega }_{{j}^{{\prime} }}-\Omega\) by down-conversion pumps. In the frequency-domain spectrum measurements, we send in a coherent state as the input signal at ωj + Ω and detect the reflected signal at the same frequency; the spectrum is mapped out as a function of the probe detuning Ω. On the other hand, in the phase-dependent transport measurements, we set Ω = 0 such that the signal and idler frequencies coincide for every mode, which generates interference between signal and idler frequencies and results in phase-dependent transport. We send in a coherent tone as the input signal at the frequency ωj while scanning the input phase, and detect the reflected signal at ωj and transported signal at \({\omega }_{{j}^{{\prime} }}\).
After the rotating wave approximation, the general quadratic Hamiltonian that can be programmed in our AQS takes the following form:
Here we work in the reference frame rotating at the frequency ωj at mode j, and the pump detunings are defined as \(\delta {\omega }_{j}={\omega }_{j}-{\omega }_{j}^{(0)}\). The second sum runs over all pairs of connected sites j and \({j}^{{\prime} }\), and the hopping and pairing strengths \({t}_{j{j}^{{\prime} }}\) and \({\Delta }_{j{j}^{{\prime} }}\) have tunable phases \({\varphi }_{j{j}^{{\prime} }}^{t}=-{\varphi }_{{j}^{{\prime} }j}^{t}\) and \({\varphi }_{j{j}^{{\prime} }}^{\Delta }={\varphi }_{{j}^{{\prime} }j}^{\Delta }\), respectively.
Under the local gauge transformation \({\hat{a}}_{j}\to U{\hat{a}}_{j}{U}^{{{{\dagger}}} }={\hat{a}}_{j}{e}^{i{\theta }_{j}}\), the link phases in the transformed Hamiltonian \({\hat{{{{{{{{\mathcal{H}}}}}}}}}}_{S}^{{\prime} }=U{\hat{{{{{{{{\mathcal{H}}}}}}}}}}_{S}{U}^{{{{\dagger}}} }\) become \({\varphi }_{j{j}^{{\prime} }}^{t}\to {\varphi }_{j{j}^{{\prime} }}^{t}-{\theta }_{j}+{\theta }_{{j}^{{\prime} }}\) and \({\varphi }_{j{j}^{{\prime} }}^{\Delta }\to {\varphi }_{j{j}^{{\prime} }}^{\Delta }+{\theta }_{j}+{\theta }_{{j}^{{\prime} }}\). Therefore, in agreement with our twisted-tubes picture, the phases \({\varphi }_{j{j}^{{\prime} }}^{\pm }=({\varphi }_{j{j}^{{\prime} }}^{t}\pm {\varphi }_{j{j}^{{\prime} }}^{\Delta })/2\) transform as \({\varphi }_{j{j}^{{\prime} }}^{+}\to {\varphi }_{j{j}^{{\prime} }}^{+}+{\theta }_{{j}^{{\prime} }}\) and \({\varphi }_{j{j}^{{\prime} }}^{-}\to {\varphi }_{j{j}^{{\prime} }}^{-}-{\theta }_{j}\); furthermore, for two links j–\({j}^{{\prime} }\) and \({j}^{{\prime} }\)–j ″, the phase \({\varphi }_{j{j}^{{\prime} }}^{+}+{\varphi }_{{j}^{{\prime} }{j}^{{\prime\prime} }}^{-}\) is fully gauge invariant.
Input signals are sent in through a measurement line, which is coupled to the system by the Hamiltonian
Here \({\kappa }_{j}^{{{{{{{{\rm{ext}}}}}}}}}\) is the external coupling rate to the input mode \({\hat{a}}_{{{{{{{{\rm{in}}}}}}}},j}\). Taking the onsite single-photon loss rate κj into account, we can solve the Heisenberg-Langevin equations of motion for the signal/idler modes \({\hat{a}}_{j}^{{{{{{{{\rm{S/I}}}}}}}}}\) in terms of the corresponding input modes \({\hat{a}}_{{{{{{{{\rm{in}}}}}}}},j}^{{{{{{{{\rm{S/I}}}}}}}}}\):
Finally, employing the input-output relation
we express the output modes in terms of the input modes,
Transport properties are now given by the scattering matrix \({S}_{j{j}^{{\prime} }}^{\alpha \beta }\). In the frequency-domain spectrum measurements, each idler frequency generally only receives a vacuum input, which is negligible compared to the coherent state input at the signal frequencies; thus the measured reflection amplitude at mode j is simply \({S}_{jj}^{{{{{{{{\rm{SS}}}}}}}}}\). On the other hand, in the phase-dependent transport measurements at zero probe detuning Ω = 0, the signal and idler frequencies always have identical inputs. The transport from \({j}^{{\prime} }\) to j is therefore characterized by the normalized phase-dependent transport coefficient with unit input,
where ϕ is the input phase at mode \({j}^{{\prime} }\).
We extract the x and p quadrature wavefunctions in chiral and trivial chains in Fig. 5c, d as follows. For an input signal with phase ϕ = 0° at site b, the x wavefunction at the site j = a, b, c is related to the x quadrature of the output mode via the input-output relation \(\langle {\hat{x}}_{j}\rangle=(\langle {\hat{x}}_{{{{{{{{\rm{out}}}}}}}},j}\rangle+\langle {\hat{x}}_{{{{{{{{\rm{in}}}}}}}},j}\rangle )/\sqrt{{\kappa }_{j}^{{{{{{{{\rm{ext}}}}}}}}}}\) (cf. Eq. (12)), where \(\langle {\hat{x}}_{{{{{{{{\rm{in}}}}}}}},j}\rangle\) is nonzero for j = b. To compare the different cases, we normalize the weights of the wavefunction such that \({\sum }_{j}{\langle {\hat{x}}_{j}\rangle }^{2}=1\). We repeat the same process to extract the p wavefunctions.
2-mode and 3-mode open chains
In this section we consider the chiral transport in 2-mode and 3-mode BKCs with open boundary conditions in the framework of the input-output theory. For a single link a–b with hopping \({t}_{ab}{e}^{i{\varphi }_{ab}^{t}}\) and pairing \({\Delta }_{ab}{e}^{i{\varphi }_{ab}^{\Delta }}\), solving Eqs. (10) and (11), we obtain the transport coefficient Eq. (14) from a to b as
and the transport coefficient from b to a as
Here, the common denominator is
and we define the sum and difference phases as before, \({\varphi }_{ab}^{\pm }=({\varphi }_{ab}^{t}\pm {\varphi }_{ab}^{\Delta })/2\). As \({\Delta }_{ab}\to \sqrt{{t}_{ab}^{2}+{\kappa }_{a}{\kappa }_{b}}\), Eq. (17) indicates the transport coefficients increase rapidly and the system approaches instability, which is consistent with the upper bound of Δ in Eq. (7) for N = 2.
Equations (15) and (16) clearly show phase-dependent chirality, as explained in Fig. 2. The input phase \(\phi={\varphi }_{ab}^{-}\) maximizes the magnitude of Eq. (15), and corresponds to the output phase \(9{0}^{\circ }-{\varphi }_{ab}^{+}\). These constitute the favored quadratures in the transport from a to b (major axes of the red ellipses). By contrast, the orthogonal quadratures with the input phase \(\phi=9{0}^{\circ }+{\varphi }_{ab}^{-}\) and the output phase \(18{0}^{\circ }-{\varphi }_{ab}^{+}\) are suppressed (minor axes of the red ellipses). The favored quadratures from a to b may be different at both sites, e.g., appearing as x at site a but as p at site b, xa → pb; it is \({\varphi }_{ab}^{-}\) (\({\varphi }_{ab}^{+}\)) that determines the favored quadrature at site a (b). On the other hand, the quadratures suppressed in the transport from a to b are favored in the transport from b to a (major axes of the blue ellipses), in our example pa ← xb. By programming the coupling phases on a single link, we can tune continuously between different transport scenarios: for instance, starting from favored quadratures xa → pb and pa ← xb, we can vary \({\varphi }_{ab}^{+}\) to arrive at xa → xb and pa ← pb, or vary \({\varphi }_{ab}^{-}\) to arrive at pa → pb and xa ← xb.
In an open chain with two links a–b and b–c, the transport coefficient from b to a is still given by Eq. (16), except that the real prefactor is replaced by \(2\sqrt{{\kappa }_{a}^{{{{{{{{\rm{ext}}}}}}}}}{\kappa }_{b}^{{{{{{{{\rm{ext}}}}}}}}}}{\kappa }_{c}/{D}_{3}^{{{{{{{{\rm{o}}}}}}}}}\), where
Similarly, making the substitutions b → c, a → b and replacing the real prefactor with \(2\sqrt{{\kappa }_{b}^{{{{{{{{\rm{ext}}}}}}}}}{\kappa }_{c}^{{{{{{{{\rm{ext}}}}}}}}}}{\kappa }_{a}/{D}_{3}^{{{{{{{{\rm{o}}}}}}}}}\) in Eq. (15), we obtain the transport coefficient from b to c. As discussed in the main text, the favored quadrature from b to a is \(\phi=-{\varphi }_{ab}^{+}\) while that from b to c is \(\phi={\varphi }_{bc}^{-}\), and the difference between the two is given by the gauge-invariant phase \(\Theta={\varphi }_{ab}^{+}+{\varphi }_{bc}^{-}\). We also mention that Eq. (18) is consistent with Eq. (7) for N = 3.
The phase dependence of the transport coefficient from a to c is slightly more complicated:
While \({\varphi }_{ab}^{-}\) (\({\varphi }_{bc}^{+}\)) determines the favored quadrature at site a (c), the other two link phases enter Eq. (19) only through the gauge-invariant linear combination Θ. When Θ = ± 90° (\(\Theta={0}^{\circ }\,{{{{{{\mathrm{mod}}}}}}}\,\,18{0}^{\circ }\)), a constructive (destructive) interference ensues, and we find a strong (weak) dependence on the input phase ϕ. All the above points are consistent with the twisted tubes picture.
3-mode closed chain
In this section we study how a closed 3-mode chain can approach dynamical instability when we tune its link phases. Here we choose to examine the determinant of the coefficient matrix of the 6 coupled Heisenberg-Langevin equations, Eqs. (10) and (11), which come from the input-output theory; solving an eigenvalue problem will lead to the same quantity. We focus on vanishing pump detuning Ω = 0, where discontinuities are seen to arise in the spectrum in Fig. 6. For a closed chain of three links, a–b, b–c and c–a, the determinant is found as
where we have defined the gauge-invariant combinations \({\Theta }_{a}={\varphi }_{ca}^{+}+{\varphi }_{ab}^{-}\), \({\Theta }_{b}={\varphi }_{ab}^{+}+{\varphi }_{bc}^{-}\) and \({\Theta }_{c}={\varphi }_{bc}^{+}+{\varphi }_{ca}^{-}\), and “perm.” represents the permutation-symmetric contributions (abc) → (bca), (cab). The last term \({C}_{3}^{{{{{{{{\rm{c}}}}}}}}}\) is independent of phases,
Note that the link phases appear exclusively through the three gauge-invariant combinations. In the open chain limit of tbc = Δbc = 0, we can explicitly verify that all phase dependence drops out. Also, if all pairing strengths vanish, only the loop phase Θa + Θb + Θc remains.
When all pairing terms are turned off, simple algebra shows that \({D}_{3}^{{{{{{{{\rm{c}}}}}}}}}\) is positive definite as a function of pump phases. Once the pairing strength exceeds a threshold value, \({D}_{3}^{{{{{{{{\rm{c}}}}}}}}}\) can become zero for some link phases and Eqs. (10) and (11) become singular: this is the point where the closed chain turns unstable and the linear theory fails. This motivates us to find the minimum of \({D}_{3}^{{{{{{{{\rm{c}}}}}}}}}\) with respect to all link phases.
Equation (20) allows an analytical minimization with respect to Θa by expanding out the cosine and sine terms,
Since cosines range between −1 and 1, it is clear that the minimum with respect to Θb and Θc is found at \(\cos 2{\Theta }_{b}=\cos 2{\Theta }_{c}=-\!1\). Consistent with the permutation symmetry, at this point we also have \(\cos 2{\Theta }_{a}=-\!1\). This is reminiscent of the chiral transport regime in the open 3-mode system, achieved at Θ = ±90°. The minimum of \({D}_{3}^{{{{{{{{\rm{c}}}}}}}}}\) with respect to link phases reads
In particular, when all hopping/pairing terms have the same strength tbc = tca = tab = t and Δbc = Δca = Δab = Δ and all loss rates are equal κa = κb = κc = κ,
Therefore, if the pairing strength satisfies Δ > κ/2, the closed chain can become unstable for some choices of link phases. This is consistent with the stability/instability condition Eq. (7).
Following the procedure in the main text, we align the a–b and b–c tubes at b by fixing \(\cos 2{\Theta }_{b}=-\!1\). Consequently
where we have used \({\Theta }_{c}+{\Theta }_{a}=-{\varphi }_{ac}^{t}+{\varphi }_{ab}^{-}+{\varphi }_{bc}^{+}\) and \({\Theta }_{c}-{\Theta }_{a}=-{\varphi }_{ac}^{\Delta }-{\varphi }_{ab}^{-}+{\varphi }_{bc}^{+}\). Equation (25) indicates that provided \(\cos 2{\Theta }_{b}=-\!1\), \({D}_{3}^{{{{{{{{\rm{c}}}}}}}}}\) has two nonequivalent minima in the \({\varphi }_{ac}^{t}\)–\({\varphi }_{ac}^{\Delta }\) plane, found at \(\cos ({\varphi }_{ac}^{t}-{\varphi }_{ab}^{-}-{\varphi }_{bc}^{+})=-\cos ({\varphi }_{ac}^{\Delta }+{\varphi }_{ab}^{-}-{\varphi }_{bc}^{+})=\pm 1\). This agrees with Fig. 6b, c, which shows both \({\varphi }_{ac}^{t}\) and \({\varphi }_{ac}^{\Delta }\) change by 180° between the two discontinuities in the spectrum.
Calibration and characterization
An arbitrary chain is calibrated by activating each link separately while turning off the remaining links. First, upon activating the hopping term of the link j–\({j}^{{\prime} }\), the single mode resonance splits into two resonances whose frequency difference gives twice the coupling strength, \(2{t}_{j{j}^{{\prime} }}\). We choose \({t}_{j{j}^{{\prime} }}\) to be in the strong-coupling regime: the resonance splitting is greater than photon decay rates, \(2{t}_{j{j}^{{\prime} }}\, > \,{\kappa }_{j},{\kappa }_{{j}^{{\prime} }}\), such that the split resonances are resolved. Furthermore, we set \({t}_{j{j}^{{\prime} }}\) to be roughly equal along the chain.
Then, we activate each link’s pairing term. In this case, there is no simple spectral feature that quantifies the pairing strength, \({\Delta }_{j{j}^{{\prime} }}\). We roughly calibrate \({\Delta }_{j{j}^{{\prime} }}\) using transport measurements in the following way. With both hopping and pairing terms applied, we sweep the input phase of the signal described above at site j and measure the contrast of the transport at site \({j}^{{\prime} }\) as a function of input phase. Here the contrast, defined as \(({A}_{\max }-{A}_{\min })/({A}_{\max }+{A}_{\min })\) where \({A}_{\max }\) (\({A}_{\min }\)) is the maximum (minimum) magnitude of the transported signal, evaluates to \(\min \{t,\Delta \}/\max \{t,\Delta \}\) according to Eqs. (15) and (16). We then vary the pairing pump power, interpreting the power where the observed contrast is maximum as \({\Delta }_{j{j}^{{\prime} }}\, \approx \, {t}_{j{j}^{{\prime} }}\). As we further increase the pairing pump power, the system eventually becomes dynamically unstable, as discussed below Eq. (17). We choose a pairing strength that satisfies the stability/instability condition Eq. (7) for N = 3; this means the 3-site open chain is dynamically stable, and the spectrum of the 3-chain closed chain shows a discontinuity at some pump phases (Fig. 6).
The phases of the signals acquire arbitrary offsets from traveling through the measurement lines between the instrument and the sample. We calibrate the input and output phases of the system as follows. First, we connect the 3-site chain with hopping and pairing terms, and perform phase-sensitive transport measurements when a signal is injected in the center site b. We then tune the chain gauge-invariant phase Θ until we realize the chiral chain. To conveniently present the results, we finally set the right-moving quadratures as I and the left-moving quadratures as Q in the chiral chain.
Data availability
The data that support the findings of this study are available from the corresponding author upon request.
Code availability
The code used for data analysis and simulation are available from the corresponding author upon request. Because the analysis code is closely integrated with instrument-control code, which depends on specialized instrument drivers, it is not meaningful to share the code in a general repository.
References
Daley, A. J. et al. Practical quantum advantage in quantum simulation. Nature 607, 667 (2022).
Chiu, C. S. et al. String patterns in the doped hubbard model. Science 365, 251 (2019).
Semeghini, G. et al. Probing topological spin liquids on a programmable quantum simulator. Science 374, 1242 (2021).
Macridin, A., Spentzouris, P., Amundson, J. & Harnik, R. Electron-phonon systems on a universal quantum computer. Phys. Rev. Lett. 121, 110504 (2018).
Sawaya, N. P. et al. Resource-efficient digital quantum simulation of d-level systems for photonic, vibrational, and spin-s hamiltonians. npj Quantum Inf. 6, 49 (2020).
Bauer, C. W. et al. Quantum simulation for high-energy physics. PRX Quantum 4, 027001 (2023).
Peng, B. et al. Quantum simulation of boson-related hamiltonians: techniques, effective hamiltonian construction, and error analysis. arXiv preprint arXiv:2307.06580 (2023).
Loh Jr, E. et al. Sign problem in the numerical simulation of many-electron systems. Phys. Rev. B 41, 9301 (1990).
Iazzi, M., Soluyanov, A. A. & Troyer, M. Topological origin of the fermion sign problem. Phys. Rev. B 93, 115102 (2016).
Smith, A., Golan, O. & Ringel, Z. Intrinsic sign problems in topological quantum field theories. Phys. Rev. Res. 2, 033515 (2020).
Ashida, Y., Gong, Z. & Ueda, M. Non-hermitian physics. Adv. Phys. 69, 249 (2020).
Bergholtz, E. J., Budich, J. C. & Kunst, F. K. Exceptional topology of non-hermitian systems. Rev. Modern Phys. 93, 015005 (2021).
Kawabata, K., Shiozaki, K., Ueda, M. & Sato, M. Symmetry and topology in non-hermitian physics. Phys. Rev. X 9, 041015 (2019).
Makris, K. G., El-Ganainy, R., Christodoulides, D. & Musslimani, Z. H. Beam dynamics in \({{{{{{{\mathcal{P}}}}}}}}{{{{{{{\mathcal{T}}}}}}}}\) symmetric optical lattices. Phys. Rev. Lett. 100, 103904 (2008).
Kawabata, K. & Sato, M. Real spectra in non-hermitian topological insulators. Phys. Rev. Res. 2, 033391 (2020).
Lin, Z. et al. Unidirectional invisibility induced by \({{{{{{{\mathcal{P}}}}}}}}{{{{{{{\mathcal{T}}}}}}}}\)-symmetric periodic structures. Phys. Rev. Lett. 106, 213901 (2011).
Chong, Y., Ge, L. & Stone, A. D. \({{{{{{{\mathcal{P}}}}}}}}{{{{{{{\mathcal{T}}}}}}}}\)-symmetry breaking and laser-absorber modes in optical scattering systems. Phys. Rev. Lett. 106, 093902 (2011).
Wiersig, J. Enhancing the sensitivity of frequency and energy splitting detection by using exceptional points: application to microcavity sensors for single-particle detection. Phys. Rev. Lett. 112, 203901 (2014).
Gong, Z. et al. Topological phases of non-hermitian systems. Phys. Rev. X 8, 031079 (2018).
Porras, D. & Fernández-Lorenzo, S. Topological amplification in photonic lattices. Phys. Rev. Lett. 122, 143901 (2019).
Herviou, L., Bardarson, J. H. & Regnault, N. Defining a bulk-edge correspondence for non-hermitian hamiltonians via singular-value decomposition. Phys. Rev. A 99, 052118 (2019).
Okuma, N. & Sato, M. Hermitian zero modes protected by nonnormality: application of pseudospectra. Phys. Rev. B 102, 014203 (2020).
Daley, A. J. Quantum trajectories and open many-body quantum systems. Adv. Phys. 63, 77 (2014).
McDonald, A., Hanai, R. & Clerk, A. A. Nonequilibrium stationary states of quantum non-hermitian lattice models. Phys. Rev. B 105, 064302 (2022).
Li, J. et al. Observation of parity-time symmetry breaking transitions in a dissipative floquet system of ultracold atoms. Nat. Commun. 10, 855 (2019).
Zhang, W. et al. Observation of non-hermitian topology with nonunitary dynamics of solid-state spins. Phys. Rev. Lett. 127, 090501 (2021).
Naghiloo, M., Abbasi, M., Joglekar, Y. N. & Murch, K. Quantum state tomography across the exceptional point in a single dissipative qubit. Nat. Phys. 15, 1232 (2019).
Abbasi, M., Chen, W., Naghiloo, M., Joglekar, Y. N. & Murch, K. W. Topological quantum state control through exceptional-point proximity. Phys. Rev. Lett. 128, 160401 (2022).
Lecocq, F. et al. Nonreciprocal microwave signal processing with a field-programmable josephson amplifier. Phys. Rev. Appl. 7, 024028 (2017).
Del Pino, J., Slim, J. J. & Verhagen, E. Non-hermitian chiral phononics through optomechanically induced squeezing. Nature 606, 82 (2022).
Wang, Y.-X. & Clerk, A. Non-hermitian dynamics without dissipation in quantum systems. Phys. Rev. A 99, 063834 (2019).
Liang, Q. et al. Dynamic signatures of non-hermitian skin effect and topology in ultracold atoms. Phys. Rev. Lett. 129, 070401 (2022).
Hung, J. S. et al. Quantum simulation of the bosonic creutz ladder with a parametric cavity. Phys. Rev. Lett. 127, 100503 (2021).
Wang, K. et al. Generating arbitrary topological windings of a non-hermitian band. Science 371, 1240 (2021).
Wang, K., Dutt, A., Wojcik, C. C. & Fan, S. Topological complex-energy braiding of non-hermitian bands. Nature 598, 59–64 (2021).
Wang, Y.-P. et al. Realizing and characterizing chiral photon flow in a circuit quantum electrodynamics necklace. Sci. Rep. 5, 8352 (2015).
Roushan, P. et al. Chiral ground-state currents of interacting photons in a synthetic magnetic field. Nat. Phys. 13, 146 (2017).
McDonald, A., Pereg-Barnea, T. & Clerk, A. Phase-dependent chiral transport and effective non-hermitian dynamics in a bosonic kitaev-majorana chain. Phys. Rev. X 8, 041031 (2018).
Xiao, L. et al. Observation of non-bloch parity-time symmetry and exceptional points. Phys. Rev. Lett. 126, 230402 (2021).
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 (2020).
Weidemann, S. et al. Topological funneling of light. Science 368, 311 (2020).
Kitaev, A. Y. Unpaired majorana fermions in quantum wires. Physics-Uspekhi 44, 131 (2001).
Yao, S. & Wang, Z. Edge states and topological invariants of non-hermitian systems. Phys. Rev. Lett. 121, 086803 (2018).
Okuma, N., Kawabata, K., Shiozaki, K. & Sato, M. Topological origin of non-hermitian skin effects. Phys. Rev. Lett. 124, 086801 (2020).
Okuma, N. & Sato, M. Quantum anomaly, non-hermitian skin effects, and entanglement entropy in open systems. Phys. Rev. B 103, 085428 (2021).
Lee, G., McDonald, A. & Clerk, A. Anomalously large relaxation times in dissipative lattice models beyond the non-hermitian skin effect. Phys. Rev. B 108, 064311 (2023).
Hatano, N. & Nelson, D. R. Localization transitions in non-hermitian quantum mechanics. Phys. Rev. Lett. 77, 570 (1996).
Flynn, V. P., Cobanera, E. & Viola, L. Deconstructing effective non-hermitian dynamics in quadratic bosonic hamiltonians. New J. Phys. 22, 083004 (2020).
De Léséleuc, S. et al. Observation of a symmetry-protected topological phase of interacting bosons with Rydberg atoms. Science 365, 775 (2019).
Kawabata, K., Shiozaki, K. & Ryu, S. Many-body topology of non-hermitian systems. Phys. Rev. B 105, 165137 (2022).
Zhang, S.-B., Denner, M. M., Bzdušek, Tcv, Sentef, M. A. & Neupert, T. Symmetry breaking and spectral structure of the interacting hatano-nelson model. Phys. Rev. B 106, L121102 (2022).
Faugno, W. N. & Ozawa, T. Interaction-induced non-hermitian topological phases from a dynamical gauge field. Phys. Rev. Lett. 129, 180401 (2022).
Wang, Y.-N., You, W.-L. & Sun, G. Quantum criticality in interacting bosonic kitaev-hubbard models. Phys. Rev. A 106, 053315 (2022).
Wanjura, C. C., Brunelli, M. & Nunnenkamp, A. Topological framework for directional amplification in driven-dissipative cavity arrays. Nat. Commun. 11, 3149 (2020).
Ramos, T., García-Ripoll, J. J. & Porras, D. Topological input-output theory for directional amplification. Phys. Rev. A 103, 033513 (2021).
Huh, J., Guerreschi, G. G., Peropadre, B., McClean, J. R. & Aspuru-Guzik, A. Boson sampling for molecular vibronic spectra. Nat. Photonics 9, 615 (2015).
McDonald, A. & Clerk, A. A. Exponentially-enhanced quantum sensing with non-hermitian lattice dynamics. Nat. Commun. 11, 5382 (2020).
Kawabata, K., Numasawa, T. & Ryu, S. Entanglement phase transition induced by the non-hermitian skin effect. Phys. Rev. X 13, 021007 (2023).
Zakka-Bajjani, E. et al. Quantum superposition of a single microwave photon in two different’colour’states. Nat. Phys. 7, 599 (2011).
Chang, C. S. et al. Generating multimode entangled microwaves with a superconducting parametric cavity. Phys. Rev. Appl. 10, 044019 (2018).
Acknowledgements
The authors wish to thank B. Plourde, J.J. Nelson and M. Hutchings at Syracuse University for invaluable help in junction fabrication for the device used to produce Fig. 6. C.M.W., J.H.B., Z.S., D.D., I.N., J.S.C.H., and C.W.S.C. acknowledge the Canada First Research Excellence Fund (CFREF), NSERC of Canada, the Canadian Foundation for Innovation, the Ontario Ministry of Research and Innovation, and Industry Canada for financial support. Device fabrication was done at the University of Waterloo’s QNFCF facility. This infrastructure would not be possible without the significant contributions of CFREF-TQT, CFI, ISED, the OMRI, and Mike & Ophelia Lazaridis. A.A.C., A.M. acknowledge support from the Air Force Office of Scientific Research under Grant No. FA9550-19-1-0362, and the Simons Foundation through a Simons Investigator Award (Grant No. 669487, A.A.C.).
Author information
Authors and Affiliations
Contributions
J.H.B., D.D., and C.M.W. contributed to the experimental work. J.S.C.H. contributed to the experimental work while a student at the University of Waterloo, prior to joining AWS Center for Quantum Computing. I.N., D.D., and C.W.S.C. fabricated the devices. Z.S., A.M., and A.A.C contributed to the theoretical work.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Busnaina, J.H., Shi, Z., McDonald, A. et al. Quantum simulation of the bosonic Kitaev chain. Nat Commun 15, 3065 (2024). https://doi.org/10.1038/s41467-024-47186-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-47186-8
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.