--- title: Two-Qubit Non-Hermitian XY Model url: https://www.emergentmind.com/topics/two-qubit-non-hermitian-xy-model type: topic --- # Two-Qubit Non-Hermitian XY Model Searching arXiv for the specified papers and closely related two-qubit non-Hermitian XY models. First search: measurement-induced non-Hermitian entanglement generation and two-qubit XY realizations. Second search: RT/PT-symmetric non-Hermitian XY dimers, geometric phases, and thermodynamic applications. The two-qubit non-Hermitian XY model denotes a class of spin-\(\tfrac12\) dimer Hamiltonians in which the XY exchange sector is rendered non-Hermitian by complex anisotropy, imaginary local \(z\)-fields, conditional decay terms, or measurement-induced postselection. In the recent literature, this label covers measurement-induced effective dimers with complex exchange couplings and imaginary local fields, RT-symmetric anisotropic XY dimers obtained as the \(N=2\) reduction of non-Hermitian chains, PT-symmetric dimers with a staggered imaginary field, and heralded no-jump XY dimers generated by spontaneous decay. Across these realizations, the common technical structures are \(2\times2\) invariant subspaces, biorthogonal spectral decompositions, exceptional points in specific parameter regimes, and Bell-state selection under normalized non-unitary dynamics [2009.10004, 1210.5613, 2203.05371, 1402.6700]. ## 1. Canonical Hamiltonian forms In the two-qubit setting, the non-Hermitian XY label does not refer to a single universal Hamiltonian. Instead, it refers to several analytically tractable dimer models that share an XY-type spin-exchange core and a non-Hermitian deformation. A general measurement-induced form is \[ H_{\mathrm{XY}} = J_x \sigma_1^x \sigma_2^x + J_y \sigma_1^y \sigma_2^y + h_1^z \sigma_1^z + h_2^z \sigma_2^z + J_z \sigma_1^z \sigma_2^z, \] with complex \(J_x,J_y\) and, in general, complex local fields \(h_1^z,h_2^z\). In the symmetric ancilla-induced construction, the identifications are \[ J_x=J_y=\gamma_{xy}-i\tau g_{xy}^2,\qquad J_z=\gamma_z,\qquad h_1^z=h_2^z=g_z+i\tau g_{xy}^2. \] A distinct RT-symmetric dimer arises from the \(N=2\) limit of a non-Hermitian anisotropic XY chain with periodic boundary conditions: \[ H_2=J\Big[(1+i\gamma)\,\sigma_1^x\sigma_2^x+(1-i\gamma)\,\sigma_1^y\sigma_2^y\Big] +J\lambda\left(\sigma_1^z+\sigma_2^z\right). \] Here the non-Hermiticity is encoded directly in the complex anisotropy, with \(J_x=J(1+i\gamma)/2\) and \(J_y=J(1-i\gamma)/2\). A PT-symmetric formulation uses a real transverse field and a staggered imaginary field, \[ H = -J\big[(1+\gamma)\,\sigma_1^x\sigma_2^x+(1-\gamma)\,\sigma_1^y\sigma_2^y\big] -h\,(\sigma_1^z+\sigma_2^z)+i\eta(-\sigma_1^z+\sigma_2^z), \] or, equivalently, in reduced parameters \(h_0=h/J\), \(\eta_0=\eta/J\), \[ H = -J\left[(1+\gamma)\,\sigma_1^x\sigma_2^x + (1-\gamma)\,\sigma_1^y\sigma_2^y\right] -J h_0 \left(\sigma_1^z + \sigma_2^z\right) + i J \eta_0 \left(-\sigma_1^z + \sigma_2^z\right). \] A further realization is the heralded no-jump XY dimer. There the effective Hamiltonian has the generic form \[ H_{\mathrm{eff}}=H-\frac{i\gamma}{2}\sum_n P_{\uparrow,n} = \sum_n \Big[2J(\sigma_n^+\sigma^-_{n+1}+\sigma_n^-\sigma^+_{n+1}) +2J'(\sigma_n^+\sigma_{n+1}^+ + \sigma_n^-\sigma_{n+1}^-) -\frac{i\gamma}{4}(\sigma_n^z+1)\Big], \] so the non-Hermiticity arises from conditioning on no fluorescence events rather than from a static gain-loss term in a closed model. In all of these formulations, the computational basis \(\{|00\rangle,|01\rangle,|10\rangle,|11\rangle\}\) or the equivalent spin basis \(\{|\uparrow\uparrow\rangle,|\uparrow\downarrow\rangle,|\downarrow\uparrow\rangle,|\downarrow\downarrow\rangle\}\) yields a block structure that makes the two-qubit problem exactly solvable [2009.10004, 1210.5613, 2203.05371, 2606.22405, 1402.6700]. ## 2. Measurement-induced construction from system–ancilla dynamics A physically explicit route to a two-qubit non-Hermitian XY model starts from a Hermitian three-spin Hamiltonian and repeated measurements on an ancilla spin. The system–ancilla unitary is \[ U_{SA}(\tau)=e^{-iH\tau}, \] and a projective measurement on the ancilla with selected outcome \(|0_A\rangle\) defines the Kraus operator \[ K(\tau)=\langle 0_A|U_{SA}(\tau)|0_A\rangle. \] After \(n\) successful repetitions, the non-normalized conditional state is \[ \rho_S^c(n\tau)=K^n(\tau)\,\rho_S(0)\,\big(K^n(\tau)\big)^\dagger. \] In the stroboscopic limit, \[ K(\tau)=\exp\!\left(-i\tau H_0^{\mathrm S}-\frac{\tau^2}{2}\Gamma_{\mathrm S}+o(\tau^2)\right), \] with \[ H_0^{\mathrm S}=\langle 0_A|H|0_A\rangle,\qquad \Gamma_{\mathrm S}=\langle 0_A|H^2|0_A\rangle-\big(H_0^{\mathrm S}\big)^2\ge 0, \] so that repeated measurements generate the effective non-Hermitian Hamiltonian \[ H_{\mathrm{eff}}=H_0^{\mathrm S}-\frac{i\tau}{2}\Gamma_{\mathrm S}. \] For the symmetric three-spin model with pairwise XXZ-like couplings, \[ H = \gamma_{xy} \left( \sigma_1^x \sigma_2^x + \sigma_1^y \sigma_2^y \right) + \gamma_z \sigma_1^z \sigma_2^z + g_{xy} \left( \sigma_1^x \sigma_3^x + \sigma_1^y \sigma_3^y \right) + g_z \sigma_1^z \sigma_3^z + g_{xy} \left( \sigma_2^x \sigma_3^x + \sigma_2^y \sigma_3^y \right) + g_z \sigma_2^z \sigma_3^z, \] preparing and repeatedly measuring the ancilla in the \(z\)-basis with outcome \(0\) yields, up to an identity shift, \[ H'_{\mathrm{eff}}= \gamma_{xy}\left(\sigma_1^x\sigma_2^x+\sigma_1^y\sigma_2^y\right) +\gamma_z\sigma_1^z\sigma_2^z +g_z(\sigma_1^z+\sigma_2^z) -i\tau g_{xy}^2\left(\sigma_1^x\sigma_2^x+\sigma_1^y\sigma_2^y-\sigma_1^z-\sigma_2^z\right). \] The imaginary terms arise from the positive semidefinite decay operator \(\Gamma_{\mathrm S}\) generated by ancilla flips due to the XY couplings \(g_{xy}\). Dropping the identity term affects the postselection probability but not the normalized conditional dynamics. The anisotropic generalization starts from \[ \tilde{H} = \sum_{\mu=x,y,z} \gamma_\mu \sigma_1^\mu \sigma_2^\mu + \sum_{\mu=x,y,z} \alpha_\mu \sigma_1^\mu \sigma_3^\mu + \sum_{\mu=x,y,z} \beta_\mu \sigma_2^\mu \sigma_3^\mu, \] which induces \[ \tilde{H}_{\mathrm{eff}}= (\alpha_z+i\tau\alpha_x\alpha_y)\sigma_1^z +(\beta_z+i\tau\beta_x\beta_y)\sigma_2^z +(\gamma_z-i\tau\alpha_z\beta_z)\sigma_1^z\sigma_2^z +(\gamma_x-i\tau\alpha_x\beta_x)\sigma_1^x\sigma_2^x +(\gamma_y-i\tau\alpha_y\beta_y)\sigma_1^y\sigma_2^y. \] This yields an anisotropic XY-like dimer with complex couplings and complex local \(z\)-fields. A key structural point is that XX/YY-type system–ancilla couplings together with ancilla measurements in \(\sigma_3^z\) generate precisely the imaginary XY exchange terms and imaginary \(z\)-fields characteristic of the non-Hermitian two-qubit XY model [2009.10004]. ## 3. Symmetry classes, invariant blocks, and exceptional points Two-qubit non-Hermitian XY dimers split naturally into symmetry sectors. In the RT-symmetric anisotropic dimer, the operator \[ \mathcal R=\exp\!\left[-i\frac{\pi}{4}(\sigma_1^z+\sigma_2^z)\right],\qquad \mathcal T=K \] satisfies \([\mathcal{RT},H_2]=0\), although the Hamiltonian is not PT-symmetric with the usual lattice parity. In the computational basis, \(H_2\) is block diagonal: \[ H_2= \begin{pmatrix} 2J\lambda & 0 & 0 & 2Ji\gamma\\ 0 & 0 & 2J & 0\\ 0 & 2J & 0 & 0\\ 2Ji\gamma & 0 & 0 & -2J\lambda \end{pmatrix}. \] The single-excitation block \(\{|01\rangle,|10\rangle\}\) has real eigenvalues \(\pm 2J\) with eigenvectors \((|01\rangle\pm|10\rangle)/\sqrt2\). The zero/double-excitation block \(\{|00\rangle,|11\rangle\}\) has \[ E_\pm=\pm 2J\sqrt{\lambda^2-\gamma^2}, \] so the unbroken RT-symmetric regime is \(\lambda^2\ge\gamma^2\), the broken regime is \(\lambda^2<\gamma^2\), and the exceptional point occurs at \(\lambda^2=\gamma^2\), where the block becomes a Jordan block. At \(\lambda=\gamma>0\), the coalesced eigenvector is \[ |\phi_{\mathrm{EP}}\rangle\propto |00\rangle-i|11\rangle. \] The generalized RT-symmetric chain with deformation parameters \(G(\xi)\) and \(A(\xi)\) reduces at \(N=2\) to a pair sector with \[ E_\pm^{(\mathrm{pair})}=\pm 2J\sqrt{\lambda^2-c^2},\qquad c=|G(\xi)|\sqrt{|A(\xi)|}, \] under periodic boundary conditions. The exceptional line is \(\lambda=\pm c\), and the corresponding unbroken region is \(\lambda^2\ge c^2\). In this formulation, RT symmetry is generated by \[ R=\exp\Big(-\frac{i\phi}{2}\sum_{j=1}^N \sigma_j^z\Big),\qquad [RT,H(\xi)]=0, \] with \(\phi=\arg G\). In the PT-symmetric staggered-field dimer, parity is the qubit-swap operator and time reversal is complex conjugation. The Hamiltonian decomposes exactly as \[ H=H_{\mathrm{work}}\oplus H_{\mathrm{idle}}, \] with \[ H_{\mathrm{work}}= \begin{pmatrix} -2Jh_0 & -2J\gamma\\ -2J\gamma & 2Jh_0 \end{pmatrix}, \qquad H_{\mathrm{idle}}= \begin{pmatrix} -2iJ\eta_0 & -2J\\ -2J & 2iJ\eta_0 \end{pmatrix}. \] Its spectrum is \[ E_{1,2}=\mp 2J\sqrt{h_0^2+\gamma^2},\qquad E_{3,4}=\mp 2J\sqrt{1-\eta_0^2}. \] The working pair is always real. The idle pair is real for \(|\eta_0|<1\), degenerate at zero for \(|\eta_0|=1\), and forms a complex-conjugate pair for \(|\eta_0|>1\). Thus \(|\eta_0|=1\) is the PT exceptional point. The heralded no-jump dimer provides a different non-Hermitian mechanism. For \(N=2\) with periodic boundary conditions, \[ H_{\mathrm{eff}}= \begin{pmatrix} -i\gamma & 0 & 0 & 4J'\\ 0 & -i\gamma/2 & 4J & 0\\ 0 & 4J & -i\gamma/2 & 0\\ 4J' & 0 & 0 & 0 \end{pmatrix}. \] The pairing block \(\{|{\uparrow\uparrow}\rangle,|{\downarrow\downarrow}\rangle\}\) has \[ \lambda_1^\pm=-i\frac{\gamma}{2}\pm \frac12\sqrt{64J'^2-\gamma^2}, \] and exhibits an exceptional point at \[ |J'|=\frac{\gamma}{8}. \] By contrast, the flip-flop block has \[ \lambda_2^\pm=-i\frac{\gamma}{2}\pm 4J, \] and its level meeting at \(J=0\) is not an exceptional point because the block remains diagonalizable. One recurring misconception is therefore avoided explicitly in this literature: eigenvalue degeneracy alone is not sufficient for an exceptional point; non-diagonalizability is essential [1210.5613, 1308.4057, 2203.05371, 2606.22405, 1402.6700]. ## 4. Entanglement generation, Bell sectors, and long-time selection In the measurement-induced symmetric dimer, the invariant \(|01\rangle/|10\rangle\) block is \[ H_{\mathrm{block}}=-\gamma_z I_2+\left(2\gamma_{xy}-i2\tau g_{xy}^2\right)\sigma_x. \] Its right eigenvalues and right eigenvectors are \[ E_\pm=-\gamma_z\pm \left(2\gamma_{xy}-i2\tau g_{xy}^2\right),\qquad |\psi_\pm^R\rangle=\frac{|01\rangle\pm|10\rangle}{\sqrt2}\equiv |\Psi^\pm\rangle. \] Because \[ \mathrm{Im}\,E_+=-2\tau g_{xy}^2<0,\qquad \mathrm{Im}\,E_-=+2\tau g_{xy}^2>0, \] postselection-normalized dynamics selects the growing eigenmode \(|\Psi^-\rangle\), so \[ |\psi_\infty\rangle=|\Psi^-\rangle,\qquad \varrho_S^c(t\to\infty)=|\Psi^-\rangle\langle\Psi^-|. \] For an initial state \(|01\rangle\), the concurrence evolves as \[ C(t)=\sqrt{\sin^2(2\gamma t)+\cos^2(2\gamma t)\tanh^2(2gt)}, \] with \(g\equiv 2\tau g_{xy}^2\) and \(\gamma\equiv 2\gamma_{xy}\), so \(\lim_{t\to\infty}C(t)=1\). In the anisotropic construction, suitable choices of \(\alpha_\mu,\beta_\mu,\gamma_\mu,\tau\) can instead select \[ |\psi_\infty\rangle=|\Phi^+\rangle=\frac{|00\rangle+|11\rangle}{\sqrt2}. \] This asymptotic Bell-state purification has a direct postselection cost. The success probability is \[ p(t)=\mathrm{Tr}\,\rho_S^c(t),\qquad \frac{dp(t)}{dt}=-\tau\,\mathrm{Tr}\!\left[\Gamma_{\mathrm S}\rho_S^c(t)\right]\le 0. \] Increasing \(g_{xy}\) while keeping \(g_{xy}\tau\ll 1\) increases the selection rate \(g=2\tau g_{xy}^2\) but decreases \(p(t)\) more rapidly. In the Hermitian limit \(\tau=0\) or \(g_{xy}=0\), the same \(|01\rangle/|10\rangle\) sector undergoes only unitary oscillations with \[ C_{\mathrm{Herm}}(t)=|\sin(2\gamma t)|, \] so no stationary Bell state emerges. The EP-based Bell mechanism in the open-chain single-excitation model is structurally different. For \(N=2\), the relevant effective Hamiltonian is \[ H_{\mathrm{eq}}^{(N=2)}= \begin{pmatrix} V+i\gamma & 1\\ 1 & V-i\gamma \end{pmatrix}, \] with eigenvalues \[ \varepsilon_\pm=V\pm \sqrt{1-\gamma^2}. \] At \(|\gamma|=1\), the exceptional point produces the coalescing eigenvector \[ |\psi_{\mathrm{EP}}\rangle\propto |01\rangle-i|10\rangle\quad (\gamma=+1), \] or \[ |\psi_{\mathrm{EP}}\rangle\propto |01\rangle+i|10\rangle\quad (\gamma=-1). \] These states are Bell states up to local phase rotations about \(z\). In the PT-broken regime \(|\gamma|>1\), normalized non-unitary evolution converges exponentially to the eigenvector with positive imaginary part, with convergence rate \(\kappa=\sqrt{\gamma^2-1}\). The staggered imaginary-field dimer yields a complementary entanglement diagnostic. In the odd sector, \[ C^{(3)}(\eta_0)= \begin{cases} 1,& |\eta_0|\le 1,\\[3pt] 1/|\eta_0|,& |\eta_0|>1, \end{cases} \] so the concurrence is exactly unity in the PT-unbroken region and develops a cusp at the exceptional point. In the even sector, \[ C_g^{(1)}=\frac{|\gamma|}{\sqrt{h_0^2+\gamma^2}}, \] which is independent of \(\eta_0\). This indicates that maximal entanglement is not tied uniquely to exceptionality: in some formulations it is selected asymptotically by imaginary spectral ordering, while in others it already characterizes an entire unbroken symmetry sector [2009.10004, 1503.05991, 2203.05371]. ## 5. Geometric phase, thermal structure, and quantum-engine applications In the generalized RT-symmetric non-Hermitian XY model, the geometric phase is formulated biorthogonally through \[ A_i(\xi)=i\,\langle\!\langle G(\xi)|\,\partial_{\eta_i}G(\xi)\rangle,\qquad F_{ij}(\xi)=\partial_{\eta_i}A_j-\partial_{\eta_j}A_i. \] For the three-parameter example in cylindrical coordinates \((\rho,\phi,\lambda)\), the curvature has only an azimuthal component, \[ Q_\phi(\xi)=4J^3\,\rho\sum_{00\) for \(\eta_0>\eta_c\). In the same example, the efficiency rises from about \(3.5\%\) at \(\eta_0=0.66\) to \(42.7\%\) at \(\eta_0=0.9\), and the limit \(\eta_0\to 1^-\) yields \[ \eta \to 1-\frac{A_C}{A_H} =1-\frac{\sqrt{h_{C,0}^{2}+\gamma^2}}{\sqrt{h_{H,0}^{2}+\gamma^2}}, \] which is the Otto bound and remains below the Carnot bound. This suggests a broader significance of the two-qubit non-Hermitian XY model: it is not only a minimal EP-bearing spin dimer, but also a microscopic platform in which non-Hermiticity controls geometry, thermal occupation, and thermodynamic mode switching [1308.4057, 2606.22405]. ## 6. Inverse engineering, experimental platforms, and limitations A general inverse-construction recipe makes the two-qubit non-Hermitian XY model physically simulable from a Hermitian enlarged system. Given a target finite-dimensional non-Hermitian Hamiltonian \(H_{\mathrm{target}}\) on the system \(S\), one decomposes it into Hermitian and anti-Hermitian parts, chooses \(\tau\) such that the maximum Bohr frequency \(f\) of \(i(H_{\mathrm{target}}-H_{\mathrm{target}}^\dagger)\) satisfies \(f\tau\ll 1\), defines \[ c=\max\!\left(0,-\lambda_{\min}\left[\frac{i}{\tau}\left(H_{\mathrm{target}}-H_{\mathrm{target}}^\dagger\right)\right]\right), \] and constructs the Hermitian Hamiltonian on \(S\otimes A\), \[ H= \frac12\left(H_{\mathrm{target}}+H_{\mathrm{target}}^\dagger\right)\otimes |0_A\rangle\langle 0_A| +\sqrt{\,cI+\frac{i}{\tau}\left(H_{\mathrm{target}}-H_{\mathrm{target}}^\dagger\right)} \otimes \left(|0_A\rangle\langle 1_A|+|1_A\rangle\langle 0_A|\right). \] Repeatedly measuring the ancilla in the \(z\)-basis and postselecting outcome \(0\) reproduces \(H_{\mathrm{target}}\) in the stroboscopic limit. For an anisotropic non-Hermitian XY target, the mapping is \[ \gamma_\mu=\mathrm{Re}\,J_\mu,\qquad \mathrm{Im}\,J_\mu=-\tau\alpha_\mu\beta_\mu, \] together with \[ \alpha_z=\mathrm{Re}\,h_1^z,\qquad \beta_z=\mathrm{Re}\,h_2^z,\qquad \mathrm{Im}\,h_1^z=\tau\alpha_x\alpha_y,\qquad \mathrm{Im}\,h_2^z=\tau\beta_x\beta_y. \] Several experimental platforms are identified as feasible. Trapped ions, superconducting qubits, and photonic qubits can implement spin-\(\tfrac12\) interactions together with fast, high-fidelity projective measurements. For the PT-symmetric Otto-engine dimer, trapped-ion realizations use two hyperfine or Zeeman levels per ion, Mølmer–Sørensen interactions to engineer anisotropic XY couplings, AC-Stark shifts for the real field, and quantum-trajectory postselection to generate the effective staggered imaginary field. For the heralded no-jump model, a three-level scheme with \(|\uparrow\rangle\), \(|\downarrow\rangle\), and an auxiliary state \(|a\rangle\) implements the effective non-Hermitian Hamiltonian by conditioning on the absence of fluorescence; trapped ions, cavity QED arrays, and Rydberg atoms in optical lattices are listed as candidate platforms. The main limitations are formulation-dependent but technically sharp. In the measurement-induced construction, the stroboscopic approximation requires \(g_{xy}\tau\ll 1\); larger \(\tau\) or explicitly time-dependent \(H\) invalidates the simple \(H_{\mathrm{eff}}\) form. Measurement errors and finite detection efficiency reduce postselection success probability and can perturb \(\Gamma_{\mathrm S}\), although moderate imperfections primarily increase postselection cost rather than alter the selected Bell state. Decoherence during the inter-measurement intervals adds Lindbladian terms not captured by the idealized effective Hamiltonian. In the PT-symmetric engine model, one must remain within \(|\eta_0|<1\) for a real spectrum, and adiabaticity degrades as the idle gap \(B\to 0\) near the exceptional point; the explicit recommendation is to truncate experimentally at about \(\eta_0\lesssim 0.95\). These caveats do not diminish the central result: two-qubit non-Hermitian XY models constitute a controlled and experimentally grounded setting in which non-Hermitian spectral theory, postselected dynamics, and entanglement generation can be analyzed exactly [2009.10004, 2606.22405, 1402.6700]. Source: https://www.emergentmind.com/topics/two-qubit-non-hermitian-xy-model