# -------------------------------------------------------------
# AKMC
# -------------------------------------------------------------
saddle_search_parameters = SaddleSearchParameters(max_neb_images=5)
markov_chain = MarkovChain(bulk_configuration, TotalEnergy(bulk_configuration).evaluate())
akmc = AdaptiveKineticMonteCarlo(markov_chain,
                                 kmc_temperature=300.0*Kelvin,
                                 md_temperature=2000*Kelvin,
                                 calculator=calculator,
                                 saddle_search_parameters=saddle_search_parameters,
                                 write_searches=True,
                                 constraints=constraints,
                                 confidence=0.95)
akmc.run(100)