Morse/Long-range potential

Last updated

The Morse/Long-range potential (MLR potential) is an interatomic interaction model for the potential energy of a diatomic molecule. Due to the simplicity of the regular Morse potential (it only has three adjustable parameters), it is very limited in its applicability in modern spectroscopy. The MLR potential is a modern version of the Morse potential which has the correct theoretical long-range form of the potential naturally built into it. [1] It has been an important tool for spectroscopists to represent experimental data, verify measurements, and make predictions. It is useful for its extrapolation capability when data for certain regions of the potential are missing, its ability to predict energies with accuracy often better than the most sophisticated ab initio techniques, and its ability to determine precise empirical values for physical parameters such as the dissociation energy, equilibrium bond length, and long-range constants. Cases of particular note include:

Contents

  1. the c-state of dilithium (Li2): where the MLR potential was successfully able to bridge a gap of more than 5000 cm−1 in experimental data. [2] Two years later it was found that the MLR potential was able to successfully predict the energies in the middle of this gap, correctly within about 1 cm−1. [3] The accuracy of these predictions was much better than the most sophisticated ab initio techniques at the time. [4]
  2. the A-state of Li2: where Le Roy et al. [1] constructed an MLR potential which determined the C3 value for atomic lithium to a higher-precision than any previously measured atomic oscillator strength, by an order of magnitude. [5] This lithium oscillator strength is related to the radiative lifetime of atomic lithium and is used as a benchmark for atomic clocks and measurements of fundamental constants.
  3. the a-state of KLi: where the MLR was used to build an analytic global potential successfully despite there only being a small amount of levels observed near the top of the potential. [6]

Historical origins

The MLR potential is based on the classic Morse potential which was first introduced in 1929 by Philip M. Morse. A primitive version of the MLR potential was first introduced in 2006 by Robert J. Le Roy and colleagues for a study on N2. [7] This primitive form was used on Ca2, [8] KLi [6] and MgH, [9] [10] [11] before the more modern version was introduced in 2009. [1] A further extension of the MLR potential referred to as the MLR3 potential was introduced in a 2010 study of Cs2, [12] and this potential has since been used on HF, [13] [14] HCl, [13] [14] HBr [13] [14] and HI. [13] [14]

Function

The Morse/Long-range potential energy function is of the form

where for large ,

so is defined according to the theoretically correct long-range behavior expected for the interatomic interaction. is the depth of the potential at equilibrium.

This long-range form of the MLR model is guaranteed because the argument of the exponent is defined to have long-range behavior:

where is the equilibrium bond length.

There are a few ways in which this long-range behavior can be achieved, the most common is to make a polynomial that is constrained to become at long-range:

where n is an integer greater than 1, which value is defined by the model chosen for the long-range potential .

It is clear to see that:

Applications

The MLR potential has successfully summarized all experimental spectroscopic data (and/or virial data) for a number of diatomic molecules, including: N2, [7] Ca2, [8] KLi, [6] MgH, [9] [10] [11] several electronic states of Li2, [1] [2] [15] [3] [10] Cs2, [16] [12] Sr2, [17] ArXe, [10] [18] LiCa, [19] LiNa, [20] Br2, [21] Mg2, [22] HF, [13] [14] HCl, [13] [14] HBr, [13] [14] HI, [13] [14] MgD, [9] Be2, [23] BeH, [24] and NaH. [25] More sophisticated versions are used for polyatomic molecules.

It has also become customary to fit ab initio points to the MLR potential, to achieve a fully analytic ab initio potential and to take advantage of the MLR's ability to incorporate the correct theoretically known short- and long-range behavior into the potential (the latter usually being of higher accuracy than the molecular ab initio points themselves because it is based on atomic ab initio calculations rather than molecular ones, and because features like spin-orbit coupling which are difficult to incorporate into molecular ab initio calculations can more easily be treated in the long-range). MLR has been used to represent ab initio points for KLi [26] and KBe. [27]

See also

Related Research Articles

Density-functional theory (DFT) is a computational quantum mechanical modelling method used in physics, chemistry and materials science to investigate the electronic structure of many-body systems, in particular atoms, molecules, and the condensed phases. Using this theory, the properties of a many-electron system can be determined by using functionals, i.e. functions of another function. In the case of DFT, these are functionals of the spatially dependent electron density. DFT is among the most popular and versatile methods available in condensed-matter physics, computational physics, and computational chemistry.

