Register to reply 
2form and dissipative systems 
Share this thread: 
#1
Feb2211, 01:46 PM

P: 307

Suppose you have a dissipative system where
[tex] \dot{q}=p/m [/tex] [tex] \dot{p}=\gamma p k\sin(q) [/tex] So there isn't a Hamiltonian for this system and Louiville's theorem doesn't hold. But the equations of motion still give us a vector field on phase space and we can still take the Lie derivative of [tex] \omega = dp\wedge dq [/tex] along it. If I do that I get [tex] L_X \omega = \gamma dp\wedge dq [/tex] (let me know if you want to see the details, I'm omitting them here to keep things brief) where X is the dynamical vector field. I'm still new to the geometric understanding of mechanics and so the above 'looks' like a differential equation, but I'm not sure what to make of this? I suppose I'd like to see how to use this to interpret [tex] \omega [/tex] as a function of time, but I'm stuck. Any ideas are welcome. 


#2
Feb2211, 03:25 PM

P: 31

For this simple dissipation your Liouville operator is:
[tex]\boldsymbol{\mathcal{L}} = \frac{p}{m} \frac{\partial}{\partial q} +k \sin(q) \frac{\partial}{\partial p} + \gamma \frac{\partial}{\partial p} p[/tex] and so in addition to the dissipative force [tex]\gamma p \frac{\partial}{\partial p}[/tex], you have the probability preserving term [tex]\gamma[/tex]. Instead of Liouville's theorem you can apply its generalization: the method of characteristic curves. I am not sure (off the top of my head) what it means to do what you are doing as dissipative mechanics is not symplectic and you are working with a symplectic form. I would also be interested if anyone knows the dissipative generalization of symplectic evolution in classical mechanics. In quantum mechanics, the group structure is much simpler and this is known. 


#3
Feb2311, 07:57 AM

P: 307

Okay, well your Louiville operator corresponds to vector field I'm differentiating [tex]\omega[/tex] along. I think I understand your objection to possibly using [tex]\omega[/tex] however I can't see why there still wouldn't be such a two form on the phase space, Hamiltonian or no Hamiltonian? clearly dq and dp are one forms, why wouldn't their wedge exist?



#4
Feb2311, 08:41 AM

P: 31

2form and dissipative systems
What you have constructed exists  I don't question that. It's the interpretation I wonder about. The fact that it doesn't vanish along your dynamical vector field is probably telling you that no symplectic transformation is volume preserving here.
Moreover, if you omit the probability preserving [tex]\gamma[/tex] term from the correct Liouville operator, as to obtain a purely differential operator, then that differential part of the Liouville operator would actually be shrinking phasespace volume. That is probably what corresponds to the Liedifferential equation you have obtained [tex]L_X \omega = \gamma \omega[/tex], which looks like exponential decay. 


#5
Feb2311, 02:05 PM

P: 307

So do you know if we can 'do' anything with such a 'differential equation'?
I think we both agree that omega exists and gives area, but that area is time dependent. I'd like to compute that. So give some time t and vectors v and w, I'd like to know how omega evaluates on v and w at time t. Then also, how that area changes with time. My first thought is to evaluate the 'differential' equation on some vectors and perhaps obtain a differential equation for the area? 


#6
Feb2311, 05:41 PM

P: 31

I believe you can use the flow map to convert it to an ODE in time.



#7
Feb2411, 05:41 AM

P: 307

Do you have a favorite reference for learning about flow maps? I'm not necessarily familiar with this, or do you mean the flow obtained by equations of motion?
By the way, I appreciate the time you've taken to engage in this thread, thank you. 


#8
Feb2411, 07:12 PM

P: 31

I come from a dissipative quantum mechanics background, so my only good reference to what you are doing is a fine book called "The Geometry of Physics" by Theodore Frankel. But, in brief the flow map is simply the map generated by your dynamical vector field. The dynamical vector field belongs to the Lie algebra, whereas the flow map belongs to the Lie group. Perhaps you know the flow map by another name. In QM I would simply call it the propagator or transition matrix (for the density matrix).
And no problem about engagement, I am weak on exterior forms, so this is making me think about things I know in a different context. 


#9
Feb2511, 05:39 AM

P: 307

I have that book! I'll have to check it out this weekend. But it does sound like yes, the flow you're talking about is the 'solution' to the dynamical vector field.



#10
Feb2511, 10:36 AM

Sci Advisor
P: 1,943




#11
Feb2511, 12:35 PM

P: 31

In quantum mechanics, closed system evolution has unitary evolution with antihermitian generators ([tex]\imath \mathbf{H}[/tex]) whereas dissipative evolution has completelypositive evolution with LindbladGKS generators ([tex]\mathcal{L}[/tex]). (And here I am specifically talking about the algebraic generators. The timetranslation generators (a.k.a. the master equation) are only Lindblad in the Markovian regime. (For unitary dynamics this distinction is irrelevant even though the generators are also not always equivalent).). For classical physics the closed system evolution is symplectic, but what characterizes the dissipative evolution? Is there, or can there be, something like the LindbladGoriniKossakowskiSudarshan theorem for the generators of classical dissipative evolution? I don't do much classical physics, so I have never attempted to derive such a thing, but I should probably sit down and try one day. 


