Open AccessArticle

Quantum Simulation of Pseudo-Hermitian-φ-Symmetric Two-Level Systems

Chao Zheng

Department of Physics, College of Science, North China University of Technology, Beijing 100144, China

Entropy 2022, 24(7), 867; https://doi.org/10.3390/e24070867

Submission received: 20 May 2022 / Revised: 17 June 2022 / Accepted: 22 June 2022 / Published: 24 June 2022

(This article belongs to the Special Issue Quantum Computing for Complex Dynamics)

Download

Browse Figures

Versions Notes

Abstract

Non-Hermitian (NH) quantum theory has been attracting increased research interest due to its featured properties, novel phenomena, and links to open and dissipative systems. Typical NH systems include PT-symmetric systems, pseudo-Hermitian systems, and their anti-symmetric counterparts. In this work, we generalize the pseudo-Hermitian systems to their complex counterparts, which we call pseudo-Hermitian-φ-symmetric systems. This complex extension adds an extra degree of freedom to the original symmetry. On the one hand, it enlarges the non-Hermitian class relevant to pseudo-Hermiticity. On the other hand, the conventional pseudo-Hermitian systems can be understood better as a subgroup of this wider class. The well-defined inner product and pseudo-inner product are still valid. Since quantum simulation provides a strong method to investigate NH systems, we mainly investigate how to simulate this novel system in a Hermitian system using the linear combination of unitaries in the scheme of duality quantum computing. We illustrate in detail how to simulate a general P-pseudo-Hermitian-

φ

-symmetric two-level system. Duality quantum algorithms have been recently successfully applied to similar types of simulations, so we look forward to the implementation of available quantum devices.

Keywords:

quantum simulation; linear combination of unitaries; non-Hermitian; pseudo-Hermitian

1. Introduction

Hermitian quantum systems are well known since closed quantum systems were focused on at the birth of conventional quantum mechanics. However, open and dissipative quantum systems are more common than closed systems in the real physical world, and cannot be described by Hermitian quantum theory. Therefore, non-Hermitian (NH) quantum theory [1,2] is attracting increased research interest. On the one hand, it extends conventional quantum theory and closely relates to open and dissipative systems [3,4,5,6,7,8,9,10,11]. On the other hand, NH systems have many novel properties and applications.

Typical NH systems include the parity-time-reversal (PT) symmetric systems, pseudo-Hermitian (PH) systems, and their anti-symmetric counterparts. PT-symmetric NH quantum systems have been focused on and investigated heavily since 1998 [12,13,14,15]. One important reason for this is that, besides Hermitian systems, PT-symmetric systems also keep the eigenvalues of H real in the exact PT phase. Due to its significance in both theory and potential applications, PT-symmetric quantum physics is developing rapidly and being investigated thoroughly from different aspects in a variety of systems [16,17,18,19,20]. In recent years, quantum simulations of PT-symmetric systems have been carried out, for example, fast and slow evolutions in the quantum brachistochrone problem [21,22,23], a generalized PT two-level system [24,25,26,27], and a PT-arbitrary-phase-symmetric system [28]. P-pseudo-Hermiticitian Hamiltonians were found to have real spectra in some conditions [29,30,31,32,33,34], while a necessary and sufficient condition is given for the reality of the spectrum of NH Hamiltonians admitting a complete set of biorthonormal eigenvectors in [31]. Properties of pseudo-Hermitian systems and their relationships with PT-symmetric systems are often investigated and discussed [35,36,37,38,39,40,41]. In this work, we extend pseudo-Hermitian Hamiltonians to the complex domain for the first time by introducing a phase factor

e^{i φ}

Inspired by Feynman, quantum simulation provides an efficient way to investigate nature [42]. It has become a strong tool to simulate novel quantum systems and discover featured properties. After a detailed quantum circuit is designed, an effective Hamiltonian is constructed, and the time-evolution of a quantum system can be simulated. Plenty of Hermitian systems and relevant phenomena have been investigated via quantum simulation methods [43,44,45,46,47,48,49,50,51,52,53]. In addition, they can also be applied to investigate NH systems in an effective way [23,24,25,26,27,28]. For example, quantum simulation of a P-pseudo Hermitian two-level system and its anti-symmetric counterpart has been proposed [54,55,56,57,58,59,60,61,62,63,64,65,66], making it possible to investigate these two NH systems in small quantum devices.

In this work, we investigate quantum simulation of the generalized pseudo-Hermitian-φ-symmetric (PH-φ) system, using the linear combination of unitaries (LCU) in the scheme of duality quantum computing [67] and the unitary-expansion (UE) techniques [8,9]. We optimize the quantum circuit and calculate the success probability. Furthermore, we discuss the implementations in NMR and quantum optics systems, expecting experimental realizations in the near future.

2. Complex Generalization of Pseudo-Hermitian Symmetry

Pseudo-Hermitian (PH) Hamiltonians

H_{P H}

satisfy

η^{- 1} {H_{P H}}^{†} η = H_{P H}

, where

η

is a linear Hermitian automorphism (invertible transformation) on the Hilbert space [33]. For example,

η

can be the parity (P) operator. Since

η

is not unique for a given

H_{P H}

[34], it is called an

η

-pseudo-Hermitian for a fixed

η

. Notice that different

η

values define different symmetries, though they are referred to as pseudo-Herimiticity in general. Similar to PT and anti-PT symmetry, an anti-symmetry of PH Hamiltonian has been introduced [35], and we call this a pseudo-Hermitian anti-symmetric (PHA) Hamiltonian

H_{P H A}

, if

η^{- 1} {H_{P H A}}^{†} η = - H_{P H A}

H_{P H}

and

H_{P H A}

can be seen as the real and imaginary counterparts of each other, since the latter can be obtained by the former Hamiltonian times i, and vice versa.

Consider a non-Hermitian Hamiltonian

H_{φ}

, which satisfies

η^{- 1} {H_{φ}}^{†} η = e^{i φ} H_{φ} .

(1)

H_{φ}

can be obtained by a phase factor

e^{- i \frac{φ}{2}}

times a relevant

H_{P H}

(i.e.,

H_{φ} = e^{- i \frac{φ}{2}} H_{P H}

). Therefore,

H_{φ}

can be seen as a complex generalization of a conventional pseudo-Hermitian symmetry

H_{φ} = e^{- i \frac{φ}{2}} H_{P H} = cos \frac{φ}{2} H_{P H} - sin \frac{φ}{2} (i H_{P H}) .

(2)

H_{φ}

can also be seen as a combination of an

H_{P H}

and a relevant

H_{P H A} = i H_{P H}

, which should have properties intermediate between them.

H_{φ - P H} = cos \frac{φ}{2} H_{P H} - sin \frac{φ}{2} H_{P H A} .

(3)

Thus, we introduce one extra degree to the conventional pseudo-Hermitian symmetry, called

H_{φ}

, which is of

η

-pseudo-Hermitian-arbitrary-phase or η-pseudo-Hermitian-φ symmetry. Both the inner product and pseudo-inner product [34] of

H_{P H}

are still well defined for

H_{φ}

. The relationships between PH, PH-anti, and PH-

φ

symmetry can be analogous to that of PT, anti-PT, and PT-arbitrary-phase symmetry, and can also be analogous to relationships between boson, fermion, and anyon. In fact, the relation in Equation (1) unifies the PH and PH-anti symmetries for

φ = 2 k π

and

φ = (2 k + 1) π

(k is integral), respectively.

3. Quantum Simulation Using LCU by Duality Quantum Computing

We now propose how to simulate the time-evolution of a pseudo-Hermitian-

φ

- symmetric (PH-

φ

) system in a conventional unitary quantum computer. The time-evolutionary operator should be an implicit function of the PH-

φ

Hamiltonian

