The application of optical heterodyne detection for transient grating spectroscopy (TGS) using a fixed, binary phase mask often relies on taking the difference between signals captured at multiple heterodyne phases. To date, this has been accomplished by manually controlling the heterodyne phase between measurements with an optical flat. In this letter, an optical configuration is presented which allows for collection of TGS measurements at two heterodyne phases concurrently through the use of two independently phase controlled interrogation paths. This arrangement allows for complete, heterodyne amplified TGS measurements to be made in a manner not constrained by a mechanical actuation time. Measurements are instead constrained only by the desired signal-to-noise ratio. A temporal resolution of between 1 and 10 s, demonstrated here on single crystal metallic samples, will allow TGS experiments to be used as an in-situ, time-resolved monitoring technique for many material processing applications.

For almost two decades, transient grating spectroscopy (TGS), also known as impulse stimulated thermal scattering (ISTS), has been implemented using an optical heterodyne detection technique in an arrangement described in Ref. 1. This methodology has been used to study the acoustic and thermal transport properties of a wide variety of systems including thin-film semiconductors,2 microsphere lattices,3 and suspended silicon films,4 among many others.5 It has also been applied to study the magnetic properties of thin films,6 the electronic properties of 2D semiconductors,7 and transport properties in high-Tc superconductors.8 Similar methods have been used for grain mapping in polycrystalline metals.9 Recently, TGS has been shown to be a useful tool in evaluating the effects of irradiation in metallic systems.10–12 

In this technique, a spatially periodic material excitation is generated by periodic thermal excitation, itself produced by the interference pattern of two crossed pumping laser beams at the material surface. In reflective materials, this excitation manifests itself as surface displacement and/or transient reflectivity change. To interrogate the dynamics of this excitation, a third, probing laser beam is diffracted from the excitation. The first order diffraction of this probe is directed into a photodetector such that the amplitude of the response oscillates with the frequency of the material excitation. In TGS experiments to date, a fourth laser beam, generated from the same source as the probe, is used as a reference oscillator, allowing the diffracted probe signal to be heterodyne amplified. In this configuration, the measured signal intensity has the form

IS(t)=IR+ID(t)+2IRID(t)cosϕ,
(1)

where IR is the reference oscillator intensity, ID(t) is the diffracted signal intensity, and ϕ is the heterodyne phase difference between the reference and diffracted beams.1 As IR is constant and usually of much greater amplitude than ID(t), the result is an overall amplification of the time-dependent signal of interest.

To recover large amplitude signals suitable for the quantitative analysis of the acoustic and thermal transport properties of materials, it is a common practice to record traces at both ϕ=0° and ϕ=180°. These two signals can be subtracted from one another to recover the total signal amplitude

IStot(t)=IS0(t)IS180(t)=4IRID(t).
(2)

In addition to further amplifying the recorded signal intensity, taking a set of measurements in this manner allows for the removal of any systematic noise from signals used for quantitative analysis. Such noise could include the impulse response of the photodetector or other electrical noise present near the experiment. In practice, the heterodyne phase is controlled by a manually adjustable, highly parallel optical flat in the path of the probe laser beam. A phase difference between probing and reference beams is generated as a function of the tilt angle of this flat due to small changes in the path length.2 Measurements relying on this type of manual phase control are time-limited by the actuation time between collections at different values of ϕ, typically in the range of tens of seconds to single minutes. This limitation must be lifted if material systems undergoing dynamic changes are to be studied using TGS. Such systems of interest may include irradiation-induced material evolution, the thermal transport properties of surfaces as oxide layer growth occurs, and the kinetics of low-temperature phase change materials.

In this letter, a modification to the optical arrangement for heterodyne amplified TGS experiments is presented which allows for the simultaneous collection of two signals with differing heterodyne phases, hereafter referred to as dual heterodyne phase collection TGS (DH-TGS). For materials that exhibit strong excitations in both displacement and reflectivity under an imposed transient grating, different dynamics can be probed by choosing the values of ϕ other than ϕ=0° and ϕ=180°.2 Of these two excitations, it is the displacement dynamics that contains information about both acoustic and thermal properties, while the reflectivity dynamics contain information about thermal transport properties alone.13,14 However, the surface displacement response cannot be interrogated independently, causing the thermal transport dynamics of measurements sensitive to acoustic propagation to be a combination of both displacement and reflectivity responses.2 In light of this complexity, the performance of the DH-TGS methodology is characterized in terms of the measured acoustic response alone. For the case considered here, fixing ϕ to 0° and 180° allows for the collection of a complete TGS measurement of the excitation of interest without the constraint of mechanical actuation.

