Lecture 15: Unraveling Open System Quantum Dynamics

• Download English-US transcript (PDF)

The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free. To make a donation, or view additional materials from hundreds of MIT courses, visit MIT OpenCourseWare at ocw.mit.edu.

PROFESSOR: I apologize for the delay. But we are ready for multimedia show with three different applications. One of them are the clickers. I hope you all took a clicker. It's just two questions today, but really good questions.

[LAUGHTER]

So today's actually very nice topic. I think almost every lecture we have a nice topic, but today's especially nice. It's called Unraveling Open Quantum System. That may not tell you so much. Well unravelling means we want to sort of untangle it, We want to understand what is inside.

And so we want to develop a better understanding what is inside a density matrix. How can we think about the quantum system in an even more microscopic way than the density matrix tell us? And again, this finds applications in atomic physics, because one application of it is spontaneous emission optical block equations. On the other hand, actually, this class will also teach you something very general about quantum physics.

So today's lecture is this part. We'll have three topics. One is I want to provide a motivation for single quantum systems. And you should appreciate this. A single quantum system is not what Schrodinger and Heisenberg told us about. Quantum mechanics was formulated as an ensemble average. But today we want to talk about, is there something in addition? About single quantum systems.

Part two is-- addresses the derivation of the quantum Monte Carlo wave function method. And this is a formalism, which is really simulating in quite general terms, quantum physics by following-- using quantum Monte Carlo simulations-- trajectories of individual quantum systems. So here the individual quantum system-- even if you have one, even if you don't have one-- becomes a computational tool.

And eventually in part two I want to give you another example for unraveling open quantum system, which differ from the models for dephasing. And if you remember that this was on one of your homework assignment that's correct. I sort of want to here take your homework into the classroom discussion.

Well I think it's nice to start this lecture with the quote by Erwin Schrodinger. Where he said, "We never experiment with just one electron, or atom, or molecule. In thought experiments, we sometimes assume that we do; this inevitably entails ridiculous consequences." Well today we want to talk about those ridiculous consequences, because many people in our field trap, address, and observe single electrons and single atoms.

Now this is related to quantum jumps. Quantum jumps has even made it into the popular realm. You can find lots and lots of references on quantum jumping, but I warn you this is about quantum jumping into different worlds. It's about broadening your horizon. And those authors, those people actually, sell you something by claiming that people like Einstein and Stephen Hawking would support the idea that you can jump into different universe, experience another version of yourself-- probably a better version of yourself-- and this will inspire you to change your life in the world you live.

But-- So this is quantum jumping. I'm not talking about quantum jumping, I will only talk about quantum chunks. So what you see here is really a classic experiment. One of three experiments which were done almost simultaneously. There have always been races for, you know, the first observation of something. And there was the [? Hamburg ?] group by [? Toshak ?], the Seattle group by Hunt [INAUDIBLE], and the Boulder Group by Dave Wylen, who observed the phenomenon of quantum jumps for the first time.

Let me just tell you what it is, because it's exactly what you're going to describe. Here in this situation you have a system. It's a trapped barium ion with a p and s state. And there is very fast transition where you can rapidly scatter light, and get lots of fluorescence, at 493 nanometer. But what happened is, every once in a while, the particle has a branching ratio to decay to d state, and, well, d to s 2 units of angular momentum.

The transition here down to the ground state is very slow. It's slow because it's forbidden. It's even slower because it happens at longer wavelengths of 1,600 nanometer. And this is slow. And it has a lifetime of about 50 second.

So you-- If you think about it naively and that's probably the correct approach-- you would say the particles in the s state, it's like a light bulb. It scatters many, many photons, but then it branches to the d state. And for 50 seconds the light bulb is switched off, until the particle spontaneously goes down to the s state, and then the light bulb is switched on. And this is exactly what people observed.

But I mean I was-- This was '86. That was the time I was a graduate student, and I talked with theorists in the field, and respectable theorists told me that there was a controversy-- probably not among experimentalists, who use the naive approach-- but theorists who believed too much in what they calculate. Because if you solve it quantum mechanically, and look what you get. It's a steady state solution. Atom in steady state. Laser in steady state. And if you calculate the quantum mechanical expectation value, the expectation value is constant as a function of time.

So a simple quantum calculation would say the fluorescence is on average. Well just the time average and his constant, and his time independent. So there were some serious doubts whether you would observe this, or something which looks more like the ensemble average.

Well made all of the theorists soon realize you cannot calculate the steady state average. You have to calculate the correlation function in steady state. And the steady state correlation function has structure exactly at 50 seconds. So in a steady state system, if you're interested in fluctuations and temporal behavior, you can't just calculate quantities, you have to calculate correlation functions of quantities.

So anyway-- But then of course, if you now average over many, many realizations, you find exactly an exponential decay of 50 seconds. And this is sort of what you would get for an ensemble average.

So this is what quantum jumps are. This is what you should always imagine when we talk about quantum jumps in the rest of this lecture. But I also want to show you more recent example from Serge Haroche, in Paris, where they observed quantum jumps. Recording a single photon in a cavity. So what happens is there's a single photon in this cavity, and they can read out whether there is a photon or not by sending Rydberg atoms through. And a single photon is never observed, but it causes a phase shift, an AC-Stark shift, on the Rydberg atom. And then they can read out with a Ramsey interferometer the AC-Stark shift.