#12
Feb2511, 01:05 PM

Sci Advisor
P: 1,943

The most general Markov processes are a combination of diffusion processes and jump processes. See the differential ChapmanKolmogorov equation in Gardiner's Handbook of stochastic methods  (3.4.22) in the second edition. The derivation there is on the level of rigor of theoretical physics; but it is very likely that there is a fully rigorous version of this result in terms of measuretheoretic stochastic processes. Should you or someone else find a reference to such a mathematical presentation, I'd be interested. 


#13
Feb2511, 02:47 PM

P: 31

So really, I am asking about the underlying Lie group and algebraic generators of dissipative classical mechanics. Choi and LindbladGKS characterize both for quantum mechanics. That is what I compare to the Hamiltonian and its unitary evolution. What characterizes dissipative classical mechanics? That is what I would compare to the symplectic evolution if I knew it. It would be something that is not (phasespace) volume preserving, but still probability preserving. When the question was posed about applying symplectic thinking to a dissipative system, my reaction was one of hesitance, because the evolution is not symplectic. That was why I made the comment that started this. I wish I knew the better algebra to think about. The quantum correspondence in application between the algebraic and Markovian dynamics of the LindbladGKS theorem could be useful classically... if it still exists. Even then, I don't remember a Markovian classical theorem like LindbladGKS which says "the Liouvillian can only take the form ... in terms of [tex]x[/tex], [tex]\frac{\partial}{\partial x}[/tex], ...". LindbladGKS is extremely robust because it doesn't care what the model is, how you introduced the stochastic process, etc. I will look back over Gardiner though. Thanks for the recommendation. 


#14
Feb2511, 03:02 PM

Sci Advisor
P: 1,943




#15
Feb2511, 04:14 PM

P: 31

[tex]\boldsymbol{\rho}(t) = \boldsymbol{\mathcal{G}}(t,0) \boldsymbol{\rho}(0)[/tex] where [tex]\boldsymbol{\rho}(t)[/tex] is the density matrix of the reduced system at time [tex]t[/tex] and we have a nonunitary theory (likely with a traced out environment) such that we can consider any initial state [tex]\boldsymbol{\rho}(0)[/tex] and also we might later include ancillary degrees of freedom (e.g. external entanglement). Then with very few assumptions, [tex]\boldsymbol{\mathcal{G}}(t,0)[/tex] are completelypositive (CP) maps or semigroup elements. Then consider Choi's theorem on CP maps which immediately characterizes them: http://en.wikipedia.org/wiki/Choi%27..._positive_maps The (algebraic) generators of these semigroup elements are given by the LindbladGKS theorem, though it is usually only useful in the Markovian regime where the algebraic and dynamical generators are equivalent. (Otherwise, one can not extract very much from LindbladGKS.) If you refer back to the original papers of Lindblad and Gorini, Kossakowski and Sudarshan, you will see that they refer back to Choi. Choi's theorem describes the semigroup, LindbladGKS then describes the algebra which generates it. It's all a very beautiful structure. (I almost finished writing a review paper wherein I try to explain these lesser discussed (and applied) details.) I think classically, one would have to think akin to nonsymplectic flow maps and their generators. (Flow would then be a misnomer.) I am unfamiliar with what kind of structure these nonsymplectic maps would be constrained to have. Maybe it is something simple that every classical physicist knows. My knowledge of symplectic manifolds is very weak. 


#16
Feb2511, 06:16 PM

Sci Advisor
P: 1,943

The associated dynamical semigroups that preserve the trace (i.e., the sum of the entries) are the Markov chains, while the more general version you are after seem to be Markov chains with arbitrarily long memory. In infinite dimensions, and assuming appropriate topologies, you get in place of a Markov chain a combined jump&diffusion process, and presumably the more general version is an arbitrary stochastic process. But since you did not specify the quantum version precisely enough, I can't tell. Symplecticity never enters. The latter is present only when one specifies a Heisenberg algebra of distinguished operators, which provide a symplectic phase space structure. 


#17
Feb2511, 07:55 PM

P: 31

You will have to excuse my ignorance. Why are we mapping between nonnegative vectors? Naively I would imagine the starting point to be normpreserving positive linear maps between positive functions of the phasespace coordinates (i.e. density functions instead of density matrices). Is there some representation that I am missing?
I also need to think about classical correlations to ancillary degrees of freedom and whatever the analog to complete positivity would be, if any. With the Markov chain that doesn't seem to matter. 


#18
Feb2611, 03:36 AM

Sci Advisor
P: 1,943

I don't know whether the analogue of Choi's theorem has been proved rigorously. Certainly the corresponding Lindblad operators are used in quantum optics. This is too cryptic to make sense to me. 


Register to reply 
Related Discussions  
Liouville's theorem for Dissipative Systems  Classical Physics  0  
Entropy increase in dissipative systems  General Physics  1  
Can Liouville's theorem apply to dissipative systems  Classical Physics  6  
Dissipative vs. conservative  General Physics  1  
Rayleigh Dissipative Function  General Physics  2 