REVIEW 4 major objections 4 minor 27 references
Optical characterisation of ilmenite by reflectance spectroscopy
T0 review · 4 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read A wavelength-dependent scattering coefficient is the key to extracting ilmenite's optical constants from reflectance spectra.
desk verdict The paper targets a real gap and ships a plausible data product, but the central validation is circular: the internal scattering coefficient is chosen to match the film emissivity, so the agreement is a fit, not a prediction. 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 argument rests on three linked analyses: a multiple-scattering bidirectional reflectance model for particulate media that ties single-scattering albedo to $n$ and $k$; a two-flux approximation for a plane-parallel slab medium whose directional-hemispherical reflectance depends mainly on $n$; and an electric-field transfer matrix for the three-phase air/pigmented-film/aluminium system used as an independent check on the extracted constants. The pivotal parameter is the internal scattering coefficient $s$, which governs how much of the extinction is assigned to scattering rather than absorption. Holding $s$ constant fails to reproduce absorption beyond 13 µm, whereas assigning $s = 50\,\mu\mathrm{m}^{-1}$ there recovers the observed emissivity.
What would settle it
Measure the directional-hemispherical reflectance of the same ilmenite pellet directly with an integrating sphere and compare it with the converted value used in the paper; a mismatch would propagate directly into the fitted $n$ and cascade into $k$. A second check would be to measure emissivity of films of several thicknesses and test whether the transfer matrix with the reported $n$ and $k$ reproduces all thicknesses.
Extended reading notes
Core claim
The central discovery is that the internal scattering coefficient $s$ in the multiple-scattering reflectance model must be allowed to depend on wavelength: $s = 50\,\mu\mathrm{m}^{-1}$ in the strong-absorption regions ($\lambda < 7\,\mu\mathrm{m}$ and $\lambda > 13\,\mu\mathrm{m}$) and $s = 10^{-7}\,\mu\mathrm{m}^{-1}$ in the high-transmission window (7–13 µm). With this band-wise choice, the combined inversion of pellet and powder reflectance determines $n$ and $k$, and the resulting constants reproduce the independently measured emissivity/absorptivity of a 50 µm ilmenite-pigmented film on aluminium. The paper reports $n$ varying weakly with wavelength and $k$ on the order of $10^{-2}$ in the atmospheric window, rising above $0.1$ beyond 13 µm.
Load-bearing premise
The pellet's directional-hemispherical reflectance is not measured; it is converted from the measured bidirectional reflectance using an empirical relation that ignores the material's volume fraction, and the refractive index $n$ is fit to that converted value, so any error in the relation shifts $n$ and then $k$.
Editorial extensions
If this is right
- For absorbing minerals, the internal scattering coefficient cannot be treated as a constant; small values bias the extracted absorptive index in strong-absorption infrared regions.
- The two-flux pellet measurement supplies an independent route to $n$, removing the need to assume a constant refractive index from external data.
- The extracted $n$ and $k$ enable forward computation of reflectance and emissivity for ilmenite-bearing mixtures and remote sensing spectra.
- A film emissivity measurement provides a validation channel for optical constants of materials that cannot be characterised by ellipsometry.
- The combined procedure generalises to other dark opaque minerals that resist uniform thin-film deposition.
Reading between the lines
- The band-wise $s$ scheme is likely a coarse proxy for a spectrally fine $s(\lambda)$; if higher-resolution transmission data for ilmenite became available, the method could be refined to resolve rapid changes in $k$ inside the strong-absorption regions.
- Because the conversion from bidirectional to directional-hemispherical reflectance is empirical and volume-fraction-free, a direct integrating-sphere measurement of pellet reflectance would provide a sharper test of the absolute scale of $n$ than the film validation alone.
- The same combination of particulate, slab, and film measurements could be applied to mineral powders of planetary regoliths, where only particulate samples are available and direct optical-constant measurements are impossible.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a combined experimental and modeling workflow to retrieve the complex refractive index (n, k) of particulate ilmenite from laboratory reflectance measurements. The approach uses the Hapke bidirectional-reflectance model for powder samples, a two-flux approximation for a pressed pellet, and a transfer-matrix model for an ilmenite-pigmented polymer film on aluminum. The authors claim that a wavelength-dependent internal scattering coefficient, s = 50 µm−1 in strong-absorption regions (<7 µm and >13 µm) and s = 10−7 µm−1 in the 7–13 µm transmission window, is the key to reproducing the measured film emissivity and thereby validating the method. They report a weakly wavelength-dependent n and a k that is about 0.01 in the 7–13 µm window and exceeds 0.1 beyond 13 µm.
Significance. The target problem—determining reliable optical constants of opaque, non-sputterable minerals—is of genuine importance for remote sensing and radiative heat transfer. The two-flux and transfer-matrix framework is physically reasonable in outline, and the paper usefully demonstrates that the assumed value of the internal scattering coefficient strongly affects the retrieved k (Fig. 8). However, the central validation step is circular: the internal scattering coefficient is selected specifically to make the model match the film-emissivity data, and the resulting agreement is then presented as an independent prediction. Consequently, the reported n, k, and s values are not established by the evidence presented. The work offers a promising methodological framework, but not a validated determination of ilmenite's optical constants.
major comments (4)
- [Sec. 4.2, Figs. 12–13] The validation of the wavelength-dependent internal scattering coefficient is circular. The manuscript states that with s1 = 10 µm−1 the emissivity underestimation in the λ > 13 µm region remains, and only when 'increasing the value to s2 = 50 µm−1' is the 'accurate prediction of emissivity/absorptivity' achieved. Since s2 is selected from exactly this emissivity comparison, the subsequent agreement shown in Fig. 13 is a fit, not an independent prediction. This does not establish the claim that a wavelength-dependent s = 50/10−7 µm−1 'paves the way in successfully predicting absorption features.'
- [Sec. 3.3 / Sec. 4.2] The transfer-matrix calculation for the film (Sec. 3.3) uses the retrieved n and k of ilmenite as input; these in turn depend on the chosen s. Hence the film-emissivity agreement is a direct consequence of the same data used to choose s. An independent test would require either a measurement of s, a comparison with literature optical constants of ilmenite, or a prediction for a film with a different pigment loading or thickness. None is provided.
- [Sec. 2.2, Sec. 3.2] The directional–hemispherical reflectance of the pellet is not measured but computed from the bi-directional reflectance via log Rd−d = 1.088 log Rd−h, with an explicit neglect of volume-fraction effects. The refractive index n is then inverted from this converted spectrum (Sec. 3.2), and k depends on n through the Hapke inversion. Any systematic error in that empirical relation—calibrated for other media—propagates directly into both n and k. Without a sensitivity estimate or a verification of the relation for ilmenite pellets, the n and k values are not robust.
- [Sec. 4.1, Fig. 8] No uncertainty or degeneracy analysis is provided for the retrieved n, k, or s. Fig. 8 shows that k is highly sensitive to s in the strong-absorption regions, yet the representative particle sizes (15 µm and 22.5 µm) and the choice to treat s as a step function with only three spectral bands are never varied in a sensitivity study. Consequently confidence intervals for the final optical constants in Fig. 14 are absent, and the possibility of alternative (s, n, k) combinations fitting the same data is not discussed.
minor comments (4)
- [Sec. 2.2] Specify the base of the logarithms in log Rd−d = 1.088 log Rd−h, and report whether the fit is in base 10 or natural log, since the numerical conversion differs.
- [Abstract and Sec. 3.1] Typos: 'electric filed transfer matrix' in the abstract should be 'electric field transfer matrix'; Sec. 3.1 contains 'for for 0–30 µm separate'.
- [Sec. 4.3] The statement that 'k > 0.1 ... demonstrated in the study [10]' appears to cite Hapke (1981), a general theory paper; if the intended reference is Roush et al. (2021) [20] or another ilmenite-specific source, correct the citation.
- [Fig. 11] The y-axis caption reads 'sλ (µm−1)' with values 10^2 to 10^-7; the caption should clarify whether the plotted curve is the step function used in the calculation or a continuous interpolation.
Circularity Check
Circular validation: s=50 µm−1 is chosen to match the film emissivity, then the same agreement is presented as an accurate prediction.
-
fitted input called prediction
[Section 4.2, Figs 12–13]
"In the spectral range of strong absorption ( λ <7 µm and λ >13 µm), values of s = 10, 50, 100, 1000 µm−1 are employed in the calculation for determining the optical properties. ... Increasing the value to s2 = 50 µm−1, this issue is resolved and accurate prediction of emissivity/absorptivity in the region is achieved (see Fig. 13)."
The measured film emissivity (Fig. 3) is the selection target for choosing s2=50 among four trial values. The transfer-matrix calculation (Sec 3.3) uses the optical constants obtained from reflectance inversion with that same s, and Fig. 8 shows that increasing s raises the retrieved k in the strong-absorption bands. Hence the agreement in Fig. 13 is not an independent prediction; it is the criterion that selected s2. Presenting that same agreement as 'accuracy and validity' makes the validation circular: the final s and the strong-absorption k values are forced by the fit to the very data used for validation.
-
fitted input called prediction
[Section 5 (Conclusions)]
"The validation of the combined spectroscopy proposed in this work and the determination of the internal scattering coefficient are achieved by scrutinizing the absorptivity/emissivity of the multi-layer medium of the material, which is modelled as a N -phase stratified medium and analyzed by investigating the propagation of electromagnetic radiation within."
This statement names the film absorptivity/emissivity dataset as the validation instrument, but that same dataset was used in Sec 4.2 to decide between s=10, 50, 100 and 1000 µm−1 and to retain s2=50 because it matched the measured emissivity. Since the match with that dataset is what selected the parameter, citing the match as validation is circular. No independent measurement of s, no withheld spectral region, and no external optical constants are used to test the prediction.
full rationale
The forward-modeling chain (Hapke/Mie for particulate reflectance, two-flux for pellet reflectance, transfer matrix for film emissivity) is internally consistent, and the n/k inversion is a standard iterative fit to the pellet and particulate reflectance data. That part is not circular. The circularity enters with the internal scattering coefficient s, which is the central new parameter of the paper. In Sec 4.2, s=50 µm−1 is selected in the strong-absorption bands specifically because the computed film emissivity matches the measured film emissivity; the same match is then presented in Sec 4.2 and Sec 5 as an 'accurate prediction' and as validation of the combined method. Because Fig. 8 shows that raising s raises the retrieved k, choosing s2=50 is effectively choosing the k that reproduces the film emission data, so the agreement is a fit rather than an independent test. The paper also explicitly states that validation is achieved by scrutinizing the same absorptivity/emissivity dataset used for the selection. No withheld data, external optical constants, or independent measurements of s are offered. The Shkuratov–Grynko conversion in Sec 2.2 is an admitted, non-circular limitation ('The effect of material volume fraction on the relation is neglected'), and the self-citations ([6], [24]) are not load-bearing for the central inversion. The central claim therefore reduces by construction: the 'prediction' of absorption at λ>13 µm and the resulting k>0.1 are consequences of choosing s2=50 to match the film data.
Assumptions & free parameters
free parameters (4)
- Internal scattering coefficient in strong absorption regions (lambda < 7 um and > 13 um) =
50 um^-1 (tested 10, 50, 100, 1000)
- Internal scattering coefficient in high transmission region (7-13 um) =
10^-7 um^-1
- Representative particle size for particulate samples =
15 um (0-30 um sieve) and 22.5 um (0-45 um sieve)
- Band edges for strong/weak absorption switch =
7 um and 13 um
assumptions (6)
- domain assumption Hapke bidirectional reflectance model with internal scattering coefficient s is valid for particulate ilmenite.
- domain assumption The two-flux approximation describes the pellet as a homogeneous, isotropic, optically thick, diffuse slab.
- domain assumption The Shkuratov-Grynko relation converts bi-directional to directional-hemispherical reflectance for the pellet.
- domain assumption Mie theory for spherical particles applies to irregular ilmenite grains.
- domain assumption The 50 um ilmenite-pigmented polyethylene film is a homogeneous stratified medium whose effective optical constants are those of the powder.
- domain assumption Kirchhoff's law holds so emissivity equals absorptivity equals 1 minus reflectance.
Cite this review
Pith. "Pith review of Optical characterisation of ilmenite by reflectance spectroscopy." pith.science (2026). https://pith.science/paper/SARNQIRK
@misc{pith2026241118132,
author = {Pith},
title = {Pith review of: Optical characterisation of ilmenite by reflectance spectroscopy},
year = {2026},
howpublished = {\url{https://pith.science/paper/SARNQIRK}},
note = {Machine review of arXiv:2411.18132}
}
read the original abstract
Bi-directional reflection spectroscopy based on multiple scattering of particulate surfaces is employed in identifying the optical properties of titanium-iron oxide mineral from laboratory reflection measurements.However, the approach suffers from issues including: i) both n and k are to be extracted from a single spectroscopy spectrum, ii) imposing constraints of the K-K correlation relating spectral n and k is weakened by its fundamental insensitivity, and iii) incapability in addressing intrinsic strong absorption features of absorbing materials. We resolve these issues by employing additional optical information of directional-hemispherical reflection and emission/absorption for a slab and a stratified multi-layer medium of the material, respectively. The accompanied analyses consist of radiation transfer in a slab medium investigated using the two-flux approximation method and electromagnetic radiation propagation in a stratified multi-layer medium investigated using the electric filed transfer matrix. We further find that an understanding of the internal scattering coefficient of grains in the multiple scattering model paves the way in successfully predicting absorption features of materials. A wavelength-dependent internal scattering coefficient of the material is then found to be 50 1/micrometer and 1/10000000 1/micrometer in regions of strong absorption and high transmission (between 7 and 13 micrometer), respectively. The value of the refractive index n varies weakly on the wavelength. A pronounced change in the determined absorptive index k with the wavelength is observed. Low values of the absorptive index k on the magnitude of 0.01 are obtained in the transmission window spectral range. In the strong absorption spectral range starting from 13 micrometer, values of the absorptive index k are higher than 0.1.
Figures
Figures from the paper (11 more)
Reference graph
Works this paper leans on
- [1]
-
[2]
K. Bell and P. Reza. Non-contact reflection-mode optical absorption spectroscopy using photoacoustic remote sensing. Optics Letters, 45(13):3427–3430, 2020
work page 2020
-
[3]
M. A. Butt, S. N. Khonina, and N. L. Kazanskiy. Recent advances in photonic crystal optical devices: A review. Optics & Laser Technology, 142:107265, 2021
work page 2021
-
[4]
W. Calvin. Additions and corrections to the absorption coefficients of co 2 Ice: Applications to the Martian south polar cap. Journal of Geophysical Research: Solid Earth, 95(B9):14743–14750, 1990
work page 1990
-
[5]
I. Carusotto, A. A. Houck, A. J. Koll´ ar, P. Roushan, D. I. Schuster, and J. Simon. Photonic materials in circuit quantum electrodynamics. Nature Physics, 16:268–279, 2020
work page 2020
-
[6]
J. Chen, A. Riaz, M. Taheri, A. Kumar, J. Coventry, and W. Lipi´ nski. Optical and radiative characterisa- tion of alumina–silica based ceramic materials for high-temperature solar thermal applications. Journal of Quantitative Spectroscopy & Radiative Transfer, 272:107754, May 2021
work page 2021
-
[7]
R. N. Clark and T. L. Roush. Reflectance spectroscopy: quantitative analysis techniques for remote sensing applications. Journal of Geophysical Research: Solid Earth, 89(B7):6329–6340, 1984
work page 1984
-
[8]
L.A. Dombrovsky, J.H. Randrianalisoa, and D. Baillis. Modified two-flux approximation for identification of radiative properties of absorbing and scattering media from directional-hemispherical measurements. Journal of the Optical Society of America A, 23(1):91–98, January 2006
work page 2006
Show all 27 references
-
[9]
W.N. Hansen. Electric fields produced by the propagation of plane coherent electromagnetic radiation in a stratified medium. Journal of the Optical Society of America, 58(3):380–390, 1968
1968
-
[10]
B. Hapke. Bidirectional reflectance spectroscopy: 1. Theory. Journal of Geophysical Research, 86(B4):3039– 3054, 1981
1981
-
[11]
Tompkins
H.G. Tompkins. A User’s Guide to Ellipsometry. Academic Press, San Diego, 1993
1993
-
[12]
Hiroi and C
T. Hiroi and C. Pieters. Estimation of grain sizes and mixing ratios of fine powder mixtures of common geologic minerals. Journal of Geophysical Research, 99(E5):10867–10879, 1994
1994
-
[13]
Hunt, M.P
J.M. Hunt, M.P. Wisherd, and L.C. Bonham. Infrared absorption spectra of minerals and other inorganic compounds. Analytical Chemistry, 22(12):1478–1497, 1950
1950
-
[14]
Howell, M.P
J.R. Howell, M.P. Meng¨ u¸ c, K. Daun, and R. Siegel. Thermal Radiation Heat Transfer. CRC Press, Boca Raton, 7th edition, 2020
2020
-
[15]
P. Lucey. Model near-infrared optical constantsof olivine and pyroxene as a function of iron content. Journal of Geophysical Research, 103(103):1703–1713, 1998
1998
-
[16]
E. J. Milton. Review article principles of field spectroscopy. Remote Sensing, 8(12):1807–1827, 1987
1987
-
[17]
Modest and S
M.F. Modest and S. Mazumder. Radiative Heat Transfer. Academic Press, New York, 4th edition, 2021
2021
-
[18]
Raman, M.A
A.P. Raman, M.A. Anoma, L. Zhu, E. Rephaeli, and S. Fan. Passive radiative cooling below ambient air temperature under direct sunlight. Nature, 515:540–544, 2014
2014
-
[19]
Roush, J
T. Roush, J. Pollack, F. Witteborn, J. Bregman, and J. Simpson. Ice and minerals on Callisto: A reassessment of the reflectance spectra. Icarus, 86(2):355–382, 1990. 11
1990
-
[20]
Roush, L
T. Roush, L. Teodoro, D. Blewett, and J. Cahill. Optical constants and diffuse reflectance of opaque minerals: A modelin study using magnetite. Icarus, 361:114331, 2021
2021
-
[21]
Roush, F
T.L. Roush, F. Esposito, G. Rossman, and L. Colangeli. Estimated optical constants of gypsum in the regions of weak absorptions: Application of scattering theories and comparisons to independent measurements. Journal of Geophysical Research, 112(B4):E10003, 2007
2007
-
[22]
Shkuratov and Y
Y. Shkuratov and Y. Grynko. Light scattering by media composed of semitransparent particles of different shapes in ray optics approximation: consequences for spectroscopy, photometry, and polarimetry of planetary regoliths. Icarus, 173:16–28, 2005
2005
-
[23]
D. Sun, Y. Zhang, D. Wang, W. Song, X. Liu, J. Pang, D. Geng, Y. Sang, and H. Liu. Microstructure and domain engineering of lithium niobate crystal films for integrated photonic applications. Light: Science & Applications, 9(1):197, 2020
2020
-
[24]
Y. Zhai, Y. Ma, S. David, D. Zhao, R. Lou, G. Tan, R. Yang, and X. Yin. Scalable-manufactured randomized glass-polymer hybrid metamaterial for daytime radiative cooling. Science, 350:1062–1066, 2017. 12 Figures 0–45 µm 0–30 µm λ (µm) Rd− d 20 13 8 4
2017
-
[25]
13 Rd− h Rd− d λ (µm) Rd− d, R d− h 20 13 8 4
05 0 Figure 1: Bi-directional reflectance spectra of particulate ilmenite from RELAB spectral library. 13 Rd− h Rd− d λ (µm) Rd− d, R d− h 20 13 8 4
-
[26]
05 0 Figure 2: Bi-directional reflectance of pellet ilmenite from RELAB spectral library and calculated directional– hemispherical reflectance employing the Shkuratov and Grynko model [22]. 14 Ilmenite-pigmented film λ (µm) ǫλ 20 13 8 4 1 0.8 0.6 0.4 0.2 0 Figure 3: Emissivity/...
-
[27]
22 s4 s3 s2 s1 λ (µm) sλ (µm− 1) 25 20 15 10 5 0 102 10− 1 10− 4 10− 7 Figure 11: Wide band internal scattering coefficient
2 0 Figure 10: Transmission spectra of ilmenite powder deposited on a conventional rock salt window. 22 s4 s3 s2 s1 λ (µm) sλ (µm− 1) 25 20 15 10 5 0 102 10− 1 10− 4 10− 7 Figure 11: Wide band internal scattering coefficient. 23 experimental theoretical λ (µm) ǫλ 20 13 8 4 1 0...
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.