Locating Conical Intersections (CIs) and Interstate Crossings (ISCs) with Firefly

 

               Conical Intersections (CIs) and Interstate Crossings (ISCs) are the specific manifolds on the Potential Energy Surfaces (PESs) of two different electronic states where, at the same geometry, the energies of these states are strictly equal each other and thus two states are exactly degenerate. The difference between CIs and ISCs is that ISCs are always formed by the states of different spatial or spin symmetry. As a result these states belong to the different blocks of the Hamiltonian matrix expressed in a block-diagonal form using symmetry-adapted (i.e. adapted according to the different possible irreducible representations) basis set. Let us consider small perturbations of the Hamiltonian matrix caused by the various deformations of the selected molecular geometry that preserve the initial molecular symmetry and hence the symmetry of the Hamiltonian. Evidently, there is no any "interaction" possible between the two states of interest, even exactly at their crossing points, as any "interaction" which could be caused by such perturbations is vanished due to symmetry reasons. As the result, the dimension of the ISC manifold is equal to the dimension of the original PES minus one, as there is in fact only one constraint (i.e. Ei = Ej) to be satisfied.

               By contract, CIs are formed by the states of exactly the same symmetry so that they are allowed to "interact". Mathematically, the second constraint arises: "interactions" of that kind must vanish at the CI manifold as otherwise the states cannot be exactly degenerate. The dimension of CI manifold is thus equal to the dimension of the original PES minus two, and the local topology on a two-dimension slice over the crossing manifold is the same as for two crossing cones with a single point in common, hence the "Conical Intersection" name. One could argue there is still only one constraint, namely: Ei = Ej at any CI point. While this is formally true, for CIs the equality Ei = Ej is in fact of the type: x+∙x + y+∙y = 0. The latter, being a single equality, implies both x = 0 and y = 0!

               Usually, there is no need to locate all the geometries belonging to ISC or CI manifolds. Quite often, one is interested in locating just the Minimum Energy Crossing Point (MECP or MECI). These are specific ISC or CI points having the minimum possible energy. Unlike CIs and ISCs, MECP and MECI are not the multidimensional manifolds of complex structure but rather are just the single point on PES. MECPs and MECIs are typically exactly what the most people calls ISCs and CIs.

               With Firefly versions prior to Firefly v. 8.0.0, it was possible to locate MECIs and MECPs only at the state-averaged MCSCF level of theory, using an efficient procedure developed for calculation of semi-numerical state-specific gradients for SA-MCSCF, which is documented elsewhere. In addition, only two intersecting states were allowed to be averaged, naturally with the equal weights, in the SA-MCSCF procedure. With Firefly v. 8.0.0, the MECIs/MECPs location code was extended to optionally use numerical gradients, thus allowing XMCQDPT2, MCQDPT2, and any configuration interaction procedure to be used for location of MECIs and MECPs - in addition to plain SA-MCSCF. Moreover, an arbitrary averaging of states can now be used for the purpose of location of MECIs and MECPs at the SA-MCSCF level provided the weights of two crossing states are equal. However, the energies of both intersecting states must be obtained as a result of a single computation (e.g. a single SA-MCSCF or XMCQDPT2 procedure). The primary reason for this is to avoid an unbalanced description of two different states of interest which could arise from the use of different computational schemes for different states.

               There are two very different ways the MECIs/MECPs can be located with Firefly. The first one is based on the use of penalty functions, the second one is based on the use of the Lagrange multiplier technique. Both of them can be used with Cartesian and internal coordinates (including DLCs), and can handle additional geometry constraints like frozen coordinates etc.

               The penalty function based methods programmed in Firefly are based on the minimization of the following expression (where ΔE = Ej-Ei):

               F(Ei, Ej) = alpha*Ei + (1-alpha)*Ej + Penalty(ΔE) = min                                                                    (1)

rather than the minimization of the Ei itself, as it would be in the case of unconstrained energy optimization. With Firefly, the possible values of alpha are 1.0 and 0.5. 1.0 corresponds to the minimization of the energy of the first state of interest, 0.5 corresponds to the minimization of the average energy of two states. Evidently, at CI point this is exactly the same but differs when away from the CI geometry. The Penalty(ΔE) is the penalty function which depends on the energy splitting between two states.

               The well designed penalty function must obey some restrictions. It should be zero if ΔE is zero and should be (very) small if ΔE is small enough. At the same time, the penalty function and its gradient should become (very) large when ΔE increases. In addition, it is convenient if Penalty(ΔE) is smooth and differentiable function as all standard geometry optimization engines assume that PES is a smooth and differentiable function. By forcing Penalty(ΔE) be smooth enough it is possible to use any standard geometry optimization code for location of MECPs and MECIs.

               Strictly speaking, the minimum of F(Ei, Ej) does not exactly correspond to MECI or MECP point and Ei and Ej are not exactly equal at the corresponding geometry. However, for well designed penalty function this point is very close to true MECI/MECP.

               There are three different penalty functions presently available in Firefly. The choice of the penalty function is controlled by the penlty variable of the $mcaver input group. Each of three penalty functions is parameterized by two parameters which can be altered by the user. See the description of the $mcaver input group below for details.

               The method based on the Lagrange multiplier technique as implemented in Firefly is the original self-consistent variant of the Sequential Quadratic Programming approach. Basically, the energy of state i (i.e. Ei) is optimized subject to constraint Ei = Ej . The constraint is taken into account by constructing the proper Lagrange function:

               L(R) = Ei(R) + λ(R)*(Ej(R) - Ei(R))                                                                                                            (2)

where λ is the Lagrange multiplier for constraint Ei = Ej and (R) denotes the dependence of all quantities on the molecular geometry. The Lagrange function is then optimized. Its minimum is located exactly at the MECI/MECP geometry while the computed value of λ contains important information on the type of located MECI/MECP.

               More precisely, if λ is less than 1.0, gradients of two degenerated states at their MECI/MECP point are always oppositely directed. Otherwise, they are directed in the same direction.

               Naturally, one can consider another form of L:

               L'(R) = (Ei(R) + Ej(R))*0.5 + λ'(R)*(Ej(R) - Ei(R))                                                                                   (3)

Optimization of both L(R) and L'(R) result in the same MECI/MECP geometry, while the values of two Lagrange multipliers at this geometry differ by 0.5 (i.e. λ = λ' + 0.5). However, the number of steps required to converge optimization is generally different. The trajectories formed by the sequential points taken by the optimization procedure are different as well.

               With Firefly, both penalty function and Lagrange multiplier based code require only state specific energy gradients for both MECI and MECP location. None of the approaches requires computation of the so-called non-adiabatic coupling vector. This vector corresponds to the second direction (the first being the difference of gradients of two states) along which the CI manifold (i.e. the manifold of the exact degeneracy, which is also called seam) is destroyed under the small variations of molecular geometry. In another words, this vector describes the second constraint which is specific to CIs but not to ISCs.

               Taken together with gradients difference non-adiabatic coupling vector forms the so-called "branching plane".

               Because non-adiabatic coupling vector is not used nor it is computed by the Firefly's routines, both MECIs and MECPs are handled by the same code and are virtually identical from the point of view of MECI/MECP geometry optimizer. In Firefly, the MECI or MECP location is not implemented as a separate runtyp. The standard $contrl runtyp=optimize is used instead together with some additional keywords controlling the process of MECI/MECP optimization. In addition to runtyp=optimize, the runtyp=rsurface is also allowed. While the penalty function based approach can be used with any of the available Firefly's geometry optimization engines such as GDIIS, QA, or NR, the Lagrange multiplier-based approach requires the dedicated geometry optimizer which is called "conic" and should be explicitly requested as follows: $statpt method=conic

               The $mcaver input group is relevant for both penalty function and Lagrange multiplier based code. In addition, the Lagrange multiplier based approach is also controlled by the dedicated $conic input group. If MECI/MECP optimization is based on fully numerical gradients (e.g. optimization of MECI at XMCQDPT2 level), the istate, jstate, and ngrads keywords of the $numgrd input group should be properly set as described in the documentation on the numerical gradient code.

               Alternatively, if semi-numerical state-specific gradients for SA-MCSCF are used, the istate keyword of the $mcscf input group must be properly set, as well as other variables that are essential for computation of state-specific gradients for SA-MCSCF.

               Finally, when locating CIs/MECPs  at SA-MCSCF level of theory, the MCSCF state tracking can be of considerable help. The latter is controlled by the $track input group as described elsewhere.

The $mcaver input group

 

jstate                   Selects the state number jstate, the second state of interest. Normally should be equal to istate + 1 unless you are interested in search of the simultaneous crossing of several states. The default value is jstate=0 i.e. no second state.

multiw                 Logical variable. Should be set for runs computing semi-numerical SS gradients for SA-MCSCF for two states (i.e. states istate and jstate) and there are other states in averaging. If not set, these runs will abort after performing some preliminary setup. If there are only states # istate and # jstate in MCSCF averaging, this flag should not be used. The default value is .false.  i.e. to assume there are no other states in averaging. Internally, this flag directs Firefly to compute SS gradients for two states independently on each other. Otherwise, the SS gradient is computed only for state # istate then gradient for state # jstate is computed based on the averaged effective gradient and on the gradient of state # istate. The second approach is approximately 5/3 times faster but is only possible if there are no additional states in state averaging. 

shift                     The value of energy difference shift during location of MECIs/MECPs. Firefly always searches for geometries such that Ej(R) - Ei(R) = shift. The default value of shift is 0.0 i.e. to search for true MECIs/MECPs. This option only affects penalty function based code.

target                  Symbolic variable which can be set to either "mixed", "pure", or "pure2" (without quotes). Any other values are ignored. "Mixed" means to set alpha in eq. (1) to 0.5 i.e. to use half-sum of state energies. This also sets eq. (3) for computation of Lagrange function. "Pure" means to use pure energy of state # istate in eq. (1) and to use eq. (2). "Pure2" means to use pure energy of state # jstate in eq. (1) while still using eq. (2) for Lagrange multiplier based method. The default value is target=mixed.

penlty                  Integer variable. Valid values are 1, 2, or 3. Selects one of the three programmed penalty functions. The default value is penlty=1. Note these functions will only be used in MECP/MECI location runs which do not use Lagrange multiplier approach.

The first penalty function has the form:

Penalty(Δ) = B*ln(1.0 + (Δ/A)2), the so-called Ciminelli penalty function.

(Note with Firefly versions prior to Firefly v. 8.0.0, the A and B parameters must be exchanged for this penalty function).

The second has the form suggested by Levine, Coe, and Martinez:

Penalty(Δ) = B*(Δ2)/(Δ+A)

Finally, the third one is specific to Firefly:

Penalty(Δ) = B*((Δ2 + A2)1/2) - B*|A|

Here A and B are the parameters defining the penalty function (see below).

a                           The parameter A of the selected penalty function.

b                           The parameter B of the selected penalty function.

               The default values of A and B depend on the particular penalty function in use.

               For penlty=1 these are A=0.008, B=0.2,

               for penlty=2 these are A=0.02, B=3.5,

               for penlty=2 these are A=0.01, B=3.5.

xgrad                   Logical variable. Should be set to request the use of Lagrange multiplier based code. Default is .true. if $statpt method=conic. Otherwise it is .false.

xhess                   Logical variable. If set, requests the use of the special optimized form of the approximate hessian used by all quasi-Newton type geometry optimizers. This option is only supported for penlty=1 penalty function and does not have effect otherwise. The default value is xhess=.false.

ssgrad                 Logical variable. If set to .false., disables all features requesting Firefly's state-specific gradients for SA-MCSCF. The jobs depending on these features will then abort. This option is now obsolete and should not be used. The default value is ssgrad=.true.

Note this input group contains some additional keywords controlling SS gradients for SA-MCSCF. These keywords are described elsewhere.

The $conic input group

 

               This group controls the detailed behavior of the $statpt method=CONIC Lagrange multiplier based CIs/ISCs optimizer called CONIC. In addition to the parameters described below, CONIC is also controlled by the parameters of the $statpt group.

               shift                     The final target value of Δ = Ej - Ei The default value is 0.0001 Hartree. For true CIs, it is almost impossible to reliably converge calculations to values of Δ less than ca. 0.0001 Hartree. For ISCs, it is possible to use tighter convergence criteria, up to shift=1.0d-7 or 1.0d-8. The difference in convergence behavior between true CIs and ISCs is related to non-zero non-adiabatic coupling between quasi-degenerate states in the case of CIs and the absence of this coupling in the case of ISCs.

               shift0                   The intermediate target value of Δ = Ej - Ei The default value is 0.0001 Hartree. For CIs, the recommended value is 0.001 Hartree. For ISCs, the recommendation is to always set shift0 to the same value as shift. The CONIC optimizer first tries to converge CIs/ISCs to Δ values equal to or below of shift0 using loosened values of various convergence thresholds. During this process (which is internally called "phase 1"), the CONIC optimizer gathers information on the behavior of the particular system of interest. Upon achieving the initial convergence, the optimizer switches to the original values of various convergence thresholds and continue optimization trying to achieve tighter convergence with the value of Δ equal to or below of shift. This is internally called "phase 2". Here, the previously gathered information is used to achieve the final convergence. In particular, the first few optimization steps of phase 2 attempt to rapidly decrease Δ using directed steps along the "most significant" modes while the rest of the steps take care of the final relaxation of the "less significant" modes. 

               hdgss                   The value to put on the diagonal of the initial diagonal guess to the Hessian matrix of the Lagrange function. The default value is 1.0. It is ignored if internal coordinates are used for optimization. In the latter case, the standard Firefly's guess for Hessian in internal coordinates is used.

               tde                       The value of the delta energy convergence threshold. It defines the first of the five tests all of which needs to be positive for the optimization process to complete successfully. This first test is satisfied if |tde| is less than |Δ - shift|. The default value of tde is 1.0d-5

               tdxmax               The value of the maximum primitive coordinate displacement convergence threshold. It defines the second of the five convergence tests. The displacement along each primitive coordinate should be less than tdxmax. The default value is 20.0*opttol.

               tdxrms                 The value of the RMS coordinate displacement convergence threshold. It defines the third of the five convergence tests. The RMS displacement should be less than tdxrms. The default value is 15.0*opttol.

               tgmax                 The threshold value for the maximum primitive component of the gradient of Lagrange function. It defines the fourth of the five convergence tests. The component of gradient along each primitive coordinate should be less than tgmax. The default value is 5.0*opttol.

               tgrms                   The threshold value for the RMS gradient of Lagrange function. It defines the last of the five convergence tests. The RMS gradient should be less than tgrms. The default value is 3.0*opttol.

               stpsiz                    The maximum permissible RMS length of the geometry update step. The default is 0.05 Bohr.

               demax                 The maximum permissible change of the Lagrange function on every geometry update step. As the change of Lagrange function can only be estimated approximately, demax is applied to this estimate, not to the actual yet unknown change. The default value is 0.03 Hartree.

               reset                    Parameter controlling whether and when to reset the current approximation to the Hessian of the Lagrange function.  Reset is performed automatically using the value of reset in some internal consistency checks. The default value for reset is 10.0. Smaller values will reset approximate Hessian more often. To disable reset-based tests, use reset=-1.0. The disabled reset may have both negative and positive effect on the convergence.

               lstart                    Provides the initial value for Lagrange multiplier λ. The default value is 0.0 i.e. to automatically compute a starting approximation to λ rather than to use user's input. This parameter can be used to facilitate restart of CONIC runs.

               lmax                    The upper bound on the absolute value of Lagrange multiplier λ: |λ| ≤ lmax. The default value of lmax is 3.0

               track                    Logical flag. If enabled, activates additional state tracking which is based on the analysis of the values of Lagrange multipliers. This additional tracking can correct mis-predictions of the base tracking and detects flip-flop of states. Default is .true. i.e. this additional tracking is enabled.

               stage2                 Logical flag. If set, tells CONIC code that CI/ISC optimization should be started directly from the "phase 2" of algorithm thus skipping the first phase. This flag is useful for restarts. Default is .false.

               ndown                Integer variable. Defines the number of directed steps performed at the beginning of "phase 2" of CONIC algorithm in order to decrease Δ to its final target value. The more steps performed the less is the "strain" on each particular step so that the molecular system of interest has more possibilities to relax during this procedure. The default is ndown=5.

               radown                              Defines the relative decrease of Δ which CONIC code attempts to achieve on each step at the beginning of "phase 2" of CONIC algorithm in order to decrease Δ to its final target value. If non-zero, triggers the alternative approach to directed decrease of Δ, namely the exponential decrease: Δnext = Δcurrent*(1.0-radown). The default value is 0.0 i.e. not to use this approach. The recommended values are in the range 0.1 to 0.2

               trdown                The multiplier for trman and trmin variables of  $statpt group which is applied when switching from "phase 1" to "phase 2" of CONIC code. The default is 0.2

               filter                     The upper limit on the absolute values of eigenvalues of the approximate Hessian matrix of Lagrange function. Any larger values are filtered off and zeroed. The default value is 5.0d+5

               diabas                                Logical flag. If set, directs CONIC code to construct approximate diabatic states based on the extremals  of gradients and to use these diabatic states rather than adiabatic ones during optimization of MECIs/MECPs. Note, with this option turned on, the value of λ at MECI/MECP point will always be equal 1.0 so it cannot be longer used to characterize MECI/MECP type. Default is .false. This option is still experimental as of Firefly version 8.0.0.

               tols                       Double precision array of undocumented length containing undocumented parameters defining the intimate details of CONIC algorithm's behavior. These parameters are implementation-dependent and thus are not documented as of now.

 

Sample inputs and outputs

 

