# Set up lattice
vector_a = [5.13498110527, 0.0, 0.0]*Angstrom
vector_b = [2.56749055263, 4.61963940331, 1.72948648768]*Angstrom
vector_c = [2.56749055263, 4.61963940331, -1.72948648768]*Angstrom
lattice = UnitCell(vector_a, vector_b, vector_c)

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

# Define coordinates
fractional_coordinates = [[ 0.002326096083,  0.014907734955,  0.014907734955],
                          [-0.330024007974,  0.352824332727,  0.352824332727],
                          [ 0.339245858849,  0.184651811908,  0.184651811908],
                          [ 0.334964646014,  0.062786363154,  0.062786363154],
                          [ 0.08886232479 ,  0.311877655258,  0.311877655258],
                          [ 0.591462796298,  0.187277838684,  0.187277838684]]

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