REVIEW 4 major objections 5 minor 32 references
Apparatus for the measurement of birefringence maps of optical materials: the case of crystalline silicon for Einstein Telescope
T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read A heterodyne polarimeter can measure optical-path differences below $10^{-12}$ m in silicon, and maps of (100) wafers show position-dependent birefringence of order $10^{-7}$.
desk verdict Solid apparatus paper with believable maps, but the 'intrinsic' (110) value is not stress-secured and needs a caveat. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the Jones-matrix description of the sample as a birefringent Fabry-Perot etalon (Equation 12). Interference among the multiple internal reflections converts part of the per-pass ellipticity $\psi$ into a measurable rotation $\Phi$; measuring both $\Psi$ and $\Phi$ gives the etalon phase $\delta$ through $\Phi/(i\Psi) = 2Z\sin\delta/(1-Z^2)$, and Equation 19 recovers $\psi$ from the transmitted power, so the birefringence $\Delta n$ follows from $\psi = (\pi/\lambda)\Delta D$. In the instrument, ellipticity and rotation are separated by frequency: the rotating half-wave plates move the sample signal to $4\nu_w$, the photoelastic modulator oscillates ellipticity at $\nu_m$, and the Faraday cell oscillates rotation at $\nu_F$, so lock-in demodulation gives $\Psi_0$ and $\Phi_0$ independently.
What would settle it
Map the same (100) wafer at both 1064 and 1550 nm, where $Z$ and the etalon finesse differ by more than a factor of two; if the per-pass $\Delta n$ values extracted from Equations 15-19 disagree by more than the stated uncertainties, the single-etalon inversion with $f=1$ is not correct.
Extended reading notes
Core claim
The paper's core discovery is that crystalline silicon, even when cut on a nominally isotropic (100) plane, is measurably birefringent at the $10^{-7}$ level, with a pattern that varies across the wafer in both magnitude and axis orientation; the same measurement on a (110) face gives an intrinsic birefringence of $(-1.50\pm0.15)\times10^{-6}$ at 1550 nm whose fast axis coincides with $[100]$. The paper further establishes that the polarimeter can separate the ellipticity from the rotation produced by the sample, and that the ratio of the two determines the etalon phase, allowing the per-pass birefringence to be extracted even though interference inside the 1-mm wafer modifies the raw signal by up to tens of percent. On this basis it argues that birefringence of order $10^{-7}$ is present in commercially produced silicon and cannot be ignored when specifying substrates for low-frequency gravitational-wave interferometry.
Load-bearing premise
The analysis assumes that each illuminated spot on the sample is a single uniform birefringent etalon with known interface reflectivity, no dichroism, and full spatial overlap of the interfering beams; if stress varies through the thickness, surfaces are rough, or the beam only partially overlaps the etalon, the extracted $\Delta n$ values are biased even though the raw ellipticity is accurate.
Editorial extensions
If this is right
- If these measurements are representative, next-generation interferometer substrates will need a birefringence budget that includes spatial gradients, not just a mean value, because direction-varying birefringence cannot be suppressed by aligning the input polarization.
- For (110)-oriented silicon, keeping birefringence below $10^{-8}$ requires aligning the crystal's $[100]$ direction to the beam within about 3 mrad, turning crystal-axis orientation into a tolerance requirement.
- The same instrument can map other substrate and coating materials with sensitivity below $10^{-12}$ m optical path, since the method is material-independent once the etalon parameters are known.
- The measured average near $10^{-7}$ is above what the authors consider acceptable for the low-frequency arm of a next-generation detector, so substrate selection and mounting procedures will have to aim at lower-stress material.
Reading between the lines
- A natural extension of the single-etalon model would be to include birefringence that varies along the beam; this could be probed by comparing maps taken with different beam diameters, since a wider beam averages over more of the transverse stress field.
- The gravity-induced birefringence estimate in the paper ($\sim 3.5\times10^{-8}$ for a 96.5-mm wafer) implies that support geometry affects the apparent birefringence of large thin samples; future measurements on mounted test masses should correct for this or suspend the sample at its midline.
- The polarimeter's ability to separate ellipticity from rotation could be turned into a dichroism microscope: any residual rotation signal after the etalon correction would be direct evidence of anisotropic absorption.
- The measured (110) intrinsic birefringence could be used to set a systematic-error floor for any future silicon-based interferometer whose beam path is not perfectly aligned to the $[100]$ axis.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper describes a high-sensitivity polarimeter, based on co-rotating half-wave plates and heterodyne detection, for two-dimensional mapping of optical birefringence, and applies it to crystalline silicon samples for the Einstein Telescope. A Jones-matrix etalon model is used to invert measured ellipticity and rotation into birefringence Δn, with calibration by the Cotton-Mouton effect in air. The authors report position-dependent birefringence of order 10^-7 for (100)-oriented 1-mm silicon wafers at 1064 nm, with average values in Figures 8-10, and an intrinsic (110) birefringence Δn(110)=(-1.50±0.15)×10^-6 at 1550 nm with the fast axis along [100]. Implications for ET substrate specifications, including a 3-mrad alignment tolerance, are discussed.
Significance. If correct, the apparatus achieves an optical-path-difference sensitivity around 10^-12 m and provides spatially resolved birefringence data for commercially available silicon, a candidate ET substrate material. Strengths include an independent Cotton-Mouton calibration in air, use of literature values for n and k rather than free parameters for the (100) samples, and a zoom repeatability check in Figure 9 that supports the reality of the observed spatial patterns. The (110) result, however, is the highest-impact quantitative claim and is also the least supported: the mechanical stress state of the strip is not characterized, and the error budget for the inversion is not given. These gaps currently prevent the paper from fully supporting its central specifications.
major comments (4)
- [§4.2, Eq. (31)] The central claim that Δn(110)=(-1.50±0.15)×10^-6 is the intrinsic spatial-dispersion birefringence of silicon is not supported by the evidence given for mechanical stress. The (110) sample is a 2×4.1×55 mm strip manufactured for channeling, and unlike the (100) samples described in Section 3, no mounting, clamping, or stress state is described. With a stress-optic coefficient C_so≈2×10^-11 Pa^-1 (Ref. [34]), a modest clamping stress of ~1 MPa would produce Δn≈2×10^-5, while the strip's own weight under simply supported conditions gives stresses of order 10^4-10^5 Pa, corresponding to Δn≈10^-7-10^-6, comparable to the quoted uncertainty. Equation (32) estimates gravity effects only for the 96.5-mm wafer, not for this strip. The authors should either characterize the stress state (for example, by finite-element modeling or by measuring the birefringence under different support conditions) or re-state the result as an apparent birefringence that includes mounting-induced contributions.
- [§2.2 and §4.1, Eqs. (15)-(19), (30)] The paper does not provide a propagated uncertainty budget for the inversion from measured Ψ and Φ to Δn. The etalon parameter is measured with the same apparatus as Z=0.128±0.018, i.e., about 14% uncertainty, and Eqs. (15)-(18) are nonlinear in Z. The text itself notes that Ψ can differ from ψ by about ±20% at 1064 nm and ±50% at 1550 nm (Section 2.2), so the Z uncertainty should dominate the final Δn uncertainty. Yet the average values quoted for Figures 8-10 and the value Δn(110)=(-1.50±0.15)×10^-6 in Eq. (31) are presented without an error budget or a covariance analysis of Z, R, f, and the measured Fourier amplitudes. Please provide a per-point uncertainty propagation and state whether the claimed ~10^-7 level and the 10% error bar on Δn(110) survive when the Z uncertainty is included.
- [Figs. 8-10] The birefringence maps show no per-point error bars or confidence intervals. The claim that birefringence is position dependent in both magnitude and axis orientation requires that the observed spatial variation exceed the measurement noise at each point. The zoom repeatability check in Figure 9 supports reproducibility of the pattern, but it does not quantify the per-vector uncertainty. Because the reported values span roughly 0.5-3.0×10^-7 with a step of 2 mm, a statement of typical per-point uncertainty—from lock-in noise, spurious-ellipticity subtraction, and Z propagation—is needed to establish that the map features are not artifacts of the analysis.
- [§4.2] For the (110) sample, the paper states 'In this case we measured Z=0.257 corresponding to f=Z/R=0.84', but no uncertainty or measurement method for this Z is given, and the interference fraction f is not independently verified. The extraction of Δn(110) therefore relies on an assumed value of f that can bias the result if the beam is not fully interfering or if surface roughness varies across the strip. Similarly, the assumption ε=0 (no dichroism) is asserted in Section 3.1 without a direct experimental test. A sensitivity analysis of the extracted Δn to f and ε, or an independent measurement of these parameters, is required to support the quantitative result in Eq. (31).
minor comments (5)
- [§3.1, Eq. (24)] The displayed equation contains an apparent typo: the text reads 'P'(t) = P_out [...] = ≈ P_0 [...]'. Please correct the equality/approximation chain, and check the dimensions of the terms.
- [Fig. 3 caption] The caption contains a typo: 'Elliticity' should be 'Ellipticity'.
- [Table 2] The table legend is partly in Italian ('Legenda'); please use English throughout, and ensure that all column headers are explicit.
- [§4.2] The sign convention leading to negative Δn and the fast axis along [100] is introduced without derivation; a sentence connecting the measured ellipticity/rotation phases to the sign of Δn would help the reader verify the result.
- [§3.2] The claimed ultimate sensitivity S_ΔD≲10^-12 m is derived indirectly as |ψ_spurious|/5; it would be useful to state the assumed noise bandwidth or integration time associated with this sensitivity.
Circularity Check
No circular derivation: silicon birefringence results are measured quantities obtained through independently calibrated optics, with self-citations confined to instrumental formalism.
full rationale
The central result, the silicon birefringence values, is not forced by the paper's inputs. The measured quantities are the ellipticity Ψ and rotation Φ, demodulated at known frequencies (Eq. 25). These are converted to ψ and then Δn using the independent etalon parameters R from literature n and k values (Eqs. 20 and 22), Z from a direct interference-swing measurement (Eq. 30), and the Cotton–Mouton calibration in air (Eqs. 28–29). No free parameter is adjusted to match a target Δn, and no prediction is statistically forced by a fit to a subset of the same silicon data. The self-citations to VMB@CERN and PVLAS ([7], [9], [10], [24], [25]) supply the modulation formalism, spurious-ellipticity analysis, and cavity-etalon language, but the silicon-specific extraction chain in Eqs. (15)–(19) is written out and used with measured inputs; the cited prior work does not itself determine the reported birefringence maps or the Δn(110) value. The stress-state concern about the (110) strip is a validity/systematics issue, not a circularity, because it does not involve a load-bearing step reducing to its own inputs. The derivation is self-contained against the measured data and independent calibration, so no circular step is identified.
Assumptions & free parameters
free parameters (3)
- Etalon parameter Z for (100) 1-mm Si at 1064 nm =
0.128 +/- 0.018
- Etalon parameter Z for (110) 2-mm Si at 1550 nm =
0.257 (f = 0.84)
- Interference fraction f for (100) samples =
1
assumptions (5)
- standard math Jones calculus and the small-angle approximation Delta phi << 1 are valid for these measurements.
- domain assumption A birefringent sample can be modeled as a single etalon with uniform phase difference Delta phi and no coupling between birefringence and interference except through Z.
- domain assumption Silicon samples have no linear dichroism at the measurement wavelengths (epsilon = 0).
- domain assumption The spurious ellipticity measured without the sample can be vectorially subtracted from the measurement with the sample in place.
- domain assumption Literature values of n and k for silicon at 1064 nm and 1550 nm apply to the tested samples at 295 +/- 1 K.
Cite this review
Pith. "Pith review of Apparatus for the measurement of birefringence maps of optical materials: the case of crystalline silicon for Einstein Telescope." pith.science (2026). https://pith.science/paper/K4XYTRNH
@misc{pith2026250609578,
author = {Pith},
title = {Pith review of: Apparatus for the measurement of birefringence maps of optical materials: the case of crystalline silicon for Einstein Telescope},
year = {2026},
howpublished = {\url{https://pith.science/paper/K4XYTRNH}},
note = {Machine review of arXiv:2506.09578}
}
abstract
Einstein Telescope (ET) is expected to achieve sensitivity improvements exceeding an order of magnitude compared to current gravitational-wave detectors. The rigorous characterization in optical birefringence of materials and coatings has become a critical task for next-generation detectors, especially since this birefringence is generally spatially non-uniform. A highly sensitive optical polarimeter has been developed at the Department of Physics and Earth Sciences of the University of Ferrara and INFN - Ferrara Section, Italy, aimed at performing two-dimensional birefringence mapping of substrates. In this paper we describe the design and working principle of the system and present results for crystalline silicon, a candidate material for substrates in the low-frequency (LF) interferometers of ET. We find that the birefringence is of order $10^{-7}$ for commercially available samples and is position dependent in the silicon (100)-oriented samples, with variations in both magnitude and axis orientation. We also measure the intrinsic birefringence of the (110) surface. Implications for the performance of gravitational-wave interferometers are discussed.
Reference graph
Works this paper leans on
-
[34]
He, S., Zheng, T., Danyluk, S.: Analysis and determination of the stress-optic coefficients of thin single crystal silicon samples. J. Appl. Phys.96, 3103 (2004) https://doi.org/10.1063/1.1774259 20 T able 2Review of birefringence measurements of non-stressed mono-crystalline silicon. The values of ∆n(110) Si found in literature have undetermined signs; i...
-
[1]
Rudiger, R.S., Strain, K.A., Danzmann, K.: Birefringence- induced loss in interferometers
Winkler, W., A. Rudiger, R.S., Strain, K.A., Danzmann, K.: Birefringence- induced loss in interferometers. Opt. Commun.112, 245 (1994) https://doi.org/ 10.1016/0030-4018(94)90626-2
-
[2]
Akutsu, T., et al., (KAGRA Collaboration): Overview of KAGRA: Detector design and construction history. Prog. Theor. Exp. Phys.2020, 05–101 (2020) https://doi.org/10.1093/ptep/ptaa125
-
[3]
Hirose, E.,et al.: Characterization of Core Optics in Gravitational-Wave Detec- tors: Case Study of KAGRA Sapphire Mirrors. Phys. Rev. Appl.14, 014021 (2020) https://doi.org/10.1103/PhysRevApplied.14.014021
-
[4]
Zeidler, S., Eisenmann, M., Bazzan, M., Li, P., Leonardi, M.: Correlation between birefringence and absorption mapping in large-size Sapphire substrates for gravitational-wave interferometry. Sci. Rep.13, 10868 (2023) https://doi.org/ 10.1038/s41598-023-45928-0
-
[5]
Lorentz, H.A.: Double refraction by regular crystals. Proc. R. Acad. Amst.24, 22 (1921)
work page 1921
-
[6]
Springer-Verlag, Berlin (1984)
Agranovich, V.M., Ginzburg, V.L.: Crystal Optics with Spatial Dispersion, and Excitons. Springer-Verlag, Berlin (1984). 2nd ed. 17
work page 1984
-
[7]
Zavattini, G., et al., (VMBCERN collaboration): Polarimetry for measuring the vacuum magnetic birefringence with quasi-static fields: a systematics study for the VMB@CERN experiment. Eur. Phys. J. C82, 159 (2022) https://doi.org/ 10.1140/epjc/s10052-022-10100-x
Show all 32 references
-
[8]
Brandi, F., et al., (PVLAS Collaboration): Measurement of the phase anisotropy of very high reflectivity interferential mirrors. Appl. Phys. B65, 351 (1997) https: //doi.org/10.1007/s003400050283
1997 doi
-
[9]
Zavattini, G., et al., (PVLAS Collaboration): On measuring birefringences and dichroisms using Fabry–P´ erot cavities. Appl. Phys. B83, 571 (2006) https://doi. org/10.1007/s00340-006-2189-y
2006 doi
-
[10]
Della Valle, F., et al., (PVLAS Collaboration): The PVLAS experiment: mea- suring vacuum magnetic birefringence and dichroism with a birefringent Fabry- Perot cavity. Eur. Phys. J. C76, 24 (2016) https://doi.org/10.1140/epjc/ s10052-015-3869-8
2016 doi
-
[11]
Phys Rev.77, 287 (1950) https://doi.org/10.1103/PhysRev.77.287
Briggs, H.B.: Optical effects in bulk silicon and germanium. Phys Rev.77, 287 (1950) https://doi.org/10.1103/PhysRev.77.287
1950 doi
-
[12]
Green, M.A., Keevers, M.J.: Optical properties of intrinsic silicon at 300 K. Prog. Photovoltaics Res. Appl.3, 189 (1995) https://doi.org/10.1002/pip.4670030303
1995 doi
-
[13]
Green, M.A.: Self-consistent optical parameters of intrinsic silicon at 300 K including temperature coefficients. Sol. Energy Mater. Sol. Cells92, 1305 (2008) https://doi.org/10.1016/j.solmat.2008.06.009
2008 doi
-
[14]
Schinke, C., Hinken, D., Schmidt, J., Bothe, K., Brendel, R.: Modeling the spec- tral luminescence emission of silicon solar cells and wafers. IEEE J. Photovoltaics 3, 1038 (2013) https://doi.org/10.1109/JPHOTOV.2013.2263985
2013
-
[17]
Franta, D.,et al.: Temperature-dependent dispersion model of float zone crys- talline silicon. Appl. Surf. Sci.421, 405 (2017) https://doi.org/10.1016/j.apsusc. 2017.02.021
2017 doi
-
[18]
Weakliem, H.A., Redfield, D.: Temperature dependence of the optical properties of silicon. J. Appl. Phys.50, 1491 (1979) https://doi.org/10.1063/1.326135 18
1979 doi
-
[19]
Sin, E.H., Ong, C.K., Tan, H.S.: Temperature dependence of interband optical absorption of silicon at 1152, 1064, 750, and 694 nm. Phys. Status Solidi A85, 199 (1984) https://doi.org/10.1002/pssa.2210850124
1984 doi
-
[20]
Nguyen, H.T., Rougieux, F.E., Mitchell, B., Macdonald, D.: Temperature depen- dence of the band-band absorption coefficient in crystalline silicon from pho- toluminescence. J. Appl. Phys.115, 043710 (2014) https://doi.org/10.1063/1. 4862912
2014 doi
-
[21]
Pergamon Press, Oxford (1980)
Born, M., Wolf, E.: Principles of Optics. Pergamon Press, Oxford (1980). 6th ed
1980
-
[22]
Salzberg, C.D., Villa, J.J.: Infrared refractive indexes of silicon germanium and modified selenium glass. J. Opt. Soc. Am.47, 244 (1957) https://doi.org/10. 1364/JOSA.47.000244
1957
-
[23]
Li, H.H.: Refractive index of silicon and germanium and its wavelength and tem- perature derivatives. J. Phys. Chem. Ref. Data9, 561 (1993) https://doi.org/10. 1063/1.555624
1993
-
[24]
Ejlli, A.,et al.: The PVLAS experiment: a 25 year effort to measure vacuum mag- netic birefringence. Phys. Rep.871, 1 (2020) https://doi.org/10.1016/j.physrep. 2020.06.001
2020 doi
-
[25]
Zavattini, G., Della Valle, F., Ejlli, A., Ruoso, G.: A polarisation modulation scheme for measuring vacuum magnetic birefringence with static fields. Eur. Phys. J. C76, 294 (2016) https://doi.org/10.1140/epjc/s10052-016-4139-0 . [Erratum: Eur. Phys. J. C 77, (2017), p. 873, d...
2016 doi
-
[26]
Rizzo, C., Rizzo, A., Bishop, D.M.: The Cotton-Mouton effect in gases: Experi- ment and theory. Int. Rev. Phys. Chem.16, 81 (1997) https://doi.org/10.1080/ 014423597230316
1997
-
[27]
Kyryllin, I.V.,et al.: Crystal assisted steering of muons and antimuons at the muon collider. Nucl. Instrum. Methods Phys. Res. A1075, 170385 (2025) https: //doi.org/10.1016/j.nima.2025.170385
2025
-
[28]
Pastrnak, J., Vedam, K.: Optical anisotropy of silicon single crystals. Phys. Rev. B3, 2567 (1971) https://doi.org/10.1103/PhysRevB.3.2567
1971 doi
-
[29]
Chu, T., Yamada, M., Donecker, J., Rossberg, M., Alex, V., Riemann, H.: Optical anisotropy and strain-induced birefringence in dislocation-free silicon single crys- tals. Mater. Sci. Eng.B91, 174 (2002) https://doi.org/10.1016/S0921-5107(01) 00986-2
2002 doi
-
[30]
Lederhandler, S.R.: Infrared studies of birefringence in silicon. J. Appl. Phys.30, 1631 (1959) https://doi.org/10.1063/1.1735026 19
1959 doi
-
[31]
Fukuzawa, M., Yamada, M.: Photoelastic characterization of si wafers by scanning infrared polariscope. J. Cryst. Growth229, 22 (2001) https://doi.org/10.1016/ S0022-0248(01)01043-0
2001
-
[32]
Classical Quantum Gravity33, 015012 (2016) https://doi.org/10.1088/0264-9381/33/1/ 015012
Kr¨ uger, C.,et al.: Birefringence measurements on crystalline silicon. Classical Quantum Gravity33, 015012 (2016) https://doi.org/10.1088/0264-9381/33/1/ 015012
2016 doi
-
[33]
Hamedan, V.J., Adam, A., Ju, C.B.L., Zhao, C.: Precision mapping of a silicon test mass birefringence. Appl. Phys. Lett.122, 064101 (2023) https://doi.org/10. 1063/5.0136869
2023
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.