## Bosonic Holes in Quadratic Bosonic Systems

Jia-Ming Hu, 1 Bo Wang, 1 and Ze-Liang Xiang 1, ∗

1 School of Physics, Sun Yat-sen University, Guangzhou 510275, China.

The concept of electron holes plays a significant role in condensed matter physics. Here we develop the concept of bosonic holes, which exhibit negative particle excitations, in quadratic bosonic systems. Unlike electron holes, the Fock states of bosonic holes are biorthogonal, and their excitation can be interpreted as removing particles from a mean-particle number with a mean-field background. Furthermore, we find that quadratic bosonic Hamiltonians related by non-unitary and local particlehole (PH) transformation possess the same locality structure and spectral properties in different spaces, reflecting the PH duality. Based on this, we study the generation of PH entanglement in two-mode cases and the PH Aharonov-Bohm effect in the three-mode case, which results in a PH chiral flow with time-reversal symmetry breaking. Our findings provide a new way to understand and explore unusual physical phenomena in particle-non-conserving and non-Hermitian systems.

In quadratic systems, Hamiltonians with pairing interactions feature particle non-conservation and play an essential role in various interesting quantum phenomena. In condensed matter physics, they are the basic form of weakly interacting Bose gases and fermionic superconductors in the mean-field approximation, giving rise to superfluidity [1] and superconductivity [2]. In quantum optics, the bosonic pairing is known as the parametric amplifier interaction [3, 4] and is regarded as a key resource for generating squeezing and entanglement [5-8].

The excitation spectrum of quadratic systems can be well described by a Bogoliubov-de-Gennes (BdG) Hamiltonian [9-11], in which the hole degrees of freedom appear. In the fermionic systems, the quasi-particles in a fermionic superconductor are superpositions of electrons and holes [12, 13]. The absence of the Pauli exclusion principle implies that bosonic holes can significantly differ from electron holes. For instance, bosonic pairing interactions in the BdG Hamiltonian can be regarded as dissipative particle-hole conversions [14]. Therefore, pairing interactions can be used as an alternative way to introduce non-Hermiticity in bosonic systems besides coupling with the dissipative baths [15, 16]. This has garnered growing attention in studying the quadratic bosonic topological models [5-8, 10, 17], in which a set of remarkable properties are predicted, such as quadraturedependent chiral transport and amplification [18-20], and unstable edge modes [21-25]. In particular, the nonHermitian Aharonov-Bohm (AB) interference effect in bosonic PH loops has been confirmed in experiments [14]. However, the bosonic hole, which is so important as a new degree of freedom, still lacks a clear physical picture and a rigorous theory to depict it.

In this letter, we characterize bosonic holes and reveal the PH duality in quadratic bosonic systems (QBSs). We find that bosonic holes , forming one-half of the basis of the bosonic BdG Hamiltonian, are pseudo-bosonic modes. A non-unitary PH transformation can be used to obtain the biorthogonal Fock states of bosonic holes, which are occupied by negative particle numbers due to the constraint of the bosonic commutation relation. Nev- ertheless, they may be observable when the excitation of holes is relative to a mean-particle number or macroscopic occupation of particles. This mean-particle number plays a role as the 'Fermi surface' and can be formed by displacing the particle or hole mode with a driving field, whose frequency can be analogous to the 'Fermi level' (see Fig. 1). In this picture, effective Hermitian and non-Hermitian QBHs can be related by the local PH transformation, while the locality structure and spectral properties are preserved, such called the PH duality. The dual QBSs show the same physics in the spaces with a single PH inversion, where their biorthogonal eigenstates can be well interpreted as the superposition states involving holes, allowing the study of novel physical phenomena in the PH space. In this setting, we explore the generation of PH entanglement in the two-mode cases, including the PH entangled ground states and the sensitivity of PH entanglement to initial states. Moreover, the Hermitian AB effect in PH loops is investigated, and we find the chiral flow in the 'single-hole trimer' with time-reversal symmetry breaking.

Hole degrees of freedom in QBHs. -We start with a general Hermitian quadratic Hamiltonian ˆ H = ∑ i,j ˆ a † i A ij ˆ a j +1 2(ˆ / a † i B ij ˆ a † j +ˆ a i B ∗ ij ˆ ). a j The Hermiticity of ˆ H imposes the condition A ij = A ∗ ji in the particleconserving terms, and the bosonic commutation relation allows us to take B ij = B ji in the particle-nonconserving pairing interactions. Defining the Nambu array ˆ = (ˆ α a 1 · · · ˆ a N , a ˆ † 1 · · · ˆ a † N ) T , the Hamiltonian can be rewritten as ˆ H = 1 2ˆ / α Hα † ˆ, while the Heisenberg equation is given by i∂ α t ˆ = H ˆ, where α

$$\underset { \underset { \varepsilon _ { - } } { \mathcal { H } } } { \mathfrak { r } _ { - } } \quad H = \left ( \begin{array} { c c } \mathcal { A } & \mathcal { B } \\ \mathcal { B } ^ { * } & \mathcal { A } ^ { * } \end{array} \right ), \quad \mathcal { H } = \tau _ { 3 } H = \left ( \begin{array} { c c } \mathcal { A } & \mathcal { B } \\ - \mathcal { B } ^ { * } & - \mathcal { A } ^ { * } \end{array} \right ), \ \ ( 1 )$$

