# Set up lattice
lattice = FaceCenteredCubic(4.04958*Angstrom)

# Define elements
elements = [Aluminium]

# Define coordinates
fractional_coordinates = [[ 0.,  0.,  0.]]

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