Unified Impact Theory for Velocity-Changing Effects and Speed Dependencies in Neutral Species Lineshapes

A dipole correlation function which incorporates velocity-changing (motional narrowing) effects and the effects of speeddependent Lorentz relaxation rates into otherwise Voigt profile correlation functions is developed, based partly upon previous work by the author. For the first time simple closed expressions, which lend themselves to elementary calculation beginning only with the relevant parts of intermolecular interaction energies, are developed for the cubic time-dependent term within the exponent describing the decay of the correlation function. This term is of first order in perturber number density, as are the Lorentz parameters, and is complex, thereby allowing for narrowing, changing in shape and asymmetry in the line profile. “Soft” and “hard” collisions play no explicit role, though both are variously present for each line. Quartic time dependencies are also discussed, though they are thought to be negligible in nonhydrogen molecular spectroscopy. Finally, some comments are added about a relevant technique for hydrogen spectra.


Introduction
For several decades, now, neutral species coherence transport and relaxation has followed a number major thrusts, above and beyond the simplest convolution of the optical coherence interruption (Lorentz) and Gaussian (Doppler) contours, as embodied in the Voigt profile for isolated spectral lines.To mention a few, five such developments come immediately to mind: (1) the incorporation of the speed dependence of coherence destruction (with speed dependent Lorentz lineshape parameters) within the Voigt convolution (SDVP effect) [1][2][3]; (2) the incorporation of collisional reduction of free-streaming molecular motions, which had previously been assumed in the initial inclusion of Doppler broadening, known as Dicke (or motional) narrowing [4]; (3) finite duration of collision effects [5] as they influence narrow lines that otherwise are regarded as having been generated by effectively instantaneous collisions; (4) line-mixing effects, in which coherent collisional mixing of the transition of immediate interest with other nearby (frequency-wise) transitions leads to an effective blending of contours in which each becomes distorted into combination Lorentziananomalous dispersion (Fano) contour [5][6][7]; (5) nonlinear spectroscopies [8], where, in two-photon absorption, the Doppler broadening may be cancelled out, due to the opposite propagation direction (and therefore compensatory momentum uptake by the excited molecule) of the two photons involved in producing the absorption signals.Also, in saturation spectroscopies, a primary interaction causes a population "hole" to be produced at a selected frequency corresponding to a specific Doppler velocity group, which can then be broadened directly by velocity-changing collisions for detection by a secondary beam.These are what might be called "special effects" which lie outside of the scope of this article, although the methods described herein are sometimes useful in dealing with these phenomena [9].Line mixing effects can be readily introduced into the formalism that will be described, though they will again be omitted from the present discussion, for simplicity.Finite duration effects fall into two categories, the first being that for a single spectroscopic line the correlation function for nearly zero time (equal to or less than the duration of typical collisions) will differ from the more typical decaying exponential time dependence seen in the longer time behavior.This gives rise to a short time correction to the correlation function which leads to the appearance of an additive (as opposed International Journal of Spectroscopy to multiplicative as in some literature; see [10], e.g.) low intensity, broad frequency contribution that will ensure finite nth moments of the intensity distributions for all n, as must result for analytic forms of the correlation functions.[11][12][13][14] The second finite-duration effect would be described as true binary collisions throughout the collision sequences which give rise to Lorentzian shapes.These would primarily lead to slight quadratic (in density) Lorentzian width and shift dependences.The finite duration-of-collision effects are indeed important for the accurate description of narrow isolated spectral features.Nonetheless, for purposes of simplicity, they again will be omitted from further discussion in the present theoretical treatment.
Of special interest in the present paper are the highly detailed experimental studies of individual spectral lines, such as those carried out by Pine et al. [15][16][17], and Wehr [18,19] for hydrogen halide spectra, by Rohart et al. [20,21] for atmospheric ozone, by Lisak et al. [22] for water vapor and by Priem et al. [23] and Wehr et al. [24,25] for CO.Self-and foreign gas-perturbed hydrogen Raman spectra [26,27] show effects related to those cases mentioned above, though the analysis of the latter lies a bit outside of the present development.Most of the articles cited contain theoretical analyses and there exist many independent theoretical analyses, too numerous to cite here.In the related analysis, there have emerged two models that seem, in varying situations, to be applicable to an understanding of the experimental results: the "soft" collision model first introduced by Galatry [10] in which, in collisions leading to internal state coherence modification, collisions gently act to reduce the free-streaming motions of molecules as assumed in Doppler broadening, together with the introduction of diffusional translational motions, and the "hard" collision model first introduced by Rautian and Sobel'man [28], in which, in addition to internal state coherence modification, molecular motion is totally randomized to the extent allowed by dynamical constraints.To a large extent these models are still in use, though some efforts have been made to bridge the gaps between the two pictures through the introduction of an additional hardness parameter.In actuality, both pictures operate separately but simultaneously, to a greater or lesser extent in producing any given spectral line, in that the distant collisions, which traditionally do not substantially limit internal state coherence but are often the primary cause of line shifts would typically be characterized as being soft, only slightly altering translational motions in any of the numerous such encounters, while the hard collisions would accompany events in which the coherence is often nearly terminated.In the present paper, we will introduce a treatment which formally makes no delineation between the soft and hard collision cases, including both types in a "seamless" unified description of all of the abovementioned processes.Also included in a natural way are all aspects of interference between internal state coherence limitation and velocity changes, which can give rise to slight asymmetries of the lines which have, on occasion, been observed.
One of the difficulties of the descriptions of Dicke narrowing is that there seem not to be any reasonably transparent first principles expressions, or series of expressions giving the relevant parameters associated with these phenomena.In the present paper we shall, for the first time, derive intuitively appealing closed expressions for the relevant quantities, in terms of relevant parts of intermolecular potential functions and intrinsic parameters such as number density and temperature.The present work follows a partial description which was first published in 2000 [29], which carried the capricious title utilizing the phrase "Theory for Everything . .." henceforth referred to as paper I.However, the detailed expressions for relevant parameters occurring in the correlation functions are newly derived in the present paper.While the closely allied SDVP effects have been variously discussed, the first description in which it was systematically related to the time dependence of the correlation function was provided by Looney [30], and that treatment will be followed here, with modifications.
The dipole-dipole (or other optically active operators, depending upon the spectroscopy being described) correlation function for an isolated transition can be described as having an exponential behavior of type [30,31]: In the first term γ = γ r + iγ i represents the standard internal state coherence interruption effect, providing a Lorentz shape half-width at half maximum, γ r , and line shift, −γ i .The second term expresses the Doppler effect which, by itself, leads to the standard Gaussian lineshape with k being the optical wavevector magnitude.The cubic term depends upon the coefficient β, which we might call the first motional narrowing coefficient.It is intrinsically complex, due to correlations of internal phase perturbation and translational velocity changes.Then δ is the second motional narrowing coefficient, and so on.We will supply closed expressions for ready calculation of β, and indicate how one may approach the calculation of δ, in the unlikely event that it would be needed.Higher order terms are required to give a convergent expression for C(t) at very large times, opposing the indefinite increase in the cubic term.But explicit calculation would only proceed to the point that C(t) becomes essentially zero, even with only the cubic term, and the higher terms would remain as being negligible, at which point the Fourier transform integration would be truncated.Higher parameters are of unlikely practical significance in virtually all nonhydrogen spectroscopy.In hydrogen spectroscopy (through Raman, quadrupole or permanent and collisionally induced dipole sharp line spectra in HD) motional narrowing is extremely important [26,27] at moderately low densities (corresponding to around one to ten atmospheres at room temperature).For these cases the above-mentioned higher terms become important.While the present theory surely applies (in a different limit) to these cases, detailed studies will not be carried out at this time.
Nonetheless, an indication of how one might approach the molecular hydrogen problem is suggested.We will assume that classical paths are sufficient for the description of collisional processes.As indicated, one can International Journal of Spectroscopy 3 presumably make the entire calculation quantum mechanically as has been done in the past for the calculation of Lorentz lineshape parameters [32].But there are other possibilities: One might accept fully quantum computation for the Lorentz parameters, accept as true in quantum mechanics as well as classical mechanics the simple Doppler terms, and use classical methods for the description of motional narrowing.If ones interest is solely to understand motional narrowing, measured values of the Lorentz parameters can be accepted, while ensuring that the correct density dependences are strictly adhered to throughout.
A note of caution should be expressed at the present time.This is, that obtaining agreement between theory and experiment on a single molecular line at a single density can lead to deceptive satisfaction.One must strive to analyze all lines that have been observed at all possible pressures, while rigorously adhering to the correct pressure dependencies of the real and imaginary parts of γ and β throughout.

