# Open Quantum Systems/Dissipative preparation of many-body states

Resource type: this resource contains a lecture or lecture notes. |

In this chapter, we will see how one can use the interplay between coherent and dissipative dynamics can give rise to the dissipative preparation of quantum many-body states with useful and interesting properties. For simplicity, we focus on Markovian master equations in Lindblad form, given by

#### Stationary states

[edit | edit source]The stationary states of a master equation are found by solving the equation . In most cases, the stationary states do not have an imaginary component of the dynamics so this is equivalent to solving the equation . The detailed nature of the stationary state of course depends on the details of the master equation, but we can make a few general statements.

For example, consider the case where all jump operators are Hermitian, i.e., . This corresponds to a pure dephasing where only the off-diagonal elements of the density matrix get damped out, but there is no reshuffling of the diagonal elements (probabilities) taking place. In this case, one can show that the maximally mixed state given by

with being the Hilbert space dimension, is a stationary state of the dynamics. To prove this, we first note that the commutator with the Hamiltonian vanishes since the maximally mixed state is proportional to the identity. The dissipative part can be written as

which also vanishes in the case of all jump operators being Hermitian.

#### Dark states

[edit | edit source]On the other hand, there can be situations where the stationary state is a pure state, i.e., . To see how this comes about, we first introduce an effective non-Hermitian Hamiltonian,

Then, we can write the Lindblad master equation as

Then, to have a pure stationary state two conditions have to be fulfilled ^{[1]}:

- is an eigenstate of the effective Hamiltonian, i.e., .
- is an eigenstate of all jump operators. One can always re-define the jump operators such that one has for all without changing the master equation.

As an example, consider a three level atom interacting with the radiation field, consisting of two electronic ground states (e.g., different hyperfine states) and one electronically excited state that can decay into both ground states. In this case, the Hamiltonian in the rotating wave approximation is given by

where we have assumed that the Rabi frequencies on both transitions are identical. The jump operators take the form

and again we assume that both decay rates have the same value, . Then, we find that the antisymmetric combination is a stationary state of the quantum optical master equation. It is an eigenstate of the Hamiltonian with eigenvalue as the two excitation paths to the excited state interfere destructively. It is also annihilated by both jump operators as it does not contain any contribution from the state . As this state does not couple to the radiation field anymore, it is often called a "dark state".

Importantly, the dark state is also the only stationary state of the master equation. Any density matrix that has a contribution from will couple to the dark state as the quantity

is nonzero. Hence, the only remaining possiblity is the symmetric state . However, is not an eigenstate of the Hamiltonian as we find

Consequently, there are no other stationary density matrices besides . This phenomenon is known as "coherent population trapping" and can also be generalized to unequal decay rates ^{[2]}.

#### Dissipative preparation of a Bose condensate

[edit | edit source]We will now see how the concept of coherent population trapping can also be generalized to many-body systems. Consider a system of bosons in a lattice. Then, in the formalism of second quantization, we can assign a bosonic creation operator, to each site . The Hilbert space of a single state is a Fock space encoding the number of bosons located at the site, i.e.,

- ,

where is the vacuum state of zero particles that is destroyed by the annihilation operator, i.e., . The total Hilbert space of the system is a product space of all the Fock spaces, with the commutation relation for bosons being

With the help of a Fourier transform, we can also define creation and annihilation operators in momentum space,

where is the volume of the system.

In such a lattice system without interactions, the kinetic part of the Hamiltonian is given by the discretization of the Laplace operator,

where denotes the strength of the tunneling between adjacent lattice sites. Its ground state is a Bose-Einstein condensate (BEC) of all particles put into the mode with zero momentum, i.e.,

This state characterizes a distinct state of matter exhibiting long-range phase coherence of the form

The construction of suitable jump operators such that is a dark state can in general be split up into a product . Here, is the "interrogation part" that checks whether the system is already in the correct state. Hence, we need to satisfy and obtain a finite value for all other states. The "pump part" then reshuffles the populations. The precise choice of is usually not important, as long as it does not create a subspace in the Hilbert space, from which a decay into the dark state is no longer possible.

It might be tempting to choose to involve the projection onto the dark state, i.e.,

However, this jump operator contains terms that involve a product of operators involving sites separated by an arbitrary distance, so such a jump operator is unphysical. Nevertheless, it is possible to achieve the desired result using jump operators involving only adjacent sites and ^{[3]}, by making the choice

This jump operator can be understood as its interrogation part probing whether the phase is constant between adjacent sites. This makes sense because a BEC is a state of matter with long-range phase coherence. The pump part ensures that a phase fluctuation with nonzero will be projected onto a phase coherent state. To show that the BEC state is indeed a dark state of the dynamics, it is instructive to look at the Fourier transform of the jump operators ^{[3]},

From this one can see that the state with (i.e., the BEC) is the only dark state of the dynamics.

