neb_configuration = nlread('ci-optimized-neb-configuration.nc',NudgedElasticBand)[-1]

# -------------------------------------------------------------
# HTST Event
# -------------------------------------------------------------
htst_event = HTSTEvent(
    neb_configuration,
    finite_difference_method=Central,
    minimum_displacement=0.05*Angstrom,
    )
nlsave('htst-event.nc', htst_event)
nlprint(htst_event)