Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- *****************
- * O R C A *
- *****************
- #,
- ###
- ####
- #####
- ######
- ########,
- ,,################,,,,,
- ,,#################################,,
- ,,##########################################,,
- ,#########################################, ''#####,
- ,#############################################,, '####,
- ,##################################################,,,,####,
- ,###########'''' ''''###############################
- ,#####'' ,,,,##########,,,, '''####''' '####
- ,##' ,,,,###########################,,, '##
- ' ,,###'''' '''############,,,
- ,,##'' '''############,,,, ,,,,,,###''
- ,#'' '''#######################'''
- ' ''''####''''
- ,#######, #######, ,#######, ##
- ,#' '#, ## ## ,#' '#, #''# ###### ,####,
- ## ## ## ,#' ## #' '# # #' '#
- ## ## ####### ## ,######, #####, # #
- '#, ,#' ## ## '#, ,#' ,# #, ## #, ,#
- '#######' ## ## '#######' #' '# #####' # '####'
- #######################################################
- # -***- #
- # Department of theory and spectroscopy #
- # Directorship and core code : Frank Neese #
- # Max Planck Institute fuer Kohlenforschung #
- # Kaiser Wilhelm Platz 1 #
- # D-45470 Muelheim/Ruhr #
- # Germany #
- # #
- # All rights reserved #
- # -***- #
- #######################################################
- Program Version 5.0.3 - RELEASE -
- With contributions from (in alphabetic order):
- Daniel Aravena : Magnetic Suceptibility
- Michael Atanasov : Ab Initio Ligand Field Theory (pilot matlab implementation)
- Alexander A. Auer : GIAO ZORA, VPT2 properties, NMR spectrum
- Ute Becker : Parallelization
- Giovanni Bistoni : ED, misc. LED, open-shell LED, HFLD
- Martin Brehm : Molecular dynamics
- Dmytro Bykov : SCF Hessian
- Vijay G. Chilkuri : MRCI spin determinant printing, contributions to CSF-ICE
- Dipayan Datta : RHF DLPNO-CCSD density
- Achintya Kumar Dutta : EOM-CC, STEOM-CC
- Dmitry Ganyushin : Spin-Orbit,Spin-Spin,Magnetic field MRCI
- Miquel Garcia : C-PCM and meta-GGA Hessian, CC/C-PCM, Gaussian charge scheme
- Yang Guo : DLPNO-NEVPT2, F12-NEVPT2, CIM, IAO-localization
- Andreas Hansen : Spin unrestricted coupled pair/coupled cluster methods
- Benjamin Helmich-Paris : MC-RPA, TRAH-SCF, COSX integrals
- Lee Huntington : MR-EOM, pCC
- Robert Izsak : Overlap fitted RIJCOSX, COSX-SCS-MP3, EOM
- Marcus Kettner : VPT2
- Christian Kollmar : KDIIS, OOCD, Brueckner-CCSD(T), CCSD density, CASPT2, CASPT2-K
- Simone Kossmann : Meta GGA functionals, TD-DFT gradient, OOMP2, MP2 Hessian
- Martin Krupicka : Initial AUTO-CI
- Lucas Lang : DCDCAS
- Marvin Lechner : AUTO-CI (C++ implementation), FIC-MRCC
- Dagmar Lenk : GEPOL surface, SMD
- Dimitrios Liakos : Extrapolation schemes; Compound Job, initial MDCI parallelization
- Dimitrios Manganas : Further ROCIS development; embedding schemes
- Dimitrios Pantazis : SARC Basis sets
- Anastasios Papadopoulos: AUTO-CI, single reference methods and gradients
- Taras Petrenko : DFT Hessian,TD-DFT gradient, ASA, ECA, R-Raman, ABS, FL, XAS/XES, NRVS
- Peter Pinski : DLPNO-MP2, DLPNO-MP2 Gradient
- Christoph Reimann : Effective Core Potentials
- Marius Retegan : Local ZFS, SOC
- Christoph Riplinger : Optimizer, TS searches, QM/MM, DLPNO-CCSD(T), (RO)-DLPNO pert. Triples
- Tobias Risthaus : Range-separated hybrids, TD-DFT gradient, RPA, STAB
- Michael Roemelt : Original ROCIS implementation
- Masaaki Saitow : Open-shell DLPNO-CCSD energy and density
- Barbara Sandhoefer : DKH picture change effects
- Avijit Sen : IP-ROCIS
- Kantharuban Sivalingam : CASSCF convergence, NEVPT2, FIC-MRCI
- Bernardo de Souza : ESD, SOC TD-DFT
- Georgi Stoychev : AutoAux, RI-MP2 NMR, DLPNO-MP2 response
- Willem Van den Heuvel : Paramagnetic NMR
- Boris Wezisla : Elementary symmetry handling
- Frank Wennmohs : Technical directorship
- We gratefully acknowledge several colleagues who have allowed us to
- interface, adapt or use parts of their codes:
- Stefan Grimme, W. Hujo, H. Kruse, P. Pracht, : VdW corrections, initial TS optimization,
- C. Bannwarth, S. Ehlert DFT functionals, gCP, sTDA/sTD-DF
- Ed Valeev, F. Pavosevic, A. Kumar : LibInt (2-el integral package), F12 methods
- Garnet Chan, S. Sharma, J. Yang, R. Olivares : DMRG
- Ulf Ekstrom : XCFun DFT Library
- Mihaly Kallay : mrcc (arbitrary order and MRCC methods)
- Jiri Pittner, Ondrej Demel : Mk-CCSD
- Frank Weinhold : gennbo (NPA and NBO analysis)
- Christopher J. Cramer and Donald G. Truhlar : smd solvation model
- Lars Goerigk : TD-DFT with DH, B97 family of functionals
- V. Asgeirsson, H. Jonsson : NEB implementation
- FAccTs GmbH : IRC, NEB, NEB-TS, DLPNO-Multilevel, CI-OPT
- MM, QMMM, 2- and 3-layer-ONIOM, Crystal-QMMM,
- LR-CPCM, SF, NACMEs, symmetry and pop. for TD-DFT,
- nearIR, NL-DFT gradient (VV10), updates on ESD,
- ML-optimized integration grids
- S Lehtola, MJT Oliveira, MAL Marques : LibXC Library
- Liviu Ungur et al : ANISO software
- Your calculation uses the libint2 library for the computation of 2-el integrals
- For citations please refer to: http://libint.valeyev.net
- Your ORCA version has been built with support for libXC version: 5.1.0
- For citations please refer to: https://tddft.org/programs/libxc/
- This ORCA versions uses:
- CBLAS interface : Fast vector & matrix operations
- LAPACKE interface : Fast linear algebra routines
- SCALAPACK package : Parallel linear algebra routines
- Shared memory : Shared parallel matrices
- BLAS/LAPACK : OpenBLAS 0.3.19 NO_AFFINITY VORTEX SINGLE_THREADED
- Core in use : VORTEX
- Copyright (c) 2011-2014, The OpenBLAS Project
- XCFun DFT library Copyright 2009-2010 Ulf Ekstrom and contributors.
- See http://admol.org/xcfun for more information.
- This is free software; see the source code for copying conditions.
- There is ABSOLUTELY NO WARRANTY; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
- For details see the documentation.
- Scientific users of this library should cite U. Ekstrom, L. Visscher, R. Bast, A. J. Thorvaldsen and K. Ruud;
- J.Chem.Theor.Comp. 2010, DOI: 10.1021/ct100117s
- XCFun Version 0.99
- ================================================================================
- ----- Orbital basis set information -----
- Your calculation utilizes the basis: STO-3G
- H-Ne : W. J. Hehre, R. F. Stewart and J. A. Pople, J. Chem. Phys. 2657 (1969).
- Na-Ar : W. J. Hehre, R. Ditchfield, R. F. Stewart and J. A. Pople, J. Chem. Phys. 2769 (1970).
- K,Ca,Ga-Kr : W. J. Pietro, B. A. Levy, W. J. Hehre and R. F. Stewart, J. Am. Chem. Soc. 19, 2225 (1980).
- Sc-Zn,Y-Cd : W. J. Pietro and W. J. Hehre, J. Comp. Chem. 4, 241 (1983).
- ----- AuxJ basis set information -----
- Your calculation utilizes the auxiliary basis: def2/J
- F. Weigend, Phys. Chem. Chem. Phys. 8, 1057 (2006).
- ================================================================================
- WARNINGS
- Please study these warnings very carefully!
- ================================================================================
- WARNING: Minnesota functionals are quite sensitive to the integration grid.
- see SE Wheeler, KN Houk, JCTC 2010, 6, 395,
- N Mardirossian, M Head-Gordon, JCTC 2016, 12, 4303.
- DEFGRID3 seems to be a minimum grid for reliable results with these functionals!
- ===> : Please increase the integration grid!
- INFO : the flag for use of the SHARK integral package has been found!
- ================================================================================
- INPUT FILE
- ================================================================================
- NAME = en.inp
- | 1> ! m062x STO-3G PAL8 PrintMOs PrintBasis Freq
- | 2>
- | 3> * xyz 0 1
- | 4> F -2.92273910742335 1.36469000000000 0.00000000000000
- | 5> H -2.00034089257665 1.36469000000000 0.00000000000000
- | 6> *
- | 7>
- | 8> ****END OF INPUT****
- ================================================================================
- ****************************
- * Single Point Calculation *
- ****************************
- ---------------------------------
- CARTESIAN COORDINATES (ANGSTROEM)
- ---------------------------------
- F -2.922739 1.364690 0.000000
- H -2.000341 1.364690 0.000000
- ----------------------------
- CARTESIAN COORDINATES (A.U.)
- ----------------------------
- NO LB ZA FRAG MASS X Y Z
- 0 F 9.0000 0 18.998 -5.523176 2.578890 0.000000
- 1 H 1.0000 0 1.008 -3.780096 2.578890 0.000000
- --------------------------------
- INTERNAL COORDINATES (ANGSTROEM)
- --------------------------------
- F 0 0 0 0.000000000000 0.00000000 0.00000000
- H 1 0 0 0.922398214847 0.00000000 0.00000000
- ---------------------------
- INTERNAL COORDINATES (A.U.)
- ---------------------------
- F 0 0 0 0.000000000000 0.00000000 0.00000000
- H 1 0 0 1.743080012478 0.00000000 0.00000000
- ---------------------
- BASIS SET INFORMATION
- ---------------------
- There are 2 groups of distinct atoms
- Group 1 Type F : 6s3p contracted to 2s1p pattern {33/3}
- Group 2 Type H : 3s contracted to 1s pattern {3}
- Atom 0F basis set group => 1
- Atom 1H basis set group => 2
- -------------------------
- BASIS SET IN INPUT FORMAT
- -------------------------
- # Basis set for element : H
- NewGTO H
- S 3
- 1 3.4252509100 0.1543289707
- 2 0.6239137300 0.5353281424
- 3 0.1688554000 0.4446345420
- end;
- # Basis set for element : F
- NewGTO F
- S 3
- 1 166.6791300000 0.1543289701
- 2 30.3608120000 0.5353281404
- 3 8.2168207000 0.4446345403
- S 3
- 1 6.4648032000 -0.0999672291
- 2 1.5022812000 0.3995128265
- 3 0.4885885000 0.7001154638
- P 3
- 1 6.4648032000 0.1559162698
- 2 1.5022812000 0.6076837191
- 3 0.4885885000 0.3919573894
- end;
- ---------------------------------
- AUXILIARY/J BASIS SET INFORMATION
- ---------------------------------
- There are 2 groups of distinct atoms
- Group 1 Type F : 12s5p4d2f1g contracted to 6s4p3d1f1g pattern {711111/2111/211/2/1}
- Group 2 Type H : 5s2p1d contracted to 3s1p1d pattern {311/2/1}
- Atom 0F basis set group => 1
- Atom 1H basis set group => 2
- -------------------------------------
- AUXILIARY/J BASIS SET IN INPUT FORMAT
- -------------------------------------
- # Auxiliary/J basis set for element : H
- NewAuxJGTO H
- S 3
- 1 15.6752927000 0.1007184237
- 2 3.6063578000 0.3404257499
- 3 1.2080016000 0.6491996172
- S 1
- 1 0.4726794000 1.0000000000
- S 1
- 1 0.2018100000 1.0000000000
- P 2
- 1 2.0281365000 0.5596620919
- 2 0.5358730000 0.5596620919
- D 1
- 1 2.2165124000 1.0000000000
- end;
- # Auxiliary/J basis set for element : F
- NewAuxJGTO F
- S 7
- 1 3514.9201549000 0.0176189594
- 2 1231.2458815000 0.0358663858
- 3 454.1367336000 0.1239867649
- 4 175.8376258000 0.2758620517
- 5 71.2081775000 0.4252764982
- 6 30.0316883000 0.2651963387
- 7 13.1259023000 -0.0078514564
- S 1
- 1 5.9126861000 -1.0000000000
- S 1
- 1 2.7284431000 1.0000000000
- S 1
- 1 1.2813993000 1.0000000000
- S 1
- 1 0.6082793000 1.0000000000
- S 1
- 1 0.2897894000 1.0000000000
- P 2
- 1 49.3740409000 -0.8404235528
- 2 13.4874956000 -0.2331250081
- P 1
- 1 3.9664900000 1.0000000000
- P 1
- 1 1.2226629000 -1.0000000000
- P 1
- 1 0.3830317000 1.0000000000
- D 2
- 1 29.9719140000 -0.2160273306
- 2 6.7860915000 -0.8915140390
- D 1
- 1 1.8137313000 -1.0000000000
- D 1
- 1 0.5157840000 -1.0000000000
- F 2
- 1 2.8301877000 0.5762078303
- 2 0.9155962000 0.5762078303
- G 1
- 1 1.6097500000 -1.0000000000
- end;
- ************************************************************
- * Program running with 8 parallel MPI-processes *
- * working on a common directory *
- ************************************************************
- ------------------------------------------------------------------------------
- ORCA GTO INTEGRAL CALCULATION
- -- RI-GTO INTEGRALS CHOSEN --
- ------------------------------------------------------------------------------
- ------------------------------------------------------------------------------
- ___
- / \ - P O W E R E D B Y -
- / \
- | | | _ _ __ _____ __ __
- | | | | | | | / \ | _ \ | | / |
- \ \/ | | | | / \ | | | | | | / /
- / \ \ | |__| | / /\ \ | |_| | | |/ /
- | | | | __ | / /__\ \ | / | \
- | | | | | | | | __ | | \ | |\ \
- \ / | | | | | | | | | |\ \ | | \ \
- \___/ |_| |_| |__| |__| |_| \__\ |__| \__/
- - O R C A' S B I G F R I E N D -
- &
- - I N T E G R A L F E E D E R -
- v1 FN, 2020, v2 2021
- ------------------------------------------------------------------------------
- Reading SHARK input file en.SHARKINP.tmp ... ok
- ----------------------
- SHARK INTEGRAL PACKAGE
- ----------------------
- Number of atoms ... 2
- Number of basis functions ... 6
- Number of shells ... 4
- Maximum angular momentum ... 1
- Integral batch strategy ... SHARK/LIBINT Hybrid
- RI-J (if used) integral strategy ... SPLIT-RIJ (Revised 2003 algorithm where possible)
- Printlevel ... 1
- Contraction scheme used ... SEGMENTED contraction
- Coulomb Range Separation ... NOT USED
- Exchange Range Separation ... NOT USED
- Finite Nucleus Model ... NOT USED
- Auxiliary Coulomb fitting basis ... AVAILABLE
- # of basis functions in Aux-J ... 60
- # of shells in Aux-J ... 20
- Maximum angular momentum in Aux-J ... 4
- Auxiliary J/K fitting basis ... NOT available
- Auxiliary Correlation fitting basis ... NOT available
- Auxiliary 'external' fitting basis ... NOT available
- Integral threshold ... 1.000000e-10
- Primitive cut-off ... 1.000000e-11
- Primitive pair pre-selection threshold ... 1.000000e-11
- Calculating pre-screening integrals ... done ( 0.0 sec) Dimension = 4
- Organizing shell pair data ... done ( 0.0 sec)
- Shell pair information
- Total number of shell pairs ... 10
- Shell pairs after pre-screening ... 10
- Total number of primitive shell pairs ... 90
- Primitive shell pairs kept ... 90
- la=0 lb=0: 6 shell pairs
- la=1 lb=0: 3 shell pairs
- la=1 lb=1: 1 shell pairs
- Calculating one electron integrals ... done ( 0.0 sec)
- Calculating RI/J V-Matrix + Cholesky decomp.... done ( 0.0 sec)
- Calculating Nuclear repulsion ... done ( 0.0 sec) ENN= 5.163274167320 Eh
- SHARK setup successfully completed in 0.0 seconds
- Maximum memory used throughout the entire GTOINT-calculation: 6.2 MB
- ************************************************************
- * Program running with 8 parallel MPI-processes *
- * working on a common directory *
- ************************************************************
- -------------------------------------------------------------------------------
- ORCA SCF
- -------------------------------------------------------------------------------
- ------------
- SCF SETTINGS
- ------------
- Hamiltonian:
- Density Functional Method .... DFT(GTOs)
- Exchange Functional Exchange .... M062X
- Correlation Functional Correlation .... M062X
- Gradients option PostSCFGGA .... off
- Hybrid DFT is turned on
- Fraction HF Exchange ScalHFX .... 0.540000
- Scaling of DF-GGA-X ScalDFX .... 0.460000
- Scaling of DF-GGA-C ScalDFC .... 1.000000
- Scaling of DF-LDA-C ScalLDAC .... 1.000000
- Perturbative correction .... 0.000000
- Density functional embedding theory .... OFF
- RI-approximation to the Coulomb term is turned on
- Number of AuxJ basis functions .... 60
- RIJ-COSX (HFX calculated with COS-X)).... on
- General Settings:
- Integral files IntName .... en
- Hartree-Fock type HFTyp .... RHF
- Total Charge Charge .... 0
- Multiplicity Mult .... 1
- Number of Electrons NEL .... 10
- Basis Dimension Dim .... 6
- Nuclear Repulsion ENuc .... 5.1632741673 Eh
- Convergence Acceleration:
- DIIS CNVDIIS .... on
- Start iteration DIISMaxIt .... 12
- Startup error DIISStart .... 0.200000
- # of expansion vecs DIISMaxEq .... 5
- Bias factor DIISBfac .... 1.050
- Max. coefficient DIISMaxC .... 10.000
- Trust-Rad. Augm. Hess. CNVTRAH .... auto
- Auto Start mean grad. ratio tolernc. .... 1.125000
- Auto Start start iteration .... 20
- Auto Start num. interpolation iter. .... 10
- Max. Number of Micro iterations .... 16
- Max. Number of Macro iterations .... Maxiter - #DIIS iter
- Number of Davidson start vectors .... 2
- Converg. threshold I (grad. norm) .... 5.000e-05
- Converg. threshold II (energy diff.) .... 1.000e-06
- Grad. Scal. Fac. for Micro threshold .... 0.100
- Minimum threshold for Micro iter. .... 0.010
- NR start threshold (gradient norm) .... 0.001
- Initial trust radius .... 0.400
- Minimum AH scaling param. (alpha) .... 1.000
- Maximum AH scaling param. (alpha) .... 1000.000
- Orbital update algorithm .... Taylor
- White noise on init. David. guess .... on
- Maximum white noise .... 0.010
- Quad. conv. algorithm .... NR
- SOSCF CNVSOSCF .... on
- Start iteration SOSCFMaxIt .... 150
- Startup grad/error SOSCFStart .... 0.003300
- Level Shifting CNVShift .... on
- Level shift para. LevelShift .... 0.2500
- Turn off err/grad. ShiftErr .... 0.0010
- Zerner damping CNVZerner .... off
- Static damping CNVDamp .... on
- Fraction old density DampFac .... 0.7000
- Max. Damping (<1) DampMax .... 0.9800
- Min. Damping (>=0) DampMin .... 0.0000
- Turn off err/grad. DampErr .... 0.1000
- Fernandez-Rico CNVRico .... off
- SCF Procedure:
- Maximum # iterations MaxIter .... 125
- SCF integral mode SCFMode .... Direct
- Integral package .... SHARK and LIBINT hybrid scheme
- Reset frequency DirectResetFreq .... 20
- Integral Threshold Thresh .... 1.000e-10 Eh
- Primitive CutOff TCut .... 1.000e-11 Eh
- Convergence Tolerance:
- Convergence Check Mode ConvCheckMode .... Total+1el-Energy
- Convergence forced ConvForced .... 0
- Energy Change TolE .... 1.000e-06 Eh
- 1-El. energy change .... 1.000e-03 Eh
- Orbital Gradient TolG .... 5.000e-05
- Orbital Rotation angle TolX .... 5.000e-05
- DIIS Error TolErr .... 1.000e-06
- Diagonalization of the overlap matrix:
- Smallest eigenvalue ... 4.168e-01
- Time for diagonalization ... 0.001 sec
- Threshold for overlap eigenvalues ... 1.000e-08
- Number of eigenvalues below threshold ... 0
- Time for construction of square roots ... 0.001 sec
- Total time needed ... 0.004 sec
- Time for model grid setup = 0.008 sec
- ------------------------------
- INITIAL GUESS: MODEL POTENTIAL
- ------------------------------
- Loading Hartree-Fock densities ... done
- Calculating cut-offs ... done
- Initializing the effective Hamiltonian ... done
- Setting up the integral package (SHARK) ... done
- Starting the Coulomb interaction ... done ( 0.0 sec)
- Reading the grid ... done
- Mapping shells ... done
- Starting the XC term evaluation ... done ( 0.0 sec)
- promolecular density results
- # of electrons = 9.999296584
- EX = -10.276158856
- EC = -0.330468699
- EX+EC = -10.606627555
- Transforming the Hamiltonian ... done ( 0.0 sec)
- Diagonalizing the Hamiltonian ... done ( 0.0 sec)
- Back transforming the eigenvectors ... done ( 0.0 sec)
- Now organizing SCF variables ... done
- ------------------
- INITIAL GUESS DONE ( 0.0 sec)
- ------------------
- -------------------
- DFT GRID GENERATION
- -------------------
- General Integration Accuracy IntAcc ... 4.388
- Radial Grid Type RadialGrid ... OptM3 with GC (2021)
- Angular Grid (max. ang.) AngularGrid ... 4 (Lebedev-302)
- Angular grid pruning method GridPruning ... 4 (adaptive)
- Weight generation scheme WeightScheme... Becke
- Basis function cutoff BFCut ... 1.0000e-10
- Integration weight cutoff WCut ... 1.0000e-14
- Angular grids for H and He will be reduced by one unit
- Partially contracted basis set ... off
- Rotationally invariant grid construction ... off
- Total number of grid points ... 9454
- Total number of batches ... 149
- Average number of points per batch ... 63
- Average number of grid points per atom ... 4727
- Time for grid setup = 0.031 sec
- --------------------
- COSX GRID GENERATION
- --------------------
- GRIDX 1
- -------
- General Integration Accuracy IntAcc ... 3.816
- Radial Grid Type RadialGrid ... OptM3 with GC (2021)
- Angular Grid (max. ang.) AngularGrid ... 1 (Lebedev-50)
- Angular grid pruning method GridPruning ... 4 (adaptive)
- Weight generation scheme WeightScheme... Becke
- Basis function cutoff BFCut ... 1.0000e-10
- Integration weight cutoff WCut ... 1.0000e-14
- Angular grids for H and He will be reduced by one unit
- Partially contracted basis set ... on
- Rotationally invariant grid construction ... off
- Total number of grid points ... 1210
- Total number of batches ... 20
- Average number of points per batch ... 60
- Average number of grid points per atom ... 605
- UseSFitting ... on
- GRIDX 2
- -------
- General Integration Accuracy IntAcc ... 4.020
- Radial Grid Type RadialGrid ... OptM3 with GC (2021)
- Angular Grid (max. ang.) AngularGrid ... 2 (Lebedev-110)
- Angular grid pruning method GridPruning ... 4 (adaptive)
- Weight generation scheme WeightScheme... Becke
- Basis function cutoff BFCut ... 1.0000e-10
- Integration weight cutoff WCut ... 1.0000e-14
- Angular grids for H and He will be reduced by one unit
- Partially contracted basis set ... on
- Rotationally invariant grid construction ... off
- Total number of grid points ... 2693
- Total number of batches ... 43
- Average number of points per batch ... 62
- Average number of grid points per atom ... 1346
- UseSFitting ... on
- GRIDX 3
- -------
- General Integration Accuracy IntAcc ... 4.338
- Radial Grid Type RadialGrid ... OptM3 with GC (2021)
- Angular Grid (max. ang.) AngularGrid ... 3 (Lebedev-194)
- Angular grid pruning method GridPruning ... 4 (adaptive)
- Weight generation scheme WeightScheme... Becke
- Basis function cutoff BFCut ... 1.0000e-10
- Integration weight cutoff WCut ... 1.0000e-14
- Angular grids for H and He will be reduced by one unit
- Partially contracted basis set ... on
- Rotationally invariant grid construction ... off
- Total number of grid points ... 6003
- Total number of batches ... 95
- Average number of points per batch ... 63
- Average number of grid points per atom ... 3002
- UseSFitting ... on
- Time for X-Grid setup = 0.022 sec
- --------------
- SCF ITERATIONS
- --------------
- ITER Energy Delta-E Max-DP RMS-DP [F,P] Damp
- *** Starting incremental Fock matrix formation ***
- 0 -98.9088085073 0.000000000000 0.07050674 0.01640221 0.1953843 0.7000
- 1 -98.9171080261 -0.008299518831 0.06089516 0.01431967 0.1389558 0.7000
- ***Turning on DIIS***
- 2 -98.9226417645 -0.005533738415 0.14024293 0.03309291 0.0902280 0.0000
- 3 -98.9335415036 -0.010899739110 0.02082397 0.00466988 0.0228956 0.0000
- 4 -98.9338134775 -0.000271973818 0.00411635 0.00094162 0.0044111 0.0000
- *** Initiating the SOSCF procedure ***
- *** Shutting down DIIS ***
- *** Re-Reading the Fockian ***
- *** Removing any level shift ***
- ITER Energy Delta-E Grad Rot Max-DP RMS-DP
- 5 -98.93382218 -0.0000087001 0.000387 0.000387 0.001039 0.000248
- *** Restarting incremental Fock matrix formation ***
- 6 -98.93372388 0.0000982962 0.000071 0.000053 0.000154 0.000035
- **** Energy Check signals convergence ****
- ***Rediagonalizing the Fockian in SOSCF/NRSCF***
- *****************************************************
- * SUCCESS *
- * SCF CONVERGED AFTER 7 CYCLES *
- *****************************************************
- Old exchange energy = -5.745884802 Eh
- New exchange energy = -5.746089493 Eh
- Exchange energy change after final integration = -0.000204691 Eh
- Total energy after final integration = -98.933928582 Eh
- Final COS-X integration done in = 0.011 sec
- ----------------
- TOTAL SCF ENERGY
- ----------------
- Total Energy : -98.93392858 Eh -2692.12906 eV
- Components:
- Nuclear Repulsion : 5.16327417 Eh 140.49983 eV
- Electronic Energy : -104.09720275 Eh -2832.62889 eV
- One Electron Energy: -149.67869597 Eh -4072.96438 eV
- Two Electron Energy: 45.58149322 Eh 1240.33549 eV
- Max COSX asymmetry : 0.00000174 Eh 0.00005 eV
- Virial components:
- Potential Energy : -196.69719938 Eh -5352.40291 eV
- Kinetic Energy : 97.76327080 Eh 2660.27384 eV
- Virial Ratio : 2.01197441
- DFT components:
- N(Alpha) : 4.999999355858 electrons
- N(Beta) : 4.999999355858 electrons
- N(Total) : 9.999998711716 electrons
- E(XC) : -5.256638287626 Eh
- DFET-embed. en. : 0.000000000000 Eh
- ---------------
- SCF CONVERGENCE
- ---------------
- Last Energy change ... -1.0028e-08 Tolerance : 1.0000e-06
- Last MAX-Density change ... 1.5543e-15 Tolerance : 1.0000e-05
- Last RMS-Density change ... 3.9422e-16 Tolerance : 1.0000e-06
- Last Orbital Gradient ... 1.5388e-05 Tolerance : 5.0000e-05
- Last Orbital Rotation ... 1.3800e-05 Tolerance : 5.0000e-05
- **** THE GBW FILE WAS UPDATED (en.gbw) ****
- **** DENSITY en.scfp WAS UPDATED ****
- **** ENERGY FILE WAS UPDATED (en.en.tmp) ****
- **** THE GBW FILE WAS UPDATED (en.gbw) ****
- **** DENSITY en.scfp WAS UPDATED ****
- ----------------
- ORBITAL ENERGIES
- ----------------
- NO OCC E(Eh) E(eV)
- 0 2.0000 -24.870168 -676.7517
- 1 2.0000 -1.200143 -32.6575
- 2 2.0000 -0.452147 -12.3036
- 3 2.0000 -0.286164 -7.7869
- 4 2.0000 -0.286164 -7.7869
- 5 0.0000 0.449015 12.2183
- ------------------
- MOLECULAR ORBITALS
- ------------------
- 0 1 2 3 4 5
- -24.87017 -1.20014 -0.45215 -0.28616 -0.28616 0.44902
- 2.00000 2.00000 2.00000 2.00000 2.00000 0.00000
- -------- -------- -------- -------- -------- --------
- 0F 1s 0.994638 -0.245080 0.095375 0.000000 -0.000000 -0.083573
- 0F 2s 0.022921 0.915097 -0.482529 -0.000000 0.000000 0.548727
- 0F 1pz 0.000000 -0.000000 0.000000 -0.259678 0.965695 0.000000
- 0F 1px 0.003069 0.126279 0.691623 0.000000 -0.000000 0.826914
- 0F 1py 0.000000 -0.000000 0.000000 -0.965695 -0.259678 0.000000
- 1H 1s -0.006088 0.182232 0.519968 0.000000 -0.000000 -1.082661
- ********************************
- * MULLIKEN POPULATION ANALYSIS *
- ********************************
- -----------------------
- MULLIKEN ATOMIC CHARGES
- -----------------------
- 0 F : -0.202111
- 1 H : 0.202111
- Sum of atomic charges: -0.0000000
- --------------------------------
- MULLIKEN REDUCED ORBITAL CHARGES
- --------------------------------
- 0 F s : 3.947688 s : 3.947688
- pz : 2.000000 p : 5.254423
- px : 1.254423
- py : 2.000000
- 1 H s : 0.797889 s : 0.797889
- *******************************
- * LOEWDIN POPULATION ANALYSIS *
- *******************************
- ----------------------
- LOEWDIN ATOMIC CHARGES
- ----------------------
- 0 F : -0.144987
- 1 H : 0.144987
- -------------------------------
- LOEWDIN REDUCED ORBITAL CHARGES
- -------------------------------
- 0 F s : 3.874758 s : 3.874758
- pz : 2.000000 p : 5.270229
- px : 1.270229
- py : 2.000000
- 1 H s : 0.855013 s : 0.855013
- *****************************
- * MAYER POPULATION ANALYSIS *
- *****************************
- NA - Mulliken gross atomic population
- ZA - Total nuclear charge
- QA - Mulliken gross atomic charge
- VA - Mayer's total valence
- BVA - Mayer's bonded valence
- FA - Mayer's free valence
- ATOM NA ZA QA VA BVA FA
- 0 F 9.2021 9.0000 -0.2021 0.9592 0.9592 -0.0000
- 1 H 0.7979 1.0000 0.2021 0.9592 0.9592 0.0000
- Mayer bond orders larger than 0.100000
- B( 0-F , 1-H ) : 0.9592
- -------
- TIMINGS
- -------
- Total SCF time: 0 days 0 hours 0 min 0 sec
- Total time .... 0.404 sec
- Sum of individual times .... 0.358 sec ( 88.7%)
- Fock matrix formation .... 0.249 sec ( 61.7%)
- Split-RI-J .... 0.125 sec ( 50.2% of F)
- Chain of spheres X .... 0.057 sec ( 22.9% of F)
- XC integration .... 0.040 sec ( 16.2% of F)
- Basis function eval. .... 0.001 sec ( 3.4% of XC)
- Density eval. .... 0.005 sec ( 13.1% of XC)
- XC-Functional eval. .... 0.009 sec ( 22.7% of XC)
- XC-Potential eval. .... 0.001 sec ( 1.3% of XC)
- Diagonalization .... 0.013 sec ( 3.2%)
- Density matrix formation .... 0.000 sec ( 0.1%)
- Population analysis .... 0.003 sec ( 0.8%)
- Initial guess .... 0.007 sec ( 1.7%)
- Orbital Transformation .... 0.000 sec ( 0.0%)
- Orbital Orthonormalization .... 0.000 sec ( 0.0%)
- DIIS solution .... 0.021 sec ( 5.1%)
- SOSCF solution .... 0.004 sec ( 1.0%)
- Grid generation .... 0.061 sec ( 15.0%)
- Maximum memory used throughout the entire SCF-calculation: 42.5 MB
- ------------------------- --------------------
- FINAL SINGLE POINT ENERGY -98.933928582316
- ------------------------- --------------------
- ***************************************
- * ORCA property calculations *
- ***************************************
- ---------------------
- Active property flags
- ---------------------
- (+) Dipole Moment
- ------------------------------------------------------------------------------
- ORCA ELECTRIC PROPERTIES CALCULATION
- ------------------------------------------------------------------------------
- Dipole Moment Calculation ... on
- Quadrupole Moment Calculation ... off
- Polarizability Calculation ... off
- GBWName ... en.gbw
- Electron density ... en.scfp
- The origin for moment calculation is the CENTER OF MASS = (-5.435352, 2.578890 0.000000)
- -------------
- DIPOLE MOMENT
- -------------
- X Y Z
- Electronic contribution: -0.37809 0.00000 0.00000
- Nuclear contribution : 0.86483 0.00000 0.00000
- -----------------------------------------
- Total Dipole Moment : 0.48674 0.00000 0.00000
- -----------------------------------------
- Magnitude (a.u.) : 0.48674
- Magnitude (Debye) : 1.23720
- --------------------
- Rotational spectrum
- --------------------
- Rotational constants in cm-1: 0.000000 20.699118 20.699118
- Rotational constants in MHz : 0.000000 620543.949424 620543.949424
- Dipole components along the rotational axes:
- x,y,z [a.u.] : 0.486743 0.000000 0.000000
- x,y,z [Debye]: 1.237202 0.000000 0.000000
- ************************************************************
- * Program running with 8 parallel MPI-processes *
- * working on a common directory *
- ************************************************************
- -------------------------------------------------------------------------------
- ORCA SCF HESSIAN
- -------------------------------------------------------------------------------
- --------- SHARK INITIALIZATION DONE ---------
- Hessian of the Kohn-Sham DFT energy:
- Kohn-Sham wavefunction type ... RKS
- Hartree-Fock exchange scaling ... 0.540
- Number of operators ... 1
- Number of atoms ... 2
- Basis set dimensions ... 6
- Integral neglect threshold ... 1.0e-10
- Integral primitive cutoff ... 1.0e-11
- Setting up DFT Hessian calculations ...
- HESS GRID
- ---------
- General Integration Accuracy IntAcc ... 4.629
- Radial Grid Type RadialGrid ... OptM3 with GC (2021)
- Angular Grid (max. ang.) AngularGrid ... 5 (Lebedev-434)
- Angular grid pruning method GridPruning ... 4 (adaptive)
- Weight generation scheme WeightScheme... Becke
- Basis function cutoff BFCut ... 1.0000e-10
- Integration weight cutoff WCut ... 1.0000e-14
- Angular grids for H and He will be reduced by one unit
- Partially contracted basis set ... off
- Rotationally invariant grid construction ... off
- Total number of grid points ... 14832
- Total number of batches ... 233
- Average number of points per batch ... 63
- Average number of grid points per atom ... 7416
- Building densities ... done ( 0.0 sec)
- Calculating rho(r) on the grid ... done ( 0.0 sec)
- Building xc-kernel on the grid ... done ( 0.0 sec)
- done ( 0.1 sec)
- Nuclear repulsion Hessian (SHARK) ... done ( 0.0 sec)
- ----------------------------------------------
- Forming right-hand sides of CP-SCF equations ...
- ----------------------------------------------
- One electron integral derivatives (SHARK) ... done ( 0.0 sec)
- Transforming the overlap derivative matrices ... done ( 0.0 sec)
- Making the Q(x) pseudodensities ... done ( 0.0 sec)
- Adding the E*S(x)*S(y) terms to the Hessian ... done ( 0.0 sec)
- Calculating energy weighted overlap derivatives ... done ( 0.0 sec)
- Two electron integral derivatives (SHARK-RI) ... done ( 0.1 sec)
- Exchange-correlation integral derivatives ... done ( 0.0 sec)
- tr(F(y)Q(x)) contribution to the Hessian ... done ( 0.0 sec)
- Response fock operator R(S(x)) (SHARK-RIJCOSX) ... done ( 0.1 sec)
- XC Response fock operator R(S(x)) ... done ( 0.0 sec)
- tr(F(y)S(x)) contribution to the Hessian ... done ( 0.0 sec)
- Transforming and finalizing RHSs ... done ( 0.0 sec)
- ----------------------------------------------
- Solving the CP-SCF equations (RIJCOSX) ...
- ----------------------------------------------
- CP-SCF ITERATION 0: 6.9389e-18 ( 0.1 sec 6/ 6 done)
- *** THE CP-SCF HAS CONVERGED ***
- ... done ( 0.1 sec)
- Forming perturbed density Hessian contributions ... done ( 0.0 sec)
- 2nd integral derivative contribs (SHARK-RI) ... done ( 0.1 sec)
- Exchange-correlation Hessian ... done ( 0.0 sec)
- Dipole derivatives ... (center of mass: -5.435352, 2.578890 0.000000)done ( 0.0 sec)
- Total SCF Hessian time: 0 days 0 hours 0 min 0 sec
- Writing the Hessian file to the disk ... done
- Maximum memory used throughout the entire SCFHESS-calculation: 48.1 MB
- -----------------------
- VIBRATIONAL FREQUENCIES
- -----------------------
- Scaling factor for frequencies = 1.000000000 (already applied!)
- 0: 0.00 cm**-1
- 1: 0.00 cm**-1
- 2: 0.00 cm**-1
- 3: 0.00 cm**-1
- 4: 0.00 cm**-1
- 5: 4829.27 cm**-1
- ------------
- NORMAL MODES
- ------------
- These modes are the Cartesian displacements weighted by the diagonal matrix
- M(i,i)=1/sqrt(m[i]) where m[i] is the mass of the displaced atom
- Thus, these vectors are normalized but *not* orthogonal
- 0 1 2 3 4 5
- 0 0.000000 0.000000 0.000000 0.000000 0.000000 -0.052984
- 1 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
- 2 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
- 3 0.000000 0.000000 0.000000 0.000000 0.000000 0.998595
- 4 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
- 5 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
- -----------
- IR SPECTRUM
- -----------
- Mode freq eps Int T**2 TX TY TZ
- cm**-1 L/(mol*cm) km/mol a.u.
- ----------------------------------------------------------------------------
- 5: 4829.27 0.018642 94.21 0.001205 (-0.034708 0.000000 0.000000)
- * The epsilon (eps) is given for a Dirac delta lineshape.
- ** The dipole moment derivative (T) already includes vibrational overlap.
- The first frequency considered to be a vibration is 5
- The total number of vibrations considered is 1
- --------------------------
- THERMOCHEMISTRY AT 298.15K
- --------------------------
- Temperature ... 298.15 K
- Pressure ... 1.00 atm
- Total Mass ... 20.01 AMU
- The molecule is recognized as being linear
- Throughout the following assumptions are being made:
- (1) The electronic state is orbitally nondegenerate
- (2) There are no thermally accessible electronically excited states
- (3) Hindered rotations indicated by low frequency modes are not
- treated as such but are treated as vibrations and this may
- cause some error
- (4) All equations used are the standard statistical mechanics
- equations for an ideal gas
- (5) All vibrations are strictly harmonic
- freq. 4829.27 E(vib) ... 0.00
- ------------
- INNER ENERGY
- ------------
- The inner energy is: U= E(el) + E(ZPE) + E(vib) + E(rot) + E(trans)
- E(el) - is the total energy from the electronic structure calculation
- = E(kin-el) + E(nuc-el) + E(el-el) + E(nuc-nuc)
- E(ZPE) - the the zero temperature vibrational energy from the frequency calculation
- E(vib) - the the finite temperature correction to E(ZPE) due to population
- of excited vibrational states
- E(rot) - is the rotational thermal energy
- E(trans)- is the translational thermal energy
- Summary of contributions to the inner energy U:
- Electronic energy ... -98.93392858 Eh
- Zero point energy ... 0.01100188 Eh 6.90 kcal/mol
- Thermal vibrational correction ... 0.00000000 Eh 0.00 kcal/mol
- Thermal rotational correction ... 0.00094418 Eh 0.59 kcal/mol
- Thermal translational correction ... 0.00141627 Eh 0.89 kcal/mol
- -----------------------------------------------------------------------
- Total thermal energy -98.92056625 Eh
- Summary of corrections to the electronic energy:
- (perhaps to be used in another calculation)
- Total thermal correction 0.00236045 Eh 1.48 kcal/mol
- Non-thermal (ZPE) correction 0.01100188 Eh 6.90 kcal/mol
- -----------------------------------------------------------------------
- Total correction 0.01336233 Eh 8.38 kcal/mol
- --------
- ENTHALPY
- --------
- The enthalpy is H = U + kB*T
- kB is Boltzmann's constant
- Total free energy ... -98.92056625 Eh
- Thermal Enthalpy correction ... 0.00094421 Eh 0.59 kcal/mol
- -----------------------------------------------------------------------
- Total Enthalpy ... -98.91962204 Eh
- Note: Rotational entropy computed according to Herzberg
- Infrared and Raman Spectra, Chapter V,1, Van Nostrand Reinhold, 1945
- Point Group: Cinfv, Symmetry Number: 1
- Rotational constants in cm-1: 0.000000 20.699118 20.699118
- Vibrational entropy computed according to the QRRHO of S. Grimme
- Chem.Eur.J. 2012 18 9955
- -------
- ENTROPY
- -------
- The entropy contributions are T*S = T*(S(el)+S(vib)+S(rot)+S(trans))
- S(el) - electronic entropy
- S(vib) - vibrational entropy
- S(rot) - rotational entropy
- S(trans)- translational entropy
- The entropies will be listed as multiplied by the temperature to get
- units of energy
- Electronic entropy ... 0.00000000 Eh 0.00 kcal/mol
- Vibrational entropy ... -0.00000000 Eh -0.00 kcal/mol
- Rotational entropy ... 0.00311931 Eh 1.96 kcal/mol
- Translational entropy ... 0.01659227 Eh 10.41 kcal/mol
- -----------------------------------------------------------------------
- Final entropy term ... 0.01971157 Eh 12.37 kcal/mol
- -------------------
- GIBBS FREE ENERGY
- -------------------
- The Gibbs free energy is G = H - T*S
- Total enthalpy ... -98.91962204 Eh
- Total entropy correction ... -0.01971157 Eh -12.37 kcal/mol
- -----------------------------------------------------------------------
- Final Gibbs free energy ... -98.93933362 Eh
- For completeness - the Gibbs free energy minus the electronic energy
- G-E(el) ... -0.00540503 Eh -3.39 kcal/mol
- Timings for individual modules:
- Sum of individual times ... 2.477 sec (= 0.041 min)
- GTO integral calculation ... 0.886 sec (= 0.015 min) 35.8 %
- SCF iterations ... 0.585 sec (= 0.010 min) 23.6 %
- Analytical frequency calculation... 1.007 sec (= 0.017 min) 40.6 %
- ****ORCA TERMINATED NORMALLY****
- TOTAL RUN TIME: 0 days 0 hours 0 minutes 2 seconds 776 msec
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement