Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Next revision
Previous revision
ab_initio_multiple_spawning_dynamics [2020/06/11 18:17] – external edit 127.0.0.1ab_initio_multiple_spawning_dynamics [2022/08/16 14:15] (current) – external edit 127.0.0.1
Line 54: Line 54:
 Control.dat is a fortran NAMELIST file storing the main parameters (in au) of an AIMS simulation. At the end of the file, one can optionally specify user-defined atom types (see below). Control.dat is a fortran NAMELIST file storing the main parameters (in au) of an AIMS simulation. At the end of the file, one can optionally specify user-defined atom types (see below).
  
-Initialization+^  Initialization  ^^^
 ^Parameter             ^Options (Default)          ^Description                                      ^ ^Parameter             ^Options (Default)          ^Description                                      ^
-|Required parameters:                                                                           |+|Required parameters:  |||
 |''NumStates''     |''Integer''            |Number of electronic states.                     | |''NumStates''     |''Integer''            |Number of electronic states.                     |
 |''InitState''     |''Integer''            |Initial electronic state.                        | |''InitState''     |''Integer''            |Initial electronic state.                        |
 |''NumParticles''  |''Integer''            |Number of atoms.                                 | |''NumParticles''  |''Integer''            |Number of atoms.                                 |
-|Common parameters:    |                                                                           |+|Common parameters:    |||
 |''IMethod''       |''0''                  |First principles dynamics (No NACV).             | |''IMethod''       |''0''                  |First principles dynamics (No NACV).             |
 |                      |''%%(1)%%''                |CASSCF dynamics, with NACV.                      | |                      |''%%(1)%%''                |CASSCF dynamics, with NACV.                      |
 |''IRestart''      |''%%(0)/1%%''              |Restart the calculation when set to 1.           | |''IRestart''      |''%%(0)/1%%''              |Restart the calculation when set to 1.           |
-|''RestartTime''   |''%%Real (-100.0)%%''      |Time from which to restart in Checkpoint.txt     | +|''RestartTime''   |''%%Real (-100.0)%%''      |Time from which to restart in Checkpoint.txtIf disabled, use Last_Bundle.txt.                |
-|                      |                           |If disabled, use Last_Bundle.txt.                |+
 |''IRndSeed''      |''%%Integer (0)%%''        |Initial condition random number seed.            | |''IRndSeed''      |''%%Integer (0)%%''        |Initial condition random number seed.            |
 |''InitialCond''   |''%%(’NoSample’)%%''       |Read initial conditions from Geometry.dat.       | |''InitialCond''   |''%%(’NoSample’)%%''       |Read initial conditions from Geometry.dat.       |
-|                      |''%%’Wigner’%%''           |Sample from the Wigner distribution of harmonic  | +|                      |''%%’Wigner’%%''           |Sample from the Wigner distribution of harmonic normal modes in Frequencies.dat.                 |
-|                      |                           |normal modes in Frequencies.dat.                 |+
 |                      |''%%’Quasiclassical’%%''   |As above, for the Quasiclassical distribution.   | |                      |''%%’Quasiclassical’%%''   |As above, for the Quasiclassical distribution.   |
 |                      |''%%’Boltzmann’%%''        |As above, for the Boltzmann distribution.        | |                      |''%%’Boltzmann’%%''        |As above, for the Boltzmann distribution.        |
 |''Temperature''   |''%%Real (0.0)%%''         |Temperature (Kelvin) of the initial conditions.  | |''Temperature''   |''%%Real (0.0)%%''         |Temperature (Kelvin) of the initial conditions.  |
-|Expert parameters:    |                                                                           +|Expert parameters:    ||| 
-|''InitBright''    |''%%Logical (.false.)%%''  |Select optically brightest of ''InitState''  | +|''InitBright''    |''%%Logical (.false.)%%''  |Select optically brightest of ''InitState'' and ''IDark'' as initial state.              |
-|                      |                           |and ''IDark'' as initial state.              |+
 |''InitDark''      |''%%Logical (.false.)%%''  |As above, but for darkest state.                 | |''InitDark''      |''%%Logical (.false.)%%''  |As above, but for darkest state.                 |
 |''IDark''         |''%%Integer (-1)%%''       |See above.                                       | |''IDark''         |''%%Integer (-1)%%''       |See above.                                       |