Theoretical Considerations
The development of a correlation function describing isolated spectral lines, omitting duration-of-collision effects, essentially follows that of paper I. Minor changes will be made, to more accurately reflect current notation.In particular, the subscripts " f i" will be dropped from expressions involving the correlation function, the choice of final and initial states in the transition being obvious.Moreover, τ and t will be replaced by t, t , and v o , v will be replaced by v, v , respectively.Vectors will be denoted by boldface type.
The spectral line profile for absorption from a lower state i to upper state f (normalized to the relevant dipole moment matrix element squared) can be represented by the form where C f i (t), the correlation function, satisfies the ensemble average: where k is the optical propagation vector, R is the molecular displacement following time t relative to the position at time t = 0, and U f f and U ii are diagonal time-development matrix elements for the respective internal molecular states resulting from whatever collisional events may occur in the time interval 0 → t.As indicated above, classical paths are envisioned in the present description, though quantum mechanical descriptions of the translational motion problem can be readily developed.For the internal states, U * f f and U ii can be written as product matrix elements for noncorrelated random collisions.However, the same factorization cannot be made for the translational factor.And it is this nonMarkovian property that has possibly inhibited the development which now is presented here.Concentrating on the internal states, for a single collision it is well known [33,34] that the product of U matrix elements can be put into the form f (v; b, u) − S (2)  i (v; b, u) +S * (4)   f + S (4)  i where the superscripts (2), ( 4) indicate the order of collisional inelastic perturbation.The phase shift term can be represented quite exactly to all orders of perturbation [35] while the inelastic terms, S (n) f ,i , often are subject to various mathematical approximations as the strength of the collisional interruption increases [36].While the ATC treatment [33,34] had envisioned the S (n)  f ,i terms as being real, they are in general complex, at times contributing significantly to line shifts, otherwise given solely by the phase shifts [37][38][39].
In the present work, we will regard the calculation of the collisional interruption quantities S(v; b, u) as being completely settled, along with the calculation of outgoing velocities following collisions.We will indicate such dynamical processes as leading to post-collisional velocities v (v; b, u).Here v is the input velocity of the optically active molecule, b being the vector impact parameter lying in the plane perpendicular to (v − u), with u being the perturber initial velocity in any particular encounter.
Following the development of paper I, we now think in terms of a collision as being one in which b lies within a somewhat large radius of collision B, in an attempt to think about probabilities for zero collisions and a time distribution for a first collision.The actual value of B is unimportant, and it will be seen to fall away from the description as the theory is developed.For now, in the interest of simplicity, let us envision the perturbers as being motionless.According to Poisson statistics, for collisions within radius B the probability for there being no collisions in the interval 0 → t is simply while the probability per unit time of a first collision is with N being the perturber number density.Of course, following the first collision there can be a multitude of collisional events, though they will be included implicitly in the appearance of correlation functions for times following the first collision.For a specific initial velocity v, in the absence of any collision, the correlation function will simply be a translational (free-streaming) exponential exp[ik • vt], so that C(v; t) will have a no-collision contribution: But now, in addition, we must include the first collision contribution, distributed throughout the interval 0 → t: where S contains the collisional interruption effects, and v represents the outgoing velocity in this first single collision.
Of course one has to specify the actual impact parameter integration within B, and the integration on u, utilizing the Maxwell distribution of perturber velocities f (u), with the result The integrations over b and u are two-and threedimensional, respectively.This equation contains the fundamental idea of this paper.As yet, however, this equation is somewhat difficult to utilize in a practical way, and the problem of how to specify B remains unsolved.These difficulties can be remedied through the transformations indicated in I and included in this paper, in the appendix.The final equation from the appendix, which has eliminated the troublesome feature of specifying B and which is capable of yielding fundamental expressions for the line shape parameters of current interest is repeated here: keeping in mind that v = v (v; b, u) in accordance with classical dynamics throughout.Note that this equation is equivalent to (10) of paper I.
Over the years, as mentioned above, the collision efficiency function S(v; bu) has been the subject of traditional line shape studies within the classical path formalism.More recently quantum mechanical treatments of the translational motions have been included [32], with considerable success.Within the present theoretical development it is difficult to see how to include quantum translational motions because of the difficulty of computing differential cross-sections for internal state disruption.