The optical arrangement for a reflective sample is shown in Fig. 1(a). The generation of the excitation grating on the sample under investigation is accomplished by passing a short pulsed laser through a volumetric diffraction grating, hereafter referred to as the “phase mask,” which is optimized to diffract most of the light into the ±1 orders. All other orders are spatially filtered, and the two pump paths are recombined at the sample using an imaging system of cemented achromatic doublets, L1 and L2. The wavelength of the intensity grating projected onto the sample is given by Λ = MΛ0/2, where Λ0 is the wavelength of the phase mask and M is the magnification ratio between L1 and L2.1 In this implementation, a matched pair of doublets is used for L1 and L2 such that M = 1.

FIG. 1.

(a) Top-view block schematic of the DH-TGS optical geometry showing paths from laser sources to detectors. BS – beam splitter, L1 and L2 – achromatic doublets, ND – neutral density filter, PA1 and PA2 – phase adjustment flats, and APD1 and APD2 – silicon avalanche photodiodes. (b) Beam arrangement in the plane of lens L2, just before the sample surface. The positions below the dotted line represent the reflection and diffraction of the reference oscillator and probe beams from the sample surface.

FIG. 1.

(a) Top-view block schematic of the DH-TGS optical geometry showing paths from laser sources to detectors. BS – beam splitter, L1 and L2 – achromatic doublets, ND – neutral density filter, PA1 and PA2 – phase adjustment flats, and APD1 and APD2 – silicon avalanche photodiodes. (b) Beam arrangement in the plane of lens L2, just before the sample surface. The positions below the dotted line represent the reflection and diffraction of the reference oscillator and probe beams from the sample surface.

Close modal

To record two heterodyne phases concurrently, two additional laser paths are added to the standard TGS implementation as a probing beam and reference oscillator pair. All of the probing beams are generated from the same continuous wave laser source. Prior to incidence onto the phase mask, the probe path is split using a non-polarizing 50:50 beam splitter (BS). The transmitted leg of this probe is focused coincident with the pump source onto the phase mask at an angle θ1 (out of the plane of the page in Fig. 1(a)) such that between L1 and L2, there is a vertical separation of the pump and first probing paths. The reflected probing beam from the beam splitter is re-directed so it, too, is coincident onto the phase mask as the same spot as the pump but with incoming angle θ2 > θ1 such that the vertical separation between the second probe paths and the pump is greater than that of the first pair. For both probes, the ±1 diffraction orders are again retained, while the transmitted beam and higher orders are blocked. The vertical placement of each of the six beam paths in the plane of L2 is shown in Fig. 1(b).

Between L1 and L2, the two reference oscillator beams are passed through a neutral density filter (ND) to avoid saturating the photodetectors downstream. Both probing paths are passed through a first phase adjust flat (PA1), which is used to uniquely control the heterodyne phase of IS180(t). The extremal probing path is passed through a vertically separated, second phase adjust flat (PA2), which adds an additional path length difference and uniquely controls the phase of IS0(t). Alternative geometries of the phase adjustment optics with the same independent control are possible. However, the configuration described here allows PA2 to be set and fixed to 180° offset from PA1, necessitating only the adjustment of PA1 when tuning the response of a new sample. Upon reflection and diffraction from the sample surface, the spatial coincidence of the reference oscillator and the diffracted signal for both probing paths is guaranteed by the optical geometry. In the reflection geometry, both signals pass back through L2. As these paths are now vertically separated below the plane of the pump paths, they can be picked off, separated, and directed into two fast silicon avalanche photodiodes (APD1 and APD2). The path lengths between the sample surface and APD1 and APD2 are managed such that the two signals are aligned in time when recorded on a digital oscilloscope.

The primary advantage of the optical arrangement described above over previous implementations of the technique is time resolution. The pumping laser in these types of experiments is commonly operated with a repetition rate in the single kHz range. Each measurement, therefore, is taken as an average over hundreds to thousands of laser pulses to reduce the system noise. Previously, a complete measurement of both heterodyne phases could be time resolution limited either by the signal collection time (in the case where many averages are taken) or by the mechanical action necessary to change the heterodyne phase (in the case where only a few averages are needed). By recording both these phases independently, the time resolution of the system becomes, in all cases, signal-to-noise ratio limited.

