Those who are used to the old SIRIUS input will not experience any dramatic changes with respect to the DALTON input, although a few minor differences exist in order to create a more user-friendly environment. The input sections are the same as before starting with a wavefunction specification, currently, including options for SCF, MP2, and MCSCF reference states. Depending on this choice of reference state the following input sections take different form, from the simplest SCF input to the more complex MCSCF input. Minimal SCF and SCF+MP2 using the new feature of generating the HF occupation on the basis of an initial Hückel calculation and then possibly change the occupation during the first DIIS iterations will be shown. There will also be an example showing how an old f21 file is used for restart in the input. In this case the Hartree-Fock occupation will be then read from the f21-file and used as initial Hartree-Fock occupation . The HF occupation is needed for an MCSCF calculation , as this is anyway determined when establishing a suitably chosen active space . An example of an MCSCF-CAS input without starting orbitals will be give, as well as an MCSCF-RAS with starting orbitals. We note that and MCSCF wave function by default will be optimized using spin-adapted configurations (CSFs) , unless the wave function optimization is followed by a calculation of molecular properties or geometry optimization, when determinants will be used. The three following examples illustrate calculations on excited states , on a core hole state with frozen core orbital and on a core hole state with relaxed core orbital. The last example shows an MCSCF calculation of non-equilibrium solvation energy .
The following input example gives a minimal SCF input with default starting orbitals (that is, Hückel guess), and automatic Hartree-Fock occupation, first based on the Hückel guess, and then updated during the DIIS iterations:
**DALTON INPUT .RUN WAVE FUNCTIONS **WAVE FUNCTIONS .HF *END OF DALTON
The following input give a minimal input for an MP2 calculation using all default settings for the Hartree-Fock calculation (see previous example):
**DALTON INPUT .RUN WAVE FUNCTIONS **WAVE FUNCTION .HF .MP2 *END OF DALTON
If we would like to calculate molecular properties in several geometries, we may take advantage of the fact the molecular orbitals at the previous geometry probably is quite close to the optimized MOs at the new geometry, and thus restart from the MOs contained in the f21 file, as indicated in the following example:
**DALTON INPUT .RUN WAVE FUNCTIONS **WAVE FUNCTION .HF *ORBITAL INPUT .MOSTART NEWORB *END OF DALTON
The next input example gives the necessary input for an Complete Active Space SCF (CASSCF) calculation where we use MP2 to provide starting orbitals for our MCSCF. The active space may for instance be chosen on the basis of an MP2 natural orbital occupation analysis as described in Ref. [5]. The input would like:
**DALTON INPUT .RUN WAVE FUNCTIONS **WAVE FUNCTIONS .HF .MP2 .MCSCF *HF INPUT .HF OCCUPATION 2 0 0 0 *CONFIGURATION INPUT .SYMMET 1 .SPIN MULT 1 .INACTIVE 0 0 0 0 .ELECTRONS 4 .CAS SPACE 6 4 4 0 *ORBITAL INPUT .NOSUPSYM *END OF DALTON
As for Hartree-Fock calculation, we may want to use available molecular orbitals on the fort.21 file from previous calculations as starting orbitals for our MCSCF as indicated in this input example:
**DALTON INPUT .RUN WAVE FUNCTIONS **WAVE FUNCTION .MCSCF *HF INPUT .HF OCCUPATION 2 0 0 0 *CONFIGURATION INPUT .SYMMET 1 .SPIN MULT 1 .INACTIVE 0 0 0 0 .ELECTRONS 4 .RAS1 SPACE 2 1 1 0 .RAS2 SPACE 2 2 2 0 .RAS3 SPACE 6 4 4 2 .RAS1 ELECTRONS 0 2 .RAS3 ELECTRONS 0 2 *OPTIMIZATION .TRACI .FOCKONLY *ORBITAL INPUT .MOSTART NEWORB *END OF SIRIUS
The next input describes the optimization of the first excited state of the same symmetry as the ground state. To speed up convergence, we employ optimal orbital trial vectors as described in Ref. [6]. Such an input would look like:
**DALTON INPUT .RUN WAVE FUNCTIONS **WAVE FUNCTION .TITLE 4 2 2 0 CAS on first excited 1A_1 state, converging to 1.D-07 .MCSCF .NSYM 4 *CONFIGURATION INPUT .SYMMETRY 1 | same symmetry as ground state .SPIN MULTIPLICITY 1 .INACTIVE 1 0 0 0 .ELECTRONS 8 .CAS SPACE 4 2 2 0 *ORBITAL INPUT .MOSTART NEWORB *CI VECTOR .STARTHDIAG | Compute start vector from Hessian CI-diagonal *OPTIMIZATION .THRESHOLD 1.D-07 .SIMULTANEOUS ROOTS 2 2 .STATE 2 | 2 since the first exited state has the same symmetry ! as the ground state .OPTIMAL ORBITAL TRIAL VECTORS *PRINT LEVELS .PRINTUNITS 6 6 .PRINTLEVELS 5 5 **END OF DALTON
The next input describes the calculation of a core-hole state of the carbon 1s orbital in carbon monoxide , the first example employing a frozen core :
**DALTON INPUT .RUN WAVE FUNCTIONS **WAVE FUNCTION .TITLE C1s core hole state of CO, 4 2 2 0 valence CAS + C1s. Frozen core orbital calculation. .MCSCF .NSYM 4 *CONFIGURATION INPUT .SYMMETRY 1 .SPIN MULTIPLICITY 2 | doublet spin symmetry because of opened core orbital .INACTIVE 2 0 0 0 | O1s and O2s orbitals are inactive while the ! opened core orbital, C1s, must be active .ELECTRONS 9 | all valence electrons plus the core hole electron ! are active .RAS1 SPACE 1 0 0 0 | the opened core orbital (OBS always only this orbital) .RAS1 ELECTRONS 1 1 | We impose single occupancy for the opened core orbital .RAS2 SPACE 4 2 2 0 | The same as the CAS space in the ground state calculation *OPTIMIZATION .COREHOLE 1 2 | symmetry of core orbital and the orbital in this symmetry ! with the core hole according to list of input orbitals. ! The same thing could be obatined by reordering ! the core orbital to the first active orbital (by .REORDER), ! and specifying .FREEZE and .NEO ALWAYS .TRACI *ORBITAL INPUT .MOSTART NEWORB | start from corresponding MCSCF ground state *CI VECTOR .STARTHDIAG **END OF DALTON
whereas we in a calculation where we would allow the core to relax only would require the following changes compared to the previous input, assuming that we start out from orbitals and CI vectors generated by the previous calculation ".STARTHDIAG" is therefore replaced by ".STARTOLDCI", and the core orbital has number 3 (instead of 1) in the list of orbitals in the first symmetry. ".CORERELAX" is specifed for relaxation of the core orbital using the NR algorithm.
*OPTIMIZATION .COREHOLE 1 3 .CORERELAX *CI VEC .STARTOLDCI
This example describes calculations for non-equilibrium solvation . Ususally one starts with a calculation of a reference state (most often the ground state) with equilibrium solvation, using keyword "INESRF". The interface file is then used (without user interference) for a non-equilibrium excited state calculation; keyword "INERSI".
**DALTON INPUT .RUN WAVE FUNCTIONS **WAVE FUNCTION .TITLE 2-RAS(2p2p') : on F+ (1^D) in Glycol .MCSCF .NSYM 8 *CONFIGURATION INPUT .SPIN MULTIPLICITY 1 .SYMMETRY 1 .INACTIVE ORBITALS 1 0 0 0 0 0 0 0 .ELECTRONS 6 .RAS1 SPACE 0 0 0 0 0 0 0 0 .RAS2 SPACE 1 2 2 0 2 0 0 0 .RAS3 SPACE 8 4 4 3 4 3 3 1 .RAS1 ELECTRONS 0 0 .RAS3 ELECTRONS 0 2 *OPTIMIZATION .NEO ALWAYS .OPTIMA *ORBITAL INPUT .MOSTART NEWORB *CI VECTOR .STARTOLDCI *SOLVENT .CAVITY 2.5133D0 .INERSI | initial state inertial polarization 37.7D0 2.050D0 | static and optic dielectric constants for Glycol .MAX L 10 .PRINT 6 *END OF DALTON