We develop a formulation of few-photon Fock-space waveguide transport that includes dissipation in the form of reservoir coupling. We develop the formalism for the case of a two-level atom and then show that our formalism leads to a simple rule that allows one to obtain the dissipative description of a system from the nondissipative case.
The topic of dissipation in quantum optics is both important and subtle. Dissipation is important because many quantum optical systems are open; therefore tracking all of the out flowing energy is often not possible, but accounting for it is often quite crucial. Dissipation is subtle because energy in quantum systems cannot simply disappear, and naively describing dissipation with phenomenological damping factors in the system’s equations of motion is formally incorrect, leading to a non-Hermitian Hamiltonian, violation of commutation relations [1], and time-reversal invariance [2].
A correct treatment of dissipation entails coupling the system under study, which usually contains discrete energy levels, to a reservoir that is modeled as having a continuum of energy states [1]. When this is done, the presence of dissipation in the system is accompanied by a source of fluctuation from the reservoir, enabling two-way energy exchange.
In quantum optics, the system–reservoir model is treated with many different approaches. In the Schrödinger picture, this model has been treated in the master equation approach, where multiple-time averages are calculated via the quantum regression theorem [2]. Alternatively in the Schrödinger picture one can use the stochastic approach of the Monte Carlo wave function [3], which involves evolving the wave function of the system with a non-Hermitian Hamiltonian term. In the Heisenberg picture, the input–output formalism [4] is widely used. In typical use of the input–output formalism, one assumes that the input is a weak coherent state. The nonlinear response is then often treated approximately in the weak excitation limit [5]. The input–output formalism is attractive because the external degrees of freedom, such as those of the reservoir, are integrated out, allowing one to focus on the dynamics of the system alone.
Sign up for Photonics Research TOC Get the latest issue of Advanced Photonics delivered right to you!Sign up now
Recently, there has been great interest in describing the properties of photons propagating in a waveguide coupled to quantum multilevel atoms either directly [6–16] or through the use of a microcavity [17,18]. The one-dimensional nature of waveguide propagation allows incident and scattered photons to interfere coherently, which can give rise to significant capabilities in controlling transport [6,7,11,16,19–23], switching [17,24], and amplification [9,15,25] processes for few-photon states.
A formulation of input–output formalism has been recently developed that enables analytical calculation of few-photon scattering matrices in these systems of waveguide coupled to quantum multilevel atoms. In this paper we extend this formalism to include dissipation. We model loss by introducing an auxiliary reservoir that transports energy away from the atom.
In applying the input–output formalism to treat two-photon transport [8,18], a key step is to insert a complete set of basis states for the single-excitation Hilbert space. Here we show that, to obtain the correct result for two-photon transport in the presence of dissipation, such a basis set needs to include the states of the reservoir.
The paper is organized as follows. In Section 2 we present the Hamiltonian describing a two-level atom with reservoir coupling. In Section 3 we obtain the one and two-photon S-matrices of the system using the input–output formalism, showing that excluding the reservoir from the system description leads to violation of outgoing particle exchange symmetry. Finally, in Section 4 we conclude with a general recipe for the inclusion of dissipation in a general class of important and relevant systems.
2. SYSTEM HAMILTONIAN
We consider a two-level atom side-coupled to a waveguide with right and left moving photons, as described by the Hamiltonian Here and create (annihilate) a right or left moving photon with group velocity , respectively. is the coupling rate between the atom and waveguide. The form of the interaction here assumes the rotating-wave approximation. raises (lowers) the state of the two-level atom whose transition frequency is . We note that the waveguide model in can describe any single-mode dielectric or plasmonic waveguide in the optical regime [13] or a transmission line in the microwave regime [26]. Since the atomic linewidth is typically quite narrow, a linear approximation of the waveguide dispersion relation, which we adopt here, is generally sufficient to describe the waveguide. We further note that the two-level atom in can describe either a real atom or an artificial atom, such as a quantum dot [13] or superconducting qubit [15].
To model dissipation we introduce reservoir–atom coupling as described by and consider the composite system as shown schematically in Fig. 1, where Here and are the reservoir photon operators and is the atom–reservoir coupling rate. We assume the reservoir has a continuum of modes centered around the atomic transition frequency . We further assume that within the spectral range of interest the dispersion of these reservoir modes is well described by a group velocity . The reservoir will initially be set to the vacuum state. As a result of interaction with the atom, photons originating in the waveguide may leak back into the waveguide, described in Eq. (1), or they may leak into the reservoir as described in Eq. (2), resulting in loss.
Figure 1.Schematic of sample system considered. A waveguide side coupled with rate to a two-level atom with transition frequency . The atom is additionally coupled with rate to a reservoir. The initial state in the waveguide is either a one- or a two-photon Fock state, whereas the auxiliary waveguide is initially in the vacuum state.
To proceed, we exploit the spatial inversion symmetry of and decompose it to even and odd subspaces , where . We transform from right- and left-moving waveguide photon operators to even and odd photon operators. We omit from here on the subscripts of the even mode operators: with the above operators, the even subspace Hamiltonian and the odd subspace Hamiltonian follow
Here describes a one-way propagating chiral waveguide mode that interacts with the atom. describes a one-way propagating waveguide mode without any interaction with the atom. Once the solutions for are known, the solution to can be straightforwardly obtained by using Eqs. (3a) and (3b). In what follows we consider only the even Hamiltonian, since it contains all of the nontrivial physics of the system.
The Hamiltonian can be recast in k-space by defining momentum-space photon operators and , We note that the total number of excitations commutes with , , allowing for the solution to proceed independently in each excitation number manifold.
3. INPUT–OUTPUT FORMALISM
In the absence of dissipation, the one- and two-photon scattering matrices of this system were solved in [10], to which we refer the reader for a detailed derivation of the formalism we use. Our focus here is on the addition of dissipation to the formalism.
By defining input and output field operators for the waveguide, and reservoir, we obtain (in the limits and ) the input–output equations, We see that the presence of dissipation in Eq. (6) in the form of damping rate is accompanied by the input field operator from the reservoir. Below we use Eqs. (4)–(6) to obtain one and two-photon S-matrices.
A. One-Photon S-Matrix
The one-photon S-matrix relates an incoming free-photon state with momentum to an outgoing free-photon state with momentum and is given by , where is the scattering operator and is an incoming one-excitation interacting eigenstate that evolved from a free-photon state in the distant past. is an outgoing one-excitation interacting eigenstate that will evolve into a free photon state in the distant future. The one-photon S-matrix can be reexpressed as follows: where we have used the input–output relation in Eq. (4) and [10]. To obtain the matrix element , we solve for its time-domain counterpart using Eq. (6), noting that and , where is an incoming reservoir eigenstate. The resulting one-photon transport matrix elements are
We note that, as a result of reservoir damping, . This seemingly innocuous difference between the dissipative case and the dissipation-free case stems from the fact that the set is not complete in the presence of dissipation. We will see the crucial importance of this fact in the two-photon calculation.
B. Two-Photon S-Matrix
The transport of two waveguide photons is given by the two-photon S-matrix, which describes the probability amplitude of incoming photons with momenta and to scatter into two outgoing photons with momenta and and can be expressed as To proceed, we insert the resolution of the identity in the single-excitation subspace which comprises of a sum of both waveguide and reservoir interacting eigenstates outer products where we have used the definitions of the one-photon scattering amplitudes and . Inserting the Fourier transform of Eq. (4) and noting that , we have We now need to solve for the matrix elements and . To do so we sandwich Eq. (6) between the appropriate states, getting the equations where we note that and . To proceed, we solve for the source terms in Eqs. (8) and (9) by using and inserting a resolution of the identity in the zero-excitation Hilbert space, Having calculated the source terms, we may now solve Eqs. (8) and (9) and Fourier transform their solution to obtain Finally, we insert Eqs. (10) and (11) into Eq. (7) to get the two-photon S-matrix, We note that the two-photon S-matrix is invariant under exchange of incoming and outgoing momenta as required by time-reversal symmetry, and also with respect to exchange of incoming or outgoing momenta and , as required by the Bose statistics of photons.
Instead of the correct approach outlined above, one might instead consider an ad hoc approach where an imaginary part is added to the frequency of the atom, resulting in the following Hamiltonian: From this ad hoc Hamiltonian, one then follows the steps in [10]. With this approach, the Hilbert space consists of only waveguide and atom states—with no reservoir. For one photon, this approach yields the correct S-matrix. For two photons, the term of Eq. (10) would have been missing, and the resulting incorrect two-photon S-matrix would have been This form of the S-matrix violates the particle exchange symmetries mentioned above.
If we examine the form of the S-matrix in Eq. (13), it might seem odd that, even though the result is expressible by single-photon scattering amplitudes, it violates exchange symmetry of outgoing particle momenta while the one-photon S-matrix does not. This is explained as follows. The physical origin of the symmetry violation is in the scattering amplitude . This amplitude is proportional to the Fourier coefficient of the atom’s spontaneous emission amplitude into the waveguide. However, in the presence of dissipation, the atom may also spontaneously emit a photon into the reservoir, corresponding to the amplitude . This amplitude is provided by the term in Eq. (10) and originates in our inclusion of the reservoir when inserting the two-excitation resolution of the identity .
4. A GENERAL RULE FOR DISSIPATION
When examining the correct form of the two-photon S-matrix in Eq. (12), we see that the presence of dissipation is seen as the addition of a term in the one-excitation transmission and atomic excitation amplitudes. It follows that the one- and two-photon S-matrices of a two-level atom with reservoir coupling rate , as described by , can be obtained from the dissipation-free S-matrices by making the replacement . However, we stress that this substitution should be made in the S-matrix and not in the Hamiltonian, since, as described above, doing so would leads to an incorrect result in the two-photon case.
This result is even more general. In the presence of reservoir dissipation channels with coupling rates , the resulting S-matrix is obtained from the dissipation-free S-matrix by making the substitution . Moreover, this result holds for a general quantum multilevel system as long as the Hamiltonian light-matter interaction term can be described in the rotating-wave approximation. Furthermore, the quantum multilevel system can be either in a cavity or directly coupled to a waveguide.
Finally, the eigenstate of the ad hoc Hamiltonian turns out to be the same as the true scattering eigenstate—other than the reservoir degree of freedom, which it is missing. This fact allows a correct real-space interacting eigenstate solution [6,7,27,28] for the nonreservoir dynamics that account for the loss of energy that accompanies reservoir coupling, in complete agreement with the solutions presented here. In particular, the existence of two-photon bound states, which is one of the important predictions made by using the wave-function method, is still valid in the presence of dissipation. In conclusion, the method described in this paper provides a straightforward and correct construction of the S-matrix in the presence of dissipation while working with the full Hermitian Hamiltonian of the system in the context of the input–output formalism.
ACKNOWLEDGMENTS
Acknowledgment. The work is supported by an AFOSR-MURI program on quantum metamaterial (grant FA9550-12-1-0488).