DH-TGS experiments using an optical arrangement as described above are carried out using a passively Q-switched, 532 nm, solid-state laser with a pulse length of 300 ps, a repetition rate of 1 kHz as a pump laser source with a spot size of 210 μm, and pulse energy at the sample surface of 2.38 μJ. The probe laser source is a 785 nm, CW diode laser modulated to a repetition rate of 1 kHz with a 20% duty cycle, a power at the sample surface of 4.7 mW for both probe/reference oscillator pairs, and a spot size of 175 μm. The probe laser is RF modulated to the frequency of the pump to reduce sample heating. The silicon avalanche photodetectors used have a 3 dB bandwidth of 1.0 GHz and manual external gain adjustment. Their outputs are concurrently recorded on a dual-band 5 GHz bandwidth digital oscilloscope.

In metallic samples investigated using this technique, IStot(t) has two features of interest, a damped acoustic oscillation overlaid on a non-exponential thermal decay curve.13 A filtered power spectrum of IStot(t) can be used to identify the dominant frequencies of the surface acoustic wave (SAW) oscillations induced by the imposed transient grating. These power spectra can be used as one measure of the quality of the recorded signal. To determine a representative minimum collection time, a series of measurements is carried out on a {001} oriented sample of 99.99% pure, single crystal tungsten, varying the number of traces collected in a single measurement from 10 to 50 000. These measurements are carried out using an imposed grating wavelength of Λ = 5.50 μm. The time-dependent signals for a 250 trace averaged measurement are shown in Fig. 2. The long-timescale oscillation in IS0(t) and IS180(t) is the systematic noise due to the photodetector impulse response present in the measurement system which is removed in the composite signal, IStot(t). The filtered power spectra for each measurement are given in Fig. 3. Given the repetition rate of the pump laser, the sampling rate of the oscilloscope, and the communication time between the oscilloscope and the acquisition computer, the total collection time of these measurements ranged from 0.550 s to 743 s (12.5 min). As the collection time is increased, the signal-to-noise ratio of the identifiable acoustic peaks in the spectrum is increased, as expected.

FIG. 2.

Example traces of IS0(t) (red), IS180(t) (blue), and IStot(t) (black) from DH-TGS measurements on {001} oriented single crystal tungsten using Λ = 5.50 μm and a collection time of 4.19 s.

FIG. 2.

Example traces of IS0(t) (red), IS180(t) (blue), and IStot(t) (black) from DH-TGS measurements on {001} oriented single crystal tungsten using Λ = 5.50 μm and a collection time of 4.19 s.

Close modal
FIG. 3.

Filtered power spectra of DH-TGS measurements collected on {001} oriented single crystal tungsten using Λ = 5.50 μm with collection times varying from 0.550 s to 743 s. Each spectrum is a representative selection from a series of three measurements at a fixed collection time.

FIG. 3.

Filtered power spectra of DH-TGS measurements collected on {001} oriented single crystal tungsten using Λ = 5.50 μm with collection times varying from 0.550 s to 743 s. Each spectrum is a representative selection from a series of three measurements at a fixed collection time.

Close modal

Table I shows statistics for each collection time, including the signal-to-noise ratio and the 95% confidence interval of the peak position parameter from a Gaussian fit to the spectrum. As there is often a DC component retained in the power spectrum from slow thermal decay, the signal-to-noise ratio is calculated over the noise floor from 0.7 to 1.3 GHz as the DC component is clearly separable from the acoustic peak of interest. The entries in Table I are the average values given from three measurements of the same record length. For measurements averaging 10 and 25 laser pulses, identifiable peaks in the power spectrum were not always apparent above the noise floor. From this series of measurements, based solely on the consistency of the appearance of peaks in the acoustic spectrum, a minimum temporal resolution of between 1 and 10 s is observed.

TABLE I.

Characteristics of power spectra from IStot(t) at varying collection times on single crystal tungsten. The signal-to-noise ratio is calculated as SNR=10log10(Ps/Pn), where Ps is the peak signal power and Pn is the maximum noise power in the range of 0.7 to 1.3 GHz. The percent peak position error is the Gaussian fit error divided by the peak position value.

