<div dir="ltr">Dear  Hongshen,<div><br></div><div>does it work when lda_plus_u = .false.? from your email it is not clear whether the error is due to the +U part or not. also, are you using the latest version of the code or an older one? at what point of the calculation does the error appear? is it in the calculation of forces, stresses or before?</div><div><br></div><div>Matteo</div><div><br><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 8, 2017 at 5:49 PM, Hongsheng Liu <span dir="ltr"><<a href="mailto:lhs.happy2007@gmail.com" target="_blank">lhs.happy2007@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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><br></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><br>I appreciate any help!<br></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></div></div>