Force field methods transition structure modelling

Ken pioneered the modeling of transition states with force field methods. Before modem tools existed for locating transition structures in all but the simplest reactions, his group used ab initio calculations to find the geometries of transition states and to determine force constants for distortions away from these preferred geometries. These force constants could then be used in standard molecular mechanics calculations, in order to predict how steric effects would affect the geometries and energies of the transition structures when substituent were present. 

With this information in hand, it seemed reasonable to attempt to use force field methods to model the transition states of more complex, chiral systems. To that end, transition state.s for the delivery of hydrogen atom from stannanes 69 71 derived from cholic acid to the 2.2,.3-trimethy 1-3-pentyl radical 72 (which was chosen as the prototypical prochiral alkyl radical) were modeled in a similar manner to that published for intramolecular free-radical addition reactions (Beckwith-Schicsscr model) and that for intramolecular homolytic substitution at selenium [32]. The array of reacting centers in each transition state 73 75 was fixed at the geometry of the transition state determined by ah initio (MP2/DZP) molecular orbital calculations for the attack of methyl radical at trimethyltin hydride (viz. rsn-n = 1 Si A rc-H = i -69 A 6 sn-H-C = 180°) [33]. The remainder of each structure 73-75 was optimized using molecular mechanics (MM2) in the usual way. In all, three transition state conformations were considered for each mode of attack (re or ) in structures 73-75 (Scheme 14). In general, the force field method described overestimates experimentally determined enantioseleclivities (Scheme 15), and the development of a flexible model is now being considered [33]. 

The rapid rise in computer speed over recent years has led to atom-based simulations of liquid crystals becoming an important new area of research. Molecular mechanics and Monte Carlo studies of isolated liquid crystal molecules are now routine. However, care must be taken to model properly the influence of a nematic mean field if information about molecular structure in a mesophase is required. The current state-of-the-art consists of studies of (in the order of) 100 molecules in the bulk, in contact with a surface, or in a bilayer in contact with a solvent. Current simulation times can extend to around 10 ns and are sufficient to observe the growth of mesophases from an isotropic liquid. The results from a number of studies look very promising, and a wealth of structural and dynamic data now exists for bulk phases, monolayers and bilayers. Continued development of force fields for liquid crystals will be particularly important in the next few years, and particular emphasis must be placed on the development of all-atom force fields that are able to reproduce liquid phase densities for small molecules. Without these it will be difficult to obtain accurate phase transition temperatures. It will also be necessary to extend atomistic models to several thousand molecules to remove major system size effects which are present in all current work. This will be greatly facilitated by modern parallel simulation methods that allow molecular dynamics simulations to be carried out in parallel on multi-processor systems [115]. 

---