Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- All connections between all procs tested: SUCCESS
- argument 1 = h2o_raman_3.nw
- ============================== echo of input deck ==============================
- echo
- start h2o_raman_3
- title "TDDFT h2o BP86/aug-cc-pVTZ"
- geometry
- symmetry c1
- O 0.00000000 0.00000000 0.12982363
- H 0.75933475 0.00000000 -0.46621158
- H -0.75933475 0.00000000 -0.46621158
- end
- basis spherical
- O library aug-cc-pVTZ
- H library aug-cc-pVTZ
- end
- dft
- xc becke88 perdew86
- convergence energy 1d-8 nolevelshifting
- direct
- end
- hessian
- thresh 1d-8
- profile
- end
- #set prop:newaoresp 1
- property
- response 1 8.8559E-2
- damping 0.007
- end
- raman
- normal
- end
- task dft optimize
- task dft raman
- ================================================================================
- Northwest Computational Chemistry Package (NWChem) 6.3
- ------------------------------------------------------
- Environmental Molecular Sciences Laboratory
- Pacific Northwest National Laboratory
- Richland, WA 99352
- Copyright (c) 1994-2013
- Pacific Northwest National Laboratory
- Battelle Memorial Institute
- NWChem is an open-source computational chemistry package
- distributed under the terms of the
- Educational Community License (ECL) 2.0
- A copy of the license is included with this distribution
- in the LICENSE.TXT file
- ACKNOWLEDGMENT
- --------------
- This software and its documentation were developed at the
- EMSL at Pacific Northwest National Laboratory, a multiprogram
- national laboratory, operated for the U.S. Department of Energy
- by Battelle under Contract Number DE-AC05-76RL01830. Support
- for this work was provided by the Department of Energy Office
- of Biological and Environmental Research, Office of Basic
- Energy Sciences, and the Office of Advanced Scientific Computing.
- Job information
- ---------------
- hostname = arcen
- program = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/../bin/LINUX64/nwchem
- date = Fri Jun 21 02:32:06 2013
- compiled = Thu_Jun_20_21:48:15_2013
- source = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17
- nwchem branch = Development
- nwchem revision = 24322
- ga revision = 10386
- input = h2o_raman_3.nw
- prefix = h2o_raman_3.
- data base = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.db
- status = startup
- nproc = 4
- time left = -1s
- Memory information
- ------------------
- heap = 13107201 doubles = 100.0 Mbytes
- stack = 13107201 doubles = 100.0 Mbytes
- global = 26214400 doubles = 200.0 Mbytes (distinct from heap & stack)
- total = 52428802 doubles = 400.0 Mbytes
- verify = yes
- hardfail = no
- Directory information
- ---------------------
- 0 permanent = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir
- 0 scratch = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir
- NWChem Input Module
- -------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- --------------------------
- Scaling coordinates for geometry "geometry" by 1.889725989
- (inverse scale = 0.529177249)
- Turning off AUTOSYM since
- SYMMETRY directive was detected!
- ------
- auto-z
- ------
- Geometry "geometry" -> ""
- -------------------------
- Output coordinates in angstroms (scale by 1.889725989 to convert to a.u.)
- No. Tag Charge X Y Z
- ---- ---------------- ---------- -------------- -------------- --------------
- 1 O 8.0000 0.00000000 0.00000000 0.11920704
- 2 H 1.0000 0.75933475 0.00000000 -0.47682817
- 3 H 1.0000 -0.75933475 0.00000000 -0.47682817
- Atomic Mass
- -----------
- O 15.994910
- H 1.007825
- Effective nuclear repulsion energy (a.u.) 9.1194413826
- Nuclear Dipole moment (a.u.)
- ----------------------------
- X Y Z
- ---------------- ---------------- ----------------
- 0.0000000000 0.0000000000 0.0000000000
- Z-matrix (autoz)
- --------
- Units are Angstrom for bonds and degrees for angles
- Type Name I J K L M Value
- ----------- -------- ----- ----- ----- ----- ----- ----------
- 1 Stretch 1 2 0.96532
- 2 Stretch 1 3 0.96532
- 3 Bend 2 1 3 103.74015
- XYZ format geometry
- -------------------
- 3
- geometry
- O 0.00000000 0.00000000 0.11920704
- H 0.75933475 0.00000000 -0.47682817
- H -0.75933475 0.00000000 -0.47682817
- ==============================================================================
- internuclear distances
- ------------------------------------------------------------------------------
- center one | center two | atomic units | angstroms
- ------------------------------------------------------------------------------
- 2 H | 1 O | 1.82419 | 0.96532
- 3 H | 1 O | 1.82419 | 0.96532
- ------------------------------------------------------------------------------
- number of included internuclear distances: 2
- ==============================================================================
- ==============================================================================
- internuclear angles
- ------------------------------------------------------------------------------
- center 1 | center 2 | center 3 | degrees
- ------------------------------------------------------------------------------
- 2 H | 1 O | 3 H | 103.74
- ------------------------------------------------------------------------------
- number of included internuclear angles: 1
- ==============================================================================
- library name resolved from: environment
- library file name is: </home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/../src/basis/libraries/>
- Basis "ao basis" -> "" (spherical)
- -----
- O (Oxygen)
- ----------
- Exponent Coefficients
- -------------- ---------------------------------------------------------
- 1 S 1.53300000E+04 0.000508
- 1 S 2.29900000E+03 0.003929
- 1 S 5.22400000E+02 0.020243
- 1 S 1.47300000E+02 0.079181
- 1 S 4.75500000E+01 0.230687
- 1 S 1.67600000E+01 0.433118
- 1 S 6.20700000E+00 0.350260
- 1 S 6.88200000E-01 -0.008154
- 2 S 1.53300000E+04 -0.000115
- 2 S 2.29900000E+03 -0.000895
- 2 S 5.22400000E+02 -0.004636
- 2 S 1.47300000E+02 -0.018724
- 2 S 4.75500000E+01 -0.058463
- 2 S 1.67600000E+01 -0.136463
- 2 S 6.20700000E+00 -0.175740
- 2 S 6.88200000E-01 0.603418
- 3 S 1.75200000E+00 1.000000
- 4 S 2.38400000E-01 1.000000
- 5 S 7.37600000E-02 1.000000
- 6 P 3.44600000E+01 0.015928
- 6 P 7.74900000E+00 0.099740
- 6 P 2.28000000E+00 0.310492
- 7 P 7.15600000E-01 1.000000
- 8 P 2.14000000E-01 1.000000
- 9 P 5.97400000E-02 1.000000
- 10 D 2.31400000E+00 1.000000
- 11 D 6.45000000E-01 1.000000
- 12 D 2.14000000E-01 1.000000
- 13 F 1.42800000E+00 1.000000
- 14 F 5.00000000E-01 1.000000
- H (Hydrogen)
- ------------
- Exponent Coefficients
- -------------- ---------------------------------------------------------
- 1 S 3.38700000E+01 0.006068
- 1 S 5.09500000E+00 0.045308
- 1 S 1.15900000E+00 0.202822
- 2 S 3.25800000E-01 1.000000
- 3 S 1.02700000E-01 1.000000
- 4 S 2.52600000E-02 1.000000
- 5 P 1.40700000E+00 1.000000
- 6 P 3.88000000E-01 1.000000
- 7 P 1.02000000E-01 1.000000
- 8 D 1.05700000E+00 1.000000
- 9 D 2.47000000E-01 1.000000
- Summary of "ao basis" -> "" (spherical)
- ------------------------------------------------------------------------------
- Tag Description Shells Functions and Types
- ---------------- ------------------------------ ------ ---------------------
- O aug-cc-pVTZ 14 46 5s4p3d2f
- H aug-cc-pVTZ 9 23 4s3p2d
- perdew86 is a nonlocal functional; adding perdew81 local functional.
- NWChem Geometry Optimization
- ----------------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- maximum gradient threshold (gmax) = 0.000450
- rms gradient threshold (grms) = 0.000300
- maximum cartesian step threshold (xmax) = 0.001800
- rms cartesian step threshold (xrms) = 0.001200
- fixed trust radius (trust) = 0.300000
- maximum step size to saddle (sadstp) = 0.100000
- energy precision (eprec) = 5.0D-06
- maximum number of steps (nptopt) = 20
- initial hessian option (inhess) = 0
- line search option (linopt) = 1
- hessian update option (modupd) = 1
- saddle point option (modsad) = 0
- initial eigen-mode to follow (moddir) = 0
- initial variable to follow (vardir) = 0
- follow first negative mode (firstneg) = T
- apply conjugacy (opcg) = F
- source of zmatrix = autoz
- -------------------
- Energy Minimization
- -------------------
- Names of Z-matrix variables
- 1 2 3
- Variables with the same non-blank name are constrained to be equal
- Using diagonal initial Hessian
- Scaling for Hessian diagonals: bonds = 1.00 angles = 0.25 torsions = 0.10
- --------
- Step 0
- --------
- Geometry "geometry" -> "geometry"
- ---------------------------------
- Output coordinates in angstroms (scale by 1.889725989 to convert to a.u.)
- No. Tag Charge X Y Z
- ---- ---------------- ---------- -------------- -------------- --------------
- 1 O 8.0000 0.00000000 0.00000000 0.11920704
- 2 H 1.0000 0.75933475 0.00000000 -0.47682817
- 3 H 1.0000 -0.75933475 0.00000000 -0.47682817
- Atomic Mass
- -----------
- O 15.994910
- H 1.007825
- Effective nuclear repulsion energy (a.u.) 9.1194413826
- Nuclear Dipole moment (a.u.)
- ----------------------------
- X Y Z
- ---------------- ---------------- ----------------
- 0.0000000000 0.0000000000 0.0000000000
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- itol2e modified to match energy
- convergence criterion.
- tol_rho modified to match energy
- convergence criterion.
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 5.00D-04
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: medium
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 49 10.0 434
- H 0.35 45 12.0 434
- Grid pruning is: on
- Number of quadrature shells: 139
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Superposition of Atomic Density Guess
- -------------------------------------
- Sum of atomic energies: -75.78090571
- Non-variational initial energy
- ------------------------------
- Total energy = -75.930721
- 1-e energy = -121.602639
- 2-e energy = 36.552476
- HOMO = -0.479146
- LUMO = 0.009494
- Time after variat. SCF: 0.2
- Time prior to 1st pass: 0.2
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 7 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 13.02 13017841
- Stack Space remaining (MW): 13.11 13106626
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4209090826 -8.55D+01 6.82D-03 4.71D-01 1.0
- d= 0,ls=0.0,diis 2 -76.3451351812 7.58D-02 3.89D-03 1.31D+00 1.6
- d= 0,ls=0.0,diis 3 -76.4635023713 -1.18D-01 5.50D-04 3.25D-02 2.4
- d= 0,ls=0.0,diis 4 -76.4659494806 -2.45D-03 8.92D-05 8.00D-04 3.0
- d= 0,ls=0.0,diis 5 -76.4660251828 -7.57D-05 1.50D-05 8.01D-06 3.8
- d= 0,ls=0.0,diis 6 -76.4660260201 -8.37D-07 2.27D-06 1.17D-07 4.7
- d= 0,ls=0.0,diis 7 -76.4660260302 -1.01D-08 1.58D-07 1.26D-09 5.5
- d= 0,ls=0.0,diis 8 -76.4660260303 -1.22D-10 1.35D-08 4.86D-12 6.3
- Total DFT energy = -76.466026030303
- One electron energy = -122.862097984700
- Coulomb energy = 46.594399688616
- Exchange-Corr. energy = -9.317769116771
- Nuclear repulsion energy = 9.119441382552
- Numeric. integr. density = 10.000001897403
- Total iterative time = 6.1s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878098D+01
- MO Center= -3.5D-17, 1.1D-17, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.974179 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.339673D-01
- MO Center= -5.6D-15, 2.0D-15, -8.7D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.481521 1 O s 4 0.277365 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.850122D-01
- MO Center= 2.9D-15, -3.3D-16, -9.6D-02, r^2= 7.8D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.324514 1 O px 6 0.301665 1 O px
- 48 0.228415 2 H s 71 -0.228415 3 H s
- Vector 4 Occ=2.000000D+00 E=-3.479937D-01
- MO Center= -2.2D-16, 9.6D-15, 1.8D-01, r^2= 7.5D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.359025 1 O pz 4 0.347780 1 O s
- 8 0.320073 1 O pz 14 0.204348 1 O pz
- 2 0.201613 1 O s 48 -0.176063 2 H s
- 71 -0.176063 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.701046D-01
- MO Center= -1.5D-15, -2.5D-15, 7.2D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428090 1 O py 7 0.364992 1 O py
- 13 0.323481 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.172497D-02
- MO Center= 1.7D-13, -6.2D-13, -7.8D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.491133 1 O s 50 -0.467854 2 H s
- 73 -0.467854 3 H s 4 0.368247 1 O s
- 49 -0.354418 2 H s 72 -0.354418 3 H s
- 14 -0.191310 1 O pz
- Vector 7 Occ=0.000000D+00 E= 1.551591D-02
- MO Center= 3.9D-13, -1.8D-13, -7.6D-01, r^2= 1.5D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.668862 2 H s 73 -2.668862 3 H s
- 49 0.814387 2 H s 72 -0.814387 3 H s
- 15 -0.648799 1 O px 12 -0.342808 1 O px
- 48 0.159539 2 H s 71 -0.159539 3 H s
- Vector 8 Occ=0.000000D+00 E= 6.650085D-02
- MO Center= 5.5D-13, -3.8D-12, -5.4D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.427495 1 O s 49 -2.151189 2 H s
- 72 -2.151189 3 H s 4 1.380804 1 O s
- 17 -1.153835 1 O pz 50 -1.018077 2 H s
- 73 -1.018077 3 H s 57 0.679097 2 H px
- 80 -0.679097 3 H px 59 -0.544366 2 H pz
- Vector 9 Occ=0.000000D+00 E= 9.413783D-02
- MO Center= -9.7D-14, 5.2D-12, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.498969 1 O py 58 -0.313568 2 H py
- 81 -0.313568 3 H py
- Vector 10 Occ=0.000000D+00 E= 9.929854D-02
- MO Center= -4.4D-12, 1.3D-12, 1.6D-01, r^2= 8.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.444642 1 O pz 49 0.507814 2 H s
- 72 0.507814 3 H s 5 0.413544 1 O s
- 50 -0.378808 2 H s 73 -0.378808 3 H s
- 59 -0.330998 2 H pz 82 -0.330998 3 H pz
- 4 0.309097 1 O s 57 0.300298 2 H px
- Vector 11 Occ=0.000000D+00 E= 1.022251D-01
- MO Center= 3.1D-12, 5.7D-14, -8.3D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.940552 2 H s 73 -2.940552 3 H s
- 49 2.206957 2 H s 72 -2.206957 3 H s
- 15 -1.598312 1 O px 57 -0.939764 2 H px
- 80 -0.939764 3 H px 59 0.927590 2 H pz
- 82 -0.927590 3 H pz 12 -0.565438 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.479623D-01
- MO Center= 1.1D-12, 4.3D-13, 6.4D-01, r^2= 8.7D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.196865 1 O px 49 -2.806830 2 H s
- 72 2.806830 3 H s 50 -2.144384 2 H s
- 73 2.144384 3 H s 12 1.398485 1 O px
- 48 -0.949065 2 H s 71 0.949065 3 H s
- 31 0.583776 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.686730D-01
- MO Center= 2.3D-13, -3.0D-13, -1.9D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.369602 1 O s 49 -1.626885 2 H s
- 72 -1.626885 3 H s 59 -0.986314 2 H pz
- 82 -0.986314 3 H pz 57 -0.634402 2 H px
- 80 0.634402 3 H px 17 0.466293 1 O pz
- 4 0.437173 1 O s 30 -0.427721 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.915048D-01
- MO Center= 3.6D-13, -3.8D-13, -5.0D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 58 1.596935 2 H py 81 -1.596935 3 H py
- 28 -0.473639 1 O d -2 55 0.235063 2 H py
- 78 -0.235063 3 H py 65 -0.186953 2 H d -2
- 88 -0.186953 3 H d -2 66 0.168551 2 H d -1
- 89 -0.168551 3 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.138476D-01
- MO Center= -2.7D-13, -2.2D-11, -2.7D-01, r^2= 6.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.986271 1 O py 58 -1.727290 2 H py
- 81 -1.727290 3 H py 13 1.509212 1 O py
- 29 -0.399224 1 O d -1 55 -0.278141 2 H py
- 78 -0.278141 3 H py 65 0.247567 2 H d -2
- 88 -0.247567 3 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.159430D-01
- MO Center= -6.5D-13, 2.0D-11, 2.3D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.572020 1 O s 49 -3.849437 2 H s
- 72 -3.849437 3 H s 17 -2.512368 1 O pz
- 14 -2.019315 1 O pz 4 1.482861 1 O s
- 48 -1.100027 2 H s 71 -1.100027 3 H s
- 59 0.640329 2 H pz 82 0.640329 3 H pz
- Vector 17 Occ=0.000000D+00 E= 3.085504D-01
- MO Center= -6.7D-12, -2.1D-14, -3.0D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 57 4.126608 2 H px 80 4.126608 3 H px
- 15 -3.422136 1 O px 49 -3.385593 2 H s
- 72 3.385593 3 H s 59 1.768722 2 H pz
- 82 -1.768722 3 H pz 12 -1.195530 1 O px
- 50 0.836354 2 H s 73 -0.836354 3 H s
- Vector 18 Occ=0.000000D+00 E= 3.141721D-01
- MO Center= -7.5D-12, 4.7D-14, -8.2D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.824770 1 O s 49 -4.955274 2 H s
- 72 -4.955274 3 H s 4 4.157586 1 O s
- 57 2.084038 2 H px 80 -2.084038 3 H px
- 59 -1.803963 2 H pz 82 -1.803963 3 H pz
- 50 -0.754804 2 H s 73 -0.754804 3 H s
- Vector 19 Occ=0.000000D+00 E= 3.264317D-01
- MO Center= 1.4D-11, 3.9D-14, -9.1D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 7.979626 2 H s 72 -7.979626 3 H s
- 57 -3.017551 2 H px 80 -3.017551 3 H px
- 59 2.190324 2 H pz 82 -2.190324 3 H pz
- 50 2.038065 2 H s 73 -2.038065 3 H s
- 15 -1.992020 1 O px 12 -1.675598 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.550070D-01
- MO Center= -1.7D-13, 2.4D-13, 2.2D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 9.787791 2 H s 72 -9.787791 3 H s
- 12 -3.657142 1 O px 15 -3.016793 1 O px
- 57 -2.481838 2 H px 80 -2.481838 3 H px
- 31 -1.976175 1 O d 1 48 1.961322 2 H s
- 71 -1.961322 3 H s 50 1.505700 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.962259D-01
- MO Center= -1.6D-12, 4.0D-14, -6.4D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 6.837980 1 O s 49 -4.358384 2 H s
- 72 -4.358384 3 H s 4 1.570924 1 O s
- 32 1.212725 1 O d 2 59 -1.104487 2 H pz
- 82 -1.104487 3 H pz 17 -0.916073 1 O pz
- 30 -0.902317 1 O d 0 14 0.872050 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.448294D-01
- MO Center= -2.4D-13, -2.1D-13, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.125830 1 O s 49 -1.979604 2 H s
- 72 -1.979604 3 H s 14 -1.963457 1 O pz
- 17 -1.892946 1 O pz 30 1.613648 1 O d 0
- 59 1.018440 2 H pz 82 1.018440 3 H pz
- 4 0.967071 1 O s 57 0.603691 2 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.09921406
- moments of inertia (a.u.)
- ------------------
- 2.270969023714 0.000000000000 0.000000000000
- 0.000000000000 6.421267634361 0.000000000000
- 0.000000000000 0.000000000000 4.150298610647
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 0.000000 0.000000 0.000000 0.000000
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.712863 -0.356432 -0.356432 0.000000
- 2 2 0 0 -3.386308 -3.752191 -3.752191 4.118075
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 0.000000 0.000000 0.000000 0.000000
- 2 0 2 0 -5.836453 -2.918226 -2.918226 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.836459 -3.433149 -3.433149 2.029839
- NWChem DFT Gradient Module
- --------------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- charge = 0.00
- wavefunction = closed shell
- DFT ENERGY GRADIENTS
- atom coordinates gradient
- x y z x y z
- 1 O 0.000000 0.000000 0.225269 0.000000 0.000000 -0.005672
- 2 H 1.434935 0.000000 -0.901075 -0.005254 0.000000 0.002836
- 3 H -1.434935 0.000000 -0.901075 0.005254 0.000000 0.002836
- ----------------------------------------
- | Time | 1-e(secs) | 2-e(secs) |
- ----------------------------------------
- | CPU | 0.00 | 1.53 |
- ----------------------------------------
- | WALL | 0.01 | 2.93 |
- ----------------------------------------
- @ Step Energy Delta E Gmax Grms Xrms Xmax Walltime
- @ ---- ---------------- -------- -------- -------- -------- -------- --------
- @ 0 -76.46602603 0.0D+00 0.00588 0.00484 0.00000 0.00000 12.5
- Z-matrix (autoz)
- --------
- Units are Angstrom for bonds and degrees for angles
- Type Name I J K L M Value Gradient
- ----------- -------- ----- ----- ----- ----- ----- ---------- ----------
- 1 Stretch 1 2 0.96532 -0.00588
- 2 Stretch 1 3 0.96532 -0.00588
- 3 Bend 2 1 3 103.74015 -0.00098
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 5.00D-04
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: medium
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 49 10.0 434
- H 0.35 45 12.0 434
- Grid pruning is: on
- Number of quadrature shells: 139
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Loading old vectors from job with title :
- TDDFT h2o BP86/aug-cc-pVTZ
- Time after variat. SCF: 8.5
- Time prior to 1st pass: 8.5
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 7 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 13.02 13017841
- Stack Space remaining (MW): 13.11 13106626
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4660413621 -8.55D+01 2.13D-04 5.52D-04 9.0
- d= 0,ls=0.0,diis 2 -76.4660632123 -2.19D-05 8.76D-05 2.38D-04 9.7
- d= 0,ls=0.0,diis 3 -76.4660648700 -1.66D-06 4.59D-05 2.22D-04 10.5
- d= 0,ls=0.0,diis 4 -76.4660831635 -1.83D-05 3.16D-06 7.75D-07 11.3
- d= 0,ls=0.0,diis 5 -76.4660832337 -7.02D-08 3.49D-07 1.06D-08 12.1
- d= 0,ls=0.0,diis 6 -76.4660832347 -1.03D-09 1.38D-08 9.63D-12 12.6
- Total DFT energy = -76.466083234744
- One electron energy = -122.715623573819
- Coulomb energy = 46.522775086465
- Exchange-Corr. energy = -9.309455874776
- Nuclear repulsion energy = 9.036221127386
- Numeric. integr. density = 10.000002457006
- Total iterative time = 4.2s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878214D+01
- MO Center= -4.3D-10, -1.4D-18, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.973964 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.299922D-01
- MO Center= -8.2D-10, 4.9D-16, -8.3D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.482989 1 O s 4 0.282226 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.831361D-01
- MO Center= 3.4D-10, -3.9D-17, -9.6D-02, r^2= 7.8D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.325417 1 O px 6 0.300717 1 O px
- 48 0.225249 2 H s 71 -0.225249 3 H s
- Vector 4 Occ=2.000000D+00 E=-3.463817D-01
- MO Center= -3.1D-10, -6.4D-16, 1.8D-01, r^2= 7.6D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.359603 1 O pz 4 0.344100 1 O s
- 8 0.320055 1 O pz 14 0.206803 1 O pz
- 2 0.199619 1 O s 48 -0.175150 2 H s
- 71 -0.175150 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.693241D-01
- MO Center= -4.4D-10, 2.2D-16, 7.3D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428288 1 O py 7 0.365089 1 O py
- 13 0.324788 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.268044D-02
- MO Center= 4.4D-09, 3.0D-15, -7.7D-01, r^2= 7.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.474686 1 O s 50 -0.460282 2 H s
- 73 -0.460282 3 H s 4 0.364473 1 O s
- 49 -0.352734 2 H s 72 -0.352734 3 H s
- 14 -0.191694 1 O pz 2 0.151681 1 O s
- Vector 7 Occ=0.000000D+00 E= 1.495628D-02
- MO Center= -1.1D-08, 2.0D-17, -7.7D-01, r^2= 1.5D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.612716 2 H s 73 -2.612716 3 H s
- 49 0.795707 2 H s 72 -0.795707 3 H s
- 15 -0.641016 1 O px 12 -0.339606 1 O px
- 48 0.157658 2 H s 71 -0.157658 3 H s
- Vector 8 Occ=0.000000D+00 E= 6.574366D-02
- MO Center= 4.3D-08, -4.9D-15, -7.2D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.370085 1 O s 49 -2.110737 2 H s
- 72 -2.110737 3 H s 4 1.359686 1 O s
- 17 -1.122883 1 O pz 50 -1.023913 2 H s
- 73 -1.023913 3 H s 57 0.688523 2 H px
- 80 -0.688523 3 H px 59 -0.541246 2 H pz
- Vector 9 Occ=0.000000D+00 E= 9.422570D-02
- MO Center= 1.5D-10, -2.2D-15, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.491757 1 O py 58 -0.307444 2 H py
- 81 -0.307444 3 H py
- Vector 10 Occ=0.000000D+00 E= 9.868144D-02
- MO Center= -5.4D-08, 1.7D-16, 1.6D-01, r^2= 8.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.467255 1 O pz 49 0.594714 2 H s
- 72 0.594714 3 H s 50 -0.354001 2 H s
- 73 -0.354001 3 H s 59 -0.301349 2 H pz
- 82 -0.301349 3 H pz 57 0.275734 2 H px
- 80 -0.275734 3 H px 4 0.251059 1 O s
- Vector 11 Occ=0.000000D+00 E= 1.002086D-01
- MO Center= 1.7D-08, -1.8D-15, -8.2D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.887615 2 H s 73 -2.887615 3 H s
- 49 2.085433 2 H s 72 -2.085433 3 H s
- 15 -1.546318 1 O px 57 -0.917097 2 H px
- 80 -0.917097 3 H px 59 0.903016 2 H pz
- 82 -0.903016 3 H pz 12 -0.538022 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.477023D-01
- MO Center= 1.8D-08, -5.9D-16, 6.4D-01, r^2= 8.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.211003 1 O px 49 -2.844879 2 H s
- 72 2.844879 3 H s 50 -2.154571 2 H s
- 73 2.154571 3 H s 12 1.384824 1 O px
- 48 -0.935707 2 H s 71 0.935707 3 H s
- 31 0.577224 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.677130D-01
- MO Center= 2.0D-08, -2.5D-15, -1.8D-01, r^2= 6.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.468321 1 O s 49 -1.655475 2 H s
- 72 -1.655475 3 H s 59 -0.988225 2 H pz
- 82 -0.988225 3 H pz 57 -0.601719 2 H px
- 80 0.601719 3 H px 17 0.454876 1 O pz
- 4 0.444475 1 O s 30 -0.422397 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.908806D-01
- MO Center= 2.7D-09, -1.1D-17, -5.0D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 58 1.575598 2 H py 81 -1.575598 3 H py
- 28 -0.470898 1 O d -2 55 0.228424 2 H py
- 78 -0.228424 3 H py 65 -0.184356 2 H d -2
- 88 -0.184356 3 H d -2 66 0.166365 2 H d -1
- 89 -0.166365 3 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.138537D-01
- MO Center= -3.3D-09, -6.6D-14, -2.7D-01, r^2= 6.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.983094 1 O py 58 -1.719757 2 H py
- 81 -1.719757 3 H py 13 1.479191 1 O py
- 29 -0.393702 1 O d -1 55 -0.268968 2 H py
- 78 -0.268968 3 H py 65 0.246688 2 H d -2
- 88 -0.246688 3 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.156112D-01
- MO Center= -2.6D-08, 7.1D-14, 2.4D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.583261 1 O s 49 -3.855146 2 H s
- 72 -3.855145 3 H s 17 -2.513751 1 O pz
- 14 -1.994179 1 O pz 4 1.469430 1 O s
- 48 -1.068222 2 H s 71 -1.068221 3 H s
- 59 0.631100 2 H pz 82 0.631100 3 H pz
- Vector 17 Occ=0.000000D+00 E= 3.077501D-01
- MO Center= -3.4D-07, 1.4D-15, -3.1D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 57 3.987596 2 H px 80 3.987600 3 H px
- 15 -3.419630 1 O px 49 -3.060455 2 H s
- 72 3.060463 3 H s 59 1.760175 2 H pz
- 82 -1.760172 3 H pz 12 -1.198502 1 O px
- 50 0.835142 2 H s 73 -0.835141 3 H s
- Vector 18 Occ=0.000000D+00 E= 3.121412D-01
- MO Center= 4.7D-07, -1.8D-15, -8.2D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.670227 1 O s 49 -4.849687 2 H s
- 72 -4.849682 3 H s 4 4.128243 1 O s
- 57 2.087233 2 H px 80 -2.087226 3 H px
- 59 -1.785236 2 H pz 82 -1.785239 3 H pz
- 50 -0.750973 2 H s 73 -0.750975 3 H s
- Vector 19 Occ=0.000000D+00 E= 3.240693D-01
- MO Center= -1.5D-07, -1.5D-17, -9.0D-01, r^2= 7.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 7.737068 2 H s 72 -7.737068 3 H s
- 57 -2.995707 2 H px 80 -2.995707 3 H px
- 59 2.143751 2 H pz 82 -2.143751 3 H pz
- 50 2.005772 2 H s 73 -2.005772 3 H s
- 15 -1.943175 1 O px 12 -1.649512 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.553791D-01
- MO Center= 4.5D-08, -3.9D-15, 2.2D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 9.618988 2 H s 72 -9.618989 3 H s
- 12 -3.613936 1 O px 15 -3.037920 1 O px
- 57 -2.452982 2 H px 80 -2.452982 3 H px
- 31 -1.953276 1 O d 1 48 1.901980 2 H s
- 71 -1.901980 3 H s 50 1.497828 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.944955D-01
- MO Center= -2.2D-08, -2.7D-15, -6.4D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 6.995230 1 O s 49 -4.460651 2 H s
- 72 -4.460650 3 H s 4 1.574263 1 O s
- 32 1.186815 1 O d 2 59 -1.128876 2 H pz
- 82 -1.128876 3 H pz 17 -0.936584 1 O pz
- 30 -0.897166 1 O d 0 14 0.865872 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.451787D-01
- MO Center= -1.5D-08, 2.4D-15, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.139134 1 O s 49 -1.992376 2 H s
- 72 -1.992376 3 H s 14 -1.951359 1 O pz
- 17 -1.890137 1 O pz 30 1.607841 1 O d 0
- 59 1.008637 2 H pz 82 1.008637 3 H pz
- 4 0.997262 1 O s 57 0.614244 2 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10059264
- moments of inertia (a.u.)
- ------------------
- 2.281000936723 0.000000000000 0.000000003368
- 0.000000000000 6.541958071871 0.000000000000
- 0.000000003368 0.000000000000 4.260957135149
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 0.000000 0.000000 0.000000 0.000000
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.710460 -0.361028 -0.361028 0.011597
- 2 2 0 0 -3.355875 -3.791874 -3.791874 4.227874
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 0.000000 0.000000 0.000000 0.000000
- 2 0 2 0 -5.852722 -2.926361 -2.926361 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.856445 -3.447632 -3.447632 2.038819
- Line search:
- step= 1.00 grad=-2.1D-04 hess= 1.6D-04 energy= -76.466083 mode=downhill
- new step= 0.68 predicted energy= -76.466099
- --------
- Step 1
- --------
- Geometry "geometry" -> "geometry"
- ---------------------------------
- Output coordinates in angstroms (scale by 1.889725989 to convert to a.u.)
- No. Tag Charge X Y Z
- ---- ---------------- ---------- -------------- -------------- --------------
- 1 O 8.0000 0.00000000 0.00000000 0.11981191
- 2 H 1.0000 0.76619470 0.00000000 -0.47713060
- 3 H 1.0000 -0.76619470 0.00000000 -0.47713060
- Atomic Mass
- -----------
- O 15.994910
- H 1.007825
- Effective nuclear repulsion energy (a.u.) 9.0624764652
- Nuclear Dipole moment (a.u.)
- ----------------------------
- X Y Z
- ---------------- ---------------- ----------------
- -0.0000000039 0.0000000000 0.0080012488
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 5.00D-04
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: medium
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 49 10.0 434
- H 0.35 45 12.0 434
- Grid pruning is: on
- Number of quadrature shells: 139
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Loading old vectors from job with title :
- TDDFT h2o BP86/aug-cc-pVTZ
- Time after variat. SCF: 12.7
- Time prior to 1st pass: 12.7
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 7 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 13.02 13017841
- Stack Space remaining (MW): 13.11 13106626
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4660947356 -8.55D+01 6.70D-05 5.39D-05 13.6
- d= 0,ls=0.0,diis 2 -76.4660966580 -1.92D-06 2.79D-05 2.54D-05 14.5
- d= 0,ls=0.0,diis 3 -76.4660970350 -3.77D-07 1.43D-05 2.14D-05 15.1
- d= 0,ls=0.0,diis 4 -76.4660988032 -1.77D-06 9.98D-07 7.87D-08 16.0
- d= 0,ls=0.0,diis 5 -76.4660987885 1.47D-08 1.10D-07 1.04D-09 16.9
- d= 0,ls=0.0,diis 6 -76.4660987886 -1.01D-10 3.68D-09 6.60D-13 17.5
- Total DFT energy = -76.466098788589
- One electron energy = -122.761897427067
- Coulomb energy = 46.545401800047
- Exchange-Corr. energy = -9.312079626795
- Nuclear repulsion energy = 9.062476465226
- Numeric. integr. density = 10.000002306187
- Total iterative time = 4.8s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878177D+01
- MO Center= -2.9D-10, -3.4D-18, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.974031 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.312436D-01
- MO Center= -5.6D-10, -1.6D-15, -8.4D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.482527 1 O s 4 0.280698 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.837315D-01
- MO Center= 2.3D-10, -1.1D-15, -9.6D-02, r^2= 7.8D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.325136 1 O px 6 0.301017 1 O px
- 48 0.226241 2 H s 71 -0.226241 3 H s
- Vector 4 Occ=2.000000D+00 E=-3.468925D-01
- MO Center= -2.1D-10, -5.7D-15, 1.8D-01, r^2= 7.5D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.359422 1 O pz 4 0.345274 1 O s
- 8 0.320061 1 O pz 14 0.206029 1 O pz
- 2 0.200251 1 O s 48 -0.175442 2 H s
- 71 -0.175442 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.695710D-01
- MO Center= -3.0D-10, 6.1D-14, 7.3D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428227 1 O py 7 0.365057 1 O py
- 13 0.324375 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.237131D-02
- MO Center= 3.0D-09, -1.1D-12, -7.7D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.479938 1 O s 50 -0.462720 2 H s
- 73 -0.462720 3 H s 4 0.365679 1 O s
- 49 -0.353263 2 H s 72 -0.353263 3 H s
- 14 -0.191584 1 O pz 2 0.150968 1 O s
- Vector 7 Occ=0.000000D+00 E= 1.513718D-02
- MO Center= -7.8D-09, -7.9D-16, -7.7D-01, r^2= 1.5D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.630530 2 H s 73 -2.630530 3 H s
- 49 0.801558 2 H s 72 -0.801558 3 H s
- 15 -0.643514 1 O px 12 -0.340624 1 O px
- 48 0.158241 2 H s 71 -0.158241 3 H s
- Vector 8 Occ=0.000000D+00 E= 6.598967D-02
- MO Center= 2.9D-08, -1.0D-10, -6.6D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.389022 1 O s 49 -2.123919 2 H s
- 72 -2.123919 3 H s 4 1.366631 1 O s
- 17 -1.132851 1 O pz 50 -1.022110 2 H s
- 73 -1.022110 3 H s 57 0.685660 2 H px
- 80 -0.685660 3 H px 59 -0.542355 2 H pz
- Vector 9 Occ=0.000000D+00 E= 9.419808D-02
- MO Center= 1.0D-10, 1.7D-10, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.494050 1 O py 58 -0.309390 2 H py
- 81 -0.309390 3 H py
- Vector 10 Occ=0.000000D+00 E= 9.887753D-02
- MO Center= -3.0D-08, -4.7D-11, 1.6D-01, r^2= 8.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.460070 1 O pz 49 0.567132 2 H s
- 72 0.567132 3 H s 50 -0.361964 2 H s
- 73 -0.361964 3 H s 59 -0.310776 2 H pz
- 82 -0.310776 3 H pz 57 0.283672 2 H px
- 80 -0.283672 3 H px 5 0.279488 1 O s
- Vector 11 Occ=0.000000D+00 E= 1.008517D-01
- MO Center= 5.3D-09, -1.6D-14, -8.2D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.903892 2 H s 73 -2.903892 3 H s
- 49 2.123153 2 H s 72 -2.123153 3 H s
- 15 -1.562342 1 O px 57 -0.924290 2 H px
- 80 -0.924290 3 H px 59 0.910770 2 H pz
- 82 -0.910770 3 H pz 12 -0.546506 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.477820D-01
- MO Center= 1.2D-08, 6.2D-14, 6.4D-01, r^2= 8.7D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.206850 1 O px 49 -2.833647 2 H s
- 72 2.833647 3 H s 50 -2.151834 2 H s
- 73 2.151834 3 H s 12 1.389336 1 O px
- 48 -0.940113 2 H s 71 0.940113 3 H s
- 31 0.579342 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.680154D-01
- MO Center= 1.4D-08, -9.4D-13, -1.8D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.437251 1 O s 49 -1.646560 2 H s
- 72 -1.646561 3 H s 59 -0.987649 2 H pz
- 82 -0.987649 3 H pz 57 -0.612086 2 H px
- 80 0.612086 3 H px 17 0.458450 1 O pz
- 4 0.442183 1 O s 30 -0.424103 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.910792D-01
- MO Center= 1.9D-09, 6.9D-15, -5.0D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 58 1.582323 2 H py 81 -1.582323 3 H py
- 28 -0.471786 1 O d -2 55 0.230512 2 H py
- 78 -0.230512 3 H py 65 -0.185180 2 H d -2
- 88 -0.185180 3 H d -2 66 0.167062 2 H d -1
- 89 -0.167062 3 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.138497D-01
- MO Center= -2.3D-09, -8.1D-10, -2.7D-01, r^2= 6.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.984100 1 O py 58 -1.722151 2 H py
- 81 -1.722151 3 H py 13 1.488733 1 O py
- 29 -0.395478 1 O d -1 55 -0.271877 2 H py
- 78 -0.271877 3 H py 65 0.246982 2 H d -2
- 88 -0.246982 3 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.157152D-01
- MO Center= -1.8D-08, 7.9D-10, 2.3D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.579816 1 O s 49 -3.853388 2 H s
- 72 -3.853387 3 H s 17 -2.513338 1 O pz
- 14 -2.002142 1 O pz 4 1.473602 1 O s
- 48 -1.078257 2 H s 71 -1.078256 3 H s
- 59 0.634004 2 H pz 82 0.634004 3 H pz
- Vector 17 Occ=0.000000D+00 E= 3.080028D-01
- MO Center= -2.2D-07, -3.2D-14, -3.1D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 57 4.032125 2 H px 80 4.032127 3 H px
- 15 -3.420015 1 O px 49 -3.163629 2 H s
- 72 3.163634 3 H s 59 1.762388 2 H pz
- 82 -1.762386 3 H pz 12 -1.197607 1 O px
- 50 0.834956 2 H s 73 -0.834956 3 H s
- Vector 18 Occ=0.000000D+00 E= 3.127859D-01
- MO Center= 3.1D-07, -6.6D-13, -8.2D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.719573 1 O s 49 -4.883405 2 H s
- 72 -4.883401 3 H s 4 4.137977 1 O s
- 57 2.086399 2 H px 80 -2.086395 3 H px
- 59 -1.791239 2 H pz 82 -1.791240 3 H pz
- 50 -0.752206 2 H s 73 -0.752207 3 H s
- Vector 19 Occ=0.000000D+00 E= 3.248154D-01
- MO Center= -9.9D-08, -8.9D-15, -9.1D-01, r^2= 7.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 7.812150 2 H s 72 -7.812150 3 H s
- 57 -3.001442 2 H px 80 -3.001442 3 H px
- 59 2.158936 2 H pz 82 -2.158937 3 H pz
- 50 2.016100 2 H s 73 -2.016100 3 H s
- 15 -1.959468 1 O px 12 -1.657962 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.552601D-01
- MO Center= 3.1D-08, -3.2D-14, 2.2D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 9.671971 2 H s 72 -9.671971 3 H s
- 12 -3.627642 1 O px 15 -3.031274 1 O px
- 57 -2.462086 2 H px 80 -2.462086 3 H px
- 31 -1.960554 1 O d 1 48 1.920849 2 H s
- 71 -1.920849 3 H s 50 1.500396 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.950389D-01
- MO Center= -1.6D-08, -1.1D-11, -6.4D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 6.945443 1 O s 49 -4.428278 2 H s
- 72 -4.428278 3 H s 4 1.572622 1 O s
- 32 1.195002 1 O d 2 59 -1.121261 2 H pz
- 82 -1.121261 3 H pz 17 -0.930047 1 O pz
- 30 -0.898956 1 O d 0 14 0.868374 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.450654D-01
- MO Center= -1.1D-08, -1.4D-11, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.135590 1 O s 49 -1.988778 2 H s
- 72 -1.988778 3 H s 14 -1.955563 1 O pz
- 17 -1.891164 1 O pz 30 1.609826 1 O d 0
- 59 1.011741 2 H pz 82 1.011741 3 H pz
- 4 0.988117 1 O s 57 0.611012 2 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10016521
- moments of inertia (a.u.)
- ------------------
- 2.277888159786 0.000000000000 0.000000002297
- 0.000000000000 6.503514452891 0.000000000000
- 0.000000002297 0.000000000000 4.225626293105
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 0.000000 0.000000 0.000000 0.000000
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.711234 -0.359618 -0.359618 0.008001
- 2 2 0 0 -3.365608 -3.779213 -3.779213 4.192817
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 0.000000 0.000000 0.000000 0.000000
- 2 0 2 0 -5.847566 -2.923783 -2.923783 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.850107 -3.443068 -3.443068 2.036029
- NWChem DFT Gradient Module
- --------------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- charge = 0.00
- wavefunction = closed shell
- DFT ENERGY GRADIENTS
- atom coordinates gradient
- x y z x y z
- 1 O 0.000000 0.000000 0.226412 0.000000 0.000000 0.000060
- 2 H 1.447898 0.000000 -0.901646 0.000025 0.000000 -0.000030
- 3 H -1.447898 0.000000 -0.901646 -0.000025 0.000000 -0.000030
- ----------------------------------------
- | Time | 1-e(secs) | 2-e(secs) |
- ----------------------------------------
- | CPU | 0.00 | 2.60 |
- ----------------------------------------
- | WALL | 0.00 | 2.61 |
- ----------------------------------------
- Step Energy Delta E Gmax Grms Xrms Xmax Walltime
- ---- ---------------- -------- -------- -------- -------- -------- --------
- @ 1 -76.46609879 -7.3D-05 0.00004 0.00003 0.00614 0.01299 28.4
- ok ok
- Z-matrix (autoz)
- --------
- Units are Angstrom for bonds and degrees for angles
- Type Name I J K L M Value Gradient
- ----------- -------- ----- ----- ----- ----- ----- ---------- ----------
- 1 Stretch 1 2 0.97128 0.00004
- 2 Stretch 1 3 0.97128 0.00004
- 3 Bend 2 1 3 104.15567 -0.00001
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 5.00D-04
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: medium
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 49 10.0 434
- H 0.35 45 12.0 434
- Grid pruning is: on
- Number of quadrature shells: 139
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Loading old vectors from job with title :
- TDDFT h2o BP86/aug-cc-pVTZ
- Time after variat. SCF: 20.8
- Time prior to 1st pass: 20.8
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 7 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 13.02 13017841
- Stack Space remaining (MW): 13.11 13106626
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4660987911 -8.55D+01 1.32D-06 1.75D-08 21.3
- d= 0,ls=0.0,diis 2 -76.4660987926 -1.51D-09 3.18D-07 2.73D-09 22.2
- Total DFT energy = -76.466098792590
- One electron energy = -122.762684266314
- Coulomb energy = 46.545862158007
- Exchange-Corr. energy = -9.312132265746
- Nuclear repulsion energy = 9.062855581463
- Numeric. integr. density = 10.000002304603
- Total iterative time = 1.4s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878175D+01
- MO Center= -1.0D-10, 5.2D-18, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.974032 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.312504D-01
- MO Center= -2.2D-10, 2.2D-15, -8.4D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.482522 1 O s 4 0.280680 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.837516D-01
- MO Center= 9.6D-11, 2.9D-16, -9.6D-02, r^2= 7.8D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.325124 1 O px 6 0.301016 1 O px
- 48 0.226262 2 H s 71 -0.226262 3 H s
- Vector 4 Occ=2.000000D+00 E=-3.468786D-01
- MO Center= -8.7D-11, 1.1D-15, 1.8D-01, r^2= 7.5D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.359427 1 O pz 4 0.345280 1 O s
- 8 0.320066 1 O pz 14 0.206034 1 O pz
- 2 0.200245 1 O s 48 -0.175440 2 H s
- 71 -0.175440 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.695671D-01
- MO Center= -1.1D-10, -3.3D-15, 7.3D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428223 1 O py 7 0.365055 1 O py
- 13 0.324371 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.236468D-02
- MO Center= 7.0D-10, -3.6D-14, -7.7D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.480099 1 O s 50 -0.462773 2 H s
- 73 -0.462773 3 H s 4 0.365725 1 O s
- 49 -0.353307 2 H s 72 -0.353307 3 H s
- 14 -0.191567 1 O pz 2 0.150961 1 O s
- Vector 7 Occ=0.000000D+00 E= 1.513959D-02
- MO Center= -2.8D-09, 5.7D-15, -7.7D-01, r^2= 1.5D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.630672 2 H s 73 -2.630672 3 H s
- 49 0.801623 2 H s 72 -0.801623 3 H s
- 15 -0.643579 1 O px 12 -0.340656 1 O px
- 48 0.158251 2 H s 71 -0.158251 3 H s
- Vector 8 Occ=0.000000D+00 E= 6.599518D-02
- MO Center= 1.0D-08, -7.3D-13, -6.6D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.389489 1 O s 49 -2.124195 2 H s
- 72 -2.124195 3 H s 4 1.366806 1 O s
- 17 -1.132979 1 O pz 50 -1.022097 2 H s
- 73 -1.022097 3 H s 57 0.685721 2 H px
- 80 -0.685721 3 H px 59 -0.542294 2 H pz
- Vector 9 Occ=0.000000D+00 E= 9.419781D-02
- MO Center= 2.3D-12, 7.5D-13, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.494151 1 O py 58 -0.309476 2 H py
- 81 -0.309476 3 H py
- Vector 10 Occ=0.000000D+00 E= 9.887826D-02
- MO Center= -1.1D-08, 6.6D-13, 1.6D-01, r^2= 8.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.460017 1 O pz 49 0.566798 2 H s
- 72 0.566798 3 H s 50 -0.362042 2 H s
- 73 -0.362042 3 H s 59 -0.310915 2 H pz
- 82 -0.310915 3 H pz 57 0.283761 2 H px
- 80 -0.283761 3 H px 5 0.280149 1 O s
- Vector 11 Occ=0.000000D+00 E= 1.008637D-01
- MO Center= 2.2D-09, -4.6D-16, -8.2D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.904173 2 H s 73 -2.904173 3 H s
- 49 2.123801 2 H s 72 -2.123801 3 H s
- 15 -1.562911 1 O px 57 -0.924350 2 H px
- 80 -0.924350 3 H px 59 0.910874 2 H pz
- 82 -0.910874 3 H pz 12 -0.546783 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.477860D-01
- MO Center= 5.1D-09, 1.7D-14, 6.4D-01, r^2= 8.7D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.206731 1 O px 49 -2.833380 2 H s
- 72 2.833380 3 H s 50 -2.151478 2 H s
- 73 2.151478 3 H s 12 1.389496 1 O px
- 48 -0.940257 2 H s 71 0.940257 3 H s
- 31 0.579424 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.680165D-01
- MO Center= 4.8D-09, -3.5D-14, -1.8D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.438152 1 O s 49 -1.647014 2 H s
- 72 -1.647014 3 H s 59 -0.987626 2 H pz
- 82 -0.987626 3 H pz 57 -0.612032 2 H px
- 80 0.612032 3 H px 17 0.458287 1 O pz
- 4 0.442469 1 O s 30 -0.424069 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.910792D-01
- MO Center= 6.3D-10, 4.2D-15, -5.0D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 58 1.582360 2 H py 81 -1.582360 3 H py
- 28 -0.471843 1 O d -2 55 0.230555 2 H py
- 78 -0.230555 3 H py 65 -0.185192 2 H d -2
- 88 -0.185192 3 H d -2 66 0.167067 2 H d -1
- 89 -0.167067 3 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.138588D-01
- MO Center= -7.5D-10, -2.8D-12, -2.7D-01, r^2= 6.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.984206 1 O py 58 -1.722279 2 H py
- 81 -1.722279 3 H py 13 1.488941 1 O py
- 29 -0.395484 1 O d -1 55 -0.271924 2 H py
- 78 -0.271924 3 H py 65 0.247016 2 H d -2
- 88 -0.247016 3 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.157235D-01
- MO Center= -6.7D-09, 2.5D-12, 2.3D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.579372 1 O s 49 -3.853227 2 H s
- 72 -3.853227 3 H s 17 -2.513391 1 O pz
- 14 -2.002301 1 O pz 4 1.473712 1 O s
- 48 -1.078325 2 H s 71 -1.078325 3 H s
- 59 0.634325 2 H pz 82 0.634326 3 H pz
- Vector 17 Occ=0.000000D+00 E= 3.080045D-01
- MO Center= -7.9D-08, 4.2D-15, -3.1D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 57 4.031830 2 H px 80 4.031831 3 H px
- 15 -3.420498 1 O px 49 -3.162555 2 H s
- 72 3.162556 3 H s 59 1.762655 2 H pz
- 82 -1.762654 3 H pz 12 -1.197847 1 O px
- 50 0.835350 2 H s 73 -0.835350 3 H s
- Vector 18 Occ=0.000000D+00 E= 3.127968D-01
- MO Center= 1.1D-07, 2.5D-15, -8.2D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.720222 1 O s 49 -4.883872 2 H s
- 72 -4.883871 3 H s 4 4.138308 1 O s
- 57 2.086532 2 H px 80 -2.086531 3 H px
- 59 -1.791161 2 H pz 82 -1.791161 3 H pz
- 50 -0.752216 2 H s 73 -0.752217 3 H s
- Vector 19 Occ=0.000000D+00 E= 3.248310D-01
- MO Center= -3.6D-08, -1.8D-15, -9.1D-01, r^2= 7.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 7.813608 2 H s 72 -7.813608 3 H s
- 57 -3.002205 2 H px 80 -3.002205 3 H px
- 59 2.158704 2 H pz 82 -2.158704 3 H pz
- 50 2.016088 2 H s 73 -2.016088 3 H s
- 15 -1.959345 1 O px 12 -1.657984 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.552613D-01
- MO Center= 1.1D-08, 5.2D-15, 2.2D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 9.671425 2 H s 72 -9.671425 3 H s
- 12 -3.628209 1 O px 15 -3.031297 1 O px
- 57 -2.461718 2 H px 80 -2.461718 3 H px
- 31 -1.960879 1 O d 1 48 1.921414 2 H s
- 71 -1.921414 3 H s 50 1.500308 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.950364D-01
- MO Center= -5.6D-09, -7.2D-14, -6.4D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 6.946074 1 O s 49 -4.428599 2 H s
- 72 -4.428599 3 H s 4 1.573213 1 O s
- 32 1.195059 1 O d 2 59 -1.121207 2 H pz
- 82 -1.121207 3 H pz 17 -0.930028 1 O pz
- 30 -0.898837 1 O d 0 14 0.868119 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.450725D-01
- MO Center= -4.2D-09, 4.0D-14, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.134191 1 O s 49 -1.987888 2 H s
- 72 -1.987889 3 H s 14 -1.955349 1 O pz
- 17 -1.890991 1 O pz 30 1.609738 1 O d 0
- 59 1.011979 2 H pz 82 1.011979 3 H pz
- 4 0.987542 1 O s 57 0.610797 2 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10009567
- moments of inertia (a.u.)
- ------------------
- 2.277381949965 0.000000000000 0.000000000830
- 0.000000000000 6.502988861392 0.000000000000
- 0.000000000830 0.000000000000 4.225606911426
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 0.000000 0.000000 0.000000 0.000000
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.711178 -0.359297 -0.359297 0.007416
- 2 2 0 0 -3.365545 -3.779172 -3.779172 4.192798
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 0.000000 0.000000 0.000000 0.000000
- 2 0 2 0 -5.847559 -2.923779 -2.923779 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.850198 -3.442887 -3.442887 2.035576
- Line search:
- step= 1.00 grad=-7.7D-09 hess= 3.7D-09 energy= -76.466099 mode=accept
- new step= 1.00 predicted energy= -76.466099
- --------
- Step 2
- --------
- Geometry "geometry" -> "geometry"
- ---------------------------------
- Output coordinates in angstroms (scale by 1.889725989 to convert to a.u.)
- No. Tag Charge X Y Z
- ---- ---------------- ---------- -------------- -------------- --------------
- 1 O 8.0000 0.00000000 0.00000000 0.11976769
- 2 H 1.0000 0.76619295 0.00000000 -0.47710849
- 3 H 1.0000 -0.76619295 0.00000000 -0.47710849
- Atomic Mass
- -----------
- O 15.994910
- H 1.007825
- Effective nuclear repulsion energy (a.u.) 9.0628555815
- Nuclear Dipole moment (a.u.)
- ----------------------------
- X Y Z
- ---------------- ---------------- ----------------
- -0.0000000014 0.0000000000 0.0074162831
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- The DFT is already converged
- Total DFT energy = -76.466098792590
- NWChem DFT Gradient Module
- --------------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- charge = 0.00
- wavefunction = closed shell
- DFT ENERGY GRADIENTS
- atom coordinates gradient
- x y z x y z
- 1 O 0.000000 0.000000 0.226328 0.000000 0.000000 0.000006
- 2 H 1.447895 0.000000 -0.901604 0.000002 0.000000 -0.000003
- 3 H -1.447895 0.000000 -0.901604 -0.000002 0.000000 -0.000003
- ----------------------------------------
- | Time | 1-e(secs) | 2-e(secs) |
- ----------------------------------------
- | CPU | 0.00 | 2.96 |
- ----------------------------------------
- | WALL | 0.00 | 2.98 |
- ----------------------------------------
- Step Energy Delta E Gmax Grms Xrms Xmax Walltime
- ---- ---------------- -------- -------- -------- -------- -------- --------
- @ 2 -76.46609879 -4.0D-09 0.00000 0.00000 0.00003 0.00008 34.2
- ok ok ok ok
- Z-matrix (autoz)
- --------
- Units are Angstrom for bonds and degrees for angles
- Type Name I J K L M Value Gradient
- ----------- -------- ----- ----- ----- ----- ----- ---------- ----------
- 1 Stretch 1 2 0.97124 0.00000
- 2 Stretch 1 3 0.97124 0.00000
- 3 Bend 2 1 3 104.16171 0.00000
- ----------------------
- Optimization converged
- ----------------------
- Step Energy Delta E Gmax Grms Xrms Xmax Walltime
- ---- ---------------- -------- -------- -------- -------- -------- --------
- @ 2 -76.46609879 -4.0D-09 0.00000 0.00000 0.00003 0.00008 34.2
- ok ok ok ok
- Z-matrix (autoz)
- --------
- Units are Angstrom for bonds and degrees for angles
- Type Name I J K L M Value Gradient
- ----------- -------- ----- ----- ----- ----- ----- ---------- ----------
- 1 Stretch 1 2 0.97124 0.00000
- 2 Stretch 1 3 0.97124 0.00000
- 3 Bend 2 1 3 104.16171 0.00000
- Geometry "geometry" -> "geometry"
- ---------------------------------
- Output coordinates in angstroms (scale by 1.889725989 to convert to a.u.)
- No. Tag Charge X Y Z
- ---- ---------------- ---------- -------------- -------------- --------------
- 1 O 8.0000 0.00000000 0.00000000 0.11976769
- 2 H 1.0000 0.76619295 0.00000000 -0.47710849
- 3 H 1.0000 -0.76619295 0.00000000 -0.47710849
- Atomic Mass
- -----------
- O 15.994910
- H 1.007825
- Effective nuclear repulsion energy (a.u.) 9.0628555815
- Nuclear Dipole moment (a.u.)
- ----------------------------
- X Y Z
- ---------------- ---------------- ----------------
- -0.0000000014 0.0000000000 0.0074162831
- Final and change from initial internal coordinates
- --------------------------------------------------
- Z-matrix (autoz)
- --------
- Units are Angstrom for bonds and degrees for angles
- Type Name I J K L M Value Change
- ----------- -------- ----- ----- ----- ----- ----- ---------- ----------
- 1 Stretch 1 2 0.97124 0.00592
- 2 Stretch 1 3 0.97124 0.00592
- 3 Bend 2 1 3 104.16171 0.42157
- ==============================================================================
- internuclear distances
- ------------------------------------------------------------------------------
- center one | center two | atomic units | angstroms
- ------------------------------------------------------------------------------
- 2 H | 1 O | 1.83538 | 0.97124
- 3 H | 1 O | 1.83538 | 0.97124
- ------------------------------------------------------------------------------
- number of included internuclear distances: 2
- ==============================================================================
- ==============================================================================
- internuclear angles
- ------------------------------------------------------------------------------
- center 1 | center 2 | center 3 | degrees
- ------------------------------------------------------------------------------
- 2 H | 1 O | 3 H | 104.16
- ------------------------------------------------------------------------------
- number of included internuclear angles: 1
- ==============================================================================
- Task times cpu: 25.8s wall: 33.0s
- NWChem Input Module
- -------------------
- Raman Analysis
- --------------
- NWChem Nuclear Hessian and Frequency Analysis
- ---------------------------------------------
- NWChem Analytic Hessian
- -----------------------
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 1.00D-08
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: fine
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 70 10.0 590
- H 0.35 60 13.0 590
- Grid pruning is: on
- Number of quadrature shells: 190
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Loading old vectors from job with title :
- TDDFT h2o BP86/aug-cc-pVTZ
- Time after variat. SCF: 25.8
- Time prior to 1st pass: 25.8
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 12 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 12.96 12956396
- Stack Space remaining (MW): 13.11 13106626
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4660924934 -8.55D+01 5.93D-07 3.41D-09 26.6
- d= 0,ls=0.0,diis 2 -76.4660924926 8.27D-10 3.47D-07 1.33D-08 27.7
- d= 0,ls=0.0,diis 3 -76.4660924937 -1.10D-09 1.59D-08 1.21D-11 29.0
- Total DFT energy = -76.466092493711
- One electron energy = -122.762568931589
- Coulomb energy = 46.545731009595
- Exchange-Corr. energy = -9.312110153179
- Nuclear repulsion energy = 9.062855581463
- Numeric. integr. density = 10.000000000532
- Total iterative time = 3.1s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878177D+01
- MO Center= -1.0D-10, 2.2D-13, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.974032 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.312562D-01
- MO Center= -4.9D-11, 5.6D-11, -8.4D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.482524 1 O s 4 0.280683 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.837557D-01
- MO Center= -3.2D-11, -2.5D-11, -9.6D-02, r^2= 7.8D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.325127 1 O px 6 0.301018 1 O px
- 48 0.226259 2 H s 71 -0.226259 3 H s
- Vector 4 Occ=2.000000D+00 E=-3.468834D-01
- MO Center= -2.8D-10, -3.7D-10, 1.8D-01, r^2= 7.5D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.359429 1 O pz 4 0.345274 1 O s
- 8 0.320068 1 O pz 14 0.206033 1 O pz
- 2 0.200244 1 O s 48 -0.175438 2 H s
- 71 -0.175438 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.695736D-01
- MO Center= -2.0D-10, -8.1D-11, 7.3D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428225 1 O py 7 0.365057 1 O py
- 13 0.324368 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.237113D-02
- MO Center= 4.7D-09, 2.2D-09, -7.7D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.481394 1 O s 50 -0.462949 2 H s
- 73 -0.462949 3 H s 4 0.366227 1 O s
- 49 -0.353909 2 H s 72 -0.353909 3 H s
- 14 -0.191736 1 O pz 2 0.150945 1 O s
- Vector 7 Occ=0.000000D+00 E= 1.509566D-02
- MO Center= -7.0D-09, 2.4D-10, -7.7D-01, r^2= 1.5D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.633144 2 H s 73 -2.633144 3 H s
- 49 0.805015 2 H s 72 -0.805015 3 H s
- 15 -0.645808 1 O px 12 -0.342008 1 O px
- 48 0.158955 2 H s 71 -0.158955 3 H s
- Vector 8 Occ=0.000000D+00 E= 6.596699D-02
- MO Center= 1.3D-08, 2.7D-09, -6.7D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.391413 1 O s 49 -2.125060 2 H s
- 72 -2.125060 3 H s 4 1.367710 1 O s
- 17 -1.132810 1 O pz 50 -1.022267 2 H s
- 73 -1.022267 3 H s 57 0.686017 2 H px
- 80 -0.686017 3 H px 59 -0.542587 2 H pz
- Vector 9 Occ=0.000000D+00 E= 9.419700D-02
- MO Center= 1.8D-09, -5.7D-09, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.494142 1 O py 58 -0.309468 2 H py
- 81 -0.309468 3 H py
- Vector 10 Occ=0.000000D+00 E= 9.887440D-02
- MO Center= 2.8D-08, 7.3D-09, 1.6D-01, r^2= 8.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.460227 1 O pz 49 0.566970 2 H s
- 72 0.566970 3 H s 50 -0.361787 2 H s
- 73 -0.361787 3 H s 59 -0.310837 2 H pz
- 82 -0.310837 3 H pz 57 0.283645 2 H px
- 80 -0.283645 3 H px 5 0.279397 1 O s
- Vector 11 Occ=0.000000D+00 E= 1.008120D-01
- MO Center= -3.1D-08, 1.9D-10, -8.2D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.904870 2 H s 73 -2.904870 3 H s
- 49 2.127572 2 H s 72 -2.127572 3 H s
- 15 -1.565803 1 O px 57 -0.924719 2 H px
- 80 -0.924719 3 H px 59 0.911515 2 H pz
- 82 -0.911515 3 H pz 12 -0.548585 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.477602D-01
- MO Center= -1.1D-09, 1.7D-09, 6.4D-01, r^2= 8.7D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.206137 1 O px 49 -2.833342 2 H s
- 72 2.833342 3 H s 50 -2.149102 2 H s
- 73 2.149102 3 H s 12 1.389795 1 O px
- 48 -0.940442 2 H s 71 0.940442 3 H s
- 31 0.579437 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.680148D-01
- MO Center= 7.6D-09, -1.6D-09, -1.8D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.438628 1 O s 49 -1.647344 2 H s
- 72 -1.647344 3 H s 59 -0.987614 2 H pz
- 82 -0.987614 3 H pz 57 -0.611981 2 H px
- 80 0.611981 3 H px 17 0.458137 1 O pz
- 4 0.442679 1 O s 30 -0.424050 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.910788D-01
- MO Center= 1.1D-10, -1.8D-09, -5.0D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 58 1.582359 2 H py 81 -1.582359 3 H py
- 28 -0.471841 1 O d -2 55 0.230555 2 H py
- 78 -0.230555 3 H py 65 -0.185192 2 H d -2
- 88 -0.185192 3 H d -2 66 0.167066 2 H d -1
- 89 -0.167066 3 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.138581D-01
- MO Center= -2.1D-09, -5.7D-09, -2.7D-01, r^2= 6.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.984210 1 O py 58 -1.722278 2 H py
- 81 -1.722278 3 H py 13 1.488935 1 O py
- 29 -0.395482 1 O d -1 55 -0.271923 2 H py
- 78 -0.271923 3 H py 65 0.247015 2 H d -2
- 88 -0.247015 3 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.157182D-01
- MO Center= -9.4D-09, 1.0D-09, 2.3D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.579737 1 O s 49 -3.853565 2 H s
- 72 -3.853565 3 H s 17 -2.513453 1 O pz
- 14 -2.002379 1 O pz 4 1.474124 1 O s
- 48 -1.078351 2 H s 71 -1.078351 3 H s
- 59 0.634250 2 H pz 82 0.634250 3 H pz
- Vector 17 Occ=0.000000D+00 E= 3.080008D-01
- MO Center= -8.0D-08, -6.9D-11, -3.1D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 57 4.030204 2 H px 80 4.030205 3 H px
- 15 -3.421382 1 O px 49 -3.158274 2 H s
- 72 3.158276 3 H s 59 1.763771 2 H pz
- 82 -1.763771 3 H pz 12 -1.198876 1 O px
- 50 0.836072 2 H s 73 -0.836072 3 H s
- Vector 18 Occ=0.000000D+00 E= 3.127818D-01
- MO Center= 1.3D-07, 2.9D-10, -8.2D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.719547 1 O s 49 -4.883678 2 H s
- 72 -4.883677 3 H s 4 4.138463 1 O s
- 57 2.086512 2 H px 80 -2.086510 3 H px
- 59 -1.791194 2 H pz 82 -1.791194 3 H pz
- 50 -0.752091 2 H s 73 -0.752091 3 H s
- Vector 19 Occ=0.000000D+00 E= 3.248071D-01
- MO Center= -5.4D-08, 5.8D-10, -9.0D-01, r^2= 7.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 7.816029 2 H s 72 -7.816029 3 H s
- 57 -3.004510 2 H px 80 -3.004510 3 H px
- 59 2.157839 2 H pz 82 -2.157839 3 H pz
- 50 2.014907 2 H s 73 -2.014907 3 H s
- 15 -1.957354 1 O px 12 -1.657851 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.552484D-01
- MO Center= 1.3D-08, 7.3D-10, 2.2D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 9.670633 2 H s 72 -9.670633 3 H s
- 12 -3.628126 1 O px 15 -3.030710 1 O px
- 57 -2.461533 2 H px 80 -2.461533 3 H px
- 31 -1.960931 1 O d 1 48 1.921420 2 H s
- 71 -1.921420 3 H s 50 1.499638 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.950339D-01
- MO Center= -5.9D-09, 3.8D-11, -6.4D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 6.945897 1 O s 49 -4.428543 2 H s
- 72 -4.428543 3 H s 4 1.573254 1 O s
- 32 1.195068 1 O d 2 59 -1.121165 2 H pz
- 82 -1.121165 3 H pz 17 -0.930025 1 O pz
- 30 -0.898818 1 O d 0 14 0.868065 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.450711D-01
- MO Center= -2.1D-09, 6.5D-10, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.134059 1 O s 49 -1.987817 2 H s
- 72 -1.987817 3 H s 14 -1.955360 1 O pz
- 17 -1.890979 1 O pz 30 1.609753 1 O d 0
- 59 1.012002 2 H pz 82 1.012002 3 H pz
- 4 0.987503 1 O s 57 0.610790 2 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10009567
- moments of inertia (a.u.)
- ------------------
- 2.277381949965 0.000000000000 0.000000000830
- 0.000000000000 6.502988861392 0.000000000000
- 0.000000000830 0.000000000000 4.225606911426
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 0.000000 0.000000 0.000000 0.000000
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.711192 -0.359304 -0.359304 0.007416
- 2 2 0 0 -3.365482 -3.779140 -3.779140 4.192798
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 0.000000 0.000000 0.000000 0.000000
- 2 0 2 0 -5.847499 -2.923750 -2.923750 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.850140 -3.442858 -3.442858 2.035576
- HESSIAN: the one electron contributions are done in 0.1s
- HESSIAN: 2-el 1st deriv. term done in 3.9s
- HESSIAN: 2-el 2nd deriv. term done in 15.5s
- stpr_wrt_fd_from_sq: overwrite of existing file:/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.hess
- stpr_wrt_fd_dipole: overwrite of existing file/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.fd_ddipole
- HESSIAN: the two electron contributions are done in 31.0s
- NWChem CPHF Module
- ------------------
- scftype = RHF
- nclosed = 5
- nopen = 0
- variables = 435
- # of vectors = 9
- tolerance = 0.10D-03
- level shift = 0.00D+00
- max iterations = 50
- max subspace = 90
- SCF residual: 1.424635967525494E-006
- Iterative solution of linear equations
- No. of variables 435
- No. of equations 9
- Maximum subspace 90
- Iterations 50
- Convergence 1.0D-04
- Start time 79.4
- iter nsub residual time
- ---- ------ -------- ---------
- 1 9 1.46D-01 86.6
- 2 18 7.06D-02 94.4
- 3 27 1.62D-03 102.4
- 4 36 7.05D-05 110.3
- HESSIAN: the CPHF contributions are done
- stpr_wrt_fd_from_sq: overwrite of existing file:/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.hess
- stpr_wrt_fd_dipole: overwrite of existing file/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.fd_ddipole
- -----------------------
- Performance information
- -----------------------
- Timer overhead = 3.00D-07 seconds/call
- Nr. of calls CPU time (s) Wall time (s) GFlops
- --------------- ------------------- ------------------------------ -------------------
- Name Min Avg Max Min Avg Max Min Avg Max Mx/calls Min Max Sum
- hess: tot 1 1 1 53. 57. 59. 76. 76. 76. 76. 0.0 0.0 0.0
- hess: int 1 1 1 5.70E-2 6.67E-2 7.30E-2 7.88E-2 7.88E-2 7.89E-2 7.89E-2 0.0 0.0 0.0
- hess: nuc 1 1 1 0.0 0.0 0.0 9.54E-7 1.73E-6 2.15E-6 2.15E-6 0.0 0.0 0.0
- hess: int 1 1 1 19. 24. 27. 31. 31. 31. 31. 0.0 0.0 0.0
- hess: cph 1 1 1 29. 30. 31. 41. 41. 41. 41. 0.0 0.0 0.0
- The average no. of pstat calls per process was 5.00D+00
- with a timing overhead of 1.50D-06s
- HESSIAN: the Hessian is done
- Vibrational analysis via the FX method
- See chapter 2 in "Molecular Vibrations" by Wilson, Decius and Cross
- Vib: Default input used
- Nuclear Hessian passed symmetry test
- ---------------------------- Atom information ----------------------------
- atom # X Y Z mass
- --------------------------------------------------------------------------
- O 1 -1.9804391D-10 0.0000000D+00 2.2632811D-01 1.5994910D+01
- H 2 1.4478947D+00 0.0000000D+00 -9.0160432D-01 1.0078250D+00
- H 3 -1.4478947D+00 0.0000000D+00 -9.0160432D-01 1.0078250D+00
- --------------------------------------------------------------------------
- ----------------------------------------------------
- MASS-WEIGHTED NUCLEAR HESSIAN (Hartree/Bohr/Bohr/Kamu)
- ----------------------------------------------------
- 1 2 3 4 5 6 7 8 9
- ----- ----- ----- ----- -----
- 1 3.93638D+01
- 2 1.75061D-06-6.29578D-02
- 3 1.97812D-06 2.45049D-06 2.63943D+01
- 4 -7.85208D+01 5.11574D-08 4.65814D+01 3.42534D+02
- 5 6.57097D-08-6.71233D-04 2.00149D-07 3.87999D-07-1.04376D-02
- 6 6.11684D+01 8.48166D-08-5.26966D+01-2.14634D+02-2.98543D-09 1.99873D+02
- 7 -7.85208D+01-1.89584D-08-4.65814D+01-2.97248D+01-1.42822D-07-2.90549D+01 3.42534D+02
- 8 1.99875D-07-6.70996D-04 3.95664D-07 3.75723D-07 1.77511D-03-2.51595D-07-4.37406D-07-1.04377D-02
- 9 -6.11684D+01 1.20706D-07-5.26966D+01 2.90567D+01-7.86112D-08 1.00493D+01 2.14634D+02-5.78446D-08 1.99873D+02
- -------------------------------------------------
- NORMAL MODE EIGENVECTORS IN CARTESIAN COORDINATES
- -------------------------------------------------
- (Frequencies expressed in cm-1)
- 1 2 3 4 5 6
- Frequency -40.79 -38.26 -37.04 -18.33 -17.96 -15.12
- 1 -0.00004 0.00440 -0.24074 -0.00117 0.00000 0.00000
- 2 0.25000 -0.00008 -0.00004 0.00000 0.00000 0.00437
- 3 -0.00007 -0.23560 -0.00430 -0.00100 0.00000 0.00000
- 4 -0.00002 0.00457 -0.14766 -0.43358 0.00017 0.00000
- 5 0.01231 0.00000 0.00000 0.00029 0.70435 -0.70426
- 6 -0.00005 -0.23534 0.11515 -0.55607 0.00022 0.00000
- 7 -0.00002 0.00450 -0.14766 -0.43358 0.00017 0.00000
- 8 0.01230 0.00000 0.00000 -0.00028 -0.70437 -0.70424
- 9 -0.00009 -0.23570 -0.12374 0.55407 -0.00022 0.00000
- 7 8 9
- Frequency 1592.72 3680.50 3784.67
- 1 0.00000 0.00000 -0.06742
- 2 0.00000 0.00000 0.00000
- 3 -0.06809 0.04855 0.00000
- 4 0.40884 0.57355 0.53507
- 5 0.00000 0.00000 0.00000
- 6 0.54054 -0.38529 -0.41682
- 7 -0.40884 -0.57356 0.53507
- 8 0.00000 0.00000 0.00000
- 9 0.54054 -0.38529 0.41683
- ----------------------------------------------------------------------------
- Normal Eigenvalue || Derivative Dipole Moments (debye/angs)
- Mode [cm**-1] || [d/dqX] [d/dqY] [d/dqZ]
- ------ ---------- || ------------------ ------------------ -----------------
- 1 -40.793 || 0.000 -0.720 0.000
- 2 -38.262 || 0.000 0.000 0.000
- 3 -37.037 || -0.282 0.000 0.000
- 4 -18.326 || -1.310 0.000 0.000
- 5 -17.964 || -0.001 0.000 0.000
- 6 -15.115 || 0.000 2.146 0.000
- 7 1592.720 || 0.000 0.000 1.260
- 8 3680.501 || 0.000 0.000 0.199
- 9 3784.675 || 1.071 0.000 0.000
- ----------------------------------------------------------------------------
- ----------------------------------------------------------------------------
- Normal Eigenvalue || Infra Red Intensities
- Mode [cm**-1] || [atomic units] [(debye/angs)**2] [(KM/mol)] [arbitrary]
- ------ ---------- || -------------- ----------------- ---------- -----------
- 1 -40.793 || 0.022462 0.518 21.897 4.812
- 2 -38.262 || 0.000000 0.000 0.000 0.000
- 3 -37.037 || 0.003444 0.079 3.358 0.738
- 4 -18.326 || 0.074337 1.715 72.467 15.926
- 5 -17.964 || 0.000000 0.000 0.000 0.000
- 6 -15.115 || 0.199621 4.605 194.600 42.768
- 7 1592.720 || 0.068761 1.586 67.032 14.732
- 8 3680.501 || 0.001711 0.039 1.668 0.367
- 9 3784.675 || 0.049744 1.148 48.493 10.657
- ----------------------------------------------------------------------------
- Vibrational analysis via the FX method
- --- with translations and rotations projected out ---
- --- via the Eckart algorithm ---
- Projected Nuclear Hessian trans-rot subspace norm:4.2370D-33
- (should be close to zero!)
- From the projected analysis
- The Zero-Point Energy (Kcal/mol) = 12.94946962
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10009567
- moments of inertia (a.u.)
- ------------------
- 2.277381949965 0.000000000000 0.000000000830
- 0.000000000000 6.502988861392 0.000000000000
- 0.000000000830 0.000000000000 4.225606911426
- Rotational Constants
- --------------------
- A= 26.433939 cm-1 ( 38.031692 K)
- B= 14.246516 cm-1 ( 20.497100 K)
- C= 9.257309 cm-1 ( 13.318905 K)
- Temperature = 298.15K
- frequency scaling parameter = 1.0000
- Zero-Point correction to Energy = 12.939 kcal/mol ( 0.020620 au)
- Thermal correction to Energy = 14.718 kcal/mol ( 0.023455 au)
- Thermal correction to Enthalpy = 15.310 kcal/mol ( 0.024399 au)
- Total Entropy = 46.508 cal/mol-K
- - Translational = 34.593 cal/mol-K (mol. weight = 18.0106)
- - Rotational = 11.908 cal/mol-K (symmetry # = 1)
- - Vibrational = 0.008 cal/mol-K
- Cv (constant volume heat capacity) = 6.013 cal/mol-K
- - Translational = 2.979 cal/mol-K
- - Rotational = 2.979 cal/mol-K
- - Vibrational = 0.054 cal/mol-K
- -------------------------------------------------
- NORMAL MODE EIGENVECTORS IN CARTESIAN COORDINATES
- -------------------------------------------------
- (Projected Frequencies expressed in cm-1)
- 1 2 3 4 5 6
- P.Frequency 0.00 0.00 0.00 0.00 0.00 0.00
- 1 0.00000 0.00000 0.00000 -0.23578 0.01756 -0.04552
- 2 0.25004 0.00000 0.00000 0.00000 0.00000 0.00000
- 3 0.00000 0.00000 0.00000 0.00000 -0.21984 -0.08481
- 4 0.00000 0.00000 0.00000 -0.23445 -0.14164 0.36714
- 5 0.00000 0.99611 0.00000 0.00000 0.00000 0.00000
- 6 0.00000 0.00000 0.00000 0.00170 -0.42420 0.44491
- 7 0.00000 0.00000 0.00000 -0.23445 -0.14164 0.36714
- 8 0.00000 0.00000 0.99611 0.00000 0.00000 0.00000
- 9 0.00000 0.00000 0.00000 -0.00171 -0.01548 -0.61454
- 7 8 9
- P.Frequency 1591.50 3680.21 3783.86
- 1 0.00000 0.00000 -0.06743
- 2 0.00000 0.00000 0.00000
- 3 -0.06813 0.04854 0.00000
- 4 0.40869 0.57366 0.53507
- 5 0.00000 0.00000 0.00000
- 6 0.54061 -0.38514 -0.41683
- 7 -0.40869 -0.57366 0.53506
- 8 0.00000 0.00000 0.00000
- 9 0.54061 -0.38515 0.41682
- ----------------------------------------------------------------------------
- Normal Eigenvalue || Projected Derivative Dipole Moments (debye/angs)
- Mode [cm**-1] || [d/dqX] [d/dqY] [d/dqZ]
- ------ ---------- || ------------------ ------------------ -----------------
- 1 0.000 || 0.000 -0.757 0.000
- 2 0.000 || 0.000 1.508 0.000
- 3 0.000 || 0.000 1.508 0.000
- 4 0.000 || -0.004 0.000 0.000
- 5 0.000 || 0.482 0.000 0.000
- 6 0.000 || -1.250 0.000 0.000
- 7 1591.504 || 0.000 0.000 1.260
- 8 3680.211 || 0.000 0.000 0.198
- 9 3783.859 || 1.071 0.000 0.000
- ----------------------------------------------------------------------------
- ----------------------------------------------------------------------------
- Normal Eigenvalue || Projected Infra Red Intensities
- Mode [cm**-1] || [atomic units] [(debye/angs)**2] [(KM/mol)] [arbitrary]
- ------ ---------- || -------------- ----------------- ---------- -----------
- 1 0.000 || 0.024856 0.573 24.231 5.325
- 2 0.000 || 0.098614 2.275 96.133 21.127
- 3 0.000 || 0.098614 2.275 96.133 21.127
- 4 0.000 || 0.000001 0.000 0.001 0.000
- 5 0.000 || 0.010076 0.232 9.823 2.159
- 6 0.000 || 0.067704 1.562 66.001 14.505
- 7 1591.504 || 0.068767 1.587 67.037 14.733
- 8 3680.211 || 0.001705 0.039 1.662 0.365
- 9 3783.859 || 0.049744 1.148 48.493 10.657
- ----------------------------------------------------------------------------
- vib:animation F
- FA nfreq= 1
- (mode_ini,mode_end)=( 1, 3)
- task_raman_doit: no rminfo found
- (begin,last,steps)=( 7, 9, 2)
- ENTER raman default
- Generating fname_raman3=/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.raman_0007-0009
- (ii,rminfo(ii,1))=( 7, 0.00000000)
- ramanpolfile=/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0001-1
- ramanpol_read: failed to open /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0001-1
- NWChem Property Module
- ----------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 5.00D-04
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: fine
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 70 10.0 590
- H 0.35 60 13.0 590
- Grid pruning is: on
- Number of quadrature shells: 190
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Loading old vectors from job with title :
- TDDFT h2o BP86/aug-cc-pVTZ
- Time after variat. SCF: 81.9
- Time prior to 1st pass: 81.9
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 12 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 12.96 12956396
- Stack Space remaining (MW): 13.11 13106484
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4660846408 -8.55D+01 4.44D-05 2.68D-05 83.0
- d= 0,ls=0.0,diis 2 -76.4660876222 -2.98D-06 1.10D-05 2.15D-06 84.2
- d= 0,ls=0.0,diis 3 -76.4660873233 2.99D-07 7.03D-06 4.56D-06 85.3
- d= 0,ls=0.0,diis 4 -76.4660877958 -4.72D-07 8.61D-07 9.05D-08 86.4
- d= 0,ls=0.0,diis 5 -76.4660878031 -7.29D-09 5.29D-08 9.90D-11 87.6
- Total DFT energy = -76.466087803064
- One electron energy = -122.766074961697
- Coulomb energy = 46.548001608885
- Exchange-Corr. energy = -9.312284467854
- Nuclear repulsion energy = 9.064270017602
- Numeric. integr. density = 9.999999990913
- Total iterative time = 5.6s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878156D+01
- MO Center= 9.1D-11, 2.0D-13, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.974039 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.310549D-01
- MO Center= 1.6D-09, 4.5D-11, -8.4D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.482679 1 O s 4 0.280803 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.846473D-01
- MO Center= -2.3D-09, -2.7D-11, -9.6D-02, r^2= 7.8D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.324835 1 O px 6 0.300895 1 O px
- 48 0.226612 2 H s 71 -0.226612 3 H s
- Vector 4 Occ=2.000000D+00 E=-3.460614D-01
- MO Center= -4.8D-10, -3.8D-10, 1.8D-01, r^2= 7.5D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.359929 1 O pz 4 0.344633 1 O s
- 8 0.320424 1 O pz 14 0.206797 1 O pz
- 2 0.199454 1 O s 48 -0.175086 2 H s
- 71 -0.175086 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.695111D-01
- MO Center= 5.6D-11, -6.6D-11, 7.2D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428202 1 O py 7 0.365037 1 O py
- 13 0.324307 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.228939D-02
- MO Center= -1.3D-08, 2.2D-09, -7.7D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.484977 1 O s 50 -0.463654 2 H s
- 73 -0.463654 3 H s 4 0.367667 1 O s
- 49 -0.355567 2 H s 72 -0.355567 3 H s
- 14 -0.190856 1 O pz 2 0.151172 1 O s
- Vector 7 Occ=0.000000D+00 E= 1.507486D-02
- MO Center= 4.8D-08, 2.4D-10, -7.6D-01, r^2= 1.5D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.627184 2 H s 73 -2.627184 3 H s
- 49 0.804206 2 H s 72 -0.804206 3 H s
- 15 -0.647113 1 O px 12 -0.342938 1 O px
- 48 0.159152 2 H s 71 -0.159152 3 H s
- Vector 8 Occ=0.000000D+00 E= 6.603930D-02
- MO Center= -2.3D-07, 2.7D-09, -6.3D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.400941 1 O s 49 -2.129439 2 H s
- 72 -2.129439 3 H s 4 1.371677 1 O s
- 17 -1.131922 1 O pz 50 -1.022908 2 H s
- 73 -1.022908 3 H s 57 0.690995 2 H px
- 80 -0.690995 3 H px 59 -0.538496 2 H pz
- Vector 9 Occ=0.000000D+00 E= 9.415274D-02
- MO Center= -1.4D-10, -5.8D-09, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.497439 1 O py 58 -0.312304 2 H py
- 81 -0.312304 3 H py
- Vector 10 Occ=0.000000D+00 E= 9.874046D-02
- MO Center= 3.0D-07, 7.4D-09, 1.6D-01, r^2= 8.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.462375 1 O pz 49 0.568888 2 H s
- 72 0.568890 3 H s 50 -0.360206 2 H s
- 73 -0.360205 3 H s 59 -0.311304 2 H pz
- 82 -0.311304 3 H pz 57 0.282599 2 H px
- 80 -0.282598 3 H px 4 0.269562 1 O s
- Vector 11 Occ=0.000000D+00 E= 1.009203D-01
- MO Center= -1.1D-07, 1.9D-10, -8.2D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.908073 2 H s 73 -2.908073 3 H s
- 49 2.133127 2 H s 72 -2.133127 3 H s
- 15 -1.583525 1 O px 57 -0.922265 2 H px
- 80 -0.922265 3 H px 59 0.910997 2 H pz
- 82 -0.910997 3 H pz 12 -0.556747 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.478961D-01
- MO Center= -9.2D-08, 1.6D-09, 6.4D-01, r^2= 8.7D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.202914 1 O px 49 -2.827581 2 H s
- 72 2.827581 3 H s 50 -2.132945 2 H s
- 73 2.132945 3 H s 12 1.395089 1 O px
- 48 -0.944943 2 H s 71 0.944943 3 H s
- 31 0.582203 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.678337D-01
- MO Center= -1.2D-07, -1.6D-09, -1.8D-01, r^2= 6.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.508819 1 O s 49 -1.677898 2 H s
- 72 -1.677897 3 H s 59 -0.986644 2 H pz
- 82 -0.986644 3 H pz 57 -0.601630 2 H px
- 80 0.601631 3 H px 4 0.459663 1 O s
- 17 0.446791 1 O pz 30 -0.420942 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.909144D-01
- MO Center= -1.4D-08, -1.8D-09, -5.0D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 58 1.579485 2 H py 81 -1.579485 3 H py
- 28 -0.474109 1 O d -2 55 0.231332 2 H py
- 78 -0.231332 3 H py 65 -0.185247 2 H d -2
- 88 -0.185247 3 H d -2 66 0.166796 2 H d -1
- 89 -0.166796 3 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.142990D-01
- MO Center= 1.3D-08, -5.9D-09, -2.7D-01, r^2= 6.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.989332 1 O py 58 -1.727310 2 H py
- 81 -1.727310 3 H py 13 1.492922 1 O py
- 29 -0.394502 1 O d -1 55 -0.272275 2 H py
- 78 -0.272275 3 H py 65 0.248557 2 H d -2
- 88 -0.248557 3 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.160347D-01
- MO Center= 1.3D-07, 1.3D-09, 2.3D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.559447 1 O s 49 -3.846654 2 H s
- 72 -3.846657 3 H s 17 -2.516796 1 O pz
- 14 -2.004780 1 O pz 4 1.477259 1 O s
- 48 -1.074624 2 H s 71 -1.074625 3 H s
- 59 0.648922 2 H pz 82 0.648921 3 H pz
- Vector 17 Occ=0.000000D+00 E= 3.079040D-01
- MO Center= 1.7D-06, -6.7D-11, -3.2D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 57 3.982523 2 H px 80 3.982505 3 H px
- 15 -3.446830 1 O px 49 -3.027298 2 H s
- 72 3.027255 3 H s 59 1.776349 2 H pz
- 82 -1.776364 3 H pz 12 -1.212204 1 O px
- 50 0.857211 2 H s 73 -0.857218 3 H s
- Vector 18 Occ=0.000000D+00 E= 3.128354D-01
- MO Center= -2.7D-06, 2.8D-10, -8.1D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.716932 1 O s 49 -4.883121 2 H s
- 72 -4.883150 3 H s 4 4.148240 1 O s
- 57 2.094049 2 H px 80 -2.094084 3 H px
- 59 -1.782811 2 H pz 82 -1.782798 3 H pz
- 50 -0.751737 2 H s 73 -0.751731 3 H s
- Vector 19 Occ=0.000000D+00 E= 3.250242D-01
- MO Center= 1.1D-06, 5.7D-10, -8.9D-01, r^2= 7.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 7.837745 2 H s 72 -7.837741 3 H s
- 57 -3.040582 2 H px 80 -3.040580 3 H px
- 59 2.134431 2 H pz 82 -2.134430 3 H pz
- 50 2.006948 2 H s 73 -2.006947 3 H s
- 15 -1.938599 1 O px 12 -1.652300 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.553489D-01
- MO Center= -2.3D-07, 7.1D-10, 2.3D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 9.604990 2 H s 72 -9.604989 3 H s
- 12 -3.648481 1 O px 15 -3.036853 1 O px
- 57 -2.435864 2 H px 80 -2.435864 3 H px
- 31 -1.972948 1 O d 1 48 1.938026 2 H s
- 71 -1.938026 3 H s 50 1.493405 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.944602D-01
- MO Center= 1.1D-07, 2.4D-11, -6.4D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.014892 1 O s 49 -4.468637 2 H s
- 72 -4.468641 3 H s 4 1.605492 1 O s
- 32 1.192218 1 O d 2 59 -1.123547 2 H pz
- 82 -1.123547 3 H pz 17 -0.933857 1 O pz
- 30 -0.891096 1 O d 0 14 0.852753 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.454550D-01
- MO Center= 7.5D-08, 6.3D-10, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.061304 1 O s 14 -1.940624 1 O pz
- 49 -1.942484 2 H s 72 -1.942483 3 H s
- 17 -1.880978 1 O pz 30 1.603618 1 O d 0
- 59 1.022395 2 H pz 82 1.022395 3 H pz
- 4 0.962826 1 O s 57 0.601602 2 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10009567
- moments of inertia (a.u.)
- ------------------
- 2.252866621448 0.000000000000 -0.000000012475
- 0.000000000000 6.502362218946 0.000000000000
- -0.000000012475 0.000000000000 4.249495597498
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 0.000000 0.000000 0.000000 0.000000
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.708562 -0.360670 -0.360670 0.012778
- 2 2 0 0 -3.351539 -3.784020 -3.784020 4.216501
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 0.000000 0.000000 0.000000 0.000000
- 2 0 2 0 -5.847697 -2.923848 -2.923848 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.861033 -3.437354 -3.437354 2.013675
- Entering AOResponse driver routine
- --------------------------------------------------------------------------------
- ****************
- *** RESPONSE ***
- ****************
- Response module for NWChem and dynamic CPKS solver
- developed by J. Autschbach and coworkers, SUNY Buffalo
- The methodology used in this program is described in
- J. Chem. Phys. 123 (2005), 114103
- J. Chem. Phys. 122 (2005), 224115
- J. Chem. Phys. 122 (2005), 074105
- Comp. Lett. 3 (2007), 131-150 (contact JA for a copy)
- Please cite this work in publications based on results
- obtained with this code. Thank you!
- For updated response module please acknowledge:
- F. Aquino, Northwestern University, Schatz Rsrch Group
- (version June 9th,2012)
- ++++++++++++++++++++++++++++++++++++++++++++++++++++++
- -----------------------------------------------
- Solving response equations for perturbing field
- -----------------------------------------------
- number of frequencies: 1
- frequency in a.u.: 0.8855900E-01
- Perturbing field: electric
- Using Dipole Length Gauge
- Applying DAMPING constant of 0.007000 a.u.
- Setting up CPKS
- with frequency omega = 0.08855900 a.u.
- and a damping of 0.00700000 a.u.
- NWChem Dynamic CPHF Module
- --------------------------
- scftype = RHF
- nclosed = 5
- nopen = 0
- variables = 1740
- num. vecs = 3
- tolerance = 1.0D-04
- max. iter = 50
- SCF residual: 1.160124836524185E-006
- Iterative solution of linear equations
- No. of variables 870
- No. of equations 3
- Maximum subspace 33
- Convergence 1.0D-04
- Start time 118.4
- iter nsub residual time
- ---- ------ -------- ---------
- dft_CPHF2_read: failed to open/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1
- File /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1 does not exist, proceed to generate (z1,Az1)
- 1 3 2.11D+00 127.4
- 2 6 1.85D+00 137.5
- 3 9 1.60D-01 147.5
- 4 12 1.62D-02 157.4
- 5 15 4.65D-04 166.8
- 6 18 6.24D-05 177.3
- Electric Dipole Response Matrix (nonzero elements):
- 1 2 3
- 1 11.0521 0.0000 0.0000
- 2 0.0000 10.7154 0.0000
- 3 0.0000 0.0000 10.7262
- DFT Linear Response polarizability / au
- Frequency = 0.0885590 / au
- Wavelength = 514.4971488 / nm
- X Y Z
- -----------------------------------------------
- X 11.0521361 0.0000000 0.0000000
- Y 0.0000000 10.7153775 0.0000000
- Z 0.0000000 0.0000000 10.7261865
- -----------------------------------------------
- Eigenvalues = 11.0521361 10.7153775 10.7261865
- Isotropic = 10.8312334
- Anisotropic = 0.3314863
- -----------------------------------------------
- task_raman_doit: no rminfo found
- task_raman_doit: no rminfo found
- task_raman_doit: no rminfo found
- Imaginary part:
- 1 2 3
- 1 0.1217 0.0000 0.0000
- 2 0.0000 0.1910 0.0000
- 3 0.0000 0.0000 0.1412
- average: 10.8312333700002 + I 0.151285913811853
- Magnetic Dipole Response Matrix (nonzero elements):
- Optical rotation tensor Beta
- 1 2 3
- 1 0.0000 0.4351 0.0000
- 2 0.7364 0.0000 0.0000
- 3 0.0000 0.0000 0.0000
- Imaginary part:
- 1 2 3
- 1 0.0000 0.0157 0.0000
- 2 0.0552 0.0000 0.0000
- 3 0.0000 0.0000 0.0000
- average: -9.218171754989208E-011 + I -5.045882332644614E-012
- Exiting AOResponse driver routine
- ramanpolfile=/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0001-2
- ramanpol_read: failed to open /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0001-2
- NWChem Property Module
- ----------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 5.00D-04
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: fine
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 70 10.0 590
- H 0.35 60 13.0 590
- Grid pruning is: on
- Number of quadrature shells: 190
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Loading old vectors from job with title :
- TDDFT h2o BP86/aug-cc-pVTZ
- Time after variat. SCF: 133.9
- Time prior to 1st pass: 133.9
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 12 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 12.96 12956396
- Stack Space remaining (MW): 13.11 13106484
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4660749991 -8.55D+01 8.86D-05 1.07D-04 135.0
- d= 0,ls=0.0,diis 2 -76.4660869410 -1.19D-05 2.17D-05 8.18D-06 136.1
- d= 0,ls=0.0,diis 3 -76.4660861423 7.99D-07 1.38D-05 1.76D-05 137.3
- d= 0,ls=0.0,diis 4 -76.4660876041 -1.46D-06 1.72D-06 3.58D-07 138.3
- d= 0,ls=0.0,diis 5 -76.4660876330 -2.89D-08 1.05D-07 3.93D-10 139.5
- d= 0,ls=0.0,diis 6 -76.4660876330 -4.36D-11 1.68D-08 2.64D-12 140.6
- Total DFT energy = -76.466087633041
- One electron energy = -122.758816379244
- Coulomb energy = 46.543339683172
- Exchange-Corr. energy = -9.311920540679
- Nuclear repulsion energy = 9.061309603711
- Numeric. integr. density = 10.000000011609
- Total iterative time = 6.8s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878197D+01
- MO Center= -3.0D-10, 9.6D-15, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.974024 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.314530D-01
- MO Center= -1.9D-09, 2.3D-12, -8.5D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.482371 1 O s 4 0.280569 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.828553D-01
- MO Center= 2.4D-09, -1.0D-12, -9.7D-02, r^2= 7.8D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.325425 1 O px 6 0.301140 1 O px
- 48 0.225898 2 H s 71 -0.225898 3 H s
- Vector 4 Occ=2.000000D+00 E=-3.477059D-01
- MO Center= 1.2D-10, -1.8D-11, 1.8D-01, r^2= 7.5D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.358928 1 O pz 4 0.345906 1 O s
- 8 0.319709 1 O pz 14 0.205273 1 O pz
- 2 0.201031 1 O s 48 -0.175787 2 H s
- 71 -0.175787 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.696351D-01
- MO Center= -3.7D-10, -5.3D-12, 7.3D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428249 1 O py 7 0.365077 1 O py
- 13 0.324430 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.245656D-02
- MO Center= 1.8D-08, 1.7D-10, -7.8D-01, r^2= 7.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.477751 1 O s 50 -0.462213 2 H s
- 73 -0.462213 3 H s 4 0.364770 1 O s
- 49 -0.352243 2 H s 72 -0.352242 3 H s
- 14 -0.192618 1 O pz 2 0.150726 1 O s
- Vector 7 Occ=0.000000D+00 E= 1.511567D-02
- MO Center= -5.7D-08, 1.6D-11, -7.7D-01, r^2= 1.5D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.639094 2 H s 73 -2.639094 3 H s
- 49 0.805832 2 H s 72 -0.805832 3 H s
- 15 -0.644484 1 O px 12 -0.341063 1 O px
- 48 0.158751 2 H s 71 -0.158751 3 H s
- Vector 8 Occ=0.000000D+00 E= 6.589044D-02
- MO Center= 2.4D-07, 2.5D-10, -7.0D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.381583 1 O s 49 -2.120499 2 H s
- 72 -2.120498 3 H s 4 1.363685 1 O s
- 17 -1.133513 1 O pz 50 -1.021645 2 H s
- 73 -1.021645 3 H s 57 0.681044 2 H px
- 80 -0.681044 3 H px 59 -0.546657 2 H pz
- Vector 9 Occ=0.000000D+00 E= 9.424099D-02
- MO Center= 1.8D-09, -4.4D-10, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.490832 1 O py 58 -0.306624 2 H py
- 81 -0.306624 3 H py
- Vector 10 Occ=0.000000D+00 E= 9.900618D-02
- MO Center= -3.4D-07, 6.1D-10, 1.6D-01, r^2= 8.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.458219 1 O pz 49 0.565393 2 H s
- 72 0.565391 3 H s 50 -0.363245 2 H s
- 73 -0.363246 3 H s 59 -0.310295 2 H pz
- 82 -0.310295 3 H pz 5 0.288233 1 O s
- 57 0.284575 2 H px 80 -0.284575 3 H px
- Vector 11 Occ=0.000000D+00 E= 1.006990D-01
- MO Center= 1.4D-07, 9.4D-12, -8.2D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.901587 2 H s 73 -2.901587 3 H s
- 49 2.121984 2 H s 72 -2.121984 3 H s
- 15 -1.547994 1 O px 57 -0.927181 2 H px
- 80 -0.927181 3 H px 59 0.911924 2 H pz
- 82 -0.911924 3 H pz 12 -0.540370 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.476263D-01
- MO Center= 9.8D-08, 1.3D-10, 6.4D-01, r^2= 8.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.209206 1 O px 49 -2.839177 2 H s
- 72 2.839177 3 H s 50 -2.165322 2 H s
- 73 2.165322 3 H s 12 1.384445 1 O px
- 48 -0.935919 2 H s 71 0.935919 3 H s
- 31 0.576688 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.681924D-01
- MO Center= 1.3D-07, -1.4D-10, -1.9D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.367759 1 O s 49 -1.616487 2 H s
- 72 -1.616489 3 H s 59 -0.988494 2 H pz
- 82 -0.988494 3 H pz 57 -0.622213 2 H px
- 80 0.622212 3 H px 17 0.469713 1 O pz
- 4 0.425781 1 O s 30 -0.427173 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.912430D-01
- MO Center= 1.5D-08, -1.3D-10, -5.1D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 58 1.585242 2 H py 81 -1.585242 3 H py
- 28 -0.469545 1 O d -2 55 0.229766 2 H py
- 78 -0.229766 3 H py 65 -0.185143 2 H d -2
- 88 -0.185143 3 H d -2 66 0.167319 2 H d -1
- 89 -0.167319 3 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.134186D-01
- MO Center= -1.6D-08, -8.6D-10, -2.7D-01, r^2= 6.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.979118 1 O py 58 -1.717255 2 H py
- 81 -1.717255 3 H py 13 1.484924 1 O py
- 29 -0.396450 1 O d -1 55 -0.271560 2 H py
- 78 -0.271560 3 H py 65 0.245474 2 H d -2
- 88 -0.245474 3 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.154064D-01
- MO Center= -1.4D-07, 4.5D-10, 2.3D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.599866 1 O s 49 -3.860412 2 H s
- 72 -3.860409 3 H s 17 -2.510010 1 O pz
- 14 -1.999753 1 O pz 4 1.470895 1 O s
- 48 -1.081948 2 H s 71 -1.081947 3 H s
- 59 0.619469 2 H pz 82 0.619471 3 H pz
- Vector 17 Occ=0.000000D+00 E= 3.080932D-01
- MO Center= -2.2D-06, -9.8D-12, -3.0D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 57 4.078617 2 H px 80 4.078638 3 H px
- 15 -3.395112 1 O px 49 -3.292449 2 H s
- 72 3.292498 3 H s 59 1.750130 2 H pz
- 82 -1.750112 3 H pz 12 -1.184533 1 O px
- 50 0.814145 2 H s 73 -0.814137 3 H s
- Vector 18 Occ=0.000000D+00 E= 3.127172D-01
- MO Center= 3.5D-06, 2.4D-11, -8.2D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.721375 1 O s 49 -4.883765 2 H s
- 72 -4.883727 3 H s 4 4.128398 1 O s
- 57 2.078890 2 H px 80 -2.078848 3 H px
- 59 -1.799527 2 H pz 82 -1.799542 3 H pz
- 50 -0.752421 2 H s 73 -0.752427 3 H s
- Vector 19 Occ=0.000000D+00 E= 3.245879D-01
- MO Center= -1.3D-06, 3.2D-11, -9.2D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 7.791530 2 H s 72 -7.791534 3 H s
- 57 -2.966445 2 H px 80 -2.966447 3 H px
- 59 2.181728 2 H pz 82 -2.181730 3 H pz
- 50 2.022944 2 H s 73 -2.022945 3 H s
- 15 -1.976988 1 O px 12 -1.663721 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.551516D-01
- MO Center= 2.5D-07, 2.4D-11, 2.2D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 9.737082 2 H s 72 -9.737083 3 H s
- 12 -3.607933 1 O px 15 -3.024516 1 O px
- 57 -2.487402 2 H px 80 -2.487402 3 H px
- 31 -1.948971 1 O d 1 48 1.904916 2 H s
- 71 -1.904916 3 H s 50 1.505867 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.956109D-01
- MO Center= -1.2D-07, -8.4D-12, -6.4D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 6.877168 1 O s 49 -4.388637 2 H s
- 72 -4.388632 3 H s 4 1.541558 1 O s
- 32 1.197919 1 O d 2 59 -1.118764 2 H pz
- 82 -1.118764 3 H pz 17 -0.926064 1 O pz
- 30 -0.906575 1 O d 0 14 0.883227 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.446723D-01
- MO Center= -8.3D-08, 2.2D-11, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.206601 1 O s 49 -2.032976 2 H s
- 72 -2.032978 3 H s 14 -1.969310 1 O pz
- 17 -1.900855 1 O pz 30 1.615461 1 O d 0
- 4 1.011736 1 O s 59 1.001249 2 H pz
- 82 1.001250 3 H pz 57 0.619761 2 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10009567
- moments of inertia (a.u.)
- ------------------
- 2.302029943550 0.000000000000 0.000000014169
- 0.000000000000 6.503815503837 0.000000000000
- 0.000000014169 0.000000000000 4.201785560287
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 0.000000 0.000000 0.000000 0.000000
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.713807 -0.357931 -0.357931 0.002054
- 2 2 0 0 -3.379424 -3.774293 -3.774293 4.169162
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 0.000000 0.000000 0.000000 0.000000
- 2 0 2 0 -5.847335 -2.923667 -2.923667 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.839239 -3.448420 -3.448420 2.057602
- Entering AOResponse driver routine
- --------------------------------------------------------------------------------
- ****************
- *** RESPONSE ***
- ****************
- Response module for NWChem and dynamic CPKS solver
- developed by J. Autschbach and coworkers, SUNY Buffalo
- The methodology used in this program is described in
- J. Chem. Phys. 123 (2005), 114103
- J. Chem. Phys. 122 (2005), 224115
- J. Chem. Phys. 122 (2005), 074105
- Comp. Lett. 3 (2007), 131-150 (contact JA for a copy)
- Please cite this work in publications based on results
- obtained with this code. Thank you!
- For updated response module please acknowledge:
- F. Aquino, Northwestern University, Schatz Rsrch Group
- (version June 9th,2012)
- ++++++++++++++++++++++++++++++++++++++++++++++++++++++
- -----------------------------------------------
- Solving response equations for perturbing field
- -----------------------------------------------
- number of frequencies: 1
- frequency in a.u.: 0.8855900E-01
- Perturbing field: electric
- Using Dipole Length Gauge
- Applying DAMPING constant of 0.007000 a.u.
- Setting up CPKS
- with frequency omega = 0.08855900 a.u.
- and a damping of 0.00700000 a.u.
- NWChem Dynamic CPHF Module
- --------------------------
- scftype = RHF
- nclosed = 5
- nopen = 0
- variables = 1740
- num. vecs = 3
- tolerance = 1.0D-04
- max. iter = 50
- SCF residual: 3.353019351982887E-007
- Iterative solution of linear equations
- No. of variables 870
- No. of equations 3
- Maximum subspace 33
- Convergence 1.0D-04
- Start time 186.6
- iter nsub residual time
- ---- ------ -------- ---------
- dft_CPHF2_read: failed to open/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1
- File /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1 does not exist, proceed to generate (z1,Az1)
- 1 3 2.14D+00 196.6
- 2 6 1.86D+00 205.2
- 3 9 1.59D-01 213.9
- 4 12 1.68D-02 225.5
- 5 15 4.68D-04 234.9
- 6 18 6.28D-05 244.4
- Electric Dipole Response Matrix (nonzero elements):
- 1 2 3
- 1 11.0364 0.0000 0.0000
- 2 0.0000 10.7148 0.0000
- 3 0.0000 0.0000 10.7400
- DFT Linear Response polarizability / au
- Frequency = 0.0885590 / au
- Wavelength = 514.4971488 / nm
- X Y Z
- -----------------------------------------------
- X 11.0363948 0.0000000 0.0000000
- Y 0.0000000 10.7148453 0.0000000
- Z 0.0000000 0.0000000 10.7399839
- -----------------------------------------------
- Eigenvalues = 11.0363948 10.7148453 10.7399839
- Isotropic = 10.8304080
- Anisotropic = 0.3097463
- -----------------------------------------------
- Imaginary part:
- 1 2 3
- 1 0.1215 0.0000 0.0000
- 2 0.0000 0.1910 0.0000
- 3 0.0000 0.0000 0.1408
- average: 10.8304080244247 + I 0.151090073986837
- Magnetic Dipole Response Matrix (nonzero elements):
- Optical rotation tensor Beta
- 1 2 3
- 1 0.0000 0.4213 0.0000
- 2 0.7469 0.0000 0.0000
- 3 0.0000 0.0000 0.0000
- Imaginary part:
- 1 2 3
- 1 0.0000 0.0152 0.0000
- 2 0.0555 0.0000 0.0000
- 3 0.0000 0.0000 0.0000
- average: 5.741073056714595E-012 + I 3.125787453284648E-013
- Exiting AOResponse driver routine
- Stepsize : 1.040753821702078E-002
- 1591.50 REAL 0.834861
- 1591.50 IMAGINARY 0.000549
- 1591.50 TOTAL 0.835409
- Writing mini:/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.raman_0007-0009_nblock
- Writing modecounter= 1
- (ii,rminfo(ii,1))=( 8, 0.00000000)
- ramanpolfile=/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0002-1
- ramanpol_read: failed to open /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0002-1
- NWChem Property Module
- ----------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 5.00D-04
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: fine
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 70 10.0 590
- H 0.35 60 13.0 590
- Grid pruning is: on
- Number of quadrature shells: 190
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Loading old vectors from job with title :
- TDDFT h2o BP86/aug-cc-pVTZ
- Time after variat. SCF: 185.3
- Time prior to 1st pass: 185.3
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 12 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 12.96 12956396
- Stack Space remaining (MW): 13.11 13106484
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4660590342 -8.55D+01 8.74D-05 1.03D-04 186.8
- d= 0,ls=0.0,diis 2 -76.4660630776 -4.04D-06 3.95D-05 4.93D-05 187.9
- d= 0,ls=0.0,diis 3 -76.4660636788 -6.01D-07 1.99D-05 4.30D-05 189.0
- d= 0,ls=0.0,diis 4 -76.4660672312 -3.55D-06 1.07D-06 8.15D-08 190.1
- d= 0,ls=0.0,diis 5 -76.4660672386 -7.39D-09 1.42D-07 1.78D-09 191.2
- Total DFT energy = -76.466067238612
- One electron energy = -122.699712463512
- Coulomb energy = 46.514738396743
- Exchange-Corr. energy = -9.308565634201
- Nuclear repulsion energy = 9.027472462358
- Numeric. integr. density = 10.000000023249
- Total iterative time = 6.0s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878236D+01
- MO Center= 1.5D-09, 2.0D-13, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.973940 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.297150D-01
- MO Center= 1.4D-08, 4.5D-11, -8.3D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.483050 1 O s 4 0.282633 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.824867D-01
- MO Center= -1.9D-08, -2.8D-11, -9.7D-02, r^2= 7.9D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.325655 1 O px 6 0.300689 1 O px
- 48 0.224764 2 H s 71 -0.224764 3 H s
- Vector 4 Occ=2.000000D+00 E=-3.466572D-01
- MO Center= -3.3D-09, -3.8D-10, 1.8D-01, r^2= 7.6D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.359396 1 O pz 4 0.344066 1 O s
- 8 0.319865 1 O pz 14 0.206641 1 O pz
- 2 0.199838 1 O s 48 -0.175242 2 H s
- 71 -0.175242 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.692823D-01
- MO Center= 2.1D-09, -7.1D-11, 7.3D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428318 1 O py 7 0.365110 1 O py
- 13 0.324945 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.282411D-02
- MO Center= -9.5D-08, 2.1D-09, -7.7D-01, r^2= 7.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.472480 1 O s 50 -0.459339 2 H s
- 73 -0.459339 3 H s 4 0.363837 1 O s
- 49 -0.352294 2 H s 72 -0.352294 3 H s
- 14 -0.192354 1 O pz 2 0.151762 1 O s
- Vector 7 Occ=0.000000D+00 E= 1.486863D-02
- MO Center= 1.6D-07, 2.4D-10, -7.7D-01, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.612785 2 H s 73 -2.612784 3 H s
- 49 0.797644 2 H s 72 -0.797644 3 H s
- 15 -0.641757 1 O px 12 -0.340145 1 O px
- 48 0.158073 2 H s 71 -0.158073 3 H s
- Vector 8 Occ=0.000000D+00 E= 6.560367D-02
- MO Center= -3.7D-07, 2.7D-09, -7.6D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.361109 1 O s 49 -2.105246 2 H s
- 72 -2.105246 3 H s 4 1.356311 1 O s
- 17 -1.120185 1 O pz 50 -1.024270 2 H s
- 73 -1.024270 3 H s 57 0.687048 2 H px
- 80 -0.687048 3 H px 59 -0.543301 2 H pz
- Vector 9 Occ=0.000000D+00 E= 9.425703D-02
- MO Center= -1.8D-08, -5.9D-09, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.489326 1 O py 58 -0.305364 2 H py
- 81 -0.305364 3 H py
- Vector 10 Occ=0.000000D+00 E= 9.868887D-02
- MO Center= -2.0D-06, 7.6D-09, 1.6D-01, r^2= 8.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.468496 1 O pz 49 0.602168 2 H s
- 72 0.602160 3 H s 50 -0.352186 2 H s
- 73 -0.352196 3 H s 59 -0.298160 2 H pz
- 82 -0.298163 3 H pz 57 0.273787 2 H px
- 80 -0.273790 3 H px 4 0.245488 1 O s
- Vector 11 Occ=0.000000D+00 E= 9.990159D-02
- MO Center= 2.3D-06, 1.7D-10, -8.2D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.881754 2 H s 73 -2.881753 3 H s
- 49 2.074797 2 H s 72 -2.074798 3 H s
- 15 -1.535126 1 O px 57 -0.916491 2 H px
- 80 -0.916490 3 H px 59 0.901497 2 H pz
- 82 -0.901496 3 H pz 12 -0.533037 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.475832D-01
- MO Center= -4.0D-08, 1.7D-09, 6.4D-01, r^2= 8.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.213223 1 O px 49 -2.851164 2 H s
- 72 2.851164 3 H s 50 -2.161405 2 H s
- 73 2.161405 3 H s 12 1.380977 1 O px
- 48 -0.932200 2 H s 71 0.932200 3 H s
- 31 0.575175 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.677140D-01
- MO Center= 4.2D-08, -1.6D-09, -1.8D-01, r^2= 6.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.441169 1 O s 49 -1.642382 2 H s
- 72 -1.642382 3 H s 59 -0.988851 2 H pz
- 82 -0.988851 3 H pz 57 -0.603966 2 H px
- 80 0.603965 3 H px 17 0.459643 1 O pz
- 4 0.436494 1 O s 30 -0.423487 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.909067D-01
- MO Center= -5.0D-08, -1.8D-09, -5.1D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 58 1.575048 2 H py 81 -1.575048 3 H py
- 28 -0.469429 1 O d -2 55 0.227376 2 H py
- 78 -0.227376 3 H py 65 -0.184073 2 H d -2
- 88 -0.184073 3 H d -2 66 0.166287 2 H d -1
- 89 -0.166287 3 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.136239D-01
- MO Center= 4.3D-08, -6.2D-09, -2.7D-01, r^2= 6.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.980112 1 O py 58 -1.716391 2 H py
- 81 -1.716391 3 H py 13 1.474192 1 O py
- 29 -0.393664 1 O d -1 55 -0.267900 2 H py
- 78 -0.267900 3 H py 65 0.245788 2 H d -2
- 88 -0.245788 3 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.154101D-01
- MO Center= 2.1D-09, 1.5D-09, 2.4D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.595202 1 O s 49 -3.859559 2 H s
- 72 -3.859559 3 H s 17 -2.512168 1 O pz
- 14 -1.990548 1 O pz 4 1.466831 1 O s
- 48 -1.067113 2 H s 71 -1.067113 3 H s
- 59 0.622426 2 H pz 82 0.622425 3 H pz
- Vector 17 Occ=0.000000D+00 E= 3.077187D-01
- MO Center= -2.0D-07, -7.2D-11, -3.1D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 57 3.998121 2 H px 80 3.998123 3 H px
- 15 -3.406526 1 O px 49 -3.095192 2 H s
- 72 3.095196 3 H s 59 1.753364 2 H pz
- 82 -1.753362 3 H pz 12 -1.192378 1 O px
- 50 0.824265 2 H s 73 -0.824265 3 H s
- Vector 18 Occ=0.000000D+00 E= 3.119028D-01
- MO Center= -2.2D-06, 2.8D-10, -8.2D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.655829 1 O s 49 -4.839445 2 H s
- 72 -4.839455 3 H s 4 4.120182 1 O s
- 57 2.083480 2 H px 80 -2.083482 3 H px
- 59 -1.787812 2 H pz 82 -1.787817 3 H pz
- 50 -0.750660 2 H s 73 -0.750664 3 H s
- Vector 19 Occ=0.000000D+00 E= 3.237056D-01
- MO Center= 2.3D-06, 5.8D-10, -9.1D-01, r^2= 7.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 7.704376 2 H s 72 -7.704368 3 H s
- 57 -2.976201 2 H px 80 -2.976198 3 H px
- 59 2.150963 2 H pz 82 -2.150960 3 H pz
- 50 2.005826 2 H s 73 -2.005824 3 H s
- 15 -1.947015 1 O px 12 -1.650182 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.553521D-01
- MO Center= -2.7D-09, 7.2D-10, 2.2D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 9.636232 2 H s 72 -9.636232 3 H s
- 12 -3.599100 1 O px 15 -3.036050 1 O px
- 57 -2.463404 2 H px 80 -2.463404 3 H px
- 31 -1.944850 1 O d 1 48 1.887588 2 H s
- 71 -1.887588 3 H s 50 1.499589 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.946298D-01
- MO Center= -2.7D-08, 2.7D-11, -6.5D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 6.974111 1 O s 49 -4.449492 2 H s
- 72 -4.449492 3 H s 4 1.558069 1 O s
- 32 1.185848 1 O d 2 59 -1.129889 2 H pz
- 82 -1.129889 3 H pz 17 -0.936582 1 O pz
- 30 -0.900654 1 O d 0 14 0.873020 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.450087D-01
- MO Center= -7.9D-09, 6.3D-10, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.178057 1 O s 49 -2.017077 2 H s
- 72 -2.017078 3 H s 14 -1.957663 1 O pz
- 17 -1.895047 1 O pz 30 1.610434 1 O d 0
- 4 1.012696 1 O s 59 1.002180 2 H pz
- 82 1.002179 3 H pz 57 0.620051 2 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10009567
- moments of inertia (a.u.)
- ------------------
- 2.294928319106 0.000000000000 -0.000000111161
- 0.000000000000 6.554085601829 0.000000000000
- -0.000000111161 0.000000000000 4.259157282723
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 0.000000 0.000000 0.000000 0.000000
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.711599 -0.357598 -0.357598 0.003596
- 2 2 0 0 -3.360265 -3.793177 -3.793177 4.226088
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 0.000000 0.000000 0.000000 0.000000
- 2 0 2 0 -5.854192 -2.927096 -2.927096 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.851744 -3.451499 -3.451499 2.051255
- Entering AOResponse driver routine
- --------------------------------------------------------------------------------
- ****************
- *** RESPONSE ***
- ****************
- Response module for NWChem and dynamic CPKS solver
- developed by J. Autschbach and coworkers, SUNY Buffalo
- The methodology used in this program is described in
- J. Chem. Phys. 123 (2005), 114103
- J. Chem. Phys. 122 (2005), 224115
- J. Chem. Phys. 122 (2005), 074105
- Comp. Lett. 3 (2007), 131-150 (contact JA for a copy)
- Please cite this work in publications based on results
- obtained with this code. Thank you!
- For updated response module please acknowledge:
- F. Aquino, Northwestern University, Schatz Rsrch Group
- (version June 9th,2012)
- ++++++++++++++++++++++++++++++++++++++++++++++++++++++
- -----------------------------------------------
- Solving response equations for perturbing field
- -----------------------------------------------
- number of frequencies: 1
- frequency in a.u.: 0.8855900E-01
- Perturbing field: electric
- Using Dipole Length Gauge
- Applying DAMPING constant of 0.007000 a.u.
- Setting up CPKS
- with frequency omega = 0.08855900 a.u.
- and a damping of 0.00700000 a.u.
- NWChem Dynamic CPHF Module
- --------------------------
- scftype = RHF
- nclosed = 5
- nopen = 0
- variables = 1740
- num. vecs = 3
- tolerance = 1.0D-04
- max. iter = 50
- SCF residual: 1.540600463862569E-006
- Iterative solution of linear equations
- No. of variables 870
- No. of equations 3
- Maximum subspace 33
- Convergence 1.0D-04
- Start time 252.7
- iter nsub residual time
- ---- ------ -------- ---------
- dft_CPHF2_read: failed to open/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1
- File /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1 does not exist, proceed to generate (z1,Az1)
- 1 3 2.15D+00 262.5
- 2 6 1.87D+00 274.4
- 3 9 1.60D-01 283.8
- 4 12 1.72D-02 293.7
- 5 15 4.73D-04 304.5
- 6 18 6.32D-05 314.0
- Electric Dipole Response Matrix (nonzero elements):
- 1 2 3
- 1 11.1280 0.0000 0.0000
- 2 0.0000 10.7431 0.0000
- 3 0.0000 0.0000 10.7877
- DFT Linear Response polarizability / au
- Frequency = 0.0885590 / au
- Wavelength = 514.4971488 / nm
- X Y Z
- -----------------------------------------------
- X 11.1279576 0.0000000 0.0000002
- Y 0.0000000 10.7430984 0.0000000
- Z 0.0000002 0.0000000 10.7876906
- -----------------------------------------------
- Eigenvalues = 11.1279576 10.7430984 10.7876906
- Isotropic = 10.8862489
- Anisotropic = 0.3646139
- -----------------------------------------------
- Imaginary part:
- 1 2 3
- 1 0.1229 0.0000 0.0000
- 2 0.0000 0.1919 0.0000
- 3 0.0000 0.0000 0.1420
- average: 10.8862488651764 + I 0.152283452185855
- Magnetic Dipole Response Matrix (nonzero elements):
- Optical rotation tensor Beta
- 1 2 3
- 1 0.0000 0.4273 0.0000
- 2 0.7787 0.0000 0.0000
- 3 0.0000 0.0000 0.0000
- Imaginary part:
- 1 2 3
- 1 0.0000 0.0156 0.0000
- 2 0.0580 0.0000 0.0000
- 3 0.0000 0.0000 0.0000
- average: -9.100574952237800E-011 + I -4.844225207359661E-012
- Exiting AOResponse driver routine
- ramanpolfile=/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0002-2
- ramanpol_read: failed to open /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0002-2
- NWChem Property Module
- ----------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 5.00D-04
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: fine
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 70 10.0 590
- H 0.35 60 13.0 590
- Grid pruning is: on
- Number of quadrature shells: 190
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Loading old vectors from job with title :
- TDDFT h2o BP86/aug-cc-pVTZ
- Time after variat. SCF: 234.2
- Time prior to 1st pass: 234.2
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 12 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 12.96 12956396
- Stack Space remaining (MW): 13.11 13106484
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4660352745 -8.56D+01 1.91D-04 3.88D-04 235.3
- d= 0,ls=0.0,diis 2 -76.4660536580 -1.84D-05 6.97D-05 1.52D-04 236.3
- d= 0,ls=0.0,diis 3 -76.4660545356 -8.78D-07 3.79D-05 1.44D-04 237.4
- d= 0,ls=0.0,diis 4 -76.4660663444 -1.18D-05 3.42D-06 9.69D-07 238.5
- d= 0,ls=0.0,diis 5 -76.4660664315 -8.71D-08 3.12D-07 8.15D-09 239.4
- d= 0,ls=0.0,diis 6 -76.4660664323 -8.02D-10 1.07D-08 6.11D-12 240.4
- Total DFT energy = -76.466066432312
- One electron energy = -122.825891104185
- Coulomb energy = 46.576999976696
- Exchange-Corr. energy = -9.315692437600
- Nuclear repulsion energy = 9.098517132777
- Numeric. integr. density = 9.999999987811
- Total iterative time = 6.3s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878117D+01
- MO Center= -1.7D-09, 2.0D-13, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.974126 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.328116D-01
- MO Center= -1.5D-08, 4.6D-11, -8.5D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.481996 1 O s 4 0.278712 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.850320D-01
- MO Center= 1.9D-08, -2.8D-11, -9.6D-02, r^2= 7.8D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.324585 1 O px 6 0.301349 1 O px
- 48 0.227785 2 H s 71 -0.227785 3 H s
- Vector 4 Occ=2.000000D+00 E=-3.471111D-01
- MO Center= 2.5D-09, -3.8D-10, 1.8D-01, r^2= 7.5D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.359461 1 O pz 4 0.346460 1 O s
- 8 0.320269 1 O pz 14 0.205423 1 O pz
- 2 0.200647 1 O s 48 -0.175627 2 H s
- 71 -0.175627 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.698670D-01
- MO Center= -2.5D-09, -6.5D-11, 7.2D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428130 1 O py 7 0.365006 1 O py
- 13 0.323786 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.192599D-02
- MO Center= 1.0D-07, 2.2D-09, -7.7D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.490310 1 O s 50 -0.466520 2 H s
- 73 -0.466520 3 H s 4 0.368620 1 O s
- 49 -0.355547 2 H s 72 -0.355547 3 H s
- 14 -0.191102 1 O pz 2 0.150132 1 O s
- Vector 7 Occ=0.000000D+00 E= 1.531706D-02
- MO Center= -1.7D-07, 2.4D-10, -7.6D-01, r^2= 1.5D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.653376 2 H s 73 -2.653376 3 H s
- 49 0.812445 2 H s 72 -0.812445 3 H s
- 15 -0.649808 1 O px 12 -0.343894 1 O px
- 48 0.159881 2 H s 71 -0.159881 3 H s
- Vector 8 Occ=0.000000D+00 E= 6.632493D-02
- MO Center= 3.6D-07, 2.8D-09, -5.7D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.420793 1 O s 49 -2.144488 2 H s
- 72 -2.144487 3 H s 4 1.378721 1 O s
- 17 -1.145413 1 O pz 50 -1.020153 2 H s
- 73 -1.020153 3 H s 57 0.684764 2 H px
- 80 -0.684764 3 H px 59 -0.541657 2 H pz
- Vector 9 Occ=0.000000D+00 E= 9.413586D-02
- MO Center= 2.2D-08, -5.4D-09, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.498961 1 O py 58 -0.313581 2 H py
- 81 -0.313581 3 H py
- Vector 10 Occ=0.000000D+00 E= 9.906081D-02
- MO Center= 1.0D-06, 6.9D-09, 1.6D-01, r^2= 8.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.451837 1 O pz 49 0.531337 2 H s
- 72 0.531341 3 H s 50 -0.371409 2 H s
- 73 -0.371404 3 H s 5 0.357192 1 O s
- 59 -0.323540 2 H pz 82 -0.323538 3 H pz
- 4 0.294113 1 O s 57 0.293466 2 H px
- Vector 11 Occ=0.000000D+00 E= 1.017117D-01
- MO Center= -1.3D-06, 2.0D-10, -8.2D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.928863 2 H s 73 -2.928864 3 H s
- 49 2.181498 2 H s 72 -2.181498 3 H s
- 15 -1.597393 1 O px 57 -0.932761 2 H px
- 80 -0.932762 3 H px 59 0.921501 2 H pz
- 82 -0.921501 3 H pz 12 -0.564659 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.479442D-01
- MO Center= 3.1D-08, 1.6D-09, 6.4D-01, r^2= 8.7D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.198339 1 O px 49 -2.814055 2 H s
- 72 2.814055 3 H s 50 -2.135815 2 H s
- 73 2.135815 3 H s 12 1.398406 1 O px
- 48 -0.948500 2 H s 71 0.948500 3 H s
- 31 0.583704 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.683194D-01
- MO Center= -4.0D-08, -1.5D-09, -1.9D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.435579 1 O s 49 -1.652053 2 H s
- 72 -1.652052 3 H s 59 -0.986283 2 H pz
- 82 -0.986283 3 H pz 57 -0.620130 2 H px
- 80 0.620131 3 H px 17 0.456667 1 O pz
- 4 0.448832 1 O s 30 -0.424578 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.912512D-01
- MO Center= 4.8D-08, -1.8D-09, -5.0D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 58 1.589710 2 H py 81 -1.589710 3 H py
- 28 -0.474229 1 O d -2 55 0.233758 2 H py
- 78 -0.233758 3 H py 65 -0.186303 2 H d -2
- 88 -0.186303 3 H d -2 66 0.167836 2 H d -1
- 89 -0.167836 3 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.140946D-01
- MO Center= -4.5D-08, -5.4D-09, -2.7D-01, r^2= 6.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.988322 1 O py 58 -1.728166 2 H py
- 81 -1.728166 3 H py 13 1.503687 1 O py
- 29 -0.397258 1 O d -1 55 -0.275962 2 H py
- 78 -0.275963 3 H py 65 0.248222 2 H d -2
- 88 -0.248222 3 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.160274D-01
- MO Center= -8.9D-09, 7.7D-10, 2.3D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.564013 1 O s 49 -3.847395 2 H s
- 72 -3.847396 3 H s 17 -2.514740 1 O pz
- 14 -2.014303 1 O pz 4 1.481419 1 O s
- 48 -1.089669 2 H s 71 -1.089669 3 H s
- 59 0.646125 2 H pz 82 0.646125 3 H pz
- Vector 17 Occ=0.000000D+00 E= 3.082837D-01
- MO Center= 1.6D-07, -6.9D-11, -3.1D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 57 4.062859 2 H px 80 4.062858 3 H px
- 15 -3.436123 1 O px 49 -3.223227 2 H s
- 72 3.223223 3 H s 59 1.774024 2 H pz
- 82 -1.774026 3 H pz 12 -1.204967 1 O px
- 50 0.847849 2 H s 73 -0.847849 3 H s
- Vector 18 Occ=0.000000D+00 E= 3.136665D-01
- MO Center= 2.2D-06, 2.8D-10, -8.2D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.783157 1 O s 49 -4.927783 2 H s
- 72 -4.927773 3 H s 4 4.156118 1 O s
- 57 2.089325 2 H px 80 -2.089323 3 H px
- 59 -1.794476 2 H pz 82 -1.794471 3 H pz
- 50 -0.753514 2 H s 73 -0.753510 3 H s
- Vector 19 Occ=0.000000D+00 E= 3.259142D-01
- MO Center= -2.3D-06, 5.8D-10, -9.0D-01, r^2= 7.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 7.929873 2 H s 72 -7.929881 3 H s
- 57 -3.033032 2 H px 80 -3.033035 3 H px
- 59 2.164834 2 H pz 82 -2.164837 3 H pz
- 50 2.024267 2 H s 73 -2.024268 3 H s
- 15 -1.968237 1 O px 12 -1.666110 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.551469D-01
- MO Center= 9.5D-09, 7.1D-10, 2.3D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 9.704769 2 H s 72 -9.704769 3 H s
- 12 -3.657291 1 O px 15 -3.025144 1 O px
- 57 -2.459435 2 H px 80 -2.459435 3 H px
- 31 -1.977066 1 O d 1 48 1.955303 2 H s
- 71 -1.955303 3 H s 50 1.499477 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.954442D-01
- MO Center= 2.3D-08, 3.0D-11, -6.4D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 6.917532 1 O s 49 -4.407506 2 H s
- 72 -4.407506 3 H s 4 1.589217 1 O s
- 32 1.204386 1 O d 2 59 -1.112248 2 H pz
- 82 -1.112248 3 H pz 17 -0.923616 1 O pz
- 30 -0.896764 1 O d 0 14 0.862375 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.451385D-01
- MO Center= 1.1D-08, 6.4D-10, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.088843 1 O s 14 -1.952695 1 O pz
- 49 -1.957809 2 H s 72 -1.957808 3 H s
- 17 -1.886780 1 O pz 30 1.609057 1 O d 0
- 59 1.021875 2 H pz 82 1.021875 3 H pz
- 4 0.961424 1 O s 57 0.601513 2 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10009567
- moments of inertia (a.u.)
- ------------------
- 2.259902915757 0.000000000000 0.000000111949
- 0.000000000000 6.452092120954 0.000000000000
- 0.000000111949 0.000000000000 4.192189205197
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 0.000000 0.000000 0.000000 0.000000
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.710774 -0.361005 -0.361005 0.011236
- 2 2 0 0 -3.370678 -3.765159 -3.765159 4.159640
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 0.000000 0.000000 0.000000 0.000000
- 2 0 2 0 -5.840782 -2.920391 -2.920391 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.848494 -3.434227 -3.434227 2.019960
- Entering AOResponse driver routine
- --------------------------------------------------------------------------------
- ****************
- *** RESPONSE ***
- ****************
- Response module for NWChem and dynamic CPKS solver
- developed by J. Autschbach and coworkers, SUNY Buffalo
- The methodology used in this program is described in
- J. Chem. Phys. 123 (2005), 114103
- J. Chem. Phys. 122 (2005), 224115
- J. Chem. Phys. 122 (2005), 074105
- Comp. Lett. 3 (2007), 131-150 (contact JA for a copy)
- Please cite this work in publications based on results
- obtained with this code. Thank you!
- For updated response module please acknowledge:
- F. Aquino, Northwestern University, Schatz Rsrch Group
- (version June 9th,2012)
- ++++++++++++++++++++++++++++++++++++++++++++++++++++++
- -----------------------------------------------
- Solving response equations for perturbing field
- -----------------------------------------------
- number of frequencies: 1
- frequency in a.u.: 0.8855900E-01
- Perturbing field: electric
- Using Dipole Length Gauge
- Applying DAMPING constant of 0.007000 a.u.
- Setting up CPKS
- with frequency omega = 0.08855900 a.u.
- and a damping of 0.00700000 a.u.
- NWChem Dynamic CPHF Module
- --------------------------
- scftype = RHF
- nclosed = 5
- nopen = 0
- variables = 1740
- num. vecs = 3
- tolerance = 1.0D-04
- max. iter = 50
- SCF residual: 1.848828471317094E-007
- Iterative solution of linear equations
- No. of variables 870
- No. of equations 3
- Maximum subspace 33
- Convergence 1.0D-04
- Start time 323.1
- iter nsub residual time
- ---- ------ -------- ---------
- dft_CPHF2_read: failed to open/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1
- File /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1 does not exist, proceed to generate (z1,Az1)
- 1 3 2.10D+00 333.4
- 2 6 1.84D+00 343.3
- 3 9 1.58D-01 353.4
- 4 12 1.59D-02 363.1
- 5 15 4.57D-04 372.9
- 6 18 6.20D-05 382.3
- Electric Dipole Response Matrix (nonzero elements):
- 1 2 3
- 1 10.9610 0.0000 0.0000
- 2 0.0000 10.6871 0.0000
- 3 0.0000 0.0000 10.6787
- DFT Linear Response polarizability / au
- Frequency = 0.0885590 / au
- Wavelength = 514.4971488 / nm
- X Y Z
- -----------------------------------------------
- X 10.9609968 0.0000000 -0.0000002
- Y 0.0000000 10.6870849 0.0000000
- Z -0.0000002 0.0000000 10.6786594
- -----------------------------------------------
- Eigenvalues = 10.9609968 10.6870849 10.6786594
- Isotropic = 10.7755804
- Anisotropic = 0.2782204
- -----------------------------------------------
- Imaginary part:
- 1 2 3
- 1 0.1203 0.0000 0.0000
- 2 0.0000 0.1901 0.0000
- 3 0.0000 0.0000 0.1400
- average: 10.7755803675108 + I 0.150109917531141
- Magnetic Dipole Response Matrix (nonzero elements):
- Optical rotation tensor Beta
- 1 2 3
- 1 0.0000 0.4292 0.0000
- 2 0.7050 0.0000 0.0000
- 3 0.0000 0.0000 0.0000
- Imaginary part:
- 1 2 3
- 1 0.0000 0.0153 0.0000
- 2 0.0527 0.0000 0.0000
- 3 0.0000 0.0000 0.0000
- average: -8.989191951370877E-011 + I -4.826458870561861E-012
- Exiting AOResponse driver routine
- Stepsize : 1.022109686160916E-002
- 3680.21 REAL 115.555937
- 3680.21 IMAGINARY 0.040682
- 3680.21 TOTAL 115.596619
- Writing mini:/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.raman_0007-0009_nblock
- Writing modecounter= 2
- (ii,rminfo(ii,1))=( 9, 0.00000000)
- ramanpolfile=/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0003-1
- ramanpol_read: failed to open /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0003-1
- NWChem Property Module
- ----------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 5.00D-04
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: fine
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 70 10.0 590
- H 0.35 60 13.0 590
- Grid pruning is: on
- Number of quadrature shells: 190
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Loading old vectors from job with title :
- TDDFT h2o BP86/aug-cc-pVTZ
- Time after variat. SCF: 286.1
- Time prior to 1st pass: 286.1
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 12 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 12.96 12956396
- Stack Space remaining (MW): 13.11 13106484
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4660494347 -8.55D+01 1.38D-04 2.00D-04 287.2
- d= 0,ls=0.0,diis 2 -76.4660616355 -1.22D-05 4.09D-05 4.35D-05 288.3
- d= 0,ls=0.0,diis 3 -76.4660608151 8.20D-07 2.40D-05 5.41D-05 289.4
- d= 0,ls=0.0,diis 4 -76.4660652308 -4.42D-06 4.04D-06 1.05D-06 290.4
- d= 0,ls=0.0,diis 5 -76.4660653284 -9.76D-08 4.45D-07 9.05D-09 291.5
- d= 0,ls=0.0,diis 6 -76.4660653294 -9.98D-10 6.07D-08 2.78D-10 292.6
- Total DFT energy = -76.466065329398
- One electron energy = -122.762769199443
- Coulomb energy = 46.545837194126
- Exchange-Corr. energy = -9.312125549415
- Nuclear repulsion energy = 9.062992225334
- Numeric. integr. density = 10.000000004434
- Total iterative time = 6.5s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878176D+01
- MO Center= -3.6D-04, -2.5D-17, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.974034 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.312679D-01
- MO Center= -3.1D-03, -2.8D-15, -8.4D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.482516 1 O s 4 0.280645 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.837544D-01
- MO Center= 4.2D-03, 2.1D-15, -9.6D-02, r^2= 7.8D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.325125 1 O px 6 0.301016 1 O px
- 71 -0.227636 3 H s 48 0.224884 2 H s
- Vector 4 Occ=2.000000D+00 E=-3.468832D-01
- MO Center= 6.4D-04, -1.5D-14, 1.8D-01, r^2= 7.5D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.359427 1 O pz 4 0.345274 1 O s
- 8 0.320067 1 O pz 14 0.206029 1 O pz
- 2 0.200243 1 O s 48 -0.175069 2 H s
- 71 -0.175809 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.695741D-01
- MO Center= -5.1D-04, 9.2D-15, 7.3D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428224 1 O py 7 0.365058 1 O py
- 13 0.324364 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.237868D-02
- MO Center= 2.2D-02, -2.0D-14, -7.7D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.481225 1 O s 50 -0.470821 2 H s
- 73 -0.454915 3 H s 4 0.366168 1 O s
- 49 -0.356203 2 H s 72 -0.351537 3 H s
- 14 -0.191714 1 O pz 2 0.150950 1 O s
- Vector 7 Occ=0.000000D+00 E= 1.509528D-02
- MO Center= -3.1D-02, 8.7D-14, -7.7D-01, r^2= 1.5D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 73 -2.641516 3 H s 50 2.624495 2 H s
- 72 -0.813631 3 H s 49 0.796412 2 H s
- 15 -0.645763 1 O px 12 -0.341999 1 O px
- 71 -0.160679 3 H s 48 0.157262 2 H s
- Vector 8 Occ=0.000000D+00 E= 6.595381D-02
- MO Center= 6.7D-02, -7.9D-14, -6.7D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.389549 1 O s 49 -2.131265 2 H s
- 72 -2.116936 3 H s 4 1.367328 1 O s
- 17 -1.132181 1 O pz 50 -1.072448 2 H s
- 73 -0.971856 3 H s 57 0.696342 2 H px
- 80 -0.675390 3 H px 59 -0.554273 2 H pz
- Vector 9 Occ=0.000000D+00 E= 9.419621D-02
- MO Center= 4.3D-03, 1.4D-13, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.494148 1 O py 81 -0.314656 3 H py
- 58 -0.304297 2 H py
- Vector 10 Occ=0.000000D+00 E= 9.877945D-02
- MO Center= 2.3D-01, -4.0D-13, 1.2D-01, r^2= 8.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.432846 1 O pz 72 0.995204 3 H s
- 50 -0.951918 2 H s 59 -0.483048 2 H pz
- 57 0.461263 2 H px 15 0.319580 1 O px
- 4 0.253737 1 O s 73 0.253986 3 H s
- 5 0.236774 1 O s 14 -0.221702 1 O pz
- Vector 11 Occ=0.000000D+00 E= 1.009127D-01
- MO Center= -2.9D-01, 5.9D-14, -7.8D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 73 -2.938221 3 H s 50 2.744448 2 H s
- 49 2.133492 2 H s 72 -2.036688 3 H s
- 15 -1.532943 1 O px 80 -0.980172 3 H px
- 82 -0.974367 3 H pz 57 -0.830662 2 H px
- 59 0.811510 2 H pz 12 -0.537905 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.477604D-01
- MO Center= 1.0D-02, -1.1D-13, 6.4D-01, r^2= 8.7D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.206024 1 O px 49 -2.863025 2 H s
- 72 2.802935 3 H s 50 -2.155694 2 H s
- 73 2.142271 3 H s 12 1.389647 1 O px
- 48 -0.939937 2 H s 71 0.940655 3 H s
- 31 0.579391 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.680166D-01
- MO Center= 1.9D-03, 1.6D-13, -1.8D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.437602 1 O s 49 -1.646870 2 H s
- 72 -1.646849 3 H s 59 -0.995789 2 H pz
- 82 -0.979315 3 H pz 80 0.619077 3 H px
- 57 -0.605104 2 H px 17 0.458309 1 O pz
- 4 0.442455 1 O s 30 -0.424054 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.910774D-01
- MO Center= 1.2D-02, -2.6D-14, -5.0D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 58 1.595641 2 H py 81 -1.568719 3 H py
- 28 -0.471762 1 O d -2 55 0.230477 2 H py
- 78 -0.230524 3 H py 65 -0.187091 2 H d -2
- 88 -0.183219 3 H d -2 66 0.167611 2 H d -1
- 89 -0.166468 3 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.138609D-01
- MO Center= -1.1D-02, 5.3D-13, -2.7D-01, r^2= 6.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.984149 1 O py 81 -1.747648 3 H py
- 58 -1.696852 2 H py 13 1.489020 1 O py
- 29 -0.395459 1 O d -1 78 -0.278853 3 H py
- 55 -0.265064 2 H py 88 -0.250696 3 H d -2
- 65 0.243309 2 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.157199D-01
- MO Center= -2.3D-03, -2.7D-13, 2.3D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.579667 1 O s 72 -3.865783 3 H s
- 49 -3.841316 2 H s 17 -2.513410 1 O pz
- 14 -2.002500 1 O pz 4 1.474339 1 O s
- 71 -1.103911 3 H s 48 -1.053110 2 H s
- 82 0.644671 3 H pz 59 0.623855 2 H pz
- Vector 17 Occ=0.000000D+00 E= 3.079962D-01
- MO Center= 1.2D-02, 3.9D-14, -3.1D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 57 4.096837 2 H px 80 3.971761 3 H px
- 15 -3.415075 1 O px 49 -3.345978 2 H s
- 72 2.998685 3 H s 82 -1.833991 3 H pz
- 59 1.681982 2 H pz 12 -1.194089 1 O px
- 73 -0.858070 3 H s 50 0.803987 2 H s
- Vector 18 Occ=0.000000D+00 E= 3.126891D-01
- MO Center= 2.6D-01, 1.4D-15, -8.2D-01, r^2= 7.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.669144 1 O s 49 -5.394235 2 H s
- 72 -4.314792 3 H s 4 4.115097 1 O s
- 57 2.197742 2 H px 59 -2.031853 2 H pz
- 80 -1.955615 3 H px 82 -1.531851 3 H pz
- 50 -0.957026 2 H s 73 -0.539514 3 H s
- Vector 19 Occ=0.000000D+00 E= 3.249113D-01
- MO Center= -2.7D-01, 3.0D-15, -9.0D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 72 8.302074 3 H s 49 -7.282222 2 H s
- 80 3.188676 3 H px 57 2.804418 2 H px
- 82 2.314665 3 H pz 73 2.077064 3 H s
- 59 -1.981323 2 H pz 15 1.947210 1 O px
- 50 -1.935335 2 H s 12 1.649465 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.552491D-01
- MO Center= 5.1D-03, -3.9D-14, 2.2D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 9.702745 2 H s 72 -9.637402 3 H s
- 12 -3.627990 1 O px 15 -3.030618 1 O px
- 57 -2.478905 2 H px 80 -2.443793 3 H px
- 31 -1.960837 1 O d 1 71 -1.942526 3 H s
- 48 1.899948 2 H s 50 1.501568 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.950400D-01
- MO Center= 3.4D-03, -4.7D-15, -6.4D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 6.945142 1 O s 49 -4.447359 2 H s
- 72 -4.408747 3 H s 4 1.573023 1 O s
- 32 1.195122 1 O d 2 59 -1.131741 2 H pz
- 82 -1.110318 3 H pz 17 -0.929974 1 O pz
- 30 -0.898788 1 O d 0 14 0.867968 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.450733D-01
- MO Center= 1.9D-03, -4.9D-15, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.133752 1 O s 49 -1.992292 2 H s
- 72 -1.982968 3 H s 14 -1.955454 1 O pz
- 17 -1.890957 1 O pz 30 1.609761 1 O d 0
- 82 1.019640 3 H pz 59 1.004467 2 H pz
- 4 0.987604 1 O s 57 0.616660 2 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10009567
- moments of inertia (a.u.)
- ------------------
- 2.277417037350 0.000000000000 0.024329657605
- 0.000000000000 6.503089094835 0.000000000000
- 0.024329657605 0.000000000000 4.225672057484
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 0.002230 -0.001538 -0.001538 0.005307
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.711194 -0.359305 -0.359305 0.007416
- 2 2 0 0 -3.365449 -3.779154 -3.779154 4.192859
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 -0.007358 0.007791 0.007791 -0.022940
- 2 0 2 0 -5.847488 -2.923744 -2.923744 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.850120 -3.442865 -3.442865 2.035611
- Entering AOResponse driver routine
- --------------------------------------------------------------------------------
- ****************
- *** RESPONSE ***
- ****************
- Response module for NWChem and dynamic CPKS solver
- developed by J. Autschbach and coworkers, SUNY Buffalo
- The methodology used in this program is described in
- J. Chem. Phys. 123 (2005), 114103
- J. Chem. Phys. 122 (2005), 224115
- J. Chem. Phys. 122 (2005), 074105
- Comp. Lett. 3 (2007), 131-150 (contact JA for a copy)
- Please cite this work in publications based on results
- obtained with this code. Thank you!
- For updated response module please acknowledge:
- F. Aquino, Northwestern University, Schatz Rsrch Group
- (version June 9th,2012)
- ++++++++++++++++++++++++++++++++++++++++++++++++++++++
- -----------------------------------------------
- Solving response equations for perturbing field
- -----------------------------------------------
- number of frequencies: 1
- frequency in a.u.: 0.8855900E-01
- Perturbing field: electric
- Using Dipole Length Gauge
- Applying DAMPING constant of 0.007000 a.u.
- Setting up CPKS
- with frequency omega = 0.08855900 a.u.
- and a damping of 0.00700000 a.u.
- NWChem Dynamic CPHF Module
- --------------------------
- scftype = RHF
- nclosed = 5
- nopen = 0
- variables = 1740
- num. vecs = 3
- tolerance = 1.0D-04
- max. iter = 50
- SCF residual: 9.594869511617160E-007
- Iterative solution of linear equations
- No. of variables 870
- No. of equations 3
- Maximum subspace 33
- Convergence 1.0D-04
- Start time 391.0
- iter nsub residual time
- ---- ------ -------- ---------
- dft_CPHF2_read: failed to open/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1
- File /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1 does not exist, proceed to generate (z1,Az1)
- 1 3 2.11D+00 400.4
- 2 6 1.85D+00 409.4
- 3 9 1.59D-01 418.8
- 4 12 1.66D-02 430.4
- 5 15 4.65D-04 439.7
- 6 18 6.26D-05 449.1
- Electric Dipole Response Matrix (nonzero elements):
- 1 2 3
- 1 11.0443 0.0000 -0.0432
- 2 0.0000 10.7151 0.0000
- 3 -0.0432 0.0000 10.7332
- DFT Linear Response polarizability / au
- Frequency = 0.0885590 / au
- Wavelength = 514.4971488 / nm
- X Y Z
- -----------------------------------------------
- X 11.0442906 0.0000000 -0.0431938
- Y 0.0000000 10.7150696 0.0000000
- Z -0.0431938 0.0000000 10.7332067
- -----------------------------------------------
- Eigenvalues = 11.0501767 10.7150696 10.7273207
- Isotropic = 10.8308557
- Anisotropic = 0.3291525
- -----------------------------------------------
- Imaginary part:
- 1 2 3
- 1 0.1216 0.0000 -0.0008
- 2 0.0000 0.1910 0.0000
- 3 -0.0008 0.0000 0.1410
- average: 10.8308556511775 + I 0.151194815164448
- Magnetic Dipole Response Matrix (nonzero elements):
- Optical rotation tensor Beta
- 1 2 3
- 1 0.0000 0.4284 0.0000
- 2 0.7417 0.0000 0.0397
- 3 0.0000 -0.0335 0.0000
- Imaginary part:
- 1 2 3
- 1 0.0000 0.0155 0.0000
- 2 0.0554 0.0000 0.0027
- 3 0.0000 -0.0015 0.0000
- average: -3.938026519322422E-016 + I -3.075505758841691E-016
- Exiting AOResponse driver routine
- ramanpolfile=/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0003-2
- ramanpol_read: failed to open /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.ramanpol_0003-2
- NWChem Property Module
- ----------------------
- TDDFT h2o BP86/aug-cc-pVTZ
- NWChem DFT Module
- -----------------
- TDDFT h2o BP86/aug-cc-pVTZ
- Caching 1-el integrals
- General Information
- -------------------
- SCF calculation type: DFT
- Wavefunction type: closed shell.
- No. of atoms : 3
- No. of electrons : 10
- Alpha electrons : 5
- Beta electrons : 5
- Charge : 0
- Spin multiplicity: 1
- Use of symmetry is: off; symmetry adaption is: off
- Maximum number of iterations: 30
- This is a Direct SCF calculation.
- AO basis - number of functions: 92
- number of shells: 32
- Convergence on energy requested: 1.00D-08
- Convergence on density requested: 1.00D-05
- Convergence on gradient requested: 5.00D-04
- XC Information
- --------------
- Becke 1988 Exchange Functional 1.000
- Perdew 1981 Correlation Functional 1.000 local
- Perdew 1986 Correlation Functional 1.000 non-local
- Grid Information
- ----------------
- Grid used for XC integration: fine
- Radial quadrature: Mura-Knowles
- Angular quadrature: Lebedev.
- Tag B.-S. Rad. Rad. Pts. Rad. Cut. Ang. Pts.
- --- ---------- --------- --------- ---------
- O 0.60 70 10.0 590
- H 0.35 60 13.0 590
- Grid pruning is: on
- Number of quadrature shells: 190
- Spatial weights used: Erf1
- Convergence Information
- -----------------------
- Convergence aids based upon iterative change in
- total energy or number of iterations.
- Levelshifting, if invoked, occurs when the
- HOMO/LUMO gap drops below (HL_TOL): 1.00D-02
- DIIS, if invoked, will attempt to extrapolate
- using up to (NFOCK): 10 stored Fock matrices.
- Damping( 0%) Levelshifting(0.0) DIIS
- --------------- ------------------- ---------------
- dE on: start N/A start
- dE off: 2 iters N/A 30 iters
- Screening Tolerance Information
- -------------------------------
- Density screening/tol_rho: 1.00D-11
- AO Gaussian exp screening on grid/accAOfunc: 18
- CD Gaussian exp screening on grid/accCDfunc: 20
- XC Gaussian exp screening on grid/accXCfunc: 20
- Schwarz screening/accCoul: 1.00D-09
- Loading old vectors from job with title :
- TDDFT h2o BP86/aug-cc-pVTZ
- Time after variat. SCF: 337.5
- Time prior to 1st pass: 337.5
- Grid_pts file = /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.gridpts.0
- Record size in doubles = 12289 No. of grid_pts per rec = 3070
- Max. records in memory = 12 Max. recs in file = 7245
- Memory utilization after 1st SCF pass:
- Heap Space remaining (MW): 12.96 12956396
- Stack Space remaining (MW): 13.11 13106484
- convergence iter energy DeltaE RMS-Dens Diis-err time
- ---------------- ----- ----------------- --------- --------- --------- ------
- d= 0,ls=0.0,diis 1 -76.4660331410 -8.55D+01 1.99D-04 4.02D-04 338.7
- d= 0,ls=0.0,diis 2 -76.4660623426 -2.92D-05 3.36D-05 3.11D-05 339.6
- d= 0,ls=0.0,diis 3 -76.4660647504 -2.41D-06 1.15D-05 5.71D-06 340.6
- d= 0,ls=0.0,diis 4 -76.4660653115 -5.61D-07 2.32D-06 2.02D-07 341.7
- d= 0,ls=0.0,diis 5 -76.4660653193 -7.76D-09 1.05D-06 1.20D-07 342.8
- Total DFT energy = -76.466065319301
- One electron energy = -122.762035161584
- Coulomb energy = 46.545003085844
- Exchange-Corr. energy = -9.312025793021
- Nuclear repulsion energy = 9.062992549459
- Numeric. integr. density = 10.000000004457
- Total iterative time = 5.2s
- DFT Final Molecular Orbital Analysis
- ------------------------------------
- Vector 1 Occ=2.000000D+00 E=-1.878176D+01
- MO Center= 3.6D-04, 7.5D-19, 1.2D-01, r^2= 1.5D-02
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 1 0.974034 1 O s
- Vector 2 Occ=2.000000D+00 E=-9.312680D-01
- MO Center= 3.1D-03, 1.0D-15, -8.4D-02, r^2= 5.0D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 2 0.482516 1 O s 4 0.280645 1 O s
- Vector 3 Occ=2.000000D+00 E=-4.837546D-01
- MO Center= -4.2D-03, 8.0D-17, -9.6D-02, r^2= 7.8D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 9 0.325125 1 O px 6 0.301016 1 O px
- 48 0.227635 2 H s 71 -0.224884 3 H s
- Vector 4 Occ=2.000000D+00 E=-3.468833D-01
- MO Center= -6.4D-04, 8.6D-16, 1.8D-01, r^2= 7.5D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 11 0.359427 1 O pz 4 0.345274 1 O s
- 8 0.320067 1 O pz 14 0.206029 1 O pz
- 2 0.200243 1 O s 48 -0.175809 2 H s
- 71 -0.175069 3 H s
- Vector 5 Occ=2.000000D+00 E=-2.695742D-01
- MO Center= 5.1D-04, -1.5D-15, 7.3D-02, r^2= 7.2D-01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 10 0.428224 1 O py 7 0.365058 1 O py
- 13 0.324364 1 O py
- Vector 6 Occ=0.000000D+00 E=-3.237871D-02
- MO Center= -2.2D-02, -5.7D-15, -7.7D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 0.481224 1 O s 73 -0.470822 3 H s
- 50 -0.454913 2 H s 4 0.366168 1 O s
- 72 -0.356204 3 H s 49 -0.351537 2 H s
- 14 -0.191714 1 O pz 2 0.150950 1 O s
- Vector 7 Occ=0.000000D+00 E= 1.509525D-02
- MO Center= 3.1D-02, 5.3D-16, -7.7D-01, r^2= 1.5D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.641515 2 H s 73 -2.624494 3 H s
- 49 0.813630 2 H s 72 -0.796411 3 H s
- 15 -0.645763 1 O px 12 -0.341998 1 O px
- 48 0.160679 2 H s 71 -0.157262 3 H s
- Vector 8 Occ=0.000000D+00 E= 6.595379D-02
- MO Center= -6.7D-02, -1.2D-14, -6.7D-02, r^2= 1.4D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 5.389547 1 O s 72 -2.131265 3 H s
- 49 -2.116934 2 H s 4 1.367328 1 O s
- 17 -1.132180 1 O pz 73 -1.072449 3 H s
- 50 -0.971855 2 H s 80 -0.696343 3 H px
- 57 0.675389 2 H px 82 -0.554273 3 H pz
- Vector 9 Occ=0.000000D+00 E= 9.419621D-02
- MO Center= -4.3D-03, 5.0D-15, 3.5D-01, r^2= 7.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.494148 1 O py 58 -0.314658 2 H py
- 81 -0.304295 3 H py
- Vector 10 Occ=0.000000D+00 E= 9.877943D-02
- MO Center= -2.3D-01, 1.2D-14, 1.2D-01, r^2= 8.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 17 1.432839 1 O pz 49 0.995255 2 H s
- 73 -0.951985 3 H s 82 -0.483066 3 H pz
- 80 -0.461282 3 H px 15 -0.319619 1 O px
- 4 0.253733 1 O s 50 0.254059 2 H s
- 5 0.236765 1 O s 14 -0.221700 1 O pz
- Vector 11 Occ=0.000000D+00 E= 1.009127D-01
- MO Center= 2.9D-01, 9.5D-16, -7.8D-01, r^2= 1.1D+01
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 50 2.938216 2 H s 73 -2.744425 3 H s
- 72 -2.133494 3 H s 49 2.036665 2 H s
- 15 -1.532934 1 O px 57 -0.980175 2 H px
- 59 0.974371 2 H pz 80 -0.830650 3 H px
- 82 -0.811497 3 H pz 12 -0.537902 1 O px
- Vector 12 Occ=0.000000D+00 E= 1.477604D-01
- MO Center= -1.0D-02, 5.9D-15, 6.4D-01, r^2= 8.7D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 15 3.206024 1 O px 72 2.863030 3 H s
- 49 -2.802932 2 H s 73 2.155695 3 H s
- 50 -2.142271 2 H s 12 1.389647 1 O px
- 48 -0.940653 2 H s 71 0.939938 3 H s
- 31 0.579391 1 O d 1
- Vector 13 Occ=0.000000D+00 E= 1.680166D-01
- MO Center= -1.9D-03, -5.6D-15, -1.8D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.437603 1 O s 49 -1.646848 2 H s
- 72 -1.646871 3 H s 82 -0.995790 3 H pz
- 59 -0.979314 2 H pz 57 -0.619076 2 H px
- 80 0.605104 3 H px 17 0.458309 1 O pz
- 4 0.442455 1 O s 30 -0.424053 1 O d 0
- Vector 14 Occ=0.000000D+00 E= 1.910774D-01
- MO Center= -1.2D-02, -7.8D-15, -5.0D-01, r^2= 5.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 81 -1.595649 3 H py 58 1.568711 2 H py
- 28 -0.471762 1 O d -2 55 0.230523 2 H py
- 78 -0.230478 3 H py 88 -0.187092 3 H d -2
- 65 -0.183217 2 H d -2 89 -0.167612 3 H d -1
- 66 0.166467 2 H d -1
- Vector 15 Occ=0.000000D+00 E= 2.138609D-01
- MO Center= 1.1D-02, 1.5D-14, -2.7D-01, r^2= 6.4D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 16 1.984149 1 O py 58 -1.747655 2 H py
- 81 -1.696846 3 H py 13 1.489020 1 O py
- 29 -0.395459 1 O d -1 55 -0.278854 2 H py
- 78 -0.265063 3 H py 65 0.250697 2 H d -2
- 88 -0.243308 3 H d -2
- Vector 16 Occ=0.000000D+00 E= 2.157199D-01
- MO Center= 2.3D-03, -8.6D-15, 2.3D-01, r^2= 6.8D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 7.579668 1 O s 49 -3.865791 2 H s
- 72 -3.841308 3 H s 17 -2.513410 1 O pz
- 14 -2.002500 1 O pz 4 1.474339 1 O s
- 48 -1.103911 2 H s 71 -1.053110 3 H s
- 59 0.644671 2 H pz 82 0.623856 3 H pz
- Vector 17 Occ=0.000000D+00 E= 3.079962D-01
- MO Center= -1.2D-02, -4.1D-16, -3.1D-01, r^2= 8.0D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 80 4.096846 3 H px 57 3.971753 2 H px
- 15 -3.415074 1 O px 72 3.346003 3 H s
- 49 -2.998668 2 H s 59 1.833997 2 H pz
- 82 -1.681973 3 H pz 12 -1.194088 1 O px
- 50 0.858073 2 H s 73 -0.803983 3 H s
- Vector 18 Occ=0.000000D+00 E= 3.126890D-01
- MO Center= -2.6D-01, -1.4D-15, -8.2D-01, r^2= 7.9D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 8.669136 1 O s 72 -5.394267 3 H s
- 49 -4.314750 2 H s 4 4.115093 1 O s
- 80 -2.197744 3 H px 82 -2.031873 3 H pz
- 57 1.955611 2 H px 59 -1.531827 2 H pz
- 73 -0.957043 3 H s 50 -0.539497 2 H s
- Vector 19 Occ=0.000000D+00 E= 3.249113D-01
- MO Center= 2.7D-01, -3.3D-16, -9.0D-01, r^2= 7.3D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 49 8.302105 2 H s 72 -7.282182 3 H s
- 57 -3.188688 2 H px 80 -2.804403 3 H px
- 59 2.314676 2 H pz 50 2.077068 2 H s
- 82 -1.981309 3 H pz 15 -1.947209 1 O px
- 73 -1.935329 3 H s 12 -1.649464 1 O px
- Vector 20 Occ=0.000000D+00 E= 4.552491D-01
- MO Center= -5.1D-03, 1.9D-15, 2.2D-01, r^2= 6.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 72 -9.702756 3 H s 49 9.637394 2 H s
- 12 -3.627990 1 O px 15 -3.030618 1 O px
- 80 -2.478907 3 H px 57 -2.443793 2 H px
- 31 -1.960837 1 O d 1 48 1.942525 2 H s
- 71 -1.899948 3 H s 50 1.497516 2 H s
- Vector 21 Occ=0.000000D+00 E= 4.950400D-01
- MO Center= -3.4D-03, -3.8D-15, -6.4D-01, r^2= 4.1D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 6.945142 1 O s 72 -4.447353 3 H s
- 49 -4.408752 2 H s 4 1.573023 1 O s
- 32 1.195122 1 O d 2 82 -1.131742 3 H pz
- 59 -1.110317 2 H pz 17 -0.929974 1 O pz
- 30 -0.898788 1 O d 0 14 0.867968 1 O pz
- Vector 22 Occ=0.000000D+00 E= 5.450733D-01
- MO Center= -1.9D-03, 2.1D-15, 2.9D-01, r^2= 4.2D+00
- Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
- ----- ------------ --------------- ----- ------------ ---------------
- 5 3.133753 1 O s 49 -1.982974 2 H s
- 72 -1.992287 3 H s 14 -1.955454 1 O pz
- 17 -1.890957 1 O pz 30 1.609761 1 O d 0
- 59 1.019639 2 H pz 82 1.004468 3 H pz
- 4 0.987604 1 O s 80 -0.616659 3 H px
- center of mass
- --------------
- x = 0.00000000 y = 0.00000000 z = 0.10009567
- moments of inertia (a.u.)
- ------------------
- 2.277416903664 0.000000000000 -0.024329655108
- 0.000000000000 6.503088627949 0.000000000000
- -0.024329655108 0.000000000000 4.225671724284
- Multipole analysis of the density
- ---------------------------------
- L x y z total alpha beta nuclear
- - - - - ----- ----- ---- -------
- 0 0 0 0 0.000000 -5.000000 -5.000000 10.000000
- 1 1 0 0 -0.002228 0.001540 0.001540 -0.005307
- 1 0 1 0 0.000000 0.000000 0.000000 0.000000
- 1 0 0 1 -0.711194 -0.359305 -0.359305 0.007416
- 2 2 0 0 -3.365448 -3.779154 -3.779154 4.192859
- 2 1 1 0 0.000000 0.000000 0.000000 0.000000
- 2 1 0 1 0.007357 -0.007791 -0.007791 0.022940
- 2 0 2 0 -5.847486 -2.923743 -2.923743 0.000000
- 2 0 1 1 0.000000 0.000000 0.000000 0.000000
- 2 0 0 2 -4.850118 -3.442865 -3.442865 2.035611
- Entering AOResponse driver routine
- --------------------------------------------------------------------------------
- ****************
- *** RESPONSE ***
- ****************
- Response module for NWChem and dynamic CPKS solver
- developed by J. Autschbach and coworkers, SUNY Buffalo
- The methodology used in this program is described in
- J. Chem. Phys. 123 (2005), 114103
- J. Chem. Phys. 122 (2005), 224115
- J. Chem. Phys. 122 (2005), 074105
- Comp. Lett. 3 (2007), 131-150 (contact JA for a copy)
- Please cite this work in publications based on results
- obtained with this code. Thank you!
- For updated response module please acknowledge:
- F. Aquino, Northwestern University, Schatz Rsrch Group
- (version June 9th,2012)
- ++++++++++++++++++++++++++++++++++++++++++++++++++++++
- -----------------------------------------------
- Solving response equations for perturbing field
- -----------------------------------------------
- number of frequencies: 1
- frequency in a.u.: 0.8855900E-01
- Perturbing field: electric
- Using Dipole Length Gauge
- Applying DAMPING constant of 0.007000 a.u.
- Setting up CPKS
- with frequency omega = 0.08855900 a.u.
- and a damping of 0.00700000 a.u.
- NWChem Dynamic CPHF Module
- --------------------------
- scftype = RHF
- nclosed = 5
- nopen = 0
- variables = 1740
- num. vecs = 3
- tolerance = 1.0D-04
- max. iter = 50
- SCF residual: 2.908535687994365E-006
- Iterative solution of linear equations
- No. of variables 870
- No. of equations 3
- Maximum subspace 33
- Convergence 1.0D-04
- Start time 456.9
- iter nsub residual time
- ---- ------ -------- ---------
- dft_CPHF2_read: failed to open/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1
- File /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.aoresp_fiao_f1 does not exist, proceed to generate (z1,Az1)
- 1 3 2.11D+00 466.8
- 2 6 1.85D+00 477.2
- 3 9 1.59D-01 487.8
- 4 12 1.66D-02 497.8
- 5 15 4.65D-04 507.5
- 6 18 6.26D-05 516.8
- Electric Dipole Response Matrix (nonzero elements):
- 1 2 3
- 1 11.0443 0.0000 0.0432
- 2 0.0000 10.7151 0.0000
- 3 0.0432 0.0000 10.7332
- DFT Linear Response polarizability / au
- Frequency = 0.0885590 / au
- Wavelength = 514.4971488 / nm
- X Y Z
- -----------------------------------------------
- X 11.0442896 0.0000000 0.0431929
- Y 0.0000000 10.7150689 0.0000000
- Z 0.0431930 0.0000000 10.7332059
- -----------------------------------------------
- Eigenvalues = 11.0501754 10.7150689 10.7273200
- Isotropic = 10.8308548
- Anisotropic = 0.3291520
- -----------------------------------------------
- Imaginary part:
- 1 2 3
- 1 0.1216 0.0000 0.0008
- 2 0.0000 0.1910 0.0000
- 3 0.0008 0.0000 0.1410
- average: 10.8308547911949 + I 0.151194883504067
- Magnetic Dipole Response Matrix (nonzero elements):
- Optical rotation tensor Beta
- 1 2 3
- 1 0.0000 0.4284 0.0000
- 2 0.7417 0.0000 -0.0397
- 3 0.0000 0.0336 0.0000
- Imaginary part:
- 1 2 3
- 1 0.0000 0.0155 0.0000
- 2 0.0554 0.0000 -0.0027
- 3 0.0000 0.0015 0.0000
- average: 7.929949419144013E-016 + I -4.577021045239255E-017
- Exiting AOResponse driver routine
- Stepsize : 1.039961178096392E-002
- 3783.86 REAL 28.406861
- 3783.86 IMAGINARY 0.010009
- 3783.86 TOTAL 28.416870
- Writing mini:/home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.raman_0007-0009_nblock
- Writing modecounter= 3
- Showing rminfo:
- (start,last,rmmodes,nfreq,plot,width)=( 1, 9, 3, 1,normal , 20.00000000)
- rminfo( 1)=( 1591.50355320, 0.83486088, 0.00054855, 0.83540943)
- rminfo( 2)=( 3680.21084179, 115.55593697, 0.04068209, 115.59661905)
- rminfo( 3)=( 3783.85869559, 28.40686131, 0.01000883, 28.41687015)
- Raman scattering data written to /home/d3y133/nwchem-dev/nwchem-dev-2013-06-17/QA/scratchdir/h2o_raman_3.normal
- Task times cpu: 360.4s wall: 482.8s
- Summary of allocated global arrays
- -----------------------------------
- No active global arrays
- GA Statistics for process 0
- ------------------------------
- create destroy get put acc scatter gather read&inc
- calls: 6797 6797 3.52e+06 1.60e+06 8.42e+05 4430 0 3766
- number of processes/call 1.15e+00 1.09e+00 1.19e+00 2.64e+00 0.00e+00
- bytes total: 4.98e+09 1.52e+08 2.10e+09 1.52e+07 0.00e+00 3.01e+04
- bytes remote: 4.06e+09 5.11e+07 1.77e+09 -3.30e+07 0.00e+00 0.00e+00
- Max memory consumed for GA by this process: 1740160 bytes
- MA_summarize_allocated_blocks: starting scan ...
- MA_summarize_allocated_blocks: scan completed: 0 heap blocks, 0 stack blocks
- MA usage statistics:
- allocation statistics:
- heap stack
- ---- -----
- current number of blocks 0 0
- maximum number of blocks 24 52
- current total bytes 0 0
- maximum total bytes 1205184 113917288
- maximum total K-bytes 1206 113918
- maximum total M-bytes 2 114
- NWChem Input Module
- -------------------
- CITATION
- --------
- Please cite the following reference when publishing
- results obtained with NWChem:
- M. Valiev, E.J. Bylaska, N. Govind, K. Kowalski,
- T.P. Straatsma, H.J.J. van Dam, D. Wang, J. Nieplocha,
- E. Apra, T.L. Windus, W.A. de Jong
- "NWChem: a comprehensive and scalable open-source
- solution for large scale molecular simulations"
- Comput. Phys. Commun. 181, 1477 (2010)
- doi:10.1016/j.cpc.2010.04.018
- AUTHORS & CONTRIBUTORS
- ----------------------
- E. Apra, E. J. Bylaska, W. A. de Jong, N. Govind, K. Kowalski,
- T. P. Straatsma, M. Valiev, H. J. J. van Dam, D. Wang, T. L. Windus,
- J. Hammond, J. Autschbach, K. Bhaskaran-Nair, J. Brabec, K. Lopata,
- F. Aquino, S. Hirata, M. T. Hackler, J. Mullin, P. Nichols, R. Peverati,
- J. Pittner, Y. Zhao, P.-D. Fan, R. J. Harrison, M. Dupuis, D. Silverstein,
- D. M. A. Smith, J. Nieplocha, V. Tipparaju, M. Krishnan, B. E. Van Kuiken,
- A. Vazquez-Mayagoitia, L. Jensen, M. Swart, Q. Wu, T. Van Voorhis,
- A. A. Auer, M. Nooijen, L. D. Crosby, E. Brown, G. Cisneros, G. I. Fann,
- H. Fruchtl, J. Garza, K. Hirao, R. Kendall, J. A. Nichols, K. Tsemekhman,
- K. Wolinski, J. Anchell, D. Bernholdt, P. Borowski, T. Clark, D. Clerc,
- H. Dachsel, M. Deegan, K. Dyall, D. Elwood, E. Glendening, M. Gutowski,
- A. Hess, J. Jaffe, B. Johnson, J. Ju, R. Kobayashi, R. Kutteh, Z. Lin,
- R. Littlefield, X. Long, B. Meng, T. Nakajima, S. Niu, L. Pollack, M. Rosing,
- K. Glaesemann, G. Sandrone, M. Stave, H. Taylor, G. Thomas, J. H. van Lenthe,
- A. Wong, Z. Zhang.
- Total times cpu: 386.2s wall: 517.0s
Advertisement
Add Comment
Please, Sign In to add comment