Morse potential

Last updated

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 (only three fitting parameters), 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.

Contents

Potential energy function

The Morse potential (blue) and harmonic oscillator potential (green). Unlike the energy levels of the harmonic oscillator potential, which are evenly spaced by ho, the Morse potential level spacing decreases as the energy approaches the dissociation energy. The dissociation energy De is larger than the true energy required for dissociation D0 due to the zero point energy of the lowest (v = 0) vibrational level. Morse-potential.png
The Morse potential (blue) and harmonic oscillator potential (green). Unlike the energy levels of the harmonic oscillator potential, which are evenly spaced by ħω, the Morse potential level spacing decreases as the energy approaches the dissociation energy. The dissociation energy De is larger than the true energy required for dissociation D0 due to the zero point energy of the lowest (v = 0) vibrational level.

The Morse potential energy function is of the form

Here is the distance between the atoms, is the equilibrium bond distance, is the well depth (defined relative to the dissociated atoms), and controls the 'width' of the potential (the smaller is, the larger the well). The dissociation energy of the bond can be calculated by subtracting the zero point energy from the depth of the well. The force constant (stiffness) of the bond can be found by Taylor expansion of around to the second derivative of the potential energy function, from which it can be shown that the parameter, , is

where is the force constant at the minimum of the well.

Since the zero of potential energy is arbitrary, the equation for the Morse potential can be rewritten any number of ways by adding or subtracting a constant value. When it is used to model the atom-surface interaction, the energy zero can be redefined so that the Morse potential becomes

which is usually written as

where is now the coordinate perpendicular to the surface. This form approaches zero at infinite and equals at its minimum, i.e. . It clearly shows that the Morse potential is the combination of a short-range repulsion term (the former) and a long-range attractive term (the latter), analogous to the Lennard-Jones potential.

Vibrational states and energies

Like the quantum harmonic oscillator, the energies and eigenstates of the Morse potential can be found using operator methods. [1] One approach involves applying the factorization method to the Hamiltonian.

To write the stationary states on the Morse potential, i.e. solutions and of the following Schrödinger equation:

it is convenient to introduce the new variables:

Then, the Schrödinger equation takes the simple form:

Its eigenvalues (reduced by ) and eigenstates can be written as: [2]

where

with denoting the largest integer smaller than , and

where (which satisfies the normalization condition ) and is a generalized Laguerre polynomial:

There also exists the following analytical expression for matrix elements of the coordinate operator: [3]

which is valid for and . The eigenenergies in the initial variables have the form:

where is the vibrational quantum number and has units of frequency. The latter is mathematically related to the particle mass, , and the Morse constants via

Whereas the energy spacing between vibrational levels in the quantum harmonic oscillator is constant at , the energy between adjacent levels decreases with increasing in the Morse oscillator. Mathematically, the spacing of Morse levels is

This trend matches the anharmonicity found in real molecules. However, this equation fails above some value of where is calculated to be zero or negative. Specifically,

integer part.

This failure is due to the finite number of bound levels in the Morse potential, and some maximum that remains bound. For energies above , all the possible energy levels are allowed and the equation for is no longer valid.

Below , is a good approximation for the true vibrational structure in non-rotating diatomic molecules. In fact, the real molecular spectra are generally fit to the form1

in which the constants and can be directly related to the parameters for the Morse potential. Specifically,

and

Note that if and are given in cm, is in cm/s (not m/s), is in kg, and is in J·s, then will be in m and will be in cm.

As is clear from dimensional analysis, for historical reasons the last equation uses spectroscopic notation in which represents a wavenumber obeying , and not an angular frequency given by .

Harmonic oscillator (grey) and Morse (black) potentials curves are shown along with their eigenfunctions (respectively green and blue for harmonic oscillator and morse) for the same vibrational levels for nitrogen. N2ground.png
Harmonic oscillator (grey) and Morse (black) potentials curves are shown along with their eigenfunctions (respectively green and blue for harmonic oscillator and morse) for the same vibrational levels for nitrogen.

Morse/Long-range potential

