Two-level Rabi Problem-numerical solution

Questions
, , , ,
1

Hello everyone! I want to numberically simulate the Rabi problem. It’s basically to solve the following equations:

屏幕截图_20230127_205446.png
屏幕截图_20230127_205446.png1750×275 20 KB

my code goes well when [CapitalDelta]= 10^-4,
[CapitalGamma] = 10^4, but breaks down when the order difference
goes higher. How to fix it?

Open in Wolfram Cloud ↗

askquestion_RWA.nb (38.1 KB)

2

Try:

   $Version
   (*"13.2.0 for Microsoft Windows (64-bit) (November 18, 2022)"*)


   timeRange =  2 Pi;
   \[CapitalDelta] = 10^-25; \[CapitalGamma] = 10^25;
   solution = 
     NDSolve[{I x'[
          t] == (Exp[-I \[CapitalDelta] t] + 
           Exp[-I \[CapitalGamma] t]) y[t], 
       I y'[t] == (Exp[I \[CapitalDelta] t] + 
           Exp[I \[CapitalGamma] t]) x[t], x[0] == 1, y[0] == 0}, {x, 
       y}, {t, 0, timeRange}, MaxSteps -> 50000, AccuracyGoal -> 4];
   Plot[{Evaluate[Norm[x[t] /. solution]^2], 
     Evaluate[Norm[y[t] /. solution]^2]}, {t, 0, timeRange}, 
    PlotTheme -> "Scientific", 
    FrameLabel -> {"t(s)", "|a(t)\!\(\*SuperscriptBox[\(|\), \(2\)]\)"}]

you get a nice plot.

Regards M.I.

3

It does. Thank you very much!