I recommend against using the function callback form of the Hamiltonian - it's not really supported anymore and I couldn't make it work for your scenario for some reason.
Instead, you can use the following lines. Your code runs with these two modifications and seems to produce results like what I think you're expecting. Hope this helps.
wp = lambda t, args: (0.132)*exp(-((t-2.5*tau)**2)/tau**2)
output = mesolve([H0, [H1, wp]], psi0, tlist, c_op_list,
[gg0*gg0.dag(), uu0*uu0.dag(), xx0*xx0.dag(), yy0*yy0.dag(),
yy1*yy1.dag(), gg1*gg1.dag(), gg2*gg2.dag(), gg3*gg3.dag()])
As an aside, I don't think the expressions for P_0/P_1/P_2 in that paper are correct - I've been working on rigorously deriving exact solutions for the generic of problem of a system scattering light under coherent laser pulses. I don't see how my results can even reduce to what is written in this paper in some limit. Rather the quantities in that paper seem like heuristics. I'll be positing my manuscript soon on the arXiv if you're interested.
Regards,
Kevin
On Wednesday, October 4, 2017 at 3:43:28 AM UTC-7, Jitendra Verma wrote:
Dear Kelvin,
I have seen modifications done by Paul. But, problem I am facing is that I want to simulate the emission probability from the cavity leakage such as 2 photon emission probability from the state |g,2> by integrating with varying limit from 0 to t. Equation is attached in file. This is for same code attached in previous files.
Looking forward for your help.
Thank You,
Kind Regards,
Jitendra