Raman scattering

Raman scattering is inelastic scattering of light in which the outgoing photon has a different frequency because matter gains or loses internal energy. Molecular vibration is a particularly useful recipient: a visible photon can reveal an internal energy spacing normally associated with infrared frequencies. Rayleigh scattering supplies the much stronger, approximately unchanged-frequency reference. The distinction concerns the transferred energy, rather than merely a change of direction or the colour of the illuminating beam.

The elementary bookkeeping is conservation of energy:Here $h$ is the Planck constant, ${\nu }_0$ and ${\nu }_s$ are the incident and scattered photon frequencies, and $E_i,E_f$ are initial and final internal molecular energies. This expression neglects translational recoil and Doppler broadening; those can matter in sufficiently precise gas measurements. A positive internal energy change gives Stokes Raman scattering, and a negative one gives Anti-Stokes Raman scattering. The energy difference determines a possible line position, but it does not by itself establish that a transition is optically allowed or predict its strength.

The spontaneous molecular process developed here has a connection to the Raman transition used in atomic measurement methods: two photon interactions transfer an internal energy difference. Coherently driving that process with two applied light fields is a different experimental regime from collecting weak spontaneous Raman scattering beside an illuminating beam.

Chandrasekhara Venkata Raman and Kariamanikkam Srinivasa Krishnan jointly reported the effect in 1928. Their letter describes light from the Sun concentrated by an $18\mathrm{cm}$-aperture, $230\mathrm{cm}$-focal length telescope objective and a second optical lens onto a liquid purified by repeated vacuum distillation. With blue-violet and yellow-green optical filters together in the incident beam, the visible track disappeared; moving the yellow-green optical filter between liquid and observer restored it. The complementary arrangement rejected unchanged-colour scattering of light while admitting newly generated light. Persistence after purification and strong wave polarization in many liquids argued against ordinary fluorescent contamination.^[ref]

In his 1930 lecture, Chandrasekhara Venkata Raman describes a cobalt-glass optical filter selecting indigo-violet light from the Sun, with a separated blue-green band in the scattered optical spectrum. Substituting a quartz mercury arc and condenser gave sharp new spectral lines. Small Hilger quartz spectrographs recorded them photographically; Kariamanikkam Srinivasa Krishnan obtained more precise records with a larger instrument and established violet-displaced lines. The optical prism disperses wavelengths for spatial readout; it does not create the Raman shift. The diagram gives functional optical paths rather than a reconstruction of a specific Hilger instrument.^[ref]

Their May 1928 report strengthens the comparison: filtered mercury illumination of toluene produced extra spectral lines. Removing the $4358\text{A˚}$ incident line removed its associated longer-wavelength group; a quinine filter removing the $4047$, $4078$ and $4109\text{A˚}$ group removed a different shifted group while preserving the one associated with $4358\text{A˚}$. These are actual selective-excitation controls, rather than inferred colours.^[ref]

The inference has limits. A complementary-filter test alone establishes newly admitted wavelengths, not their precise spacing or a unique microscopic mechanism. Fluorescence can also survive purification when intrinsic to the liquid, and fluorescence can be polarized. Useful distinctions therefore combine resolved spectral lines, excitation dependence, wave polarization and sample preparation. A proposed modern extension is to change a single incident frequency: a molecular Raman shift should follow the same internal energy gap, whereas ordinary relaxed molecular fluorescence usually retains approximately its emission-band positions. This is a diagnostic under stated conditions, not a universal theorem about every fluorescent material or resonant scattering process. Raman scattering need not involve a populated, relaxed electronic excited state as ordinary molecular fluorescence does; neither a redder photon nor the word Stokes identifies the mechanism by itself.

The 1930 Nobel Prize in Physics went to Chandrasekhara Venkata Raman alone for his work on scattering of light and the discovery of the effect. That award allocation should be distinguished from the joint experimental report with Kariamanikkam Srinivasa Krishnan.^[ref]

The signed Raman shift subtracts scattered from incident spectroscopic wavenumber:Here $c$ is the speed of light, and both wavelengths are in vacuum. Positive shifts denote Stokes Raman scattering; negative shifts denote Anti-Stokes Raman scattering. Some plots reverse the horizontal direction or report separate unsigned magnitudes, so check the convention. Spectroscopic wavenumber has no $2\pi$ factor, unlike the common angular wavenumber $k=2\pi /\lambda$. A shift of $1000\mathrm{c}{\mathrm{m}}^{-1}$ means $100000{\mathrm{m}}^{-1}$, not a frequency of $1000\mathrm{Hz}$.

For a one-quantum change of molecular vibration with ordinary frequency ${\nu }_v$, the magnitude is ${\tilde{\nu }}_v={\nu }_v/c$. This makes Raman spectroscopy attractive: the same molecular energy spacing can be read using different visible excitation wavelengths. Equal and opposite Raman shifts are not equal and opposite changes measured in nanometres, because wavelength is the reciprocal of spectroscopic wavenumber.

A Raman shift as a probe of bond stiffness interprets a measured vibrational frequency through a mechanical model of the assigned mode. For one stretching coordinate modeled as a harmonic oscillator,where $k_{\mathrm{eff}}$ is an effective spring constant and ${\mu }_{\mathrm{eff}}$ an effective mass. For an isolated diatomic stretch the latter is the reduced mass. A collective normal mode can involve many atoms, so its effective stiffness is not automatically the force constant of one selected bond.

An original same-mode example shifts from $1500$ to $1530\mathrm{c}{\mathrm{m}}^{-1}$. If the effective mass is unchanged, the harmonic model givesor a $4.04\%$ increase in effective stiffness. The small-change version isA ten-percent mass increase at fixed stiffness would instead reduce the frequency by $1/\sqrt{1.10}\approx 0.9535$. This reduced-mass ambiguity in a Raman stiffness inference is why a frequency change alone does not uniquely measure a stiffness change.