Sample 1. Optimization of ISC in 1,1-difluoro-1,3-butadiene molecule using Lagrange multiplier based code and semi-numerical gradients for SA-MCSCF.

 $CONTRL SCFTYP=mcscf RUNTYP=optimize inttyp=hondo icut=11 d5=.t. nzvar=1 $END                                                                 

 $SYSTEM TIMLIM=6000 MWORDS=10 nojac=1 kdiag=0 $END                             

 $moorth nostf=1 nozero=1 tole=0 tolz=0 syms=1 symden=1 symvec=1 $end                                   

 $GUESS GUESS=MOREAD NORB=26 $END   

 $trans cuttrf=1d-12 $end                           

 $DRT GROUP=C1 FORS=.TRUE. NMCC=18 NDOC=5 NVAL=3 $END                          

 $GUGDIA NSTATE=2  ITERMX=200 cvgtol=1d-7 $END                                 

 $GUGDM2 cutoff=1d-12 WSTATE(1)=1,1 $END                                       

 $GUGEM  cutoff=1d-12 $END                                                      

 $MCSCF CISTEP=GUGA MAXIT=350 ISTATE=1 acurcy=5d-8 ENGTOL=5.0d-13 ntrack=2 $END                                                     

 $mcaver jstate=2 conic=2 $end                                    

 $track tol=1.4 freeze=.t. sticky=.f. update=.t. reset=.f. delciv=.f. $end     

 $statpt nstep=1000 opttol=1.0d-4 method=conic $end         

 $conic $end

 $zmat dlc=1 auto=1 $end

 $DATA

 optimization of ISC in CF2=CH-CH=CH2 by state-averaged CAS

CS

 

 H           1.0   0.7702498843  -0.0264069111   0.0000000000

 S   3

  1        13.01000000         0.1968500000E-01

  2        1.962000000         0.1379770000   

  3       0.4446000000         0.4781480000   

 S   1

  1       0.1220000000          1.000000000   

 P   1

  1       0.7270000000          1.000000000   

 

 C           6.0  -0.3178094497  -0.0863191256   0.0000000000

 S   8

  1        6665.000000         0.6920000000E-03

  2        1000.000000         0.5329000000E-02

  3        228.0000000         0.2707700000E-01

  4        64.71000000         0.1017180000   

  5        21.06000000         0.2747400000   

  6        7.495000000         0.4485640000   

  7        2.797000000         0.2850740000   

  8       0.5215000000         0.1520400000E-01

 S   8

  1        6665.000000        -0.1460000000E-03

  2        1000.000000        -0.1154000000E-02

  3        228.0000000        -0.5725000000E-02

  4        64.71000000        -0.2331200000E-01

  5        21.06000000        -0.6395500000E-01

  6        7.495000000        -0.1499810000   

  7        2.797000000        -0.1272620000   

  8       0.5215000000         0.5445290000   

 S   1

  1       0.1596000000          1.000000000   

 P   3

  1        9.439000000         0.3810900000E-01

  2        2.002000000         0.2094800000    

  3       0.5456000000         0.5085570000   

 P   1

  1       0.1517000000          1.000000000   

 D   1

  1       0.5500000000          1.000000000   

 

 C           6.0  -0.8695993539  -1.3160814708   0.0000000000

 S   8

  1        6665.000000         0.6920000000E-03

  2        1000.000000         0.5329000000E-02

  3        228.0000000         0.2707700000E-01

  4        64.71000000         0.1017180000   

  5        21.06000000         0.2747400000   

  6        7.495000000         0.4485640000   

  7        2.797000000         0.2850740000   

  8       0.5215000000         0.1520400000E-01

 S   8

  1        6665.000000        -0.1460000000E-03

  2        1000.000000        -0.1154000000E-02

  3        228.0000000        -0.5725000000E-02

  4        64.71000000        -0.2331200000E-01

  5        21.06000000        -0.6395500000E-01

  6        7.495000000        -0.1499810000   

  7        2.797000000        -0.1272620000   

  8       0.5215000000         0.5445290000   

 S   1

  1       0.1596000000          1.000000000   

 P   3

  1        9.439000000         0.3810900000E-01

  2        2.002000000         0.2094800000   

  3       0.5456000000         0.5085570000   

 P   1

  1       0.1517000000          1.000000000   

 D   1

  1       0.5500000000          1.000000000   

 

 C           6.0  -1.0185012979   1.1687688440   0.0000000000

 S   8

  1        6665.000000         0.6920000000E-03

  2        1000.000000         0.5329000000E-02

  3        228.0000000         0.2707700000E-01

  4        64.71000000         0.1017180000   

  5        21.06000000         0.2747400000   

  6        7.495000000         0.4485640000   

  7        2.797000000         0.2850740000   

  8       0.5215000000         0.1520400000E-01

 S   8

  1        6665.000000        -0.1460000000E-03

  2        1000.000000        -0.1154000000E-02

  3        228.0000000        -0.5725000000E-02

  4        64.71000000        -0.2331200000E-01

  5        21.06000000        -0.6395500000E-01

  6        7.495000000        -0.1499810000   

  7        2.797000000        -0.1272620000   

  8       0.5215000000         0.5445290000   

 S   1

  1       0.1596000000          1.000000000   

 P   3

  1        9.439000000         0.3810900000E-01

  2        2.002000000         0.2094800000   

  3       0.5456000000         0.5085570000   

 P   1

  1       0.1517000000          1.000000000   

 D   1

  1       0.5500000000          1.000000000   

 

 H           1.0  -0.2496902066  -2.2018838869   0.0000000000

 S   3

  1        13.01000000         0.1968500000E-01

  2        1.962000000         0.1379770000   

  3       0.4446000000         0.4781480000   

 S   1

  1       0.1220000000          1.000000000   

 P   1

  1       0.7270000000          1.000000000   

 

 C           6.0  -0.3881093034   2.4164000813   0.0000000000

 S   8

  1        6665.000000         0.6920000000E-03

  2        1000.000000         0.5329000000E-02

  3        228.0000000         0.2707700000E-01

  4        64.71000000         0.1017180000   

  5        21.06000000         0.2747400000   

  6        7.495000000         0.4485640000   

  7        2.797000000         0.2850740000   

  8       0.5215000000         0.1520400000E-01

 S   8

  1        6665.000000        -0.1460000000E-03

  2        1000.000000        -0.1154000000E-02

  3        228.0000000        -0.5725000000E-02

  4        64.71000000        -0.2331200000E-01

  5        21.06000000        -0.6395500000E-01

  6        7.495000000        -0.1499810000   

  7        2.797000000        -0.1272620000   

  8       0.5215000000         0.5445290000   

 S   1

  1       0.1596000000          1.000000000   

 P   3

  1        9.439000000         0.3810900000E-01

  2        2.002000000         0.2094800000   

  3       0.5456000000         0.5085570000   

 P   1

  1       0.1517000000          1.000000000   

 D   1

  1       0.5500000000          1.000000000   

 

 H           1.0  -1.9445187301  -1.4525604073   0.0000000000

 S   3

  1        13.01000000         0.1968500000E-01

  2        1.962000000         0.1379770000   

  3       0.4446000000         0.4781480000   

 S   1

  1       0.1220000000          1.000000000   

 P   1

  1       0.7270000000          1.000000000   

 

 H           1.0  -2.1013245134   1.2299991265   0.0000000000

 S   3

  1        13.01000000         0.1968500000E-01

  2        1.962000000         0.1379770000   

  3       0.4446000000         0.4781480000   

 S   1

  1       0.1220000000          1.000000000   

 P   1

  1       0.7270000000          1.000000000   

 

 F           9.0  -0.0719399830   3.1332504676   1.0527749090

 S   8

  1        14710.00000         0.7210000000E-03

  2        2207.000000         0.5553000000E-02

  3        502.8000000         0.2826700000E-01

  4        142.6000000         0.1064440000   

  5        46.47000000         0.2868140000   

  6        16.70000000         0.4486410000   

  7        6.356000000         0.2647610000   

  8        1.316000000         0.1533300000E-01

 S   8

  1        14710.00000        -0.1650000000E-03

  2        2207.000000        -0.1308000000E-02

  3        502.8000000        -0.6495000000E-02

  4        142.6000000        -0.2669100000E-01

  5        46.47000000        -0.7369000000E-01

  6        16.70000000        -0.1707760000   

  7        6.356000000        -0.1123270000   

  8        1.316000000         0.5628140000   

 S   1

  1       0.3897000000          1.000000000   

 P   3

  1        22.67000000         0.4487800000E-01

  2        4.977000000         0.2357180000   

  3        1.347000000         0.5085210000   

 P   1

  1       0.3471000000          1.000000000   

 D   1

  1        1.640000000          1.000000000   

 

 $END

--- NATURAL ORBITALS OF MCSCF --- GENERATED AT 13:56:45 LT  29-FEB-2008

 optimization of first CI in CF2=CH-CH=CH2 by state-averaged CAS

