VALBOND

Last updated January 15, 2021

In molecular mechanics, VALBOND is a method for computing the angle bending energy that is based on valence bond theory.^[1] It is based on orbital strength functions, which are maximized when the hybrid orbitals on the atom are orthogonal. The hybridization of the bonding orbitals are obtained from empirical formulas based on Bent's rule, which relates the preference towards p character with electronegativity.

The VALBOND functions are suitable for describing the energy of bond angle distortion not only around the equilibrium angles, but also at very large distortions. This represents an advantage over the simpler harmonic oscillator approximation used by many force fields, and allows the VALBOND method to handle hypervalent molecules ^[2] and transition metal complexes.^[3]^[4] The VALBOND energy term has been combined with force fields such as CHARMM ^[5] and UFF to provide a complete functional form that includes also bond stretching, torsions, and non-bonded interactions.

Functional form

Non-hypervalent molecules

For an angle α between normal (non-hypervalent) bonds involving an sp^mdⁿ hybrid orbital, the energy contribution is

E(\alpha )=k(S^{max}-S(\alpha ))

,

where k is an empirical scaling factor that depends on the elements involved in the bond, S^max, the maximum strength function, is

S^{max}={\sqrt {\frac {1}{1+m+n}}}(1+{\sqrt {3m}}+{\sqrt {5n}})

and S(α) is the strength function

S(\alpha )=S^{max}{\sqrt {1-{\frac {1-{\sqrt {1-\Delta ^{2}}}}{2}}}}

which depends on the nonorthogonality integral Δ:

\Delta ={\frac {1}{1+m+n}}\left[1+m\cos \alpha +{\frac {n}{2}}(3\cos ^{2}\alpha -1)\right]

The energy contribution is added twice, once per each of the bonding orbitals involved in the angle (which may have different hybridizations and different values for k).

For non-hypervalent p-block atoms, the hybridization value n is zero (no d-orbital contribution), and m is obtained as %p(1-%p), where %p is the p character of the orbital obtained from

\%p_{i}={\frac {n_{p}wt_{i}}{\sum _{j}wt_{j}}}

where the sum over j includes all ligands, lone pairs, and radicals on the atom, n_p is the "gross hybridization" (for example, for an "sp²" atom, n_p = 2). The weight wt_i depends on the two elements involved in the bond (or just one for lone pair or radicals), and represents the preference for p character of different elements. The values of the weights are empirical, but can be rationalized in terms of Bent's rule.

Hypervalent molecules

For hypervalent molecules, the energy is represented as a combination of VALBOND configurations, which are akin to resonance structures that place three-center four-electron bonds (3c4e) in different ways. For example, ClF₃ is represented as having one "normal" two-center bond and one 3c4e bond. There are three different configurations for ClF₃, each one using a different Cl-F bond as the two-center bond. For more complicated systems the number of combinations increases rapidly; SF₆ has 45 configurations.

E_{tot}=\sum _{j}c_{j}E_{j}

where the sum is over all configurations j, and the coefficient c_j is defined by the function

c_{j}={\frac {\displaystyle \prod _{i=1}^{hype}\cos ^{2}\alpha _{i}}{\displaystyle \sum _{j=1}^{config}\prod _{i=1}^{hype}\cos ^{2}\alpha _{i}}}

where "hype" refers to the 3c4e bonds. This function ensures that the configurations where the 3c4e bonds are linear are favored.

The energy terms are modified by multiplying them by a bond order factor, BOF, which is the product of the formal bond orders of the two bonds involved in the angle (for 3c4e bonds, the bond order is 0.5). For 3c4e bonds, the energy is calculated as

E(\alpha )=BOF\times k_{\alpha }[1-\Delta (\alpha +\pi )^{2}]

where Δ is again the non-orthogonality function, but here the angle α is offset by 180 degrees (π radians).

Finally, to ensure that the axial vs equatorial preference of different ligands in hypervalent compounds is reproduced, an "offset energy" term is subtracted. It has the form

E_{offset}=\sum _{i=1}^{config}c_{i}\sum _{j=1}^{hype}{\frac {EN_{ija}+EN_{ijb}}{2}}

where the EN terms depend on the electronegativity difference between the ligand and the central atom as follows:

EN_{ija}=30\times (en_{lig}-en_{c.a.})\times ss

where ss is 1 if the electronegativity difference is positive and 2 if it is negative.

For p-block hypervalent molecules, d orbitals are not used, so n = 0. The p contribution m is estimated from ab initio quantum chemistry methods and a natural bond orbital (NBO) analysis.

