<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<pre wrap="">Dear Paolo:
Thank you for your response.
I change the input according to your suggestion.
The QE code runs but produces results that may not
be correct.
The complete QE input and part of output are
attached by the end of this message.
The crystal is determined by QE to be "ibrav=-13"
(see the output), which is an undocumented feature
in QE 5.3.0. (I assume it is like "ibrav=13" with
a different unqiue axis.)
The cif file from which I took the structure data
shows that the crystal has symmetry "C 1 2/c 1"
(monoclinic P) with unique axis along b, which should
correspond to "ibrav=-12" rather than "ibrav=13".
By the way, I also used cif2qe.sh to convert the
cif file to QE input, the results are correct.
The cif file is attached at the very end of this
message.
Could you please point out if I make a mistake
in the input?
Thank you very much. I truly appeciate your work.
Tsung-Lung Li
------------------------------------------------------
QE input file
------------------------------------------------------
&control
title = 'Gypsum (CaSO4.2H20 Z=4)'
prefix = 's-12'
outdir = '/work2/u2equa02/tempdir/s-12'
pseudo_dir = '/home/u2equa02/quantumhome/pseudo'
calculation = 'scf'
restart_mode = 'from_scratch'
disk_io = 'low'
verbosity = 'high'
wf_collect = .true.
tstress = .true.
tprnfor = .true.
etot_conv_thr = 1.0D-4
forc_conv_thr = 1.0D-3
/
&system
celldm(1) = 11.86181088389489 ! bohr
celldm(2) = 2.41851202803887
celldm(3) = 0.903616377250276
celldm(5) = -0.408489773780522
ntyp = 4
nat = 7
space_group = 15
uniqueb = .true.
origin_choice = 1
ecutwfc = 10
ecutrho = 80
/
&electrons
electron_maxstep = 200
conv_thr = 1.0D-6
startingpot = 'atomic'
startingwfc = 'atomic+random'
mixing_mode = 'plain'
mixing_beta = 0.5
mixing_ndim= 8
diagonalization = 'david'
/
ATOMIC_SPECIES
Ca 40.0780000000 Ca_pbe_v1.uspp.F.UPF
S 32.0650000000 S_pbe_v1.2.uspp.F.UPF
O 15.9994000000 O.pbe-n-kjpaw_psl.0.1.UPF
H 1.0079000000 H.pbe-rrkjus_psl.0.1.UPF
ATOMIC_POSITIONS crystal_sg
Ca 0.00000 0.17050 0.25000
S 0.00000 0.32727 0.75000
O 0.08319 0.27218 0.59103
O 0.19997 0.38195 0.91298
O -0.20823 0.06826 -0.07831
H -0.25800 0.08700 -0.23400
H -0.24400 0.02000 -0.07700
K_POINTS automatic
2 1 2 0 0 0
------------------------------------------------------
QE output file (partial)
------------------------------------------------------
bravais-lattice index = -13
lattice parameter (alat) = 11.8618 a.u.
unit-cell volume = 1664.6131 (a.u.)<sup class="moz-txt-sup"><span style="display:inline-block;width:0;height:0;overflow:hidden">^</span>3</sup>
number of atoms/cell = 24
number of atomic types = 4
number of electrons = 112.00
number of Kohn-Sham states= 56
kinetic-energy cutoff = 10.0000 Ry
charge density cutoff = 80.0000 Ry
convergence threshold = 1.0E-06
mixing beta = 0.5000
number of iterations used = 8 plain mixing
Exchange-correlation = PBE ( 1 4 3 4 0 0)
celldm(1)= 11.861811 celldm(2)= 2.418512 celldm(3)= 0.903616
celldm(4)= 0.000000 celldm(5)= -0.408490 celldm(6)= 0.000000
crystal axes: (cart. coord. in units of alat)
a(1) = ( 0.500000 -1.209256 0.000000 )
a(2) = ( 0.500000 1.209256 0.000000 )
a(3) = ( -0.369118 0.000000 0.824788 )
reciprocal axes: (cart. coord. in units 2 pi/alat)
b(1) = ( 1.000000 -0.413477 0.447531 )
b(2) = ( 1.000000 0.413477 0.447531 )
b(3) = ( 0.000000 0.000000 1.212434 )
------------------------------------------------------
cif file
------------------------------------------------------
data_global
_chemical_name_mineral 'Gypsum'
loop_
_publ_author_name
'Comodi P'
'Nazzareni S'
'Zanazzi P F'
'Speziale S'
_journal_name_full 'American Mineralogist'
_journal_volume 93
_journal_year 2008
_journal_page_first 1530
_journal_page_last 1537
_publ_section_title
;
High-pressure behavior of gypsum: A single-crystal X-ray study
Note: P = 0.0001 GPa
;
_database_code_amcsd 0004651
_chemical_compound_source 'Valle di Caramanico, Abruzzo, Italy'
_chemical_formula_sum 'Ca S O6 H4'
_cell_length_a 6.277
_cell_length_b 15.181
_cell_length_c 5.672
_cell_angle_alpha 90
_cell_angle_beta 114.11
_cell_angle_gamma 90
_cell_volume 493.340
_exptl_crystal_density_diffrn 2.318
_symmetry_space_group_name_H-M 'C 1 2/c 1'
loop_
_space_group_symop_operation_xyz
'x,y,z'
'1/2+x,1/2+y,z'
'x,-y,1/2+z'
'1/2+x,1/2-y,1/2+z'
'-x,y,1/2-z'
'1/2-x,1/2+y,1/2-z'
'-x,-y,-z'
'1/2-x,1/2-y,-z'
loop_
_atom_site_type_symbol
_atom_site_fract_x
_atom_site_fract_y
_atom_site_fract_z
_atom_site_U_iso_or_equiv
Ca 0.00000 0.17050 0.25000 0.01170
S 0.00000 0.32727 0.75000 0.00990
O 0.08319 0.27218 0.59103 0.01690
O 0.19997 0.38195 0.91298 0.01690
O -0.20823 0.06826 -0.07831 0.02410
H -0.25800 0.08700 -0.23400 0.03300
H -0.24400 0.02000 -0.07700 0.04400
</pre>
<div class="moz-cite-prefix">On 10/13/2016 09:35 PM, Paolo Giannozzi
wrote:<br>
</div>
<blockquote
cite="mid:CAPMgbCuw+BVAgGKiHhrEXj8Qtie-xB6RtXn4qKQVB6FvmpcVgA@mail.gmail.com"
type="cite">
<div dir="ltr">
<div>I don't think you need to specify "ibrav" (the code should
selecvt the proper one), but you need to specify celldm(2) and
celldm(3).<br>
<br>
</div>
Paolo</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Thu, Oct 13, 2016 at 2:39 PM,
Tsung-Lung Li <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:quantum@mail.ncyu.edu.tw" target="_blank">quantum@mail.ncyu.edu.tw</a>></span>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">Dear QE
Friends:<br>
<br>
I encounter a problem with the space_group option. A part of
my QE input is<br>
attached at the end of this message.<br>
<br>
The crystal I am working on has symmetry "C 1 2/c 1" (as
suggested by the<br>
cif file from which I take the structure data).<br>
<br>
Suggested by the following two webpages, I set "space_group
= 15".<br>
<a moz-do-not-send="true"
href="http://www.cryst.ehu.es/cgi-bin/cryst/programs/nph-wp-list?gnum=15&gua=b"
rel="noreferrer" target="_blank">http://www.cryst.ehu.es/cgi-<wbr>bin/cryst/programs/nph-wp-<wbr>list?gnum=15&gua=b</a><br>
<a moz-do-not-send="true"
href="http://img.chem.ucl.ac.uk/sgp/large/015ay1.htm"
rel="noreferrer" target="_blank">http://img.chem.ucl.ac.uk/sgp/<wbr>large/015ay1.htm</a><br>
<br>
Because the unique axis is along b in my problem,<br>
I take "ibrav = -12" (monoclinic).<br>
<br>
Running QE, I got the following error message.<br>
<br>
Error in routine input (1):<br>
Input ibrav not compatible with space group number<br>
<br>
There is not much information about space_group option in QE
manual.<br>
Could ould you please point out what is wrong with the
input?<br>
<br>
Thbak you.<br>
<br>
Tsung-Lung Li<br>
<br>
------------------------------<wbr>---------------------<br>
&system<br>
ibrav = -12<br>
celldm(1) = 11.86181088389489 ! bohr<br>
ntyp = 4<br>
nat = 7<br>
space_group = 15<br>
uniqueb = .true.<br>
origin_choice = 1<br>
....<br>
/<br>
<br>
ATOMIC_SPECIES<br>
Ca 40.0780000000 Ca_pbe_v1.uspp.F.UPF<br>
S 32.0650000000 S_pbe_v1.2.uspp.F.UPF<br>
O 15.9994000000 O.pbe-n-kjpaw_psl.0.1.UPF<br>
H 1.0079000000 H.pbe-rrkjus_psl.0.1.UPF<br>
<br>
ATOMIC_POSITIONS crystal_sg<br>
Ca 0.00000 0.17050 0.25000<br>
S 0.00000 0.32727 0.75000<br>
O 0.08319 0.27218 0.59103<br>
O 0.19997 0.38195 0.91298<br>
O -0.20823 0.06826 -0.07831<br>
H -0.25800 0.08700 -0.23400<br>
H -0.24400 0.02000 -0.07700<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
<br>
______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a moz-do-not-send="true"
href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a moz-do-not-send="true"
href="http://pwscf.org/mailman/listinfo/pw_forum"
rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
</font></span></blockquote>
</div>
<br>
<br clear="all">
<br>
-- <br>
<div class="gmail_signature" data-smartmail="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>Paolo Giannozzi, Dip. Scienze Matematiche
Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone +39-0432-558216, fax +39-0432-558222<br>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Pw_forum mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>
<a class="moz-txt-link-freetext" href="http://pwscf.org/mailman/listinfo/pw_forum">http://pwscf.org/mailman/listinfo/pw_forum</a></pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
Tsung-Lung Li, Ph. D.
Professor
Department of Electrophysics
National Chia-Yi University
300 Hsueh-Fu Road, Chiayi 60004, Taiwan
Phone: 886-5-2717904. FAX: 886-5-2717909.
<a class="moz-txt-link-abbreviated" href="mailto:E-mail:quantum@mail.ncyu.edu.tw">E-mail:quantum@mail.ncyu.edu.tw</a> </pre>
</body>
</html>