Lock-in piezoreflectance (PzR) yields derivative-like signatures at direct optical transitions across a broad set of van der Waals (vdW) crystals and semiconductor microstructures. Figure 6 shows typical raw output from a single spot in WS₂ acquired with the workflow in Sections 4–6: strain-modulated AC reflectance ΔR, reference DC component proportional to reflectance R (Figure 6A), and the baseline-corrected ratio ΔR/R (Figure 6B). With a piezoceramic disk (thickness TH = 1 mm, OD = 30 mm) driven at 280 Hz and 100 V, ΔR/R peak magnitudes on the order of 10⁻5–10⁻4 are routinely obtained using lock-in time constants of 3 s to 30 s, allowing clear visualization of weak excitonic resonances that are often ambiguous or invisible in reflectance contrast (RC).

Figure 6: Representative outputs of the piezoreflectance workflow: strain-modulated reflectance and calculated PzR. (A) Strain-modulated AC reflectance ΔR (red, left axis) and chopper-referenced DC reflectance R (blue, right axis) from a single spot, measured at room temperature using a piezoceramic disk (TH = 1 mm, OD = 30 mm) driven at 280 Hz and 100 V. (B) Baseline-corrected piezoreflectance spectra (PzR) obtained by calculating ΔR/R. Please click here to view a larger version of this figure.
The ΔR/R baseline is typically flat (<1 x 10⁻5) in non-resonant regions when the spectral range extends beyond all strong features (Sections 4.11, 6.3). Residual curvature usually indicates an insufficient spectral window or slow lamp performance drift; both are corrected by extending the scan and/or subtracting a linear baseline from ΔR prior to forming ΔR/R (Section 6.3–6.4).
To quantify optical transitions in WS₂, the ΔR/R spectrum in Figure 7A (black) was fit with the Aspnes third-derivative functional form (Section 6.6, Equation 1). The resulting fit (cyan) reproduces all features; Figure 7B decomposes the model into four components (colored curves) with the dashed line showing their sum. The gray lines at the bottom of Figure 7A depict the phase-independent moduli Δρi(E) for direct transitions computed from Equation 2 using the parameters returned by the Aspnes fit, providing a convenient measure of the relative strength of each direct transition. The transition energies represented by peaks 1-4 are, respectively, 1.829 eV, 1.991 eV, 2.291 eV, and 2.421 eV. The energies obtained for the A and B transitions are consistent with prior modulation-spectroscopy reports on WS229.

Figure 7: Aspnes-model fitting of WS₂ piezoreflectance and decomposition into critical-point components. (A) Room-temperature PzR spectrum of bulk WS₂ (black; plotted as 104 ΔR/R vs. photon energy) with best fit (cyan) obtained using the Aspnes third-derivative functional form. Thin gray curves show the decomposition into the individual modulus components of ΔR/R used in the fit. (B) Component line-shapes from the Aspnes model corresponding to panel A: Fit Peak 1 (green), Fit Peak 2 (red), Fit Peak 3 (blue), and Fit Peak 4 (magenta). The dashed black line represents the cumulative fit (sum of all components). Please click here to view a larger version of this figure.
Dependence on measurement parameters
Here, we assess how the measurement parameters specified in the protocol influence the experimental outcome. The parameter choices are not universal: each optical setup must be calibrated first on well-characterized reference materials such as common transition-metal dichalcogenides (TMDs, namely WS₂, WSe₂, MoS₂, and MoSe₂), whose excitonic transitions lie in the measurement window for the silicon diode utilized as a detector and provide known resonance energies. The absolute magnitude of the detected ΔR is mainly determined by fixed optical factors – the halogen-lamp output, transmission losses through lenses/beam splitters/mirrors, and the illumination spot size defined by the aperture, objective, and monochromator slit. After alignment, the electronic settings determine the signal-to-noise ratio, in particular the preamplifier/amplifier gain and bandwidth, the modulation frequency and alternating voltage amplitude applied to the piezoceramic and the lock-in time constant. When preparing a new setup or starting measurements for a new material family, these parameters need to be optimized with signal-to-noise ratio (SNR) as the primary figure of merit, ensuring the possibility of observing resonances in measured spectra.
Two practical dependencies are illustrated in the figures. First, the strain-drive scaling (Figure 8): with τ = 1 s, the ΔR line shape in WS₂ is preserved while the amplitude increases monotonically as the AC voltage is increased from 100 to 1200 V (Figure 8A). The peak magnitude at selected energies scales linearly with the voltage amplitude (Figure 8B), confirming the operation in the elastic regime of the transducer and the coherent strain coupling at the resonances. Second, the lock-in time constant (Figure 9A): increasing τ reduces high-frequency noise and reveals weaker features, as seen by comparing MoS₂ spectra at 3 s, 10 s, and 30 s measured with the same AC amplitude of 100 V. The trade-off is a greater susceptibility to baseline drift and longer scan durations; in practice, τ ≈ 3-10 s offers a robust compromise that preserves derivative-like line shapes while delivering adequate signal-to-noise ratio. However, for some materials, it can become impossible to obtain satisfying PzR spectra without setting the time constant to 30 s. These calibrations define an operating window that can be reused across samples and revisited after any optical realignment.

Figure 8: Drive-voltage dependence of ΔR: amplitude scaling and elastic-regime verification. (A) Room temperature spectra acquired with lock-in time constant τ = 1 s as the AC drive is increased from 100 V to 1200 V (light to dark blue). The ΔR line shape is preserved while the amplitude grows monotonically. (B) Peak magnitudes at selected energies (1.29 eV, 1.58 eV, 1.88 eV, 1.93 eV) vs. applied voltage with linear fits (red), confirming elastic-regime transducer operation and coherent strain coupling at the resonances. Please click here to view a larger version of this figure.

