***************** * O R C A * ***************** #, ### #### ##### ###### ########, ,,################,,,,, ,,#################################,, ,,##########################################,, ,#########################################, ''#####, ,#############################################,, '####, ,##################################################,,,,####, ,###########'''' ''''############################### ,#####'' ,,,,##########,,,, '''####''' '#### ,##' ,,,,###########################,,, '## ' ,,###'''' '''############,,, ,,##'' '''############,,,, ,,,,,,###'' ,#'' '''#######################''' ' ''''####'''' ,#######, #######, ,#######, ## ,#' '#, ## ## ,#' '#, #''# ,####, ,#, ## ## ## ,#' ## #' '# #' ,# # ## ## ####### ## ,######, #####, # '#, ,#' ## ## '#, ,#' ,# #, #, # # '#######' ## ## '#######' #' '# '####' # # ######################################################### # -***- # # Department of theory and spectroscopy # # # # Frank Neese # # # # Directorship, Architecture, Infrastructure # # SHARK, DRIVERS # # Core code/Algorithms in most modules # # # # Max Planck Institute fuer Kohlenforschung # # Kaiser Wilhelm Platz 1 # # D-45470 Muelheim/Ruhr # # Germany # # # # All rights reserved # # -***- # ######################################################### Program Version 6.1.1 - RELEASE - With contributions from (in alphabetic order): [Max-Planck-Institut fuer Kohlenforschung] 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 : All parallelization in ORCA, NUMFREQ, NUMCALC Giovanni Bistoni : ED, misc. LED, open-shell LED, HFLD Dmytro Bykov : pre 5.0 version of the SCF Hessian Marcos Casanova-Páez : Triplet and SCS-CIS(D). UHF-(DLPNO)-IP/EA/STEOM-CCSD. UHF-CVS-IP/STEOM-CCSD Vijay G. Chilkuri : MRCI spin determinant printing, contributions to CSF-ICE Pauline Colinet : FMM embedding Dipayan Datta : RHF DLPNO-CCSD density Achintya Kumar Dutta : EOM-CC, STEOM-CC Nicolas Foglia : Exact transition moments, OPA infrastructure, MCD improvements Dmitry Ganyushin : Spin-Orbit,Spin-Spin,Magnetic field MRCI Miquel Garcia-Rates : C-PCM and meta-GGA Hessian, CCSD/C-PCM, Gaussian charge scheme Tiago L. C. Gouveia : GS-ROHF, GS-ROCIS Yang Guo : DLPNO-NEVPT2, F12-NEVPT2, CIM, IAO-localization Andreas Hansen : Spin unrestricted coupled pair/coupled cluster methods Ingolf Harden : AUTO-CI MPn and infrastructure Benjamin Helmich-Paris : MC-RPA, TRAH-(SCF,CASSCF), AVAS, COSX integrals, SCF dyn. polar., MC-PDFT, srDFT Lee Huntington : MR-EOM, pCC Robert Izsak : Overlap fitted RIJCOSX, COSX-SCS-MP3, EOM Riya Kayal : Wick's Theorem for AUTO-CI, AUTO-CI UHF-CCSDT Emily Kempfer : AUTO-CI RHF CISDT and CCSDT, approximate NEVPT4 Christian Kollmar : KDIIS, OOCD, Brueckner-CCSD(T), CCSD density, CASPT2, CASPT2-K, improved NEVPT2 Axel Koslowski : Symmetry handling Simone Kossmann : meta-GGA functionals, TD-DFT gradient, OOMP2, (MP2 Hessian; deprecated post 5.0) Lucas Lang : DCDCAS, Hyperfine gauge corrections, ICE-SOC+SSC Marvin Lechner : AUTO-CI (C++ implementation), FIC-MRCC Spencer Leger : CASSCF response Dagmar Lenk : GEPOL surface, SMD, ORCA-2-JSON Dimitrios Liakos : Extrapolation schemes; Compound Job, Property file Dimitrios Manganas : Further ROCIS development; embedding schemes. LFT, Crystal Embedding Dimitrios Pantazis : SARC Basis sets Anastasios Papadopoulos: AUTO-CI, single reference methods and gradients Taras Petrenko : pre 6.0 DFT Hessian and TD-DFT gradient, ECA, NRVS Petra Pikulova : Analytic Raman intensities Peter Pinski : DLPNO-MP2, DLPNO-MP2 Gradient Shashank Vittal Rao : ES-AILFT, MagRelax Christoph Reimann : Effective Core Potentials Marius Retegan : Local ZFS, SOC Christoph Riplinger : Optimizer, TS searches, QM/MM, DLPNO-CCSD(T), (RO)-DLPNO pert. Triples Michael Roemelt : Original ROCIS implementation, recursive CI coupling coefficients Masaaki Saitow : Open-shell DLPNO-CCSD energy and density Barbara Sandhoefer : DKH picture change effects Yorick L. A. Schmerwitz: GMF and freeze-and-release deltaSCF, NEB S-IDPP initial path Kantharuban Sivalingam : CASSCF convergence/infrastructure, NEVPT2, NEVPT3, NEVPT4(SD), FIC-MRCI and CEPA variants Bernardo de Souza : ESD, SOC TD-DFT Georgi L. Stoychev : AutoAux, RI-MP2 NMR, DLPNO-MP2 response, X2C Van Anh Tran : RI-MP2 g-tensors Willem Van den Heuvel : Paramagnetic NMR Zikuan Wang : NOTCH, Electric field optimization Frank Wennmohs : Technical directorship and infrastructure Hang Xu : AUTO-CI-Response properties [FACCTs GmbH] Markus Bursch, Nicolas Foglia, Miquel Garcia-Rates, Ingolf Harden, Hagen Neugebauer, Anastasios Papadopoulos, Christoph Riplinger, Bernardo de Souza, Georgi L. Stoychev APM, various basis sets, CI-OPT, improved COSX, DLPNO-Multilevel, DOCKER, DRACO, updates on ESD, Fragmentator, GOAT, IRC, LR-CPCM, L-BFGS, MBIS, meta-GGA TD-DFT gradient, ML-optimized integration grids, MM, NACMEs, nearIR, NEB, NEB-TS, NL-DFT gradient (VV10), 2- and 3-layer-ONIOM, interface openCOSMO-RS, QMMM, Crystal-QMMM, RESP, rigid body optimization, SF, symmetry and pop. for TD-DFT, various functionals, SOLVATOR [Other institutions] V. Asgeirsson : NEB Christoph Bannwarth : sTDA-DFT, sTD-DFT, PBEh-3c, B97-3c, D3 Giovanni Bistoni : ETS/NOCV, ADLD/ADEX, COVALED Martin Brehm : Molecular dynamics Ronald Cardenas : ETS/NOCV Martina Colucci : COVALED Sebastian Ehlert : rSCAN, r2SCAN, r2SCAN-3c, D4, dhf basis sets Marvin Friede : D4 for Fr, Ra, Ac-Lr Lars Goerigk : TD-DFT with DH, B97 family of functionals Stefan Grimme : VdW corrections, initial TS optimization, DFT functionals, gCP, sTDA/sTD-DF Waldemar Hujo : DFT-NL H. Jonsson : NEB Holger Kruse : gCP Marcel Mueller : wB97X-3c, vDZP basis set Hagen Neugebauer : wr2SCAN, Native XTB Gianluca Regni : ADLD/ADEX Tobias Risthaus : pre 6.0 range-separated hybrid DFT and stability analysis Lukas Wittmann : regularized MP2, r2SCAN double-hybrids, wr2SCAN We gratefully acknowledge several colleagues who have allowed us to interface, adapt or use parts of their codes: 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) Frank Weinhold : gennbo (NPA and NBO analysis) Simon Mueller : openCOSMO-RS Christopher J. Cramer and Donald G. Truhlar : smd solvation model 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: 7.0.0 For citations please refer to: https://libxc.gitlab.io This ORCA versions uses: CBLAS interface : Fast vector & matrix operations LAPACKE interface : Fast linear algebra routines Shared memory : Shared parallel matrices *********************************** * Starting time: Fri Jun 5 17:21:46 2026 * Host name: KDH02 * Process ID: 19984 * Working dir.: C:\Users\denny\orca_work\examples\01_sp *********************************** Warning: RI is on but no J-basis has been assigned. Assigning Def2/J (nothing to worry about!) ================================================================================ ----- Orbital basis set information ----- Your calculation utilizes the basis: def2-SVP F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys. 7, 3297 (2005). ----- AuxJ basis set information ----- Your calculation utilizes the auxiliary basis: def2/J H-Rn: F. Weigend, Phys. Chem. Chem. Phys. 8, 1057 (2006). Fr-Lr: K. Eichkorn, F. Weigend, O. Treutler, R. Ahlrichs; Theor. Chem. Acc. 97, 119 (1997). ================================================================================ WARNINGS Please study these warnings very carefully! ================================================================================ ================================================================================ INPUT FILE ================================================================================ NAME = water.inp | 1> # Single-point energy - water, B3LYP/def2-SVP | 2> ! B3LYP def2-SVP | 3> * xyz 0 1 | 4> O 0.000000 0.000000 0.119262 | 5> H 0.000000 0.763239 -0.477047 | 6> H 0.000000 -0.763239 -0.477047 | 7> * | 8> | 9> ****END OF INPUT**** ================================================================================ **************************** * Single Point Calculation * **************************** --------------------------------- CARTESIAN COORDINATES (ANGSTROEM) --------------------------------- O 0.000000 0.000000 0.119262 H 0.000000 0.763239 -0.477047 H 0.000000 -0.763239 -0.477047 ---------------------------- CARTESIAN COORDINATES (A.U.) ---------------------------- NO LB ZA FRAG MASS X Y Z 0 O 8.0000 0 15.999 0.000000 0.000000 0.225373 1 H 1.0000 0 1.008 0.000000 1.442313 -0.901488 2 H 1.0000 0 1.008 0.000000 -1.442313 -0.901488 -------------------------------- INTERNAL COORDINATES (ANGSTROEM) -------------------------------- O 0 0 0 0.000000000000 0.00000000 0.00000000 H 1 0 0 0.968565018263 0.00000000 0.00000000 H 1 2 0 0.968565018263 103.99987510 0.00000000 --------------------------- INTERNAL COORDINATES (A.U.) --------------------------- O 0 0 0 0.000000000000 0.00000000 0.00000000 H 1 0 0 1.830322627413 0.00000000 0.00000000 H 1 2 0 1.830322627413 103.99987510 0.00000000 --------------------- BASIS SET INFORMATION --------------------- There are 2 groups of distinct atoms Group 1 Type O : 7s4p1d contracted to 3s2p1d pattern {511/31/1} Group 2 Type H : 4s1p contracted to 2s1p pattern {31/1} Atom 0O basis set group => 1 Atom 1H basis set group => 2 Atom 2H basis set group => 2 --------------------------------- AUXILIARY/J BASIS SET INFORMATION --------------------------------- There are 2 groups of distinct atoms Group 1 Type O : 12s5p4d2f1g contracted to 6s4p3d1f1g pattern {711111/2111/211/2/1} Group 2 Type H : 5s2p1d contracted to 3s1p1d pattern {311/2/1} Atom 0O basis set group => 1 Atom 1H basis set group => 2 Atom 2H basis set group => 2 ------------------------------------------------------------------------------ ORCA STARTUP CALCULATIONS -- 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, v3 2022-2024 ------------------------------------------------------------------------------ ---------------------- SHARK INTEGRAL PACKAGE ---------------------- Number of atoms ... 3 Number of basis functions ... 24 Number of shells ... 12 Maximum angular momentum ... 2 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 Prescreening option ... SCHWARTZ Thresh ... 1.000e-10 Tcut ... 1.000e-11 Tpresel ... 1.000e-11 Coulomb Range Separation ... NOT USED Exchange Range Separation ... NOT USED Multipole approximations ... NOT USED Finite Nucleus Model ... NOT USED CABS basis ... NOT available Auxiliary Coulomb fitting basis ... AVAILABLE # of basis functions in Aux-J ... 71 # of shells in Aux-J ... 25 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 Checking pre-screening integrals ... done ( 0.0 sec) Dimension = 12 => SHARK Basis and OBASIS are compatible. Storing Pre-screening Shell pair information Shell pair cut-off parameter TPreSel ... 1.0e-11 Total number of shell pairs ... 78 Shell pairs after pre-screening ... 78 Total number of primitive shell pairs ... 272 Primitive shell pairs kept ... 265 la=0 lb=0: 28 shell pairs la=1 lb=0: 28 shell pairs la=1 lb=1: 10 shell pairs la=2 lb=0: 7 shell pairs la=2 lb=1: 4 shell pairs la=2 lb=2: 1 shell pairs Checking whether 4 symmetric matrices of dimension 24 fit in memory :Max Core in MB = 4096.00 MB in use = 3.20 MB left = 4092.80 MB needed = 0.01 Data fit in memory = YES Calculating RI/J V-Matrix + Cholesky decomp.... done ( 0.0 sec) Calculating Nuclear repulsion ... done ( 0.0 sec) ENN= 9.088293724142 Eh Diagonalization of the overlap matrix: Smallest eigenvalue ... 3.828e-02 Time for diagonalization ... 0.002 sec Threshold for overlap eigenvalues ... 1.000e-07 Number of eigenvalues below threshold ... 0 Time for construction of square roots ... 0.003 sec Total time needed ... 0.006 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... mBecke (2022) Basis function cutoff BFCut ... 1.0000e-10 Integration weight cutoff WCut ... 1.0000e-14 Partially contracted basis set ... off Rotationally invariant grid construction ... off Angular grids for H and He will be reduced by one unit Total number of grid points ... 12754 Total number of batches ... 201 Average number of points per batch ... 63 Average number of grid points per atom ... 4251 -------------------- 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... mBecke (2022) Basis function cutoff BFCut ... 1.0000e-10 Integration weight cutoff WCut ... 1.0000e-14 Partially contracted basis set ... on Rotationally invariant grid construction ... off Angular grids for H and He will be reduced by one unit Total number of grid points ... 1650 Total number of batches ... 14 Average number of points per batch ... 117 Average number of grid points per atom ... 550 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... mBecke (2022) Basis function cutoff BFCut ... 1.0000e-10 Integration weight cutoff WCut ... 1.0000e-14 Partially contracted basis set ... on Rotationally invariant grid construction ... off Angular grids for H and He will be reduced by one unit Total number of grid points ... 3496 Total number of batches ... 29 Average number of points per batch ... 120 Average number of grid points per atom ... 1165 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... mBecke (2022) Basis function cutoff BFCut ... 1.0000e-10 Integration weight cutoff WCut ... 1.0000e-14 Partially contracted basis set ... on Rotationally invariant grid construction ... off Angular grids for H and He will be reduced by one unit Total number of grid points ... 7784 Total number of batches ... 62 Average number of points per batch ... 125 Average number of grid points per atom ... 2595 UseSFitting ... on Grids setup in 0.1 sec Initializing property integral containers ... done ( 0.0 sec) SHARK setup successfully completed in 0.2 seconds Maximum memory used throughout the entire STARTUP-calculation: 8.8 MB ------------------------------------------------------------------------------- ORCA GUESS Start orbitals & Density for SCF / CASSCF ------------------------------------------------------------------------------- ------------ SCF SETTINGS ------------ Hamiltonian: Density Functional Method .... DFT(GTOs) Exchange Functional Exchange .... B88 X-Alpha parameter XAlpha .... 0.666667 Becke's b parameter XBeta .... 0.004200 Correlation Functional Correlation .... LYP LDA part of GGA corr. LDAOpt .... VWN-5 Gradients option PostSCFGGA .... off Hybrid DFT is turned on Fraction HF Exchange ScalHFX .... 0.200000 Scaling of DF-GGA-X ScalDFX .... 0.720000 Scaling of DF-GGA-C ScalDFC .... 0.810000 Scaling of DF-LDA-C ScalLDAC .... 1.000000 Perturbative correction .... 0.000000 NL short-range parameter .... 4.800000 RI-approximation to the Coulomb term is turned on Number of AuxJ basis functions .... 71 RIJ-COSX (HFX calculated with COS-X)).... on General Settings: Integral files IntName .... water Hartree-Fock type HFTyp .... RHF Total Charge Charge .... 0 Multiplicity Mult .... 1 Number of Electrons NEL .... 10 Basis Dimension Dim .... 24 Nuclear Repulsion ENuc .... 9.0882937241 Eh Convergence Acceleration: AO-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 MO-DIIS CNVKDIIS .... off Trust-Rad. Augm. Hess. CNVTRAH .... auto Auto Start mean grad. ratio tolernc. .... 1.125000 Auto Start start iteration .... 50 Auto Start num. interpolation iter. .... 10 Max. Number of Micro iterations .... 24 Max. Number of Macro iterations .... Maxiter - #DIIS iter Number of Davidson start vectors .... 2 Converg. threshold (grad. norm) .... 5.000e-05 Grad. Scal. Fac. for Micro threshold .... 0.100 Minimum threshold for Micro iter. .... 1.000e-02 NR start threshold (gradient norm) .... 1.000e-04 Initial trust radius .... 0.400 Minimum AH scaling param. (alpha) .... 1.000 Maximum AH scaling param. (alpha) .... 1000.000 Quad. conv. algorithm .... NR White noise on init. David. guess .... on Maximum white noise .... 0.010 Pseudo random numbers .... off Inactive MOs .... canonical Orbital update algorithm .... Taylor Preconditioner .... Diag Full preconditioner red. dimension .... 250 SOSCF CNVSOSCF .... on Start iteration SOSCFMaxIt .... 150 Startup grad/error SOSCFStart .... 0.003300 Hessian update SOSCFHessUp .... L-BFGS Autom. constraints SOSCFAutoConstrain .... off 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 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 ------------------------------ 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) Making the grid ... done ( 0.0 sec) Mapping shells ... done Starting the XC term evaluation ... done ( 0.0 sec) promolecular density results # of electrons = 9.999259806 EX = -8.779131706 EC = -0.334384305 EX+EC = -9.113516011 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.1 sec) ------------------ **** ENERGY FILE WAS UPDATED (water.en.tmp) **** Finished Guess after 0.1 sec Maximum memory used throughout the entire GUESS-calculation: 5.4 MB ------------------------------------------------------------------------------------------- ORCA LEAN-SCF memory conserving SCF solver ------------------------------------------------------------------------------------------- ----------------------------------------D-I-I-S-------------------------------------------- Iteration Energy (Eh) Delta-E RMSDP MaxDP DIISErr Damp Time(sec) ------------------------------------------------------------------------------------------- *** Starting incremental Fock matrix formation *** 1 -76.2268238847872937 0.00e+00 1.32e-02 8.80e-02 3.06e-01 0.700 0.1 2 -76.2733081135723978 -4.65e-02 7.46e-03 5.02e-02 1.29e-01 0.700 0.1 ***Turning on AO-DIIS*** 3 -76.2890197169736552 -1.57e-02 1.90e-03 1.26e-02 3.41e-02 0.700 0.1 4 -76.2988422159047985 -9.82e-03 2.41e-03 1.59e-02 1.41e-02 0.000 0.1 5 -76.3213131728440857 -2.25e-02 5.56e-04 3.21e-03 4.70e-03 0.000 0.1 *** Initializing SOSCF *** ---------------------------------------S-O-S-C-F-------------------------------------- Iteration Energy (Eh) Delta-E RMSDP MaxDP MaxGrad Time(sec) -------------------------------------------------------------------------------------- 6 -76.3213541719220245 -4.10e-05 1.67e-04 1.12e-03 7.30e-04 0.2 *** Restarting incremental Fock matrix formation *** 7 -76.3212858050746661 6.84e-05 1.06e-04 8.97e-04 2.23e-04 0.1 8 -76.3212860247906377 -2.20e-07 3.28e-05 2.32e-04 1.49e-04 0.1 **** Energy Check signals convergence **** ***************************************************** * SUCCESS * * SCF CONVERGED AFTER 8 CYCLES * ***************************************************** Recomputing exchange energy using gridx3 ... done ( 0.038 sec) Old exchange energy : -1.789659509 Eh New exchange energy : -1.789647862 Eh Exchange energy change after final integration : 0.000011647 Eh Total energy after final integration : -76.321274411 Eh **** ENERGY FILE WAS UPDATED (water.en.tmp) **** ---------------- TOTAL SCF ENERGY ---------------- Total Energy : -76.32127441114484 Eh -2076.80746 eV Components: Nuclear Repulsion : 9.08829372414239 Eh 247.30504 eV Electronic Energy : -85.40957978265600 Eh -2324.11282 eV One Electron Energy: -122.92827980607993 Eh -3345.04855 eV Two Electron Energy: 37.51870002342393 Eh 1020.93573 eV Virial components: Potential Energy : -152.14561447236935 Eh -4140.09265 eV Kinetic Energy : 75.82434006122452 Eh 2063.28519 eV Virial Ratio : 2.00655375766567 DFT components: N(Alpha) : 4.999998186273 electrons N(Beta) : 4.999998186273 electrons N(Total) : 9.999996372546 electrons E(X) : -7.120829933981 Eh E(C) : -0.402136766460 Eh E(XC) : -7.522966700442 Eh --------------- SCF CONVERGENCE --------------- Last Energy change ... 2.1972e-07 Tolerance : 1.0000e-06 Last MAX-Density change ... 2.3157e-04 Tolerance : 1.0000e-05 Last RMS-Density change ... 3.2831e-05 Tolerance : 1.0000e-06 Last DIIS Error ... 7.2997e-04 Tolerance : 1.0000e-06 Last Orbital Gradient ... 1.4869e-04 Tolerance : 5.0000e-05 Last Orbital Rotation ... 9.0713e-05 Tolerance : 5.0000e-05 ---------------- ORBITAL ENERGIES ---------------- NO OCC E(Eh) E(eV) 0 2.0000 -19.116064 -520.1746 1 2.0000 -0.974696 -26.5228 2 2.0000 -0.509052 -13.8520 3 2.0000 -0.362494 -9.8640 4 2.0000 -0.287790 -7.8312 5 0.0000 0.047531 1.2934 6 0.0000 0.126225 3.4348 7 0.0000 0.559920 15.2362 8 0.0000 0.619174 16.8486 9 0.0000 0.920193 25.0397 10 0.0000 0.923716 25.1356 11 0.0000 1.009683 27.4749 12 0.0000 1.096404 29.8347 13 0.0000 1.348157 36.6852 14 0.0000 1.410908 38.3928 15 0.0000 1.550285 42.1854 *Only the first 10 virtual orbitals were printed. ******************************** * MULLIKEN POPULATION ANALYSIS * ******************************** ----------------------- MULLIKEN ATOMIC CHARGES ----------------------- 0 O : -0.297627 1 H : 0.148813 2 H : 0.148813 Sum of atomic charges: -0.0000000 -------------------------------- MULLIKEN REDUCED ORBITAL CHARGES -------------------------------- 0 O s : 3.687048 s : 3.687048 pz : 1.499228 p : 4.599836 px : 1.959651 py : 1.140957 dz2 : 0.001493 d : 0.010743 dxz : 0.001199 dyz : 0.007297 dx2y2 : 0.000753 dxy : 0.000000 1 H s : 0.781083 s : 0.781083 pz : 0.021862 p : 0.070104 px : 0.019575 py : 0.028667 2 H s : 0.781083 s : 0.781083 pz : 0.021862 p : 0.070104 px : 0.019575 py : 0.028667 ******************************* * LOEWDIN POPULATION ANALYSIS * ******************************* ---------------------- LOEWDIN ATOMIC CHARGES ---------------------- 0 O : -0.151941 1 H : 0.075970 2 H : 0.075970 ------------------------------- LOEWDIN REDUCED ORBITAL CHARGES ------------------------------- 0 O s : 3.439936 s : 3.439936 pz : 1.571684 p : 4.692137 px : 1.893086 py : 1.227367 dz2 : 0.000937 d : 0.019868 dxz : 0.000369 dyz : 0.015320 dx2y2 : 0.003241 dxy : 0.000000 1 H s : 0.749916 s : 0.749916 pz : 0.060743 p : 0.174113 px : 0.053272 py : 0.060097 2 H s : 0.749916 s : 0.749916 pz : 0.060743 p : 0.174113 px : 0.053272 py : 0.060097 ***************************** * 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 O 8.2976 8.0000 -0.2976 2.0142 2.0142 -0.0000 1 H 0.8512 1.0000 0.1488 1.0127 1.0127 0.0000 2 H 0.8512 1.0000 0.1488 1.0127 1.0127 0.0000 Mayer bond orders larger than 0.100000 B( 0-O , 1-H ) : 1.0071 B( 0-O , 2-H ) : 1.0071 ------- TIMINGS ------- Total SCF time: 0 days 0 hours 0 min 1 sec Total time .... 1.225 sec Sum of individual times .... 1.127 sec ( 92.0%) SCF preparation .... 0.069 sec ( 5.6%) Fock matrix formation .... 0.483 sec ( 39.4%) Startup .... 0.005 sec ( 1.0% of F) Split-RI-J .... 0.034 sec ( 7.0% of F) Chain of spheres X .... 0.213 sec ( 44.1% of F) XC integration .... 0.162 sec ( 33.5% of F) Basis function eval. .... 0.045 sec ( 27.8% of XC) Density eval. .... 0.020 sec ( 12.3% of XC) XC-Functional eval. .... 0.016 sec ( 9.9% of XC) XC-Potential eval. .... 0.018 sec ( 11.1% of XC) Diagonalization .... 0.000 sec ( 0.0%) Density matrix formation .... 0.119 sec ( 9.7%) Total Energy calculation .... 0.047 sec ( 3.8%) Population analysis .... 0.032 sec ( 2.6%) Orbital Transformation .... 0.019 sec ( 1.6%) Orbital Orthonormalization .... 0.000 sec ( 0.0%) DIIS solution .... 0.185 sec ( 15.1%) SOSCF solution .... 0.173 sec ( 14.1%) Finished LeanSCF after 1.3 sec Maximum memory used throughout the entire LEANSCF-calculation: 4.9 MB ------------------------- -------------------- FINAL SINGLE POINT ENERGY -76.321274411145 ------------------------- -------------------- ------------------------------------------------------------------------------ ORCA PROPERTY CALCULATIONS ------------------------------------------------------------------------------ GBWName ... water.gbw Number of atoms ... 3 Number of basis functions ... 24 Max core memory ... 4096 MB Electric properties: Dipole moment ... YES Quadrupole moment ... NO Static polarizability (Dipole/Dipole) ... NO Static polarizability (Dipole/Quad.) ... NO Static polarizability (Quad./Quad.) ... NO Static polarizability (Velocity) ... NO Static hyperpolarizability ... NO Atomic electric properties: Dipole moment ... NO Quadrupole moment ... NO Static polarizability ... NO Choice of electric origin ... Center of mass Position of electric origin ... 0.000000 0.000000 0.099269 General magnetic properties: Magnetizability ... NO EPR properties: g-Tensor (aka g-matrix) ... NO Zero-Field splitting spin-orbit ... NO Zero-field splitting spin-spin ... NO Hyperfine couplings ... NO ( 0 nuclei) Quadrupole couplings ... NO ( 0 nuclei) Contact density ... NO ( 0 nuclei) NMR properties: Chemical shifts ... NO ( 0 nuclei) Spin-rotation constants ... NO ( 0 nuclei) Spin-spin couplings ... NO ( 0 nuclei, 0 pairs) Choice of magnetic origin ... GIAO Position of magnetic origin ... 0.000000 0.000000 0.000000 Properties with geometric perturbations: SCF Hessian ... NO IR spectrum ... NO VCD spectrum ... NO X-ray spectroscopy properties: SCF XES/XAS/RIXS spectra ... NO SCF SOC stabilization energy ... NO Diagonal Born-Oppenheimer correction ... NO ------------- DIPOLE MOMENT ------------- Method : SCF Type of density : Electron Density Multiplicity : 1 Irrep : 0 Energy : -76.3212744111448416 Eh Basis : AO X Y Z Electronic contribution: -0.000000001 -0.000000005 0.208533992 Nuclear contribution : 0.000000000 0.000000000 -0.992688278 ----------------------------------------- Total Dipole Moment : -0.000000001 -0.000000005 -0.784154286 ----------------------------------------- Magnitude (a.u.) : 0.784154286 Magnitude (Debye) : 1.993161796 -------------------- Rotational spectrum -------------------- Rotational constants in cm-1: 26.479190 14.354395 9.308336 Rotational constants in MHz : 793826.144330 430333.930912 279056.907719 Dipole components along the rotational axes: x,y,z [a.u.] : -0.000000 -0.784154 -0.000000 x,y,z [Debye]: -0.000000 -1.993162 -0.000000 Dipole moment calculation done in 0.0 sec Maximum memory used throughout the entire PROP-calculation: 3.3 MB -------------------------------- SUGGESTED CITATIONS FOR THIS RUN -------------------------------- Below you find a list of papers that are relevant to this ORCA run We neither can nor want to force you to cite these papers, but we appreciate if you do You receive ORCA, which is the product of decades of hard work by many enthusiastic individuals, for free The only thing we kindly ask in return is that you cite our papers, We deeply appreciate it, if you show your appreciation for ORCA by not just citing the generic ORCA reference. Please note that relegating all ORCA citations to the supporting information does *not* help us. SI sections are not indexed - citations you put there will not count into any citation statistics But we need these citations in order to attract the funding resources that allow us to do what we are doing Therefore, if you are a happy ORCA user, please consider citing a few of the papers listed below in the main body of your paper In addition to the list printed below, the program has created the file water.bibtex that contains the list in bibtex format You can import this file easily into all common literature databanks and citation aid programs It goes without saying that in many instances, there are alternative algorithms to achieve similar results as the ones you have gotten from ORCA. It is, of course, also the case that in some instances ORCA just re-implements algorithms worked out by others. We are fully aware of that and we are also fully appreciative of our colleagues work. Hence this citation list should not be read as indicating that the listed papers, which are focused on our own work, are the only ones worth citing. It simply meant to make it easier for users to cite ORCA specific papers. It is not a substitute for doing your own literature research and citing the relevant literature in a scientifically appropriate manner. List of essential papers. We consider these as the minimum necessary citations 1. Neese, F. Software update: the ORCA program system, version 6.0 WIRES Comput. Molec. Sci. 2025 15(1), e70019 doi.org/10.1002/wcms.70019 List of papers to cite with high priority. The work reported in these papers was absolutely necessary for this run to complete. Our perspective: the developers of density functionals and basis sets usually get cited in chemistry papers Good! But without the algorithms to do something with them, the functionals or basis sets would not do anything. Hence, in our opinion, the algorithm design and method developments papers are equally worthy of getting cited 1. Neese, F. An improvement of the resolution of the identity approximation for the formation of the Coulomb matrix J. Comp. Chem. 2003 24(14), 1740-1747 doi.org/10.1002/jcc.10318 2. Neese, F.; Wennmohs, F.; Hansen, A.; Becker, U. Efficient, approximate and parallel Hartree-Fock and hybrid DFT calculations. A 'chain-of-spheres' algorithm for the Hartree-Fock exchange Chem. Phys. 2009 356(1-3), 98-109 doi.org/10.1016/j.chemphys.2008.10.036 3. Helmich-Paris, B.; de Souza, B.; Neese, F.; Izsák, R. An improved chain of spheres for exchange algorithm J. Chem. Phys. 2021 155(10), 104109 doi.org/10.1063/5.0058766 4. Neese, F. The SHARK Integral Generation and Digestion System J. Comp. Chem. 2022 44(3), 381 doi.org/10.1002/jcc.26942 List of suggested additional citations. These are papers that are important in the 'surrounding' of of this run, or papers that preceded the highly important papers. If you like your results we are grateful for a citation. 1. Izsak, R.; Neese, F. An overlap fitted chain of spheres exchange method J. Chem. Phys. 2011 135 , 144105 doi.org/10.1063/1.3646921 2. Izsak, R.; Hansen, A.; Neese, F. The resolution of identity and chain of spheres approximations for the LPNO-CCSD singles Fock term Molec. Phys. 2012 110 , 2413-2417 doi.org/10.1080/00268976.2012.687466 3. Neese, F. The ORCA program system WIRES Comput. Molec. Sci. 2012 2(1), 73-78 doi.org/10.1002/wcms.81 4. Izsak, R.; Neese, F.; Klopper, W. Robust fitting techniques in the chain of spheres approximation to the Fock exchange: The role of the complementary space J. Chem. Phys. 2013 139 , doi.org/10.1063/1.4819264 5. Neese, F. Software update: the ORCA program system, version 4.0 WIRES Comput. Molec. Sci. 2018 8(1), 1-6 doi.org/10.1002/wcms.1327 6. Neese, F.; Wennmohs, F.; Becker, U.; Riplinger, C. The ORCA quantum chemistry program package J. Chem. Phys. 2020 152(22), 224108 doi.org/10.1063/5.0004608 7. Neese, F. Software update: The ORCA program system—Version 5.0 WIRES Comput. Molec. Sci. 2022 12(1), e1606 doi.org/10.1002/wcms.1606 List of optional additional citations 1. Neese, F. Approximate second-order SCF convergence for spin unrestricted wavefunctions Chem. Phys. Lett. 2000 325(1-3), 93-98 doi.org/10.1016/s0009-2614(00)00662-x Timings for individual modules: Sum of individual times ... 2.234 sec (= 0.037 min) Startup calculation ... 0.326 sec (= 0.005 min) 14.6 % SCF iterations ... 1.728 sec (= 0.029 min) 77.4 % Property calculations ... 0.180 sec (= 0.003 min) 8.1 % ****ORCA TERMINATED NORMALLY**** TOTAL RUN TIME: 0 days 0 hours 0 minutes 2 seconds 657 msec