So every atom which passes through the cavity will tell you whether there is a photon or not. And-- And this is shown here. Blue means no, there is no photon. Red means there is a photon. So here you repeatedly figure out again there is no photon. And suddenly there is a quantum jump when the cavity has absorbed a photon, or when a thermal photon was released.

You may say, what happened here? Well no measurement is perfect. There is some noise. But by seeing that these are outliers, you clearly see that here is a photon in the cavity, and here where the red read out of the atom overwhelms the blue readout here you have no photon. And then as time goes by, you can observe the death and birth of micro-cavity photons. So these are now quantum jumps in the-- not in atomic population, but quantum jumps in the photon number in a cavity.

AUDIENCE:What's the time scale on that? Seconds?

PROFESSOR: Yeah, time is second. So what is shown here is-- I think what is shown is here are 300 milliseconds, and this is 2.5 second element.

AUDIENCE: What was the density of each line mean?

PROFESSOR: Each line is a Rydberg atom, which has been readout with a Ramsey interferometer. And--

AUDIENCE: So that means that they're leaving no, like between--

PROFESSOR: I think that means there was a Rydberg atom in state e, which means there is a photon. And blue means, the blue bar means, there was a Rydberg atom in state g. These are two highlighting states of the rubidium atom. Sort of each spike here is one readout of the microchannel plate. So these are detector clicks. You can say, yes. You can sort of say, blue is yes, and red is no. And these are the two clicks.

AUDIENCE: Yeah, but, like the separation between two clicks is not uniform.

PROFESSOR: No, because if you're an atomic beam the atoms come not equally spaced out. There's a randomness when the atom comes. I think it's posed only in statistics in this experiment. It's a real experiment. It's amazing. Other questions? Team One.

AUDIENCE: So one more question. So the, for example, the first bit of blue is all measurements for the same state of the cavity. I know blue means photon or no photon, but if it does mean photon, that means it's the same photon interacting with--

PROFESSOR: Yes.

AUDIENCE: OK.

PROFESSOR: Yes, because, I mean, this is something people know that the photon number in this cavity only changes over very long time scales. And you sort of see, kind of, exactly hear how-- You see here two quantum jumps. One from blue to red, and one from red to blue in this cavity. Yes?

AUDIENCE: I don't if this is maybe just my accent's not very good, but in the bottom left graph-- Right? Or on the bottom graph, in sort of the left corner, it looks like there's some overlap between the dark blue and the dark red. Yeah.

PROFESSOR: You mean here?

AUDIENCE: Yeah, like right there. Is that like sort of there is a photon in that [INAUDIBLE] simultaneously?

PROFESSOR: I think, to the best of my knowledge, I haven't really looked-- I've heard wonderful talks about it.You can go back and read the paper, but this probably means in the atom number in the cavity, and the algorithm they have used to analyze the data says that here the photon number was 0, jumped to 1, and jumped back again. So you-- If you see here, for instance, there is some overlap, but you would still say, since either there is a photon or not a photon, as long as you have a-- as you have a higher number of red clicks than blue clicks, you would say there is a photon.

And the little overlap, there is maybe a little bit of time where you don't know whether there is a photon or not. But this really-- This is real, because every measurement takes a certain amount of time. Since you have a few random clicks, if they are a few random clicks, you do not know yet if the system has jumped, you have to wait a little bit. But this is really also showing how our knowledge about the quantum system comes with time. And therefore, the wave function, or the statistical operator changes with time, because of wave function and statistical operator are simply a description of our knowledge about the quantum system.

AUDIENCE: [INAUDIBLE] is too dense

PROFESSOR: You could have a situation, for instance, you could prepare the cavity, that it has an equal probability foreseeable in one. And then what you would read out may just be completely random. So that may be possible, but in their experiment, I think, the nature here is that you have most of the time stable photons, because the cavity is cooled to very low temperature, and only occasionally does the thermal distribution, the Bose-Einstein distribution of photons in the cavity have a one. And what we observe are the thermal fluctuations between mainly 0 and occasionally a 1.

If you would make the cavity and little bit hotter, it would have maybe half of the time a 1 and half of the time 0. But that wouldn't mean that you have completely randomness in your results, because there is a thermal relaxation time over which a photon is created, and a photon goes back. The photon would not jump back and forth from the cavity to the walls of the cavity in microseconds. It would probably have an exchange time, and you could sort of measure the thermal correlation time of photons in a cavity.

It's a really good experiment. OK. So in that sense I think I showed you that what Schrodinger said is ridiculous is real, is the subject of current research, and since for quantum information processing in quantum computers, we need single quantum system is at the heart of quantum information processing.

So therefore, we want to now describe single quantum systems.

But when I do that, I-- come very -- It will become clear at the end of this unit, but I want to give you the warning right away. The way how I describe now a single quantum system may not be unique. There will be, there is actually, an infinite number of ways.

And I've already mentioned that to you, that there is an infinite number of ways to unravel a density matrix. I mentioned that you can write every density matrix as a probability of pure states. Weighted probability of being in pure states, but there's always different choices. You can get an infinite number of different combinations of pure states, and construct your density matrix. So, so this was sort of in at the level of the density matrix, but similarly-- and this will come out of this discussion today-- there is also an infinite number of ways to unravel the dynamics-- the dynamics involved in the evolution of an open quantum system.

So-- so with that I want to describe the quantum Monte Carlo wave function method. So in the quantum Monte Carlo wave function method, we perform-- We don't need an experiment. We don't need sophisticated experiments to observe a single quantum system. It's a theoretical method, but we perform thought, or gedanken experiments.

So our gedanken experiments is that we perform measurements. So we assume-- we have-- When we want to describe an ensemble of atoms, we assume we have a single atom. We assume it emits photon, and then we assume that the photon is detected. And the mom-- So the moment the photon is detected we know, for instance, that the atom is in the ground state.

And so our quantum Monte Carlo simulation follows now the trajectory of the single quantum system by saying, I toss a coin, there's a certain probability for detecting the photon. If I detect the photon, the system I know is in the ground state. If I don't detect the photon, it's not in the ground state, and such. And these probabilistic approach is just simulated in the quantum Monte Carlo sense.

So-- So what is in essence-- what is the essence of the quantum Monte Carlo method is that we allow a small time step delta t. It's a real time simulation. And there is a small probability that the system has decayed in that time. So we may now toss a coin-- use a random function generator-- and ask, has spontaneous emission taken place?

And then the two possibilities are yes, and we continue whatever the outcome, whatever the quantum state is. Or we say no, and then with the probability-- with a-- so with a probability, p, we say yes. With a probability 1 minus p we say no, and then we continue the simulation accordingly.

So this will become much, much clearer when I write down the formalism for you, but there is one important aspect, which I really want to discuss with you. And this is the essence of the method. Namely, how does the wave function change by your measurement?

So-- so this is the essential, the only conceptional part for the quantum Monte Carlo wave function. Afterwards I just show you. Once you have this, once we have cleared up this conceptional part, the rest is just a few equations. And we immediately show that the quantum Monte Carlo wave function method is equivalent to the optical block equation. But let's really feature and take some time for discussing the concept, namely, the change over wave function by a quantum measurement.

And maybe what this carries home to you is that the nature of an open system-- we have a quantum system, which interacts with the environment-- is, actually, that the environment permanently performs measurements. And you will actually see that the idea of performing measurements, formalized in the quantum Monte Carlo sense, gives us an equation for the density matrix, which is identical to the optical block equation.

So in other words, you can say I am re-deriving for you master equation now, but from a very, very different perspective. From the perspective of an environment which constantly performs measurement.

OK. So let's go through that. How-- What happens in a measurement? Well if our initial state is the ground state-- starting very in trivial ways-- well there is no measurement, because the particle will stay in the ground state forever.

OK. Now let's go to the next situation that our atom starts out in the excited state. Now in a time step delta t, we have two possibilities. A photon can be observed, or nor photon is observed. When the photon is observed, we know for sure that the system is now projected into the ground state. If no photon is observed, well the system is still in the excited state.

OK. So we've discussed what happens when the system is initially in the ground, and when the system is initially in the excited state. But now comes the following question. We prepare the system in a superposition of ground and excited state. And we want to discuss what happens when the photon is observed? What happens when no photon is observed? I think if you observe the photon it's pretty clear we have to assume after photon emission the system is in the ground state.

But if no photon is observed, I want to give you three choices, and you should address them with the clicker. One choice is that the particle stays in the same state. The second choice is that the particle stays in a superposition state, but the probability to be in the excited state decays. Also no photon has being emitted, or no photon has been observed. And the third possibility is something else.

So in other words, assume you're an observer, and you know with 100% detection probability-- let's not discuss technical issues-- you have completely surrounded your system with perfect detectors. And you know in the time step delta t, in the time t, no photon was observed. And now you should make a prediction what happens next to the quantum system. And you make the prediction by assuming that the system is described by a wave function or density matrix. And I'm asking you, which of the three choices correctly describes what the system will do afterwards? What the system will do next?

Is the question very clear? Good.

[TAPPING]

OK. All right. Pretty good. A lot of people when I taught the class before said it's A, the wave function is not changing, because no photon has been emitted. Now in-- For those people, but there only few among those, I would have now asked the next question.

If you have-- If you have a 50-50 superposition, of ground and excited state-- And let's say the spot time for spontaneous emission is 10 nanoseconds-- But now you wait one second, and after 1 second, no photon has been observed. What is your prediction? In what state is the system at this point? Is the system in the ground state? Or is the system in, still in, a superposition of ground and excited state? Well since we have to clicker, you can tell me. Is it-- Is it still a superposition state, or is the system in the ground state?

OK. So now you have to help me out, because the people who said B, if they would go back to here-- And most people said it is B. If you put in gamma 10 nanoseconds, the inverse of 10 nanosecond, and the time is 1 second, the excited state amplitude, the excited stated mixture, is 10 to the minus nothing. So-- OK. If you want to be a mathematician and say--

[LAUGHING]

--there is epsilon 10 to the minus 10, yes. Then you can say there is still a small admixture to the excited state. But by making it so extreme-- 10 nanoseconds versus a second-- at least in all practical terms, the system is in the ground state. And I think this should teach you what it means if you have a system which is in the superposition of ground and excited state, and after 1 second it has not emitted a photon. I mean--

Then you would say, it had its chance. It had plenty of chances to emit, but it decided it doesn't want to emit. And that would mean your 50/50 superposition of ground and excited state means the atom decided that it's in the ground state. You should also realize a 50% superposition of ground and excited state means in half, in 50% of the cases there will never be a photon emitted.

So in other words, what you should realize is if your 50-50% superposition of ground and excited state, that the quantum evolution of this system is-- after a long time the system is definitely positively, absolutely in the ground state. It's always in the ground state after a long time. In 50% of the cases with emission of a photon, in 50% of the cases without emission of a photon.

