ABSTRACT
Spontaneous Rayleigh-Brillouin (RB) scattering experiments have been performed in air for pressures in the range 0.25–3 bar and temperatures in the range 273–333 K. The functional behaviour of the RB-spectral profile as a function of experimental parameters, such as the incident wavelength, scattering angle, pressure and temperature is analysed, as well as the dependence on thermodynamic properties of the gas, as the shear viscosity, the thermal conductivity, the internal heat capacity and the bulk viscosity. Measurements are performed in a scattering geometry detecting at a scattering angle and an incident wavelength of , at which the Brillouin features become more pronounced than in a right angles geometry and for ultraviolet light. For pressure conditions of 1–3 bar the RB-spectra, measured at high signal-to-noise ratio, are compared to Tenti-S6 model calculations and values for the bulk viscosity of air are extracted. Values of are found to exhibit a linear dependence on temperature over the measurement interval in the range . A temperature dependent value is deduced from a collection of experiments to yield: . These results are implemented in model calculations that were verified for the low pressure conditions (p<1 bar) relevant for the Earth's atmosphere. As a result we demonstrate that the RB-scattering spectral profiles for air under sub-atmospheric conditions can be generated via the Tenti-S6 model, for given gas-phase and detection conditions (p, T, , and θ), and for values for the gas transport coefficients. Spectral profiles for coherent RB-scattering in air are also computed, based on the Tenti-S6 formalism, and the predictions are compared with profiles of spontaneous RB-scattering. Finally data on RB-scattering in air, obtained under a variety of pressure, temperature, wavelength and scattering angles, are analysed in terms of universal scaling, involving the dimensionless uniformity parameter y and the dimensionless frequency x. Such scaling behaviour is shown to be well behaved for a wide parameter space and implies that RB-scattering spectra can be generated for a wide range of atmospheric applications of RB-scattering. The verification of this dimensionless scaling also shows that air can be treated as an ideal gas in the atmospheric regime, where .
GRAPHICAL ABSTRACT
1. Introduction
Light scattering phenomena and the measurement of spectral profiles of Rayleigh-Brillouin (RB) scattered light is of great relevance in a variety of applications [Citation1]. Foremost, studies of the Earth's atmosphere are often based on the analysis of light scattering, e.g. via satellite remote sensing, where the Doppler-wind lidar of the ADM-Aeolus mission is a marked example [Citation2–7]. Temperature profiles in the atmosphere can be probed from RB-scattered light in Lidar configurations [Citation8]. In industrial environments and in combustion light scattering is applied for thermometry [Citation9], for probing gas phase velocity fields [Citation10–12], complex flows [Citation13] and jets [Citation14].
Information on the measurement conditions, such as the density and temperature of the gas, but also on the intrinsic properties of the gas, such as the mass of the molecules, the thermal conductivity, the viscosity and the internal heat capacity is contained in the spectral line shape of the RB-scattered light. That means that some of these parameters can be extracted from dedicated light scattering experiments. So can the temperature of air be determined from the RB-line shape [Citation15], but so can also the viscosity of gases be determined [Citation16, Citation17].
The overall behaviour of RB-spectra can be cast into a generalised description to depend on the uniformity parameter y, a dimensionless parameter dividing the physics of light scattering of a gas into different regimes. For , the scattering is in the Knudsen regime and governed by the Doppler effect, for which the profile can be expressed by a Gaussian function. For the domain of hydrodynamics is entered where a gas is treated as a continuum fluid where a threefold resonance structure is obtained with two fully distinct Brillouin side-peaks and the spectrum represented by three Lorentzian functions [Citation18]. The hydrodynamic approach can be extrapolated to cover the regime of somewhat lower pressures [Citation19]. In the intermediate regime, of , Brillouin side peaks marginally overlap with the central Rayleigh peak, making the profile complex and no exact solutions are available.
Various models have been developed to describe the light scattering spectral line shape in this intermediate, or kinetic regime, which is of relevance for the Earth's atmosphere. Models are expressed in terms of gas transport coefficients, and are generally based on linearization or approximate treatment of the Boltzmann equation representing the collisions in the medium. One approach is the Tenti-model [Citation20, Citation21] providing a rather straightforward representation of the scattering profile in terms of six moments connected to transport coefficients. The Tenti approach has been shown to work well in the case of diatomic and linear molecules [Citation16, Citation22–24], where usually the bulk viscosity is included in the model as a fit parameter. Another approach for the kinetic regime dealt with collisions of molecules treated as rough spheres, which well described the RB-spectrum of SF [Citation25].
Laboratory experiments on RB-scattering of gases typically serve the purpose of verifying and optimising the models and determine the thermodynamic gas transport coefficients that are essential building blocks for models of the light scattering spectral profiles. For this reason a number of RB-scattering studies were performed on air, in view of its relevance for atmospheric science. RB-profiles of air were determined over a variety of pressures and temperatures as well as incident wavelengths (366 nm and 403 nm), where a comparison was made with scattering in pure nitrogen and pure oxygen, while a model description of the Tenti-model was pursued [Citation15, Citation26–28]. Shang et al. [Citation29] measured the spectra of air at an incident wavelength of 532 nm at elevated pressures of 4.0–7.0 bar for varying temperatures, thus extending the parameter space of RB-scattering studies.
Most of the laboratory studies were performed under scattering angles of 90. In the present study the parameter space for RB-scattering in air at 532 nm is widened to the study at a sharper scattering angle, in the forward direction, which bears the advantage that Brillouin side peaks become more pronounced. Also pressure and temperature were varied in the measurements of spectral profiles. The study verifies the validity of the Tenti description for RB scattering in air, and results in values for the bulk viscosity of air.
2. The Rayleigh-Brillouin lineshape
At low pressures in the gas phase, in the regime of mbars or the Knudsen regime, collective effects are absent, light scattering is elastic and can be well approximated as pure Rayleigh scattering. The elastic peak depends on the angle under which the scattered radiation is detected, as illustrated in Figure .Since in light scattering experiments the magnitudes of scattered and incident wave vectors are similar, , with n the refractive index and the incident wavelength, the magnitude of the scattering wave vector is represented by: (1) (1) Rayleigh scattering of a gas involves the thermal motion of the molecules. Depending on the velocity v the detected photon frequency is shifted from the emitted central frequency via: (2) (2) where c is speed of light in vacuum. Under thermal conditions the molecules in a gas exhibit a velocity distribution corresponding to the Maxwell distribution: (3) (3) with the thermal velocity (), the Boltzmann constant, T the temperature, and m the molecular mass. Insertion of this velocity distribution function into Equation (Equation2(2) (2) ) results in a spectral function centred at , with a Gaussian functional form of width: (4) (4) This is the Doppler effect associated with the molecular motion, turning the elastic Rayleigh scattering peak into a spectral profile of Gaussian nature. While the Doppler width depends on the incident wavelength, the molecular mass, and the temperature, it also strongly depends on the scattering angle θ, as illustrated in Figure (a). At smaller scattering angles the scattering profile becomes narrower, and for scattering in the exact forward direction only a coherent wave remains [Citation1]. This overall width of the scattering profile is retained under conditions of collisions in the kinetic regime.
At somewhat higher pressures, including the regime of the Earth's atmosphere at various altitudes, collisions set in to play a role in the light scattering process, giving rise to a fine structure of the Rayleigh line, predicted by Brillouin [Citation30], and independently by Mandel'shtam [Citation31], which will be referred to in the following as the Rayleigh-Brillouin lineshape. Fluctuations in the pressure or concentration of the medium are the source of elastic waves that propagate in all directions. By taking a fixed geometry with an incident vector and detection of a scattered vector in the direction of angle θ, scattering vectors associated with certain elastic waves fulfilling the Bragg condition will be selected for maximum scattering intensity [Citation32]. The magnitude q of the scattering vector is similar to that given in Equation (Equation1(1) (1) ). This form of scattering can also be considered as a Doppler effect, where the change in light frequency is not a result of scattering from a molecule at velocity v, but rather a reflection from a travelling elastic wave propagating through the medium at the sound velocity . This scattering process gives rise to positive and negative frequency shifts in units of angular frequency: (5) (5) Hence the Brillouin shift scales with the angle of observation as shown in Figure .
Under such conditions of increased density or pressure the RB-light scattering spectral profile is determined not only by (T, , m, θ), but by additional thermodynamic properties of the gas, such as the pressure p, the thermal conductivity κ, the shear viscosity , the internal specific heat capacity , and a multitude of relaxation processes occurring in the gas, usually expressed in terms of a bulk viscosity . The latter are known as the macroscopic gas transport coefficients. Air is treated as consisting of a hypothetical single component gas of diatomic molecules of mass , the geometric mean of the masses of nitrogen and oxygen accounting for their abundances, with its experimental gas transport parameters [Citation28]. So the overall width of the scattering feature, mainly determined by , depending on the molecular mass m, is rather similar. It is noted that the minor constituents of air do not have a significant influence on the RB-scattering spectrum of air, as was experimentally verified for water vapour saturated air samples [Citation26].
Collisions in gaseous media are treated by the Boltzmann kinetic equation, which was extended to the case of gases with internal degrees of freedom by Wang Chang and Uhlenbeck [Citation33] to the full WCU kinetic equation for the phase-space distribution function : (6) (6) where is the space position and is the velocity of the particle at time t. Indexes i, j, k and l are internal states of the molecules (v and velocities of first and second molecules) exhibiting elastic and inelastic collisions via state-to-state cross sections σ. As collisions between molecules only cause small deviations from equilibrium, the distribution function of the internal level can be linearised via: (7) (7) with the number density, the Maxwell distribution function, the fraction of molecules in state i and the deviation from equilibrium.
The thermodynamic properties of a gas can then be expressed in a perturbative manner, as deviations from equilibrium in terms of h. Boley et al. [Citation20] developed a model, based on the WCU-equation, involving non-degenerate eigenvectors, that were written as such linearised approximations. These seven eigenvectors are related to seven physical moments (1) the fraction of particles in different internal states, (2) the momentum, (3) the translational energy, (4) the translational heat flux, (5) the internal energy, (6) the internal heat flux and (7) the traceless pressure tensor. This model, later referred to as the Tenti-S7 model, was subsequently truncated into a 6-dimensional model (Tenti-S6) by Tenti et al. [Citation21] by neglecting this seventh parameter. The truncated S6-version of the Tenti-model is considered to be the standard for treating RBS profiles in the kinetic regime. The description in terms of six independent eigenvectors and some detail of the mathematical evaluation is presented in the Supplementary Material, and also a numerical model in the form of a Matlab-code is provided in the Supplementary Material.
With these approximations the WCU-equation can be simplified into a vectorial equation: (8) (8) where is the collision operator, in the form of a matrix (of dimension 6 for the Tenti-S6 model), and each matrix element represents a collision integral. This equation is used to generate a linear system of equations, by expanding h into eigenfunctions of . The six eigenvectors of the Tenti-S6 model can be expressed in terms of three transport coefficients (the shear viscosity , the bulk viscosity and the thermal conductivity κ), the atomic mass m of the particles and the internal specific heat capacity per molecule . In this way the model calculations are intimately connected to well-known and measurable thermodynamic properties of the gas.
Based on this framework Tenti-model codes were developed to straightforwardly calculate the spectral line shape of Rayleigh-Brillouin scattering [Citation20, Citation34–36]. These calculations provide insight into the RB-spectra as a function of relevant variables. In Figure (b) the dependence on the scattering angle θ is again displayed, now also for elevated pressures. While the overall width follows the pattern of the pure Rayleigh scattering, the effect of the Brillouin side peaks is apparent mainly in the forward scattering geometry, for small angles θ, for the highest pressure and the lowest temperatures values. This implies that the Brillouin phenomenon and the testing of models can be most sensitively conducted in a forward scattering geometry. This forms the rationale of choosing a smaller scattering angle, in this study , than in most other studies performed thusfar [Citation15, Citation26–29].
In Figure RB-spectra calculated within this Tenti-S6 formalism are displayed, as a function of incident wavelength, temperature and pressure. The results show that for longer wavelengths, for lower temperatures, and for higher pressures, the Brillouin side peaks become more pronounced. The Brillouin side peaks and their distinctness are, in general terms, associated with the collisional and internal relaxation behaviour of gases. So in order to extract information on these relaxation phenomena occurring in gases, e.g. extracting information on the bulk viscosity , conditions must be chosen under which the RB-profiles undergo a decisive effect. The simulations show that such information can best be obtained from gases at low temperatures and at high pressures. As for the wavelength it would in principle be favourable to use a long wavelength, but a consideration is the strong decrease of the scattering cross section [Citation37, Citation38], scaling with , that deteriorates the signal-to-noise ratio in light scattering experiments. Moreover at longer wavelengths () the quantum efficiency of photon detectors decreases severely. Here a wavelength of 532 nm is chosen since it is well detectable and an often used wavelength in atmospheric lidar experiments.
In Figure simulations of the spectral lineshapes are plotted, for standard instrumental settings of and and for a pressure of 1 bar and room temperature. In the simulations, where a particle mass of m = 29 amu is adopted for air, the values of the thermal conductivity κ, shear viscosity , and the internal heat capacity as well as the ratio of bulk viscosity over shear viscosity are varied over a range around the known literature values. These simulations show that the RB-spectral profiles only marginally depend on the gas transport coefficients. Even for large deviations of the parameters, only small shifts in the line profiles are found. Specifically for the case of the bulk viscosity a variation of 50% in causes a maximum variation of 2% in the Brillouin side peaks.
The thermal conductivity and shear viscosity of air under the given standard conditions are known to good accuracy according to the Sutherland formulas or Sutherland laws [Citation39]: (9) (9) and (10) (10) where , , , , . This leads to a value for the thermal conductivity for air of and a value for the shear viscosity of , both at room temperature. For the bulk viscosity a value of is adopted in the simulations [Citation28].
As was noticed in previous studies on light scattering the values of some transport coefficients depend on the excitation frequency. This holds in particular for the bulk viscosity, which was found to deviate by some four orders of magnitude for the case of CO between measurements performed at acoustic frequencies, i.e. in the MHz range, or light scattering experiments for GHz scattering frequencies [Citation17, Citation22, Citation40]. Also for the case of the thermal conductivity such a frequency dependence was found, specifically for NO [Citation23]. This phenomenon is associated with relaxation of the internal degrees of freedom of the molecule, i.e. rotations and vibrations [Citation17, Citation41]. When the product of the sound frequency and the relaxation time , the corresponding internal degree of freedom will become frozen. Generally, the rotational relaxation time is about seconds, that is, and the rotation will be excited. For N and O, the vibrational relaxation time is larger than seconds at room temperature [Citation42], which means the vibrational degrees remain frozen, both for acoustic (MHz) and light scattering frequencies (GHz). For air, which mainly includes N and O, we suppose that the vibrational degree of freedom is not excited and the thermal conductivities are the same for low and high frequencies. Hence we adopt the values as given by the Sutherland relations and the values reported in Ref. [Citation39]. These values are in agreement for other reported values for the shear viscosity and thermal conductivity of air [Citation43, Citation44].
For completeness also the dependence of the RB-spectra profile on the internal heat capacity is displayed in Figure (c). For air, composed of diatomic molecules with two degrees of rotational relaxation, , where R is the gas constant. A calculation is performed for a triatomic molecular species, with , and otherwise the same gas transport coefficients. The calculation shows a small dependency on molecular composition via its dependence on .
As illustrated in Figures , , and the RB-spectral line shapes exhibit more pronounced Brillouin side peaks for certain conditions where the effects of collisional and internal relaxation phenomena become apparent. For those conditions the information content of the measured spectra is enhanced, allowing for extraction of the gas transport coefficients related to collisional phenomena. The scaling of the sharpness and the pronouncedness of the Brillouin side peaks, here calculated via the Tenti-model, can be approached in a more direct manner by considering that the width of the side peaks is physically connected to damping of elastic waves in the medium, as was analysed by Leontovich [Citation45]. Based on a model for damping of elastic waves in a liquid medium the magnitude of the broadening of the Brillouin peak due to damping was derived (in units of angular frequency): (11) (11) where is the Brillouin shift as defined in Equation (Equation5(5) (5) ) and Γ a damping parameter in the Navier-Stokes equation, characterising the classical dissipation of the acoustic mode [Citation32, Citation46]: (12) (12) with the density of the medium, , and the other thermodynamic parameters as defined before. This approach leads to a ratio of the Brillouin shift over the Brillouin linewidth representing sharpness and the degree to which the Brillouin side peaks are resolved [Citation32]: (13) (13) Even though this model derives its validity from hydrodynamics in essence it will be applicable to RB-scattering in the kinetic regime, thus providing an explanation and a scaling law for the separation of Brillouin side peaks at forward scattering angles θ. In principle the side peaks are most pronounced at the smallest value of θ, but as Fabilinskii already discussed, a further experimental limitation lies in the instrument width of the measurement apparatus [Citation32]. That should be sufficiently small to resolve the Brillouin side peaks at gradually smaller shifts when approaching forward scattering.
Note that the dependence of the RB-spectral lineshape on the measurement conditions (θ, , T and p) is much more pronounced than for the gas transport coefficients (, κ and ). In particular the bulk viscosity , cannot be determined easily from gas-transport experiments and its values are not well known, while Figure (d) shows that the RB-profile is only marginally dependent on the exact structure of the Brillouin side lobes in the scattering spectrum. This implies that spectra of high quality are needed to extract a reasonably accurate value of the bulk viscosity, while at the same time the measurement conditions have to be set or measured at high accuracy.
Concluding from these simulations a strategy can be formulated to measure at conditions of pronounced Brillouin side bands, and then extrapolate its value to those conditions where the side bands are less pronounced, possibly taking into account the scaling of as a function of values for some of the other transport coefficients. That is the rationale for a measurement of RB-spectra at a scattering angle of than scattering at , cf. Figure (b).
The present study focuses on spontaneous Rayleigh-Brillouin (RB) scattering, and most of the experimental studies on RBS have employed this form of light scattering. However, through the availability of narrowband pulsed lasers coherent RBS has been explored as well [Citation47–55]. In such studies an experimental configuration is setup where laser beams drive the density fluctuations. Dipole forces, proportional to the molecular polarisability α, and induced by the pump beams (with amplitudes and ) create a standing-wave lattice at wave vector = ( the wave vectors of the pump fields). This effect produces an acceleration: (14) (14) which adds as an additional term [Citation51] to the kinetic WCU equation, which then reads: (15) (15) These driven, coherent, density fluctuations add to the thermal fluctuations probed in spontaneous RBS, and give rise to a different line profile under conditions of coherent RBS. The line profiles can be calculated via an extended version of the Tenti-S6 model, which was developed by Pan and coworkers [Citation34, Citation51]. Via this formalism RB-spectral profiles are calculated for typical conditions of air (p = 1 bar, ), for a typical wavelength used in cases of pulsed lasers, and for scattering angles and , while the transport coefficients are set at the nominal values for these conditions. In Figure a comparison is made between line shapes for coherent RBS and spontaneous RBS for the same conditions, including the same scattering angle θ. The predicted spectra show that the Brillouin side peaks become more pronounced in the case of coherent RBS. So in principle the underlying collisional and relaxation parameters may be measured at better accuracy compared to spontaneous RBS. However, it is noted that typical experimental configurations for coherent RBS employ a geometry with a scattering angle of close to ; in fact is used in many studies [Citation48, Citation49, Citation53, Citation54]. Under such angles distinctiveness of the side-peaks deteriorates.
3. Experimental
Rayleigh-Brillouin scattering spectra of air are measured in an experimental setup, which was used in a previous study [Citation17] and is schematically shown in Figure . The laser light with wavelength of 532.22 nm, power of 5 Watts, and bandwidth less than 5 MHz travels through a gas cell. For this scattering cell, two Brewster-angled windows are mounted at entrance and exit ports to reduce the loss of incident light and the inside walls were painted black to reduce stray light. A pressure gauge is connected to the cell to monitor the pressure change and a temperature control system consisting of PT-100 sensor, Peltier elements as well as water cooling are used to keep the cell at a constant temperature with uncertainty less than .
Data were recorded in a pressure regime of 0.25–3 bar and for temperatures of 273.2–333.2 K. The forward RB-scattering was recorded under an angle , measured by a home-built goniometer rotation stage. The opening angle is less than determined by the geometry of the scattering gas cell and a slit positioned behind the gas cell. The scattered light propagates through a bandpass filter (Materion, at , bandwidth ) onto a half-confocal Fabry-Perot interferometer (FPI).
The FPI is aligned using a reference laser beam. The radius of the curved mirror is with reflectivity of 99%. The FPI has a nominal free spectral range (FSR) of some 12 GHz, but this turns into an effective FSR if the 4-mode transmission pattern in considered [Citation56]. The precise value of the FSR is determined through frequency-scanning a laser (a narrowband tunable cw-ring dye laser) over more than 1000 modes of the FPI, while measuring the laser wavelength by a wavelength meter (Toptica HighFinesse WSU-30), and yields an uncertainty in the FSR below 1 MHz and an .
The instrument width, yielding a value of (FWHM), is determined by using the reference beam while scanning the piezo-actuated FPI, following methods discussed by Gu et al. [Citation56]. The instrument function is verified to exhibit a functional form of an Airy function, which may be well approximated by a Lorentzian function during data analysis.
RB-scattering spectral profiles were recorded by piezo-scanning the FPI at integration times of 1 s for each step, usually over 18 MHz. A full spectrum covering a large number of consecutive RB-peaks and 10,000 data points were obtained in about 3 h. The piezo-voltage scans were linearised and converted to a frequency scale by fitting the RB-peak separations to the calibrated FSR-value [Citation36, Citation56].
4. Results
The experimental data of spontaneous RBS spectra of air were measured at a wavelength of , a scattering angle of and at various pressures and temperatures. The resulting spectra are split into two subsets, one set obtained at the higher pressures bar, for deducing and verifying the values of the gas-transport coefficients, in particular to extract a value for the bulk viscosity , based on the Tenti-S6 model. The data for these measurements are displayed in Figures , and . A second subset, for data obtained at sub-atmospheric pressure conditions, bar, are then used to verify the obtained model description in the realm of atmospheric applications, to be discussed in Section 5. The conditions under which all data were obtained are listed in Table . The data underlying the Figures – and Figures – are included in the Supplementary Material.
Table 1. Determination of bulk viscosities from RBS measurements on air at and (p, T) experimental conditions, corresponding to uniformity parameters y as indicated.
The data of the first subset were incorporated in a Tenti-S6 model description, for which a Matlab-code was written [Citation36]. Here, the values for the shear viscosity and thermal conductivity κ were taken for the specific measurement temperatures, interpolated from the Sutherland formulas given in Equations (Equation9(9) (9) ) and (Equation10(10) (10) ). A fit was made to extract an optimum value for the bulk viscosity via a least-squares procedure for each (p, T) combination and plotted in Figure as a function of temperature. Also data from previous investigations are included [Citation29, Citation57].
Table 2. Values of .
Inspection of the deduced values of indicates that all values fall within the range , with some pressure dependence, although not systematic. However, the temperature dependencies follow an approximately linear temperature dependence for all individual data sets. To describe this a functional form is chosen to represent this behaviour: (16) (16) where the reference temperature is chosen such that the results for various pressures, and also for various previous experiments exhibit a common origin point, at = 0.86 (0.01). The resulting parameters representing the monotonous increase of the bulk viscosity versus temperature are listed in Table , while the slopes are included in Figure .
In Figure also the resulting values for from previous studies on RBS in air, all measured at a scattering angle of , are included [Citation28, Citation29, Citation57]. The result of Figure shows the spread in results for the bulk viscosity obtained between the present study and three previous studies [Citation28, Citation29, Citation57] The simulations performed in Section 2 and displayed in Figures and show how insensitive the RBS line shape in fact is for the exact value of the bulk viscosity , and its much greater sensitivity to experimental parameters as the scattering angle θ, and the readings of pressure and temperature. As for the temperature reading it may be noted that for all present and previous results a rather smooth temperature functionality is found, indicating that there should be no source of error connected to temperature. In order to deduce a set of values for the bulk viscosity of air an average is taken over the concatenated data sets. Since all data follow the behaviour of Equation (Equation16(16) (16) ), crossing at a common point at the entirety of results can be represented by a single value of . From a weighted average over all data a mean value of /K results, which can be considered as a generic value for the calculation of RBS profiles in air, in conjunction with the other gas transport coefficients.
5. Atmospheric RB spectral profiles
In the previous sections a description is given, based on the Tenti-S6 model on how to compute Rayleigh-Brillouin spectra profiles of air under various conditions and settings of an experiment. The scattering profile depends on some experimental settings as the scattering angle θ and the incident wavelength , as well as on some conditions of the gas, such as temperature and pressure, where air is treated as a hypothetical molecule with mass and . Furthermore the RBS line shape depends on the thermodynamic gas properties of air, the shear viscosity and the thermal conductivity κ, values of which can be determined from generic measurements on the gas [Citation39, Citation43, Citation44]. Here it must be considered that these gas transport coefficients may be temperature dependent, as modelled in the Sutherland formulas, and that in most cases the vibrational relaxation is frozen in air under regular atmospheric conditions. Finally, the last thermodynamic gas property of air, the bulk viscosity can be described by . This combined set of parameters, in combination with the Tenti-S6 model, should provide a universal description for RB-scattering in air, under all experimental settings and gas conditions, the latter with restriction to the kinetic regime.
This model description is now used to compute the RB-spectral profiles under conditions measured in the present experiment for low pressures. Experimental and computed spectra for pressures of 0.25 bar, 0.50 bar and 0.75 bar, and in all cases for a variety of temperatures in the range 273–333 K, are compared in Figures , , and . The spectra show that good agreement is found, well within an error margin of 1% of the peak intensity. In a few cases a small residual in the form of a positive photon signal is found exactly at the centre frequency. This phenomenon was observed in some of our previous studies [Citation15], and can be safely ascribed to spurious Mie reflection of dust particles. Hence, it may be concluded that the set of parameters finds very good agreement with the observed RB-scattering profiles of air measured under atmospheric conditions of p<1 bar.
6. Universal scaling of RB-scattering in air
As is shown in the above the Tenti-S6 model provides a representation of the RB-light scattering function based on the conditions of the scattering process and the specific conditions of the gas. The Tenti-S6 model provides also a means to transfer the scattering profile from one set of conditions to another, a transfer based on the knowledge of the relevant underlying parameters. Such scaling can be connected to the concept of uniformity in kinetic gas theory, where universal scaling parameters are represented in dimensionless coordinates (y, x) to describe the light scattering response of the gas. The dimensionless uniformity parameter y is proportional to the ratio of scattering wavelength Λ (connected to the scattering wave vector as ) to the mean-free-path between collisions: (17) (17) Note that in the community of rarefied gas dynamics the rarefaction parameter is also used, proportional to the uniformity parameter via . With , the thermal velocity of the molecules, and the collisional frequency , p and being the pressure and the shear viscosity of the gas, the uniformity parameter can be written as: (18) (18) It may be verified (see Figures –) that the isolated prominence and the distinctiveness of the Brillouin side peaks is highest for the largest values of the uniformity parameter y. Large y corresponds with higher pressure p, lower temperature T, smaller scattering angle θ and longer wavelength .
Also the frequency scale of the RB scattering profiles can be converted to a reduced and dimensionless frequency scale defined as (19) (19) where q is the scattering vector and , with f the frequency of the incident light beam.
This description in terms of uniformity in a rarefied gas dynamics approach considers the gas as independent of the gas transport coefficients, which relate to internal relaxation mechanisms, and to the intermolecular dynamics connected to a potential describing the interactions. In this context the shear viscosity is a parameter connected to the kinetics describing the molecular collisions in terms of a mean-free-path. Such scaling in terms of dimensionless parameters was discussed and tested for gases with atomic constituents [Citation58–60]. The scaling is expected to hold for densities for which the Boltzmann equation may be applied, thus where the mean field or average interaction between molecules may be neglected.
Even though the description of the light scattering process and its resulting spectral profile can be connected to dimensionless (y, x) coordinates, it is still needed to exploit an underlying model for generating spectra and for demonstrating the scaling of spectra. Such procedure is here undertaken with application of the Tenti-S6 model, where it must be considered that the Tenti-model involves aspects of gas transport. To avoid the effects of transport coefficients on the RBS-line shape their values are taken at fixed nominal values (, , and ), even though these coefficients may be temperature dependent, or be subject to freezing out of internal relaxation modes. Fixing the values of these thermodynamic properties implies that they do not affect the dimensionsless or universal scaling.
The outcome of such scaling is shown in Figure , where calculated spectra are displayed for a single value of the dimensionless uniformity parameter y = 1.10. Note that these RB-scattering spectra, plotted on a common dimensionless x-axis, are calculated for different combinations of , θ, p and T. While wavelengths (366 nm, 403 nm, and 532 nm) and scattering angles (, ) are taken to coincide with those used in experiments, values for p and T are chosen such as to match y = 1.10. In these calculations no convolution with an instrument function is included. The results plotted in Figure for four different experimental conditions indeed show identical RB-scattering profiles for a common uniformity parameter. For a pair of dimensionless (y, x) values the RB-scattering spectra are not influenced by the incident wavelength and scattering angle, nor of a set of p and T. Hence, the RB-spectral can be expressed as a function of y and x alone. This demonstration of dimensionless scaling of calculated spectra verifies that spectra, recorded in one set of conditions can be transferred to another set of conditions.
This universal scaling was investigated and tested for experimental RB-spectra of noble gases [Citation58–60]. The scaling in those studies was done for constant scattering conditions (θ, ), while combinations of gaseous settings (p, T) were chosen as to match a similar value of y. In these studies spectra for different noble gases He, Ne and Ar were compared, which made it necessary to invoke a slight rescaling of the reduced frequency axis for the Brillouin shift . As a result of the experimental comparisons it was concluded that the RB-spectra in the kinetic regime grossly follow ideal gas conditions and that they are rather insensitive to details of the intermolecular potential. However, small deviations were found, in particular at higher densities with gradual increase of deviations at larger values of y, that were attributed to variations in the sound velocity and thermal conductivity [Citation59].
In the present study we extend this comparison of dimensionless scaling laws for RB-scattering to the case of air. In this specific case the difference between atoms and molecules needs to be considered, where molecules undergo internal relaxation between the modes of translation, vibration and rotation, which is expressed in terms of a bulk viscosity parameter . Experimentally the present data set obtained at and , and various values of p and T, leading to uniformity parameters y as listed in Table , are compared with previously obtained data for RB-scattering in air at [Citation57] and [Citation28]. Hence the scaling comparison covers all (, θ, p, T) parameters.
The data sets of RB-spectra contain examples where the uniformity parameters y are near-coincident, while the underlying physical conditions are vastly different, thus allowing for an experimental comparison of scaling, while remaining in the kinetic regime of . In Figure a comparison is made between spectra recorded for and , and other parameters as indicated, resulting in a uniformity parameter of . In Figure a comparison is made for two cases also for a near-coincident uniformity parameter for wavelengths and .
The comparisons of Figures and demonstrate that the dimensionless scaling of the spectra at the two values of y produces excellent agreement. The residuals indicate agreement better than 1.5% for the maximum differences in both cases, where it should be noted that the y-parameters do not match exactly, in both cases differing by 0.02. Moreover the instrument widths used in the different experiments were different. When converted to the reduced x-axis scale the instrument widths were , , and for the present experimental study .
The aptness of the universal scaling approach for the case of air demonstrates that under the conditions of the included experiments, mainly set by the range of the uniformity parameter , air very much behaves as an ideal gas. Real gas effects as a result of molecular interactions, relaxation phenomena, freezing out of certain modes, and represented in terms of a bulk viscosity parameters, which may be dependent on temperature, only marginally affect the scaling and the shape of the Rayleigh-Brillouin spectra. Indeed, from the a priori simulations shown in Figures – it is evident that the effect of small changes at the few %-level in the bulk viscosity or thermal conductivity is only marginal and not detectable.
7. Conclusion
We have investigated Rayleigh-Brillouin scattering in air over a range of pressures and temperatures, measuring high quality data at a signal-to-noise ratio well within the 1% level of peak intensities. Based on an a priori modelling a smaller scattering angle () was chosen for the experimental geometry, because under such condition the Brillouin side peaks become more pronounced. From the higher pressure data bar values of the bulk viscosity parameter were derived. From these measurements, and including information from previous studies, a temperature dependent functional form is established for this elusive gas transport coefficient .
Based on the Tenti-S6 model description and a set of reliable gas transport coefficients for air a universal framework is established to predict Rayleigh-Brillouin scattering profiles under various density and temperature conditions of atmospheric air, and for different settings of the scattering geometry and incident wavelength. This framework holds for various scattering geometries and experimental parameters, such as the incident wavelength. The available data thus are at the basis of a predictive tool for delivering RB-spectra of air, for application in atmospheric sensing.
Finally from the now available experimental data sets for RB-scattering a universal scaling is tested relating spectra to a dimensionless uniformity parameter and a reduced frequency scale. Very good agreement is found for the regime of , demonstrating that in the range of atmospheric pressures and temperatures air very much behaves like an ideal gas.
Supplemental Material
Download Zip (3 MB)Acknowledgments
This research was supported by the China Exchange Program jointly run by the Netherlands Royal Academy of Sciences (KNAW) and the Chinese Ministry of Education. YW acknowledges support from the Chinese Scholarship Council (CSC) for his stay at VU Amsterdam. The authors wish to thank Willem van de Water (Delft University) for fruitful discussions.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- R.B. Miles, W.R. Lempert and J.N. Forkey, Meas. Sci. Technol. 12, R33 (2001). doi:10.1088/0957-0233/12/5/201.
- A. Stoffelen, G.J. Marseille, F. Bouttier, D. Vasiljevic, S. de Haan and C. Cardinali, Q. J. R. Meteorol. Soc. 132, 1927 (2006). doi:10.1256/qj.05.83.
- O. Reitebuch, C. Lemmerz, E. Nagel, U. Paffrath, Y. Durand, M. Endemann, F. Fabre and M. Chaloupy, J. Atmos. Ocean. Tech. 26, 2501 (2009). doi:10.1175/2009JTECHA1309.1.
- U. Marksteiner, C. Lemmerz, O. Lux, S. Rahm, A. Schäfler, B. Witschas and O. Reitebuch, Remote Sens. 10, 2056 (2018). doi:10.3390/rs10122056.
- M. Šavli, J. de Kloe, G.-J. Marseille, M. Rennie, N. Žagar and N. Wedi, Q. J. R. Meteorol. Soc. 145, 3499 (2019). doi:10.1002/qj.v145.725. doi: 10.1002/qj.3634
- O. Lux, C. Lemmerz, F. Weiler, U. Marksteiner, B. Witschas, S. Rahm, A. Geiß and O. Reitebuch, Atmos. Meas. Techn. 13, 2075 (2020). doi:10.5194/amt-13-2075-2020.
- B. Witschas, C. Lemmerz, A. Geiß, O. Lux, U. Marksteiner, S. Rahm, O. Reitebuch and F. Weiler, Atmos. Meas. Techn. 13, 2381 (2020). doi:10.5194/amt-13-2381-2020.
- B. Witschas, C. Lemmerz and O. Reitebuch, Opt. Lett. 39, 1972 (2014). doi:10.1364/OL.39.001972.
- J.A. Sutton and J.F. Driscoll, Opt. Lett. 29, 2620 (2004). doi:10.1364/OL.29.002620.
- J. Forkey, W. Lempert and R. Miles, Exp. Fluids 24, 151 (1998). doi:10.1007/s003480050162.
- M. Boguszko and G.S. Elliot, Exp. Fluids 38, 33 (2005). doi:10.1007/s00348-004-0881-4.
- U. Doll, G. Stockhausen and C. Willert, Opt. Lett. 42, 3773 (2017). doi:10.1364/OL.42.003773.
- R.B. Miles and and W.R. Lempert, Ann. Rev. Fluid Mech. 29, 285 (1997). doi:10.1146/annurev.fluid.29.1.285.
- F. Benhassen, M.D. Polanka and M.F. Reeder, J. Aerosp. Eng. 30, 04016067 (2017). doi:10.1061/(ASCE)AS.1943-5525.0000655.
- B. Witschas, Z.Y. Gu and W. Ubachs, Opt. Express 22, 29655 (2014). doi:10.1364/OE.22.029655.
- Z. Gu and W. Ubachs, Opt. Lett. 38, 1110 (2013). doi:10.1364/OL.38.001110.
- Y. Wang, W. Ubachs and W. van de Water, J. Chem. Phys. 150, 154502 (2019). doi:10.1063/1.5093541.
- R.D. Mountain, Rev. Mod. Phys. 38, 205 (1966). doi:10.1103/RevModPhys.38.205.
- C.M. Hammond and T.A. Wiggins, J. Chem. Phys. 65, 2788 (1976). doi:10.1063/1.433425.
- C.D. Boley, R.C. Desai and G. Tenti, Can. J. Phys. 50, 2158 (1972). doi:10.1139/p72-286.
- G. Tenti, C.D. Boley and R.C. Desai, Can. J. Phys. 52, 285 (1974). doi:10.1139/p74-041.
- Z. Gu, W. Ubachs and W. van de Water, Opt. Lett. 39, 3301 (2014). doi:10.1364/OL.39.003301.
- Y. Wang, K. Liang, W. van de Water, W. Marques and W. Ubachs, J. Quant. Spectros. Ra. 206, 63 (2018). doi:10.1016/j.jqsrt.2017.10.029.
- J. Shang, T. Wu, C. Yang, R. Hu, X. He and Z. Chen, Opt. Comm. 436, 127 (2019). doi:10.1016/j.optcom.2018.12.017.
- Y. Wang, Y. Yu, K. Liang, W. Marques, W. van de Water and W. Ubachs, Chem. Phys. Lett. 669, 137 (2017). doi:10.1016/j.cplett.2016.12.033.
- B. Witschas, M.O. Vieitez, E.J. van Duijn, O. Reitebuch, W. van de Water and W. Ubachs, Appl. Opt. 49, 4217 (2010). doi:10.1364/AO.49.004217.
- B. Witschas, Appl. Opt. 50, 267 (2011). doi:10.1364/AO.50.000267.
- Z. Gu and W. Ubachs, J. Chem. Phys. 141, 104320 (2014). doi:10.1063/1.4895130.
- J. Shang, T. Wu, H. Wang, C. Yang, C. Ye, R. Hu, J. Tao and X. He, IEEE Access 7, 136439 (2019). doi:10.1109/Access.6287639.
- L. Brillouin, Ann. d. Phys. (Paris) 17, 88 (1922). doi:10.1051/anphys/192209170088.
- L.I. Mandelstam, Zh. Russ. Fiz-Khim. Ova. 58, 381 (1926).
- I.L. Fabelinskii, Molecular Scattering of Light (Springer Science, Heidelberg, 2012).
- C.S. Wang Chang and G.E. Uhlenbeck, Tech. Rep. Res. Rep. No. CM-681, Univ. of Michigan, 1951.
- X. Pan, Ph. D. dissertation, Princeton University, 2003.
- Z. Gu, Ph. D. dissertation, Vrije Universiteit Amsterdam, 2015.
- Y. Wang, Ph. D. dissertation, Vrije Universiteit Amsterdam, 2019.
- J.W. Strutt (Lord Rayleigh), Philos. Mag. 47, 375 (1899). doi:10.1080/14786449908621276.
- M. Sneep and W. Ubachs, J. Quant. Spectrosc. Ra. 92, 293 (2005). doi:10.1016/j.jqsrt.2004.07.025.
- F.M. White and I. Corfield, Viscous Fluid Flow (McGraw-Hill, New York, 2006).
- X. Pan, M.N. Shneider and R.B. Miles, Phys. Rev. A 71, 045801 (2005). doi:10.1103/PhysRevA.71.045801.
- N.A. Clark, G.R. Mellman and T.J. Greytak, Phys. Rev. Lett. 29, 150 (1972). doi:10.1103/PhysRevLett.29.150.
- B. Taylor, D. Kessler and E. Oran, in 24th International Colloquium on the Dynamics of Explosive and Reactive Systems, Taipei Taiwan, 2013.
- E.W. Lemmon and R.T. Jacobsen, Int. J. Thermophys. 25, 21 (2004). doi:10.1023/B:IJOT.0000022327.04529.f3.
- R. Montgomery, J. Meteorology 4, 193 (1947). doi:10.1175/1520-0469(1947)004<0193:VATCOA>2.0.CO;2.
- M. Leontovich, J. Phys. USSR 4, 499 (1941). https://mathscinet.ams.org/mathscinet-getitem?mr=4808
- T.D. Rossing, Springer Handbook of Acoustics (Springer, New York, 2014).
- C.Y. She, G.C. Herring, H. Moosmüller and S.A. Lee, Phys. Rev. Lett. 51, 1648 (1983). doi:10.1103/PhysRevLett.51.1648.
- J.H. Grinstead and P.F. Barker, Phys. Rev. Lett. 85, 1222 (2000). doi:10.1103/PhysRevLett.85.1222.
- X. Pan, M.N. Shneider and R.B. Miles, Phys. Rev. Lett. 89, 183001 (2002). doi:10.1103/PhysRevLett.89.183001.
- X. Pan, P.F. Barker, A. Meschanov, J.H. Grinstead, M.N. Shneider and R.B. Miles, Opt. Lett. 27, 161 (2002). doi:10.1364/OL.27.000161.
- X. Pan, M.N. Shneider and R.B. Miles, Phys. Rev. A 69, 033814 (2004). doi:10.1103/PhysRevA.69.033814.
- M.O. Vieitez, E.J. van Duijn, W. Ubachs, B. Witschas, A. Meijer, A.S. de Wijn, N.J. Dam and W. van de Water, Phys. Rev. A 82, 043836 (2010). doi:10.1103/PhysRevA.82.043836.
- A.S. Meijer, A.S. de Wijn, M.F.E. Peters, N.J. Dam and W. van de Water, J. Chem. Phys. 133, 164315 (2010). doi:10.1063/1.3491513.
- A. Manteghi, N.J. Dam, A.S. Meijer, A.S. de Wijn and W. van de Water, Phys. Rev. Lett. 107, 173903 (2011). doi:10.1103/PhysRevLett.107.173903.
- J. Graul and T. Lilly, Opt. Express 22, 20117 (2014). doi:10.1364/OE.22.020117.
- Z. Gu, M.O. Vieitez, E.J. van Duijn and W. Ubachs, Rev. Scient. Instrum. 83, 053112 (2012). doi:10.1063/1.4728575.
- Z. Gu, B. Witschas, W. van de Water and W. Ubachs, Appl. Opt. 52, 4640 (2013). doi:10.1364/AO.52.004640.
- J.E. Fookson, W.S. Gornall and H.D. Cohen, J. Chem. Phys. 65, 350 (1976). doi:10.1063/1.432774.
- V. Ghaem-Maghami and A.D. May, Phys. Rev. A 22, 692 (1980). doi:10.1103/PhysRevA.22.692.
- V. Ghaem-Maghami and A.D. May, Phys. Rev. A 22, 698 (1980). doi:10.1103/PhysRevA.22.698.