and τ i = σ i ⊗ I N from the Pauli matrices σ i ( i = 1 2 3). , , Instead of the Hermitian H , the Heisenberg dynamical matrix H with pseudo Hermiticity, τ 3 H τ 3 = H † , inherently respects PH symmetry, τ 1 H τ 1 = -H ∗ . Furthermore, the eigenvalues of H are consistent with the system spectrum, so that it can be treated as the bosonic BdG

̸

Hamiltonian [10, 11]. Substituting the BdG Hamiltonian into the second-quantized Hamiltonian ˆ H = 1 2ˆ ¯ / α † H ˆ α with ˆ ¯ α † = ˆ α τ † 3 = (ˆ a † 1 · · · ˆ a N † , -ˆ a 1 · · · -ˆ a N ), the annihilation operators of bosonic holes can be regarded as ˆ a † i in ˆ [14, 20], while their creation operators are given α by -ˆ a i in ˆ , ¯ α † instead of ˆ a i in ˆ , α † to fulfill the bosonic commutation relation, [ˆ a , † i -ˆ ] a j = [ ˆ h , h i ˆ ¯ † j ] = δ ij . Similar to the situation in electron systems, bosonic particle and hole operators are redundant, i.e., the bosonic holes are equivalent to the bosonic particles with opposite frequencies. However, bosonic holes feature negative particle excitations due to the bosonic commutation relation, ˆ ˆ a a † i i + ˆ ¯ h h † i ˆ i = -1, and they are pseudo-bosonic modes because ¯ ˆ h † i = ˆ h † i [26-28]. These introduce a redundancy in the Nambu notation with ˆ ¯ α † = ˆ α τ τ T 1 3 , leading to the combination of PH symmetry and pseudo-Hermiticity of H , i.e., τ 3 τ 1 H τ 1 τ 3 = -H T .

Hence, the eigenvalues of H can be divided into two groups { ϵ n , -ϵ n } . Like particles and holes, the eigenmodes with the frequencies ϵ n are equivalent to those with the frequencies -ϵ n . However, they cannot simultaneously hold a positive excitation on the eigenstates of ˆ . H For this reason, ˆ H can be rewritten in terms of two sets of eigenmodes with two sets of eigenstates. In the real spectrum regime ( ϵ n ∈ R ),

$$\hat { H } = \sum _ { n = 1 } ^ { N } \epsilon _ { n } ( \hat { \psi } _ { n } ^ { \dagger } \hat { \psi } _ { n } + \frac { 1 } { 2 } ) = \sum _ { n = 1 } ^ { N } - \epsilon _ { n } ( \hat { \phi } _ { n } ^ { \dagger } \hat { \phi } _ { n } + \frac { 1 } { 2 } ), \quad ( 2 ) \quad \text{PH tta} \\ \text{frequ} \\ \text{Fi} \sigma$$

where ˆ ψ † n = ˆ α τ † 3 | ψ n ⟩ = ˆ ϕ n , ˆ ¯ ϕ † n = -ˆ ψ n , and | ψ n ⟩ are the eigenvectors of H with the eigenvalues ϵ n . It is clear that ˆ ψ n are normal bosonic modes only composed of particles, while ˆ ϕ n are pseudo-bosonic modes only composed of bosonic holes due to PH symmetry, thus called quasi-particles and quasi-holes. In the complex spectrum regime ( ϵ n / ∈ R ),

$$\hat { H } = \sum _ { n = 1 } ^ { N } \epsilon _ { n } ( \hat { \psi } _ { n } ^ { \dagger } \hat { \psi } _ { n ^ { * } } + \frac { 1 } { 2 } ) = \sum _ { n = 1 } ^ { N } - \epsilon _ { n } ( \hat { \bar { \phi } } _ { n ^ { * } } ^ { \dagger } \hat { \phi } _ { n } + \frac { 1 } { 2 } ), \ \ ( 3 ) \quad \begin{array} { c } \text{states} \\ \hat { \Omega } ^ { - 1 } | \gamma \\ \hat { \iota }, \i \cap \_ \_ \_ \end{array} \end{array}$$

where ˆ ψ n ∗ = ⟨ ψ n ∗ | τ 3 ˆ = α -ˆ ¯ ϕ † n ∗ and | ψ n ∗ ⟩ are the eigenvectors of H with the eigenvalues ϵ ∗ n . In this case, both sets of eigenmodes are pseudo-bosonic modes and, therefore, composed of particles and holes. The composition of eigenmodes in terms of particles and holes can be examined by the expectation of each particle number operator in eigenstates, where a positive (negative) value indicates the particle (hole) components and the complex value indicates PH mixed ones [29]. The spectrum of these eigenmodes with hole components has been observed experimentally [14, 30], confirming the presence of hole degrees of freedom in QBHs. Then, let us turn to find a rigorous formalism of bosonic holes to characterize the negative particle occupation space.

FIG. 1. Bosonic excitation spectrum in QBHs with respect to the frequency of a driving field ω L ('Fermi level'), where ± ∆ are the relative frequencies of particles and holes, respectively, with ω a being the intrinsic frequency of the boson.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01059v1-figure-751c6784-f6dc-4749-9876-a70ab0b15449.jpg)

Particle-hole transformation. -To obtain the Fock states of bosonic holes, we here introduce the PH transformation with a non-unitary but Hermitian operator,

