Perfectly capturing traveling single photons of arbitrary
temporal wavepackets with a single tunable device

arXiv:1609.05643v1 [quant-ph] 19 Sep 2016

Hendra I. Nurdin∗, Matthew R. James†, and Naoki Yamamoto‡

Abstract
We derive the explicit analytical form of the time-dependent coupling parameter to an external field for perfect absorption of traveling single photon fields with
arbitrary temporal profiles by a tunable single input-output open quantum system,
which can be realized as either a single qubit or single resonator system. However, the
time-dependent coupling parameter for perfect absorption has a singularity at t = 0
and constraints on real systems prohibit a faithful physical realization of the perfect
absorber. A numerical example is included to illustrate the absorber’s performance
under practical limitations on the coupling strength.

1

Introduction

Quantum networks composed of open quantum systems as network nodes that are interconnected by traveling quantum optical fields are of much interest for quantum information
applications, including quantum communication and quantum computing; see, e.g., [1, 2]
and the references therein. Photons in optical fields are ideal particles for transmitting information between nodes in a quantum network as they can propagate through free space
and various physical media. Quantum systems at the nodes are typically matter systems
such as an atom, atomic ensembles, superconducting circuits, amongst many possibilities.
Each node can process or store classical or quantum information and exchange its information content with an optical field, either to receive information contained in the field or to
transmit information into the field. Thus an important problem for quantum networks is
the problem of state transfer between an optical field impinging on a node and the matter
system at the node. That is, the transfer of the quantum state of the optical field to the
matter and vice-versa.
∗

H. I. Nurdin is with the School of Electrical Engineering and Telecommunincations, UNSW Australia,
Sydney NSW 2052, Australia (email: h.nurdin@unsw.edu.au). Research was supported by the Australian
Research Council.
†
M. R. James is with the School of Engineering, The Australian National University, Canberra ACT
0200, Australia (email: matthew.james@anu.edu.au). Research was supported by the Australian Research
Council.
‡
N. Yamamoto is with the Department of Applied Physics and Physico-Informatics, Keio University,
Hiyoshi 3-14-1, Kohoku, Yokohama, Japan (email: yamamoto@appi.keio.ac.jp). Research was supported
by JSPS Grant-in-Aid No. 15K06151.

1

The problem of state transfer in a quantum network with two nodes was first considered
in the seminal work [3]. In this work, the authors consider the perfect transfer of an
arbitrary superposition state α|0i + β|1i from a qubit on one node to a qubit on the other
via a one way optical field connecting the two nodes. Each node is a cavity QED system
that can be modelled as a qubit coupled to an optical cavity which is in turn coupled to the
optical field. It was shown that perfect transfer can be achieved by suitably modulating the
Rabi frequency and phase of a Raman laser driving the qubits at each node. They employ
the quantum filtering equation or quantum trajectory equation [4, 5] for the network and
apply the dark state principle to derive differential equations characterising the required
modulation. The dark state principle states that there should be no photon counted at the
output optical field reflected from the receiving node when a photon counting measurement
is performed on that output field.
In this paper, we are interested in the problem of the perfect transfer of a traveling
single photon field into a system that can function as a perfect single photon absorber.
The state of a travelling single photon field is characterized
by a temporal wavepacket
R∞
ξ, a complex-valued function of time, that satisfies 0 |ξ(s)|2 ds = 1. The wavepacket
gives the detection probability of the single photon state, that is, the probability that the
photon will be detected (say, byR registering a click in a photo detector) in a time interval
t
[t1 , t2 ] (0 ≤ t1 < t2 ) is given by t12 |ξ(s)|2 ds. We are interested in a tunable single photon
absorber with a coupling parameter that can be modulated to absorb single photons of any
wavepacket shape.
In [6] it was shown how single photon fields with symmetric temporal profiles can be
mode-matched to be absorbed by a coupled cavity-oscillator system into a single photon
state of the oscillator, by suitably modulating the coupling parameter between the cavity
and oscillator modes. In the work [7], it was shown that an absorber can be implemented
with a cavity QED system composed of a three level atom coupled to an optical cavity.
Perfect absorption of a photon with an arbitrary temporal shape could be achieved by appropriately modulating the amplitude of a laser beam driving a Raman transition process
in the atom, under the assumption of no spontaneous emission from the atom. The modulation of the laser beam is tailored to the temporal profile of the incoming single photon field.
In this paper, we consider a different class of single photon absorbers, a single two-level
system (qubit) or, equivalently, a single resonator with a tunable coupling parameter to an
external optical field. This type of coupling has been theoretically proposed for the transfer
of state between two microwave resonators connected by a one-way travelling optical field
between the resonators [8]. Importantly, such a tunable coupling parameter has already
been demonstrated experimentally in a microwave superconducting resonator [9], where
the tunability is realized using an externally controlled variable inductance. This coupling
is analogous to a mirror with tunable transmissivity on an optical cavity. The contribution of this work is to analytically derive the exact form for the time-dependent coupling
parameter for perfect absorption of a single photon with an arbitrary wavepacket shape.
We describe the system with a QSDE and derive the optimal evolving coupling parameter
by two methods: by explicitly solving the QSDE, and by application of the zero-dynamics
principle from [10]. We note that a form of this principle had also been employed in earlier
works [6, 8, 7] without being referred to as such. That is, these works require that the
incoming field and the field reflected from the system interfere destructively, resulting in
no photon in the output field (zero output dynamics).

2

