3
introduced as an artificially engineered environment to the test qubit, but whose dynamics was not included in the
effective NH Hamiltonian dynamics. In distinct contrast, in the present system the degree of freedom of Rconstitutes
a part of the Hamiltonian, whose NH term is induced by the natural environment.
The experimental demonstration of the NH physics is performed in a circuit quantum electrodynamics architecture,
where the qubit-photon model is realized with a superconducting qubit Qand its resonator Rwith a fixed frequency
ωr/2π= 6.66 GHz (see Supplemental Material, Sec. S2 [45]). The decaying rates of Qand Rare κq≃0.07 MHz
and κf= 5 MHz, respectively. The accessible maximum frequency of Q,ωmax = 2π×6.01 GHz, is lower than ωrby
an amount much larger than their on-resonance swapping coupling gr= 2π×41 MHz (Fig. 1c). To observe the EP
physics, an ac flux is applied to Q, modulating its frequency as ωq=ω0+εcos(νt), where ω0is the mean |e⟩-|g⟩energy
difference, and εand νdenote the modulating amplitude and frequency. This modulation enables Qto interact with
Rat a preset sideband, with the photon swapping rate Ω tunable by ε(see Supplemental Material, Sec. S3 [45]).
Our experiment focuses on the single-excitation case (n= 1). Before the experiment, the system is initialized to
the ground state |g, 0⟩. The experiment starts by transforming Qfrom the ground state |g⟩to excited state |e⟩with
aπpulse, following which the parametric modulation is applied to Qto initiate the Q-Rinteraction (see Fig. S4 of
Supplementary Material for the pulse sequence). This interaction, together with the natural dissipations, realizes the
NH Hamiltonian of Eq. (1). After the modulating pulse, the Q-Rstate is measured with the assistance of an ancilla
qubit (Qa) and a bus resonator (Rb), which is coupled to both Qand Qa. The subsequent Q→Rb,Rb→Qa, and
R→Qquantum state transferring operations map the Q-Routput state to the Qa-Qsystem, whose state can be
read out by quantum state tomography (see Supplemental Material, Sec. S5 [45]).
A defining feature of our system is the conservation of the excitation number under the NH Hamiltonian. This
Hamiltonian evolves the system within the subspace {|e, 0⟩,|g, 1⟩}. A quantum jump would disrupt this conversion,
moving the system out of this subspace. This property in turn enables us to post-select the output state governed the
NH Hamiltonian simply by discarding the joint Qa-Qoutcome |g, g⟩after the state mapping. With a correction of the
quantum state distortion caused by the decoherence occurring during the state mapping, we can infer the quantum
Rabi oscillatory signal, characterized by the evolution of the joint probability |e, 0⟩, denoted as Pe,0. Fig. 2a shows
thus-obtained Pe,0as a function of the rescaled coupling ηand the Q-Rinteraction time t. The results clearly show
that both the shapes and periods of vacuum Rabi oscillations are significantly modulated by the NH term around the
EP η= 1, where incoherent dissipation is comparable to the coherent interaction. These experimental results are in
well agreement with numerical simulations (see Supplemental Material, Sec. S4 [45]).
The Rabi signal does not unambiguously reveal the system’s quantum behavior. To extract full information of
the two-qubit entangled state, it is necessary to individually measure all three Bloch vectors for each qubit, and
then correlate the results for the two qubits. The z-component can be directly measured by state readout, while
measurements of the x- and y-components require y-rotations and x-rotations before state readout, which breaks
down excitation-number conservation, and renders it impossible to distinguish individual jump events from no-jump
ones. We circumvent this problem by first reconstructing the two-qubit density matrix using all the measurement
outcomes, and then discarding the matrix elements associated with |g, g⟩(see Supplemental Material, Sec. S5 [45]).
This effectively post-casts the two-qubit state to the subspace {|e, g⟩,|g, e⟩}. We note that such a technique is in
sharp contrast with the conventional post-selection method, where some auxiliary degree of freedom (e.g., propagation
direction of a photon) enables reconstruction of the relevant conditional output state, but which is unavailable in the
NH qubit-photon system. With a proper correction for the state mapping error, we obtain the Q-Routput state
governed by the non-Hermitian Hamiltonian. In Fig. 2b, we present the resulting Q-Rconcurrence, as a function of
ηand t. To show the entanglement behaviors more clearly, we present the concurrence evolutions for η= 5 and 0.5
in Fig. 2c and d, respectively. The results demonstrate that the entanglement exhibits distinct evolution patterns in
the regimes above and below the EP.
To reveal the close relation between the exceptional entanglement behavior and the EP, we infer the eigenvalues and
eigenstates of the NH Hamiltonian from the output states, measured for different evolution times (see Supplemental
Material, Sec. S7 [45]). The concurrences associated with the two eigenstates, obtained for different values of η, are
shown in Fig. 3a and b (dots), respectively. As expected, each of these entanglements exhibits a linear scaling below
the EP, but is saturated at the EP, and no longer depends upon ηafter crossing the EP. The singular features around
the EP are manifested by their derivatives to η, which are shown in the insets. The discontinuity of these derivatives
indicates the occurrence of an entanglement transition at the EP. Such a nonclassical behavior, as a consequence of
the competition between the coherent coupling and incoherent dissipation, represents a unique character of strongly
correlated NH quantum systems, but has not been reported so far. These results demonstrate a longitudinal merging
of the two entangled eigenstates at the EP (left panel of Fig. 2b). Accompanying this is the onset of the transverse
splitting, which can be characterized by the relative phase difference between |Φ1,±⟩, defined as φ=φ+−φ−, with
φ±representing the relative phases between |g, 1⟩and |e, 0⟩in the eigenstates |Φ1,±⟩. Such a phase difference, inferred