$$\hat { \Omega } = e ^ { i \frac { \pi } { 4 } ( e ^ { - i \theta } \hat { a } ^ { 2 } - e ^ { i \theta } \hat { a } ^ { \dagger 2 } ) },$$

such that the bosonic operators ˆ and ˆ a a † can now be transformed into the hole operators, Ω ˆ -1 ˆ ˆ a Ω = -ia e ˆ † iθ = ˆ h and Ω ˆ -1 ˆ Ω a † ˆ = -iae ˆ -iθ = ˆ ¯ h † , where θ is an arbitrary gauge phase. Different from PH symmetry, the PH transformation maps particles to holes with the same frequency, i.e., Ω ˆ -1 ∆(ˆ ˆ + 1 a a † / 2) ˆ Ω = ∆( ˆ ¯ ˆ + 1 h h † / 2) (see Fig. 1), which plays a role as the charge-conjugation transformation in fermionic systems [31].

Then, the vacuum state of the hole can be defined as | 0 ⟩ h = ˆ Ω -1 | 0 ⟩ a to make -ia e ˆ † iθ | 0 ⟩ h = ˆ 0 h | ⟩ h = 0. Notably, this state is biorthogonal because of the nonunitarity of Ω. ˆ We take ˆ and ¯ h ˆ h † as the annihilation and creation operators of the right vector, while their transposed conjugations ˆ h † and ˆ ¯ as the creation and annihih lation operators of the left vector. The biorthogonal Fock states of the hole are given by | n ⟩ h = ( ˆ ¯ h † ) n / √ n ! 0 | ⟩ h = ˆ Ω -1 | n ⟩ a and | m ⟩ ¯ h = ( ˆ h † ) m / √ m ! 0 | ⟩ ¯ h = ˆ Ω | m ⟩ a , where ˆ ¯ 0 h | ⟩ ¯ h = 0 and ¯ h ⟨ mn | ⟩ h = δ mn . Therefore, we obtain the hole Fock states that occupy by negative particle number ¯ h ⟨ n a a n | ˆ ˆ † | ⟩ h = -( n +1).

Such negative particle occupation can be understood by introducing the mean-field theory, wherein the positive and negative particle excitation is relative to the mean particle number, which plays a role as the 'Fermi surface.' We bring it with a displacement operator ˆ (¯) = D a e ¯ˆ aa † -¯ a a ∗ ˆ . In the displaced space, n bosonic holes can be regarded as removing n + 1 particles from the mean particle number ¯ = n | ¯ a | 2 , i.e.,

$$_ { \bar { h } } \langle n | \hat { D } ^ { - 1 } ( \bar { a } ) \hat { a } ^ { \dagger } \hat { a } \hat { D } ( \bar { a } ) | n \rangle _ { h } = | \bar { a } | ^ { 2 } - ( n + 1 ). \quad ( 5 )$$

This displacement can be realized by considering a system driven by a coherent field, described by the Hamil-

tonian ˆ H D = ∆ˆ ˆ + a a † λ a (ˆ † + ˆ). a Here, ∆ denotes the detuning of the system mode and driving field, and λ represents the driving strength. The Lindblad master equation of the system can be obtained,

$$\frac { d \hat { \rho } } { d t } = - i [ \hat { H } _ { \text{D} }, \hat { \rho } ] + \frac { \gamma } { 2 } \{ 2 \hat { z } \hat { \rho } \hat { z } ^ { \dagger } - \hat { z } ^ { \dagger } \hat { z } \hat { \rho } - \hat { \rho } \hat { z } ^ { \dagger } \hat { z } \}, \quad ( 6 ) \quad \Big |$$

where ˆ = ˆ (ˆ ) denote the z a a † incoherent loss (pump). Their steady-state solutions, i.e., dρ ˆ ss /dt = 0, are [29, 32]

$$\hat { \rho } _ { s s } = \hat { D } ( \bar { a } ) | 0 \rangle _ { a a } \langle 0 | \hat { D } ^ { - 1 } ( \bar { a } ), \ \bar { a } = \frac { - i \lambda } { \gamma / 2 + i \Delta }, \ \ ( 7 ) \quad \Big |$$

and

$$\hat { \rho } _ { s s } = \hat { D } ( \bar { a } ) | 0 \rangle _ { h h } \langle 0 | \hat { D } ^ { - 1 } ( \bar { a } ), \ \bar { a } = \frac { i \lambda } { \gamma / 2 - i \Delta }, \ \ ( 8 ) \quad \begin{matrix} \hat { H } _ { \text{DBS} } \\ \text{herein} \end{matrix}$$

for incoherent loss and pump, respectively. Thus, the displaced hole vacuum may be generated because of the competition between the coherent driving and incoherent pump processes. Further, the frequency of particles, ∆, is now relative to the frequency of the driving field, which can be analogized to the 'Fermi level.' When the driving field is sufficiently strong, the quadratic bosonic interaction can be regarded as a perturbation so that eigenmodes of effective QBHs, composed of particles and/or holes, can be regarded as excitation near the mean field.

Particle-hole duality and entanglement. -Introducing hole Fock states allows us to explore the ground-state properties and state evolution of QBSs in the PH space. In particular, the local PH transformation builds a relation between different QBHs. They can surprisingly possess the same BdG Hamiltonian in different PH representations with only a single PH inversion, i.e.,

