Raman mapping

Raman mapping records Raman spectra at known sample positions and constructs spatial maps from their spectral features. A map can show composition or structural variation that is absent from an ordinary reflected-light image. A point-scanning apparatus combines a focused laser, microscope objective, rejection of unshifted light, a spectrometer and controlled sample positioning. A primary university facility description identifies single-point acquisition and two-dimensional mapping as measurement modes.^[ref]

The acquired data can be represented as $I(x,y,\tilde{\nu })$, with position coordinates $x,y$ and Raman shift $\tilde{\nu }$. It is a form of hyperspectral imaging. The displayed image is a derived summary of those spectra. For example, a baseline-corrected Raman band-area map useswhile a Raman peak-position map records a fitted band's position. Different summaries can reveal different properties of the same sample. Manufacturer examples include a silicon-carbide wafer whose spectral map separates inclusions and structural variation.^[ref]

For an illustrative20µm-by20µm point map sampled every1µm with both boundaries included, there are $21\times 21=441$ positions. At0.10s exposure per position, exposure time alone is44.1s; motion and readout add overhead. Reducing step size increases sampling density, while optical resolution remains constrained by the excitation/detection point spread function and instrument configuration. A finely interpolated display cannot recover an unsampled small inclusion.

Intensity also depends on illumination, focus and the spectral background. A single-frequency brightness map can mix a desired Raman band with broad background, whereas fitted or integrated features can answer a more specific question. Quantitative composition or stress estimates require an appropriate measurement calibration; a colourful spatial contrast is not by itself a calibrated concentration or stress measurement.^[ref]

A Raman peak-position map assigns to each sample position the fitted Raman shift of a selected spectral band. Starting from position-resolved spectra, the analysis estimates ${\tilde{\nu }}_{\mathrm{peak}}(x,y)$ and displays its spatial variation. It answers a different question from a Raman band-area map, which records an integrated spectral intensity. Manufacturer documentation uses fitted band position to reveal structural variation, including strain-related changes in a silicon-carbide wafer.^[ref]

For a simple isolated Gaussian function band above background $B$,the fitted center ${\tilde{\nu }}_0$ is the peak position. A maximum chosen from discrete channels can instead jump in steps of the channel spacing. Fitting uses the neighboring samples too, although its result still depends on the line-shape model, background and noise. Raman peak-fit uncertainty belongs with the resulting map.

For an explicitly hypothetical band, suppose a reference region is at1000.0cm⁻¹ and another region at1000.5cm⁻¹. The shift difference is0.5cm⁻¹. If an independently established local stress calibration were $\Delta \tilde{\nu }=K\sigma$ with $K=2\mathrm{c}{\mathrm{m}}^{-1}/GPa$, that model would imply0.25GPa. The coefficient is illustrative, not a material constant asserted for the manufacturer's wafer. A quantitative stress map requires the applicable measurement calibration and control of other causes of spectral change.

Overlapping bands can move the maximum of their sum when only their relative intensities change. A stable center estimate must therefore track a specified band rather than interpret every maximum shift as strain. Instrument drift, temperature and laser-induced heating can also affect position. Retaining the underlying spectra and fit results makes it possible to test those interpretations rather than rely only on the map's colours.

Raman peak-fit uncertainty characterizes the uncertainty of a spectral-band position estimated by fitting a model to a Raman spectrum. It accompanies a Raman peak-position map: a colour difference between pixels need not establish a physical shift if the estimated positions are insufficiently precise. In Raman mapping, the instrument collects spectra at successive sample positions, and fitting extracts band parameters such as position and width.^[ref]

For an isolated band, a least squares fit might estimate the center, height, width and background of a Gaussian function. The uncertainty assigned to its center is conditional on the chosen model and description of spectral noise. Repeating the complete acquisition and fit at one unchanged location provides another assessment: if $n$ independent center estimates have standard deviation $s$, the estimated standard uncertainty of their sample mean is $s/\sqrt{n}$. This reduction applies to independent variation; it does not eliminate an offset shared by every acquisition.^[ref]

