|
1 | 1 | # -*- coding: utf-8 -*- |
| 2 | +# ruff: noqa: E741 |
2 | 3 |
|
3 | 4 | """ |
4 | | - Pair-state basis level diagram calculations |
5 | | -
|
6 | | - Calculates Rydberg spaghetti of level diagrams, as well as pertubative C6 |
7 | | - and similar properties. It also allows calculation of Foster resonances |
8 | | - tuned by DC electric fields. |
9 | | -
|
10 | | - Example: |
11 | | - Calculation of the Rydberg eigenstates in pair-state basis for Rubidium |
12 | | - in the vicinity of the |
13 | | - :math:`|60~S_{1/2}~m_j=1/2,~60~S_{1/2}~m_j=1/2\\rangle` state. Colour |
14 | | - highlights coupling strength from state :math:`6~P_{1/2}~m_j=1/2` with |
15 | | - :math:`\\pi` (:math:`q=0`) polarized light. |
16 | | - eigenstates:: |
17 | | -
|
18 | | - from arc import * |
19 | | - calc1 = PairStateInteractions(Rubidium(), 60, 0, 0.5, 60, 0, 0.5,0.5, 0.5) |
20 | | - calc1.defineBasis( 0., 0., 4, 5,10e9) |
21 | | - # optionally we can save now results of calculation for future use |
22 | | - saveCalculation(calc1,"mycalculation.pkl") |
23 | | - calculation1.diagonalise(linspace(1,10.0,30),250,progressOutput = True,drivingFromState=[6,1,0.5,0.5,0]) |
24 | | - calc1.plotLevelDiagram() |
25 | | - calc1.ax.set_xlim(1,10) |
26 | | - calc1.ax.set_ylim(-2,2) |
27 | | - calc1.showPlot() |
| 5 | +Pair-state basis level diagram calculations |
| 6 | +
|
| 7 | +Calculates Rydberg spaghetti of level diagrams, as well as pertubative C6 |
| 8 | +and similar properties. It also allows calculation of Foster resonances |
| 9 | +tuned by DC electric fields. |
| 10 | +
|
| 11 | +Example: |
| 12 | + Calculation of the Rydberg eigenstates in pair-state basis for Rubidium |
| 13 | + in the vicinity of the |
| 14 | + :math:`|60~S_{1/2}~m_j=1/2,~60~S_{1/2}~m_j=1/2\\rangle` state. Colour |
| 15 | + highlights coupling strength from state :math:`6~P_{1/2}~m_j=1/2` with |
| 16 | + :math:`\\pi` (:math:`q=0`) polarized light. |
| 17 | + eigenstates:: |
| 18 | +
|
| 19 | + from arc import * |
| 20 | + calc1 = PairStateInteractions(Rubidium(), 60, 0, 0.5, 60, 0, 0.5,0.5, 0.5) |
| 21 | + calc1.defineBasis( 0., 0., 4, 5,10e9) |
| 22 | + # optionally we can save now results of calculation for future use |
| 23 | + saveCalculation(calc1,"mycalculation.pkl") |
| 24 | + calculation1.diagonalise(linspace(1,10.0,30),250,progressOutput = True,drivingFromState=[6,1,0.5,0.5,0]) |
| 25 | + calc1.plotLevelDiagram() |
| 26 | + calc1.ax.set_xlim(1,10) |
| 27 | + calc1.ax.set_ylim(-2,2) |
| 28 | + calc1.showPlot() |
28 | 29 |
|
29 | 30 | """ |
30 | 31 |
|
|
0 commit comments