An extension of the Morse potential that made the Morse form useful for modern (high-resolution) spectroscopy is the MLR (Morse/Long-range) potential. [4] The MLR potential is used as a standard for representing spectroscopic and/or virial data of diatomic molecules by a potential energy curve. It has been used on N2, [5] Ca2, [6] KLi, [7] MgH, [8] [9] [10] several electronic states of Li2, [4] [11] [12] [13] [9] Cs2, [14] [15] Sr2, [16] ArXe, [9] [17] LiCa, [18] LiNa, [19] Br2, [20] Mg2, [21] HF, [22] [23] HCl, [22] [23] HBr, [22] [23] HI, [22] [23] MgD, [8] Be2, [24] BeH, [25] and NaH. [26] More sophisticated versions are used for polyatomic molecules.

See also

Related Research Articles

In particle physics, the Dirac equation is a relativistic wave equation derived by British physicist Paul Dirac in 1928. In its free form, or including electromagnetic interactions, it describes all spin-1/2 massive particles, called "Dirac particles", such as electrons and quarks for which parity is a symmetry. It is consistent with both the principles of quantum mechanics and the theory of special relativity, and was the first theory to account fully for special relativity in the context of quantum mechanics. It was validated by accounting for the fine structure of the hydrogen spectrum in a completely rigorous way. It has become vital in the building of the Standard Model.

<span class="mw-page-title-main">Schrödinger equation</span> Description of a quantum-mechanical system

The Schrödinger equation is a partial differential equation that governs the wave function of a quantum-mechanical system. Its discovery was a significant landmark in the development of quantum mechanics. It is named after Erwin Schrödinger, who postulated the equation in 1925 and published it in 1926, forming the basis for the work that resulted in his Nobel Prize in Physics in 1933.

<span class="mw-page-title-main">Magnetic circular dichroism</span>

Magnetic circular dichroism (MCD) is the differential absorption of left and right circularly polarized light, induced in a sample by a strong magnetic field oriented parallel to the direction of light propagation. MCD measurements can detect transitions which are too weak to be seen in conventional optical absorption spectra, and it can be used to distinguish between overlapping transitions. Paramagnetic systems are common analytes, as their near-degenerate magnetic sublevels provide strong MCD intensity that varies with both field strength and sample temperature. The MCD signal also provides insight into the symmetry of the electronic levels of the studied systems, such as metal ion sites.

A bound state is a composite of two or more fundamental building blocks, such as particles, atoms, or bodies, that behaves as a single object and in which energy is required to split them.

<span class="mw-page-title-main">Nonlinear Schrödinger equation</span> Nonlinear form of the Schrödinger equation

In theoretical physics, the (one-dimensional) nonlinear Schrödinger equation (NLSE) is a nonlinear variation of the Schrödinger equation. It is a classical field equation whose principal applications are to the propagation of light in nonlinear optical fibers and planar waveguides and to Bose–Einstein condensates confined to highly anisotropic, cigar-shaped traps, in the mean-field regime. Additionally, the equation appears in the studies of small-amplitude gravity waves on the surface of deep inviscid (zero-viscosity) water; the Langmuir waves in hot plasmas; the propagation of plane-diffracted wave beams in the focusing regions of the ionosphere; the propagation of Davydov's alpha-helix solitons, which are responsible for energy transport along molecular chains; and many others. More generally, the NLSE appears as one of universal equations that describe the evolution of slowly varying packets of quasi-monochromatic waves in weakly nonlinear media that have dispersion. Unlike the linear Schrödinger equation, the NLSE never describes the time evolution of a quantum state. The 1D NLSE is an example of an integrable model.

The Compton wavelength is a quantum mechanical property of a particle, defined as the wavelength of a photon whose energy is the same as the rest energy of that particle. It was introduced by Arthur Compton in 1923 in his explanation of the scattering of photons by electrons.

In probability theory and mathematical physics, a random matrix is a matrix-valued random variable—that is, a matrix in which some or all of its entries are sampled randomly from a probability distribution. Random matrix theory (RMT) is the study of properties of random matrices, often as they become large. RMT provides techniques like mean-field theory, diagrammatic methods, the cavity method, or the replica method to compute quantities like traces, spectral densities, or scalar products between eigenvectors. Many physical phenomena, such as the spectrum of nuclei of heavy atoms, the thermal conductivity of a lattice, or the emergence of quantum chaos, can be modeled mathematically as problems concerning large, random matrices.

