• Register
PhysicsOverflow is a next-generation academic platform for physicists and astronomers, including a community peer review system and a postgraduate-level discussion forum analogous to MathOverflow.

Welcome to PhysicsOverflow! PhysicsOverflow is an open platform for community peer review and graduate-level Physics discussion.

Please help promote PhysicsOverflow ads elsewhere if you like it.


New printer friendly PO pages!

Migration to Bielefeld University was successful!

Please vote for this year's PhysicsOverflow ads!

Please do help out in categorising submissions. Submit a paper to PhysicsOverflow!

... see more

Tools for paper authors

Submit paper
Claim Paper Authorship

Tools for SE users

Search User
Reclaim SE Account
Request Account Merger
Nativise imported posts
Claim post (deleted users)
Import SE post

Users whose questions have been imported from Physics Stack Exchange, Theoretical Physics Stack Exchange, or any other Stack Exchange site are kindly requested to reclaim their account and not to register as a new user.

Public \(\beta\) tools

Report a bug with a feature
Request a new functionality
404 page design
Send feedback


(propose a free ad)

Site Statistics

204 submissions , 162 unreviewed
5,030 questions , 2,184 unanswered
5,344 answers , 22,705 comments
1,470 users with positive rep
816 active unimported users
More ...

  How to handle nonmarkovian dynamics in open quantum system?

+ 3 like - 0 dislike
1229 views asked May 26, 2014 in Theoretical Physics by sciencestrider (15 points) [ no revision ]

Is it possible to clarify the question to be more specific? Do you mean a quantum system interacting with a bath which has a long relaxation time? There are many possible approaches.

If the bath has a long relaxation time, then the dynamics is non-Markovian i.e. the Markovian approximation is not valid. For a quantum system interacting a bath or a light pulse, I want to know a general answer without Markovian approximation.

4 Answers

+ 4 like - 0 dislike

Usually equations of motion containing a time convolution arise when a subsystem is projected out. A standard procedure is to use the Zwanzig projection operator method to obtain a generalised master equation. But it can also happen that such an equation of motion is directly given. This happens with the phenomenological Maxwell's equations for dispersive and absorptive systems. Then a "reverse" procedure can be followed. A new, auxiliary, field is introduced and the combined set of fields obey a unitary time evolution which can then be quantised. A posteriory the new field can here be identified as a polarisation field. References are

A. Tip: Phys. Rev. A 57, 1998, p. 4818.

B. Gralak and A. Tip: j.Math.Phys 51, 2010, p. 052902.

The method should also work for other linear cases with a time convolution term.

It is interesting to note that one can start from a unitary time evolution, use the Zwanzig method to obtain a time evolution with a convolution term and then use the reversed method to obtain a new unitary time evolution. The latter is not the same as the original one. It essentially picks up only stuff that is needed to describe the projected system. Note further that the method is exact. No approximations are made.

answered May 27, 2014 by Urgje [ no revision ]

Zwanzig projection operator method works when the environment (bath) stays in equilibrium. However, there are cases where the environment is can change due to the interaction with the system. For example, if we have a light pulse as input and let it interact with a qubit, the method does not work anymore.  

But the equation I displayed can be handled. In case we encounter non-linear terms the situation becomes more complicated and I have only made partial progress in that case.. There still exists a set of equations without time convolutions but I was not able to obtain a Lagrange-Hamilton formalism in that situation.

In that case the difficulty lies in the fact we may need infinite number of equations to know the dynamics. This is really a headache.

It is not that bad. For instance if you have in the macroscopic Maxwell's equations n non-linear susceptibilities then n additional equations are needed to convert to a set of equations without time convolutions.

+ 2 like - 0 dislike

The question is rather general. Are you thinking of an equation of motion of
the type

\(\begin{equation*} \partial _{t}f(t)=Af(t)+\int_{0}^{t}ds\rho (t-s)f(s)? \end{equation*}\)

Then I know of a procedure to convert it into a set of two equations without time-convolutions. Let me know. I used it to quantise the phenomenological Maxwell's equations for an absorptive system.

answered May 27, 2014 by Urgje [ no revision ]

This is the form of the master equation with a memory kernel. Can you convert it to time local equations?

+ 2 like - 0 dislike

You may either add a retarded memory term to your differential equation, giving an integro-differential equation of the kind appearing in the Zwanzig projection method before taking the Markovian limit, and use some phenomenological term for the memory. (This is treated, e.g., in the book by Honerkamp). 

Or you may extend the number of variables to include variables that act as memory, and then work with a Markovian master equation in the augmented set of variables. E.g., going from a hydrodynamic description with memory to a kinetic description.

answered May 28, 2014 by Arnold Neumaier (15,787 points) [ no revision ]
+ 2 like - 0 dislike

If you're interested in a time-local prescription, then you can use the projection operator technique to derive a time-convolutionless (TCL) master equation (see the book by Breuer & Petruccione for example). This technique works through a perturbative expansion in the system-bath coupling. At lowest order the TCL equation gives results of comparable (and often much better) accuracy than the standard Nakajima-Zwanzig equation.  If your system-bath coupling is in fact strong, then the perturbative expansion will fail. In this case it is often a good idea to redefine your system to include the slow variables of your bath, and treat only the coupling to the residual fast variables perturbatively. The polaron transformed spin-boson model and its generalisations furnish a well-known example of this idea.

It would really help if you describe the Hamiltonian properly, since the question is far too general for us to give you a useful answer at present. In essence you are currently asking for the general solution for an arbitrary dynamical system: almost all physical problems are non-Markovian if you choose your variables poorly.

answered May 29, 2014 by Mark Mitchison (270 points) [ revision history ]

Your answer

Please use answers only to (at least partly) answer questions. To comment, discuss, or ask for clarification, leave a comment instead.
To mask links under text, please type your text, highlight it, and click the "link" button. You can then enter your link URL.
Please consult the FAQ for as to how to format your post.
This is the answer box; if you want to write a comment instead, please use the 'add comment' button.
Live preview (may slow down editor)   Preview
Your name to display (optional):
Privacy: Your email address will only be used for sending these notifications.
Anti-spam verification:
If you are a human please identify the position of the character covered by the symbol $\varnothing$ in the following word:
Then drag the red bullet below over the corresponding character of our banner. When you drop it there, the bullet changes to green (on slow internet connections after a few seconds).
Please complete the anti-spam verification

user contributions licensed under cc by-sa 3.0 with attribution required

Your rights