Rotation-time symmetry in bosonic systems and the existence of exceptional points in the absence of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{PT}}$$\end{document}PT symmetry

We study symmetries of open bosonic systems in the presence of laser pumping. Non-Hermitian Hamiltonians describing these systems can be parity-time (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathscr{PT}}}$$\end{document}PT) symmetric in special cases only. Systems exhibiting this symmetry are characterised by real-valued energy spectra and can display exceptional points, where a symmetry-breaking transition occurs. We demonstrate that there is a more general type of symmetry, i.e., rotation-time (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{RT}}$$\end{document}RT) symmetry. We observe that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{RT}}$$\end{document}RT-symmetric non-Hermitian Hamiltonians exhibit real-valued energy spectra which can be made singular by symmetry breaking. To calculate the spectra of the studied bosonic non-diagonalisable Hamiltonians we apply diagonalisation methods based on bosonic algebra. Finally, we list a versatile set rules allowing to immediately identifying or constructing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{RT}}$$\end{document}RT-symmetric Hamiltonians. We believe that our results on the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr{RT}}$$\end{document}RT-symmetric class of bosonic systems and their spectral singularities can lead to new applications inspired by those of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathscr{PT}}}$$\end{document}PT-symmetric systems.

We also find general rules for constructing R T-symmetric symmetric non-Hermitian Hamiltonians. This provides a framework for exploring the physics of singular energy spectra in terms of symmetries in a range of bosonic systems including non-linear interactions between modes, associated, e.g., with unconventional photon blockades or second-and third-harmonic generation processes 51,[53][54][55][56] . Consequently, future experiments could explore properties of spectral singularities in areas of quantum and atom optics that have never been studied before.

