The Story:
Reaction pathways can be difficult things to deal with using electronic structure methods. Intrinsic reaction coordinates (IRCs) can be determined a number of ways (a discussion of these methods lie way outside the scope of this post). Suffice to say, the reaction pathway is an 'energy' pathway of reaction, connecting reactants and products, usually through some intermediate and/or transition state (TS).
For the purposes of this question, imagine implementing some sort of IRC technique which begins with a previously resolved transition state (Figure 1). The computation follows the reaction coordinate in both the forward and reverse directions from the TS energetically downhill until it reaches a (for the sake of simplicity) a minimum energy point with features which satisfy the method used. In general, the structure and energy of this final IRC geometry should be very representative of the structure which is connected to the TS via the pathway determined. If you were to optimize this final geometry, you should get good agreement between the optimized structure and the IRC configuration.
Now consider this scenario. An optimization was performed on the IRC resolved structure B and the energy relaxation was significant and the structure significantly changed, giving rise to a new (and lower) energy. One could no longer (within the formalism of IRC methods) say that the optimized structure B lies on the reaction pathway, and therefore is not a representative configuration at point B on the reaction pathway associated with the TS.
The Question:
Does anyone know of an example in the literature where this type of behavior occurs? How does one 'deal' with a situation like this?
Other Information:
Let us assume the following:
That this scenario occurs with a few IRC procedures (not just one).
That the IRC computations were given a ridiculously large amount of cycles and steps
That various stepsizes were sampled in parallel IRC jobs with no significant difference in the results
(I will add more as I think of them)
Potentially Helpful Imagery
Imagine two somewhat complex molecules coming together to form a single covalent bond. You run an IRC from the TS (which describes the bond formation) in the reverse direction and the intermolecular separation begins to increase. This separation continues to increase until the IRC job terminates properly. Now imagine running an optimization on these two (slightly complex) molecules which are separated by maybe 4.0 Å. Obviously these molecules are going to interact because they are not infinitely separated. You observe the (slightly complex) molecules moving toward each other during the optimization procedure (rotating, twisting, gyrating along the way) until you now have a minimum energy structure which is much lower in energy than the resulting IRC computation. These two molecules interact via intermolecular interactions (no covalent bonds are formed between the molecules). This optimized structure is analogous to a reactant complex (as opposed to separated reactants) and has an energy a few kcal mol$^{-1}$ lower than the corresponding IRC geometry. How on earth does this fit in to the reaction pathway?