Jump to content

SJL1218 H2 OPT POP.LOG

From ChemWiki
Default is to use a total of   8 processors:
                               8 via shared-memory
                               1 via Linda
Entering Link 1 = C:\G09W\l1.exe PID=      9964.
 
Copyright (c) 1988,1990,1992,1993,1995,1998,2003,2009,2013,
           Gaussian, Inc.  All Rights Reserved.
 
This is part of the Gaussian(R) 09 program.  It is based on
the Gaussian(R) 03 system (copyright 2003, Gaussian, Inc.),
the Gaussian(R) 98 system (copyright 1998, Gaussian, Inc.),
the Gaussian(R) 94 system (copyright 1995, Gaussian, Inc.),
the Gaussian 92(TM) system (copyright 1992, Gaussian, Inc.),
the Gaussian 90(TM) system (copyright 1990, Gaussian, Inc.),
the Gaussian 88(TM) system (copyright 1988, Gaussian, Inc.),
the Gaussian 86(TM) system (copyright 1986, Carnegie Mellon
University), and the Gaussian 82(TM) system (copyright 1983,
Carnegie Mellon University). Gaussian is a federally registered
trademark of Gaussian, Inc.
 
This software contains proprietary and confidential information,
including trade secrets, belonging to Gaussian, Inc.
 
This software is provided under written license and may be
used, copied, transmitted, or stored only in accord with that
written license.
 
The following legend is applicable only to US Government
contracts under FAR:
 
                   RESTRICTED RIGHTS LEGEND
 
Use, reproduction and disclosure by the US Government is
subject to restrictions as set forth in subparagraphs (a)
and (c) of the Commercial Computer Software - Restricted
Rights clause in FAR 52.227-19.
 
Gaussian, Inc.
340 Quinnipiac St., Bldg. 40, Wallingford CT 06492
 
 
---------------------------------------------------------------
Warning -- This program may not be used in any manner that
competes with the business of Gaussian, Inc. or will provide
assistance to any competitor of Gaussian, Inc.  The licensee
of this program is prohibited from giving any competitor of
Gaussian, Inc. access to this program.  By using this program,
the user acknowledges that Gaussian, Inc. is engaged in the
business of creating and licensing software in the field of
computational chemistry and represents and warrants to the
licensee that it is not a competitor of Gaussian, Inc. and that
it will not use this program in any manner prohibited above.
---------------------------------------------------------------
 
Cite this work as:
Gaussian 09, Revision D.01,
M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, 
M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, 
G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, 
A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, 
M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, 
Y. Honda, O. Kitao, H. Nakai, T. Vreven, J. A. Montgomery, Jr., 
J. E. Peralta, F. Ogliaro, M. Bearpark, J. J. Heyd, E. Brothers, 
K. N. Kudin, V. N. Staroverov, T. Keith, R. Kobayashi, J. Normand, 
K. Raghavachari, A. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, 
M. Cossi, N. Rega, J. M. Millam, M. Klene, J. E. Knox, J. B. Cross, 
V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, 
O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. W. Ochterski, 
R. L. Martin, K. Morokuma, V. G. Zakrzewski, G. A. Voth, 
P. Salvador, J. J. Dannenberg, S. Dapprich, A. D. Daniels, 
O. Farkas, J. B. Foresman, J. V. Ortiz, J. Cioslowski, 
and D. J. Fox, Gaussian, Inc., Wallingford CT, 2013.

******************************************
Gaussian 09:  EM64W-G09RevD.01 13-Apr-2013
               12-Mar-2019 
******************************************
%chk=H:\justin\1styearlab\sjl1218_h2_opt2_pop.chk
Default route:  MaxDisk=10GB
----------------------------------------------------------------------
# opt freq b3lyp/6-31g(d,p) geom=connectivity integral=grid=ultrafine 
pop=(full,nbo)
----------------------------------------------------------------------
1/14=-1,18=20,19=15,26=4,38=1,57=2/1,3;
2/9=110,12=2,17=6,18=5,40=1/2;
3/5=1,6=6,7=101,11=2,16=1,25=1,30=1,71=1,74=-5,75=-5/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=3,28=1,40=1/1,7;
7//1,2,3,16;
1/14=-1,18=20,19=15,26=4/3(2);
2/9=110/2;
99//99;
2/9=110/2;
3/5=1,6=6,7=101,11=2,16=1,25=1,30=1,71=1,74=-5,75=-5/1,2,3;
4/5=5,16=3,69=1/1;
5/5=2,38=5/2;
7//1,2,3,16;
1/14=-1,18=20,19=15,26=4/3(-5);
2/9=110/2;
6/7=3,19=2,28=1,40=1/1,7;
99/9=1/99;
---------------
h2 optimisation
---------------
Symbolic Z-matrix:
Charge =  0 Multiplicity = 1
H                     0.        0.        0.37 
H                     0.        0.       -0.37 

GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Initialization pass.
                          ----------------------------
                          !    Initial Parameters    !
                          ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,2)                  0.74           estimate D2E/DX2                !
--------------------------------------------------------------------------------
Trust Radius=3.00D-01 FncErr=1.00D-07 GrdErr=1.00D-06
Number of steps in this run=     20 maximum allowed number of steps=    100.
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
                         Input orientation:                          
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
     1          1           0        0.000000    0.000000    0.370000
     2          1           0        0.000000    0.000000   -0.370000
---------------------------------------------------------------------
Stoichiometry    H2
Framework group  D*H[C*(H.H)]
Deg. of freedom     1
Full point group                 D*H     NOp   8
Largest Abelian subgroup         D2H     NOp   8
Largest concise Abelian subgroup C2      NOp   2
                        Standard orientation:                         
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
     1          1           0        0.000000    0.000000    0.370000
     2          1           0        0.000000    0.000000   -0.370000
---------------------------------------------------------------------
Rotational constants (GHZ):      0.0000000   1831.4649389   1831.4649389
Standard basis: 6-31G(d,p) (6D, 7F)
There are     3 symmetry adapted cartesian basis functions of AG  symmetry.
There are     0 symmetry adapted cartesian basis functions of B1G symmetry.
There are     1 symmetry adapted cartesian basis functions of B2G symmetry.
There are     1 symmetry adapted cartesian basis functions of B3G symmetry.
There are     0 symmetry adapted cartesian basis functions of AU  symmetry.
There are     3 symmetry adapted cartesian basis functions of B1U symmetry.
There are     1 symmetry adapted cartesian basis functions of B2U symmetry.
There are     1 symmetry adapted cartesian basis functions of B3U symmetry.
There are     3 symmetry adapted basis functions of AG  symmetry.
There are     0 symmetry adapted basis functions of B1G symmetry.
There are     1 symmetry adapted basis functions of B2G symmetry.
There are     1 symmetry adapted basis functions of B3G symmetry.
There are     0 symmetry adapted basis functions of AU  symmetry.
There are     3 symmetry adapted basis functions of B1U symmetry.
There are     1 symmetry adapted basis functions of B2U symmetry.
There are     1 symmetry adapted basis functions of B3U symmetry.
   10 basis functions,    14 primitive gaussians,    10 cartesian basis functions
    1 alpha electrons        1 beta electrons
      nuclear repulsion energy         0.7151043359 Hartrees.
NAtoms=    2 NActive=    2 NUniq=    1 SFac= 4.00D+00 NAtFMM=   60 NAOKFM=F Big=F
Integral buffers will be    131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
One-electron integrals computed using PRISM.
NBasis=    10 RedAO= T EigKep=  1.36D-01  NBF=     3     0     1     1     0     3     1     1
NBsUse=    10 1.00D-06 EigRej= -1.00D+00 NBFU=     3     0     1     1     0     3     1     1
ExpMin= 1.61D-01 ExpMax= 1.87D+01 ExpMxC= 1.87D+01 IAcc=3 IRadAn=         5 AccDes= 0.00D+00
Harris functional with IExCor=  402 and IRadAn=       5 diagonalized for initial guess.
HarFok:  IExCor=  402 AccDes= 0.00D+00 IRadAn=         5 IDoV= 1 UseB2=F ITyADJ=14
ICtDFT=  3500011 ScaDFX=  1.000000  1.000000  1.000000  1.000000
FoFCou: FMM=F IPFlag=           0 FMFlag=      100000 FMFlg1=           0
        NFxFlg=           0 DoJE=T BraDBF=F KetDBF=T FulRan=T
        wScrn=  0.000000 ICntrl=     500 IOpCl=  0 I1Cent=   200000004 NGrid=           0
        NMat0=    1 NMatS0=      1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
Petite list used in FoFCou.
Initial guess orbital symmetries:
      Occupied  (SGG)
      Virtual   (SGU) (SGG) (SGU) (PIU) (PIU) (SGG) (PIG) (PIG)
                (SGU)
The electronic state of the initial guess is 1-SGG.
Keep R1 ints in memory in symmetry-blocked form, NReq=883589.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Integral accuracy reduced to 1.0D-05 until final iterations.
Initial convergence to 1.0D-05 achieved.  Increase integral accuracy.
SCF Done:  E(RB3LYP) =  -1.17853402690     A.U. after    6 cycles
           NFock=  6  Conv=0.56D-09     -V/T= 2.0290
