Engineering PapersSearch

Engineering topics

Tipping, R. H.

Publications and source records attributed to Tipping, R. H..

At least 19 records

Relaxation Matrix for Symmetric Tops with Inversion Symmetry: Line Coupling and Line Mixing Effects on NH3 Lines in the V4 Band

Line shape parameters including the half-widths and the off-diagonal elements of the relaxation matrix have been calculated for self-broadened NH3 lines in the perpendicular v4 band. As in the pure rotational and the parallel v1 bands, the small inversion splitting in this band causes a complete failure of the isolated line approximation. As a result, one has to use formalisms not relying on this approximation. However, due to differences between parallel and perpendicular bands of NH3, the applicability of the formalism used in our previous studies of the v1 band and other parallel bands must be carefully verified. We have found that, as long as potential models only contain components with K1 equals K2 equals 0, whose matrix elements require the selection rule delta k equals 0, the formalism is applicable for the v4 band with some minor adjustments. Based on both theoretical considerations and results from numerical calculations, the non-diagonality of the relaxation matrices in all the PP, RP, PQ, RQ, PR, and RR branches is discussed. Theoretically calculated self-broadened half-widths are compared with measurements and the values listed in HITRAN 2012. With respect to line coupling effects, we have compared our calculated intra-doublet off-diagonal elements of the relaxation matrix with reliable measurements carried out in the PP branch where the spectral environment is favorable. The agreement is rather good since our results do well reproduce the observed k and j dependences of these elements, thus validating our formalism.

Absorption coefficient

Theoretical Studies of Spectroscopic Line Mixing in Remote Sensing Applications

The phenomenon of collisional transfer of intensity due to line mixing has an increasing importance for atmospheric monitoring. From a theoretical point of view, all relevant information about the collisional processes is contained in the relaxation matrix where the diagonal elements give half-widths and shifts, and the off-diagonal elements correspond to line interferences. For simple systems such as those consisting of diatom-atom or diatom-diatom, accurate fully quantum calculations based on interaction potentials are feasible. However, fully quantum calculations become unrealistic for more complex systems. On the other hand, the semi-classical Robert-Bonamy (RB) formalism, which has been widely used to calculate half-widths and shifts for decades, fails in calculating the off-diagonal matrix elements. As a result, in order to simulate atmospheric spectra where the effects from line mixing are important, semi-empirical fitting or scaling laws such as the ECS (Energy-Corrected Sudden) and IOS (Infinite-Order Sudden) models are commonly used. Recently, while scrutinizing the development of the RB formalism, we have found that these authors applied the isolated line approximation in their evaluating matrix elements of the Liouville scattering operator given in exponential form. Since the criterion of this assumption is so stringent, it is not valid for many systems of interest in atmospheric applications. Furthermore, it is this assumption that blocks the possibility to calculate the whole relaxation matrix at all. By eliminating this unjustified application, and accurately evaluating matrix elements of the exponential operators, we have developed a more capable formalism. With this new formalism, we are now able not only to reduce uncertainties for calculated half-widths and shifts, but also to remove a once insurmountable obstacle to calculate the whole relaxation matrix. This implies that we can address the line mixing with the semi-classical theory based on interaction potentials between molecular absorber and molecular perturber. We have applied this formalism to address the line mixing for Raman and infrared spectra of molecules such as N2, C2H2, CO2, NH3, and H2O. By carrying out rigorous calculations, our calculated relaxation matrices are in good agreement with both experimental data and results derived from the ECS model.

Formalism

Line Mixing in Parallel and Perpendicular Bands of CO2: A Further Test of the Refined Robert-Bonamy Formalism

Starting from the refined Robert-Bonamy formalism [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)], we propose here an extension of line mixing studies to infrared absorptions of linear polyatomic molecules having stretching and bending modes. The present formalism does not neglect the internal degrees of freedom of the perturbing molecules, contrary to the energy corrected sudden (ECS) modeling, and enables one to calculate the whole relaxation matrix starting from the potential energy surface. Meanwhile, similar to the ECS modeling, the present formalism properly accounts for roles played by all the internal angular momenta in the coupling process, including the vibrational angular momentum. The formalism has been applied to the important case of CO2 broadened by N2. Applications to two kinds of vibrational bands (sigma yields sigma and sigma yields pi) have shown that the present results are in good agreement with both experimental data and results derived from the ECS model.