Collection time (s)  0.550  0.770  1.14  1.85  4.19  7.92  15.5  38.5  75.7  149  374  743 
Traces averaged  10  25  50  100  250  500  1000  2500  5000  10 000  25 000  50 000 
Signal-to-noise ratio (dB)  …  …  3.73  6.55  9.97  11.8  14.1  15.0  20.7  22.9  24.8  26.3 
Peak position error (×105 Hz)  …  …  5.3  3.7  3.0  2.5  2.2  2.1  1.7  1.8  1.6  1.6 
% Peak position error  …  …  0.109  0.076  0.061  0.051  0.043  0.042  0.035  0.037  0.034  0.034 
Collection time (s)  0.550  0.770  1.14  1.85  4.19  7.92  15.5  38.5  75.7  149  374  743 
Traces averaged  10  25  50  100  250  500  1000  2500  5000  10 000  25 000  50 000 
Signal-to-noise ratio (dB)  …  …  3.73  6.55  9.97  11.8  14.1  15.0  20.7  22.9  24.8  26.3 
Peak position error (×105 Hz)  …  …  5.3  3.7  3.0  2.5  2.2  2.1  1.7  1.8  1.6  1.6 
% Peak position error  …  …  0.109  0.076  0.061  0.051  0.043  0.042  0.035  0.037  0.034  0.034 

Using this implementation of a TGS experiment, the thermal and acoustic properties of specimens subject to some external driving force can be monitored. Drivers such as temperature, voltage, or irradiation exposure can cause drastic changes in the acoustic and thermal transport properties of materials under investigation. As a simple proof of the principle experiment, a series of DH-TGS measurements is carried out on a {111} oriented sample of single crystal aluminum subject to a temperature ramp between 50 and 250 °C. Measurements were taken with Λ = 4.78 μm, using 1000 trace collections on a 20 s interval over a period of about 30 min. If acoustic peaks are identified in the power spectrum at dominant frequencies fTG, then the surface acoustic wave (SAW) speed for each measurement can be calculated as cTG = ΛfTG since the imposed grating wavelength is calibrated before each experiment, as in Ref. 16. Figure 4 shows the variation in SAW speed with temperature as measured on single crystal aluminum. The apparent linear decrease in cTG over this temperature range is consistent with SAW speeds calculated from temperature-dependent elastic constants of aluminum.15,17

FIG. 4.

Change in SAW speed, cTG, as a function of temperature for {111} oriented single crystal aluminum. Each data point indicates a collection time of 15.5 s recorded at 20 s intervals. Red circles indicate points taken on the heating cycle, blue triangles indicate the cooling cycle, and the dashed line is the expected SAW speed calculated using the method of Ref. 15.

FIG. 4.

Change in SAW speed, cTG, as a function of temperature for {111} oriented single crystal aluminum. Each data point indicates a collection time of 15.5 s recorded at 20 s intervals. Red circles indicate points taken on the heating cycle, blue triangles indicate the cooling cycle, and the dashed line is the expected SAW speed calculated using the method of Ref. 15.

Close modal

In general, collecting traces of IStot(t) with a signal-to-noise ratio in the time domain sufficiently high enough to accurately fit thermal transport parameters requires a longer collection time than that required to identify the primary excited acoustic oscillations as discussed here. However, the difference in these time scales varies widely between different material systems. For this reason and due to contributions of the reflectivity response in the thermal decay profile, the discussion thus far has been restricted to the minimum time scale necessary to analyze the acoustic component of the DH-TGS measurement.

A comment is required concerning the difference in the angles θ1 and θ2 which control the vertical separation of the two probe paths in Fig. 1(b). In practice, these two angles should be set as close as possible to each other while maintaining the ability to independently control the heterodyne phase of each path. The lenses L1 and L2 used in the current iteration of this experiment are achromatically, but not aspherically, corrected doublets. A small spherical aberration present in the lens system causes a focal length shift between the two probing paths. The overlap volume of the two pump beams is large enough with the spot sizes used here to allow for the simultaneous collection of both IS0(t) and IS180(t) with careful positioning of the sample surface in the beam path direction. Upon initial alignment, the photodetectors must be tuned to compensate for differences in signal amplitude occurring due to lens aberration and differences in intrinsic gain. The response amplitudes can be matched by setting the heterodyne phase to ϕ=0° for both probe paths, taking the real-time Fourier transform of the photodetector response on the acquisition oscilloscope, and manually tuning the photodetector gain such that the SAW peak in the transform of each detector has the same power. Future implementations of this technique should utilize both achromatically and aspherically corrected lens systems for L1 and L2 to reduce the sensitivity to this aberration.