Mechanical strain can shift Raman bands by changing the local vibrational dynamics. HORIBA's application documentation illustrates stress mapping in microfabricated silicon, including shallow-trench regions, and describes a correlation between applied stress and peak shift.^[ref] That calibrated material response is more specific than treating every shift as a direct scalar bond-stiffness measurement.

A motivating experiment would track the same assigned band while varying a known load, holding temperature, excitation and polarization conditions fixed and measuring the wavelength reference. Heating, isotopic composition, mode mixing and structural changes can also alter frequencies. The Raman shift gives an energy spacing; mechanical interpretation additionally needs the mode assignment, mass model and appropriate stress/strain calibration.

Reduced-mass ambiguity in a Raman stiffness inference means that a measured frequency ratio does not determine a stiffness ratio without a mass model. In the assigned harmonic oscillator approximation used in Raman shift as a probe of bond stiffness, the vibrational wavenumber obeys $\tilde{\nu }\propto \sqrt{k_{\mathrm{eff}}/{\mu }_{\mathrm{eff}}}$. ThusThe reduced mass of an isolated two-body oscillator is $m_1{m}_2/(m_1+m_2)$. An effective mass for a collective normal mode additionally depends on which atoms move and their relative displacements. Keeping the same peak label does not establish that this effective mass stayed fixed.

For an original example, a band moves from $1000$ to $980\mathrm{c}{\mathrm{m}}^{-1}$, giving $q=0.98$. With unchanged effective mass, the inferred stiffness ratio is $0.9604$, a $3.96\%$ reduction. With an effective-mass ratio $1.06$, the same observed shift instead gives $k_2/k_1=1.018024$: stiffness has increased by about $1.80\%$. The two interpretations satisfy the same measured frequency ratio.

This ambiguity matters when comparing compositions, isotopes or changing mode assignments. Identifiability of stiffness from isotope-resolved Raman shifts requires additional constraints on the mass model and assignment; a frequency shift alone supplies one equation for two ratios. A controlled experiment can compare a known composition under load, stabilize temperature and optical conditions, and measure the spectral reference while checking that the same normal mode is followed.

HORIBA's silicon microstructure application maps stress using a material-specific relationship between stress and the Raman peak shift.^[ref] That experiment motivates an independently calibrated response. It does not make every redshift a universal direct measurement of a reduced spring constant.

Identifiability of stiffness from isotope-resolved Raman shifts asks which stiffness parameters become determined once the moving masses are independently constrained. In the assigned harmonic oscillator model of reduced-mass ambiguity in a Raman stiffness inference,Here $c$ is the speed of light, ${\tilde{\nu }}_i$ the vibrational wavenumber and ${\mu }_i$ the mode's known effective mass. Known masses and identified modes make each effective stiffness identifiable within this model. Unknown masses still permit the simultaneous rescaling $(k_i,{\mu }_i)\mapsto (a{k}_i,a{\mu }_i)$ without changing a frequency; knowing only mass ratios determines only stiffness ratios.

Controlled isotope substitution is useful because it changes nuclear mass while approximately preserving the interatomic force model. In a pure-carbon mode for which every moving mass changes by the same ratio, unchanged stiffness predictsThe mass numbers here are an explicit approximation. They replace a full normal mode calculation only for the stipulated uniform mass scaling.

For an original numerical test, a chosen carbon-12 reference band at $1580\mathrm{c}{\mathrm{m}}^{-1}$ predicts about $1518.01\mathrm{c}{\mathrm{m}}^{-1}$ after uniform carbon-13 replacement with unchanged stiffness. If a hypothetical carbon-13 band is instead $1520\mathrm{c}{\mathrm{m}}^{-1}$, its isotope-corrected Raman stiffness ratio isThe residual inferred change is about $0.262\%$, not the several-percent softening one would infer by ignoring mass. This is a chosen calculation, not a quoted experimental peak pair.

The primary graphene study measured several Raman bands at different carbon-13 concentrations and laser excitation energies; its measured mass-dependent frequencies followed the harmonic-oscillator relation.^[ref] Such a controlled series tests the mass model. It does not by itself rule out changes in strain, temperature or mode mixing in a different sample. Uncertainty in isotope composition, spectral calibration and mode assignment must accompany any small residual stiffness inference.

An isotope-corrected Raman stiffness ratio removes the stipulated mass contribution before interpreting a vibrational frequency change as a stiffness change. For corresponding harmonic modes,The identifiability of stiffness from isotope-resolved Raman shifts article explains why the effective masses and mode assignment are required. The ratio does not by itself validate those assumptions.

For an original calculation, take hypothetical pure-carbon bands at $1600\mathrm{c}{\mathrm{m}}^{-1}$ for carbon-12 and $1530\mathrm{c}{\mathrm{m}}^{-1}$ for carbon-13. Uniform isotope replacement with the mass-number approximation gives ${\mu }_2/{\mu }_1=13/12$ and $K\approx 0.99062$. The inferred stiffness reduction is about $0.938\%$, despite the $4.375\%$ decrease in wavenumber. These chosen values are not a published experimental pair.

Uncertainty in the ratio follows fromWith exact masses in the model and independent wavenumber standard uncertainties $0.5\mathrm{c}{\mathrm{m}}^{-1}$, the relative standard uncertainty of $K$ is approximately $0.0904\%$. A single common multiplicative spectral-scale error cancels exactly in the frequency ratio. A common additive offset does not: its first-order effect is $2\delta (1/1530-1/1600)$, giving only about $0.00286\%$ relative uncertainty for a $0.5\mathrm{c}{\mathrm{m}}^{-1}$ offset uncertainty. The relative calibration error in an isotope-corrected stiffness ratio therefore depends on the spectral calibration model and correlations, not just each peak's marginal error bar.^[ref]