**********************************************************************
           Population analysis using the SCF density.
**********************************************************************
Orbital symmetries:
      Occupied  (SGG)
      Virtual   (SGU) (SGG) (SGU) (PIU) (PIU) (SGG) (PIG) (PIG)
                (SGU)
The electronic state is 1-SGG.
Alpha  occ. eigenvalues --   -0.43216
Alpha virt. eigenvalues --    0.10176   0.57599   1.03840   1.67444   1.67444
Alpha virt. eigenvalues --    2.38417   2.61878   2.61878   4.16664
    Molecular Orbital Coefficients:
                          1         2         3         4         5
                      (SGG)--O  (SGU)--V  (SGG)--V  (SGU)--V  (PIU)--V
    Eigenvalues --    -0.43216   0.10176   0.57599   1.03840   1.67444
  1 1   H  1S          0.32507   0.17164   0.74093  -0.90042   0.00000
  2        2S          0.26989   1.64077  -0.67779   1.48333   0.00000
  3        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.00000   0.61059
  5        3PZ        -0.01537   0.01055  -0.04015  -0.25232   0.00000
  6 2   H  1S          0.32507  -0.17164   0.74093   0.90042   0.00000
  7        2S          0.26989  -1.64077  -0.67779  -1.48333   0.00000
  8        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000   0.61059
 10        3PZ         0.01537   0.01055   0.04015  -0.25232   0.00000
                          6         7         8         9        10
                      (PIU)--V  (SGG)--V  (PIG)--V  (PIG)--V  (SGU)--V
    Eigenvalues --     1.67444   2.38417   2.61878   2.61878   4.16664
  1 1   H  1S          0.00000   0.39005   0.00000   0.00000  -1.43461
  2        2S          0.00000  -0.19286   0.00000   0.00000  -0.18637
  3        3PX         0.61059   0.00000   0.87112   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.87112   0.00000
  5        3PZ         0.00000   0.64761   0.00000   0.00000   1.58406
  6 2   H  1S          0.00000   0.39005   0.00000   0.00000   1.43461
  7        2S          0.00000  -0.19286   0.00000   0.00000   0.18637
  8        3PX         0.61059   0.00000  -0.87112   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000  -0.87112   0.00000
 10        3PZ         0.00000  -0.64761   0.00000   0.00000   1.58406
    Density Matrix:
                          1         2         3         4         5
  1 1   H  1S          0.21135
  2        2S          0.17547   0.14569
  3        3PX         0.00000   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.00000
  5        3PZ        -0.00999  -0.00830   0.00000   0.00000   0.00047
  6 2   H  1S          0.21135   0.17547   0.00000   0.00000  -0.00999
  7        2S          0.17547   0.14569   0.00000   0.00000  -0.00830
  8        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00999   0.00830   0.00000   0.00000  -0.00047
                          6         7         8         9        10
  6 2   H  1S          0.21135
  7        2S          0.17547   0.14569
  8        3PX         0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00999   0.00830   0.00000   0.00000   0.00047
   Full Mulliken population analysis:
                          1         2         3         4         5
  1 1   H  1S          0.21135
  2        2S          0.11551   0.14569
  3        3PX         0.00000   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.00000
  5        3PZ         0.00000   0.00000   0.00000   0.00000   0.00047
  6 2   H  1S          0.09623   0.08933   0.00000   0.00000   0.00470
  7        2S          0.08933   0.12443   0.00000   0.00000   0.00129
  8        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00470   0.00129   0.00000   0.00000   0.00019
                          6         7         8         9        10
  6 2   H  1S          0.21135
  7        2S          0.11551   0.14569
  8        3PX         0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00000   0.00000   0.00000   0.00000   0.00047
    Gross orbital populations:
                          1
  1 1   H  1S          0.51711
  2        2S          0.47624
  3        3PX         0.00000
  4        3PY         0.00000
  5        3PZ         0.00665
  6 2   H  1S          0.51711
  7        2S          0.47624
  8        3PX         0.00000
  9        3PY         0.00000
 10        3PZ         0.00665
         Condensed to atoms (all electrons):
              1          2
    1  H    0.588525   0.411475
    2  H    0.411475   0.588525
Mulliken charges:
              1
    1  H    0.000000
    2  H    0.000000
Sum of Mulliken charges =   0.00000
Mulliken charges with hydrogens summed into heavy atoms:
              1
Electronic spatial extent (au):  <R**2>=              5.1101
Charge=              0.0000 electrons
Dipole moment (field-independent basis, Debye):
   X=              0.0000    Y=              0.0000    Z=              0.0000  Tot=              0.0000
Quadrupole moment (field-independent basis, Debye-Ang):
  XX=             -2.0255   YY=             -2.0255   ZZ=             -1.5071
  XY=              0.0000   XZ=              0.0000   YZ=              0.0000
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
  XX=             -0.1728   YY=             -0.1728   ZZ=              0.3456
  XY=              0.0000   XZ=              0.0000   YZ=              0.0000
Octapole moment (field-independent basis, Debye-Ang**2):
 XXX=              0.0000  YYY=              0.0000  ZZZ=              0.0000  XYY=              0.0000
 XXY=              0.0000  XXZ=              0.0000  XZZ=              0.0000  YZZ=              0.0000
 YYZ=              0.0000  XYZ=              0.0000
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=             -1.8704 YYYY=             -1.8704 ZZZZ=             -2.7539 XXXY=              0.0000
XXXZ=              0.0000 YYYX=              0.0000 YYYZ=              0.0000 ZZZX=              0.0000
ZZZY=              0.0000 XXYY=             -0.6235 XXZZ=             -0.7854 YYZZ=             -0.7854
XXYZ=              0.0000 YYXZ=              0.0000 ZZXY=              0.0000
N-N= 7.151043359324D-01 E-N=-3.651704855992D+00  KE= 1.145276304114D+00
Symmetry AG   KE= 1.145276304114D+00
Symmetry B1G  KE= 0.000000000000D+00
Symmetry B2G  KE= 2.272325598493D-34
Symmetry B3G  KE= 2.272325598493D-34
Symmetry AU   KE= 0.000000000000D+00
Symmetry B1U  KE= 8.594876516240D-33
Symmetry B2U  KE= 6.501636522888D-35
Symmetry B3U  KE= 6.501636522888D-35
Orbital energies and kinetic energies (alpha):
                                1                 2
  1         (SGG)--O         -0.432158          0.572638
  2         (SGU)--V          0.101759          0.492676
  3         (SGG)--V          0.575986          1.406130
  4         (SGU)--V          1.038396          2.259771
  5         (PIU)--V          1.674442          2.449079
  6         (PIU)--V          1.674442          2.449079
  7         (SGG)--V          2.384170          3.280561
  8         (PIG)--V          2.618778          3.362504
  9         (PIG)--V          2.618778          3.362504
 10         (SGU)--V          4.166643          5.437219
Total kinetic energy from orbitals= 1.145276304114D+00
******************************Gaussian NBO Version 3.1******************************
            N A T U R A L   A T O M I C   O R B I T A L   A N D
         N A T U R A L   B O N D   O R B I T A L   A N A L Y S I S
******************************Gaussian NBO Version 3.1******************************
      /RESON  / : Allow strongly delocalized NBO set
Analyzing the SCF density
Job title: h2 optimisation                                                 
Storage needed:       370 in NPA,       439 in NBO ( 805306340 available)


NATURAL POPULATIONS:  Natural atomic orbital occupancies 
                                                         
  NAO  Atom  No  lang   Type(AO)    Occupancy      Energy
----------------------------------------------------------
    1    H    1  S      Val( 1S)     0.99944      -0.03191
    2    H    1  S      Ryd( 2S)     0.00000       0.71694
    3    H    1  px     Ryd( 2p)     0.00000       2.14661
    4    H    1  py     Ryd( 2p)     0.00000       2.14661
    5    H    1  pz     Ryd( 2p)     0.00056       3.23237
    6    H    2  S      Val( 1S)     0.99944      -0.03191
    7    H    2  S      Ryd( 2S)     0.00000       0.71694
    8    H    2  px     Ryd( 2p)     0.00000       2.14661
    9    H    2  py     Ryd( 2p)     0.00000       2.14661
   10    H    2  pz     Ryd( 2p)     0.00056       3.23237


Summary of Natural Population Analysis:                  
                                                         
                                      Natural Population 
               Natural  -----------------------------------------------
   Atom  No    Charge         Core      Valence    Rydberg      Total
-----------------------------------------------------------------------
     H    1    0.00000      0.00000     0.99944    0.00056     1.00000
     H    2    0.00000      0.00000     0.99944    0.00056     1.00000
=======================================================================
  * Total *    0.00000      0.00000     1.99888    0.00112     2.00000
                                Natural Population      
--------------------------------------------------------
  Valence                    1.99888 ( 99.9439% of   2)
  Natural Minimal Basis      1.99888 ( 99.9439% of   2)
  Natural Rydberg Basis      0.00112 (  0.0561% of   2)
--------------------------------------------------------
   Atom  No          Natural Electron Configuration
