Hi,
I still not understand the tutorial. Do I need to set it in bulk configuration ? Can you please check it and tell me where i did my mistakes
?
# Set up lattice
vector_a = [10.0, 0.0, 0.0]*Angstrom
vector_b = [0.0, 6.60121396209, 0.0]*Angstrom
vector_c = [0.0, 0.0, 20.9820399218]*Angstrom
lattice = UnitCell(vector_a, vector_b, vector_c)
# Define elements
elements = [Carbon, Carbon, Carbon, Carbon, Carbon, Carbon, Carbon, Carbon,
Carbon, Carbon, Carbon, Carbon, Carbon, Carbon, Carbon, Carbon,
Carbon, Carbon, Carbon, Carbon, Carbon, Carbon, Carbon, Carbon,
Carbon, Carbon, Carbon, Carbon, Carbon, Carbon, Carbon, Carbon,
Carbon, Carbon, Carbon, Carbon, Carbon, Carbon, Carbon, Carbon,
Hydrogen, Hydrogen, Hydrogen, Hydrogen, Hydrogen, Hydrogen,
Hydrogen, Hydrogen, Hydrogen, Hydrogen, Hydrogen, Hydrogen,
Hydrogen, Hydrogen, Hydrogen, Hydrogen, Hydrogen, Hydrogen,
Hydrogen, Hydrogen]
# Define coordinates
fractional_coordinates = [[ 0.5 , 0.220392144016, 0.127551466439],
[ 0.5 , 0.220392144016, 0.330705212018],
[ 0.5 , 0.220392144016, 0.533858957596],
[ 0.5 , 0.220392144016, 0.737012703175],
[ 0.5 , 0.220392144016, 0.940166448754],
[ 0.5 , 0.406797381339, 0.161410424035],
[ 0.5 , 0.406797381339, 0.364564169614],
[ 0.5 , 0.406797381339, 0.567717915193],
[ 0.5 , 0.406797381339, 0.770871660772],
[ 0.5 , 0.406797381339, 0.974025406351],
[ 0.5 , 0.220392144016, 0.059833551246],
[ 0.5 , 0.220392144016, 0.262987296825],
[ 0.5 , 0.220392144016, 0.466141042404],
[ 0.5 , 0.220392144016, 0.669294787982],
[ 0.5 , 0.220392144016, 0.872448533561],
[ 0.5 , 0.593202618661, 0.127551466439],
[ 0.5 , 0.593202618661, 0.330705212018],
[ 0.5 , 0.593202618661, 0.533858957596],
[ 0.5 , 0.593202618661, 0.737012703175],
[ 0.5 , 0.593202618661, 0.940166448754],
[ 0.5 , 0.779607855984, 0.161410424035],
[ 0.5 , 0.779607855984, 0.364564169614],
[ 0.5 , 0.779607855984, 0.567717915193],
[ 0.5 , 0.779607855984, 0.770871660772],
[ 0.5 , 0.779607855984, 0.974025406351],
[ 0.5 , 0.406797381339, 0.025974593649],
[ 0.5 , 0.406797381339, 0.229128339228],
[ 0.5 , 0.406797381339, 0.432282084807],
[ 0.5 , 0.406797381339, 0.635435830386],
[ 0.5 , 0.406797381339, 0.838589575965],
[ 0.5 , 0.593202618661, 0.059833551246],
[ 0.5 , 0.593202618661, 0.262987296825],
[ 0.5 , 0.593202618661, 0.466141042404],
[ 0.5 , 0.593202618661, 0.669294787982],
[ 0.5 , 0.593202618661, 0.872448533561],
[ 0.5 , 0.779607855984, 0.025974593649],
[ 0.5 , 0.779607855984, 0.229128339228],
[ 0.5 , 0.779607855984, 0.432282084807],
[ 0.5 , 0.779607855984, 0.635435830386],
[ 0.5 , 0.779607855984, 0.838589575965],
[ 0.5 , 0.077393049027, 0.153526060088],
[ 0.5 , 0.077393049027, 0.356679805667],
[ 0.5 , 0.077393049027, 0.559833551246],
[ 0.5 , 0.077393049027, 0.762987296825],
[ 0.5 , 0.077393049027, 0.966141042404],
[ 0.5 , 0.077393049027, 0.033858957596],
[ 0.5 , 0.077393049027, 0.237012703175],
[ 0.5 , 0.077393049027, 0.440166448754],
[ 0.5 , 0.077393049027, 0.643320194333],
[ 0.5 , 0.077393049027, 0.846473939912],
[ 0.5 , 0.922606950973, 0.135435830386],
[ 0.5 , 0.922606950973, 0.338589575965],
[ 0.5 , 0.922606950973, 0.541743321544],
[ 0.5 , 0.922606950973, 0.744897067122],
[ 0.5 , 0.922606950973, 0.948050812701],
[ 0.5 , 0.922606950973, 0.051949187299],
[ 0.5 , 0.922606950973, 0.255102932878],
[ 0.5 , 0.922606950973, 0.458256678456],
[ 0.5 , 0.922606950973, 0.661410424035],
[ 0.5 , 0.922606950973, 0.864564169614]]
# Set up configuration
bulk_configuration = BulkConfiguration(
bravais_lattice=lattice,
elements=elements,
fractional_coordinates=fractional_coordinates
)
#----------------------------------------
# Basis Set
#----------------------------------------
basis_set = [
LDABasis.Hydrogen_SingleZeta,
LDABasis.Carbon_SingleZeta,
]
#----------------------------------------
# Exchange-Correlation
#----------------------------------------
exchange_correlation = LSDA.PZ
k_point_sampling = MonkhorstPackGrid(
na=1,
nb=1,
nc=100,
symmetries=[
([[ 1., 0., 0.],
[ 0., 1., 0.],
[ 0., 0., 1.]], [ 0., 0., 0.]),
([[-1., 0., 0.],
[ 0.,-1., 0.],
[ 0., 0.,-1.]], [ 0., 0., 0.]),
],
force_timereversal=True,
shift_to_gamma=[True, True, True],
)
numerical_accuracy_parameters = NumericalAccuracyParameters(
k_point_sampling=k_point_sampling,
density_mesh_cutoff=150.0*Rydberg,
)
calculator = LCAOCalculator(
basis_set=basis_set,
exchange_correlation=exchange_correlation,
numerical_accuracy_parameters=numerical_accuracy_parameters,
)
bulk_configuration.setCalculator(calculator)