$$\hat { \Omega } _ { i } ^ { - 1 } \hat { H } \hat { \Omega } _ { i } = \frac { 1 } { 2 } \hat { \Omega } _ { i } ^ { - 1 } \hat { \bar { \alpha } } _ { i } ^ { \dagger } \mathcal { H } \hat { \alpha } _ { i } \hat { \Omega } _ { i } = \frac { 1 } { 2 } \hat { \bar { \beta } } _ { i } ^ { \dagger } \mathcal { H } \hat { \beta } _ { i } = \hat { \bar { H } }, \quad ( 9 ) \quad \stackrel { \wedge } { \theta } = \\ \text{As ssl}$$

where ˆ β i = ˆ Ω -1 i ˆ ˆ α Ω i = (ˆ a 1 · · · -ia e ˆ † i iθ · · · ˆ a N , a ˆ † 1 · · · -iae ˆ -iθ · · · ˆ a † N ) T and ˆ ¯ β † i = ˆ Ω -1 i ˆ Ω ¯ α † ˆ i with [ ˆ β , β i ˆ ¯ ] † i = I 2 N . For any eigenmode ˆ ψ j in ˆ H with Fock states | n ⟩ ψ j , there exists an eigenmode ˆ ˜ ψ j in ˜ , denoted by Ω H ˆ -1 i ˆ ˆ ψ j Ω . This i eigenmode has the same frequency as ˆ ψ j and its Fock states can be given by | n ⟩ ˜ ψ j = ˆ Ω -1 i | n ⟩ ψ j . The PH transformation preserves the locality structure and spectral properties of the original Hamiltonian, and the transformed states and the original states are different representations for essentially the same physics, implying the PH duality [33, 34]. Further, the information about the evolution of dual states under the dual Hamiltonian is encoded in the initial system. For example, if | ψ (0) ⟩ is the initial state of the original system, we have

$$e ^ { - i \hat { \tilde { H } } t } \hat { \Omega } _ { i } ^ { - 1 } | \psi ( 0 ) \rangle = \hat { \Omega } _ { i } ^ { - 1 } e ^ { - i \hat { H } t } | \psi ( 0 ) \rangle = \hat { \Omega } _ { i } ^ { - 1 } | \psi ( t ) \rangle. \ \ ( 1 0 ) \quad t e r a c t$$

This PH duality is nontrivial because the solution can be well-known before the PH transformation but still hidden

FIG. 2. Basic forms of two-mode QBHs, (a) ˆ H BS , (b) ˆ H DP , (c) ˆ H DBS , and (d) ˆ H P , where the solid and dashed lines denote coherent and dissipative couplings, respectively. (a)[(c)] shows ˆ H BS ( ˆ H DBS ) with coherent (dissipative) particle-particle and hole-hole conversions, while (b)[(d)] shows ˆ H DP ( ˆ H P ) with coherent (dissipative) PH conversions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01059v1-figure-22fe4009-50d9-47b0-a2df-6060c07a9697.jpg)

after. In addition, such a non-unitary duality can relate the Hermitian and non-Hermitian QBHs and apply to any regime beyond the unitary duality [35-37].

To demonstrate the PH duality, it is instructive to consider the basic forms of two-mode-interacting QBHs and their non-Hermitian dual forms, i.e., the beamsplitter (pairing) Hamiltonian and the dissipative pairing [38] (dissipative beamsplitter [39-41]) Hamiltonian, given by

$$\underset { \text{reg} } { \text{acc.} } \quad \hat { H } _ { \text{P/DP} } = - \Delta _ { 1 } \hat { a } _ { 1 } ^ { \dagger } \hat { a } _ { 1 } + \Delta _ { 2 } \hat { a } _ { 2 } ^ { \dagger } \hat { a } _ { 2 } + g ( \hat { a } _ { 1 } ^ { \dagger } \hat { a } _ { 2 } ^ { \dagger } \pm \hat { a } _ { 1 } \hat { a } _ { 2 } ), \\ \text{igly} \quad \hat { H } _ { \text{DBS/BS} } = & \Delta _ { 1 } \hat { a } _ { 1 } ^ { \dagger } \hat { a } _ { 1 } + \Delta _ { 2 } \hat { a } _ { 2 } ^ { \dagger } \hat { a } _ { 2 } + i g ( \hat { a } _ { 1 } \hat { a } _ { 2 } ^ { \dagger } \pm \hat { a } _ { 1 } ^ { \dagger } \hat { a } _ { 2 } ),$$

where ˆ a 1 (ˆ a 2 ) is the annihilation operator and ∆ 1 (∆ ) 2 is the frequency relative to the driving field. Setting θ = 0 in Eq. (4), we have ˆ H P DP / = ˆ Ω -1 1 ˆ H DBS BS / ˆ Ω . 1 As shown in Fig. 2 (a) and (b), the duality between ˆ H BS and ˆ H DP shows that the non-Hermitian ˆ H DP can be regarded as the coherent PH conversion and therefore has stable PH eigenstates with a completely real spectrum. This can serve as the generator of the PH transformation in the single-mode case [29]. In addition, biorthogonal PH eigenstates also give the Hermitian ˆ H P a complex spectrum, in which ˆ H P should be regarded as the dissipative PH conversion, as shown in Fig. 2 (d). Conversely, in the real spectrum regime, the particle-particle (holehole) eigenstates of ˆ H P show that biorthogonal eigenstates of the non-Hermitian ˆ H DBS are PH eigenstates. In this case, ˆ H DBS should be inversely regarded as the PH pairing Hamiltonian and, therefore, can work as a resource for generating PH entanglement.

