Nonadiabatic histidine dissociation of hexacoordinate heme in neuroglobin protein.

In the present work, density functional theory and canonical nonadiabatic Monte Carlo transition state theory have been used to investigate the histidine dissociation process from hexacoordinate heme in Ngb protein. The potential energy surfaces (PES) of the lowest singlet, triplet, and quintet states are calculated by stepwise optimization along with the histidine dissociation pathway. Based on the calculated two-dimensional PES, the histidine dissociation rates for the spin-forbidden processes via singlet to triplet and singlet to quintet transitions have been calculated by the nonadiabatic Monte Carlo transition state theory in canonical ensemble. The present study provides a quantitative description on spin-forbidden histidine dissociation processes.