Extension

More recent extensions, available in the CHARMM suite of codes, include the trans-influence (or trans effect) within VALBOND-TRANS^[5] and the possibility to run reactive molecular dynamics^[6] with "Multi-state VALBOND".^[7]

Related Research Articles

Tidal acceleration is an effect of the tidal forces between an orbiting natural satellite, and the primary planet that it orbits. The acceleration causes a gradual recession of a satellite in a prograde orbit away from the primary, and a corresponding slowdown of the primary's rotation. The process eventually leads to tidal locking, usually of the smaller first, and later the larger body. The Earth–Moon system is the best-studied case.

In physics, a phonon is a collective excitation in a periodic, elastic arrangement of atoms or molecules in condensed matter, specifically in solids and some liquids. Often referred to as a quasiparticle, it is an excited state in the quantum mechanical quantization of the modes of vibrations for elastic structures of interacting particles. Phonons can be thought of as quantized sound waves, similar to photons as quantized light waves.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space $under the operation of composition. By definition, a rotation about the origin is a transformation that preserves the origin, Euclidean distance, and orientation. Every non-trivial rotation is determined by its axis of rotation and its angle of rotation. Composing two rotations results in another rotation; every rotation has a unique inverse rotation; and the identity map satisfies the definition of a rotation. Owing to the above properties, the set of all rotations is a group under composition. Rotations are not commutative, making it a nonabelian group. Moreover, the rotation group has a natural structure as a manifold for which the group operations are smoothly differentiable; so it is in fact a Lie group. It is compact and has dimension 3.$

Assisted Model Building with Energy Refinement (AMBER) is a family of force fields for molecular dynamics of biomolecules originally developed by Peter Kollman's group at the University of California, San Francisco. AMBER is also the name for the molecular dynamics software package that simulates these force fields. It is maintained by an active collaboration between David Case at Rutgers University, Tom Cheatham at the University of Utah, Adrian Roitberg at University of Florida, Ken Merz at Michigan State University, Carlos Simmerling at Stony Brook University, Ray Luo at UC Irvine, and Junmei Wang at Encysive Pharmaceuticals.

In mathematics, a Catalan solid, or Archimedean dual, is a dual polyhedron to an Archimedean solid. There are 13 Catalan solids. They are named for the Belgian mathematician, Eugène Catalan, who first described them in 1865.

The rigid rotor is a mechanical model of rotating systems. An arbitrary rigid rotor is a 3-dimensional rigid object, such as a top. To orient such an object in space requires three angles, known as Euler angles. A special rigid rotor is the linear rotor requiring only two angles to describe, for example of a diatomic molecule. More general molecules are 3-dimensional, such as water, ammonia, or methane. The rigid-rotor Schroedinger equation is discussed in Section 11.2 on pages 240-253 of the textbook by Bunker and Jensen.

In chemistry, orbital hybridisation is the concept of mixing atomic orbitals into new hybrid orbitals suitable for the pairing of electrons to form chemical bonds in valence bond theory. For example, in a carbon atom which forms four single bonds the valence-shell s orbital combines with three valence-shell p orbitals to form four equivalent sp³ mixtures which are arranged in a tetrahedral arrangement around the carbon to bond to four different atoms. Hybrid orbitals are useful in the explanation of molecular geometry and atomic bonding properties and are symmetrically disposed in space. Usually hybrid orbitals are formed by mixing atomic orbitals of comparable energies.

The 3-center 4-electron (3c–4e) bond is a model used to explain bonding in certain hypervalent molecules such as tetratomic and hexatomic interhalogen compounds, sulfur tetrafluoride, the xenon fluorides, and the bifluoride ion. It is also known as the Pimentel–Rundle three-center model after the work published by George C. Pimentel in 1951, which built on concepts developed earlier by Robert E. Rundle for electron-deficient bonding. An extended version of this model is used to describe the whole class of hypervalent molecules such as phosphorus pentafluoride and sulfur hexafluoride as well as multi-center π-bonding such as ozone and sulfur trioxide.

