Example: COSMO: HCl
*******************

"Download Solv_HCl.run"

   #! /bin/sh


   # Computing solvent effects, with the COSMO model, is illustrated in the HCl
   # example.

   # After a non-solvent (reference) calculation, which is omitted here, two
   # solvent runs are presented, with somewhat different settings for a few input
   # parameters. The block key Solvation controls all solvent-related input.

   # All subkeys in the SOLVATION block are discussed in the User's Guide. Most of
   # them are rather technical and should not severely affect the outcome.
   # Physically relevant is the specification of the solute properties, by the
   # SOLVENT subkey: the dielectric constant and the effective radius of the
   # solvent molecule.

   # Note that a non-electrostatic terms as a function of surface area is included
   # in the COSMO calculation, by setting the values for CAV0 and CAV1 in the
   # subkey SOLVENT of the key SOLVATION. In ADF2010 one should explicitly include
   # such values for CAV0 and CAV1, otherwise this non-electrostatic term will be
   # taken to be zero, since the defaults have changed in ADF2010.

   # A rather strong impact on the computation times has the method of treating the
   # 'C-matrix'. There are 3 options (see the User's Guide): EXACT is the most
   # expensive, but presumably most accurate. POTENTIAL is the cheapest alternative
   # and is usually quite adequate. EXACT uses the exact charge density for the
   # Coulomb interaction between the molecular charge distribution and the point
   # charges (on the Van der Waals type molecular surface) which model the effects
   # of the solvent. The alternatives, notably 'POTENTIAL', use the fitted charge
   # density instead. Assuming that the fit is a fairly accurate approximation to
   # the exact charge density, the difference in outcome should be marginal.


   AMS_JOBNAME=Gas $AMSBIN/ams <<eor
   System
     symmetrize
     atoms
        H    0.000000     0.000000    0.000000
        Cl   1.304188     0.000000    0.000000
     end
   end

   Task SinglePoint

   Engine ADF
     title HCl(0) reference run (gas phase)
     eprint
       scf NoEigvec
     end
     basis
       type DZP
       CreateOutput Yes
     end
     noprint Bas EigSFO EKin SFO, frag, functions
   EndEngine
   eor



   AMS_JOBNAME=Exact $AMSBIN/ams <<eor
   System
     atoms
        H    0.000000     0.000000    0.000000    adf.R=1.18
        Cl   1.304188     0.000000    0.000000    adf.R=1.75
     end
   end

   Task SinglePoint

   Engine ADF
     title HCl(1) Solv-excl surface; Gauss-Seidel (old std options)
     eprint
       scf NoEigvec
     end
     noprint Bas EigSFO EKin SFO, frag, functions
     basis
       type DZP
     end
     solvation
       c-mat Exact
       charged Method=Gauss-Seidel
       disc SCale=0.01  LEGendre=10 TOLerance=1.0e-2
       div ND=4  min=0.5  Ofac=0.8
       scf Variational
       solv epsilon=78.8 radius=1.4 cav0=1.321 cav1=0.0067639
       surf delley
     end
     symmetry NOSYM
   EndEngine
   eor


   # In the second solvent run, another (technical) method is used for determining
   # the charge distribution on the cavity surface (conjugate-gradient versus
   # Gauss-Seidel in the previous calculation), and the POTENTIAL variety is used
   # for the C-matrix handling. The results show that it makes little difference in
   # outcome, but quite a bit in computation times.

   AMS_JOBNAME=Potential $AMSBIN/ams <<eor
   System
     symmetrize
     atoms
        H    0.000000     0.000000    0.000000    adf.R=1.18
        Cl   1.304188     0.000000    0.000000    adf.R=1.75
     end
   end

   Task SinglePoint

   Engine ADF
     title HCl(9) NoDisk and Cmatrix potential
     eprint
       scf NoEigvec
     end
     noprint Bas EigSFO EKin SFO, frag, functions
     basis
       type DZP
     end
     solvation
       c-mat POTENTIAL
       charged Method=conjugate-gradient
       div ND=4  min=0.5  Ofac=0.8
       scf Variational
       solv epsilon=78.8 radius=1.4 cav0=1.321 cav1=0.0067639
       surf delley
     end
   EndEngine
   eor