Applications
The beauties of (10) are twofold.The first is purely conceptual, in that C(v; t) follows all possible collisional histories, with many possible original and intermediate velocities, for ending at time t having velocity v.But also, by redefining the various time intervals and utilizing the principle of detailed balance, by which the (v, u) → (v , u ) process is equally likely as its reverse process, one can see that C(v; t) also can be viewed as applying to all histories such that one begins at t = 0 with velocity v and ends up with the reverse array of intermediate and final velocity v .The second is the more practical aspect of how to apply this equation to finding correlation functions that are useful in describing the wealth of spectra that can be illustrative of the effects of collisional narrowing.To accomplish this, we will explicitly find the complex t 3 -dependent term in the correlation function in a closed form that represents such an expression for the first time.In addition, we will point the reader in a direction such as to be able to recover the t 4 terms, if the need should arise.We propose using an iterative technique, in which (10) will be solved in successive orders of refinement, such that the first nonVoigt term will be linear in number density of perturbers, the second will be quadratic, and so on.To begin with, let us rewrite (10) in the alternate form But the (1 − S(v; b, u)) term in this equation, by itself, is simply the ordinary pressure broadening, thus The third term in the exponent represents the complete effect of motional narrowing.Denoting C (1) (v; t) as a first iteration correlation function, we now use the collisionless form of the ratio of C's shown in (12) to be exp[−ik following integration on t , as the motional narrowing exponent.In Taylor expansion form, this is representing the lowest order motional narrowing terms which should, in all cases, be computer friendly and which, once computed, should allow one to find C(t) through the t 3 motional narrowing.There are t 3 terms in C(t) arising from speed dependent Voigt profile effects, however, as will be seen in the next section.It is not so easy to proceed further from this point without (elementary) computer analysis.Nonetheless, we can make arguments which may indicate some of the basic properties.Beginning with the first Taylor term we note that integrating only over the angles for b and u, the average velocity difference (v − v ) must take the direction of v itself.
Therefore, this part of the exponent can be written in the form exp[−ik • vXt 2 ], with which is first order in perturber number density N.Although X in the above expression probably varies with speed, in general this will be a small effect compared r.m.s.value of speed effect (and, indeed, with the speed dependence of γ itself, which will be addressed below).At this stage, therefore, From ( 15) it is seen that X is complex, with the imaginary part of β in C(t) ultimately arising from an interference between internal collisional phase shifts and velocity (zcomponent) changes.
The second term in (14), is less obvious in its effects, and the averaging over b and v angles is more subtle.The average will depend upon v x and v y , as well as v z .As such, this term explicitly handles Doppler effects following collision arising from the production of v z components generated from all velocity components present at t = 0.The overall effect of such terms might be quite small, however, as can be seen by reference to typical situations.For collisions between hydrogen halide and other molecules of somewhat similar mass, for the longer range encounters, v will hardly change, leaving (v − v ) as a very small quantity; for the closer collisions that can indeed change v through scattering more toward right angles, for which (v − v ) may be quite large, the collisional interruption given by S(v; b, u) may be so large as to effectively screen the dipole moments from being sensitive to these velocity changes.Of course, for He being the foreign gas, large angle scattering of the optically active molecule will not be much of a factor, leaving (v − v ) as a very small term.Similar arguments might hold for the term linear in (v − v ) analyzed above, but in view of its first order appearance this term might remain the principal source of the t 3 term.To address this question for either case, let us imagine that the smaller angle scattering is most likely to survive the coherence loss dictated by S(v; b, u).Thus, for polar and azimuthal scattering angles ψ, χ, measured form the incident velocity v direction, the term (1/2)k • (v − v )t 2 would be proportional to (1/2)kv z (1 − cos ψ)t 2 , and for small angles, (1/4)kv z ψ 2 t 2 prior to averaging over ψ.Relative to this, the term x +v 2 y )sin 2 ψcos 2 χ).One immediately averages over χ; the first term is quartic in the small scattering angle ψ, and can therefore be ignored.(This shows that the major effect in these terms comes from the v x , v y components of v being collisionally diverted toward the ± z-directions.)Therefore this term is approximately equal to (−k 2 /6)((v 2 x + v 2 y )/2)ψ 2 t 3 .This suggests that the entire exponent for this term can be written as with the factor 2 somewhat compensates the cos 2 χ averaging.In this form, Y will tend to be insensitive to temperature, and is proportional to N.
The overall correlation function C(t) is now found upon averaging over v z , in the form noting that v 2 x + v 2 y = 2k B T/m.The t 4 -term is quadratic in N. By itself it gives an incomplete accounting of the entire t 4 dependence so we will, for present purposes, work only to order t 3 , for which Comparison of small angle approximations reveals that Y is about equal to 2X/3, tending to diminish the motional narrowing.It actually represents a motional-broadening effect associated with collisional production of v z velocity components (from previously v y and v y components), as noted.It is seen that both X and Y are complex quantities, the imaginary parts leading to asymmetries in line shapes; some of which have been observed.
For many spectra the fourth order, t 4 terms in the exponent of C(t) are entirely negligible.If they are needed they will arise from the neglected t 4 terms as noted above, together with the second iteration C (2) in which the results of the first iteration for C (1) (v; t) are included.The calculation will follow the lines already set out, here, for C (1) (v; t), leading to further integrations, for pairs of collisions, in which not only (v, u) → v , but also (v , u ) → v .To work this out in detail is straightforward, though this lies beyond the scope of the present paper.In addition, there is a (neglected) t 4 term found in the next section.
Finally, as previously noted, hydrogen spectra are extremely collisionally narrowed, such that an expansion in terms of individually computed powers of t would seem pointless.In this case, it might be possible to use a model for the collisionally narrowed C(v; t) for all v, then for purposes of calculating a next iteration result, insert the model forms into (10) for use in the ratio C(v ; t)/C(v; t).In this manner one can recover a once-refined result for C(v; t) for all v and t.To the extent that the initial guess was reliable, the corrections within the first refinement result will indeed be small.This would then lead to correlation functions which can then be averaged over v leading to a refined expression for C(t) which in turn can be used directly to obtain the line profile.