To see this more explicitly, we set -∆ 1 = ∆ 2 = ∆ &gt; g , where ˆ H P is the basic form of weakly interacting Bose gases [1] in the mean-field approximation, while ˆ H DBS becomes an anti-parity-time ( APT ) symmetric Hamiltonian [42, 43] in the real spectrum regime. In this case, the quasi-particles of ˆ H P are

Bogoliubov modes ˆ ψ i = ˆ S a a 1 2 ( r a S )ˆ i ˆ -1 a a 1 2 ( r ), while the quasi-holes, ˆ ϕ i = ˆ S h h 1 2 ( r h S ) ˆ i ˆ -1 h h 1 2 ( r ), can be regarded as removing quasi-particles from the mean-field, where i = 1, 2 and ˆ S i 1 2 i are the squeezing operators [29]. The ground states of ˆ H P in terms of quasi-particles and quasi-holes are the two-mode squeezed vacuum ˆ S a a 1 2 ( r ) 00 | ⟩ a a 1 2 and ˆ S h h 1 2 ( r ) 00 | ⟩ h h 1 2 , respectively, with entanglement E N = 2 | r | [44]. Due to the PH duality, the ground states of the APT symmetric Hamiltonian in terms of two sets of redundant eigenmodes are given by Ω ˆ 1 ˆ S a a 1 2 ( r ) 00 | ⟩ a a 1 2 = ˆ S ¯ h a 1 2 ( r ) 00 | ⟩ ¯ h a 1 2 and ˆ Ω ˆ 1 S h h 1 2 ( r ) 00 | ⟩ h h 1 2 = ˆ S a h 1 2 ( r ) 00 | ⟩ a h 1 2 , which are indeed the PH two-mode squeezed vacuum [29]. Therefore, the PH duality reveals the ground-state PH entanglement of the APT symmetric Hamiltonian, which might be detectable in the driven-dissipative system.

Another way to generate entanglement through ˆ H P is to let the particle vacuum or coherent states evolve under ˆ H P in the complex spectrum regime, especially at the resonant point ∆ 1 = ∆ 2 = ∆ [8, 45]. In the interaction picture, the evolution can be described as e -iH ˆ P t | 00 ⟩ a a 1 2 = ˆ S a a 1 2 ( -igt ) 00 | ⟩ a a 1 2 , with entanglement E N = 2 | r | = 2 gt . The PH duality of Eq. (10) shows that the PH vacuum evolving under ˆ H DBS will result in the infinite PH entanglement, i.e., ˆ Ω 1 e -iH ˆ P t | 00 ⟩ a a 1 2 = e -iH ˆ DBS t | 00 ⟩ ¯ h a 1 2 = ˆ S ¯ h a 1 2 ( -igt ) 00 | ⟩ ¯ h a 1 2 . These unbounded evolutions arise from the fact that the initial state and the system's eigenstates are not in the same space. In this case, the eigenstates of ˆ H DBS are particleparticle states or hole-hole states, which are the Fock states of the quasi-particles ˆ ψ ± = ( ± ˆ a 1 + ˆ ) a 2 / √ 2 or the quasi-holes ˆ ϕ ± = ( ± ˆ h 1 + ˆ h 2 ) / √ 2 with frequencies ϵ ± = ∆ ± ig and -ϵ ± , respectively. Conversely, when the initial state and the eigenstates are in the same space, the evolution will be confined to the excitation conservation subspace, leading to a steady state that is the eigenstate with the largest imaginary frequency. For example, when | 01 ⟩ a a 1 2 evolves under ˆ H DBS , the steady state is a Bell state ˆ ψ † + | 00 ⟩ a a 1 2 = ( 10 | ⟩ a a 1 2 + 01 | ⟩ a a 1 2 ) / √ 2 [46]. With the PH duality, the single-excited PH state Ω ˆ -1 1 | 01 ⟩ a a 1 2 will evolve to a steady PH Bell state under ˆ H P , i.e., ˆ Ω -1 1 ˆ ψ + | 00 ⟩ a a 1 2 = ( 10 | ⟩ h a 1 2 + 01 | ⟩ h a 1 2 ) / √ 2. Obviously, the space of initial states determines the representation of Hamiltonians and the following evolution. For example, pairing interactions describe the coherent creation or annihilation of excitations in pairs on particle or hole states while describing dissipative PH conversions on PH states. However, this is invisible from the Heisenberg equations of motion because particle and hole operators are redundant.

Chiral flows of particle-hole excitations. -The story becomes even more interesting when the hole degrees of freedom meet the network, where particle and hole nodes are connected by the four basic interactions in Fig. 2. Particularly, when interactions form a closed loop, along

FIG. 3. Schematics of chiral flows in the trimers with the Hamiltonians (a) ˆ H BST and (b) ˆ H SHT , when Φ = -π 2 . (c) and (d) plot the population evolutions of single-excited states in each loop when Φ = π/ 2 and Φ = -π/ 2, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01059v1-figure-73e19779-dc54-49ee-b678-6fb6d40a39b3.jpg)

