To my understanding the nucleophilic attack of H2C=O is a two step process:
1) The slight positive charge on the carbon is open to attack by a nucleophile. The C=O loses its pi bond and reduces down to a sigma bond.
2) This enables the large negative charge on the pi orbital of O to attract a proton from its environment.
How would one best model the reaction pathway of this mechanism in Gaussian? I am working on a larger reaction pathway which includes a couple of these reactions. Should I attempt to search for a transition state of both steps happening simultaneously, or as two transition states with NuHC-(O-) as an intermediate? My nucleophile is lysine and the carbonyl molecule is linear glucose.
Very new to Gaussian, I would appreciate as much detail as possible.