This paper is organized as follows. Section 2 introduces the notation of the paper
and gives a brief overview of quantum stochastic calculus, quantum stochastic differential
equations, single photon generators, and systems driven by a single photon field. Section
3 introduces the model of the single photon absorber of interest, analytically derives the
modulating function for the coupling parameter for perfect absorption of any single photon
field, and develops a numerical example showing the effects of practical limitations on the
absorber’s performance. Finally, Section 4 draws a conclusion for the paper.

2

Preliminaries

√
Notation. We will use ı = −1, ∗ to denote the adjoint of a linear operator as well as
the conjugate of a complex number. If A = [ajk ] then A# = [a∗jk ], and A† = (A# )> , where
(·)> denotes matrix transposition. <{A} = (A + A# )/2 and ={A} = 2ı1 (A − A# ). We
use R+ to denote the set of non-negative real numbers and L2 (R+ ; C) to denote the set of
square-integrable functions on R+ . hXi denotes the quantum expectation of an operator
X, and Tr(X) denotes the trace of X.

2.1

Quantum stochastic calculus and quantum stochastic differential equations

We will be working with open Markov quantum systems that are coupled to n continuousmode boson fields indexed by j = 1, 2, . . . , m with annihilation field operators ηj (t) satisfying the field commutation relations [ηj (t), ηk (t0 )∗ ] = δjk δ(t − t0 ) and [ηj (t), ηk (t0 )] = 0. For
our purposes, we can focus on fields
R t in a vacuum state. Let us introduce the integrated
field annihilation process Aj (t) = 0 ηj (s)ds and its adjoint process, the integrated field
Rt
creation process, A∗j (t) = 0 ηj∗ (s)ds. In the vacuum representation, their future-pointing
Itō increments dAj (t) = Aj (t + dt) − Aj (t) and dA∗j (t) = A∗j (t + dt) − A∗j (t) satisfy the
quantum Itō table [11, 12, 13]
×
dA∗k dAk
dAj δjk dt 0
dA∗j 0
0
We may also define the counting process (or gauge process)
Z t
Λjk (t) =
b∗j (r)bk (r)dr,
0

which may be included in the Itō table [11]. The additional non-trivial products of differentials are
dΛjk dA∗l = δkl dA∗j , dAj dΛkl = δjk dAl , dΛjk AΛli = δkl dΛji .
Using the processes A = (A1 , A2 , . . . , Am )> , A# = (A∗1 , A∗2 , . . . , A∗m )> and Λ = [Λjk ]j,k=1,...,m ,
one may define quantum stochastic integrals of adapted processes on the tensor product
of the system and joint boson (symmetric) Fock space of the fields. The system is the
quantum mechanical object that is being coupled to the fields, and adapted means that
at time t the process acts trivially on the portion of the boson Fock space after time t ,
see, e.g., [11, 12, 13] for details. An adapted process commutes at time t with all of the
3

future pointing differentials. The product of two adapted processes X(t) and Y (t) is again
adapted, and the increment of the product obeys the quantum Itō rule
d(X(t)Y (t)) = (dX(t))Y (t) + X(t)dY (t) + dX(t)dY (t),
Based on these quantum stochastic integrals, one may define quantum stochastic differential equations (QSDEs).
For the remainder of the paper, we will only consider a system coupled to a single
field, m = 1 and A(t) = A1 (t), A∗ (t) = A∗1 (t), Λ(t) = Λ11 (t). An important QSDE that
describes the joint unitary evolution of an open Markov system coupled to vacuum boson
fields, common in quantum optics and related fields, is the Hudson-Parthasarathy QSDE
given by
dU (t) = (−(ıH(t) + 1/2L(t)∗ L(t))dt + dA(t)∗ L(t) − L(t)∗ S(t)dA(t) + (S(t) − I)dΛ(t))U (t),
(1)
with initial condition U (0) = I. We adopt a very general setting [14] where H(t) = H(t)∗
can be a general adapted process representing the time-dependent Hamiltonian of the
system, L(t) an adapted process representing the (possibly time-dependent) coupling of
the system to the creation process A∗ (t), and S(t) a unitary adapted process, S(t)∗ S(t) =
S(t)S(t)∗ = I, representing the (possibly time-dependent) coupling of the system to the
process Λ of the field. This QSDE has a unique solution whenever S(t), L(t), H(t) are
bounded operators for all t ≥ 0 [14]. Moreover, in that case the solution is guaranteed to
be unitary. The input field A(t) after the interaction is transformed into the output field
Y (t) = U (t)∗ A(t)U (t). Let jt (X) = U (t)∗ XU (t) denote the Heisenberg picture evolution
of a system operator X, then we have the following QSDEs for the Heisenberg picture
evolution of jt (X) and Y (t) [15],
djt (X) = jt (Lt (X)) + dA(t)∗ jt (S(t)∗ [X, L(t)]) + jt ([L(t)∗ , X]S(t))dA(t)
+jt (S(t)∗ XS(t) − X)dΛ(t),
dY (t) = jt (L(t)) + jt (S(t))dA(t),
with
Lt (X) = 1/2L(t)∗ [X, L(t)] + 1/2[X, L(t)]L(t) − ı[X, H(t)].
We will denote an input-output open quantum system G with parameters S(t), L(t), H(t)
with the notation G = (S(t), L(t), H(t)) [15]. The output of system G1 = (S1 (t), L1 (t), H1 (t))
can be passed as the input to a second system G2 = (S2 (t), L2 (t), H2 (t)) forming a cascaded
system which is again an input-output open quantum system G3 , the parameters of which
is given by the series product rule /,
G3 = G2 / G1 = (S2 (t)S1 (t), L2 (t) + S2 (t)L1 (t), H1 (t) + H2 (t) + ={L2 (t)∗ S2 (t)L1 (t)}).
We note that G1 and G2 can be subsystems of one system sharing the same Hilbert space,
see [15] for details. The series product is associative, thus if Gn = (Sn (t), Ln (t), Hn (t))
then the series product Gn / Gn−1 / . . . / G1 is unambiguously defined.

4

2.2

Single photon generator

Our main interest is in systems driven by a traveling single photon field with wavepacket
ξ, defined by a continuous-mode single photon state |1ξ i on the boson Fock space, given by
Rt
Rt
|1ξ i = 0 ξ(s)η(s)∗ ds |Ωi = 0 ξ(s)dA∗ (s)|Ωi [16, 17, 18, 19], where |Ωi denotes the vacuum
state of the field.
In [20, 19] it has been shown that |1ξ i can be generated as the output of a two-level
system (qubit) that is coupled to single vacuum input field through a modulated coupling
coefficient. The generator G1 has the Hilbert space H1 = C2 , where we take as basis
vectors |0i1 = (0, 1)> and |1i1 = (1, 0)> , and is given by G1 = (I, λ(t)σ1,− , 0), where σ1,−
is the qubit’s lowering operator,


0 0
σ1,− = |0i1 h1| =
1 0
and λ(t) is a time-dependent complex coupling coefficient of the system to the field given
by
ξ(t)
λ(t) = p
,
w(t)

(2)

R∞
with w(t) = t |ξ(s)|2 ds. The coupling coefficient depends on the wavepacket ξ of the
single photon field that one wishes to generate. Let U1 (t) denote the unitary solution to
the generator’s QSDE. When the generator is initialized in the state |1i1 then we have
[20, 19],
Z t
p
ξ(s)dA∗ (s)|Ωi
U1 (t)|1i1 |Ωi = w(t)|1i1 ⊗ |Ωi + |0i1 ⊗
0

and so limt→∞ U1 (t)|1i|Ωi = |0i ⊗ |1ξ i. Thus, if a photon detection measurement is performed at the output of the generatorR at any time t ≥ 0, the probability of detecting a
t
photon in the interval [0, t] is exactly 0 |ξ(s)|2 ds, as expected.

2.3

Systems driven by a single photon

Systems driven by a traveling single photon field with wavepacket ξ can be viewed equivalently as being driven by an ancillary generator that outputs the single photon field, as
described in the preceding subsection. Suppose that a system G2 = (S2 (t), L2 (t), H2 (t))
on the Hilbert space H2 is driven by a single photon field with wavepacket ξ. Then the
dynamics of the system can be analyzed by studying the equivalent cascaded system
G2 / G1 = (S2 (t), L2 (t) + λ(t)S2 (t)∗ σ1,− , H(t) + ={λ(t)L2 (t)∗ S2 (t)σ1,− })
that is driven by a vacuum field. The cascaded system is defined on the composite Hilbert
space C2 ⊗ H2 . Note that for notational simplicity for an arbitrary operator X1 on the
generator and X2 on the system, we will often denote the ampliations X1 ⊗ I, I ⊗ X2 , and
X1 ⊗ X2 simply as X1 , X2 , and X1 X2 , etc, with the tensor product being omitted.

5

3

A tunable single photon absorber with time-dependent
coupling to an external field

We now describe our tunable photon absorber that can perfectly absorb a single photon
field with an arbitrary wavepacket. The Hilbert space of the absorber is that of a qubit,
but since only a single excitation will be involved (i.e., a single photon) it can actually be
implemented by a quantum harmonic oscillator with a tunable coupling parameter such as
the microwave system realized in [7]. We will say more about the physical realization of
the absorber near the end of this section.
We choose |0i2 = (0, 1)> and |1i2 = (1, 0)> as basis vectors for the absorber Hilbert
space H2 = C2 . Our perfect single photon absorber is then a system G2 = (I, γ(t)σ2,− , 0),
with σ2,− the qubit lowering operator on H2 . Note that the absorber Hamiltonian has been
set to 0, so we are implicitly working in a rotating frame with the respect to the absorber’s
original Hamiltonian. We only require that the original Hamiltonian commutes with σ2,−
or only rotates the latter in time as e−ıωt σ2,− for some real constant ω ( e−ıωt can then
be absorbed into γ(t)), but is otherwise arbitrary. In principle, the analysis to follow can
be adapted to accommodate more general Hamiltonians but we do not pursue this here.
The time-dependent parameter γ will be chosen to enable perfect absorption of any single
photon field, with γ dependent on the wavepacket ξ of the latter. Let G1 be the ancillary
generator for a single photon as given in Section 2.2 for a given ξ. The cascade of G1 onto
G2 is given by:
G2 / G1 = (I, λ(t)σ1,− + γ(t)σ2,− , ={γ(t)∗ λ(t)σ2,+ σ1,− }),

(3)

∗
where σj,+ = σj,−
denotes the raising operator for the system Gj . The generator G1
is initialized in the excited state |1i1 while the absorber G2 is initialized in the ground
state |0i2 . The field to which the cascade is coupled to is in the vacuum state |Ωi. We
shall show that for any wavepacket ξ(t), γ(t) can always be chosen such that in the limit
t → ∞, the joint system and field state converges to |0i1 |1i2 |Ωi. That is, a single photon
is asymptotically perfectly absorbed by G2 . We shall derive the form for γ explicitly in
two ways: (i) by explictly solving the QSDE for the cascaded system in Section 3.1, and
(ii) by application of the zero-dynamics principle in Section 3.2.
It will turn out that the coupling parameter γ will have a singularity at t = 0 for
any wavepacket ξ. In actual practical implementation this singularity must be truncated
to some finite value. We propose a truncation in Section 3.3 and apply it in a numerical
example for an exponentially decaying wavepacket to give an illustration of the effect of this
truncation on the absorber’s performance in capturing a single photon with this wavepacket
shape.
Coming back to the discussion at the beginning of this section, we note that since the
absorber can have no more than a single quanta of excitation when driven by a single
photon field, the dynamics of G2 when initialized at |0i2 can be embedded in the dynamics
of an equivalent quantum harmonic oscillator system G02 = (I, λ(t)a, 0), where a is the
annihilation operator of the oscillator, satisfying the commutation relation [a, a∗ ] = 1.
This is because H2 is isomorphic to span{|0io , |1io }, where |0io and |1io is the vacuum
and 1-photon Fock state of the oscillator, respectively, and σ2,− can be identified with the
restriction of a to span{|0io , |1io }. Indeed, this kind of embedding of the dynamics of a

6

qubit system into a quantum harmonic oscillator under single photon driving has already
been exploited in [21].

3.1

Determining γ by solving the QSDE for G2 / G1

The wavefunction of the cascaded system (generator and absorber) and the external field is
initially |1i|0i|Ωi and at time t ≥ 0 is of the form ψ(t) = (ψ1 (t), ψ2 (t), ψ3 (t), ψ4 (t))> ⊗|Ω[t i,
where ψ(t) is adapted on the joint cascaded system and Fock space, in the sense that ψ(t)
lives in the tensor product of the system and the portion of the Fock space up to time
t. That is, U (t)|1i|0i|Ωi = ψ(t) ⊗ |Ω[t i, where U (t) is a solution of the QSDE (1) with
coefficients given by (3). Here |Ωt] i and |Ω[t i denote the portion of Fock vacuum from time
0 up to time t, and from time t onwards, respectively, see [11, 12]. We have that the vector
ψ(t) satisfies the QSDE:


−(|λ(t)|2 + |γ(t)|2 )/2
0
0
0

0
−|λ(t)|2 /2
0
0 
 ψ(t) ⊗ |Ω[t idt
d(ψ(t) ⊗ |Ω[t i) = 
∗
2

0
−γ(t) λ(t) −|γ(t)| /2 0 
0
0
0
0


0
0
0 0
 γ(t) 0
0 0 
 ψ(t) ⊗ dA(t)∗ |Ω[t i,
+
 λ(t) 0
0 0 
0 λ(t) γ(t) 0
with initial condition ψ(0) = |1i1 |0i2 = (0, 1, 0, 0)> . From the above equation one readily
gets that


Z t
2
2
1
ψ1 (t) = exp − /2 (|λ(s)| + |γ(s)| )ds ψ1 (0)|Ωt] i
0

= 0,
for all t ≥ 0. Using this fact and recalling from [19] (using (2)) that

 sR ∞
Z t
|ξ(τ )|2 dτ
,
exp −1/2
|λ(τ )|2 dτ = Rt∞
|ξ(τ )|2 dτ
s
s
for 0 ≤ s ≤ t, we can also solve for ψ2 (t), obtaining


Z t
2
1
|λ(τ )| dτ ψ2 (0)|Ωt] i
ψ2 (t) = exp − /2
0
sZ
∞

|ξ(τ )|2 dτ |Ωt] i.

=
t

We want that as t → ∞, ψ(t) → |0i1 |1i2 |Ωi = (0, 0, |Ωi, 0)> . That is, G2 completely
absorbs the single photon emitted by G1 , so that the output field state goes to the vacuum.
From the equation for ψ4 (t) we have that
Z t
(λ(s)ψ2 (s)dA(s)∗ + γ(t)ψ3 (s)dA(s)∗ )|Ωt] i,
ψ4 (t) = ψ4 (0)|Ωt] i +
0
Z t
=
(ξ(s)dA(s)∗ + γ(s)ψ3 (s)dA(s)∗ )|Ωt] i,
0

