nwout file for Id=66558 bylaska@archive.emsl.pnl.gov:chemdb2/74/97/nwchemarrows-2021-11-20-11-13-134211.out-267365-2021-11-20-11:40:38 argument 1 = /home/bylaska/Projects/Work/RUNARROWS0/nwchemarrows-2021-11-20-11-13-134211.nw ============================== echo of input deck ============================== permanent_dir /home/bylaska/Projects/Work/RUNARROWS0 scratch_dir /home/bylaska/Projects/Work/RUNARROWS0 ######################### START NWCHEM INPUT DECK - NWJOB 134211 ######################## # # NWChemJobId: 6198935409fd0929a2452554 # # NWChem Input Generation (tnt_submit5) - The current time is Fri Nov 19 22:18:54 2021 # - adding tag osmiles:[C]:osmiles to input deck. # # - pubchem_synonyms = ['(11c)methane', '(13c)methane', '02329_FLUKA', '0CB689EE-132E-4559-A597-C79A40192203', '14493-06-2', '150036-83-2', '295477_ALDRICH', '3B4-2254', '463035_ALDRICH', '490210_ALDRICH', '74-82-8', '8006-14-2', 'AC1L18XA', 'AC1Q2825', # # - queue_number = 134211 # - mformula = C1 # - name = [C] # - smiles = [C] # - csmiles = [C] # - InChI = InChI=1S/C # - InChIKey = OKTJSMMVPCPJKN-UHFFFAOYSA-N # - pubchem_cid = 297 # - pubchem_smiles = C # - pubchem_iupac = methane # - pubchem_synonym0 = (11c)methane # - theory = pspw # - pspw4 = False # - paw = False # - xc = pbe0 # - basis = default # - basisHZ = default # - theory_property = pspw # - property_pspw4 = False # - property_paw = False # - xc_property = pbe0 # - basis_property = default # - basisHZ_property = default # - type = ovb # - solvation_type = # - charge = 0 # - mult = 3 # - babel gen. xyz = True # - cactus gen. xyz = False # - bonds rotated = False # - machine = Shirky # - emailresults = # # - twirl webpage = TwirlMol Link # - image webpage = GIF Image Link # - nmrdb webpage = 1H NMR prediction # - nmrdb webpage = 13C NMR prediction # - nmrdb webpage = COSY prediction # - nmrdb webpage = HSQC/HMBC prediction # # # # # # # # # # # # # # # # # # # # # # # # # # # # .. # C # # # # # title "swnc: ovb theory=pspw xc=pbe0 formula=C1 charge=0 mult=3" #machinejob:Shirky #vtag= osmiles:[C]:osmiles echo start pspw-pbe0-134211 memory 1900 mb charge 0 geometry units angstroms print xyz noautosym C 0.000000 -0.000000 0.000000 end nwpw cutoff 50.0 mult 3 xc pbe0 lmbfgs end nwpw simulation_cell fcc 38.0 end end driver; default; maxiter 50; clear; end task pspw optimize ignore task pspw freq numerical ### Generating HOMO and LUMO Gaussian cube files ### nwpw virtual 8 dplot orbital 1 homo-alpha.cube orbital 123456789 homo-beta.cube density total density.cube end end task pspw energy ignore task pspw pspw_dplot nwpw dplot vectors pspw-pbe0-134211.emovecs orbital 1 lumo-alpha.cube orbital 123456789 lumo-beta.cube end end task pspw pspw_dplot ######################### END NWCHEM INPUT DECK - NWJOB 134211 ######################## ================================================================================ Northwest Computational Chemistry Package (NWChem) 7.0.1 -------------------------------------------------------- Environmental Molecular Sciences Laboratory Pacific Northwest National Laboratory Richland, WA 99352 Copyright (c) 1994-2020 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 = arrow13 program = /home/bylaska/bin/nwchem date = Sat Nov 20 11:13:06 2021 compiled = Wed_Sep_01_21:43:00_2021 source = /home/bylaska/nwchem-releases/nwchem nwchem branch = 7.0.0 nwchem revision = nwchem_on_git-2849-g9a2f8ae ga revision = 5.8.0 use scalapack = F input = /home/bylaska/Projects/Work/RUNARROWS0/nwchemarrows-2021-11-20-11-13-134211.nw prefix = pspw-pbe0-134211. data base = /home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.db status = startup nproc = 32 time left = -1s Memory information ------------------ heap = 62259198 doubles = 475.0 Mbytes stack = 62259195 doubles = 475.0 Mbytes global = 124518400 doubles = 950.0 Mbytes (distinct from heap & stack) total = 249036793 doubles = 1900.0 Mbytes verify = yes hardfail = no Directory information --------------------- 0 permanent = /home/bylaska/Projects/Work/RUNARROWS0 0 scratch = /home/bylaska/Projects/Work/RUNARROWS0 NWChem Input Module ------------------- swnc: ovb theory=pspw xc=pbe0 formula=C1 charge=0 mult=3 -------------------------------------------------------- Scaling coordinates for geometry "geometry" by 1.889725989 (inverse scale = 0.529177249) Geometry "geometry" -> "" ------------------------- Output coordinates in angstroms (scale by 1.889725989 to convert to a.u.) No. Tag Charge X Y Z ---- ---------------- ---------- -------------- -------------- -------------- 1 C 6.0000 0.00000000 0.00000000 0.00000000 Atomic Mass ----------- C 12.000000 Effective nuclear repulsion energy (a.u.) 0.0000000000 Nuclear Dipole moment (a.u.) ---------------------------- X Y Z ---------------- ---------------- ---------------- 0.0000000000 0.0000000000 0.0000000000 XYZ format geometry ------------------- 1 geometry C 0.00000000 0.00000000 0.00000000 Deleted DRIVER restart files NWChem Geometry Optimization ---------------------------- swnc: ovb theory=pspw xc=pbe0 formula=C1 charge=0 mult=3 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) = 1.0D-07 maximum number of steps (nptopt) = 50 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 = ------------------- Energy Minimization ------------------- Using diagonal initial Hessian -------- 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 C 6.0000 0.00000000 0.00000000 0.00000000 Atomic Mass ----------- C 12.000000 Effective nuclear repulsion energy (a.u.) 0.0000000000 Nuclear Dipole moment (a.u.) ---------------------------- X Y Z ---------------- ---------------- ---------------- 0.0000000000 0.0000000000 0.0000000000 **************************************************** * * * NWPW PSPW Calculation * * * * [ (Grassmann/Stiefel manifold implementation) ] * * * * [ NorthWest Chemistry implementation ] * * * * version #5.10 06/12/02 * * * * This code was developed by Eric J. Bylaska, * * and was based upon algorithms and code * * developed by the group of Prof. John H. Weare * * * **************************************************** >>> JOB STARTED AT Sat Nov 20 11:13:10 2021 <<< ================ input data ======================== library name resolved from: compiled reference NWCHEM_NWPW_LIBRARY set to: Generating 1d pseudopotential for C Generated formatted_filename: /home/bylaska/Projects/Work/RUNARROWS0/C.vpp random planewave guess, initial psi:pspw-pbe0-134211.movecs - spin, nalpha, nbeta: 2 3 1 input psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs Warning - K.S. orbitals are not orthonormal. Applying Gram-Schmidt orthonormalization. - exact norm=0.400000E+01 norm=0.288088E+01 corrected norm=0.400000E+01 (error=0.111912E+01) number of processors used: 32 processor grid : 32 x 1 parallel mapping :2d-hilbert parallel mapping : balanced number of threads : 1 parallel io : off options: boundary conditions = periodic (version3) electron spin = unrestricted exchange-correlation = PBE0 (White and Bird) parameterization - HFX relaxed - HFX alpha orbitals: 1 2 3 - HFX beta orbitals : 4 - HFX screened coulomb solver - HFX screening radius(pspw:HFX_screening_radius): 0.800E+01 - HFX screening power (pspw:HFX_screening_power) : 0.800E+01 - HFX screening type (pspw:HFX_screening_type) : 0 - HFX scaling parameter (pspw:HFX_parameter) : 0.250E+00 elements involved in the cluster: 1: C valence charge: 4.0000 lmax= 2 comment : Parameterized (J.Phys.Chem., vol 100, page 6966) Hamman psp pseudpotential type : 0 highest angular component : 2 local potential used : 2 number of non-local projections: 4 cutoff = 0.800 0.850 0.850 total charge: 0.000 atomic composition: C : 1 number of electrons: spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) number of orbitals : spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) supercell: cell_name: cell_default lattice: a1=< 19.000 19.000 0.000 > a2=< 19.000 0.000 19.000 > a3=< 0.000 19.000 19.000 > reciprocal: b1=< 0.165 0.165 -0.165 > b2=< 0.165 -0.165 0.165 > b3=< -0.165 0.165 0.165 > lattice: a= 26.870 b= 26.870 c= 26.870 alpha= 60.000 beta= 60.000 gamma= 60.000 omega= 13718.0 density cutoff= 25.000 fft= 50x 50x 50( 34894 waves 1090 per task) wavefnc cutoff= 12.500 fft= 50x 50x 50( 14406 waves 450 per task) Ewald summation: cut radius= 8.55 and 1 Madelung Wigner-Seitz= 1.79174723 (alpha= 2.88828212 rs=14.85026019) technical parameters: time step= 5.80 fictitious mass= 400000.0 tolerance=0.100E-06 (energy) 0.100E-06 (density) maximum iterations = 1000 ( 10 inner 100 outer ) == Energy Calculation == ============ Grassmann lmbfgs iteration ============ >>> ITERATION STARTED AT Sat Nov 20 11:13:12 2021 <<< iter. Energy DeltaE DeltaRho ------------------------------------------------------ - 15 steepest descent iterations performed 10 -0.5001666091E+01 -0.76797E-01 0.29358E+00 - 10 steepest descent iterations performed 20 -0.5358782898E+01 -0.12966E-01 0.21530E-01 - 10 steepest descent iterations performed 30 -0.5395961871E+01 -0.74937E-04 0.27991E-03 40 -0.5396070335E+01 -0.13233E-05 0.17309E-05 50 -0.5396267981E+01 -0.71840E-04 0.31828E-02 - 10 steepest descent iterations performed 60 -0.5396450490E+01 -0.16331E-05 0.15712E-05 70 -0.5396490299E+01 -0.12567E-04 0.13119E-03 - 10 steepest descent iterations performed 80 -0.5396528269E+01 -0.13799E-05 0.17561E-05 90 -0.5396604511E+01 -0.15326E-04 0.55707E-03 - 10 steepest descent iterations performed 100 -0.5396671644E+01 -0.25814E-05 0.21652E-05 110 -0.5396686981E+01 -0.17811E-05 0.11044E-04 120 -0.5396927841E+01 -0.51167E-04 0.68551E-02 - 10 steepest descent iterations performed 130 -0.5397024062E+01 -0.89297E-06 0.48435E-06 140 -0.5397026730E+01 -0.18322E-06 0.51744E-06 150 -0.5397030697E+01 -0.82759E-06 0.13843E-04 - 10 steepest descent iterations performed 160 -0.5397033960E+01 -0.13868E-06 0.16091E-06 170 -0.5397035552E+01 -0.24913E-06 0.23632E-05 - 10 steepest descent iterations performed 180 -0.5397036973E+01 -0.10897E-06 0.26676E-06 190 -0.5397037197E+01 -0.99300E-07 0.41008E-07 *** tolerance ok. iteration terminated >>> ITERATION ENDED AT Sat Nov 20 11:13:35 2021 <<< == Summary Of Results == number of electrons: spin up= 3.00000 down= 1.00000 (real space) total energy : -0.5397037197E+01 ( -0.53970E+01/ion) total orbital energy: -0.1656467820E+01 ( -0.41412E+00/electron) hartree energy : 0.3467809872E+01 ( 0.86695E+00/electron) exc-corr energy : -0.1148985624E+01 ( -0.28725E+00/electron) HF exchange energy : -0.3538776217E+00 ( -0.88469E-01/electron) ion-ion energy : -0.9652341209E+00 ( -0.96523E+00/ion) kinetic (planewave) : 0.3082731002E+01 ( 0.77068E+00/electron) V_local (planewave) : -0.9852490124E+01 ( -0.24631E+01/electron) V_nl (planewave) : 0.3730094187E+00 ( 0.93252E-01/electron) V_Coul (planewave) : 0.6935619744E+01 ( 0.17339E+01/electron) V_xc. (planewave) : -0.1487582618E+01 ( -0.37190E+00/electron) K.S. HFX energy : -0.7077552435E+00 ( -0.17694E+00/electron) Virial Coefficient : -0.1307750685E+01 orbital energies: -0.2761263E+00 ( -7.514eV) -0.2761471E+00 ( -7.514eV) -0.6162685E+00 ( -16.770eV) -0.4879259E+00 ( -13.277eV) Total PSPW energy : -0.5397037197E+01 === Spin Contamination === = 2.0000000000000000 = 2.0042949888431107 == Center of Charge == spin up ( -0.0003, -0.0003, -0.0002 ) spin down ( 0.0008, 0.0009, 0.0010 ) total ( -0.0000, -0.0000, 0.0001 ) ionic ( 0.0000, 0.0000, 0.0000 ) == Molecular Dipole wrt Center of Mass == mu = ( 0.0000, 0.0000, -0.0005 ) au |mu| = 0.0005 au, 0.0013 Debye output psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs == Timing == cputime in seconds prologue : 0.137855E+01 main loop : 0.232981E+02 epilogue : 0.587840E-01 total : 0.247354E+02 cputime/step: 0.464105E-01 ( 502 evalulations, 182 linesearches) Time spent doing total step percent total time : 0.247687E+02 0.493400E-01 100.0 % i/o time : 0.727819E+00 0.144984E-02 2.9 % FFTs : 0.810942E+01 0.161542E-01 32.7 % dot products : 0.535253E+00 0.106624E-02 2.2 % geodesic : 0.743701E+00 0.148148E-02 3.0 % ffm_dgemm : 0.110182E+00 0.219486E-03 0.4 % fmf_dgemm : 0.312348E+00 0.622207E-03 1.3 % mmm_dgemm : 0.144452E-01 0.287753E-04 0.1 % m_diagonalize : 0.202607E-01 0.403599E-04 0.1 % exchange correlation : 0.113470E+02 0.226036E-01 45.8 % local pseudopotentials : 0.128984E-03 0.256941E-06 0.0 % non-local pseudopotentials : 0.203748E+00 0.405873E-03 0.8 % hartree potentials : 0.249357E+00 0.496726E-03 1.0 % ion-ion interaction : 0.542307E-02 0.108029E-04 0.0 % structure factors : 0.580072E-02 0.115552E-04 0.0 % phase factors : 0.691314E-05 0.137712E-07 0.0 % masking and packing : 0.242567E+01 0.483201E-02 9.8 % queue fft : 0.789381E+01 0.157247E-01 31.9 % queue fft (serial) : 0.193594E+01 0.385645E-02 7.8 % queue fft (message passing): 0.586812E+01 0.116895E-01 23.7 % HFX potential : 0.630804E+01 0.125658E-01 25.5 % non-local psp FFM : 0.141020E+00 0.280916E-03 0.6 % non-local psp FMF : 0.323672E-01 0.644766E-04 0.1 % non-local psp FFM A : 0.684035E-02 0.136262E-04 0.0 % non-local psp FFM B : 0.120428E-01 0.239896E-04 0.0 % == HFX Counters == HFX Evalulations = 502 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:13:35 2021 <<< **************************************************** * * * NWPW PSPW Calculation * * * * [ (Grassmann/Stiefel manifold implementation) ] * * * * [ NorthWest Chemistry implementation ] * * * * version #5.10 06/12/02 * * * * This code was developed by Eric J. Bylaska, * * and was based upon algorithms and code * * developed by the group of Prof. John H. Weare * * * **************************************************** >>> JOB STARTED AT Sat Nov 20 11:13:35 2021 <<< ================ input data ======================== pseudopotential is not correctly formatted:C.vpp Generated formatted_filename: /home/bylaska/Projects/Work/RUNARROWS0/C.vpp Error reading psi - bad grid nfft : 50 50 50 ngrid: 72 72 72 Error reading psi - bad grid nfft : 50 50 50 ngrid: 72 72 72 Error reading psi - bad grid Grid is being converted: ------------------------ To turn off automatic grid conversion: set nwpw:psi_nogrid .false. old_filename: pspw-pbe0-134211.movecs new_filename: pspw-pbe0-134211.movecs converting : 50x 50x 50 --> 72x 72x 72 converting .... psi: 1 spin: 1 converting .... psi: 2 spin: 1 converting .... psi: 3 spin: 1 converting .... psi: 1 spin: 2 input psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs number of processors used: 32 processor grid : 32 x 1 parallel mapping :2d-hilbert parallel mapping : balanced number of threads : 1 parallel io : off options: boundary conditions = periodic (version3) electron spin = unrestricted exchange-correlation = PBE0 (White and Bird) parameterization - HFX relaxed - HFX alpha orbitals: 1 2 3 - HFX beta orbitals : 4 - HFX screened coulomb solver - HFX screening radius(pspw:HFX_screening_radius): 0.800E+01 - HFX screening power (pspw:HFX_screening_power) : 0.800E+01 - HFX screening type (pspw:HFX_screening_type) : 0 - HFX scaling parameter (pspw:HFX_parameter) : 0.250E+00 elements involved in the cluster: 1: C valence charge: 4.0000 lmax= 2 comment : Parameterized (J.Phys.Chem., vol 100, page 6966) Hamman psp pseudpotential type : 0 highest angular component : 2 local potential used : 2 number of non-local projections: 4 cutoff = 0.800 0.850 0.850 total charge: 0.000 atomic composition: C : 1 number of electrons: spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) number of orbitals : spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) supercell: cell_name: cell_default lattice: a1=< 19.000 19.000 0.000 > a2=< 19.000 0.000 19.000 > a3=< 0.000 19.000 19.000 > reciprocal: b1=< 0.165 0.165 -0.165 > b2=< 0.165 -0.165 0.165 > b3=< -0.165 0.165 0.165 > lattice: a= 26.870 b= 26.870 c= 26.870 alpha= 60.000 beta= 60.000 gamma= 60.000 omega= 13718.0 density cutoff= 50.000 fft= 72x 72x 72( 101771 waves 3180 per task) wavefnc cutoff= 25.000 fft= 72x 72x 72( 40966 waves 1280 per task) Ewald summation: cut radius= 8.55 and 1 Madelung Wigner-Seitz= 1.79174723 (alpha= 2.88828212 rs=14.85026019) technical parameters: time step= 5.80 fictitious mass= 400000.0 tolerance=0.100E-06 (energy) 0.100E-06 (density) maximum iterations = 1000 ( 10 inner 100 outer ) == Energy Calculation == ============ Grassmann lmbfgs iteration ============ >>> ITERATION STARTED AT Sat Nov 20 11:13:38 2021 <<< iter. Energy DeltaE DeltaRho ------------------------------------------------------ 10 -0.5421180352E+01 -0.77241E-05 0.45017E-03 20 -0.5421196586E+01 -0.14330E-06 0.10528E-06 30 -0.5421196683E+01 -0.96363E-07 0.24421E-10 *** tolerance ok. iteration terminated >>> ITERATION ENDED AT Sat Nov 20 11:13:44 2021 <<< == Summary Of Results == number of electrons: spin up= 3.00000 down= 1.00000 (real space) total energy : -0.5421196683E+01 ( -0.54212E+01/ion) total orbital energy: -0.1651936480E+01 ( -0.41298E+00/electron) hartree energy : 0.3501193253E+01 ( 0.87530E+00/electron) exc-corr energy : -0.1157253137E+01 ( -0.28931E+00/electron) HF exchange energy : -0.3565210285E+00 ( -0.89130E-01/electron) ion-ion energy : -0.9652341209E+00 ( -0.96523E+00/ion) kinetic (planewave) : 0.3183809469E+01 ( 0.79595E+00/electron) V_local (planewave) : -0.9948520910E+01 ( -0.24871E+01/electron) V_nl (planewave) : 0.3213297918E+00 ( 0.80332E-01/electron) V_Coul (planewave) : 0.7002386506E+01 ( 0.17506E+01/electron) V_xc. (planewave) : -0.1497899279E+01 ( -0.37447E+00/electron) K.S. HFX energy : -0.7130420569E+00 ( -0.17826E+00/electron) Virial Coefficient : -0.1294896548E+01 orbital energies: -0.2802694E+00 ( -7.627eV) -0.2802704E+00 ( -7.627eV) -0.6107618E+00 ( -16.620eV) -0.4806348E+00 ( -13.079eV) Total PSPW energy : -0.5421196683E+01 === Spin Contamination === = 2.0000000000000000 = 2.0045854389723328 == Center of Charge == spin up ( -0.0001, -0.0001, -0.0001 ) spin down ( 0.0004, 0.0005, 0.0004 ) total ( 0.0000, 0.0000, 0.0000 ) ionic ( 0.0000, 0.0000, 0.0000 ) == Molecular Dipole wrt Center of Mass == mu = ( -0.0001, -0.0001, -0.0001 ) au |mu| = 0.0002 au, 0.0004 Debye output psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs == Timing == cputime in seconds prologue : 0.310912E+01 main loop : 0.610300E+01 epilogue : 0.124020E+00 total : 0.933614E+01 cputime/step: 0.127146E+00 ( 48 evalulations, 21 linesearches) Time spent doing total step percent total time : 0.936936E+01 0.195195E+00 100.0 % i/o time : 0.836695E+00 0.174311E-01 8.9 % FFTs : 0.183984E+01 0.383299E-01 19.6 % dot products : 0.959685E-01 0.199934E-02 1.0 % geodesic : 0.155963E+00 0.324923E-02 1.7 % ffm_dgemm : 0.147412E-01 0.307109E-03 0.2 % fmf_dgemm : 0.944006E-01 0.196668E-02 1.0 % mmm_dgemm : 0.137234E-02 0.285904E-04 0.0 % m_diagonalize : 0.230049E-02 0.479269E-04 0.0 % exchange correlation : 0.272009E+01 0.566684E-01 29.0 % local pseudopotentials : 0.212908E-03 0.443558E-05 0.0 % non-local pseudopotentials : 0.428238E-01 0.892162E-03 0.5 % hartree potentials : 0.721615E-01 0.150336E-02 0.8 % ion-ion interaction : 0.111699E-02 0.232706E-04 0.0 % structure factors : 0.185441E-02 0.386335E-04 0.0 % phase factors : 0.595946E-05 0.124156E-06 0.0 % masking and packing : 0.400161E+00 0.833668E-02 4.3 % queue fft : 0.243487E+01 0.507264E-01 26.0 % queue fft (serial) : 0.554765E+00 0.115576E-01 5.9 % queue fft (message passing): 0.185157E+01 0.385744E-01 19.8 % HFX potential : 0.196438E+01 0.409245E-01 21.0 % non-local psp FFM : 0.297883E-01 0.620589E-03 0.3 % non-local psp FMF : 0.895238E-02 0.186508E-03 0.1 % non-local psp FFM A : 0.246832E-02 0.514234E-04 0.0 % non-local psp FFM B : 0.302697E-02 0.630618E-04 0.0 % == HFX Counters == HFX Evalulations = 48 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:13:44 2021 <<< **************************************************** * * * NWPW PSPW Calculation * * * * [ (Grassmann/Stiefel manifold implementation) ] * * * * [ NorthWest Chemistry implementation ] * * * * version #5.10 06/12/02 * * * * This code was developed by Eric J. Bylaska, * * and was based upon algorithms and code * * developed by the group of Prof. John H. Weare * * * **************************************************** >>> JOB STARTED AT Sat Nov 20 11:13:44 2021 <<< ================ input data ======================== pseudopotential is not correctly formatted:C.vpp Generated formatted_filename: /home/bylaska/Projects/Work/RUNARROWS0/C.vpp Error reading psi - bad grid nfft : 72 72 72 ngrid: 90 90 90 Error reading psi - bad grid nfft : 72 72 72 ngrid: 90 90 90 Error reading psi - bad grid Grid is being converted: ------------------------ To turn off automatic grid conversion: set nwpw:psi_nogrid .false. old_filename: pspw-pbe0-134211.movecs new_filename: pspw-pbe0-134211.movecs converting : 72x 72x 72 --> 90x 90x 90 converting .... psi: 1 spin: 1 converting .... psi: 2 spin: 1 converting .... psi: 3 spin: 1 converting .... psi: 1 spin: 2 input psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs number of processors used: 32 processor grid : 32 x 1 parallel mapping :2d-hilbert parallel mapping : balanced number of threads : 1 parallel io : off options: boundary conditions = periodic (version3) electron spin = unrestricted exchange-correlation = PBE0 (White and Bird) parameterization - HFX relaxed - HFX alpha orbitals: 1 2 3 - HFX beta orbitals : 4 - HFX screened coulomb solver - HFX screening radius(pspw:HFX_screening_radius): 0.800E+01 - HFX screening power (pspw:HFX_screening_power) : 0.800E+01 - HFX screening type (pspw:HFX_screening_type) : 0 - HFX scaling parameter (pspw:HFX_parameter) : 0.250E+00 elements involved in the cluster: 1: C valence charge: 4.0000 lmax= 2 comment : Parameterized (J.Phys.Chem., vol 100, page 6966) Hamman psp pseudpotential type : 0 highest angular component : 2 local potential used : 2 number of non-local projections: 4 cutoff = 0.800 0.850 0.850 total charge: 0.000 atomic composition: C : 1 number of electrons: spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) number of orbitals : spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) supercell: cell_name: cell_default lattice: a1=< 19.000 19.000 0.000 > a2=< 19.000 0.000 19.000 > a3=< 0.000 19.000 19.000 > reciprocal: b1=< 0.165 0.165 -0.165 > b2=< 0.165 -0.165 0.165 > b3=< -0.165 0.165 0.165 > lattice: a= 26.870 b= 26.870 c= 26.870 alpha= 60.000 beta= 60.000 gamma= 60.000 omega= 13718.0 density cutoff= 75.000 fft= 90x 90x 90( 192337 waves 6010 per task) wavefnc cutoff= 37.500 fft= 90x 90x 90( 75079 waves 2346 per task) Ewald summation: cut radius= 8.55 and 1 Madelung Wigner-Seitz= 1.79174723 (alpha= 2.88828212 rs=14.85026019) technical parameters: time step= 5.80 fictitious mass= 400000.0 tolerance=0.100E-06 (energy) 0.100E-06 (density) maximum iterations = 1000 ( 10 inner 100 outer ) == Energy Calculation == ============ Grassmann lmbfgs iteration ============ >>> ITERATION STARTED AT Sat Nov 20 11:13:49 2021 <<< iter. Energy DeltaE DeltaRho ------------------------------------------------------ 10 -0.5424239426E+01 -0.10623E-06 0.23648E-04 20 -0.5424239505E+01 -0.78997E-07 0.81399E-09 *** tolerance ok. iteration terminated >>> ITERATION ENDED AT Sat Nov 20 11:13:54 2021 <<< == Summary Of Results == number of electrons: spin up= 3.00000 down= 1.00000 (real space) total energy : -0.5424239505E+01 ( -0.54242E+01/ion) total orbital energy: -0.1652276646E+01 ( -0.41307E+00/electron) hartree energy : 0.3504490698E+01 ( 0.87612E+00/electron) exc-corr energy : -0.1158267143E+01 ( -0.28957E+00/electron) HF exchange energy : -0.3568035381E+00 ( -0.89201E-01/electron) ion-ion energy : -0.9652341209E+00 ( -0.96523E+00/ion) kinetic (planewave) : 0.3193607923E+01 ( 0.79840E+00/electron) V_local (planewave) : -0.9957912866E+01 ( -0.24895E+01/electron) V_nl (planewave) : 0.3158795420E+00 ( 0.78970E-01/electron) V_Coul (planewave) : 0.7008981396E+01 ( 0.17522E+01/electron) V_xc. (planewave) : -0.1499225565E+01 ( -0.37481E+00/electron) K.S. HFX energy : -0.7136070762E+00 ( -0.17840E+00/electron) Virial Coefficient : -0.1293921356E+01 orbital energies: -0.2807901E+00 ( -7.641eV) -0.2807979E+00 ( -7.641eV) -0.6103835E+00 ( -16.610eV) -0.4803051E+00 ( -13.070eV) Total PSPW energy : -0.5424239505E+01 === Spin Contamination === = 2.0000000000000000 = 2.0045752944359538 == Center of Charge == spin up ( -0.0001, -0.0001, -0.0001 ) spin down ( 0.0004, 0.0004, 0.0004 ) total ( 0.0000, 0.0000, 0.0000 ) ionic ( 0.0000, 0.0000, 0.0000 ) == Molecular Dipole wrt Center of Mass == mu = ( -0.0001, -0.0001, -0.0001 ) au |mu| = 0.0001 au, 0.0003 Debye output psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs == Timing == cputime in seconds prologue : 0.451895E+01 main loop : 0.505579E+01 epilogue : 0.209893E+00 total : 0.978463E+01 cputime/step: 0.187251E+00 ( 27 evalulations, 11 linesearches) Time spent doing total step percent total time : 0.981804E+01 0.363631E+00 100.0 % i/o time : 0.154948E+01 0.573882E-01 15.8 % FFTs : 0.148999E+01 0.551847E-01 15.2 % dot products : 0.593508E-01 0.219818E-02 0.6 % geodesic : 0.121346E+00 0.449429E-02 1.2 % ffm_dgemm : 0.893691E-02 0.330997E-03 0.1 % fmf_dgemm : 0.930833E-01 0.344753E-02 0.9 % mmm_dgemm : 0.756502E-03 0.280186E-04 0.0 % m_diagonalize : 0.128078E-02 0.474361E-04 0.0 % exchange correlation : 0.251847E+01 0.932768E-01 25.7 % local pseudopotentials : 0.353813E-03 0.131042E-04 0.0 % non-local pseudopotentials : 0.365186E-01 0.135254E-02 0.4 % hartree potentials : 0.108395E+00 0.401463E-02 1.1 % ion-ion interaction : 0.136685E-02 0.506242E-04 0.0 % structure factors : 0.230312E-02 0.853008E-04 0.0 % phase factors : 0.691414E-05 0.256079E-06 0.0 % masking and packing : 0.387169E+00 0.143396E-01 3.9 % queue fft : 0.164149E+01 0.607958E-01 16.7 % queue fft (serial) : 0.607998E+00 0.225185E-01 6.2 % queue fft (message passing): 0.986414E+00 0.365338E-01 10.0 % HFX potential : 0.139751E+01 0.517597E-01 14.2 % non-local psp FFM : 0.243190E-01 0.900702E-03 0.2 % non-local psp FMF : 0.845909E-02 0.313300E-03 0.1 % non-local psp FFM A : 0.344562E-02 0.127615E-03 0.0 % non-local psp FFM B : 0.296427E-02 0.109788E-03 0.0 % == HFX Counters == HFX Evalulations = 27 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:13:54 2021 <<< **************************************************** * * * NWPW PSPW Calculation * * * * [ (Grassmann/Stiefel manifold implementation) ] * * * * [ NorthWest Chemistry implementation ] * * * * version #5.10 06/12/02 * * * * This code was developed by Eric J. Bylaska, * * and was based upon algorithms and code * * developed by the group of Prof. John H. Weare * * * **************************************************** >>> JOB STARTED AT Sat Nov 20 11:13:54 2021 <<< ================ input data ======================== pseudopotential is not correctly formatted:C.vpp Generated formatted_filename: /home/bylaska/Projects/Work/RUNARROWS0/C.vpp Error reading psi - bad grid nfft : 90 90 90 ngrid: 100 100 100 Error reading psi - bad grid nfft : 90 90 90 ngrid: 100 100 100 Error reading psi - bad grid Grid is being converted: ------------------------ To turn off automatic grid conversion: set nwpw:psi_nogrid .false. old_filename: pspw-pbe0-134211.movecs new_filename: pspw-pbe0-134211.movecs converting : 90x 90x 90 --> 100x100x100 converting .... psi: 1 spin: 1 converting .... psi: 2 spin: 1 converting .... psi: 3 spin: 1 converting .... psi: 1 spin: 2 input psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs number of processors used: 32 processor grid : 32 x 1 parallel mapping :2d-hilbert parallel mapping : balanced number of threads : 1 parallel io : off options: boundary conditions = periodic (version3) electron spin = unrestricted exchange-correlation = PBE0 (White and Bird) parameterization - HFX relaxed - HFX alpha orbitals: 1 2 3 - HFX beta orbitals : 4 - HFX screened coulomb solver - HFX screening radius(pspw:HFX_screening_radius): 0.800E+01 - HFX screening power (pspw:HFX_screening_power) : 0.800E+01 - HFX screening type (pspw:HFX_screening_type) : 0 - HFX scaling parameter (pspw:HFX_parameter) : 0.250E+00 elements involved in the cluster: 1: C valence charge: 4.0000 lmax= 2 comment : Parameterized (J.Phys.Chem., vol 100, page 6966) Hamman psp pseudpotential type : 0 highest angular component : 2 local potential used : 2 number of non-local projections: 4 cutoff = 0.800 0.850 0.850 total charge: 0.000 atomic composition: C : 1 number of electrons: spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) number of orbitals : spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) supercell: cell_name: cell_default lattice: a1=< 19.000 19.000 0.000 > a2=< 19.000 0.000 19.000 > a3=< 0.000 19.000 19.000 > reciprocal: b1=< 0.165 0.165 -0.165 > b2=< 0.165 -0.165 0.165 > b3=< -0.165 0.165 0.165 > lattice: a= 26.870 b= 26.870 c= 26.870 alpha= 60.000 beta= 60.000 gamma= 60.000 omega= 13718.0 density cutoff=100.000 fft=100x100x100( 282860 waves 8839 per task) wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 3616 per task) Ewald summation: cut radius= 8.55 and 1 Madelung Wigner-Seitz= 1.79174723 (alpha= 2.88828212 rs=14.85026019) technical parameters: time step= 5.80 fictitious mass= 400000.0 tolerance=0.100E-06 (energy) 0.100E-06 (density) maximum iterations = 1000 ( 10 inner 100 outer ) == Energy Calculation == ============ Grassmann lmbfgs iteration ============ >>> ITERATION STARTED AT Sat Nov 20 11:14:00 2021 <<< iter. Energy DeltaE DeltaRho ------------------------------------------------------ 10 -0.5424770381E+01 -0.36270E-07 0.22844E-05 20 -0.5424770404E+01 -0.23083E-07 0.15339E-09 *** tolerance ok. iteration terminated >>> ITERATION ENDED AT Sat Nov 20 11:14:07 2021 <<< == Summary Of Results == number of electrons: spin up= 3.00000 down= 1.00000 (real space) total energy : -0.5424770404E+01 ( -0.54248E+01/ion) total orbital energy: -0.1652099566E+01 ( -0.41302E+00/electron) hartree energy : 0.3505362470E+01 ( 0.87634E+00/electron) exc-corr energy : -0.1158538277E+01 ( -0.28963E+00/electron) HF exchange energy : -0.3568837030E+00 ( -0.89221E-01/electron) ion-ion energy : -0.9652341209E+00 ( -0.96523E+00/ion) kinetic (planewave) : 0.3196306508E+01 ( 0.79908E+00/electron) V_local (planewave) : -0.9960290399E+01 ( -0.24901E+01/electron) V_nl (planewave) : 0.3145071174E+00 ( 0.78627E-01/electron) V_Coul (planewave) : 0.7010724940E+01 ( 0.17527E+01/electron) V_xc. (planewave) : -0.1499580326E+01 ( -0.37490E+00/electron) K.S. HFX energy : -0.7137674061E+00 ( -0.17844E+00/electron) Virial Coefficient : -0.1293567640E+01 orbital energies: -0.2808744E+00 ( -7.643eV) -0.2808803E+00 ( -7.643eV) -0.6102180E+00 ( -16.605eV) -0.4801269E+00 ( -13.065eV) Total PSPW energy : -0.5424770404E+01 === Spin Contamination === = 2.0000000000000000 = 2.0045822031296376 == Center of Charge == spin up ( -0.0001, -0.0001, -0.0001 ) spin down ( 0.0003, 0.0003, 0.0003 ) total ( 0.0000, 0.0000, 0.0000 ) ionic ( 0.0000, 0.0000, 0.0000 ) == Molecular Dipole wrt Center of Mass == mu = ( -0.0000, -0.0001, -0.0001 ) au |mu| = 0.0001 au, 0.0002 Debye Translation force removed: ( -0.00000 -0.00000 -0.00000) ============= Ion Gradients ================= Ion Forces: 1 C ( 0.000000 0.000000 0.000000 ) C.O.M. ( 0.000000 0.000000 0.000000 ) =============================================== |F| = 0.000000E+00 |F|/nion = 0.000000E+00 max|Fatom|= 0.000000E+00 ( 0.000eV/Angstrom) output psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs == Timing == cputime in seconds prologue : 0.544189E+01 main loop : 0.699241E+01 epilogue : 0.276857E+00 total : 0.127112E+02 cputime/step: 0.258978E+00 ( 27 evalulations, 11 linesearches) Time spent doing total step percent total time : 0.127459E+02 0.472070E+00 100.0 % i/o time : 0.199101E+01 0.737410E-01 15.6 % FFTs : 0.198610E+01 0.735592E-01 15.6 % dot products : 0.891641E-01 0.330237E-02 0.7 % geodesic : 0.181052E+00 0.670562E-02 1.4 % ffm_dgemm : 0.130902E-01 0.484821E-03 0.1 % fmf_dgemm : 0.147182E+00 0.545117E-02 1.2 % mmm_dgemm : 0.834703E-03 0.309149E-04 0.0 % m_diagonalize : 0.126385E-02 0.468092E-04 0.0 % exchange correlation : 0.345522E+01 0.127971E+00 27.1 % local pseudopotentials : 0.140691E-02 0.521077E-04 0.0 % non-local pseudopotentials : 0.586622E-01 0.217267E-02 0.5 % hartree potentials : 0.153229E+00 0.567514E-02 1.2 % ion-ion interaction : 0.228214E-02 0.845238E-04 0.0 % structure factors : 0.419497E-02 0.155369E-03 0.0 % phase factors : 0.905891E-05 0.335515E-06 0.0 % masking and packing : 0.599690E+00 0.222107E-01 4.7 % queue fft : 0.228392E+01 0.845897E-01 17.9 % queue fft (serial) : 0.816878E+00 0.302547E-01 6.4 % queue fft (message passing): 0.138830E+01 0.514185E-01 10.9 % HFX potential : 0.194079E+01 0.718812E-01 15.2 % non-local psp FFM : 0.386806E-01 0.143261E-02 0.3 % non-local psp FMF : 0.124009E-01 0.459291E-03 0.1 % non-local psp FFM A : 0.622772E-02 0.230656E-03 0.0 % non-local psp FFM B : 0.491858E-02 0.182170E-03 0.0 % == HFX Counters == HFX Evalulations = 27 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:14:07 2021 <<< @ Step Energy Delta E Gmax Grms Xrms Xmax Walltime @ ---- ---------------- -------- -------- -------- -------- -------- -------- @ 0 -5.42477040 0.0D+00 0.00000 0.00000 0.00000 0.00000 61.6 ok ok ok ok Warning ... line search gradient +ve 1.0000000000000000 0.0000000000000000 **************************************************** * * * NWPW PSPW Calculation * * * * [ (Grassmann/Stiefel manifold implementation) ] * * * * [ NorthWest Chemistry implementation ] * * * * version #5.10 06/12/02 * * * * This code was developed by Eric J. Bylaska, * * and was based upon algorithms and code * * developed by the group of Prof. John H. Weare * * * **************************************************** >>> JOB STARTED AT Sat Nov 20 11:14:07 2021 <<< ================ input data ======================== input psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs number of processors used: 32 processor grid : 32 x 1 parallel mapping :2d-hilbert parallel mapping : balanced number of threads : 1 parallel io : off options: boundary conditions = periodic (version3) electron spin = unrestricted exchange-correlation = PBE0 (White and Bird) parameterization - HFX relaxed - HFX alpha orbitals: 1 2 3 - HFX beta orbitals : 4 - HFX screened coulomb solver - HFX screening radius(pspw:HFX_screening_radius): 0.800E+01 - HFX screening power (pspw:HFX_screening_power) : 0.800E+01 - HFX screening type (pspw:HFX_screening_type) : 0 - HFX scaling parameter (pspw:HFX_parameter) : 0.250E+00 elements involved in the cluster: 1: C valence charge: 4.0000 lmax= 2 comment : Parameterized (J.Phys.Chem., vol 100, page 6966) Hamman psp pseudpotential type : 0 highest angular component : 2 local potential used : 2 number of non-local projections: 4 cutoff = 0.800 0.850 0.850 total charge: 0.000 atomic composition: C : 1 number of electrons: spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) number of orbitals : spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) supercell: cell_name: cell_default lattice: a1=< 19.000 19.000 0.000 > a2=< 19.000 0.000 19.000 > a3=< 0.000 19.000 19.000 > reciprocal: b1=< 0.165 0.165 -0.165 > b2=< 0.165 -0.165 0.165 > b3=< -0.165 0.165 0.165 > lattice: a= 26.870 b= 26.870 c= 26.870 alpha= 60.000 beta= 60.000 gamma= 60.000 omega= 13718.0 density cutoff=100.000 fft=100x100x100( 282860 waves 8839 per task) wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 3616 per task) Ewald summation: cut radius= 8.55 and 1 Madelung Wigner-Seitz= 1.79174723 (alpha= 2.88828212 rs=14.85026019) technical parameters: time step= 5.80 fictitious mass= 400000.0 tolerance=0.100E-06 (energy) 0.100E-06 (density) maximum iterations = 1000 ( 10 inner 100 outer ) == Energy Calculation == ============ Grassmann lmbfgs iteration ============ >>> ITERATION STARTED AT Sat Nov 20 11:14:11 2021 <<< iter. Energy DeltaE DeltaRho ------------------------------------------------------ 10 -0.5424770414E+01 -0.93663E-08 0.12345E-10 *** tolerance ok. iteration terminated >>> ITERATION ENDED AT Sat Nov 20 11:14:12 2021 <<< == Summary Of Results == number of electrons: spin up= 3.00000 down= 1.00000 (real space) total energy : -0.5424770414E+01 ( -0.54248E+01/ion) total orbital energy: -0.1652102287E+01 ( -0.41303E+00/electron) hartree energy : 0.3505359397E+01 ( 0.87634E+00/electron) exc-corr energy : -0.1158537600E+01 ( -0.28963E+00/electron) HF exchange energy : -0.3568835187E+00 ( -0.89221E-01/electron) ion-ion energy : -0.9652341209E+00 ( -0.96523E+00/ion) kinetic (planewave) : 0.3196296086E+01 ( 0.79907E+00/electron) V_local (planewave) : -0.9960279222E+01 ( -0.24901E+01/electron) V_nl (planewave) : 0.3145085647E+00 ( 0.78627E-01/electron) V_Coul (planewave) : 0.7010718794E+01 ( 0.17527E+01/electron) V_xc. (planewave) : -0.1499579473E+01 ( -0.37489E+00/electron) K.S. HFX energy : -0.7137670375E+00 ( -0.17844E+00/electron) Virial Coefficient : -0.1293569564E+01 orbital energies: -0.2808750E+00 ( -7.643eV) -0.2808809E+00 ( -7.643eV) -0.6102188E+00 ( -16.605eV) -0.4801275E+00 ( -13.065eV) Total PSPW energy : -0.5424770414E+01 === Spin Contamination === = 2.0000000000000000 = 2.0045823741933062 == Center of Charge == spin up ( -0.0001, -0.0001, -0.0001 ) spin down ( 0.0003, 0.0003, 0.0003 ) total ( 0.0000, 0.0000, 0.0000 ) ionic ( 0.0000, 0.0000, 0.0000 ) == Molecular Dipole wrt Center of Mass == mu = ( -0.0000, -0.0001, -0.0001 ) au |mu| = 0.0001 au, 0.0002 Debye output psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs == Timing == cputime in seconds prologue : 0.376158E+01 main loop : 0.130750E+01 epilogue : 0.276157E+00 total : 0.534524E+01 cputime/step: 0.261500E+00 ( 5 evalulations, 1 linesearches) Time spent doing total step percent total time : 0.537846E+01 0.107569E+01 100.0 % i/o time : 0.512237E+00 0.102447E+00 9.5 % FFTs : 0.373209E+00 0.746419E-01 6.9 % dot products : 0.140345E-01 0.280690E-02 0.3 % geodesic : 0.161095E-01 0.322189E-02 0.3 % ffm_dgemm : 0.152397E-02 0.304795E-03 0.0 % fmf_dgemm : 0.140948E-01 0.281896E-02 0.3 % mmm_dgemm : 0.967979E-04 0.193596E-04 0.0 % m_diagonalize : 0.215053E-03 0.430105E-04 0.0 % exchange correlation : 0.636979E+00 0.127396E+00 11.8 % local pseudopotentials : 0.472069E-03 0.944138E-04 0.0 % non-local pseudopotentials : 0.105569E-01 0.211139E-02 0.2 % hartree potentials : 0.972211E-01 0.194442E-01 1.8 % ion-ion interaction : 0.783920E-03 0.156784E-03 0.0 % structure factors : 0.939846E-03 0.187969E-03 0.0 % phase factors : 0.905891E-05 0.181178E-05 0.0 % masking and packing : 0.114644E+00 0.229288E-01 2.1 % queue fft : 0.419447E+00 0.838893E-01 7.8 % queue fft (serial) : 0.152171E+00 0.304341E-01 2.8 % queue fft (message passing): 0.252774E+00 0.505548E-01 4.7 % HFX potential : 0.356195E+00 0.712390E-01 6.6 % non-local psp FFM : 0.724292E-02 0.144858E-02 0.1 % non-local psp FMF : 0.225830E-02 0.451660E-03 0.0 % non-local psp FFM A : 0.113296E-02 0.226593E-03 0.0 % non-local psp FFM B : 0.866891E-03 0.173378E-03 0.0 % == HFX Counters == HFX Evalulations = 5 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:14:12 2021 <<< Line search: step=-1.00 grad= 0.0D+00 hess=-9.4D-09 energy= -5.424770 mode=accept new step=-1.00 predicted energy= -5.424770 -------- 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 C 6.0000 0.00000000 0.00000000 0.00000000 Atomic Mass ----------- C 12.000000 Effective nuclear repulsion energy (a.u.) 0.0000000000 Nuclear Dipole moment (a.u.) ---------------------------- X Y Z ---------------- ---------------- ---------------- 0.0000000000 0.0000000000 0.0000000000 **************************************************** * * * NWPW PSPW Calculation * * * * [ (Grassmann/Stiefel manifold implementation) ] * * * * [ NorthWest Chemistry implementation ] * * * * version #5.10 06/12/02 * * * * This code was developed by Eric J. Bylaska, * * and was based upon algorithms and code * * developed by the group of Prof. John H. Weare * * * **************************************************** >>> JOB STARTED AT Sat Nov 20 11:14:12 2021 <<< ================ input data ======================== input psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs number of processors used: 32 processor grid : 32 x 1 parallel mapping :2d-hilbert parallel mapping : balanced number of threads : 1 parallel io : off options: boundary conditions = periodic (version3) electron spin = unrestricted exchange-correlation = PBE0 (White and Bird) parameterization - HFX relaxed - HFX alpha orbitals: 1 2 3 - HFX beta orbitals : 4 - HFX screened coulomb solver - HFX screening radius(pspw:HFX_screening_radius): 0.800E+01 - HFX screening power (pspw:HFX_screening_power) : 0.800E+01 - HFX screening type (pspw:HFX_screening_type) : 0 - HFX scaling parameter (pspw:HFX_parameter) : 0.250E+00 elements involved in the cluster: 1: C valence charge: 4.0000 lmax= 2 comment : Parameterized (J.Phys.Chem., vol 100, page 6966) Hamman psp pseudpotential type : 0 highest angular component : 2 local potential used : 2 number of non-local projections: 4 cutoff = 0.800 0.850 0.850 total charge: 0.000 atomic composition: C : 1 number of electrons: spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) number of orbitals : spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) supercell: cell_name: cell_default lattice: a1=< 19.000 19.000 0.000 > a2=< 19.000 0.000 19.000 > a3=< 0.000 19.000 19.000 > reciprocal: b1=< 0.165 0.165 -0.165 > b2=< 0.165 -0.165 0.165 > b3=< -0.165 0.165 0.165 > lattice: a= 26.870 b= 26.870 c= 26.870 alpha= 60.000 beta= 60.000 gamma= 60.000 omega= 13718.0 density cutoff=100.000 fft=100x100x100( 282860 waves 8839 per task) wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 3616 per task) Ewald summation: cut radius= 8.55 and 1 Madelung Wigner-Seitz= 1.79174723 (alpha= 2.88828212 rs=14.85026019) technical parameters: time step= 5.80 fictitious mass= 400000.0 tolerance=0.100E-06 (energy) 0.100E-06 (density) maximum iterations = 1000 ( 10 inner 100 outer ) == Energy Calculation == ============ Grassmann lmbfgs iteration ============ >>> ITERATION STARTED AT Sat Nov 20 11:14:16 2021 <<< iter. Energy DeltaE DeltaRho ------------------------------------------------------ 10 -0.5424770419E+01 -0.48307E-08 0.37640E-11 *** tolerance ok. iteration terminated >>> ITERATION ENDED AT Sat Nov 20 11:14:18 2021 <<< == Summary Of Results == number of electrons: spin up= 3.00000 down= 1.00000 (real space) total energy : -0.5424770419E+01 ( -0.54248E+01/ion) total orbital energy: -0.1652104422E+01 ( -0.41303E+00/electron) hartree energy : 0.3505356986E+01 ( 0.87634E+00/electron) exc-corr energy : -0.1158536985E+01 ( -0.28963E+00/electron) HF exchange energy : -0.3568833781E+00 ( -0.89221E-01/electron) ion-ion energy : -0.9652341209E+00 ( -0.96523E+00/ion) kinetic (planewave) : 0.3196287131E+01 ( 0.79907E+00/electron) V_local (planewave) : -0.9960270773E+01 ( -0.24901E+01/electron) V_nl (planewave) : 0.3145107209E+00 ( 0.78628E-01/electron) V_Coul (planewave) : 0.7010713971E+01 ( 0.17527E+01/electron) V_xc. (planewave) : -0.1499578717E+01 ( -0.37489E+00/electron) K.S. HFX energy : -0.7137667562E+00 ( -0.17844E+00/electron) Virial Coefficient : -0.1293571143E+01 orbital energies: -0.2808755E+00 ( -7.643eV) -0.2808814E+00 ( -7.643eV) -0.6102195E+00 ( -16.605eV) -0.4801280E+00 ( -13.065eV) Total PSPW energy : -0.5424770419E+01 === Spin Contamination === = 2.0000000000000000 = 2.0045825228655865 == Center of Charge == spin up ( -0.0001, -0.0001, -0.0001 ) spin down ( 0.0003, 0.0003, 0.0003 ) total ( 0.0000, 0.0000, 0.0000 ) ionic ( 0.0000, 0.0000, 0.0000 ) == Molecular Dipole wrt Center of Mass == mu = ( -0.0000, -0.0001, -0.0001 ) au |mu| = 0.0001 au, 0.0002 Debye Translation force removed: ( -0.00000 -0.00000 -0.00000) ============= Ion Gradients ================= Ion Forces: 1 C ( 0.000000 0.000000 0.000000 ) C.O.M. ( 0.000000 0.000000 0.000000 ) =============================================== |F| = 0.000000E+00 |F|/nion = 0.000000E+00 max|Fatom|= 0.000000E+00 ( 0.000eV/Angstrom) output psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs == Timing == cputime in seconds prologue : 0.384504E+01 main loop : 0.132655E+01 epilogue : 0.273739E+00 total : 0.544532E+01 cputime/step: 0.265310E+00 ( 5 evalulations, 1 linesearches) Time spent doing total step percent total time : 0.547867E+01 0.109573E+01 100.0 % i/o time : 0.565587E+00 0.113117E+00 10.3 % FFTs : 0.375983E+00 0.751966E-01 6.9 % dot products : 0.168976E-01 0.337953E-02 0.3 % geodesic : 0.164201E-01 0.328403E-02 0.3 % ffm_dgemm : 0.129772E-02 0.259543E-03 0.0 % fmf_dgemm : 0.140109E-01 0.280218E-02 0.3 % mmm_dgemm : 0.891685E-04 0.178337E-04 0.0 % m_diagonalize : 0.250100E-03 0.500200E-04 0.0 % exchange correlation : 0.638562E+00 0.127712E+00 11.7 % local pseudopotentials : 0.134706E-02 0.269413E-03 0.0 % non-local pseudopotentials : 0.140457E-01 0.280914E-02 0.3 % hartree potentials : 0.979433E-01 0.195887E-01 1.8 % ion-ion interaction : 0.163913E-02 0.327826E-03 0.0 % structure factors : 0.154972E-02 0.309944E-03 0.0 % phase factors : 0.905891E-05 0.181178E-05 0.0 % masking and packing : 0.127768E+00 0.255537E-01 2.3 % queue fft : 0.426532E+00 0.853065E-01 7.8 % queue fft (serial) : 0.151115E+00 0.302229E-01 2.8 % queue fft (message passing): 0.260945E+00 0.521890E-01 4.8 % HFX potential : 0.362205E+00 0.724410E-01 6.6 % non-local psp FFM : 0.860238E-02 0.172048E-02 0.2 % non-local psp FMF : 0.224733E-02 0.449467E-03 0.0 % non-local psp FFM A : 0.121212E-02 0.242423E-03 0.0 % non-local psp FFM B : 0.996113E-03 0.199223E-03 0.0 % == HFX Counters == HFX Evalulations = 5 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:14:18 2021 <<< Step Energy Delta E Gmax Grms Xrms Xmax Walltime ---- ---------------- -------- -------- -------- -------- -------- -------- @ 1 -5.42477042 -1.4D-08 0.00000 0.00000 0.00000 0.00000 72.6 ok ok ok ok ---------------------- Optimization converged ---------------------- Step Energy Delta E Gmax Grms Xrms Xmax Walltime ---- ---------------- -------- -------- -------- -------- -------- -------- @ 1 -5.42477042 -1.4D-08 0.00000 0.00000 0.00000 0.00000 72.6 ok ok ok ok Geometry "geometry" -> "geometry" --------------------------------- Output coordinates in angstroms (scale by 1.889725989 to convert to a.u.) No. Tag Charge X Y Z ---- ---------------- ---------- -------------- -------------- -------------- 1 C 6.0000 0.00000000 0.00000000 0.00000000 Atomic Mass ----------- C 12.000000 Effective nuclear repulsion energy (a.u.) 0.0000000000 Nuclear Dipole moment (a.u.) ---------------------------- X Y Z ---------------- ---------------- ---------------- 0.0000000000 0.0000000000 0.0000000000 Task times cpu: 70.5s wall: 71.8s NWChem Input Module ------------------- NWChem Nuclear Hessian and Frequency Analysis --------------------------------------------- NWChem Finite-difference Hessian -------------------------------- **************************************************** * * * NWPW PSPW Calculation * * * * [ (Grassmann/Stiefel manifold implementation) ] * * * * [ NorthWest Chemistry implementation ] * * * * version #5.10 06/12/02 * * * * This code was developed by Eric J. Bylaska, * * and was based upon algorithms and code * * developed by the group of Prof. John H. Weare * * * **************************************************** >>> JOB STARTED AT Sat Nov 20 11:14:18 2021 <<< ================ input data ======================== input psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs number of processors used: 32 processor grid : 32 x 1 parallel mapping :2d-hilbert parallel mapping : balanced number of threads : 1 parallel io : off options: boundary conditions = periodic (version3) electron spin = unrestricted exchange-correlation = PBE0 (White and Bird) parameterization - HFX relaxed - HFX alpha orbitals: 1 2 3 - HFX beta orbitals : 4 - HFX screened coulomb solver - HFX screening radius(pspw:HFX_screening_radius): 0.800E+01 - HFX screening power (pspw:HFX_screening_power) : 0.800E+01 - HFX screening type (pspw:HFX_screening_type) : 0 - HFX scaling parameter (pspw:HFX_parameter) : 0.250E+00 elements involved in the cluster: 1: C valence charge: 4.0000 lmax= 2 comment : Parameterized (J.Phys.Chem., vol 100, page 6966) Hamman psp pseudpotential type : 0 highest angular component : 2 local potential used : 2 number of non-local projections: 4 cutoff = 0.800 0.850 0.850 total charge: 0.000 atomic composition: C : 1 number of electrons: spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) number of orbitals : spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) supercell: cell_name: cell_default lattice: a1=< 19.000 19.000 0.000 > a2=< 19.000 0.000 19.000 > a3=< 0.000 19.000 19.000 > reciprocal: b1=< 0.165 0.165 -0.165 > b2=< 0.165 -0.165 0.165 > b3=< -0.165 0.165 0.165 > lattice: a= 26.870 b= 26.870 c= 26.870 alpha= 60.000 beta= 60.000 gamma= 60.000 omega= 13718.0 density cutoff=100.000 fft=100x100x100( 282860 waves 8839 per task) wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 3616 per task) Ewald summation: cut radius= 8.55 and 1 Madelung Wigner-Seitz= 1.79174723 (alpha= 2.88828212 rs=14.85026019) technical parameters: time step= 5.80 fictitious mass= 400000.0 tolerance=0.100E-06 (energy) 0.100E-06 (density) maximum iterations = 1000 ( 10 inner 100 outer ) == Energy Calculation == ============ Grassmann lmbfgs iteration ============ >>> ITERATION STARTED AT Sat Nov 20 11:14:22 2021 <<< iter. Energy DeltaE DeltaRho ------------------------------------------------------ 10 -0.5424770422E+01 -0.34752E-08 0.31211E-11 *** tolerance ok. iteration terminated >>> ITERATION ENDED AT Sat Nov 20 11:14:23 2021 <<< == Summary Of Results == number of electrons: spin up= 3.00000 down= 1.00000 (real space) total energy : -0.5424770422E+01 ( -0.54248E+01/ion) total orbital energy: -0.1652106483E+01 ( -0.41303E+00/electron) hartree energy : 0.3505354662E+01 ( 0.87634E+00/electron) exc-corr energy : -0.1158536399E+01 ( -0.28963E+00/electron) HF exchange energy : -0.3568832482E+00 ( -0.89221E-01/electron) ion-ion energy : -0.9652341209E+00 ( -0.96523E+00/ion) kinetic (planewave) : 0.3196278536E+01 ( 0.79907E+00/electron) V_local (planewave) : -0.9960262579E+01 ( -0.24901E+01/electron) V_nl (planewave) : 0.3145127276E+00 ( 0.78628E-01/electron) V_Coul (planewave) : 0.7010709324E+01 ( 0.17527E+01/electron) V_xc. (planewave) : -0.1499577994E+01 ( -0.37489E+00/electron) K.S. HFX energy : -0.7137664964E+00 ( -0.17844E+00/electron) Virial Coefficient : -0.1293572658E+01 orbital energies: -0.2808760E+00 ( -7.643eV) -0.2808820E+00 ( -7.643eV) -0.6102201E+00 ( -16.605eV) -0.4801284E+00 ( -13.065eV) Total PSPW energy : -0.5424770422E+01 === Spin Contamination === = 2.0000000000000000 = 2.0045826789685925 == Center of Charge == spin up ( -0.0001, -0.0001, -0.0001 ) spin down ( 0.0003, 0.0003, 0.0003 ) total ( 0.0000, 0.0000, 0.0000 ) ionic ( 0.0000, 0.0000, 0.0000 ) == Molecular Dipole wrt Center of Mass == mu = ( -0.0000, -0.0001, -0.0001 ) au |mu| = 0.0001 au, 0.0002 Debye output psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs == Timing == cputime in seconds prologue : 0.373890E+01 main loop : 0.131647E+01 epilogue : 0.278618E+00 total : 0.533398E+01 cputime/step: 0.263294E+00 ( 5 evalulations, 1 linesearches) Time spent doing total step percent total time : 0.536761E+01 0.107352E+01 100.0 % i/o time : 0.515451E+00 0.103090E+00 9.6 % FFTs : 0.373968E+00 0.747937E-01 7.0 % dot products : 0.139279E-01 0.278559E-02 0.3 % geodesic : 0.158570E-01 0.317140E-02 0.3 % ffm_dgemm : 0.124407E-02 0.248814E-03 0.0 % fmf_dgemm : 0.143151E-01 0.286302E-02 0.3 % mmm_dgemm : 0.932217E-04 0.186443E-04 0.0 % m_diagonalize : 0.216006E-03 0.432012E-04 0.0 % exchange correlation : 0.642178E+00 0.128436E+00 12.0 % local pseudopotentials : 0.491142E-03 0.982285E-04 0.0 % non-local pseudopotentials : 0.115800E-01 0.231600E-02 0.2 % hartree potentials : 0.974018E-01 0.194804E-01 1.8 % ion-ion interaction : 0.717878E-03 0.143576E-03 0.0 % structure factors : 0.951767E-03 0.190353E-03 0.0 % phase factors : 0.905891E-05 0.181178E-05 0.0 % masking and packing : 0.116353E+00 0.232705E-01 2.2 % queue fft : 0.420149E+00 0.840298E-01 7.8 % queue fft (serial) : 0.152391E+00 0.304782E-01 2.8 % queue fft (message passing): 0.253164E+00 0.506328E-01 4.7 % HFX potential : 0.357959E+00 0.715918E-01 6.7 % non-local psp FFM : 0.824619E-02 0.164924E-02 0.2 % non-local psp FMF : 0.232506E-02 0.465012E-03 0.0 % non-local psp FFM A : 0.130009E-02 0.260019E-03 0.0 % non-local psp FFM B : 0.891922E-03 0.178384E-03 0.0 % == HFX Counters == HFX Evalulations = 5 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:14:23 2021 <<< Saving state for pspw with suffix hess /home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs initial hessian zero matrix atom: 1 xyz: 1(+) wall time: 78.1 date: Sat Nov 20 11:14:23 2021 Total PSPW energy : -0.5424770242E+01 == Timing == cputime in seconds prologue : 0.373923E+01 main loop : 0.113446E+02 epilogue : 0.274254E+00 total : 0.153581E+02 cputime/step: 0.257832E+00 ( 44 evalulations, 20 linesearches) == HFX Counters == HFX Evalulations = 44 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:14:39 2021 <<< atom: 1 xyz: 1(-) wall time: 93.8 date: Sat Nov 20 11:14:39 2021 Total PSPW energy : -0.5424770258E+01 == Timing == cputime in seconds prologue : 0.375943E+01 main loop : 0.123633E+02 epilogue : 0.274694E+00 total : 0.163974E+02 cputime/step: 0.280983E+00 ( 44 evalulations, 20 linesearches) == HFX Counters == HFX Evalulations = 44 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:14:56 2021 <<< atom: 1 xyz: 2(+) wall time: 110.6 date: Sat Nov 20 11:14:56 2021 Total PSPW energy : -0.5424770251E+01 == Timing == cputime in seconds prologue : 0.373778E+01 main loop : 0.113867E+02 epilogue : 0.276789E+00 total : 0.154013E+02 cputime/step: 0.258790E+00 ( 44 evalulations, 20 linesearches) == HFX Counters == HFX Evalulations = 44 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:15:14 2021 <<< atom: 1 xyz: 2(-) wall time: 128.2 date: Sat Nov 20 11:15:14 2021 Total PSPW energy : -0.5424770267E+01 == Timing == cputime in seconds prologue : 0.375183E+01 main loop : 0.113556E+02 epilogue : 0.278070E+00 total : 0.153855E+02 cputime/step: 0.258083E+00 ( 44 evalulations, 20 linesearches) == HFX Counters == HFX Evalulations = 44 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:15:29 2021 <<< atom: 1 xyz: 3(+) wall time: 143.9 date: Sat Nov 20 11:15:29 2021 Total PSPW energy : -0.5424770233E+01 == Timing == cputime in seconds prologue : 0.374044E+01 main loop : 0.114193E+02 epilogue : 0.275751E+00 total : 0.154355E+02 cputime/step: 0.259529E+00 ( 44 evalulations, 20 linesearches) == HFX Counters == HFX Evalulations = 44 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:15:45 2021 <<< atom: 1 xyz: 3(-) wall time: 159.7 date: Sat Nov 20 11:15:45 2021 Total PSPW energy : -0.5424770243E+01 == Timing == cputime in seconds prologue : 0.374128E+01 main loop : 0.113841E+02 epilogue : 0.275818E+00 total : 0.154012E+02 cputime/step: 0.258729E+00 ( 44 evalulations, 20 linesearches) == HFX Counters == HFX Evalulations = 44 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:16:01 2021 <<< finite difference hessian delta = 1.00000000000000002E-002 1 2 3 1 0.0043 -0.0072 -0.0051 2 -0.0072 0.0041 -0.0068 3 -0.0051 -0.0068 0.0045 finite difference derivative dipole; delta = 1.00000000000000002E-002 X vector of derivative dipole (au) [debye/angstrom] d_dipole_x/ = 0.1290 [ 0.6199] d_dipole_x/ = -0.0698 [ -0.3354] d_dipole_x/ = -0.0592 [ -0.2844] Y vector of derivative dipole (au) [debye/angstrom] d_dipole_y/ = -0.0702 [ -0.3372] d_dipole_y/ = 0.1200 [ 0.5766] d_dipole_y/ = -0.0670 [ -0.3218] Z vector of derivative dipole (au) [debye/angstrom] d_dipole_z/ = -0.0591 [ -0.2838] d_dipole_z/ = -0.0665 [ -0.3196] d_dipole_z/ = 0.1328 [ 0.6379] triangle hessian written to /home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.hess derivative dipole written to /home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.fd_ddipole Deleting state for pspw with suffix hess /home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs 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 -------------------------------------------------------------------------- C 1 0.0000000D+00 0.0000000D+00 0.0000000D+00 1.2000000D+01 -------------------------------------------------------------------------- ---------------------------------------------------- MASS-WEIGHTED NUCLEAR HESSIAN (Hartree/Bohr/Bohr/Kamu) ---------------------------------------------------- 1 2 3 ----- ----- ----- ----- ----- 1 3.54761D-01 2 -6.00512D-01 3.40295D-01 3 -4.21028D-01 -5.68755D-01 3.75315D-01 ------------------------------------------------- NORMAL MODE EIGENVECTORS IN CARTESIAN COORDINATES ------------------------------------------------- (Frequencies expressed in cm-1) 1 2 3 Frequency 0.00 0.00 0.00 1 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 ---------------------------------------------------------------------------- Normal Eigenvalue || Derivative Dipole Moments (debye/angs) Mode [cm**-1] || [d/dqX] [d/dqY] [d/dqZ] ------ ---------- || ------------------ ------------------ ----------------- 1 0.000 || 0.000 0.000 0.000 2 0.000 || 0.000 0.000 0.000 3 0.000 || 0.000 0.000 0.000 ---------------------------------------------------------------------------- ---------------------------------------------------------------------------- Normal Eigenvalue || Infra Red Intensities Mode [cm**-1] || [atomic units] [(debye/angs)**2] [(KM/mol)] [arbitrary] ------ ---------- || -------------- ----------------- ---------- ----------- 1 0.000 || 0.000000 0.000 0.000 NaN 2 0.000 || 0.000000 0.000 0.000 NaN 3 0.000 || 0.000000 0.000 0.000 NaN ---------------------------------------------------------------------------- Vibrational analysis via the FX method --- with translations and rotations projected out --- --- via the Eckart algorithm --- Dependent rotation vector no. 1 found in ECKART; assuming linear geometry Dependent rotation vector no. 2 found in ECKART; assuming linear geometry Dependent rotation vector no. 3 found in ECKART; assuming linear geometry Projected Nuclear Hessian trans-rot subspace norm:0.0000D+00 (should be close to zero!) -------------------------------------------------------- MASS-WEIGHTED PROJECTED HESSIAN (Hartree/Bohr/Bohr/Kamu) -------------------------------------------------------- 1 2 3 ----- ----- ----- ----- ----- 1 0.00000D+00 2 0.00000D+00 0.00000D+00 3 0.00000D+00 0.00000D+00 0.00000D+00 center of mass -------------- x = 0.00000000 y = 0.00000000 z = 0.00000000 moments of inertia (a.u.) ------------------ 0.000000000000 0.000000000000 0.000000000000 0.000000000000 0.000000000000 0.000000000000 0.000000000000 0.000000000000 0.000000000000 Rotational Constants -------------------- A= 0.000000 cm-1 ( 0.000000 K) B= 0.000000 cm-1 ( 0.000000 K) C= 0.000000 cm-1 ( 0.000000 K) Temperature = 298.15K frequency scaling parameter = 1.0000 Atom Zero-Point correction to Energy = 0.000 kcal/mol ( 0.000000 au) Thermal correction to Energy = 0.888 kcal/mol ( 0.001416 au) Thermal correction to Enthalpy = 1.481 kcal/mol ( 0.002359 au) Total Entropy = 33.383 cal/mol-K - Translational = 33.383 cal/mol-K (mol. weight = 12.0000) - Rotational = 0.000 cal/mol-K (symmetry # = 1) - Vibrational = 0.000 cal/mol-K Cv (constant volume heat capacity) = 2.979 cal/mol-K - Translational = 2.979 cal/mol-K - Rotational = 0.000 cal/mol-K - Vibrational = 0.000 cal/mol-K ------------------------------------------------- NORMAL MODE EIGENVECTORS IN CARTESIAN COORDINATES ------------------------------------------------- (Projected Frequencies expressed in cm-1) 1 2 3 P.Frequency 0.00 0.00 0.00 1 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 ---------------------------------------------------------------------------- Normal Eigenvalue || Projected Derivative Dipole Moments (debye/angs) Mode [cm**-1] || [d/dqX] [d/dqY] [d/dqZ] ------ ---------- || ------------------ ------------------ ----------------- 1 0.000 || 0.000 0.000 0.000 2 0.000 || 0.000 0.000 0.000 3 0.000 || 0.000 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.000000 0.000 0.000 NaN 2 0.000 || 0.000000 0.000 0.000 NaN 3 0.000 || 0.000000 0.000 0.000 NaN ---------------------------------------------------------------------------- vib:animation F Task times cpu: 98.6s wall: 103.2s NWChem Input Module ------------------- **************************************************** * * * NWPW PSPW Calculation * * * * [ (Grassmann/Stiefel manifold implementation) ] * * * * [ NorthWest Chemistry implementation ] * * * * version #5.10 06/12/02 * * * * This code was developed by Eric J. Bylaska, * * and was based upon algorithms and code * * developed by the group of Prof. John H. Weare * * * **************************************************** >>> JOB STARTED AT Sat Nov 20 11:16:01 2021 <<< ================ input data ======================== input psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs number of processors used: 32 processor grid : 32 x 1 parallel mapping :2d-hilbert parallel mapping : balanced number of threads : 1 parallel io : off options: boundary conditions = periodic (version3) electron spin = unrestricted exchange-correlation = PBE0 (White and Bird) parameterization - HFX relaxed - HFX alpha orbitals: 1 2 3 - HFX beta orbitals : 4 - HFX screened coulomb solver - HFX screening radius(pspw:HFX_screening_radius): 0.800E+01 - HFX screening power (pspw:HFX_screening_power) : 0.800E+01 - HFX screening type (pspw:HFX_screening_type) : 0 - HFX scaling parameter (pspw:HFX_parameter) : 0.250E+00 elements involved in the cluster: 1: C valence charge: 4.0000 lmax= 2 comment : Parameterized (J.Phys.Chem., vol 100, page 6966) Hamman psp pseudpotential type : 0 highest angular component : 2 local potential used : 2 number of non-local projections: 4 cutoff = 0.800 0.850 0.850 total charge: 0.000 atomic composition: C : 1 number of electrons: spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) number of orbitals : spin up= 3 ( 3 per task) down= 1 ( 1 per task) (Fourier space) supercell: cell_name: cell_default lattice: a1=< 19.000 19.000 0.000 > a2=< 19.000 0.000 19.000 > a3=< 0.000 19.000 19.000 > reciprocal: b1=< 0.165 0.165 -0.165 > b2=< 0.165 -0.165 0.165 > b3=< -0.165 0.165 0.165 > lattice: a= 26.870 b= 26.870 c= 26.870 alpha= 60.000 beta= 60.000 gamma= 60.000 omega= 13718.0 density cutoff=100.000 fft=100x100x100( 282860 waves 8839 per task) wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 3616 per task) Ewald summation: cut radius= 8.55 and 1 Madelung Wigner-Seitz= 1.79174723 (alpha= 2.88828212 rs=14.85026019) technical parameters: time step= 5.80 fictitious mass= 400000.0 tolerance=0.100E-06 (energy) 0.100E-06 (density) maximum iterations = 1000 ( 10 inner 100 outer ) == Energy Calculation == ============ Grassmann lmbfgs iteration ============ >>> ITERATION STARTED AT Sat Nov 20 11:16:05 2021 <<< iter. Energy DeltaE DeltaRho ------------------------------------------------------ 10 -0.5424770425E+01 -0.27853E-08 0.22266E-11 *** tolerance ok. iteration terminated >>> ITERATION ENDED AT Sat Nov 20 11:16:06 2021 <<< == Summary Of Results == number of electrons: spin up= 3.00000 down= 1.00000 (real space) total energy : -0.5424770425E+01 ( -0.54248E+01/ion) total orbital energy: -0.1652108162E+01 ( -0.41303E+00/electron) hartree energy : 0.3505352777E+01 ( 0.87634E+00/electron) exc-corr energy : -0.1158535906E+01 ( -0.28963E+00/electron) HF exchange energy : -0.3568831490E+00 ( -0.89221E-01/electron) ion-ion energy : -0.9652341209E+00 ( -0.96523E+00/ion) kinetic (planewave) : 0.3196271168E+01 ( 0.79907E+00/electron) V_local (planewave) : -0.9960255795E+01 ( -0.24901E+01/electron) V_nl (planewave) : 0.3145146007E+00 ( 0.78629E-01/electron) V_Coul (planewave) : 0.7010705554E+01 ( 0.17527E+01/electron) V_xc. (planewave) : -0.1499577392E+01 ( -0.37489E+00/electron) K.S. HFX energy : -0.7137662980E+00 ( -0.17844E+00/electron) Virial Coefficient : -0.1293573922E+01 orbital energies: -0.2808764E+00 ( -7.643eV) -0.2808824E+00 ( -7.643eV) -0.6102206E+00 ( -16.605eV) -0.4801288E+00 ( -13.065eV) == Virtual Orbital Calculation == random planewave guess, initial psi:pspw-pbe0-134211.emovecs - spin, nalpha, nbeta: 2 8 8 input epsi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.emovecs orbital 1 current e=-0.146E+00 (error=0.728E-06) iterations 121( 104 preconditioned, Ep,Sp= 20.0 200.0) orbital 1 current e=-0.146E+00 (error=0.439E-07) iterations 10( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 2 current e=-0.523E-02 (error=0.144E-06) iterations 121( 78 preconditioned, Ep,Sp= 20.0 200.0) orbital 2 current e=-0.523E-02 (error=0.415E-07) iterations 2( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 3 current e= 0.336E-01 (error=0.145E-04) iterations 121( 77 preconditioned, Ep,Sp= 20.0 200.0) orbital 3 current e= 0.330E-01 (error=0.459E-05) iterations 121( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 3 current e= 0.324E-01 (error=0.973E-07) iterations 87( 28 preconditioned, Ep,Sp= 20.0 200.0) orbital 4 current e= 0.345E-01 (error=0.321E-04) iterations 121( 120 preconditioned, Ep,Sp= 20.0 200.0) orbital 4 current e= 0.337E-01 (error=0.131E-05) iterations 121( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 4 current e= 0.396E-01 (error=0.230E-04) iterations 121( 102 preconditioned, Ep,Sp= 20.0 200.0) orbital 4 current e= 0.363E-01 (error=0.503E-04) iterations 121( 39 preconditioned, Ep,Sp= 20.0 200.0) orbital 4 current e= 0.336E-01 (error=0.314E-06) iterations 121( 30 preconditioned, Ep,Sp= 20.0 200.0) orbital 4 current e= 0.336E-01 (error=0.947E-07) iterations 65( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 5 current e= 0.372E-01 (error=0.590E-04) iterations 121( 107 preconditioned, Ep,Sp= 20.0 200.0) orbital 5 current e= 0.336E-01 (error=0.146E-05) iterations 121( 68 preconditioned, Ep,Sp= 20.0 200.0) orbital 5 current e= 0.396E-01 (error=0.229E-04) iterations 121( 102 preconditioned, Ep,Sp= 20.0 200.0) orbital 5 current e= 0.394E-01 (error=0.626E-05) iterations 121( 4 preconditioned, Ep,Sp= 20.0 200.0) orbital 5 current e= 0.373E-01 (error=0.177E-04) iterations 121( 56 preconditioned, Ep,Sp= 20.0 200.0) orbital 5 current e= 0.363E-01 (error=0.185E-06) iterations 121( 9 preconditioned, Ep,Sp= 20.0 200.0) orbital 5 current e= 0.362E-01 (error=0.935E-07) iterations 107( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 6 current e= 0.387E-01 (error=0.114E-04) iterations 121( 88 preconditioned, Ep,Sp= 20.0 200.0) orbital 6 current e= 0.338E-01 (error=0.590E-05) iterations 121( 53 preconditioned, Ep,Sp= 20.0 200.0) orbital 6 current e= 0.396E-01 (error=0.229E-04) iterations 121( 102 preconditioned, Ep,Sp= 20.0 200.0) orbital 6 current e= 0.394E-01 (error=0.207E-05) iterations 121( 4 preconditioned, Ep,Sp= 20.0 200.0) orbital 6 current e= 0.389E-01 (error=0.970E-05) iterations 121( 6 preconditioned, Ep,Sp= 20.0 200.0) orbital 6 current e= 0.339E-01 (error=0.712E-05) iterations 121( 96 preconditioned, Ep,Sp= 20.0 200.0) orbital 6 current e= 0.336E-01 (error=0.919E-07) iterations 104( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 7 current e= 0.391E-01 (error=0.473E-06) iterations 121( 84 preconditioned, Ep,Sp= 20.0 200.0) orbital 7 current e= 0.391E-01 (error=0.943E-07) iterations 45( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 8 current e= 0.401E-01 (error=0.756E-05) iterations 121( 91 preconditioned, Ep,Sp= 20.0 200.0) orbital 8 current e= 0.392E-01 (error=0.479E-06) iterations 121( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 8 current e= 0.397E-01 (error=0.251E-04) iterations 121( 100 preconditioned, Ep,Sp= 20.0 200.0) orbital 8 current e= 0.395E-01 (error=0.985E-07) iterations 74( 4 preconditioned, Ep,Sp= 20.0 200.0) orbital 9 current e=-0.108E+00 (error=0.540E-03) iterations 121( 120 preconditioned, Ep,Sp= 20.0 200.0) orbital 9 current e=-0.114E+00 (error=0.781E-07) iterations 38( 13 preconditioned, Ep,Sp= 20.0 200.0) orbital 10 current e=-0.538E-01 (error=0.950E-07) iterations 113( 76 preconditioned, Ep,Sp= 20.0 200.0) orbital 11 current e=-0.533E-01 (error=0.624E-04) iterations 121( 120 preconditioned, Ep,Sp= 20.0 200.0) orbital 11 current e=-0.538E-01 (error=0.950E-07) iterations 37( 6 preconditioned, Ep,Sp= 20.0 200.0) orbital 12 current e= 0.439E-03 (error=0.882E-07) iterations 106( 70 preconditioned, Ep,Sp= 20.0 200.0) orbital 13 current e= 0.422E-01 (error=0.638E-05) iterations 121( 75 preconditioned, Ep,Sp= 20.0 200.0) orbital 13 current e= 0.415E-01 (error=0.844E-06) iterations 121( 4 preconditioned, Ep,Sp= 20.0 200.0) orbital 13 current e= 0.418E-01 (error=0.101E-04) iterations 121( 71 preconditioned, Ep,Sp= 20.0 200.0) orbital 13 current e= 0.412E-01 (error=0.191E-05) iterations 121( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 13 current e= 0.412E-01 (error=0.958E-07) iterations 25( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 14 current e= 0.425E-01 (error=0.572E-05) iterations 121( 83 preconditioned, Ep,Sp= 20.0 200.0) orbital 14 current e= 0.417E-01 (error=0.741E-06) iterations 121( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 14 current e= 0.430E-01 (error=0.945E-07) iterations 115( 78 preconditioned, Ep,Sp= 20.0 200.0) orbital 15 current e= 0.424E-01 (error=0.119E-04) iterations 121( 55 preconditioned, Ep,Sp= 20.0 200.0) orbital 15 current e= 0.418E-01 (error=0.145E-05) iterations 121( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 15 current e= 0.566E-01 (error=0.325E-03) iterations 121( 54 preconditioned, Ep,Sp= 20.0 200.0) orbital 15 current e= 0.436E-01 (error=0.474E-06) iterations 121( 77 preconditioned, Ep,Sp= 20.0 200.0) orbital 15 current e= 0.419E-01 (error=0.345E-05) iterations 121( 27 preconditioned, Ep,Sp= 20.0 200.0) orbital 15 current e= 0.415E-01 (error=0.488E-06) iterations 121( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 15 current e= 0.413E-01 (error=0.281E-06) iterations 121( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 16 current e= 0.424E-01 (error=0.432E-05) iterations 121( 76 preconditioned, Ep,Sp= 20.0 200.0) orbital 16 current e= 0.418E-01 (error=0.854E-06) iterations 121( 1 preconditioned, Ep,Sp= 20.0 200.0) orbital 16 current e= 0.565E-01 (error=0.181E-03) iterations 121( 55 preconditioned, Ep,Sp= 20.0 200.0) orbital 16 current e= 0.436E-01 (error=0.159E-06) iterations 121( 76 preconditioned, Ep,Sp= 20.0 200.0) orbital 16 current e= 0.436E-01 (error=0.599E-07) iterations 2( 1 preconditioned, Ep,Sp= 20.0 200.0) virtual orbital energies: 0.3948473E-01 ( 1.074eV) 0.4363889E-01 ( 1.187eV) 0.3907676E-01 ( 1.063eV) 0.4295568E-01 ( 1.169eV) 0.3621855E-01 ( 0.986eV) 0.4128119E-01 ( 1.123eV) 0.3362617E-01 ( 0.915eV) 0.4122574E-01 ( 1.122eV) 0.3362365E-01 ( 0.915eV) 0.4388081E-03 ( 0.012eV) 0.3237430E-01 ( 0.881eV) -0.5378498E-01 ( -1.464eV) -0.5230856E-02 ( -0.142eV) -0.5378801E-01 ( -1.464eV) -0.1459188E+00 ( -3.971eV) -0.1139015E+00 ( -3.099eV) output epsi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.emovecs Total PSPW energy : -0.5424770425E+01 === Spin Contamination === = 2.0000000000000000 = 2.0045828181603347 == Center of Charge == spin up ( -0.0001, -0.0001, -0.0001 ) spin down ( 0.0003, 0.0003, 0.0003 ) total ( 0.0000, 0.0000, 0.0000 ) ionic ( 0.0000, 0.0000, 0.0000 ) == Molecular Dipole wrt Center of Mass == mu = ( -0.0000, -0.0001, -0.0001 ) au |mu| = 0.0001 au, 0.0002 Debye output psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs == Timing == cputime in seconds prologue : 0.376686E+01 main loop : 0.432094E+03 epilogue : 0.449321E+00 total : 0.436310E+03 cputime/step: 0.864188E+02 ( 5 evalulations, 1 linesearches) Time spent doing total step percent total time : 0.436345E+03 0.872690E+02 100.0 % i/o time : 0.242030E+01 0.484060E+00 0.6 % FFTs : 0.274242E+03 0.548484E+02 62.8 % dot products : 0.249236E+02 0.498473E+01 5.7 % geodesic : 0.158331E-01 0.316663E-02 0.0 % ffm_dgemm : 0.139118E-02 0.278235E-03 0.0 % fmf_dgemm : 0.166142E-01 0.332284E-02 0.0 % mmm_dgemm : 0.917912E-04 0.183582E-04 0.0 % m_diagonalize : 0.404595E-03 0.809191E-04 0.0 % exchange correlation : 0.769211E+00 0.153842E+00 0.2 % local pseudopotentials : 0.471830E-03 0.943661E-04 0.0 % non-local pseudopotentials : 0.124701E+02 0.249402E+01 2.9 % hartree potentials : 0.171281E+02 0.342563E+01 3.9 % ion-ion interaction : 0.735044E-03 0.147009E-03 0.0 % structure factors : 0.114521E+01 0.229041E+00 0.3 % phase factors : 0.715256E-05 0.143051E-05 0.0 % masking and packing : 0.521538E+02 0.104308E+02 12.0 % queue fft : 0.422278E+00 0.844556E-01 0.1 % queue fft (serial) : 0.149233E+00 0.298467E-01 0.0 % queue fft (message passing): 0.258212E+00 0.516425E-01 0.1 % HFX potential : 0.297981E+03 0.595963E+02 68.3 % non-local psp FFM : 0.672862E+01 0.134572E+01 1.5 % non-local psp FMF : 0.392055E+01 0.784111E+00 0.9 % non-local psp FFM A : 0.175180E+01 0.350360E+00 0.4 % non-local psp FFM B : 0.753658E+00 0.150732E+00 0.2 % == HFX Counters == HFX Evalulations = 5 Lin ACE HFX Evalulations = 0 >>> JOB COMPLETED AT Sat Nov 20 11:23:17 2021 <<< Task times cpu: 433.9s wall: 436.4s NWChem Input Module ------------------- >>>> PSPW Parallel Module - pspw_dplot <<<< **************************************************** * * * pspw DPLOT * * * * [ Generates density and orbital grids ] * * * * [ NorthWest Chemistry implementation ] * * * * version #1.00 08/22/01 * * * **************************************************** >>> JOB STARTED AT Sat Nov 20 11:23:18 2021 <<< ============ PSPW DPLOT input data ================= input psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.movecs number of processors used: 32 parallel mapping :2d hilbert parallel mapping : balanced options: boundary conditions = periodic (version3) number of electrons: spin up= 3 spin down= 1 ( Fourier space) ncell = 0 0 0 position tolerance = 0.100000E-02 origin=< 0.000 0.000 0.000 > supercell: lattice: a1=< 19.000 19.000 0.000 > a2=< 19.000 0.000 19.000 > a3=< 0.000 19.000 19.000 > b1=< 0.165 0.165 -0.165 > b2=< 0.165 -0.165 0.165 > b3=< -0.165 0.165 0.165 > volume : 13718.0 density cutoff=100.000 fft=100x100x100( 282860 waves 8839 per task) wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 3616 per task) writing orbital 1 to filename: homo-alpha.cube writing orbital 4 to filename: homo-beta.cube writing total density to filename: density.cube ----------------- cputime in seconds total : 14.772498846054077 >>> JOB COMPLETED AT Sat Nov 20 11:23:32 2021 <<< Task times cpu: 14.6s wall: 14.8s NWChem Input Module ------------------- >>>> PSPW Parallel Module - pspw_dplot <<<< **************************************************** * * * pspw DPLOT * * * * [ Generates density and orbital grids ] * * * * [ NorthWest Chemistry implementation ] * * * * version #1.00 08/22/01 * * * **************************************************** >>> JOB STARTED AT Sat Nov 20 11:23:32 2021 <<< ============ PSPW DPLOT input data ================= input psi filename:/home/bylaska/Projects/Work/RUNARROWS0/pspw-pbe0-134211.emovecs number of processors used: 32 parallel mapping :2d hilbert parallel mapping : balanced options: boundary conditions = periodic (version3) number of electrons: spin up= 8 spin down= 8 ( Fourier space) ncell = 0 0 0 position tolerance = 0.100000E-02 origin=< 0.000 0.000 0.000 > supercell: lattice: a1=< 19.000 19.000 0.000 > a2=< 19.000 0.000 19.000 > a3=< 0.000 19.000 19.000 > b1=< 0.165 0.165 -0.165 > b2=< 0.165 -0.165 0.165 > b3=< -0.165 0.165 0.165 > volume : 13718.0 density cutoff=100.000 fft=100x100x100( 282860 waves 8839 per task) wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 3616 per task) writing orbital 1 to filename: lumo-alpha.cube writing orbital 9 to filename: lumo-beta.cube ----------------- cputime in seconds total : 10.762134075164795 >>> JOB COMPLETED AT Sat Nov 20 11:23:43 2021 <<< Task times cpu: 10.6s wall: 10.8s NWChem Input Module ------------------- Summary of allocated global arrays ----------------------------------- No active global arrays GA Statistics for process 0 ------------------------------ create destroy get put acc scatter gather read&inc calls: 0 0 0 0 0 0 0 0 number of processes/call 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 bytes total: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 bytes remote: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Max memory consumed for GA by this process: 0 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 284 25 current total bytes 0 0 maximum total bytes 31410776 8161272 maximum total K-bytes 31411 8162 maximum total M-bytes 32 9 CITATION -------- Please cite the following reference when publishing results obtained with NWChem: E. Apra, E. J. Bylaska, W. A. de Jong, N. Govind, K. Kowalski, T. P. Straatsma, M. Valiev, H. J. J. van Dam, Y. Alexeev, J. Anchell, V. Anisimov, F. W. Aquino, R. Atta-Fynn, J. Autschbach, N. P. Bauman, J. C. Becca, D. E. Bernholdt, K. Bhaskaran-Nair, S. Bogatko, P. Borowski, J. Boschen, J. Brabec, A. Bruner, E. Cauet, Y. Chen, G. N. Chuev, C. J. Cramer, J. Daily, M. J. O. Deegan, T. H. Dunning Jr., M. Dupuis, K. G. Dyall, G. I. Fann, S. A. Fischer, A. Fonari, H. Fruchtl, L. Gagliardi, J. Garza, N. Gawande, S. Ghosh, K. Glaesemann, A. W. Gotz, J. Hammond, V. Helms, E. D. Hermes, K. Hirao, S. Hirata, M. Jacquelin, L. Jensen, B. G. Johnson, H. Jonsson, R. A. Kendall, M. Klemm, R. Kobayashi, V. Konkov, S. Krishnamoorthy, M. Krishnan, Z. Lin, R. D. Lins, R. J. Littlefield, A. J. Logsdail, K. Lopata, W. Ma, A. V. Marenich, J. Martin del Campo, D. Mejia-Rodriguez, J. E. Moore, J. M. Mullin, T. Nakajima, D. R. Nascimento, J. A. Nichols, P. J. Nichols, J. Nieplocha, A. Otero-de-la-Roza, B. Palmer, A. Panyala, T. Pirojsirikul, B. Peng, R. Peverati, J. Pittner, L. Pollack, R. M. Richard, P. Sadayappan, G. C. Schatz, W. A. Shelton, D. W. Silverstein, D. M. A. Smith, T. A. Soares, D. Song, M. Swart, H. L. Taylor, G. S. Thomas, V. Tipparaju, D. G. Truhlar, K. Tsemekhman, T. Van Voorhis, A. Vazquez-Mayagoitia, P. Verma, O. Villa, A. Vishnu, K. D. Vogiatzis, D. Wang, J. H. Weare, M. J. Williamson, T. L. Windus, K. Wolinski, A. T. Wong, Q. Wu, C. Yang, Q. Yu, M. Zacharias, Z. Zhang, Y. Zhao, and R. J. Harrison "NWChem: Past, present, and future J. Chem. Phys. 152, 184102 (2020) doi:10.1063/5.0004997 AUTHORS ------- E. Apra, E. J. Bylaska, N. Govind, K. Kowalski, M. Valiev, W. A. de Jong, T. P. Straatsma, H. J. J. van Dam, D. Wang, T. L. Windus, N. P. Bauman, A. Panyala, J. Hammond, J. Autschbach, K. Bhaskaran-Nair, J. Brabec, K. Lopata, S. A. Fischer, S. Krishnamoorthy, M. Jacquelin, W. Ma, M. Klemm, O. Villa, Y. Chen, V. Anisimov, F. Aquino, S. Hirata, M. T. Hackler, Eric Hermes, L. Jensen, J. E. Moore, J. C. Becca, V. Konjkov, D. Mejia-Rodriguez, T. Risthaus, M. Malagoli, A. Marenich, A. Otero-de-la-Roza, J. Mullin, P. Nichols, R. Peverati, J. Pittner, Y. Zhao, P.-D. Fan, A. Fonari, M. J. Williamson, R. J. Harrison, J. R. Rehr, M. Dupuis, D. Silverstein, D. M. A. Smith, J. Nieplocha, V. Tipparaju, M. Krishnan, B. E. Van Kuiken, A. Vazquez-Mayagoitia, 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. A. Kendall, J. A. Nichols, K. Tsemekhman, K. Wolinski, J. Anchell, D. E. Bernholdt, P. Borowski, T. Clark, D. Clerc, H. Dachsel, M. J. O. Deegan, K. Dyall, D. Elwood, E. Glendening, M. Gutowski, A. C. Hess, J. Jaffe, B. G. 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. T. Wong, Z. Zhang. Total times cpu: 628.9s wall: 637.7s # MYMACHINENAME: Eric Bylaska - arrow13.emsl.pnl.gov :MYMACHINENAME