which excitations experience a geometric phase, the AB effect will occur. Recently, the non-Hermitian AB effect has been observed in PH loops connected by pairing and beamsplitter interactions [14]. Here, we will show that hole states can successfully describe the excitation space and predict the Hermitian AB effect in PH loops.

̸

We consider the ring networks composed of three bosonic modes, as shown in Fig. 3 (a) and (b). In the particle-conserving case, it is named 'beamsplitter trimer,' [14] described by the Hamiltonian ˆ H BST = ∑ 3 i =1 ,j = i ge -iφ ij ˆ ˆ , a a † i j with φ ji = -φ ij . The dual version can be obtained by replacing ˆ a 1 (ˆ a † 1 ) with ˆ h 1 ( ¯ ˆ h † 1 ) through the PH transformation, called a 'single-hole trimer,' which composed of one beamsplitter and two dissipative pairing interactions, described by

̸

$$\hat { H } _ { \text{SHT} } = & \sum _ { i = 2, j \neq i } ^ { 3 } \Delta \hat { \bar { h } } _ { 1 } ^ { \dagger } \hat { h } _ { 1 } + \Delta \hat { a } _ { i } ^ { \dagger } \hat { a } _ { i } + g e ^ { i \varphi _ { j i } } \hat { a } _ { i } ^ { \dagger } \hat { a } _ { j } \\, & + g ( e ^ { i \varphi _ { 1 i } } \hat { a } _ { i } ^ { \dagger } \hat { h } _ { 1 } + e ^ { i \varphi _ { i 1 } } \hat { h } _ { 1 } ^ { \dagger } \hat { a } _ { i } ).$$

The restoration of excitation-conservation of ˆ H SHT reflects its invariance under the U(1) gauge transformation in the PH space, i.e., ˆ U PH = e iφ h h ( ˆ ¯ † 1 ˆ 1 + ∑ 3 i =2 ˆ a a † i ˆ ) i , which is the product of the local gauge transformation for each node. The phase sum around the particle and single-hole loops is Φ = φ 12 + φ 23 + φ 31 , which is invariant under the local gauge transformations on loop nodes, i.e., ˆ a † i → ˆ a e † i iφ and ¯ ˆ h † i → ˆ ¯ h e † i iφ for particle and hole nodes, respectively. This represents the synthetic flux threading loops and is a useful notation for determining whether time-reversal symmetry holds [47]. At Φ = zπ ( z ∈ Z ), ˆ H BST and ˆ H SHT hold time-reversal symmetry, while any

other flux breaks time-reversal symmetry, enabling chiral energy transport at each loop when Φ = π/ 2 + zπ .

At Φ = -π/ 2, the Heisenberg equation of nodes in the particle and single-hole loop exhibits a counterclockwise circulation [47, 48], ˆ ( ˆ a † 1 ¯ h † 1 ) → ˆ a † 2 → ˆ a † 3 → ˆ ( ˆ a † 1 ¯ h † 1 ), and then repeating the pattern. Note that, the ground state in each loop representation is the direct product of the vacuum states of its nodes. Applying the circulation of creation operators on the corresponding vacuum of each loop, we can obtain the chiral flow of the single-excited states, e.g., | 100 ⟩ a a a 1 2 3 → | 010 ⟩ a a a 1 2 3 → | 001 ⟩ a a a 1 2 3 → | 100 ⟩ a a a 1 2 3 in the particle loop and | 100 ⟩ h a a 1 2 3 → | 010 ⟩ h a a 1 2 3 → | 001 ⟩ h a a 1 2 3 → | 100 ⟩ h a a 1 2 3 in the single-hole loop. When Φ = π/ 2, the single-excited states in these loops will circulate in the opposite direction, as shown in Fig. 3 (c) and (d).

It is worth noting that the concept of 'holes' in hardcore bosonic systems differs from that in QBSs. In the former systems, ˆ H BST is accompanied by strong on-site interactions, and a 'hole' state is an unoccupied state for particles. In this case, the chiral flows exhibit the opposite chirality when exchanging single-excited particle and 'hole' states, e.g., | 100 ⟩ a a a 1 2 3 ↔ | 011 ⟩ a a a 1 2 3 . This can be interpreted as a 'hole' possessing the opposite 'charge' to the particle [48]. Conversely, without strong on-site interactions, a hole state is occupied by negative particle number. In the opposite excitation space, ˆ H BST and ˆ H SHT describe the hole and single-particle loop, respectively, with the synthetic flux Φ = ′ -Φ+ , as shown in Fig. 3 (a) and (b). Exchanging π particle and hole states in the single-excited case, e.g., | 100 ⟩ a a a /h a a 1 2 3 1 2 3 ↔ | 100 ⟩ h h h /a h h 1 2 3 1 2 3 in ˆ H BST /H ˆ SHT , does not change the chirality. This is because the excitation of holes is equivalent to removing particles from the mean field, and their energy chiral flow is symmetric about the 'Fermi level.'

Conclusions. -In summary, we present the concept of bosonic holes , which features the negative particle excitation near the mean particle number in QBSs. The key property, namely, the PH duality, not only reveals the duality between the Hermitian and non-Hermitian QBHs, it also enables various quantum applications with QBSs, such as PH entanglement and chiral flows. This paves a potential way to explore novel phenomena in particlenon-conserving and non-Hermitian systems, benefiting future quantum technologies.

