next up previous contents index
Next: Geometry optimization Up: Solvent calculations Previous: General considerations

Input description

 

tex2html_wrap9579

The necessary input for a solvent calculation is given in the **INTEGRALS   and **WAVE FUNCTIONS   input modules. A typical input file for an SCF calculation of the nuclear shielding constants of a molecule in a dielectric medium will look like  :

**DALTON INPUT
.RUN PROPERTIES
**INTEGRALS
*ONEINT
.MAX L
 10
**WAVE FUNCTIONS
.HF
*SOLVENT
.DIELECTRIC
 78.5
.MAX L
 10
.CAVITY
 3.98
**PROPERTIES
.SHIELD
*END OF INPUT

In **INTEGRALS   we request the evaluation of the undifferentiated solvent multipole integrals  as given in for instance Ref. [56] by the keyword .MAX L   in the *ONEINT   submodule. We request all integrals up to L=10 to be evaluated. This is needed if static or dynamic (response) properties calculations are to be done, but is not needed for a run of the wave function only (**WAVE FUNCTIONS  ).

In **WAVE FUNCTIONS   there is a separate input module for the solvent input, headed by the name *SOLVENT  . We refer to Sec. gif for a presentation of all possible keywords in this submodule. The interaction between the solute and the dielectric medium  is characterized by three parameters; the dielectric constant , the cavity  size, and the order of the multipole expansion . In the above input we have requested a dielectric constant of 78.5 (corresponding to water ) through the keyword .DIELECTRIC  , a cavity radius of 3.98 atomic units with the keyword .CAVITY  , and the multipole expansion is to include all terms up to L=10, as can be seen from the keyword .MAX L  . Note that this number cannot be larger than the number given for .MAX L   in the *ONEINT   input module.





Kenneth Ruud
Sat Apr 5 10:26:29 MET DST 1997