The Hückel method or Hückel molecular orbital theory, proposed by Erich Hückel in 1930, is a very simple linear combination of atomic orbitals molecular orbitals method for the determination of energies of molecular orbitals of π-electrons in π-delocalized molecules, such as ethylene, benzene, butadiene, and pyridine. It is the theoretical basis for Hückel's rule for the aromaticity of $π-electron cyclic, planar systems. It was later extended to conjugated molecules such as pyridine, pyrrole and furan that contain atoms other than carbon, known in this context as heteroatoms. A more dramatic extension of the method to include σ-electrons, known as the extended Hückel method (EHM), was developed by Roald Hoffmann. The extended Hückel method gives some degree of quantitative accuracy for organic molecules in general and was used to provide computational justification for the Woodward-Hoffmann rules. To distinguish the original approach from Hoffmann's extension, the Hückel method is also known as the simple Hückel method (SHM). An elementary description of the application of the simple Huckel method to benzene is given in Sections 3.4.3 and 10.2 of a student level text book.$

In chemistry, Bent's rule describes and explains the relationship between the orbital hybridization of central atoms in molecules and the electronegativities of substituents. The rule was stated by Henry A. Bent as follows:

Atomic s character concentrates in orbitals directed toward electropositive substituents.

The OPLS force field was developed by Prof. William L. Jorgensen at Purdue University and later at Yale University.

The Jaynes–Cummings model is a theoretical model in quantum optics. It describes the system of a two-level atom interacting with a quantized mode of an optical cavity, with or without the presence of light. It was originally developed to study the interaction of atoms with the quantized electromagnetic field in order to investigate the phenomena of spontaneous emission and absorption of photons in a cavity.

A pendulum is a body suspended from a fixed support so that it swings freely back and forth under the influence of gravity. When a pendulum is displaced sideways from its resting, equilibrium position, it is subject to a restoring force due to gravity that will accelerate it back toward the equilibrium position. When released, the restoring force acting on the pendulum's mass causes it to oscillate about the equilibrium position, swinging it back and forth. The mathematics of pendulums are in general quite complicated. Simplifying assumptions can be made, which in the case of a simple pendulum allow the equations of motion to be solved analytically for small-angle oscillations.

Photon polarization is the quantum mechanical description of the classical polarized sinusoidal plane electromagnetic wave. An individual photon can be described as having right or left circular polarization, or a superposition of the two. Equivalently, a photon can be described as having horizontal or vertical linear polarization, or a superposition of the two.

In atomic, molecular, and optical physics and quantum chemistry, the molecular Hamiltonian is the Hamiltonian operator representing the energy of the electrons and nuclei in a molecule. This operator and the associated Schrödinger equation play a central role in computational chemistry and physics for computing properties of molecules and aggregates of molecules, such as thermal conductivity, specific heat, electrical conductivity, optical, and magnetic properties, and reactivity.

The Wigner D-matrix is a unitary matrix in an irreducible representation of the groups SU(2) and SO(3). The complex conjugate of the D-matrix is an eigenfunction of the Hamiltonian of spherical and symmetric rigid rotors. The matrix was introduced in 1927 by Eugene Wigner. $D$ stands for Darstellung, which means "representation" in German.

In organic chemistry, Möbius aromaticity is a special type of aromaticity believed to exist in a number of organic molecules. In terms of molecular orbital theory these compounds have in common a monocyclic array of molecular orbitals in which there is an odd number of out-of-phase overlaps, the opposite pattern compared to the aromatic character to Hückel systems. The nodal plane of the orbitals, viewed as a ribbon, is a Möbius strip, rather than a cylinder, hence the name. The pattern of orbital energies is given by a rotated Frost circle (with the edge of the polygon on the bottom instead of a vertex), so systems with 4n electrons are aromatic, while those with 4n + 2 electrons are anti-aromatic/non-aromatic. Due to incrementally twisted nature of the orbitals of a Möbius aromatic system, stable Möbius aromatic molecules need to contain at least 8 electrons, although 4 electron Möbius aromatic transition states are well known in the context of the Dewar-Zimmerman framework for pericyclic reactions. Möbius molecular systems were considered in 1964 by Edgar Heilbronner by application of the Hückel method, but the first such isolable compound was not synthesized until 2003 by the group of Rainer Herges. However, the fleeting trans-C₉H₉⁺ cation, one conformation of which is shown on the right, was proposed to be a Möbius aromatic reactive intermediate in 1998 based on computational and experimental evidence.

In pure and applied mathematics, quantum mechanics and computer graphics, a tensor operator generalizes the notion of operators which are scalars and vectors. A special class of these are spherical tensor operators which apply the notion of the spherical basis and spherical harmonics. The spherical basis closely relates to the description of angular momentum in quantum mechanics and spherical harmonic functions. The coordinate-free generalization of a tensor operator is known as a representation operator.

In chemistry, isovalent or second order hybridization is an extension of orbital hybridization, the mixing of atomic orbitals into hybrid orbitals which can form chemical bonds, to include fractional numbers of atomic orbitals of each type. It allows for a quantitative depiction of bond formation when the molecular geometry deviates from ideal bond angles.

References

↑ Root, D. M.; Landis, C. R.; Cleveland, T. Valence Bond Concepts Applied to the Molecular Mechanics Description of Molecular Shapes. 1. Application to Nonhypervalent Molecules of the P-Block. J. Am. Chem. Soc.1993, 115, 4201-4209.
↑ Cleveland, T.; Landis, C. R. Valence Bond Concepts Applied to the Molecular Mechanics Description of Molecular Shapes. 2. Application to Hypervalent Molecules of the P-Block. J. Am. Chem. Soc.1996, 118, 6020-6030. doi : 10.1021/ja9506521
↑ Landis, C. R.; Cleveland, T.; Firman; T. K. Valence Bond Concepts Applied to the Molecular Mechanics Description of Molecular Shapes. 3. Application to Transition Metal Alkyls and Hydrides. J. Am. Chem. Soc.1998, 120, 2641-2649. doi : 10.1021/ja9734859
↑ Firman; T. K.; Landis, C. R. Valence Bond Concepts Applied to the Molecular Mechanics Description of Molecular Shapes. 4. Transition Metals with π-Bonds. J. Am. Chem. Soc.2001, 123, 11728-11742. doi : 10.1021/ja002586v
1 2 I. Tubert-Brohman, M. Schmid, and M. Meuwly (2009). "Molecular Mechanics Force Field for Octahedral Organometallic Compounds with Inclusion of the Trans Influence". J. Chem. Theo. Chem. 5: 530–539.CS1 maint: multiple names: authors list (link)
↑ T. Nagy, J. Yosa, and M. Meuwly. "Multisurface Adiabatic Reactive Molecular Dynamics". J. Chem. Theo. Chem. 10: 1366–1375.CS1 maint: multiple names: authors list (link)
↑ M. Schmid, A. K. Das, C. R. Landis and M. Meuwly (2018). "Multi-State VALBOND for Atomistic Simulations of Hypervalent Molecules, Metal Complexes, and Reactions". J. Chem. Theo. Chem. 14: 3565–3578.CS1 maint: multiple names: authors list (link)

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Root, D. M.; Landis, C. R.; Cleveland, T. Valence Bond Concepts Applied to the Molecular Mechanics Description of Molecular Shapes. 1. Application to Nonhypervalent Molecules of the P-Block. J. Am. Chem. Soc.1993, 115, 4201-4209.

[2] Cleveland, T.; Landis, C. R. Valence Bond Concepts Applied to the Molecular Mechanics Description of Molecular Shapes. 2. Application to Hypervalent Molecules of the P-Block. J. Am. Chem. Soc.1996, 118, 6020-6030. doi : 10.1021/ja9506521

[3] Landis, C. R.; Cleveland, T.; Firman; T. K. Valence Bond Concepts Applied to the Molecular Mechanics Description of Molecular Shapes. 3. Application to Transition Metal Alkyls and Hydrides. J. Am. Chem. Soc.1998, 120, 2641-2649. doi : 10.1021/ja9734859

[4] Firman; T. K.; Landis, C. R. Valence Bond Concepts Applied to the Molecular Mechanics Description of Molecular Shapes. 4. Transition Metals with π-Bonds. J. Am. Chem. Soc.2001, 123, 11728-11742. doi : 10.1021/ja002586v

[:0-5] 1 2 I. Tubert-Brohman, M. Schmid, and M. Meuwly (2009). "Molecular Mechanics Force Field for Octahedral Organometallic Compounds with Inclusion of the Trans Influence". J. Chem. Theo. Chem. 5: 530–539.CS1 maint: multiple names: authors list (link)

[6] T. Nagy, J. Yosa, and M. Meuwly. "Multisurface Adiabatic Reactive Molecular Dynamics". J. Chem. Theo. Chem. 10: 1366–1375.CS1 maint: multiple names: authors list (link)

[7] M. Schmid, A. K. Das, C. R. Landis and M. Meuwly (2018). "Multi-State VALBOND for Atomistic Simulations of Hypervalent Molecules, Metal Complexes, and Reactions". J. Chem. Theo. Chem. 14: 3565–3578.CS1 maint: multiple names: authors list (link)

[1]

[2]

[3]

[4]

[5]

[6]

[7]