7

which we want to go to 0 as t → ∞. Also, since
dψ3 (t) = −1/2|γ(t)|2 ψ3 (t)dt − γ(t)∗ λ(t)ψ2 (t)dt,
we can solve this to give:


Z t
2
1
ψ3 (t) = exp − /2
|γ(s)| ds ψ3 (0)|Ωt] i
0


Z t
Z t
2
exp −1/2
|γ(τ )| dτ γ(s)∗ λ(s)ψ2 (s)ds)|Ωt] i,
−
0
s


Z t
Z t
2
exp −1/2
|γ(τ )| dτ γ(s)∗ λ(s)ψ2 (s)ds)|Ωt] i,
= −
0

s

the second line following from the fact that ψ3 (0) = 0. The expressions for ψ4 (t) and ψ3 (t)
just given, and the form of λ(t) for the generator, suggest the following as an educated
guess for γ(t) to achieve limt→∞ ψ4 (t) = 0 and limt→∞ ψ3 (t) = |Ωi,
γ(t) = −eıφ0 qR

t
0

ξ(t)

,

|ξ(s)|2 ds

for some arbitrary real constant φ0 . We will verify below that this form of the modulation
function γ(t) indeed achieves perfect absorption. Note that although γ is singular at t = 0,
γ(0) = −∞, it is square integrable on [s, t] for any 0 < s ≤ t ≤ ∞. Indeed, by direct
integration,

