Molecular potential force constants

Z-matriccs arc commonly used as input to quantum mechanical ab initio and serai-empirical) calculations as they properly describe the spatial arrangement of the atoms of a molecule. Note that there is no explicit information on the connectivity present in the Z-matrix, as there is, c.g., in a connection table, but quantum mechanics derives the bonding and non-bonding intramolecular interactions from the molecular electronic wavefunction, starting from atomic wavefiinctions and a crude 3D structure. In contrast to that, most of the molecular mechanics packages require the initial molecular geometry as 3D Cartesian coordinates plus the connection table, as they have to assign appropriate force constants and potentials to each atom and each bond in order to relax and optimi-/e the molecular structure. Furthermore, Cartesian coordinates are preferable to internal coordinates if the spatial situations of ensembles of different molecules have to be compared. Of course, both representations are interconvertible. 

Once HyperChem calculates potential energy, it can obtain all of the forces on the nuclei at negligible additional expense. This allows for rapid optimization of equilibrium and transition-state geometries and the possibility of computing force constants, vibrational modes, and molecular dynamics trajectories. 

Each molecular mechanics method has its own functional form MM+, AMBER, OPES, and BIO+. The functional form describes the analytic form of each of the terms in the potential. Eor example, MM+ has both a quadratic and a cubic stretch term in the potential whereas AMBER, OPES, and BIO+ have only quadratic stretch terms. The functional form is referred to here as the force field. Eor example, the functional form of a quadratic stretch with force constant Kj. and equilibrium distance rQ is ... 

At the same time, many lattice dynamics models have been constructed from force-constant models or ab-initio methods. Recently, the technique of molecular dynamics (MD) simulation has been widely used" " to study vibrations, surface melting, roughening and disordering. In particular, it has been demonstrated " " " that the presence of adatoms modifies drastically the vibrational properties of surfaces. Lately, the dynamical properties of Cu adatoms on Cu(lOO) " and Cu(lll) faces have been calculated using MD simulations and a many-body potential based on the tight-binding (TB) second-moment aproximation (SMA). " ... 

Vibrational spectroscopy has played a very important role in the development of potential functions for molecular mechanics studies of proteins. Force constants which appear in the energy expressions are heavily parameterized from infrared and Raman studies of small model compounds. One approach to the interpretation of vibrational spectra for biopolymers has been a harmonic analysis whereby spectra are fit by geometry and/or force constant changes. There are a number of reasons for developing other approaches. The consistent force field (CFF) type potentials used in computer simulations are meant to model the motions of the atoms over a large ranee of conformations and, implicitly temperatures, without reparameterization. It is also desirable to develop a formalism for interpreting vibrational spectra which takes into account the variation in the conformations of the chromophore and surroundings which occur due to thermal motions. 

Terms representing these interactions essentially make up the difference between the traditional force fields of vibrational spectroscopy and those described here. They are therefore responsible for the fact that in many cases spectroscopic force constants cannot be transferred to the calculation of geometries and enthalpies (Section 2.3.). As an example, angle deformation potential constants derived for force fields which involve nonbonded interactions often deviate considerably from the respective spectroscopic constants (7, 7 9, 21, 22). Nonbonded interactions strongly influence molecular geometries, vibrational frequencies, and enthalpies. They are a decisive factor for the transferability of force fields between systems of different strain (Section 2.3.). 

Molecular mechanics force fields rest on four fundamental principles. The first principle is derived from the Bom-Oppenheimer approximation. Electrons have much lower mass than nuclei and move at much greater velocity. The velocity is sufficiently different that the nuclei can be considered stationary on a relative scale. In effect, the electronic and nuclear motions are uncoupled, and they can be treated separately. Unlike quantum mechanics, which is involved in determining the probability of electron distribution, molecular mechanics focuses instead on the location of the nuclei. Based on both theory and experiment, a set of equations are used to account for the electronic-nuclear attraction, nuclear-nuclear repulsion, and covalent bonding. Electrons are not directly taken into account, but they are considered indirectly or implicitly through the use of potential energy equations. This approach creates a mathematical model of molecular structures which is intuitively clear and readily available for fast computations. The set of equations and constants is defined as the force... 

It is finally assumed that with all force constants and potential functions correctly specified in terms of the electronic configuration of the molecule, the nuclear arrangement that minimizes the steric strain corresponds to the observed structure of the isolated (gas phase) molecule. In practice however, the adjustable parameters, in virtually all cases, are chosen to reproduce molecular structures observed by solid-state diffraction methods. The parameters are therefore conditioned by the crystal environment and the minimized structure corresponds to neither gas phase nor isolated molecule [109],... 

A potentially much more adaptable technique is force-field vibrational modeling. In this method, the effective force constants related to distortions of a molecule (such as bond stretching) are used to estimate unknown vibrahonal frequencies. The great advantage of this approach is that it can be applied to any material, provided a suitable set of force constants is known. For small molecules and complexes, approximate force constants can often be determined using known (if incomplete) vibrational specha. These empirical force-field models, in effect, represent a more sophisticated way of exhapolating known frequencies than the rule-based method. A simple type of empirical molecular force field, the modified Urey-Bradley force field (MUBFF), is introduced below. 

While Eq. (9.49) has a well-defined potential energy function, it is quite difficult to solve in the indicated coordinates. However, by a clever transfonnation into a unique set of mass-dependent spatial coordinates q, it is possible to separate the 3 Ai-dirncnsional Eq. (9.49) into 3N one-dimensional Schrodinger equations. These equations are identical to Eq. (9.46) in form, but have force constants and reduced masses that are defined by the action of the transformation process on the original coordinates. Each component of q corresponding to a molecular vibration is referred to as a normal mode for the system, and with each component there is an associated set of harmonic oscillator wave functions and eigenvalues that can be written entirely in terms of square roots of the force constants found in the Hessian matrix and the atomic masses. 

---