Dear all<br>I performed CP molecular dynamics for a system that is metal with LDA and insulator with LDA+U.It contains Fe ions so I must use LDA+U for it .Now I want to obtain external pressure affect on it.These are my plan to BO MD:<br>
(dt=6)<br><span style="color: rgb(0, 0, 153);">1.electron minimization </span><br style="color: rgb(0, 0, 153);"><span style="color: rgb(0, 0, 153);">&nbsp;&amp;electrons</span><br>&nbsp;&nbsp; emass = 200.d0,<br>&nbsp;&nbsp; emass_cutoff = 3.0d0,<br>
&nbsp;&nbsp; orthogonalization = &#39;ortho&#39;,<br>&nbsp;&nbsp; electron_dynamics = &#39;damp&#39;,<br>&nbsp;/<br>&nbsp;&amp;ions<br>&nbsp;ion_dynamics = &#39;none&#39;,<br>/<br>&nbsp;&amp;cell<br>&nbsp;&nbsp; cell_dynamics = &#39;none&#39;,<br>&nbsp;&nbsp; press = 0.0d0,<br>-----------------------------------------<br>
&nbsp;&amp;electrons<br>&nbsp;&nbsp; emass = 200.d0,<br>&nbsp;&nbsp; emass_cutoff = 3.0d0,<br>&nbsp;&nbsp; orthogonalization = &#39;ortho&#39;,<br>&nbsp;&nbsp; electron_dynamics = &#39;damp&#39;,<br>&nbsp;/<br>&nbsp;&amp;ions<br>&nbsp;ion_dynamics = &#39;none&#39;,<br>/<br>&nbsp;&amp;cell<br>
&nbsp;&nbsp; cell_dynamics = &#39;none&#39;,<br>&nbsp;&nbsp; press = 0.0d0,<br>-----------------------------------------<br><span style="color: rgb(0, 0, 153);">2.randomization</span><br>&nbsp;&amp;electrons<br>&nbsp;&nbsp; emass = 200.d0,<br>&nbsp; ortho_max=40<br>
&nbsp;&nbsp; emass_cutoff = 3.0d0,<br>&nbsp;&nbsp; orthogonalization = &#39;ortho&#39;,<br>electron_dynamics = &#39;sd&#39;, <br>/<br>&nbsp;&amp;ions<br>&nbsp;ion_dynamics = &#39;none&#39;,<br>&nbsp; tranp(2) = .TRUE.<br>&nbsp; amprp(2) = 0.05<br>/<br>&nbsp;&amp;cell<br>
&nbsp;&nbsp; cell_dynamics = &#39;none&#39;,<br>&nbsp;&nbsp; press = 0.0d0,<br>------------------------------------------<br><span style="color: rgb(51, 51, 153);">3. dynamics</span><br>&nbsp;&amp;electrons<br>&nbsp;&nbsp; emass = 200.d0,<br>&nbsp;&nbsp; emass_cutoff = 3.0d0,<br>
&nbsp;&nbsp; orthogonalization = &#39;ortho&#39;,<br>electron_dynamics = &#39;cg&#39;, <br>electron_temperature = &#39;nose&#39;,<br>ekincw=0.0001,<br>fnosee=0.5,<br>/<br>&nbsp;&amp;ions<br>ion_dynamics = &#39;verlet&#39;,<br>&nbsp; ion_temperature = &#39;nose&#39;,<br>
&nbsp; tempw = 300.0,<br>&nbsp; fnosep = 70<br>/<br>&nbsp;&amp;cell<br>&nbsp;&nbsp; cell_dynamics = &#39;pr&#39;,<br>&nbsp;&nbsp; press = 0.0d0,<br><br><span style="color: rgb(204, 0, 0);">Is this way correct ?<br><span style="color: rgb(0, 0, 0);">I look forward to hearing from you?</span><br style="color: rgb(0, 0, 0);">
<span style="color: rgb(0, 0, 0);">Thanks</span><br style="color: rgb(0, 0, 0);"><span style="color: rgb(0, 0, 0);">Mansoureh</span><br style="color: rgb(204, 0, 0);"></span><br><br>