----------------------------------------------------------------------------
     H    1            1S( 1.00)
     H    2            1S( 1.00)


NATURAL BOND ORBITAL ANALYSIS:
                      Occupancies       Lewis Structure    Low   High
          Occ.    -------------------  -----------------   occ   occ
 Cycle   Thresh.   Lewis   Non-Lewis     CR  BD  3C  LP    (L)   (NL)   Dev
=============================================================================
  1(1)    1.90     2.00000   0.00000      0   1   0   0     0      0    0.00
-----------------------------------------------------------------------------
Structure accepted: No low occupancy Lewis orbitals
--------------------------------------------------------
  Valence Lewis             2.00000 (100.000% of   2)
 ==================       ============================
  Total Lewis               2.00000 (100.000% of   2)
 -----------------------------------------------------
  Valence non-Lewis         0.00000 (  0.000% of   2)
  Rydberg non-Lewis         0.00000 (  0.000% of   2)
 ==================       ============================
  Total non-Lewis           0.00000 (  0.000% of   2)
--------------------------------------------------------


      (Occupancy)   Bond orbital/ Coefficients/ Hybrids
---------------------------------------------------------------------------------
    1. (2.00000) BD ( 1) H   1 - H   2  
               ( 50.00%)   0.7071* H   1 s( 99.94%)p 0.00(  0.06%)
                                           0.9997  0.0000  0.0000  0.0000 -0.0237
               ( 50.00%)   0.7071* H   2 s( 99.94%)p 0.00(  0.06%)
                                           0.9997  0.0000  0.0000  0.0000  0.0237
    2. (0.00000) RY*( 1) H   1           s(100.00%)
    3. (0.00000) RY*( 2) H   1           s(  0.00%)p 1.00(100.00%)
    4. (0.00000) RY*( 3) H   1           s(  0.00%)p 1.00(100.00%)
    5. (0.00000) RY*( 4) H   1           s(  0.06%)p99.99( 99.94%)
    6. (0.00000) RY*( 1) H   2           s(100.00%)
    7. (0.00000) RY*( 2) H   2           s(  0.00%)p 1.00(100.00%)
    8. (0.00000) RY*( 3) H   2           s(  0.00%)p 1.00(100.00%)
    9. (0.00000) RY*( 4) H   2           s(  0.06%)p99.99( 99.94%)
   10. (0.00000) BD*( 1) H   1 - H   2  
               ( 50.00%)   0.7071* H   1 s( 99.94%)p 0.00(  0.06%)
               ( 50.00%)  -0.7071* H   2 s( 99.94%)p 0.00(  0.06%)


NHO Directionality and "Bond Bending" (deviations from line of nuclear centers)
        [Thresholds for printing:  angular deviation  >  1.0 degree]
                                   hybrid p-character > 25.0%
                                   orbital occupancy  >  0.10e
                               Line of Centers        Hybrid 1              Hybrid 2
                               ---------------  -------------------   ------------------
               NBO               Theta   Phi    Theta   Phi    Dev    Theta   Phi    Dev
========================================================================================
   None exceeding thresholds


Second Order Perturbation Theory Analysis of Fock Matrix in NBO Basis
    Threshold for printing:   0.50 kcal/mol
                                                                             E(2)  E(j)-E(i) F(i,j)
        Donor NBO (i)                     Acceptor NBO (j)                 kcal/mol   a.u.    a.u. 
===================================================================================================
within unit  1
      None above threshold


Natural Bond Orbitals (Summary):
                                                           Principal Delocalizations
          NBO                        Occupancy    Energy   (geminal,vicinal,remote)
====================================================================================
Molecular unit  1  (H2)
    1. BD (   1) H   1 - H   2          2.00000    -0.43216   
    2. RY*(   1) H   1                  0.00000     0.71694   
    3. RY*(   2) H   1                  0.00000     2.14661   
    4. RY*(   3) H   1                  0.00000     2.14661   
    5. RY*(   4) H   1                  0.00000     3.22573   
    6. RY*(   1) H   2                  0.00000     0.71694   
    7. RY*(   2) H   2                  0.00000     2.14661   
    8. RY*(   3) H   2                  0.00000     2.14661   
    9. RY*(   4) H   2                  0.00000     3.22573   
   10. BD*(   1) H   1 - H   2          0.00000     0.38161   
      -------------------------------
             Total Lewis    2.00000  (100.0000%)
       Valence non-Lewis    0.00000  (  0.0000%)
       Rydberg non-Lewis    0.00000  (  0.0000%)
      -------------------------------
           Total unit  1    2.00000  (100.0000%)
          Charge unit  1    0.00000
Calling FoFJK, ICntrl=      2127 FMM=F ISym2X=1 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0.
***** Axes restored to original set *****
-------------------------------------------------------------------
Center     Atomic                   Forces (Hartrees/Bohr)
Number     Number              X              Y              Z
-------------------------------------------------------------------
     1        1           0.000000000    0.000000000    0.002025594
     2        1           0.000000000    0.000000000   -0.002025594
-------------------------------------------------------------------
Cartesian Forces:  Max     0.002025594 RMS     0.001169477
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4.
Internal  Forces:  Max     0.002025594 RMS     0.002025594
Search for a local minimum.
Step number   1 out of a maximum of   20
All quantities printed in internal units (Hartrees-Bohrs-Radians)
Mixed Optimization -- RFO/linear search
Second derivative matrix not updated -- first step.
The second derivative matrix:
                         R1
          R1           0.39139
ITU=  0
    Eigenvalues ---    0.39139
RFO step:  Lambda=-1.04828338D-05 EMin= 3.91394209D-01
Linear search not attempted -- first point.
Iteration  1 RMS(Cart)=  0.00365941 RMS(Int)=  0.00000000
Iteration  2 RMS(Cart)=  0.00000000 RMS(Int)=  0.00000000
ClnCor:  largest displacement from symmetrization is 1.58D-19 for atom     1.
Variable       Old X    -DE/DX   Delta X   Delta X   Delta X     New X
                                (Linear)    (Quad)   (Total)
   R1        1.39840   0.00203   0.00000   0.00518   0.00518   1.40357
        Item               Value     Threshold  Converged?
Maximum Force            0.002026     0.000450     NO 
RMS     Force            0.002026     0.000300     NO 
Maximum Displacement     0.002588     0.001800     NO 
RMS     Displacement     0.003659     0.001200     NO 
Predicted change in Energy=-5.241557D-06
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
                         Input orientation:                          
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
     1          1           0        0.000000    0.000000    0.371369
     2          1           0        0.000000    0.000000   -0.371369
---------------------------------------------------------------------
Stoichiometry    H2
Framework group  D*H[C*(H.H)]
Deg. of freedom     1
Full point group                 D*H     NOp   8
Largest Abelian subgroup         D2H     NOp   8
Largest concise Abelian subgroup C2      NOp   2
                        Standard orientation:                         
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
     1          1           0        0.000000    0.000000    0.371369
     2          1           0        0.000000    0.000000   -0.371369
---------------------------------------------------------------------
Rotational constants (GHZ):      0.0000000   1817.9840451   1817.9840451
Standard basis: 6-31G(d,p) (6D, 7F)
There are     3 symmetry adapted cartesian basis functions of AG  symmetry.
There are     0 symmetry adapted cartesian basis functions of B1G symmetry.
There are     1 symmetry adapted cartesian basis functions of B2G symmetry.
There are     1 symmetry adapted cartesian basis functions of B3G symmetry.
There are     0 symmetry adapted cartesian basis functions of AU  symmetry.
There are     3 symmetry adapted cartesian basis functions of B1U symmetry.
There are     1 symmetry adapted cartesian basis functions of B2U symmetry.
There are     1 symmetry adapted cartesian basis functions of B3U symmetry.
There are     3 symmetry adapted basis functions of AG  symmetry.
There are     0 symmetry adapted basis functions of B1G symmetry.
There are     1 symmetry adapted basis functions of B2G symmetry.
There are     1 symmetry adapted basis functions of B3G symmetry.
There are     0 symmetry adapted basis functions of AU  symmetry.
There are     3 symmetry adapted basis functions of B1U symmetry.
There are     1 symmetry adapted basis functions of B2U symmetry.
There are     1 symmetry adapted basis functions of B3U symmetry.
   10 basis functions,    14 primitive gaussians,    10 cartesian basis functions
    1 alpha electrons        1 beta electrons
      nuclear repulsion energy         0.7124676349 Hartrees.
NAtoms=    2 NActive=    2 NUniq=    1 SFac= 4.00D+00 NAtFMM=   60 NAOKFM=F Big=F
Integral buffers will be    131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
One-electron integrals computed using PRISM.
NBasis=    10 RedAO= T EigKep=  1.37D-01  NBF=     3     0     1     1     0     3     1     1
NBsUse=    10 1.00D-06 EigRej= -1.00D+00 NBFU=     3     0     1     1     0     3     1     1
Initial guess from the checkpoint file:  "H:\justin\1styearlab\sjl1218_h2_opt2_pop.chk"
B after Tr=     0.000000    0.000000    0.000000
        Rot=    1.000000    0.000000    0.000000    0.000000 Ang=   0.00 deg.