The removal of any long-timescale systematic error by subtracting IS0(t) and IS180(t) to find IStot(t) provides several improvements in the analysis of elastic properties through the measured acoustic response. As the signal intensity measured on the photodetectors is intrinsically shot-noise limited, the method proposed here relying on the upstream splitting of the same probe laser source to create both probing paths should, in theory, not result in an improvement in the overall acoustic signal intensity compared to a single-phase collection using the same probe laser source. However, as the reference oscillator needs to be attenuated to avoid saturating the photodetector, splitting the probe laser source allows for the ratio of the reference oscillator intensity to the diffracted signal intensity, IR/ID(t), to be increased across two detectors compared to the single detector in traditional TGS experiments. Additionally, the power spectra of short time scale collections with traditional TGS experiments which do not correct the systematic component of IStot(t) may be more complicated to analyze in the low frequency region. This complication becomes detrimental for measurements at long acoustic wavelengths or for particularly compliant materials.

In contrast to the acoustic analysis, the determination of the thermal transport properties from IStot(t) relies on model fitting, for which the removal of the systematic error by subtracting IS0(t) and IS180(t) is a necessary step.2,11 Given this requirement and the presence of a systematic component to the captured signals, which is common, DH-TGS will allow for thermal transport properties to be measured on much shorter time scales than traditional experiments. Converting existing TGS experiments to DH-TGS experiments is possible as the optical arrangement described here requires only the addition of a small number of optics and a second photodetector. The true advantage using the methodology described in this letter will be the ability to near-continuously characterize both the thermal and elastic properties of dynamic material systems undergoing some external forcing. In contrast, traditional TGS experiments have largely been limited to static material systems. This static-to-dynamic shift will allow the transient grating method to be used in in-situ applications which have not yet been practical with single heterodyne phase collection.

The authors would like to thank A. Maznev and F. Hofmann for useful discussions and M. Rieth for providing the tungsten single crystal material. This paper is based upon work supported by the NSF Graduate Research Fellowship under Grant No. 1122374, the DOE NNSA Stewardship Science Graduate Fellowship under cooperative agreement No. DE-NA0002135, and the MIT-SUTD International Design Center (IDC). M.P.S. acknowledges funding from the U.S. Nuclear Regulatory Commission's MIT Nuclear Education Faculty Development Program under Grant No. NRC-HQ-84-15-G-0045.

