On Parity in Diatomic Molecules and Application of a Rigorous Algorithm for the Prediction of Nitric Oxide Spectra

On Parity in Diatomic Molecules and Application of a Rigorous Algorithm for the Prediction of Nitric Oxide Spectra James O Hornkohl1 and Christian G Parigger2* 1Hornkohl Consulting, 344 Turkey Creek Road, Tullahoma, TN, USA 2University of Tennessee, University of Tennessee Space Institute, Center for Laser Applications, Tullahoma, TN, USA International Journal of Molecular and Theoretical Physics Open Access Research Article


Introduction
The diatomic Hamiltonian matrix is historically paritypartitioned, thereby giving parity a more important role in diatomic spectroscopy than in atomic spectroscopy, e.g., see Zare, et al [1].Typically, each spectral line in a highly resolved diatomic spectrum is labeled with the e/f parity designation of the lower level [2] in addition to the branch designation and the value of the J quantum number for the lower level.Several authors, for instance Hougen [3], Røeggen [4], Judd [5] and Larsson [6], have presented treatments of diatomic parity using the Born-Oppenheimer approximate separation of the diatomic eigenfunction into rotational, vibrational and electronic factors.They also used atomic or molecular orbital models of the electronic eigenfunction to determine its parity.An exact formula relating the parity eigenvalue, p = ± 1, to the total angular momentum quantum number, J, is needed.
In this work, the standard framework of parity and angular momentum methods are applied.The effect of parity on the prediction of heteronuclear diatomic molecular spectra is presented.Subsequently, the algorithm for the computation of molecular spectra is described.As a specific example, results for nitric oxide spectra are compared with experimental data.

Methods
Operation of the parity operator on the general Wigner-Witmer [7] diatomic eigenfunction is used to yield the parity eigenvalues that are composed of a constant and an angular momentum dependent part.The computation of diatomic molecular spectra is accomplished without the need of explicitly including parity selection rules.The fundamental Wigner-Witmer diatomic eigenfunction simplifies the determination of rotational line strengths.

Parity operator
The discrete parity operation can be accomplished with a rotation and a reflection.The parity operator, P can be written as a product, The determinant of the matrix representations ( , ) v y z σ and 2

C ( )
x in laboratory xyz-coordinates are -1 and +1, respectively.The 2 C operator is a proper rotation that can be expressed as a discrete transformation of Euler angles.The Euler angles are the arguments of the Wigner D -function used to formulate rotational symmetry, consequently, the eigenvalues of C 2 are controlled by the angular momentum, J .The ( , )   v y z σ operation results in a constant factor, and the 2 C operation yields the angular momentum dependent part of the parity eigenvalue.

Aplastic anemia in systemic lupus erythematosus: A better prognosis acquired aplastic anemia
Copyright: © 2016 Pannu and Varma Angular momentum is the generator of rotations, therefore, one can expect that application of the discrete operator 2 C to the arguments of the function would yield a relationship between angular momentum and parity.In terms of spatial and angular coordinates of diatomic molecules, Equation (2) can be written as Here ρ is the distance of one electron (the electron arbitrarily labeled 1 but it could be any one of the electrons), ζ is the distance of that electron above or below the plane that passes through the center of mass of the two nuclei (the coordinate origin), and χ is the angle of rotation of that electron about the internuclear vector (r, , )  θ ϕ  , , , , , r, , | = , , ' , , ' , r| ( , , ). ( 5) The values of the quantum numbers, J and Ω influence the electronic-vibrational eigenfunction but the electronic-vibrational eigenfunction is not an angular momentum state vector.

Parity of diatomic states
Parity is rotationally invariant.Inversion of the signs of all rotated coordinates inverts the signs of all un-rotated coordinates, and vice versa.Therefore, the parity operator can be represented by 2 C ( ) ( , ).The constant part of the parity eigenvalue, p Σ is labeled in accord with standard spectroscopic notation.The imaginary values of (-1) J occurring when J is half-integer can be avoided if one adopts the convention [2] to always subtract 1/ 2 from J when J is half-integer.With this convention, Eq. ( 6) is replaced by