Initial guess orbital symmetries:
      Occupied  (SGG)
      Virtual   (SGG) (SGG) (PIG) (PIG) (SGU) (SGU) (SGU) (PIU)
                (PIU)
Keep R1 ints in memory in symmetry-blocked form, NReq=883589.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
SCF Done:  E(RB3LYP) =  -1.17853935531     A.U. after    4 cycles
           NFock=  4  Conv=0.10D-10     -V/T= 2.0313
Calling FoFJK, ICntrl=      2127 FMM=F ISym2X=1 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0.
***** Axes restored to original set *****
-------------------------------------------------------------------
Center     Atomic                   Forces (Hartrees/Bohr)
Number     Number              X              Y              Z
-------------------------------------------------------------------
     1        1           0.000000000    0.000000000    0.000039427
     2        1           0.000000000    0.000000000   -0.000039427
-------------------------------------------------------------------
Cartesian Forces:  Max     0.000039427 RMS     0.000022763
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Using GEDIIS/GDIIS optimizer.
Internal  Forces:  Max     0.000039427 RMS     0.000039427
Search for a local minimum.
Step number   2 out of a maximum of   20
All quantities printed in internal units (Hartrees-Bohrs-Radians)
Mixed Optimization -- En-DIIS/RFO-DIIS
Update second derivatives using D2CorX and points    1    2
DE= -5.33D-06 DEPred=-5.24D-06 R= 1.02D+00
TightC=F SS=  1.41D+00  RLast= 5.18D-03 DXNew= 5.0454D-01 1.5526D-02
Trust test= 1.02D+00 RLast= 5.18D-03 DXMaxT set to 3.00D-01
The second derivative matrix:
                         R1
          R1           0.38379
ITU=  1  0
Use linear search instead of GDIIS.
    Eigenvalues ---    0.38379
RFO step:  Lambda= 0.00000000D+00 EMin= 3.83786282D-01
Quartic linear search produced a step of  0.02003.
Iteration  1 RMS(Cart)=  0.00007330 RMS(Int)=  0.00000000
Iteration  2 RMS(Cart)=  0.00000000 RMS(Int)=  0.00000000
ClnCor:  largest displacement from symmetrization is 3.17D-21 for atom     1.
Variable       Old X    -DE/DX   Delta X   Delta X   Delta X     New X
                                (Linear)    (Quad)   (Total)
   R1        1.40357   0.00004   0.00010   0.00000   0.00010   1.40368
        Item               Value     Threshold  Converged?
Maximum Force            0.000039     0.000450     YES
RMS     Force            0.000039     0.000300     YES
Maximum Displacement     0.000052     0.001800     YES
RMS     Displacement     0.000073     0.001200     YES
Predicted change in Energy=-2.025007D-09
Optimization completed.
   -- Stationary point found.
                          ----------------------------
                          !   Optimized Parameters   !
                          ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,2)                  0.7427         -DE/DX =    0.0                 !
--------------------------------------------------------------------------------
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
                         Input orientation:                          
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
     1          1           0        0.000000    0.000000    0.371369
     2          1           0        0.000000    0.000000   -0.371369
---------------------------------------------------------------------
Stoichiometry    H2
Framework group  D*H[C*(H.H)]
Deg. of freedom     1
Full point group                 D*H     NOp   8
Largest Abelian subgroup         D2H     NOp   8
Largest concise Abelian subgroup C2      NOp   2
                        Standard orientation:                         
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
     1          1           0        0.000000    0.000000    0.371369
     2          1           0        0.000000    0.000000   -0.371369
---------------------------------------------------------------------
Rotational constants (GHZ):      0.0000000   1817.9840451   1817.9840451
**********************************************************************
           Population analysis using the SCF density.
**********************************************************************
Orbital symmetries:
      Occupied  (SGG)
      Virtual   (SGU) (SGG) (SGU) (PIU) (PIU) (SGG) (PIG) (PIG)
                (SGU)
The electronic state is 1-SGG.
Alpha  occ. eigenvalues --   -0.43159
Alpha virt. eigenvalues --    0.10098   0.57698   1.03614   1.67470   1.67470
Alpha virt. eigenvalues --    2.37898   2.61573   2.61573   4.15313
    Molecular Orbital Coefficients:
                          1         2         3         4         5
                      (SGG)--O  (SGU)--V  (SGG)--V  (SGU)--V  (PIU)--V
    Eigenvalues --    -0.43159   0.10098   0.57698   1.03614   1.67470
  1 1   H  1S          0.32476   0.17232   0.74181  -0.90214   0.00000
  2        2S          0.27047   1.63371  -0.67787   1.48132   0.00000
  3        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.00000   0.61121
  5        3PZ        -0.01536   0.01056  -0.04072  -0.24982   0.00000
  6 2   H  1S          0.32476  -0.17232   0.74181   0.90214   0.00000
  7        2S          0.27047  -1.63371  -0.67787  -1.48132   0.00000
  8        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000   0.61121
 10        3PZ         0.01536   0.01056   0.04072  -0.24982   0.00000
                          6         7         8         9        10
                      (PIU)--V  (SGG)--V  (PIG)--V  (PIG)--V  (SGU)--V
    Eigenvalues --     1.67470   2.37898   2.61573   2.61573   4.15313
  1 1   H  1S          0.00000   0.38958   0.00000   0.00000  -1.41514
  2        2S          0.00000  -0.19262   0.00000   0.00000  -0.18970
  3        3PX         0.61121   0.00000   0.00000   0.86934   0.00000
  4        3PY         0.00000   0.00000   0.86934   0.00000   0.00000
  5        3PZ         0.00000   0.64661   0.00000   0.00000   1.57600
  6 2   H  1S          0.00000   0.38958   0.00000   0.00000   1.41514
  7        2S          0.00000  -0.19262   0.00000   0.00000   0.18970
  8        3PX         0.61121   0.00000   0.00000  -0.86934   0.00000
  9        3PY         0.00000   0.00000  -0.86934   0.00000   0.00000
 10        3PZ         0.00000  -0.64661   0.00000   0.00000   1.57600
    Density Matrix:
                          1         2         3         4         5
  1 1   H  1S          0.21093
  2        2S          0.17568   0.14631
  3        3PX         0.00000   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.00000
  5        3PZ        -0.00998  -0.00831   0.00000   0.00000   0.00047
  6 2   H  1S          0.21093   0.17568   0.00000   0.00000  -0.00998
  7        2S          0.17568   0.14631   0.00000   0.00000  -0.00831
  8        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00998   0.00831   0.00000   0.00000  -0.00047
                          6         7         8         9        10
  6 2   H  1S          0.21093
  7        2S          0.17568   0.14631
  8        3PX         0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00998   0.00831   0.00000   0.00000   0.00047
   Full Mulliken population analysis:
                          1         2         3         4         5
  1 1   H  1S          0.21093
  2        2S          0.11565   0.14631
  3        3PX         0.00000   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.00000
  5        3PZ         0.00000   0.00000   0.00000   0.00000   0.00047
  6 2   H  1S          0.09552   0.08926   0.00000   0.00000   0.00468
  7        2S          0.08926   0.12482   0.00000   0.00000   0.00129
  8        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00468   0.00129   0.00000   0.00000   0.00019
                          6         7         8         9        10
  6 2   H  1S          0.21093
  7        2S          0.11565   0.14631
  8        3PX         0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00000   0.00000   0.00000   0.00000   0.00047
    Gross orbital populations:
                          1
  1 1   H  1S          0.51604
  2        2S          0.47733
  3        3PX         0.00000
  4        3PY         0.00000
  5        3PZ         0.00663
  6 2   H  1S          0.51604
  7        2S          0.47733
  8        3PX         0.00000
  9        3PY         0.00000
 10        3PZ         0.00663
         Condensed to atoms (all electrons):
              1          2
    1  H    0.589008   0.410992
    2  H    0.410992   0.589008
Mulliken charges:
              1
    1  H    0.000000
    2  H    0.000000
Sum of Mulliken charges =   0.00000
Mulliken charges with hydrogens summed into heavy atoms:
              1
Electronic spatial extent (au):  <R**2>=              5.1229
Charge=              0.0000 electrons
Dipole moment (field-independent basis, Debye):
   X=              0.0000    Y=              0.0000    Z=              0.0000  Tot=              0.0000
Quadrupole moment (field-independent basis, Debye-Ang):
  XX=             -2.0290   YY=             -2.0290   ZZ=             -1.5076
  XY=              0.0000   XZ=              0.0000   YZ=              0.0000
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
  XX=             -0.1738   YY=             -0.1738   ZZ=              0.3476
  XY=              0.0000   XZ=              0.0000   YZ=              0.0000