-|''InitGap''       |''%%Real (0.0)%%''         |Select initial ground-excited energy gap         | +|''InitGap''       |''%%Real (0.0)%%''         |Select initial ground-excited energy gap to within $\pm 0.5\times$''InitGapWidth''  |
-|                      |                           |to within $\pm 0.5\times$''InitGapWidth''  |+
 |''InitGapWidth''  |''%%Real (0.0)%%''         |See above.                                       | |''InitGapWidth''  |''%%Real (0.0)%%''         |See above.                                       |
-|''IRestartTraj''  |''%%Integer array (0)%%''  |Subset of trajectories to restart from. Can      | +|''IRestartTraj''  |''%%Integer array (0)%%''  |Subset of trajectories to restart from. Can restart a single dead trajectory.                |
-|                      |                           |restart a single dead trajectory.                |+
 |''NumInitBasis''  |''%%Integer (1)%%''        |Number of initial coupled basis functions.       | |''NumInitBasis''  |''%%Integer (1)%%''        |Number of initial coupled basis functions.       |
-|''FirstGauss''    |''%%Logical (.false.)%%''  |Place the first basis function at the            | +|''FirstGauss''    |''%%Logical (.false.)%%''  |Place the first basis function at the molecular origin (that of Geometry.dat), with zero momentum.                              | 
-|                      |                           |molecular origin (that of Geometry.dat),         | +|''MirrorBasis''   |''%%Logical (.false.)%%''  |Mirror the initial basis on another electronic state, ''MirrorState''          |
-|                      |                           |with zero momentum.                              | +
-|''MirrorBasis''   |''%%Logical (.false.)%%''  |Mirror the initial basis on another              | +
-|                      |                           |electronic state, ''MirrorState''          |+
 |''MirrorState''   |''%%Integer (1)%%''        |See above.                                       | |''MirrorState''   |''%%Integer (1)%%''        |See above.                                       |
-|''EnergyAdjust''  |''%%Logical (.false.)%%''  |Adjust momenta of mirrored basis (along          | +|''EnergyAdjust''  |''%%Logical (.false.)%%''  |Adjust momenta of mirrored basis (along NACV) to match original basis energy.            | 
-|                      |                           |NACV) to match original basis energy.            | +|''SharpEnergy''   |''%%Logical (.false.)%%''  |Adjust momenta of all initial basis functions to match first basis function energy.            |
-|''SharpEnergy''   |''%%Logical (.false.)%%''  |Adjust momenta of all initial basis functions    | +
-|                      |                           |to match first basis function energy.            |+
  
  
-Integration and numerics+^Integration and numerics  ^^^
 ^Parameter                   ^Options (Default)              ^Description                                 ^ ^Parameter                   ^Options (Default)              ^Description                                 ^
-|Required parameters:        |                                                                          |+|Required parameters:        |||
 |''TimeStep''            |''Real''                   |Simulation timestep.                        | |''TimeStep''            |''Real''                   |Simulation timestep.                        |
 |''SimulationTime''      |''Real''                   |Simulation length.                          | |''SimulationTime''      |''Real''                   |Simulation length.                          |
