Bonds breaking force

Ea = Arrhenius activation energy Es = excess stress energy AEr = potential barrier for bond rotation Eel = molecular elastic energy F = mean force potential f = average force on the chain fb = bond breaking force H0 = Hookean spring constant kB = Boltzmann constant... 

Xsuda Y, Yasutake H, Ishijima A, Yanagida X (1996) Xorsional rigidity of single actin filaments and actin-actin bond breaking force under torsion measured directly by in vitro micromanipulation. Proc Natl Acad Sci USA 93 12937... 

In this chapter, we wiU review electrochemical electron transfer theory on metal electrodes, starting from the theories of Marcus [1956] and Hush [1958] and ending with the catalysis of bond-breaking reactions. On this route, we will explore the relation to ion transfer reactions, and also cover the earlier models for noncatalytic bond breaking. Obviously, this will be a tour de force, and many interesting side-issues win be left unexplored. However, we hope that the unifying view that we present, based on a framework of model Hamiltonians, will clarify the various aspects of this most important class of electrochemical reactions. 

Despite advent of theoretical methods and techniques and faster computers, no single theoretical method seems to be capable of reliable computational studies of reactivities of biocatalysts. Ab initio quantum mechanical (QM) methods may be accurate but are still too expensive to apply to large systems like biocatalysts. Semi-empirical quantum methods are not as accurate but are faster, but may not be fast enough for long time simulation of large molecular systems. Molecular mechanics (MM) force field methods are not usually capable of dealing with bond-breaking and formation... 

Estimates of the ultimate shear strength r0 can be obtained from molecular mechanics calculations that are applied to perfect polymer crystals, employing accurate force fields for the secondary bonds between the chains. When the crystal structure of the polymer is known, the increase in the energy can be calculated as a function of the shear displacement of a chain. The derivative of this function is the attracting force between the chains. Its maximum value represents the breaking force, and the corresponding displacement allows the calculation of the maximum allowable shear strain. In Sect. 4 we will present a model for the dependence of the strength on time and temperature. In this model a constant shear modulus g is used, thus r0=gyb. 

When addressing problems in computational chemistry, the choice of computational scheme depends on the applicability of the method (i.e. the types of atoms and/or molecules, and the type of property, that can be treated satisfactorily) and the size of the system to be investigated. In biochemical applications the method of choice - if we are interested in the dynamics and effects of temperature on an entire protein with, say, 10,000 atoms - will be to run a classical molecular dynamics (MD) simulation. The key problem then becomes that of choosing a relevant force field in which the different atomic interactions are described. If, on the other hand, we are interested in electronic and/or spectroscopic properties or explicit bond breaking and bond formation in an enzymatic active site, we must resort to a quantum chemical methodology in which electrons are treated explicitly. These phenomena are usually highly localized, and thus only involve a small number of chemical groups compared with the complete macromolecule. 

Molecular dynamics simulations, with quantum-mechanically derived energy and forces, can provide valuable insights into the dynamics and structure of systems in which electronic excitations or bond breaking processes are important. In these cases, conventional techniques with classical analytical potentials, are not appropriate. Since the quantum mechanical calculation has to be performed many times, one at each time step, the choice of a computationally fast method is crucial. Moreover, the method should be able to simulate electronic excitations and breaking or forming of bonds, in order to provide a proper treatment of those properties for which classical potentials fail. 

The dynamic NMR technique allows investigations on the rate of exchange between 3-substituted quinuclidinium ions and water. The rate of dissociation of amine/water (or amine/alcohol) complexes is determined71 by the free energy contribution from the pKa-dependent hydrogen bond breaking, and from dispersion forces between acceptor and donor which may be at the most 40% of the activation energy of the dissociation of the complex. Similar importance may be attributed to a term for the formation of a cavity prior to the dissociation of the complexes. 

---