E(MCSCF)=     -352.5912216598, 22 ITERS, E(NUC)=  215.2380986153

 $VEC

 1  1 3.94108091E-05-1.36899842E-04 1.53350671E-05 6.92804971E-05 0.00000000E+00

 1  2 1.92770369E-05 3.62947163E-05-2.23328986E-04 1.09488689E-05 1.33918231E-05

 1  3 0.00000000E+00 1.72644900E-04-1.07495094E-03 0.00000000E+00 7.36293651E-06

 1  4-1.93874994E-05 1.20245629E-05 2.63917967E-05 0.00000000E+00 0.00000000E+00

 1  5 3.19574911E-05 5.42264019E-05-7.37465162E-04 6.37842521E-06-1.07304036E-05

 1  6 0.00000000E+00-3.70436277E-04-2.37571546E-04 0.00000000E+00-1.76600598E-05

 1  7-9.60079326E-06 2.72608530E-05-7.55281671E-06 0.00000000E+00 0.00000000E+00

 1  8-9.67536584E-05-2.03363389E-04 1.43814497E-03-7.00491791E-05-1.21969283E-04

 1  9 0.00000000E+00 5.68092261E-04-4.16075356E-04 0.00000000E+00 2.38670423E-05

 1 10-4.57693966E-05 2.19023543E-05-8.01100726E-05 0.00000000E+00 0.00000000E+00

 1 11 5.76649296E-05 2.94221434E-05 6.81668680E-06 5.94672176E-06 0.00000000E+00

 1 12 1.26830791E-04 2.79952235E-04 2.22791422E-03 1.32920883E-04 3.32207227E-04

 1 13 0.00000000E+00 6.25530182E-04 1.55886736E-03 0.00000000E+00-1.63570860E-04

 1 14 8.77118821E-05 7.58589781E-05 1.96593025E-04 0.00000000E+00 0.00000000E+00

 1 15 1.68657987E-06-1.34240922E-04 7.12103757E-06 1.16800087E-05 0.00000000E+00

 1 16 3.26044818E-05 1.01525430E-04-2.13729272E-05 1.67733296E-05 0.00000000E+00

 1 17 7.07834125E-01 8.80434641E-04-3.69624410E-03-1.46244408E-04-3.67551714E-04

 1 18 6.03302491E-04 1.71470878E-04 3.64293719E-04-9.89627010E-04 3.83166019E-05

 1 19 3.50642366E-05-7.33808385E-05 8.27817536E-07 4.59882776E-05 1.05723768E-04

 1 20 7.07834125E-01 8.80434641E-04-3.69624410E-03-1.46244408E-04-3.67551714E-04

 1 21-6.03302491E-04 1.71470878E-04 3.64293719E-04 9.89627010E-04 3.83166019E-05

 1 22 3.50642366E-05-7.33808385E-05 8.27817536E-07-4.59882776E-05-1.05723768E-04

 2  1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.36794388E-05

 2  2 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 2  3-1.52563709E-06 0.00000000E+00 0.00000000E+00-6.48036620E-05 0.00000000E+00

 2  4 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.28457203E-05 2.66079495E-06

 2  5 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 2  6 3.15004727E-05 0.00000000E+00 0.00000000E+00-3.21389116E-05 0.00000000E+00

 2  7 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.48711871E-05-8.11064292E-06

 2  8 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 2  9-3.53229101E-05 0.00000000E+00 0.00000000E+00 6.83105411E-04 0.00000000E+00

 2 10 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.16844474E-04 2.13306345E-04

 2 11 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 4.28737175E-06

 2 12 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 2 13-1.31782864E-03 0.00000000E+00 0.00000000E+00-1.99072650E-03 0.00000000E+00

 2 14 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.98566733E-06-4.20403684E-05

 2 15 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 8.87159022E-06

 2 16 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00-3.07817688E-05

 2 17 7.07800399E-01 9.46139750E-04-3.32226393E-03-5.06239289E-04-1.09014652E-03

 2 18 1.28944608E-03 8.57718115E-05 2.39969313E-04-2.06747714E-04 7.81643142E-06

 2 19-5.22329978E-06-2.59313164E-06-8.38256630E-06 2.12156746E-05 4.52245242E-05

 2 20-7.07800399E-01-9.46139750E-04 3.32226393E-03 5.06239289E-04 1.09014652E-03

 2 21 1.28944608E-03-8.57718115E-05-2.39969313E-04-2.06747714E-04-7.81643142E-06

 2 22 5.22329978E-06 2.59313164E-06 8.38256630E-06 2.12156746E-05 4.52245242E-05

 3  1 1.79285524E-05 3.68477748E-04-6.64403562E-05-1.26816283E-04 0.00000000E+00

 3  2-8.00164585E-05-1.43054272E-04-2.60779627E-04 8.83810691E-06-9.28321470E-05

 3  3 0.00000000E+00-3.33596963E-04 1.88801940E-03 0.00000000E+00-9.89187327E-05

 3  4 1.17799373E-04-1.88806398E-05-2.70410237E-05 0.00000000E+00 0.00000000E+00

 3  5-1.21616120E-04-2.08387605E-04 1.76321208E-03 6.57472525E-06-8.43175450E-06

 3  6 0.00000000E+00 5.15563894E-04 9.12449675E-04 0.00000000E+00 2.17945234E-05

 3  7 7.60809310E-05-9.78754544E-05-2.49729018E-05 0.00000000E+00 0.00000000E+00

 3  8 2.09590509E-03-3.04654893E-04 3.49442756E-03 7.32079394E-05 1.92911229E-04

 3  9 0.00000000E+00 8.12994206E-04 3.32654136E-03 0.00000000E+00-1.19734778E-04

 3 10 3.00924684E-04-1.81189906E-04 2.97085262E-04 0.00000000E+00 0.00000000E+00

 3 11-1.04035039E-04 1.27726854E-04-5.17490664E-06-3.35341325E-05 0.00000000E+00

 3 12 1.00169467E+00 3.18350296E-03-7.28721641E-03 7.45013644E-04 2.09901000E-03

 3 13 0.00000000E+00 1.56472190E-03 3.13581930E-03 0.00000000E+00-6.32763307E-04

 3 14 9.32205918E-05 5.39542715E-04 3.91699884E-04 0.00000000E+00 0.00000000E+00

 3 15-7.58640645E-05 2.52168961E-04 8.46998987E-07-3.43948079E-05 0.00000000E+00

 3 16-2.17480551E-05 3.29924677E-04 1.25444070E-04-2.50432976E-05 0.00000000E+00

 3 17 1.94778457E-04 6.42346188E-04-1.50727239E-03-2.76780372E-04-1.94405684E-05

 3 18 2.96990258E-04-4.76243354E-05 3.50692141E-04-1.25934865E-03 1.10573576E-04

 3 19 4.21735983E-05-1.52747174E-04-4.08942991E-05 4.40354837E-05 1.17172789E-04

 3 20 1.94778457E-04 6.42346188E-04-1.50727239E-03-2.76780372E-04-1.94405684E-05

 3 21-2.96990258E-04-4.76243354E-05 3.50692141E-04 1.25934865E-03 1.10573576E-04

 3 22 4.21735983E-05-1.52747174E-04-4.08942991E-05-4.40354837E-05-1.17172789E-04

 4  1 3.41655260E-04-9.32530878E-04-4.82373673E-04-8.90013819E-06 0.00000000E+00

 4  2-1.00079397E+00-4.32048268E-03 9.16602643E-03-8.44189614E-05-2.58194776E-04

 4  3 0.00000000E+00-1.37124717E-03 2.06802504E-05 0.00000000E+00-1.65836720E-04

 4  4-4.26659439E-04 5.92496159E-04-8.01160624E-05 0.00000000E+00 0.00000000E+00

 4  5-2.95071021E-02 3.36729732E-04-2.25649839E-03 1.87342874E-05 7.26288149E-05

 4  6 0.00000000E+00-7.22467791E-05-1.35051118E-03 0.00000000E+00 2.47911476E-04

 4  7 4.66451907E-05-2.94556667E-04 4.70902145E-05 0.00000000E+00 0.00000000E+00

 4  8-4.06411710E-02 5.80092899E-05-8.88332577E-04 3.75026487E-05-1.48568747E-04

 4  9 0.00000000E+00-7.34652037E-05 1.80409010E-03 0.00000000E+00 1.80046551E-04

 4 10 1.96904236E-04-3.76950787E-04-3.45541257E-04 0.00000000E+00 0.00000000E+00

 4 11-8.04249721E-05-3.60373817E-04 1.63266031E-05-1.13091969E-04 0.00000000E+00

 4 12 2.21390630E-04 2.25964965E-04-1.42392405E-03-9.36897972E-05-5.19320634E-06

 4 13 0.00000000E+00 4.07905458E-04 6.44244651E-04 0.00000000E+00-3.33643302E-05

 4 14-2.15048426E-05 5.48691728E-05-1.92044869E-05 0.00000000E+00 0.00000000E+00

 4 15-4.73610351E-05-5.50189850E-05-1.62872885E-04-4.65478005E-05 0.00000000E+00

 4 16-6.22194002E-05-1.31010699E-04-1.42160607E-04-9.87959720E-05 0.00000000E+00

 4 17-1.18864327E-08-4.91146297E-07-5.10260400E-06-1.00287174E-04 6.83104649E-05

 4 18-1.96976956E-06-1.20858801E-04-4.05983247E-05-1.02737366E-05 2.93461840E-06

 4 19-2.13830172E-06-7.96316684E-07-2.06195567E-06-4.06804887E-06 3.06662027E-06

 4 20-1.18864327E-08-4.91146297E-07-5.10260400E-06-1.00287174E-04 6.83104649E-05

 4 21 1.96976956E-06-1.20858801E-04-4.05983247E-05 1.02737366E-05 2.93461840E-06

 4 22-2.13830172E-06-7.96316684E-07-2.06195567E-06 4.06804887E-06-3.06662027E-06

 5  1 2.45856820E-05-5.69558225E-05 1.41487023E-04 6.83535563E-05 0.00000000E+00

 5  2 4.10458597E-02 4.81790357E-04-2.06548922E-03-8.86602279E-06 1.02445109E-04

 5  3 0.00000000E+00 2.80675013E-04-2.61198275E-03 0.00000000E+00 1.13263624E-04

 5  4 1.98154115E-04-3.11417739E-04-1.08598745E-04 0.00000000E+00 0.00000000E+00

 5  5 2.79553802E-03 1.06938288E-04-1.57972271E-03-1.13965229E-05-1.85415261E-05

 5  6 0.00000000E+00-6.08232928E-04-4.98707633E-04 0.00000000E+00-8.13524248E-05

 5  7-6.87889742E-05 1.50141399E-04-8.45249630E-06 0.00000000E+00 0.00000000E+00

 5  8-1.00139871E+00-4.74331097E-03 9.47883792E-03-2.86840933E-04-7.57262038E-04

 5  9 0.00000000E+00 1.72268662E-03-2.66128398E-04 0.00000000E+00-1.37929185E-04

 5 10-2.79066768E-04 4.16995954E-04-2.39714343E-05 0.00000000E+00 0.00000000E+00

 5 11 1.34579290E-04-1.23285101E-05-5.97072726E-06 2.74214483E-05 0.00000000E+00

 5 12 2.69115616E-03 2.46099815E-04-2.02617384E-03-1.81223607E-04-1.82044484E-04

 5 13 0.00000000E+00 4.12954685E-05 6.61724893E-04 0.00000000E+00 6.26762630E-05

 5 14 3.18917431E-04-3.81593694E-04 3.76274042E-04 0.00000000E+00 0.00000000E+00

 5 15 6.63529585E-05-1.17871902E-04 4.06436482E-05 4.36040810E-05 0.00000000E+00

 5 16 3.37323047E-04-8.77726577E-04 4.16319866E-04 3.58150480E-05 0.00000000E+00

 5 17 5.08222011E-05 1.16407936E-04-4.91765859E-05-6.62136852E-05-7.73383088E-05

 5 18 1.83295583E-05 3.42728195E-05 7.73464810E-05-2.31540273E-05 1.77648751E-05

 5 19 1.38514664E-05-3.16163415E-05 2.25021619E-06-1.48719376E-05-1.69481360E-05

 5 20 5.08222011E-05 1.16407936E-04-4.91765859E-05-6.62136852E-05-7.73383088E-05

 5 21-1.83295583E-05 3.42728195E-05 7.73464810E-05 2.31540273E-05 1.77648751E-05

 5 22 1.38514664E-05-3.16163415E-05 2.25021619E-06 1.48719376E-05 1.69481360E-05

 6  1 1.69212382E-04 1.27171464E-04-2.13242141E-04-4.71801256E-05 0.00000000E+00

 6  2-3.03984191E-02-4.56706803E-04 1.24590582E-03 5.76203897E-06 1.53731242E-04

 6  3 0.00000000E+00-1.19569285E-04-1.29133570E-03 0.00000000E+00-2.29439107E-04

 6  4-1.89717311E-04 4.19156418E-04-1.40038111E-06 0.00000000E+00 0.00000000E+00

 6  5 1.00164231E+00 4.37367334E-03-7.63971778E-03-2.96735981E-05-7.07571451E-05

 6  6 0.00000000E+00-2.91946581E-04-5.64244938E-04 0.00000000E+00 3.23599816E-04

 6  7 3.93692226E-04-7.17292042E-04 2.39462370E-05 0.00000000E+00 0.00000000E+00

 6  8 1.47349753E-03-7.29033334E-05 4.86072270E-04 3.38155801E-05-2.76341833E-05

 6  9 0.00000000E+00 6.85708840E-05-9.15656460E-04 0.00000000E+00-1.64888411E-05

 6 10 1.51831242E-04-1.35342401E-04-2.71537660E-05 0.00000000E+00 0.00000000E+00

 6 11-4.03964495E-04 1.20268102E-03 3.73934382E-04-4.32518048E-04 0.00000000E+00

 6 12-7.31833030E-05-1.19820232E-04 7.53503379E-04 4.11192244E-05 2.49642040E-05

 6 13 0.00000000E+00-1.89176782E-04-2.29971412E-04 0.00000000E+00-2.00985971E-05

 6 14-3.76768811E-06 2.38662852E-05-8.33459883E-06 0.00000000E+00 0.00000000E+00

 6 15-3.78943476E-04 1.18530950E-03-5.60973521E-04-8.79133324E-06 0.00000000E+00

 6 16-5.37422772E-05-4.31391513E-05-8.73008694E-06 7.58203613E-05 0.00000000E+00

 6 17 6.50517230E-06 1.15984520E-05-5.35343385E-05 1.26583557E-05-9.51844901E-06

 6 18-1.41949500E-06 4.12458846E-05 2.81170807E-05-1.30133367E-05-1.34498901E-08

 6 19-1.47452227E-06 1.48797216E-06-1.29988057E-06 3.81340927E-06 1.91587976E-06

 6 20 6.50517230E-06 1.15984520E-05-5.35343385E-05 1.26583557E-05-9.51844901E-06

 6 21 1.41949500E-06 4.12458846E-05 2.81170807E-05 1.30133367E-05-1.34498901E-08

 6 22-1.47452227E-06 1.48797216E-06-1.29988057E-06-3.81340927E-06-1.91587976E-06

 7  1-9.85415128E-05-1.26745007E-03 1.03991712E-04 9.02592803E-04 0.00000000E+00

 7  2 3.42056013E-04-5.09788076E-04 1.77868166E-03-2.17278542E-04 8.49965966E-06

 7  3 0.00000000E+00 1.07687226E-03-1.21651431E-02 0.00000000E+00 3.58214808E-05

 7  4 4.00972001E-04-4.36793482E-04 2.31979318E-04 0.00000000E+00 0.00000000E+00

 7  5 3.71951583E-04 4.27445821E-04-1.04173694E-02-8.82074015E-05-5.45639769E-04

 7  6 0.00000000E+00-4.72650997E-03-4.50721012E-03 0.00000000E+00-2.09218644E-04

 7  7-2.56244719E-04 4.65463363E-04-3.22495683E-04 0.00000000E+00 0.00000000E+00

 7  8-8.72182534E-05-1.50003326E-02 1.28621637E-03-5.42852255E-03-1.32050316E-02

 7  9 0.00000000E+00 3.59194012E-03-8.63271164E-03 0.00000000E+00-1.05598077E-04

 7 10-3.81588706E-03 3.92148514E-03-2.12332821E-03 0.00000000E+00 0.00000000E+00

 7 11 6.74421484E-04 9.67874455E-05 1.16868567E-04 1.30375474E-04 0.00000000E+00

 7 12 1.08377983E-02-1.56357464E-01-3.24473354E-02-3.45684984E-02-8.60895349E-02

 7 13 0.00000000E+00 8.77863878E-04 5.18855883E-04 0.00000000E+00 1.53090626E-02

 7 14 6.58017644E-04-1.59670802E-02-7.77568934E-03 0.00000000E+00 0.00000000E+00

 7 15 1.63098057E-04-2.13266755E-03 9.16357096E-05 8.16905325E-05 0.00000000E+00

 7 16-7.37079905E-04 1.57902234E-03 1.90209329E-04 5.56068694E-04 0.00000000E+00

 7 17 3.46123646E-03-3.21230030E-01-3.41785066E-01 1.81920433E-02 2.69643988E-02

 7 18-5.77583092E-02 8.75797728E-03 1.04842402E-02-2.75408633E-02 3.42939602E-03

 7 19 1.44781740E-03-4.87721342E-03-1.40956567E-03 2.78104259E-03 5.90749218E-03

 7 20 3.46123646E-03-3.21230030E-01-3.41785066E-01 1.81920433E-02 2.69643988E-02

 7 21 5.77583092E-02 8.75797728E-03 1.04842402E-02 2.75408633E-02 3.42939602E-03

 7 22 1.44781740E-03-4.87721342E-03-1.40956567E-03-2.78104259E-03-5.90749218E-03

 8  1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00-1.92365264E-04

 8  2 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 8  3-1.30275030E-03 0.00000000E+00 0.00000000E+00-1.56820781E-03 0.00000000E+00

 8  4 0.00000000E+00 0.00000000E+00 0.00000000E+00 6.14076169E-04 1.66896291E-04

 8  5 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 8  6-1.03201437E-03 0.00000000E+00 0.00000000E+00-7.21380469E-04 0.00000000E+00

 8  7 0.00000000E+00 0.00000000E+00 0.00000000E+00-5.10138168E-05 2.71331318E-05

 8  8 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 8  9 2.78360792E-03 0.00000000E+00 0.00000000E+00 2.29257199E-03 0.00000000E+00

 8 10 0.00000000E+00 0.00000000E+00 0.00000000E+00-1.31474273E-03-2.49173835E-03

 8 11 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.15599948E-05

 8 12 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 8 13-2.07411045E-02 0.00000000E+00 0.00000000E+00-4.79082817E-03 0.00000000E+00

 8 14 0.00000000E+00 0.00000000E+00 0.00000000E+00-9.97233647E-03-2.29025865E-02

 8 15 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00-3.09408452E-05

 8 16 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00-6.50803710E-05

 8 17-1.74097094E-03 3.53885421E-01 3.89156391E-01 2.54655087E-02 6.33379245E-02

 8 18-5.27913978E-02 1.86299621E-02 4.62938400E-02-3.97028722E-02-6.28565130E-05

 8 19-3.95795233E-04 4.58651746E-04-1.49518058E-04-5.62818699E-05 1.05788150E-05

 8 20 1.74097094E-03-3.53885421E-01-3.89156391E-01-2.54655087E-02-6.33379245E-02

 8 21-5.27913978E-02-1.86299621E-02-4.62938400E-02-3.97028722E-02 6.28565130E-05

 8 22 3.95795233E-04-4.58651746E-04 1.49518058E-04-5.62818699E-05 1.05788150E-05

 9  1 8.84608336E-02 1.06417120E-02-1.37685425E-02-1.71323177E-05 0.00000000E+00

 9  2-1.77399261E-02 2.36716332E-01 1.67147314E-01-4.65630326E-02 1.90245790E-02

 9  3 0.00000000E+00-3.91306443E-02-1.52389935E-02 0.00000000E+00-6.97779805E-03

 9  4 1.95817687E-02-1.26039706E-02-2.58346295E-03 0.00000000E+00 0.00000000E+00

 9  5-1.08228638E-02 1.40410893E-01 4.88766917E-02 2.30266879E-02 6.66508654E-02

 9  6 0.00000000E+00 8.70897688E-03-5.71876808E-03 0.00000000E+00-6.66156547E-03

 9  7 9.79320356E-03-3.13163809E-03 5.32306133E-03 0.00000000E+00 0.00000000E+00

 9  8-1.85569872E-02 2.69157163E-01 1.49334921E-01 5.78811834E-02 9.81329648E-03

 9  9 0.00000000E+00 3.24559157E-02 2.14230030E-03 0.00000000E+00-4.56587145E-03

 9 10 1.74056512E-02-1.28397798E-02-6.60167029E-04 0.00000000E+00 0.00000000E+00

 9 11 4.32157109E-02-2.60882748E-03-2.71993437E-03 7.54002623E-03 0.00000000E+00

 9 12 1.73727820E-03 2.06241788E-01 9.58433464E-02-5.26648300E-02-1.11057183E-01

 9 13 0.00000000E+00-2.41676239E-02-4.14485777E-02 0.00000000E+00-5.44474667E-03

 9 14 9.72183736E-03-4.27709069E-03 7.63863907E-03 0.00000000E+00 0.00000000E+00

 9 15 5.18987932E-02 5.45926182E-03 9.02432265E-03 2.28694669E-03 0.00000000E+00

 9 16 9.71121536E-02 1.04987558E-02 1.47198701E-02-4.17625190E-04 0.00000000E+00

 9 17-9.31914116E-04-6.15019265E-02-8.04602725E-02-3.50424644E-02-5.98218354E-02

 9 18 4.31582088E-02-1.83674495E-02-3.02373525E-02 2.74607338E-02-1.23079717E-03

 9 19 9.33130920E-04 2.97666247E-04 1.52273868E-03-1.76776289E-03-3.61347215E-03

 9 20-9.31914116E-04-6.15019265E-02-8.04602725E-02-3.50424644E-02-5.98218354E-02

 9 21-4.31582088E-02-1.83674495E-02-3.02373525E-02-2.74607338E-02-1.23079717E-03

 9 22 9.33130920E-04 2.97666247E-04 1.52273868E-03 1.76776289E-03 3.61347215E-03

10  1-7.63256773E-02-1.79581658E-02 1.11067970E-02 2.00957101E-03 0.00000000E+00

10  2 1.44334651E-02-1.97669306E-01-1.57739225E-01 3.17698444E-02 1.17576319E-01

10  3 0.00000000E+00 3.54192996E-02 3.25990977E-02 0.00000000E+00 5.03332131E-04

10  4-8.36906140E-03 7.86572927E-03-1.10233544E-02 0.00000000E+00 0.00000000E+00

10  5 1.83356170E-02-2.37162777E-01-1.33181820E-01-3.77233002E-02-7.62428706E-02

10  6 0.00000000E+00-9.21703876E-03-6.04863009E-03 0.00000000E+00 3.10100796E-03

10  7-8.99368776E-03 5.89267980E-03-1.04926655E-02 0.00000000E+00 0.00000000E+00

10  8-4.12761961E-03 1.05739090E-01 2.24273129E-02 1.88762610E-02 1.36178541E-01

10  9 0.00000000E+00-2.03322966E-02-1.65495869E-02 0.00000000E+00 1.33744822E-03

10 10 2.54551608E-03-3.88296430E-03 1.26414716E-02 0.00000000E+00 0.00000000E+00

10 11-9.57794509E-02-5.50651607E-03 7.68966293E-03-1.35879874E-02 0.00000000E+00

10 12 4.12279257E-03 2.58473639E-01 1.71383465E-01-3.78023192E-02-5.38195380E-02

10 13 0.00000000E+00-3.98299428E-02-6.91181357E-02 0.00000000E+00-5.50561712E-03

10 14 3.97355915E-03 1.53205797E-03 1.09217048E-02 0.00000000E+00 0.00000000E+00

10 15-9.19076629E-02-9.16992394E-03-1.40573095E-02-3.38777755E-03 0.00000000E+00

10 16 4.00377681E-02-6.66857642E-03 5.77751802E-03 3.96853969E-03 0.00000000E+00

10 17-1.75205827E-03-8.52610796E-02-1.15390384E-01-4.85438536E-02-9.92288343E-02

10 18 1.11019643E-01-2.46470542E-02-5.24457782E-02 6.85642609E-02-2.70392182E-03

10 19 1.05659138E-04 2.59826268E-03 2.05505411E-03-2.97137031E-03-6.09744882E-03

10 20-1.75205827E-03-8.52610796E-02-1.15390384E-01-4.85438536E-02-9.92288343E-02

10 21-1.11019643E-01-2.46470542E-02-5.24457782E-02-6.85642609E-02-2.70392182E-03

10 22 1.05659138E-04 2.59826268E-03 2.05505411E-03 2.97137031E-03 6.09744882E-03

11  1 4.69071521E-03 7.69423367E-03-6.94491517E-04 3.20849059E-03 0.00000000E+00

11  2 1.02330838E-03 3.54651344E-03-2.37997878E-02-4.59422404E-03 1.28269613E-01

11  3 0.00000000E+00-1.01806068E-02 1.21261766E-02 0.00000000E+00 8.48340250E-04

11  4-2.92660601E-06-8.45413644E-04-1.32024574E-02 0.00000000E+00 0.00000000E+00

11  5 9.88630181E-03-1.42688641E-01-1.09581752E-01-1.51760411E-02-6.25840869E-03

11  6 0.00000000E+00-7.52928483E-03 3.47917074E-03 0.00000000E+00-1.71172123E-03

11  7-3.53654766E-04 2.06537599E-03-3.56617818E-03 0.00000000E+00 0.00000000E+00

11  8-1.62137628E-02 1.59706512E-01 2.29666230E-01 1.00949303E-02-1.09196750E-03

11  9 0.00000000E+00 5.35377274E-02 3.45773663E-02 0.00000000E+00 1.24176030E-03

11 10 3.79045991E-03-5.03222021E-03 2.13553055E-04 0.00000000E+00 0.00000000E+00

11 11-7.81519101E-02-6.69627134E-03 6.35468910E-03-8.90510518E-03 0.00000000E+00

11 12-4.76751832E-03-1.48534486E-01-1.09746443E-01-1.02822440E-01-2.34270994E-01

11 13 0.00000000E+00 1.18786540E-02 3.52841928E-02 0.00000000E+00 1.28167225E-02