Z t
Z τ
Z t
2
2
2
|ξ(y)| dy dτ,
|γ(τ )| dτ =
|ξ(τ )| /
s
0
s

Z t  Z τ
2
=
d ln
|ξ(y)| dy ,
s
0
!
Rt
2
|ξ(y)|
dy
= ln R 0s
.
|ξ(y)|2 dy
0
Moreover, from this it follows immediately that
 rR
s

|ξ(y)|2 dy



Z t
R0t
s>0

2
0 |ξ(y)| dy
exp −1/2
|γ(τ )|2 dτ
=
,
0
s = 0, t > 0

s


1
s=t=0
which is thus well-defined for any 0 ≤ s ≤ t ≤ ∞. Substituting the expression for γ(t) and
solving for ψ3 (t) gives
(
Rt
e−ıφ0 √R t 1 2
|ξ(s)|2 ds|Ωt] i t > 0
0
|ξ(τ
)|
dτ
ψ3 (t) =
.
0
0
t=0

8

Notice that limt→∞ ψ3 (t) = |Ωi, as desired. Substituting all of the above into the expression
for ψ4 (t) gives,
Z t
(λ(s)ψ2 (s)dA(s)∗ + γ(s)ψ3 (s)dA(s)∗ )|Ωt] i,
ψ4 (t) =
Z0 t
=
(ξ(s)dA(s)∗ − ξ(s)dA(s)∗ )|Ωt] i,
0

= 0,
again exactly as expected. Thus, weq
conclude that a two-level system of the form G =
Rt
(I, γ(t)σ− , 0) with γ(t) = −eıφ0 ξ(t)/ 0 |ξ(τ )|2 dτ can perfectly absorb a single photon
with an arbitrary wavepacket ξ(t).

3.2

Determining γ by application of the zero-dynamics principle

We now show that the modulation signal γ(t) can also be found by application of the
zero-dynamics principle proposed in [10]. Recall that the cascaded generator and absorber
has the coupling operator L(t) = λ(t)σ1,− + γ(t)σ2,− and the Hamiltonian is H(t) =
={γ(t)∗ λ(t)σ2,+ σ1,− }. The initial state of the composite generator and absorber system is
|1i1 |0i2 . The formulation of the zero-dynamics
R t principle given in [10] is in the Heisenberg
picture. It states that the output field Y (t) = 0 js (L(s))ds+A(t) should be a vacuum field.
This is necessary as no photon should be present in the output field if perfect absorption
is taking place. This means that one must have [22],

 Z ∞



Z ∞
∗
∗
2
1
exp ı
(u(t) dY (t) + u(t)dY (t) )
= exp − /2
|u(s)| ds ,
0

0

for any u ∈ L2 (R+ ; C), where the right hand side is the characteristic function of the
vacuum state of a bosonic field. Evaluating the left hand side is in general difficult, if
not intractable, except in special instances like the linear quantum networks analysed
in [10]. However, the zero-dynamics principle can still be applied, but by looking at it
from another equivalent viewpoint. Zero-dynamics corresponds to the state of the output
field always being vacuum at all times t ≥ 0, U (t)|1i1 |0i2 |Ωi = |ϕ(t)i|Ωi for some pure
state vector |ϕ(t)i on H1 ⊗ H2 for all t ≥ 0. This is similar in spirit to the form of
the principle employed in [8, 7]. Thus we may state the ansatz that the wavefunction is
|ϕ(t)i = α(t)|1i1 |0i2 + β(t)|0i1 |1i2 for t ≥ 0 for some scalar complex functions α(t) and
β(t) that will be sought, satisfying |α(t)|2 + |β(t)|2 = 1. We also have the initial condition
α0 = 1 and β0 = 0. Substituting this ansatz into the joint QSDE of the system and field
gives:
d|ϕ(t)i ⊗ |Ωi = (−(ıH(t) + 1/2L(t)∗ L(t))dt + L(t)dA(t)∗ − L(t)∗ dA(t)) |ϕ(t)i ⊗ |Ωi.
Since the input to the generator-absorber is vacuum, for zero-dynamics to hold the term
L(t)dA(t)∗ in the QSDE should be 0 for all times (otherwise dA∗ (t) will create a photon
in the field after time t and the field will no longer be in the vacuum state), meaning that
L(t)|ϕ(t)i = 0 for all t ≥ 0. This gives the condition
(λ(t)α(t) + β(t)γ(t))|0i1 |0i2 = 0, ∀t ≥ 0
9

and thus,
λ(t)α(t) + β(t)γ(t) = 0, ∀t ≥ 0.

(4)

Moreover, noting that L(t)∗ dA(t)|ϕ(t)i|Ωi = 0 (since dA(t) annihilates the portion of the
vacuum after time t), the QSDE reduces to the deterministic equations,
d|ϕ(t)i = −ıH(t)|ϕ(t)idt,
α̇(t)|1i1 |0i2 + β̇(t)|0i1 |1i2 = 1/2γ(t)λ(t)∗ β(t)|1i1 |0i2 − 1/2γ(t)∗ λ(t)α(t)|0i1 |1i2
Thus we arrive at the following set of differential equations,
α̇(t) = γ(t)λ(t)∗ β(t),
β̇(t) = −γ(t)∗ λ(t)α(t),
λ(t)α(t) + β(t)γ(t) = 0
Using the constraint (4), subtituting β(t)γ(t) = −λ(t)α(t) gives,
α̇(t) = −1/2|λ(t)|2 α(t),
β̇(t) = 1/2|γ(t)|2 β(t),
λ(t)α(t) + β(t)γ(t) = 0
The solution for α with initial condition α0 = 1 is:
sZ
α(t) = e−

