You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: README.md
+8-4Lines changed: 8 additions & 4 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -8,22 +8,26 @@ Please note this repo is for documentation. Usual developement will NOT be made
8
8
9
9
### Issues
10
10
2023.9.11<br>
11
-
(1) The missing scaling factor sqrt(2E) in the analysis script plot_part.m has been corrected.
11
+
(1) The missing scaling factor sqrt(2E) in the analysis script plot_part.m has been corrected.
12
12
13
-
(2) The format of output for the real-time projection onto bound states has been corrected.
13
+
(2) The format of output "proj_val" for the real-time projection onto bound states has been corrected. Now it follows the format in the user guide.
14
14
15
15
### Important
16
16
2023.9.11<br>
17
-
(1) The variable "lmd" in an output file differs a phase factor to the P_{lm}(E) defined in eq.(58):
17
+
(1) The variable "lmd1" in an output file actually differs a phase factor to the P_{lm}(E) as
18
18
19
19
LMD(k,m,l)= (-i)^l \exp(i\Delta_{l}(k)) P_{lm}(E)
20
20
21
-
This will not affect the use of eq.(65) but could cause some confusion when calculating physical quantities such as the time-delay.
21
+
Here P_{lm}(E) is defined in eq.(58). This feature will not affect the use of eq.(65), but could cause confusion when calculating physical quantities such as the time-delay.
22
22
23
23
(2) The gaussian envelope in the original release is really "gaussian", which means it is non-zero everywhere. The simulation starts and ends at 2FWHM far from its central peak. If the PCS method is applied for this case, some unphysical results would occur since the instantaneous A(t_f) is non-zero.
24
24
25
25
A modified gaussian envelope with its tail cut should be customized instead. One may refer to include/structure/field.h for its format.
26
26
27
+
(3) The projection onto field-free eigenstates could differ a sign between two runs using different box size parameters in the paper version code, since the LAPACK diagonalization routines randomly choose the sign of the eigenvectors. This may cause some confusion.
28
+
29
+
In the latest version we always fix the sign of eigen states to be positive near the origin.
30
+
27
31
### References
28
32
[1] Zhao-Han Zhang, Yang Li, Yi-Jia Mao, Feng He,
29
33
*QPC-TDSE: A parallel TDSE solver for atoms and small molecules in strong lasers*
0 commit comments