Acknowledgments. -We thank Z. Yan for stimulating discussions. This work was supported by the National Natural Science Foundation of China (Grant No. 12375025 and No. 11874432), the National Key R&amp;D Program of China (Grant No. 2019YFA0308200).

∗ xiangzliang@mail.sysu.edu.cn

- [1] N. Bogoliubov, ON THE THEORY OF SUPERFLUIDITY, J. Phys. (Moscowd). 11 , 23 (1947).
- [2] J. Bardeen, L. N. Cooper, and, J. R. Schrierfer, Theory of Superconductivity, Phys. Rev. 108 , 1175 (1957).
- [3] W. H. Louisell, A. Yariv, and A. E. Siegman, Quantum fluctuations and noise in parametric processes. I., Phys. Rev. 124 , 1646 (1961).
- [4] J. P. Gordon, W. H. Louisell, and L. R. Walker, Quantum Fluctuations and Noise in Parametric Processes. II, Phys. Rev. 129 , 481 (1963).
- [5] B. R. Mollow and R. J. Glauber, Quantum Theory of Parametric Amplification. I., Phys. Rev. 160 , 1076 (1967).
- [6] C. M. Caves and B. L. Schumaker, New formalism for two-photon quantum optics. I. Quadrature phases and squeezed states, Phys. Rev. A. 31 , 3068 (1985).
- [7] B. Kraus and J. I. Cirac, Discrete Entanglement Distribution with Squeezed Light, Phy. Rev. Lett. 92 , 013602 (2004).
- [8] C. C. Gerry and P. L. Knight, Introductory Quantum Optics (Cambridge University Press, Cambridge, UK, 2005).
- [9] S. Sachdev, Quantum Phase of Matter (Cambridge University Press, Cambridge, UK, 2023).
- [10] S. Lieu, Topological symmetry classes for non-Hermitian models and connections to the bosonic Bogoliubov-de Gennes equation, Phys. Rev. B. 98 , 115135 (2018).
- [11] V. P. Flynn, E. Cobanera, and L. Viola, Deconstructing effective non-Hermitian dynamics in quadratic bosonic Hamiltonians, New J. Phys. 22 , 083004 (2020).
- [12] C. W. J. Beenakker, Search for Majorana Fermions in Superconductors, Annu. Rev. Condens. Matter Phys. 4 , 113 (2013).
- [13] K. Flensberg, F. v. Oppen, and A. Stern, Engineered platforms for topological superconductivity and Majorana zero modes. Nat Rev Mater 6, 944-958 (2021).
- [14] J. D. Pino, J. J. Slim, and E. Verhagen, NonHermitian chiral phononics through optomechanically induced squeezing, Nature (London). 606 , 82-87 (2022).
- [15] Y. X. Wang, and A. A. Clerk, Non-Hermitian dynamics without disspation in quantum systems, Phys. Rev. A. 99 , 063834 (2019).
- [16] N. Okuma, Boundary-dependent dynamical instability of bosonic Green's function: Dissipative Bogoliubov-de Gennes Hamiltonian and its application to nonHermitian skin effect, Phys. Rev. B. 105 , 224301 (2022).
- [17] L.-L Wan, and X.-Y. L¨, Quantum-Squeezing Induced u Point-Gap Topology and Skin Effect, Phys. Rev. Lett. 130 , 203605 (2023).
- [18] A. McDonald, T. P.-Barnea, and A. A. Clerk, Phase-Dependent Chiral Transport and Effective NonHermitian Dynamics in a Bosonic Kitaev-Majorana Chain, Phys. Rev. X. 8 , 041031 (2018).
- [19] J. J. Slim, C. C. Wanjura, M. Brunelli, J. d. Pino, A. Nunnenkamp, and E. Verhagen, Optomechanical realization of the bosonic Kitaev chain, Nature (London). 627 , 767-771 (2024).
- [20] C. C. Wanjura, J. J. Slim, J. d. Pino, M. Brunelli, E. Verhagen, and A. Nunnenkamp, Quadrature nonreciprocity in bosonic networks without breaking time-reversal symmetry, Nat. Phys. 19 , 1429-1436 (2023).
- [21] V. Peano, M. Houde, C. Brendel, F. Marquardt, and A. A. Clerk, Topological phase transitions and chiral inelastic transport induced by the squeezing of light, Nat. Com-

mun. 7 , 10779 (2016).