H_{φ}

, say

E (t) = E (H_{φ} (t))

. Given that

H_{φ}

is not Hermitian,

E (t)

is not unitary and cannot be simulated directly by a quantum computer governed by conventional quantum mechanics, of which the time-evolution is unitary. However,

E (t)

can be extended to a summation of unitary operators. Therefore, we are able to simulate the non-unitary time-evolution of our PH-

φ

NH systems using LCU based on the duality quantum algorithm [67].

LCU and duality quantum computing were proposed in 2002 [67], and they developed rapidly [67,68,69,70,71,72,73], becoming some of the strongest tools for designing quantum algorithms [74]. Recently, we developed LCU to simulate NH systems [23,24,25,26,27,28,35,64,65], other novel systems, and time-dependent non-unitary operators [8].

Assume that the time-evolutionary operator can be extended to m terms as

E (t) = \sum_{k = 0}^{m - 1} c_{k} U_{k},

(4)

where each

U_{k}

is a unitary operator, and

c_{k}

are complex UE parameters (

k = 0, \dots, m - 1

). The unitary expansion of

E (t)

is not unique, and we only show the schematic strategy to simulate the time-evolution of a general PH-

φ

non-Hermitian system in this section. We do not discuss details of how to extend the general non-unitary operator by our UE techniques, because this is a significant question and deserves to be investigated alone. However, it can save qubits and reduce the complexity of quantum simulation if the non-unitary operators can be expressed by fewer UE terms. We will illustrate in detail the quantum simulation of a P-pseudo-Hermitian-

φ

-symmetric two-level system in the next section.

Quantum simulation of the time-evolution of a PH-

φ

system in Equation (4) can be achieved using either qudits or qubits as the quantum circuit, as shown in Figure 1. The whole system is composed of an ancillary subsystem a and an evolutionary subsystem e. The simulation can be achieved in either a qudit system or a qubit system. The dimensions of qudits or the total number of the qubits are decided by the dimensions d of the PH-

φ

non-Hermitian system and the number m of the UE terms in Equation (4). In detail, an m-dimensional qudit and a d-dimensional qudit are able to be the ancillary and evolutionary subsystems, respectively. If we simulate using qubits, the qubit numbers of the ancillary and evolutionary subsystems should not be less than

n_{1} = l o g_{2} m

and

n_{2} = l o g_{2} d

, respectively. In general cases, for different dimensional NH systems, m has a maximum value of four [8,9]. In detail, the maximum of m is three in general and two in special cases for

d = 2

[8], while the maximum of m is four for higher dimensions [9]. At the beginning, the ancillary subsystem is initialized to a logic state

{| 0 〉}_{a}

, and the PH-

φ

non-Hermitian subsystem is initialized to an arbitrary state

{| ψ 〉}_{e}

. In the middle part, the operator

U_{E 1}

is applied on the ancillary subsystem to assign the UE parameters, and the k-controlled gates (

k = 0, 1, \dots, m - 1

) together with the operator

U_{E 2}

achieve the UE terms’ generation [8], constructing the PH-

φ

system. At the end of the circuit, quantum measurements will be performed on the ancillary system to complete the time-evolution governed by the PH-

φ

Hamiltonian in an indeterministic way when the ancillary qubits collapse into the logic state

{| 0 〉}_{a}

. In this case, the evolutionary subsystem will evolve as the PH-

φ

system requires. Otherwise, this simulation will be terminated and started over until success.

4. Quantum Simulation of $P$ -Pseudo-Hermitian- $φ$ -Symmetric Two-Level Systems

We take a P-pseudo-Hermitian-

φ

-symmetric system as an example to illustrate how to simulate it using LCU and duality quantum algorithms. P-pseudo-Hermitian systems are one class of typical NH systems, and a series of theoretical investigations have been conducted on them. Recently, pseudo-Hermiticity was found to be able to protect unitary scattering [41]. P-PH systems have close relationships with PT-symmetric systems, though they are different. We have generalized the PT symmetry to PT-arbitrary-phase symmetry [28], which also has close relationships with P-PH-

φ

symmetry. We show the sets of the P-PH-

φ

, PT-

φ

-symmetric, P-PH, PT-symmetric, and Hermitian systems in Figure 2. They have intersections in the P-PH-related sets, PT-symmetry-related sets, and Hermitian set.

4.1. P-PH- $φ$ Two-Level Systems

In the two-dimensional case, the most general form of a P-PH-

φ

Hamiltonian can be obtained by the relevant P-pseudo-Hermitian Hamiltonian [14,24] times

e^{- i \frac{φ}{2}}

, as follows:

H_{φ} = e^{- i \frac{φ}{2}} (\begin{matrix} r e^{i θ} & u \\ v & r e^{- i θ} \end{matrix}),

(5)

where

φ

, r, u, v, and

θ

are real parameters, the parity operator

P = (\begin{matrix} 0 & 1 \\ 1 & 0 \end{matrix})

, and

φ

is a fixed parameter that links to the symmetry.

The eigenvalues of

H_{φ}

are

ε_{\pm} = e^{- i \frac{φ}{2}} (r cos θ \pm \frac{ω_{0}}{2})

with respect to the two eigenvectors

|ε_{\pm}〉

, respectively, where

ω_{0} = e^{i \frac{φ}{2}} (ε_{+} - ε_{-}) = 2 \sqrt{u v - r^{2} {sin}^{2} θ}

is the energy difference of

H_{P H} : = e^{i \frac{φ}{2}} H_{φ}

. The exceptional points (EPs) of

H_{φ}

in the parametric space [75] are composed of the points leading

ω_{0}

to be zero. The EPs of

H_{φ}

, forming the boundary of real (

ω_{0}^{2} > 0

) and imaginary (

ω_{0}^{2} < 0

) phases of

H_{P H}

, are also EPs of

H_{φ}

Investigations into novel phenomena using a controllable, currently available quantum device is one of the main tasks of quantum simulation. Because of the Hermiticity, the time-evolution of a Hermitian system is unitary and can be simulated by a conventional quantum system directly in the Hilbert space of the same dimensions.

We will simulate the two-dimensional time-evolutionary operator

e^{- i \frac{t}{ℏ} H_{φ}},

(6)

which is not unitary. Therefore, instead of the Hermitian case using one qubit, we will construct a general P-PH-

φ

subsystem in a larger Hilbert space and simulate the time-evolution in the scheme of duality quantum computing using the LCU method.

Our quantum simulation method is applicable to the whole parametric space (except the EPs), including the neighborhoods of EPs. We still use the Hilbert–Schmidt inner product of the conventional quantum mechanics because this NH system will be simulated in a Hermitian system, while the pseudo-inner product introduced by A. Mostafazadeh [34] is well defined in this system.

4.2. UE of the Time-Evolutionary Operator

First of all, the UE techniques [8] will be applied to the non-unitary time-evolutionary operator in Equation (4). We calculate the UE terms of

e^{- i \frac{t}{ℏ} H_{φ}}

in detail, which can be expanded by four or three UE terms in general, as follows:

e^{- i \frac{t}{ℏ} H_{φ}} = f_{0} σ_{0} + f_{1} σ_{1} + f_{2} (i σ_{2}) + f_{3} σ_{3},

(7)

where the Pauli matrices

σ_{0} = [\begin{matrix} 1 & 0 \\ 0 & 1 \end{matrix}]

σ_{1} = [\begin{matrix} 0 & 1 \\ 1 & 0 \end{matrix}]

σ_{2} = [\begin{matrix} 0 & - i \\ i & 0 \end{matrix}]

, and

σ_{3} = [\begin{matrix} 1 & 0 \\ 0 & - 1 \end{matrix}]

, and

f_{k} = |f_{k}| e^{i θ_{k}}

