1 |
chuckv |
4 |
!---------------------------------------------------------------------------- |
2 |
|
|
! Module Constants |
3 |
|
|
! |
4 |
|
|
!---------------------------------------------------------------------------- |
5 |
|
|
! Module supplies various constants to simulation. |
6 |
|
|
! |
7 |
|
|
! . Mathematical constants: |
8 |
|
|
! PI Value for pi from atan |
9 |
|
|
! . Defined Physical constants: |
10 |
|
|
! |
11 |
|
|
! amu_to_kg a.m.u. to kg |
12 |
|
|
! atm_to_pascals Atmospheres to pascals |
13 |
|
|
! ev_to_kj electron volts to kJ/mol |
14 |
|
|
! k_boltz Boltzmans constant |
15 |
|
|
! kcal_to_kj kcal/mol to kJ/mol |
16 |
|
|
! navogadro Avogadro's number |
17 |
|
|
! |
18 |
|
|
! . Derived Physical constants: |
19 |
|
|
! |
20 |
|
|
! R Gas constant |
21 |
|
|
! TO_DEGREES Converts radians to degrees |
22 |
|
|
! TO_HZ cm^-1 to fs^-1 |
23 |
|
|
! TO_RADIANS Degrees to radians |
24 |
|
|
! TO_WAVENUMBERS To cm^-1 |
25 |
|
|
! |
26 |
|
|
! . Constants used for simulation |
27 |
|
|
! |
28 |
|
|
! amua2fs2_to_K From amu (A fs^-1)^2 to Kelvin |
29 |
|
|
! amua3fs2_to_kjmol From amu (A fs^-1)^2 to kJ/mol |
30 |
|
|
! ms_to_afs From m/s to A/fs |
31 |
|
|
! kJmol_to_amua2fs2 From kJ/mol to amu (A/fs)^2 |
32 |
|
|
! amua2fs2_to_kcalmol From amu (A fs^-1)^2 to kcal/mol |
33 |
|
|
! kcalmol_to_amua2fs2 From Kcal/mol to amu (A/fs)^2 |
34 |
|
|
! . Notes |
35 |
|
|
! Fundamental units for this simulations are |
36 |
|
|
! time: femtoseconds |
37 |
|
|
! distance: angstroms |
38 |
|
|
! energy: kcal/mol |
39 |
|
|
! temperature: Kelvin |
40 |
|
|
!--------------------------------------------------------------------------------- |
41 |
|
|
|
42 |
|
|
module constants |
43 |
|
|
use definitions, ONLY : DP |
44 |
|
|
implicit none |
45 |
|
|
public |
46 |
|
|
|
47 |
|
|
!-------------------------------------------------------------------------------- |
48 |
|
|
! . Mathematical and Physical constants |
49 |
|
|
!-------------------------------------------------------------------------------- |
50 |
|
|
|
51 |
|
|
! . Atomic mass units to Kg. |
52 |
|
|
real( kind = DP ), parameter :: amu_to_kg = 1.6605402E-27_DP |
53 |
|
|
|
54 |
|
|
! . Atmospheres to Pascals. |
55 |
|
|
REAL ( KIND = DP ), PARAMETER :: ATM_TO_PASCALS = 1.013250E+5_DP |
56 |
|
|
|
57 |
|
|
! . The conversion factor from atomic units to Debyes. |
58 |
|
|
REAL ( KIND = DP ), PARAMETER :: AU_TO_DB = 2.54176568_DP |
59 |
|
|
|
60 |
|
|
! . The conversion factor from atomic units to electron-volts. |
61 |
|
|
REAL ( KIND = DP ), PARAMETER :: AU_TO_EV = 27.21_DP |
62 |
|
|
|
63 |
|
|
! . The conversion factor from atomic units to kJ/mole. |
64 |
|
|
REAL ( KIND = DP ), PARAMETER :: AU_TO_KJ = 27.21_DP * 4.184_DP * 23.061_DP |
65 |
|
|
|
66 |
|
|
! . The conversion factor from Bohrs to Angstroms. |
67 |
|
|
REAL ( KIND = DP ), PARAMETER :: BOHRS_TO_ANGSTROMS = 0.529177249_DP |
68 |
|
|
|
69 |
|
|
! . The speed of light. |
70 |
|
|
REAL ( KIND = DP ), PARAMETER :: C = 299792458.0_DP |
71 |
|
|
|
72 |
|
|
! . The conversion factor from electron-volts to kJ/mole. |
73 |
|
|
REAL ( KIND = DP ), PARAMETER :: EV_TO_KJ = 4.184_DP * 23.061_DP |
74 |
|
|
|
75 |
|
|
! . The permittivity of the vacuum (F m^-1). |
76 |
|
|
REAL ( KIND = DP ), PARAMETER :: E0 = 8.854187817E-12_DP |
77 |
|
|
|
78 |
|
|
! . Boltzmann's constant. |
79 |
|
|
REAL ( KIND = DP ), PARAMETER :: K_BOLTZ = 1.380658E-23_DP |
80 |
|
|
|
81 |
|
|
! . kcal mol^-1 to kJ mol^-1. |
82 |
|
|
REAL ( KIND = DP ), PARAMETER :: KCAL_TO_KJ = 4.184_DP |
83 |
|
|
|
84 |
|
|
! . Avogadro's number. |
85 |
|
|
REAL ( KIND = DP ), PARAMETER :: NAVOGADRO = 6.0221367E+23_DP |
86 |
|
|
|
87 |
|
|
! . Pi. Calculated to machine precision |
88 |
|
|
REAL ( KIND = DP ), PARAMETER :: PI = 3.141592653589793_DP |
89 |
|
|
!REAL ( KIND = DP ) :: pi_tmp = 4.0_DP*atan(1.0_DP) |
90 |
|
|
!REAL ( KIND = DP ), PARAMETER :: PI = pi_tmp |
91 |
|
|
REAL ( KIND = DP ), PARAMETER :: TWOPI = 2.0E0_DP * PI |
92 |
|
|
|
93 |
|
|
! . Planck's constant. |
94 |
|
|
REAL ( KIND = DP ), PARAMETER :: PLANCK = 6.6260755E-34_DP |
95 |
|
|
|
96 |
|
|
! . Elementary charge in coulombs. |
97 |
|
|
REAL ( KIND = DP ), PARAMETER :: TO_COULOMBS = 1.60217733E-19_DP |
98 |
|
|
|
99 |
|
|
!----------------------------------------------------------------------------------- |
100 |
|
|
! . Derived constants |
101 |
|
|
!----------------------------------------------------------------------------------- |
102 |
|
|
|
103 |
|
|
! . The conversion factor from Angstroms to Bohrs. |
104 |
|
|
REAL ( KIND = DP ), PARAMETER :: ANGSTROMS_TO_BOHRS = 1.0_DP / BOHRS_TO_ANGSTROMS |
105 |
|
|
|
106 |
|
|
! . The conversion factor from e^2/Angstroms to kJ mol^-1. |
107 |
|
|
REAL ( KIND = DP ), PARAMETER :: ELECT_CONST = ( 1.0E+7_DP * NAVOGADRO * TO_COULOMBS**2 ) / ( 4.0_DP * PI * E0 ) |
108 |
|
|
|
109 |
|
|
! . The conversion factor from atm. Angstroms^3 to kJ mol^-1. |
110 |
|
|
REAL ( KIND = DP ), PARAMETER :: PV_TO_KJ_MOLE = ATM_TO_PASCALS * NAVOGADRO * 1.0E-33_DP |
111 |
|
|
|
112 |
|
|
! . The gas constant (kJ mol^-1 K^-1). |
113 |
|
|
REAL ( KIND = DP ), PARAMETER :: R = K_BOLTZ * NAVOGADRO * 1.0E-3_DP |
114 |
|
|
|
115 |
|
|
! . The conversion factor from cm^-1 to ps^-1. |
116 |
|
|
REAL ( KIND = DP ), PARAMETER :: TO_HZ = C / 1.0E+10_DP |
117 |
|
|
|
118 |
|
|
! . Radians to degrees. |
119 |
|
|
REAL ( KIND = DP ), PARAMETER :: TO_DEGREES = 180.0_DP / PI |
120 |
|
|
|
121 |
|
|
! . The conversion factor from internal units to km/mole. |
122 |
|
|
REAL ( KIND = DP ), PARAMETER :: TO_KM_PER_MOLE = ( NAVOGADRO * TO_COULOMBS**2 ) / ( 12.0_DP * 1.0E+3_DP * E0 * C**2 * AMU_TO_KG ) |
123 |
|
|
|
124 |
|
|
! . Degrees to radians. |
125 |
|
|
REAL ( KIND = DP ), PARAMETER :: TO_RADIANS = PI / 180.0_DP |
126 |
|
|
|
127 |
|
|
! . The conversion factor from internal units to cm^-1. |
128 |
|
|
REAL ( KIND = DP ), PARAMETER :: TO_WAVENUMBERS = 1.0E+11_DP / ( 2.0_DP * PI * C ) |
129 |
|
|
|
130 |
|
|
!-------------------------------------------------------------------------------------------- |
131 |
|
|
! . Simulation constants |
132 |
|
|
|
133 |
|
|
! . The conversion factor from amu A^2 fs^-2 to kJ mol^-1 (equivalent to |
134 |
|
|
! . AMU_TO_KG * NAVOGADRO * 10^-3 / MS_TO_AFS^2 ). |
135 |
|
|
REAL ( KIND = DP ), PARAMETER :: AMUA2FS2_TO_KJMOL = 1.0E+4_DP |
136 |
|
|
|
137 |
|
|
! . The conversion factor from m s^-1 to A fs^-1. |
138 |
|
|
REAL ( KIND = DP ), PARAMETER :: MS_TO_AFS = 1.0E-5_DP |
139 |
|
|
|
140 |
|
|
! . The conversion factor from amu A^2 fs^-2 to Kelvin. |
141 |
|
|
REAL ( KIND = DP ), PARAMETER :: AMUA2FS2_TO_K = AMU_TO_KG / ( K_BOLTZ * MS_TO_AFS**2 ) |
142 |
|
|
|
143 |
|
|
! . The conversion factor from kJ mol^-1 to amu A^2 fs^-2. |
144 |
|
|
REAL ( KIND = DP ), PARAMETER :: KJMOL_TO_AMUA2FS2 = 1.0_DP / AMUA2FS2_TO_KJMOL |
145 |
|
|
|
146 |
|
|
! . The conversion factor from amu A^2 fs to Kcal/mol |
147 |
|
|
REAL( KIND = DP ), PARAMETER :: AMUA2FS2_TO_KCALMOL = AMUA2FS2_TO_KJMOL/KCAL_TO_KJ |
148 |
|
|
! . The conversion factor from kcal mol^-1 to amu A^2 fs^-2. |
149 |
|
|
REAL ( KIND = DP ), PARAMETER :: KCALMOL_TO_AMUA2FS2 = 1.0_DP / AMUA2FS2_TO_KCALMOL |
150 |
|
|
|
151 |
|
|
|
152 |
|
|
|
153 |
|
|
|
154 |
|
|
end module constants |