11 14-7.15252324E-03-5.66419931E-03-9.08482151E-03 0.00000000E+00 0.00000000E+00

11 15-6.70629648E-02-9.05630482E-03-8.98356597E-03-1.60345213E-03 0.00000000E+00

11 16 8.47258294E-02 2.72008030E-02 1.11918239E-02-1.13301402E-03 0.00000000E+00

11 17 2.20316843E-03 9.68851468E-02 1.39596730E-01 3.25211509E-02 5.75412042E-02

11 18-3.09270579E-01 1.75369174E-02 2.73191653E-02-1.89913429E-01 5.72587384E-03

11 19 4.65665017E-03-1.03825240E-02-8.18662051E-04 3.67405821E-03 7.93295215E-03

11 20 2.20316843E-03 9.68851468E-02 1.39596730E-01 3.25211509E-02 5.75412042E-02

11 21 3.09270579E-01 1.75369174E-02 2.73191653E-02 1.89913429E-01 5.72587384E-03

11 22 4.65665017E-03-1.03825240E-02-8.18662051E-04-3.67405821E-03-7.93295215E-03

12  1-1.21084098E-01-2.39360092E-02 1.30896984E-02-3.21915570E-03 0.00000000E+00

12  2 6.80732699E-03-1.68479636E-01-1.15882215E-01-2.72211695E-02-1.43007980E-01

12  3 0.00000000E+00-6.36376727E-03-5.34758436E-02 0.00000000E+00-7.08201112E-03

12  4 3.65401677E-03 3.42799435E-03 1.36940759E-02 0.00000000E+00 0.00000000E+00

12  5-1.32072404E-02 1.90523485E-01 1.47301355E-01-3.11126380E-02-9.38795947E-02

12  6 0.00000000E+00-2.04371236E-02-3.07964182E-02 0.00000000E+00 5.20183973E-03

12  7-4.03551267E-03-1.16632707E-03-5.45340999E-03 0.00000000E+00 0.00000000E+00

12  8 4.57644209E-03-5.13004303E-02-6.46985225E-02-2.17082002E-02 1.91191357E-01

12  9 0.00000000E+00-1.77422432E-02 5.09945953E-02 0.00000000E+00-5.91251133E-03

12 10 4.29539000E-03 1.61712133E-03 2.21866690E-02 0.00000000E+00 0.00000000E+00

12 11 1.53985288E-01 3.17245931E-02-1.17664397E-02 1.29647823E-02 0.00000000E+00

12 12 1.19698254E-03 8.16913135E-02 7.23101571E-02-7.41998549E-02-1.31146871E-01

12 13 0.00000000E+00-2.18708353E-02-5.58316895E-02 0.00000000E+00 2.44100284E-03

12 14-2.55574194E-02 2.31164165E-02-9.30033036E-03 0.00000000E+00 0.00000000E+00

12 15 1.41104810E-01 2.49284487E-02 1.69524363E-02 6.71602308E-05 0.00000000E+00

12 16-1.28113321E-02-4.42186443E-03-4.14549521E-04 4.93021345E-03 0.00000000E+00

12 17-5.10083472E-04 5.87096228E-03 1.52476607E-02-6.72630350E-02-1.88145088E-01

12 18-1.77647168E-01-4.22295579E-02-1.20754668E-01-1.11423670E-01 1.44541498E-03

12 19 4.78439273E-03-6.22980771E-03 2.14709468E-03-7.23961866E-04-1.89806009E-03

12 20-5.10083472E-04 5.87096228E-03 1.52476607E-02-6.72630350E-02-1.88145088E-01

12 21 1.77647168E-01-4.22295579E-02-1.20754668E-01 1.11423670E-01 1.44541498E-03

12 22 4.78439273E-03-6.22980771E-03 2.14709468E-03 7.23961866E-04 1.89806009E-03

13  1 1.99593953E-01 4.19981895E-02-1.93372184E-02-2.48672558E-03 0.00000000E+00

13  2-5.55987196E-03 1.47949689E-01 9.08231539E-02 1.54149717E-01-1.90145420E-02

13  3 0.00000000E+00 6.46627644E-02 2.54033219E-02 0.00000000E+00 3.07823739E-03

13  4-4.03782437E-03 9.59586978E-04 4.88674233E-03 0.00000000E+00 0.00000000E+00

13  5 4.55489149E-03-9.19180329E-02-4.03082689E-02 1.06517663E-01 7.53835147E-02

13  6 0.00000000E+00 4.58998179E-02 3.96607924E-02 0.00000000E+00-1.46051771E-03

13  7 2.67266066E-03-1.21214295E-03 9.76962618E-03 0.00000000E+00 0.00000000E+00

13  8 1.04494698E-02-1.39579994E-01-1.49494378E-01 1.43872762E-01-1.82314198E-02

13  9 0.00000000E+00 4.25221825E-02-1.24786999E-02 0.00000000E+00-1.86845250E-03

13 10 3.48027274E-03-1.61182024E-03-1.62902728E-03 0.00000000E+00 0.00000000E+00

13 11-6.01910084E-02-1.60340870E-02 5.90824401E-03-4.32886451E-03 0.00000000E+00

13 12 2.94551864E-03 3.47140066E-02 3.75246344E-03 2.93058941E-02 2.23502975E-02

13 13 0.00000000E+00-1.01376349E-02-2.67161249E-02 0.00000000E+00-1.27402317E-02

13 14-2.11815214E-02 3.39217531E-02-2.15973350E-02 0.00000000E+00 0.00000000E+00

13 15-1.40106809E-01-2.95768760E-02-1.41103506E-02-2.72894841E-04 0.00000000E+00

13 16-1.96376845E-01-4.38351059E-02-1.94737000E-02-5.29587365E-04 0.00000000E+00

13 17-1.14452059E-03-6.13258320E-03 1.13605656E-04-1.39441655E-01-2.21572146E-01

13 18-1.00096813E-01-9.45391521E-02-1.45788197E-01-6.47539035E-02-9.82742753E-05

13 19 3.09305633E-03-2.99478205E-03 2.34932422E-03-2.32764812E-03-3.93290712E-03

13 20-1.14452059E-03-6.13258320E-03 1.13605656E-04-1.39441655E-01-2.21572146E-01

13 21 1.00096813E-01-9.45391521E-02-1.45788197E-01 6.47539035E-02-9.82742753E-05

13 22 3.09305633E-03-2.99478205E-03 2.34932422E-03 2.32764812E-03 3.93290712E-03

14  1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00-3.53101825E-04

14  2 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

14  3 4.43486039E-05 0.00000000E+00 0.00000000E+00-4.06759122E-04 0.00000000E+00

14  4 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.08581518E-03-1.66315081E-04

14  5 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

14  6-1.31441698E-03 0.00000000E+00 0.00000000E+00-2.10510772E-03 0.00000000E+00

14  7 0.00000000E+00 0.00000000E+00 0.00000000E+00-1.49728465E-04-2.11754300E-05

14  8 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

14  9 9.87087563E-03 0.00000000E+00 0.00000000E+00 7.00481017E-03 0.00000000E+00

14 10 0.00000000E+00 0.00000000E+00 0.00000000E+00-3.02075894E-03-6.17603264E-03

14 11 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 5.54517287E-05

14 12 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

14 13-2.34790836E-02 0.00000000E+00 0.00000000E+00 3.44652326E-02 0.00000000E+00

14 14 0.00000000E+00 0.00000000E+00 0.00000000E+00-1.67594981E-02-3.67785077E-02

14 15 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00-3.45610410E-05

14 16 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00-5.36093545E-04

14 17-1.39474833E-03 3.46831251E-03 1.22730081E-02 1.29798693E-01 2.88187147E-01

14 18 3.78396151E-01 8.80253720E-02 1.95726707E-01 2.57251516E-01-2.35750676E-03

14 19-7.21294784E-03 9.57045459E-03-3.11414785E-03 6.18838350E-05 1.00332790E-04

14 20 1.39474833E-03-3.46831251E-03-1.22730081E-02-1.29798693E-01-2.88187147E-01

14 21 3.78396151E-01-8.80253720E-02-1.95726707E-01 2.57251516E-01 2.35750676E-03

14 22 7.21294784E-03-9.57045459E-03 3.11414785E-03 6.18838350E-05 1.00332790E-04

15  1-1.60740608E-01-3.83454457E-02 1.21546358E-02 3.82532660E-04 0.00000000E+00

15  2 2.47861656E-04 3.28632821E-02 4.22339254E-02-2.30310148E-01 2.93738624E-02

15  3 0.00000000E+00-9.89694047E-02 3.58256518E-02 0.00000000E+00 4.19457185E-03

15  4-2.82380208E-03-1.37076977E-03-5.54769938E-03 0.00000000E+00 0.00000000E+00

15  5 3.87439971E-03-2.25055211E-02-3.15616651E-02-1.21794066E-01 1.59565199E-01

15  6 0.00000000E+00-2.98490742E-02 4.55384192E-02 0.00000000E+00-4.88600142E-03

15  7 4.80017077E-03 8.58306475E-05 4.10259793E-03 0.00000000E+00 0.00000000E+00

15  8 3.31224784E-03 2.33304237E-02-4.33498096E-02-1.29408549E-01-2.30254626E-01

15  9 0.00000000E+00-7.14053284E-02-7.97497325E-02 0.00000000E+00-1.31396763E-04

15 10-2.50001254E-03 2.63140930E-03-1.34344442E-02 0.00000000E+00 0.00000000E+00

15 11-1.73979804E-01-5.57264496E-02 7.62970593E-03-1.05216297E-02 0.00000000E+00

15 12-3.25853106E-03-8.90223615E-02-7.37865393E-02 9.21860946E-02 2.07272840E-01

15 13 0.00000000E+00-1.37201320E-03 2.86509415E-03 0.00000000E+00-9.78931936E-03

15 14-2.50722058E-02 3.48615252E-02-9.13878599E-03 0.00000000E+00 0.00000000E+00

15 15 6.23987850E-02 2.03582976E-02 3.04897014E-03 4.90607915E-03 0.00000000E+00

15 16 1.07624574E-01 2.69502844E-02 7.64874539E-03-6.00300191E-03 0.00000000E+00

15 17-9.37437938E-04-5.59284347E-04 8.80033196E-03-9.53367592E-02-2.30630022E-01

15 18-2.41301599E-02-6.39918968E-02-1.58104768E-01-1.77732523E-02-1.32829794E-03

15 19 9.62703265E-04 3.65594672E-04 1.38539429E-03-1.72360902E-03-4.24187927E-03

15 20-9.37437938E-04-5.59284347E-04 8.80033196E-03-9.53367592E-02-2.30630022E-01

15 21 2.41301599E-02-6.39918968E-02-1.58104768E-01 1.77732523E-02-1.32829794E-03

15 22 9.62703265E-04 3.65594672E-04 1.38539429E-03 1.72360902E-03 4.24187927E-03

16  1 2.64889981E-02 1.34541574E-02-2.42013073E-03-6.44618697E-03 0.00000000E+00

16  2-1.93211002E-03 1.07562191E-02 2.16072632E-02 3.86872422E-02-2.92813856E-01

16  3 0.00000000E+00 7.56275475E-03-9.13556533E-02 0.00000000E+00 9.10658918E-04

16  4-3.26451807E-03 2.35385915E-03 3.85646806E-03 0.00000000E+00 0.00000000E+00

16  5 3.51001868E-03-1.14674700E-02-3.52340189E-02 2.64884869E-01 1.50024282E-01

16  6 0.00000000E+00 1.00134936E-01 5.67803690E-02 0.00000000E+00-8.62887340E-03

16  7 9.62133193E-03-9.92458536E-04 1.08168180E-02 0.00000000E+00 0.00000000E+00

16  8-3.25870051E-03-4.88565547E-03-7.63428363E-03-2.86228711E-01 1.28781380E-01

16  9 0.00000000E+00-1.31811802E-01 6.95736831E-02 0.00000000E+00 3.65980054E-03

16 10-6.16980017E-03 2.50999963E-03 4.98613740E-03 0.00000000E+00 0.00000000E+00

16 11 1.24057822E-02 3.15637688E-03 5.62705728E-03 5.20788434E-03 0.00000000E+00

16 12-1.66172112E-03-1.42217295E-02-1.05485710E-02-3.09248431E-02 3.89608035E-02

16 13 0.00000000E+00 1.13045204E-04 1.00004173E-03 0.00000000E+00 1.48344190E-02

16 14-1.64872417E-02 1.65282267E-03 1.26512673E-02 0.00000000E+00 0.00000000E+00

16 15-2.45589950E-01-8.05787464E-02-1.77185859E-02 1.28488910E-03 0.00000000E+00

16 16 2.47027124E-01 7.49164453E-02 1.86665730E-02 2.79114631E-03 0.00000000E+00

16 17 6.41127217E-05 4.40201943E-04 5.21004976E-04 7.01442246E-02-4.69091168E-02

16 18-9.77319529E-04 5.08610595E-02-3.28963322E-02-7.14055924E-04-3.34863550E-04

16 19 1.90543437E-04 1.44320113E-04-2.43460532E-04 8.96566102E-04-6.46510877E-04

16 20 6.41127217E-05 4.40201943E-04 5.21004976E-04 7.01442246E-02-4.69091168E-02

16 21 9.77319529E-04 5.08610595E-02-3.28963322E-02 7.14055924E-04-3.34863550E-04

16 22 1.90543437E-04 1.44320113E-04-2.43460532E-04-8.96566102E-04 6.46510877E-04

17  1-1.20752660E-01-5.44061776E-02 6.98080232E-03-7.61571859E-03 0.00000000E+00

17  2 1.26467084E-03 3.77441857E-02 1.55298284E-02-1.50142357E-01-2.41013235E-01

17  3 0.00000000E+00-5.54107417E-02-8.00560552E-02 0.00000000E+00-1.00701926E-02

17  4 8.26145286E-03 1.80873979E-03 2.62978114E-02 0.00000000E+00 0.00000000E+00

17  5 3.20222811E-04 2.51049044E-02 1.56883153E-02-6.23036807E-02 3.45663940E-01

17  6 0.00000000E+00-2.90936594E-02 1.37554053E-01 0.00000000E+00-5.19097312E-03

17  7 5.48683955E-03-2.95866423E-04 1.75670993E-02 0.00000000E+00 0.00000000E+00

17  8 4.09849882E-03-5.60930914E-02-6.46289273E-02 1.41747221E-01 1.81777360E-01

17  9 0.00000000E+00 5.33383724E-02 7.56778348E-02 0.00000000E+00-1.78410039E-02

17 10 2.00677407E-02-2.22673673E-03 2.07655490E-02 0.00000000E+00 0.00000000E+00

17 11-2.67105399E-01-1.02731145E-01 1.23716767E-02-1.17190701E-02 0.00000000E+00

17 12 5.33997386E-03 1.12027945E-01 1.10516073E-01-7.52533218E-02-1.80453549E-01

17 13 0.00000000E+00-7.24696579E-03-2.85257871E-02 0.00000000E+00 2.99252696E-03

17 14 1.26302715E-02-1.56227984E-02 9.62268044E-04 0.00000000E+00 0.00000000E+00

17 15 2.84478666E-02 7.29646638E-03 1.67002909E-03 9.49898176E-03 0.00000000E+00

17 16-1.69389546E-01-7.57774877E-02-1.05093138E-02 6.22612080E-03 0.00000000E+00

17 17 3.44259485E-04-5.29035124E-03-1.56275481E-02 4.57332538E-02 1.13158249E-01

17 18-8.14910767E-03 3.23191855E-02 8.16375426E-02-4.37626150E-03 8.72017424E-04

17 19 3.72497502E-04-1.24451493E-03-3.14104612E-04 7.94688315E-04 1.79890171E-03

17 20 3.44259485E-04-5.29035124E-03-1.56275481E-02 4.57332538E-02 1.13158249E-01

17 21 8.14910767E-03 3.23191855E-02 8.16375426E-02 4.37626150E-03 8.72017424E-04

17 22 3.72497502E-04-1.24451493E-03-3.14104612E-04-7.94688315E-04-1.79890171E-03

18  1 2.88191151E-01 1.54676074E-01-1.43961864E-02-2.85832647E-03 0.00000000E+00

18  2-4.35306867E-03-1.21476054E-02 1.88943156E-02 2.89686659E-01-4.16824129E-02

18  3 0.00000000E+00 1.04721848E-01 2.21056468E-02 0.00000000E+00 3.37246528E-02

18  4-3.44147507E-02 6.90097894E-04 1.67745182E-02 0.00000000E+00 0.00000000E+00

18  5-2.52059495E-04-2.23502028E-02-1.95601315E-03-2.83383543E-01 3.42020078E-02

18  6 0.00000000E+00-1.05648271E-01 2.40941068E-02 0.00000000E+00 2.68769729E-02

18  7-2.54153947E-02-1.46157824E-03 2.11369665E-02 0.00000000E+00 0.00000000E+00

18  8 1.75554309E-03-4.99298957E-02-1.20951199E-01-1.69613413E-01 1.28863738E-01

18  9 0.00000000E+00-1.02323374E-01 6.94922170E-02 0.00000000E+00 1.61363025E-02

18 10-1.61183993E-02-1.79031975E-05-5.41575893E-03 0.00000000E+00 0.00000000E+00

18 11-2.06189947E-01-1.02754861E-01-1.05451751E-03-1.30370480E-02 0.00000000E+00

18 12-1.52962534E-03-1.13102300E-02 1.11771295E-02-2.40074794E-02 2.55443494E-02

18 13 0.00000000E+00-3.82630368E-03-8.69665703E-03 0.00000000E+00 1.05724490E-02

18 14-9.97384948E-03-5.98599556E-04 9.15911139E-03 0.00000000E+00 0.00000000E+00

