<div dir="ltr"><div dir="ltr"><div style="color:rgb(0,0,0);box-sizing:border-box;font-family:"Segoe UI",system-ui,"Apple Color Emoji","Segoe UI Emoji",sans-serif;font-size:14px;font-variant-ligatures:normal"><div>Hello,</div><div><br>I can't find a way to make ph.x 6.6 run correctly on my system.</div><div>It is a MAPbI3 orthorhombic cell containing 48 atoms (Pb, I, C, N, H), with a 2x2x4 grid of k-points. ph.x correctly lists all representations (last line printed is "Representation 144 1 modes - To be done"), it writes the .wfcXX files and the _ph0 directory in a couple of minutes, then gets stuck indefinitely (I tested up to one hour). High verbosity does not give further clues.</div><div>A similar system (half the size, 24 atoms) runs fine, starting the self-consistent calculation right after it lists the representations.</div><div>I tried changing the following:</div><div>1) crystal setting both specifying the space group and describing the unit cell with lattice vectors.</div><div>2) PPs: PAW/US/full relativistic/PBE/PBEsol</div><div>3) the q-grid, or just the gamma point</div><div>4) computer: the behavior is the same on a local cluster and Marconi100. In particular, I ran these tests on M100 with -np 64 and -npool 16</div><div>Can you give me any advice on a way to work around this issue, or on some other tests I could run to help pinpoint what the problem may be?</div><div>I attach the pw.x and ph.x input files.</div><div>&CONTROL<br>title ='prova'<br>calculation ='scf'<br>restart_mode='from_scratch'<br>nstep = 200<br>prefix = 'mapi'<br>pseudo_dir ='../pseudo/'<br>outdir ='./'<br>wf_collect = .TRUE.<br>/</div><div>&SYSTEM<br>a = 1.26230e+01<br>c = 6.31612e+00<br>ibrav = 6<br>nat = 48<br>ntyp = 5<br>ecutwfc = 40<br>ecutrho = 360<br>occupations ='smearing'<br>smearing ='marzari-vanderbilt'<br>degauss = 0.02<br>/</div><div>&ELECTRONS<br>conv_thr = 1.0d-6<br>mixing_beta = 0.7<br>electron_maxstep = 200<br>startingpot = 'atomic'<br>startingwfc = 'atomic+random'<br>/</div><div>ATOMIC_SPECIES<br>Pb 207.21 Pb.pbe-dn-kjpaw_psl.1.0.0.UPF<br>I 126.90 I.pbe-n-kjpaw_psl.1.0.0.UPF<br>C 12.01 C.pbe-n-kjpaw_psl.1.0.0.UPF<br>N 14.00674 N.pbe-n-kjpaw_psl.0.1.UPF<br>H 1.00794 H.pbe-kjpaw_psl.0.1.UPF</div><div>K_POINTS {automatic}<br>2 2 4 0 0 0</div><div>ATOMIC_POSITIONS {angstrom}<br>Pb 3.261162 3.037754 3.394201<br>I 0.105163 3.037754 3.420198<br>I 3.261162 3.037754 0.181197<br>C 0.105163 6.193274 5.863317<br>N 0.105163 6.193274 0.930200<br>H 0.105163 5.241276 1.267197<br>H 0.930160 0.357750 1.267197<br>H 5.591685 0.357750 1.267197<br>H 5.689684 5.466275 5.499319<br>H 1.098160 5.927275 5.500317<br>H 6.150684 0.874751 5.500317<br>I 3.261162 0.000000 3.394201<br>Pb 3.261162 9.349274 3.394201<br>I 0.105163 9.349274 3.420198<br>I 3.261162 9.349274 0.181197<br>C 0.105163 12.504794 5.863317<br>N 0.105163 12.504794 0.930200<br>H 0.105163 11.552796 1.267197<br>H 0.930160 6.669270 1.267197<br>H 5.591685 6.669270 1.267197<br>H 5.689684 11.777795 5.499319<br>H 1.098160 12.238795 5.500317<br>H 6.150684 7.186271 5.500317<br>I 3.261162 6.311520 3.394201<br>Pb 9.572682 3.037754 3.394201<br>I 6.416683 3.037754 3.420198<br>I 9.572682 3.037754 0.181197<br>C 6.416683 6.193274 5.863317<br>N 6.416683 6.193274 0.930200<br>H 6.416683 5.241276 1.267197<br>H 7.241680 0.357750 1.267197<br>H 11.903205 0.357750 1.267197<br>H 12.001204 5.466275 5.499319<br>H 7.409680 5.927275 5.500317<br>H 12.462204 0.874751 5.500317<br>I 9.572682 0.000000 3.394201<br>Pb 9.572682 9.349274 3.394201<br>I 6.416683 9.349274 3.420198<br>I 9.572682 9.349274 0.181197<br>C 6.416683 12.504794 5.863317<br>N 6.416683 12.504794 0.930200<br>H 6.416683 11.552796 1.267197<br>H 7.241680 6.669270 1.267197<br>H 11.903205 6.669270 1.267197<br>H 12.001204 11.777795 5.499319<br>H 7.409680 12.238795 5.500317<br>H 12.462204 7.186271 5.500317<br>I 9.572682 6.311520 3.394201</div><div><br>&INPUTPH<br>tr2_ph=1.0d-14,<br>prefix='mapi'<br>ldisp = .true.<br>nq1 = 8,<br>nq2 = 8,<br>nq3 = 8,<br>outdir='./'<br>fildyn='mapi.dynG',<br>/</div><div><br></div><div>----------------</div><div>Candida Pipitone</div><div>phD student</div><div>University of Palermo</div></div></div></div>