Results
PT symmetry in bosonic system. Let us begin by considering a P T-symmetric Hamiltonian describing a damped bosonic linear system with a classical laser drive.
Thus, let us choose a system depicted in Fig. 1 and described with Hamiltonian where a and b are the annihilation operators for two modes, respectively, g is the coupling strength between these two modes, κ ≡ κ a = κ b is the field loss/gain rate of the a/b mode and ε is the strength of the modes drive. We assume for the sake of simplicity that g and ε are real and positive. It is known that a Hamiltonian H is regarded as P T-symmetric if it commutes with the operator P T . Since P T is a reflection, i.e., P T = (P T ) −1 , we can rewrite the condition for H to be P T-symmetric as (P T )H(P T ) = H . So, if the transformation does not change H, then H is P T-symmetric. Now we want to check whether H 1 has P T symmetry properties. For that purpose, it is convenient to know how P T transforms bosonic operators. For a single field mode, one can easily derive all needed formulas knowing the effect of the space-reflection operator (also known as the parity operator) P and the time-reversal operator T , where [P , T ] = 0 , on the position operator x and the momentum operator p : PxP = −x , TxT =x , PpP = −p , and TpT = −p . From the form of the position operator x = (a + a † )/ √ 2 and the momentum operator p = i(a † − a)/ √ 2 , we can infer the result of applying the symmetry transformation to the bosonic field, i.e., (P T )a(P T ) = −a , (P T )a † (P T ) = −a † and (P T )i(P T ) = −i 37,38 . We have used [P , T ] = 0 for derivations. In the case of two-mode fields, these transformations of both field-mode annihilation operators are sometimes sufficient to reveal P T symmetry manifesting itself as real-valued energy spectra (e.g., see Refs. 37,38 ). Unfortunately, it is not the case of Hamiltonian (1). One can easily check using P = exp[iπ(a † a + b † b)] 37,38 that (P T )H 1 (P T ) = H 1 because of the last two terms of H 1 , which describe loss in the mode A and gain in the mode B. The gain-loss terms are, however, crucial for the existence of EPs, because non-Hermiticity is necessary for the emergence of EPs 9,57 . Moreover, modelling losses is unavoidable as they are present in all real systems. Therefore, to properly account for the symmetry between loss and gain we can, e.g., modify the spacereflection operator by multiplying it by the exchange operator P S 58 , which can be interpreted as exchanging the modes spatially (i.e., a ↔ b ). Since there are only two modes, the exchange operator acts in the same way as a permutation operator or a cyclic-shift operator, and thus, we can use all these terms interchangeably here. A matrix representation of P S is given by a perfect shuffle 59 . Hence, we use here the modified parity operator given by P = P S P and one can easily check that ( P T )H 1 ( P T ) = H 1 , so H 1 is P T-symmetric. Note that P S does not change any important features of the parity operator because P S is also a reflection operator satisfying P 2 S = 1 , [P S , P ] = 0 , and [P S , T ] = 0 . Thus, P T is also a reflection and [ P , T ] = 0 . Using these features of P T and applying the Cayley-Hamilton theorem, one can prove that the parameters of the characteristic equation of P T -symmetric Hamiltonians are real 38 . Consequently, the eigenvalues of H 1 may only appear as complex-conjugate pairs or are real. Moreover, one can show that an eigenvalue of H 1 is real if the corresponding eigenstate of H 1 is also an eigenstate of P T . If all the eigenfunctions are identical for both operators, then the entire spectrum of H 1 is real. In such a case H 1 is said to have unbroken P T symmetry. Otherwise, the P T symmetry is broken and then, complex-conjugate pairs of eigenvalues appear. It is obvious that values of various system parameters determine whether the symmetry of H 1 is broken. For example, for κ = 0 , the entire spectrum is real because www.nature.com/scientificreports/ then H 1 is Hermitian. However, one can expect that large enough κ values result in symmetry breaking. We can easily illustrate this feature of H 1 in the case when ε = 0 , i.e., when lasers are turned off, because the solution for such a Hamiltonian is well known 5 . For ε = 0 , we can diagonalise the Hamiltonian by introducing new bosonic operators, which are superpositions of the old ones, and which describe eigenmodes. Eigenvalues are then given by n , where = g 2 − κ 2 and n is the difference between the numbers of excitations in the first and second eigenmode.
In Fig. 2, we plot real and imaginary parts of three eigenvalues (corresponding to n = −1, 0, 1 ) as functions of κ for g = 1 . One can see that for κ < g all these eigenvalues are real. However, for κ > g , one eigenvalue is real, while the two other make a complex-conjugate pair. One can also see that exactly at the point κ = g a symmetry-breaking transition takes place. All eigenvalues have there the same real and imaginary parts -they are identical. In such points not only the eigenvalues are the same, but also the corresponding eigenvectors are parallel, and therefore, the Hamiltonian is not diagonalisable 3,41 . Such points in the parameter space are known as EPs. Since some of eigenvectors are parallel in EP, Hermitian Hamiltonians cannot exhibit any EP 9,57 . The Hamiltonian (1) is non-Hermitian, so in Fig. 2 we can observe an EP.
Spectral singularities are of great importance because non-Hermitian systems in their vicinity can exhibit phenomena, which cannot be observed in Hermitian systems. It is worth mentioning that in case of two or more dimensional parameter spaces it possible to observe not single points but exceptional lines or even exceptional surfaces 60,61 . Exceptional surfaces have similar properties as EPs, but are more stable and easier achievable in experiments 62 .
The effect of laser pumping on a spectral singularity. We are going to investigate the influence of a non-zero ε on the spectrum of the Hamiltonian H 1 and on its spectral singularity. For this purpose, we diagonalise H 1 by expressing it in terms of the bosonic operators c ε , c + ε , d ε , and d + ε (see Methods), i.e., where 0 = −2gε 2 / 2 . Let us restrict ourselves to considering only three different eigenvalues -one of them corresponding to the cases, in which there are equal numbers of excitations in both c ε and d ε modes: and the two others corresponding to the cases, in which there is one excitation more in one mode than in the other: (2)  www.nature.com/scientificreports/ As one can see from Fig. 3, the effect of laser pumping on the spectrum is to shift real parts of all the eigenvalues by the same value, i.e., 0 . The absolute value of the shift increases with the increase of κ/g ratio. As a result we observe an unusual spectral singularity. All three eigenvalues tend to the same value in their real and imaginary parts as κ → g , but unlike typical EP, in the case of this spectral singularity, all eigenvalues go to −∞ . It is also seen from Fig. 3 that a non-zero ε does not influence the imaginary parts of the eigenvalues.
However, this is true only for the exact resonant case considered here. If we introduce a detuning of the two modes by the substitution κ → κ + i δ/2 then we can observe that also the imaginary parts of the eigenvalues are shifted. In order to show this we present Riemann surfaces of the eigenvalues of H 1 in Fig. 4. Note that despite of the presence of this unusual spectral singularity, eigenvalues behave as we expect from a P T-symmetric system. There are ranges of κ , for which eigenvalues are real, and such ranges that a complex-conjugate pair appear. Spectral singularity marks the boundary between these two ranges and in its vicinity we can observe phenomena, which cannot be observed in Hermitian systems. For instance, one can check that we can obtain the enhancement of sensing, since eigenvalue splitting is given by �E = E + − E 0 ≈ 2κ�g , where �g = g − κ can be interpreted as a perturbation strength.
We conclude this section by noting that the lasers pumping process does not change any typical feature of P T-symmetric systems except that eigenvalues tend to ±∞ while going to coalescence. To our knowledge, this kind of behaviour was not predicted yet even in other systems.
Spectral singularities beyond the PT symmetry. Remarkably, spectral singularities can be found also in non-Hermitian systems that are not P T-symmetric 9,57 . Therefore now the question arises: can a given bosonic non-P T-symmetric systems display any singularity and how to engineer such systems? In order to answer this question let us consider the following Hamiltonian It is seen that H 2 is not P T-symmetric. Nevertheless, eigenvalues of H 2 are exactly the same as those of H 1 , because U −1 H 2 U = H 1 and U = exp[iπ/2(a † a + b † b)] is a unitary transformation. We calculate eigenvalues of H 2 in the section Methods. Therefore it is reasonable to expect that H 2 has some other type of symmetry. In order to find this symmetry, we note that the reflection operator P = exp[iπ(a † a + b † b)] describes, in fact, the rotation of the frame by an angle π . Zhang et al. have shown that in fermionic systems the rotation-time www.nature.com/scientificreports/ symmetry is a generalisation of the parity-time symmetry 43,44 . Therefore, we can expect that the symmetry we are looking for in our bosonic system is also the rotation-time symmetry. Let us introduce the rotation operator given by where θ is an arbitrary angle and let us see whether the R T symmetry leads to the same conclusions as those for the P T symmetry. We need to apply again the permutation operator P S because of the gain-loss terms. So, we define R = P S R.
Note that R has different properties than P , i.e., it is not a reflection operator, i.e., R −1 � = R , and it does not commute with the time-reflection operator T , i.e., [R , T ] � = 0 . It seems to be an obstacle. Fortunately, the operator R T is a reflection operator: Using the above property of R T , we are going to show now that a non-Hermitian Hamiltonian H, which commutes with R T , has real energy spectra, if eigenstates of R T are also the eigenstates of H. To this end, let us first note that, from [H, R T ] = 0 and ( R T ) 2 = 1 , we can infer that the characteristic equation of H is real 38 , and thus, eigenvalues of H can only be real or appear as complex-conjugate pairs.
Let us demonstrate now that eigenvalues of the operator R T , which we denote by χ , cannot be equal to zero. In fact, they are just phase factors, i.e., the absolute value of each eigenvalue is equal to unity. The eigenvalue equation reads Multiplying Eq. (8) on the left by R T and using the fact that ( R T ) 2 = 1 we obtain Now we can use Eqs. (8) and (7), and T = T −1 to make further transformation: and, then, using T iT = −i , we finally obtain This concludes the proof that χ cannot be equal to zero.
Next, let us show that the eigenvalue E of a Hamiltonian H is real if | � is an eigenstate of both operators, H and R T . We begin with the time-independent Schrödinger equation Multiplying it on the left by R T , once more recalling ( R T ) 2 = 1 and using assumption that | � is also an eigenstate of R T we obtain   Since χ is nonzero, pure phase, we deduce that E = E * , which means that the eigenvalue E is real. Hence, the symmetry described by the rotation-time operator R T leads to the same conclusions as those for P T-symmetry: (i) eigenvalues of R T-symmetric Hamiltonians are real if corresponding eigenstates are the same for both R T and H; (ii) or eigenvalues appear as complex-conjugate pairs, when the symmetry is broken. Most importantly, we can find EPs at the boundary between the broken and unbroken R T symmetry regions in the parameter space.
We conclude this section noting that H 2 is R T-symmetric in the rotation operator (6) with θ = 0.