18 15 2.61250898E-01 1.39249204E-01 1.11403787E-02 6.00153009E-03 0.00000000E+00

18 16 1.41382132E-01 6.88131850E-02 6.10879798E-03 1.34959434E-03 0.00000000E+00

18 17 1.63302759E-04 4.90545813E-05-1.15890701E-03 4.49442434E-02-2.48757464E-02

18 18-9.67282398E-04 3.28011216E-02-1.56869078E-02-9.64723396E-04-7.03452799E-05

18 19 3.62067654E-05 3.41385144E-05-2.24080148E-04 5.13254884E-04-2.11867360E-04

18 20 1.63302759E-04 4.90545813E-05-1.15890701E-03 4.49442434E-02-2.48757464E-02

18 21 9.67282398E-04 3.28011216E-02-1.56869078E-02 9.64723396E-04-7.03452799E-05

18 22 3.62067654E-05 3.41385144E-05-2.24080148E-04-5.13254884E-04 2.11867360E-04

19  1 1.46430508E-02 5.11254556E-03-7.91082541E-04-1.06359914E-03 1.52430005E-12

19  2-2.82472895E-04-5.89713357E-03-1.12548578E-02 1.79218173E-02-3.47405183E-02

19  3-6.48352513E-13 1.85457967E-02-5.90452979E-03-2.34961884E-11 1.48867017E-03

19  4-2.73606692E-03 1.24739675E-03 4.22105718E-03-3.02871003E-12 4.54325262E-12

19  5-6.06762904E-04 2.30694119E-03 9.07033279E-03 2.25813647E-03 1.31157922E-02

19  6-1.07792087E-11-1.44349921E-03 1.01126087E-02-7.52769273E-12 1.29995019E-03

19  7-9.00974018E-04-3.98976176E-04 1.85429573E-03 1.19279714E-12-6.94374771E-13

19  8-7.15338631E-04-6.61547812E-04 8.32038227E-04-6.34576260E-02 2.95266419E-02

19  9 1.09463491E-10-1.41216774E-02-4.23398752E-03 8.66364465E-11-9.31836041E-03

19 10 8.92782568E-03 3.90534730E-04-8.00032204E-03 2.72013759E-11 5.07778124E-11

19 11-8.65959246E-03-3.80568363E-03 3.42374123E-05-3.22528301E-04-6.01463714E-13

19 12-1.51510366E-03-1.90250280E-02-1.82572122E-03-2.07611766E-01 9.42631099E-02

19 13 9.88253315E-10-1.03001345E-01 4.68888555E-02 4.25831028E-11-1.96981066E-02

19 14 2.07710100E-02-1.07290338E-03-2.20848497E-02 6.90368776E-11 1.75443903E-10

19 15-2.14981653E-04 1.23817296E-03-3.50858266E-04 3.66678836E-04-2.68131417E-14

19 16 4.63408834E-02 8.79963680E-03 5.47910851E-03-1.60369824E-03 1.85166513E-12

19 17 2.05421736E-04 1.66922588E-04 1.34612319E-03-3.76351825E-01 1.78484433E-01

19 18-5.05632145E-03-2.41764989E-01 1.14143836E-01-2.87642808E-03 4.89169691E-03

19 19-4.75981916E-03-1.31877753E-04 4.73512029E-03-9.87736695E-03 4.97410675E-03

19 20 2.05421714E-04 1.66922762E-04 1.34612309E-03-3.76351826E-01 1.78484432E-01

19 21 5.05632007E-03-2.41764990E-01 1.14143835E-01 2.87642727E-03 4.89169689E-03

19 22-4.75981915E-03-1.31877738E-04 4.73512032E-03 9.87736698E-03-4.97410668E-03

20  1 3.83231405E-11 1.36279135E-11-2.03315641E-12-2.77744025E-12-6.04209876E-04

20  2-6.78276284E-13-1.61181178E-11-3.22430070E-11 4.72925440E-11-9.22770413E-11

20  3-1.62475027E-03 4.93392774E-11-1.62836230E-11 6.85428500E-03 3.96457659E-12

20  4-7.19645100E-12 3.23187441E-12 1.10783433E-11 1.23956488E-03-1.73826486E-03

20  5-1.58352094E-12 6.23358315E-12 2.43985989E-11 6.14267941E-12 3.47843067E-11

20  6 1.96957187E-03-3.68720703E-12 2.70128207E-11 1.75069229E-03 3.44102217E-12

20  7-2.37941671E-12-1.06160545E-12 4.91110969E-12-5.00246940E-04 1.69489615E-04

20  8-1.88988911E-12-1.79428683E-12 2.06171725E-12-1.66950755E-10 7.74814559E-11

20  9-4.24301704E-02-3.69994744E-11-1.21122178E-11-3.28048854E-02-2.39363629E-11

20 10 2.29120641E-11 1.02429880E-12-2.06235048E-11-1.04379748E-02-1.93001846E-02

20 11-2.28460769E-11-1.01504632E-11 9.53003318E-14-8.43552924E-13 1.41756921E-04

20 12-4.00652445E-12-4.98406343E-11-3.73929706E-12-5.36864898E-10 2.43789298E-10

20 13-3.76460778E-01-2.65319431E-10 1.20426954E-10-1.67276983E-02-5.17119445E-11

20 14 5.45363781E-11-2.82443352E-12-5.79850100E-11-3.00051030E-02-6.50673842E-02

20 15-6.05143562E-13 3.28549149E-12-9.31836147E-13 9.79541628E-13-7.85761496E-05

20 16 1.23124218E-10 2.50835705E-11 1.43707025E-11-4.19109297E-12-7.06600211E-04

20 17-4.26087269E-03 3.32548972E-02-1.92663528E-02-1.10634007E-01-2.76715693E-01

20 18 2.62848572E-01-6.53483159E-02-1.66577414E-01 1.53262651E-01-5.28107849E-03

20 19 2.27579700E-03 3.00528149E-03 4.60374194E-03-5.12169310E-03-1.20094351E-02

20 20 4.26087269E-03-3.32548972E-02 1.92663528E-02 1.10634005E-01 2.76715694E-01

20 21 2.62848572E-01 6.53483146E-02 1.66577414E-01 1.53262651E-01 5.28107852E-03

20 22-2.27579703E-03-3.00528149E-03-4.60374192E-03-5.12169304E-03-1.20094351E-02

21  1 3.29938950E-12 7.98625152E-13-2.31379262E-13-2.48448352E-13-5.15287961E-04

21  2-1.50701156E-13-3.52015128E-13 1.44188465E-12 3.48392074E-12-5.89300377E-12

21  3 3.22377993E-06 3.03413638E-12-1.45108662E-13 1.06202938E-02 2.37676943E-13

21  4-5.65490219E-13 3.27813276E-13 9.06594992E-13-1.64807095E-04 9.06477039E-04

21  5-1.43098548E-13 2.40142389E-13 1.10164729E-12 1.76590301E-13 2.30129695E-12

21  6 8.30944329E-04-4.54627523E-13 1.49858274E-12-3.76398823E-03 2.37461578E-13

21  7-1.72428427E-13-6.50331507E-14 3.34853562E-13-1.52920648E-04-3.05901925E-04

21  8-1.35868305E-13-5.92749517E-14 3.52818086E-13-1.30534840E-11 6.35949686E-12

21  9-1.42500167E-03-3.12361572E-12 5.17673829E-13-8.70339787E-03-2.74359676E-12

21 10 2.65875311E-12 8.48436496E-14-2.25159241E-12 1.21274079E-05-1.16043200E-04

21 11-1.69062762E-12-5.85778502E-13-4.14114182E-16-7.34627800E-14 2.34947889E-04

21 12-2.82497676E-13-4.21655802E-12-1.89402202E-12-5.60376745E-11 2.53951695E-11

21 13-8.31009956E-03-2.92752770E-11 1.38310068E-11 4.04538681E-03-4.21157602E-12

21 14 4.42980673E-12-2.18230714E-13-4.71141551E-12 4.41930086E-02-2.08793491E-02

21 15 1.22177724E-14 2.14765042E-13-5.96674457E-14 5.42403703E-14 2.41694882E-04

21 16 7.81148673E-12-9.47950601E-13 1.19028419E-12-3.69959424E-13-4.22004102E-04

21 17-3.10815921E-04-1.23622024E-03-4.23190934E-04-4.36457243E-01 1.86920747E-01

21 18 7.87053854E-03-3.06561099E-01 1.31712176E-01 5.00869672E-03 4.05302040E-03

21 19-4.17159488E-03 1.18574473E-04 4.48451797E-03-7.97610777E-03 3.18644247E-03

21 20 3.10815921E-04 1.23622024E-03 4.23190935E-04 4.36457243E-01-1.86920747E-01

21 21 7.87053854E-03 3.06561099E-01-1.31712176E-01 5.00869672E-03-4.05302040E-03

21 22 4.17159488E-03-1.18574473E-04-4.48451797E-03-7.97610777E-03 3.18644247E-03

22  1 5.07457448E-13-1.02344580E-12-2.07952372E-13-6.67505517E-14 1.30433149E-02

22  2-3.05592285E-13 3.11450345E-12 1.31263749E-11-1.26219738E-12 5.37026927E-12

22  3 3.84348264E-01-3.24664121E-12 3.82296270E-12 2.48931079E-01-2.80563856E-13

22  4 7.84876200E-14 2.02076236E-13-4.92146726E-15-1.44248617E-02-6.09710462E-03

22  5-4.27130829E-14-8.70317603E-13-2.88682077E-12-1.05058812E-12-1.76760764E-12

22  6 3.08402746E-01-4.89636728E-13-2.29993545E-12 2.07227841E-01-1.43347102E-13

22  7 7.26855123E-14 7.06615898E-14-2.17626040E-13 1.07321243E-02 2.29796610E-02

22  8 6.12202903E-14 2.82159380E-13 5.90655135E-13 2.03081808E-12 2.60243066E-14

22  9 2.52566940E-01-2.91299231E-13 4.85393302E-12 1.54172851E-01-2.51110119E-12

22 10 2.50827773E-12 2.82345650E-15-1.80275271E-12 1.23642211E-02-1.36014436E-02

22 11 5.85270640E-13 7.91406729E-13-2.64336333E-14-1.38627302E-14 1.04455471E-02

22 12 1.47588403E-13-4.20836372E-13-5.10119991E-12-3.86539958E-11 1.73874454E-11

22 13 2.87901850E-02-2.41843796E-11 1.27224128E-11 1.97873000E-02 8.80919659E-14

22 14-1.30819820E-13 4.27278545E-14 1.34310838E-13-1.11824928E-02-2.43651664E-02

22 15 1.98660116E-13-1.75310971E-13 5.37211331E-14-8.37262871E-14 1.07426436E-02

22 16-7.33091945E-12-9.65333196E-12 3.94362462E-14-8.48402023E-14 9.46426272E-03

22 17-9.41193426E-05 1.09611964E-02 1.97151737E-02-1.58715966E-02-4.08019248E-02

22 18 1.13027774E-02-1.33247652E-02-3.20992177E-02 5.58410790E-03 1.80255780E-04

22 19-6.62646516E-05-1.13991128E-04-9.14693673E-05-7.38468365E-05-2.19198263E-04

22 20 9.41193425E-05-1.09611964E-02-1.97151737E-02 1.58715966E-02 4.08019247E-02

22 21 1.13027774E-02 1.33247653E-02 3.20992177E-02 5.58410790E-03-1.80255779E-04

22 22 6.62646510E-05 1.13991128E-04 9.14693679E-05-7.38468356E-05-2.19198264E-04

23  1-6.40135887E-14 3.65980330E-13 6.18514503E-14 1.43177166E-14 3.07284106E-03

23  2 9.69097175E-14-1.04767854E-12-4.32256047E-12 5.30783224E-13-1.97449657E-12

23  3 8.88568454E-02 1.17690003E-12-1.27687387E-12 9.53083860E-02 1.00924336E-13

23  4-4.40934266E-14-5.68309096E-14 3.04524433E-14 2.32897065E-03-4.19961240E-02

23  5 9.66654088E-15 2.97269634E-13 9.95732750E-13 3.55242817E-13 6.61393376E-13

23  6 3.13814990E-01 1.48498947E-13 8.13036338E-13 2.78293996E-01 5.52524788E-14

23  7-2.96697044E-14-2.55827745E-14 8.30675591E-14 4.40289080E-03 7.15872997E-03

23  8-2.46921366E-14-9.57851093E-14-1.85353183E-13-1.09072870E-12 1.93465215E-13

23  9-4.66244644E-01-2.33605084E-15-1.59890992E-12-4.15980920E-01 7.48480604E-13

23 10-7.50237586E-13 1.75698201E-15 5.28383501E-13-7.40965204E-03-1.66953006E-02

23 11-2.48508754E-13-2.81993977E-13 8.78381190E-15 2.27855198E-15 1.15828321E-02

23 12-5.80953585E-14 6.03588446E-15 1.63745412E-12 1.10827823E-11-4.97930718E-12

23 13-1.17853061E-01 7.11796836E-12-3.79435650E-12-1.02357222E-01-1.63174439E-13

23 14 1.84330142E-13-2.11557028E-14-1.94442398E-13 2.72732228E-02 5.67749899E-02

23 15-6.57247179E-14 6.51685143E-14-1.97745398E-14 2.95876064E-14 1.02699998E-02

23 16 2.68796638E-12 3.18244896E-12 2.47070166E-14 1.64758308E-14-1.77060621E-02

23 17 6.93182761E-04-2.87168345E-02-7.64535795E-02 4.57052663E-02 9.53094384E-02

23 18-3.33176548E-02 3.87749327E-02 8.17461952E-02-2.26152013E-02-4.43240820E-04

23 19 8.18971100E-04-3.75730280E-04 8.93862431E-04 8.45487642E-05-1.36827353E-05

23 20-6.93182761E-04 2.87168345E-02 7.64535795E-02-4.57052663E-02-9.53094384E-02

23 21-3.33176548E-02-3.87749327E-02-8.17461952E-02-2.26152013E-02 4.43240820E-04

23 22-8.18971100E-04 3.75730280E-04-8.93862431E-04 8.45487640E-05-1.36827352E-05

24  1-5.38436405E-03 1.89362782E-02 3.42101501E-03 9.09345225E-04 6.24618730E-13

24  2 5.23245741E-03-5.53736359E-02-2.30207034E-01 2.60620488E-02-1.01208260E-01

24  3 1.60746680E-11 6.06135083E-02-6.76585192E-02 8.71520497E-12 5.21143961E-03

24  4-1.99868847E-03-3.21275114E-03 1.06941780E-03-7.97440101E-13 4.90549949E-13

24  5 5.99093268E-04 1.56272709E-02 5.21659025E-02 1.87413956E-02 3.37043066E-02

24  6 1.70514328E-11 8.15271231E-03 4.22335683E-02 9.05770295E-12 2.78839850E-03

24  7-1.47026194E-03-1.31813655E-03 4.20547464E-03 5.15559622E-13 1.16231315E-12

24  8-1.22820095E-03-5.04626865E-03-1.00467091E-02-5.00044195E-02 6.43258602E-03

24  9 1.78983585E-11 1.75819781E-03-8.51438588E-02 2.11383806E-11 4.13646245E-02

24 10-4.14067204E-02 4.20959481E-05 2.93906713E-02 1.23170365E-12 9.11021616E-13

24 11-1.21662163E-02-1.46085404E-02 4.66284376E-04 1.65156208E-04 4.54643342E-13

24 12-2.91209405E-03 2.86093419E-03 8.80181843E-02 6.21774856E-01-2.79484547E-01

24 13 1.21205112E-11 3.95291544E-01-2.09647674E-01 1.04536188E-11-6.11970849E-03

24 14 7.11065243E-03-9.90943941E-04-7.47748553E-03-2.47254790E-12-3.34250194E-12

24 15-3.49445941E-03 3.32818056E-03-1.01320736E-03 1.53712054E-03 5.15398687E-13

24 16 1.37906204E-01 1.69419537E-01 5.93650466E-04 1.09706305E-03 1.06872901E-12

24 17 2.93898442E-04-1.86032875E-03-8.96950364E-03-2.15571584E-01 1.01623010E-01

24 18-3.35825493E-03-1.87700509E-01 9.20826638E-02-4.26682719E-03-3.49987295E-03

24 19 3.64657425E-03-1.46701301E-04-3.89898240E-03 6.26141588E-03-2.62435011E-03

24 20 2.93898442E-04-1.86032876E-03-8.96950365E-03-2.15571584E-01 1.01623010E-01

24 21 3.35825493E-03-1.87700509E-01 9.20826638E-02 4.26682719E-03-3.49987294E-03

24 22 3.64657425E-03-1.46701301E-04-3.89898240E-03-6.26141588E-03 2.62435011E-03

25  1 6.23036358E-14-1.67130110E-13-3.17682312E-14-9.22796645E-15-1.18743691E-02

25  2-4.77379329E-14 4.97036771E-13 2.07852894E-12-2.20025530E-13 8.86446415E-13

25  3-5.68368491E-01-5.32967319E-13 6.08501177E-13-4.66907000E-01-4.59207791E-14

25  4 1.56232696E-14 3.02975095E-14-5.82753902E-15 1.33624015E-04 3.59852620E-04

25  5-5.99327357E-15-1.39681060E-13-4.65016326E-13-1.67982980E-13-2.93783612E-13

25  6 5.48607612E-01-7.53020611E-14-3.73934685E-13 3.51167814E-01-2.41077672E-14

25  7 1.25135727E-14 1.15941944E-14-3.64570489E-14-3.13982329E-03-1.41585109E-02

