"Pro question": There is a strange idea about the *OH radical. According to some ideas, those are organized in "layers". Do you know some method that could simulate different ways how *OH "approach" to molecules?

I didnt get you. If I understand you correctly, are you asking if we can incorporate OH radical in MD simulations? Maybe QM based MD techniques can incorporate reactions. Not in classical MD.