Octapole moment (field-independent basis, Debye-Ang**2):
 XXX=              0.0000  YYY=              0.0000  ZZZ=              0.0000  XYY=              0.0000
 XXY=              0.0000  XXZ=              0.0000  XZZ=              0.0000  YZZ=              0.0000
 YYZ=              0.0000  XYZ=              0.0000
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=             -1.8758 YYYY=             -1.8758 ZZZZ=             -2.7677 XXXY=              0.0000
XXXZ=              0.0000 YYYX=              0.0000 YYYZ=              0.0000 ZZZX=              0.0000
ZZZY=              0.0000 XXYY=             -0.6253 XXZZ=             -0.7887 YYZZ=             -0.7887
XXYZ=              0.0000 YYXZ=              0.0000 ZZXY=              0.0000
N-N= 7.124676349026D-01 E-N=-3.645572510598D+00  KE= 1.142762958424D+00
Symmetry AG   KE= 1.142762958424D+00
Symmetry B1G  KE= 0.000000000000D+00
Symmetry B2G  KE= 2.251668286257D-34
Symmetry B3G  KE= 2.251668286257D-34
Symmetry AU   KE= 0.000000000000D+00
Symmetry B1U  KE= 1.137098012505D-32
Symmetry B2U  KE= 6.480579541981D-35
Symmetry B3U  KE= 6.480579541981D-35
Orbital energies and kinetic energies (alpha):
                                1                 2
  1         (SGG)--O         -0.431588          0.571381
  2         (SGU)--V          0.100982          0.493348
  3         (SGG)--V          0.576985          1.407549
  4         (SGU)--V          1.036138          2.257827
  5         (PIU)--V          1.674704          2.448647
  6         (PIU)--V          1.674704          2.448647
  7         (SGG)--V          2.378979          3.275210
  8         (PIG)--V          2.615726          3.359638
  9         (PIG)--V          2.615726          3.359638
 10         (SGU)--V          4.153129          5.416951
Total kinetic energy from orbitals= 1.142762958424D+00
******************************Gaussian NBO Version 3.1******************************
            N A T U R A L   A T O M I C   O R B I T A L   A N D
         N A T U R A L   B O N D   O R B I T A L   A N A L Y S I S
******************************Gaussian NBO Version 3.1******************************
      /RESON  / : Allow strongly delocalized NBO set
Analyzing the SCF density
Job title: h2 optimisation                                                 
Storage needed:       370 in NPA,       439 in NBO ( 805306340 available)


NATURAL POPULATIONS:  Natural atomic orbital occupancies 
                                                         
  NAO  Atom  No  lang   Type(AO)    Occupancy      Energy
----------------------------------------------------------
    1    H    1  S      Val( 1S)     0.99944      -0.03358
    2    H    1  S      Ryd( 2S)     0.00000       0.71747
    3    H    1  px     Ryd( 2p)     0.00000       2.14521
    4    H    1  py     Ryd( 2p)     0.00000       2.14521
    5    H    1  pz     Ryd( 2p)     0.00056       3.22342
    6    H    2  S      Val( 1S)     0.99944      -0.03358
    7    H    2  S      Ryd( 2S)     0.00000       0.71747
    8    H    2  px     Ryd( 2p)     0.00000       2.14521
    9    H    2  py     Ryd( 2p)     0.00000       2.14521
   10    H    2  pz     Ryd( 2p)     0.00056       3.22342


Summary of Natural Population Analysis:                  
                                                         
                                      Natural Population 
               Natural  -----------------------------------------------
   Atom  No    Charge         Core      Valence    Rydberg      Total
-----------------------------------------------------------------------
     H    1    0.00000      0.00000     0.99944    0.00056     1.00000
     H    2    0.00000      0.00000     0.99944    0.00056     1.00000
=======================================================================
  * Total *    0.00000      0.00000     1.99888    0.00112     2.00000
                                Natural Population      
--------------------------------------------------------
  Valence                    1.99888 ( 99.9438% of   2)
  Natural Minimal Basis      1.99888 ( 99.9438% of   2)
  Natural Rydberg Basis      0.00112 (  0.0562% of   2)
--------------------------------------------------------
   Atom  No          Natural Electron Configuration
----------------------------------------------------------------------------
     H    1            1S( 1.00)
     H    2            1S( 1.00)


NATURAL BOND ORBITAL ANALYSIS:
                      Occupancies       Lewis Structure    Low   High
          Occ.    -------------------  -----------------   occ   occ
 Cycle   Thresh.   Lewis   Non-Lewis     CR  BD  3C  LP    (L)   (NL)   Dev
=============================================================================
  1(1)    1.90     2.00000   0.00000      0   1   0   0     0      0    0.00
-----------------------------------------------------------------------------
Structure accepted: No low occupancy Lewis orbitals
--------------------------------------------------------
  Valence Lewis             2.00000 (100.000% of   2)
 ==================       ============================
  Total Lewis               2.00000 (100.000% of   2)
 -----------------------------------------------------
  Valence non-Lewis         0.00000 (  0.000% of   2)
  Rydberg non-Lewis         0.00000 (  0.000% of   2)
 ==================       ============================
  Total non-Lewis           0.00000 (  0.000% of   2)
--------------------------------------------------------


      (Occupancy)   Bond orbital/ Coefficients/ Hybrids
---------------------------------------------------------------------------------
    1. (2.00000) BD ( 1) H   1 - H   2  
               ( 50.00%)   0.7071* H   1 s( 99.94%)p 0.00(  0.06%)
                                           0.9997  0.0000  0.0000  0.0000 -0.0237
               ( 50.00%)   0.7071* H   2 s( 99.94%)p 0.00(  0.06%)
                                           0.9997  0.0000  0.0000  0.0000  0.0237
    2. (0.00000) RY*( 1) H   1           s(100.00%)
    3. (0.00000) RY*( 2) H   1           s(  0.00%)p 1.00(100.00%)
    4. (0.00000) RY*( 3) H   1           s(  0.00%)p 1.00(100.00%)
    5. (0.00000) RY*( 4) H   1           s(  0.06%)p99.99( 99.94%)
    6. (0.00000) RY*( 1) H   2           s(100.00%)
    7. (0.00000) RY*( 2) H   2           s(  0.00%)p 1.00(100.00%)
    8. (0.00000) RY*( 3) H   2           s(  0.00%)p 1.00(100.00%)
    9. (0.00000) RY*( 4) H   2           s(  0.06%)p99.99( 99.94%)
   10. (0.00000) BD*( 1) H   1 - H   2  
               ( 50.00%)   0.7071* H   1 s( 99.94%)p 0.00(  0.06%)
               ( 50.00%)  -0.7071* H   2 s( 99.94%)p 0.00(  0.06%)


NHO Directionality and "Bond Bending" (deviations from line of nuclear centers)
        [Thresholds for printing:  angular deviation  >  1.0 degree]
                                   hybrid p-character > 25.0%
                                   orbital occupancy  >  0.10e
                               Line of Centers        Hybrid 1              Hybrid 2
                               ---------------  -------------------   ------------------
               NBO               Theta   Phi    Theta   Phi    Dev    Theta   Phi    Dev
========================================================================================
   None exceeding thresholds


Second Order Perturbation Theory Analysis of Fock Matrix in NBO Basis
    Threshold for printing:   0.50 kcal/mol
                                                                             E(2)  E(j)-E(i) F(i,j)
        Donor NBO (i)                     Acceptor NBO (j)                 kcal/mol   a.u.    a.u. 
===================================================================================================
within unit  1
      None above threshold


Natural Bond Orbitals (Summary):
                                                           Principal Delocalizations
          NBO                        Occupancy    Energy   (geminal,vicinal,remote)
====================================================================================
Molecular unit  1  (H2)
    1. BD (   1) H   1 - H   2          2.00000    -0.43159   
    2. RY*(   1) H   1                  0.00000     0.71747   
    3. RY*(   2) H   1                  0.00000     2.14521   
    4. RY*(   3) H   1                  0.00000     2.14521   
    5. RY*(   4) H   1                  0.00000     3.21685   
    6. RY*(   1) H   2                  0.00000     0.71747   
    7. RY*(   2) H   2                  0.00000     2.14521   
    8. RY*(   3) H   2                  0.00000     2.14521   
    9. RY*(   4) H   2                  0.00000     3.21685   
   10. BD*(   1) H   1 - H   2          0.00000     0.37758   
      -------------------------------
             Total Lewis    2.00000  (100.0000%)
       Valence non-Lewis    0.00000  (  0.0000%)
       Rydberg non-Lewis    0.00000  (  0.0000%)
      -------------------------------
           Total unit  1    2.00000  (100.0000%)
          Charge unit  1    0.00000
1|1| IMPERIAL COLLEGE-SKCH-232A-010|FOpt|RB3LYP|6-31G(d,p)|H2|SJL1218|
12-Mar-2019|0||# opt freq b3lyp/6-31g(d,p) geom=connectivity integral=
grid=ultrafine pop=(full,nbo)||h2 optimisation||0,1|H,0.,0.,0.37136929
64|H,0.,0.,-0.3713692964||Version=EM64W-G09RevD.01|State=1-SGG|HF=-1.1
785394|RMSD=9.952e-012|RMSF=2.276e-005|Dipole=0.,0.,0.|Quadrupole=-0.1
292199,-0.1292199,0.2584398,0.,0.,0.|PG=D*H [C*(H1.H1)]||@


The juvenile sea squirt wanders through the sea searching
for a suitable rock or hunk of coral to cling to and make
its home for life.  For this task it has a rudimentary
nervous system.  When it finds its spot and takes root, it
doesn't need its brain any more so it eats it.
It's rather like getting tenure.
                              -- source unknown
