CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace:
exercises:2019_conexs_newcastle:ex2

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
exercises:2019_conexs_newcastle:ex2 [2019/09/11 13:08] – [Making a XAS_TDP compatible ground state calculation] abussyexercises:2019_conexs_newcastle:ex2 [2020/08/21 10:15] (current) – external edit 127.0.0.1
Line 11: Line 11:
 You can find documentation for the XAS_TDP method [[ http://htmlpreview.github.io/?https://github.com/abussy/CONEXS/blob/master/xas_tdp_docs/CP2K_INPUT/FORCE_EVAL/DFT/XAS_TDP.html | here]] and refer to the {{:exercises:2019_conexs_newcastle:xas_tdp_conexs.pdf|slides}}. You can find documentation for the XAS_TDP method [[ http://htmlpreview.github.io/?https://github.com/abussy/CONEXS/blob/master/xas_tdp_docs/CP2K_INPUT/FORCE_EVAL/DFT/XAS_TDP.html | here]] and refer to the {{:exercises:2019_conexs_newcastle:xas_tdp_conexs.pdf|slides}}.
  
-**Remember** that the LR-TDDFT code is not yet included in the main CP2K code, hence you have to modify the `module load` line in the ''cp2k.sh'' file in order to load the correct version:+**Remember** that the LR-TDDFT code is not yet included in the main CP2K code, hence you have to modify the ''module load'' line in the ''cp2k.sh'' file in order to load the correct version:
  
 <code> <code>
Line 17: Line 17:
 </code> </code>
  
-To make sure you are using the proper version, type `module unload CP2K/6.1-foss-2017b` before launching any calculation.+To make sure you are using the proper version, type ''module unload CP2K/6.1-foss-2017b'' before launching any calculation.
 ===== Part 1: the standard workflow ===== ===== Part 1: the standard workflow =====
  
Line 252: Line 252:
 Once the calculation is over, open ''part1_check.out'' and look for ''XAS_TDP''. This will lead you to the relevant part of the output file, in which you can check the quality of the donor MOs. **Remember**: for the theory to apply correctly, the donor states need to be local in space and of 1s character. Look at the ''overlap'' and ''charge'' indicators, are we good to proceed ? Once the calculation is over, open ''part1_check.out'' and look for ''XAS_TDP''. This will lead you to the relevant part of the output file, in which you can check the quality of the donor MOs. **Remember**: for the theory to apply correctly, the donor states need to be local in space and of 1s character. Look at the ''overlap'' and ''charge'' indicators, are we good to proceed ?
  
-We are not, because the Aluminium atoms are equivalent under symmetry and the 2 lowest energy MOs will be arbitrary linear combinations of the 1s states. To get the needed donor state quality, uncomment the LOCALIZE keyword in the input file and rerun.+**No**, we are not, because the Aluminium atoms are equivalent under symmetry and the 2 lowest energy MOs will be arbitrary linear combinations of the 1s states. To get the needed donor state quality, uncomment the LOCALIZE keyword in the input file and rerun.
  
  
Line 366: Line 366:
 Note that again, the program reuses the previous ground state calculation. The spectral data can be found in the ''part1.spectrum'' file. The excitation energies are given in eV and the oscillator strength are within the dipole approximation (arbitrary units). There are two blocks of data, one for each Al 1s. Note that the values are practically the same since the Al atoms are indistinguishable.  Note that again, the program reuses the previous ground state calculation. The spectral data can be found in the ''part1.spectrum'' file. The excitation energies are given in eV and the oscillator strength are within the dipole approximation (arbitrary units). There are two blocks of data, one for each Al 1s. Note that the values are practically the same since the Al atoms are indistinguishable. 
  
-If you wish, you can plot the XANES spectrum using this [[https://raw.githubusercontent.com/abussy/CONEXS/master/plot_spectrum.py | python script]]. Make sure the script is located in the same directory as your ''part1.spectrum'' file and type: **Note: ** before using this script, you have to load a proper python installation by typing: ''module load part1.speAnaconda3/''+If you wish, you can plot the XANES spectrum using this [[https://raw.githubusercontent.com/abussy/CONEXS/master/plot_spectrum.py | python script]]. Make sure the script is located in the same directory as your ''part1.spectrum'' file and type: **Note: ** before using this script, you have to load a proper python installation by typing: ''module load Anaconda3/''
  
 <code> <code>
exercises/2019_conexs_newcastle/ex2.1568207290.txt.gz · Last modified: 2020/08/21 10:15 (external edit)