The primary graphene experiment varied carbon-13 concentration and laser excitation energy and observed Raman-frequency dependence following the harmonic-oscillator mass relation.^[ref] A residual stiffness comparison additionally needs controlled strain, temperature and composition. Mode mixing or uncertain isotope distribution can dominate the small residual even when fitted peak centres are precise.

Relative calibration error in an isotope-corrected stiffness ratio is the fractional error remaining when spectral calibration affects the two compared bands differently. Write the model of isotope-corrected Raman stiffness ratio as $K=m({\nu }_2/{\nu }_1{)}^2$, where $m={\mu }_2/{\mu }_1$ and the positive ${\nu }_i$ denote vibrational wavenumbers. For small errors,The mass correction does not remove differential spectral calibration error.

A useful calibration model is $\delta {\nu }_i=(g_c+g_i){\nu }_i+b_c+b_i$, with common and separate gain errors $g_c,g_i$ and offsets $b_c,b_i$. Its spectral contribution becomeswith the mass term added if uncertain. A single multiplicative scale cancels exactly when no additive error is present. Separate gains do not: comparing two spectra calibrated independently can lose that cancellation.

For an original example, take ${\nu }_1=1000\mathrm{c}{\mathrm{m}}^{-1}$ and ${\nu }_2=900\mathrm{c}{\mathrm{m}}^{-1}$. A differential gain $g_2-g_1=0.001$ creates about $0.20\%$ relative error in $K$. A common offset $b_c=1\mathrm{c}{\mathrm{m}}^{-1}$ instead contributes about $0.0222\%$. If the separate gains each have independent $0.1\%$ standard uncertainty, their combined contribution to the relative standard uncertainty of $K$ is $2\sqrt{2}\times 0.1\%\approx 0.2828\%$.

The separate spectral-calibration gains in isotope stiffness inference must therefore be distinguished from a shared reference error. Covariance propagation keeps precisely those shared and separate influences.^[ref] Recording both samples with the same validated spectral scale can help, but drift or changed dispersion can still produce differential errors.

The primary graphene isotope experiment measured Raman bands over isotope concentrations and laser excitation energies and found the harmonic-oscillator mass dependence.^[ref] Small residual stiffness changes require the additional calibration and mode-model checks above; the numerical errors here are chosen examples, not quoted experimental uncertainties.

Separate spectral-calibration gains in isotope stiffness inference arise when the two Raman spectra retain different multiplicative calibration errors. Let the corrected wavenumbers satisfy ${\hat{\nu }}_i=(1+g_i){\nu }_i$. For $K=({\mu }_2/{\mu }_1)({\nu }_2/{\nu }_1{)}^2$, small residual gains giveThe relative calibration error in an isotope-corrected stiffness ratio article treats the accompanying offset terms. Here the important quantity is the gain difference, not either gain alone.

If the two residual gains have standard uncertainties $u_1,u_2$ and correlation $\rho$, their contribution isThis follows by applying the full covariance propagation rule.^[ref] A shared scale component cancels; independent session errors do not. The covariance of differential spectral-calibration gains therefore determines whether combining calibrations improves the ratio.

For an original example, let $u_1=u_2=0.2\%$. Independent calibrations give relative uncertainty about $0.566\%$ in $K$. A justified correlation $\rho =0.9$ reduces it to about $0.179\%$; a correlation $-0.9$ instead gives about $0.780\%$. These are different calibration models, not alternative ways of adding the same independent errors.

A practical comparison can measure the same known reference line during each sample session. If raw readings are $y_i=c_i{\nu }_i$ and reference readings $r_i=c_i{\nu }_{\mathrm{ref}}$, a purely multiplicative model givesBoth session gains cancel, but reference-reading noise, offsets and drift between reference and sample readings remain. Repeated readings do not automatically reveal an unknown common bias.

The primary graphene isotope experiment compared Raman bands over isotope concentrations and excitation energies, finding the harmonic-oscillator mass dependence.^[ref] The calibration model above concerns precision of the inferred residual stiffness ratio; it does not validate mode assignment or the effective-mass approximation.

The covariance of differential spectral-calibration gains describes the joint errors of gain differences used in several spectral ratios. Let the calibration-error vector be $\mathbf{g}$ with covariance ${\Sigma }_g$, and form differences $\mathbf{d}=B\mathbf{g}$. ThenThis covariance identity is exact for the linear differences; no normal-distribution assumption is required. The first-order relation between gain differences and isotope stiffness ratios is given in separate spectral-calibration gains in isotope stiffness inference.^[ref]

For an original shared-reference model, take three gain errors $g_i=c+e_i$, where $c$ is a common scale error and the $e_i$ are independent residual errors with variances $s_i^2$. Two sample spectra are compared with the same reference spectrum:The common scale cancels, but the residual reference error remains in both differences:Cancellation of one common influence therefore does not establish independent comparisons.

Choose all three residual standard uncertainties as $0.1\%$. Then ${\Sigma }_d=10^{-6}\left(\begin{array}{cc}2& 1\\ 1& 2\end{array}\right)$, with correlation $1/2$ between the differences. Even a much larger common scale uncertainty cancels in this stipulated model. Each inferred log stiffness ratio has error approximately twice its gain difference. Averaging the two log-ratio estimates gives variance $6\times 10^{-6}$ and standard uncertainty about $0.245\%$. Incorrectly treating them as independent would give only $0.200\%$.

This shared-reference floor in averaged isotope stiffness ratios matters when many sample spectra reuse one calibrated reference. Independent extra sample readings reduce their own errors, while a reference component shared by all comparisons remains. Assess covariance from the actual reference reuse, calibration drift and repeated acquisition protocol. The matrix calculation concerns multiplicative calibration errors; uncertain peak assignment, isotope composition and additive spectral offsets need their own contributions.