-|Common parameters:          |                                                                          +|Common parameters:          ||| 
-|''ExShift''             |''%%Real (0.0)%%''             |Apply a shift of ''%%+ExShift%%'' to        | +|''ExShift''             |''%%Real (0.0)%%''             |Apply a shift of ''%%+ExShift%%'' to diagonals of Hamiltonian to improve amplitude/phase propagation. Good choice is average potential energy of ground and excited states.               | 
-|                            |                               |diagonals of Hamiltonian to improve         | +|''CoupTimeStep''        |''%%Real (TimeStep/4)%%''      |Timestep to use in regions of large nonadiabatic coupling.                      | 
-|                            |                               |amplitude/phase propagation. Good           | +|Expert parameters:          ||| 
-|                            |                               |choice is average potential energy          | +|''TimeStepRejection''   |''%%Logical (.true.)%%''       |Enables adaptive time-stepping when encountering integration problems.          | 
-|                            |                               |of ground and excited states.               | +|''MinTimeStep''         |''%%Real (CoupTimeStep/4)%%''  |Minimum timestep to use when adaptive time-stepping.                     | 
-|''CoupTimeStep''        |''%%Real (TimeStep/4)%%''      |Timestep to use in regions of large         | +|''DieOnMinStep''        |''%%Logical (.true.)%%''       |Stop calculation if minimum timestep is reached.                                 | 
-|                            |                               |nonadiabatic coupling.                      | +|''EnergyStepCons''      |''%%Real (0.005)%%''           |Allowed violation of classical energy conservation before timestep is rejected.                                   | 
-|Expert parameters:          |                                                                          +|''NormCons''            |''%%Real (0.25)%%''            |Permitted deviation in Norm from beginning of simulation.                    | 
-|''TimeStepRejection''   |''%%Logical (.true.)%%''       |Enables adaptive time-stepping when         | +|''RejectAllStateFlip''  |''%%Logical (.true.)%%''       |Reject a timestep if any two electronic states flip. If ''.false.'' then only reject if state is coupled to occupied state.                          | 
-|                            |                               |encountering integration problems.          | +|''OlapThresh''          |''%%Real (0.001)%%''           |Threshold overlap between two trajectories below which matrix elements not calculated.                    | 
-|''MinTimeStep''         |''%%Real (CoupTimeStep/4)%%''  |Minimum timestep to use when                | +|''RegThresh''           |''%%Real (0.0001)%%''          |Threshold to regularize the overlap matrix before inversion.                    | 
-|                            |                               |adaptive time-stepping.                     | +|''Constrain''           |''%%Logical (.false.)%%''      |Apply geometric constraints? If yes, supply Constraints.dat                      |
-|''DieOnMinStep''        |''%%Logical (.true.)%%''       |Stop calculation if minimum timestep        | +
-|                            |                               |is reached.                                 | +
-|''EnergyStepCons''      |''%%Real (0.005)%%''           |Allowed violation of classical energy       | +
-|                            |                               |conservation before timestep is             | +
-|                            |                               |rejected.                                   | +
-|''NormCons''            |''%%Real (0.25)%%''            |Permitted deviation in Norm from            | +
-|                            |                               |beginning of simulation.                    | +
-|''RejectAllStateFlip''  |''%%Logical (.true.)%%''       |Reject a timestep if any two                | +
-|                            |                               |electronic states flip. If ''.false.''  | +
-|                            |                               |then only reject if state is coupled        | +
-|                            |                               |to occupied state.                          | +
-|''OlapThresh''          |''%%Real (0.001)%%''           |Threshold overlap between two               | +
-|                            |                               |trajectories below which matrix             | +
-|                            |                               |elements not calculated.                    | +
-|''RegThresh''           |''%%Real (0.0001)%%''          |Threshold to regularize the overlap         | +
-|                            |                               |matrix before inversion.                    | +
-|''Constrain''           |''%%Logical (.false.)%%''      |Apply geometric constraints? If yes,        | +
-|                            |                               |supply Constraints.dat                      |+
  
-Spawning+^Spawning  ^^^
 ^Parameter                ^Options (Default)          ^Description                                           ^ ^Parameter                ^Options (Default)          ^Description                                           ^
-|Common parameters:                                                                                      |+|Common parameters:       |||
 |''CSThresh''         |''%%Real (0.1)%%''         |Threshold Coupling, above which to spawn.             | |''CSThresh''         |''%%Real (0.1)%%''         |Threshold Coupling, above which to spawn.             |
-|''PopToSpawn''       |''%%Real (0.1)%%''         |Minimum population a trajectory must have             | +|''PopToSpawn''       |''%%Real (0.1)%%''         |Minimum population a trajectory must have before it can spawn.                                  | 
-|                                                   |before it can spawn.                                  | +|''OMax''             |''%%Real (0.8)%%''         |Prevent spawning if overlap of spawned trajectory with existing basis greater than ''Omax''                                   | 
-|''OMax''             |''%%Real (0.8)%%''         |Prevent spawning if overlap of spawned                | +|''MaxTraj''          |''%%Real (100)%%''         |Maximum number of allowed trajectories. Above which, spawning is disabled.                    | 
-|                                                   |trajectory with existing basis greater                | +|''IgnoreState''      |''%%Integer (0)%%''        |Stop dynamics of trajectories on ''IgnoreState''if not coupled for at least ''DecoherenceTime''. Choose 1 if only interested in excited states.        | 
-|                                                   |than ''Omax''                                   | +|''MaxEDiff''         |''%%Real (0.03)%%''        |Energy gap threshold for calculating nonadiabatic couplings, above which, coupling taken to be zero.| 
-|''MaxTraj''          |''%%Real (100)%%''         |Maximum number of allowed trajectories.               | +|Expert parameters:       |||
-|                                                   |Above which, spawning is disabled.                    | +
-|''IgnoreState''      |''%%Integer (0)%%''        |Stop dynamics of trajectories on ''IgnoreState''  | +
-|                                                   |if not coupled for at least ''DecoherenceTime'' | +
-|                                                   |Choose 1 if only interested in excited states.        | +
-|''MaxEDiff''         |''%%Real (0.03)%%''        |Energy gap threshold for calculating                  | +
-|                                                   |nonadiabatic couplings, above which,                  | +
-|                                                   |coupling taken to be zero.                            +
-|Expert parameters:                                                                                      |+
 |''DecoherenceTime''  |''%%Real (200.0)%%''       |See above.                                            | |''DecoherenceTime''  |''%%Real (200.0)%%''       |See above.                                            |