RT-symmetric bosonic systems. The R T symmetry can play an important role in searching for
bosonic systems with EPs. Note that the R T symmetry includes much larger class of bosonic systems than the P T symmetry because operator P T is just one specific instance (corresponding to θ = π ) of all R T operators. To give an example of a class of R T-symmetric Hamiltonians with potential EPs, let us consider a Hamiltonian, which is a generalisation of both previously analysed Hamiltonians H 1 and H 2 : The Hamiltonian H 3 has the same eigenvalues as H 1 , since they are related by a unitary transformation. The Hamiltonians H 1 , H 2 and H 3 describe the same quantum system but in different frames. The effect of an action of the R T operator on the Hamiltonian can be quickly inferred using the formulas: = H 3 , if we choose θ = −2φ. It is straightforward to give expressions in their most general form, which can be a part of the Hamiltonian that governs a bosonic R T-symmetric system. Each of them guarantees itself the R T invariance.
The most relevant Hamiltonian terms having this property are presented in Table 1. By utilising these R T -invariant terms, one can compose a wide class of R T-symmetric models, which can exhibit exceptional points.

Discussion
In summary, we have studied the R T symmetry and its applications, and conditions for its presence in a laser pumped bosonic system with losses and gain. To our knowledge, this is the first proposal of adopting this type of symmetry to quantum optics. We have shown that the P T symmetry is only a special case of a rotation-time symmetry. Our study of the effect of laser pumping on symmetries and symmetry phase transitions resulted in discovering a new type of spectral singularity (as shown in Fig. 3). This was possible by creating a new types of operators to diagonalise a Hamiltonian describing a bosonic system with a classical laser drive (see Methods). There, we have obtained real energy values without applying requirements of the P T symmetry, which turned out to be too restrictive for many types of photonic systems. We have created a versatile framework for certifying the R T symmetry based on a set of expressions that guarantee the R T invariance. We believe that this work can significantly contribute to quantum optics and we hope that can open opportunities to study new classes of systems and related physical effects. Table 1. Different types of R T-invariant terms for a specific choice of R , i.e., θ = −2φ . Parameters n, m, j, and l are arbitrary natural numbers.
The Hamiltonian H 1 can be diagonalised in the same way but with one additional step. H 1 in terms of the operators c, c + , d, and d + is given by In this additional step we introduce new operators: ic → c , −ic + → c + , id → d , and −id + → d + . These new operators satisfy the same commutation relations as original operators. In this way we transform H 1 to the form given by Eq. (19). Then we apply transformations given by the displacement operators to obtain the same eigenvalues as for H 2 .