For an original numerical example, model two fitted positions aswhere $c$ represents a shared uncertain additive correction from measurement calibration, while $e_0,e_1$ are independent fitting errors, also independent of $c$. Suppose their standard uncertainties are $u(c)=0.10\mathrm{c}{\mathrm{m}}^{-1}$ and $u(e_i)=0.05\mathrm{c}{\mathrm{m}}^{-1}$. Each position then has combined standard uncertainty $\sqrt{0.10^2+0.05^2}\approx 0.112\mathrm{c}{\mathrm{m}}^{-1}$. Their difference cancels the same $c$, givingHere $u$ denotes standard uncertainty, so the difference has uncertainty about $0.0707\mathrm{c}{\mathrm{m}}^{-1}$. The covariance term is essential: treating the two positions as entirely independent would instead give about $0.158\mathrm{c}{\mathrm{m}}^{-1}$. A common additive correction cancels only insofar as it really is shared; instrument drift between measurements or a scale error need not cancel.^[ref]

For a hypothetical stress model $\sigma =({\hat{\nu }}_1-{\hat{\nu }}_0)/K$ with an exactly known illustrative $K=2\mathrm{c}{\mathrm{m}}^{-1}/GPa$, this difference uncertainty contributes $0.0354\mathrm{GPa}$ to the standard uncertainty of stress. Uncertainty in $K$ would require its own contribution. Raman peak-model bias, overlapping bands, an inadequate background and laser-induced sample changes can invalidate an otherwise small fit error. Retaining spectra and repeating fits with physically defensible models therefore matters alongside quoting a numerical uncertainty.

Raman peak-model bias is a systematic displacement of an estimated Raman shift or another fitted band parameter caused by an inadequate spectral model. The discrepancy can persist in repeated low-noise measurements. It differs from the random scatter contributing to Raman peak-fit uncertainty. Manufacturer documentation explains that fitted band parameters depend on theoretical curves and warns that a broad spectral background can contaminate simpler intensity measurements.^[ref]

An original three-channel example isolates a Raman baseline-induced center shift. Write $x$ for displacement from the true band center and let the true sampled signal be a Gaussian function plus a small residual baseline slope:Suppose an equal-weight least squares fit omits that slope and estimates only the center $\mu$, keeping $H$ and $s$ fixed:Near $\mu =0$, its derivative with respect to the center is $j_i=(H{x}_i/s^2)e^{-x_i^2/(2s^2)}$, so linearizing gives $m(x_i;\mu )\approx m(x_i;0)+\mu j_i$. Minimizing the resulting squared residuals givesThe odd slope in the unmodeled background therefore projects onto the odd center-shift response of the symmetric band.

For $H=10$ intensity units, $s=2\mathrm{c}{\mathrm{m}}^{-1}$ and $b=0.05$ intensity units per $\mathrm{c}{\mathrm{m}}^{-1}$, the first-order center bias is about $0.0330\mathrm{c}{\mathrm{m}}^{-1}$. This is a deliberately restricted toy fit, not a universal bias formula: fitting the background, width and height too, or using other sample positions and statistical weights, changes the result. The calculation nevertheless demonstrates why reducing random noise alone cannot remove a wrong-model displacement.

In Raman mapping, a spatially varying background or unresolved neighboring band can consequently create an apparent Raman peak-position map contrast. Comparing physically defensible background/line-shape models and examining their residuals tests that interpretation. A good residual plot is useful evidence but does not prove uniqueness of a spectral decomposition. The GUM requires significant corrections and inadequately represented influences to be included in the measurement model; its discussion of fit-parameter uncertainties does not make an incomplete model adequate.^[ref]

A Raman baseline-induced center shift is a displacement of an estimated band center caused by the spectral background used or omitted in the analysis. It is one mechanism of Raman peak-model bias. A background slope can move the local maximum of a band-plus-background signal even if the underlying band has not shifted. A fitted center can also move when an incomplete background model transfers part of that slope into the fitted band position.

For an original local model, write $x$ for Raman shift relative to the true center, and consider a Gaussian function band plus a small residual slope:Its stationary-point equation near the band center is $-(Hx/s^2)e^{-x^2/(2s^2)}+b=0$. When the displacement is small compared with $s$, the local maximum is displaced byFor $H=10$ intensity units, $s=2\mathrm{c}{\mathrm{m}}^{-1}$ and $b=0.05$ intensity units per $\mathrm{c}{\mathrm{m}}^{-1}$, this is approximately $0.020\mathrm{c}{\mathrm{m}}^{-1}$. A positive residual slope moves this local maximum toward increasing Raman shift. This local spectral-window model does not assert a physically unbounded linear background.

