# -------------------------------------------------------------
# Bulk Configuration
# -------------------------------------------------------------

# Set up lattice
vector_a = [2.84772899324, 0.0, 0.0]*Angstrom
vector_b = [0.0, 2.84772899324, 0.0]*Angstrom
vector_c = [0.0, 0.0, 23.0545939284]*Angstrom
lattice = UnitCell(vector_a, vector_b, vector_c)

# Define elements
elements = [Palladium, Palladium, Palladium, Palladium, Oxygen, Carbon]

# Define coordinates
fractional_coordinates = [[ 0.5           ,  0.5           ,  0.304219128903],
                          [ 0.            ,  0.            ,  0.391561742195],
                          [ 0.5           ,  0.5           ,  0.478552505581],
                          [ 0.            ,  0.            ,  0.565044244745],
                          [ 0.348570877238,  0.526851983339,  0.650146521047],
                          [ 0.351157010507,  0.52673551576 ,  0.594241652772]]

# Set up configuration
bulk_configuration = BulkConfiguration(
    bravais_lattice=lattice,
    elements=elements,
    fractional_coordinates=fractional_coordinates
    )