sample qmd input from nwchem official page does not work


Just Got Here
Hi!
I tried to run the following input (which you can also find in this link, http://nwchemgit.github.io/index.php/Release66:Gaussian_Basis_AIMD#Sample_input),

start qmd_tddft_h2o_svr
echo
print low
geometry noautosym noautoz
 O   0.00000000    -0.01681748     0.11334792
H 0.00000000 0.81325914 -0.34310308
H 0.00000000 -0.67863597 -0.56441201
end
basis
 * library 6-31G*
end
dft
 xc pbe0
end
tddft
 nroots 5
notriplet
target 1
civecs
grad
root 1
end
end
qmd
 nstep_nucl  200
dt_nucl 10.0
com_step 10
thermostat none
print_xyz 5
end
task tddft qmd


However, I got this error,


Summary of "ao basis" -> "" (cartesian)


     Tag                 Description            Shells   Functions and Types

------------------------------ ------ ---------------------
  • 6-31G* on all atoms


tddft_input: filename not found; default name will be used


qmd_input: unknown directive 0




current input line :
27: com_step 10




There is an error in the input file


For more information see the NWChem manual at http://nwchemgit.github.io/index.php/NWChem_Documentation


For further details see manual section:


In principle, this input should work but does not. What is happening??

Best,

J

Forum Regular
Hi, I am unable to reproduce your issue with the input from our online manual.
Can you provide more information ?

Thanks.

Best,
-Niri

niri.govind@pnnl.gov

Just Got Here
Hi,
Find below the rest of the output,


          Memory information
------------------

   heap     =   13107198 doubles =    100.0 Mbytes
stack = 13107195 doubles = 100.0 Mbytes
global = 26214400 doubles = 200.0 Mbytes (distinct from heap & stack)
total = 52428793 doubles = 400.0 Mbytes
verify = yes
hardfail = no


          Directory information
---------------------

 0 permanent = .
0 scratch = .




                               NWChem Input Module
-------------------



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 O 8.0000 0.00000000 -0.01682581 0.11342109
2 H 1.0000 0.00000000 0.81325081 -0.34302991
3 H 1.0000 0.00000000 -0.67864430 -0.56433884

     Atomic Mass 
-----------

     O                 15.994910
H 1.007825


Effective nuclear repulsion energy (a.u.)       9.2887672039

           Nuclear Dipole moment (a.u.) 
----------------------------
X Y Z
---------------- ---------------- ----------------
0.0000000000 0.0000000000 0.0000000000


           XYZ format geometry
-------------------
3
geometry
O 0.00000000 -0.01682581 0.11342109
H 0.00000000 0.81325081 -0.34302991
H 0.00000000 -0.67864430 -0.56433884

==============================================================================
                               internuclear distances
------------------------------------------------------------------------------
center one | center two | atomic units | angstroms
------------------------------------------------------------------------------
2 H | 1 O | 1.79013 | 0.94730
3 H | 1 O | 1.79012 | 0.94729
------------------------------------------------------------------------------
number of included internuclear distances: 2
==============================================================================



==============================================================================
                                internuclear angles
------------------------------------------------------------------------------
center 1 | center 2 | center 3 | degrees
------------------------------------------------------------------------------
2 H | 1 O | 3 H | 105.51
------------------------------------------------------------------------------
number of included internuclear angles: 1
==============================================================================



 library name resolved from: environment
library file name is: </aplic/nwchem/6.5_ompi1.6.5_ifort12.1/data/libraries/>



Summary of "ao basis" -> "" (cartesian)
------------------------------------------------------------------------------
Tag Description Shells Functions and Types
---------------- ------------------------------ ------ ---------------------
* 6-31G* on all atoms


tddft_input: filename not found; default name will be used
------------------------------------------------------------------------
qmd_input: unknown directive 0
------------------------------------------------------------------------
------------------------------------------------------------------------
current input line :
27: com_step 10
------------------------------------------------------------------------
------------------------------------------------------------------------
There is an error in the input file
------------------------------------------------------------------------
For more information see the NWChem manual at http://nwchemgit.github.io/index.php/NWChem_Documentation


For further details see manual section:                                                                                                                                                                                    


I got an error file from the node,
0:qmd_input: unknown directive:Received an Error in Communication


MPI_ABORT was invoked on rank 0 in communicator MPI COMMUNICATOR 3 DUP FROM 0
with errorcode -1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.




mpirun has exited due to process rank 0 with PID 10405 on
node gatchan29 exiting improperly. There are two reasons this could occur:

1. this process did not call "init" before exiting, but others in
the job did. This can cause a job to hang indefinitely while it waits
for all processes to call "init". By rule, if one process calls "init",
then ALL processes must call "init" prior to termination.

2. this process called "init", but exited without calling "finalize".
By rule, all processes that call "init" MUST call "finalize" prior to
exiting or it will be considered an "abnormal termination"

This may have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).


Forum Vet
Jrt2018
I can see that you are using NWChem 6.5 from this output line

/aplic/nwchem/6.5_ompi1.6.5_ifort12.1


Unfortunately, the com_step input option was introduced in NWChem 6.6.

Therefore, the only solution to your problem is to move to a more recent NWChem version (6.8 being the latest so far)


Forum >> NWChem's corner >> Running NWChem