Semiclassical theories

Effects on Calculated Half-Widths and Shifts from the Line Coupling for Asymmetric-Top Molecules

The refinement of the Robert-Bonamy formalism by considering the line coupling for linear molecules developed in our previous studies [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013); 140, 104304 (2014)] have been extended to asymmetric-top molecules. For H2O immersed in N2 bath, the line coupling selection rules applicable for the pure rotational band to determine whether two specified lines are coupled or not are established. Meanwhile, because the coupling strengths are determined by relative importance of off-diagonal matrix elements versus diagonal elements of the operator -iS1 -S2, quantitative tools are developed with which one is able to remove weakly coupled lines from consideration. By applying these tools, we have found that within reasonable tolerances, most of the H2O lines in the pure rotational band are not coupled. This reflects the fact that differences of energy levels of the H2O states are pretty large. But, there are several dozen strongly coupled lines and they can be categorized into different groups such that the line couplings occur only within the same groups. In practice, to identify those strongly coupled lines and to confine them into sub-linespaces are crucial steps in considering the line coupling. We have calculated half-widths and shifts for some groups, including the line coupling. Based on these calculations, one can conclude that for most of the H2O lines, it is unnecessary to consider the line coupling. However, for several dozens of lines, effects on the calculated half-widths from the line coupling are small, but remain noticeable and reductions of calculated half-widths due to including the line coupling could reach to 5%. Meanwhile, effects on the calculated shifts are very significant and variations of calculated shifts could be as large as 25%.

formalism

Two Dimensional Symmetric Correlation Functions of the S Operator and Two Dimensional Fourier Transforms: Considering the Line Coupling for P and R Lines of Linear Molecules

The refinement of the Robert-Bonamy (RB) formalism by considering the line coupling for isotropic Raman Q lines of linear molecules developed in our previous study [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)] has been extended to infrared P and R lines. In these calculations, the main task is to derive diagonal and off-diagonal matrix elements of the Liouville operator iS1 − S2 introduced in the formalism. When one considers the line coupling for isotropic Raman Q lines where their initial and final rotational quantum numbers are identical, the derivations of off-diagonal elements do not require extra correlation functions of the ^S operator and their Fourier transforms except for those used in deriving diagonal elements. In contrast, the derivations for infrared P and R lines become more difficult because they require a lot of new correlation functions and their Fourier transforms. By introducing two dimensional correlation functions labeled by two tensor ranks and making variable changes to become even functions, the derivations only require the latters' two dimensional Fourier transforms evaluated at two modulation frequencies characterizing the averaged energy gap and the frequency detuning between the two coupled transitions. With the coordinate representation, it is easy to accurately derive these two dimensional correlation functions. Meanwhile, by using the sampling theory one is able to effectively evaluate their two dimensional Fourier transforms. Thus, the obstacles in considering the line coupling for P and R lines have been overcome. Numerical calculations have been carried out for the half-widths of both the isotropic Raman Q lines and the infrared P and R lines of C2H2 broadened by N2. In comparison with values derived from the RB formalism, new calculated values are significantly reduced and become closer to measurements.

infrared radiation

Refinement of the Robert-Bonamy Formalism: Considering Effects from the Line Coupling

Since it was developed in 1979, the Robert-Bonamy (RB) formalism has been widely used in calculating pressure broadened half-widths and induced shifts for many molecular systems. However, this formalism contains several approximations whose applicability has not been thoroughly justified. One of them is that lines of interest are well isolated. When these authors developed the formalism, they have relied on this assumption twice. First, in calculating the spectral density F(ω), they have only considered the diagonal matrix elements of the relaxation operator. Due to this simplification, effects from the line mixing are ignored. Second, when they applied the linked cluster theorem to remove the cutoff, they have assumed the matrix elements of the operator exp(-iS(sub 1) - S(sub 2)) can be replaced by the exponential of the matrix elements of -iS(sub 1) - S(sub 2). With this replacement, effects from the line coupling are also ignored. Although both these two simplifications relied on the same approximation, their validity criteria are completely different and the latter is more stringent than the former. As a result, in many cases where the line mixing becomes negligible, significant effects from the line coupling have been completely missed. In the present study, we have developed a new method to evaluate the matrix elements of exp(-iS(sub 1) - S(sub 2)) and have refined the RB formalism such that line coupling can be taken into account. Our numerical calculations of the half-widths for Raman Q lines of the N(sub 2)-N(sub 2) pair have demonstrated that effects from the line coupling are important. In comparison with values derived from the RB formalism, new calculated values for these lines are significantly reduced. A recent study has shown that in comparison with the measurements and the most accurate close coupling calculations, the RB formalism overestimates the half-widths by a large amount. As a result, the refinement of the RB formalism goes in the right direction and these new calculated half-widths become closer to the "true" values.

