# Define lattice
lattice = FaceCenteredCubic(4.08*Angstrom)

# Define elements
elements = [Gold]

# Define coordinates
coordinates = [[2.0, 2.0, 2.0]]*Angstrom

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

# Define a a calculator
bulk_configuration.setCalculator(LCAOCalculator())

# Calculate and save the Bloch state with quantum number 5
bloch_state = BlochState(bulk_configuration, quantum_number=5,
                         k_point=[0.0, 0.5, 0.5])
nlsave('bloch_state.nc', bloch_state)