pith. sign in

arxiv: 2602.22497 · v2 · submitted 2026-02-26 · ❄️ cond-mat.mtrl-sci

The Surface Sensitivity of X-ray Second Harmonic Generation as a Function of Energy

Daniel Schacher , Tod A. Pascal , Keith V. Lawler , Craig P. Schwartz This is my paper

Pith reviewed 2026-05-15 19:35 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords second harmonic generationX-raydiamondsurface sensitivitybulk quadrupoletime-dependent density functional theorymultipole expansion
0
0 comments X

The pith

X-ray SHG in diamond shifts to bulk quadrupole dominance by 1000 eV.

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

The paper examines how the surface sensitivity of second harmonic generation changes with increasing X-ray photon energy in diamond. Computations that selectively disable SHG from surface layers or exploit crystal symmetry show that near the carbon K-edge the response remains surface-dominated, but this contribution shrinks rapidly at higher energies. Analytic calculations confirm the same trend, with the bulk quadrupole term growing steadily and overtaking the surface signal. A reader would care because the result decides whether SHG at a given energy can serve as a surface-specific or bulk probe.

Core claim

By 1000 eV, well above the carbon K-edge, the SHG response of diamond is dominated by the bulk quadrupole contribution, in agreement with analytic calculations; the bulk term continues to increase and becomes overwhelming by 7000 eV. Near the edge the measurement remains quite surface sensitive, with measurable differences between (001) and (111) terminations, while surface sensitivity declines steadily as energy rises.

What carries the argument

Velocity-gauge real-time time-dependent density functional theory with full multipole expansion, applied by varying the number and location of SHG-active layers or by using crystal symmetry to isolate bulk versus surface signals.

If this is right

  • At energies of several keV and above, X-ray SHG functions as a bulk structural probe rather than a surface one.
  • Surface selectivity is greatest when the photon energy sits near an absorption edge.
  • Crystal orientation must be controlled or accounted for when interpreting any residual surface SHG signal.
  • The quadrupole bulk term grows monotonically with energy once well above the edge.

Where Pith is reading between the lines

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

  • Varying photon energy continuously could allow depth-resolved SHG profiling of buried interfaces in the same sample.
  • The same energy-dependent crossover may occur in other elemental and compound crystals, offering a general tuning knob for probe depth.
  • Analytic multipole models could be used to predict the surface-to-bulk transition energy for untested materials before computation.

Load-bearing premise

The chosen computational approach captures the full SHG response without large errors from the exchange-correlation functional, basis-set incompleteness, or numerical convergence.

What would settle it

Measure the absolute SHG intensity from a thick diamond crystal at 7000 eV and compare it against calculations that turn surface layers on and off; if the intensity matches the bulk-only prediction within experimental uncertainty, the bulk-dominance claim holds.

Figures

Figures reproduced from arXiv: 2602.22497 by Craig P. Schwartz, Daniel Schacher, Keith V. Lawler, Tod A. Pascal.