Job cpu time:       0 days  0 hours  0 minutes 34.0 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1
Normal termination of Gaussian 09 at Tue Mar 12 16:35:39 2019.
Link1:  Proceeding to internal job step number  2.
----------------------------------------------------------------------
#N Geom=AllCheck Guess=TCheck SCRF=Check GenChk RB3LYP/6-31G(d,p) Freq
----------------------------------------------------------------------
1/10=4,29=7,30=1,38=1,40=1/1,3;
2/12=2,40=1/2;
3/5=1,6=6,7=101,11=2,14=-4,16=1,25=1,30=1,70=2,71=2,74=-5,75=-5,116=1,140=1/1,2,3;
4/5=101/1;
5/5=2,98=1/2;
8/6=4,10=90,11=11/1;
11/6=1,8=1,9=11,15=111,16=1/1,2,10;
10/6=1/2;
6/7=3,18=1,28=1,40=1/1,7;
7/8=1,10=1,25=1/1,2,3,16;
1/10=4,30=1/3;
99//99;
Structure from the checkpoint file:  "H:\justin\1styearlab\sjl1218_h2_opt2_pop.chk"
---------------
h2 optimisation
---------------
Charge =  0 Multiplicity = 1
Redundant internal coordinates found in file.
H,0,0.,0.,0.3713692964
H,0,0.,0.,-0.3713692964
Recover connectivity data from disk.
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Initialization pass.
                          ----------------------------
                          !    Initial Parameters    !
                          ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,2)                  0.7427         calculate D2E/DX2 analytically  !
--------------------------------------------------------------------------------
Trust Radius=3.00D-01 FncErr=1.00D-07 GrdErr=1.00D-07
Number of steps in this run=      2 maximum allowed number of steps=      2.
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
                         Input orientation:                          
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
     1          1           0        0.000000    0.000000    0.371369
     2          1           0        0.000000    0.000000   -0.371369
---------------------------------------------------------------------
Stoichiometry    H2
Framework group  D*H[C*(H.H)]
Deg. of freedom     1
Full point group                 D*H     NOp   8
Largest Abelian subgroup         D2H     NOp   8
Largest concise Abelian subgroup C2      NOp   2
                        Standard orientation:                         
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
     1          1           0        0.000000    0.000000    0.371369
     2          1           0        0.000000    0.000000   -0.371369
---------------------------------------------------------------------
Rotational constants (GHZ):      0.0000000   1817.9840451   1817.9840451
Standard basis: 6-31G(d,p) (6D, 7F)
There are     3 symmetry adapted cartesian basis functions of AG  symmetry.
There are     0 symmetry adapted cartesian basis functions of B1G symmetry.
There are     1 symmetry adapted cartesian basis functions of B2G symmetry.
There are     1 symmetry adapted cartesian basis functions of B3G symmetry.
There are     0 symmetry adapted cartesian basis functions of AU  symmetry.
There are     3 symmetry adapted cartesian basis functions of B1U symmetry.
There are     1 symmetry adapted cartesian basis functions of B2U symmetry.
There are     1 symmetry adapted cartesian basis functions of B3U symmetry.
There are     3 symmetry adapted basis functions of AG  symmetry.
There are     0 symmetry adapted basis functions of B1G symmetry.
There are     1 symmetry adapted basis functions of B2G symmetry.
There are     1 symmetry adapted basis functions of B3G symmetry.
There are     0 symmetry adapted basis functions of AU  symmetry.
There are     3 symmetry adapted basis functions of B1U symmetry.
There are     1 symmetry adapted basis functions of B2U symmetry.
There are     1 symmetry adapted basis functions of B3U symmetry.
   10 basis functions,    14 primitive gaussians,    10 cartesian basis functions
    1 alpha electrons        1 beta electrons
      nuclear repulsion energy         0.7124676349 Hartrees.
NAtoms=    2 NActive=    2 NUniq=    1 SFac= 4.00D+00 NAtFMM=   60 NAOKFM=F Big=F
Integral buffers will be    131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
One-electron integrals computed using PRISM.
NBasis=    10 RedAO= T EigKep=  1.37D-01  NBF=     3     0     1     1     0     3     1     1
NBsUse=    10 1.00D-06 EigRej= -1.00D+00 NBFU=     3     0     1     1     0     3     1     1
Initial guess from the checkpoint file:  "H:\justin\1styearlab\sjl1218_h2_opt2_pop.chk"
B after Tr=     0.000000    0.000000    0.000000
        Rot=    1.000000    0.000000    0.000000    0.000000 Ang=   0.00 deg.
Initial guess orbital symmetries:
      Occupied  (SGG)
      Virtual   (SGU) (SGG) (SGU) (PIU) (PIU) (SGG) (PIG) (PIG)
                (SGU)
Keep R1 ints in memory in symmetry-blocked form, NReq=883589.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Skip diagonalization as Alpha Fock matrix is already diagonal.
SCF Done:  E(RB3LYP) =  -1.17853935531     A.U. after    1 cycles
           NFock=  1  Conv=0.00D+00     -V/T= 2.0313
DoSCS=F DFT=T ScalE2(SS,OS)=  1.000000  1.000000
Range of M.O.s used for correlation:     1    10
NBasis=    10 NAE=     1 NBE=     1 NFC=     0 NFV=     0
NROrb=     10 NOA=     1 NOB=     1 NVA=     9 NVB=     9
Symmetrizing basis deriv contribution to polar:
IMax=3 JMax=2 DiffMx= 0.00D+00
G2DrvN: will do     3 centers at a time, making    1 passes.
Calling FoFCou, ICntrl=  3107 FMM=F I1Cent=   0 AccDes= 0.00D+00.
End of G2Drv F.D. properties file   721 does not exist.
End of G2Drv F.D. properties file   722 does not exist.
End of G2Drv F.D. properties file   788 does not exist.
         IDoAtm=11
         Differentiating once with respect to electric field.
               with respect to dipole field.
         Differentiating once with respect to nuclear coordinates.
Keep R1 ints in memory in symmetry-blocked form, NReq=860931.
         There are     6 degrees of freedom in the 1st order CPHF.  IDoFFX=4 NUNeed=     6.
     6 vectors produced by pass  0 Test12= 1.50D-16 1.67D-08 XBig12= 3.62D+00 1.90D+00.
AX will form     6 AO Fock derivatives at one time.
     2 vectors produced by pass  1 Test12= 1.50D-16 1.67D-08 XBig12= 4.65D-02 2.14D-01.
     1 vectors produced by pass  2 Test12= 1.50D-16 1.67D-08 XBig12= 1.15D-06 1.07D-03.
InvSVY:  IOpt=1 It=  1 EMax= 3.64D-18
Solved reduced A of dimension     9 with     6 vectors.
Isotropic polarizability for W=    0.000000        2.51 Bohr**3.
End of Minotr F.D. properties file   721 does not exist.
End of Minotr F.D. properties file   722 does not exist.
End of Minotr F.D. properties file   788 does not exist.
**********************************************************************
           Population analysis using the SCF density.
**********************************************************************
Orbital symmetries:
      Occupied  (SGG)
      Virtual   (SGU) (SGG) (SGU) (PIU) (PIU) (SGG) (PIG) (PIG)
                (SGU)