-|''SpawnCoupV''       |''%%Logical (.false.)%%''  |Use Coup.velocity as criterion for spawning           | +|''SpawnCoupV''       |''%%Logical (.false.)%%''  |Use Coup.velocity as criterion for spawning rather than Coupling magnitude alone.                 | 
-|                                                   |rather than Coupling magnitude alone.                 | +|''EShellOnly''       |''%%Logical (.true.)%%''   |Require spawned trajectories to have same classical energy as parent by adjusting momentum along NACV.                                  |
-|''EShellOnly''       |''%%Logical (.true.)%%''   |Require spawned trajectories to have same             | +
-|                                                   |classical energy as parent by adjusting               | +
-|                                                   |momentum along NACV.                                  |+
  
-Output+^  Output  ^^^^
 ^Parameter                        ^Options (Default)          ^Description                                           ^ ^Parameter                        ^Options (Default)          ^Description                                           ^
-|Common parameters:                                                                                              |+|Common parameters:               |||
 |''zAllText''                 |''%%Logical (.false.)%%''  |Write all output files?                               | |''zAllText''                 |''%%Logical (.false.)%%''  |Write all output files?                               |
 |''zEDatFile''                |''%%Logical (.true.)%%''   |Write Energy file?                                    | |''zEDatFile''                |''%%Logical (.true.)%%''   |Write Energy file?                                    |
Line 173: Line 134:
 |''zTDipoleFile''             |''%%Logical (.false.)%%''  |Write Transition Dipole files?                        | |''zTDipoleFile''             |''%%Logical (.false.)%%''  |Write Transition Dipole files?                        |
 |''zQMRRFile''                |''%%Logical (.false.)%%''  |Write $\langle r^2 \rangle$ quadrupole moment files?  | |''zQMRRFile''                |''%%Logical (.false.)%%''  |Write $\langle r^2 \rangle$ quadrupole moment files?  |
-|''zCIVecFile''               |''%%Logical (.false.)%%''  |Write CI Vector files? (Not                           | +|''zCIVecFile''               |''%%Logical (.false.)%%''  |Write CI Vector files? (Not recommended for large CI)                             |
-|                                                           |recommended for large CI)                             |+
 |''zCorrfile''                |''%%Logical (.false.)%%''  |Write Correlation function file?                      | |''zCorrfile''                |''%%Logical (.false.)%%''  |Write Correlation function file?                      |
 |''zBundMatFile''             |''%%Logical (.false.)%%''  |Write Matrices (H, S, etc)?                           | |''zBundMatFile''             |''%%Logical (.false.)%%''  |Write Matrices (H, S, etc)?                           |
 |''WriteMolden''              |''%%Logical (.false.)%%''  |Write Molden files?                                   | |''WriteMolden''              |''%%Logical (.false.)%%''  |Write Molden files?                                   |