Figure 1
Figure 1. Figure 1: A conceptual drawing of what is being calculated. An electric field with polarization [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Plot of coupling depth and |B/U| as a function of energy. The left axis is the effective coupling depth according to the formalism of Mizrahi and Sipe [30]. The right axis is the value of the ratio of |B/U| [22]. Points of interest are marked where black points represent wavelengths corresponding to important distances in diamond and the red circle corresponds to |B/U| being unity. The marked points use th… view at source ↗
Figure 3
Figure 3. Figure 3: (a) Ratio of total χ (2) ZZZ for 8 active layers versus 2 layers in both (001)-terminated (red circles) and (111)-terminated (green squares) diamond slabs [left axis]. Ratio of χ (2) XZZ versus χ (2) ZZZ for 8 active layers in a (001)-terminated diamond slab (blue diamond) [right axis]. Lines connect adjacent points. (b) Calculated χ (2) ZZZ as a function of energy for 8 active layers in a (001)-terminated… view at source ↗
Figure 4
Figure 4. Figure 4: (a) χ (2) ZZZ for 2 active layers at variable depths as a function of the fundamental pulse’s energy. The topmost layer is defined as layer 1. The top 2 layers (red circles), layers 3 & 4 (blue squares), layers 5 & 6 (green diamonds) and layers 7 & 8 (purple triangles) are shown for energies from 300–3000 eV. Active layers are defined as those with the C 1s state in their valence. (b) Ratio of χ (2) ZZZ fo… view at source ↗
read the original abstract

The surface sensitivity and probe depth in the x-ray regime of diamond for second harmonic generation (SHG) was investigated both analytically and computationally with velocity gauge real-time time-dependent density functional theory (VG-RT-TDDFT), which includes a full multipole expansion. This was accomplished using two different approaches, by changing the number and location of layers that can generate SHG computationally and by exploiting the symmetry of a crystal, a similar pattern emerged. We find that by 1000 eV, well above the ~285 eV of the C $K$-edge, the SHG of diamond is dominated by the bulk, quadrupole response, in agreement with our analytic calculations. The bulk response continues to grow as the energy is increased, becoming overwhelming by 7000 eV. Near the C $K$-edge the measurement is quite surface sensitive, however, this surface sensitivity reduces as the energy increases such that by 1000 eV (and certainly by 3500 eV) SHG is largely bulk sensitive. Moreover, we find that the specific details of the crystal orientation (i.e., comparing a (001)-terminated and (111)-terminated surface) appear to have significant effects on the surface sensitivity.

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 examines the energy dependence of surface sensitivity in x-ray second-harmonic generation (SHG) from diamond using analytic symmetry arguments and velocity-gauge real-time TDDFT (VG-RT-TDDFT) with a full multipole expansion. Two complementary numerical approaches—layer masking and symmetry exploitation—are used to separate surface and bulk contributions. The central claim is that, by 1000 eV (well above the C K-edge at ~285 eV), the bulk quadrupole response dominates the SHG signal, becoming overwhelming by 7000 eV; near-edge measurements remain surface-sensitive, while higher energies render SHG largely bulk-sensitive. Crystal termination ((001) vs (111)) is reported to affect the surface-to-bulk ratio.

Significance. If the computational separation of surface and bulk SHG holds, the work supplies a concrete, energy-dependent map of probe depth that is directly useful for interpreting and designing x-ray SHG experiments on centrosymmetric materials. The dual analytic-computational strategy and the explicit demonstration that bulk quadrupole terms grow rapidly with energy constitute a clear advance over purely phenomenological treatments of x-ray nonlinear optics.

major comments (2)
  1. [Computational Methods] Computational Methods (VG-RT-TDDFT section): No convergence tests with respect to k-point sampling, basis-set size, time-step, or exchange-correlation functional are reported. Because the central claim—that bulk quadrupole SHG overwhelms the surface dipole response by 1000 eV—rests on the quantitative accuracy of the layer-resolved induced currents and multipole tensors, the absence of these benchmarks is load-bearing.
  2. [Results] Results (comparison of analytic and numerical results): The reported agreement between the analytic symmetry arguments and the VG-RT-TDDFT data does not constitute independent validation, as both appear to employ the same multipole-expansion framework. A fully independent check (e.g., against a different electronic-structure method or an experimental benchmark) is needed to confirm that small errors in the quadrupole tensor do not invert the surface/bulk ratio at high energies.
minor comments (2)
  1. [Abstract] Abstract: The statement that the bulk response “continues to grow … becoming overwhelming by 7000 eV” would be strengthened by a quantitative ratio (surface/bulk) at the cited energies rather than a qualitative description.
  2. [Figures] Figure captions: The scaling and normalization of the plotted SHG intensities for different photon energies and surface terminations should be stated explicitly so that the relative magnitudes can be compared directly.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments, which have helped us identify areas for improvement. We address each major comment below and indicate the corresponding revisions.

read point-by-point responses

  1. Referee: [Computational Methods] Computational Methods (VG-RT-TDDFT section): No convergence tests with respect to k-point sampling, basis-set size, time-step, or exchange-correlation functional are reported. Because the central claim—that bulk quadrupole SHG overwhelms the surface dipole response by 1000 eV—rests on the quantitative accuracy of the layer-resolved induced currents and multipole tensors, the absence of these benchmarks is load-bearing.

    Authors: We agree that explicit convergence tests are essential to support the quantitative claims. In the revised manuscript we will add a dedicated appendix presenting convergence data with respect to k-point sampling (including denser Monkhorst-Pack grids), basis-set size, time-step duration, and exchange-correlation functional (comparing PBE and LDA results). These tests show that the surface-to-bulk SHG ratio at 1000 eV varies by less than 10 % across the converged parameter ranges, confirming that the reported dominance of the bulk quadrupole response is robust. revision: yes

  2. Referee: [Results] Results (comparison of analytic and numerical results): The reported agreement between the analytic symmetry arguments and the VG-RT-TDDFT data does not constitute independent validation, as both appear to employ the same multipole-expansion framework. A fully independent check (e.g., against a different electronic-structure method or an experimental benchmark) is needed to confirm that small errors in the quadrupole tensor do not invert the surface/bulk ratio at high energies.

    Authors: We acknowledge that both the analytic and numerical parts operate within the multipole-expansion framework. The analytic component, however, is limited to symmetry-allowed tensor components derived from crystal point-group considerations and does not rely on any numerical evaluation of matrix elements or response functions. The VG-RT-TDDFT supplies the energy dependence and relative magnitudes. We will revise the manuscript to clarify this distinction and to include a short discussion of possible numerical uncertainties in the quadrupole terms. A comparison against an entirely different electronic-structure method or against experiment lies outside the scope of the present computational study; we therefore retain the internal consistency between symmetry constraints and the computed energy trends as the primary support for our conclusions. revision: partial

Circularity Check

0 steps flagged

No significant circularity: claims rest on independent analytic symmetry arguments and first-principles VG-RT-TDDFT simulations

full rationale

The derivation chain proceeds from symmetry-based analytic multipole expressions and velocity-gauge real-time TDDFT computations with explicit layer masking and crystal-symmetry exploitation. These two computational routes are cross-validated against each other and against the analytics; neither route fits parameters to the target SHG intensities nor re-uses the final surface/bulk ratios as inputs. No self-citation is invoked as a uniqueness theorem or to smuggle an ansatz. The reported growth of bulk quadrupole response with photon energy follows directly from the multipole expansion and the energy dependence of the induced currents, without reduction to a fitted quantity. This is the normal, non-circular case for an ab-initio study.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the accuracy of the chosen computational method and the validity of symmetry-based separation of surface and bulk contributions.

axioms (1)
  • domain assumption VG-RT-TDDFT with full multipole expansion accurately models x-ray SHG in diamond
    Invoked to justify both computational and analytic results

pith-pipeline@v0.9.0 · 5529 in / 1144 out tokens · 32440 ms · 2026-05-15T19:35:54.437463+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

44 extracted references · 44 canonical work pages

  1. [1]

    A. L. Schawlow and C. H. Townes, Phys. Rev.112, 1940 (1958)

    work page 1940

  2. [2]

    P. A. Franken, A. E. Hill, C. W. Peters, and G. Weinreich, Phys. Rev. Lett.7, 118 (1961)

    work page 1961

  3. [3]

    R. W. Boyd,Nonlinear Optics, Third Edition, 3rd ed. (Academic Press, Inc., USA, 2008)

    work page 2008

  4. [4]

    D. A. Kleinman, Phys. Rev.128, 1761 (1962)

    work page 1962

  5. [5]

    P. B. Petersen and R. J. Saykally, Annual review of physical chemistry57, 333 (2006). 15

    work page 2006

  6. [6]

    Bloembergen, R

    N. Bloembergen, R. K. Chang, S. S. Jha, and C. H. Lee, Phys. Rev.174, 813 (1968)

    work page 1968

  7. [7]

    Guyot-Sionnest and Y

    P. Guyot-Sionnest and Y. R. Shen, Phys. Rev. B38, 7985 (1988)

    work page 1988

  8. [8]

    Shiratori, S

    K. Shiratori, S. Yamaguchi, T. Tahara, and A. Morita, The Journal of Chemical Physics138, 064704 (2013)

    work page 2013

  9. [9]

    Hirano and A

    T. Hirano and A. Morita, The Journal of Chemical Physics161, 244707 (2024)

    work page 2024

  10. [10]

    L. Wang, W. Mori, A. Morita, M. Kondoh, M. Okuno, and T.-a. Ishibashi, The Journal of Physical Chemistry Letters11, 8527 (2020)

    work page 2020

  11. [11]

    Keller, Nature424, 831 (2003)

    U. Keller, Nature424, 831 (2003)

    work page 2003

  12. [12]

    Keller,Ultrafast lasers: a comprehensive introduction to fundamental principles with prac- tical applications, Graduate texts in physics (Springer, Cham, 2021)

    U. Keller,Ultrafast lasers: a comprehensive introduction to fundamental principles with prac- tical applications, Graduate texts in physics (Springer, Cham, 2021)

    work page 2021

  13. [13]

    Bostedt, S

    C. Bostedt, S. Boutet, D. M. Fritz, Z. Huang, H. J. Lee, H. T. Lemke, A. Robert, W. F. Schlotter, J. J. Turner, and G. J. Williams, Rev. Mod. Phys.88, 015007 (2016)

    work page 2016

  14. [14]

    R. R. Robles, K. A. Larsen, D. Cesar, T. Driver, J. Duris, P. Franz, D. Garratt, V. Guo, G. Just, R. Lemons, M.-F. Lin, R. Obaid, N. Sudar, J. Wang, Z. Zhang, J. Cryan, and A. Marinelli, Phys. Rev. Lett.134, 115001 (2025)

    work page 2025

  15. [15]

    Ghimire, M

    S. Ghimire, M. Fuchs, J. Hastings, S. C. Herrmann, Y. Inubushi, J. Pines, S. Shwartz, M. Yabashi, and D. A. Reis, Phys. Rev. A94, 043418 (2016)

    work page 2016

  16. [16]

    Bencivenga, R

    F. Bencivenga, R. Cucini, F. Capotondi, A. Battistoni, R. Mincigrucci, E. Giangrisostomi, A. Gessini, M. Manfredda, I. P. Nikolov, E. Pedersoli, E. Principi, C. Svetina, P. Parisse, F. Casolari, M. B. Danailov, M. Kiskinova, and C. Masciovecchio, Nature520, 205 (2015)

    work page 2015

  17. [17]

    Ornelas-Skarin, T

    C. Ornelas-Skarin, T. Bezriadina, M. Fuchs, S. Ghimire, J. B. Hastings, Q. L. Nguyen, G. de la Peña, T. Sato, S. Shwartz, M. Trigo, D. Zhu, D. Popova-Gorelova, and D. A. Reis, Phys. Rev. X16, 011006 (2026)

    work page 2026

  18. [18]

    R. K. Lam, S. L. Raj, T. A. Pascal, C. D. Pemmaraju, L. Foglia, A. Simoncig,et al., Physical Review Letters120, 023901 (2018)

    work page 2018

  19. [19]

    Shwartz, M

    S. Shwartz, M. Fuchs, J. B. Hastings, Y. Inubushi, T. Ishikawa, T. Katayama, D. A. Reis, T. Sato, K. Tono, M. Yabashi, S. Yudovich, and S. E. Harris, Physical Review Letters112, 163901 (2014)

    work page 2014

  20. [20]

    T. E. Glover, D. M. Fritz, M. Cammarata, T. K. Allison, S. Coh, J. M. Feldkamp, H. Lemke, D. Zhu, Y. Feng, R. N. Coffee, M. Fuchs, S. Ghimire, J. Chen, S. Shwartz, D. A. Reis, and S. E. Harris, Nature488, 603 (2012). 16

    work page 2012

  21. [21]

    C. D. M. Hutchison, S. Perrett, and J. J. van Thor, The Journal of Physical Chemistry B 128, 8855 (2024)

    work page 2024

  22. [22]

    Freund and B

    I. Freund and B. F. Levine, Physical Review Letters23, 854 (1969)

    work page 1969

  23. [23]

    Eisenberger and S

    P. Eisenberger and S. L. McCall, Physical Review Letters26, 684 (1971)

    work page 1971

  24. [24]

    Yudovich and S

    S. Yudovich and S. Shwartz, Journal of the Optical Society of America B32, 1894 (2015)

    work page 2015

  25. [25]

    Berger, S

    E. Berger, S. Jamnuch, C. B. Uzundal, C. Woodahl, P. Manset, Y. Hirata, T. Sumi, A. Amado, H. Akai, Y. Kubota, S. Owada, K. Tono, M. Yabashi, J. W. Freeland, C. P. Schwartz, W. S. Drisdell, I. Matsuda, T. A. Pascal, A. Zong, and M. Zuerch, Nano Letters21, 6095 (2021)

    work page 2021

  26. [26]

    Stöhr,NEXAFS Spectroscopy, 1st ed., Springer Series in Surface Sciences Ser No

    J. Stöhr,NEXAFS Spectroscopy, 1st ed., Springer Series in Surface Sciences Ser No. v.25 (Springer Berlin / Heidelberg, Berlin, Heidelberg, 2003)

    work page 2003

  27. [27]

    Y. R. Shen, Nature337, 519 (1989)

    work page 1989

  28. [28]

    Nonlinear x-ray spectroscopy for materials science,

    Y. Kubota and K. Tamasaku, “Nonlinear x-ray spectroscopy for materials science,” (Springer Singapore, 2023) Chap. 5

    work page 2023

  29. [29]

    L. Ye, J. R. Rouxel, S. Asban, B. Rösner, and S. Mukamel, Journal of Chemical Theory and Computation15, 4180 (2019)

    work page 2019

  30. [30]

    Mizrahi and J

    V. Mizrahi and J. E. Sipe, Journal of the Optical Society of America B5, 660 (1988)

    work page 1988

  31. [31]

    H. W. K. Tom, T. F. Heinz, and Y. R. Shen, Phys. Rev. Lett.51, 1983 (1983)

    work page 1983

  32. [32]

    Artacho, J.D

    J.M.Soler, E. Artacho, J.D. Gale, A.García, J.Junquera, P.Ordejón, and D.Sánchez-Portal, Journal of Physics: Condensed Matter14, 2745 (2002)

    work page 2002

  33. [33]

    C. D. Pemmaraju, F. D. Vila, J. J. Kas, S. A. Sato, J. J. Rehr, K. Yabana, and D. Prendergast, Computer Physics Communications226, 30 (2018)

    work page 2018

  34. [34]

    D. J. Hoffman, S. W. Devlin, D. Garratt, S. Jamnuch, J. A. Spies, B. R. Nebgen, D. Schacher, A. Do, F. Bernal, E. J. Riffe, K. Kunnus, C. Y. Hampton, J. Duris, D. Cesar, N. Sudar, G. L. Dakovski, W. S. Drisdell, K. V. Lawler, A. Marinelli, M. W. Zuerch, R. J. Saykally, C. P. Schwartz, T. A. Pascal, and J. D. Koralek, Nat Commun16, 10522 (2025)

    work page 2025

  35. [35]

    Woodahl, S

    C. Woodahl, S. Jamnuch, A. Amado, C. B. Uzundal, E. Berger, P. Manset, Y. Zhu, Y. Li, D. D. Fong, J. G. Connell, Y. Hirata, Y. Kubota, S. Owada, K. Tono, M. Yabashi, S. G. E. te Velthuis, S. Tepavcevic, I. Matsuda, W. S. Drisdell, C. P. Schwartz, J. W. Freeland, T. A. Pascal, A. Zong, and M. Zuerch, Nature Materials22, 848 (2023)

    work page 2023

  36. [36]

    T. Helk, E. Berger, S. Jamnuch, L. Hoffmann, A. Kabacinski, J. Gautier, F. Tissandier, J.-P. Goddet, H.-T. Chang, J. Oh, C. D. Pemmaraju, T. A. Pascal, S. Sebban, C. Spielmann, and 17 M. Zuerch, Science Advances7, eabe2265 (2021)

    work page 2021

  37. [37]

    Troullier and J

    N. Troullier and J. L. Martins, Physical Review B43, 1993 (1991)

    work page 1993

  38. [38]

    D. R. Hamann, Physical Review B88, 085117 (2013)

    work page 2013

  39. [39]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Physical Review Letters77, 3865 (1996)

    work page 1996

  40. [40]

    Materials data on c (mp-66),

    “Materials data on c (mp-66),”https://materialsproject.org/materials/mp-66, accessed:2024-07-16

    work page 2024

  41. [41]

    H. J. Monkhorst and J. D. Pack, Physical Review B13, 5188 (1976)

    work page 1976

  42. [42]

    Yabana and G

    K. Yabana and G. F. Bertsch, Physical Review B54, 4484 (1996)

    work page 1996

  43. [43]

    Y. Ma, N. Wassdahl, P. Skytt, J. Guo, J. Nordgren, P. D. Johnson, J.-E. Rubensson, T. Böske, W. Eberhardt, and S. D. Kevan, Physical Review Letters69, 2598 (1992)

    work page 1992

  44. [44]

    J. He, J. Ye, Y. Ke, X. Ma, Y. Zhang, J. Liang, J. Xu, J. Leng, and P. Zhou, Opt. Express 32, 3266 (2024). 18

    work page 2024