lattice = BodyCenteredCubic( 3.509 * Ang)
li_bcc_bulk = BulkConfiguration(
    lattice,
    [Lithium],
    [ ( 0.0, 0.0, 0.0 ) * Ang ]
    )