A shared-reference floor in averaged isotope stiffness ratios is the calibration uncertainty that remains when many sample measurements reuse one uncertain reference. In the first-order gain model of covariance of differential spectral-calibration gains, let the error in log stiffness-ratio estimate $i$ be $2(e_i-e_r)$. Assume independent sample errors with variance $s^2$ and one independent reference error with variance $s_r^2$. Averaging $m$ estimates givesThe sample contribution decreases with replication, while the shared reference contribution is unchanged. This is the covariance propagation rule applied to the actual reuse of the reference, rather than treating every ratio as independent.^[ref]

For an original example, choose sample calibration standard uncertainty $s=0.2\%$ and reference uncertainty $s_r=0.1\%$. One ratio has approximately $0.447\%$ standard uncertainty on the log-ratio scale. Sixteen sample measurements reduce it to about $0.224\%$, but the limiting floor is $2s_r=0.200\%$. Reaching within ten percent of that floor requiresso $m\ge 20$ for these chosen numbers. Replication beyond that point has progressively less effect unless the reference uncertainty is reduced.

If a new independent reference is measured for each sample, the model changes: both error contributions can average down, giving $4(s^2+s_r^2)/m$ under the corresponding independence assumptions. Alternatively, averaging several independent reference measurements reduces their random component, while a systematic component shared by all references still remains. A reference-versus-sample allocation in repeated isotope-ratio measurements must distinguish those resources and common influences.

The floor is conditional on the measurement model. Additive spectral calibration errors, reference drift, isotope-composition uncertainty and mode-assignment error are separate contributions. A useful acquisition protocol records which sample spectra share which reference observations; numerical averaging alone cannot remove an uncertainty component reused unchanged in every estimate.

Reference-versus-sample allocation in repeated isotope-ratio measurements distributes a finite measurement budget between sample readings and the reference used to calibrate them. Suppose independent sample and reference readings have first-order gain-error variances $s^2$ and $r^2$. Averaging $n_s$ sample readings and $n_r$ reference readings gives log stiffness-ratio variancewhere $s_c^2$ is the variance of a differential calibration error shared by all ratio estimates. The shared-reference floor in averaged isotope stiffness ratios explains why that last component does not average down. Covariance propagation must reflect which observations and calibration influences are actually shared.^[ref]

Let one sample reading cost $c_s$ budget units and one reference reading cost $c_r$, with $c_s{n}_s+c_r{n}_r=T$. Treating counts continuously for the initial optimization, a constraint multiplier givesThe counts are proportional to standard uncertainty divided by the square root of cost, rather than simply to each variance. A budget-independent common component changes the attainable floor but not this optimum ratio under the stated model.

For an original example, take equal costs, $s=0.3\%$, $r=0.1\%$ and 40 total readings. The optimum is 30 sample and ten reference readings. Without a common component, the predicted log-ratio standard uncertainty is about $0.1265\%$, compared with $0.1414\%$ for twenty of each. A common standard uncertainty $s_c=0.05\%$ raises those values to about $0.1612\%$ and $0.1732\%$, respectively.

An integer allocation near a continuous isotope-reference optimum must compare feasible neighboring count choices and any acquisition-order constraints. More reference readings help only to the extent that their errors are independent; repeatedly measuring an unchanged unknown calibration bias does not create independent references. If reference drift depends on elapsed time, a single constant cost and common-variance term may be inadequate.

In practice, record which reference observations calibrate which sample spectra and account for any correlated peak fitting, isotope-composition uncertainty and spectral offsets. The variance and cost model above is a chosen allocation problem, not an experimentally established precision specification for an isotope measurement system.

An integer allocation near a continuous isotope-reference optimum converts a real-valued budget solution into feasible counts of sample and reference observations. For the model in reference-versus-sample allocation in repeated isotope-ratio measurements, minimizewith positive integer counts. The shared component $s_c^2$ does not affect the choice when it is genuinely independent of the allocation. Counting new reference readings does not make a shared calibration error independent.^[ref]

For equal unit costs and a budget of $N$ readings, put $n_r=N-n_s$. With $s,r>0$, the count-dependent objective $f(n)=s^2/n+r^2/(N-n)$ hasIts continuous minimum is $n_{\ast }=Ns/(s+r)$. Therefore an integer minimum lies among the feasible floor and ceiling of $n_{\ast }$, clipped to $1,\dots ,N-1$. Comparing those candidates is sufficient because the function decreases before its minimum and increases after it. Rounding both continuous counts independently can violate the budget.

For an original example, let $N=17$, $s=0.003$ and $r=0.001$. The continuous sample count is12.75. Counts $(n_s,n_r)=(12,5)$ give variance $3.8000\times 10^{-6}$ before a common component, while $(13,4)$ give approximately $3.7692\times 10^{-6}$. The latter is preferable, with log-ratio standard uncertainty about0.1941%. The nearby $(14,3)$ gives approximately $3.9048\times 10^{-6}$.

Unequal costs introduce a stepped feasible boundary. One exact finite search enumerates positive $n_r$ and, for each, useswhen that value is positive. For fixed $n_r$, using the largest feasible $n_s$ improves this stated objective. For an original $T=20$, $c_s=3$, $c_r=2$, $s=3$ and $r=1$ in consistent planning units, the best pair is $(5,2)$, consuming19 units. Its count-dependent expression $s^2/n_s+r^2/n_r$ is2.3; $(6,1)$ and $(4,4)$ both give2.5. An unused budget unit can be a legitimate integer remainder.

Acquisition ordering, paired exposures, minimum replicate counts or drift constraints can restrict the feasible set further. The continuous optimum remains a useful guide, but only comparisons among genuinely feasible integer designs establish the chosen discrete optimum.

