summaryrefslogtreecommitdiff
path: root/xmds2/Genas_system
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2011-08-24 23:01:34 -0400
committerEugeniy Mikhailov <evgmik@gmail.com>2011-08-24 23:01:34 -0400
commita5e94f0c56f2388b4edbd8b73250ccb41e0bf674 (patch)
treeea74071878e8490e7c22574cbc28021b9b2cae98 /xmds2/Genas_system
parentdad3107ccbbf9a231d343b98c80f5e2f5fade86e (diff)
downloadNresonances-a5e94f0c56f2388b4edbd8b73250ccb41e0bf674.tar.gz
Nresonances-a5e94f0c56f2388b4edbd8b73250ccb41e0bf674.zip
Rho evolves according to Simon's equations
Diffstat (limited to 'xmds2/Genas_system')
-rw-r--r--xmds2/Genas_system/Genas_system.xmds10
1 files changed, 5 insertions, 5 deletions
diff --git a/xmds2/Genas_system/Genas_system.xmds b/xmds2/Genas_system/Genas_system.xmds
index cccee85..24cc2b4 100644
--- a/xmds2/Genas_system/Genas_system.xmds
+++ b/xmds2/Genas_system/Genas_system.xmds
@@ -178,14 +178,14 @@
// Equations of motions according to Simon's mathematica code
dr11_dt = gt - gt*r11 + (-Ey - Ex*i)*r12 + (-Ey + Ex*i)*r14 + (-Eyc + Exc*i)*r21 + g0*r22 + g0*r33 + (-Eyc - Exc*i)*r41 + g0*r44;
- dr12_dt = (Eyc - Exc*i)*r11 + (-g0/2. - gt - delta*i + i*OL)*r12 + i*rf*r13 + (-Eyc + Exc*i)*r22 + (-Eyc - Exc*i)*r42;
+ dr12_dt = (Eyc - Exc*i)*r11 + (-g0/2. - gt - delta*i + deltaL*i)*r12 + i*rf*r13 + (-Eyc + Exc*i)*r22 + (-Eyc - Exc*i)*r42;
dr13_dt = i*rf*r12 + (-g0/2. - gt - delta*i + i*delta0z(z))*r13 + i*rf*r14 + (-Eyc + Exc*i)*r23 + (-Eyc - Exc*i)*r43;
- dr14_dt = (Eyc + Exc*i)*r11 + i*rf*r13 + (-g0/2. - gt - delta*i - i*OL)*r14 + (-Eyc + Exc*i)*r24 + (-Eyc - Exc*i)*r44;
+ dr14_dt = (Eyc + Exc*i)*r11 + i*rf*r13 + (-g0/2. - gt - delta*i - deltaL*i)*r14 + (-Eyc + Exc*i)*r24 + (-Eyc - Exc*i)*r44;
dr22_dt = (Ey + Ex*i)*r12 + (Eyc - Exc*i)*r21 + (-g0 - gt)*r22 + i*rf*r23 - i*rf*r32;
- dr23_dt = (Ey + Ex*i)*r13 + i*rf*r22 + (-g0 - gt - i*OL + i*delta0z(z))*r23 + i*rf*r24 - i*rf*r33;
- dr24_dt = (Ey + Ex*i)*r14 + (Eyc + Exc*i)*r21 + i*rf*r23 + (-g0 - gt - 2*i*OL)*r24 - i*rf*r34;
+ dr23_dt = (Ey + Ex*i)*r13 + i*rf*r22 + (-g0 - gt - deltaL*i + i*delta0z(z))*r23 + i*rf*r24 - i*rf*r33;
+ dr24_dt = (Ey + Ex*i)*r14 + (Eyc + Exc*i)*r21 + i*rf*r23 + (-g0 - gt - 2*deltaL*i)*r24 - i*rf*r34;
dr33_dt = -(i*rf*r23) + i*rf*r32 + (-g0 - gt)*r33 + i*rf*r34 - i*rf*r43;
- dr34_dt = -(i*rf*r24) + (Eyc + Exc*i)*r31 + i*rf*r33 + (-g0 - gt - i*OL - i*delta0z(z))*r34 - i*rf*r44;
+ dr34_dt = -(i*rf*r24) + (Eyc + Exc*i)*r31 + i*rf*r33 - i*(deltaL - (g0 + gt)*i + delta0z(z))*r34 - i*rf*r44;
dr44_dt = (Ey - Ex*i)*r14 - i*rf*r34 + (Eyc + Exc*i)*r41 + i*rf*r43 + (-g0 - gt)*r44;
]]>
</operator>