If you have a radioactive-- ensemble of radioactive nuclei-- and let's say half of the atoms are in the excited state and [INAUDIBLE] and half are in the ground state-- and you send the sample to your friend, and your friend waits many, many, many half times, common sense tells you that all the particles are now in the ground state. Also only half of them have decayed, because the other half were in the ground state to begin with. I mean this is what quantum physics tells us. This is what a 50/50 superposition is.

But quantitatively this is, of course, included in this result-- which magically that almost all of you got right-- namely there is a superposition of ground and excited state, but the excited state amplitude decays. And to the limit, in the limit that the time is much longer than the decay time, no matter what the coefficient initially was in front of the excited state, this coefficient has decayed to 0.

Any question about that? Let me just emphasize that, again in giving you also an example in the reference, a 0 measurement.

A 0 measurement modifies the wave function. You can also put it like this, what the wave function is is the best knowledge you would have about this system. The wave function is a way to predict what happens next, and it's the most precise way to predict it. More accurate predictions than the wave function always cannot be made, because of Heisenberg's uncertainty relation.

But what happens is, if you have a 50/50% superposition of ground and excited state, and the system has not emitted, then you would say initially my guess was half of the atoms were excited, half of the atoms are in the ground state. But if none of them has decayed for awhile, you would say I revise my estimate now. Now I have to assume that more are in the ground state. Because having not decayed means there is a higher probability that the system has actually-- is in the ground state to begin with. So this is how this way of dealing with a wave function automatically adjusts for the knowledge you have gained about the system.

Let me give you another example which carries home the same message, and I take it from a pedagogical paper from Dickey, from American Journal of Physics, 49, 926, 1981. And what he discussed is the following. You have a box and originally your atomic wave function is completely localized. Your Bose-Einstein condensate if you want in this box, or let's just assume one atom to be more precise, and it's completely delocalized.

But now you focus a laser beam into the system. And the laser beam would ionize your atom, and you could count the ion with 100% probability. And for a short moment, you flash on the laser and the result is yes, you count. You get a count.

In that moment you would say my wave function is more localized. I revise my estimate. Well maybe I shouldn't say ionized. If it's ionized you destroy the atom. Let's just say you observe fluorescence so the atom is still alive. But now you have to say that the atom is localized in that region.

However, if you don't observe anything, you would say, well now it's actually more probable that the atom is outside. Because if it had been in the laser beam, with a certain probability, I would have detected a photon. But the non-detection of a photon means that I revise my estimate. And I say it's much more probable that the atom is outside the laser beam than inside the laser beam.

So this here is exactly how a non-observation of anything, a 0 measurement, modifies your wave function. And the formalism to incorporate that is what we have discussed before.

Yes? Cody?

AUDIENCE: So your arguments here so far work the same way as if it's in a wave function, you have just an entire Pascal probability, and we're updating our understanding of probability. Like you haven't included anything about the relative phase between these two--

PROFESSOR: I want to do that now. I mean, I-- This is sort of just to address the basic concepts about what does a measurement mean, what does a 0 measurement mean. But now I want to write down everything for you in amplitudes. We want to do exactly the time evolution of the system. Nancy?

AUDIENCE: Is this, like, inconsistent with [INAUDIBLE] quantum equation? Like especially the [INAUDIBLE], because when you are saying that it is admixture of the ground and excited state, and we did not observe--

PROFESSOR: Then the superposition say that ground and excited state. We are not starting out in an energy heightened state, so therefore, we-- what we have is not a sharp value of the energy, but an expectation value. And also if you repeat the measurement many, many times, the average energy is conserved, and this is exactly what energy conservation in quantum physics says.

AUDIENCE: [INAUDIBLE]

PROFESSOR: Well if your particle, if your system, is not in an eigenstate of the Hamiltonian. It doesn't have a sharp energy. And that would mean that when you measure the energy now, you have fluctuations. Sometimes you measure higher, sometimes you measure lower. It's only when you're in eigenstate that you measure a sharp value.

And therefore you know from the beginning that you will measure an energy distribution. And sometimes you measure higher, and sometimes you measure lower than the average value. There's nothing wrong about it, but your question is a very good one. You should-- I mean those examples are really-- On the one hand they are trivial, but on the other hand, it's very profound what they tell us about quantum physics and how to apply conservation laws and such.

OK. I think-- Let's-- Let me now formalize exactly how this is done with all the bells and whistles. We assume we have an initial wave function, which is now an arbitrary superposition state of ground and excited state. And we have the environment, which, in this case, we assume is the 0 photon state. It's a vacuum.

And we want to exactly solve the Schrodinger equation for that system. But we restrict ourselves too much, to a very small time step. The small time step is smaller than, you know, anything else. Than the natural decay time. Maybe the inverse [INAUDIBLE] frequency if you drive the system. You're not doing it here, but I will later show how you can add a laser beam and drive the system. Or it should also be smaller than the inverse tuning.

So what is very important here is that we want to deal only with simple cases, namely one photon has been emitted or not photon. So you want to make sure that there is, at most, one spontaneous emission event during the time delta t. So what I'm writing down now is the total wave function of the system plus the environment.

And now we can evolve it. We can ask what happens? What is the wave function time later? This can be exactly done by time-dependent perturbation theory. And the result is that the system will be in a superposition of ground and excited state. And we want to calculate the coefficients alpha prime and beta prime.

And we still have the vacuum state. But then we have the possibility that a photon has been emitted. In that case we are in the ground state. And the direct product with the environment involves now the photon emitted into a certain direction with wave vector k, with the polarization epsilon. And we have coefficients beta, k epsilon, and we have to sum over all possibilities for the photon to be emitted.

So we call this-- So the wave function of the total system has now two parts. It's a wave function, which I call psi 0. This is a wave function which involves no photon in the environment. And the wave function psi 1 involves one photon in the environment.

So this is the [INAUDIBLE], which seems very natural, but it's also-- you could immediately prove that these are the only possibilities, how the system can evolve. And we can verify this by using time-dependent perturbation theory.

It's actually almost everyone of you has seen it in either 8.21 or in a more basic course on quantum physics. It's time-dependent perturbation theory for the emission of a photon. But usually when you see those treatments there reservoir, the environment, is not treated as explicitly as we do it here. The theory where this is treated exactly in this way, how we need it, is the Viegener Biscoff theory.

Which is nothing else in the perturbative approach, but it's really a perturbation theory. Not just for the atomic system-- how we sometimes present it in a simplified version-- it's perturbation theory for the total wave function of the complete system. Philmore?

AUDIENCE: So you mentioned the time interval [INAUDIBLE] is much smaller than a lot of things. But you were ignoring, I take it, the counter rotating terms? In which case we're still at a large enough time step that's average out, so to speak? I simply say this, because I expected to be an e with a photon term. From the g going up, so to speak.

PROFESSOR: Yes. Well-- The counter rotating term-- The first answer I wanted to give you, no this is exact perturbation theory. But I think I get myself in trouble if I would allow-- I would get myself in trouble if I would allow the time step to be extremely short, because then I'm in shorter than 1 over omega. Because during the time 1 over omega, we have-- you can say we have counter rotating terms. Or in other words, this particle in the ground state during a time 1 over omega, particle in the ground state can emit a photon and reabsorb it. These are where those weird diagrams with virtual states, which we discussed earlier in this course.

So I think I want to assume here, which I haven't done, that the time step is larger than omega to the minus 1. And as you know the counter rotating term has a detuning delta, which is 2 omega, so therefore it is excluded. Or in other words, when we talk about photons sent into the environment, we want to talk about real photons, and not virtual photons. Yes, good point.

OK. So if you do simple lowest order time dependent perturbation theory, in the Viegener Biscoff approach, you get an exact result for beta prime. You find that beta prime in this superposition is the original amplitude beta, but it has decayed with e to the minus gamma over 2dt. And since we are only interested in small time steps, we can do a linear expansion of that.

The probability, dp, that a photon has been emitted is the norm of this wave function psi 1. And this, using Viegener Biscoff perturbation theory, is gamma dt times the amplitude squared, that the system was excited to begin with. And because of the conservation of the norm, the wave function, the norm of the wave function psi 0-- which is a wave function without a photon being emitted-- is 1 minus dp.

So what I'm telling you here is that this is an exact result of time-dependent perturbation theory. Often when this-- when textbooks treat spontaneous emission, they're more interested in getting Fermi's gold rule at this rate. But the same approach-- if you just write it down-- tells you how the amplitude in the excited state evolves. So therefore, what we learn from perturbation theory, that alpha prime and beta prime, the wave, the coefficients of the wave function, without emission of a photon, this evolution of the wave function occurs with a non-emission Hamiltonian.

And this non-emission Hamiltonian is our Hamiltonian for the atomic system. But then we have to account for spontaneous decay, and this is done by this non-emission part. So occasionally you may have heard that people wave their hands and say, your system is described by a Hamiltonian, which has an imaginary part for decay. And this is sort of phenomenological.

This is definitely not the case, because, even if you have a Hamiltonian like this, a pure state would simply decay and remain a pure state. What we are doing here explicitly is we deal with probability in the correct way. Every time the system can branch out into different possibilities, two different dimensions of the density matrix-- one has probability delta p, one has probability 1 minus delta p-- and it is only the wave function associated with probability 1 minus delta p, which evolves with this Hamiltonian.

So what I'm telling you is the exact solution for the evolution of the total system, in terms of a density matrix for the atomic system. And this is exact. So that means we can now-- write down the-- Just get some extra space. Oops. How does it do it?

Nope So we have to insert new page. Yep. OK.

OK. So the procedure, how we implement this exact solution of perturbation theory, is the following. We have a time step delta t. We compute what is the probability that a photon will be emitted. Then we need a ran-- Then we need a random number generator. So we need a number epsilon, which is a random number chosen in the interval 0 and 1.

If this random number turns out to be smaller than delta p, then we continue our time evolution on the computer. That psi is now in the ground state, and maybe there is a laser beam which excites it again. And such I will add a few bells and whistles later. Otherwise psi, our wave function, is now the wave function which is the time evolution with a non-emission operator of the original wave function psi. And since we have a real wave function now with probability 1, we have to re-normalize the wave function by this denominator.

And then we execute the next time step. This means go to 2 and do the next time step. And then you have, so to speak, if you do it many, many times, you get a time direct trajectory of, so to speak, one experiment. And then you start again with your system in a wave function psi, and you accumulate a second experiment. And maybe you do it 10,000 times to get enough statistics. And then you sum up, you know, everything you want and calculate all the expectation values you want to know about your quantum system.

So the claim is that this method, called quantum Monte Carlo wave function method, is fully equivalent to the optical block equations. And I want to prove it to you by showing that, if I take a density matrix, which is an ensemble average, over many realizations of those quantum Monte Carlo trajectories. That then this density matrix follows the differential equation, which is your optical block equation.

And the proof is shown here. So what I told you is the density matrix after time delta t, has now two matrix elements. One with probability delta p, and this is-- photon has been emitted. System is in the ground state, with probability 1 minus delta p. We evolve the quantum state with a non-emission Hamiltonian.

And then all what is done in the next few steps we assume that delta t is small. We do a Taylor expansion of the exponent. We neglect quadratic terms in delta t, and then we come to this line. And if we write this as rho of t plus delta t, minus the original rho of t, we find that this follows differential equation, which is exactly the optical block equations.

So therefore what we sort of implemented as a form of doing many, many quantum measurements in the environment is a procedure, which is rigorously the same as the optical block equations, which we derived using a master equation approach.

OK. This method is very powerful for the following reasons. If you simulate a density matrix with, you know, an internal states, or external, internal states-- I should've said with in quantum states, you need n times n matrix elements, which can become quite a memory hog for your computer. The wave function, at any given moment, as only n components. So therefore there is a computational advantage in using a stochastic wave function approach over simulation of the density matrix.

The second thing is that a lot of people, especially experimentalists, like sort of this approach, because it reflects directly what they do in the experiment. And I think it's pretty clear that with this approach you can deal with a great variety of situations. Let me just mention two obvious extensions. One is polarization.

A photon is detected. That means your random number produced a number epsilon, which was smaller than delta p. At that point you can create a second random number, which determines if your pull-- if the photon has been detected with sigma plus, sigma minus, or pi polarization.

Or you can discuss recoil. If the photon is-- If the photon is detected, you use-- you throw on another random function generator, which determines what the k vector of the photon is. So what direction each part k the photon has taken, and this determines now what is the recoil, the change in momentum of your atom. So you see, kind of, you can start with an atom at 0 momentum. You can see it emits a photon, and then your computer always makes a choice based on the random number. And then you say, OK my photon has now received a recoil kick at 45 degrees.

Now your wave function of the photon is such and such, and if you add a laser to the Hamiltonian, the recoil may have kicked-- may have Doppler shifted the resonance, but everything is easily taken into account. So you see this quantum Monte Carlo wave function method can easily be extended to describe external degrees of freedom, multiple laser fields, and all that. Any questions?

OK. Let me now just generalize the thought, but I think you know already everything about it. We talked about the master equation. We talked about the most general master equation in the [INAUDIBLE] platform, where those operators L are the jump operators. And the prominent example for jump operator was the signal minus operator, which takes a particle from the excited state to the ground state. And this is the jump operator for spontaneous emission.

But you may have many jump operators. May be spontaneous emission of different photons with different polarization and such. Or the cavity may lose a photon, and then we have a jump operator for the cavity. We went through that.

Well that means now in an exact way that the non-emission Hamiltonian is the original Hamiltonian with an imaginary part, which comes because of those jump operators. And the general quantum Monte Carlo wave function procedure is that first you ask, has something happened? Has any jump happened? And this gives you the probability delta p.

And if-- If no jump has happened, you just do a time evolution with a non-emission Hamiltonian. But if a jump has happened, then the jump is now-- the wave function is projected by the jump operator. If you have a sigma minus operator, it takes a particle to the ground state. If you have several operators, a different jump operator may take your particle to another state. And you-- And you then have a branching ratio that you know first a jump has happened.

And then the question is, which jump has happened? And you just play the probabilities game. So that you see that the quantum Monte Carlo wave function method can be immediately be generalized.

So let me come back to the special case of spontaneous emission. In that case, we have only one jump operator, which is sigma minus. And the normalization is the square root of gamma. And that actually means that our quantum trajectory ground state, excited state, is very simple.

If you start with the system 100% in the excited state, it's very trivial what happens. Nothing happens. But if then the jump occurs, the particle is in the ground state. Then nothing happens then until infinite time.

You then have to start-- The next trajectory you start in the excited state. And this time, by chance of the random number generator, the jump happens later. Another time the jump happens earlier. And what you then have to do is, you have to average, over all of those realizations, what is the probability for the particle to be in the excited state. And this is now your estimator for the excited state diagonal matrix element of the density matrix. And if you-- If you've written your code correctly, you will find that you get a wonderful exponential decay exactly as you have expected.

Questions? Now-- Yes? Let-- Let me now talk about dephasing. What I've shown you so far is how we can unravel the density matrix in many microscopic realizations. And I've given you a specific example. We assume probability for spontaneous emission, and with that we propagate our wave function.

But now I come to, sort of, the nitty gritty details, or the dirty truth, that what we have assumed is by no means unique. And I want to explain it to you first by reminding you of a very, very nice homework assignment you have solved. And this was about we have optical block equations for the density matrix, and the dense -- and the solution of the optical block equation is that we have population damping with a damping time t1. And the off diagonal matrix elements are damped with a time t2.

Remember if you've only spontaneous emission t2 is 2 times t1. But we can have a lot of other processes, which can lead to much, much faster defacing time than the population, than the time of population changes. And in your homework, you have discussed three possibilities.

One is, well, spontaneous emission is energy loss. But you have discussed three different possibilities for t2. One-- Three different possibilities how dephasing, loss of coherence, and phase damping can happen. And in this homework assignment, you assumed either that there is a random phase, that the elastic collision which project onto the excited state, and then you got the most crazy and artificial model. That every time, randomly, the phase of the excited state flips form plus to minus.

