#acl +All:read apw185:read,write,delete Known:read All: <> [[AJMPublic/potentials/water-methane|Back to Water-Methane Potentials page]] = Water-Methane : Isotropic model with SAPT-5s charge arrangement on water and non-polarisable methane = Summary: 1. 5 site charge model on water and ISA charges on methane for electrostatics. No damping. 1. L1, isotropic undamped polarization model on water. 1. $C_{12}$, isotropic damped dispersion model. 1. Isotropic exchange-repulsion terms only. Units: Atomic Units == Citations == * A. J. Misquitta & A. J. Stone, //Ab initio atom–atom potentials using CamCASP: Theory// * Mas, E. M., Bukowski, R., Szalewicz, K., Groenenboom, G. C., Wormer, P. E., & van der Avoird, A. (2000). Water pair potential of near spectroscopic accuracy. I. Analysis of potential surface and virial coefficients. The Journal of Chemical Physics, 113(16), 6687-6701. * == Geometry and charges == Source for water: Mas et al. 2000 Source for methane: BS-ISA, A+DF $\zeta=0.1$ aug-cc-pVTZ main basis with aQZ/set2 aux basis. PBE0/AC. CamCASP 5.8. [[attachment:H2O_SAPT5s.mom]] {{{ O 0.000000 0.000000 0.000000 Type O Rank 0 0.2585046 H1 -1.453652 0.000000 -1.121687 Type HW Rank 0 0.5640504 H2 1.453652 0.000000 -1.121687 Type HW Rank 0 0.5640504 L1 0.000000 0.206721 -0.371792 Type L Rank 0 -0.6933027 L2 0.000000 -0.206721 -0.371792 Type L Rank 0 -0.6933027 }}} [[attachment:CH4-ISA.mom]] {{{ ! Based on DF-type : ISA C 0.00000000 0.00000000 0.00000000 Type C Rank 0 -0.413530 H1 0.00000000 2.07703578 0.00000000 Type H Rank 0 0.103598 H2 1.69589261 -0.69234525 0.97912406 Type H Rank 0 0.103383 H3 -1.69589261 -0.69234525 0.97912406 Type H Rank 0 0.103383 H4 0.00000000 -0.69234525 -1.95824810 Type H Rank 0 0.103166 }}} == Polarizability == Source: Averaged terms taken from WSM L3,aniso; then multiplied by 0.9 (see [[ajm/potentials/water/isowater5stet| water model]]. Main basis: d-aug-cc-pVTZ; Aux-basis: aug-cc-pVQZ. PBE0/AC with ALDA+CHF kernel. CamCASP 5.8. [[attachment:H2O_daTZ_ref_wt4_L3_static-0_9averaged.pol]] {{{ # Static polarizabilities O O 0.000 0.0000000 0.0000000 0.0000000 0.000 5.40801522 0.0000000 0.0000000 0.000 0.0000000 5.40801522 0.0000000 0.000 0.0000000 0.0000000 5.40801522 H1 H1 0.000 0.0000000 0.0000000 0.0000000 0.000 1.60741755 0.0000000 0.0000000 0.000 0.0000000 1.60741755 0.0000000 0.000 0.0000000 0.0000000 1.60741755 H2 H2 0.000 0.0000000 0.0000000 0.0000000 0.000 1.60741755 0.0000000 0.0000000 0.000 0.0000000 1.60741755 0.0000000 0.000 0.0000000 0.0000000 1.60741755 }}} == Potential & Damping == Dispersion coefficients taken from $C_{12}$ anisotropic model. Tang-Toennies damping. Single dispersion damping coefficient, no induction damping. [[attachment:CH4-H2O-iso-sapt5s-nonpol.pot]] {{{ C O rho alpha 00 00 0 6.326747 1.838297 END C HW rho alpha 00 00 0 4.314904 2.005177 END H O rho alpha 00 00 0 4.706678 1.895820 END H HW rho alpha 00 00 0 3.880222 1.719991 END ! Dispersion model C O Dispersion damping factor 1.88 End C HW Dispersion damping factor 1.88 End H O Dispersion damping factor 1.88 End H HW Dispersion damping factor 1.88 End C O C6 C7 C8 C9 C10 C11 C12 00 00 0 16.88311 0.0 497.2595 0.0 24469.86 0.0 668674.5 End C HW C6 C7 C8 C9 C10 C11 C12 00 00 0 4.753501 0.0 87.49391 0.0 4100.564 End H O C6 C7 C8 C9 C10 C11 C12 00 00 0 7.994712 0.0 85.56545 0.0 1658.857 End H HW C6 C7 C8 C9 C10 C11 C12 00 00 0 2.285648 End }}} == Sample Orient command file == You will need the above files present in the directory as they are '''included''' into this Orient file. Orient version: 4.7/4.8 [[attachment:CH4-H2O-iso-sapt5s-nonpol.ornt]] {{{ UNITS BOHR Parameters Sites 20 polarizable 20 S-functions 50000 Alphas 50000 Parameter-sets 50000 Pairs 100000 End Types C Z 6 H Z 1 O Z 8 HW Z 1 L Z 0 End Variables Rx 0.0 B Ry 0.0 B Rz 0.0 B alpha 0.0 D Nx 0.0 Ny 0.0 Nz 0.0 Index 0 End Molecule CH4 at 0.0 0.0 0.0 rotated by 0.0 about 1.0 0.0 0.0 #include CH4-ISA.mom End Molecule H2O at Rx Ry Rz ROTATED BY alpha ABOUT Nx Ny Nz #include H2O_SAPT5s.mom End Polarizabilities for H2O Read rank 1 #include H2O_daTZ_ref_wt4_L3_static-0_9averaged.pol End Limit rank 1 for O H1 H2 End Pairs #include CH4-H2O-iso-sapt5s-nonpol.pot End Units Bohr kJ/mol Switch Induce On Iterate On Options Induction Iterations 100 Convergence 1e-12 End Energy Table Format f12.6 Print es er ind disp total Variables Index Rx Ry Rz alpha Nx Ny Nz 3 0.00000000 -5.23172000 0.10250400 152.84713000 0.68617877 -0.68617877 0.24148992 4 0.00000000 -5.55870250 0.10891050 152.84713000 0.68617877 -0.68617877 0.24148992 5 0.00000000 -5.88568500 0.11531700 152.84713000 0.68617877 -0.68617877 0.24148992 6 0.00000000 -6.21266750 0.12172350 152.84713000 0.68617877 -0.68617877 0.24148992 7 0.00000000 -6.40885700 0.12556740 152.84713000 0.68617877 -0.68617877 0.24148992 8 0.00000000 -6.53965000 0.12813000 152.84713000 0.68617877 -0.68617877 0.24148992 9 0.00000000 -6.67044300 0.13069260 152.84713000 0.68617877 -0.68617877 0.24148992 10 0.00000000 -6.86663250 0.13453650 152.84713000 0.68617877 -0.68617877 0.24148992 11 0.00000000 -7.19361500 0.14094300 152.84713000 0.68617877 -0.68617877 0.24148992 12 0.00000000 -7.52059750 0.14734950 152.84713000 0.68617877 -0.68617877 0.24148992 13 0.00000000 -7.84758000 0.15375600 152.84713000 0.68617877 -0.68617877 0.24148992 14 0.00000000 -8.50154500 0.16656900 152.84713000 0.68617877 -0.68617877 0.24148992 15 0.00000000 -9.15551000 0.17938200 152.84713000 0.68617877 -0.68617877 0.24148992 End Finish }}} Output should end with [[attachment:end of output]] {{{ Rx Ry Rz alpha Nx Ny Nz Index es er ind disp total 3.00000 0.00000 -5.23172 0.10250 152.84713 0.68618 -0.68618 0.24149 -8.784589 61.964036 -0.222941 -27.479352 25.477154 4.00000 0.00000 -5.55870 0.10891 152.84713 0.68618 -0.68618 0.24149 -6.124523 35.196338 -0.106422 -19.249492 9.715901 5.00000 0.00000 -5.88568 0.11532 152.84713 0.68618 -0.68618 0.24149 -4.356535 19.985452 -0.053192 -13.540474 2.035252 6.00000 0.00000 -6.21267 0.12172 152.84713 0.68618 -0.68618 0.24149 -3.156715 11.344097 -0.027735 -9.579946 -1.420299 7.00000 0.00000 -6.40886 0.12557 152.84713 0.68618 -0.68618 0.24149 -2.623425 8.074759 -0.019121 -7.810300 -2.378088 8.00000 0.00000 -6.53965 0.12813 152.84713 0.68618 -0.68618 0.24149 -2.326471 6.436868 -0.015033 -6.826881 -2.731516 9.00000 0.00000 -6.67044 0.13069 152.84713 0.68618 -0.68618 0.24149 -2.068271 5.130966 -0.011885 -5.975322 -2.924512 10.00000 0.00000 -6.86663 0.13454 152.84713 0.68618 -0.68618 0.24149 -1.741479 3.651335 -0.008440 -4.905996 -3.004581 11.00000 0.00000 -7.19362 0.14094 152.84713 0.68618 -0.68618 0.24149 -1.322327 2.070761 -0.004894 -3.558695 -2.815154 12.00000 0.00000 -7.52060 0.14735 152.84713 0.68618 -0.68618 0.24149 -1.017308 1.174194 -0.002921 -2.607489 -2.453524 13.00000 0.00000 -7.84758 0.15376 152.84713 0.68618 -0.68618 0.24149 -0.792139 0.665747 -0.001791 -1.930770 -2.058952 14.00000 0.00000 -8.50155 0.16657 152.84713 0.68618 -0.68618 0.24149 -0.496130 0.214006 -0.000722 -1.093547 -1.376394 15.00000 0.00000 -9.15551 0.17938 152.84713 0.68618 -0.68618 0.24149 -0.322766 0.068808 -0.000316 -0.646459 -0.900733 }}}