Figure 9: Measurement optimization and cross-validation: lock-in time constant effects and comparison to photoreflectance. (A) Room temperature MoS₂ PzR spectra acquired at 100 V with lock-in time constants τ = 3 s, 10 s, and 30 s (blue to red to gray). (B) Room temperature results obtained for bulk WS₂: PzR data (black) with Aspnes fit (cyan) compared to photoreflectance (PR, orange). Please click here to view a larger version of this figure.
Comparison with photoreflectance
Photoreflectance (PR) is a pump-probe variant of modulation spectroscopy in which a spectrally scanned probe beam measures the relative reflectance change ΔR/R induced by a periodically modulated pump. Moreover, PR is contactless electroreflectance: the chopped pump creates excess carriers that modulate built-in electric fields (typically in surface/depletion regions), thereby perturbing the complex dielectric function and, as a result, the sample reflectance. Because the detection is phase-sensitive (lock-in at the pump modulation frequency), PR yields derivative-like line shapes with high sensitivity to interband critical points while suppressing a large, slowly varying background12,13,30.
For this section, in addition to the PzR signal, the PR signal was measured using the lock-in technique, which also allows extracting: (i) weak AC signals proportional to ΔR from the background and (ii) the DC component (proportional to reflectance, R). Both components were recorded with the same Si PIN photodiode and lock-in detection used for PzR. The relative changes in the reflection spectrum (ΔR/R) were evoked by modulation of the built-in surface electric field of the investigated sample, generated by illuminating the sample with a CW 405 nm laser that was mechanically chopped at 280 Hz. All measurements were performed in a dark configuration30.
Figure 9B presents the results for WS₂, a representative of the material class called transition-metal dichalcogenides with the general formula MX2 (where M = Mo, W, and X = S, Se, Te). In WS₂, PR produces clear derivative-like features at the A/B excitons, and the transition energies agree with those obtained from PzR (cyan Aspnes fits over the black PzR traces). This reflects that both modulation spectroscopies probe the same optical critical points in WS₂. The match between PR features and the Aspnes fits to PzR confirms that the PzR response is governed by the same excitonic dielectric-function derivative, validating PzR as a promising alternative to PR for extracting excitonic transition energies in TMDs.
Piezoceramic Actuator Selection
To generate a measurable piezoreflectance signal, the primary actuator-selection criterion is the achievable strain amplitude at the sample. In practice, we therefore prioritize piezoceramic elements with a large piezoelectric charge/strain coefficient (piezo modulus) dij, which quantifies either the induced charge density per applied stress or, the mechanical strain per applied electric field. For a thickness-poled disk actuator driven by a voltage (U) across its thickness (TH), the in-plane deformation is conveniently estimated in the small-signal regime by the radial (effective) coefficient d31: the disk outer diameter changes by ΔOD ≈ d31(OD/TH)U, corresponding to an average in-plane strain ε ≈ ΔOD/OD ≈ d31U/TH. Because thin disks satisfy OD >> TH, this geometry is efficient for producing in-plane strain at modest voltages. When the sample is placed on the top electrode, a thickness-driven disk typically produces an approximately biaxial in-plane (radial) strain (Figure 10), which is sufficient for lock-in detection of ΔR/R as the optical probe measures a spot-averaged modulation. The strain profile becomes important primarily when symmetry is crucial to the scientific question (e.g., uniaxial strain to deliberately break in-plane rotational symmetry or split anisotropic features), whereas for routine enhancement/derivative detection of weak optical transitions, the repeatable strain amplitude is usually the dominant requirement. For completeness, note that the lowest mechanical resonance of thin disks is commonly the radial extensional mode, with a characteristic scaling fs ≈ Np/OD, where Np is the frequency coefficient of the planar oscillation of a round disk; while resonant driving can increase strain, we generally operate well below resonance for linear response and stable phase, and practical resonance matching is often constrained by available actuator geometries/materials. Finally, electrode design can matter for reliability and mounting: thin-film electrodes deposited by PVD (e.g., sputtering) are typically ~1µm thick; we use silver electrodes but encourage testing alternative metallizations where needed. For an overview of available actuator shapes (disks, plates, shear elements), polarization directions, and electrode options, the PI Ceramic catalog provides a useful reference31.

Figure 10: Thin piezoceramic disk actuator geometry and voltage-induced in-plane deformation. (A) Schematic of a thin, thickness-poled piezoceramic disk (outer diameter OD, thickness TH, with OD ≫ TH) equipped with electrodes on the top and bottom faces and driven by an applied voltage U across the thickness. On the left, the orthogonal coordinate system (1-2-3) used to describe poled piezoelectric ceramics is shown; the polarization vector P (right side of panel A) is parallel to the 3 (Z) axis. (B) Illustration of the dominant in-plane (radial) deformation produced under thickness excitation, yielding an approximately biaxial lateral strain at the top surface relevant for strain-modulated optical measurements. Please click here to view a larger version of this figure.
Supplementary Figure 1: User interface of the control and acquisition software. Representative screenshot of the UI of the control/acquisition application supporting experimental setup.Please click here to download this file.
Supplementary Figure 2: Preamplifier configuration for piezoreflectance measurements. Sample settings.Please click here to download this file.
Supplementary Figure 3: Lock-in amplifier configurations for AC (ΔR) and DC (R) reflectance detection. Sample settings (A) of the lock-in amplifier for the measurement of lock-in AC component proportional to ΔR and DC component proportional to R, and (B) of the chopper-referenced lock-in amplifier for measurement of the additional DC reflectance component (R).Please click here to download this file.
Supplementary Figure 4: Instrument settings for AC drive and optical chopping. Sample settings (A) of the optical chopper controller and (B) of the custom-made AC voltage generator.Please click here to download this file.