geometry optimization output coordinates

From NWChem

Viewed 2690 times, With a total of 4 Posts
Jump to: navigation, search

Clicked A Few Times
Threads 2
Posts 5
I am attempting to do a geometry optimization and output the optimized coordinates. I use the "xyz" flag, but the coordinate file is not created. The job runs through succesfully; it just doesn't create a coordinate file as I expect it to. My input file is:

start
md
 system test
noshake solute
end
basis
 * library 6-311G*
end
dft
 xc b3lyp
end
qmmm
 region all
density espfit
xyz foo_output_coords
end
task qmmm dft optimize

In the end, there are no files with the tag "foo_output_coords". Do I need an additional command to direct it to write out the coordinate file using that tag?

  • Marat Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
    Profile
    Send PM
Clicked A Few Times
Threads 2
Posts 43
It will not print xyz coordinates when you do all region optimization. Given, typically, large size of the system in QM/MM runs the xyz option was only meant for region qm. You can always extract PDB file at the end of the run.

Clicked A Few Times
Threads 2
Posts 5
Thank you for your help! What is the best way to extract the PDB file and which file should I extract it from? I know for output from the Molecular Dynamics module, PDB files can be extracted from the trajectory file (*.trj) using the Analysis module. However, no trajectory file is created from the qmmm dft geometry optimization calculation, and it's unclear to me in which output file the final optimized coordinates can be found. If necessary, I can create a PDB formatted file from raw xyz coordinates, if I just knew how to access the coordinates.

Clicked A Few Times
Threads 5
Posts 19
I don't know if it is the faster way, but to extract the optimized geometry I use a prepare module at the end of the input:

prepare
read file_name_ref.qrs
write file_name_opt.pdb
end
task prepare

Ale
Edited On 1:13:08 AM PDT - Wed, Oct 23rd 2013 by Varnon

Clicked A Few Times
Threads 2
Posts 5
That's great; thank you for your help!


Forum >> NWChem's corner >> QMMM



Who's here now Members 0 Guests 0 Bots/Crawler 1


AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC