Advances in Mathematical Physics

Volume 2018, Article ID 3653851, 9 pages

https://doi.org/10.1155/2018/3653851

## Loss, Gain, and Singular Points in Open Quantum Systems

^{1}Institute for Quantum Science and Engineering, Texas A&M University, College Station, Texas 77843, USA^{2}Department of Applied Science and Mathematics, College of Arts and Sciences, Abu Dhabi University, Abu Dhabi, UAE^{3}Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany

Correspondence should be addressed to Ingrid Rotter; ed.gpm.skp@rettor

Received 27 June 2017; Accepted 19 April 2018; Published 2 September 2018

Academic Editor: Soheil Salahshour

Copyright © 2018 Hichem Eleuch and Ingrid Rotter. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Non-Hermitian quantum physics is used successfully for the description of different puzzling experimental results, which are observed in open quantum systems. Mostly, the influence of exceptional points on the dynamical properties of the system is studied. At these points, two complex eigenvalues of the non-Hermitian Hamiltonian coalesce (where is the energy and is the inverse lifetime of the state ). We show that also the eigenfunctions of the two states play an important role, sometimes even the dominant one. Besides exceptional points, other critical points exist in non-Hermitian quantum physics. At these points in the parameter space, the biorthogonal eigenfunctions of become orthogonal. For illustration, we show characteristic numerical results.

#### 1. Introduction

The properties of open quantum systems are described usually by averaging over the spectroscopic features of the individual states of the system. Open quantum systems are characterized by, e.g., their Markovian and non-Markovian behavior or by equilibrium and nonequilibrium properties. The present-day high-resolution experimental studies provide however more information, namely, concrete information on the spectroscopic properties of the individual states. It is a challenge for the theory to describe these results.

In many current theoretical studies, a non-Hermitian formalism is used for the description of individual states of an open quantum system. In these papers, the part of the system, which is localized in a certain finite space region, is considered to be embedded into an infinitely extended environment. Mathematically, the total function space consists of two parts: the localized part of the system and the extended environment. Mostly, the environment is assumed to be the continuum of scattering wavefunctions^{1}.

This formalism is, from the point of view of mathematics, much more complicated than the familiar Hermitian formalism. For example, the eigenfunctions of a non-Hermitian operator are biorthogonal with nontrivial consequences for their normalization. Most important is however the existence of singularities.

In mathematics, the so-called exceptional points (EPs) are known for a long time; see the book [1] on the perturbation theory for linear operators. At the EPs, two eigenvalues of coalesce. This coalescence is, on principle, different from the well-known degeneration of two eigenstates of a Hermitian operator . In the last case, the eigenfunctions of the two states are different and orthogonal while those of two coalescing eigenstates of are biorthogonal and differ only by a phase from one another at the EP. The related interesting topological features are studied also experimentally in detail [2–4].

Due to the important role of EPs in non-Hermitian systems, in most theoretical studies of non-Hermitian quantum physics only the eigenvalues of the non-Hermitian Hamiltonian are considered in detail. The eigenfunctions are often assumed to be of minor interest. Examples are the method of complex scaling [5] and the review [6] on superradiance. The meaning of the eigenvalues as well as of the eigenfunctions of is discussed theoretically in the review [7] and compared to experimental results in the review [8]. The eigenfunctions around EPs play an important role for the explanation of some experimental results that are puzzling from the point of view of Hermitian quantum physics. Examples are the so-called phase lapses [9] observed experimentally in mesoscopic systems and the high efficiency of the photosynthesis [10]. The role of EPs in quantum physics is discussed in many other papers [11–16].

The meaning of the eigenvalues and that of the eigenfunctions of are less different from one another than it seems to be. On the one hand, the eigenfunctions of play a role also in those cases in which they are usually not considered [7, 8]. Among others, they are responsible for the fact that the influence of the EPs occurs in a larger parameter range around the EP than expected from the behavior of the eigenvalues. Moreover, higher-order EPs are shielded from observation in physical systems [17] because of the mixing of the eigenfunctions of via the environment. On the other hand, the eigenvalues of determine the position of the EPs and influence therefore also the behavior of the eigenfunctions of .

Further studies in non-Hermitian quantum physics have shown the following nontrivial results. The coupling between system and environment may occur via an exchange of particles or of information (mostly in terms of excitons). We have the following in detail.

(i) Every state of the system is related to different decay channels according to the values of the so-called partial widths of the state relative to the channels . The coupling of the whole system to the environment is characterized by the sum of partial widths. This case is realized in many-body systems. In standard Hermitian quantum physics, the exchange of particles between system and environment is described as tunneling of particles.

(ii) Every state of the system is coupled exclusively to only one channel with, respectively, some gain from the environment and some loss to the environment. This case is realized, e.g., in the photosynthesis in which visible light is captured in the light-harvesting complex [10]. The exciton transfer occurs without excitation of any eigenstates of . It cannot be described in standard theory.

