Exploring solvent effects upon the Menshutkin reaction using a polarizable force field

Area of Science:

• Computational Chemistry
• Chemical Kinetics
• Physical Chemistry

Background:

• The Menshutkin reaction is a fundamental SN2 reaction involving alkyl halides and tertiary amines.
• Accurate computation of reaction energetics, especially transition states, is crucial for understanding reaction mechanisms.
• Solvent effects play a significant role in reaction rates and mechanisms, particularly for polar transition states.

Purpose of the Study:

• To compute the energetics of the Menshutkin reaction between triethylamine and ethyl iodide.
• To evaluate the performance of different computational methods and basis sets.
• To investigate the impact of solvent models on the accuracy of reaction energetics.

Main Methods:

• Density functional theory (B3LYP) and Møller–Plesset perturbation theory (MP2) were employed.
• Various basis sets (LANL2DZ, MIDI!, 6-311G(d,p), aug-cc-PVTZ) and pseudopotentials were utilized.
• Quantum mechanics/molecular mechanics (QM/MM) Monte Carlo simulations with free-energy perturbation theory were used for solvent effects, employing PDDG/PM3 and OPLS force fields (nonpolarizable and polarizable OPLS-AAP).

Main Results:

• The B3LYP/MIDI! level of theory yielded the best activation free energies (DeltaG(++)), with a mean absolute error (MAE) of 4.9 kcal/mol.
• Nonpolarizable force fields underestimated relative rates in low-dielectric solvents like cyclohexane.
• Using a polarizable OPLS-AAP force field reduced the MAE to 3.1 kcal/mol, highlighting its importance for dipolar transition states.
• PDDG/PM3/OPLS and OPLS-AAP methods offered comparable accuracy to B3LYP at lower computational cost.

Conclusions:

• A fully polarizable force field is essential for accurately computing solvent effects on highly dipolar transition structures in low-dielectric media.
• Solute-solvent interactions significantly stabilize the emerging charge separation at the transition state, leading to rate accelerations.
• Computational methods like PDDG/PM3/OPLS-AAP provide a cost-effective alternative for studying Menshutkin reaction energetics.