Results from an EMSL Arrows Calculation
 |
EMSL Arrows is a revolutionary approach to materials and chemical simulations that uses NWChem and chemical computational databases to make materials and chemical modeling accessible via a broad spectrum of digital communications including posts to web APIs, social networks, and traditional email. |
Molecular modeling software has previously been extremely complex, making it prohibative to all but experts in the field, yet even experts can struggle to perform calculations. This service is designed to be used by experts and non-experts alike. Experts can carry out and keep track of large numbers of complex calculations with diverse levels of theories present in their workflows. Additionally, due to a streamlined and easy-to-use input, non-experts can carry out a wide variety of molecular modeling calculations previously not accessible to them.
Link back to EMSL Arrows API
##################### start nwoutput #######################
nwout file for Id=44766
bylaska@archive.emsl.pnl.gov:chemdb2/80/32/k+.out-2017-11-21-13:7:38
argument 1 = k+.nw
NWChem w/ OpenMP: maximum threads = 8
============================== echo of input deck ==============================
# NWChem Input Generation (tnt_submit5) - The current time is Tue Nov 21 12:01:30 2017
# - adding tag osmiles:[K+]:osmiles to input deck.
#
# - pubchem_synonyms = ['potassium ion', 'Potassium cation', 'potassium(1+)', 'Potassium(+)', 'Potassium (ion)', 'Potassium ion(+)', 'Potassium monocation', 'Potassium (K+)', 'potassium(I) cation', 'potassium(1+) ion', 'Potassium ion(1+)', 'Potassium ion
#
# - queue_number = 79089
# - mformula = K1
# - name = [K+]
# - smiles = [K+]
# - csmiles = [K+]
# - InChI = InChI=1S/K/q+1
# - InChIKey = NPYPAHLBTDXSSS-UHFFFAOYSA-N
# - pubchem_cid = 813
# - pubchem_smiles = [K+]
# - pubchem_iupac = potassium(1+)
# - pubchem_synonym0 = potassium ion
# - theory = pspw
# - pspw4 = True
# - paw = False
# - xc = pbe
# - basis = default
# - basisHZ = default
# - theory_property = pspw
# - property_pspw4 = True
# - property_paw = False
# - xc_property = pbe
# - basis_property = default
# - basisHZ_property = default
# - type = ovb
# - solvation_type =
# - charge = 1
# - mult = 1
# - 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
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
#
# +
# K
#
#
#
#
#
title "swnc: ovb theory=pspw xc=pbe formula=K1 charge=1 mult=1 machinejob:Shirky "
#vtag= osmiles:[K+]:osmiles
echo
start pspw-pbe-K1-79089
memory 1900 mb
charge 1
geometry units angstroms print xyz noautosym
K 0.93527 0.06775 0.03493
end
nwpw
cutoff 50.0
mult 1
xc pbe96
lmbfgs
end
nwpw
simulation_cell
fcc 38.0
boundary_conditions aperiodic
end
end
driver; 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-restricted.cube
density total density.cube
end
end
#task pspw energy ignore
#task pspw pspw_dplot
nwpw
dplot
vectors pspw-pbe-K1-79089.emovecs
orbital 1 lumo-restricted.cube
end
end
#task pspw pspw_dplot
================================================================================
Northwest Computational Chemistry Package (NWChem) 6.6
------------------------------------------------------
Environmental Molecular Sciences Laboratory
Pacific Northwest National Laboratory
Richland, WA 99352
Copyright (c) 1994-2015
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 = WE34304
program = nwchem
date = Tue Nov 21 13:02:17 2017
compiled = Tue_Nov_21_13:02:01_2017
source = /Users/bylaska/nwchem-releases/nwchem
nwchem branch = Development
nwchem revision = 29547
ga revision = N/A
input = k+.nw
prefix = pspw-pbe-K1-79089.
data base = ./pspw-pbe-K1-79089.db
status = startup
nproc = 4
time left = -1s
Memory information
------------------
heap = 62259196 doubles = 475.0 Mbytes
stack = 62259201 doubles = 475.0 Mbytes
global = 124518400 doubles = 950.0 Mbytes (distinct from heap & stack)
total = 249036797 doubles = 1900.0 Mbytes
verify = yes
hardfail = no
Directory information
---------------------
0 permanent = .
0 scratch = .
NWChem Input Module
-------------------
swnc: ovb theory=pspw xc=pbe formula=K1 charge=1 mult=1 machinejob:Shirky
-------------------------------------------------------------------------
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 K 19.0000 0.00000000 0.00000000 0.00000000
Atomic Mass
-----------
K 38.963710
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
K 0.00000000 0.00000000 0.00000000
Deleted DRIVER restart files
NWChem Geometry Optimization
----------------------------
swnc: ovb theory=pspw xc=pbe formula=K1 charge=1 mult=1 machinejob:Shirky
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 K 19.0000 0.00000000 0.00000000 0.00000000
Atomic Mass
-----------
K 38.963710
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 Tue Nov 21 13:02:17 2017 <<<
================ input data ========================
pseudopotential is not correctly formatted:K.vpp
Generated formatted_filename: ./K.vpp
number of processors used: 4
processor grid : 4 x 1
parallel mapping :2d hilbert
parallel mapping : balanced
number of threads : 8
parallel io : off
options:
boundary conditions = aperiodic (version4)
exchange-correlation = PBE96 (White and Bird) parameterization
elements involved in the cluster:
1: K valence charge: 1.0000 lmax= 1
comment : Parameterized on March 19,2009 TM psp
pseudpotential type : 0
highest angular component : 1
local potential used : 0
number of non-local projections: 3
aperiodic cutoff radius : 1.000
semicore corrections included : 2.000 (radius) 4.376 (charge)
cutoff = 4.200 4.050
total charge: 1.000
number of electrons: spin up= 0 ( 0 per task) down= 0 ( 0 per task) (Fourier space)
number of orbitals : spin up= 0 ( 0 per task) down= 0 ( 0 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 8723 per task)
wavefnc cutoff= 12.500 fft= 50x 50x 50( 14406 waves 3601 per task)
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.665641E-01
main loop : -0.103879-316
epilogue : 0.568627E+01
total : 0.561970E+01
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
Time spent doing total step percent
total time : 0.562260E+01 Infinity 100.0 %
i/o time : 0.419612E+01 Infinity 74.6 %
FFTs : 0.149860E-01 Infinity 0.3 %
dot products : 0.500679E-04 Infinity 0.0 %
phase factors : 0.405212E-05 Infinity 0.0 %
masking and packing : 0.257635E-02 Infinity 0.0 %
>>> JOB COMPLETED AT Tue Nov 21 13:02:22 2017 <<<
****************************************************
* *
* 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 Tue Nov 21 13:02:22 2017 <<<
================ input data ========================
pseudopotential is not correctly formatted:K.vpp
Generated formatted_filename: ./K.vpp
number of processors used: 4
processor grid : 4 x 1
parallel mapping :2d hilbert
parallel mapping : balanced
number of threads : 8
parallel io : off
options:
boundary conditions = aperiodic (version4)
exchange-correlation = PBE96 (White and Bird) parameterization
elements involved in the cluster:
1: K valence charge: 1.0000 lmax= 1
comment : Parameterized on March 19,2009 TM psp
pseudpotential type : 0
highest angular component : 1
local potential used : 0
number of non-local projections: 3
aperiodic cutoff radius : 1.000
semicore corrections included : 2.000 (radius) 4.376 (charge)
cutoff = 4.200 4.050
total charge: 1.000
number of electrons: spin up= 0 ( 0 per task) down= 0 ( 0 per task) (Fourier space)
number of orbitals : spin up= 0 ( 0 per task) down= 0 ( 0 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 25442 per task)
wavefnc cutoff= 25.000 fft= 72x 72x 72( 40966 waves 10241 per task)
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.568804E+01
main loop : -0.103879-316
epilogue : 0.345145E+02
total : 0.288264E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
Time spent doing total step percent
total time : 0.288281E+02 Infinity 100.0 %
i/o time : 0.119659E+02 Infinity 41.5 %
FFTs : 0.398478E-01 Infinity 0.1 %
dot products : 0.197907E-04 Infinity 0.0 %
phase factors : 0.190635E-05 Infinity 0.0 %
masking and packing : 0.404882E-02 Infinity 0.0 %
>>> JOB COMPLETED AT Tue Nov 21 13:02:51 2017 <<<
****************************************************
* *
* 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 Tue Nov 21 13:02:51 2017 <<<
================ input data ========================
pseudopotential is not correctly formatted:K.vpp
Generated formatted_filename: ./K.vpp
number of processors used: 4
processor grid : 4 x 1
parallel mapping :2d hilbert
parallel mapping : balanced
number of threads : 8
parallel io : off
options:
boundary conditions = aperiodic (version4)
exchange-correlation = PBE96 (White and Bird) parameterization
elements involved in the cluster:
1: K valence charge: 1.0000 lmax= 1
comment : Parameterized on March 19,2009 TM psp
pseudpotential type : 0
highest angular component : 1
local potential used : 0
number of non-local projections: 3
aperiodic cutoff radius : 1.000
semicore corrections included : 2.000 (radius) 4.376 (charge)
cutoff = 4.200 4.050
total charge: 1.000
number of electrons: spin up= 0 ( 0 per task) down= 0 ( 0 per task) (Fourier space)
number of orbitals : spin up= 0 ( 0 per task) down= 0 ( 0 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 48084 per task)
wavefnc cutoff= 37.500 fft= 90x 90x 90( 75079 waves 18769 per task)
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.345162E+02
main loop : -0.103879-316
epilogue : 0.751874E+02
total : 0.406712E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
Time spent doing total step percent
total time : 0.406729E+02 Infinity 100.0 %
i/o time : 0.238023E+02 Infinity 58.5 %
FFTs : 0.906601E-01 Infinity 0.2 %
dot products : 0.267029E-04 Infinity 0.0 %
phase factors : 0.286002E-05 Infinity 0.0 %
masking and packing : 0.165279E-01 Infinity 0.0 %
>>> JOB COMPLETED AT Tue Nov 21 13:03:32 2017 <<<
****************************************************
* *
* 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 Tue Nov 21 13:03:32 2017 <<<
================ input data ========================
pseudopotential is not correctly formatted:K.vpp
Generated formatted_filename: ./K.vpp
number of processors used: 4
processor grid : 4 x 1
parallel mapping :2d hilbert
parallel mapping : balanced
number of threads : 8
parallel io : off
options:
boundary conditions = aperiodic (version4)
exchange-correlation = PBE96 (White and Bird) parameterization
elements involved in the cluster:
1: K valence charge: 1.0000 lmax= 1
comment : Parameterized on March 19,2009 TM psp
pseudpotential type : 0
highest angular component : 1
local potential used : 0
number of non-local projections: 3
aperiodic cutoff radius : 1.000
semicore corrections included : 2.000 (radius) 4.376 (charge)
cutoff = 4.200 4.050
total charge: 1.000
number of electrons: spin up= 0 ( 0 per task) down= 0 ( 0 per task) (Fourier space)
number of orbitals : spin up= 0 ( 0 per task) down= 0 ( 0 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 70715 per task)
wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 28933 per task)
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.751898E+02
main loop : -0.841473-233
epilogue : 0.124823E+03
total : 0.496337E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
Time spent doing total step percent
total time : 0.496358E+02 Infinity 100.0 %
i/o time : 0.322795E+02 Infinity 65.0 %
FFTs : 0.132630E+00 Infinity 0.3 %
dot products : 0.379086E-04 Infinity 0.0 %
phase factors : 0.190635E-05 Infinity 0.0 %
masking and packing : 0.290196E-01 Infinity 0.1 %
>>> JOB COMPLETED AT Tue Nov 21 13:04:21 2017 <<<
@ Step Energy Delta E Gmax Grms Xrms Xmax Walltime
@ ---- ---------------- -------- -------- -------- -------- -------- --------
@ 0 0.00000000 0.0D+00 0.00000 0.00000 0.00000 0.00000 124.8
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 Tue Nov 21 13:04:21 2017 <<<
================ input data ========================
number of processors used: 4
processor grid : 4 x 1
parallel mapping :2d hilbert
parallel mapping : balanced
number of threads : 8
parallel io : off
options:
boundary conditions = aperiodic (version4)
exchange-correlation = PBE96 (White and Bird) parameterization
elements involved in the cluster:
1: K valence charge: 1.0000 lmax= 1
comment : Parameterized on March 19,2009 TM psp
pseudpotential type : 0
highest angular component : 1
local potential used : 0
number of non-local projections: 3
aperiodic cutoff radius : 1.000
semicore corrections included : 2.000 (radius) 4.376 (charge)
cutoff = 4.200 4.050
total charge: 1.000
number of electrons: spin up= 0 ( 0 per task) down= 0 ( 0 per task) (Fourier space)
number of orbitals : spin up= 0 ( 0 per task) down= 0 ( 0 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 70715 per task)
wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 28933 per task)
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.124834E+03
main loop : -0.693235-309
epilogue : 0.142099E+03
total : 0.172649E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
Time spent doing total step percent
total time : 0.172670E+02 Infinity 100.0 %
i/o time : 0.589540E-01 Infinity 0.3 %
FFTs : 0.127377E+00 Infinity 0.7 %
dot products : 0.293016E-03 Infinity 0.0 %
phase factors : 0.309944E-05 Infinity 0.0 %
masking and packing : 0.196469E-01 Infinity 0.1 %
>>> JOB COMPLETED AT Tue Nov 21 13:04:39 2017 <<<
Line search:
step=-1.00 grad= 0.0D+00 hess= 0.0D+00 energy= 0.000000 mode=accept
new step=-1.00 predicted energy= 0.000000
--------
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 K 19.0000 0.00000000 0.00000000 0.00000000
Atomic Mass
-----------
K 38.963710
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 Tue Nov 21 13:04:39 2017 <<<
================ input data ========================
number of processors used: 4
processor grid : 4 x 1
parallel mapping :2d hilbert
parallel mapping : balanced
number of threads : 8
parallel io : off
options:
boundary conditions = aperiodic (version4)
exchange-correlation = PBE96 (White and Bird) parameterization
elements involved in the cluster:
1: K valence charge: 1.0000 lmax= 1
comment : Parameterized on March 19,2009 TM psp
pseudpotential type : 0
highest angular component : 1
local potential used : 0
number of non-local projections: 3
aperiodic cutoff radius : 1.000
semicore corrections included : 2.000 (radius) 4.376 (charge)
cutoff = 4.200 4.050
total charge: 1.000
number of electrons: spin up= 0 ( 0 per task) down= 0 ( 0 per task) (Fourier space)
number of orbitals : spin up= 0 ( 0 per task) down= 0 ( 0 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 70715 per task)
wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 28933 per task)
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.142106E+03
main loop : -0.841473-233
epilogue : 0.160162E+03
total : 0.180563E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
Time spent doing total step percent
total time : 0.180585E+02 Infinity 100.0 %
i/o time : 0.675161E-01 Infinity 0.4 %
FFTs : 0.127299E+00 Infinity 0.7 %
dot products : 0.288486E-04 Infinity 0.0 %
phase factors : 0.214577E-05 Infinity 0.0 %
masking and packing : 0.212531E-01 Infinity 0.1 %
>>> JOB COMPLETED AT Tue Nov 21 13:04:57 2017 <<<
Step Energy Delta E Gmax Grms Xrms Xmax Walltime
---- ---------------- -------- -------- -------- -------- -------- --------
@ 1 0.00000000 0.0D+00 0.00000 0.00000 0.00000 0.00000 160.2
ok ok ok ok
----------------------
Optimization converged
----------------------
Step Energy Delta E Gmax Grms Xrms Xmax Walltime
---- ---------------- -------- -------- -------- -------- -------- --------
@ 1 0.00000000 0.0D+00 0.00000 0.00000 0.00000 0.00000 160.2
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 K 19.0000 0.00000000 0.00000000 0.00000000
Atomic Mass
-----------
K 38.963710
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: 160.1s wall: 160.1s
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 Tue Nov 21 13:04:57 2017 <<<
================ input data ========================
number of processors used: 4
processor grid : 4 x 1
parallel mapping :2d hilbert
parallel mapping : balanced
number of threads : 8
parallel io : off
options:
boundary conditions = aperiodic (version4)
exchange-correlation = PBE96 (White and Bird) parameterization
elements involved in the cluster:
1: K valence charge: 1.0000 lmax= 1
comment : Parameterized on March 19,2009 TM psp
pseudpotential type : 0
highest angular component : 1
local potential used : 0
number of non-local projections: 3
aperiodic cutoff radius : 1.000
semicore corrections included : 2.000 (radius) 4.376 (charge)
cutoff = 4.200 4.050
total charge: 1.000
number of electrons: spin up= 0 ( 0 per task) down= 0 ( 0 per task) (Fourier space)
number of orbitals : spin up= 0 ( 0 per task) down= 0 ( 0 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 70715 per task)
wavefnc cutoff= 50.000 fft=100x100x100( 115733 waves 28933 per task)
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.160169E+03
main loop : -0.693235-309
epilogue : 0.178067E+03
total : 0.178984E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
Time spent doing total step percent
total time : 0.179006E+02 Infinity 100.0 %
i/o time : 0.771141E-01 Infinity 0.4 %
FFTs : 0.125451E+00 Infinity 0.7 %
dot products : 0.330186E-02 Infinity 0.0 %
phase factors : 0.309844E-05 Infinity 0.0 %
masking and packing : 0.239525E-01 Infinity 0.1 %
>>> JOB COMPLETED AT Tue Nov 21 13:05:15 2017 <<<
Saving state for pspw with suffix hess
initial hessian
zero matrix
atom: 1 xyz: 1(+) wall time: 178.1 date: Tue Nov 21 13:05:15 2017
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.178080E+03
main loop : -0.841473-233
epilogue : 0.195455E+03
total : 0.173743E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
>>> JOB COMPLETED AT Tue Nov 21 13:05:32 2017 <<<
atom: 1 xyz: 1(-) wall time: 195.5 date: Tue Nov 21 13:05:32 2017
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.195460E+03
main loop : -0.841473-233
epilogue : 0.212835E+03
total : 0.173747E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
>>> JOB COMPLETED AT Tue Nov 21 13:05:49 2017 <<<
atom: 1 xyz: 2(+) wall time: 212.8 date: Tue Nov 21 13:05:49 2017
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.212841E+03
main loop : -0.841473-233
epilogue : 0.230654E+03
total : 0.178130E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
>>> JOB COMPLETED AT Tue Nov 21 13:06:07 2017 <<<
atom: 1 xyz: 2(-) wall time: 230.7 date: Tue Nov 21 13:06:07 2017
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.230661E+03
main loop : -0.841473-233
epilogue : 0.248261E+03
total : 0.176007E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
>>> JOB COMPLETED AT Tue Nov 21 13:06:25 2017 <<<
atom: 1 xyz: 3(+) wall time: 248.3 date: Tue Nov 21 13:06:25 2017
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.248267E+03
main loop : -0.841473-233
epilogue : 0.266144E+03
total : 0.178767E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
>>> JOB COMPLETED AT Tue Nov 21 13:06:43 2017 <<<
atom: 1 xyz: 3(-) wall time: 266.1 date: Tue Nov 21 13:06:43 2017
Total PSPW energy : 0.0000000000E+00
== Timing ==
cputime in seconds
prologue : -0.266149E+03
main loop : -0.841473-233
epilogue : 0.283454E+03
total : 0.173048E+02
cputime/step: -Infinity ( 0 evalulations, 0 linesearches)
>>> JOB COMPLETED AT Tue Nov 21 13:07:00 2017 <<<
finite difference hessian delta = 1.0000000000000000E-002
zero matrix
finite difference derivative dipole; delta = 1.0000000000000000E-002
X vector of derivative dipole (au) [debye/angstrom]
d_dipole_x/ = 0.0000 [ 0.0000]
d_dipole_x/ = 0.0000 [ 0.0000]
d_dipole_x/ = 0.0000 [ 0.0000]
Y vector of derivative dipole (au) [debye/angstrom]
d_dipole_y/ = 0.0000 [ 0.0000]
d_dipole_y/ = 0.0000 [ 0.0000]
d_dipole_y/ = 0.0000 [ 0.0000]
Z vector of derivative dipole (au) [debye/angstrom]
d_dipole_z/ = 0.0000 [ 0.0000]
d_dipole_z/ = 0.0000 [ 0.0000]
d_dipole_z/ = 0.0000 [ 0.0000]
triangle hessian written to ./pspw-pbe-K1-79089.hess
derivative dipole written to ./pspw-pbe-K1-79089.fd_ddipole
Deleting state for pspw with suffix hess
./pspw-pbe-K1-79089.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
--------------------------------------------------------------------------
K 1 2.2204460D-16 0.0000000D+00 0.0000000D+00 3.8963710D+01
--------------------------------------------------------------------------
----------------------------------------------------
MASS-WEIGHTED NUCLEAR 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
-------------------------------------------------
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 = 36.892 cal/mol-K
- Translational = 36.892 cal/mol-K (mol. weight = 38.9637)
- 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: 123.3s wall: 123.3s
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 267 12
current total bytes 0 0
maximum total bytes 213040040 48482024
maximum total K-bytes 213041 48483
maximum total M-bytes 214 49
CITATION
--------
Please cite the following reference when publishing
results obtained with NWChem:
M. Valiev, E.J. Bylaska, N. Govind, K. Kowalski,
T.P. Straatsma, H.J.J. van Dam, D. Wang, J. Nieplocha,
E. Apra, T.L. Windus, W.A. de Jong
"NWChem: a comprehensive and scalable open-source
solution for large scale molecular simulations"
Comput. Phys. Commun. 181, 1477 (2010)
doi:10.1016/j.cpc.2010.04.018
AUTHORS
-------
E. Apra, E. J. Bylaska, W. A. de Jong, N. Govind, K. Kowalski,
T. P. Straatsma, M. Valiev, H. J. J. van Dam, D. Wang, T. L. Windus,
J. Hammond, J. Autschbach, K. Bhaskaran-Nair, J. Brabec, K. Lopata,
S. A. Fischer, S. Krishnamoorthy, W. Ma, M. Klemm, O. Villa, Y. Chen,
V. Anisimov, F. Aquino, S. Hirata, M. T. Hackler, V. Konjkov, 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. 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, L. Jensen, M. Swart, Q. Wu, T. Van Voorhis,
A. A. Auer, M. Nooijen, L. D. Crosby, E. Brown, G. Cisneros, G. I. Fann,
H. Fruchtl, J. Garza, K. Hirao, R. 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: 283.5s wall: 283.5s
##################### end nwoutput #######################
All requests to Arrows were successful.
KEYWORDs -
reaction: :reaction
chinese_room: :chinese_room
molecule: :molecule
nmr: :nmr
predict: :predict
submitesmiles: :submitesmiles
nosubmitmissingesmiles
resubmitmissingesmiles
submitmachines: :submitmachines
useallentries
nomodelcorrect
eigenvalues: :eigenvalues
frequencies: :frequencies
nwoutput: :nwoutput
xyzfile: :xyzfile
alleigs: :alleigs
allfreqs: :allfreqs
reactionenumerate:
energytype:[erxn(gas) hrxn(gas) grxn(gas) delta_solvation grxn(aq)] :energytype
energytype:[kcal/mol kj/mol ev cm-1 ry hartree au] :energytype
tablereactions:
reaction: ... :reaction
reaction: ... :reaction
...
:tablereactions
tablemethods:
method: ... :method
method: ... :method
...
:tablemethods
:reactionenumerate
rotatebonds
xyzinput:
label: :label
xyzdata:
... xyz data ...
:xyzdata
:xyzinput
submitHf: :submitHf
nmrexp: :nmrexp
findreplace: old text | new text :findreplace
listnwjobs
fetchnwjob: :fetchnwjob
pushnwjob: :pushnwjob
printcsv: :printcsv
printeig: :printeig
printfreq: :printfreq
printxyz: :printxyz
printjobinfo: :printjobinfo
printnwout: :printnwout
badids: :badids
hup_string:
database:
table:
request_table:
listallesmiles
queuecheck
This software service and its documentation were developed at the Environmental Molecular Sciences Laboratory (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, and Department of Defense environmental science and technology program (SERDP). THE SOFTWARE SERVICE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE SERVICE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE SERVICE.