Rt
0

|λ(s)|2 ds

∞

|ξ(s)|2 ds,

α0 =
t

Rt
But we also have that |α(t)|2 + |β(t)|2 = 1, so |β(t)|2 = 0 |ξ(s)|2 ds. It follows from the
constraint (4) that
|ξ(t)|
.
|γ(t)| = |λ(t)α(t)/β(t)| = qR
t
2 ds
|ξ(s)|
0
So, we arrive at
sZ
β(t) =
0

t

|ξ(s)|2 ds e−iφ(t) , γ(t) = − qR

ξ(t)
t
0

eiφ(t) .

|ξ(s)|2 ds

for some arbitrary real-valued differentiable function φ(t). However, to satisfy the differential equation for β(t), φ(t) must in fact be a real constant, say, φ0 . Whence,
sZ
t
ξ(t)
β(t) =
|ξ(s)|2 ds e−iφ0 , γ(t) = − qR
eiφ0 .
t
0
|ξ(s)|2 ds
0
Thus we recover the form of the modulating function γ(t) that was obtained in Section
3.1 by explicitly solving the QSDE, indicating the power of the zero-dynamics principle
for state-transfer problems such as this. Since one does not need to solve a QSDE, this
approach would be more widely applicable. Moreover, even if the QSDE is explicitly
solvable, it allows the differential equations characterizing γ(t) to be derived thus avoiding
having to make an educated guess about γ(t).
10

3.3

Effect of imperfect realization of γ

We have seen that the analytical form of γ for perfect absorption has a singularity that
grows as O(t−1/2 ) as t → 0 for any wavepacket ξ. This poses a practical challenge as
an infinite coupling magnitude cannot be physically realized. Moreover, in the potential
implementation of the absorber using a single mode resonator, the coupling parameter
should be much smaller than the resonator’s free spectral range as to not excite higher
frequency modes of the resonator. A sub-optimal implementation would be to truncate
the magnitude of γ(t) for small values of t to a finite value. However, this will mean that
the absorber will no longer perfectly absorb a single photon. In this section, we numerically
evaluate the effect of bounding the magnitude of γ(t) for small t on a particular example.
Let nj = σj,+ σj,− denote the number operator for the generator (j = 1) and absorber
(j = 2). To ease the notation, in the following we will often not explicitly write the
time dependence on operators, with the time dependence being implicitly understood. We
will be interested in the evolution of the mean number operator for the absorber, hn2 i,
in the Heisenberg picture, as this gives the probability of excitation of a single photon
in the absorber. In the perfect absorption case we have already discussed, we have that
limt→∞ hn2 i = 1. Since the combined generator-absorber is driven by a vacuum field, the
evolution of the mean hX1 X2 i is given by the ordinary differential equation (the general
expression for Lt (X1 X2 ) is derived in the Appendix),
d
hX1 X2 i = hLt (X1 X2 )i,
dt
= γ(t)∗ λ(t) X1 σ1,− [σ2,+ , X2 ]
+1/2|λ(t)|2 (σ1,+ [X1 , σ1,− ] + [σ1,+ , X1 ]σ1,− )X2
+λ(t)∗ γ(t) σ1,+ X1 [X2 , σ2,− ]
+1/2|γ(t)|2 (σ2,+ [X2 , σ2,− ] + [σ2,+ , X2 ]σ2,− )X1 .
When X2 = n2 and X1 = 1, we get that:
d
hn2 i = hLt (n2 )i,
dt
= γ(t)∗ λ(t) σ1,− [σ2,+ , n2 ] + λ(t)∗ γ(t) σ1,+ [n2 , σ2,− ]
+1/2|γ(t)|2 (σ2,+ [n2 , σ2,− ] + [σ2,+ , n2 ]σ2,− ) ,
= −|γ(t)|2 hn2 i − γ(t)∗ λ(t) σ1,+ σ2,−

∗

− γ(t)λ(t)∗ σ1,+ σ2,− .

(5)

(6)

Let σ2,z = [σ2,+ , σ2,− ]. Then when X2 = σ2,− and X1 = σ1,+ the equation is:
d
hσ1,+ σ2,− i = hL(σ1,+ σ2,− )i,
dt
= −1/2(|λ(t)|2 + |γ(t)|2 )hσ1,+ σ2,− i + γ(t)∗ λ(t) n1 σ2,z i.

(7)

When X1 = n1 and X2 = σ2,z , the equation is:
d
hn1 σ2,z i = hLt (n1 σ2,z )i,
dt
= −|λ(t)|2 hn1 σ2,z i + 2|γ(t)|2 n1 n2 i.
11

(8)

Finally, when X1 = n1 and X2 = n2 the equation is:
d
hn1 n2 i = hLt (n1 n2 )i,
dt
= −(|λ(t)|2 + |γ(t)|2 )hn1 n2 i.

(9)

The initial conditions are hn2 (0)i = 0, hσ1,+ (0)σ2,− (0)i = 0, hn1 (0)σ2,z (0)i = −1, and
hn1 (0)n2 (0)i = 0. From the last initial condition, we can solve for (9) to get that
hn1 (t)n2 (t)i = 0 for all t ≥ 0. Thus there are actually only three coupled equations
that can have a non-trivial solution. The remaining initial conditions can then be used to
sequentially solve equations (8), (7), and (6) (in that order), yielding the explicit solutions.
However, here we are interested in evaluating the effect of approximating the function γ(t),
which is singular at t = 0, with a function γT (t) given by

 − √R Tξ(T ) 2 , t ≤ T,