A least squares fitted center answers another question. If $j_i$ is the derivative of a chosen band model with respect to its center, omitted residual background $r_i$ produces the first-order fitted displacementfor fixed other parameters and positive weights $w_i$. An omitted background component matters when it projects onto the center-shift response. For symmetric samples and weights, a constant residual is orthogonal to the odd derivative of a symmetric band, whereas a residual slope is not. Asymmetric sampling or weights can remove even that constant-background cancellation.

The three-channel fit in Raman peak-model bias gives approximately $0.0330\mathrm{c}{\mathrm{m}}^{-1}$ for the same slope and band parameters, because estimating a center from those samples is not the same operation as finding the continuous local maximum. At sample positions $-s,0,s$, the center derivative is proportional to $x$; allowing both a center change and an unknown linear background therefore introduces a first-order identifiability problem in that restricted model. Additional spectral samples or independent background information can help separate their responses.

Examining residuals and changing physically justified background models are useful tests before interpreting a Raman peak-position map as a map of material changes. Manufacturer guidance distinguishes spectral intensity, background contamination and fitted curve parameters.^[ref] A small conditional Raman peak-fit uncertainty alone does not establish that the assumed background is adequate.^[ref]

A Raman band-area map assigns each spatial position an integrated, usually background-corrected spectral intensity over a stated Raman shift range. If $I(x,y,\tilde{\nu })$ is the spectrum and $B$ is its estimated baseline,This is one way to summarize the spectra collected in Raman mapping. It is useful for following a band's spatial distribution while preserving a different interpretation from a Raman peak-position map.^[ref]

Peak height alone does not determine area. For an isolated Gaussian function band of height $H$ and width parameter $s$, integrating the complete band givesTwo hypothetical bands of equal height10 intensity units and $s=2$ or4cm⁻¹ have full areas about50.13 and100.27 intensity-units·cm⁻¹. Thus an area map can change while a height map stays constant. A finite integration window captures only part of the complete band, so the actual window must accompany the value.

Raman integration-window bias arises if shifting or broadening moves part of the chosen band outside that fixed range. Neighboring bands and an incorrectly subtracted background can also change the result. For a sampled spectrum, a quadrature uses the actual channel spacing, for example $A\approx {\sum }_i(I_i-B_i)\Delta {\tilde{\nu }}_i$. A bare sum without the spacing factor has a different numerical convention when sampling changes.

Comparing areas across positions also requires consistent exposure, illumination and focus, or an appropriate normalization and measurement calibration. A larger area is therefore not automatically a larger concentration of one material. The manufacturer's discussion of spectral analysis explicitly cautions that broad background can contaminate simple intensity-based images.^[ref]

Raman integration-window bias is a systematic difference between the spectral-band area one intends to compare and the area captured by a selected finite Raman shift interval. It can create contrast in a Raman band-area map when the complete band area is unchanged: shifting or broadening redistributes intensity across the fixed interval's edges. The target must be stated clearly, because an area over that interval can be correctly measured even when it is a biased proxy for the complete band.

For an original calculation, take an isolated background-corrected Gaussian function bandIts complete area is $A_{\mathrm{full}}=Hs\sqrt{2\pi }$. Integrating only from $a$ to $b$ gives the Raman captured-area fractionwhere $\Phi$ is the cumulative distribution function of the standard normal distribution. The substitution $z=(\tilde{\nu }-{\tilde{\nu }}_0)/s$ gives this expression directly. The complete area depends on height and width; the captured fraction also depends on center position relative to the window.

Suppose $s=2\mathrm{c}{\mathrm{m}}^{-1}$ and the fixed window runs from996 to1004cm⁻¹. A center at1000cm⁻¹ places its edges at two standard widths on each side, giving $f=\Phi (2)-\Phi (-2)\approx 0.95450$. Shift the center to1004cm⁻¹ without changing height or width: now $f=\Phi (0)-\Phi (-4)\approx 0.49997$. The measured fixed-window area falls to about52.38% of its previous value although the complete band area is identical. This contrast is therefore insufficient evidence for a decrease in the amount of material.