The electronic state is 1-SGG.
Alpha  occ. eigenvalues --   -0.43159
Alpha virt. eigenvalues --    0.10098   0.57698   1.03614   1.67470   1.67470
Alpha virt. eigenvalues --    2.37898   2.61573   2.61573   4.15313
    Molecular Orbital Coefficients:
                          1         2         3         4         5
                      (SGG)--O  (SGU)--V  (SGG)--V  (SGU)--V  (PIU)--V
    Eigenvalues --    -0.43159   0.10098   0.57698   1.03614   1.67470
  1 1   H  1S          0.32476   0.17232   0.74181  -0.90214   0.00000
  2        2S          0.27047   1.63371  -0.67787   1.48132   0.00000
  3        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.00000   0.61121
  5        3PZ        -0.01536   0.01056  -0.04072  -0.24982   0.00000
  6 2   H  1S          0.32476  -0.17232   0.74181   0.90214   0.00000
  7        2S          0.27047  -1.63371  -0.67787  -1.48132   0.00000
  8        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000   0.61121
 10        3PZ         0.01536   0.01056   0.04072  -0.24982   0.00000
                          6         7         8         9        10
                      (PIU)--V  (SGG)--V  (PIG)--V  (PIG)--V  (SGU)--V
    Eigenvalues --     1.67470   2.37898   2.61573   2.61573   4.15313
  1 1   H  1S          0.00000   0.38958   0.00000   0.00000  -1.41514
  2        2S          0.00000  -0.19262   0.00000   0.00000  -0.18970
  3        3PX         0.61121   0.00000   0.86934   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.86934   0.00000
  5        3PZ         0.00000   0.64661   0.00000   0.00000   1.57600
  6 2   H  1S          0.00000   0.38958   0.00000   0.00000   1.41514
  7        2S          0.00000  -0.19262   0.00000   0.00000   0.18970
  8        3PX         0.61121   0.00000  -0.86934   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000  -0.86934   0.00000
 10        3PZ         0.00000  -0.64661   0.00000   0.00000   1.57600
    Density Matrix:
                          1         2         3         4         5
  1 1   H  1S          0.21093
  2        2S          0.17568   0.14631
  3        3PX         0.00000   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.00000
  5        3PZ        -0.00998  -0.00831   0.00000   0.00000   0.00047
  6 2   H  1S          0.21093   0.17568   0.00000   0.00000  -0.00998
  7        2S          0.17568   0.14631   0.00000   0.00000  -0.00831
  8        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00998   0.00831   0.00000   0.00000  -0.00047
                          6         7         8         9        10
  6 2   H  1S          0.21093
  7        2S          0.17568   0.14631
  8        3PX         0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00998   0.00831   0.00000   0.00000   0.00047
   Full Mulliken population analysis:
                          1         2         3         4         5
  1 1   H  1S          0.21093
  2        2S          0.11565   0.14631
  3        3PX         0.00000   0.00000   0.00000
  4        3PY         0.00000   0.00000   0.00000   0.00000
  5        3PZ         0.00000   0.00000   0.00000   0.00000   0.00047
  6 2   H  1S          0.09552   0.08926   0.00000   0.00000   0.00468
  7        2S          0.08926   0.12482   0.00000   0.00000   0.00129
  8        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00468   0.00129   0.00000   0.00000   0.00019
                          6         7         8         9        10
  6 2   H  1S          0.21093
  7        2S          0.11565   0.14631
  8        3PX         0.00000   0.00000   0.00000
  9        3PY         0.00000   0.00000   0.00000   0.00000
 10        3PZ         0.00000   0.00000   0.00000   0.00000   0.00047
    Gross orbital populations:
                          1
  1 1   H  1S          0.51604
  2        2S          0.47733
  3        3PX         0.00000
  4        3PY         0.00000
  5        3PZ         0.00663
  6 2   H  1S          0.51604
  7        2S          0.47733
  8        3PX         0.00000
  9        3PY         0.00000
 10        3PZ         0.00663
         Condensed to atoms (all electrons):
              1          2
    1  H    0.589008   0.410992
    2  H    0.410992   0.589008
Mulliken charges:
              1
    1  H    0.000000
    2  H    0.000000
Sum of Mulliken charges =   0.00000
Mulliken charges with hydrogens summed into heavy atoms:
              1
APT charges:
              1
    1  H    0.000000
    2  H    0.000000
Sum of APT charges =   0.00000
APT charges with hydrogens summed into heavy atoms:
              1
Electronic spatial extent (au):  <R**2>=              5.1229
Charge=              0.0000 electrons
Dipole moment (field-independent basis, Debye):
   X=              0.0000    Y=              0.0000    Z=              0.0000  Tot=              0.0000
Quadrupole moment (field-independent basis, Debye-Ang):
  XX=             -2.0290   YY=             -2.0290   ZZ=             -1.5076
  XY=              0.0000   XZ=              0.0000   YZ=              0.0000
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
  XX=             -0.1738   YY=             -0.1738   ZZ=              0.3476
  XY=              0.0000   XZ=              0.0000   YZ=              0.0000
Octapole moment (field-independent basis, Debye-Ang**2):
 XXX=              0.0000  YYY=              0.0000  ZZZ=              0.0000  XYY=              0.0000
 XXY=              0.0000  XXZ=              0.0000  XZZ=              0.0000  YZZ=              0.0000
 YYZ=              0.0000  XYZ=              0.0000
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=             -1.8758 YYYY=             -1.8758 ZZZZ=             -2.7677 XXXY=              0.0000
XXXZ=              0.0000 YYYX=              0.0000 YYYZ=              0.0000 ZZZX=              0.0000
ZZZY=              0.0000 XXYY=             -0.6253 XXZZ=             -0.7887 YYZZ=             -0.7887
XXYZ=              0.0000 YYXZ=              0.0000 ZZXY=              0.0000
N-N= 7.124676349026D-01 E-N=-3.645572510598D+00  KE= 1.142762958424D+00
Symmetry AG   KE= 1.142762958424D+00
Symmetry B1G  KE= 0.000000000000D+00
Symmetry B2G  KE= 2.251668286257D-34
Symmetry B3G  KE= 2.251668286257D-34
Symmetry AU   KE= 0.000000000000D+00
Symmetry B1U  KE= 1.715305121897D-32
Symmetry B2U  KE= 6.480579541981D-35
Symmetry B3U  KE= 6.480579541981D-35
Orbital energies and kinetic energies (alpha):
                                1                 2
  1         (SGG)--O         -0.431588          0.571381
  2         (SGU)--V          0.100982          0.493348
  3         (SGG)--V          0.576985          1.407549
  4         (SGU)--V          1.036138          2.257827
  5         (PIU)--V          1.674704          2.448647
  6         (PIU)--V          1.674704          2.448647
  7         (SGG)--V          2.378979          3.275210
  8         (PIG)--V          2.615726          3.359638
  9         (PIG)--V          2.615726          3.359638
 10         (SGU)--V          4.153129          5.416951
Total kinetic energy from orbitals= 1.142762958424D+00
 Exact polarizability:   0.582   0.000   0.582   0.000   0.000   6.364
Approx polarizability:   0.599   0.000   0.599   0.000   0.000   7.763
******************************Gaussian NBO Version 3.1******************************
            N A T U R A L   A T O M I C   O R B I T A L   A N D
         N A T U R A L   B O N D   O R B I T A L   A N A L Y S I S
******************************Gaussian NBO Version 3.1******************************
      /RESON  / : Allow strongly delocalized NBO set
Analyzing the SCF density
Job title: h2 optimisation                                                 
Storage needed:       370 in NPA,       439 in NBO ( 805306340 available)


NATURAL POPULATIONS:  Natural atomic orbital occupancies 
                                                         
  NAO  Atom  No  lang   Type(AO)    Occupancy      Energy
----------------------------------------------------------
    1    H    1  S      Val( 1S)     0.99944      -0.03358
    2    H    1  S      Ryd( 2S)     0.00000       0.71747
    3    H    1  px     Ryd( 2p)     0.00000       2.14521
    4    H    1  py     Ryd( 2p)     0.00000       2.14521
    5    H    1  pz     Ryd( 2p)     0.00056       3.22342
    6    H    2  S      Val( 1S)     0.99944      -0.03358
    7    H    2  S      Ryd( 2S)     0.00000       0.71747
    8    H    2  px     Ryd( 2p)     0.00000       2.14521
    9    H    2  py     Ryd( 2p)     0.00000       2.14521
   10    H    2  pz     Ryd( 2p)     0.00056       3.22342


Summary of Natural Population Analysis:                  
                                                         
                                      Natural Population 
               Natural  -----------------------------------------------
   Atom  No    Charge         Core      Valence    Rydberg      Total
-----------------------------------------------------------------------
     H    1    0.00000      0.00000     0.99944    0.00056     1.00000
     H    2    0.00000      0.00000     0.99944    0.00056     1.00000
=======================================================================
  * Total *    0.00000      0.00000     1.99888    0.00112     2.00000
                                Natural Population      
--------------------------------------------------------
  Valence                    1.99888 ( 99.9438% of   2)
  Natural Minimal Basis      1.99888 ( 99.9438% of   2)
  Natural Rydberg Basis      0.00112 (  0.0562% of   2)
--------------------------------------------------------
   Atom  No          Natural Electron Configuration
----------------------------------------------------------------------------
     H    1            1S( 1.00)
     H    2            1S( 1.00)


NATURAL BOND ORBITAL ANALYSIS:
                      Occupancies       Lewis Structure    Low   High
          Occ.    -------------------  -----------------   occ   occ
 Cycle   Thresh.   Lewis   Non-Lewis     CR  BD  3C  LP    (L)   (NL)   Dev
=============================================================================
  1(1)    1.90     2.00000   0.00000      0   1   0   0     0      0    0.00
-----------------------------------------------------------------------------
Structure accepted: No low occupancy Lewis orbitals
--------------------------------------------------------
  Valence Lewis             2.00000 (100.000% of   2)
 ==================       ============================
  Total Lewis               2.00000 (100.000% of   2)
 -----------------------------------------------------
  Valence non-Lewis         0.00000 (  0.000% of   2)
  Rydberg non-Lewis         0.00000 (  0.000% of   2)
 ==================       ============================
  Total non-Lewis           0.00000 (  0.000% of   2)
--------------------------------------------------------


      (Occupancy)   Bond orbital/ Coefficients/ Hybrids
