Jump to content

Mod:Hunt Research Group/atomic radii

From ChemWiki

The atomic Radii

General comments

  • atomic radii are not well defined
  • most van der Waals radii refer back to an article by Bondi in 1964 doi:10.1021/j100785a001
  • UFF radii are defined in an article by Rappe in 1992 doi:10.1021/jp8111556
  • Batsnov compares radii in 2001 doi:10.1023/A:1011625728803
  • Truhlar has a redefinition based on the density in 2009 doi:10.1021/jp8111556
  • PCM and SMD use different radii, go to the end of this page to see more information
  • PCM do not have internal radii for some atoms, mainly TMs!

modify all radii

If you want to manually change ALL the VDW radii use the read command and specify an alpha, this changes all radii

# opt b3lyp/6-311g(d,p) int=ultrafine scf=conver=9
  geom=cartesian scrf(smd,solvent=water,read)

Title

-2 1
   6           0        0.000142    0.000521    0.000000
   8           0        0.007404   -1.299082    0.000000
   8           0        1.122188    0.655527    0.000000
   8           0       -1.129734    0.643021   -0.000000

alpha=1.15

without you will get:

 Atomic radii         : SMD-Coulomb.
 Cavity type          : VdW (van der Waals Surface) (Alpha=1.000).

using the above you will get:

 Atomic radii         : SMD-Coulomb.
 Cavity type          : Scaled VdW (van der Waals Surface) (Alpha=1.150).

modify radii definition

If you want to change the VDW radii definition used the read command and specify the type of radii, this changes all radii to the method identified

# opt b3lyp/6-311g(d,p) int=ultrafine scf=conver=9
  geom=cartesian scrf(pcm,solvent=water,read)

Title

-2 1
   6           0        0.000142    0.000521    0.000000
   8           0        0.007404   -1.299082    0.000000
   8           0        1.122188    0.655527    0.000000
   8           0       -1.129734    0.643021   -0.000000

Radii=Bondi

without you will get:

 Atomic radii         : UFF (Universal Force Field).
 Cavity type          : Scaled VdW (van der Waals Surface) (Alpha=1.100).

using the above you will get:

 Atomic radii         : Bondi.
 Cavity type          : Scaled VdW (van der Waals Surface) (Alpha=1.100).

modify specific radii

this currently does not work, the instructions on the gaussian web-site result in an error

ModifySph: Alters parameters for one or more spheres. The modified spheres 
can be indicated in the PCM input in lines following this keyword having
the following format:
   atom radius [alpha]
where atom is the atom number or element type.

include

ModifySph
    O   1.6000

and get

Using the following non-standard input for PCM:
 ModifySph
     O   1.6000
 ------------------------------------------------------------------------------
 The following wrong keyword(s) have been found:

 1.6000
 ------------------------------------------------------------------------------
 Error termination via Lnk1e in /Applications/g16/l301.exe at Sat Dec 23 16:41:53 2017.

The correct way to approach this is

  • use atom number and not the symbol (and not the atomic number!)
  • this option must be terminated with a blank line

Below the oxygen which happens to be the 4th atom in the z-matrix is given the radius 1.6 angstroms

ModifySph

    4   1.6000