Dilithium, Li2, is a strongly electrophilic, diatomic molecule comprising two lithium atoms covalently bonded together. Li2 is known in the gas phase. It has a bond order of 1, an internuclear separation of 267.3 pm and a bond energy of 102 kJ/mol or 1.06 eV in each bond. The electron configuration of Li2 may be written as σ2.

<span class="mw-page-title-main">Morse potential</span> Model for the potential energy of a diatomic molecule

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.

Koopmans' theorem states that in closed-shell Hartree–Fock theory (HF), the first ionization energy of a molecular system is equal to the negative of the orbital energy of the highest occupied molecular orbital (HOMO). This theorem is named after Tjalling Koopmans, who published this result in 1934.

The isothermal–isobaric ensemble is a statistical mechanical ensemble that maintains constant temperature and constant pressure applied. It is also called the -ensemble, where the number of particles is also kept as a constant. This ensemble plays an important role in chemistry as chemical reactions are usually carried out under constant pressure condition. The NPT ensemble is also useful for measuring the equation of state of model systems whose virial expansion for pressure cannot be evaluated, or systems near first-order phase transitions.

<span class="mw-page-title-main">Barium fluoride</span> Chemical compound

Barium fluoride is an inorganic compound with the formula BaF2. It is a colorless solid that occurs in nature as the rare mineral frankdicksonite. Under standard conditions it adopts the fluorite structure and at high pressure the PbCl2 structure. Like CaF2, it is resilient to and insoluble in water.

<span class="mw-page-title-main">Spartan (chemistry software)</span>

Spartan is a molecular modelling and computational chemistry application from Wavefunction. It contains code for molecular mechanics, semi-empirical methods, ab initio models, density functional models, post-Hartree–Fock models, and thermochemical recipes including G3(MP2) and T1. Quantum chemistry calculations in Spartan are powered by Q-Chem.

<span class="mw-page-title-main">Stretched exponential function</span>

The stretched exponential function

Car–Parrinello molecular dynamics or CPMD refers to either a method used in molecular dynamics or the computational chemistry software package used to implement this method.

This page deals with the electron affinity as a property of isolated atoms or molecules. Solid state electron affinities are not listed here.

<span class="mw-page-title-main">Helium hydride ion</span> Chemical compound

The helium hydride ion, hydridohelium(1+) ion, or helonium is a cation (positively charged ion) with chemical formula HeH+. It consists of a helium atom bonded to a hydrogen atom, with one electron removed. It can also be viewed as protonated helium. It is the lightest heteronuclear ion, and is believed to be the first compound formed in the Universe after the Big Bang.

<span class="mw-page-title-main">Water cluster</span> Cluster of water molecules held together through hydrogen bonding

In chemistry, a water cluster is a discrete hydrogen bonded assembly or cluster of molecules of water. Many such clusters have been predicted by theoretical models (in silico), and some have been detected experimentally in various contexts such as ice, bulk liquid water, in the gas phase, in dilute mixtures with non-polar solvents, and as water of hydration in crystal lattices. The simplest example is the water dimer (H2O)2.

<span class="mw-page-title-main">Two-dimensional infrared spectroscopy</span> Nonlinear infrared spectroscopy technique

Two-dimensional infrared spectroscopy is a nonlinear infrared spectroscopy technique that has the ability to correlate vibrational modes in condensed-phase systems. This technique provides information beyond linear infrared spectra, by spreading the vibrational information along multiple axes, yielding a frequency correlation spectrum. A frequency correlation spectrum can offer structural information such as vibrational mode coupling, anharmonicities, along with chemical dynamics such as energy transfer rates and molecular dynamics with femtosecond time resolution. 2DIR experiments have only become possible with the development of ultrafast lasers and the ability to generate femtosecond infrared pulses.

In computational chemistry, spin contamination is the artificial mixing of different electronic spin-states. This can occur when an approximate orbital-based wave function is represented in an unrestricted form – that is, when the spatial parts of α and β spin-orbitals are permitted to differ. Approximate wave functions with a high degree of spin contamination are undesirable. In particular, they are not eigenfunctions of the total spin-squared operator, Ŝ2, but can formally be expanded in terms of pure spin states of higher multiplicities.

CONQUEST is a linear scaling, or O(N), density functional theory (DFT) electronic structure open-source code. The code is designed to perform DFT calculations on very large systems containing many thousands of atoms. It can be run at different levels of precision ranging from ab initio tight binding up to full DFT with plane wave accuracy. It has been applied to the study of three-dimensional reconstructions formed by Ge on Si(001), containing over 20,000 atoms. Tests on the UK's national supercomputer HECToR in 2009 demonstrated the capability of the code to perform ground-state calculations on systems of over 1,000,000 atoms.