0 |ξ(s)| ds
,
γT (t) =
 − √R t ξ(t) 2 , t > T
0

|ξ(s)| ds

where T is a real parameter taking a value in the open interval (0, ∞). By construction,
γT (t) is continuous for all t ≥ 0. To evaluate the probability of the absorber being excited
to the single photon state |1i2 when γ(t) is replaced by γT (t) for some chosen value of T ,
we can solve for (8) to obtain
Z ∞
hn1 (t)σ2,z (t)i = −
|ξ(τ )|2 dτ,
t

and numerically integrate Eqs. (6) and (7), with the given initial conditions.
Let us consider the case√where the input wavepacket ξ(t) is a decaying exponential
function of the form ξ(t) = ce−ct/2 for some positive real constant c. This is the form of
the wavepacket that would be produced at the output of an optical cavity that is initialized
in the 1-photon Fock state. Let us take c = 7.2 × 107 (this is a value that can be realized
in table-top quantum optical experiments, see, e.g., [23]) and allow the absorber to run up
to time t1 = 10/c = 1.3889 × 10−7 . Fig. 1 shows the evolution of hn2 (t)i for T = 0.001t1 ,
T = 0.01t1 , and T = 0.1t1 . It can be seen that for larger T (wider truncation) the
excitation probabillity of the absorber is lower for all t ≥ 0, as can be expected. The steadystate excitation probability is approximately 0.9957, 0.9575, and 0.6037 for T = 0.001t1 ,
T = 0.01t1 , and T = 0.1t1 , respectively.

4

Conclusion

In this work we have considered a single photon absorber with a tunable coupling parameter
to an external travelling single photon field. We analytically derived the exact form of the
time-dependent coupling parameter for perfect absorption of a single photon field of any
temporal wavepacket shape. The ideal modulation function has a singularity at t = 0 which
cannot be attained in real devices, therefore it is approximated with a continuous function
that is truncated to a finite value at t = 0. In a numerical example, we illustrate the effect
of this truncation on the ability of the absorber to perfectly absorb a single photon for a
particular truncation scheme applied to an exponentially decaying wavepacket.
12

Evolution of the probability of single photon excitation in the absorber

1
0.9
0.8
0.7

< n2(t) >

0.6
0.5
0.4
0.3
0.2
0.1
0

0

0.2

0.4

0.6

0.8

1

1.2

Time (s)

1.4
−7

x 10

Figure 1: Time evolution of the excitation probability hn2 i. with the absorber initalized in
the state |0i2 , from t = 0 to t = t1 = 1.3889 × 10−7 , for T = 0.001t1 (top solid black curve),
T = 0.01t1 (middle dashed blue curve), and T = 0.1t1 (bottom dash dotted red curve).

Appendix: General expression for Lt(X1X2)
LL(t) (X1 X2 ) =
=
=
−ı[X1 X2 , H(t)] =
=
=
=