(

k = 0, 1, 2, 3

) are the UE parameters, being time-dependent complex-functions of

H_{φ}

in Equation (5), with no limitations on their norms. The explicit forms of

f_{k}

’s (

k = 0, 1, 2, 3

) are shown in Appendix A (i). To simplify the simulation, the number of UE terms should be further reduced. Thus, we merge the four UE terms in Equation (7) to three, as

e^{- i \frac{t}{ℏ} H_{φ}} = e^{i θ_{0}} |f_{0}| σ_{0} + e^{i θ_{3}} |g_{1}| V_{1} + e^{i θ_{3}} |g_{2}| V_{2} .

(8)

where the UE parameters

g_{1}

and

g_{2}

are complex functions of

f_{k}

’s (

k = 0, 1, 2, 3

) in Equation (7). The explicit forms of

g_{k}

and

V_{k}

(

k = 1, 2

) are presented in Appendix A (ii) and (iii), respectively.

The above UEs are valid for a general P-PH-

φ

Hamiltonian. If

f_{k}

’s (

k = 0, 1, 2, 3

) in Equation (7) satisfy one of the phase-matching conditions, which have been introduced and investigated in detail in [8] to judge the minimum numbers of terms in LCU for an arbitrary two-dimensional nonunitary operator, the UE terms in Equation (8) can be merged to fewer terms, as follows:

e^{- i \frac{t}{ℏ} H} = h_{0} W_{0} + h_{1} W_{1} .

(9)

The explicit forms of

h_{0}

h_{1}

W_{0}

, and

W_{1}

vary in different phase-matching conditions [8], and they are complex functions with respect to time t.

For convenience, we set a normalizing factor f as

f = |f| = \sum_{k = 0}^{3} \sqrt{{|f_{k}|}^{2}} = \sqrt{{|f_{0}|}^{2} + {|g_{1}|}^{2} + {|g_{2}|}^{2}} = \sqrt{{|h_{0}|}^{2} + {|h_{1}|}^{2}} .

(10)

Now, we are able to simulate the non-unitary evolution in Equation (6) using LCU in the scheme of duality quantum computing. A qudit- or a qudit–qubit-hybrid device is able to achieve the simulation, and qudits take advantage over qubits in some quantum algorithms [76] (e.g., they reach a higher accuracy when solving the eigenvalue problem using quantum phase estimation algorithms with qudits than with qubits [77]). However, we focus on quantum simulations using qubits here, since qubit-quantum computers are currently available technologies.

4.3. Qubit Simulation

Three qubits or fewer are able to simulate the time-evolution of a P-PH-

φ

system by our theory. In a general case, the three qubits are divided into an evolutionary qubit e and an ancillary subsystem a of the remaining two qubits. The evolutionary qubit will evolve as per Equation (6) with the assistance of the ancillary subsystem in a probabilistic way. Only a six-dimensional Hilbert subspace, extended by

{| 00 〉}_{a} {| k 〉}_{e}

{| 01 〉}_{a} {k 〉}_{e}

, and

{| 10 〉}_{a} {k 〉}_{e}

(

k = 0, 1

), is needed, while the remaining two dimensions are spared. The success probability using six dimensions is larger than that using the full eight dimensions [8].

The quantum circuit used to achieve the quantum simulation is shown in Figure 3. At the beginning, the whole system is initialized to a pure state

{| 00 〉}_{a} {| 0 〉}_{e}

, and the evolutionary qubit e will be prepared in an arbitrary state

{| ψ 〉}_{e}

, as needed by a single-qubit rotation

R_{ψ}

. The two ancillary qubits will assist the evolutionary qubit to evolve, governed by the P-PH-

φ

Hamiltonian. The first block aims at, on the one hand, deleting the basis

{| 11 〉}_{a}

of the ancillary subsystem, so that the rest of the bases

{| 00 〉}_{a}

{| 01 〉}_{a}

and

{| 10 〉}_{a}

together with the two bases

{| 0 〉}_{e}

and

{| 1 〉}_{e}

of the evolutionary qubit are used to construct a six-dimensional subspace; on the other hand, it aims at assigning the three UE parameters in Equation (8) to

{| 00 〉}_{a} {| ψ 〉}_{e}

{| 01 〉}_{a} {| ψ 〉}_{e}

, and

{| 10 〉}_{a} {| ψ 〉}_{e}

. This is the first key step to simulating the time-evolution. In detail, the first and second ancillary qubits are swapped, and two single-qubit rotations

S_{1}

and

- σ_{3}

are applied to them. Then, a controlled-NOT gate is applied, in which the first and the second qubits take roles as the target and control qubit, respectively. After the first qubit is rotated by

S_{2}

, the two ancillary qubits are swapped again. The explicit forms of two single-qubit operators are

S_{1} = \frac{1}{f} [\begin{matrix} g_{1} & - \sqrt{{|f|}^{2} - {|g_{1}|}^{2}} \\ \sqrt{{|f|}^{2} - {|g_{1}|}^{2}} & g_{1}^{*} \end{matrix}]

(11)

and

S_{2} = \frac{1}{\sqrt{{|f|}^{2} - {|g_{1}|}^{2}}} [\begin{matrix} g_{2}^{*} & f_{0} \\ - f_{0}^{*} & g_{2} \end{matrix}],

(12)

where

f_{0}

g_{1}

g_{2}

, and f are as in Equations (8) and (10).

In the second block of Figure 3, the three UE terms in Equation (8) will be generated. Notice that the unit matrix

C_{00 - σ_{0}}

is a trivial operation that can be removed in practice. The other two jointly controlled gates are necessary, and their matrix forms are

C_{01 - V_{1}} = [\begin{matrix} σ_{0} & 0 & 0 & 0 \\ 0 & V_{1} & 0 & 0 \\ 0 & 0 & σ_{0} & 0 \\ 0 & 0 & 0 & σ_{0} \end{matrix}]

(13)

and

C_{10 - V_{2}} = [\begin{matrix} σ_{0} & 0 & 0 & 0 \\ 0 & σ_{0} & 0 & 0 \\ 0 & 0 & V_{2} & 0 \\ 0 & 0 & 0 & σ_{0} \end{matrix}],

(14)

respectively. For the explicit expressions of

V_{1}

and

V_{2}

, refer to Equation (8) and Appendix A (iii). Now, the three UE terms are generated and entangled with the three bases of the ancillary subspace.

The third block aims at superposing the three UE terms by swapping the two ancillary qubits three times and applying

H_{2}

and

S_{3}

in between them, as shown in Figure 3, where

S_{3} = \frac{1}{\sqrt{3}} [\begin{matrix} \sqrt{2} & 1 \\ 1 & - \sqrt{2} \end{matrix}] .

(15)

Now, the whole system evolves to a superposition state

\frac{1}{\sqrt{3} f} [{| 00 〉}_{a} e^{- i \frac{t}{ℏ} H_{φ}} {| ψ 〉}_{e} + f \sum_{k = 01, 10} {| k 〉}_{a} {| s_{k} 〉}_{e}],

(16)

where values for

| s_{k} 〉_{e}

are not given explicitly because the relevant terms will be discarded after quantum measurements. Notice that the three UE terms are superposed as the time-evolution in the first term relevant to

{| 00 〉}_{a}

Finally, quantum measurements are performed on the ancillary subsystem. If it outputs the state

{| 00 〉}_{a}

, the evolutionary qubit will evolve to

e^{- i \frac{t}{ℏ} H_{φ}} {| ψ 〉}_{e}

, which is governed by the NH Hamiltonian in Equation (5), with a success probability of

\frac{1}{3 f^{2}}_{e} 〈 ψ | e^{i \frac{t}{ℏ} H_{φ}^{†}} e^{- i \frac{t}{ℏ} H_{φ}} {| ψ 〉}_{e} .

