## Coherent trapping in small quantum networks

T.J. Pope, , J. Rajendran 1 2 3 , , A. Ridolfo, 2 3 , , E. Paladino 2 3 4 , , , F.M.D. Pellegrino 2 3 , , and G. Falci 2 3 4 , ,

- 1 School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne NE1 7RU
- 2 Dipartimento di Fisica e Astronomia 'Ettore Majorana', Via S. Sofia 64, Ctanaia (I)
- 3 Istituto Nazionale di Fisica Nucleare, Sezione di Catania, Via S. Sofia 64, Ctanaia (I)
- 4 Consiglio Nazionale delle Ricerche, Istituto di Microelettronica e Microsistemi, UoS Universit`, Via S. Sofia 64, Ctanaia (I) a

E-mail: gfalci@dmfci.unict.it

Abstract. We consider a three-node fully connected network (Delta network) showing that a coherent population trapping phenomenon occurs, generalizing results for the Lambda network known to support a dark state. Transport in such structures provides signatures of detrapping, which can be triggered by external controls. In the presence of an environment it turns out to be sensitive to its Markovianity. Adiabatic modulation of the system's parameters may yield coherent population transfer, analogous to the stimulated Raman adiabatic passage phenomenon. Robustness of this protocol against non-adiabatic transitions is studied. Coherent nanostructures where these phenomena are relevant for quantum transport and quantum protocols are suggested.

## 1. Introduction

Coherent trapping is a well known phenomenon in atomic physics, occurring when various ground states are coupled to a common upper level by multiple resonant laser beams [1, 2], the simplest instance being a three-level system called Lambda (Λ) network (see Fig. 1). As a result population is trapped in a subspace, the upper level being never populated. Since this prevents radiative decay the trapped state is often called dark state. Adiabatic manipulations of the dark state may allow faithful and robust population transfer by techniques as the stimulated emission Raman adiabatic passage (STIRAP) [3] or superadiabatic extensions [4], and are nowadays of interest for applications in solid-state quantum information [5-11].

In this work we study a Delta (∆) systems, a three-site quantum network depicted in Fig. 1. Each site defines a basis state {| 0 〉 , | 1 〉 , | 2 〉} . We will show that under suitable conditions also this more general network admits a stable state, trapped in the subspace

Coherent trapping in small quantum networks

|

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01145v1-figure-01ce9b9d-f322-454f-ac6e-25ed7af53f0b.jpg)

0

〉

glyph[negationslash]

Figure 1. Left panel schematic of the ∆ system consisting of three levels and three couplings. One of the levels may dacay to a sink with a rate Γ 2 S . This network may admit a stable 'trapped' eigenstate | D 〉 , the population being confined in the subspace {| 0 〉 , | 1 〉} , if the condition ¯ = 0 is fullfilled (see text). δ This is an effective model for a simple quantum network, or for a three-level atom driven by three coherent fields. In the Λ system one of the couplings vanishes, Ω 0 = 0. Right panel: the same network in a 'dressed' basis, emphasizing that detrapping occurs only for ¯ = 0. δ

{| 0 〉 , | 1 〉} . The Hamiltonian of the ∆ systems is written as

$$\cdot \cdot \cdot \cdot \cdot \\ H = \left ( \begin{array} { c c c } 0 & \Omega _ { 0 } & \Omega _ { p } \\ \Omega _ { 0 } ^ { * } & \delta & \Omega _ { s } \\ \Omega _ { p } & \Omega _ { s } & \delta _ { p } \end{array} \right )$$

where we have chosen the arbitrary phases of the basis states in such a way that the two off-diagonal elements Ω p and Ω s are real. Here we use the notation of quantum optics where Eq.(1) is the Hamiltonian of a multilevel atom driven by three nearly resonant lasers expressed in a multiple rotating frame, in the rotating wave approximation. Diagonal elements represent detunings of the fields while off-diagonal elements are the Rabi energies, related to the amplitudes (see 4) of the driving fields. When the system is driven by only two lasers special configurations are obtained: the Λ scheme for Ω 0 = 0, the ladder scheme (Ω p = 0) and the Vee scheme (Ω s = 0).

Seen as a Hamiltonian in the laboratory frame Eq.(1) describes a quantum network where for instance three localized states with on-site energies δ and δ p (relative to | 0 ) and with Ω 〉 i being the tunneling amplitudes connecting the sites. The physics of population trapping also plays an important in solid-state quantum networks [12, 13] and has been invoked to explain peculiar quantum transport properties in nanostructures, from architectures of quantum dots [14] to light-harvesting systems in biological complexes [15, 16].

We will show in Sec. 2 that also in the ∆ system an 'unconventional' trapping phenomenon occurs. We discuss in Sec. 3 the implications for transport in solid-state quantum networks and in Sec. 4 applications to the dynamics of driven three-level atoms, who also provide insight on the robustness of unconventional trapping.

## 2. Trapped states

We now look for trapped states of the network, namely if eigenstates of the form

$$\left | D \right \rangle = c _ { 0 } \left | 0 \right \rangle + c _ { 1 } \left | 1 \right \rangle$$

exist. The probability of detecting the system in state | 2 〉 vanishes, 〈 2 | D 〉 = 0, therefore the system is trapped in a one-dimensional subspace, despite the non vanishing transition amplitudes. In the quantum optical version the system is trapped in | D 〉 despite the fields tend to trigger transitions to the state | 2 . 〉 Since this latter often decays radiatively, when the system is trapped no fluorescence is observed, therefore | D 〉 is called a dark state.

By imposing that | D 〉 is an eigenstate we find the equations

$$\begin{cases} \ - z _ { D } c _ { 0 } + \Omega _ { 0 } c _ { 1 } = 0 \\ \ \Omega _ { 0 } ^ { * } c _ { 0 } + ( \delta - z _ { D } ) c _ { 1 } = 0 \\ \ \Omega _ { p } c _ { 0 } + \Omega _ { s } c _ { 1 } = 0 \end{cases}$$

z D being the relative eigenvalue z . The last equation gives c 1 = -(Ω p / Ω ) s c 0 , yielding the form of the trapped state

$$| D \rangle = \frac { \Omega _ { s } \, | 0 \rangle - \Omega _ { p } \, | 1 \rangle } { \Omega _ { R M S } } = \colon \cos \theta \, | 0 \rangle - \sin \theta \, | 1 \rangle$$

where Ω RMS = √ Ω +Ω and the mixing angle is given by tan 2 p 2 s θ = Ω p / Ω . This trapped s state is an eigenstate only if the other two equations are satisfied, namely

$$\delta = \frac { \Omega _ { 0 } ^ { * } \Omega _ { s } ^ { 2 } - \Omega _ { 0 } \Omega _ { p } ^ { 2 } } { \Omega _ { p } \Omega _ { s } } = \colon \delta _ { 0 }$$

$$z = - \frac { \Omega _ { 0 } \dot { \Omega _ { p } } } { \Omega _ { s } } = \colon z _ { D }$$

Notice that besides Ω , also 0 δ and δ p may be complex. For simplicity here we focus on the case of real matrix elements. The network can be digonalized and we find for the eigenvalues the compact form

$$z _ { \pm } = \frac { 1 } { 2 } \left [ \delta _ { p } + \frac { \Omega _ { 0 } \Omega _ { s } } { \Omega _ { p } } \pm \Omega _ { A T } \right ]$$

where a quantity analogous to the Autler-Townes splitting appears

$$\Omega _ { A T } = \sqrt { \tilde { \delta } _ { p } ^ { 2 } + 4 \Omega _ { R M S } ^ { 2 } } \quad ; \quad \tilde { \delta } _ { p } = \delta _ { p } - \frac { \Omega _ { 0 } \Omega _ { s } } { \Omega _ { p } }$$

Eigenvectors can be conveniently written as

$$\left | + \right \rangle = \sin \Phi \left | B \right \rangle + \cos \Phi \left | 2 \right \rangle \ \ ; \ \left | - \right \rangle = \cos \Phi \left | B \right \rangle - \sin \Phi \left | 2 \right \rangle$$

where | B 〉 = sin θ | 0 〉 +cos θ | 1 〉 is called bright state, and forms a basis together with | D 〉 and | 2 . 〉 Finally the mixing angle is given by

$$\tan 2 \Phi = \frac { 2 \Omega _ { R M S } } { \tilde { \delta } _ { p } }$$

Figure 2. Numerically calculated final sink population for the Λ network. Here Ω = Ω s p = 1 thus Ω RMS = √ 2. Left: sink population for for a range of δ and δ p , with Γ 2 S = 1, compared to the analytical prediction (10). Right: Numerically calculated final sink population for the Λ network for a range of δ and Γ 2 S , where Ω = 1 and δ p = Ω 0 = 0, compared to the analytical prediction for sink.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01145v1-figure-f3a6caec-fb9d-4da3-927d-7a649429c9c3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01145v1-figure-f0fbbb95-9932-4f1d-bd88-f147909be895.jpg)

glyph[negationslash]

A well known special solution is obtained for Ω 0 = 0, which describes a Lambda system. The dark state is obtained at zero two-photon detuning, δ = 0, and the eigenvalue is z D = 0, independent on all the other parameters. This latter is a strong indication of the robustness of the dark state for fluctuations of the parameters. We will see that robustness is kept also for Ω 0 = 0, despite of the fact that z D = 0. Another special case is the symmetric Delta system, where Ω s = Ω , the dark state being obtained for p δ = 0, with eigenvalue z d = -Ω . 0

glyph[negationslash]

Notice that eigenvalues Eq. (5) and eigenvectors Eq. (7) have the same structure known for Lambda systems [3], with redefined parameters given by Eqs.(6,8).

## 3. Trapped states in nanostructures and detrapping

## 3.1. Transport in small quantum networks

Dark states in Lambda and Delta systems play a crucial role in the dynamics of quantum networks. The Hamiltonian (1) describes three sites where an electron or an exciton can be trapped, the off-diagonal elements describing tunneling. We also allow decay from the site | 2 〉 outside of the system. To understand the dynamics it is convenient to represent the Hamiltonian in the basis {| D , B 〉 | 〉 , | 2 〉}

$$\begin{array} {$$

where ¯ = δ δ -δ 0 is the deviation from the trapping condition.

Weuse this representation to study detrapping-induced transport in small networks. The system is initialized in the trapped state | D 〉 and eventually decays to a sink with a rate Γ 2 S (see Fig.). The sink is accounted for by adding an extra level to the network

or by adding an imaginary part -i Γ 2 s / 2 to δ p . As it is clear from Eq.(9) a nonzero ¯ triggers transitions δ | D 〉 → | B 〉 . As a consequence also | 2 〉 populates and decays to the sink. Therefore detrapping is detected by populating the sink. In Fig.2 we show population of the sink P sink ( ) t = 1 -Tr[ ρ t ( )] where ρ t ( ) is the density matrix of the network, at a large enough time t = t f , for a symmetric network Ω s = Ω p in Lambda configuration Ω 0 = 0, where the trapping condition is δ 0 = 0. In quantum optics the minimum in the curves Eq.(9) at δ = δ 0 is referred as the 'dark resonance'.

In the limit δ glyph[lessmuch] Ω RMS glyph[lessorsimilar] Γ 2 S and a simple approximate analytic form may be derived by perturbation theory (see Appendix Appendix A)

$$\text{Tr} [ \rho ( t ) ] \approx \rho _ { D$$

where ρ DD ( ) is the population of the trapped state. t The formula is derived assuming that coupling between the dark and the bright state is small, whereas the coupling between the bright state and site 2 is relatively large as well as the decay to the sink. Decay is triggered by the level splitting δ and as a result the population of the network decreases.

The left panel of Fig.2 this result is plotted against the numerical calculation, and we see good agreement for range of | δ | &lt; | Ω p,s | whereas for larger δ the processes | D 〉 → | B 〉 and | B 〉 → | 〉 2 are no longer separable. In the right panel the effect of changing the decay rate Γ 2 S is shown. Here, the prediction in Eq. (10) still holds and we see that by increasing Γ 2 S , the width of the dark resonance is reduced.

The right panel of Fig.3 shows the effect of changing the decay rate Γ 2 S . Here, the prediction in Eq. 2.19 still holds and we see that by increasing Γ 2 S , the width of the dark resonance is reduced.

Results a Delta system initialized in | 0 , are shown in the left panel of Fig.(2). 〉 It is apparent that dark resonances in the population appear for finite δ , since the trapping condition δ = δ 0 now depends on Ω . 0 Notice that population at the minima δ = δ 0 is now nonzero since only a fraction |〈 0 | D 〉| 2 = Ω 2 s / Ω 2 RMS of the population is initially trapped.

## 3.2. Effect of Noise

In this section we discuss the effect of noise, focusing on the most relevant process destroying coherence, namely dephasing in the trapped state. In addition we consider decay of population from | 2 〉 to a sink, described by a fourth state | S 〉 . The effect of a Markovian environment is accounted for by the following Lindblad equation

$$\dot { \rho } \left ( t \right ) = - \frac { i } {$$

where decay to the sink and the associated dephasing is described by the dissipator

$$\mathcal { D } _ { s i n k } \left [ \rho \left ($$

Figure 3. Left0: Numerically calculated final sink population for the Λ network for a range of δ and Γ 2 S , where Ω = 1 and δ p = Ω 0 = 0, compared to the analytical prediction for sink. Left: Right: Dark state shifting with level splitting δ in the asymmetric ∆ network.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01145v1-figure-abfb32ad-7def-4815-917c-ab411400399c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01145v1-figure-b610070b-9bcb-4c80-ad3b-09af686a342d.jpg)

where ρ is the density matrix of the four state network and curly brackets indicate the anticommutator. Pure dephasing [17] of the trapped state is accounted for by the dissipator

$$\mathcal { D } _ { d e p h } \left [ \rho \left ($$

We also consider detrapping due to non-Markovian calssical noise, which is accounted for by a stochastic process x t ( ) coupled to the Hamiltonian

$$H \, \rightarrow \, H + x \, ( t ) \, | 1 \rangle \langle 1 | \,.$$

the evolution of ρ t ( ) being obtained by averagin over all realizations of x t ( ). here we consider the Ornstein-Uhlenbeck (OU) classical stochastic process [18, 19], characterized by the un-normalized autocorrelation function; 〈 x t ( ) x t ( + t ′ ) 〉 = σ 2 ou e -t /τ ′ , where τ is the correlation time. The power spectrum of the noise is given by

$$\text{ne correlation time.                                                                                                                                                                                                        
                                                                                                                                                                                                        \mathcal { S } ( \omega ) \, = \, & \int _ { - \infty } ^ { \infty } d t ^ { \prime } \, e ^ { i \omega t ^ { \prime } } \, \left \langle x \left ( t \right ) x \left ( t + t ^ { \prime } \right ) \right \rangle \\ = \, & 2 \, \sigma _ { o u } ^ { 2 } \int _ { 0 } ^ { \infty } d t \, e ^ { - t / \tau } \, \cos \left ( \omega t \right ) = \frac { 2 \, \sigma _ { o u } ^ { 2 } \tau } { 1 + \omega ^ { 2 } \tau ^ { 2 } }. \\ \text{vhoro wo have a command that the noise is atianary } \, / \, r \left ( t \right \rangle \, r \left ( t \right \rangle \right \rangle - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, /\, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ) \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t ^ { \prime } \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( \Omega \right \rangle \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, \right \rangle \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \,. \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \rll \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rigs \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \r ) \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \r | \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rparable \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \ragle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \r]) \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \ralg \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rlash \, \ll - \, / \, r \left ( t \right \rangle \, \ll - \, / \, r \left ( t \right \rangle \, \ll -$$

where we have assumed that the noise is stationary, 〈 x t ( ) x t ( + t ′ ) 〉 = 〈 x (0) x t ( ′ ) 〉 . In order to understand the role of non-Markovianity, we will compare noises which would give the same decay time of coherence if treated in by the standard Bloch-Redfield master equation theory [20], since this latter approach yields correct decay rates in the Markovian limit. A scale for the decay time is γ ∼ S (Ω ), obtained by arguing that 0 local dephasing triggers transitions between 'dressed' eigenstates | φ n 〉 , whose splitting is Ω . ˜ 0 To this end, for a given γ we consider a OU process with width σ 2 ou given by

$$\sigma _ { o u } ^ { 2 } = \frac { 1 + \Omega _ { 0 } ^ { 2 } \tau ^ { 2 } } { \tau } \gamma, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

In the simplest case, where Ω s = Ω p = τ = 1 and Ω 0 = δ p = 0 we identify σ 2 ou = γ .

Figure 4. Effects of noise on the ∆ network, where Ω 0 = δ p = 1. The left-hand panel shows the effect of white noise and the right-hand panel shows the effect of the Ornstein Uhlenbeck process with a correlation time of τ = 1 ps. In both panels, the population of the sink at time t = 100 ps is plotted with solid lines and the sum of the absolute values of the coherences at time t = 10 ps is plotted with dashed line.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01145v1-figure-629ed09e-cab6-4a3d-8bba-a02cf0bc2aa0.jpg)

Figure 4 shows the effect of each noise source on the sink population and on the sum of the coherences as a function of the level splitting δ for the ∆ network, where Ω 0 = δ p = 1. The left-hand panel shows the effect of white noise characterized by a range of γ 1 and the right-hand panel shows the effect of the OU process characterized by a range of noise widths, σ 2 ou . The measurement for the coherence is taken at 10 ps, and the measurement for the sink population is taken at 100 ps. It is seen that while decoherence in the trapped subspace in general weakens the trapping phenomenon an environment with memory is less effective in this respect, although in both cases coherences are suppressed in a comparable fashion. Non-Markovianity is expected to produce more pronounced signatures if also fluctuations of the splitting δ p are accounted for [21] together with their correlations, a problem already explited in superconducting quantum architectures [22], where this analysis allows to determine or design optimal operating points where decoherence is minimized [23-25].

## 4. Coherent population transfer in Delta systems

Adiabatically modulated trapped states can be used for implementing quantum operation, the simplest being is population transfer [3]. This can be achieved by

Figure 5. Population transfer by STIRAP in the Λ ( α = 0 in Eq. (17)) and in the ∆ ( α = 1) systems for Ω T = 15, τ = 0 7 . T and γ = 0. Top left panel: the pulses Ω s,p ( ) t for standard STIRAP (dashed lines) and the additional pulses Ω ( ) 0 t and δ t ( ) for unconventional STIRAP (full lines). Top right panel: instantaneous eigenvalues for ( α, δ, δ p ) = (0 0 0) (dashed lines) and for ( , , α, δ, δ p ) = (1 , δ 0 ( ) t , 0) (full lines), and the eigenvalue z -( ) t for ( α, δ, δ p ) = (1 , δ 0 ( ) t , -2 Ω) (dot dashed line), as given by Eqs.(4,5). Bottom left panel: population histories for for ( α, δ, δ p ) = (0 , 0 0) , (dashed gray lines) and for ( α, δ, δ p ) = (1 , δ 0 ( ) t , 0) (full lines); these latter showing imperfect population transfer due to Zener tunneling occurring when eigenvalues z D ( ) t and z -( ) t are degenerate. Bottom right panel: population histories for ( α, δ, δ p ) = (1 , δ 0 ( ) t , -2 Ω) (full lines), for ( α, δ, δ p ) = (0 , δ 0 ( ) t , -2 Ω) (gray dashed lines) and ρ 22 ( ) t for ( α, δ, δ p ) = (0 0 0); here the variuos , , ρ 22 ( ) are magnified by a factor 30, showing t that the unconventional pattern for STIRAP is more robust against transitions to the intermediale level.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01145v1-figure-9660a0ce-439d-4e63-b749-11058336a8bf.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01145v1-figure-a1844327-b96f-4314-88fd-2067c0919884.jpg)

suitable modulation of the parameters of the Hamiltonian (1). Physically it describes a three-level atom driven by external fields W t k ( ) = 2Ω k cos( ω t k ), for k = 0 , s, p . The angular frequencies ω i of the fields are chosen to approximately match the level splittings of the atom, therefore they address specific transitions. In particular we take ω p = glyph[epsilon1] 2 -glyph[epsilon1] 0 -δ p , where | δ p | glyph[lessmuch] ω p is the pump field detuning. We also take slightly detuned, ω s = glyph[epsilon1] 2 -glyph[epsilon1] 1 -δ s , whereas ω 0 = glyph[epsilon1] 1 -glyph[epsilon1] 0 is not detuned. We define the twophoton detuning δ = δ p -δ s . Then Eq.(1) is the Hamiltonian in the rotating wave approximation as seen from a rotating frame.

We allow for a time dependence of the Rabi frequencies Ω ( ) on a time scale much k t longer than the field period 2 π/ω k . Therefore the Hamiltonian (1) may determine an adiabatic dynamics of the system. In the Λ configuration at two-photon resonance, δ = 0 this allows to achieve complete population transfer | 0 〉 → | 1 〉 by adiabatically following

Figure 6. Stability against fluctuations of detunings of population transfer by STIRAP, for Ω T = 15, τ = 0 7 . T and γ = 1 /T and various α = 0 0 3 0 5 1 0. , . , . , . We plot lines corresponding to P 1 ( t f ) = 0 95. . The stability region enlarges for increasing α showing that that the ∆ system protocol is more robust. A negative δ p increases the stability for larger α since it removes the degeneracy of z D ( ) and t z -( ) enforcing t the adiabatic pattern | 0 〉 → | 1 . 〉 Also a positive δ p may increase population transfer efficiency, due to pattens combining adiabatic dynamics an diabatic transitions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01145v1-figure-0a8b570d-2619-4ad1-b30b-a6b2b318ff88.jpg)

the trapped state Eq.(2). To this end the fields' amplitudes are adiabatically modulated in the 'counterintuitive' sequence, i.e. Ω ( ) is shined before Ω ( ), a protocol named s t p t STIRAP. In our work we consider for simplicity Gaussian pulses of equal amplitude

$$\Omega _ { s } ( t ) = \Omega \, e ^ { - \left ( \frac { t + \tau } { T } \right ) ^ { 2 } } \ \ ; \ \ \Omega _ { p } ( t ) = \Omega \, e ^ { - \left ( \frac { t - \tau } { T } \right ) ^ { 2 } }$$

where T is the width of the pulse and the delay is τ glyph[lessorsimilar] T , shown in Fig. 5. The protocol is very selective and robust, being almost insensitive to fluctuation of all parameters, except δ (see Fig. 6). Population transfer by STIRAP [3] in solid-state devices has been proposed in systems of artificial atoms based on superconductors [5, 6, 26, 27] and semiconductors [28] and has been recently demonstated in Josephson systems [29, 30].

STIRAP can also be implemented in the ∆ system, as shown in Fig. 6 (bottom panels). We have chosen the following form for the coupling

$$\Omega _ { 0 } ( t ) = \alpha \, \frac { \Omega _ { s } ( t ) \, \Omega _ { p } ( t ) } { \Omega } = \alpha \, \Omega \, e ^ { - \frac { 2 } { T ^ { 2 } } ( t ^ { 2 } + \tau ^ { 2 } ) }$$

population transfer being achieved by modulation of the two-photon detuning δ t ( ) = δ 0 ( ) where from Eq.(3) t

$$\delta _ { 0 } ( t ) = \frac { \Omega _ { 0 } ( t ) } { \Omega _ { s } ( t ) \Omega _ { p } ( t ) } \left [ \Omega _ { s } ^ { 2 } ( t ) - \Omega _ { p } ^ { 2 } ( t ) \right ] = \frac { \alpha } { \Omega } \left [ \Omega _ { s } ^ { 2 } ( t ) - \Omega _ { p } ^ { 2 } ( t ) \right ]$$

The modulated fields are shown in Fig. 5. Also the instantaneous eigenvalues are plotted, showing that an 'unconventional' adiabatic pattern for population transfer

exists also in the Delta system with modulated δ t ( ). It is interesting to notice that for δ p = 0 instantaneous eigenstates may be degenerate at some time ( z D and z -in the top right panel of Fig. 5), therefore efficiency of the protocol is lowered by Zener tunneling (see Fig. 5, lower left panel). Efficiency is recovered by using a constant δ p &lt; 0 (see Fig. 5, lower right panel) which eliminates the crossing. In this case not only the final population of the target state is larger, but population of the intermediate state (which in many cases should be avoided) is lowered enforcing trapping during the whole protocol.

Lower population of the intermediate state is a fingerprint of the fact that 'unconventional' STIRAP in ∆ configuration is more robust than the conventional Λ-STIRAP. This is clearly shown in Fig. 6, where we study sensitivity of the efficiency P 1 ( t f ) to parametric fluctuations of the detunings. Thus we allow for constant deviations of the detunings from their reference values, namely δ t ( ) = δ 0 ( ) + ¯ and t δ δ p = ¯ . δ p In order to check that population is trapped during the whole procedure, we introduce decay of the intermediate level δ p → δ p -iγ/ 2 which lowers P 1 ( t f ) if | 2 〉 is always occupied. As it is apparent the presence of the field Ω 0 enlarges the stability region. This latter is moved towards negative values of ¯ δ p which ensure that crossings of adiabatic eigenvalues are eliminated. Actually efficiency may increase also for ¯ δ p &gt; 0 (effect not shown in the figure), due to the combination of adiabatic evolution and diabatic transition where Zener tunneling restores the pattern leading to population transfer. An analogous phenomenon is known to occurs in Λ STIRAP for δ nonvanishing [3, 31], which in general guarantees population transfer also for δ sufficiently close to the trapping condition Eq.(3).

A ∆ network can be implemented in superconducting architectures by working away from symmetry points [8] but decoherence is minimized at optimal points, where one of the three couplings has to be implemented by a two-photon process [32]. Other protocols combining slow modulation of detuning with resonant ac pulses have also been proposed [33], which in artificial atoms are expected to be robust against phase noise since control is operated at microwave frequencies.

## 5. Conclusions

We have studied coherent trapping in three-site quantum networks. We have shown that under suitable conditions ∆ networks have trapped eigenstates analogous to what is found in Λ (and also Ladder and Vee) networks.

Trapping and detrapping phenomena could be relevant in several physical systems, as artificial networks of semiconducting quantum dots [14] or molecular complexes in certain biological systems [16], where the ∆ system describes either part or the whole network, and for small circuit-QED quantum networks of artificial atoms where the δ system may provide an effective low-energy description.

In driven ∆ network the physics of trapping may allow adiabatic coherent population transfer. Delta systems in superconductors have been studied [8] and recently

implemented to demonstrate superadiabatic protocols. This configuration of couplings may also be relevant for adiabatic dynamics in circuit-QED architectures [34, 35].

The rich physics of the dark resonance in ∆ systems can be used as a tool for sensing the properties of a noisy environment [22] yielding for instance direct information on dephasing and on its non-Markovianity [17] or to find optimal operating point in quantum devices [23-25].

## Appendix A. Derivation of the current in the sink

To derive Eq.(10) we use perturbation theory in the two-photon detuning ¯ , δ by writing the Hamiltonian Eq.(1) using the eigenbasis for δ = δ 0

$$H = \sum _ { k = D, \pm } z _ { k } \, | k \rangle \langle k | + \bar { \delta } \ | 1 \rangle \langle 1 |$$

In the limit δ glyph[lessmuch] Ω RMS glyph[lessorsimilar] Γ 2 S we can suppose that |±〉 is almost depopulated and Tr[ ρ t ( )] ≈ |〈 D ψ t | ( ) 〉| 2 . We prepare the system in | D 〉 . Then for small ¯ in leading order δ in ¯ δ/ Ω RMS the only nonvanishing amplitude is 〈 D ψ t | ( ) 〉 ≈ 1 and

$$\langle D | \psi ( t ) \rangle \approx e ^ { - i ( z _ { D } + \delta z _ { D } ) t }$$

where δz D is the correction to the eigenvalue z D . In the illustrative case of the symmetric Λ system, Ω 0 = 0 and Ω p = Ω s = Ω RMS / √ 2 using the eigenstates (2,7) and the eigenvalues (4,5) we easily evaluate

$$\delta z _ { D } = \frac { \bar { \delta } } { 2 } - \frac { \bar { \delta } ^ { 2 } } { 4 } \left [ \frac { \sin ^ { 2 } \Phi } { z _ { + } } + \frac { \cos ^ { 2 } \Phi } { z _ { - } } \right ] = \frac { \bar { \delta } } { 2 } - \frac { \bar { \delta } ^ { 2 } \delta _ { p } } { 4 \Omega _ { R M S } ^ { 2 } }$$

Finally letting δ p → δ p -i Γ 2 S / 2, we obtain

$$| \langle D | \psi ( t ) \rangle | \approx e ^ { 2 \, \Im ( \delta z _ { D } ) \, t } = e ^ { - \frac { \tilde { \delta } ^ { 2 } \Gamma _ { 2 S } } { 4 \Omega _ { R M S } ^ { 2 } } t }$$

which is the desired result. Eq.(10) may be derived more rigorously and extended beyond the small ¯ by adiabatic elimination of the bright state. δ

## References

- [1] Arimondo E O G 1976 Nuovo Cimento
- [2] H R Gray H, Whitley R and Stroud C 1978 Opt. Lett. 3 218-220
- [3] Vitanov N V, Rangelov A A, Shore B W and Bergmann K 2017 Rev. Mod. Phys. 89 (1) 015006
- [4] Giannelli L and Arimondo E 2014 Phys. Rev. A 89 (3) 033419
- [5] Liu Y x, You J Q, Wei L F, Sun C P and Nori F 2005 Phys. Rev. Lett. 95 (8) 087001
- [6] Siewert J, Brandes T and Falci G 2009 Phys. Rev. B 79 024504

- [7] Veps¨ al¨inen A, Danilin S, Paladino E, Falci G and Paraoanu G S 2016 a Photonics 3 62
- [8] You J Q and Nori F 2011 Nature 474 589
- [9] Zhou L, Yang L P, Li Y and Sun C P 2013 Phys. Rev. Lett. 111 (10) 103604
- [10] Falci G, Ridolfo A, Di Stefano P and Paladino E 2019 Sci. Rep. 9 9249
- [11] Ridolfo A, Falci G, Pellegrino F and Paladino E 2019 Eur. Phys. J. Special Topics 227 21832188
- [12] M¨ ulken O and Blumen A 2011 Phys. Rep. 502 37?87
- [13] Bianconi G 2015 EPL 111 56001
- [14] Contreras-Pulido L, Bruderer M, Huelga S and Plenio M 2014 New Jour. Phys., Volume 16, November 2014 16 113061
- [15] Lambert N, Chen Y N, Cheng Y C, Li C M, Chen G Y and Nori F 2013 Nature Physics 9 10-18
- [16] Mohseni M, Omar Y, Engel G and Plenio M (eds) 2104 Quantum Effects in Biology (Cambridge University Press)
- [17] Paladino E, Galperin Y, Falci G and Altshuler B 2014 Rev. Mod. Phys. 86 361-418
- [18] Gardiner C 1985 Stochastic methods (Springer-Verlag, Berlin-Heidelberg-New York-Tokyo)
- [19] Bartosch L 2001 International Journal of Modern Physics C 12 851
- [20] Cohen-Tannoudji C, Dupont-Roc J and Gryndberg G 1998 Atom-Photon Interactions: Basic Processes and Applications (Wiley-VCH,)
- [21] Falci G, La Cognata A, Berritta M, D'Arrigo A, Paladino E and Spagnolo B 2013 Phys. Rev. B 87 (21) 214515-1,214515-13
- [22] D'Arrigo A, Mastellone A, Paladino E and Falci G 2008 New Jour. Phys. 10
- [23] Chiarello F, Paladino E, Castellano M G, Cosmelli C, D'Arrigo A, Torrioli G and Falci G 2012 New Jour. Phys, 14
- [24] Paladino E, Mastellone A, D'Arrigo A and Falci G 2010 Phys. Rev. B 81
- [25] Paladino E, D'Arrigo A, Mastellone A and Falci G 2011 New Jour. Phys. 13
- [26] Siewert J, Brandes T and Falci G 2006 Optics Communications 264 435 - 440
- [27] Falci G, Di Stefano P, Ridolfo A, D'Arrigo A, Paraoanu G and Paladino E 2017 Fort. Phys. 65 1600077
- [28] Menchon-Enrich R, Benseny A, Ahufinger V, Greentree A, Busch T and Mompart J 2016 Rep. Prog. Phys. 79 074401
- [29] Kumar K, Veps¨ al¨inen A, Danilin S and Paraoanu G 2016 a Nat. Comm. 7 10628
- [30] Xu H K, Song C, Liu W Y, Xue G M, Su F F, Deng H, Tian Y, Zheng D N, Han S, Zhong Y P, Wang H, Liu Y x and Zhao S P 2016 Nature Communications 7 11018
- [31] Falci G, Berritta M, Russo A, D'Arrigo A and Paladino E 2012 Physica Scripta T151

- [32] Di Stefano P G, Paladino E, Pope T J and Falci G 2016 Phys. Rev. A 93 (5) 051801
- [33] Di Stefano P G, Paladino E, D'Arrigo A and Falci G 2015 Phys. Rev. B 91 (22) 224506
- [34] Wallraff A, Schuster D I, Blais A, Frunzio L, Huang R S, Majer J, Kumar S, Girvin S and Schoelkopf R 2004 Nature 421 162-167
- [35] Di Stefano P G, Alonso J J, Lutz E, Falci G and Paternostro M 2018 Phys. Rev. B 98 (14) 144514
T. J. Pope J. Rajendran A. Ridolfo E. Paladino F. M. D. Pellegrino G. Falci
Coherent trapping in small quantum networks
We consider a three-node fully connected network (Delta network) showing that
a coherent population trapping phenomenon occurs, generalizing results for the
Lambda network known to support a dark state. Transport in such structures
provides signatures of detrapping, which can be triggered by external controls.
In the presence of an environment it turns out to be sensitive to its
Markovianity. Adiabatic modulation of the system's parameters may yield
coherent population transfer, analogous to the stimulated Raman adiabatic
passage phenomenon. Robustness of this protocol against non-adiabatic
transitions is studied. Coherent nanostructures where these phenomena are
relevant for quantum transport and quantum protocols are suggested.