Name | Modified | Size | Downloads / Week |
---|---|---|---|
Parent folder | |||
README | 2016-05-04 | 1.6 kB | |
awk_process_spectra | 2016-05-04 | 2.6 kB | |
awk_get_dipoles_cp2k | 2016-05-04 | 66 Bytes | |
awk_get_dipoles_gamess_us | 2016-05-04 | 127 Bytes | |
awk_jdx_converter | 2015-11-06 | 817 Bytes | |
awk_autocorr_function | 2015-11-06 | 434 Bytes | |
awk_calc_spectra | 2015-11-06 | 998 Bytes | |
Totals: 7 Items | 6.6 kB | 0 |
********************************************** * AWK scripts for * * IR abscorbtion spectra calculation * * on the basis of MD simulations * ********************************************** Usage (in general): awk -f <name_of_script> <name_of_input_file> [> <output_file] The following scripts should be used to extract IR absorbtion spectra from MD simulation trajectory in the following order: ?0) awk_get_dipoles* (extracts the dipole moments on each MD simulation step from GAMESS US (*=gamess_us) or CP2K (*=cp2l) output files); 1) awk_autocorr_function (calculates the autocorrelation function from the dipole MD trajectory, it dumps the results into the file autocorr.dat); 2) awk_calc_spectra (performes Fourier-transformation of autorcorr.dat file with Gaussian-like window function yelding raw absorbtion spectra); 3) awk_process_spectra (extracts the baseline from the raw spectra and multiplies it by the quantum form factor function yelding a normed absorbtion IR spectra from MD). For NIST experimental data use awk_jdx_converter (converts the JDX spectroscopic files into two-column file frequency vs. intensity). Versions log: YYYY-MM-DD 2015-11-06 - first version of scripts is now aviliable. 2016-05-04 - a method for the correction of finite time step frequency shift was added [D. Tikhonov // J. Chem. Phys. 144, 17 (2016), p. 174108], also -- new quantum form-factors (+ the Kubo one was corrected)