<div dir="ltr"><div><div><div>Hi,<br></div>  Thank you for your reply. It failed even with kpoints of 4*4*4. It seems to be not the problem with memory. I will check with gamma point only.<br><br></div>Best regards,<br></div>Hongsheng<br></div><div class="gmail_extra"><br><div class="gmail_quote">2017-02-08 17:59 GMT+01:00 Federico Iori <span dir="ltr"><<a href="mailto:federico.iori@u-psud.fr" target="_blank">federico.iori@u-psud.fr</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:arial,helvetica,sans-serif;font-size:12pt;color:#000000"><div>Hi.<br></div><div>If it succeed for the smaller system, I would guess it is a memory problem maybe.<br></div><div>For a vc-relax of the conventional larger cell, do you really need so many kpt?<br></div><div><br></div><div><span name="x"></span><table class="m_5708662016925063169mceItemTable"><tbody><tr><td><img alt="Université Paris-Sud" src="https://zimbra.u-psud.fr/public/z_logo_signature.png"></td><td valign="top"><div><span style="font-size:10pt"><b>Federico IORI</b></span></div><div><span style="font-size:10pt">Marie Curie Fellow<br>Laboratoire de Physique des Solides<br>Bâtiment 510 - Rue André Rivière<br>91400 Orsay</span></div></td></tr></tbody></table><div><br></div><span name="x"></span><br></div><hr id="m_5708662016925063169zwchr"><div style="color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt"><b>De: </b>"Hongsheng Liu" <<a href="mailto:lhs.happy2007@gmail.com" target="_blank">lhs.happy2007@gmail.com</a>><br><b>Para: </b>"PWSCF Forum" <<a href="mailto:pw_forum@pwscf.org" target="_blank">pw_forum@pwscf.org</a>><br><b>Enviados: </b>Miércoles, 8 de Febrero 2017 17:49:26<br><b>Asunto: </b>[Pw_forum] problem with DFT+U<div><div class="h5"><br><div><br></div><div dir="ltr"><div><div>Dear all,<br></div>  I'm trying to do a geometry optimization for bulk magnetite with PBE+U. The job succeeded with primitive cell (including 14 atoms), but failed with conventional cell (including 56 atoms). The error message is shown below,<br>  'Fatal error in PMPI_Bcast: Other MPI error, error stack:<br>  PMPI_Bcast(2434)........: MPI_Bcast(buf=0x998c5a0, count=3366, MPI_DOUBLE_PRECISION,   root=4, comm=0x84000002) failed<br>  MPIR_Bcast_impl(1807)...:<br>  MPIR_Bcast(1835)........:'<br><div><br></div></div><div>Similar question was asked before on pw_forum, but it seems no solution till now. <br></div><div> My input is listed below,<br> &CONTROL<br>  calculation    = "vc-relax",<br>......<br>/<br>&SYSTEM<br>  ibrav = 1<br>  A = 8.522631 <br>  nat = 56,<br>  ntyp = 3,<br>  ecutwfc = 42.0, <br>  ecutrho = 340.0,<br>  nspin = 2,<br>  starting_magnetization(1) = 0,<br>  starting_magnetization(2) = -0.25,<br>  starting_magnetization(3) = 0.25,<br>  occupations = 'smearing',<br>  smearing = 'gaussian',<br>  degauss = 0.0008,<br>  lda_plus_u = .TRUE. <br>  lda_plus_u_kind = 0<br>  Hubbard_U(1) = 0,<br>  Hubbard_U(2) = 3,<br>  Hubbard_U(3) = 3,<br>/<br>.........<br>/<br>&CELL<br>  cell_dofree = 'volume'<br>/<br>......<br>K_POINTS automatic                     <wbr>               <br>6 6 6 0 0 0     <br><div><br></div>I appreciate any help!<br></div></div><br></div></div>______________________________<wbr>_________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br><a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a></div><div><br></div></div></div><br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br></div>