High-Precision Lunar Corner-Cube Retroreflectors: A Wave-Optics Perspective
Pith reviewed 2026-05-22 20:42 UTC · model grok-4.3
The pith
Hollow silicon-carbide corner-cube retroreflectors achieve competitive photon return with nearly an order-of-magnitude mass reduction for lunar laser ranging.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The author establishes that hollow silicon-carbide corner-cube retroreflectors not only match or exceed the photon return of solid fused-silica designs for the same aperture sizes but also provide nearly an order-of-magnitude mass reduction, with particular advantages at 1064 nm where phase errors are reduced, based on a two-dimensional Fourier-optics analysis that incorporates thermal-mechanical errors and velocity aberration.
What carries the argument
The two-dimensional Fourier-optics model that couples aperture size to aberration angular offset through diffraction lobe width and includes wavefront errors from thermal-mechanical sources.
Load-bearing premise
The modeled thermal-mechanical wavefront errors and the two-dimensional Fourier-optics treatment are sufficiently representative of actual three-dimensional lunar surface conditions and velocity-aberration geometry for the 80-110 mm apertures considered.
What would settle it
Measurement of the actual returned photon counts from a hollow silicon-carbide corner-cube retroreflector on the Moon at known velocity aberration offsets compared to the model's predictions.
Figures
read the original abstract
High-precision corner-cube retroreflectors (CCRs) are critical for advanced lunar laser ranging (LLR) because they enable sub-millimeter-scale measurements of the Earth-Moon distance -- a level of precision essential for rigorous tests of relativistic gravitation and for advancing our understanding of lunar geophysics. In this work, we develop a comprehensive two-dimensional Fourier-optics model for single CCRs with apertures ranging from 80-110 mm. Our model incorporates realistic thermal-mechanical wavefront errors, detailed diffraction effects, and velocity aberration offsets. Our analysis reveals a strong coupling between aperture size and aberration angular offset: while larger CCRs deliver high on-axis flux under near-ideal conditions, their narrow diffraction lobes suffer significant flux loss at moderate aberration offsets, thereby favoring smaller apertures with broader main lobes. Furthermore, comparisons between solid fused-silica and hollow silicon-carbide (SiC) CCRs show that hollow designs not only achieve competitive or superior photon return -- particularly at 1064 nm, where phase errors are relatively reduced -- but also offer nearly an order-of-magnitude mass reduction for the same aperture sizes. These results establish a robust quantitative framework for optimizing CCR designs to perform at the sub-millimeter level under realistic lunar conditions and underscore the advantages of precision hollow SiC CCRs for next-generation LLR operations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a two-dimensional Fourier-optics model for single corner-cube retroreflectors (CCRs) with apertures ranging from 80-110 mm. It incorporates realistic thermal-mechanical wavefront errors, diffraction effects, and velocity aberration offsets to compare solid fused-silica and hollow silicon-carbide (SiC) designs. The analysis identifies a coupling between aperture size and aberration angular offset that favors smaller apertures with broader diffraction lobes, and concludes that hollow SiC CCRs achieve competitive or superior photon return (especially at 1064 nm) while providing nearly an order-of-magnitude mass reduction for equivalent apertures.
Significance. If the modeling assumptions hold, the work supplies a quantitative framework for optimizing CCRs for sub-millimeter lunar laser ranging, with the mass-reduction advantage being independent of the optics simulation and therefore robust. The photon-return comparison, however, rests directly on the fidelity of the incorporated wavefront-error spectrum and the 2D diffraction treatment; confirmation of these assumptions would make the hollow-SiC recommendation practically useful for next-generation LLR arrays.
major comments (1)
- [wave-optics model and flux calculations] The central claim that hollow SiC CCRs deliver competitive or superior photon return (particularly at 1064 nm) is load-bearing on the two-dimensional Fourier-optics model and the assumed thermal-mechanical wavefront errors (as described in the abstract and the modeling framework). A 2D slice through the corner-cube geometry cannot capture the full three-dimensional diffraction, polarization-dependent phase shifts on the three orthogonal faces, or the vector nature of velocity aberration (including out-of-plane components), all of which may differ between hollow designs (with air paths and distinct mounting distortions) and solid fused-silica CCRs. If these effects alter the relative flux by more than the reported margin, the optical-performance advantage disappears. A direct comparison to 3D vector simulations or existing LLR return data for similar apertures would be required to support
minor comments (2)
- The abstract states that the model incorporates realistic wavefront errors but supplies neither a full parameter list nor an error budget; adding these (e.g., explicit RMS values, spatial-frequency content, and sensitivity analysis) would allow readers to assess robustness.
- Presentation of the aperture-aberration coupling would be clearer if the main-lobe width versus offset curves were shown for both wavelengths and both materials in a single figure with consistent scaling.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review of our manuscript. The comments on the dimensionality and fidelity of the wave-optics model are well taken. We have revised the manuscript to add an explicit discussion of the 2D approximation's limitations and their bearing on the reported performance margins. Our point-by-point responses follow.
read point-by-point responses
-
Referee: The central claim that hollow SiC CCRs deliver competitive or superior photon return (particularly at 1064 nm) is load-bearing on the two-dimensional Fourier-optics model and the assumed thermal-mechanical wavefront errors (as described in the abstract and the modeling framework). A 2D slice through the corner-cube geometry cannot capture the full three-dimensional diffraction, polarization-dependent phase shifts on the three orthogonal faces, or the vector nature of velocity aberration (including out-of-plane components), all of which may differ between hollow designs (with air paths and distinct mounting distortions) and solid fused-silica CCRs. If these effects alter the relative flux by more than the reported margin, the optical-performance advantage disappears. A direct comparison to 3D vector simulations or existing LLR return data for similar apertures would be required to support
Authors: We agree that a 2D Fourier-optics slice is an approximation and cannot fully replicate three-dimensional vector diffraction, polarization interactions across all three orthogonal faces, or out-of-plane aberration components. Our model is intentionally restricted to the principal plane containing the velocity-aberration vector for lunar ranging, where the dominant flux loss arises from the mismatch between the narrow diffraction lobe of larger apertures and the angular offset. The thermal-mechanical wavefront errors are taken from finite-element analyses performed separately for the solid fused-silica and hollow SiC geometries, thereby incorporating design-specific mounting distortions and air-path effects in the hollow case. Polarization phase shifts upon reflection are largely common to both designs and are folded into the effective wavefront-error spectrum used in the propagation. We have added a new paragraph in Section 2.3 that quantifies the expected magnitude of the neglected out-of-plane and polarization contributions and shows that they remain smaller than the reported performance margin at 1064 nm. While a full 3D vector simulation would provide valuable cross-checks, the present 2D framework already isolates the aperture-size versus aberration coupling that drives the design recommendation and is computationally tractable for the parameter study presented. revision: partial
- A complete 3D vector diffraction simulation or a quantitative comparison against existing LLR return data for 80-110 mm apertures, both of which would require substantial new computational development or access to proprietary flight data beyond the scope of the current study.
Circularity Check
Forward 2D Fourier-optics simulation with external inputs; no circular reduction
full rationale
The paper develops and applies a two-dimensional Fourier-optics model that takes as inputs external physical quantities (thermal-mechanical wavefront errors, diffraction, velocity-aberration offsets) and produces photon-return comparisons between solid and hollow CCRs. No parameter is fitted to the model's own output and then relabeled as a prediction; no uniqueness theorem or ansatz is imported via self-citation; the central claims rest on the forward propagation rather than on any definitional equivalence or self-referential loop. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (2)
- aperture diameter
- wavefront error amplitude and spatial spectrum
axioms (2)
- domain assumption Two-dimensional Fourier-optics propagation adequately captures the three-dimensional diffraction and far-field behavior of the CCR
- domain assumption The chosen thermal-mechanical wavefront error model represents actual lunar surface conditions
Forward citations
Cited by 2 Pith papers
-
High-Precision Amplitude-Modulated Continuous-Wave Lunar Laser Ranging
The paper outlines metrology for AM-CW lunar laser ranging that projects 0.08 mm absolute range precision under favorable conditions with dedicated high-power stations.
-
Fundamental Physics in 2025: Status, Decisive Targets, and Path Forward
The review summarizes the baseline SM+GR+Lambda-CDM framework, lists major anomalies and missing pieces, surveys theoretical and experimental approaches, and outlines a staged roadmap organized by decision points.
Reference graph
Works this paper leans on
-
[1]
62 m/ s at its equator, corresponding to a negligible one-way aberration o f ∼ 0. 015 µ rad. Lunar librations (up to ± 7. 9◦ longitude, ± 6. 7◦ latitude) and orbital inclination (5 . 145◦) introduce slight variations in the effective orientation of lunar CCRs. Although instantaneous aberrations du e to librations are typically below 0 . 1 µ rad, their cumu...
-
[2]
with D = 100 mm. These are central intensities of the laterally shifted Airy pattern (no receiver-apertur e integration). Velocity Aberration Retained Flux α (µ rad) λ = 532 nm λ = 1064 nm 0 100.0% 100.0% 2 69.8% 91.6% 4 20.0% 69.8% 6 0.45% 43.0% 8 1.45% 20.0% One way to mitigate the impact of velocity aberration at 532 nm is to int entionally offset one o...
-
[3]
Solid CCRs (Fused Silica) Wavefront distortions in solid fused-silica CCRs arise from internal t hermal gradients due to lunar diurnal temper- ature fluctuations [ 18, 21, 22]. Although these monolithic prisms provide high thermal and mechanic al stability, they are susceptible to refractive distortions caused by internal temp erature gradients. Over a ful...
-
[4]
Hollow CCRs (Silicon Carbide, SiC) Hollow CCRs do not experience refractive distortions but are more s usceptible to mechanically induced WFEs. Although, hollow CCR designs eliminate bulk refractive distortions but are susceptible to misalignment due to me- chanical flexure, mounting stress, and differential thermal expa nsion effects. Unlike solid CCRs, whi...
-
[5]
Technical Summary Table V summarizes the estimated RMS WFEs for both solid and hollow CCRs at 5 32 nm and 1064 nm. Table VI presents the median achievable WFE values the corresponding Stre hl ratios (7), which quantify the impact of wavefront distortions on the diffraction-limited performance of solid and hollow C CRs. High-quality fused silica CCRs typica...
-
[6]
Aperture Discretization and Grid Resolution Each CCR is modeled as a circular aperture of diameter D, sampled onto a uniform Cartesian grid: x, y ∈ [ − 1 2 D, + 1 2 D ] . (10) The aperture field is discretized using a grid resolution of N × N points, where N ≥ 512 (typically 1024), ensuring accurate representation of both the aperture function and WFE dist...
-
[7]
Wavefront Error Modeling and Phase Distortion WFEs are introduced as a spatially varying phase distortion φ(x, y ), which represents optical path deviations relative to an ideal retroreflecting wavefront. These distortions are dec omposed using a Zernike polynomial expansion: φ(x, y ) = ∑ n,m an,mZn,m(x, y ), (12) where the coefficients an,m are determined t...
-
[8]
Velocity Aberration and Beam Misalignment Velocity aberration arises from the relative motion between the Ear th-based laser station and the lunar CCR array. The Moon’s orbital velocity ( ∼ 1 km/s), Earth’s rotation, and lunar librations induce an angular offs et α between the transmitted and returned beam directions, see Table II. We implement the round-t...
-
[9]
(15) The far-field intensity distribution is given by: Ifar(kx, k y) = ⏐ ⏐Efar(kx, k y) ⏐ ⏐2
Far-Field Diffraction Computation The aperture field is propagated into the far-field domain using the F raunhofer diffraction integral, computed via a two-dimensional fast Fourier transform (FFT): Efar(kx, k y) = F {Eap(x, y )}. (15) The far-field intensity distribution is given by: Ifar(kx, k y) = ⏐ ⏐Efar(kx, k y) ⏐ ⏐2 . (16) Here, kx and ky represent angula...
-
[10]
Scaling of Photon Return Efficiency When the returned beam from a CCR is perfectly aligned ( α = 0), the total photon flux F (D) collected by an Earth-based telescope can be approximated by F (D) ∝ D4 λ 2 ρ Ifar(0, 0) Ifar,ideal(0, 0) , (17) where: D2 represents the geometric collecting area of the CCR; ( D/λ ) 2 captures diffraction-limited beam collimation ...
-
[11]
Flux Normalization and Comparative Analysis To enable direct performance comparison across different CCR con figurations, photon return flux values are nor- malized relative to the maximum observed flux: Fnorm(D, α ) = F (D, α )/ maxD,design[F (D, α )], yielding Fnorm(D, α ) ∝ ( D Dref ) 4 ρ S (λ) [ 2 J1 (πDα/λ ) πDα/λ ]2 . (20) This normalization facilitates...
-
[12]
Zerodur Mirror Bonding Zerodur (lithium aluminosilicate glass-ceramic), with ultra-low therma l expansion ( α CTE = 0 . 05 × 10− 6 K− 1) and moderate stiffness ( E ≈ 90 GPa), demands bonding methods that minimally influence dimensiona l stability and wavefront accuracy: • Optical Contacting: Surfaces polished to < 1 nm RMS roughness achieve molecular-level ...
-
[13]
Fused Silica Mirror Bonding Fused silica, exhibiting moderate stiffness ( E = 72 GPa) and low thermal expansion ( α CTE = 0 . 55 × 10− 6 K− 1), benefits from chemically stable, minimally invasive bonding techniques: • Hydroxide-Catalysis (Silicate) Bonding: Employing sodium silicate solutions, this method achieves chemically stable bond thicknesses around 1...
-
[14]
Silicon Carbide (SiC) Mirror Bonding SiC, characterized by high stiffness ( E ≈ 410 GPa), superior thermal conductivity (120–270 W m − 1K− 1), and moderate thermal expansion ( α CTE = 2 . 2 × 10− 6 K− 1), requires bonding methods capable of withstanding thermal and mechanical stresses: • Active-Metal Brazing (Ti-Ag-Cu alloys): This robust metallurgical bon...
-
[15]
Ultra-Precision Surface Polishing and Metrology Mirror substrates (SiC, fused silica, Zerodur) undergo determinis tic polishing processes utilizing magnetorheological finishing (MRF) or ion-beam figuring (IBF), achieving surface figure accuracy < 1 nm RMS and micro-roughness < 0. 5 nm RMS. Surface quality and wavefront accuracy are verified inte rferometrical...
-
[16]
Precision Dihedral Angle Alignment and Verification CCR optical performance critically depends on maintaining precise dih edral angles (90 ◦ ± 0. 2′′ or ± 1 µ rad). Align- ment is executed via robotic hexapod positioning stages integrated with real-time interferometric feedback. Achieved positional accuracies routinely reach < 0. 5 µ m translationally and ...
-
[17]
Optimized Bonding Techniques and Thermal Management Substrate-appropriate bonding methods are meticulously selecte d to minimize induced wavefront distortion and thermal mismatch stresses: • Optical contacting (Zerodur, fused silica): Ultra-thin, adhesive-free bonding ( < 10 nm), resulting in negligible wavefront distortion ( < 1 nm RMS), bond strengths 0...
-
[18]
Comprehensive Thermal and Mechanical Qualification CCRs must undergo rigorous environmental qualification testing to verify structural and optical robustness under lunar and launch conditions: • Thermal vacuum cycling : > 500 cycles spanning − 170◦C to +120 ◦C, demonstrating WFEs within 5 nm RMS. • Vibration and shock testing : Random vibration tests confor...
-
[19]
Long-term Environmental Stability and Outgassing Contr ol Extended vacuum bake-outs ( >500 hours at pressures < 10− 8 Torr) ensure ultra-low outgassing rates ( < 10− 6 g cm− 2s− 1). Reflective coatings, comprising Ag or Au layers with dielectric prot ective multilayers, undergo radi- ation exposure (proton/electron flux) and accelerated thermal aging tests,...
-
[20]
Quantitative Technical Summary Technical capabilities achieved using current precision assembly and qualification techniques for hollow CCRs are quantitatively summarized in Table IX. Collectively, these state-of-the-art assembly practices, mate rial-specific bond- ing techniques, and comprehensive qualification strategies ensure the CCRs achieve the stringe...
-
[21]
Table X summarizes the CTEs for key materials used in hollow CCRs
Thermal Expansion Effects on Optical Stability Differential thermal expansion across bonded mirror assemblies ca n lead to dihedral angle distortions, introducing wavefront errors exceeding 50 nm RMS. Table X summarizes the CTEs for key materials used in hollow CCRs. TABLE X: Coefficient of thermal expansion (CTE) for materials used in hollow CCR mirrors. Mat...
-
[22]
Results are summariz ed in Table XI
Finite Element Analysis (FEA) of Hollow CCR Stability An FEA study was conducted for a 100-mm hollow CCR with SiC mirrors b onded to a titanium support structure to analyze thermal stress-induced flexure. Results are summariz ed in Table XI. TABLE XI: Thermally induced flexure in a 100-mm hollow CCR und er lunar diurnal cycling. Material Combination Max Fle...
-
[23]
Mitigation Strategies To reduce the thermal distortion effects described above, sever al approaches can be employed: • Use Zerodur or fused-silica mounts rather than titanium. These ultra-low expansion ceramics and glasses exhibit minimal CTE mismatch when bonded to SiC or other low-CTE mirror subs trates. Consequently, the overall mount–mirror assembly is...
-
[24]
Aperture Diameter and Velocity Aberration Tolerance Numerical simulations (Secs. V–VI) indicate that CCR apertures below ∼ 80 mm generally fail to gather enough photons for sub-millimeter LLR precision. On the other hand, larger diameters (exceeding ∼ 110 mm) produce relatively narrow main diffraction lobes, making them highly susceptible to modest velocit...
-
[25]
Lunar Thermal Cycle and Wavefront Deformations Fused silica possesses a modest CTE and a relatively low thermal cond uctivity, see Table I. During a diurnal lunar temperature swing of ∼ 290◦C, these properties cause significant gradients within a solid fused- silica prism, giving rise to refractive-index variations that elevate WFEs to about 50– 70 nm RMS ...
-
[26]
Reflectivity: TIR vs. Mirror-Based Reflection Although TIR in a fused-silica prism can provide per-bounce reflectiv ities of ≳ 0. 99, Fresnel losses at the entrance and exit faces typically reduce the overall reflectivity to around ∼ 0. 92–0. 93. In contrast, hollow SiC CCRs rely on three external mirror reflections, each with ρ ≈ 0. 98 when using modern prote...
-
[27]
Mass Budget and Manufacturing Feasibility For a 100 mm-aperture CCR, monolithic fused-silica designs typically w eigh between 1 . 8 and 2 . 4 kg, depending on prism height and mounting details. In contrast, hollow SiC CCRs of c omparable aperture weigh only about
-
[28]
5 kg, representing an approximate 80% mass reduction
4–0. 5 kg, representing an approximate 80% mass reduction. This lighter mass is especially advantageous in lunar robotic missions, where stringent payload constraints heavily influe nce both mission design and capabilities. Fused-silica CCRs can be polished to high optical standards and asse mbled using hydroxide-catalysis bonding or optical contacting. Th...
-
[29]
Quantitative Summary Table XII summarizes key performance metrics for a 100 mm-aperture CCR o perating under the demanding lunar thermal environment of ∼ 290◦C. The hollow SiC CCR consistently exhibits superior thermo-mechanic al stability, as evidenced by minimal wavefront distortion ( ≤ 5 nm RMS), and significantly reduced mass (0 . 4–0. 5 kg). Although ...
-
[30]
Dual-Reflector Configuration and Co-Boresighting The dual-CCR layout serves three goals: (i) improve azimuthal robustness to velocity-aberration/libration by “clocking” the two far-field patterns so at least one strong lobe re mains near boresight; (ii) enable differential ranging across a short, thermally monitored baseline to estimate and remov e lander-ind...
-
[31]
Diffraction-Limited Performance and Acceptance Angles For a CCR aperture of diameter D = 100 mm operating at a wavelength λ = 1064 nm, the far-field diffraction pattern follows the Airy function, with a first-null angular radius giv en by θnull ≈ 1. 22λ/D ≈ 12. 97 µ rad. Velocity aberrations induced by lunar orbital motion ( ∼ 1 km/ s) and Earth’s rotation (2...
-
[32]
H ere the fixed baseline is L = 0
Differential Ranging and Thermal Expansion Compensation In earlier work [ 10] km-scale baselines between widely separated CCRs were considere d primarily to enable the sensitivity to deep lunar interior with the differential LLR paradigm. H ere the fixed baseline is L = 0. 50 m and targets a different goal: in-situ compensation of local thermo-mechanical expa...
-
[33]
Aperture Selection, Wavelength Considerations, and Mas s Benefits Wave-optics simulations indicate that CCR apertures larger than ap proximately 110 mm become increasingly sen- sitive to modest velocity aberration offsets ( ∼ 4–7 µ rad). These larger apertures produce narrower diffraction lobes , leading to substantial reductions (often exceeding 50%) in on-...
work page 2025
-
[34]
J. O. Dickey, P. L. Bender, J. E. Faller, X. X. Newhall, R. R icklefs, J. G. Ries, P. J. Shelus, C. Veillet, A. L. Whipple, J. R. Wiant, et al., Science 265, 482 (1994)
work page 1994
-
[35]
T. W. Murphy, Jr., E. G. Adelberger, J. B. R. Battat, L. N. C arey, C. D. Hoyle, P. LeBlanc, E. L. Michelsen, K. Nordtvedt, A. E. Orin, J. D. Strasburg, et al., PASP 120, 20 (2008)
work page 2008
-
[36]
J. G. Williams, S. G. Turyshev, and D. H. Boggs, IJMPD 18, 1129 (2009)
work page 2009
-
[37]
T. W. Murphy, Jr., E. G. Adelberger, J. B. R. Battat, C. D. H oyle, R. J. McMillan, E. L. Michelsen, R. L. Samad, C. W. Stubbs, and H. E. Swanson, Icarus 208, 31 (2010)
work page 2010
- [38]
- [39]
-
[40]
S. G. Turyshev, J. G. Williams, W. M. Folkner, G. M. Gutt, R . T. Baran, R. C. Hein, R. P. Somawardhana, J. A. Lipa, and S. Wang, Experimental Astronomy 36, 105 (2013)
work page 2013
- [41]
-
[42]
J. J. Degnan, Photonics 10, 1215 (2023)
work page 2023
- [43]
-
[44]
J. G. Williams, S. G. Turyshev, and D. H. Boggs, Phys. Rev . Lett. 93, 261101 (2004)
work page 2004
-
[45]
J. G. Williams and D. H. Boggs, Interoffice Memorandum IOM 335N-20-01, Jet Propulsion Laboratory, California Instit ute of Technology, Pasadena, California (2020)
work page 2020
-
[46]
J. G. Williams, L. Porcelli, S. Dell’Agnello, L. Mauro, M. Muccino, D. G. Currie, D. Wellnitz, C. Wu, D. H. Boggs, and N. H. Johnson, Planet. Sci. J. 4, 89 (2023)
work page 2023
-
[47]
J. E. Faller and E. J. Wampler, Scientific American 223, 38 (1970)
work page 1970
-
[48]
D. A. Arnold, Cross section of the apollo lunar retroreflector arrays (2005), NASA/GSFC Report, URL https://ilrs.gsfc.nasa.gov/docs/apollo_arrays.pdf
work page 2005
-
[49]
I. H. Malitson, JOSA 55, 1205 (1965)
work page 1965
- [50]
-
[51]
S. D. Goodrow and T. W. Murphy, Applied Optics 51, 8793 (2012)
work page 2012
-
[52]
M. Born and E. Wolf, Principles of Optics: Electromagnetic Theory of Propagati on, Interference and Diffraction of Light (Cambridge University Press, 1999), 7th ed
work page 1999
-
[53]
J. W. Goodman, Introduction to Fourier Optics (Roberts and Company Publishers, 2017), 4th ed
work page 2017
-
[54]
T. W. Murphy, Jr., Rep. Progr. Phys. 76, 076901 (2013). 26
work page 2013
-
[55]
T. W. Murphy and S. D. Goodrow, Applied Optics 52, 117 (2013)
work page 2013
-
[56]
J. B. R. Battat, T. W. Murphy, E. G. Adelberger, B. Gilles pie, C. D. Hoyle, R. J. McMillan, E. L. Michelsen, K. Nordtved t, A. E. Orin, C. W. Stubbs, et al., PASP 121, 29 (2009)
work page 2009
- [57]
- [58]
- [59]
-
[60]
S. G. Turyshev, Ann. Rev. Nucl. Part. Sci. 58, 207 (2008)
work page 2008
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.