25  8 1.04876164E-14 4.51815892E-14 9.19264006E-14 3.95459478E-13-3.12539982E-14

25  9 1.74537194E-01-2.89215669E-14 7.68694599E-13 2.01171569E-01-3.83893485E-13

25 10 3.83916830E-13-2.33441456E-17-2.74033266E-13-1.57258355E-04 3.67672706E-02

25 11 1.02442324E-13 1.29060955E-13-4.19955463E-15-1.79478067E-15 1.03180941E-02

25 12 2.50325866E-14-4.34303813E-14-8.00340434E-13-5.83254225E-12 2.62256293E-12

25 13 1.38328656E-01-3.68140268E-12 1.94536977E-12 1.20493405E-01 3.74405422E-14

25 14-4.54469721E-14 8.00642986E-15 4.75690614E-14-1.83699065E-02-3.64721736E-02

25 15 3.15107856E-14-2.90623991E-14 8.87149911E-15-1.36108812E-14 1.04824933E-02

25 16-1.20878547E-12-1.52920689E-12-3.55228108E-16-1.14272632E-14 6.65974169E-03

25 17-1.08092871E-03 2.22164977E-02 7.78801984E-02-3.33432351E-02-6.17461419E-02

25 18 3.36515166E-02-2.40744017E-02-4.93976695E-02 1.82746714E-02 4.82931011E-04

25 19-9.88574535E-04 5.05643525E-04-1.16270474E-03 1.18761313E-04 4.84124073E-04

25 20 1.08092871E-03-2.22164977E-02-7.78801984E-02 3.33432351E-02 6.17461419E-02

25 21 3.36515166E-02 2.40744017E-02 4.93976695E-02 1.82746714E-02-4.82931011E-04

25 22 9.88574535E-04-5.05643525E-04 1.16270474E-03 1.18761314E-04 4.84124073E-04

26  1-4.75262415E-14 1.42474135E-13 2.64865515E-14 7.41230811E-15-2.38664163E-03

26  2 4.01070952E-14-4.20569465E-13-1.75422966E-12 1.91343656E-13-7.58246269E-13

26  3-1.08600207E-01 4.55097826E-13-5.14441576E-13-8.36203229E-02 3.91746220E-14

26  4-1.40806636E-14-2.50939583E-14 6.33245213E-15-7.59791303E-03 4.12512581E-03

26  5 4.84241824E-15 1.18410932E-13 3.94674584E-13 1.42228886E-13 2.51836686E-13

26  6 1.99391209E-02 6.29282937E-14 3.18322105E-13 7.51275163E-02 2.07411188E-14

26  7-1.08424095E-14-9.89870929E-15 3.13281681E-14-2.02224220E-04-7.65653880E-03

26  8-9.07359251E-15-3.82729230E-14-7.71348568E-14-3.54455532E-13 3.62870902E-14

26  9 6.67212833E-01 1.95895595E-14-6.48784390E-13-5.90892915E-02 3.20150068E-13

26 10-3.20301550E-13 1.51482330E-16 2.28076314E-13-2.29639223E-02-6.89734246E-02

26 11-8.91946827E-14-1.09973118E-13 3.54820653E-15 1.40260452E-15 1.50954868E-03

26 12-2.15925784E-14 3.01530381E-14 6.73389466E-13 4.84179008E-12-2.17677015E-12

26 13-6.53779793E-01 3.06552122E-12-1.62246370E-12-5.05290802E-01-3.81795068E-14

26 14 4.52842653E-14-7.10475852E-15-4.75148604E-14 3.19178331E-02 7.30628386E-02

26 15-2.66093360E-14 2.48928035E-14-7.58955275E-15 1.15861578E-14 4.82696761E-04

26 16 1.03361936E-12 1.29079001E-12 2.14860117E-15 9.08213541E-15-5.74132157E-03

26 17-4.21466882E-03-1.01633535E-01-3.02102573E-01 1.02432945E-01 2.38523041E-01

26 18-1.59942491E-01 5.20891788E-02 1.15645352E-01-5.39727172E-02-4.68977408E-03

26 19 7.89856408E-03-3.20879000E-03 8.04731115E-03-2.45366871E-03-6.20345588E-03

26 20 4.21466882E-03 1.01633535E-01 3.02102573E-01-1.02432945E-01-2.38523041E-01

26 21-1.59942491E-01-5.20891788E-02-1.15645352E-01-5.39727172E-02 4.68977408E-03

26 22-7.89856408E-03 3.20879000E-03-8.04731115E-03-2.45366871E-03-6.20345588E-03

 $end

 

 

Sample 2. Optimization of CI in ethene molecule at XMCQDPT2 level using penalty function based code and numerical gradients.

 

$CONTRL                                                                        

! set nzvar to non-zero to allow automatic DLCs

  nzvar=1                                

  SCFTYP=mcscf mplevl=2

! Some standard stuff for XMCQDPT2 numerical derivatives.

  RUNTYP=optimize inttyp=hondo icut=13 numder=1                                                                     

! We want to use spherical basis set in this particular run

  d5=.t.

 $END                                                                          

 $SYSTEM

  TIMLIM=6000 MWORDS=40

! To speed up things a bit

  nojac=1 kdiag=0

 $END                            

! This group is required for extreme precision needed for numerical derivatives

 $moorth nostf=1 nozero=1 tole=0 tolz=0 $end          

 

 $BASIS  GBASIS=n311 ngauss=6 NDFUNC=1 POLAR=DUNNING $END                      

 $guess guess=moread norb=9 $end                                               

! we need to increase precision here a bit

 $trans altpar=.t. mptran=2 dirtrf=.t. mode=100 cuttrf=1d-13 $end

! we need to increase precision here a bit

 $MCSCF CISTEP=ALDET soscf=.t. focas=.f. MAXIT=100 acurcy=1d-8 $END

 $DET

   group=c1 NCORE=7 NACT=2 NELS=2 NSTATE=4 WSTATE(1)=1,1

! we need to increase precision here a bit

   itermx=1000 cvgtol=1d-8

 $END

! This input is specific to penalty-based CI search algorithm

! (using penalty function due to work of Todd Martinez)

 $mcaver target=mixed penlty=2 xgrad=.f. $end                                                     

 $xmcqdpt                                                                       

   nstate=2 edshft=0.02 istsym=1 wstate(1)=1,1

! we need to increase precision a bit

   thrgen=1d-20 threri=1d-20 genzro=1d-20 thrcon=1d-8             

 $end                                                                           

! To compute two (ngrads=2) gradients for istate=1 and jstate=2

! Using special fitting for CIs (spline=4 confit=1 praxes=1) to improve

! precision of numerical gradients near CI point

 $numgrd                                                                        

  order=6 delta=0.01 istate=1 jstate=2 ngrads=2 spline=4 confit=1 praxes=1

 $end                                                                          

! Using qa optimizer

 $STATPT NSTEP=400 method=qa opttol=3d-5 fmaxt=1000 $END                       

! To allow DLCs

 $zmat dlc=1 auto=1 $end                                                       

 $DATA

 C2H4

C1

 CARBON      6.0  -0.7011340567  -0.0539373538  -0.0388703175

 CARBON      6.0   0.7442178559   0.0435817054   0.0094540087

 HYDROGEN    1.0  -1.1838841441   0.8659009018   0.2781695799

 HYDROGEN    1.0  -1.4854226635  -0.8335637860  -0.0770387381

 HYDROGEN    1.0   0.9848952054  -0.1578508595  -1.0643856054

 HYDROGEN    1.0   0.1756554144  -0.8166558976   0.6576508518

 $end

--- OPTIMIZED MCSCF MO-S --- GENERATED AT 16:47:16 30-OCT-2011   

 C2H4                                                                          

E(MCSCF)=      -77.7107004274, 14 ITERS, E(NUC)=   32.6937128999

 $VEC  

 1  1-7.14795677E-02-5.91203307E-02-1.85285756E-04 9.09327041E-05-1.27695598E-04

 1  2-8.58902665E-04-3.13010663E-04-3.47842651E-05 1.70829653E-04 1.72345239E-04

 1  3-1.25284817E-05 7.12738059E-06 1.22125691E-04-1.81489875E-04 9.24988078E-05

 1  4 8.89910670E-05-2.09705033E-06 4.14078680E-05 1.60096666E-06-5.58987249E-01

 1  5-4.62937380E-01 1.36100583E-03 1.63536729E-03 5.18428728E-04-3.10680461E-03

 1  6-6.02780586E-04-7.81871785E-04-1.20924519E-04-1.11501703E-04 2.72005938E-04

 1  7 2.07180112E-04-1.51456852E-04 3.46002515E-05-1.43268711E-05-2.02733804E-05

 1  8 2.16810574E-05 1.91798281E-04-7.64143630E-05-2.71968506E-05 8.33247144E-05

 1  9-6.60529605E-05-4.54412382E-05-7.44805470E-05-5.90966009E-05-2.65459230E-04

 1 10 8.74392705E-05-2.26824378E-04-2.71324604E-04-5.40500830E-06 1.13304517E-04

 2  1 5.58794991E-01 4.62835829E-01 6.74567532E-04-5.64088954E-04 1.31539123E-03

 2  2 4.06624131E-03-1.13002204E-04 1.05513989E-04-3.30915316E-04-4.58197576E-04

 2  3 1.51855311E-05-1.87215936E-04 1.21407172E-04 5.71252637E-05 1.45730933E-04

 2  4-2.02856197E-04-4.93355927E-05-1.31603127E-04-5.58335825E-05-7.16024497E-02

 2  5-5.93508612E-02 3.92976427E-05 2.16568890E-04 6.86941782E-05-1.75365672E-04

 2  6-1.86639119E-04-1.54761721E-04 4.76411883E-05 9.61655553E-05-9.19445831E-05

 2  7 4.02644123E-05-2.23633931E-05 2.99112841E-05-2.04944698E-05-9.41681426E-06

 2  8 6.00106629E-05-3.55014290E-06-2.70883716E-05 3.08836208E-04-4.36490040E-06

 2  9 2.31340427E-04 3.51150653E-04-5.89672836E-05 1.12037076E-04-5.29243854E-05

 2 10-1.46766205E-05 1.27104036E-05 2.65229142E-04-2.21832361E-04 9.86333663E-05

 3  1-8.31504372E-02-1.39721680E-01 4.81228440E-02-2.14996439E-02 1.75066190E-02

 3  2 3.65764221E-01 8.31264234E-02-4.17797828E-02 4.01229619E-02 2.25131433E-01

 3  3 1.45299599E-02-5.46265093E-03 1.24849146E-02 1.49502628E-02 5.33686731E-03

 3  4-2.02871301E-02-1.08382432E-02 3.58005524E-03-8.91100661E-03-6.20567421E-02

 3  5-1.03677332E-01-4.85641677E-02-2.17776105E-02-2.83607282E-03 2.70519291E-01

 3  6-8.49961925E-02-4.84267454E-02-6.54998680E-03 1.46782523E-01-7.91160723E-03

 3  7-1.32432128E-02-5.82372253E-03 1.49860930E-02-4.61810559E-03-1.03679874E-02

 3  8 1.77980719E-02-7.15206205E-03-1.91273736E-03 4.77872841E-02 4.21407919E-02

 3  9-1.34484059E-02 5.58348951E-02 4.75733769E-02-4.40160009E-03 3.42103930E-02

 3 10 2.85898341E-02-1.04490569E-02 1.32945085E-01 1.08518315E-01-1.16833159E-02

 4  1-5.20169029E-02-8.87473752E-02-8.29051857E-02-1.59966210E-02 3.83053268E-02

 4  2 2.33744234E-01-1.44313672E-01-3.61482034E-02 7.61632498E-02 2.79084547E-01

 4  3-2.04841928E-02-1.08093295E-02 3.51457244E-02-1.43509661E-02 1.85573276E-02

 4  4-4.20636151E-03 1.00082897E-02-1.12446151E-02-9.08791870E-03 7.29828702E-02

 4  5 1.23369790E-01-3.18350095E-02 1.09316529E-02 4.73350230E-02-3.28717733E-01

 4  6-5.32667700E-02 3.01766710E-02 8.79885166E-02-3.41170060E-01 1.48311128E-02

 4  7 2.70598592E-02 3.25554260E-02 3.98436885E-03 6.19640768E-03-1.01807765E-02

 4  8-3.41302448E-03 1.65734068E-02-3.59385821E-03 5.34847611E-02 5.80813265E-02

 4  9 6.63534518E-03 1.13257277E-01 1.40199738E-01 2.32201016E-02-1.13208189E-01

 4 10-1.31193441E-01-2.38573153E-02-7.41141931E-03-9.70666304E-03 6.91520293E-03

 5  1-2.22469839E-02-3.81468608E-02-4.26815515E-02 1.62063091E-01-7.12268177E-03

 5  2 1.00116858E-01-7.01016068E-02 2.88213346E-01-8.07412152E-03 1.47369231E-01

 5  3-3.51430313E-02 1.29868790E-01 8.29862547E-03-1.39568931E-02 1.62508856E-02

 5  4-2.29399248E-03 1.30835516E-02-1.02165080E-02 2.75633210E-02-1.38553607E-02

 5  5-2.33336398E-02 2.23617186E-02 6.09994802E-02-2.95863919E-02 6.10277638E-02

 5  6 3.64939768E-02 1.03214348E-01-5.29423955E-02 7.77749784E-02 8.77178788E-03

 5  7 3.17989999E-02-2.72330284E-02 3.26849409E-03-1.30856285E-02 9.81713445E-03

 5  8-2.23279622E-02-2.56382973E-03 3.00759800E-03 1.48837206E-01 1.99219184E-01

 5  9 5.61960804E-02-2.95966960E-02-5.44937484E-02-1.81999065E-02 4.38435037E-02

 5 10 5.50172841E-02 1.95751632E-02-1.21513256E-01-1.84109416E-01-4.99541141E-02

 6  1-3.53349381E-03-6.23207118E-03-1.30462105E-01-7.16322226E-02 8.85690754E-03

 6  2 9.41294304E-03-2.24803331E-01-1.29061934E-01 1.41973227E-02 4.38695801E-02

 6  3-9.82107012E-02-8.13883710E-02 1.00259994E-02-9.14366445E-03 5.19805182E-03

 6  4 3.94561263E-03 2.22882311E-02-2.50092396E-02-1.53964855E-02-1.32117079E-02

 6  5-2.27097063E-02 1.22408717E-01 3.84992580E-03-7.53131541E-02 5.99635482E-02

 6  6 2.12802905E-01 1.47532072E-02-1.35119138E-01 8.51736764E-02 9.47094178E-02

 6  7 5.96491589E-03-7.25385835E-02-1.95915087E-02-7.02094290E-03 2.66124516E-02

 6  8-1.47774471E-02-1.48141216E-02 8.15609584E-03-3.56757531E-02-4.93732949E-02

 6  9-9.63945608E-03 1.27958136E-01 1.91983582E-01 6.79742843E-02 1.26164532E-01

 6 10 1.76149352E-01 6.14639382E-02-5.15842454E-02-6.76830385E-02-1.27522719E-02

 7  1-3.36465217E-02-5.86573748E-02 3.30857132E-02-2.13513234E-02-1.61662942E-01

 7  2 1.78064641E-01 7.09201557E-02-4.59707539E-02-2.69297602E-01 2.11010829E-01

 7  3 3.04918255E-02-8.44988496E-02-3.03519729E-01 1.24490341E-02 4.38216085E-03

 7  4-1.68311949E-02-5.75232261E-03-1.89645701E-02 1.46762889E-03 3.25242794E-02

 7  5 5.57926275E-02-6.24008785E-02 8.91456655E-04-1.19842832E-01-1.58945633E-01

 7  6-9.73570164E-02 7.28335586E-03-2.02826753E-01-1.70026386E-01-9.85874326E-02

 7  7 1.61267447E-03-1.91127132E-01 6.26109876E-03-1.05349133E-02 4.27381449E-03

 7  8 3.56521550E-03 2.93436688E-02 1.16101365E-02-3.20782425E-02-2.39680205E-02

 7  9-2.38058977E-02-2.56653904E-02-2.79775085E-02-6.94034445E-03 1.61306081E-02

 7 10-1.19809284E-03 3.88112829E-03-1.57700511E-02-6.17772283E-02-9.07825204E-02

 8  1 2.13033949E-03 4.12921500E-03-2.88722296E-02 5.93356210E-02 1.21900883E-02

 8  2-9.96696627E-03-6.69339315E-02 1.14100288E-01 1.13188537E-02-1.22503066E-01

 8  3-1.33109762E-01 1.32887805E-01 3.23555318E-02-1.60929995E-02 1.21697051E-02

 8  4 3.92329447E-03-3.45019299E-02-2.83442404E-03 1.34033575E-02 3.20473065E-02

 8  5 5.45387801E-02 5.62719028E-02-2.15902334E-01-3.66698553E-02-1.63517759E-01

 8  6 8.80506465E-02-3.49848273E-01-5.22839180E-02-1.92366868E-01 6.87564464E-02

 8  7-4.43801753E-01-8.29050940E-02-1.08466915E-02 1.40654059E-02-3.21871439E-03

 8  8 9.88217400E-03-5.82676126E-03 8.10566210E-03 8.71952826E-02 1.71854819E-01

 8  9 1.50407148E-01-4.96094155E-02-1.22280775E-01-1.41939012E-01 3.75904698E-02

 8 10 5.11022917E-02 3.46805179E-02 6.20229047E-02 1.60638142E-01 1.77763971E-01

 9  1-8.93330037E-03-1.46694330E-02-6.75859323E-02 1.55631391E-03-1.66079821E-01

 9  2 4.56972414E-02-1.24918431E-01 8.84031312E-03-2.77335135E-01 2.28869804E-02

 9  3-1.54525264E-01 1.79741897E-02-4.36261107E-01-1.88777737E-03 5.44944041E-03

 9  4-3.56166304E-03 2.87785749E-03-2.11850693E-03 8.54845905E-03-2.98451460E-02

 9  5-5.09313044E-02 8.33214734E-02-2.96086168E-02 1.63299116E-01 1.40748286E-01

 9  6 1.19670267E-01-5.51940698E-02 2.56094423E-01 3.70903210E-01 1.31119178E-01

 9  7-8.64022400E-02 2.99040966E-01-1.68085722E-03 9.67950010E-03-7.99864288E-03

 9  8-6.32329100E-03 1.95423287E-02-1.39544444E-02-1.57436788E-02-1.16314196E-02

 9  9 2.04025498E-02-1.80527699E-02-1.84263845E-02-1.60320921E-02-7.57055466E-02

 9 10-1.68246279E-01-8.97894260E-02-2.50868305E-02-1.88534568E-02 7.83963584E-02

 $END  

 

