Thank you for your help. Can you please explain how to get the “periodic equivalent”? for example, for the “basis 0.666 0.333 -0.166”, is the "basis 0.666 0.333 1-0.166" its periodic equivalent? sorry, my background is not physics or materials sceince.
The 12 atom basis you were using does not look like wurtzite. Below is
from an old LAMMPS script, but you should be able to modify it to
match the latest LAMMPS syntax:
# Set up 8 atom CdSe wurtzite unit cell,
variable 1_3 equal div(1,3)
variable 2_3 equal div(2,3)
variable 1_6 equal div(1,6)
variable 5_6 equal div(5,6)
variable 1_12 equal div(1,12)
variable 5_12 equal div(5,12)
lattice user 1.0 &
a1 $a 0.0 0.0 &
a2 0.0 $b 0.0 &
a3 0.0 0.0 c &
basis 0\.0 0\.0 0\.0 &
basis 0\.5 0\.5 0\.0 &
basis {1_3} 0.0 0.5 &
basis \{5\_6\} 0\.5 0\.5 &
basis 0\.0 0\.0 0\.625 &
basis 0\.5 0\.5 0\.625 &
basis {1_3} 0.0 0.125 &
basis ${5_6} 0.5 0.125
Another questiuon is that " is it nessesary to use a prism region to build a
suppercell, say 20*10*10 unit cell, for a wurtzite structure?"
not automatically, usually a specific multiplicity is required. please
note that this is not really a question about LAMMPS but about
crystallography. there are *lots* of text books introducing concepts
like lattices and primitive vs. conventional vs. supercells. the
purpose of this mailing list is not to make up for your lack of
training in this matter. if you are serious about doing simulations
where this kind of knowledge is required, then it should go without
saying that you'll have to acquire the suitable basic training first.