pith. sign in

arxiv: 2512.14167 · v2 · submitted 2025-12-16 · ❄️ cond-mat.mtrl-sci

Study of the acoustic and thermal response of an elastically anisotropic solid to a sub-nanosecond laser pulse in transient grating spectroscopy

Pith reviewed 2026-05-16 22:20 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords transient grating spectroscopyfinite element modelinganisotropic solidsthermoelastic couplingacoustic responsethermal grating decaysub-nanosecond laser pulse
0
0 comments X

The pith

A finite element model of transient grating spectroscopy fully captures thermoelastic coupling and optical detection in anisotropic solids.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper builds a numerical simulation of laser-induced transient grating spectroscopy that tracks both heat flow and mechanical deformation inside opaque anisotropic samples. The model incorporates the full interaction between temperature changes and elastic strains together with the heterodyning optical readout of surface motion. It reproduces the expected directional relaxation of the thermal grating and isolates several acoustic signals that appear only in the first few nanoseconds. These short-time features carry independent information on the material's elastic constants. The resulting tool lets researchers test new experimental configurations on a computer before any laboratory work.

Core claim

The central claim is that a custom two-dimensional finite-element formulation accurately represents the coupled thermal and mechanical fields generated by a sub-nanosecond laser pulse together with the optical heterodyning detection of surface displacement, thereby reproducing both the anisotropic thermal decay and the ultra-transient acoustic responses that appear at sub-nanosecond scales.

What carries the argument

Custom-designed two-dimensional finite elements that embed the complete thermoelastic coupling equations and the heterodyning optical detection scheme for anisotropic media.

If this is right

  • The simulation supplies independent elastic constants from the ultra-transient acoustic signals in addition to thermal diffusivity from the longer-time decay.
  • Researchers can explore the effects of different laser pulse durations, grating periods, or crystal orientations entirely inside the model before conducting experiments.
  • The framework extends quantitative transient grating spectroscopy to materials whose elastic anisotropy prevents simple analytical solutions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same numerical approach could be adapted to study time-dependent or viscoelastic response by adding frequency-dependent material parameters.
  • Coupling the forward model to an optimization routine would allow direct fitting of unknown elastic and thermal constants from measured signals.
  • Extending the mesh to three dimensions would reveal out-of-plane wave components that remain invisible in the current two-dimensional treatment.

Load-bearing premise

The custom finite elements correctly include every relevant thermoelastic interaction and detection step at sub-nanosecond scales without numerical artifacts.

What would settle it

A direct comparison between simulated surface-displacement waveforms and experimental transient grating signals recorded on a calibrated anisotropic crystal, checking whether predicted acoustic frequencies and decay rates match the measured traces within experimental uncertainty.

Figures

