NWChem and COSMO: custom radii

To contact us Click HERE
There are two approaches described in the nwchem manual for using custom radii in COSMO calculations:

geometry  H  0 0 0  H  0 1 0  O  1 0 0endcosmo       radius 1.1              1.1              1.8end

and

cosmoendset cosmo:map custom.par

where custom.par looks like this:

H 1.1O 1.8

The downside to the first example is that it's a PITA to use -- you need to enter each vdw radius in the order you are listing the atoms in the geometry section. It means that for a 50 atom geometry you need to enter 50 values, even if all 50 atoms are the same element.

The downside with the second example is that you need to first create the run folder, put cosmo.par there and first then can you submit.

An easier approach is to create the custom.par on the fly using task shell:

task shell "echo -e 'H 1.1 \n O 1.8' > mycosmopars.par"cosmoendset cosmo:map mycosmopars.par