Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Variational MPO approach for finding NESS #32

Open
4psireal2 opened this issue Jun 1, 2024 · 6 comments
Open

Variational MPO approach for finding NESS #32

4psireal2 opened this issue Jun 1, 2024 · 6 comments

Comments

@4psireal2
Copy link

Would there soon be an implementation of Variational MPO approach for finding NESS? See ref

@amartyabose
Copy link
Owner

This would be good to put in. However I don't have the time to immediately work on this. Would you be willing to contribute? I can help you out.

@4psireal2
Copy link
Author

I'd be! However, I've encountered great difficulties with finding NESS for my system of interest using this approach. I'm afraid that this approach is not numerically stable... I'll share my experience more and perhaps contribute once I get better convincing results!

@4psireal2
Copy link
Author

Also what's your ansatz for the density matrix in the Lindblad master equation?

@amartyabose
Copy link
Owner

I'd be! However, I've encountered great difficulties with finding NESS for my system of interest using this approach. I'm afraid that this approach is not numerically stable... I'll share my experience more and perhaps contribute once I get better convincing results!

Hmmm... By the way, a naive question: how are you defining steady state if these open quantum systems? Are you talking about the equilibrium reduced density matrix?

@amartyabose
Copy link
Owner

Also what's your ansatz for the density matrix in the Lindblad master equation?

In the GQME module? It's represented just as a matrix. See here: GQME documentation

@4psireal2
Copy link
Author

So I'm interested in non-equilibrium steady state (NESS) . Given dρ/dt = Lρ, steady state is a fixed point of the evolution dρ_s/dt = 0 and hence, ρ_s is a zero eigenvector of L. I hope that answers your question!

I see. I was thinking in a many-body setup, where the state is represented with a tensor network.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants