Differences

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

--- exercises:2018_ethz_mmm:infrared_2018 [2018/04/20 10:38] – dpasserone
+++ exercises:2018_ethz_mmm:infrared_2018 [2020/08/21 10:15] (current) – external edit 127.0.0.1
@@ Line 29: / Line 29: @@
 </code>
-<note warning>The ** .mol ** file for c6h6 is already there since the job is quite long (**C6H6-VIBRATIONS-1.ref.mol**)
+<note warning>The ** .mol ** file for c6h6 and methanol obtained with better precision (basis set) is already in the directory. (**C6H6-VIBRATIONS-1.ref.mol**) Using the command ** diff vib.c6h6.inp vib.c6h6.ref ** you can see the difference in the input parameters.
 </note>
@@ Line 68: / Line 68: @@
 ===== 2. Task: Computing vibrational spectra using DFTB molecular dynamics =====
-You will find a fortran program in the repository, called ** dipole_correlation.f90 **
+You will find a fortran program in the repository, called ** dipole_correlation.f90 ** . This is already compiled and the executable is dipole.x.
-Compile it (module load gcc; gfortran -o dipole.x dipole_correlation.f90 ). This program computes the correlation function of the (derivative of) the dipole moment and performs the Fourier transform.
+ This program computes the correlation function of the (derivative of) the dipole moment and performs the Fourier transform.
+Run ** cp2k ** with the ** md*.inp ** input files (for the two molecules). Note that the dipole moment and derivatives are extracted from simulation and saved in a file dip*traj (check the input). Run first 5000 steps, then edit the file dipole.in  and run ** ./dipole.x < dipole.in **.
-Run ** cp2k ** with the ** md*.inp ** input files (for the two molecules). Note that the dipole moment and derivatives are extracted from simulation and saved in a file dip*traj (check the input). Run first 5000 steps, then edit the file dipole.in  and run ** dipole.x < dipole.in **.
 This will generate the autocorrelation function of the dipole derivative (why?) and its Fourier transform (frequency domain).
-<note important>  - Check the result. Is it satisfactory with respect to the DFT static calculation and literature? Why?
+<note important>  - Check the result. **gnuplot** with the command **plot "file" u 1:2 w l** will help. Is it satisfactory with respect to the DFT static calculation and literature? Why?
   - Run 40000 more steps. Check the new results. Discuss what you obtained. Discuss the behavior of the autocorrelation in the time domain.
+  - ** WEB SITE [[https://www.cfa.harvard.edu/hitran/vibrational.html|Vibrational modes]] **
 </note>