spectra

Comparison of Trajectory Models in Calculations of N2-broadened Half-widths and N2-induced Line Shifts for the Rotational Band of H2O-16 and Comparison with Measurements

In this work, Complex Robert-Bonamy calculations of half-widths and line shifts were done for N2-broadening of water for 1639 transitions in the rotational band using two models for the trajectories. The first is a model correct to second order in time, the Robert-Bonamy parabolic approximation. The second is the solution of Hamilton's equations. Both models use the isotropic part of the atom-atom potential to determine the trajectories. The present calculations used an intermolecular potential expanded to 20th order to assure the convergence of the half-widths and line shifts. The aim of the study is to assess if the difference in the half-widths and line shifts determined from the two trajectory models is greater than the accuracy requirements of the spectroscopic and remote sensing communities. The results of the calculations are compared with measurements of the half-widths and line shifts. It is shown that the effects of the trajectory model greatly exceed the needs of current remote sensing measurements and that line shape parameters calculated using trajectories determined by solving Hamilton's equations agree better with measurement.

Lamouroux, J.

Causal Correlation Functions and Fourier Transforms: Application in Calculating Pressure Induced Shifts

By adopting a concept from signal processing, instead of starting from the correlation functions which are even, one considers the causal correlation functions whose Fourier transforms become complex. Their real and imaginary parts multiplied by 2 are the Fourier transforms of the original correlations and the subsequent Hilbert transforms, respectively. Thus, by taking this step one can complete the two previously needed transforms. However, to obviate performing the Cauchy principal integrations required in the Hilbert transforms is the greatest advantage. Meanwhile, because the causal correlations are well-bounded within the time domain and band limited in the frequency domain, one can replace their Fourier transforms by the discrete Fourier transforms and the latter can be carried out with the FFT algorithm. This replacement is justified by sampling theory because the Fourier transforms can be derived from the discrete Fourier transforms with the Nyquis rate without any distortions. We apply this method in calculating pressure induced shifts of H2O lines and obtain more reliable values. By comparing the calculated shifts with those in HITRAN 2008 and by screening both of them with the pair identity and the smooth variation rules, one can conclude many of shift values in HITRAN are not correct.

correlation function

Theoretical Studies of N2-broadened Half-widths of H2O Lines Involving High j States