(17)

If one of the remaining two results of

{| 01 〉}_{a}

{| 01 〉}_{a}

is measured, the simulation will be terminated. The whole process will then be started over until

{| 00 〉}_{a}

is obtained. Therefore, it is an indeterministic protocol to simulate the time-evolution of P-PH-

φ

two-level systems.

The number of qubits can be reduced to two if the time-evolutionary operator in Equation (8) can be united into two UE terms, as in Equation (9), and the quantum circuit is shown in Figure 4.

The two-qubit system is initialized to

{| 0 〉}_{a} {| 0 〉}_{e}

, and then qubit e is rotated to

{| ψ 〉}_{e}

as needed. A single-qubit unitary

W = \frac{1}{f} [\begin{matrix} h_{0} & - h_{1}^{*} \\ h_{1} & h_{0}^{*} \end{matrix}],

(18)

is applied on the ancillary qubit to assign the UE parameters. Notice that

h_{0}

and

h_{1}

always satisfy Equation (10), while their explicit forms change with different phase-matching conditions, as seen in [8].

Then, two controlled gates follow, which are

C_{0 - W_{0}} = [\begin{matrix} W_{0} & 0 \\ 0 & σ_{0} \end{matrix}]

(19)

and

C_{1 - W_{1}} = [\begin{matrix} σ_{0} & 0 \\ 0 & W_{1} \end{matrix}],

(20)

where the explicit forms of

W_{0}

and

W_{1}

are decided by the specific phase-matching conditions [8]. After a Hadamard

H_{2}

is performed on the ancillary qubit, the two-qubit system evolves to

\frac{1}{\sqrt{2} f} [{| 0 〉}_{a} e^{- i \frac{t}{ℏ} H_{φ}} {| ψ 〉}_{e} + {| 1 〉}_{a} (h_{0} W_{0} - h_{1} W_{1}) {| ψ 〉}_{e}] .

(21)

Similarly, the first term links to the time-evolution governed by the P-PH-

φ

Hamiltonian, while the second term will be discarded after quantum measurements.

Finally, a quantum measurement is performed on the ancillary qubit. If

{| 0 〉}_{a}

is output, the evolutionary qubit e will evolve as

e^{- i \frac{t}{ℏ} H_{φ}} {| ψ 〉}_{e}

with a success probability of

\frac{1}{2 f^{2}}_{e} 〈 ψ | e^{i \frac{t}{ℏ} H_{φ}^{†}} e^{- i \frac{t}{ℏ} H_{φ}} {| ψ 〉}_{e} .

(22)

This is decided by not only by the initial state but also

H_{φ}

. If the ancillary qubit is observed in state

{| 1 〉}_{a}

, the process will be terminated and the result will be discarded. We start the quantum simulation again, and it continues until

{| 0 〉}_{a}

is measured. From Equations (17) and (22), it can be seen that the success probability using two qubits is increased further than that using three qubits. Therefore, it is valuable to merge the UE terms before quantum simulation to increase efficiency. On one hand, this will save qubits, decreasing the complexities of the quantum circuit. On the other hand, it will enlarge the probability of simulating the P-PH-

φ

NH system successfully.

5. Experimental Proposals

Given that LCU and duality quantum algorithms have recently been successfully applied to similar types of quantum simulations in different experimental systems [23,25,26,27], we look forward to the implementation of available quantum devices. Since qubit-quantum devices have become available technologies, it is appropriate to simulate the pseudo-Hermitian-

φ

-symmetric systems using qubits for experimental implementations. Candidate qubit-systems include nuclear-magnetic-resonance (NMR) quantum simulators, quantum optics systems, superconductor quantum systems, two energy levels of ultracold atoms, and ion-trap systems. The operations in the quantum circuit can be realized by the related controlling methods and techniques.

Take an NMR quantum simulator as an example. The evolutionary and ancillary qubits are realized by the nuclei of spin-

1 / 2

. To initialize the pseudo-pure state, the spatial-averaging method [78] can be applied at the beginning of the experiment, and then a series of magnetic pulse sequences will be applied to realize the quantum gates in the quantum circuit. In detail, hard pulses are used to realize single-qubit rotations directly, while free evolutions of the two nuclei of spin-

\frac{1}{2}

in a period [23] are necessary to realize controlled two-qubit gates.

By quantum optics, two orthogonal polarized directions of a photon take the role of a qubit. Single-qubit gates can be realized by a series of quarter-wave and half-wave plates [79]. The efficiency is too low in practice to realize a two-polarization-qubit gate, though the task can be achieved using measurement-induced nonlinearity [80]. A more practical method is to improve efficiency with the assistance of the degrees of freedom of photon locations utilized as the qubit basis (i.e., the location qubit), which can be prepared and operated by beam-splitters and Mach–Zehnder interferometers [81].

In addition, quantum processors (such as the IBM QE 5-qubit [82]) can be used to realize the quantum simulation in this work.

6. Conclusions

We complexly generalize the conventional pseudo-Hermitian system to a pseudo-Hermitian-

φ

-symmetric system by adding an extra freedom of symmetry to the original. Therefore, the pseudo-Hermitian system and its anti-symmetric counterpart can be seen as two special cases when the phase-angles

φ

are set to be

2 k π

and

(2 k + 1) π

(k is integral), respectively. This can be analogous to the case of the PT, anti-PT, and PT-arbitrary-phase symmetries, or that of bosons, fermions, and anyons. We believe that more novel properties can be found in the generalized systems, and the conventional PH systems can be better investigated through extended freedom of symmetry (as we can understand a two-dimensional plane better when we are in a three-dimensional space). The well-defined inner product and pseudo-inner product are still valid for the PH-

φ

systems.

We mainly investigate quantum simulation of the PH-

φ

system. We show a schematic proposal for a general system using LCU in the scheme of duality quantum computing, and we propose in detail a general P-PH-

φ

two-level system. A minimum six-dimensional Hilbert space is necessary to simulate the time-evolution of an arbitrary P-PH-

φ

two-level system by our unitary expansion techniques, while four-dimensional Hilbert space is enough in some special cases. The simulation is achieved in an indeterministic way, and the success probability is decided by the initial state, the Hamiltonian, and the dimensions of the used Hilbert space. The fewer dimensions there are, the higher the success probability is. Therefore, it is meaningful to merge the UE terms based on our UE techniques and phase-matching conditions to reduce the dimensions required before quantum simulation. With this, the qubit source will be saved, and the success probability will be increased.

Finally, we discuss experimental implementations of available quantum devices, such as NMR, quantum optics systems, and IBM QE. Given that LCU and duality quantum algorithms have been recently successfully applied to similar types of quantum simulations experimentally, quantum simulation of a general P-PH-

φ

two-level system can soon be implemented on small quantum devices.

Funding

This work was funded by the National Natural Science Foundation of China Grant No. 12175002, Beijing Natural Science Foundation Grant No. 1222020, and NCUT Talents Project and Special Fund for C.Z.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:

NH	non-Hermitian
PT	parity-time-reversal
PH	pseudo-Hermitian
PHA	pseudo-Hermitian anti-symmetric
PH- $φ$	pseudo-Hermitian- $φ$ -symmetric
P	parity
EP	exceptional point
LCU	linear combination of unitaries
NMR	nuclear magnetic resonance

Appendix A

(i) Explicit forms of

f_{k}