Heat transfer physics describes the kinetics of energy storage, transport, and energy transformation by principal energy carriers: phonons, electrons, fluid particles, and photons. Heat is energy stored in temperature-dependent motion of particles including electrons, atomic nuclei, individual atoms, and molecules. Heat is transferred to and from matter by the principal energy carriers. The state of energy stored within matter, or transported by the carriers, is described by a combination of classical and quantum statistical mechanics. The energy is different made (converted) among various carriers. The heat transfer processes are governed by the rates at which various related physical phenomena occur, such as the rate of particle collisions in classical mechanics. These various states and kinetics determine the heat transfer, i.e., the net rate of energy storage or transport. Governing these process from the atomic level to macroscale are the laws of thermodynamics, including conservation of energy.

In computational chemistry, a solvent model is a computational method that accounts for the behavior of solvated condensed phases. Solvent models enable simulations and thermodynamic calculations applicable to reactions and processes which take place in solution. These include biological, chemical and environmental processes. Such calculations can lead to new predictions about the physical processes occurring by improved understanding.

<span class="mw-page-title-main">Dual graviton</span> Hypothetical particle found in supergravity

In theoretical physics, the dual graviton is a hypothetical elementary particle that is a dual of the graviton under electric-magnetic duality, as an S-duality, predicted by some formulations of supergravity in eleven dimensions.

<span class="mw-page-title-main">Diargon</span> Chemical compound

Diargon or the argon dimer is a molecule containing two argon atoms. Normally, this is only very weakly bound together by van der Waals forces. However, in an excited state, or ionised state, the two atoms can be more tightly bound together, with significant spectral features. At cryogenic temperatures, argon gas can have a few percent of diargon molecules.

Adiabatic electron-transfer is a type of oxidation-reduction processes. The mechanism is ubiquitous in nature in both the inorganic and biological spheres. Adiabatic electron-transfers proceed without making or breaking chemical bonds. Adiabatic electron-transfer can occur by either optical or thermal mechanisms. Electron transfer during a collision between an oxidant and a reductant occurs adiabatically on a continuous potential-energy surface.