Speed Dependent Voigt Profile Effects
For some time now it has been known that the speeddependence of internal coherence relaxation has had importance in determining the shape of spectral lines, shifting them away from the simple Voigt profile [1][2][3]23].Because this lies outside of the major thrust of the present paper we will only sketch the manner in which these effects can provide time-dependent terms in the correlation function C(t).This, again, follows the development of Looney [30].Beginning with the form of C(v; t) leading to the voigt profile, one can make simple changes that will lead to the appropriate time dependences in the exponential decay rate for C(v; t).Let us begin with the simple form of C(v; t), modified so as to include velocity dependences of γ.Let γ(v) depend upon the speed of the optically active molecule in the following manner: as a general expression, v 2 being the value for the mean squared speed, about which we will expand.Subscript N is used to remind us that this derivative is to be considered for constant number density.For present purposes we work within the quadratic approximation, for which we can write with γ = (dγ/dv 2 ) N .The averages over velocity components are easily accomplished, ignoring, for the moment, the exp[−γ o t] term.We begin with the average over v x , using the Maxwell distribution, such that k B being the Boltzmann constant and m the mass of the optically active molecule.We define the constant for notational simplicity, so that C x (t) can be readily obtained, in the form By noting that v 2 x = k B T/m, one can find through expansion and subsequent identification of terms which is exact to terms in t 3 , provided that mod |Zt| < 1.This will always be true for short enough times, and our assumption is that it will be possible to do this for any physically realistic times.
The same steps can be followed for the y-integration, as well.With the z integration, we begin with the same integral, but now with the Doppler (translational) phase factor, exp[ikv z ], leading directly to the result Within terms to order t 3 , this is equivalent to the form Note that Z is proportional to number density, so that the simple Doppler distribution will be altered (in width, e.g.) by terms quadratic in density, while the motional narrowing term will contain parts linear, and varying as the third power in density.While the change in the simple Doppler correlation function will probably be beyond detection, we will retain that term for reference, while dropping the Z 3 term.By now multiplying the x-, yand z-results, and recovering the exp[−γ o t] factor, we find This represents the Voigt correlation function, whose width is supplemented the small (usually negligible) Z 2 t 2 term and by a further exponential t 3 term.While γ , and therefore Z are presumably complex, the speed dependence of line shifts is thought to be minimal, because the increase in collision rate will be roughly matched by corresponding decreases in phase shifts, tending toward cancelation of the speed dependence.Therefore, we will assume Z to be real.As described by Looney [30], the value of γ can be inferred from the temperature dependence of γ, provided that it is known.To see this, note that This holds only if the perturbers are motionless (i.e., have infinite mass).For finite mass perturbers, however, the relative collision speeds will increase more rapidly with temperature, and m must be replaced by the reduced mass of the collision pair, with the result with m x being the mass of the perturbing molecule.This will obviously be most valuable if the temperature dependence of γ is known from separate measurements.