’s (

k = 0, 1, 2, 3

f_{0} = m [a_{0} b_{0} - a_{2} b_{2} + i (- a_{1} b_{1} + a_{3} b_{3})],

(A1)

f_{1} = m [(a_{0} b_{1} - a_{2} b_{3} + a_{3} b_{2}) - i a_{1} b_{0}],

(A2)

f_{2} = m [(a_{0} b_{2} + a_{3} b_{1} - a_{1} b_{3}) - i a_{2} b_{0}]

(A3)

and

f_{3} = m [a_{3} b_{0} + i (a_{1} b_{2} - a_{2} b_{1} + a_{0} b_{3})], where

(A4)

m = e^{- i \frac{t}{ℏ} r e^{(- φ / 2)} cos θ},

(A5)

a_{0} = cos α = cos (ω_{a} t / 2 ℏ),

(A6)

a_{1} = (u + v) cos (φ / 2) sin α / ω_{a},

(A7)

a_{2} = (u - v) cos (φ / 2) sin α / ω_{a},

(A8)

a_{3} = 2 r cos (φ / 2) sin θ sin α / ω_{a},

(A9)

ω_{a} = 2 |cos (φ / 2)| \sqrt{u v - r^{2} {sin}^{2} θ};

(A10)

and

b_{0} = cos β = cos (ω_{b} t / 2 ℏ),

(A11)

b_{1} = - (u + v) sin (φ / 2) sin β / ω_{b},

(A12)

b_{2} = - (u - v) sin (φ / 2) sin β / ω_{b},

(A13)

b_{3} = - 2 r sin (φ / 2) sin θ sin β / ω_{b},

(A14)

ω_{b} = 2 |sin (φ / 2)| \sqrt{r^{2} {sin}^{2} θ - u v} .

(A15)

(ii) Explicit forms of

g_{1}

and

g_{2}

g_{1}

and

g_{2}

are obtained from

f_{k} = | f_{k} | e^{i θ_{k}}

(

k = 0, 1, 2, 3

), that is,

g_{1} = e^{i θ_{3}} \sqrt{{|f_{3}|}^{2} + {|f_{1}|}^{2} {cos}^{2} φ_{1} + {|f_{2}|}^{2} {sin}^{2} φ_{2}}

(A16)

and

g_{2} = e^{i θ_{3}} \sqrt{{|f_{1}|}^{2} {sin}^{2} φ_{1} + {|f_{2}|}^{2} {cos}^{2} φ_{2}},

(A17)

where

φ_{1} = θ_{1} - θ_{3} and φ_{2} = θ_{2} - θ_{3} .

(A18)

(iii) Explicit forms of

V_{1}

and

V_{2}

[\begin{matrix} cos ζ_{1} & e^{i ϕ_{1}} sin ζ_{1} \\ e^{- i ϕ_{1}} sin ζ_{1} & - cos ζ_{1} \end{matrix}] and [\begin{matrix} 0 & e^{i ϕ_{2}} \\ - e^{- i ϕ_{2}} & 0 \end{matrix}],

(A19)

where

ϕ_{1}

is decided by

cos ϕ_{1} = \frac{|f_{1}| cos φ_{1}}{\sqrt{{|f_{1}|}^{2} {cos}^{2} φ_{1} + {|f_{2}|}^{2} {sin}^{2} φ_{2}}}

(A20)

and

sin ϕ_{1} = \frac{|f_{2}| sin φ_{2}}{\sqrt{{|f_{1}|}^{2} {cos}^{2} φ_{1} + {|f_{2}|}^{2} {sin}^{2} φ_{2}}};

(A21)

and

ϕ_{2}

is decided by

cos ϕ_{2} = \frac{|f_{2}| cos φ_{2}}{|g_{2}|} and sin ϕ_{2} = \frac{|f_{1}| sin φ_{1}}{|g_{2}|};

(A22)

ζ_{1}

is decided by

cos ζ_{1} = \frac{|f_{3}|}{|g_{1}|} and sin ζ_{1} = \frac{\sqrt{{|f_{1}|}^{2} {cos}^{2} φ_{1} + {|f_{2}|}^{2} {sin}^{2} φ_{2}}}{|g_{1}|} .

(A23)

References

Gamow, G. Quantum Theory at Nucleus. Z. Phys. 1928, 51, 204. [Google Scholar] [CrossRef]
Moiseyev, N. Non-Hermitian Quantum Mechanics; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar]
Breuer, H.-P.; Petruccione, F. The Theory of Open Quantum Systems, 10th Anniversary ed.; Oxford University Press: Oxford, UK, 2002. [Google Scholar]
Barreiro, J.T.; Müller, M.; Schindler, P.; Nigg, D.; Monz, T.; Chwalla, M.; Hennrich, M.; Roos, C.F.; Zoller, P.; Blatt, R. An Open-system Quantum Simulator with Trapped Ions. Nature 2011, 470, 486–491. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Hu, Z.; Xia, R.; Kais, S. A Quantum Algorithm for Evolving Open Quantum Dynamics on Quantum Computing Devices. Sci. Rep. 2020, 10, 3301. [Google Scholar] [CrossRef]
Del Re, L.; Rost, B.; Kemper, A.F.; Freericks, J.K. Driven-Dissipative Quantum Mechanics on a Lattice: Simulating a Fermionic Reservoir on a Quantum Computer. Phys. Rev. B 2020, 102, 125112. [Google Scholar] [CrossRef]
Viyuela, O.; Rivas, A.; Gasparinetti, S.; Wallraff, A.; Filipp, S.; Martin-Delgado, M.A. Observation of Topological Uhlmann Phases with Superconducting Qubits. Njp Quantum Inf. 2018, 4, 10. [Google Scholar] [CrossRef] [Green Version]
Zheng, C. Universal Quantum Simulation of Single-Qubit Nonunitary Operators using Duality Quantum Algorithm. Sci. Rep. 2021, 11, 3960. [Google Scholar] [CrossRef]
Schlimgen, A.W.; Head-Marsden, K.; Sager, L.M.; Narang, P.; Mazziotti, D.A. Quantum Simulation of Open Quantum Systems Using a Unitary Decomposition of Operators. Phys. Rev. Lett. 2021, 127, 270503. [Google Scholar] [CrossRef]
Del Re, L.; Rost, B.; Foss-Feig, M.; Kemper, A.F.; Freericks, J.K. Robust Measurements of N-Point Correlation Functions of Driven-Dissipative Quantum Systems on a Digital Quantum Computer. arXiv 2022, arXiv:2204.12400. [Google Scholar]
Ding, P.Z.; Yi, W. Two-body exceptional points in open dissipative systems. Chin. Phys. B 2022, 31, 010309. [Google Scholar] [CrossRef]
Bender, C.M.; Boettcher, S. Real Spectra in Non-Hermitian Hamiltonians having PT Symmetry. Phys. Rev. Lett. 1998, 80, 5243–5246. [Google Scholar] [CrossRef] [Green Version]
Bender, C.M.; Boettcher, S.; Meisinger, P.N. PT-Symmetric Quantum Mechanics. J. Math. Phys. 1999, 40, 2201–2229. [Google Scholar] [CrossRef] [Green Version]
Bender, C.M.; Brody, D.C.; Jones, H.F. Complex Extension of Quantum Mechanics. Phys. Rev. Lett. 2002, 89, 270401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Bender, C.M.; Brody, D.C.; Jones, H.F. Must a Hamiltonian be Hermitian? Am. J. Phys. 2003, 71, 1095–1102. [Google Scholar] [CrossRef] [Green Version]
Bender, C.M. PT-symmetric Quantum Theory. J. Phys. Conf. Ser. 2015, 631, 012002. [Google Scholar] [CrossRef]
Zhang, S.; Jin, L.; Song, Z. Topology of a parity-time symmetric non-Hermitian rhombic lattice. Chin. Phys. B 2022, 31, 010312. [Google Scholar] [CrossRef]
Hu, Z.; Jin, L.; Zeng, Z.-Y.; Tang, J.; Luo, X.-B. Quasi-parity-time symmetric dynamics in periodically driven two-level non-Hermitian system. Acta Phys. Sin. 2022, 71, 074207. [Google Scholar] [CrossRef]
Wang, K.; Gao, Y.-P.; Jiao, R.; Wang, C. Recent progress on optomagnetic coupling and optical manipulation based on cavity-optomagnonics. Front. Phys. 2021, 14, 42201. [Google Scholar] [CrossRef]
Xu, W.-L.; Liu, X.-F.; Sun, Y.; Gao, Y.-P.; Wang, T.-J.; Wang, C. Magnon-induced chaos in an optical PT-symmetric resonator. Phys. Rev. E 2020, 101, 012205. [Google Scholar] [CrossRef]
Bender, C.M.; Brody, D.C.; Jones, H.F.; Meister, B.K. Faster than Hermitian Quantum Mechanics. Phys. Rev. Lett. 2007, 98, 040403. [Google Scholar] [CrossRef] [Green Version]
Günther, U.; Samsonov, B.F. Naimark-dilated PT-Symmetric Brachistochrone. Phys. Rev. Lett. 2008, 101, 230404. [Google Scholar] [CrossRef] [Green Version]
Zheng, C.; Hao, L.; Long, G.L. Observation of a Fast Evolution in a Parity-Time-Symmetric System. Philos. Trans. R. Soc. A 2013, 371, 20120053. [Google Scholar] [CrossRef] [PubMed]
Zheng, C. Duality Quantum Simulation of a General Parity-Time-Symmetric Two-level System. EPL 2018, 123, 40002. [Google Scholar] [CrossRef]
Wen, J.; Zheng, C.; Kong, X.; Wei, S.; Xin, T.; Long, G.L. Experimental Demonstration of a Digital Quantum Simulation of a General PT-symmetric System. Phys. Rev. A 2019, 99, 062122. [Google Scholar] [CrossRef]
Gao, W.-C.; Zheng, C.; Liu, L.; Wang, T.-J.; Wang, C. Experimental simulation of the parity-time symmetric dynamics using photonic qubits. Opt. Exp. 2021, 29, 517–526. [Google Scholar] [CrossRef]
Wen, J.; Zheng, C.; Ye, Z.; Xin, T.; Long, G.L. Stable states with nonzero entropy under broken PT-symmetry. Phys. Rev. Res. 2021, 3, 013256. [Google Scholar] [CrossRef]
Zheng, C. Quantum simulation of PT-arbitrary-phase–symmetric systems. EPL 2021, 136, 30002. [Google Scholar] [CrossRef]
Lee, T.D.; Wick, G.C. Negative Metric and the Unitarity of the S Matrix. Nucl. Phys. B 1969, 9, 209–243. [Google Scholar] [CrossRef]
Mostafazadeh, A. Pseudo-Hermiticity versus PT Symmetry: The Necessary Condition for the Reality of the Spectrum of a Non-Hermitian Hamiltonian. J. Math. Phys. 2002, 43, 205–243. [Google Scholar] [CrossRef]
Mostafazadeh, A. Pseudo-Hermiticity versus PT-symmetry. II. A complete characterization of non-Hermitian Hamiltonians with a real spectrum. J. Math. Phys. 2002, 43, 2814–2816. [Google Scholar] [CrossRef]
Konotop, V.V.; Yang, Z.; Zezyulin, D.A. Nonlinear Waves in PT-Symmetric Systems. Rev. Mod. Phys. 2016, 88, 035002. [Google Scholar] [CrossRef] [Green Version]
Mostafazadeh, A. Pseudo-Hermiticity versus PT-Symmetry III: Equivalence of Pseudo-Hermiticity and the Presence of Antilinear Symmetries. J. Math. Phys. 2002, 43, 3944–3951. [Google Scholar] [CrossRef]
Mostafazadeh, A. Pseudo-Hermiticity and generalized PT- and CPT-symmetries. J. Math. Phys. 2003, 44, 974–989. [Google Scholar] [CrossRef] [Green Version]
Zheng, C.; Tian, J.; Li, D.; Wen, J.; Wei, S.; Li, Y.-S. Efficient quantum simulation of an anti-P-pseudo-Hermitian two-level system. Entropy 2020, 22, 812. [Google Scholar] [CrossRef] [PubMed]
Solombrino, L. Weak Pseudo-Hermiticity and Antilinear Commutant. J. Math. Phys. 2002, 43, 5439–5445. [Google Scholar] [CrossRef] [Green Version]
Nixon, S.; Yang, J. All-real spectra in optical systems with arbitrary gain-and-loss distributions. Phys. Rev. A 2016, 93, 031802(R). [Google Scholar] [CrossRef] [Green Version]
Mostafazadeh, A. Time-Dependent Pseudo-Hermitian Hamiltonians and a Hidden Geometric Aspect of Quantum Mechanics. Entropy 2020, 22, 471. [Google Scholar] [CrossRef] [Green Version]
Pinske, J.; Teuber, L.; Scheel, S. Holonomic Gates in Pseudo-Hermitian Quantum Systems. Phys. Rev. A 2019, 100, 042316. [Google Scholar] [CrossRef] [Green Version]
Chu, Y.; Liu, Y.; Liu, H.; Cai, J. Quantum Sensing with a Single-Qubit Pseudo-Hermitian System. Phys. Rev. Lett. 2020, 124, 020501. [Google Scholar] [CrossRef] [Green Version]
Jin, L. Unitary Scattering Protected by Pseudo-Hermiticity. Chin. Phys. Lett. 2022, 39, 037302. [Google Scholar] [CrossRef]
Feynman, R. Simulating Physics with Computers. Int. J. Theor. Phys. 1982, 21, 467–488. [Google Scholar] [CrossRef]
Greiner, M.; Mandel, O.; Esslinger, T.; Hansch, T.W.; Bloch, I. Quantum Phase Transition from a Superfluid to a Mott Insulator in a Gas of Ultracold Atoms. Nature 2002, 415, 39–44. [Google Scholar] [CrossRef] [PubMed]
Leibfried, D.; DeMarco, B.; Meyer, V.; Rowe, M.; Ben-Kish, A.; Britton, J.; Itano, W.M.; Jelenkovic, B.; Langer, C.; Rosenband, T.; et al. Trapped-Ion Quantum Simulator: Experimental Application to Nonlinear Interferometers. Phys. Rev. Lett. 2002, 89, 247901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Friedenauer, A.; Schmitz, H.; Glueckert, J.T.; Porras, D.; Schaetz, T. Simulating a Quantum Magnet with Trapped Ions. Nat. Phys. 2008, 4, 757–761. [Google Scholar] [CrossRef]
Kim, K.; Chang, M.-S.; Korenblit, S.; Islam, R.; Edwards, E.E.; Freericks, J.K.; Lin, G.-D.; Duan, L.-M.; Monroe, C. Quantum Simulation of Frustrated Ising Spins with Trapped Ions. Nature 2010, 465, 590–593. [Google Scholar] [CrossRef] [PubMed]
Lanyon, B.P.; Whitfield, J.D.; Gillett, G.G.; Goggin, M.E.; Almeida, M.P.; Kassal, I.; Biamonte, J.D.; Mohseni, M.; Powell, B.J.; Barbieri, M.; et al. Towards Quantum Chemistry on a Quantum Computer. Nat. Chem. 2010, 2, 106–111. [Google Scholar] [CrossRef] [Green Version]
Gerritsma, R.; Kirchmair, G.; Zahringer, F.; Solano, E.; Blatt, R.; Roos, C.F. Quantum Simulation of the Dirac Equation. Nature 2010, 463, 68–71. [Google Scholar] [CrossRef] [Green Version]
Georgescu, I.M.; Ashhab, S.; Nori, F. Quantum Simulation. Nature 2014, 86, 153–185. [Google Scholar] [CrossRef] [Green Version]
Setia, K.; Bravyi, S.; Mezzacapo, A.; Whitfield, J.D. Superfast Encodings for Fermionic Quantum Simulation. Phys. Rev. Res. 2019, 1, 033033. [Google Scholar] [CrossRef] [Green Version]
Aspuru-Guzik, A.; Walther, P. Photonic Quantum Simulators. Nat. Phys. 2012, 8, 285–291. [Google Scholar] [CrossRef] [Green Version]
Sheng, Y.B.; Zhou, L. Distributed Secure Quantum Machine Learning. Sci. Bull. 2017, 62, 1025–1029. [Google Scholar] [CrossRef] [Green Version]
Tranter, A.; Love, P.J.; Mintert, F.; Wiebe, N.; Coveney, P.V. Ordering of Trotterization: Impact on Errors in Quantum Simulation of Electronic Structure. Entropy 2019, 21, 1218. [Google Scholar] [CrossRef] [Green Version]
Ge, L.; Tureci, H.E. Antisymmetric PT-Photonic Structures with Balanced Positive- and Negative-index Materials. Phys. Rev. A 2013, 88, 053810. [Google Scholar] [CrossRef] [Green Version]
Hang, C.; Huang, G.; Konotop, V.V. PT Symmetry with a System of Three-Level Atoms. Phys. Rev. Lett. 2013, 110, 083604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Antonosyan, D.A.; Solntsev, A.S.; Sukhorukov, A.A. Parity-time anti-symmetric parametric amplifier. Opt. Lett. 2015, 40, 4575–4582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Wu, J.-H.; Artoni, M.; La Rocca, G.C. Parity-Time-Antisymmetric Atomic Lattices without Gain. Phys. Rev. A 2015, 91, 033811. [Google Scholar] [CrossRef] [Green Version]
Peng, P.; Cao, W.; Shen, C.; Qu, W.; Wen, J.; Jiang, L.; Xiao, Y. Anti-Parity-Time Symmetry with Flying Atoms. Nat. Phys. 2016, 12, 1139–1145. [Google Scholar] [CrossRef]
Yang, F.; Liu, Y.C.; You, L. Anti-PT Symmetry in Dissipatively Coupled Optical Systems. Phys. Rev. A 2017, 96, 053845. [Google Scholar] [CrossRef]
Choi, Y.; Hahn, C.; Yoon, J.W.; Song, H.S. Observation of an Anti-PT-Symmetric Exceptional Point and Energy-Difference Conserving Dynamics in Electrical Circuit Resonators. Nat. Commun. 2018, 9, 2182. [Google Scholar] [CrossRef]
Konotop, V.V.; Zezyulin, D.A. Odd-Time Reversal PT Symmetry Induced by an Anti-PT-Symmetric Medium. Phys. Rev. Lett. 2018, 120, 123902. [Google Scholar] [CrossRef] [Green Version]
Chuang, Y.-L.; Ziauddin; Lee, R.-K. Realization of Simultaneously Parity-Time-Symmetric and Parity-Time-Antisymmetric Susceptibilities along the Longitudinal Direction in Atomic Systems with all Optical Controls. Opt. Express 2018, 26, 21969–21978. [Google Scholar] [CrossRef]
Li, Y.; Peng, Y.-G.; Han, L.; Miri, M.-A.; Li, W.; Xiao, M.; Zhu, X.-F.; Zhao, J.; Alu, A.; Fan, S.; et al. Odd-Time Reversal PT Symmetry Induced by an Anti-PT-Symmetric Medium. Science 2019, 364, 170–173. [Google Scholar] [CrossRef] [PubMed]
Zheng, C. Duality quantum simulation of a generalized anti-PT-symmetric two-level system. EPL 2019, 126, 30005. [Google Scholar] [CrossRef]
Wen, J.; Qin, G.; Zheng, C.; Wei, S.; Kong, X.; Xin, T.; Long, G.L. Observation of information flow in the anti-PT-symmetric system with nuclear spins. Npj Quantum Inf. 2020, 6, 28. [Google Scholar] [CrossRef]
Zhang, H.; Peng, M.; Xu, X.-W.; Jing, H. Anti-PT-symmetric Kerr gyroscope. Chin. Phys. B 2022, 31, 014215. [Google Scholar] [CrossRef]
Long, G.-L. General Quantum Interference Principle and Duality Computer. Commun. Theor. Phys. 2006, 45, 825–843. [Google Scholar]
Long, G.-L.; Liu, Y. Duality Quantum Computing. Front. Comput. Sci. 2008, 2, 167. [Google Scholar] [CrossRef]
Long, G.-L.; Liu, Y.; Wang, C. Allowable Generalized Quantum Gates. Commun. Theor. Phys. 2009, 51, 65–67. [Google Scholar]
Long, G.L. Duality Quantum Computing and Duality Quantum Information Processing. Int. J. Theor. Phys. 2011, 50, 1305–1318. [Google Scholar] [CrossRef]
Cui, J.; Zhou, T.; Long, G.L. Density Matrix Formalism of Duality Quantum Computer and the Solution of Zero-Wave-Function Paradox. Quantum Inf. Process. 2012, 11, 317–323. [Google Scholar] [CrossRef]
Qiang, X.; Zhou, X.; Wang, J.; Wilkes, C.M.; Loke, T.; O’Gara, S.; Kling, L.; Marshall, G.D.; Santagati, R.; Ralph, T.C.; et al. Large-Scale Silicon Quantum Photonics Implementing Arbitrary Two-Qubit Processing. Nat. Photon. 2018, 12, 534–539. [Google Scholar] [CrossRef] [Green Version]
Wei, S.-J.; Li, H.; Long, G.-L. A Full Quantum Eigensolver for Quantum Chemistry Simulations. Research 2020, 2020, 1486935. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Shao, C.P.; Li, Y.; Li, H.B. Quantum Algorithm Design: Techniques and Applications. J. Syst. Sci. Complex. 2019, 32, 375–495. [Google Scholar] [CrossRef]
Kato, T. Perturbation Theory for Linear Operators; Springer: Berlin, Germany, 1966. [Google Scholar]
Neeley, M.; Ansmann, M.; Bialczak, R.C.; Hofheinz, M.; Lucero, E.; O’Connell, A.D.; Sank, D.; Wang, H.; Wenner, J.; Cleland, A.N.; et al. Emulation of a Quantum Spin with a Superconducting Phase Qudit. Science 2009, 325, 722–725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Nielsen, M.A.; Chuang, I.L. Quantum Computation and Quantum Information, 10th Anniversary ed.; Cambridge University Press: Cambridge, UK, 2010; pp. 217–226. [Google Scholar]
Cory, D.G.; Price, M.D.; Havel, T.F. Nuclear Magnetic Resonance Spectroscopy: An Experimentally Accessible Paradigm for Quantum Computing. Physics D 2008, 120, 82–101. [Google Scholar] [CrossRef] [Green Version]
Hu, S.W.; Xue, K.; Ge, M.L. Optical Simulation of the Yang-Baxter Equation. Phys. Rev. A 2008, 78, 022319. [Google Scholar] [CrossRef] [Green Version]
Knill, E.; Laflamme, R.; Milburn, G. A Scheme for Efficient Quantum Computation with Linear Optics. Nature 2001, 409, 46–52. [Google Scholar] [CrossRef]
Cerf, N.J.; Adami, C.; Kwiat, P.G. Optical Simulation of Quantum Logic. Phys. Rev. A 1998, 57, R1477(R). [Google Scholar] [CrossRef] [Green Version]
Fan, C.-R.; Lu, B.; Feng, X.-T.; Gao, W.-C.; Wang, C. Efficient multi-qubit quantum data compression. Quantum Eng. 2021, 3, e67. [Google Scholar] [CrossRef]

Figure 1. Schematic circuit for the quantum simulation of a PH-

φ

system based on LCU. The whole system consists of

n_{1}

ancillary qubits and

n_{2}

evolutionary qubits, and it will pass the quantum circuit from the left to the right. The system is initialized in

{| 0 〉}_{a} {| 0 〉}_{e}

at first, and then the evolutionary qubits are prepared in arbitrary state

{| ψ 〉}_{e}

as demanded. After being operated by a unitary rotation

U_{E 1}

, m-controlled operations (i.e., 0-controlled

U_{0}

, 1-controlled

U_{1}

, …, and

(m - 1)

-controlled

U_{m - 1}

), and a single-qudit rotation

U_{E 2}

, the evolutionary subsystem will evolve as per Equation (4) if the ancillary subsystem is measured in state

{| 0 〉}_{a}

Figure 1. Schematic circuit for the quantum simulation of a PH-

φ

system based on LCU. The whole system consists of

n_{1}

ancillary qubits and

n_{2}

evolutionary qubits, and it will pass the quantum circuit from the left to the right. The system is initialized in

{| 0 〉}_{a} {| 0 〉}_{e}

at first, and then the evolutionary qubits are prepared in arbitrary state

{| ψ 〉}_{e}

as demanded. After being operated by a unitary rotation

U_{E 1}

, m-controlled operations (i.e., 0-controlled

U_{0}

, 1-controlled

U_{1}

, …, and

(m - 1)

-controlled

U_{m - 1}

), and a single-qudit rotation

U_{E 2}

, the evolutionary subsystem will evolve as per Equation (4) if the ancillary subsystem is measured in state

{| 0 〉}_{a}

Figure 2. Sets of the P-pseudo-Hermitian-

φ

-symmetric, PT-

φ

-symmetric, P-pseudo-Hermitian, and Hermitian systems. The green ellipse of the P-PH set is in the blue ellipse of the general P-PH-

φ

set because the former can be seen as a special case where

φ = 2 k π

(k is integral) of the latter. The P-PH anti-symmetric set is also a subset of the P-PH-

φ

set when

φ = (2 k + 1) π

(k is integral). The sets relevant to the P-pseudo-Hermiticity and the PT symmetry are different, though they have intersections with each other. Notice that the Hermitian sets are not in the P-PH-

φ

set only. Other NH sets may include various

η

-PH-

φ

sets (

η

is other than P), sets relevant to unknown symmetries, and so on.

Figure 2. Sets of the P-pseudo-Hermitian-

φ

-symmetric, PT-

φ

-symmetric, P-pseudo-Hermitian, and Hermitian systems. The green ellipse of the P-PH set is in the blue ellipse of the general P-PH-