Figures reproduced from arXiv: 2512.14167 by 2), (2) Faculty of Nuclear Sciences, Czech Academy of Sciences, Czech Technical University in Prague), Hanu\v{s} Seiner (1) ((1) Institute of Thermomechanics, Jakub Ku\v{s}n\'ir (1, Pavla Stoklasov\'a (1), Petr Sedl\'ak (1), Physical Engineering, Prague, Tom\'a\v{s} Grabec (1).

Figure 1
Figure 1. Figure 1: Simplified schematics of the optical paths in the TGS method (a), illustrating the excited grating at the surface (b), and the correspondence with the computational domain (c). The computational domain is shown with illustrated calculated (largely magnified) displacement, where the shading suggests the y-axis displacement magnitude. displacement decreases, and the system becomes sensitive only to thermoref… view at source ↗
Figure 2
Figure 2. Figure 2: Simulated displacement fields (ux,uy, and uz) for the direction 30◦ off the [001] in Ni(110) at various time steps: reaching the maximum surface temperature (2.2 ns, note that the thermal pulse, Eq. (3), is centered around tP = 2 ns), followed by specific times of the SAW oscillation. Below are the corresponding temperature fields showing the gradual homogenization at the surface and its spread towards the… view at source ↗
Figure 3
Figure 3. Figure 3: (a) Comparison of experimental (blue) and simulated (black) time-domain signals obtained from measurement on Ni(110)[001] at the acoustic wavelength λ = 10 µm calculated for the element size of 25 nm (b) Close-up on the first 10 nanoseconds after the pump pulse. (c) Simulated surface displacement profiles at given time steps marked in (b). 3.3 Frequency-domain signal comparison Although the time-domain sig… view at source ↗
Figure 4
Figure 4. Figure 4: Simulated time-domain signals comparing attenuation (lifetime) of the acoustic oscillations with mesh sizes of 25, 50, and 100 nm shown in black, red, and yellow, respectively. While the acoustic lifetime differs significantly, thermal profile agrees with the quasi-static simulation omitting the acoustic dynamics (denoted in green). 4 Conclusion We developed and validated a thermomechanical finite element … view at source ↗
Figure 5
Figure 5. Figure 5: (a) Frequency spectra of measured (blue) signal and signal obtained with the mesh size of 25 nm (black) and 50 nm (red), obtained for a direction 0 ◦ , 30◦ , 60◦ , and 90◦ off the direction [001]. The amplitudes were adjusted for the best fit of the low-amplitude features. (b) Frequency-angular dispersion maps measured by TGS and obtained by FEM, respectively. Note that the simulated map was obtained with … view at source ↗
read the original abstract

Transient grating spectroscopy (TGS) is a material characterization technique based on laser-induced thermoelastic excitation of thermal and acoustic gratings. On opaque samples, these gratings are dynamic surface displacements that reflect the sample's elastic and thermal properties, enabling both types of parameters to be determined from a single experiment. Here, we develop a detailed finite element model (FEM) of the TGS experiment that fully captures the coupling between the thermal and mechanical fields, as well as the optical detection of surface displacement using a heterodyning approach. Using custom-designed two-dimensional elements, the model is particularly suitable for analyzing TGS measurements on anisotropic media, for which analytical theory is insufficient. The simulation captures not only the anisotropic relaxation of the thermoelastic field but also several acoustic features that arise at very short (ultra-transient) timescales and provide additional information about the elastic properties of the examined material. The model offers new opportunities for the in silico testing of various modifications of TGS experiments and their applications to a broad class of materials.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a two-dimensional finite-element model of transient grating spectroscopy (TGS) on elastically anisotropic solids. Custom elements are used to couple the thermal and mechanical fields and to simulate optical heterodyning detection of surface displacements. The central claim is that the simulation reproduces anisotropic thermoelastic relaxation and additional acoustic signatures at ultra-transient timescales that supply extra elastic-property information unavailable from analytical theory.

Significance. If the numerical implementation is shown to be free of artifacts, the work would supply a practical forward-simulation tool for TGS on anisotropic media, where closed-form solutions are unavailable. It would also enable systematic in-silico exploration of experimental variants and could improve extraction of elastic constants from short-time acoustic features.

major comments (2)
  1. [Numerical model / FEM implementation] The section describing the custom two-dimensional finite elements provides no verification against known analytical TGS solutions in the isotropic limit, no mesh-refinement or time-step convergence studies, and no error estimates for the ultra-transient acoustic modes. Without these checks the reported additional elastic information cannot be distinguished from possible numerical artifacts arising from element formulation, boundary conditions, or integration scheme.
  2. [Results] Results section: the assertion that the ultra-transient acoustic features “provide additional information about the elastic properties” is stated qualitatively but is not supported by any sensitivity analysis, parameter-recovery test, or comparison against independently measured elastic constants.
minor comments (2)
  1. [Figures] Figure captions should explicitly state the material symmetry class, the values of the elastic constants used, and the laser-pulse parameters so that readers can reproduce the limiting cases.
  2. [Abstract / Introduction] The abstract and introduction would benefit from a single sentence clarifying whether the simulations are performed on a specific material or on a generic anisotropic medium.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We agree that the numerical implementation requires explicit verification and that the claims regarding ultra-transient acoustic features need quantitative support. We will revise the manuscript accordingly to address these points.

read point-by-point responses
  1. Referee: [Numerical model / FEM implementation] The section describing the custom two-dimensional finite elements provides no verification against known analytical TGS solutions in the isotropic limit, no mesh-refinement or time-step convergence studies, and no error estimates for the ultra-transient acoustic modes. Without these checks the reported additional elastic information cannot be distinguished from possible numerical artifacts arising from element formulation, boundary conditions, or integration scheme.

    Authors: We agree that verification against analytical solutions and convergence studies are essential. In the revised manuscript we will add a new subsection that (i) compares FEM results directly to known closed-form isotropic TGS solutions for both thermal decay and acoustic oscillations, (ii) presents mesh-refinement and time-step convergence studies with quantitative error metrics, and (iii) provides error estimates specifically for the ultra-transient acoustic modes. These additions will demonstrate that the reported features are not numerical artifacts. revision: yes

  2. Referee: [Results] Results section: the assertion that the ultra-transient acoustic features “provide additional information about the elastic properties” is stated qualitatively but is not supported by any sensitivity analysis, parameter-recovery test, or comparison against independently measured elastic constants.

    Authors: We acknowledge that the current claim is qualitative. In the revision we will add (i) a sensitivity analysis showing how individual elastic constants affect the amplitude, frequency, and decay of the ultra-transient features, (ii) parameter-recovery tests in which elastic constants are extracted from simulated ultra-transient data and compared with the input values, and (iii) where literature data exist, a comparison against independently measured elastic constants for the same material. These quantitative results will substantiate the additional information content. revision: yes

Circularity Check

0 steps flagged

No circularity: forward FEM simulation independent of fitted data

full rationale

The paper develops a custom 2D finite-element forward model for thermoelastic response and optical detection in TGS on anisotropic media. The derivation proceeds from the standard coupled thermoelastic equations and optical heterodyning formulas implemented numerically; no parameter is fitted to the target dataset and then re-used as a 'prediction,' no self-citation supplies a uniqueness theorem or ansatz that the present work relies upon, and no known empirical pattern is merely renamed. The central claim that the model captures ultra-transient acoustic features is therefore a direct numerical consequence of the implemented physics rather than a tautology. Lack of benchmark verification against analytical limits is a verification concern, not a circularity issue.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Based on abstract only: the model rests on standard assumptions of linear thermoelasticity and finite element discretization for coupled fields; no free parameters or invented entities are explicitly introduced in the provided text.

axioms (2)
  • domain assumption Linear thermoelastic coupling holds for the sub-nanosecond regime in anisotropic solids
    Invoked implicitly when stating the model captures thermal and mechanical field coupling.
  • domain assumption Custom 2D finite elements accurately discretize the surface displacement and optical detection
    Central to the claim that the model is suitable for anisotropic media where analytical theory fails.

pith-pipeline@v0.9.0 · 5560 in / 1301 out tokens · 29159 ms · 2026-05-16T22:20:58.319242+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

35 extracted references · 35 canonical work pages

  1. [1]

    R. M. White. Generation of elastic waves by transient surface heating.J. Appl. Phys., 34(12):3559–3567, dec 1963

  2. [2]

    H.J. Eichler. Laser-induced Grating Phenomena.Optica Acta: International Journal of Optics, 24(6):631–642, June 1977

  3. [3]

    Rogers, Martin Fuchs, Matthew J

    John A. Rogers, Martin Fuchs, Matthew J. Banet, John B. Hanselman, Randy Logan, and Keith A. Nel- son. Optical system for rapid materials characterization with the transient grating technique: Application to nondestructive evaluation of thin films used in microelectronics.Applied Physics Letters, 71(2):225–227, July 1997

  4. [4]

    A. A. Maznev, K. A. Nelson, and J. A. Rogers. Optical heterodyne detection of laser-induced gratings.Opt. Lett., 23(16):1319–1321, August 1998. Publisher: OSA

  5. [5]

    Stoklasová, T

    P. Stoklasová, T. Grabec, K. Zoubková, P. Sedlák, S. Krátký, and H. Seiner. Laser-Ultrasonic Characterization of Strongly Anisotropic Materials by Transient Grating Spectroscopy.Experimental Mechanics, February 2021

  6. [6]

    A. A. Maznev, Keith A. Nelson, and T. Yagi. Surface phonon spectroscopy with frequency-domain impulsive stimulated light scattering.Solid State Communications, 100(12):807–811, December 1996

  7. [7]

    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.Applied Physics Letters, 109(15):151906, October 2016

  8. [8]

    Grabec, Z

    T. Grabec, Z. Soudná, K. Repček, K. Lünser, S. Fähler, P. Stoklasová, P. Sedlák, and H. Seiner. Guided acoustic waves in thin epitaxial films: Experiment and inverse problem solution for NiTi.Ultrasonics, 138, 2024

  9. [9]

    O. W. Käding, H. Skurk, A. A. Maznev, and E. Matthias. Transient thermal gratings at surfaces for thermal characterization of bulk materials and thin films.Applied Physics A, 61(3):253–261, September 1995

  10. [10]

    Johnson, A

    Jeremy A. Johnson, A. A. Maznev, John Cuffe, Jeffrey K. Eliason, Austin J. Minnich, Timothy Kehoe, Clivia M. Sotomayor Torres, Gang Chen, and Keith A. Nelson. Direct Measurement of Room-Temperature Non- diffusive Thermal Transport Over Micron Distances in a Silicon Membrane.Phys. Rev. Lett., 110(2):025901, January 2013. Publisher: American Physical Society

  11. [11]

    Dennett and Michael P

    Cody A. Dennett and Michael P. Short. Thermal diffusivity determination using heterodyne phase insensitive transientgratingspectroscopy.Journal of Applied Physics, 123(21):215109, 2018. Publisher: AmericanInstitute of Physics

  12. [12]

    Char- acterizing microscale energy transport in materials with transient grating spectroscopy.Journal of Applied Physics, 130(23):231101, December 2021

    Usama Choudhry, Taeyong Kim, Melanie Adams, Jeewan Ranasinghe, Runqing Yang, and Bolin Liao. Char- acterizing microscale energy transport in materials with transient grating spectroscopy.Journal of Applied Physics, 130(23):231101, December 2021

  13. [13]

    Sermeus, B

    J. Sermeus, B. Verstraeten, R. Salenbien, P. Pobedinskas, K. Haenen, and C. Glorieux. Determination of elastic and thermal properties of a thin nanocrystalline diamond coating using all-optical methods.Thin Solid Films, 590:284–292, September 2015

  14. [14]

    Dennett, Michael P

    Abdallah Reza, Cody A. Dennett, Michael P. Short, John Waite, Yevhen Zayachuk, Christopher M. Magazzeni, Simon Hills, and Felix Hofmann. Non-contact, non-destructive mapping of thermal diffusivity and surface acoustic wave speed using transient grating spectroscopy.Review of Scientific Instruments, 91(5):054902, May

  15. [15]

    Publisher: American Institute of Physics. 11

  16. [16]

    M. J. Simmonds, A. Založnik, M. I. Patino, M. J. Baldwin, and N. Boechler. An increased accuracy laser- induced transient grating spectroscopy analysis method for probing near surface thermal diffusivity with giga- hertz frequency instrumentation.AIP Advances, 14(10):105226, October 2024

  17. [17]

    Apparent anisotropic thermal diffusivity measured in cubic single crystals by transient grating spectroscopy.Journal of Applied Physics, 133(12):125108, March 2023

    Jakub Kušnír, Tomáš Grabec, Kristýna Zoubková, Pavla Stoklasová, Petr Sedlák, and Hanuš Seiner. Apparent anisotropic thermal diffusivity measured in cubic single crystals by transient grating spectroscopy.Journal of Applied Physics, 133(12):125108, March 2023

  18. [18]

    Dennett, Penghui Cao, Sara E

    Cody A. Dennett, Penghui Cao, Sara E. Ferry, Alejandro Vega-Flick, Alexei A. Maznev, Keith A. Nelson, Arthur G. Every, and Michael P. Short. Bridging the gap to mesoscale radiation materials science with transient grating spectroscopy.Phys. Rev. B, 94(21):214106, December 2016. Publisher: American Physical Society

  19. [19]

    Hofmann, D

    F. Hofmann, D. R. Mason, J. K. Eliason, A. A. Maznev, K. A. Nelson, and S. L. Dudarev. Non-Contact Mea- surement of Thermal Diffusivity in Ion-Implanted Nuclear Materials.Scientific Reports, 5(1):16042, November 2015

  20. [20]

    Dennett, Hongbing Yu, Kenichiro Mizohata, and Felix Hofmann

    Abdallah Reza, Guanze He, Cody A. Dennett, Hongbing Yu, Kenichiro Mizohata, and Felix Hofmann. Thermal diffusivity recovery and defect annealing kinetics of self-ion implanted tungsten probed by insitu transient grating spectroscopy.Acta Materialia, 232:117926, June 2022

  21. [21]

    Transient Grating Spectroscopy for Complete Elastic Anisotropy: Beyond the Measurement of Surface Acoustic Waves

    Kristyna Zoubkova, Pavla Stoklasova, Tomas Grabec, Petr Sedlak, and Hanus Seiner. Transient Grating Spectroscopy for Complete Elastic Anisotropy: Beyond the Measurement of Surface Acoustic Waves. In2021 IEEE International Ultrasonics Symposium (IUS), pages 1–3, Xi’an, China, September 2021. IEEE

  22. [22]

    Ultra-transient grating spectroscopy for visualization of surface acoustics, 2025

    Tomáš Grabec, Pavla Stoklasová, Kristýna Repček, Jakub Kušnír, David Mareš, Martin Ševčík, Petr Sedlák, and Hanuš Seiner. Ultra-transient grating spectroscopy for visualization of surface acoustics, 2025

  23. [23]

    Johnson, Alexei A

    Jeremy A. Johnson, Alexei A. Maznev, Mayank T. Bulsara, Eugene A. Fitzgerald, T. C. Harman, S. Calawa, C. J. Vineis, G. Turner, and Keith A. Nelson. Phase-controlled, heterodyne laser-induced transient grating measurements of thermal transport properties in opaque material.Journal of Applied Physics, 111(2):023503, January 2012. Publisher: American Instit...

  24. [24]

    Verstraeten, J

    B. Verstraeten, J. Sermeus, R. Salenbien, J. Fivez, G. Shkerdin, and C. Glorieux. Determination of ther- moelastic material properties by differential heterodyne detection of impulsive stimulated thermal scattering. Photoacoustics, 3(2):64–77, June 2015

  25. [25]

    Forward and inverse problems for surface acoustic waves in anisotropic media: a Ritz-Rayleigh method based approach.Ultrasonics, 56:381–389, Febru- ary 2015

    Pavla Stoklasová, Petr Sedlák, Hanus Seiner, and Michal Landa. Forward and inverse problems for surface acoustic waves in anisotropic media: a Ritz-Rayleigh method based approach.Ultrasonics, 56:381–389, Febru- ary 2015

  26. [26]

    Application of the Ritz–Rayleigh method for Lamb waves in extremely anisotropic media.Wave Motion, 96:102567, July 2020

    Tomáš Grabec, Petr Sedlák, and Hanuš Seiner. Application of the Ritz–Rayleigh method for Lamb waves in extremely anisotropic media.Wave Motion, 96:102567, July 2020

  27. [27]

    Simmonds, Brandon D

    Anže Založnik, Michael J. Simmonds, Brandon D. Schwendeman, Nicholas Boechler, Matthew J. Baldwin, and George R. Tynan. Analytical model for laser-induced transient grating measurements of thermal diffusivity in non-opaque materials.Journal of Applied Physics, 135(12):125108, March 2024

  28. [28]

    Effective non-reflective boundary for Lamb waves: Theory, finite element implementation, and applications.Wave Motion, 58:22–41, November 2015

    Yanfeng Shen and Victor Giurgiutiu. Effective non-reflective boundary for Lamb waves: Theory, finite element implementation, and applications.Wave Motion, 58:22–41, November 2015

  29. [29]

    Turkel, R

    E. Turkel, R. Gordon, and D. Gordon. Local absorbing boundary conditions for the elastic wave equation. Wave Motion, 118:103109, April 2023. 12

  30. [30]

    Courant, K

    R. Courant, K. Friedrichs, and H. Lewy. On the Partial Difference Equations of Mathematical Physics.IBM Journal of Research and Development, 11(2):215–234, March 1967

  31. [31]

    Emsley.The Elements

    J. Emsley.The Elements. Oxford Chemistry Guides. Clarendon Press, 1998

  32. [32]

    Wolfgang S. M. Werner, Kathrin Glantschnig, and Claudia Ambrosch-Draxl. Optical Constants and Inelas- tic Electron-Scattering Data for 17 Elemental Metals.Journal of Physical and Chemical Reference Data, 38(4):1013–1092, December 2009. Publisher: American Institute of Physics

  33. [33]

    Veres, and Martin Ryzy

    Tomáš Grabec, István A. Veres, and Martin Ryzy. Surface acoustic wave attenuation in polycrystals: Numerical modeling using a statistical digital twin of an actual sample.Ultrasonics, 119:106585, February 2022

  34. [34]

    Haslinger, Peter Huthwaite, Peter B

    Georgios Sarris, Stewart G. Haslinger, Peter Huthwaite, Peter B. Nagy, and Michael J. S. Lowe. Attenuation of Rayleigh waves due to surface roughness.J. Acoust. Soc. Am., 149(6):4298–4308, June 2021

  35. [35]

    Kurt J. Marfurt. Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations.Geophysics, 49(5):533–549, May 1984. 13