References

  1. 1 2 3 4 Le Roy, Robert J.; N. S. Dattani; J. A. Coxon; A. J. Ross; Patrick Crozet; C. Linton (2009). "Accurate analytic potentials for Li2(X) and Li2(A) from 2 to 90 Angstroms, and the radiative lifetime of Li(2p)". Journal of Chemical Physics. 131 (20): 204309. Bibcode:2009JChPh.131t4309L. doi:10.1063/1.3264688. PMID   19947682.
  2. 1 2 Dattani, N. S.; R. J. Le Roy (8 May 2011). "A DPF data analysis yields accurate analytic potentials for Li2(a) and Li2(c) that incorporate 3-state mixing near the c-state asymptote". Journal of Molecular Spectroscopy. 268 (1–2): 199–210. arXiv: 1101.1361 . Bibcode:2011JMoSp.268..199D. doi:10.1016/j.jms.2011.03.030. S2CID   119266866.
  3. 1 2 Semczuk, M.; Li, X.; Gunton, W.; Haw, M.; Dattani, N. S.; Witz, J.; Mills, A. K.; Jones, D. J.; Madison, K. W. (2013). "High-resolution photoassociation spectroscopy of the 6Li2 13Σ+ state". Phys. Rev. A. 87 (5): 052505. arXiv: 1309.6662 . Bibcode:2013PhRvA..87e2505S. doi:10.1103/PhysRevA.87.052505. S2CID   119263860.
  4. ↑ Halls, M. S.; H. B. Schlegal; M. J. DeWitt; G. F. W. Drake (18 May 2001). "Ab initio calculation of the a-state interaction potential and vibrational levels of 7Li2" (PDF). Chemical Physics Letters. 339 (5–6): 427–432. Bibcode:2001CPL...339..427H. doi:10.1016/s0009-2614(01)00403-1.
  5. ↑ L-Y. Tang; Z-C. Yan; T-Y. Shi; J. Mitroy (30 November 2011). "Third-order perturbation theory for van der Waals interaction coefficients". Physical Review A. 84 (5): 052502. Bibcode:2011PhRvA..84e2502T. doi:10.1103/PhysRevA.84.052502.
  6. 1 2 3 Salami, H.; A. J. Ross; P. Crozet; W. Jastrzebski; P. Kowalczyk; R. J. Le Roy (2007). "A full analytic potential energy curve for the a3Σ+ state of KLi from a limited vibrational data set". Journal of Chemical Physics. 126 (19): 194313. Bibcode:2007JChPh.126s4313S. doi: 10.1063/1.2734973 . PMID   17523810.
  7. 1 2 Le Roy, R. J.; Y. Huang; C. Jary (2006). "An accurate analytic potential function for ground-state N2 from a direct-potential-fit analysis of spectroscopic data". Journal of Chemical Physics. 125 (16): 164310. Bibcode:2006JChPh.125p4310L. doi:10.1063/1.2354502. PMID   17092076.
  8. 1 2 Le Roy, Robert J.; R. D. E. Henderson (2007). "A new potential function form incorporating extended long-range behaviour: application to ground-state Ca2". Molecular Physics. 105 (5–7): 663–677. Bibcode:2007MolPh.105..663L. doi:10.1080/00268970701241656. S2CID   94174485.
  9. 1 2 3 Henderson, R. D. E.; A. Shayesteh; J. Tao; C. Haugen; P. F. Bernath; R. J. Le Roy (4 October 2013). "Accurate Analytic Potential and Born–Oppenheimer Breakdown Functions for MgH and MgD from a Direct-Potential-Fit Data Analysis". The Journal of Physical Chemistry A. 117 (50): 13373–87. Bibcode:2013JPCA..11713373H. doi:10.1021/jp406680r. PMID   24093511.
  10. 1 2 3 4 Le Roy, R. J.; C. C. Haugen; J. Tao; H. Li (February 2011). "Long-range damping functions improve the short-range behaviour of 'MLR' potential energy functions" (PDF). Molecular Physics. 109 (3): 435–446. Bibcode:2011MolPh.109..435L. doi:10.1080/00268976.2010.527304. S2CID   97119318.
  11. 1 2 Shayesteh, A.; R. D. E. Henderson; R. J. Le Roy; P. F. Bernath (2007). "Ground State Potential Energy Curve and Dissociation Energy of MgH". The Journal of Physical Chemistry A. 111 (49): 12495–12505. Bibcode:2007JPCA..11112495S. CiteSeerX   10.1.1.584.8808 . doi:10.1021/jp075704a. PMID   18020428.
  12. 1 2 Coxon, J. A.; P. G. Hajigeorgiou (2010). "The ground X 1Σ+g electronic state of the cesium dimer: Application of a direct potential fitting procedure". Journal of Chemical Physics. 132 (9): 094105. Bibcode:2010JChPh.132i4105C. doi:10.1063/1.3319739. PMID   20210387.
  13. 1 2 3 4 5 6 7 8 Li, Gang; I. E. Gordon; P. G. Hajigeorgiou; J. A. Coxon; L. S. Rothman (2013). "Reference spectroscopic data for hydrogen halides, Part II: The line lists". Journal of Quantitative Spectroscopy & Radiative Transfer. 130: 284–295. Bibcode:2013JQSRT.130..284L. doi:10.1016/j.jqsrt.2013.07.019.
  14. 1 2 3 4 5 6 7 8 Coxon, John A.; Hajigeorgiou, Photos G. (2015). "Improved direct potential fit analyses for the ground electronic states of the hydrogen halides: HF/DF/TF, HCl/DCl/TCl, HBr/DBr/TBr and HI/DI/TI". Journal of Quantitative Spectroscopy and Radiative Transfer. 151: 133–154. Bibcode:2015JQSRT.151..133C. doi:10.1016/j.jqsrt.2014.08.028.
  15. ↑ Gunton, Will; Semczuk, Mariusz; Dattani, Nikesh S.; Madison, Kirk W. (2013). "High resolution photoassociation spectroscopy of the 6Li2 A(11Σu+) state". Physical Review A. 88 (6): 062510. arXiv: 1309.5870 . Bibcode:2013PhRvA..88f2510G. doi:10.1103/PhysRevA.88.062510. S2CID   119268157.
  16. ↑ Xie, F.; L. Li; D. Li; V. B. Sovkov; K. V. Minaev; V. S. Ivanov; A. M. Lyyra; S. Magnier (2011). "Joint analysis of the Cs2 a-state and 1g(33Π11g) states". Journal of Chemical Physics. 135 (2): 02403. Bibcode:2011JChPh.135b4303X. doi:10.1063/1.3606397. PMID   21766938.
  17. ↑ Stein, A.; H. Knockel; E. Tiemann (April 2010). "The 1S+1S asymptote of Sr2 studied by Fourier-transform spectroscopy". The European Physical Journal D. 57 (2): 171–177. arXiv: 1001.2741 . Bibcode:2010EPJD...57..171S. doi:10.1140/epjd/e2010-00058-y. S2CID   119243162.
  18. ↑ Piticco, Lorena; F. Merkt; A. A. Cholewinski; F. R. W. McCourt; R. J. Le Roy (December 2010). "Rovibrational structure and potential energy function of the ground electronic state of ArXe". Journal of Molecular Spectroscopy. 264 (2): 83–93. Bibcode:2010JMoSp.264...83P. doi:10.1016/j.jms.2010.08.007. hdl: 20.500.11850/210096 .
  19. ↑ Ivanova, Milena; A. Stein; A. Pashov; A. V. Stolyarov; H. Knockel; E. Tiemann (2011). "The X2Σ+ state of LiCa studied by Fourier-transform spectroscopy". Journal of Chemical Physics. 135 (17): 174303. Bibcode:2011JChPh.135q4303I. doi:10.1063/1.3652755. PMID   22070298.
  20. ↑ Steinke, M.; H. Knockel; E. Tiemann (27 April 2012). "X-state of LiNa studied by Fourier-transform spectroscopy". Physical Review A. 85 (4): 042720. Bibcode:2012PhRvA..85d2720S. doi:10.1103/PhysRevA.85.042720.
  21. ↑ Yukiya, T.; N. Nishimiya; Y. Samejima; K. Yamaguchi; M. Suzuki; C. D. Boonec; I. Ozier; R. J. Le Roy (January 2013). "Direct-potential-fit analysis for the system of Br2". Journal of Molecular Spectroscopy. 283: 32–43. Bibcode:2013JMoSp.283...32Y. doi:10.1016/j.jms.2012.12.006.
  22. ↑ Knockel, H.; S. Ruhmann; E. Tiemann (2013). "The X-state of Mg2 studied by Fourier-transform spectroscopy". Journal of Chemical Physics. 138 (9): 094303. Bibcode:2013JChPh.138i4303K. doi:10.1063/1.4792725. PMID   23485290.
  23. ↑ Meshkov, Vladimir V.; Stolyarov, Andrey V.; Heaven, Michael C.; Haugen, Carl; Leroy, Robert J. (2014). "Direct-potential-fit analyses yield improved empirical potentials for the ground X1Σg+ state of Be2". The Journal of Chemical Physics. 140 (6): 064315. Bibcode:2014JChPh.140f4315M. doi:10.1063/1.4864355. PMID   24527923.
  24. ↑ Dattani, Nikesh S. (2015). "Beryllium monohydride (BeH): Where we are now, after 86 years of spectroscopy". Journal of Molecular Spectroscopy. 311: 76–83. arXiv: 1408.3301 . Bibcode:2015JMoSp.311...76D. doi:10.1016/j.jms.2014.09.005. S2CID   118542048.
  25. ↑ Walji, Sadru-Dean; Sentjens, Katherine M.; Le Roy, Robert J. (2015). "Dissociation energies and potential energy functions for the ground X 1Σ+ and "avoided-crossing" A 1Σ+ states of NaH". The Journal of Chemical Physics. 142 (4): 044305. Bibcode:2015JChPh.142d4305W. doi:10.1063/1.4906086. PMID   25637985.
  26. ↑ Xiao, Ke-La; Yang, Chuan-Lu; Wang, Mei-Shan; Ma, Xiao-Guang; Liu, Wen-Wang (2013). "The effect of inner-shell electrons on the ground and low-lying excited states of KLi: Ab initio study with all-electron basis sets". Journal of Quantitative Spectroscopy and Radiative Transfer. 129: 8–14. Bibcode:2013JQSRT.129....8X. doi:10.1016/j.jqsrt.2013.05.025.
  27. ↑ Xiao, Ke-La; Yang, Chuan-Lu; Wang, Mei-Shan; Ma, Xiao-Guang; Liu, Wen-Wang (2013). "An ab initio study of the ground and low-lying excited states of KBe with the effect of inner-shell electrons". The Journal of Chemical Physics. 139 (7): 074305. Bibcode:2013JChPh.139g4305X. doi:10.1063/1.4818452. PMID   23968090.
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.