In Stokes Raman scattering, matter gains internal energy and the photon loses the same amount. For a quantum harmonic oscillator representing one molecular vibration, $E_v=h{\nu }_v(v+1/2)$ with $v=0,1,\dots$. The transition $v\to v+1$ gives ${\nu }_s={\nu }_0-{\nu }_v$, requiring ${\nu }_0>{\nu }_v$ for an outgoing photon. Even a molecule initially in its vibrational ground state can undergo this process. The shared name does not turn it into the absorption/emission-band Stokes shift used for fluorescence.

In Anti-Stokes Raman scattering, matter loses internal energy and the outgoing photon gains it. The harmonic transition $v\to v-1$ requires $v\ge 1$, so an initially excited vibrational population must be available. At low temperature, high-frequency vibrational excitations are scarce and the corresponding anti-Stokes line is weak. The extra photon energy comes from the molecule; there is no conflict with conservation of energy.

The simplest off-resonant explanation is a vibrating polarizability. Let $Q(t)=Q_0\cos (2\pi {\nu }_{v}t)$ be a small normal mode coordinate of molecular vibration, and use one scalar component of polarizability, $\alpha (Q)\simeq {\alpha }_0+{\alpha }_0^{\prime }Q$. For incident electric field $E(t)=E_0\cos (2\pi {\nu }_{0}t)$, the induced electric dipole moment isThe product-to-sum formula produces the two displaced frequencies. A nonzero derivative ${\alpha }_0^{\prime }$ gives vibrational Raman scattering in this linear model. Real polarizability is generally a tensor, and the selected components depend on incident and observed wave polarization. This condition differs from vibrational infrared activity, which requires a changing molecular electric dipole moment without an applied optical electric field. A vibration can therefore appear in one measurement and be absent in the other.

This classical modulation argument explains line locations but supplies equal sideband wave amplitudes; it does not explain their thermal asymmetry or the nonzero Stokes Raman scattering from a vibrational ground state. Those require the quantum harmonic oscillator population and transition strengths. The far-off-resonance drawing uses dashed virtual intermediate levels: they represent terms in a scattering probability amplitude, not a long-lived occupied electronic energy eigenstate. Near an electronic resonance, a more detailed treatment of the intermediate response is needed.

Assume one harmonic normal mode at thermal equilibrium, weak spontaneous excitation far from electronic resonances, the same collection solid angle and wave polarization, and matched response after instrumental correction. The Boltzmann distribution givesHere $k_B$ is the Boltzmann constant, $T>0$ is temperature, and $P_v$ is the probability of vibrational quantum number $v$. Thus $P_1/P_0=r$. Linear coupling to $Q$ has squared oscillator matrix elements proportional to $v+1$ for Stokes Raman scattering and $v$ for Anti-Stokes Raman scattering. Summing over initial levels gives strengths proportional to $\overline{v}+1$ and $\overline{v}$ respectively; their ratio is $r$, even when levels above $v=1$ are occupied.

Let $I_{\mathrm{S}}$ and $I_{\mathrm{AS}}$ denote line-integrated radiant intensities, meaning power per collection solid angle. The off-resonant electric dipole radiation law adds the fourth power of emitted frequency:The photon rate divides each line's power by its own $h{\nu }_s$. Consequently, the ideal response-corrected anti-Stokes/Stokes count-rate ratio has the third power of the same frequency ratio. Peak heights alone do not give line-integrated radiant intensities if linewidths differ. Uncorrected detector counts can also include unequal detector response and unequal optical filter transmission. Resonant enhancement, overlapping spectral lines, nonthermal populations and local optical heating can invalidate a temperature inferred from this simple formula.

Choose a vibrational Raman shift magnitude ${\tilde{\nu }}_v=1000\mathrm{c}{\mathrm{m}}^{-1}$, excitation vacuum wavelength ${\lambda }_0=532\mathrm{nm}$ and temperature $T=300\mathrm{K}$. These are original illustrative inputs, not a measured 1928 optical spectrum or a claim about a particular liquid. They giveTake reciprocal lengths in consistent units before converting to nanometres. At $300\mathrm{K}$, $r\simeq 0.00826337$ and $\overline{v}\simeq 0.00833222$. The predicted line-integrated radiant intensity ratio is $0.0126522$, whereas the ideal photon count-rate ratio is $0.0113740$. Thus the radiometric ratio is about $1.27$ percent, rather than the $0.826$ percent obtained by using the Boltzmann factor alone.

Raman spectroscopy measures shifted scattering of light to infer internal energy spacings and the molecular response that makes them visible. In a modern illustrative arrangement, narrowband laser illumination replaces the historical arc, an optical filter rejects the strong excitation line, and a spectrograph resolves the weak shifted spectral lines. The illumination, rejection and spectral dispersion are separate jobs. Increasing rejection can conceal small Raman shifts close to the excitation line.

A proposed measurement should record the excitation wavelength, a dark signal, a blank container and the liquid, then compare calibrated line positions and integrated signals under the same collection conditions. Measure optical filter transmission and detector response before using the Thermal Raman line ratio for temperature. Repeat at lower incident power to look for heating or altered sample response. These are suggested checks, not an assertion that the 1928 apparatus implemented modern calibration procedures.

Its appeal is the conversion of internal molecular energy differences into visible spectral lines: changing the excitation can move the whole pattern while preserving characteristic Raman shifts. Combined with vibrational infrared activity, Raman vibrational activity supplies complementary information about molecular vibration. A line pattern can constrain a material's identity or environment, but overlapping bands, fluorescence backgrounds and the finite spectral resolution element prevent an unqualified identification from one peak alone.

A Raman spectrum records the strength of Raman scattering as a function of scattered optical frequency or, usually, Raman shift relative to the excitation. For vacuum wavelengths in nanometres, a common wavenumber convention isHere the numerical wavelengths are in nanometres. Stokes Raman scattering has positive shift under this convention; Anti-Stokes Raman scattering has negative shift. The strong unshifted Rayleigh scattering line is not a vibrational Raman band.