Based on the properties of the energy levels and wave functions of H2O states, one can categorize H2O lines into individually defined groups such that within the same group, the energy levels and the wave functions associated with two paired lines have an identity property while those associated with different pairs have a similarity property. Meanwhile, by thoroughly analyzing processes used to calculate N2-broadened half-widths, it was found that the 'Fourier series' of W(sup a)(sub L(sub 1))(sub K(sub 1))(sub K(sub 1)) (t; j(sub f) T(sub f) and W(sup a)(sub L(sub 1))(sub K(sub 1))(sub K(sub 1)) (t; j(sub i) T(sub i), and a factor P(sub 222) (j(sub f) T(sub f) j(sub i) T(sub i)) are the key items in the Robert-Bonamy formalism to distinguish contributions to ReS2(r(sub c)) among different transitions of j(sub f) T(sub f) - j(sub i). However, these items are completely determined by the energy levels and the wave functions associated with their initial and final states and they must bear the latter's features as well. Thus, it becomes obvious that for two paired lines in the same group, their calculated half-widths must be almost identical and the values associated with different pairs must vary smoothly as their ji values vary. Thus, the pair identity and the smooth variation rules are established within individual groups of lines. One can use these rules to screen half-width data listed in HITRAN and to improve the data accuracies.

water

Uncertainties Associated with Theoretically Calculated N2-Broadened Half-Widths of H2O Lines

With different choices of the cut-offs used in theoretical calculations, we have carried out extensive numerical calculations of the N2-broadend Lorentzian half-widths of the H2O lines using the modified Robert-Bonamy formalism. Based on these results, we are able to thoroughly check for convergence. We find that, with the low-order cut-offs commonly used in the literature, one is able to obtain converged values only for lines with large half-widths. Conversely, for lines with small half-widths, much higher cut-offs are necessary to guarantee convergence. We also analyse the uncertainties associated with calculated half-widths, and these are correlated as above. In general, the smaller the half-widths, the poorer the convergence and the larger the uncertainty associated with them. For convenience, one can divide all H2O lines into three categories, large, intermediate, and small, according to their half-width values. One can use this division to judge whether the calculated half-widths are converged or not, based on the cut-offs used, and also to estimate how large their uncertainties are. We conclude that with the current Robert- Bonamy formalism, for lines in category lone can achieve the accuracy requirement set by HITRAN, whereas for lines in category 3, it 'is impossible to meet this goal.

Ma, Q.

Pair Identity and Smooth Variation Rules Applicable for the Spectroscopic Parameters of H2O Transitions Involving High-J States

Two basic rules (i.e. the pair identity and the smooth variation) applicable for H2O transitions involving high-J states have been discovered. The origins of these rules are the properties of the energy levels and wavefunctions of H2O states with the quantum number J above certain boundaries. As a result, for lines involving high-J states in individually defined groups, all their spectroscopic parameters (i.e. the transition wavenumber, intensity, pressure-broadened half-width, pressure-induced shift, and temperature exponent) must follow these rules. One can use these rules to screen spectroscopic data provided by databases and to identify possible errors. In addition, by using extrapolation methods within the individual groups, one is able to predict the spectroscopic parameters for lines in this group involving very high-J states. The latter are required in developing high-temperature molecular spectroscopic databases such as HITEMP.

Ma, Q.

The Frequency Detuning Correction and the Asymmetry of Line Shapes: The Far Wings of H2O-H2O

A far-wing line shape theory which satisfies the detailed balance principle is applied to the H2O-H2O system. Within this formalism, two line shapes are introduced, corresponding to band-averages over the positive and negative resonance lines, respectively. Using the coordinate representation, the two line shapes can be obtained by evaluating 11-dimensional integrations whose integrands are a product of two factors. One depends on the interaction between the two molecules and is easy to evaluate. The other contains the density matrix of the system and is expressed as a product of two 3-dimensional distributions associated with the density matrices of the absorber and the perturber molecule, respectively. If most of the populated states are included in the averaging process, to obtain these distributions requires extensive computer CPU time, but only have to be computed once for a given temperature. The 11-dimensional integrations are evaluated using the Monte Carlo method, and in order to reduce the variance, the integration variables are chosen such that the sensitivity of the integrands on them is clearly distinguished.

Ma, Q.

The Density Matrix of H20 - N2 In the Coordinate Representation: A Monte Carlo Calculation of the Far-Wing Line Shape

The far-wing line shape theory within the binary collision and quasistatic framework has been developed using the coordinate representation. Within this formalism, the main computational task is the evaluation of multidimensional integrals whose variables are the orientational angles needed to specify the initial and final positions of the system during transition processes. Using standard methods, one is able to evaluate the 7-dimensional integrations required for linear molecular systems, or the 7-dimensional integrations for more complicated asymmetric-top (or symmetric-top) molecular systems whose interaction potential contains cyclic coordinates. In order to obviate this latter restriction on the form of the interaction potential, a Monte Carlo method is used to evaluate the 9-dimensional integrations required for systems consisting of one asymmetric-top (or symmetric-top) and one linear molecule, such as H20-N2. Combined with techniques developed previously to deal with sophisticated potential models, one is able to implement realistic potentials for these systems and derive accurate, converged results for the far-wing line shapes and the corresponding absorption coefficients. Conversely, comparison of the far-wing absorption with experimental data can serve as a sensitive diagnostic tool in order to obtain detailed information on the short-range anisotropic dependence of interaction potentials.

Ma, Q.

Line Coupling in Atmospheric Spectra

The theoretical modeling of atmospheric spectra is important for a number of different applications: for instance, in the determination of minor atmospheric constituents such as ozone, carbon dioxide, CFC's etc.; in monitoring the temperature profile for climate studies; and in measuring the incoming and outgoing radiation to input into global climate models. In order to accomplish the above mentioned goal, one needs to know the spectral parameters characterizing the individual spectral lines (frequency, width, strength, and shape) as well as the physical parameters of the atmosphere (temperature, abundances, and pressure). When all these parameters are known, it is usually assumed that the resultant spectra and concomitant absorption coefficient can then be calculated by a superposition of individual profiles of appropriate frequency, strength and shape. However, this is not true if the lines are 'coupled'. Line coupling is a subtle effect that takes place when lines of a particular molecule overlap in frequency. In this case when the initial states and the final states of two transitions are connected by collisions, there is a quantum interference resulting in perturbed shapes. In general, this results in the narrowing of Q-branches (those in which the rotational quantum number does not change), and vibration-rotational R- and P branches (those in which the rotational quantum number changes by +/- 1), and in the spectral region beyond band heads (regions where the spectral lines pile up due to centrifugal distortion). Because these features and spectral regions are often those of interest in the determination of the abundances and pressure-temperature profiles, one must take this effect into account in atmospheric models.

Tipping, R. H.

A far-wing line shape theory which satisfies the detailed balance principle

A far-wing theory in which the validity of the detailed balance principle is maintained in each step of the derivation is presented. The role of the total density matrix including the initial correlations is analyzed rigorously. By factoring out the rapidly varying terms in the complex-time development operator in the interaction representation, better approximate expressions can be obtained. As a result, the spectral density can be expressed in terms of the line-coupling functions in which two coupled lines are arranged symmetrically and whose frequency detunings are omega - 1/2(omega(sub ji) + omega (sub j'i'). Using the approximate values omega - omega(sub ji) results in expressions that do not satisfy the detailed balance principle. However, this principle remains satisfied for the symmetrized spectral density in which not only the coupled lines are arranged symmetrically, but also the initial and final states belonging to the same lines are arranged symmetrically as well.

Ma, Q.

Calculation of far wing of allowed spectra: The water continuum

A far-wing line shape theory based on the binary collision and quasistatic approximations that is applicable for both the low- and high-frequency wings of allowed vibrational-rotational lines has been developed. This theory has been applied in order to calculate the frequency and temperature dependence of the continuous absorption coefficient for frequencies up to 10,000 cm(exp -1) for pure H2O and for H2O-N2 mixtures. The calculations are made assuming an interaction potential consisting of an isotropic Lennard-Jones part and the leading long-range anisotropic part, and utilizing the measured line strengths and transition frequencies. The results compare well with existing data, both in magnitude and in temperature dependence. This leads us to the conclusion that although dimer and collision-induced absorptions are present, the primary mechanism responsible for the observed water continuum is the far-wing absorption of allowed lines. Recent progress on near-wing corrections to the theory and validations with recent laboratory measurements are discussed briefly.

Tipping, R. H.

Extension of the quasistatic far-wing line shape theory to multicomponent anisotropic potentials

The formalism developed previously for the calculation of the far-wing line shape function and the corresponding absorption coefficient using a single-component anisotropic interaction term and the binary collision and quasistatic approximations is generalized to multicomponent anisotropic potential functions. Explicit expressions are presented for several common cases, including the long-range dipole-dipole plus dipole-quadrupole interaction and a linear molecule interacting with a perturber atom. After determining the multicomponent functional representation for the interaction between the CO2 and Ar from previously published data, we calculate the theoretical line shape function and the corresponding absorption due to the nu(sub 3) band of CO2 in the frequency range 2400-2580 cm(exp -1) and compare our results with previous calculations carried out using a single-component anisotropic interaction, and with the results obtained assuming Lorentzian line shapes. The principal uncertainties in the present results, possible refinements of the theoretical formalism, and the applicability to other systems are discussed briefly.

Ma, Q.

The detailed balance requirement and general empirical formalisms for continuum absorption

Two general empirical formalisms are presented for the spectral density which take into account the deviations from the Lorentz line shape in the wing regions of resonance lines. These formalisms satisfy the detailed balance requirement. Empirical line shape functions, which are essential to provide the continuum absorption at different temperatures in various frequency regions for atmospheric transmission codes, can be obtained by fitting to experimental data.

Ma, Q.