bulk_configuration = BulkConfiguration(
    bravais_lattice=BodyCenteredCubic(2.8665*Angstrom),
    elements=[Iron],
    cartesian_coordinates=[[ 0.,  0.,  0.]]*Angstrom
    )

numerical_accuracy_parameters = NumericalAccuracyParameters(
    k_point_sampling=(10, 10, 10),
    )

calculator = LCAOCalculator(
    exchange_correlation=LSDA.PZ,
    numerical_accuracy_parameters=numerical_accuracy_parameters,
    )

initial_spin = InitialSpin(scaled_spins=[1.0])

bulk_configuration.setCalculator(
    calculator,
    initial_spin=initial_spin,
)

mulliken_population = MullikenPopulation(bulk_configuration)
nlprint(mulliken_population)