Sample 3. Optimization of CI in allene molecule using Lagrange function based code and semi-numerical gradients for SA-MCSCF.

 

 $CONTRL SCFTYP=mcscf RUNTYP=optimize inttyp=hondo icut=11 d5=.t. nzvar=1 $END

 $SYSTEM TIMLIM=6000 MWORDS=10 nojac=1 kdiag=0 $END

 $moorth nostf=1 nozero=1 tole=0 tolz=0 syms=1 symden=1 symvec=1 $end

 $GUESS  GUESS=MOREAD  NORB=13  $END                                            

 $BASIS  GBASIS=DH NDFUNC=1 POLAR=DUNNING $END                                 

 $trans  cuttrf=1d-12 $end

 $DRT    GROUP=Cs FORS=.TRUE. NMCC=9 NDOC=2 NVAL=2 istsym=1 $END               

 $GUGEM  cutoff=1d-12 $END

 $GUGDIA NSTATE=4  ITERMX=200 cvgtol=1d-7 memmax=999999 $END

 $GUGDM2 cutoff=1d-12 WSTATE(1)=1,1 $END

 $mcscf  CISTEP=GUGA maxit=100 ISTATE=1 acurcy=1d-7 ENGTOL=1.0d-12 ntrack=4 $end                                                    

 $mcaver jstate=2 conic=2 $end

 $conic  shift0=1d-2 $end

 $track  tol=1.4 freeze=.t. sticky=.f. update=.t. reset=.f. delciv=.f. $end

 $statpt nstep=1000 opttol=3.0d-5 method=conic $end

 $DATA

 C3H4

Cs

 

 CARBON      6.0   0.2215746811  -0.0198271535   0.0000000000

 CARBON      6.0  -0.3058327517   0.1108888601   1.2870548725

 HYDROGEN    1.0   0.3005544890   0.4461493143   2.1106814061

 HYDROGEN    1.0  -1.2147135738  -0.4223697576   1.7357021469

 $END

--- NATURAL ORBITALS OF MCSCF --- GENERATED AT 23:19:08 LT  30-SEP-2007

 C3H4                                                                           

E(MCSCF)=     -115.7329336740,  0 ITERS, E(NUC)=   57.4137045913

 $VEC  

 1  1 9.98268169E-01 4.07567343E-03 1.17660574E-03-4.95636722E-04 1.72436525E-04

 1  2 0.00000000E+00 4.22767469E-04-9.80774906E-05 0.00000000E+00-4.43764843E-04

 1  3-7.98998843E-04 2.11890876E-05-1.56019463E-04 0.00000000E+00 0.00000000E+00

 1  4 2.46828339E-02 3.93083778E-04-3.17429087E-04 1.00543138E-04-3.98900944E-05

 1  5 2.89858946E-04-1.54283488E-04 7.39523129E-05-1.83425196E-04-5.65602149E-05

 1  6-4.49010611E-05 6.93086458E-05-2.94299392E-05 1.01380170E-04-7.54371357E-06

 1  7 2.46828339E-02 3.93083778E-04-3.17429087E-04 1.00543138E-04-3.98900944E-05

 1  8-2.89858946E-04-1.54283488E-04 7.39523129E-05 1.83425196E-04-5.65602149E-05

 1  9-4.49010611E-05 6.93086458E-05-2.94299392E-05-1.01380170E-04 7.54371357E-06

 1 10 2.32195384E-05-7.38766046E-05 2.32195384E-05-7.38766046E-05 7.97318262E-05

 1 11 7.52447517E-05 7.97318262E-05 7.52447517E-05

 2  1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 2  2-6.30464401E-04 0.00000000E+00 0.00000000E+00 1.43667068E-04 0.00000000E+00

 2  3 0.00000000E+00 0.00000000E+00 0.00000000E+00 2.09601457E-04-4.93843686E-05

 2  4 7.06187046E-01 3.68560070E-03 5.16459749E-04 7.68557246E-05-1.48919694E-04

 2  5 3.71903349E-05-2.72209704E-05 6.96355736E-05-8.86768427E-05-4.86612979E-04

 2  6-6.20561198E-04-4.12156413E-04 1.25188701E-04 1.70715773E-05-3.29447222E-05

 2  7-7.06187046E-01-3.68560070E-03-5.16459749E-04-7.68557246E-05 1.48919694E-04

 2  8 3.71903349E-05 2.72209704E-05-6.96355736E-05-8.86768427E-05 4.86612979E-04

 2  9 6.20561198E-04 4.12156413E-04-1.25188701E-04 1.70715773E-05-3.29447222E-05

 2 10 3.57502302E-04-1.95872709E-04-3.57502302E-04 1.95872709E-04 3.48221963E-04

 2 11-1.53771701E-04-3.48221963E-04 1.53771701E-04

 3  1 3.54142752E-02-3.41095443E-04 3.10375922E-04 9.12441314E-05 1.29218582E-04

 3  2 0.00000000E+00-1.72752611E-04 1.18659132E-04 0.00000000E+00 3.22872210E-05

 3  3 1.52151884E-05-2.80028155E-04 3.31891421E-05 0.00000000E+00 0.00000000E+00

 3  4-7.05779619E-01-3.65047054E-03-5.53888899E-04-1.43655856E-04 3.23329203E-04

 3  5 2.61794600E-05 1.18585906E-05 1.28484245E-05 1.89638569E-06 5.03190921E-04

 3  6 6.05674046E-04 4.05536733E-04-1.38318050E-04 1.49763777E-05 3.86495877E-05

 3  7-7.05779619E-01-3.65047054E-03-5.53888899E-04-1.43655856E-04 3.23329203E-04

 3  8-2.61794600E-05 1.18585906E-05 1.28484245E-05-1.89638569E-06 5.03190921E-04

 3  9 6.05674046E-04 4.05536733E-04-1.38318050E-04-1.49763777E-05-3.86495877E-05

 3 10-3.46540620E-04 1.49563462E-04-3.46540620E-04 1.49563462E-04-3.38198094E-04

 3 11 1.37705770E-04-3.38198094E-04 1.37705770E-04

 4  1 1.58355072E-01-3.19421897E-01-1.73518613E-01 7.32816146E-02-1.37907571E-02

 4  2 0.00000000E+00-1.87601373E-02 4.37495472E-04 0.00000000E+00-1.10721521E-03

 4  3 5.91128890E-03-2.62143030E-02 1.91099214E-03 0.00000000E+00 0.00000000E+00

 4  4 1.18961075E-01-2.34921207E-01-1.60013320E-01-2.81929051E-02 2.49254256E-02

 4  5-8.57512161E-02 6.44394112E-03 5.83715136E-03 5.08121393E-03-7.68310794E-03

 4  6 1.24870236E-03-2.00163280E-02-6.94586577E-03-5.76025418E-03 3.84713386E-03

 4  7 1.18961075E-01-2.34921207E-01-1.60013320E-01-2.81929051E-02 2.49254256E-02

 4  8 8.57512161E-02 6.44394112E-03 5.83715136E-03-5.08121393E-03-7.68310794E-03

 4  9 1.24870236E-03-2.00163280E-02-6.94586577E-03 5.76025418E-03-3.84713386E-03

 4 10-5.75988297E-02-6.43523791E-03-5.75988297E-02-6.43523791E-03-6.78605947E-02

 4 11-1.22974589E-02-6.78605947E-02-1.22974589E-02

 5  1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 5  2 2.60593206E-01 0.00000000E+00 0.00000000E+00-5.19060378E-02 0.00000000E+00

 5  3 0.00000000E+00 0.00000000E+00 0.00000000E+00-1.78433797E-02 3.30323150E-03

 5  4 1.46338150E-01-2.85484581E-01-2.57983950E-01-1.76490541E-02 1.96419635E-02

 5  5 3.54102766E-03-1.52112970E-02 6.74793653E-04-1.99168027E-03-1.30868202E-02

 5  6 4.95396654E-04-1.56257696E-02-1.13314976E-02-1.38849011E-03 4.76771030E-03

 5  7-1.46338150E-01 2.85484581E-01 2.57983950E-01 1.76490541E-02-1.96419635E-02

 5  8 3.54102766E-03 1.52112970E-02-6.74793653E-04-1.99168027E-03 1.30868202E-02

 5  9-4.95396654E-04 1.56257696E-02 1.13314976E-02-1.38849011E-03 4.76771030E-03

 5 10-1.05715734E-01-3.31921242E-02 1.05715734E-01 3.31921242E-02-1.00586590E-01

 5 11-4.34985146E-02 1.00586590E-01 4.34985146E-02

 6  1 1.25394722E-01-2.56643050E-01-2.27243474E-01 4.21644085E-02-3.18447405E-02

 6  2 0.00000000E+00-1.14347828E-02-1.49785755E-03 0.00000000E+00-1.99094956E-03

 6  3 1.71414517E-03-1.24413859E-02 2.65221846E-03 0.00000000E+00 0.00000000E+00

 6  4-6.47188060E-02 1.34950498E-01 1.39901479E-01-8.12270938E-02-1.51459932E-02

 6  5-2.94834847E-01-2.60146182E-03-8.12847013E-03-3.12733944E-02 1.69380752E-02

 6  6 4.83889211E-03-3.16175841E-03 1.49300791E-02-1.32277756E-02-1.69362038E-03

 6  7-6.47188060E-02 1.34950498E-01 1.39901479E-01-8.12270938E-02-1.51459932E-02

 6  8 2.94834847E-01-2.60146182E-03-8.12847013E-03 3.12733944E-02 1.69380752E-02

 6  9 4.83889211E-03-3.16175841E-03 1.49300791E-02 1.32277756E-02 1.69362038E-03

 6 10 1.35091339E-01 7.99885559E-02 1.35091339E-01 7.99885559E-02 1.21823545E-01

 6 11 7.23332112E-02 1.21823545E-01 7.23332112E-02

 7  1-2.00463614E-02 4.53714651E-02 5.06303336E-02 1.40284512E-01 7.21324055E-02

 7  2 0.00000000E+00 1.88962237E-02 8.56565396E-03 0.00000000E+00-1.13044727E-02

 7  3 3.35191116E-03 1.10346593E-02-3.57454646E-03 0.00000000E+00 0.00000000E+00

 7  4-2.39517701E-03 5.67232660E-03 1.19952465E-02 2.99869651E-01 1.78964920E-01

 7  5-1.11820265E-01 5.24356314E-02 3.60448322E-02-2.39206434E-02-5.82739905E-03

 7  6-5.07109146E-03 9.12652924E-03-6.75899385E-03-1.31163091E-02-5.25604823E-03

 7  7-2.39517701E-03 5.67232660E-03 1.19952465E-02 2.99869651E-01 1.78964920E-01

 7  8 1.11820265E-01 5.24356314E-02 3.60448322E-02 2.39206434E-02-5.82739905E-03

 7  9-5.07109146E-03 9.12652924E-03-6.75899385E-03 1.31163091E-02 5.25604823E-03

 7 10 1.52208147E-01 1.05346442E-01 1.52208147E-01 1.05346442E-01-1.53405153E-01

 7 11-1.18723611E-01-1.53405153E-01-1.18723611E-01

 8  1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 8  2 3.80457292E-01 0.00000000E+00 0.00000000E+00 1.23925894E-02 0.00000000E+00

 8  3 0.00000000E+00 0.00000000E+00 0.00000000E+00-1.41341016E-02 1.03010254E-02

 8  4 1.47197164E-03 1.25424968E-03 4.45846837E-03-2.81171090E-01-9.62968297E-02

 8  5-2.36454902E-01-5.39628280E-02-2.46633455E-02-3.46723278E-02 1.32108042E-02

 8  6 6.24353320E-03-1.06847866E-02 1.51600709E-02-4.71314144E-03 6.32793684E-03

 8  7-1.47197164E-03-1.25424968E-03-4.45846837E-03 2.81171090E-01 9.62968297E-02

 8  8-2.36454902E-01 5.39628280E-02 2.46633455E-02-3.46723278E-02-1.32108042E-02

 8  9-6.24353320E-03 1.06847866E-02-1.51600709E-02-4.71314144E-03 6.32793684E-03

 8 10-8.08073434E-04 2.18739521E-03 8.08073434E-04-2.18739521E-03 1.74133111E-01

 8 11 1.29767619E-01-1.74133111E-01-1.29767619E-01

 9  1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

 9  2-2.44894961E-01 0.00000000E+00 0.00000000E+00-1.52654620E-02 0.00000000E+00

 9  3 0.00000000E+00 0.00000000E+00 0.00000000E+00 2.55710536E-02 3.99535243E-03

 9  4-7.72306266E-03 1.20926255E-02 1.34428458E-02-1.87033067E-01-1.64449013E-01

 9  5 3.00806336E-01-3.38396979E-02-3.01619603E-02 5.88372178E-02 1.42841609E-03

 9  6 2.21996466E-03-9.67189089E-03 2.66617472E-05 2.38737614E-02 5.95485448E-03

 9  7 7.72306266E-03-1.20926255E-02-1.34428458E-02 1.87033067E-01 1.64449013E-01

 9  8 3.00806336E-01 3.38396979E-02 3.01619603E-02 5.88372178E-02-1.42841609E-03

 9  9-2.21996466E-03 9.67189089E-03-2.66617472E-05 2.38737614E-02 5.95485448E-03

 9 10-1.96656976E-01-1.56766868E-01 1.96656976E-01 1.56766868E-01 9.31817735E-02

 9 11 8.09154062E-02-9.31817735E-02-8.09154062E-02

10  1 4.92986215E-02-1.07684263E-01-1.86033270E-01-3.74529951E-01 3.33683214E-01

10  2 0.00000000E+00-1.53403911E-01 1.01287781E-01 0.00000000E+00 8.37715146E-03

10  3-1.62964976E-03-1.21685988E-02-1.38942701E-02 0.00000000E+00 0.00000000E+00

10  4-2.64292360E-02 6.34757011E-02 5.28952305E-02-8.61913573E-02 2.61313258E-01

10  5 1.05323220E-01-4.65480708E-02 8.79309032E-02 2.15107572E-02-6.62313258E-03

10  6-2.38125049E-03 1.68920593E-02 4.07721177E-03-8.84208039E-03 8.17481923E-03

10  7-2.64292360E-02 6.34757011E-02 5.28952305E-02-8.61913573E-02 2.61313258E-01

10  8-1.05323220E-01-4.65480708E-02 8.79309032E-02-2.15107572E-02-6.62313258E-03

10  9-2.38125049E-03 1.68920593E-02 4.07721177E-03 8.84208039E-03-8.17481923E-03

10 10 1.64021720E-02 2.78490305E-02 1.64021720E-02 2.78490305E-02-2.94343672E-02

10 11-7.12747553E-02-2.94343672E-02-7.12747553E-02

11  1 3.41423318E-02-6.86209113E-02-1.93472922E-01-5.20237716E-01-2.32966727E-01

11  2 0.00000000E+00-3.53807022E-01-5.65609435E-02 0.00000000E+00 1.09189845E-02

11  3-1.61160561E-02 4.70803872E-03 1.48654267E-02 0.00000000E+00 0.00000000E+00

11  4-1.34987060E-02 4.43066302E-02-8.15030212E-03 1.87337119E-01-1.81305704E-01

11  5 4.94611542E-02 3.38687378E-02-1.56979959E-01-1.49306654E-02-1.83899197E-02

11  6 1.27753162E-03 2.57693522E-02-5.67627815E-03-2.36741562E-02-1.88099963E-02

11  7-1.34987060E-02 4.43066302E-02-8.15030212E-03 1.87337119E-01-1.81305704E-01

11  8-4.94611542E-02 3.38687378E-02-1.56979959E-01 1.49306654E-02-1.83899197E-02

11  9 1.27753162E-03 2.57693522E-02-5.67627815E-03 2.36741562E-02 1.88099963E-02

11 10 7.88968944E-02 1.61310145E-01 7.88968944E-02 1.61310145E-01-8.69785254E-02

11 11-2.48820055E-01-8.69785254E-02-2.48820055E-01

12  1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00

12  2-3.76055105E-02 0.00000000E+00 0.00000000E+00-4.56325900E-02 0.00000000E+00

12  3 0.00000000E+00 0.00000000E+00 0.00000000E+00-1.70950040E-02 4.17396082E-02

12  4 2.34059726E-02-5.40604169E-02-1.26208302E-01 2.11404221E-01-4.33807964E-01

12  5-1.17619545E-01 1.02110163E-01-3.01522110E-01-6.95739297E-02-1.42586903E-03

12  6 1.24698259E-03-9.84279760E-04-4.71089020E-03 2.87311627E-04-4.37566456E-03

12  7-2.34059726E-02 5.40604169E-02 1.26208302E-01-2.11404221E-01 4.33807964E-01

