<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                                    <br>6 6 6 0 0 0     <br><br>I appreciate any help!<br></div></div>