强场多电子原子的经典动力学研究 及辛算法计算

1. 强场多电子原子的经典动力学研究及辛算法计算强场多电子原子的经典动力学研究及辛算法计算 刘学深 吉林大学原子与分子物理研究所 2011.12

2. 研究工作： • 量子辛算法及其在激光原子物理中的应用 • 强场多电子原子的电离(辛算法) • 高次谐波及阿秒脉冲的产生 • 玻色-爱因斯坦凝聚理论研究 (辛算法)

3. 一个简单的例子: A2B 模型分子的经典轨迹 (C2v 对称性) Atom B moves along Z axis, two atoms A oscillate periodically. A2B模型分子的经典轨迹。 左：辛算法，右：Runge-Kutta方法 能量演化的比较

4. 物理背景 强场多电子原子的电离 • Non-sequential double ionization (NSDI) of atoms in intense laser fields is a typical topic in strong-field physics. 实验结果：

5. 非次序双电离 • “shake off” 模型 (Fittinghoff et al., Phys. Rev. Lett.69, 2642(1992)) one electron is thought to ionize very fast and then the second electron ionizes due to the sudden change of the binding potential and consequently ``shake-off'' the atom. • “correlated energy sharing” 模型 (Becker and Faisal, Phys. Rev. Lett.84, 3546(2000)) • “rescattering” 模型 (Corkum, Phys. Rev. Lett.71, 1994(1993)) one electron ionizes first and revisit the core, and then frees the second electron by collision.

6. 中红外长波长条件下的强场电离: Observation of the low-energy structure The LES structure can be characterized by a peak energy and high energy limit defined by the break in the slope. It is for the linearly polarized case. We can see the behavior of the ‘direct’ and ‘rescattered’ electrons. W. Quan, Z. Lin, M. Wu, H. Kang, H. Liu, X. Liu, J. Chen, J. Liu, X. T. He, S. G. Chen, H. Xiong, L. Guo, H. Xu, Y. Fu, Y. Cheng, and Z. Z. Xu, Classical aspects in above-threshold with a midinfrared strong laser field, Phys. Rev. Lett. 103, 093001 (2009). C. I. Blaga, F.Catoire, P. Colosimo, G. G. Paulus, H. G. Muller, P. Agostini and L. F. Dimauro, Strong-field photoionization revisited, Nature Physics, 5, 335, 2009.

7. Three-step model • P. B. Corkum, Phys. Rev. Lett. 71, 1994 (1993) • C. Winterfeldt, Colloquium: Reviews of Modern Physics, 2008, 80, 117-140 • M. Lewenstein, Ph. Balcou, M. Yu et al., Phys. Rev. A 49, 2117 (1994). The HHG and isolated attosecond pulse • Phys. Rev. A, 78(2008), 041402(R); 78(2008), 013401; 80(2009), 055404; 82(2010), 023402; 82(2010), 013411; 81(2010), 043420; (2011), revised. • J. Phys. B: At. Mol. Opt. Phys.,41(2008), 125602; 41 (2008), 015602; 42(2009), 125603.

8. 强场多电子原子的电离 • 理论方法 • 量子理论 : 求解含时薛定谔方程 • 经典理论: 将原子系统看作Coulomb势和外场作用下的质点系，系统运动由Newton方程或Hamilton正则方程描述，随机选取大量初态组成大量初值问题，求解方程的初值问题得到质点运动的轨迹，利用统计平均的方法，研究原子的动力学性质。 • 半经典理论： • 耦合干涉态（CCS）方法：The coupled coherent-state approach • 由 M. S. Child 等人发展并改进(Chemical Physics 347 (2008) 257–262) • CCS方法采用一系列的冻结高斯干涉态作为基轨道，从蒙特卡罗分布中取样，作为量子含时Schrödinger方程的基态。

9. 强场多电子原子的电离：经典理论方法(以He原子为例)强场多电子原子的电离：经典理论方法(以He原子为例) 牛顿方程 Legendre 变换 Hamilton正则方程 He的示意图

10. 为了求解Hamilton正则方程，需要随机的选取大量初值来描述真实的动力学过程。为了求解Hamilton正则方程，需要随机的选取大量初值来描述真实的动力学过程。 • 系综法随机选取初值 。(???) 采用辛算法计算的无外场时一维氦原子的经典轨迹随时间的演化。 未加激光场时分别采用辛算法与R-K法求得的系统总能量随时间的演化 • 对每一个初始状态，用辛算法数值求解激光场中He原子系统正则方程的初值问题即可得到电子的经典轨迹。 • 应用统计平均方法计算电离几率。

11. The classical ensemble method The Property of classical ensemble method it can be extended to more degrees of freedom such as multi-electron atoms and molecules; The calculation process is simple and time-saving, It can’t be used to consider the quantum effect.

12. 激光脉冲持续时间对He原子电离的影响 梯形激光脉冲中不同平台周期数（N）对应的“knee”结构。

13. 四种不同状态： • 初态 • 重碰撞态 • 电离态 • 自由运动态 P. J. Ho, R. Panfili, S. L. Hann and J.H. Eberly, Phys. Rev. Lett.94, 093002(2005). 在线性极化场中非次序双电离的能量分布。实线和虚线分别表示两电子的能量。激光场波长为390 nm. 首先，两个电子均处于束缚态上，并且均在平衡位置附近振荡，这是初态。接下来一个电子先电离，并且多次返回到核附近与第二个电子发生多重碰撞，这是重碰撞态。接下来，经过数次碰撞过程后，第二个电子终于电离（电离态）。在两个电子都电离后，它们都表现出类似于自由电子的振荡运动，相位相同或相反，这就是振动态。

