**A calibration load case must be known independently of the spectrum.** Four-point bending creates a nominally uniform uniaxial surface strain between inner loading points and is useful for bars or wafers, but thickness, support spacing, anisotropic elasticity, anticlastic curvature, and load alignment matter. Strain gauges, digital image correlation, displacement metrology, finite-element analysis, or diffraction should verify the strain actually present in the Raman sampling region. Hydrostatic pressure in a pressure cell provides a different stress state and can determine pressure coefficients over a broad range. Pressure medium hydrostaticity, pressure marker, phase stability, pressure gradients, and optical access limit accuracy. A hydrostatic coefficient cannot be substituted for an in-plane biaxial coefficient merely because both use gigapascals; their tensor contractions and mode splitting differ. Biaxial calibration can use membrane bulging, pressure-loaded windows, epitaxial standards, thermal-mismatch structures, or calibrated wafer curvature with a verified film model. Each introduces assumptions about adhesion, thickness, elastic anisotropy, edge effects, plasticity, and stress uniformity. An epitaxial layer may provide a well-defined in-plane strain from x-ray diffraction, but composition, relaxation, defects, and thermal history must be measured. Nanoindentation and patterned test structures create rich multiaxial fields valuable for validating spatial maps. Their stress state is not known from force alone; contact mechanics or finite-element models and independent deformation measurements are required. Near edges, cracks, interfaces, and free surfaces, continuum assumptions and optical averaging become especially important. Such structures are better validation artifacts than primary scalar calibrators unless the mechanics are tightly constrained. |Calibration route|Best-established quantity|Main advantage|Dominant limitation|Essential validation| |---|---|---|---|---| |Four-point bending|Surface uniaxial strain or stress in a central region|Reversible loading and multiple calibration points|Alignment, anisotropy, thickness, anticlastic bending|Strain gauge or DIC plus elastic model| |Hydrostatic pressure cell|Pressure coefficient|Broad, symmetric loading range|Hydrostaticity and mismatch to device stress state|Independent pressure marker and phase check| |Biaxial membrane or bulge|In-plane biaxial stress/strain|Closer to many thin-film boundary conditions|Geometry, edge effects, thickness, nonlinear deflection|Profile metrology and membrane mechanics| |Epitaxial reference series|Composition- and orientation-specific lattice strain|Process-relevant material stack|Composition–strain covariance and partial relaxation|Reciprocal-space x-ray mapping| |Patterned or indented validation artifact|Spatially varying multiaxial field|Tests mapping and tensor reconstruction|Model dependence and gradients below optical resolution|Finite-element model plus independent displacement or diffraction| **Temperature, composition, carriers, and phase must be separated from mechanics.** A practical peak-shift model is $$ \Delta\omega_m=\mathbf{P}_m:\boldsymbol{\epsilon}+\chi_{mT}\Delta T+\chi_{mc}\Delta c+\chi_{mn}\Delta n_c+\Delta\omega_{phase}+\cdots $$ The deformation-potential term is only one contribution. Laser heating, device self-heating, alloy fraction, doping, free carriers, isotope content, phase transformation, damage, and resonance can move or reshape the same band. Calibration specimens and unknowns should match these variables or include independently measured corrections. Temperature compensation should use a low-stress, composition-matched specimen over the relevant temperature range and optical conditions. A linear coefficient may be adequate over a narrow interval, but anharmonicity and thermal expansion can create curvature. In a powered device, temperature and stress change together; using a single peak cannot generally solve both. Multiple phonons with distinct temperature and strain sensitivities, a Stokes/anti-Stokes ratio, or an orthogonal thermometer can make the system identifiable. Alloy calibration needs at least enough independent observables to separate composition and strain. SiGe, III–V alloys, nitrides, and ternary or quaternary systems can show multiple bond-related modes, local ordering, clustering, and composition-dependent deformation potentials. X-ray diffraction, composition metrology, and relaxed reference films anchor the model. A coefficient trained on one growth method may not transfer when ordering or defect content changes. Carrier density can cause phonon self-energy shifts, linewidth changes, and asymmetric Fano coupling; polar materials can exhibit longitudinal-optical phonon–plasmon coupled modes. Electric fields can also produce inverse piezoelectric strain or modify phonon frequencies through additional coupling. Bias-dependent Raman maps therefore need electrical, thermal, and electromechanical controls before a shift is labeled mechanical stress. Phase and damage checks precede quantitative conversion. High pressure, indentation, machining, laser exposure, or process excursions can transform crystal structure or amorphize a region. Applying the original phase’s coefficient to a transformed peak is meaningless. Peak inventory, polarization, linewidth, and an orthogonal structural measurement should confirm that the calibration phase remains intact throughout loading. **Diffraction measures lattice strain and requires its own reference and geometry.** Bragg’s law is $$ 2d\sin\theta=m\lambda $$ and small changes at fixed wavelength give $$ \frac{\Delta d}{d}\approx-\cot\theta\,\Delta\theta $$ when $\Delta\theta$ is expressed in radians and peak-angle conventions are consistent. This returns the lattice-strain projection normal to the diffracting planes. Converting it to a stress tensor requires elastic constants, grain interaction assumptions, specimen orientation, and enough independent diffraction vectors. The stress-free lattice spacing $d_0$ is often the dominant uncertainty. Composition, temperature, defect concentration, chemistry, and ordering change $d_0$. In thin films, conventional symmetric scans may provide only out-of-plane strain, while device performance depends on in-plane strain. Reciprocal-space maps, asymmetric reflections, grazing incidence, or multiple specimen tilts can add components, but penetration depth and spatial resolution differ from Raman. Cross-calibration should compare compatible spatial and tensor averages. A micron-scale Raman spot, millimeter-scale x-ray beam, wafer-curvature average, and nanometer-scale electron-diffraction measurement do not observe the same field. Agreement may be accidental if tensile and compressive regions average differently. Register coordinates, model each point-spread or gauge volume, and compare the forward-predicted observable rather than raw “stress” maps. Wafer curvature can estimate average film stress when a uniform film is much thinner than its substrate, curvature is small, and the biaxial modulus is known. Patterned films, multilayers, anisotropy, or stress gradients require generalized models. Curvature is useful for wafer averages, while Raman resolves local departures. **Spatial resolution and sampling depth turn local stress into an optical average.** A confocal Raman voxel has finite lateral and axial weighting set by wavelength, numerical aperture, refractive index, absorption, pinhole, aberration, and the layered stack. If stress varies within that volume, the spectrum is an integral over shifted local responses: $$ I(\omega,\mathbf{r}_0)=\int W(\mathbf{r}-\mathbf{r}_0)\,L[\omega-\omega_0-\Delta\omega(\mathbf{r})],d\mathbf{r} $$ where $W$ is the optical weighting and $L$ is the local line shape. A fitted peak center is a weighted statistic of the distribution; it is not necessarily the stress at the voxel center. Broadening and asymmetry can contain gradient information but are also affected by defects, temperature, and resolution. Mapping with a step smaller than the spot size oversamples the optical field; it does not create independent nanoscale resolution. Deconvolution can improve localization only with a measured point-spread function, adequate signal, and regularization whose bias is quantified. Tip-enhanced Raman can shrink the near-field sampling region, but enhancement variation, tip stress, heating, polarization, and far-field background introduce a new calibration problem. At free surfaces and patterned edges, mechanical relaxation changes the field, while optical focus and collection also change. Topography can correlate with apparent Raman shift through defocus, aberration, or mixed material signal. Co-registered height, reflectance, phase, and fit-quality maps help distinguish mechanics from optics. Changing laser wavelength or focus changes depth weighting, absorption, and resonance. Differences are not direct depth derivatives; they require an optical and layered-stress model. **Regression and uncertainty determine whether calibration transfers.** A calibration should include multiple loading and unloading points, repeats, independently verified zero, and coverage of the intended operating range. Plot residuals against load, time, position, temperature, and signal level. Hysteresis or drift can reveal slip, plasticity, mounting change, heating, phase evolution, or instrumental motion. Both axes have uncertainty: the reference stress or strain is not exact, and the spectral shift has fit and calibration error. Ordinary least squares can bias the slope when reference uncertainty is material. Orthogonal-distance, generalized least-squares, hierarchical, or errors-in-variables models may be appropriate. Correlated uncertainties—such as one thickness value used for every load point—must not be treated as independent random noise. The uncertainty budget can be expressed schematically as $$ u_y^2=\mathbf{J}\mathbf{U}_x\mathbf{J}^{T}+u_{model}^2+u_{repeat}^2 $$ where $\mathbf{J}$ contains sensitivities of the reported stress or strain to inputs, $\mathbf{U}_x$ is their covariance matrix, and the remaining terms represent model inadequacy and repeatability. Inputs can include peak center, spectral calibration, temperature, composition, coefficient, elastic constants, orientation, thickness, load, geometry, and reference state. Precision is not accuracy. A spectral center repeatable to a small fraction of a wavenumber can still produce biased stress through a wrong coefficient, temperature drift, reference offset, or boundary condition. Report repeatability, calibration uncertainty, spatial reproducibility, and model uncertainty separately. Validation on a withheld specimen or geometry tests transfer better than a high coefficient of determination on the calibration data. Calibration validity should be bounded by material, phase, orientation, stress state, temperature, composition, optical configuration, and load range. Extrapolation needs new validation, and coefficients should retain versioned provenance. ```flowchart Define the required strain or stress components and coordinate system -> Choose a material-, orientation-, and geometry-matched reference series -> Apply reversible load while independently measuring strain or stress -> Control temperature, composition, carriers, phase, and optical configuration -> Acquire polarized spectra and fit components with fixed quality rules -> Regress shifts against verified tensors with errors on both axes -> Build uncertainty, hysteresis, gradient, and transfer-validity budgets -> Test the calibration on a withheld structure and orthogonal method -> Deploy only within the validated material and state domain ``` **A production calibration is a versioned measurement model, not a coefficient lookup.** Store the specimen identity, crystal and device coordinates, phase, composition, thickness, elastic constants, deformation potentials or empirical slopes, load geometry, reference state, temperature, optical recipe, peak model, regression code, covariance, residuals, validity limits, and approval history. Raw spectra and reference-load data must remain recoverable. For each unknown, report the measured shift and linewidth, selected mode component, temperature and composition corrections, inferred strain or stress components, expanded uncertainty, fit quality, and whether the point lies inside the calibration domain. Reject pixels or specimens with phase mismatch, unresolved splitting, excessive gradients, saturation, low signal, or extrapolation unless a separate model handles them. The most defensible workflow predicts what every instrument should observe from one mechanical state. Raman, x-ray diffraction, curvature, microscopy, and device simulation are then compared at their native spatial weighting and tensor projection. Disagreement becomes diagnostic evidence about references, gradients, material properties, or missing physics rather than something hidden by adjusting a scalar conversion factor. The durable way to use stress–strain calibration is through a reference-state-tensor-deformation-potential-elasticity-confounder-spatial-weighting-regression-and-traceability lens.