-|''MoldenStep''               |''%%Integer (200)%%''      |How often (every number of                            | +|''MoldenStep''               |''%%Integer (200)%%''      |How often (every number of timesteps) Molden files are written.                  | 
-|                                                           |timesteps) Molden files are written.                  | +|''RestartStep''              |''%%Integer (4)%%''        |How often (every number of timesteps) Restart files are written.                 | 
-|''RestartStep''              |''%%Integer (4)%%''        |How often (every number of                            | +|Expert parameters:               ||| 
-|                                                           |timesteps) Restart files are written.                 | +|''NStepToPrint''             |''%%Integer (1)%%''        |How often (every number of timesteps) output files are written.                  | 
-|Expert parameters:                                                                                              +|''WriteEveryStep''           |''%%Logical (.true.)%%''   |Write output every timestep, including subtimesteps of adaptive propagation. Overrides ''NStepToPrint''. | 
-|''NStepToPrint''             |''%%Integer (1)%%''        |How often (every number of                            | +|''NBonds''                   |''%%Integer (0)%%''        |Number of bond lengths to print. Define IBond below.                            | 
-|                                                           |timesteps) output files are written.                  | +|''%%IBond(2,NBonds)%%''          |''%%Integer (0)%%''        |List of atom indices that are bonded. See NBonds above.                         | 
-|''WriteEveryStep''           |''%%Logical (.true.)%%''   |Write output every timestep,                          | +|''NAngles''                  |''%%Integer (0)%%''        |Number of bond angles to print. Define IAngle below.                           | 
-|                                                           |including subtimesteps of                             | +|''%%IAngle(3,NAngles)%%''        |''%%Integer (0)%%''        |List of atom indices that define bond angles.                                   | 
-|                                                           |adaptive propagation. Overrides                       | +|''NDihedrals''               |''%%Integer (0)%%''        |Number of dihedral angles to print. Define IDihedral below.                        | 
-|                                                           |''NStepToPrint''                                +|''%%IDihedral(4,NDihedrals)%%''  |''%%Integer (0)%%''        |List of atom indices that define dihedral angles.                               | 
-|''NBonds''                   |''%%Integer (0)%%''        |Number of bond lengths to                             | +|''NPyrams''                  |''%%Integer (0)%%''        |Number of pyramidalization angles to print. Define IPyram below.                        | 
-|                                                           |print. Define IBond below.                            | +|''%%IPyram(4,NPyrams)%%''        |''%%Integer (0)%%''        |List of atom indices that define pyramidalization angles.                       |
-|''%%IBond(2,NBonds)%%''          |''%%Integer (0)%%''        |List of atom indices that                             | +
-|                                                           |are bonded. See NBonds above.                         | +
-|''NAngles''                  |''%%Integer (0)%%''        |Number of bond angles to                              | +
-|                                                           |print. Define IAngle below.                           | +
-|''%%IAngle(3,NAngles)%%''        |''%%Integer (0)%%''        |List of atom indices that                             | +
-|                                                           |define bond angles.                                   | +
-|''NDihedrals''               |''%%Integer (0)%%''        |Number of dihedral angles to                          | +
-|                                                           |print. Define IDihedral below.                        | +
-|''%%IDihedral(4,NDihedrals)%%''  |''%%Integer (0)%%''        |List of atom indices that                             | +
-|                                                           |define dihedral angles.                               | +
-|''NPyrams''                  |''%%Integer (0)%%''        |Number of pyramidalization angles                     | +
-|                                                           |to print. Define IPyram below.                        | +
-|''%%IPyram(4,NPyrams)%%''        |''%%Integer (0)%%''        |List of atom indices that                             | +
-|                                                           |define pyramidalization angles.                       |+
  
 **User-defined atom types** **User-defined atom types**
Line 296: Line 242:
 If ''Constraint=.true.'' then Constraints.dat must be supplied. This is a fortran NAMELIST file that defines the geometric constraints in the system. These constraints are enforced with the RATTLE algorithm. The parameters contained in this file are: If ''Constraint=.true.'' then Constraints.dat must be supplied. This is a fortran NAMELIST file that defines the geometric constraints in the system. These constraints are enforced with the RATTLE algorithm. The parameters contained in this file are:
  
-Constraint parameters+^Constraint parameters  ^^^
 ^Parameter                          ^Options (Default)        ^Description                           ^ ^Parameter                          ^Options (Default)        ^Description                           ^
