Photoelectron spectra with Qprop and t-SURFF

详细信息    查看全文

• 作者：Volker Mosert ; Dieter Bauer ; ^{dieter.bauer@uni-rostock.de" class="auth_mail" title="E-mail the corresponding author}
• 关键词：Photoelectron spectrum ; Time-dependent surface flux method ; Strong laser fields ; Time-dependent Schrö ; dinger equation
• 刊名：Computer Physics Communications
• 出版年：2016
• 出版时间：October 2016
• 年：2016
• 卷：207
• 期：Complete
• 页码：452-463
• 全文大小：966 K

文摘

Calculating strong-field, momentum-resolved photoelectron spectra (PES) from numerical solutions of the time-dependent Schrödinger equation (TDSE) is a very demanding task due to the large spatial excursions and drifts of electrons in intense laser fields. The time-dependent surface flux (t-SURFF) method for the calculation of PES [Tao and Scrinzi (2012)] allows to keep the numerical grid much smaller than the space over which the wavefunction would be spread at the end of the laser pulse. We present an implementation of the t-SURFF method in the well established TDSE-solver Qprop [Bauer and Koval (2006)]. Qprop efficiently propagates wavefunctions for single-active electron systems with spherically symmetric binding potentials in classical, linearly (along method=retrieve&_eid=1-s2.0-S0010465516301801&_mathId=si29.gif&_user=111111111&_pii=S0010465516301801&_rdoc=1&_issn=00104655&md5=2b95d8935b19779c629ad7a661491c1d" title="Click to view the MathML source">z$z$) or elliptically (in the method=retrieve&_eid=1-s2.0-S0010465516301801&_mathId=si30.gif&_user=111111111&_pii=S0010465516301801&_rdoc=1&_issn=00104655&md5=7ee8d8a32504c3ad46cb220bca73fca1" title="Click to view the MathML source">xy$xy$-plane) polarized laser fields in dipole approximation. Its combination with t-SURFF makes the simulation of PES feasible in cases where it is just too expensive to keep the entire wavefunction on the numerical grid, e.g., in the long-wavelength or long-pulse regime.

Program summary

Program title: Qprop

Catalogue identifier: ADXB_v2_0

Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADXB_v2_0.html

Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland

Licensing provisions: GNU General Public License, version 3

No. of lines in distributed program, including test data, etc.: 12458

No. of bytes in distributed program, including test data, etc.: 86258

Distribution format: at.gz

Programming language: C++.

Computer: x86_64.

Operating system: Linux.

RAM:   The memory requirements for calculating PES are determined by the maximum method=retrieve&_eid=1-s2.0-S0010465516301801&_mathId=si14.gif&_user=111111111&_pii=S0010465516301801&_rdoc=1&_issn=00104655&md5=da1dd11253d7f764d57a233b243670f0" title="Click to view the MathML source">ℓ$\ell$ in the spherical harmonics expansion of the wave function and the number of momentum (or energy) values for which the PES are to be calculated. The example with the largest memory demand (large-clubs) uses approximately 6GB of RAM. The size of the numerical representation of a wavefunction during propagation is modest for the examples included (53 MB for the large-club example).

Number of processors used: The evaluation of the PES can be distributed over up to method=retrieve&_eid=1-s2.0-S0010465516301801&_mathId=si32.gif&_user=111111111&_pii=S0010465516301801&_rdoc=1&_issn=00104655&md5=454e9659aaa68ce8ee8247a16e9c5214" title="Click to view the MathML source">N_k$N_k$ MPI processes (method=retrieve&_eid=1-s2.0-S0010465516301801&_mathId=si32.gif&_user=111111111&_pii=S0010465516301801&_rdoc=1&_issn=00104655&md5=454e9659aaa68ce8ee8247a16e9c5214" title="Click to view the MathML source">N_k$N_k$ is the number of momentum values).

Catalogue identifier of previous version: ADXB_v1_0

Journal reference of previous version: Comput. Phys. Comm. 174(2006)396

Classification: 2.5.

External routines: GNU Scientific Library, Open MPI (optional), BOOST (optional)

Does the new version supersede the previous version?: For TDDFT calculations the previous version should be used.

Nature of problem:

When atoms are ionized by intense laser fields electrons may escape with large momenta (especially when rescattering is involved). This translates to a rapidly spreading wavefunction in numerical simulations of these systems thus rendering the calculation of PES very costly for increasing wave lengths and peak intensities.

Solution method:

The TDSE is solved by propagating the wavefunction using a Crank–Nicolson propagator. The wavefunction is represented by an expansion in spherical harmonics. In order to reduce the requirements with respect to the grid size the t-SURFF method is used to calculate PES.

Reasons for new version:

Using the window operator method to calculate PES is increasingly costly with increasing ponderomotive energies and pulse durations. The new version of Qprop provides an implementation of the t-SURFF method which allows the use of much smaller numerical grids.

Summary of revisions:

An implementation of the t-SURFF method and examples for calculating PES are provided in the new release.

Restrictions:

The dipole approximation for the laser interaction has to be applicable. t-SURFF is only implemented for velocity gauge. Furthermore a finite cutoff for long range binding potentials has to be used in the implemented t-SURFF method.

Additional comments:

For additional information see www.qprop.de

Running time:

Depends strongly on the laser interaction studied. The examples given in this paper have run times from a few minutes to 12.5 hours.