<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JMP</journal-id><journal-title-group><journal-title>Journal of Modern Physics</journal-title></journal-title-group><issn pub-type="epub">2153-1196</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmp.2018.99109</article-id><article-id pub-id-type="publisher-id">JMP-86415</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Theory of the Condon Locus
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jeremy</surname><given-names>B. Tatum</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Physics and Astronomy, University of Victoria, Victoria, Canada</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>jtatum@uvic.ca</email></corresp></author-notes><pub-date pub-type="epub"><day>31</day><month>07</month><year>2018</year></pub-date><volume>09</volume><issue>09</issue><fpage>1735</fpage><lpage>1760</lpage><history><date date-type="received"><day>8,</day>	<month>May</month>	<year>2018</year></date><date date-type="rev-recd"><day>30,</day>	<month>July</month>	<year>2018</year>	</date><date date-type="accepted"><day>2,</day>	<month>August</month>	<year>2018</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p><html>
 <head></head>
 
  The Condon locus for a diatomic molecule is the locus, in the 
  <img src="Edit_d6a8e432-29f6-4de5-b57f-5ffa88a3edfe.bmp" alt="" /> plane, of the strongest bands in an electronic band system. The form of the locus depends upon the form of the potential energy function of the electronic states involved. We show how the locus depends on the potential energy function for simple harmonic and anharmonic oscillators, first from a classical point of view, and then from a quantum mechanical point of view. One phenomenon of interest is that, in the case of anharmonic oscillators, the upper branch of the Condon locus traces much stronger bands than the lower branch. Another phenomenon, predicted by quantum mechanics but not by classical mechanics, is the existence of secondary nested Condon loci.
 
</html></p></abstract><kwd-group><kwd>Diatomic Molecules</kwd><kwd> Franck-Condon Factors</kwd><kwd> Condon Locus</kwd><kwd> Vibrational Constants</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The Franck-Condon factors are factors that give the relative strengths of the many bands in an electronic band system of a diatomic molecule. A band is a transition between the vibrational quantum number v ′ of an upper electronic state and the vibrational quantum number v ″ of a lower electronic state. Thus a band is defined by the number pair ( v ′ , v ″ ) . There is no rigorous selection rule governing such transitions between different electronic states. However, not all such transitions are equally strong. If the Franck-Condon factors are displayed in a table with rows of constant v ′ and columns of constant v ″ , the strongest Franck-Condon factors are often seen to lie on a locus (called the Condon locus) which is roughly a parabola whose symmetry axis makes an angle of about 45˚ to the v ′ and v ″ axes. This angle may be a little more than 45˚, or a little less. The parabola may be quite narrow (small latus rectum) or quite open (large latus rectum).</p><p>A qualitative interpretation for the Franck-Condon loci was originally described by Franck [<xref ref-type="bibr" rid="scirp.86415-ref1">1</xref>] , and a quantum mechanical explanation was developed by Condon [<xref ref-type="bibr" rid="scirp.86415-ref2">2</xref>] . Their physical explanation was described with customary clarity by Herzberg [<xref ref-type="bibr" rid="scirp.86415-ref3">3</xref>] . Roughly it is as follows. At any given instant of time, a molecule is likely to be at its condition of maximum extension or compression, since near these positions in harmonic motion, the speed is slowest. Therefore transitions are most likely to take place from a condition of maximum extension or compression in one electronic state, to a similar condition in the other electronic state. The time taken for an electronic transition is assumed to be very much smaller than the period of vibration of the molecule.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref> represents the potential energy curves of two electronic states and their quantized total energy levels. It will be seen that a vertical line can easily be drawn from v ′ = 1 in the upper state to v ″ = 3 in the lower state. Thus the ( 1 , 3 ) band is likely to be strong. Likewise the (2, 0) band is also likely to be strong.</p><p>In this paper, we examine the form of the Condon locus at four levels of sophistication:</p><p>1) Simple harmonic oscillations, classical treatment.</p><p>2) Anharmonic oscillations, classical treatment.</p><p>3) Simple harmonic oscillations, wave mechanical treatment.</p><p>4) Anharmonic oscillations, wave mechanical treatment.</p><p>In addition to the accounts by Condon and Herzberg cited above, one of the most illuminating papers on this topic of which I am aware is that of Nicholls [<xref ref-type="bibr" rid="scirp.86415-ref4">4</xref>] , and there will be some inevitable overlap between that paper and this. However, I have not seen it treated systematically and quantitatively with quite the approach that I am adopting here. Nicholls described the phenomenon of nested Condon loci, which are not predicted classically but which are a result of the wave-mechanical nature of the molecular system.</p></sec><sec id="s2"><title>2. Simple Harmonic Oscillations, Classical Treatment</title><p>In the simple harmonic treatment, the potential energy V as a function of internuclear separation r is given by an equation of the form</p><p>V = V e + 1 2 k ( r − r e ) 2 . (1)</p><p>It is customary in spectroscopic practice to express energies as “term values” T, which are the energies divided by hc, so that the term values are then dimensionally similar (L<sup>−</sup><sup>1</sup>) to wavenumbers. Thus Equation (1) would customarily be written as</p><p>T = T e + k 2 h c ( r − r e ) 2 . (2)</p><p>Here T e is the electronic contribution to the potential energy. The second term on the right hand side is the elastic contribution (of the vibrating molecule) to the potential energy, r is the internuclear distance, and r<sub>e</sub> is its equilibrium value. The symbol k is the force constant, related to the molecular vibrational constant ω<sub>e</sub> by</p><p>k = 4 π 2 m c 2 ω e 2 (3)</p><p>where m is the “reduced mass” m 1 m 2 m 1 + m 2 of the molecule.</p><p>Mathematically the problem is to draw a horizontal line to intersect the upper curve of <xref ref-type="fig" rid="fig1">Figure 1</xref>; then drop vertical lines from the two points of intersection; and find the two places where these vertical lines intersect the lower curve. This will result in a relation for the strongest Franck-Condon factors in the form of an equation relating the term values T ′ in the upper electronic state to the term values T ″ in the lower electronic state.</p><p>This relation can then be transferred to the ( v ′ , v ″ ) plane via the relations</p><p>T ″ = T ″ e + ( v ″ + 1 2 ) ω ″ e and T ′ = T ′ e + ( v ′ + 1 2 ) ω ′ e (4)</p><p>The derivation and other details of the analysis are given in Hefferlin et al. [<xref ref-type="bibr" rid="scirp.86415-ref5">5</xref>] and are not repeated here. Suffice it to say that the Condon locus in the ( v ′ , v ″ ) plane is a parabola whose equation in a form that is convenient to compute can be written in the form</p><p>Ω ″ ( v ″ + 1 2 ) = ( Δ &#177; Ω ′ ( v ′ + 1 2 ) ) 2 . (5)</p><p>In this equation I make use of a quantity</p><p>L = ℏ π m c (6)</p><p>having the dimensions of a length. This is merely to avoid having to repeat ℏ / ( π m c ) in subsequent equations. If m is expressed in amu, L has the numerical value 6.743052 &#215; 10 − 17 m metres. Then</p><p>Ω ″ = 1 L ω ″ e ,   Ω ′ = 1 L ω ′ e ,   Δ = r ″ e − r ′ e L . (7)</p><p>The properties of the Condon parabola can be traced using any good text on the conic sections (I used my trusty Loney [<xref ref-type="bibr" rid="scirp.86415-ref6">6</xref>] ). The axis of the parabola makes an angle θ with the v ″ axis given by</p><p>tan θ = ω ′ ω ″ (8)</p><p>and the length of the latus rectum is</p><p>2 l = 4 ( r ″ e − r ′ e ) 2 ω ″ e 2 ω ′ e 2 L ( ω ″ e 2 + ω ′ e 2 ) 3 / 2 . (9)</p><p>The parabola is tangent to the lines v ′ = 1 2 and v ″ = 1 2 .</p><p>If the equilibrium internuclear separations in the two electronic states are equal, the parabola degenerates into a straight line. A wide Condon parabola indicates that the internuclear separations in the two electronic states are rather different. I show, in <xref ref-type="table" rid="table1">Table 1</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref>, two examples, in one of which, CN B 2 Σ + − X 2 Σ + , the internuclear distances in the two states are not very different and the Condon parabola is consequently rather narrow; and in the other, AlO A 2 Π r − X 2 Σ + , the internuclear distances are rather different, and the Condon parabola is broad.</p><p>In the classical model, a simple harmonic oscillator is most likely to be found, at some instant of time, at one of the extrema of its motion. The probability that its distance ξ from its equilibrium position will be between ξ and ξ + dξ is proportional to the reciprocal of its speed (by which I mean d r / d t rather than the speed of one of the atoms), which I shall call its “slowness”, s. The time spent in traversing a distance dξ is sdξ. Indeed the probability that the position of the system will, at some instant of time, be in the interval dξ is 2 s d ξ / P = ω s d ξ / π . Here P is the period of the motion, and ω (not to be confused with the vibrational constant ω<sub>e</sub>) is 2π/P. The factor 2 on the left hand side of the equation arises because the displacement from equilibrium passes through ξ twice per period. And, since the speed in simple harmonic motion of amplitude a is ω a 2 − ξ 2 , the probability d Φ , at some instant, that the displacement will be between ξ and ξ + dξ is</p><p>d Φ = 2 d ξ π a 2 − ξ 2 . (10)</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Molecular constants for CN and AlO used in the calculations for <xref ref-type="fig" rid="fig2">Figure 2</xref></title></caption><table><tbody><thead><tr><th align="center" valign="middle" >CN B 2 Σ + − X 2 Σ +</th><th align="center" valign="middle" >AlO A 2 Π r − X 2 Σ +</th></tr></thead><tr><td align="center" valign="middle" >m = 6.46427 amu</td><td align="center" valign="middle" >10.041 amu</td></tr><tr><td align="center" valign="middle" >L = 1.04331 &#215; 10 − 17 m</td><td align="center" valign="middle" >0.35537 &#215; 10 − 17   m</td></tr><tr><td align="center" valign="middle" >r ′ e = 1.1506 &#215; 10 − 10   m</td><td align="center" valign="middle" >1.7708 &#215; 10 − 10   m</td></tr><tr><td align="center" valign="middle" >r ″ e = 1.1718 &#215; 10 − 10   m</td><td align="center" valign="middle" >1.6179 &#215; 10 − 10   m</td></tr><tr><td align="center" valign="middle" >ω ′ e = 2.16413 &#215; 10 5   m − 1</td><td align="center" valign="middle" >7.285 &#215; 10 4   m − 1</td></tr><tr><td align="center" valign="middle" >ω ″ e = 2.068705 &#215; 10 5   m − 1</td><td align="center" valign="middle" >9.7923 &#215; 10 4   m − 1</td></tr><tr><td align="center" valign="middle" >Δ = 2.0320 &#215; 10 5</td><td align="center" valign="middle" >− 4.303 &#215; 10 6</td></tr><tr><td align="center" valign="middle" >Ω ′ = 4.42898 &#215; 10 11</td><td align="center" valign="middle" >  6.85268 &#215; 10 10</td></tr><tr><td align="center" valign="middle" >Ω ″ = 4.63328 &#215; 10 11</td><td align="center" valign="middle" >7.16877 &#215; 10 10</td></tr><tr><td align="center" valign="middle" >2l = 0.1361</td><td align="center" valign="middle" >3.8982</td></tr><tr><td align="center" valign="middle" >θ = 46.29˚</td><td align="center" valign="middle" >36.65˚</td></tr><tr><td align="center" valign="middle" >vertex: ( v ′ , v ″ ) = ( − 0.4775 , − 0.4743 )</td><td align="center" valign="middle" >(0.5510, −0.0672)</td></tr><tr><td align="center" valign="middle" >focus: ( v ′ , v ″ ) = ( − 0.4529 , − 0.4508 )</td><td align="center" valign="middle" >(1.1327, 0.7147)</td></tr></tbody></table></table-wrap><p>The coefficient of dξ, which I may call the probability density, is shown graphically in <xref ref-type="fig" rid="fig3">Figure 3</xref> for a = 1. The integral between ξ = 0 and ξ = a is unity, as befits a probability. The function goes to infinity at ξ = &#177; a , but of course d Φ is everywhere finite (zero) as d ξ → 0 .</p></sec><sec id="s3"><title>3. Anharmonic Oscillations, Classical Treatment</title><p>The vibration of a real molecule will not be simple harmonic, and the curve representing its potential energy as a function of internuclear distance will not be a parabola. For small internuclear distances, when the molecule is compressed, there will be a strong Coulomb repulsion between the nuclei, so the potential energy curve there is steep and negative. For large internuclear distances, the molecule will tend to dissociate, so the potential energy curve asymptotically approaches a dissociation limit. Qualitatively the potential energy function would be expected to look somewhat similar to one of the curves shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>. (The small difference between these curves will be described later, following Equation (15)).</p><p>The same principles apply in forming the Condon locus as in the simple harmonic case, except that the Condon locus will no longer be a parabola. Transitions are most likely to take place near the stationary points (greatest extension or compression) of the vibration as before. However, if one were to imagine a particle sliding without friction to and fro in one of the potential wells of <xref ref-type="fig" rid="fig4">Figure 4</xref>, it is easy to conclude that the particle will spend more time at large r than at</p><p>small r. That means that transitions are more likely to occur when the molecule is at greatest extension than at greatest compression. This leads to the conclusion, by qualitative argument alone and without any numerical calculation, that the upper branch of the Condon locus traces much stronger bands than the lower branch. We shall see later that the wavemechanical treatment leads to the same conclusion.</p><p>Various attempts can be made to devise empirical equations that mimic the expected qualitative potential curve. Two of them are the Morse potential [<xref ref-type="bibr" rid="scirp.86415-ref7">7</xref>] and the Lennard-Jones potential [<xref ref-type="bibr" rid="scirp.86415-ref8">8</xref>] .</p><p>The Morse potential is</p><p>V = D e [ 1 − e − ( r − r e ) / a ] 2 . (11)</p><p>V = potential energy as a function of the internuclear distance r. r<sub>e</sub> = equilibrium internuclear distance. D<sub>e</sub> = dissociation energy. The parameter a has the dimensions of length. What little geometric meaning we can give to a is such that when r − r e = a ,</p><p>V = ( 1 − e − 1 ) 2 D e = 0.399576 D e ≈ 0.4 D e . (12)</p><p>That is to say the extension of the molecule from its equilibrium separation is a when the potential energy is about 40 percent (i.e. 0.399576) of the dissociation energy. Alternatively, when the potential energy is half of the dissociation energy, the extension or compression is</p><p>r − r e = ln ( 2 &#177; 2 ) = 1.2279 a       o r       − 0.5348 a . (13)</p><p>There is another formal solution, namely when r e − r = a ,</p><p>V = ( 1 − e ) 2 D e = 2.952 D e , (14)</p><p>but this is not a physically interesting solution, because if V &gt; D e , the molecule is unstable. If the molecule is compressed by an amount r e − r = a , it will bounce back and dissociate.</p><p>The Lennard-Jones potential is</p><p>V = D e [ 1 + ( r e r ) m − 2 ( r e r ) n ] . (15)</p><p>In <xref ref-type="fig" rid="fig4">Figure 4</xref>, the dashed line is a Lennard-Jones potential with m = 12 and n = 6, while the continuous line is a Morse potential with a = 0.1772. Some slightly tedious algebra will show that the Morse and Lennard-Jones potentials (with m = 12 and n = 6) will have the same full widths at half minimum (FWHm) for</p><p>a = ( 2 + 2 ) 1 / 6 − ( 2 − 2 ) 1 / 6 ln ( 3 + 8 ) = 0.177212908 , (16)</p><p>which is the reason why I chose that value for the Morse parameter in preparing <xref ref-type="fig" rid="fig4">Figure 4</xref>.</p><p>Although the two curves look somewhat similar, only the Lennard-Jones function has the physically desirable characteristic of going to infinity as r → 0 . However, the Morse function is very steep for small r, and, in the example of <xref ref-type="fig" rid="fig4">Figure 4</xref>, it reaches the respectable value of 7.9 &#215; 10 4 at r = 0, so that the function does not give a wholly unreasonable representation of a real potential function. Another attractive feature of the Morse function is that, when it is inserted into the Schr&#246;dinger equation, the eigenfunctions can be written explicitly in terms of algebraic functions (Laguerre polynomials), and, especially, the eigenvalues (vibrational energy levels) are given as</p><p>T = T e + G ( v ) = T e + ω e ( v + 1 2 ) − ω e x e ( v + 1 2 ) 2 , (17)</p><p>with no higher powers of ( v + 1 2 ) . For these reasons I use the Morse function</p><p>in the present analysis of anharmonicity. (In Equation (17), T is the energy divided by hc. That is to say, it is the term value (in m<sup>−</sup><sup>1</sup>) of the level.)</p><p>Further comparisons between these two potential functions can be found in Lim [<xref ref-type="bibr" rid="scirp.86415-ref9">9</xref>] .</p><p>Let us introduce the dimensionless variables</p><p>ξ = ( r − r e ) / a . (18)</p><p>and</p><p>U = V / D e (19)</p><p>Then the Morse function is</p><p>U = ( 1 − e − ξ ) 2 . (20)</p><p>The Taylor expansion of this to ξ<sup>2</sup> is just U = ξ 2 . Thus, to order ξ<sup>2</sup>, the Morse potential is the same as the simple harmonic oscillator potential for which</p><p>V = 1 2 k ( r − r e ) 2</p><p>where k is the force constant, and hence</p><p>a = 2 D e / k ,         k = 2 D e / a 2 . (21)</p><p>The fundamental frequency is</p><p>ν 0 ​ ​ ​ = 1 2π k m = 1 2 π a 2 D e m , (22)</p><p>and the following relations are also of interest:</p><p>ω e = ν 0 c ,             ω e x e = h ν 0 2 4 D e c . (23)</p><p>For large a (small force constant), the graph of the potential energy versus r − r<sub>e</sub> has a wide and shallow minimum. For small a, the graph of the potential energy versus r − r<sub>e</sub> has a sharp and steep minimum. <xref ref-type="fig" rid="fig5">Figure 5</xref> shows a Morse curve with its energy levels.</p><p>In <xref ref-type="fig" rid="fig5">Figure 5</xref>, ξ is r − r e in units of a, and, for reference, the energies of the v th levels and the limits of the motion are given in <xref ref-type="table" rid="table2">Table 2</xref>.</p><p>The Condon locus resulting from transitions between two electronic states whose potential energies are given by Morse functions can be found by the same</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Limits of the motion for the vibrational levels of <xref ref-type="fig" rid="fig5">Figure 5</xref></title></caption><table><tbody><thead><tr><th align="center" valign="middle" >v</th><th align="center" valign="middle" >E<sub>v</sub>/D<sub>e</sub></th><th align="center" valign="middle" >ξ<sub>left</sub></th><th align="center" valign="middle" >ξ<sub>right</sub></th></tr></thead><tr><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0.121093750</td><td align="center" valign="middle" >−0.298611087</td><td align="center" valign="middle" >0.427688129</td></tr><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.339843750</td><td align="center" valign="middle" >−0.459297264</td><td align="center" valign="middle" >0.874575994</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.527343750</td><td align="center" valign="middle" >−0.545913410</td><td align="center" valign="middle" >1.295300309</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0.683593750</td><td align="center" valign="middle" >−0.602564316</td><td align="center" valign="middle" >1.753292606</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >0.808593750</td><td align="center" valign="middle" >−0.641442439</td><td align="center" valign="middle" >2.294799585</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >0.902343750</td><td align="center" valign="middle" >−0.667787197</td><td align="center" valign="middle" >2.994088817</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >0.964843750</td><td align="center" valign="middle" >−0.684239930</td><td align="center" valign="middle" >4.032192797</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >0.996093750</td><td align="center" valign="middle" >−0.692169184</td><td align="center" valign="middle" >6.237346629</td></tr></tbody></table></table-wrap><p>procedure as for simple harmonic functions, as described following Equation (3) in Section 2. Because of the transcendental nature of the equations, it is not possible to arrive at a simple, explicit equation for the Condon locus (which is no longer a parabola). It is straightforward, in any particular case, to carry out the procedure numerically by computer. Transformation from the ( T ′ , T ″ ) plane to the ( v ′ , v ″ ) plane is performed by inversion of Equation (17). The electronic contributions to the energies of the two states do not come into the Condon locus in the ( v ′ , v ″ ) plane, and can both conveniently be taken to be zero.</p><p>I have done the calculations for several cases below, in <xref ref-type="fig" rid="fig6">Figure 6</xref>, where we can see how the nature of the Condon locus varies with the molecular constants of the two states. Except for the anharmonicity constants ω ′ e x ′ e and ω ″ e x ″ e , I have used the molecular constants for the CN B 2 Σ + − X 2 Σ + tabulated in Section 2, <xref ref-type="table" rid="table1">Table 1</xref>. For illustrative purposes I have added various purely fictional values of ω ′ e x ′ e and ω ″ e x ″ e , shown in <xref ref-type="table" rid="table3">Table 3</xref> in order to see their effect on the Condon locus. The dashed curves show the Condon parabola in the simple harmonic approximation, in which ω ′ e x ′ e and ω ″ e x ″ e are both zero. The full lines are the Condon loci when anharmonicity is added. I have drawn the upper arm of the Condon locus in the anharmonic case with a thicker line than the lower arm, to reflect the fact that, as explained in the second paragraph of this Section, the upper arm delineates stronger Franck-Condon factors than the lower arm. The values chosen for the anharmonic constants are indicated in <xref ref-type="table" rid="table3">Table 3</xref>. Since a picture is worth a thousand words, I leave it to the reader to discern the trends that arise from various choices of the anharmonicity constants.</p><p>As in the simple harmonic case, the probability that the position of the system will, at some instant of time, be in the interval dξ is 2 s d ξ / P , where s is the slowness and P is the period, though determining this quantity is slightly less easy that in the simple harmonic case. For illustrative purposes I shall consider a Morse potential of the form given by Equation (11) and I shall determine expressions for the period P of the motion and the slowness s as a function of ξ.</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Anharmonicity constants used for the calculations of <xref ref-type="fig" rid="fig6">Figure 6</xref></title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >ω ″ e x ″ e (m<sup>−</sup><sup>1</sup>)</th><th align="center" valign="middle" >ω ′ e x ′ e (m<sup>−</sup><sup>1</sup>)</th></tr></thead><tr><td align="center" valign="middle" >(a)</td><td align="center" valign="middle" >500</td><td align="center" valign="middle" >500</td></tr><tr><td align="center" valign="middle" >(b)</td><td align="center" valign="middle" >5000</td><td align="center" valign="middle" >5000</td></tr><tr><td align="center" valign="middle" >(c)</td><td align="center" valign="middle" >600</td><td align="center" valign="middle" >300</td></tr><tr><td align="center" valign="middle" >(d)</td><td align="center" valign="middle" >300</td><td align="center" valign="middle" >600</td></tr><tr><td align="center" valign="middle" >(e)</td><td align="center" valign="middle" >6000</td><td align="center" valign="middle" >5000</td></tr><tr><td align="center" valign="middle" >(f)</td><td align="center" valign="middle" >5000</td><td align="center" valign="middle" >6000</td></tr></tbody></table></table-wrap><p>The period P of the motion, in units of a m / D e (a-Morse parameter, m = reduced mass of the molecule, D<sub>e</sub> = dissociation energy), is</p><p>P = 2 [ ∫ ξ left 0 d ξ E v − ( 1 − e − ξ ) 2 + ∫ 0 ξ r i g h t d ξ E v − ( 1 − e − ξ ) 2 ] , (24)</p><p>where E<sub>v</sub> is the total energy in the vth vibrational level. The first integral pertains to the time when the molecule is compressed; the second integral pertains to the time when the molecule is extended. The integration results in</p><p>P = 8 1 − E v [ tan − 1 ( 1 − E v 1 − E v ) + tan − 1 ( 1 + E v 1 − E v ) ] . (25)</p><p><xref ref-type="table" rid="table4">Table 4</xref> shows, for the first eight vibrational levels, the time P<sub>1</sub>, in units of a m / D e , during which the molecule is compressed; the time P<sub>2</sub> during which it is extended; the total period P; and the ratio P<sub>2</sub>/P<sub>1</sub>. It will be seen that this ratio increases with vibrational quantum number. This means that our prediction that the upper arm of the Condon locus delineates stronger Franck-Condon factors than the lower arm is more pronounced at larger quantum numbers, and less pronounced at lower quantum numbers.</p><p>The speed (i.e. d r / d t ) in units of D e / m as a function of internuclear distance is given by (nonrelativistic) energy considerations to be</p><p>d r d t = 2 [ E v − ( 1 − e − ξ ) 2 ] , (26)</p><p>where ξ = ( r − r e ) / a , and the slowness s is the reciprocal of this. The probability density 2s/P is shown in <xref ref-type="fig" rid="fig7">Figure 7</xref> as a function of ξ for the ν = 7 level. The area under the curve is unity. <xref ref-type="table" rid="table4">Table 4</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref> show that the molecule spends more time in extension ( ξ &gt; 1 ) than in compression ( ξ &lt; 1 ) , with the consequence that the upper arm of the Condon locus is stronger than the lower arm; and that for large v the molecule spends much more time in extension than in compression, with the consequence that in practice the lower arm of the Condon locus is likely to be observed only for low ν.</p><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Comparison of the time spent in compression (P<sub>1</sub>) and expansion (P<sub>2</sub>) in each level of <xref ref-type="fig" rid="fig5">Figure 5</xref></title></caption><table><tbody><thead><tr><th align="center" valign="middle" >v</th><th align="center" valign="middle" >P<sub>1</sub></th><th align="center" valign="middle" >P<sub>2</sub></th><th align="center" valign="middle" >P</th><th align="center" valign="middle" >P<sub>2</sub>/P<sub>1</sub></th></tr></thead><tr><td align="center" valign="middle" >0</td><td align="center" valign="middle" >1.833386651</td><td align="center" valign="middle" >2.905688483</td><td align="center" valign="middle" >4.739075134</td><td align="center" valign="middle" >1.584874899</td></tr><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1.650805553</td><td align="center" valign="middle" >3.817358063</td><td align="center" valign="middle" >5.468163616</td><td align="center" valign="middle" >2.312421385</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >1.559318591</td><td align="center" valign="middle" >4.903056592</td><td align="center" valign="middle" >6.462375183</td><td align="center" valign="middle" >3.144358453</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >1.501973790</td><td align="center" valign="middle" >6.396484767</td><td align="center" valign="middle" >7.898458557</td><td align="center" valign="middle" >4.258719301</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >1.463724273</td><td align="center" valign="middle" >8.691436728</td><td align="center" valign="middle" >10.155161002</td><td align="center" valign="middle" >5.937892052</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >1.438305896</td><td align="center" valign="middle" >12.778919506</td><td align="center" valign="middle" >14.217225402</td><td align="center" valign="middle" >8.884702163</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >1.422633874</td><td align="center" valign="middle" >22.272741796</td><td align="center" valign="middle" >23.695375670</td><td align="center" valign="middle" >15.655990062</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >1.415135897</td><td align="center" valign="middle" >69.670991114</td><td align="center" valign="middle" >71.086127011</td><td align="center" valign="middle" >49.232721242</td></tr></tbody></table></table-wrap></sec><sec id="s4"><title>4. Simple Harmonic Oscillations, Wave Mechanical Treatment</title><p>In the wave mechanical model, the probability density, denoted by the symbol ψ<sup>2</sup>, is a wavefunction with several nodes (zeroes) and antinodes (maxima), and which extends slightly beyond the classical limits of the motion. At a given instant of time, the extension of a molecule from its equilibrium position is most likely to be at one of the maxima of ψ<sup>2</sup>. And, since there are several maxima, this gives rise to the possibility that there will be, in the ( v ′ , v ″ ) plane, several nested Condon loci delineating strong Franck-Condon factors. Examples of such nested Frank-Condon factors will be found in Nicholls [<xref ref-type="bibr" rid="scirp.86415-ref4">4</xref>] .</p><p>For simple harmonic potentials the wavefunctions are given in many standard texts, such as that of Eyring, Walter and Kimbal [<xref ref-type="bibr" rid="scirp.86415-ref10">10</xref>] . Normalized to unit area they are</p><p>ψ 2 = 1 2 v v ! π H v 2 ( ξ ) e − ξ 2 l (27)</p><p>where</p><p>l = ℏ m k (28)</p><p>and the H v are the Hermite polynomials. The constant l has dimensions of length. (Recall that ψ 2 d r must be dimensionless, which verifies that Equation (27) balances dimensionally.) Samples of these, for v = 0 , 5 , 10 are shown in Figures 8-10, where the vertical lines are the classical limits of the motion.</p><p>Several features are worthy of comment. For large v the locus of the maxima closely follow the probability density for the classical case illustrated in <xref ref-type="fig" rid="fig3">Figure 3</xref>, and this is often cited as an example of the Bohr correspondence principle. The molecule, at some instant, is most likely to be found near (but not exactly at) the extrema of the motion. The subsidiary maxima are rather lower than the main maxima, which means that the Franck-Condon factors that they delineate will not be as strong as those delineated by the principal maxima. For intermediate ν these characteristic are not at all as pronounced, while for ν = 0 (no vibration)</p><p>the most likely condition of the molecule is not at an extremum of the potential well, but it will most likely be at its equilibrium position, which will surprise no one.</p><p>It is evident that, in order to calculate the Condon loci, we need to know the positions of the maxima of these wavefunctions. These are given in <xref ref-type="table" rid="table5">Table 5</xref>.</p><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> Classical limits, and positions and heights of the maxima of the probability densities (ψ<sup>2</sup>) for the first 11 vibrational levels of a simple harmonic oscillator</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >Classical limits</th><th align="center" valign="middle" >ψ<sup>2</sup> max abscissa</th><th align="center" valign="middle" >ψ<sup>2</sup> max ordinate</th></tr></thead><tr><td align="center" valign="middle" >ν = 0</td><td align="center" valign="middle" >&#177;1.000000</td><td align="center" valign="middle" >0.000000</td><td align="center" valign="middle" >0.564190</td></tr><tr><td align="center" valign="middle" >ν = 1</td><td align="center" valign="middle" >&#177;1.732051</td><td align="center" valign="middle" >&#177;1.000000</td><td align="center" valign="middle" >0.415107</td></tr><tr><td align="center" valign="middle" >ν = 2</td><td align="center" valign="middle" >&#177;2.236068</td><td align="center" valign="middle" >&#177;1.581139</td><td align="center" valign="middle" >0.370492</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.000000</td><td align="center" valign="middle" >0.282095</td></tr><tr><td align="center" valign="middle" >ν = 3</td><td align="center" valign="middle" >&#177;2.645751</td><td align="center" valign="middle" >&#177;2.034075</td><td align="center" valign="middle" >0.345616</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;0.602114</td><td align="center" valign="middle" >0.245552</td></tr><tr><td align="center" valign="middle" >ν = 4</td><td align="center" valign="middle" >&#177;3.000000</td><td align="center" valign="middle" >&#177;2.417686</td><td align="center" valign="middle" >0.328801</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;1.074613</td><td align="center" valign="middle" >0.225993</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.000000</td><td align="center" valign="middle" >0.211571</td></tr><tr><td align="center" valign="middle" >ν = 5</td><td align="center" valign="middle" >&#177;3.316625</td><td align="center" valign="middle" >&#177;2.756849</td><td align="center" valign="middle" >0.316282</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;1.475241</td><td align="center" valign="middle" >0.213132</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;0.476251</td><td align="center" valign="middle" >0.193521</td></tr><tr><td align="center" valign="middle" >ν = 6</td><td align="center" valign="middle" >&#177;3.605551</td><td align="center" valign="middle" >&#177;3.062508</td><td align="center" valign="middle" >0.306401</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;1.828611</td><td align="center" valign="middle" >0.203748</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;0.881604</td><td align="center" valign="middle" >0.181731</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.000000</td><td align="center" valign="middle" >0.176309</td></tr><tr><td align="center" valign="middle" >ν = 7</td><td align="center" valign="middle" >&#177;3.872983</td><td align="center" valign="middle" >&#177;3.344197</td><td align="center" valign="middle" >0.298291</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;2.147928</td><td align="center" valign="middle" >0.196458</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;1.239870</td><td align="center" valign="middle" >0.173179</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;0.406782</td><td align="center" valign="middle" >0.165096</td></tr><tr><td align="center" valign="middle" >ν = 8</td><td align="center" valign="middle" >&#177;4.123106</td><td align="center" valign="middle" >&#177;3.606369</td><td align="center" valign="middle" >0.291445</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;2.441238</td><td align="center" valign="middle" >0.190550</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;1.563978</td><td align="center" valign="middle" >0.166570</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;0.767093</td><td align="center" valign="middle" >0.156984</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.000000</td><td align="center" valign="middle" >0.154271</td></tr><tr><td align="center" valign="middle" >ν = 9</td><td align="center" valign="middle" >&#177;4.358899</td><td align="center" valign="middle" >&#177;3.852560</td><td align="center" valign="middle" >0.285543</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;2.713869</td><td align="center" valign="middle" >0.185615</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;1.861876</td><td align="center" valign="middle" >0.161237</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;1.093513</td><td align="center" valign="middle" >0.150730</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;0.361030</td><td align="center" valign="middle" >0.146450</td></tr><tr><td align="center" valign="middle" >ν = 10</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;4.085357</td><td align="center" valign="middle" >0.280369</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;2.969559</td><td align="center" valign="middle" >0.181399</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;2.138862</td><td align="center" valign="middle" >0.156801</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;1.393823</td><td align="center" valign="middle" >0.145697</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >&#177;0.688554</td><td align="center" valign="middle" >0.140428</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.000000</td><td align="center" valign="middle" >0.138844</td></tr></tbody></table></table-wrap><p>By way of example of the principles involved, I draw in <xref ref-type="fig" rid="fig1">Figure 1</xref>1 two simple harmonic potential energy curves for a fictitious molecule. I am taking the potential energy of the two electronic states involved to be</p><p>V = V e + 1 2 k ( r − r e ) 2 , (29)</p><p>where for the lower state,</p><p>V e = 0 ,       k 2 h c = 2.447 &#215; 10 25 m − 3 ,       r e = 3.9 &#215; 10 − 10 m (30)</p><p>and for the upper state,</p><p>V e h c = 2 &#215; 10 5 m − 1 ,         k 2 h c = 1.625 &#215; 10 25 m − 3 ,       r e = 4.0 &#215; 10 − 10 m . (31)</p><p>In <xref ref-type="fig" rid="fig1">Figure 1</xref>1 I have drawn the first eleven vibrational energy levels in each electronic state, and I have indicated by dots the positions of the maxima of the squares of the wavefunctions. The heights of these maxima can be found from <xref ref-type="table" rid="table5">Table 5</xref> and from Figures 8-10. At a given instant of time the most probable condition of the molecule is not that of greatest extension or compression, but at a separation corresponding to the position of one of the dots. A likely transition (i.e. a large Franck-Condon factor) is one for which a dot on a level in the lower state is vertically beneath a dot on a level in the upper state (on the supposition that the time taken for an electronic transition is much shorter than the vibrational period of the molecule).</p><p>I have also indicated, by dashed curves in the figure the approximate loci of the largest and the second largest maxima. These loci are not exact parabolas, although in the drawing I have indicated the “best” parabolas (quadratic least squares regressions of r upon V) through the dots.</p><p>These loci can now be used to calculate Condon loci in the ( v ′ , v ″ ) plane, which are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>2. The full curve is the Condon locus for simple harmonic oscillations, classical treatment. The dashed curves are the loci of the most prominent bands according to the wavemechanical calculations. In addition to large Franck-Condon factors delineated along the Condon loci, there may also be a few strong bands randomly distributed. For example, reference to <xref ref-type="fig" rid="fig1">Figure 1</xref>1 shows that there will be a strong band joining a secondary maximum at v ′ = 8 to a principal maximum at v ″ = 7 .</p></sec><sec id="s5"><title>5. Anharmonic Oscillations, Wave Mechanical Treatment</title><p>The Morse potential is given by Equation (11). When this is inserted into the Schr&#246;dinger equation, it is well known that the eigenvalues (energy levels) can be written as a series in v + 1 2 containing no powers higher than the second. The energy levels are given by</p><p>E v = h ν 0 [ ( v + 1 2 ) − h ν 0 4 D e ( v + 1 2 ) 2 ] (32)</p><p>In spectroscopic practice, this is usually written as vibrational term values G ( v ) (energy divided by hc) and the vibrational and anharmonicity constants ω e and ω e x e :</p><p>G ( v ) = ω e ( v + 1 2 ) − ω e x e ( v + 1 2 ) 2 (33)</p><p>from which we see that</p><p>ω e = ν 0 c (34)</p><p>and</p><p>ω e x e = h ν 0 2 4 D e c (35)</p><p>Another way of writing equations (32) or (33) is</p><p>E v D e = b [ ( v + 1 2 ) − 1 4 ( v + 1 2 ) 2 ] (36)</p><p>where b (dimensionless) is</p><p>b = h ν 0 D e (37)</p><p><xref ref-type="fig" rid="fig5">Figure 5</xref> was drawn with an arbitrary value of b = 1 4 .</p><p>I also introduce, for convenience, some more quantities as follows:</p><p>E ^ = h 2 2 m a 2 = π 2 c 2 ℏ 2 ω e 2 D e = c 2 h 2 ω e 2 4 D e = 1 4 b 2 D e . (38)</p><p>β = 1 b = D e h ν 0 (39)</p><p>z = 4 β e − ξ (40)</p><p>P = v ! [ 4 β − v − 1 Γ ( v + 1 ) Γ ( 4 β − v ) ] 1 / 2 (41)</p><p>Q = z 2 β − ( v + 1 2 ) e − z / 2 (42)</p><p>R = L v 4 β − 2 v − 1 ( z ) (42)</p><p>where the L are the Laguerre polynomials, generated by</p><p>L v n ( z ) = z − n e z v ! d v d z v ( z n + v e − z ) . (43)</p><p>Of these quantities, E ^ has dimensions of energy, while the others are dimensionless.</p><p>The wavefunctions for the Morse potential are then given by</p><p>ψ v ( z ) = P Q R a − 1 / 2 . (44)</p><p>The quantity P is the normalization factor chosen to ensure that ∫ 0 ∞ ψ 2 d r = 1 . In this paper I am concerned only with the Condon loci and hence only with the positions and relative heights of the maxima of the wavefunctions. Thus I shall be concerned only with the ξ-dependent part of the wavefunctions, namely the product QR. The Morse function of <xref ref-type="fig" rid="fig5">Figure 5</xref> was drawn with b = 1 4 , β = 4 , and I tabulate below, in perhaps a more comprehensible form than Equation (44), the unnormalized wave functions (i.e. the product QR) for the first eight vibrational levels.</p><p>ψ 0 = z 7.5 e − z / 2 (45)</p><p>ψ 1 = ( 14 − z ) z 6.5 e − z / 2 (46)</p><p>ψ 2 = 1 2 ( 156 − 26 z + z 2 ) z 5.5 e − z / 2 (47)</p><p>ψ 3 = 1 6 ( 1320 − 396 z + 36 z 2 − z 3 ) z 4.5 e − z / 2 (48)</p><p>ψ 4 = 1 24 ( 7920 − 3960 z + 660 z 2 − 44 z 3 + z 4 ) z 3.5 e − z / 2 (49)</p><p>ψ 5 = 1 120 ( 30240 − 25200 z + 7200 z 2 − 900 z 3 + 50 z 4 − z 5 ) z 2.5 e − z / 2 (50)</p><p>ψ 6 = 1 720 ( 60480 − 90720 z + 45360 z 2 − 10080 z 3     + 1080 z 4 − 54 z 5 + z 6 ) z 1.5 e − z / 2 (51)</p><p>ψ 7 = 1 5040 ( 40320 − 141120 z + 141120 z 2 − 58800 z 3     + 11760 z 4 − 1176 z 5 + 56 z 6 − z 7 ) z 0.5 e − z / 2 (52)</p><p>By way of example I draw, in Figures 13-15, the squares of the unnormalized wavefunctions ( Q 2 R 2 ) for ν = 0, 4 and 7. The vertical dashed lines indicate the classical limits of the motion.</p><p>As in the classical model, we see that, for large v, the molecule spends much more time in extension than in compression, so that the lower arm of the Condon locus is likely to be observed only for very small v. And, as for the simple harmonic case, for large v, the locus of the maxima becomes more and more similar to the classical “slowness” curve of <xref ref-type="fig" rid="fig7">Figure 7</xref>, providing another example of the correspondence principle.</p><p>The strongest bands in a band system depend immediately on the positions of the maxima of the squares of the wavefunctions. Accordingly I provide, for the first eight vibrational levels, these positions in <xref ref-type="table" rid="table6">Table 6</xref>, in which the entries in normal font are the values of ξ at which ψ<sup>2</sup> is a maximum, and boldface entries are the classical limits of the motion.</p><p>As with the simple harmonic case of Section 4, I draw as <xref ref-type="fig" rid="fig1">Figure 1</xref>6, by way of example of the principles involved, two anharmonic potential energy curves for</p><p>a fictitious molecule. The potential energy curves are calculated with the (arbitrary but realistic) values given in <xref ref-type="table" rid="table7">Table 7</xref>.</p><p>The dots in <xref ref-type="fig" rid="fig1">Figure 1</xref>6 indicate the positions of the maxima of the squares of the wavefunctions. These dots give the most probable extension or compression</p><table-wrap id="table6" ><label><xref ref-type="table" rid="table6">Table 6</xref></label><caption><title> Classical limits, and positions of the maxima of the probability densities (ψ<sup>2</sup>) for the first eight vibrational levels of an anharmonic (Morse) oscillator</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >ν = 0</th><th align="center" valign="middle" >ν = 1</th><th align="center" valign="middle" >ν = 2</th></tr></thead><tr><td align="center" valign="middle" >−0.298611087</td><td align="center" valign="middle" >−0.459297264</td><td align="center" valign="middle" >−0.545913410</td></tr><tr><td align="center" valign="middle" >−0.064538521</td><td align="center" valign="middle" >−0.213854184</td><td align="center" valign="middle" >−0.331747632</td></tr><tr><td align="center" valign="middle" >+0.427688129</td><td align="center" valign="middle" >+0.555024941</td><td align="center" valign="middle" >+0.239376541</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >+0.874575994</td><td align="center" valign="middle" >+0.962385978</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >+1.295300309</td></tr><tr><td align="center" valign="middle" >ν = 3</td><td align="center" valign="middle" >ν = 4</td><td align="center" valign="middle" >ν = 5</td></tr><tr><td align="center" valign="middle" >−0.602564316</td><td align="center" valign="middle" >−0.641442439</td><td align="center" valign="middle" >−0.667787197</td></tr><tr><td align="center" valign="middle" >−0.404447021</td><td align="center" valign="middle" >−0.452953895</td><td align="center" valign="middle" >−0.485304849</td></tr><tr><td align="center" valign="middle" >+0.103244511</td><td align="center" valign="middle" >+0.021148216</td><td align="center" valign="middle" >−0.030959327</td></tr><tr><td align="center" valign="middle" >+0.620563556</td><td align="center" valign="middle" >+0.475860785</td><td align="center" valign="middle" >+0.392400221</td></tr><tr><td align="center" valign="middle" >+1.388382250</td><td align="center" valign="middle" >+1.018545638</td><td align="center" valign="middle" >+0.869710180</td></tr><tr><td align="center" valign="middle" >+1.753292606</td><td align="center" valign="middle" >+1.877286233</td><td align="center" valign="middle" >+1.474643142</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >+2.294799585</td><td align="center" valign="middle" >+2.488684223</td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >+2.994088817</td></tr><tr><td align="center" valign="middle" >ν = 6</td><td align="center" valign="middle" >ν = 7</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >−0.684239930</td><td align="center" valign="middle" >−0.692169184</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >−0.505321775</td><td align="center" valign="middle" >−0.514921216</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >−0.062338833</td><td align="center" valign="middle" >−0.077179042</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >+0.344304585</td><td align="center" valign="middle" >+0.322030283</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >+0.791542138</td><td align="center" valign="middle" >+0.756684695</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >+1.328614250</td><td align="center" valign="middle" >+1.268855283</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >+2.047561579</td><td align="center" valign="middle" >+1.924187023</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >+3.355078812</td><td align="center" valign="middle" >+2.880975731</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >+4.032192797</td><td align="center" valign="middle" >+5.015474119</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >+6.237346629</td><td align="center" valign="middle" ></td></tr></tbody></table></table-wrap><table-wrap id="table7" ><label><xref ref-type="table" rid="table7">Table 7</xref></label><caption><title> Molecular constants used for the calculation of <xref ref-type="fig" rid="fig1">Figure 1</xref>6</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >T<sub>e</sub> (m<sup>−</sup><sup>1</sup>)</th><th align="center" valign="middle" >D<sub>e</sub> (m<sup>−</sup><sup>1</sup>)</th><th align="center" valign="middle" >r<sub>e</sub> (m)</th><th align="center" valign="middle" >a (m)</th><th align="center" valign="middle" >ω<sub>e</sub> (m<sup>−</sup><sup>1</sup>)</th><th align="center" valign="middle" >ω<sub>e</sub>x<sub>e</sub> (m<sup>−</sup><sup>1</sup>)</th></tr></thead><tr><td align="center" valign="middle" >Upper</td><td align="center" valign="middle" >2 &#215; 10 6</td><td align="center" valign="middle" >2.82 &#215; 10 6</td><td align="center" valign="middle" >4.00 &#215; 10 − 10</td><td align="center" valign="middle" >1.546 &#215; 10 − 11</td><td align="center" valign="middle" >3.5 &#215; 10 5</td><td align="center" valign="middle" >1.085 &#215; 10 4</td></tr><tr><td align="center" valign="middle" >Lower</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >2.42 &#215; 10 6</td><td align="center" valign="middle" >4.03 &#215; 10 − 10</td><td align="center" valign="middle" >1.728 &#215; 10 − 11</td><td align="center" valign="middle" >2.9 &#215; 10 5</td><td align="center" valign="middle" >8.689 &#215; 10 3</td></tr></tbody></table></table-wrap><p>at any instant of time. The magnitudes of ψ<sup>2</sup> at each of these position are given in <xref ref-type="table" rid="table6">Table 6</xref> and Figures 13-15. The loci of the strongest maxima are shown by dashed curves.</p><p>The Condon loci have to be determined numerically by the following procedure. We start with some value of the upper vibrational quantum number v ′ and determine the corresponding term value T ′ = V ′ / h c from Equation (17). We then determine the two values of the internuclear distance corresponding to this term value from the dashed loci in the upper state of <xref ref-type="fig" rid="fig1">Figure 1</xref>6. Then determine the term values of the lower state corresponding to these internuclear distances from the dashed curves of the lower state. Finally calculate the corresponding values of v ″ from the converse of Equation (17).</p><p>In <xref ref-type="fig" rid="fig1">Figure 1</xref>7 I show the principal Condon locus as a dashed curve. The full curve shows the “classical” locus for Morse (anharmonic) potential functions-i.e. the locus calculated on the assumption that the most likely internuclear separation at any given instant is that of full compression or full extension, rather than at the maxima of the eigenfunctions.</p><p>I have not shown any secondary loci, because the possibilities are almost endless, quantum numbers of the strong and weak bands being almost random. Strong bands occur wherever the internuclear separation corresponding to a maximum in the eigenfunction in the upper state corresponds with a maximum in the lower state at the same internuclear separation. In terms of <xref ref-type="fig" rid="fig1">Figure 1</xref>6, and expressed more simply, there is a strong band wherever a dot in the upper part of <xref ref-type="fig" rid="fig1">Figure 1</xref>6 is vertically above a dot in the lower part.</p><p>It will be understood from this that small differences in the shapes and positions of the potential curves (i.e. in the equilibrium internuclear distances and the vibrational constants) will result in differences in the positions of the dots</p><p>and hence in the designations of the bands that are likely to be strong or weak. It will also be understood that if the two potential curves are quite similar, secondary and tertiary Condon loci will be much more likely. Indeed, it could be advanced that, if a secondary Condon locus is very evident, it is likely that the internuclear separations of the two electronic states are not very different. Unfortunately in that case the entire electronic band system is likely to be weak because of the small difference in electric dipole moments of the two electronic states, and also, as noted in Section 2, the Condon locus is likely to be narrow and the strongest bands are those in which Δ v is small.</p></sec><sec id="s6"><title>6. Comparison with Observation</title><p>The thrust of this paper has been primarily theoretical. However, calculations of the predicted Condon loci, including in the simple harmonic case the latera recta of the Condon parabolas and the inclinations of their axes, have been carried out for 47 electronic band systems and compared with observations. The results for these are given in Hefferlin et al. [<xref ref-type="bibr" rid="scirp.86415-ref5">5</xref>] .</p></sec><sec id="s7"><title>Acknowledgements</title><p>I would like to thank Rev. Dr Ray Hefferlin of Southern Adventist University, Tennessee, for his encouragement and for exciting my interest in this problem during his continuing studies of the systematics of molecular spectroscopic constants.</p></sec><sec id="s8"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s9"><title>Cite this paper</title><p>Tatum, J.B. (2018) Theory of the Condon Locus. Journal of Modern Physics, 9, 1735-1760. https://doi.org/10.4236/jmp.2018.99109</p></sec></body><back><ref-list><title>References</title><ref id="scirp.86415-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Franck, J. (1925) Transactions of the Faraday Society, 21, 536-542. https://doi.org/10.1039/tf9262100536</mixed-citation></ref><ref id="scirp.86415-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Condon, E.U. (1928) Physical Review, 32, 858-872. https://doi.org/10.1103/PhysRev.32.858</mixed-citation></ref><ref id="scirp.86415-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Herzberg, G. (1950) Molecular Spectra and Molecular Structure. I. Spectra of Diatomic Molecules, Van Nostrand, Princeton, 194-204.</mixed-citation></ref><ref id="scirp.86415-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Nicholls, R.W. (1982) Journal of Quantitative Spectroscopy &amp; Radiative Transfer, 28, 481-492. https://doi.org/10.1016/0022-4073(82)90014-0</mixed-citation></ref><ref id="scirp.86415-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Hefferlin, R., Sackett, J. and Tatum, J.B. (2012) Systematics and Prediction in Franck-Condon Factors, in Quantum Systems in Chemistry and Physics. Springer, Dordrecht, 179-191. https://doi.org/10.1007/978-94-007-5297-9_8</mixed-citation></ref><ref id="scirp.86415-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Loney, S.L. (1895) The Elements of Coordinate Geometry. MacMillan, London, 359-363.</mixed-citation></ref><ref id="scirp.86415-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Morse, P.M. (1929) Physical Review, A34, 57-64. https://doi.org/10.1103/PhysRev.34.57</mixed-citation></ref><ref id="scirp.86415-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Lennard-Jones, J.E. and Dent, B.M. (1926) Proceedings of the Royal Society, A112, 230-234. ttps://doi.org/10.1098/rspa.1926.0107</mixed-citation></ref><ref id="scirp.86415-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Lim, T.-C. (2003) Zeitschrift für Naturforschung, 58a, 615-617. https://doi.org/10.1515/zna-2003-1104</mixed-citation></ref><ref id="scirp.86415-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Eyring, H., Walter, J. and Kimball, G.E. (1944) Quantum Chemistry. Wiley, New York, 75-79.</mixed-citation></ref></ref-list></back></article>