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);"> &electrons</span><br>   emass = 200.d0,<br>   emass_cutoff = 3.0d0,<br>
   orthogonalization = 'ortho',<br>   electron_dynamics = 'damp',<br> /<br> &ions<br> ion_dynamics = 'none',<br>/<br> &cell<br>   cell_dynamics = 'none',<br>   press = 0.0d0,<br>-----------------------------------------<br>
 &electrons<br>   emass = 200.d0,<br>   emass_cutoff = 3.0d0,<br>   orthogonalization = 'ortho',<br>   electron_dynamics = 'damp',<br> /<br> &ions<br> ion_dynamics = 'none',<br>/<br> &cell<br>
   cell_dynamics = 'none',<br>   press = 0.0d0,<br>-----------------------------------------<br><span style="color: rgb(0, 0, 153);">2.randomization</span><br> &electrons<br>   emass = 200.d0,<br>  ortho_max=40<br>
   emass_cutoff = 3.0d0,<br>   orthogonalization = 'ortho',<br>electron_dynamics = 'sd', <br>/<br> &ions<br> ion_dynamics = 'none',<br>  tranp(2) = .TRUE.<br>  amprp(2) = 0.05<br>/<br> &cell<br>
   cell_dynamics = 'none',<br>   press = 0.0d0,<br>------------------------------------------<br><span style="color: rgb(51, 51, 153);">3. dynamics</span><br> &electrons<br>   emass = 200.d0,<br>   emass_cutoff = 3.0d0,<br>
   orthogonalization = 'ortho',<br>electron_dynamics = 'cg', <br>electron_temperature = 'nose',<br>ekincw=0.0001,<br>fnosee=0.5,<br>/<br> &ions<br>ion_dynamics = 'verlet',<br>  ion_temperature = 'nose',<br>
  tempw = 300.0,<br>  fnosep = 70<br>/<br> &cell<br>   cell_dynamics = 'pr',<br>   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>