VALBOND

Last updated July 30, 2024

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 ^{[ broken anchor ]}. 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

Bessel functions, first defined by the mathematician Daniel Bernoulli and then generalized by Friedrich Bessel, are canonical solutions $y (x)$ of Bessel's differential equation $for an arbitrary complex number, which represents the order of the Bessel function. Although and produce the same differential equation, it is conventional to define different Bessel functions for these two values in such a way that the Bessel functions are mostly smooth functions of .$

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.$

Assisted Model Building with Energy Refinement (AMBER) is the name of a widely-used molecular dynamics software package originally developed by Peter Kollman's group at the University of California, San Francisco. It has also, subsequently, come to designate a family of force fields for molecular dynamics of biomolecules that can be used both within the AMBER software suite and with many modern computational platforms.

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

The Morse potential, named after physicist Philip M. Morse, is a convenient interatomic interaction model for the potential energy of a diatomic molecule. It is a better approximation for the vibrational structure of the molecule than the quantum harmonic oscillator because it explicitly includes the effects of bond breaking, such as the existence of unbound states. It also accounts for the anharmonicity of real bonds and the non-zero transition probability for overtone and combination bands. The Morse potential can also be used to model other interactions such as the interaction between an atom and a surface. Due to its simplicity, it is not used in modern spectroscopy. However, its mathematical form inspired the MLR (Morse/Long-range) potential, which is the most popular potential energy function used for fitting spectroscopic data.

In chemistry, a hypervalent molecule is a molecule that contains one or more main group elements apparently bearing more than eight electrons in their valence shells. Phosphorus pentachloride, sulfur hexafluoride, chlorine trifluoride, the chlorite ion in chlorous acid and the triiodide ion are examples of hypervalent molecules.

In rotordynamics, 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.

In chemistry, orbital hybridisation is the concept of mixing atomic orbitals to form 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 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.

Wetting is the ability of a liquid to displace gas to maintain contact with a solid surface, resulting from intermolecular interactions when the two are brought together. This happens in presence of a gaseous phase or another liquid phase not miscible with the first one. The degree of wetting (wettability) is determined by a force balance between adhesive and cohesive forces. There are two types of wetting: non-reactive wetting and reactive wetting.

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 simple method for calculating molecular orbitals as linear combinations of atomic orbitals. The theory predicts the molecular orbitals for π-electrons in π-delocalized molecules, such as ethylene, benzene, butadiene, and pyridine. It provides the theoretical basis for Hückel's rule that cyclic, planar molecules or ions with $π-electrons are aromatic. It was later extended to conjugated molecules such as pyridine, pyrrole and furan that contain atoms other than carbon and hydrogen (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).$

A theoretical motivation for general relativity, including the motivation for the geodesic equation and the Einstein field equation, can be obtained from special relativity by examining the dynamics of particles in circular orbits about the Earth. A key advantage in examining circular orbits is that it is possible to know the solution of the Einstein Field Equation a priori. This provides a means to inform and verify the formalism.

In chemistry, Bent's rule describes and explains the relationship between the orbital hybridization 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, and is being further developed commercially by Schrödinger, Inc.

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.

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.

A hydrogen-like atom (or hydrogenic atom) is any atom or ion with a single valence electron. These atoms are isoelectronic with hydrogen. Examples of hydrogen-like atoms include, but are not limited to, hydrogen itself, all alkali metals such as Rb and Cs, singly ionized alkaline earth metals such as Ca⁺ and Sr⁺ and other ions such as He⁺, Li²⁺, and Be³⁺ and isotopes of any of the above. A hydrogen-like atom includes a positively charged core consisting of the atomic nucleus and any core electrons as well as a single valence electron. Because helium is common in the universe, the spectroscopy of singly ionized helium is important in EUV astronomy, for example, of DO white dwarf stars.

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.

In computational chemistry, natural resonance theory (NRT) is an iterative, variational functional embedded into the natural bond orbital (NBO) program, commonly run in Gaussian, GAMESS, ORCA, Ampac and other software packages. NRT was developed in 1997 by Frank A. Weinhold and Eric D. Glendening, chemistry professors at University of Wisconsin-Madison and Indiana State University, respectively. Given a list of NBOs for an idealized natural Lewis structure, the NRT functional creates a list of Lewis resonance structures and calculates the resonance weights of each contributing resonance structure. Structural and chemical properties, such as bond order, valency, and bond polarity, may be calculated from resonance weights. Specifically, bond orders may be divided into their covalent and ionic contributions, while valency is the sum of bond orders of a given atom. This aims to provide quantitative results that agree with qualitative notions of chemical resonance. In contrast to the "wavefunction resonance theory" (i.e., the superposition of wavefunctions), NRT uses the density matrix resonance theory, performing a superposition of density matrices to realize resonance. NRT has applications in ab initio calculations, including calculating the bond orders of intra- and intermolecular interactions and the resonance weights of radical isomers.

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 (3): 530–539. doi:10.1021/ct800392n. PMID 26610220.{{cite journal}}: CS1 maint: multiple names: authors list (link)
↑ T. Nagy, J. Yosa, and M. Meuwly (2014). "Multisurface Adiabatic Reactive Molecular Dynamics". J. Chem. Theo. Chem. 10 (4): 1366–1375. doi:10.1021/ct400953f. PMID 26580356.{{cite journal}}: 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 (7): 3565–3578. doi:10.1021/acs.jctc.7b01210. PMID 29490133. S2CID 3632002.{{cite journal}}: 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 (3): 530–539. doi:10.1021/ct800392n. PMID 26610220.{{cite journal}}: CS1 maint: multiple names: authors list (link)

[6] T. Nagy, J. Yosa, and M. Meuwly (2014). "Multisurface Adiabatic Reactive Molecular Dynamics". J. Chem. Theo. Chem. 10 (4): 1366–1375. doi:10.1021/ct400953f. PMID 26580356.{{cite journal}}: 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 (7): 3565–3578. doi:10.1021/acs.jctc.7b01210. PMID 29490133. S2CID 3632002.{{cite journal}}: CS1 maint: multiple names: authors list (link)

[1]

[2]

[3]

[4]

[5]

[6]

[7]