For an original chosen excitation ${\lambda }_0=500\mathrm{nm}$ and a $1500\mathrm{c}{\mathrm{m}}^{-1}$ vibration, the Stokes wavenumber is $20000-1500=18500\mathrm{c}{\mathrm{m}}^{-1}$, giving ${\lambda }_s\approx 540.54\mathrm{nm}$. The corresponding anti-Stokes wavelength is $10^7/21500\approx 465.12\mathrm{nm}$. Changing the excitation changes these optical wavelengths while the same assumed transition retains its Raman shift. The One-mode Raman spectrum example gives an existing level diagram and a separate numerical spectrum.

The vertical quantity must be stated too: detector counts per spectral bin, calibrated scattered power density and a stimulated gain ratio are different observables. Femtosecond stimulated Raman spectroscopy can plot the Raman pump-on pump-off gain ratio rather than spontaneous photon counts. Its primary apparatus disperses probe and reference beams onto an array and calibrates shifts with known solvent lines.^[ref]

Spectral-density conversion between wavelength and Raman shift also changes the density's units: rebadging the horizontal axis alone need not preserve a spectral integral. Background, excitation rejection, instrument response and spectral resolution affect observed band heights and widths. An integrated band area is not determined by peak height alone. A chemical inference should compare calibrated positions and responses under specified conditions, rather than identify a material from one unqualified peak.

Spectral-density conversion between wavelength and Raman shift preserves the integrated physical quantity while changing its horizontal coordinate. Fix the excitation wavelength ${\lambda }_0$, and let $s$ be the Raman shift. Using numerical vacuum wavelengths in nanometres and shifts in $\mathrm{c}{\mathrm{m}}^{-1}$,Let $D_{\lambda }$ denote detected counts per nanometre and $D_s$ counts per unit Raman shift. Conservation of counts under a change of variables givesThe conversion factor carries the corresponding density-unit change. The same argument applies to calibrated power densities, provided both sides describe the same quantity.

For an original narrow bin centred at $600\mathrm{nm}$, a count density of $1200$ per nanometre becomes $43.2$ per $\mathrm{c}{\mathrm{m}}^{-1}$. A one-nanometre bin there spans approximately $27.78\mathrm{c}{\mathrm{m}}^{-1}$, so its integrated count remains about $1200$. At $800\mathrm{nm}$, the same wavelength density becomes $76.8$ per $\mathrm{c}{\mathrm{m}}^{-1}$. Relabeling bin positions without adjusting the density would conceal this coordinate dependence.

For finite bins, use their exact transformed edges and preserve each bin's integral. Raman spectral rebinning with conserved integrated counts distinguishes this operation from interpolating peak heights onto a new grid. Instrument response correction is another step; the Jacobian alone does not calibrate detector sensitivity or remove background.

The primary FSRS apparatus calibrated Raman shifts using known solvent lines and dispersed the probe spectrum onto an array.^[ref] Its dimensionless Raman pump-on pump-off gain ratio is a different observable: a ratio of matched spectral densities can cancel their common coordinate factor. It should not be treated as a count density merely because it is plotted against Raman shift.

Raman spectral rebinning with conserved integrated counts redistributes finite-bin measurements onto a new spectral grid while preserving their summed quantity over the same covered range. It is different from relabeling bin centres or interpolating peak heights. The parent spectral-density conversion between wavelength and Raman shift derives the continuous Jacobian; finite bins require a stated assumption about unresolved variation inside each bin.

For an original piecewise-constant wavelength-density model, let old bin $I_i$ have width $\Delta {\lambda }_i$ and count $N_i$. Invert the new Raman-shift edges to wavelength intervals $J_j$. Define the overlap matrixWhen the new bins partition the old covered range, ${\sum }_j{A}_{ji}=1$, hence ${\sum }_j{N}_j^{\prime }={\sum }_i{N}_i$. Clipping the new range instead excludes the corresponding counts and must be reported.

For original old wavelength bins $[600,602]$ and $[602,604]\mathrm{nm}$ with counts $100$ and $300$, choose new edges corresponding to wavelengths $600$, $601$, $603$ and $604\mathrm{nm}$. ThenThe total remains $400$. A different within-bin model can distribute that total differently; the original measurements do not resolve every sub-bin feature.

Covariance of count-conserving rebinned Raman channels transforms as $C^{\prime }=AC{A}^T$ for this fixed linear operation. Fractional rebinned estimates are therefore not automatically independent Poisson photon counts. Splitting one measured bin among neighbours creates shared uncertainty; finer output pixels do not create additional spectral resolution.

The primary FSRS instrument records a dispersed spectrum on an array and calibrates its Raman-shift axis.^[ref] The overlap algorithm above is an original illustrative reconstruction, not a claim about that instrument's published processing. Apply radiometric response and background corrections consistently with the quantity whose integral is being conserved.

The covariance of count-conserving rebinned Raman channels records the shared noise created when one measured channel contributes to several output channels. For a fixed deterministic rebinning matrix $A$ and an input count vector $\mathbf{N}$ with covariance matrix $C$, the output is ${\mathbf{N}}^{\prime }=A\mathbf{N}$ andThis identity needs finite second moments, not a Gaussian approximation. Raman spectral rebinning with conserved integrated counts describes the count-preserving map; retaining its full covariance matrix is necessary for uncertainties of sums and fitted spectral features.^[ref]

Consider independent input counts with means $100$ and $300$ and Poisson distribution variances equal to those means. Split them deterministically usingThe three output means are $50,200,150$. Their total has variance $1^T{C}^{\prime }1=400$, exactly the original total variance, and standard deviation $20$. Keeping only diagonal entries would give variance $200$ and standard deviation about $14.14$, understating the uncertainty of the conserved total.