Parity in an algorithm for computing diatomic spectra
The following describes an algorithm in which Eqs. ( 7) and ( 8) become practical equations for computing diatomic parity.
Consider the algorithm for computation of the wavelengths and intensities in the spectrum of a molecule from the first principles of quantum mechanics.The upper H ′ and lower H Hamiltonian matrices are computed and numerically diagonalized by unitary matrices, U ′ and U .The upper n v J F ′ ′ ′ ′ and lower nvJ F terms are the eigenvalues of the Hamiltonians, † = ( 9) ( : half integer,

Copyright: © 2016 Pannu and Varma
The algorithm for computation of the vacuum wavenumbers, ν ãnd line strengths, ( , ), S nvJ n v J ′ ′ ′ , of diatomic spectral lines is a straightforward application of quantum mechanics, but except for the very simplest molecules is also very far removed from the realm of the possible.However, with two very stringent caveats, the algorithm can be implemented for the diatomic molecule.The first caveat is that the vacuum wavenumbers, ν ~, for many spectral lines in many bands of a band system must have been experimentally measured with high accuracy such as that provided by Fourier transform spectroscopy.Secondly, using semi-empirical molecular constants one must be able to build upper and lower Hamiltonian matrices whose eigenvalue differences accurately predict the measured vacuum wavenumbers.A fitting process is required [1].One assumes trial values for the molecular constants, computes the spectral lines positions, ν ~, and from the differences between finds the corrections to the molecular parameters.The difference between computed and measured line positions will typically equal the measurement error margins.
The Hund's case (a) basis is mathematically complete.A sum of basis functions, | a〉 ,can be quantitatively very accurate.The parity operator, P , commutes with the Hamiltonian.Thus, the orthogonal matrix that diagonalizes the case (a) representation of the Hamiltonian will also diagonalize the case (a) representation of P .

Application to the prediction of nitric oxide spectra
As a specific example, the line position data of Faris and Cosby [10] are used for the NO beta (3,0) band for then purpose of creating a complete line list for the band with line strengths.Figure 1 illustrates a spectrum generated from the NO line list.
A multiphoton 1+1 excitation was used to observe 10 of the 12 possible branches [10], with particular attention to the parity designations of the numerous Λ doublets.These data are particularly suited for testing applications of the algorithm for the calculation of diatomic spectra.A total of 428 lines were fitted with a standard deviation of 0.030 cm -1 cm , and a line list having no missing lines for the range of upper and lower J values was computed.Table 1 provides details of the lines displayed in Figure 1.The computed parity eigenvalues agree with those assigned by Faris and Cosby [10].

Discussion
The exact separation of the coordinates of the total angular momentum in the Wigner D -function greatly simplifies implementation of the algorithm which uses non-vanishing line strengths to determine if a computed term difference represents an allowed spectral line.The Hönl-London factors are computed from the Hund's case (a) transition moment and the matrices,U and U ′ , which diagonalize the upper and lower Hamiltonians.
A single selection rule handles all types of diatomic spectra.If the Hönl-London factor, ( , ) S J J′ is non-vanishing, then the transition is allowed.Parity plays no part in the fitting process which determines the molecular parameters, but the parity eigenvalues are computed from the finalized values of the molecular parameters.The presented algorithm can be used to predict molecular spectra for the purpose of fitting measured data [9].

Conclusions
The essence of this work clearly is the theoretical simplification for the prediction of diatomic molecular spectra.Extensive parity considerations that are presented in classic texts on molecular spectroscopy appear to not be necessary, rather the straightforward application of the Wigner-Witmer function allows one to conclude that a transition is allowed if the rotational line strength is non-vanishing.Parity is not utilized in the determination of the molecular parameters, consequently, the computation and analysis of diatomic molecular spectra is significantly alleviated.

On Parity in Diatomic Molecules and Application of a Rigorous
Algorithm for the Prediction of Nitric Oxide Spectra (3,0) band.Also see Table 1.

r.
The vibrational quantum number v , has been extracted from the quantum numbers collection, n , which represents all required quantum numbers except J , and .The variables ρ , ζ and r are scalars, unaffected by rotations.The physical rotation and the angle of coordinate rotation are also about the same axis, namely the first intermediate y axis of the full coordinate rotation.The angles χ and γ are both rotations about the z′ axis.Thus, In coordinate rotations, one is at liberty to choose , and γ .If one chooses for the angles = The application to the right side of the Wigner-Witmer diatomic eigenfunction (5) yields the parity eigenvalues, wavenumbers,ν ~,of the predicted spectral lines, are term differences.Of the very large number of computed term differences only those for which the Condon and Shortley[8] line strength does not vanish are spectral lines.The line strength, ( , ), S nvJ n v J ′ ′ ′ is the sum over all M and M ′ of the irreducible tensor ( ) separation of the total angular momentum in the Wigner-Witmer diatomic eigenfunction results in a diatomic line strength composed of two parts, the unitless rotational line strength or Hönl-London factor, ( , )S J J′ .The Born-Oppenheimer approximation separates the electronicvibrational strength into electronic and vibrational parts.In the Hund's case (a) basis built from the Wigner-Witmer eigenfunction, the third Euler angle, = , χ γ appears in the Wigner D -function, Parigger, James Hornkoh (2017) On Parity in Diatomic Molecules and Application of a Rigorous Algorithm for the Prediction of Nitric Oxide Spectra.Int J Mol Theor Phy.1(1):1 Int J Hematol Blo Dis 1(1).

Figure 1 :
Figure 1: Section of computed spectrum of the NO Π − Π 2 2 X B Parigger, James Hornkoh (2017) On Parity in Diatomic Molecules and Application of a Rigorous Algorithm for the Prediction of Nitric Oxide Spectra.Int J Mol Theor Phy.1Int J Hematol Blo Dis 1(1).
The value of p Σ does not depend upon quantum numbers.It is a global value applying to all states of a given molecule.If the diatomic molecule can be said to have an intrinsic parity, then it is clearly p Σ .One would expect the product of the intrinsic parities of the fundamental particles composing the molecule to equal p Σ .

Table 1 :
Lines of the NO