howto:rtp_field_xas
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionLast revisionBoth sides next revision | ||
howto:rtp_field_xas [2023/10/11 09:43] – [CP2K input] glebreton | howto:rtp_field_xas [2023/10/16 16:15] – oschuett | ||
---|---|---|---|
Line 5: | Line 5: | ||
In this tutorial, we will present a simulation of resonant X-Ray excitation of an isolated carbon monoxide in real-time using a time-dependent field. | In this tutorial, we will present a simulation of resonant X-Ray excitation of an isolated carbon monoxide in real-time using a time-dependent field. | ||
On this page, you will find an overview of the method, some equations, and the CP2K input file. | On this page, you will find an overview of the method, some equations, and the CP2K input file. | ||
- | A longer version is available in the form of a jupyter notebook file in {{: | + | A longer version is available in the form of a jupyter notebook file in {{ : |
This kind of calculation is not easy to grasp: do not hesitate to have a first look before diving into the equations and details! | This kind of calculation is not easy to grasp: do not hesitate to have a first look before diving into the equations and details! | ||
This tutorial is connected to this article REF where you can find complementary information. | This tutorial is connected to this article REF where you can find complementary information. | ||
Line 273: | Line 273: | ||
</ | </ | ||
- | We will not detail all the parameters and instead focus on the one related to the dynamics itself, the field, and the time-dependent projection part. | + | We will not detail all the parameters and instead focus on the one related to the real-time propagation, the field, and the time-dependent projection part. |
== Real-Time Propagation parameters == | == Real-Time Propagation parameters == | ||
Line 281: | Line 281: | ||
INITIAL_WFN SCF_WFN | INITIAL_WFN SCF_WFN | ||
</ | </ | ||
- | Therefore, before | + | Before |
We use a Molecular Orbital-based description of the wave function for the propagation along with the Arnoldi approach to compute the exponential: | We use a Molecular Orbital-based description of the wave function for the propagation along with the Arnoldi approach to compute the exponential: | ||
Line 288: | Line 288: | ||
</ | </ | ||
- | Note that for very large systems, the density-based method can be used for linear scaling (see the DENSITY_PROPAGATION keyword). | + | Note that for extensive |
For each time step, the AERTS algorithm is used with a convergence threshold defined by | For each time step, the AERTS algorithm is used with a convergence threshold defined by | ||
Line 313: | Line 313: | ||
The time-dependent electric field is defined for both gauges in the EFIELD section. It should be noted that one can define several EFIELD sections to apply several fields within the same simulation. | The time-dependent electric field is defined for both gauges in the EFIELD section. It should be noted that one can define several EFIELD sections to apply several fields within the same simulation. | ||
- | The field is defined by its envelope, its intensity | + | The field is defined by its envelope, its intensity, its polarization along the laboratory x, y, and z-axis, its wavelength, and the original phase. Several types of field envelopes can be used. Here we use a Gaussian one with a width of $\sigma=0.3073$ fs and centered at $T0=1.3190$ fs. The intensity used is 4.08E+13 W.cm$^{-2}$. Along with a carrying frequency of 529 ev (approximately 2.34374655955 nm), it should promote about $10^{-3}$ from the Oxygen 1s to the first available excited state. See {{: |
< | < | ||
Line 355: | Line 355: | ||
All the time-dependent Molecular Orbitals are projected (TD_MO_INDEX = -1) and these projections are stored separately (SUM_ON_ALL_TD .FALSE.). This calculation is spin-independent so one does not have to define the spin of the MO to project. | All the time-dependent Molecular Orbitals are projected (TD_MO_INDEX = -1) and these projections are stored separately (SUM_ON_ALL_TD .FALSE.). This calculation is spin-independent so one does not have to define the spin of the MO to project. | ||
- | The reference to projected | + | The reference to projected |
- | There are $N_e/2 = 7$ molecular orbitals for carbon monoxide. There are thus 7 time dependent MOs (the $i$) and 7 reference MO (the $j$), so that there will be $7x7=49$ projection per time step: | + | There are $N_e/2 = 7$ molecular orbitals for carbon monoxide. There are thus 7 time-dependent MOs (the $i$s) and 7 reference MO (the $j$s), so that there will be $7x7=49$ projection per time step: |
$$ | $$ | ||
Line 370: | Line 370: | ||
$$ | $$ | ||
- | Now, let us focus on the second and third projections which can be viewed as excited-state projections: | + | Now, let us focus on the second and third projections which can be viewed as projections toward |
< | < | ||
Line 386: | Line 386: | ||
</ | </ | ||
- | In this case, all the time-dependent MO are involved in the projection and stored separately. This time, the reference wave function is supposed to be from an XAS_TDP calculation, | + | In this case, all the time-dependent MO are involved in the projection and stored separately. This time, the reference wave function is supposed to be from an XAS_TDP calculation, |
$$ | $$ | ||
Line 393: | Line 393: | ||
Where $c_\omega^a$ is the $a^{\text{th}}$ atomic coefficient of the excited state found in the XAS_TDP module. | Where $c_\omega^a$ is the $a^{\text{th}}$ atomic coefficient of the excited state found in the XAS_TDP module. | ||
- | There will thus be 7 projections per time step in this case. The excited state population associated | + | There will be 7 projections per time step in this case: one for each time-dependent MO. The excited state population associated |
$$ | $$ | ||
Line 399: | Line 399: | ||
$$ | $$ | ||
- | The carbon monoxide molecule has a rotational symmetry along its CO bond: if one notes this axis $z$, then the $x$ and $y$-axis are equivalent by symmetry. It happens that the first available excited state for the Oxygen 1s is degenerate: there are two available states orthogonal in the $xy$-plane. | + | The carbon monoxide molecule has a rotational symmetry along its CO bond. If one notes this axis $z$, then the $x$ and $y$-axis are equivalent by symmetry. It happens that the first available excited state for the Oxygen 1s is degenerate: there are two available |
- | Therefore, the excited state should be understood as the some over these two frequency | + | Therefore, the excited state should be understood as the sum over the two equivalent |
$$ | $$ | ||
Line 407: | Line 407: | ||
$$ | $$ | ||
- | Where $\omega' | + | Where $\omega' |
howto/rtp_field_xas.txt · Last modified: 2024/02/24 10:02 by oschuett