So if I would implement that now, with a quantum Monte Carlo method, it would have the following affect. And I hope you enjoy, sort of, the graphical representation. You can really think about it, that this is what happens. That this is what inside-- what is inside the density matrices described by this process.

So let's assume we have a system, which would be in a superposition of ground and excited state. And that would mean that in the left frame the dipole moment would just rotate at the resonance frequency. It rotates, and it is a rotating dipole moment, which emits coherent light.

But now you assume that you have a random phase. So if you add random phases to it, maybe because you're fluctuating magnetic fields, then the line becomes sort of-- That doesn't look random, but you know what I mean. It becomes sort of jagged. The phase distribution is still sort of there, but there is a jitter on top of it. And as a result, the light emitted by this dipole is spectrally broadened

Well in your second model you assumed that-- with a certain randomness, with a certain time constant-- there's an elastic collision in the excited state, which projects the system into the excited state. At this moment there is no superposition state anymore, and the dipole moment is 0. So that's what you assume. Or when you have random phase flips, you assumed that suddenly this kind of data to data ministic sine function suddenly jumps corresponding to a minus sign in the excited state. And depending where you are in the cycle, it creates jumps at random places.

And the question of course is, if all those three processes-- and that's what you showed in your homework-- lead to the same density matrix, which one is correct? Which one is real-- is really realized in an experiment on a system which is described by this kind of damping, or by the optical block equation?

So you can say all or none. Well you can assume what you want, it doesn't make a difference. And the reason is actually profound. The reason-- The reason is profound in that sense that the way-- our goal was to have an open quantum system, which is the atomic system, interacting with a reservoir.

But the way, how we phrase the question, we are only interested in what the atomic system does. We didn't do extra measurements on the environment. The environment was just a dump for photons; a dump for energy; a dump for whatever we assumed in our dephasing mechanism.

So let me just take the example of putting photons into the environment. And this is described by a certain damping term, but there is an ambiguity. And this is shown here. If you evolve the system-- our, our atomic system, the density matrix evolves, but the environment sort of also evolves. And for instance here, the unitary time evolution has taken the excitation from the atom, and we have emitted a photon into the environment. And you remember we got the optical block equation by doing the partial trace here, and just focusing on the atomic part.

But now wait a moment. If we trace out the environment, we can trace it out in a different basis set. Remember, I showed you the quantum Monte Carlo wave function, we emit a photon. But we can, for instance, detect the photon in our gedanken experiment. We can detect it with linear polarization or with circular polarization.

If you have a situation where you have a-- an s state, which decays to a p state, sigma plus takes you to m equals plus 1, sigma minus takes you to m equals minus 1. So therefore, if photon is emitted in you quantum Monte Carlo procedure, you would say, now the atom mean is in the plus 1 state, or now the atom is in the minus 1 state.

But if you detect linearly polarized light, by just putting a polarizer here-- and this is a unitary transformation in front of your detector-- you would no longer project the atomic system on plus 1 at minus 1, you would project it on m equals 0 or whatever is connected to the linear polarization of the measurement.

So therefore, if you simply assume that something has been dumped in the environment, and can be used for measurement, there are many ways, many unitary transformations, what you can do to the information, to the energy, to the photons, which have been dumped into the environment. And each of them will lead to a very different trajectory in your quantum Monte Carlo system.

So therefore if you just dump the photons and not measure them, you have equal rights to assume that your quantum Monte Carlo wave function jumps to states which correspond to a circular basis, or jump to states which correspond to linear basis. And there's many, many possibilities. But each of those possibilities is 100% correct unravelling of the density matrix. And if you just do it right, assume that the measurement is done in a certain basis, and you're consistent with it, you will 100% correctly describe the time evolution of your atomic system.

Any questions? Philmore?

AUDIENCE: Again with the dt. I'm just curious that we take a [INAUDIBLE] very short, but why don't you run into something like the quantum Zeno effect? Where if every evolution is initially quadratic, like the [INAUDIBLE] frequency, you take very short dt-- short measuring system-- Shouldn't for a choice of sufficiently short dt you, you know the simulation would give you strange results, because of--

PROFESSOR: [SIGHS]

AUDIENCE: --some sort of quantum Zeno effect there?

PROFESSOR: That's a very deep question, Philmore. And I, I want to think about it more, but I think I've excluded that by saying that the time step, dt, is larger than the correlation time of the environment. So I'm doing some kind of Markov approximation with the environment, where s-- the quadratic part of the behavior is a coherent time evolution, for very short time steps a wave function evolves quadratically, but this part is a coherent evolution.

And that's related to the fact that an atom can emit a photon, but in the very first moment, before the vacuum has transported it away, it can take the photon back. And the result of that is that at very short times the exponential decay doesn't start out exponentially, it starts out a little bit flatter. We get the exponential decay, we get Fermi's golden rule, and we get optical blocks equations only if your time step is larger than the correlation time of the reservoir.

And the fact, you remember when we derived the master equation we had to say we do a step, which is sufficiently small for the dynamics of the atomic system, but sufficiently large, that we are not getting into any memory affects of the reservoir. And we've done the same assumption here. Other questions? Nicky?

AUDIENCE: Just confused. We are not always say the environment measures the system in range at time intervals, delta t, which are large compared to the correlation time? What's more [INAUDIBLE] the evolution of the atom. But now everyone isn't that defected the Markov approximation we've made when we derived the master equation? And I was just wondering for just-- [INAUDIBLE] did you derive it by making the master equation? Maybe we don't actually need-- Maybe we don't need the idea of constantly measuring. Maybe a physical interpretation of the Markov approximation.

PROFESSOR: Oh yes, exactly, Nicky. I still think you said it really very, very nicely. We derived a master equation just by assuming that the system, you know, takes a photon. And past the correlation time, the photon is taken for good by the environment. And with a Markov approximation, we formalized that there is no memory affect, the photon is not stored like in a cavity. The photon has disappeared. And by saying that this time is very short, this was a Markov approximation in the derivation of the master equation.

But that also means, if the photon has disappeared, has separated from the atomic system, we are now free to make a measurement. And of course it shouldn't matter whether we make the measurement or not. But what I was able to do today is, I was making now the assumption, let's assume we make a measurement. We measure all the photons which have been emitted, and the measurement is probabilistic. And we fold this probabilistic evolution into our quantum system. And what we obtained was exactly the same time evolution for the density matrix as we got from the master equation.

So in other words, this should you that interactions with an open quantum system, where irreversibly energy, photons, angular momentum, or whatever-- flows into the environment. Once it has flown into the environment, you can measure it. And you get all the stochastics from the measurement. But even if you don't measure it, you get exactly the same stochastics as if you had measured it. It's all the same. Jenny?

AUDIENCE: Now that I think about it, it seems sort of weird to me that we've-- that we're making these measurements, or lack of measurements, these interactions, at regular intervals. And not--

PROFESSOR: It doesn't really matter. All we-- All we have to make is-- You could actually make delta t a random variable. It wouldn't change anything. The only thing we have to make sure is we have to make the time interval short enough that we don't have, maybe, two photons emitted in that time. We just want to make sure that it's a simple probabilistic branching, yes or no. And that's the only requirement here.

AUDIENCE: Is this a formula or is it based on time? You say that would be the wave function evolving and then collapse. And that's-- This is what-- Where the measurement comes in. The equivalent, so now say, the equivalent is a wave function evolved in a collapse, when you make a measurement you say that after it evolves. So is it possible that it could also have-- like evolves some kind of formalism? Like this? Where there is an operator? Because still like even though the density matrix operator, maybe it is kind of like-- [INAUDIBLE] describes a state.

PROFESSOR: I haven't seen it, Mark, but I'm absolutely certain that you could describe the same physics in the Heisenberg picture, where the time evolution is with operators. Because what we have here is, we have a time step, delta t, where the system evolves as an isolated quantum system, and then we measure again. And so, I think if you would use operator equation, you would have an evolution of the operator with the same Hamiltonian, which is non-emission Hamiltonian. And then the measurement would do some form of projection.

I have to think about it, what the projection would be in terms of operators. But my gut feeling is, you can always take the transformation where you put the time dependence into the operators, and not in the wave function.

On the other hand, I have some little bit misgivings about that, because the quantum Monte Carlo wave function was really developed in order not to deal with matrices, not to deal with something which is dimension n times in, if n is the number of components of the wave function. It was specifically formulated to have the simpler description with n coefficients for the wave function. But conceptually, I think, it is evolution with a Hamiltonian measurement. Evolution with a Hamiltonian measurement. And I assume this could also be done with operators. Nancy?

AUDIENCE: Suddenly like [INAUDIBLE] but when we've-- these simulations are actively done, how important is the randomness? Like up to epsilon [INAUDIBLE] because no computer is actually emitting random numbers.

PROFESSOR: Well I think how random a random number has to be is really dealt with in computer science and mathematics.

AUDIENCE: But due to our reserves we always get this exponential behavior?

PROFESSOR: I don't know the answer. I mean here, conceptionally, it should be completely random number. And I think that algorithms, which even which produce pseudo random numbers, but if the pseudo is close enough, if there is a recurrence time of your random number-- if your random number-- your series of random number repeat itself after several billion random numbers, this gives you a small error bar. And if you don't-- If you'd only do a limited sampling, a limited number of time trajectories, I don't think it matters. Of course if you want to have ultimate precision, then also the precision of the random number may enter through the backdoor at some point.

AUDIENCE: [INAUDIBLE] lack of simulation of the photons emitted? But if we know the photons emitted, if we have that information, then we can decide [INAUDIBLE].

PROFESSOR: That's correct. If we-- If we would know-- If we would know more about the environment, if we would say we measured the photons in a certain basis, or when energy is dumped in the environment by a dephasing mechanism, we look at the environment. Then, of course, we would add extra information to it, and then certain quantum trajectories would reflect the extra knowledge we have. And of course, we would then describe everything in this basis.

On the other hand, what you should learn from this is once the information, once the photon has escaped, it doesn't matter for the time evolution of your atomic system, in which basis you measure the photon. So the evolution of the atomic system itself-- at least when you averaged-- is unaffected by the basis.

However if you make a coincidence measurement, that you say your atoms which has emitted, is still flying through your vacuum chamber, and now you figure out the photon is emitted as circularly polarized light. Then of course you know that this atom, which is still available, has emitted circularly polarized light. And then you have actually obtained extra information about your atom.

In other words, the photon in the atom was entangled, like a Bell pair. And if you now do a measurement on one part of the Bell pair, you know more about the other part. But this is just how quantum physics works. If you don't use this information for anything, then you could have-- you could have as well not measured the polarization. And it would have no effect on the atomic density matrix.

OK that's it for today. We have class at the usual time on Wednesday.