<meta http-equiv="Content-Type" content="text/html; charset=GB18030"><div>Dear Giuseppe,</div><div><br></div><div>Thank you for your detailed answer which would be definitely helpful.</div><div><br></div><div>Jibiao Li</div><div style="position: relative;"><div><br></div><div><br></div><div style="font-size: 12px;font-family: Arial Narrow;padding:2px 0 2px 0;">------------------ Original ------------------</div><div style="font-size: 12px;background:#efefef;padding:8px;"><div><b>From:</b>                                                                                                                        "Quantum ESPRESSO users Forum"                                                                                    <giuseppe.mattioli@ism.cnr.it>;</div><div><b>Date:</b> Mon, Sep 19, 2022 08:39 PM</div><div><b>To:</b> "Quantum ESPRESSO users Forum"<users@lists.quantum-espresso.org>;<wbr></div><div></div><div><b>Subject:</b> Re: [QE-users] Error about incorrect atomic position for ESM</div></div><div><br></div><br>Dear Jibiao Li<br>Please read carefully the documentation, particularly when performing  <br>calculations very far from "plain vanilla" dft.<br><br> From the PW manual<br><br>                    'esm' :<br>                         Effective Screening Medium Method.<br>                         For polarized or charged slab calculation, embeds<br>                         the simulation cell within an effective semi-<br>                         infinite medium in the perpendicular direction<br>                         (along z). Embedding regions can be vacuum or<br>                         semi-infinite metal electrodes (use "esm_bc" to<br>                         choose boundary conditions). If between two<br>                         electrodes, an optional electric field<br>                         ("esm_efield") may be applied. Method described in<br>                         M. Otani and O. Sugino, "First-principles calculations<br>                         of charged surfaces and interfaces: A plane-wave<br>                         nonrepeated slab approach", PRB 73, 115407 (2006).<br><br>                         NB:<br>                            - Two dimensional (xy plane) average charge density<br>                              and electrostatic potentials are printed out to<br>                              'prefix.esm1'.<br><br>                            - Requires cell with a_3 lattice vector along z,<br>                              normal to the xy plane, with the slab centered<br>                              around z=0.<br><br>                            - For bc2 with an electric field and bc3 boundary<br>                              conditions, the inversion symmetry along  <br>z-direction<br>                              is automatically eliminated.<br><br>                            - In case of calculation='vc-relax', use<br>                              "cell_dofree"='2Dxy' or other parameters so that<br>                              c-vector along z-axis should not be  <br>moved. <------------------<br><br>and<br><br>    Variable:       esm_w<br><br>    Type:           REAL<br>    See:            assume_isolated<br>    Default:        0.d0<br>    Description:    If "assume_isolated" = 'esm', determines the  <br>position offset<br>                    [in a.u.] of the start of the effective screening region,<br>                    measured relative to the cell edge. (ESM region begins at<br>                    z = +/- [L_z/2 + esm_w] ).<br><br>In other words, you should first of all increase celldm(3); the  <br>*charge density* of your system should be contained inside a region  <br>between -c/4 and +c/4, being c the lattice parameter parallel to the z  <br>axis. You may provide an offset if your system is not placed across 0,  <br>but in my past experience you may sometimes encounter problems with  <br>the offset, thus I sugget (but maybe it is useless) to shift your  <br>system in order to place it centered around z=0.<br><br>If you are studying intercalation between graphene planes instead, and  <br>your system must be considered as strictly periodic along z, then you  <br>should not use esm at all.<br>HTH<br>Giuseppe<br><br>Quoting Jibiao Li <jibiaoli@foxmail.com>:<br><br>> Dear all,<br>><br>><br>> I try to relax a structure with ESM, but I got an error below<br>><br>><br>> &nbsp; &nbsp; &nbsp;number of k points=&nbsp; &nbsp; 12&nbsp;  <br>> Methfessel-Paxton smearing, width (Ry)=&nbsp; 0.0200<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; &nbsp;cart. coord. in units 2pi/alat<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp; 1) = (&nbsp;  <br>> &nbsp;0.0000000&nbsp; &nbsp;0.0000000&nbsp; &nbsp;0.0000000), wk  <br>> =&nbsp; &nbsp;0.1250000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp; 2) = (&nbsp;  <br>> &nbsp;0.0000000&nbsp; &nbsp;0.0000000&nbsp; &nbsp;0.4993936), wk  <br>> =&nbsp; &nbsp;0.2500000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp; 3) = (&nbsp;  <br>> &nbsp;0.0000000&nbsp; &nbsp;0.0000000&nbsp; -0.9987873), wk =&nbsp;  <br>> &nbsp;0.1250000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp; 4) = (&nbsp;  <br>> &nbsp;0.0000000&nbsp; -0.5773503&nbsp; &nbsp;0.0000000), wk =&nbsp;  <br>> &nbsp;0.1250000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp; 5) = (&nbsp;  <br>> &nbsp;0.0000000&nbsp; -0.5773503&nbsp; &nbsp;0.4993936), wk =&nbsp;  <br>> &nbsp;0.2500000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp; 6) = (&nbsp;  <br>> &nbsp;0.0000000&nbsp; -0.5773503&nbsp; -0.9987873), wk =&nbsp;  <br>> &nbsp;0.1250000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp; 7) = (&nbsp;  <br>> &nbsp;0.5000000&nbsp; -0.2886751&nbsp; &nbsp;0.0000000), wk =&nbsp;  <br>> &nbsp;0.1250000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp; 8) = (&nbsp;  <br>> -0.5000000&nbsp; -0.2886751&nbsp; &nbsp;0.0000000), wk =&nbsp;  <br>> &nbsp;0.1250000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp; 9) = (&nbsp;  <br>> &nbsp;0.5000000&nbsp; -0.2886751&nbsp; &nbsp;0.4993936), wk =&nbsp;  <br>> &nbsp;0.2500000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp;10) = (&nbsp;  <br>> -0.5000000&nbsp; -0.2886751&nbsp; &nbsp;0.4993936), wk =&nbsp;  <br>> &nbsp;0.2500000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp;11) = (&nbsp;  <br>> &nbsp;0.5000000&nbsp; -0.2886751&nbsp; -0.9987873), wk =&nbsp;  <br>> &nbsp;0.1250000<br>> &nbsp; &nbsp; &nbsp; &nbsp; k(&nbsp; &nbsp;12) = (&nbsp;  <br>> -0.5000000&nbsp; -0.2886751&nbsp; -0.9987873), wk =&nbsp;  <br>> &nbsp;0.1250000<br>><br>><br>> &nbsp; &nbsp; &nbsp;Dense&nbsp; grid:&nbsp; 1344013 G-vectors&nbsp;  <br>> &nbsp; &nbsp;FFT dimensions: ( 192, 192,&nbsp; 96)<br>><br>><br>> &nbsp; &nbsp; &nbsp;Smooth grid:&nbsp; &nbsp;389147 G-vectors&nbsp;  <br>> &nbsp; &nbsp;FFT dimensions: ( 120, 120,&nbsp; 60)<br>><br>><br>> &nbsp; &nbsp; &nbsp;Estimated max dynamical RAM per process  <br>> &gt;&nbsp; &nbsp; &nbsp;377.41 MB<br>><br>><br>> &nbsp; &nbsp; &nbsp;Estimated total dynamical RAM &gt;&nbsp; &nbsp;  <br>> &nbsp; 11.79 GB<br>><br>><br>> &nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>> &nbsp; &nbsp; &nbsp;Error in routine esm_check (1):<br>> &nbsp; &nbsp; &nbsp;incorrect atomic position for ESM<br>> &nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>><br>><br>> &nbsp; &nbsp; &nbsp;stopping ...<br>><br>><br>><br>> I examined the whole input file, but found no mistakes. The input  <br>> file is shown below and included as the attachment.&nbsp;<br>><br>><br>> Any idea to crack this problem?<br>><br>><br>><br>><br>> &nbsp;&amp;CONTROL<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; calculation  <br>> = 'vc-relax' ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; restart_mode  <br>> = 'from_scratch' ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; outdir = './' ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> pseudo_dir = '/home/jibiaoli/pseudo/PAW' ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; prefix = 'bt' ,<br>>            &nbsp; &nbsp; &nbsp; &nbsp;nstep = 199,&nbsp;<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp;tstress = .true. ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp;tprnfor = .true. ,<br>>             &nbsp; &nbsp; &nbsp; &nbsp; lfcp = .TRUE.,&nbsp; &nbsp; &nbsp;<br>> &nbsp;/<br>> &nbsp;&amp;SYSTEM<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; &nbsp;ibrav = 4,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp;celldm(1) = 28.03332990783,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp;celldm(3) = 0.5006071,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; &nbsp; &nbsp;nat = 80,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; &nbsp; ntyp = 4,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp;ecutwfc = 45 ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp;ecutrho = 411 ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp;input_dft = 'sla+pw+ggx+vdw1' ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp;occupations = 'smearing' ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp;degauss = 0.02D0 ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; smearing = 'methfessel-paxton' ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;assume_isolated = 'esm',<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; esm_bc = 'bc3',<br>> /<br>> &nbsp;&amp;ELECTRONS<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; electron_maxstep = 299,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp;mixing_beta = 0.2D0 ,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;diagonalization = 'david' ,<br>> &nbsp;/<br>> &nbsp;&amp;IONS<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp;ion_dynamics = 'bfgs' ,<br>> &nbsp;/<br>> &nbsp;&amp;CELL<br>> &nbsp;/<br>> &nbsp;&amp;FCP<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; fcp_mu = -5.219576,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; fcp_dynamics = 'bfgs',<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; fcp_tempw = 1,<br>> &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <br>> &nbsp; &nbsp; freeze_all_atoms= .FALSE.,<br>> &nbsp;/<br>> ATOMIC_SPECIES<br>> &nbsp; &nbsp; Br&nbsp; &nbsp;79.9100&nbsp; Br.pbe-dn-kjpaw_psl.1.0.0.UPF<br>> &nbsp; &nbsp; F&nbsp; &nbsp; 18.9980&nbsp; F.pbe-n-kjpaw_psl.1.0.0.UPF<br>> &nbsp; &nbsp; P&nbsp; &nbsp; 30.9740&nbsp; P.pbe-nl-kjpaw_psl.1.0.0.UPF&nbsp;<br>> &nbsp; &nbsp; C&nbsp; &nbsp; 12.0100&nbsp; C.pbe-n-kjpaw_psl.1.0.0.UPF<br>> ATOMIC_POSITIONS angstrom<br>> Br&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 4.8857064708&nbsp;  <br>> &nbsp; &nbsp; &nbsp;11.7768394079&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 3.9407468281<br>> F&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;3.8197174274&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.4597110794&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 4.7181716817<br>> F&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;6.1744957101&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.3394270012&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 4.7152037304<br>> F&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;4.8931069224&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.3705875468&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 4.5561838920<br>> P&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;4.9650890846&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 3.8015572354&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 3.7192840419<br>> F&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;3.7566375831&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 3.2595433439&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 2.7242934995<br>> F&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;6.1097595794&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 3.1430454823&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 2.7198097728<br>> F&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;5.0367554331&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 5.2311497524&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 2.8808884101<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0.0000000000&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.0000000000&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0000000000&nbsp; &nbsp; 0&nbsp; &nbsp;0&nbsp; &nbsp;0<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1.2361900607&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.7170517096&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0138171839<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2.4745523828&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.0073581771&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0448272197<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;3.7101907474&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.7241136134&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0572176935<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;4.9461994089&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.0101040361&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0863797199<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;6.1820578624&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.7251457619&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0828109698<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;7.4176619616&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.0089705107&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0883946780<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;8.6551903712&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.7202023957&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0603186182<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;9.8909505105&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.0053753049&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0427731456<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 11.1279529839&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.7166862952&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0121655509<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 12.3628679472&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.0029247779&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0025228560<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 13.5985969881&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 0.7147863623&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0110001690<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -1.2342764689&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.1429530459&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0181043859<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0.0024003497&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.8577643822&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0204839124<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1.2401865535&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.1458894037&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0013599653<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2.4761630563&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.8595809474&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0015451766<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;3.7122363408&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.1496027948&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0273663690<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;4.9464737546&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.8630747660&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0146793099<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;6.1806980931&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.1502829260&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0463756610<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;7.4161500403&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.8612203653&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0218290368<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;8.6518314190&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.1481555526&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0336305043<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;9.8901607936&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.8587098803&nbsp; &nbsp; &nbsp;  <br>> &nbsp;-0.0048730476<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 11.1267472325&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.1438196754&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0019971875<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 12.3636445079&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 2.8569763878&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0200815823<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -2.4711937241&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.2853488039&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0249202360<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -1.2338369649&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.9991066858&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0365309494<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0.0032481085&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.2854945862&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0352579279<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1.2415393747&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.9980020840&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0529725484<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2.4769033437&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.2857833200&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0429444112<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;3.7119127277&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.9978396952&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0675414101<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;4.9462507117&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.2856397140&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0404226256<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;6.1801544254&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.9984549476&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0632981383<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;7.4150498301&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.2864649840&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0325796465<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;8.6507650586&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.9981049219&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0451196154<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;9.8889676933&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.2857170240&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0212735505<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 11.1261803200&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 4.9991355997&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0312193338<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -3.7077605676&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 6.4273122556&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0502024657<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -2.4707707909&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 7.1413755217&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0579266693<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -1.2338320314&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 6.4271914837&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0527292585<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0.0035270392&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 7.1403619794&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0751485318<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1.2400808211&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 6.4260958758&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0786314435<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2.4766161323&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 7.1389471757&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1093129145<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;3.7110817006&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 6.4224940502&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1055736985<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;4.9461572831&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 7.1379282005&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1291761472<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;6.1814691185&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 6.4227956617&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1076923856<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;7.4164285695&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 7.1381818651&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1180161778<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;8.6530876314&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 6.4255094881&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0820559688<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;9.8895104371&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 7.1400398116&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0808612137<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -4.9443461491&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 8.5682115418&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1112338937<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -3.7077546465&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 9.2834166091&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1054738699<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -2.4708750334&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 8.5694590592&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0812215343<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -1.2340998158&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 9.2837141630&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0935732776<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0.0026467008&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 8.5687481228&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.0961884153<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1.2390053120&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 9.2823526931&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1257256553<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2.4752913601&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 8.5672098249&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1377526017<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;3.7104092166&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 9.2814910408&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1683633344<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;4.9455242075&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 8.5646106959&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1639022560<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;6.1810320547&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 9.2797968930&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1812835417<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;7.4170311281&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 8.5658219451&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1560512784<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;8.6534519955&nbsp;  <br>> &nbsp; &nbsp; &nbsp; 9.2807854615&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1528772627<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -6.1813449556&nbsp;  <br>> &nbsp; &nbsp; &nbsp;10.7086618913&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1767895913<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -4.9446616408&nbsp;  <br>> &nbsp; &nbsp; &nbsp;11.4244280124&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1633881046<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -3.7077392754&nbsp;  <br>> &nbsp; &nbsp; &nbsp;10.7114996002&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1300302244<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -2.4712557593&nbsp;  <br>> &nbsp; &nbsp; &nbsp;11.4275338898&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1223106464<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -1.2345466857&nbsp;  <br>> &nbsp; &nbsp; &nbsp;10.7120855486&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1099089414<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0.0019510066&nbsp;  <br>> &nbsp; &nbsp; &nbsp;11.4259908842&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1241121465<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1.2386816088&nbsp;  <br>> &nbsp; &nbsp; &nbsp;10.7107204043&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1391473407<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2.4748267037&nbsp;  <br>> &nbsp; &nbsp; &nbsp;11.4228996714&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1720290525<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;3.7096081115&nbsp;  <br>> &nbsp; &nbsp; &nbsp;10.7080418335&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.1894854994<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;4.9453339628&nbsp;  <br>> &nbsp; &nbsp; &nbsp;11.4206956373&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.2114298160<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;6.1802465576&nbsp;  <br>> &nbsp; &nbsp; &nbsp;10.7059532510&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.2075739383<br>> C&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;7.4159381224&nbsp;  <br>> &nbsp; &nbsp; &nbsp;11.4189623464&nbsp; &nbsp; &nbsp; &nbsp;  <br>> 0.2098980144<br>> K_POINTS automatic&nbsp;<br>> &nbsp; 2 2 4&nbsp; &nbsp;0 0 0&nbsp;<br>><br>><br>><br>><br>><br>><br>><br>> Jibiao Li<br>><br>> Department of Materials Science and Engineering<br>><br>> Yangtze Normal University<br>><br>> Juxian Avenue 16, Fuling, Chongqing, China 408100<br>><br>><br>><br>><br>><br>> &nbsp;<br><br><br><br>GIUSEPPE MATTIOLI<br>CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>Via Salaria Km 29,300 - C.P. 10<br>I-00015 - Monterotondo Scalo (RM)<br>Mob (*preferred*) +39 373 7305625<br>Tel + 39 06 90672342 - Fax +39 06 90672316<br>E-mail: <giuseppe.mattioli@ism.cnr.it><br><br>_______________________________________________<br>The Quantum ESPRESSO community stands by the Ukrainian<br>people and expresses its concerns about the devastating<br>effects that the Russian military offensive has on their<br>country and on the free and peaceful scientific, cultural,<br>and economic cooperation amongst peoples<br>_______________________________________________<br>Quantum ESPRESSO is supported by MaX (www.max-centre.eu)<br>users mailing list users@lists.quantum-espresso.org<br>https://lists.quantum-espresso.org/mailman/listinfo/users</div>