Dear QuTiP community,
I am currently simulating the HEOM dynamics of a three-site Frenkel-exciton model, where each site interacts with an independent bosonic bath characterized by an Ohmic spectral density. To model the environment, I used the QuTiP function env.approximate("prony", tlist_corr, separate=True, Nr=8, Ni=8) to fit the exact bath correlation function for use with HEOMSolver.
However, I’m encountering divergence in the population dynamics at steady state. I have also tested other fitting methods such as "cf" and "espira-I", and increased both the number of exponential terms and the hierarchy depth (up to 4). Unfortunately, going beyond this depth is computationally infeasible due to time constraints.
I would greatly appreciate it if someone could help verify whether there is an error in my simulation setup. For reference, I’ve attached the Python script I am using.
Thank you for your time and support.
Best regards,
Amit Kumar Upadhyay
Let me know if you'd like a shorter or more technical version.
--
You received this message because you are subscribed to the Google Groups "QuTiP: Quantum Toolbox in Python" group.
To unsubscribe from this group and stop receiving emails from it, send an email to qutip+un...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/qutip/a1e4bc9b-c69b-45e8-abc6-c9f59baef0c6n%40googlegroups.com.