<div dir="ltr"><div>Dear Developer,</div><div><br></div><div>I was running vc-relax of Hf3N4. The input parameters are for body centered cubic crystal: space group I-43d (group 220 of the international table of crystallography). I used this space group for the positions of Hf and N atoms. The main thing I want to ask is: after vc-relax the CELL_PARAMETERS values which are for ibrav= -3 changed such that all diagonal elements are same but different to all non-diagonal elements ,however, the non-diagonal elements are same to each other. I think, if they all were same then that would mean the unit cell is either shrink or expanded which would make it easy to calculated the final (relaxed) lattice constant. But, in this case where it has different numbers for diagonal and non-diagonal, how do I interpret it? Is it still a cubic lattice? or it is trigonal now? I think I can not change the atomic positions otherwise I would lose the symmetry. What do I do now? Is it true that vc-relax may change the atomic positions but it will not change the crystal symmetry or the space group?<br></div><div><br></div><div>Input: <br></div><div>  <b><i>&control<br>    calculation = 'vc-relax',<br>    restart_mode = 'from_scratch',<br>    forc_conv_thr = 1.0D-5,<br>    tstress = .true.,<br>    tprnfor = .true.,<br>    prefix = 'Hf3N4_vr3' ,<br>    pseudo_dir = '/home/z5126106/QEspresso_bin/Pseudopotentials',<br>    outdir = '/share/scratch/z5126106/TempFiles',<br>    verbosity = 'high'<br>    !disk_io = "none"<br>  /<br>  &system<br>    <span style="color:rgb(255,0,0)">ibrav = 0,  <br>    celldm(1) = 12.69,              !materialsproject</span><br>    nat = 14,<br>    ntyp = 2,<br>    ecutwfc = 60,<br>    ecutrho = 720      <br>    nbnd = 80<br>    occupations = 'fixed'<br>  /<br>  &electrons<br>    conv_thr = 1.0d-11,<br>    mixing_mode = 'plain',<br>    mixing_ndim = 10,<br>    mixing_beta = 0.3<br>  /<br>  &ions<br>    ion_dynamics =  "bfgs",<br>    bfgs_ndim    = 3<br>  /<br>  &cell<br>    cell_dynamics = "bfgs",<br>    press_conv_thr = 0.1D0,<br>    cell_dofree = "all"<br>    cell_factor = 5<br>  /<br>CELL_PARAMETERS alat  ! ibrav = -3<br><span style="color:rgb(255,0,0)"> -0.5000 0.5000 0.5000<br> 0.5000 -0.5000 0.5000<br> 0.5000 0.5000 -0.5000</span><br>ATOMIC_SPECIES<br> Hf 178.49 Hf.pbe-spn-kjpaw_psl.1.0.0.UPF<br> N 14.0067 N.pbe-n-kjpaw_psl.1.0.0.UPF <br>ATOMIC_POSITIONS alat      !crystal<br> Hf       0.375000000   0.000000000   0.250000000 0 0 0   <br>Hf       0.125000000   0.000000000   0.750000000 0 0 0<br>Hf       0.250000000   0.375000000   0.000000000 0 0 0 <br>Hf       0.750000000   0.125000000   0.000000000 0 0 0<br>Hf       0.000000000   0.250000000   0.375000000 0 0 0<br>Hf       0.000000000   0.750000000   0.125000000 0 0 0<br>N        0.099800000   0.099800000   0.099800000 <br>N        0.400200000   -0.09980000   0.599800000 <br>N        -0.09980000   0.599800000   0.400200000 <br>N        0.599800000   0.400200000  -0.099800000 <br>N        0.349800000   0.349800000   0.349800000 <br>N        0.150200000   0.650200000   0.849800000 <br>N        0.849800000   0.150200000   0.650200000 <br>N        0.650200000   0.849800000   0.150200000 <br>K_POINTS automatic<br> 6 6 6 1 1 1</i></b></div><div><br></div><div>Output:</div><div><br></div><div><b><i>Begin final coordinates<br>     new unit-cell volume =   1718.69533 a.u.^3 (   254.68437 Ang^3 )<br><br>CELL_PARAMETERS (alat= 12.69000000)<br>  <span style="color:rgb(255,0,0)">-0.426131582   0.607025249   0.607025249<br>   0.607025249  -0.426131582   0.607025249<br>   0.607025249   0.607025249  -0.426131582</span><br><br>ATOMIC_POSITIONS (crystal)<br>Hf       0.375000000  -0.000000000   0.250000000    0   0   0<br>Hf       0.125000000   0.000000000   0.750000000    0   0   0<br>Hf       0.250000000   0.375000000   0.000000000    0   0   0<br>Hf       0.750000000   0.125000000  -0.000000000    0   0   0<br>Hf      -0.000000000   0.250000000   0.375000000    0   0   0<br>Hf       0.000000000   0.750000000   0.125000000    0   0   0<br>N        0.048570739   0.048570739   0.048570739<br>N        0.146441100  -0.104982446   0.427330131<br>N       -0.104982446   0.427330131   0.146441100<br>N        0.427330131   0.146441100  -0.104982446<br>N        0.309957583   0.309957583   0.309957583<br>N        0.262463403   0.707381916   1.020329856<br>N        1.020329856   0.262463403   0.707381916<br>N        0.707381916   1.020329856   0.262463403<br>End final coordinates<br><br><br><br>     A final scf calculation at the relaxed structure.<br>     The G-vectors are recalculated for the final unit cell<br>     Results may differ from those at the preceding step.<br> </i></b></div><div>Regards,<b><i><br></i></b></div><div><br></div><div></div><div><br></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><b><font size="2">Bharat Thapa</font></b></div></div></div></div></div></div></div></div>