Release66:1D-RISM
From NWChem
Line 1: | Line 1: | ||
__NOTITLE__ | __NOTITLE__ | ||
- | 1D-RISM module in NWChem provides description of solvated systems following one-dimensional reference interaction site of model of Chandler and Anderson. Similar to ab-initio density-functional theory, 1D-RISM can be thought of as an approach where discrete particle representation of solvent degrees of freedom is replaced by average density field. Unlike traditional continuum solvation model, this density based representation is inherently inhomogenous and incorporates specific molecular features of the solvent. In the current implementation, 1D-RISM is not directly coupled to QM calculations but presumed to be used as a post processing step after QM calculations which provide ESP point charges for a given solute geometry. | + | 1D-RISM module in NWChem provides description of solvated systems following one-dimensional reference interaction site of model of Chandler and Anderson. Similar to ab-initio density-functional theory, 1D-RISM can be thought of as an approach where discrete particle representation of solvent degrees of freedom is replaced by average density field. Unlike traditional continuum solvation model, this density based representation is inherently inhomogenous and incorporates specific molecular features of the solvent. In the current implementation, 1D-RISM is not directly coupled to QM calculations but presumed to be used as a post processing step after QM calculations which provide ESP point charges for a given solute geometry. |
Then parameters for 1D-RISM calculations are defined in the rism input block | Then parameters for 1D-RISM calculations are defined in the rism input block | ||
Line 41: | Line 41: | ||
The computed chemical potentials in both HNC and gaussian approximations are written in the output file. | The computed chemical potentials in both HNC and gaussian approximations are written in the output file. | ||
+ | |||
+ | Here is the complete example input file for solvated calculation of acetic acid. | ||
+ | |||
+ | echo | ||
+ | start rism | ||
+ | |||
+ | memory global 40 mb stack 23 mb heap 5 mb | ||
+ | |||
+ | rism | ||
+ | closure kh | ||
+ | temp 298 | ||
+ | vdw rule arithmetic parameters vdw.par | ||
+ | solute configuration solute2.data | ||
+ | end | ||
+ | |||
+ | task energy rism |
Revision as of 14:33, 19 October 2015
1D-RISM module in NWChem provides description of solvated systems following one-dimensional reference interaction site of model of Chandler and Anderson. Similar to ab-initio density-functional theory, 1D-RISM can be thought of as an approach where discrete particle representation of solvent degrees of freedom is replaced by average density field. Unlike traditional continuum solvation model, this density based representation is inherently inhomogenous and incorporates specific molecular features of the solvent. In the current implementation, 1D-RISM is not directly coupled to QM calculations but presumed to be used as a post processing step after QM calculations which provide ESP point charges for a given solute geometry.
Then parameters for 1D-RISM calculations are defined in the rism input block
rism solute configuration <filename> vdw [rule <arithmetic|geometric> ] parameters <filename> [temp <float default 298.15>] [closure <hnc|kh>] end
At this point energy task is supported, which is invoked using standard directive
task rism energy
- solute configuration - points to the file that contains information about the solute geometry. charges, and atom type mapping. The format is similar to xyz style with additional fields that specify charge and atom type. The atom type maps back to the vdw parameter file (see below). The example file is shown below
7 O1 -1.092111 0.733461 1.237573 -1.104415 O O2 0.758765 -0.201687 0.473908 -1.043019 O C1 -0.212954 1.568653 -0.833617 -0.474263 C1 C2 -0.174205 0.630432 0.357135 1.276672 C2 H1 0.360636 1.160405 -1.668859 0.102898 H H2 0.242419 2.521128 -0.531952 0.118979 H H3 -1.243967 1.772778 -1.139547 0.123148 H
- vdw - defines van der waals parameters
- rule - optional setting that specifies combination rule, defaults to "arithmetic"
- parameters - points to the file that contains vdw parameters for the system. Example file is shown below (note comments in the file)
#Van der Waals parameters file for RISM # type sigma(Angstrom) epsilon (kj/mol) C 0.3400E+01 0.3601E+00 H 0.2600E+01 0.0628E-00
- temp - defines temperature of the system with default value of 298.15
- closure - specifies choice of closure
Upon completion of the run, the resulting radial distribution functions are saved into rdf_out.data file.
The computed chemical potentials in both HNC and gaussian approximations are written in the output file.
Here is the complete example input file for solvated calculation of acetic acid.
echo start rism
memory global 40 mb stack 23 mb heap 5 mb
rism closure kh temp 298 vdw rule arithmetic parameters vdw.par solute configuration solute2.data end
task energy rism