#### Frustration-free Hamiltonians

[edit | edit source]It is now very interesting, of course, to ask what kind of states we can prepare dissipatively. It should be noted that it is not possible to efficiently prepare arbitrary states because the preparation process might take a time that grows exponentially with the system size (e.g., in glassy systems) ^{[4]}. A particularly important class of models whose ground states can be prepared efficiently, are frustration-free Hamiltonians ^{[5]}. Such Hamiltonian are of the form

where the parts are projectors, i.e, , and they act only on a local subspace of the system. Additionally, while the parts do not have to commute in general, they have to commute within the ground state manifold, i.e.,

where is the ground state that can be found from the simultaneous minimization of all sub-parts, i.e., . If the parts commute with each other, the Hamiltonian belongs to a specific subclass of frustration-free Hamiltonians known as "stabilizer Hamiltonians" ^{[6]}. An example for a stabilizer Hamiltonian is Kitaev's toric code Hamiltonian ^{[7]},

where the "plaquette" operators and the "site" operators contain Pauli matrices representing four-body spin interactions. The ground state of the toric code is given by the constraints

for all plaquettes and sites, respectively, which allows to solve the system exactly. For periodic boundary conditions on a torus the stabilizers satisfy the relations and . for a system of atoms there are independent stabilizers. Consequently, the ground state manifold of
the system will be four-fold degenerate. Excitations of the toric code Hamiltonian are violations of the
stabilizer constraints and have energy gap of as every violation will affect two sites or plaquettes, respectively. Interestingly, these quasi-particles are neither bosonic nor fermionic, but will pick up minus sign when moved around each other ^{[7]}.

To show that it is possible to efficiently cool into the ground state, we consider a dynamical map of the form

where are probabilities and are projectors onto the high-energy subspace with ^{[5]}. This dynamics can be generated using local Hamiltonians and jump operators and any state in the ground state manifold, , is a stationary state of the dynamics satisfying . It can also be shown that this dissipative preparation process requires only a time that is polynomially increasing with the size of the system ^{[5]}.

#### Wave-function Monte-Carlo method

[edit | edit source]Finally, we want to conclude the course by taking a look at the practical problem of the numerical simulation of such dissipative quantum-many body systems. In general this is a very hard problem as the Hilbert space of a spin-1/2 many-body system grows exponentially like for spins. The size density matrix grows with the square of that, and including the cost for diagonalization of matrix to calculate the dynamics ^{[8]}, the whole problem requires steps, which becomes prohibitive for anything but very few spins. A way out of this difficulty is the reduction of the problem to work only on the level of wave functions and sample over many of them.

Consider an statistical ensemble of wave functions,

We now choose the size of the ensemble, to be much smaller than the Hilbert space dimension . For each member of the ensemble, , and at each timestep , we perform the following procedure, i ^{[9]}:

- Evolve the state under the non-Hermitian Hamiltonian
- This can be done approximately by

- For each jump operator , calculate the probability for a quantum jump to occur, which is given by
- Choose a uniformly distributed random number between 0 and 1.
- Apply the jump operator to the quantum state, , with chosen such that
- If is larger than the sum of all jump probabilities, no quantum jump occurs.

- Normalize the resulting state.

If we now average observables over all states of the ensemble, we will reproduce the dynamics of the observables under the full master equation according to ^{[9]}

Note that this method is only a first-order approximation to the dynamics and higher-order schemes exist ^{[10]}. Alternatively, one can approximate the dynamics in the same way by a stochastic version of the Schrödinger equation ^{[11]}.

## References

[edit | edit source]- ↑ Kraus 2008
- ↑ Cite error: Invalid
`<ref>`

tag; no text was provided for refs named`Gray1978`

- ↑
^{3.0}^{3.1}Cite error: Invalid`<ref>`

tag; no text was provided for refs named`Diehl2008`

- ↑ Cite error: Invalid
`<ref>`

tag; no text was provided for refs named`Kliesch2011`

- ↑
^{5.0}^{5.1}^{5.2}Cite error: Invalid`<ref>`

tag; no text was provided for refs named`Verstraete2009`

- ↑ Cite error: Invalid
`<ref>`

tag; no text was provided for refs named`Gottesman1996`

- ↑
^{7.0}^{7.1}Cite error: Invalid`<ref>`

tag; no text was provided for refs named`Kitaev2003`

- ↑ Cite error: Invalid
`<ref>`

tag; no text was provided for refs named`Press1992`

- ↑
^{9.0}^{9.1}Cite error: Invalid`<ref>`

tag; no text was provided for refs named`Molmer1993`

- ↑ Cite error: Invalid
`<ref>`

tag; no text was provided for refs named`Daley2014`

- ↑ Cite error: Invalid
`<ref>`

tag; no text was provided for refs named`Breuer2002`