In quantum mechanics, the Hellmann–Feynman theorem relates the derivative of the total energy with respect to a parameter to the expectation value of the derivative of the Hamiltonian with respect to that same parameter. According to the theorem, once the spatial distribution of the electrons has been determined by solving the Schrödinger equation, all the forces in the system can be calculated using classical electrostatics.

<span class="mw-page-title-main">Ponderomotive force</span> Nonlinear force experienced by a charged particle

In physics, a ponderomotive force is a nonlinear force that a charged particle experiences in an inhomogeneous oscillating electromagnetic field. It causes the particle to move towards the area of the weaker field strength, rather than oscillating around an initial point as happens in a homogeneous field. This occurs because the particle sees a greater magnitude of force during the half of the oscillation period while it is in the area with the stronger field. The net force during its period in the weaker area in the second half of the oscillation does not offset the net force of the first half, and so over a complete cycle this makes the particle move towards the area of lesser force.

<span class="mw-page-title-main">Inelastic mean free path</span> Index of how far electrons travel through a solid before losing energy

The inelastic mean free path (IMFP) is an index of how far an electron on average travels through a solid before losing energy.

<span class="mw-page-title-main">Jaynes–Cummings model</span> Model in quantum optics

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.

In quantum mechanics the delta potential is a potential well mathematically described by the Dirac delta function - a generalized function. Qualitatively, it corresponds to a potential which is zero everywhere, except at a single point, where it takes an infinite value. This can be used to simulate situations where a particle is free to move in two regions of space with a barrier between the two regions. For example, an electron can move almost freely in a conducting material, but if two conducting surfaces are put close together, the interface between them acts as a barrier for the electron that can be approximated by a delta potential.

The Gross–Pitaevskii equation describes the ground state of a quantum system of identical bosons using the Hartree–Fock approximation and the pseudopotential interaction model.

In mathematics, the spectral theory of ordinary differential equations is the part of spectral theory concerned with the determination of the spectrum and eigenfunction expansion associated with a linear ordinary differential equation. In his dissertation, Hermann Weyl generalized the classical Sturm–Liouville theory on a finite closed interval to second order differential operators with singularities at the endpoints of the interval, possibly semi-infinite or infinite. Unlike the classical case, the spectrum may no longer consist of just a countable set of eigenvalues, but may also contain a continuous part. In this case the eigenfunction expansion involves an integral over the continuous part with respect to a spectral measure, given by the Titchmarsh–Kodaira formula. The theory was put in its final simplified form for singular differential equations of even degree by Kodaira and others, using von Neumann's spectral theorem. It has had important applications in quantum mechanics, operator theory and harmonic analysis on semisimple Lie groups.

The Kapitza–Dirac effect is a quantum mechanical effect consisting of the diffraction of matter by a standing wave of light. The effect was first predicted as the diffraction of electrons from a standing wave of light by Paul Dirac and Pyotr Kapitsa in 1933. The effect relies on the wave–particle duality of matter as stated by the de Broglie hypothesis in 1924.

In linear algebra, a raising or lowering operator is an operator that increases or decreases the eigenvalue of another operator. In quantum mechanics, the raising operator is sometimes called the creation operator, and the lowering operator the annihilation operator. Well-known applications of ladder operators in quantum mechanics are in the formalisms of the quantum harmonic oscillator and angular momentum.

<span class="mw-page-title-main">Two-photon circular dichroism</span>

Two-photon circular dichroism (TPCD), the nonlinear counterpart of electronic circular dichroism (ECD), is defined as the differences between the two-photon absorption (TPA) cross-sections obtained using left circular polarized light and right circular polarized light.

In quantum probability, the Belavkin equation, also known as Belavkin-Schrödinger equation, quantum filtering equation, stochastic master equation, is a quantum stochastic differential equation describing the dynamics of a quantum system undergoing observation in continuous time. It was derived and henceforth studied by Viacheslav Belavkin in 1988.