14. Double ionization of Helium with classical ensemble simulations He2+/He+ 的比率随场强的变化是判断非次序双电离的一个重要标志。532 nm, 780 nm 和 1024 nm时，可以很清楚的看到“knee”结构. 表明随着场强的增大He的双电离从非次序向次序转变。在波长为248 nm 时观察不到“knee”结构，这说明在这一波长下次序双电离占主导地位。 Guo, Yu, and Liu, Phys. Lett. A, 372(2008), 5799

15. Nonsequential double ionization of Helium with elliptically polarized laser pulse 双电离与单电离之比随椭圆率变化 左上图：经典理论结果 左下图：量子计算结果 (Watson et al. JPB, 2000, 33, L103) 右 图：实验结果 (Dietrich P et al. PRA, 1994, 50, R3585) a）圆偏振光（b）线偏振光 Yu, Guo, Liu, Chinese Physics B, 19(2010), 023201.

16. Li原子一维模型 Li原子在外加激光场中运动，我们取质心坐标系，以Li原子核为坐标原点。假定激光场是线性极化的，核和电子的运动均沿着电场的方向。与He的计算相同,我们可得出电子的哈密顿量为 激光场与Li原子系统的相互作用势为

17. 系统哈密顿正则方程 (6) 数值求解哈密顿方程（6）的初值问题即可得到强场下Li原子在相空间中的经典运动轨迹。

18. 以上四种过程分别被定义为：存活、一次电离、二次电离和三次电离。以上四种过程分别被定义为：存活、一次电离、二次电离和三次电离。 电离判据 其中

19. Guo and Liu, Phys. Rev. A, 78(2008), 013401. 强场多电子原子的电离 Lithium ionization by intense laser fieldswith classical ensemble simulations

20. Lithium ionization by intense laser fieldswith classical ensemble simulations Guo and Liu, Phys. Rev. A, 78(2008), 013401.

21. Lithium ionization by intense laser fieldswith classical ensemble simulations 波长为780 nm 时典型的非次序双电离能量随时间的变化曲线. (初态，重碰撞阶段，电离态，自由振荡阶段) 波长为780 nm 时典型的非次序三电离能量随时间的变化曲线. (更加复杂)

22. 耦合干涉态理论在强场物理中的应用 Developed by M. S. Child et al (Chemical Physics 347 (2008) 257–262) • The advantage of CCS method • it can be extended to more degrees of freedom such as multielectron atoms and molecules; • the initial state can be chosen randomly; • the classical mechanics can be used to guide the motion of electrons.

23. 耦合干涉态理论 在坐标表象下，耦合干涉态由一系列的波包构成 • 经典轨迹由哈密顿量计算得出，这里我们考虑了量子修正 。 • 有序哈密顿量可以由将哈密顿量表示成关于产生和湮灭算符的形式而推导出， • 对氦原子，有:

24. 辛算法的应用 • 考虑量子修正后, 求解方程 • 上式为某种意义上的哈密顿正则方程,可用辛算法或R-K方法求解,从而得到电子的轨迹分布,以及经典作用量S.

25. 干涉态下的薛定谔方程的微分积分形式为 . • 其次，我们计算经典作用量 S=S(z, t) • 再次，我们计算次序哈密顿量的二阶变分 • 然后，计算 格点 上的系数积分矩阵 • 最后，用劈裂算符方法计算系数，从而得到波函数 首先，我们从上式中计算了 当给定一组初态 的值。

26. He的低能量结构(LES)区域 图2.1 低能量结构能由峰值能量和高能极限域构成。这是线性极化的情况。我们可以看到“直接”和“散射”电子的行为. J. Guo, X.S. Liu, S.-I Chu, Exploration of strong-field multiphoton double ionization, rescattering, and electron angular distribution of He atoms in intense long-wavelength laser fields: The Coupled Coherent-states Approach, Phys. Rev. A, 82(2010), 023402. L. F. DiMauro et al, Nature Physics, 5, 335(2009)

27. 低能量结构随激光参数的变化 图3 波长为2000 nm 时LES区域随激光强度的变化 图2 激光强度为100 TW/cm2时He原子的 LES区域随波长的变化

28. 非次序双电离机制 FIG. 4. (Color online) Position (a) and energy (b) distribution of a typical trajectory as a function of time for a linearly polarized laser field. The red (dashed) and black (solid) lines track the two electrons. The laser intensity is 200 TW/cm2.

29. 双势阱中玻色-爱因斯坦凝聚动力学性质研究------------几个计算问题双势阱中玻色-爱因斯坦凝聚动力学性质研究------------几个计算问题 单个玻色-爱因斯坦凝聚体隧穿动力学 GP方程为 标准双模近似 布居数差: 相位差: 系统的哈密顿量 这个方程的求解

31. 两个玻色-爱因斯坦凝聚体的自囚禁现象和对称性两个玻色-爱因斯坦凝聚体的自囚禁现象和对称性 标准双模模型: 哈密顿正则方程: 布居数差和相关相位 这个方程的求解

33. Spontaneous symmetry breaking of Bose-Fermi mixture in a two-dimensional double-well potential [25] Muruganandam P and Adhikari S K 2009 Comput. Phys. Commun. 1801888

35. Thank You