Dissociative force field

Another important experimental observable is the rate of a chemical reaction. With a validated PES at hand, such rates can be determined in a variety of ways, ranging from transition state theory (TST) to computationally demanding fully QM rates based on flux correlation functions [99]. As MMPT is a dissociable force field, it can be used in all such formalisms. Furthermore, because the environmental modes are described with a conventional force field and therefore computationally inexpensive to evaluate, and because only the motion along the PT motif is more highly parameterized, MMPT is a particularly efficient empirical representation of the intermolecular interactions. 

Conformational Adjustments The conformations of protein and ligand in the free state may differ from those in the complex. The conformation in the complex may be different from the most stable conformation in solution, and/or a broader range of conformations may be sampled in solution than in the complex. In the former case, the required adjustment raises the energy, in the latter it lowers the entropy in either case this effect favors the dissociated state (although exceptional instances in which the flexibility increases as a result of complex formation seem possible). With current models based on two-body potentials (but not with force fields based on polarizable atoms, currently under development), separate intra-molecular energies of protein and ligand in the complex are, in fact, definable. However, it is impossible to assign separate entropies to the two parts of the complex. 

Pig. 4. Photo dissociation of ArHCl. Left hand side the number of force field evaluations per unit time. Right hand side the number of Fast-Fourier-transforms per unit time. Dotted line adaptive Verlet with the Chebyshev approximation for the quantum propagation. Dash-dotted line with the Lanczos iteration. Solid line stepsize controlling scheme based on PICKABACK. If the FFTs are the most expensive operations, PiCKABACK-like schemes are competitive, and the Lanczos iteration is significantly cheaper than the Chebyshev approximation. 

For each pair of interacting atoms (/r is their reduced mass), three parameters are needed D, (depth of the potential energy minimum, k (force constant of the par-tictilar bond), and l(, (reference bond length). The Morse ftinction will correctly allow the bond to dissociate, but has the disadvantage that it is computationally very expensive. Moreover, force fields arc normally not parameterized to handle bond dissociation. To circumvent these disadvantages, the Morse function is replaced by a simple harmonic potential, which describes bond stretching by Hooke s law (Eq. (20)). 

Fluorine ions form the first coordination sphere around the tantalum ion, which is the central atom of the complex. Potassium ions form the second coordination sphere, which significantly affects the geometry and force field of the first coordination sphere. The melting of K2TaF7 leads to the dissociation of the compound into ions, as follows ... 

After the first unsuccessful attempts to record a matrix IR spectrum of the methyl radical, reliable data were obtained by the use of the vacuum pyrolysis method. IR spectra of the radicals CH3 and CD3 frozen in neon matrices were measured among the products of dissociation of CH3I, (CH3)2Hg and CD3I (Snelson, 1970a). The spectra contained three absorptions at 3162 (1 3), 1396 V2) and 617 cm (I l) belonging to the radical CH3 and three bands 2381, 1026 and 463 cm assigned to the radical CD3. Normal coordinate analysis of these intermediates was performed and a valence force field calculated. In accordance with the calculations, methyl radical is a planar species having symmetry >31,. 

In a different approach to this problem, Brenner and Garrison used molecular dynamics to examine the chemical mechanisms which lead to reordering of the atom-pairing reconstruction during atom deposition . This simulation incorporated a dissociative valence-force field potentiaF and consisted essentially of a high-temperature anneal of monolayers of silicon atoms which had been deposited on a silicon (001) reconstructed surface. 

Another reactive force field that is dependent on bond-order was developed by van Duin, Dasgupta, Loran, and Goddard [183] for hydrocarbons. The configurational energy is described as the sum of energy contributions from internal modes as well as non-bonding van der Waals and Coulombic interactions, but the parameters of the functions that describe each contribution is dependent upon the bond order of atoms involved in each description. It is assumed that the bond order between an atom pair is dependent on the interatomic separation. While this model has been used to predict bond dissociation energies, heats of formation and structures of simple hydrocarbons, it was not applied to predict condensed phase properties. However, the form of the potential should allow for condensed phase studies. 

---