Stimulated Raman spectroscopy, also referred to as stimulated Raman scattering (SRS), is a form of spectroscopy employed in physics, chemistry, biology, and other fields. The basic mechanism resembles that of spontaneous Raman spectroscopy: a pump photon, of the angular frequency , which is scattered by a molecule has some small probability of inducing some vibrational transition, as opposed to inducing a simple Rayleigh transition. This makes the molecule emit a photon at a shifted frequency. However, SRS, as opposed to spontaneous Raman spectroscopy, is a third-order non-linear phenomenon involving a second photon—the Stokes photon of angular frequency —which stimulates a specific transition. When the difference in frequency between both photons resembles that of a specific vibrational transition the occurrence of this transition is resonantly enhanced. In SRS, the signal is equivalent to changes in the intensity of the pump and Stokes beams. The signals are typically rather low, of the order of a part in 10^5, thus calling for modulation-transfer techniques: one beam is modulated in amplitude, and the signal is detected on the other beam via a lock-in amplifier. Employing a pump laser beam of a constant frequency and a Stokes laser beam of a scanned frequency allows for unraveling the molecule's spectral fingerprint. This spectral fingerprint differs from those obtained by other spectroscopy methods, such as Rayleigh scattering, as the Raman transitions confer different exclusion rules than those that apply to Rayleigh transitions.

Hamiltonian truncation is a numerical method used to study quantum field theories (QFTs) in spacetime dimensions. Hamiltonian truncation is an adaptation of the Rayleigh–Ritz method from quantum mechanics. It is closely related to the exact diagonalization method used to treat spin systems in condensed matter physics. The method is typically used to study QFTs on spacetimes of the form , specifically to compute the spectrum of the Hamiltonian along . A key feature of Hamiltonian truncation is that an explicit ultraviolet cutoff is introduced, akin to the lattice spacing a in lattice Monte Carlo methods. Since Hamiltonian truncation is a nonperturbative method, it can be used to study strong-coupling phenomena like spontaneous symmetry breaking.

References

  1. F. Cooper, A. Khare, U. Sukhatme, Supersymmetry in Quantum Mechanics, World Scientific, 2001, Table 4.1
  2. Dahl, J.P.; Springborg, M. (1988). "The Morse Oscillator in Position Space, Momentum Space, and Phase Space" (PDF). The Journal of Chemical Physics. 88 (7): 4535. Bibcode:1988JChPh..88.4535D. doi:10.1063/1.453761. S2CID   97262147.
  3. Lima, Emanuel F de; Hornos, José E M. (2005). "Matrix elements for the Morse potential under an external field". Journal of Physics B. 38 (7): 815–825. Bibcode:2005JPhB...38..815D. doi:10.1088/0953-4075/38/7/004. S2CID   119976840.
  4. 1 2 Le Roy, Robert J.; N. S. Dattani; J. A. Coxon; A. J. Ross; Patrick Crozet; C. Linton (25 November 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.
  5. 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. S2CID   32249407.
  6. 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.
  7. 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. S2CID   26105905.
  8. 1 2 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. S2CID   23016118.
  9. 1 2 3 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. Archived from the original (PDF) on 2019-01-08. Retrieved 2013-11-30.
  10. 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.
  11. Dattani, N. S.; R. J. Le Roy (8 May 2013). "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.
  12. Gunton, Will; Semczuk, Mariusz; Dattani, Nikesh S.; Madison, Kirk W. (2013). "High-resolution photoassociation spectroscopy of the 6Li2A(11Σ+
    u
    ) state". Physical Review A. 88 (6): 062510. arXiv: 1309.5870 . Bibcode:2013PhRvA..88f2510G. doi:10.1103/PhysRevA.88.062510. S2CID   119268157.
  13. 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 c-state". Phys. Rev. A. 87 (5): 052505. arXiv: 1309.6662 . Bibcode:2013PhRvA..87e2505S. doi:10.1103/PhysRevA.87.052505. S2CID   119263860.
  14. 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 1 g (33Π1g ) states". Journal of Chemical Physics. 135 (2): 02403. Bibcode:2011JChPh.135b4303X. doi:10.1063/1.3606397. PMID   21766938.
  15. 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.
  16. 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.
  17. 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 .
  18. 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.
  19. 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.
  20. 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.
  21. 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.
  22. 1 2 3 4 Li, Gang; I. E. Gordon; P. G. Hajigeorgiou; J. A. Coxon; L. S. Rothman (July 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.
  23. 1 2 3 4 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.
  24. 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 XΣg+1 state of Be2". The Journal of Chemical Physics. 140 (6): 064315. Bibcode:2014JChPh.140f4315M. doi:10.1063/1.4864355. PMID   24527923.
  25. 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.
  26. 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. S2CID   2481313.