---------------------------------------------------------------------------------
    1. (2.00000) BD ( 1) H   1 - H   2  
               ( 50.00%)   0.7071* H   1 s( 99.94%)p 0.00(  0.06%)
                                           0.9997  0.0000  0.0000  0.0000 -0.0237
               ( 50.00%)   0.7071* H   2 s( 99.94%)p 0.00(  0.06%)
                                           0.9997  0.0000  0.0000  0.0000  0.0237
    2. (0.00000) RY*( 1) H   1           s(100.00%)
    3. (0.00000) RY*( 2) H   1           s(  0.00%)p 1.00(100.00%)
    4. (0.00000) RY*( 3) H   1           s(  0.00%)p 1.00(100.00%)
    5. (0.00000) RY*( 4) H   1           s(  0.06%)p99.99( 99.94%)
    6. (0.00000) RY*( 1) H   2           s(100.00%)
    7. (0.00000) RY*( 2) H   2           s(  0.00%)p 1.00(100.00%)
    8. (0.00000) RY*( 3) H   2           s(  0.00%)p 1.00(100.00%)
    9. (0.00000) RY*( 4) H   2           s(  0.06%)p99.99( 99.94%)
   10. (0.00000) BD*( 1) H   1 - H   2  
               ( 50.00%)   0.7071* H   1 s( 99.94%)p 0.00(  0.06%)
               ( 50.00%)  -0.7071* H   2 s( 99.94%)p 0.00(  0.06%)


NHO Directionality and "Bond Bending" (deviations from line of nuclear centers)
        [Thresholds for printing:  angular deviation  >  1.0 degree]
                                   hybrid p-character > 25.0%
                                   orbital occupancy  >  0.10e
                               Line of Centers        Hybrid 1              Hybrid 2
                               ---------------  -------------------   ------------------
               NBO               Theta   Phi    Theta   Phi    Dev    Theta   Phi    Dev
========================================================================================
   None exceeding thresholds


Second Order Perturbation Theory Analysis of Fock Matrix in NBO Basis
    Threshold for printing:   0.50 kcal/mol
                                                                             E(2)  E(j)-E(i) F(i,j)
        Donor NBO (i)                     Acceptor NBO (j)                 kcal/mol   a.u.    a.u. 
===================================================================================================
within unit  1
      None above threshold


Natural Bond Orbitals (Summary):
                                                           Principal Delocalizations
          NBO                        Occupancy    Energy   (geminal,vicinal,remote)
====================================================================================
Molecular unit  1  (H2)
    1. BD (   1) H   1 - H   2          2.00000    -0.43159   
    2. RY*(   1) H   1                  0.00000     0.71747   
    3. RY*(   2) H   1                  0.00000     2.14521   
    4. RY*(   3) H   1                  0.00000     2.14521   
    5. RY*(   4) H   1                  0.00000     3.21685   
    6. RY*(   1) H   2                  0.00000     0.71747   
    7. RY*(   2) H   2                  0.00000     2.14521   
    8. RY*(   3) H   2                  0.00000     2.14521   
    9. RY*(   4) H   2                  0.00000     3.21685   
   10. BD*(   1) H   1 - H   2          0.00000     0.37758   
      -------------------------------
             Total Lewis    2.00000  (100.0000%)
       Valence non-Lewis    0.00000  (  0.0000%)
       Rydberg non-Lewis    0.00000  (  0.0000%)
      -------------------------------
           Total unit  1    2.00000  (100.0000%)
          Charge unit  1    0.00000
Calling FoFJK, ICntrl=    100127 FMM=F ISym2X=1 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0.
Full mass-weighted force constant matrix:
Low frequencies ---  -38.3800  -38.3800   -0.0001    0.0000    0.0000 4466.4614
Diagonal vibrational polarizability:
       0.0000000       0.0000000       0.0000000
Harmonic frequencies (cm**-1), IR intensities (KM/Mole), Raman scattering
activities (A**4/AMU), depolarization ratios for plane and unpolarized
incident light, reduced masses (AMU), force constants (mDyne/A),
and normal coordinates:
                     1
                    SGG
Frequencies --   4466.4614
Red. masses --      1.0078
Frc consts  --     11.8457
IR Inten    --      0.0000
 Atom  AN      X      Y      Z
    1   1     0.00   0.00   0.71
    2   1     0.00   0.00  -0.71
-------------------
- Thermochemistry -
-------------------
Temperature   298.150 Kelvin.  Pressure   1.00000 Atm.
Atom     1 has atomic number  1 and mass   1.00783
Atom     2 has atomic number  1 and mass   1.00783
Molecular mass:     2.01565 amu.
Principal axes and moments of inertia in atomic units:
                          1         2         3
    Eigenvalues --     0.00000   0.99272   0.99272
          X            0.00000   1.00000   0.00000
          Y            0.00000   0.00000   1.00000
          Z            1.00000   0.00000   0.00000
This molecule is a prolate symmetric top.
Rotational symmetry number  2.
Rotational temperature (Kelvin)     87.24937
Rotational constant (GHZ):        1817.984045
Zero-point vibrational energy      26715.4 (Joules/Mol)
                                   6.38513 (Kcal/Mol)
Vibrational temperatures:   6426.23
         (Kelvin)

Zero-point correction=                           0.010175 (Hartree/Particle)
Thermal correction to Energy=                    0.012536
Thermal correction to Enthalpy=                  0.013480
Thermal correction to Gibbs Free Energy=        -0.001312
Sum of electronic and zero-point Energies=             -1.168364
Sum of electronic and thermal Energies=                -1.166004
Sum of electronic and thermal Enthalpies=              -1.165059
Sum of electronic and thermal Free Energies=           -1.179851

                    E (Thermal)             CV                S
                     KCal/Mol        Cal/Mol-Kelvin    Cal/Mol-Kelvin
Total                    7.866              4.968             31.132
Electronic               0.000              0.000              0.000
Translational            0.889              2.981             28.080
Rotational               0.592              1.987              3.052
Vibrational              6.385              0.000              0.000
                      Q            Log10(Q)             Ln(Q)
Total Bot       0.401232D+01          0.603395          1.389369
Total V=0       0.192185D+06          5.283720         12.166216
Vib (Bot)       0.208773D-04         -4.680325        -10.776847
Vib (V=0)       0.100000D+01          0.000000          0.000000
Electronic      0.100000D+01          0.000000          0.000000
Translational   0.112481D+06          5.051078         11.630537
Rotational      0.170861D+01          0.232643          0.535679
***** Axes restored to original set *****
-------------------------------------------------------------------
Center     Atomic                   Forces (Hartrees/Bohr)
Number     Number              X              Y              Z
-------------------------------------------------------------------
     1        1           0.000000000    0.000000000    0.000039427
     2        1           0.000000000    0.000000000   -0.000039427
-------------------------------------------------------------------
Cartesian Forces:  Max     0.000039427 RMS     0.000022763
FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4.
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Internal  Forces:  Max     0.000039427 RMS     0.000039427
Search for a local minimum.
Step number   1 out of a maximum of    2
All quantities printed in internal units (Hartrees-Bohrs-Radians)
Second derivative matrix not updated -- analytic derivatives used.
The second derivative matrix:
                         R1
          R1           0.38043
ITU=  0
    Eigenvalues ---    0.38043
Angle between quadratic step and forces=   0.00 degrees.
Linear search not attempted -- first point.
Iteration  1 RMS(Cart)=  0.00007328 RMS(Int)=  0.00000000
Iteration  2 RMS(Cart)=  0.00000000 RMS(Int)=  0.00000000
ClnCor:  largest displacement from symmetrization is 3.17D-21 for atom     1.
Variable       Old X    -DE/DX   Delta X   Delta X   Delta X     New X
                                (Linear)    (Quad)   (Total)
   R1        1.40357   0.00004   0.00000   0.00010   0.00010   1.40368
        Item               Value     Threshold  Converged?
Maximum Force            0.000039     0.000450     YES
RMS     Force            0.000039     0.000300     YES
Maximum Displacement     0.000052     0.001800     YES
RMS     Displacement     0.000073     0.001200     YES
Predicted change in Energy=-2.043043D-09
Optimization completed.
   -- Stationary point found.
                          ----------------------------
                          !   Optimized Parameters   !
                          ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,2)                  0.7427         -DE/DX =    0.0                 !
--------------------------------------------------------------------------------
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
1|1| IMPERIAL COLLEGE-SKCH-232A-010|Freq|RB3LYP|6-31G(d,p)|H2|SJL1218|
12-Mar-2019|0||#N Geom=AllCheck Guess=TCheck SCRF=Check GenChk RB3LYP/
6-31G(d,p) Freq||h2 optimisation||0,1|H,0.,0.,0.3713692964|H,0.,0.,-0.
3713692964||Version=EM64W-G09RevD.01|State=1-SGG|HF=-1.1785394|RMSD=0.
000e+000|RMSF=2.276e-005|ZeroPoint=0.0101753|Thermal=0.0125358|Dipole=
0.,0.,0.|DipoleDeriv=0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0
.,0.|Polar=0.5819176,0.,0.5819176,0.,0.,6.3644219|PG=D*H [C*(H1.H1)]|N
Imag=0||-0.00002809,0.,-0.00002809,0.,0.,0.38042913,0.00002809,0.,0.,-
0.00002809,0.,0.00002809,0.,0.,-0.00002809,0.,0.,-0.38042913,0.,0.,0.3
8042913||0.,0.,-0.00003943,0.,0.,0.00003943|||@


It is also a good rule not to put too much confidence in
experimental results until they have been confirmed by
Theory.
                            -- Sir Arthur Eddington
Job cpu time:       0 days  0 hours  0 minutes 23.0 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1
Normal termination of Gaussian 09 at Tue Mar 12 16:36:03 2019.