NWChem/web/support/faq/md.html
2006-01-12 18:52:52 +00:00

122 lines
3.7 KiB
HTML

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"
"http://www.w3.org/TR/html4/loose.dtd">
<HTML>
<HEAD>
<TITLE>NWChem Frequently Asked Questions</TITLE>
<link rel="stylesheet" type="text/css" href="../../shared/nwchem_basic.css">
</HEAD>
<body text="#000000" bgcolor="#FFFFFF" link="blue" alink="blue" vlink="blue">
<table width="650"><tr><td>
<br>
<table>
<tbody>
<tr>
<td align="left" width="325"><img src="../../images/nwchem_logo_dark.gif" border=0 width=200 alt="NWChem - computational chemistry on parallel computers"></td>
<td align="right" width="325" valign="bottom">
<a href="../../nwchem_main.html">NWChem Home</a> |
<a href="../../disclaimer.html">Security & Privacy</a> |
<a href="http://www.pnl.gov" target="_blank">PNNL</a>
</td>
</tr>
</tbody>
</table>
<hr>
<table width="650">
<tbody>
<tr>
<td width="350" align="left"></td>
<td width="300" align="right"><a href="dirdyvtst.html">Next FAQ Catagory</a> | <a href="NWChem_FAQ.html">Return to Main FAQ</a></td>
</tr>
</tbody>
</table>
<br><br>
<h1>NWChem FAQ</h1>
<p><h2>Molecular Dynamics</h2></p>
<!--DO NOT REMOVE Begin Question and End Question Comments.-->
<ul>
<!--Begin Question-->
<li><a href="#0">How do I restart a Molecular Dynamics calculation?</a></li>
<!--End Question-->
</ul>
<hr><p>
<a name="0"></a>
<font color="purple">How do I restart a Molecular Dynamics calculation?</font>
<p>The restart file that was used to run the "task md dynamics" will contain the updated coordinates if you have a "record rest ###"
where ### is some number of steps (like 500). The restart file will have the prefix of the system name used in the md input block.
For example,
<pre>
# Start warming solute: initially at 50 K
# This file already exists: QJDa_md1.rst"
md
system QJDa_md1
cutoff 0.9
pme grid 64 order 4
leapfrog equil 0 data 10000 step 0.001
print step 1000 stat 10000
isotherm 50 trelax 0.1 0.1
update pairs 10 center 100 motion 100
record prop 50
record rest 500
record scoor 100
load pairs 20 size 0.75
end
task md dynamics
# Continue warming solute: at 100 K
task shell "cp QJDa_md1.rst QJDa_md2.rst"
md
system QJDa_md2
cutoff 0.9
pme grid 64 order 4
leapfrog equil 0 data 10000 step 0.001
print step 1000 stat 10000
isotherm 100 trelax 0.1 0.1
update pairs 10 center 100 motion 100
record prop 50
record rest 500
record scoor 100
load pairs 20 size 0.75
end
task md dynamics
</pre>
The QJDa_md2.rst file will have the final restart information.
</p>
<p>
<hr>
<font size="0">
<center>
<a href="../../nwchem_main.html">NWChem</a> |
<a href="../../capabilities/nwchem_capab.html">Capabilities</a> |
<a href="../../platforms/platforms.html">Platforms</a> |
<a href="../../download.html">Download</a> |
<a href="../../doc/user/index.html">User's Manual</a> |
<a href="../../doc/prog/index.html">Programmer's Manual</a> |
<a href="../../release-notes/index.html">Release Notes</a> |
<a href="NWChem_FAQ.html">FAQ</a>
</center>
<br>
<center>
<a href="../known_bugs.html">Known Bugs</a> |
<a href="../support.html">Support</a> |
<a href="../../training/training.html">Tutorial</a> |
<a href="../../contributors.html">Contributors</a> |
<a href="../../benchmarks/index.html">Benchmarks</a> |
<a href="http://www.emsl.pnl.gov/forms/search_nwchem.html">Search</a> |
<a href="http://www.emsl.pnl.gov/docs/mssg/index.html" target="_blank">Mol Sci. Soft. Group</a> |
<a href="../../citation.html">Citation</a>
</center>
</font>
<hr>Contact: <a href="mailto:ms3distribution@emsl.pnl.gov">NWChem Support</a><br>
Updated: February 22, 2005
</td></tr>
</table>
</body>
</html>