Another problem of non-Hermitian quantum physics is the following. In mathematics, the EPs are defined in relation to one environment (mostly called channel). Physical systems are related however generally to more than one channel. For example, transmission of particles through a system needs at least two channels, entrance and exit channel. In this case, the different channels are parts of the total environment and are orthogonal to one another. Another example is the above-mentioned processes with gain and loss each of which occurs relative to another channel. The role played by EPs in physical systems is therefore not at all clear.

Consideration of the eigenfunctions of the non-Hermitian Hamiltonian in detail has shown further nontrivial properties of non-Hermitian quantum physics. For example, the EPs are not the only singular points in non-Hermitian quantum systems. Other singularities occur at maximum width bifurcation or level repulsion where the eigenfunctions of are orthogonal (instead of biorthogonal). They are related to obviously nonlinear irreversible processes [17].

In the present paper, we will illustrate the relation between the eigenvalues ^{2} and eigenfunctions of the non-Hermitian Hamiltonian , on the one hand, as well as the appearance of gain and loss and different types of singularities in non-Hermitian quantum physics, on the other hand.

For this aim, we start from a system with two states the complex energies of which are . We then consider the coupling of the system to two channels under the assumption that each state is coupled, above all, to only one channel: . We vary parametrically the energies of the states and determine the eigenvalues of the non-Hermitian Hamiltonian . In such a case, width bifurcation causes automatically a transition from loss (determined by and ) to gain (determined by and ).

The paper is organized in the following manner. In Section 2, we sketch the formalism which is used in performing the calculations. Some numerical results are shown in Section 3. First we show a few eigenvalue trajectories together with the corresponding eigenfunction trajectories in Section 3.1. The resonance structure of the transmission and the transmission probability are given in Section 3.2. In the last Section 4 we summarize and discuss the results and draw some conclusions on non-Hermitian quantum physics of open quantum systems.

#### 2. Non-Hermitian Hamiltonian of an Open Quantum System

##### 2.1. Eigenvalues and Eigenfunctions of

To begin with, we sketch the features typical for an open quantum system embedded in one common continuum. Details can be found in [7] and, above all, in [17]. They can be discussed by means of the genuine non-Hermitian matrix:Here, are the complex energies of the two states which are coupled to the environment (called also channel ). are the energies of the states and are their widths. stand for the coupling matrix elements of the two states via the common environment . They are complex where Re() arises from the principal value integral and Im() from the residuum [7]. The complex eigenvalues of give energies and widths of the states of the localized part of the system.

In (1), a singular point may appear, the so-called EP, at which two eigenvalues of coalesce [1]. At these points, the two corresponding eigenfunctions are not orthogonal. Insteadaccording to analytical and numerical results [18–20]. An EP is, according to its definition, related to the common environment in which the system is embedded. In other words, it is well defined under the condition that the system is embedded in only one continuum.

The eigenfunctions of a non-Hermitian Hamilton operator are biorthogonalIn the case of the symmetric Hamiltonian (1), it isand the eigenfunctions should be normalized according toin order to smoothly describe the transition from a closed system with discrete states to a weakly open one with narrow resonance states. As a consequence of (5), the values of the standard expressions are changed:Furthermore, the phase rigidity, which is a quantitative measure for the biorthogonality of the eigenfunctions,is smaller than 1. Far from an EP, while it approaches the value when an EP is approached.

Additionally to the Hamiltonian (1), we will consider the non-Hermitian matrixwhich describes the system without any mixing of its states via the environment. In other words, corresponds to vanishing EM of the eigenstates. The eigenfunctions of can be represented in the set of eigenfunctions of ,under the condition that are normalized by . The coefficients differ from the . They contain the information on the strength of EM.

The main features characteristic of open quantum systems are described well by the eigenvalues and eigenfunctions of (1). Typical values related to the eigenfunctions are the phase rigidity (9) and the contribution of EM (11) to their purity. All these values contain the influence of the environment. They are proven experimentally; for details see [7, 8, 17].

##### 2.2. Schrödinger Equation with

The Schrödinger equation may be rewritten into a Schrödinger equation with source term [7, 17],In this representation, the coupling of the states and of the localized system via the common environment of scattering wavefunctions (EM) is contained in the source term of the Schrödinger equation; for details see [7].

Far from EPs, the coupling of the localized system to the environment influences the spectroscopic properties of the system only marginally [7, 17]. The influence of the environment is however nonvanishing also in this case; see, e.g., the experimental results [21]. That means, the eigenstates of will never be the same as those of a Hermitian operator.

In the neighborhood of EPs, the coupling between system and environment cause, according to mathematical studies, nonlinear effects in the Schrödinger equation (12). Among others, these effects lead to a conservation of the resonance structure of the cross section in the one-channel case which is finally unaffected by EM and by the existence of EPs [17]. Thus, the one-channel case cannot be used in order to test the results of the non-Hermitian formalism.

##### 2.3. Two-Channel Case with

Let us consider now the genuine non-Hermitian matrixHere, and are the complex energies of the two states relative to the two different channels . Usually . and stand for the coupling matrix elements of the two states via environments and , respectively. The Hamiltonian (13) includes the fact that the complex energy of the localized state may be different relative to the different channels .

Also the eigenvalues and eigenfunctions of (13) are characterized by two numbers: the number of the state () of the localized part of the system and the number of the channel (), called environment, in which the system is embedded. Usually, and . Further, the wave functions and differ from one another due to the EM of the eigenstates via environments and , respectively. That means the system has four states, from a mathematical point of view.

The two environments are different from and orthogonal to one another. Further, the two eigenstates with equal and different arise from the same state of the localized part of the system. The zeros in the matrix (13) express the corresponding fact that the two eigenstates relative to the two channels 1 and 2 of cannot interact with one another.

We mention here that the Hamiltonian (13) is formally the same as the Hamiltonian (24) or rather in [10]. There is however a fundamental difference: in in [10], the transition of excitons is considered while (13) describes primarily the transition of particles. Furthermore, in in [10] the whole system is fully embedded into both environments () and () which are both of completely different nature and exist independently of one another. In contrast to this, the Hamiltonian (13) describes a system, the states of which are primarily embedded partially in each of the two different environments. These two different environments are also independent of one another. They are, however, nothing but parts of the total environment.

Without an EP in the considered parameter range in relation to both channels, we have , , and . Under the influence of an EP relative to and/or relative to , the eigenvalues and eigenfunctions will be, however, different from one another; , , and in the corresponding parameter range.

In analogy to (10), we will consider also the non-Hermitian Hamiltonian:which describes the system without any mixing of its states via any environment. In other words, corresponds to vanishing EM of the eigenstates via an environment. The mixing of the eigenstates of (13) can be represented in a set of eigenfunctions of (14) in complete analogy to the relation (11) for two states coupled to one common environment.

##### 2.4. Schrödinger Equation with

Using (14), we can write down the Schrödinger equation with source term for the two-channel case in analogy to (12) for the one-channel case. The corresponding equation readsThe source depends on the coupling of the system to both channels, i.e., on and on . It does not depend on the energies .

We repeat here that, according to their definition [1], EPs occur only in the one-channel case, i.e., only in the two submatrices related to channel and channel , respectively. They are not defined in the matrix (13). However, the EPs of the two submatrices in (13) may influence the dynamics of the open two-channel system.

##### 2.5. Non-Hermitian Hamiltonian and Resonance Structure of the Matrix

The resonance structure of the matrix containing two resonance states can be obtained from the expressionwhich is derived in [22] by rewriting the Breit-Wigner expression for one or more isolated resonances that are coupled to a common continuum of scattering wavefunctions; see also [7, 17]. The expression (16) is unitary. According to (16), the resonance structure of the cross section is determined exclusively by the spectroscopic values of the localized part of the system, i.e., by the eigenvalues of the non-Hermitian Hamiltonian . The expression (16) allows us therefore to receive reliable results also in the two-channel case and, moreover, when the phase rigidity of the eigenfunctions of is reduced, and when the eigenfunctions of contain EM, i.e., when they are mixed in the set of eigenfunctions of according to (11).

According to the results obtained in [17], the evolution of the system near to an EP is driven exclusively by the nonlinear source term of the Schrödinger equation (12) which describes the coupling of the localized part of the system to the common environment and is characteristic of the open quantum system embedded in one environment. The calculations in [17] are performed for fixed ; i.e., can not be responsible for the width bifurcation occurring in these calculations under the influence of an EP. The only conclusion is therefore that the nonlinear source term is the driving force. It is able, in the one-channel case, to largely conserve the resonance structure of the cross section.

The one-channel case does, therefore, not allow us to prove the existence of the nonlinear effects and of EM, since the resonance structure of the cross section calculated with and without EM is the same in this case [17]. This result agrees, on the one hand, with the experience received from many different numerical studies in realistic cases which are performed without taking into account EM. On the other hand, it is not in contradiction with the conclusions received from the study of non-Hermitian physics of open quantum systems.

The conservation of the resonance structure of the cross section, which is possible in the one-channel case, is expected to be impossible, generally, in the two-channel (or more-channel) case.

#### 3. Numerical Results

##### 3.1. Eigenvalues and Eigenfunctions of

The aim of our numerical studies is first to show the influence of a singularity onto the eigenvalues and eigenfunctions of the non-Hermitian Hamilton operator . Knowing these results, we study then, in the following Section 3.2, the influence of the singularity onto observable values such as resonance structure and contour plot of the transmission probability.

In Figure 1, we provide numerical results for the parameter dependence of eigenvalues and eigenfunctions of the Hamiltonian . We consider the case that one of the states is coupled to a special channel much stronger than the other one, . Without loss of generality, we choose the extreme case ; ; .