A window that follows the fitted center can keep the captured fraction constant for this ideal fixed-width model. Broadening, neighboring bands, errors in a Raman peak-position map and an imperfect background still matter; a larger window can admit unwanted signal. Manufacturer documentation explicitly distinguishes simple intensity images from fitted-band analysis and warns that broad background can contaminate intensity contrast.^[ref] Reporting the actual interval, background treatment and whether the interval moves with the band makes the Raman band-area map scientifically interpretable.

The Raman captured-area fraction is the part of a specified complete spectral band's area that lies inside a chosen integration interval. For a nonnegative isolated background-corrected band $S(\tilde{\nu })$ with finite positive complete area,It lies between zero and one. It quantifies one mechanism of Raman integration-window bias: changing the center or width changes the fraction even when the complete area stays fixed. Background errors or overlapping bands can prevent measured intensity from representing the assumed isolated band.

For a Gaussian function centered at ${\tilde{\nu }}_0$ with width parameter $s$, substitution of $z=(\tilde{\nu }-{\tilde{\nu }}_0)/s$ giveswhere $\Phi$ is the cumulative distribution function of the standard normal distribution. A centered window with half-width $ks$ captures $2\Phi (k)-1$. Half-widths of one, two and three $s$ capture approximately68.27%,95.45% and99.73% of the complete Gaussian area. These percentages refer to this line shape, not every Raman spectrum.

For an original example, a window extending4cm⁻¹ either side of the fitted center captures about0.95450 when $s=2\mathrm{c}{\mathrm{m}}^{-1}$, but about0.81758 when $s=3\mathrm{c}{\mathrm{m}}^{-1}$. Even if the height decreases so the complete area is unchanged, the finite-window area falls. If one incorrectly uses the first fraction to correct the broader band, the estimated complete area is only $0.81758/0.95450\approx 0.85655$ of the true value.

Raman band-area extrapolation uses a justified fraction to estimate complete area as $A_{\mathrm{full}}=A_{[a,b]}/f$. For example, measured area10 in a centered two-width Gaussian window implies about10.477 in the same area units. The correction depends on the model and its parameters; it is not extra measured signal. At fixed $f$, an absolute area error is multiplied by $1/f$, while uncertainty in $f$ supplies another contribution. Large extrapolations from a narrow interval therefore need especially convincing evidence about the unobserved tails.

Collecting and retaining the full local spectra supports checks of center, width, background and neighboring bands instead of interpreting a single intensity summary alone. Primary manufacturer guidance makes this distinction between simple spectral-intensity images and fitted band parameters explicit.^[ref]

Raman band-area extrapolation estimates the complete area of a spectral band from a measured finite-window area and a model for the fraction outside that window. It uses the Raman captured-area fraction $f>0$:The extra area is inferred rather than observed. Its validity depends on the assumed line shape, width, center, background and separation from neighboring bands. Keeping these assumptions explicit prevents a corrected Raman band-area map from concealing the model dependence of its values.

For an original uncertainty example, let $A_{\mathrm{win}}=10.0$ in arbitrary area units with standard uncertainty0.20, and let $f=0.950$ with standard uncertainty0.010. The extrapolated area is about10.526. For independent inputs, the first-order combined standard uncertainty isin the same area units. Treating $f$ as exact would retain only the first term and understate the uncertainty in this example.

The fitted fraction and window area may be inferred from the same spectrum, so independence needs justification. The general first-order expression retains their covariance:The signs follow from the sensitivity coefficients $\partial A_{\mathrm{full}}/\partial A_{\mathrm{win}}=1/f$ and $\partial A_{\mathrm{full}}/\partial f=-A_{\mathrm{win}}/f^2$. The GUM's correlated-input propagation is the applicable measurement framework.^[ref]

A deterministic bound answers a different question. If $A_{\mathrm{win}}\in [9.8,10.2]$ and $f\in [0.94,0.97]$, positivity givesThis interval follows from allowed ranges and is not automatically a probabilistic confidence interval. Small $f$ makes both sensitivity coefficients large, so extrapolating from a small central portion can be unstable even when that portion is measured precisely.

Model inadequacy needs attention alongside parameter uncertainty. A background mistaken for a band tail, or an unresolved neighboring band, changes what is being extrapolated. Retaining the local Raman spectrum allows the fitted band parameters to be checked instead of relying solely on an intensity summary.^[ref]