This module controls the main features of the HF, SOPPA, and MCSCF property calculations, that is, which properties is to be calculated. In addition it includes directives affecting the performance of several of the program sections. This includes HF and MCSCF molecular gradients and Hessians. It should be noted, however, that the specification of what kind of walk (minimization, location of transition states, dynamical walks) is given in the *WALK or *OPTIMIZE submodules in the general input module. See also Chapter 7.
See Chapter 27 for specification of CC property calculations.
Note that RESPONSECHAPTER )
IS THE MOST GENERAL PART OF THE CODE FOR CALCULATING
MANY DIFFERENT ELECTRONIC LINEAR, QUADRATIC, OR CUBIC MOLECULAR
RESPONSE PROPERTIES BASED ON SCF, MCSCF, OR CI WAVE FUNCTIONS OR KOHN-SHAM DFT.
SOME OF THESE SCF/MCSCF PROPERTIES CAN ALSO BE REQUESTED
FROM THE **PROPERTIES INPUT MODULES DESCRIBED HERE.
NOTE: FOR SUCH PROPERTIES YOU SHOULD REQUEST THEM EITHER HERE OR
IN **RESPONSE, OTHERWISE YOU WILL CALCULATE THEM TWICE!
USUALLY THE OUTPUT IS NICEST HERE IN
THE **PROPERTIES MODULE (E.G. COLLECTED IN TABLES AND IN
OFTEN USED UNITS, MOST PROPERTIES ARE ONLY GIVEN IN ATOMIC
UNITS IN RESPONSE), AND NUCLEAR CONTRIBUTIONS ARE INCLUDED IF RELEVANT.
NO NUCLEAR CONTRIBUTIONS ARE ADDED IN RESPONSE.
SOME SPECIFIC PROPERTIES, ESPECIALLY THOSE INVOLVING NUCLEAR DERIVATIVES,
CAN ONLY BE CALCULATED VIA **PROPERTIES.
OTHER PROPERTIES, FOR EXAMPLE QUADRATIC AND CUBIC MOLECULAR RESPONSE FUNCTIONS
CAN ONLY BE CALCULATED IN THE **RESPONSE MODULE.
READ (LUCMD,*) (CAVORG(ICOOR), ICOOR = 1, 3)
READS THE ORIGIN TO BE USED FOR THE CAVITY
DURING A SOLVENT CALCULATION. BY DEFAULT THIS IS CHOSEN TO BE THE
CENTER OF MASS. SHOULD BY USED WITH CARE, AS IT
HAS TO CORRESPOND TO THE CENTER USED IN THE EVALUATION OF THE
UNDIFFERENTIATED SOLVENT INTEGRALS IN THE HERMIT SECTION, SEE CHAPTER .
READ (LUCMD, *) (DIPORG(ICOOR), ICOOR = 1, 3)
READS IN A USER DEFINED DIPOLE ORIGIN. THIS MAY
AFFECT PROPERTIES IN
WHICH CHANGES IN THE DIPOLE ORIGIN IS CANCELED BY
SIMILAR CHANGES IN THE
NUCLEAR PART. IT SHOULD ALSO BE USED WITH CARE, AS THE SAME DIPOLE
ORIGIN MUST BE USED DURING THE INTEGRAL EVALUATION SECTIONS, IN
PARTICULAR IF ONE IS DOING NUMERICAL
DIFFERENTIATION WITH RESPECT TO ELECTRIC FIELD PERTURBATIONS. FOR SUCH
FINITE-FIELD CALCULATIONS, WE REFER TO CHAPTER
, WHICH DEALS WITH FINITE FIELD CALCULATIONS. IT IS
PRIMARILY USED FOR DEBUGGING.
READ (LUCMD, *) (GAGORG(ICOOR), ICOOR = 1, 3)
READS IN A USER DEFINED GAUGE ORIGIN AND OVERWRITE BOTH THE .NOCMC OPTION, AS WELL AS THE DEFAULT VALUE OF CENTER-OF-MASS COORDINATES. NOTE THAT AN UNSYMMETRIC POSITION OF THE GAUGE ORIGIN WILL LEAD TO WRONG RESULTS IN CALCULATIONS EMPLOYING SYMMETRY, AS THE PROGRAM WILL NOT BE ABLE TO DETECT THAT SUCH A CHOICE OF GAUGE ORIGIN BREAKS THE SYMMETRY OF THE MOLECULE.
.VIBANA
ALSO IS
ACTIVATED.
FURTHERMORE, THE NATURAL CONNECTION (REF. [,]) IS THE DEFAULT IN ORDER TO ENSURE NUMERICALLY STABLE RESULTS. THE NATURAL CONNECTION CAN BE TURNED OFF BY THE KEYWORD .NODIFC, IN WHICH CASE THE SYMMETRIC CONNECTION WILL BE USED.
THE GAUGE ORIGIN IS CHOSEN TO BE THE CENTER OF MASS OF THE MOLECULE. THIS ORIGIN CAN BE CHANGED BY THE TWO KEYWORDS .GAUGEO AND .NOCMC. THIS WILL OF COURSE NOT AFFECT THE TOTAL MAGNETIZABILITY, ONLY THE MAGNITUDE OF THE DIA- AND PARAMAGNETIC TERMS.
COMBINED WITH THE KEYWORD .SOPPA OR .SOPPA(CCSD) IT INVOKES A SOPPA OR SOPPA(CCSD) CALCULATION OF THE MAGNETIZABILITY AND THE ROTATIONAL G TENSOR (REF. [61]). LONDON ORBITALS ARE AUTOMATICALLY DISABLED IN SOPPA OR SOPPA(CCSD) CALCULATIONS.
.MAGNET IN COMBINATION WITH THE KEYWORD .CTOCD INVOKES A CALCULATION WITHOUT THE USE OF LONDON ORBITALS BOTH WITH THE CTOCD-DZ METHOD (REF [,]) AND WITH THE COMMON ORIGIN METHOD. CHANGING THE DEFAULT VALUE OF THE GAUGE ORIGIN COULD GIVE WRONG RESULTS!
BY DEFINITION THE GAUGE ORIGIN OF THE MOLECULAR
G-FACTOR IS TO BE THE
CENTER OF MASS OF THE MOLECULE, AND ALTHOUGH THE
GAUGE ORIGIN CAN BE
CHANGED THROUGH THE KEYWORDS .NOCMC
AND .GAUGEO
, THIS
IS NOT RECOMMENDED, AND MAY GIVE ERRONEOUS RESULTS.
NOTE THAT IF THE ISOTOPIC CONSTITUTION OF THE MOLECULE IS SUCH THAT THE VIBRATIONAL WAVE FUNCTION HAS LOWER SYMMETRY THAN THE ELECTRONIC WAVE FUNCTION, CARE MUST BE TAKEN TO ENSURE THE SYMMETRY CORRESPONDS TO THE SYMMETRY OF THE NUCLEAR FRAMEWORK. THE AUTOMATIC SYMMETRY DETECTION ROUTINES WILL IN GENERAL ENSURE THAT THIS IS THE CASE.
COMBINED WITH THE KEYWORD .SOPPA OR .SOPPA(CCSD) IT INVOKES A SOPPA OR SOPPA(CCSD) CALCULATION OF THE MAGNETIZABILITY AND THE ROTATIONAL G TENSOR (REF. [61]). LONDON ORBITALS ARE AUTOMATICALLY DISABLED IN SOPPA OR SOPPA(CCSD) CALCULATIONS.
AS THE SYMMETRIC CONNECTION MAY GIVE NUMERICALLY INACCURATE RESULTS, IT'S USE IS NOT RECOMMENDED FOR OTHER THAN COMPARISONS WITH OTHER PROGRAMS.
*ABALNR
INPUT
MODULE. SEE ALSO DESCRIPTION IN CHAPTER
READ (LUCMD, *) (ORIGIN(ICOOR), ICOOR = 1, 3)
CHANGES THE ORIGIN OF THE PHASE-FACTORS ENTERING THE LONDON ATOMIC ORBITALS. THIS WILL CHANGE THE VALUE OF ALL OF THE CONTRIBUTIONS TO THE DIFFERENT MAGNETIC FIELD DEPENDENT PROPERTIES WHEN USING LONDON ATOMIC ORBITALS, BUT THE TOTAL MAGNETIC PROPERTIES WILL REMAIN UNCHANGED. TO BE USED FOR DEBUGGING PURPOSES ONLY.
READ (LUCMD, *) IPRDEF
SET DEFAULT PRINT LEVEL FOR THE CALCULATION. READ ONE
MORE LINE CONTAINING PRINT LEVEL. DEFAULT PRINT LEVEL IS THE
VALUE OF IPRDEF
FROM THE GENERAL INPUT MODULE.
READ (LUCMD, *) NREPS
READ (LUCMD, *) (IDOSYM(I),I = 1, NREPS)
CONSIDER PERTURBATIONS OF SELECTED SYMMETRIES ONLY. READ ONE MORE LINE SPECIFYING HOW MANY SYMMETRIES, THEN ONE LINE LISTING THE DESIRED SYMMETRIES. THIS OPTION IS CURRENTLY ONLY IMPLEMENTED FOR GEOMETRIC PERTURBATIONS.
READ (LUCMD,*) NPERT
READ (LUCMD, *) (IPOINT(I),I=1,NPERT)
SELECT WHICH NUCLEAR GEOMETRIC PERTURBATIONS ARE TO BE CONSIDERED. READ ONE MORE LINE SPECIFYING HOW MANY PERTURBATIONS, THEN ON A NEW LINE THE LIST OF PERTURBATIONS TO BE CONSIDERED. BY DEFAULT, ALL PERTURBATIONS ARE TO BE CONSIDERED, BUT BY INVOKING THIS KEYWORD, ONLY THOSE PERTURBATIONS SPECIFIED IN THE SEQUENCE WILL BE CONSIDERED.
THE PERTURBATION ORDERING FOLLOWS THE ORDERING OF THE SYMMETRIZED
NUCLEAR COORDINATES. THIS ORDERING CAN BE OBTAINED BY SETTING THE
PRINT LEVEL IN THE *READIN
MODULE TO 20 OR HIGHER.
.NOLOND
.
FURTHERMORE, THE NATURAL CONNECTION
(REF. [,]) IS THE DEFAULT IN ORDER TO ENSURE
NUMERICALLY STABLE RESULTS AS WELL AS PHYSICALLY INTERPRETABLE
RESULTS FOR THE PARAMAGNETIC AND DIAMAGNETIC TERMS. THE NATURAL
CONNECTION CAN BE TURNED OFF BY THE KEYWORD .NODIFC
IN WHICH
CASE THE SYMMETRIC CONNECTION IS USED INSTEAD.
THE GAUGE ORIGIN IS CHOSEN TO BE THE CENTER OF
MASS OF THE MOLECULE.
THIS ORIGIN CAN BE CHANGED BY THE TWO KEYWORDS .GAUGEO
AND
.NOCMC
. THIS CHOICE OF GAUGE ORIGIN WILL NOT AFFECT
THE FINAL SHIELDINGS IF LONDON ORBITALS ARE USED, ONLY THE SIZE OF THE
DIA- AND PARAMAGNETIC CONTRIBUTIONS.
COMBINED WITH THE KEYWORD .SOPPA OR .SOPPA(CCSD) IT INVOKES A SOPPA OR SOPPA(CCSD) CALCULATION OF THE NUCLEAR MAGNETIC SHIELDINGS (REF. [,,]). LONDON ORBITALS ARE AUTOMATICALLY DISABLED IN SOPPA OR SOPPA(CCSD) CALCULATIONS. GAUGE ORIGIN INDEPENDENT SOPPA OR SOPPA(CCSD) CALCULATIONS OF NUCLEAR MAGNETIC SHIELDINGS CAN BE CARRIED OUT WITH THE CTOCD-DZ METHOD (SEE REFS. [,,]) USING THE KEYWORD .CTOCD.
IN COMBINATION WITH THE KEYWORD .CTOCD THIS INVOKES A CALCULATION OF THE NUCLEAR MAGNETIC SHIELDINGS WITHOUT THE USE OF LONDON ORBITALS BUT WITH BOTH THE CTOCD-DZ METHOD (REF. [,,]) AND WITH THE COMMON ORIGIN METHOD. FOR THE CTOCD-DZ METHOD THE NUCLEAR MAGNETIC SHIELDINGS ARE GIVEN IN THE OUTPUT FILE FOR BOTH THE ORIGIN AT THE CENTER OF MASS AND AT THE RESPECTIVE ATOMS. CHANGING THE DEFAULT VALUE OF THE GAUGE ORIGIN COULD GIVE WRONG RESULTS!
.NOLOND
.
BY DEFINITION THE GAUGE ORIGIN OF THE
SPIN-ROTATION CONSTANT IS TO BE THE
CENTER OF MASS OF THE MOLECULE, AND ALTHOUGH THE
GAUGE ORIGIN CAN BE
CHANGED THROUGH THE KEYWORDS .NOCMC
AND .GAUGEO
, THIS
IS NOT RECOMMENDED, AND MAY GIVE ERRONEOUS RESULTS.
IN THE CURRENT IMPLEMENTATION, SYMMETRY DEPENDENT NUCLEI CANNOT BE USED DURING THE CALCULATION OF SPIN-ROTATION CONSTANTS.
AS THIS IS A VERY TIME CONSUMING PROPERTY, IT IS RECOMMENDED TO
CONSULT THE CHAPTER DESCRIBING THE CALCULATION OF NMR-PARAMETERS
(CH. ). THE MAIN CONTROL OF WHICH
CONTRIBUTIONS AND WHICH NUCLEI TO CALCULATE SPIN-SPIN COUPLINGS
BETWEEN IS DONE IN THE
*SPIN-S
MODULE.
.NOLOND
.
FURTHERMORE, THE NATURAL CONNECTION
(REF. [,]) IS DEFAULT IN ORDER TO ENSURE
NUMERICALLY STABLE RESULTS. THE NATURAL
CONNECTION CAN BE TURNED OFF BY THE KEYWORD .NODIFC
IN WHICH
CASE THE SYMMETRIC CONNECTION WILL BE USED.
IN THE CURRENT IMPLEMENTATION, THE KEYWORD .NOCMC WILL BE SET TRUE IN CALCULATIONS OF VIBRATIONAL CIRCULAR DICHROISM, THAT IS, THE COORDINATE SYSTEM ORIGIN WILL BE USED AS GAUGE ORIGIN. CHANGING THIS DEFAULT VALUE WILL GIVE INCORRECT RESULTS FOR VCD.
NOTE THAT IN THE CURRENT RELEASE, VCD IS NOT IMPLEMENTED FOR DENSITY FUNCTIONAL THEORY CALCULATIONS, AND THE PROGRAM WILL STOP IF VCD IS REQUESTED FOR A DFT CALCULATION.
.DIPGRA
THE IR INTENSITIES
WILL BE CALCULATED AS WELL.