1.
A. A.
Maznev
,
K. A.
Nelson
, and
J. A.
Rogers
, “
Optical heterodyne detection of laser-induced gratings
,”
Opt. Lett.
23
,
1319
1321
(
1998
).
2.
J. A.
Johnson
,
A. A.
Maznev
,
M. T.
Bulsara
,
E. A.
Fitzgerald
,
T. C.
Harman
,
S.
Calawa
,
C. J.
Vineis
,
G.
Turner
, and
K. A.
Nelson
, “
Phase-controlled, heterodyne laser-induced transient grating measurements of thermal transport properties in opaque material
,”
J. Appl. Phys.
111
,
023503
(
2012
).
3.
N.
Boechler
,
J. K.
Eliason
,
A.
Kumar
,
A. A.
Maznev
,
K. A.
Nelson
, and
N.
Fang
, “
Interaction of a contact resonance of microspheres with surface acoustic waves
,”
Phys. Rev. Lett.
111
,
036103
(
2013
).
4.
J. A.
Johnson
,
A. A.
Maznev
,
J.
Cuffe
,
J. K.
Eliason
,
A. J.
Minnich
,
T.
Kehoe
,
C. M. S.
Torres
,
G.
Chen
, and
K. A.
Nelson
, “
Direct measurement of room-temperature nondiffusive thermal transport over micron distances in a silicon membrane
,”
Phys. Rev. Lett.
110
,
025901
(
2013
).
5.
J. A.
Rogers
,
A. A.
Maznev
,
M. J.
Banet
, and
K. A.
Nelson
, “
Optical generation and characterization of acoustic waves in thin films: Fundamentals and applications
,”
Annu. Rev. Mater. Sci.
30
,
117
157
(
2000
).
6.
J.
Janusonis
,
T.
Jansma
,
C. L.
Chang
,
Q.
Liu
,
A.
Gatilova
,
A. M.
Lomonosov
,
V.
Shalagatskyi
,
T.
Pezeril
,
V. V.
Temnov
, and
R. I.
Tobey
, “
Transient grating spectroscopy in magnetic thin films: Simultaneous detection of elastic and magnetic dynamics
,”
Sci. Rep.
6
,
29143
(
2016
).
7.
L.
Yang
,
J. D.
Koralek
,
J.
Orenstein
,
D. R.
Tibbetts
,
J. L.
Reno
, and
M. P.
Lilly
, “
Measurement of electron-hole friction in an n-doped GaAs/AlGaAs quantum well using optical transient grating spectroscopy
,”
Phys. Rev. Lett.
106
,
247401
(
2011
).
8.
N.
Gedik
,
J.
Orenstein
,
R.
Liang
,
D. A.
Bonn
, and
W. N.
Hardy
, “
Diffusion of nonequilibrium quasi-particles in a cuprate superconductor
,”
Science
300
,
1410
1412
(
2003
).
9.
R. J.
Smith
,
W.
Li
,
J.
Coulson
,
M.
Clark
,
M. G.
Somekh
, and
S. D.
Sharples
, “
Spatially resolved acoustic spectroscopy for rapid imaging of material microstructure and grain orientation
,”
Meas. Sci. Technol.
25
,
055902
(
2014
).
10.
F.
Hofmann
,
D.
Nguyen-Manh
,
M. R.
Gilbert
,
C. E.
Beck
,
J. K.
Eliason
,
A. A.
Maznev
,
W.
Liu
,
D. E. J.
Armstrong
,
K. A.
Nelsond
, and
S. L.
Dudarev
, “
Lattice swelling and modulus change in a helium-implanted tungsten alloy: X-ray micro-diffraction, surface acoustic wave measurements, and multiscale modelling
,”
Acta Mater.
89
,
352
363
(
2015
).
11.
F.
Hofmann
,
D. R.
Mason
,
J. K.
Eliason
,
A. A.
Maznev
,
K. A.
Nelson
, and
S. L.
Dudarev
, “
Non-contact measurement of thermal diffusivity in ion-implanted nuclear materials
,”
Sci. Rep.
5
,
16042
(
2015
).
12.
R. A.
Duncan
,
F.
Hofmann
,
A.
Vega-Flick
,
J. K.
Eliason
,
A. A.
Maznev
,
A. G.
Every
, and
K. A.
Nelson
, “
Increase in elastic anisotropy of single crystal tungsten upon he-ion implantation measured with laser-generated surface acoustic waves
,”
Appl. Phys. Let.
109
,
151906
(
2016
).
13.
O. W.
Kading
,
H.
Skurk
,
A. A.
Maznev
, and
E.
Matthias
, “
Transient thermal gratings at surfaces for thermal characterization of bulk materials and thin films
,”
Appl. Phys. A
61
,
253
261
(
1995
).
14.
J.
Fivez
, “
The determination of the elastic constants of isotropic solids by means of transient thermal surface gratings
,”
J. Appl. Phys.
119
,
015301
(
2016
).
15.
D.
Gerlich
and
E. S.
Fisher
, “
The high temperature elastic moduli of aluminum
,”
J. Phys. Chem. Solids
30
,
1197
1205
(
1969
).
16.
C. A.
Dennett
,
P.
Cao
,
S. E.
Ferry
,
A.
Vega-Flick
,
A. A.
Maznev
,
K. A.
Nelson
,
A. G.
Every
, and
M. P.
Short
, “
Bridging the gap to mesoscale radiation materials science with transient grating spectroscopy
,”
Phys. Rev. B
94
,
214106
(
2016
).
17.
A. G.
Every
,
A. A.
Maznev
,
W.
Grill
,
M.
Pluta
,
J. D.
Comins
,
O. B.
Wright
,
O.
Matsuda
,
W.
Sachse
, and
J. P.
Wolfe
, “
Bulk and surface acoustic wave phenomena in crystals: Observation and interpretation
,”
Wave Motion
50
,
1197
1217
(
2013
).