The rank deficiency from deterministic spectral upsampling is visible without matrix inversion: the vector $(1,-1,1{)}^T$ has zero output contrast for every possible input. Three displayed channels therefore contain at most two independent degrees of freedom in this model. Fractional redistribution creates correlated numerical channels; it is different from physically assigning each detected photon randomly to one new channel.

The primary stimulated-Raman apparatus records spectra with array detectors and calibrates their spectral coordinates.^[ref] When those data are placed on a new spectral grid, the original detector counts remain the source of statistical information. Additional interpolation points do not supply independent photon measurements.

Rank deficiency from deterministic spectral upsampling occurs because a deterministic map cannot manufacture independent measurements. For $\mathbf{y}=A\mathbf{n}$ and input covariance matrix $C$,When the output has more channels than the input, its covariance matrix is necessarily singular. Covariance of count-conserving rebinned Raman channels gives the conserved-count example and its full propagated uncertainty.^[ref]

Using that parent's two-to-three channel map, the output always satisfies $y_1-y_2+y_3=0$. If the input readings are $20$ and $30$, the rebinned vector is $(10,25,15{)}^T$. Recovering the original readings as $n_1=2y_1$ and $n_2=2y_3$ shows that the third displayed number is not a third independent observation.

For an original fitting example, suppose both input means equal an unknown $\theta$ and their independent measurement-error variances are known to be $25$. Fitting in the original channels gives $\hat{\theta }=(20+30)/2=25$, with variance $25/2=12.5$. A correct generalized least squares fit after rebinning gives the same result and uncertainty. It must use the singular covariance on its supported range, for example through the Moore-Penrose inverse, together with the exact linear constraints.

The range restriction in singular-covariance spectral fitting matters if a proposed output model fails those constraints. A residual along $(1,-1,1{)}^T$ cannot arise from this rebinning model at all. Merely inserting a pseudoinverse into a quadratic objective would assign zero penalty to that impossible direction. Either enforce the range condition or supply a physically justified additional noise model; adding a small diagonal term solely to make the matrix invertible changes the statistical assumptions.

The motivating Raman apparatus measures detector-array spectra before any subsequent numerical grid change.^[ref] Working on the original grid often avoids this singular fitting problem while retaining the actual measurement information.

Range restriction in singular-covariance spectral fitting requires a residual to lie in the subspace where the measurement model allows random error. Let $\mathbf{y}=\mathbf{m}(\theta )+\mathbf{ϵ}$ with zero-mean error covariance $C$. Every vector $\mathbf{v}$ in the kernel of $C$ satisfiesso ${\mathbf{v}}^T\mathbf{ϵ}=0$ almost surely. Thus a permissible residual obeyswhere $C^{+}$ is the Moore-Penrose inverse. The matrix $C{C}^{+}$ is the orthogonal projector onto the range of the covariance.

Only after enforcing this condition should a generalized least squares objective use $(\mathbf{y}-\mathbf{m}{)}^T{C}^{+}(\mathbf{y}-\mathbf{m})$. The quadratic form alone assigns zero penalty to the forbidden directions and therefore cannot impose the condition. Rank deficiency from deterministic spectral upsampling gives its physical spectral-grid origin through $C=A{C}_{\mathrm{in}}{A}^T$. Covariance propagation preserves those exact dependencies.^[ref]

For an original idealized two-channel example, take $C=\mathrm{diag}(1,0)$, mean $\mathbf{m}(\theta )=(\theta ,2\theta {)}^T$ and observation $(3,4{)}^T$. The second channel has no random error in this stipulated model, so the range condition forces $\theta =2$. The remaining squared residual is one. Minimizing only the pseudoinverse quadratic would instead choose $\theta =3$ and objective zero, while predicting the impossible second-channel value six.

Real measurements need not contain a perfect channel: deterministic rebinning produces exact noise-free combinations of otherwise noisy channels. Additional physical noise in such a combination must be modelled if present. If covariance ranges change with the parameters, pseudoinverse fitting with a parameter-dependent covariance range requires enforcing the appropriate support for each candidate; a single fixed projection need not apply.

Raman detector-array spectra are recorded before numerical redistribution onto another spectral grid.^[ref] Fitting the original independent measurement coordinates often makes these support restrictions explicit and avoids treating interpolation points as new observations.

Pseudoinverse fitting with a parameter-dependent covariance range requires checking the allowed residual subspace separately for each fitted parameter. If the covariance is $C_{\theta }$ and the residual is ${\mathbf{r}}_{\theta }=\mathbf{y}-\mathbf{m}(\theta )$, the range restriction in singular-covariance spectral fitting becomesThe Moore-Penrose inverse quadratic ${\mathbf{r}}_{\theta }^T{C}_{\theta }^{+}{\mathbf{r}}_{\theta }$ alone does not enforce this condition. A projector computed at one trial parameter need not describe another trial's allowed noise.

For an original example, let ${\mathbf{u}}_{\theta }=(\cos \theta ,\sin \theta {)}^T$, zero mean and $C_{\theta }={\sigma }^2{\mathbf{u}}_{\theta }{\mathbf{u}}_{\theta }^T$. Noise lies on a line that rotates with $\theta$. With observation $\mathbf{y}=(2,0{)}^T$ and $\sigma =1$, the model at $\theta =0$ permits the observation and gives quadratic value four. At $\theta =\pi /2$, the pseudoinverse quadratic is zero, yet the observation is outside the allowed vertical line. Minimizing that quadratic without its range constraint would reward an impossible candidate.