- [22] R. Barnett, Edge-state instabilities of bosons in a topological band, Phys. Rev. A 88 , 063631 (2013).
- [23] R. Shindou, R. Matsumoto, S. Murakami, and J.-i. Ohe, Topological chiral magnonic edge mode in a magnonic crystal, Phys. Rev. B 87 , 174427 (2013).
- [24] B. Galilo, D. K. K. Lee, and R. Barnett, Selective Population of Edge States in a 2D Topological Band System, Phys. Rev. Lett. 115 , 245302 (2015).
- [25] V. Peano, M. Houde, F. Marquardt, and A. A. Clerk, Topological Quantum Fluctuations and Traveling Wave Amplifiers, Phys. Rev. X. 6 , 041026 (2016).
- [26] D.A. Trifonov, Pseudo-Boson Coherent and Fock States, arXiv:0902.3744 (2009).
- [27] R. Rossignoli, and A. M. Kowalski, Complex modes in unstable quadratic bosonic forms, Phys. Rev. A. 72 , 032101 (2005).
- [28] J. Garcia, and R. Rossignoli, Spectrum and normal modes of non-Hermitian quadratic boson operators, Phys. Rev. A. 96 , 062130 (2017).
- [29] See Supplemental Material for more details.
- [30] N. R. Bernier, L. D. T´th, A. K. Feofanov, and T. J. Kipo penberg, Level attraction in a microwave optomechanical circuit, Phys. Rev. A. 98 , 023841 (2018).
- [31] X.-G. Wen, Symmetry-protected topological phases in noninteracting fermion systems, Phys. Rev. B 85 , 085103 (2012).
- [32] G. S. Agarwal, Nonclassical statistics of fields in pair coherent states, J. Opt. Soc. Am. B, 5 , 1940 (1988).
- [33] D. X. Nguyen, T. Can, and A. Gromov, Particle-Hole Duality in the Lowest Landau Level, Phys. Rev. Lett. 118 , 206602 (2017).
- [34] E. Cobanera, G. Ortiz, and Z. Nussinov, Unified Approach to Quantum and Classical Dualities, Phys. Rev. Lett. 104 , 020402 (2010).
- [35] F. J. Dyson, General Theory of Spin-Wave Interactions. Phys. Rev. 102 , 1217 (1956).
- [36] M. Znojil, Hybrid form of quantum theory with nonHermitian Hamiltonians, arXiv:2211.10633 (2022).
- [37] V. P. Flynn, E. Cobanera, and L. Viola, Restoring num-
- ber conservation in quadratic bosonic Hamiltonians with dualities, EPL. 131 , 40006 (2020).
- [38] A. Pocklington, Y. X. Wang, and A. A. Clerk, Dissipative Pairing Interactions: Quantum Instabilities, Topological Light, and Volume-Law Entanglement, Phys. Rev. Lett. 130 , 123602 (2023).
- [39] M. Harder, Y. Yang, B. M. Yao, C. H. Yu, J. W. Rao, Y. S. Gui, R. L. Stamps, and C.-M. Hu, Level Attraction Due to Dissipative Magnon-Photon Coupling, Phys. Rev. Lett. 121 , 137203 (2018).
- [40] W. C. Yu, J. J. Wang, H. Y. Yuan, and J. Xiao, prediction of attracitve level crossing via a dissipative mode, Phys. Rev. Lett. 123 , 227201 (2019).
- [41] Y.-P. Wang, and C.-M. Hu, Dissipative couplings in cavity magnonics, J. Appl. Phys. 127 , 130901 (2020).
- [42] P. Peng, W. Cao, C. Shen, W. Qu, J. Wen, L. Jiang, and Y. Xiao, Anti-parity-time symmetry with flying atoms, Nat. Phys. 12, 1139 (2016).
- [43] F. Yang, Y.-C. Liu, and L. You, Anti-PT symmetry in dissipatively coupled optical systems, Phys. Rev. A 96, 053845 (2017).
- [44] Y. D. Wang and A. A. Clerk, Reservoir-Engineered Entanglement in Optomechanical Systems, Phys. Rev. Lett. 110 , 253601 (2013).
- [45] X. W. Luo, C. W. Zhang, and S. W. Du, Quantum Squeezing and Sensing with Pseudo-Anti-Parity-Time Symmetry, Phys. Rev. Lett. 128 , 173602 (2022).
- [46] H. Y. Yuan, P. Yan, S. S. Zheng, Q. Y. He, K. Xia, and M. H. Yung, Steady Bell state Generation via MagnonPhoton Couping, Phys. Rev. Lett. 124 , 053602 (2020).
- [47] J. Koch, A. A. Houck, K. L. Hur, and S. M. Girvin, Time-reversal-symmetry breaking in circuit-QED-based photon lattices, Phys. Rev. A. 82 , 043811 (2010).
- [48] P. Roushan, C. Neill, A. Megrant, Y. Chen, R. Babbush, R. Barends, B. Campbell, Z. Chen, B. Chiaro, A. Dunsworth, A. Fowler, E. Jeffrey, J. Kelly, E. Lucero, J. Mutus, P. J. J. O'Malley, M. Nelley, C. Quintana, D. Sank, A. Vainsencher, J. Wenner, T. White, E. Kapit, H. Neven, and J. Martinis, Chiral ground-state currents of interacting photons in a synthetic magnetic field, Nat. Phys. 13 , 146-151 (2017).
Jia-Ming Hu Bo Wang Ze-Liang Xiang
Bosonic Holes in Quadratic Bosonic Systems
The concept of electron holes plays a significant role in condensed matter
physics. Here we develop the concept of bosonic holes, which exhibit negative
particle excitations, in quadratic bosonic systems. Unlike electron holes, the
Fock states of bosonic holes are biorthogonal, and their excitation can be
interpreted as removing particles from a mean-particle number with a mean field
background. Furthermore, we find that quadratic bosonic Hamiltonians related by
non unitary and local particle hole transformation possess the same locality
structure and spectral properties in different spaces, reflecting the PH
duality. Based on this, we study the generation of PH entanglement in two mode
cases and the PH Aharonov Bohm effect in the three mode case, which results in
a PH chiral flow with time reversal symmetry breaking. Our findings provide a
new way to understand and explore unusual physical phenomena in particle non
conserving and non Hermitian systems.