12  8-1.17619545E-01-1.02110163E-01 3.01522110E-01-6.95739297E-02 1.42586903E-03

12  9-1.24698259E-03 9.84279760E-04 4.71089020E-03 2.87311627E-04-4.37566456E-03

12 10 2.27559714E-02 4.45340937E-02-2.27559714E-02-4.45340937E-02 1.06368058E-02

12 11-1.80991431E-02-1.06368058E-02 1.80991431E-02

13  1-9.68257699E-03 3.81602041E-02 1.41569126E-02 1.14155067E-01-6.86101209E-01

13  2 0.00000000E+00 4.16369113E-02-4.24399085E-01 0.00000000E+00 2.66463493E-02

13  3-1.48491722E-02-1.15059337E-02-2.58965157E-03 0.00000000E+00 0.00000000E+00

13  4-7.98458342E-03 2.46834268E-02 1.24291037E-01-1.87380537E-01 3.75682068E-01

13  5 4.25704843E-02-1.42511912E-01 2.16301194E-01 4.60749653E-02-4.64859882E-03

13  6 1.08429285E-02-9.93447317E-03-2.21012137E-02-8.67072903E-03-2.45262435E-02

13  7-7.98458342E-03 2.46834268E-02 1.24291037E-01-1.87380537E-01 3.75682068E-01

13  8-4.25704843E-02-1.42511912E-01 2.16301194E-01-4.60749653E-02-4.64859882E-03

13  9 1.08429285E-02-9.93447317E-03-2.21012137E-02 8.67072903E-03 2.45262435E-02

13 10 1.66203477E-02 5.98139921E-02 1.66203477E-02 5.98139921E-02-5.50608325E-02

13 11-1.54805163E-01-5.50608325E-02-1.54805163E-01

 $END  

 

Sample 4. Optimization of CI in ethene molecule at XMCQDPT2 level using Lagrange function based code and numerical gradients.

 

$CONTRL                                                                        

! set nzvar to non-zero to allow automatic DLCs

  nzvar=1                                

  SCFTYP=mcscf mplevl=2

! Some standard stuff for XMCQDPT2 numerical derivatives.

  RUNTYP=optimize inttyp=hondo icut=13 numder=1                                                                      

! We want to use spherical basis set in this particular run

  d5=.t.

 $END                                                                          

 $SYSTEM

  TIMLIM=6000 MWORDS=40

! To speed up things a bit

  nojac=1 kdiag=0

 $END                            

! This group is required for extreme precision needed for numerical derivatives

 $moorth nostf=1 nozero=1 tole=0 tolz=0 $end          

 $BASIS  GBASIS=n311 ngauss=6 NDFUNC=1 POLAR=DUNNING $END                       

 $guess guess=moread norb=9 $end                                               

! we need to increase precision here a bit

 $trans altpar=.t. mptran=2 dirtrf=.t. mode=100 cuttrf=1d-13 $end

! we need to increase precision here a bit

 $MCSCF CISTEP=ALDET soscf=.t. focas=.f. MAXIT=100 acurcy=1d-8 $END

 $DET

   group=c1 NCORE=7 NACT=2 NELS=2 NSTATE=4 WSTATE(1)=1,1

! we need to increase precision here a bit

   itermx=1000 cvgtol=1d-8

 $END

 $mcaver target=pure xgrad=.t. $end                                                      

 

 $xmcqdpt                                                                      

   nstate=2 edshft=0.02 istsym=1 wstate(1)=1,1

! we use this option to reduce printout

  lpout=10

! we need to increase precision a bit

   thrgen=1d-20 threri=1d-20 genzro=1d-20 thrcon=1d-8             

 $end                                                                          

! To compute two (ngrads=2) gradients of istate=1 and jstate=2

! Using special fitting for CIs (spline=4 confit=1 praxes=1) to improve

! precision of numerical gradients near CI point

 $numgrd                                                                       

  order=4 delta=0.01 istate=1 jstate=2 ngrads=2 spline=4 confit=1

 $end                                                                           

! Using conic optimizer

 $STATPT NSTEP=400 method=conic opttol=3d-5 $END                       

 $conic shift0=0.003 shift=0.001 diabas=.t. $end

! To allow DLCs

 $zmat dlc=1 auto=1 $end                                                        

 $DATA

 C2H4

C1

 CARBON      6.0  -0.7011340567  -0.0539373538  -0.0388703175

 CARBON      6.0   0.7442178559   0.0435817054   0.0094540087

 HYDROGEN    1.0  -1.1838841441   0.8659009018   0.2781695799

 HYDROGEN    1.0  -1.4854226635  -0.8335637860  -0.0770387381

 HYDROGEN    1.0   0.9848952054  -0.1578508595  -1.0643856054

 HYDROGEN    1.0   0.1756554144  -0.8166558976   0.6576508518

 $end

--- OPTIMIZED MCSCF MO-S --- GENERATED AT 16:47:16 30-OCT-2011   

 C2H4                                                                          

E(MCSCF)=      -77.7107004274, 14 ITERS, E(NUC)=   32.6937128999

 $VEC  

 1  1-7.14795677E-02-5.91203307E-02-1.85285756E-04 9.09327041E-05-1.27695598E-04

 1  2-8.58902665E-04-3.13010663E-04-3.47842651E-05 1.70829653E-04 1.72345239E-04

 1  3-1.25284817E-05 7.12738059E-06 1.22125691E-04-1.81489875E-04 9.24988078E-05

 1  4 8.89910670E-05-2.09705033E-06 4.14078680E-05 1.60096666E-06-5.58987249E-01

 1  5-4.62937380E-01 1.36100583E-03 1.63536729E-03 5.18428728E-04-3.10680461E-03

 1  6-6.02780586E-04-7.81871785E-04-1.20924519E-04-1.11501703E-04 2.72005938E-04

 1  7 2.07180112E-04-1.51456852E-04 3.46002515E-05-1.43268711E-05-2.02733804E-05

 1  8 2.16810574E-05 1.91798281E-04-7.64143630E-05-2.71968506E-05 8.33247144E-05

 1  9-6.60529605E-05-4.54412382E-05-7.44805470E-05-5.90966009E-05-2.65459230E-04

 1 10 8.74392705E-05-2.26824378E-04-2.71324604E-04-5.40500830E-06 1.13304517E-04

 2  1 5.58794991E-01 4.62835829E-01 6.74567532E-04-5.64088954E-04 1.31539123E-03

 2  2 4.06624131E-03-1.13002204E-04 1.05513989E-04-3.30915316E-04-4.58197576E-04

 2  3 1.51855311E-05-1.87215936E-04 1.21407172E-04 5.71252637E-05 1.45730933E-04

 2  4-2.02856197E-04-4.93355927E-05-1.31603127E-04-5.58335825E-05-7.16024497E-02

 2  5-5.93508612E-02 3.92976427E-05 2.16568890E-04 6.86941782E-05-1.75365672E-04

 2  6-1.86639119E-04-1.54761721E-04 4.76411883E-05 9.61655553E-05-9.19445831E-05

 2  7 4.02644123E-05-2.23633931E-05 2.99112841E-05-2.04944698E-05-9.41681426E-06

 2  8 6.00106629E-05-3.55014290E-06-2.70883716E-05 3.08836208E-04-4.36490040E-06

 2  9 2.31340427E-04 3.51150653E-04-5.89672836E-05 1.12037076E-04-5.29243854E-05

 2 10-1.46766205E-05 1.27104036E-05 2.65229142E-04-2.21832361E-04 9.86333663E-05

 3  1-8.31504372E-02-1.39721680E-01 4.81228440E-02-2.14996439E-02 1.75066190E-02

 3  2 3.65764221E-01 8.31264234E-02-4.17797828E-02 4.01229619E-02 2.25131433E-01

 3  3 1.45299599E-02-5.46265093E-03 1.24849146E-02 1.49502628E-02 5.33686731E-03

 3  4-2.02871301E-02-1.08382432E-02 3.58005524E-03-8.91100661E-03-6.20567421E-02

 3  5-1.03677332E-01-4.85641677E-02-2.17776105E-02-2.83607282E-03 2.70519291E-01

 3  6-8.49961925E-02-4.84267454E-02-6.54998680E-03 1.46782523E-01-7.91160723E-03

 3  7-1.32432128E-02-5.82372253E-03 1.49860930E-02-4.61810559E-03-1.03679874E-02

 3  8 1.77980719E-02-7.15206205E-03-1.91273736E-03 4.77872841E-02 4.21407919E-02

 3  9-1.34484059E-02 5.58348951E-02 4.75733769E-02-4.40160009E-03 3.42103930E-02

 3 10 2.85898341E-02-1.04490569E-02 1.32945085E-01 1.08518315E-01-1.16833159E-02

 4  1-5.20169029E-02-8.87473752E-02-8.29051857E-02-1.59966210E-02 3.83053268E-02

 4  2 2.33744234E-01-1.44313672E-01-3.61482034E-02 7.61632498E-02 2.79084547E-01

 4  3-2.04841928E-02-1.08093295E-02 3.51457244E-02-1.43509661E-02 1.85573276E-02

 4  4-4.20636151E-03 1.00082897E-02-1.12446151E-02-9.08791870E-03 7.29828702E-02

 4  5 1.23369790E-01-3.18350095E-02 1.09316529E-02 4.73350230E-02-3.28717733E-01

 4  6-5.32667700E-02 3.01766710E-02 8.79885166E-02-3.41170060E-01 1.48311128E-02

 4  7 2.70598592E-02 3.25554260E-02 3.98436885E-03 6.19640768E-03-1.01807765E-02

 4  8-3.41302448E-03 1.65734068E-02-3.59385821E-03 5.34847611E-02 5.80813265E-02

 4  9 6.63534518E-03 1.13257277E-01 1.40199738E-01 2.32201016E-02-1.13208189E-01

 4 10-1.31193441E-01-2.38573153E-02-7.41141931E-03-9.70666304E-03 6.91520293E-03

 5  1-2.22469839E-02-3.81468608E-02-4.26815515E-02 1.62063091E-01-7.12268177E-03

 5  2 1.00116858E-01-7.01016068E-02 2.88213346E-01-8.07412152E-03 1.47369231E-01

 5  3-3.51430313E-02 1.29868790E-01 8.29862547E-03-1.39568931E-02 1.62508856E-02

 5  4-2.29399248E-03 1.30835516E-02-1.02165080E-02 2.75633210E-02-1.38553607E-02

 5  5-2.33336398E-02 2.23617186E-02 6.09994802E-02-2.95863919E-02 6.10277638E-02

 5  6 3.64939768E-02 1.03214348E-01-5.29423955E-02 7.77749784E-02 8.77178788E-03

 5  7 3.17989999E-02-2.72330284E-02 3.26849409E-03-1.30856285E-02 9.81713445E-03

 5  8-2.23279622E-02-2.56382973E-03 3.00759800E-03 1.48837206E-01 1.99219184E-01

 5  9 5.61960804E-02-2.95966960E-02-5.44937484E-02-1.81999065E-02 4.38435037E-02

 5 10 5.50172841E-02 1.95751632E-02-1.21513256E-01-1.84109416E-01-4.99541141E-02

 6  1-3.53349381E-03-6.23207118E-03-1.30462105E-01-7.16322226E-02 8.85690754E-03

 6  2 9.41294304E-03-2.24803331E-01-1.29061934E-01 1.41973227E-02 4.38695801E-02

 6  3-9.82107012E-02-8.13883710E-02 1.00259994E-02-9.14366445E-03 5.19805182E-03

 6  4 3.94561263E-03 2.22882311E-02-2.50092396E-02-1.53964855E-02-1.32117079E-02

 6  5-2.27097063E-02 1.22408717E-01 3.84992580E-03-7.53131541E-02 5.99635482E-02

 6  6 2.12802905E-01 1.47532072E-02-1.35119138E-01 8.51736764E-02 9.47094178E-02

 6  7 5.96491589E-03-7.25385835E-02-1.95915087E-02-7.02094290E-03 2.66124516E-02

 6  8-1.47774471E-02-1.48141216E-02 8.15609584E-03-3.56757531E-02-4.93732949E-02

 6  9-9.63945608E-03 1.27958136E-01 1.91983582E-01 6.79742843E-02 1.26164532E-01

 6 10 1.76149352E-01 6.14639382E-02-5.15842454E-02-6.76830385E-02-1.27522719E-02

 7  1-3.36465217E-02-5.86573748E-02 3.30857132E-02-2.13513234E-02-1.61662942E-01

 7  2 1.78064641E-01 7.09201557E-02-4.59707539E-02-2.69297602E-01 2.11010829E-01

 7  3 3.04918255E-02-8.44988496E-02-3.03519729E-01 1.24490341E-02 4.38216085E-03

 7  4-1.68311949E-02-5.75232261E-03-1.89645701E-02 1.46762889E-03 3.25242794E-02

 7  5 5.57926275E-02-6.24008785E-02 8.91456655E-04-1.19842832E-01-1.58945633E-01

 7  6-9.73570164E-02 7.28335586E-03-2.02826753E-01-1.70026386E-01-9.85874326E-02

 7  7 1.61267447E-03-1.91127132E-01 6.26109876E-03-1.05349133E-02 4.27381449E-03

 7  8 3.56521550E-03 2.93436688E-02 1.16101365E-02-3.20782425E-02-2.39680205E-02

 7  9-2.38058977E-02-2.56653904E-02-2.79775085E-02-6.94034445E-03 1.61306081E-02

 7 10-1.19809284E-03 3.88112829E-03-1.57700511E-02-6.17772283E-02-9.07825204E-02

 8  1 2.13033949E-03 4.12921500E-03-2.88722296E-02 5.93356210E-02 1.21900883E-02

 8  2-9.96696627E-03-6.69339315E-02 1.14100288E-01 1.13188537E-02-1.22503066E-01

 8  3-1.33109762E-01 1.32887805E-01 3.23555318E-02-1.60929995E-02 1.21697051E-02

 8  4 3.92329447E-03-3.45019299E-02-2.83442404E-03 1.34033575E-02 3.20473065E-02

 8  5 5.45387801E-02 5.62719028E-02-2.15902334E-01-3.66698553E-02-1.63517759E-01

 8  6 8.80506465E-02-3.49848273E-01-5.22839180E-02-1.92366868E-01 6.87564464E-02

 8  7-4.43801753E-01-8.29050940E-02-1.08466915E-02 1.40654059E-02-3.21871439E-03

 8  8 9.88217400E-03-5.82676126E-03 8.10566210E-03 8.71952826E-02 1.71854819E-01

 8  9 1.50407148E-01-4.96094155E-02-1.22280775E-01-1.41939012E-01 3.75904698E-02

 8 10 5.11022917E-02 3.46805179E-02 6.20229047E-02 1.60638142E-01 1.77763971E-01

 9  1-8.93330037E-03-1.46694330E-02-6.75859323E-02 1.55631391E-03-1.66079821E-01

 9  2 4.56972414E-02-1.24918431E-01 8.84031312E-03-2.77335135E-01 2.28869804E-02

 9  3-1.54525264E-01 1.79741897E-02-4.36261107E-01-1.88777737E-03 5.44944041E-03

 9  4-3.56166304E-03 2.87785749E-03-2.11850693E-03 8.54845905E-03-2.98451460E-02

 9  5-5.09313044E-02 8.33214734E-02-2.96086168E-02 1.63299116E-01 1.40748286E-01

 9  6 1.19670267E-01-5.51940698E-02 2.56094423E-01 3.70903210E-01 1.31119178E-01

 9  7-8.64022400E-02 2.99040966E-01-1.68085722E-03 9.67950010E-03-7.99864288E-03

 9  8-6.32329100E-03 1.95423287E-02-1.39544444E-02-1.57436788E-02-1.16314196E-02

 9  9 2.04025498E-02-1.80527699E-02-1.84263845E-02-1.60320921E-02-7.57055466E-02

 9 10-1.68246279E-01-8.97894260E-02-2.50868305E-02-1.88534568E-02 7.83963584E-02

 $END  

 

Sample 5. Optimization of singlet-triplet ISC in methylene molecule at SA-MCSCF level using Lagrange function based code and semi-numerical gradients.

 

 

 $CONTRL

  SCFTYP=MCSCF RUNTYP=OPTIMIZE inttyp=hondo icut=11 NZVAR=3 COORD=ZMT

 $END

 $SYSTEM TIMLIM=6000 MWORDS=10 nojac=1 kdiag=0 $END

 $BASIS  GBASIS=n31 ngauss=6 ndfunc=1 npfunc=1 $END

 $GUESS  GUESS=HUCKEL $END

 $moorth nostf=1 nozero=1 tole=0 tolz=0 syms=1 symden=1 symvec=1 $end

 $trans  cuttrf=1d-13 $end

 $DET   

  NCORE=3 NACT=2 NELS=2 NSTATE=2 PURES=.F. CVGTOL=1D-8 GROUP=C1

  wstate(1)=1,1

 $END

 $MCSCF  CISTEP=ALDET ISTATE=1 acurcy=1d-7 ENGTOL=1.0d-12 NTRACK=2 $END

 $mcaver jstate=2 conic=2 $end

 $track  tol=1.4 freeze=.t. sticky=.f. update=.t. reset=.f. delciv=.f. $end

 $ZMAT   IZMAT(1)=1,1,2,   1,1,3,   2,2,1,3  $END

 $statpt nstep=50 opttol=1.0d-5 method=conic $end

 $conic  SHIFT0=1d-7 SHIFT=1D-7 $end

 $DATA

Methylene Singlet/Triplet ISC optimization

Cnv  2

 

C

H 1 rCH

H 1 rCH 2 aHOH

 

rCH=1.09

aHOH=99.0

 $END

 

An archive with output files can be found here.



See also:



 

Last updated: September 05, 2012