φ

set because the former can be seen as a special case where

φ = 2 k π

(k is integral) of the latter. The P-PH anti-symmetric set is also a subset of the P-PH-

φ

set when

φ = (2 k + 1) π

(k is integral). The sets relevant to the P-pseudo-Hermiticity and the PT symmetry are different, though they have intersections with each other. Notice that the Hermitian sets are not in the P-PH-

φ

set only. Other NH sets may include various

η

-PH-

φ

sets (

η

is other than P), sets relevant to unknown symmetries, and so on.

Figure 3. Three-qubit quantum circuit. The system is initialized to

{| 00 〉}_{a} {| 0 〉}_{e}

, and the evolutionary qubit e can be rotated to

{| ψ 〉}_{e}

R_{ψ}

as needed. In the first block, operations prepare the six-dimensional subspace, and assign the three UE parameters. In the second block, three controlled-controlled operators (the first dashed one can be removed) generate the UE terms. In the third block, operations are applied on the ancillary system to superpose the three UE terms in Equation (8). Finally, quantum measurements are performed on the ancillary system to evolve the qubit e as the P-PH-

φ

system in an indeterministic way if

{| 00 〉}_{a}

is output.

Figure 3. Three-qubit quantum circuit. The system is initialized to

{| 00 〉}_{a} {| 0 〉}_{e}

, and the evolutionary qubit e can be rotated to

{| ψ 〉}_{e}

R_{ψ}

φ

system in an indeterministic way if

{| 00 〉}_{a}

is output.

Figure 4. Two-qubit quantum circuit. The system includes an ancillary qubit a and an evolutionary qubit e, and is initialized to

{| 0 〉}_{a} {| 0 〉}_{e}

at first. Then, the qubit e is rotated to

{| ψ 〉}_{e}

as needed by

R_{ψ}

. In the main part, operators are applied in series (i.e., a single-qubit rotation W, two controlled operators, and a Hadamard). Finally, the evolutionary qubit e will evolve as

e^{- i \frac{t}{ℏ} H_{φ}} {| ψ 〉}_{e}

, if the ancillary qubit is measured in state

{| 0 〉}_{a}

Figure 4. Two-qubit quantum circuit. The system includes an ancillary qubit a and an evolutionary qubit e, and is initialized to

{| 0 〉}_{a} {| 0 〉}_{e}

at first. Then, the qubit e is rotated to

{| ψ 〉}_{e}

as needed by

R_{ψ}

. In the main part, operators are applied in series (i.e., a single-qubit rotation W, two controlled operators, and a Hadamard). Finally, the evolutionary qubit e will evolve as

e^{- i \frac{t}{ℏ} H_{φ}} {| ψ 〉}_{e}

, if the ancillary qubit is measured in state

{| 0 〉}_{a}

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Zheng, C. Quantum Simulation of Pseudo-Hermitian-φ-Symmetric Two-Level Systems. Entropy 2022, 24, 867. https://doi.org/10.3390/e24070867

AMA Style

Zheng C. Quantum Simulation of Pseudo-Hermitian-φ-Symmetric Two-Level Systems. Entropy. 2022; 24(7):867. https://doi.org/10.3390/e24070867

Chicago/Turabian Style

Zheng, Chao. 2022. "Quantum Simulation of Pseudo-Hermitian-φ-Symmetric Two-Level Systems" Entropy 24, no. 7: 867. https://doi.org/10.3390/e24070867

APA Style

Zheng, C. (2022). Quantum Simulation of Pseudo-Hermitian-φ-Symmetric Two-Level Systems. Entropy, 24(7), 867. https://doi.org/10.3390/e24070867

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Quantum Simulation of Pseudo-Hermitian-φ-Symmetric Two-Level Systems

Abstract

1. Introduction

2. Complex Generalization of Pseudo-Hermitian Symmetry

3. Quantum Simulation Using LCU by Duality Quantum Computing

4. Quantum Simulation of $P$ -Pseudo-Hermitian- $φ$ -Symmetric Two-Level Systems

4.1. P-PH- $φ$ Two-Level Systems

4.2. UE of the Time-Evolutionary Operator

4.3. Qubit Simulation

5. Experimental Proposals

6. Conclusions

Funding

Data Availability Statement

Conflicts of Interest

Abbreviations

Appendix A

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI

Article Menu

Quantum Simulation of Pseudo-Hermitian-φ-Symmetric Two-Level Systems

Abstract

1. Introduction

2. Complex Generalization of Pseudo-Hermitian Symmetry

3. Quantum Simulation Using LCU by Duality Quantum Computing

4. Quantum Simulation of P -Pseudo-Hermitian- φ -Symmetric Two-Level Systems

4.1. P-PH- φ Two-Level Systems

4.2. UE of the Time-Evolutionary Operator

4.3. Qubit Simulation

5. Experimental Proposals

6. Conclusions

Funding

Data Availability Statement

Conflicts of Interest

Abbreviations

Appendix A

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI

4. Quantum Simulation of $P$ -Pseudo-Hermitian- $φ$ -Symmetric Two-Level Systems

4.1. P-PH- $φ$ Two-Level Systems