The pinned SciPy documentation describes singular-covariance evaluation using pseudo-determinants, pseudoinverses and covariance rank; its implementation additionally applies a numerical support mask and returns zero density, or negative-infinite log density, off support.^[ref] Parameter-dependent supports also require care about which measure defines a density. Values defined on different subspaces are not automatically a common likelihood merely because each uses the same-looking formula.

A physically justified transverse noise resolves the idealized singularity. With perpendicular unit vector ${\mathbf{v}}_{\theta }$, use $C_{\theta ,ϵ}={\sigma }^2{\mathbf{u}}_{\theta }{\mathbf{u}}_{\theta }^T+{ϵ}^2{\mathbf{v}}_{\theta }{\mathbf{v}}_{\theta }^T$. For $ϵ=0.1$, the two quadratic values become four and four hundred. The determinant is the same at these orientations, so this full-rank Gaussian model correctly disfavors the transverse observation.

The transverse-noise limit of a rotating rank-one covariance model is an inference about a stated physical noise model, not permission to add arbitrary diagonal noise for numerical convenience. In spectral work, distinguish a parameter-dependent physical noise direction from a fixed deterministic rebinning matrix; the latter has a fixed output range when its input covariance is positive definite.

The transverse-noise limit of a rotating rank-one covariance model describes how an ordinary two-dimensional noise distribution concentrates onto a line as one physical noise scale tends to zero. For orthonormal vectors ${\mathbf{u}}_{\theta },{\mathbf{v}}_{\theta }$, setA zero-mean Gaussian observation has densityThe pseudoinverse fitting with a parameter-dependent covariance range article explains the range restriction at the exactly singular endpoint.

For a fixed point off the line, ${\mathbf{v}}_{\theta }^T\mathbf{y}\ne 0$, the exponential drives this density to zero as $ϵ\downarrow 0$. On the line the two-dimensional density grows like $1/ϵ$. There is no finite pointwise two-dimensional density limit: probability concentrates onto a one-dimensional Gaussian measure. A density on that line uses a different reference measure. Simply substituting zero for $ϵ$ in the full-rank density formula is invalid.

For an original angular example, let an observed vector have length $R=3$ along orientation ${\theta }_0$, with $\sigma =1$ and $ϵ=0.06$. At angular displacement $\delta =\theta -{\theta }_0$,At $\delta =0.02$ radians the density ratio is about $0.608$. The angular resolution set by transverse spectral noise has characteristic small-angle scale $ϵ/R\approx 0.02$ radians, or $1.15^{\circ }$, when the transverse noise is much smaller than the longitudinal noise.

SciPy's pinned implementation treats singular covariance using its rank and a separate numerical support check, rather than only discarding zero eigenvalues in a quadratic.^[ref] Numerical support tolerance and actual transverse measurement noise are different quantities. Choose the latter from the detector and processing model; a tiny arbitrary regularization can otherwise imply unjustified angular precision.

Angular resolution set by transverse spectral noise quantifies how sharply an observed spectral vector constrains the orientation of an anisotropic two-coordinate noise model. Use the covariance from transverse-noise limit of a rotating rank-one covariance model, with known longitudinal standard deviation $\sigma$ and transverse standard deviation $0<ϵ<\sigma$. For an observation of length $R>0$ aligned with orientation ${\theta }_0$, the Gaussian density givesThe determinant is independent of orientation, so the angular dependence comes entirely from the covariance-weighted residual.^[ref]

Define a local half-width by a chosen log-density drop $c>0$. For $2c\le K$, the lobe around ${\theta }_0$ has half-widthWhen this angle is small and $ϵ\ll \sigma$, it is approximately $\sqrt{2c}ϵ/R$. If $2c>K$, even an orthogonal orientation fails to reach that drop, so this local-width criterion does not exclude any orientation. At $R=0$ or equal noise scales, the observation has no angular preference in this model.

For an original example, choose $R=2$, $\sigma =0.8$ and $ϵ=0.04$ in consistent signal units. Then $K=2493.75$. A drop $c=1/2$ gives half-width about $0.0200$ radians, or $1.15^{\circ }$. Reducing only the transverse standard deviation to $0.02$ approximately halves the width. Increasing a numerical regularization parameter is not an experimental improvement in this noise scale.

The global orientation ambiguity in a rank-one covariance fit remains: orientations separated by $\pi$ give the same covariance. A narrow local lobe does not distinguish those two orientations. Nor is a fixed log-density drop automatically a confidence interval with a stated coverage probability. Multiple observations, unknown noise scales and experimental nuisance parameters require their corresponding inference model.

Repeated reference spectra under stable conditions can estimate the longitudinal and transverse noise components. Their covariance must include shared calibration effects, as required by the measurement-propagation rule.^[ref] The numerical example is a model calculation in spectral coordinates, rather than a specification of a particular optical instrument.

Global orientation ambiguity in a rank-one covariance fit means that covariance identifies an axis without identifying its positive direction. For ${\mathbf{u}}_{\theta }=(\cos \theta ,\sin \theta {)}^T$ and $\lambda >0$,because ${\mathbf{u}}_{\theta +\pi }=-{\mathbf{u}}_{\theta }$. Adding an isotropic transverse-noise term does not remove this equality. The angular resolution set by transverse spectral noise can make each local orientation lobe narrow while leaving the two directed orientations indistinguishable.

For an original example, let $\lambda =4$ and $\theta =30^{\circ }$. Both $30^{\circ }$ and $210^{\circ }$ giveThe eigenvalues are four and zero, and the principal eigenvector can be written as either $(\sqrt{3}/2,1/2{)}^T$ or its negative. More samples can improve estimation of this matrix without choosing between those signs.

In a zero-mean Gaussian covariance model, equal covariance gives equal probability laws, including the supported line in the singular case.^[ref] No fitting procedure using only that law can distinguish the two parameter labels. A fitting range $0\le \theta <\pi$ removes duplicate labels by convention; it does not measure the sign of a physical direction.