-|''NumBonds''                   |''%%Integer (0)%%''      |Number of bond lengths to constrain   | +|''NumBonds''                   |''%%Integer (0)%%''      |Number of bond lengths to constrainDefine IBondPtcles below.             | 
-|                                                           |Define IBondPtcles below.             | +|''%%IBondPtcles(2,NumBonds)%%''    |''%%Integer (0)%%''      |List of atom indices that are bonded. See NumBonds above.           | 
-|''%%IBondPtcles(2,NumBonds)%%''    |''%%Integer (0)%%''      |List of atom indices that are         | +|''NumAngles''                  |''%%Integer (0)%%''      |Number of bond angles to constrain. Define IAnglePtcles below.            | 
-|                                                           |bonded. See NumBonds above.           | +|''%%IAnglePtcles(3,NumAngles)%%''  |''%%Integer (0)%%''      |List of atom indices that define bond angles.                   | 
-|''NumAngles''                  |''%%Integer (0)%%''      |Number of bond angles to constrain.   | +|''NComAtoms''                  |''%%Integer (0)%%''      |Number of atoms that define the center-of-mass constraint. Define ICOMPtcles below.                     | 
-|                                                           |Define IAnglePtcles below.            | +|''%%ICOMPtcles(NComAtoms)%%''      |''%%Integer (0)%%''      |List of atom indices that define the center of mass.                       |
-|''%%IAnglePtcles(3,NumAngles)%%''  |''%%Integer (0)%%''      |List of atom indices that             | +
-|                                                           |define bond angles.                   | +
-|''NComAtoms''                  |''%%Integer (0)%%''      |Number of atoms that define the       | +
-|                                                           |center-of-mass constraint. Define     | +
-|                                                           |ICOMPtcles below.                     | +
-|''%%ICOMPtcles(NComAtoms)%%''      |''%%Integer (0)%%''      |List of atom indices that define the  | +
-|                                                           |center of mass.                       |+
 |''Toler''                      |''%%Real (1E-4)%%''      |Rattle tolerance                      | |''Toler''                      |''%%Real (1E-4)%%''      |Rattle tolerance                      |
 |''MaxIts''                     |''%%Integer (10000)%%''  |Maximum number of Rattle iterations   | |''MaxIts''                     |''%%Integer (10000)%%''  |Maximum number of Rattle iterations   |
Line 520: Line 459:
   orbital,orbrec   orbital,orbrec
   ciguess,cirec   ciguess,cirec
-  save,ci=cirec+  save,det=cirec
   if (status.lt.0) then   if (status.lt.0) then
     text, MCSCF failed in MPINIT.     text, MCSCF failed in MPINIT.
     text, Attempting uncoupled iterations.     text, Attempting uncoupled iterations.
-    {multi+    {multi,so
     pspace,10.0     pspace,10.0
     wf;state,nstate     wf;state,nstate
     orbital,orbrec     orbital,orbrec
     ciguess,cirec     ciguess,cirec
-    save,ci=cirec +    save,det=cirec
-    {iterations +
-      do,uncouple,1,to,15}+
     maxiter,50}     maxiter,50}
   endif   endif
Line 561: Line 498:
 ciguess,cirec ciguess,cirec
 dm,orbrec dm,orbrec
-save,ci=cirec+save,det=cirec
 diab,diabrec,save=orbrec diab,diabrec,save=orbrec
 } }
Line 573: Line 510:
   ciguess,cirec   ciguess,cirec
   dm,orbrec   dm,orbrec
-  save,ci=cirec+  save,det=cirec
   {iterations   {iterations
     do,uncouple,1,to,10}     do,uncouple,1,to,10}
Line 596: Line 533:
   orbital,orbrec   orbital,orbrec
   ciguess,cirec   ciguess,cirec
-  save,ci=cirec+  save,det=cirec
   diab,diabrec,save=orbrec   diab,diabrec,save=orbrec
   cpmcscf,grad,stateind,save=gradrec0;   cpmcscf,grad,stateind,save=gradrec0;
Line 607: Line 544:
   orbital,orbrec   orbital,orbrec
   ciguess,cirec   ciguess,cirec
-  save,ci=cirec+  save,det=cirec
   diab,diabrec,save=orbrec   diab,diabrec,save=orbrec
   cpmcscf,grad,stateind,save=gradrec0;   cpmcscf,grad,stateind,save=gradrec0;
Line 619: Line 556:
   orbital,orbrec   orbital,orbrec
   ciguess,cirec   ciguess,cirec
-  save,ci=cirec+  save,det=cirec
   diab,diabrec,save=orbrec   diab,diabrec,save=orbrec
   cpmcscf,grad,stateind,save=gradrec0;   cpmcscf,grad,stateind,save=gradrec0;
Line 632: Line 569:
   orbital,orbrec   orbital,orbrec
   ciguess,cirec   ciguess,cirec
-  save,ci=cirec+  save,det=cirec
   diab,diabrec,save=orbrec   diab,diabrec,save=orbrec
   cpmcscf,grad,stateind,save=gradrec0;   cpmcscf,grad,stateind,save=gradrec0;