diff options
| -rw-r--r-- | xmds2/lambda_no_dopler_no_z.example/lambda_no_dopler_no_z.xmds | 216 |
1 files changed, 110 insertions, 106 deletions
diff --git a/xmds2/lambda_no_dopler_no_z.example/lambda_no_dopler_no_z.xmds b/xmds2/lambda_no_dopler_no_z.example/lambda_no_dopler_no_z.xmds index fb3be35..c6eb904 100644 --- a/xmds2/lambda_no_dopler_no_z.example/lambda_no_dopler_no_z.xmds +++ b/xmds2/lambda_no_dopler_no_z.example/lambda_no_dopler_no_z.xmds @@ -1,109 +1,113 @@ <?xml version="1.0"?> <simulation xmds-version="2"> - <name>lambda_no_dopler_no_z</name> - - <author>Eugeniy Mikhailov</author> - <description> - Solving 3 level atom in lambda field configuration, - no field propagation included - at time = tset drive field is changed in step-like manner - </description> - - <features> - <benchmark /> - <bing /> - <fftw plan="patient" /> - <openmp /> - <auto_vectorise /> - <globals> - <![CDATA[ - const complex g=10; - const complex gbc=.001; - const double tset=5; - - complex Gab, Gca, Gcb; - complex Ep,Ed,Epc,Edc; - double delta_2photon, delta_d; - complex rba,rac,rbc; - complex drive_field0=10, probe_field0=0.01; - ]]> - </globals> - </features> - - <geometry> - <propagation_dimension> t </propagation_dimension> - <transverse_dimensions> - <dimension name="delta_2photon" lattice="128" domain="(-25, 25)" /> - </transverse_dimensions> - </geometry> - - <vector name="density_matrix" type="complex" dimensions="delta_2photon"> - <components>raa rbb rcc rab rca rcb</components> - <initialisation> - <![CDATA[ - raa = 0; - rbb = 1; - rcc = 0; - rab = 0; - rca = 0; - rcb= 0; - ]]> - </initialisation> - </vector> - - <sequence> - <integrate algorithm="ARK45" interval="10.0" tolerance="1e-7"> - <samples>128</samples> - <operators> - <integration_vectors>density_matrix</integration_vectors> - <![CDATA[ - delta_d=0; - - Ed=drive_field0; - Ep=probe_field0; - if (t<tset) { Ed=drive_field0;} else {Ed=drive_field0/2.0;} - Epc = conj(Ep); - Edc = conj(Ed); - - Gab=g+i*(delta_d+delta_2photon); - Gca=g-i*delta_d; - Gcb=gbc+i*delta_2photon; - - - rba=conj(rab); - rac=conj(rca); - rbc=conj(rcb); - - draa_dt = -i*Epc*rab+i*Ep*rba-i*Edc*rac+i*Ed*rca-2*g*raa; - drbb_dt = i*Epc*rab-i*Ep*rba+g*raa-gbc*rbb+gbc*rcc; - drcc_dt = i*Edc*rac-i*Ed*rca+g*raa-gbc*rcc+gbc*rbb; - drab_dt = -Gab*rab+i*Ep*(rbb-raa)+i*Ed*rcb; - drca_dt = -Gca*rca+i*Edc*(raa-rcc)-i*Epc*rcb; - drcb_dt = -Gcb*rcb-i*Ep*rca+i*Edc*rab; - ]]> - </operators> - </integrate> - </sequence> - - - - - <!-- The output to generate --> - <output format="ascii" filename="lambda_no_dopler_no_z.xsil"> - <group> - <sampling basis="delta_2photon" initial_sample="yes"> - <dependencies>density_matrix</dependencies> - <moments>rab_rlOut rab_imOut raaOut rbbOut rccOut</moments> - <![CDATA[ - rab_rlOut = rab.Re(); - rab_imOut = rab.Im(); - raaOut = raa.Re(); - rbbOut = rbb.Re(); - rccOut = rcc.Re(); - ]]> - </sampling> - </group> - </output> - -</simulation> + <name>lambda_no_dopler_no_z</name> + + <author>Eugeniy Mikhailov</author> + <description> + Solving 3 level atom in lambda field configuration, + no field propagation included + at time = tset drive field is changed in step-like manner + </description> + + <features> + <benchmark /> + <bing /> + <fftw plan="patient" /> + <openmp /> + <auto_vectorise /> + <globals> + <![CDATA[ + const complex g=10; + const complex gbc=.001; + const double tset=5; + + complex Gab, Gca, Gcb; + complex Ep,Ed,Epc,Edc; + double delta_2photon, delta_d; + complex rba,rac,rbc; + complex drive_field0=10, probe_field0=0.01; + ]]> + </globals> + </features> + + <geometry> + <propagation_dimension> t </propagation_dimension> + <transverse_dimensions> + <dimension name="delta_2photon" lattice="128" domain="(-25, 25)" /> + </transverse_dimensions> + </geometry> + + <vector name="density_matrix" type="complex" dimensions="delta_2photon"> + <components>raa rbb rcc rab rca rcb</components> + <initialisation> + <![CDATA[ + raa = 0; + rbb = 1; + rcc = 0; + rab = 0; + rca = 0; + rcb= 0; + ]]> + </initialisation> + </vector> + + <sequence> + <integrate algorithm="ARK45" interval="10.0" tolerance="1e-7"> + <samples>128</samples> + <operators> + <integration_vectors>density_matrix</integration_vectors> + <![CDATA[ + delta_d=0; + + Ed=drive_field0; + Ep=probe_field0; + if (t<tset) { Ed=drive_field0;} else {Ed=drive_field0/2.0;} + Epc = conj(Ep); + Edc = conj(Ed); + + Gab=g+i*(delta_d+delta_2photon); + Gca=g-i*delta_d; + Gcb=gbc+i*delta_2photon; + + + rba=conj(rab); + rac=conj(rca); + rbc=conj(rcb); + + draa_dt = -i*Epc*rab+i*Ep*rba-i*Edc*rac+i*Ed*rca-2*g*raa; + drbb_dt = i*Epc*rab-i*Ep*rba+g*raa-gbc*rbb+gbc*rcc; + drcc_dt = i*Edc*rac-i*Ed*rca+g*raa-gbc*rcc+gbc*rbb; + drab_dt = -Gab*rab+i*Ep*(rbb-raa)+i*Ed*rcb; + drca_dt = -Gca*rca+i*Edc*(raa-rcc)-i*Epc*rcb; + drcb_dt = -Gcb*rcb-i*Ep*rca+i*Edc*rab; + ]]> + </operators> + </integrate> + </sequence> + + + + + <!-- The output to generate --> + <output format="ascii" filename="lambda_no_dopler_no_z.xsil"> + <group> + <sampling basis="delta_2photon" initial_sample="yes"> + <dependencies>density_matrix</dependencies> + <moments>rab_rlOut rab_imOut raaOut rbbOut rccOut</moments> + <![CDATA[ + rab_rlOut = rab.Re(); + rab_imOut = rab.Im(); + raaOut = raa.Re(); + rbbOut = rbb.Re(); + rccOut = rcc.Re(); + ]]> + </sampling> + </group> + </output> + + </simulation> + +<!-- +vim: ts=2 sw=2 foldmethod=indent: +--> |