Conclusions
As a result of the above analysis the t 3 terms in the correlation function are given by the contributions coming from velocity-changing collisions (cf.( 20)) together with the above SDVP effect, equation ( 29), giving as a final result with X, Y , and Z being given by ( 15), (18), and (24).Each term is linear in perturber number density N.
To summarize, a formalism has been presented whereby one can find the relevant terms for the dipole-dipole correlation function needed to characterize the line shapes in the impact approximation for isolated (atomic or) molecular lines perturbed by foreign gases.Through standard techniques this formalism could be extended to a fully quantum description of the entirety of the collision problem treated here.In addition, one could treat the problem of line-mixing with neighboring transitions within the present formalism.For purposes of clarity, neither of these steps is presented here.The cubic time-dependent terms of the exponent describing the exponential decay of the correlation function have been explicitly found in terms of averages utilizing only equilibrium velocity distributions and the variously important terms of the perturber interaction with the optically active molecule.Although both "soft" and "hard" types of collisions occur for various impact parameters for almost every line, the present theory makes no fundamental distinction between these limits, and makes a seamless transition from one to the other type of behavior.
For the velocity changing t 3 effects there exist two competing influences: A line narrowing due to the fact that higher speed components tend to be reduced in collisions and a line broadening associated with collisional production of velocity components parallel to the optical propagation vector, derived from pre-collision components perpendicular to this direction.The two effects oppose one another, with the broadening being about two thirds of the narrowing.This means that the velocity-change broadening will not dominate the narrowing, but at the same time it is substantial.In addition, speed-dependent Voigt profile effects are active in giving rise to a further t 3 narrowing exponential term (presuming that the speed dependence of the ordinary Lorentz lineshape parameters shows an increase with increasing velocity, which is the standard case).
In hydrogen molecular spectra the coherence limiting effects are very small.Accordingly one can, at some pressures, see a pronounced line narrowing.One might approach the theoretical problem by using the best model form of C(v; t) (based upon a Dicke type approach), then compute C(v; t) using (10), with the ratio C(v; t )/C(v ; t ) replaced by the model result.C(v; t) will then show small correction terms, if the initial guess is close to being correct. is achieved, which is shown as (10) of the text.Importantly, the actual specification of B is no longer necessary with integration on all physically relevant b values being implied.