1/2L(t)∗ [X

∗
1
1 X2 , L(t)] + /2[L(t) , X1 X2 ]L(t),
∗
1/2(λ(t)σ
1,+ + γ(t) σ2,+ )[X1 X2 , λ(t)σ1,− + γ(t)σ2,− ]
+1/2[λ(t)∗ σ1,+ + γ(t)∗ σ2,+ , X1 X2 ](λ(t)σ1,− + γ(t)σ2,− )
∗
1/2(λ(t)∗ σ
1,+ + γ(t) σ2,+ )(λ(t)[X1 , σ1,− ]X2 + γ(t)X1 [X2 , σ2,− ])
+1/2(λ(t)∗ [σ1,+ , X1 ] + γ(t)∗ σ2,+ X1 [σ2,+ , X2 ](λ(t)σ1,− + γ(t)σ2,− ).
−ı[X1 X2 , ={γ(t)∗ λ(t)σ2,+ σ1,− }],
−ı[X1 X2 , 1/2ı(γ(t)∗ λ(t)σ2,+ σ1,− − γ(t)λ(t)∗ σ2,− σ1,+ )],
[X1 X2 , 1/2(γ(t)λ(t)∗ σ2,− σ1,+ − γ(t)∗ λ(t)σ2,+ σ1,− )],
∗
1/2γ(t)λ(t)∗ [X X , σ
1
1 2
2,− σ1,+ ] − /2γ(t) λ(t)[X1 X2 , σ2,+ σ1,− ].

13

Combining yields,
Lt (X1 X2 ) = −ı[X1 X2 , H(t)] + LL(t) (X1 X2 ),
= 1/2γ(t)∗ λ(t)[X1 , σ1,− ]σ2,+ X2 + 1/2γ(t)∗ λ(t)X1 [σ2,+ , X2 ]σ1,−
−1/2γ(t)∗ λ(t)(X1 σ1,− X2 σ2,+ − σ1,− X1 σ2,+ X2 ),
+1/2|λ(t)|2 (σ1,+ [X1 , σ1,− ]X2 + [σ1,+ , X1 ]σ1,− X2 ]
+1/2λ(t)∗ γ(t)σ1,+ X1 [X2 , σ2,− ] + 1/2λ(t)∗ γ(t)[σ1,+ , X1 ]X2 σ2,−
+1/2λ(t)∗ γ(t)(X1 X2 σ2,− σ1,+ − σ1,+ σ2,− X1 X2 ),
+1/2|γ(t)|2 (σ2,+ [X2 , σ2,− ]X1 + X1 [σ2,+ , X2 ]σ2,− ),
= γ(t)∗ λ(t)X1 σ1,− [σ2,+ , X2 ] + 1/2|λ(t)|2 (σ1,+ [X1 , σ1,− ] + [σ1,+ , X1 ]σ1,− )X2
+λ(t)∗ γ(t)σ1,+ X1 [X2 , σ2,− ] + 1/2|γ(t)|2 (σ2,+ [X2 , σ2,− ] + [σ2,+ , X2 ]σ2,− )X1 .

References
[1] L.-M. Duan, M. D. Lukin, J. I. Cirac, and P. Zoller, “Long-distance quantum communication with atomic ensembles and linear optics,” Nature, vol. 414, pp. 413–418,
2001.
[2] H. J. Kimble, “The quantum internet,” Nature, vol. 453, pp. 1023–1030, 2008.
[3] J. I. Cirac, P. Zoller, H. J. Kimble, and H. Mabuchi, “Quantum state transfer and
entanglement distribution among distant nodes in a quantum network,” Phys. Rev.
Lett., vol. 78, no. 3221, 1997.
[4] L. Bouten, R. van Handel, and M. R. James, “An introduction to quantum filtering,”
SIAM J. Control Optim., vol. 46, pp. 2199–2241, 2007.
[5] H. Carmichael, An Open Systems Approach to Quantum Optics.
1993.

Berlin: Springer,

[6] Q. Y. He, M. D. Reid, and P. D. Drummond, “Digital quantum memories with symmetric pulses,” Opt. Express, vol. 17, no. 12, pp. 9662–9668, 2009.
[7] J. Dilley, P. Nisbet-Jones, B. W. Shore, and A. Kuhn, “Single-photon absorption in
coupled atom-cavity systems,” Phys. Rev. A, vol. 85, no. 023834, 2012.
[8] A. N. Korotkov, “Flying microwave qubits with nearly perfect transfer efficiency,”
Phys. Rev. B, vol. 84, no. 014510, 2011.
[9] Y. Yin, Y. Chen, D. Sank, P. J. J. O’Malley, T. C. White, R. Barends, J. Kelly,
E. Lucero, M. Mariantoni, A. Megrant, C. Neill, A. Vainsencher, J. Wenner, A. N. .
Korotkov, A. N. Cleland, and J. M. Martinis, “Catch and release of microwave photon
states,” Phys. Rev. Lett., vol. 110, no. 107001, 2013.
[10] N. Yamamoto and M. R. James, “Zero-dynamics principle for perfect quantum memory in linear networks,” New J. Phys., vol. 16, no. 073032, pp. 1–30, 2014.

14

[11] R. L. Hudson and K. R. Parthasarathy, “Quantum Ito’s formula and stochastic evolution,” Commun. Math. Phys., vol. 93, pp. 301–323, 1984.
[12] K. Parthasarathy, An Introduction to Quantum Stochastic Calculus.
Birkhauser, 1992.
[13] P.-A. Meyer, Quantum Probability for Probabilists, 2nd ed.
Springer-Verlag, 1995.

Berlin:

Berlin-Heidelberg:

[14] L. Bouten and R. van Handel, “On the separation principle of quantum control,” in
Quantum Stochastics and Information: Statistics, Filtering and Control (University
of Nottingham, UK, 15 - 22 July 2006), V. P. Belavkin and M. Guta, Eds. Singapore:
World Scientific, 2008, pp. 206–238.
[15] J. Gough and M. R. James, “The series product and its application to quantum
feedforward and feedback networks,” IEEE Trans. Autom. Control, vol. 54, no. 11,
pp. 2530–2544, 2009.
[16] R. Loudon, The Quantum Theory of Light, 3rd ed. Oxford University Press, 2000.
[17] G. J. Milburn, “Coherent control of single photon states,” EPJ ST, vol. 159, no. 1,
pp. 113–117, 2008.
[18] J. E. Gough, M. R. James, and H. I. Nurdin, “Quantum filtering for systems driven by
fields in single photon states and superposition of coherent states using non-markovian
embeddings,” Quantum Inf Process, vol. 12, pp. 1469–1499, 2013.
[19] J. Gough, M. R. James, H. I. Nurdin, and J. Combes, “Quantum filtering for systems
driven by fields in single-photon states or superposition of coherent states,” Phys. Rev.
A, vol. 86, p. 043819, 2012.
[20] J. E. Gough, M. R. James, and H. I. Nurdin, “Quantum master equation and filter
for systems driven by fields in a single photon state,” in Proceedings of the 50th IEEE
Conference on Decision and Control (CDC), 2011, pp. 5570–5576.
[21] Y. Pan, G. Zhang, and M. R. James, “Analysis and control of quantum finite-level
systems driven by single-photon input states,” Automatica, vol. 69, pp. 18–23, 2016.
[22] J. E. Gough, M. R. James, and H. I. Nurdin, “Squeezing components in linear quantum
feedback networks,” Phys. Rev. A, vol. 81, pp. 023 804–1– 023 804–15, 2010.
[23] S. Iida, M. Yukawa, H. Yonezawa, N. Yamamoto, and A. Furusawa, “Experimental
demonstration of coherent feedback control on optical field squeezing,” IEEE Trans.
Autom. Control, vol. 57, no. 8, pp. 2045–2050, 2012.

15

