CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace:
howto:rtp_field_xas

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
Last revisionBoth sides next revision
howto:rtp_field_xas [2023/10/11 09:51] – [CP2K input] glebretonhowto: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 {{:howto:rtp_field_xas.zip | this zip file}} along with the simulated data so that you do not need to run this calculation yourself to perform the analysis.+A longer version is available in the form of a jupyter notebook file in {{ :howto:rtp_field_xas.zip | this zip file}} along with the simulated data so that you do not need to run this calculation yourself to perform the analysis.
 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 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 on is loaded from the file RTP-RESTART.wfn, which is the ground state obtained after the SCF cycles. Note that you can define a wave-function that is not the ground state as long as the wave-function description (basis set, number of electrons...) is the same as the one used for the real-time propagation. The 'type' of reference is the default one, REFERENCE_TYPE SCF. All the molecular orbitals available in this reference wave-function will be used as a reference for the projection (REF_MO_INDEX -1 ), and stored in separate files (SUM_ON_ALL_REF .FALSE.). The projection will be computed at each time step.+The reference to projected to is loaded from the file RTP-RESTART.wfn, which is the ground state obtained after the SCF cycles. Note that you can define a wave-function that is not the ground state as long as the wave-function description (basis set, number of electrons...) is the same as the one used for the real-time propagation. The 'type' of the reference file is the default one, REFERENCE_TYPE SCF. All the molecular orbitals available in this reference wave-function will be used as a reference for the projection (REF_MO_INDEX -1 ), and stored in separate files (SUM_ON_ALL_REF .FALSE.). The projection will be computed at each time step.
  
-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 excited-states: 
  
 <code> <code>
Line 386: Line 386:
 </code> </code>
  
-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, see XAS_TDP%PRINT%RESTART_WFN section or the Linear Response input file in the {{:howto:rtp_field_xas.zip |tutorial archive}}. Running with the proper parameters, this XAS_TDP run produces a .wfn file per excited state: it is the $|\omega>$ presented before. Thus, using  REFERENCE_TYPE XAS_TDP automatically looks for the  $|\omega>$ saved in the .wfn file and performs the projection toward it: +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, see XAS_TDP%PRINT%RESTART_WFN section or the Linear Response input file in the {{:howto:rtp_field_xas.zip |tutorial archive}}. Running with the proper parameters, this XAS_TDP run saves one .wfn file per excited state: it is the $|\omega>$ presented before. Using REFERENCE_TYPE XAS_TDP, the projection uses automatically the $|\omega>$ saved in the .wfn file as the reference file: 
  
 $$ $$
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 to this state is:+There will be 7 projections per time step in this case: one for each time-dependent MO. The excited state population associated with this excited state is:
  
 $$ $$
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.  That is why ta third projection is requested: it uses the other excited state proposed by the XAD_TDP module which has the exact same frequency as the first one. +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 excited states orthogonal in the $xy$-plane.  That is why a third projection is requested which uses the other excited state proposed by the XAD_TDP module.
  
-Therefore, the excited state should be understood as the some over these two frequency equivalent state: +Therefore, the excited state should be understood as the sum over the two equivalent excited states:
  
 $$ $$
Line 407: Line 407:
 $$ $$
  
-Where $\omega'$ stands for the other equivalent excited state. +Where $\omega'$ stands for the other equivalent excited state, with $\omega=\omega'=529$ ev. 
  
  
howto/rtp_field_xas.txt · Last modified: 2024/02/24 10:02 by oschuett