REVIEW 3 major objections 5 minor 39 references
Mode-selective Raman imaging of metal-organic frameworks reveals surface heterogeneities of single HKUST-1 crystals
T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Raman mapping of a single HKUST-1 crystal shows more spectral variation across its surface than across 100 crystals of the same batch.
desk verdict Useful reference-grade Raman imaging of single HKUST-1 crystals, but the headline within-vs-between-crystal comparison is not yet statistically supported. 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 central mechanism is mode-selective Raman imaging: a confocal Raman microscope raster-scans a diffraction-limited laser spot across the crystal, and the scattered intensity is integrated within narrow windows centered on specific vibrational modes (Cu-O stretches, ring breathing, C-H stretches), producing spatial maps of each mode. The quantitative comparison rests on Lorentzian curve fitting of 100 spectra per condition, from which peak positions and full widths at half maximum are extracted and their standard deviations compared between the one-crystal and the 100-crystal data sets. Density-functional-theory phonon calculations provide the mode assignments that justify which bands carry defect, contaminant, or linker information, and principal-component analysis is applied to the image data cube to expose features not visible in any single band.
What would settle it
Repeat the 100-position within-crystal Raman scan on at least five additional crystals from the same batch and compare the within-crystal standard deviations of peak positions and line widths with the between-crystal standard deviations. If any of the additional crystals shows within-crystal standard deviations not consistently larger than the between-crystal values, the claim that heterogeneity is high within the same crystal fails; a simple F-test on the variances would also settle whether the reported difference is statistically meaningful.
Extended reading notes
Core claim
The core discovery is that spatial heterogeneity within a single HKUST-1 crystal can exceed the heterogeneity measured across a batch of crystals. When the authors average Raman spectra from 100 positions on one crystal and from the centers of 100 crystals of the same batch, the band positions match, but the standard deviations of the peak positions and line widths are larger in the within-crystal set. They interpret this as direct evidence that single-crystal surfaces are spectroscopically nonuniform, with defect- and contaminant-related bands (notably the low-frequency Cu-O modes and the 2918 cm$^{-1}$ band) showing the largest variability. Mode-selective images localize the variation: defect-related bands show intensity enhancements or reductions at the same spots, the contaminant band appears as a connected low-intensity region, and principal-component analysis uncovers a triangular roughly 10 micrometer region that no single band shows. The 3090 cm$^{-1}$ band, previously unreported, is assigned to the out-of-phase C-H stretch of the benzene ring in the trimesic linker, while the 2918 cm$^{-1}$ band is attributed to surface contaminants of unconfirmed identity.
Load-bearing premise
The key result assumes that the single crystal chosen for the 100-position scan is representative of the batch, because the paper compares one crystal's internal spread with the spread of 100 crystal centers and reports no significance test or replicate intra-crystal scans.
Editorial extensions
If this is right
- Within a single HKUST-1 crystal, Raman spectra vary more across the surface than they vary between crystals of the same batch, so single-point measurements underestimate true local heterogeneity.
- The reported peak positions and line widths can serve as a spectral reference for batch-level characterization of HKUST-1 in the lab.
- The previously unreported 3090 cm$^{-1}$ band is assigned to the out-of-phase C-H stretch of the trimesic linker and can be used as an additional fingerprint mode.
- The 2918 cm$^{-1}$ band, assigned to surface contaminants, is absent in some crystals of the same batch, indicating irregular contamination across crystals.
- Principal component analysis of the Raman image data reveals a triangular roughly 10 micrometer surface region that is not visible in any single mode-selective image.
Reading between the lines
- If the within-crystal greater-than-between-crystal variance pattern holds for other MOFs, single-point quality control systematically misses the dominant source of spectral scatter; a multi-point-per-crystal sampling protocol would be the natural fix.
- The triangular region that appears only in PCA components 3-6 may mark a growth sector or facet-dependent defect distribution; correlating Raman PCA maps with atomic-force microscopy on the same crystal could test this.
- If the 2918 cm$^{-1}$ contaminant's chemical identity is confirmed, that band could become a non-destructive, spatially resolved cleanliness indicator for HKUST-1 surfaces.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a Raman micro-spectroscopy study of single HKUST-1 crystals. It combines DFT-based simulations with experimental spectra to assign vibrational modes, including a previously unreported band at 3090 rel.cm−1, and uses mode-selective Raman imaging to visualize spatial heterogeneity across a single crystal. The authors then compare spectral peak positions and line widths obtained from 100 positions on one crystal with those from the center positions of 100 different crystals, concluding that within-crystal spectral variability can exceed between-crystal variability. They also apply PCA to reveal features not obvious in the mode-selective images and make their data and simulation code publicly available.
Significance. If the central claim is supported, the work is a useful contribution: it provides the first spatially resolved, mode-selective Raman images of single HKUST-1 crystals, offers reference values for batch-level spectral variability, and demonstrates that single-point or bulk measurements may underestimate local heterogeneity. The open data and code strengthen reproducibility. However, the headline quantitative claim—that within-crystal variability exceeds between-crystal variability—currently rests on a confounded comparison and lacks statistical testing, so the significance of the paper's main message is not yet established.
major comments (3)
- [Spectral Variability of Single-Crystal HKUST-1, Figure 3 and Table 2] The central comparison between case (a), 100 positions equally distributed across the surface of a single crystal, and case (b), center positions of 100 different crystals, confounds two variables: spatial sampling location and crystal identity. Case (a) includes edges, vertices, and the defect-rich bottom vertex visible in Figure 1b–d, while case (b) samples only central regions. The larger standard deviations in case (a) could therefore arise from position-dependent effects such as surface morphology, focus variations, or edge scattering, rather than from a genuine within-crystal versus between-crystal property. Only one crystal contributes to the entire within-crystal distribution, so the ordering cannot be generalized without replicate crystals. Additionally, the 100 intra-crystal spectra are spatially autocorrelated, making any simple independent-sample test inappropriate. Please provide replicate within-crystal measurements, match the spatial sampling (e.g., compare centers only or use identical grids), and report a significance test such as Levene or a bootstrap that accounts for spatial correlation.
- [Spectral Variability of Single-Crystal HKUST-1, Table 2] The statement that 'spectral standard deviations are larger in case (a) than in case (b)' is not supported by a quantitative table. Table 2 reports peak positions, line widths, and standard deviations only for case (b). Without the corresponding within-crystal standard deviations, the reader cannot evaluate the claim band by band. Please add a table or figure giving the within-crystal standard deviations for the same bands, including confidence intervals or a measure of uncertainty, so the comparison is transparent and reproducible.
- [Mode-selective Raman Imaging, Figure 1f and the assignment in the section 'Single-crystal Raman spectrum of HKUST-1'] The interpretation of Figure 1f as showing 'the presence of contaminants' depends entirely on the assignment of the 2918 rel.cm−1 band to surface contaminants. The authors explicitly note that 'confirming the chemical identity of the adsorbant species would require further research which is beyond the scope of this paper.' This is an acknowledged limitation, yet the later text states that the image 'attests to the presence of contaminants present at the surface of the crystal,' which goes beyond the evidence. Please either soften the language to 'unidentified surface species' or provide independent confirmation, for example control spectra of the bare substrate, comparative measurements on cleaned crystals, or a complementary surface-sensitive technique.
minor comments (5)
- [Mode-selective Raman Imaging, figure cross-reference] The text refers to 'Figures 1 (b–g)' when describing the Raman images, but the figure caption lists only panels (b–e). Correct the cross-reference or extend the figure to include the missing panels.
- [Methods, Spectral Data Analysis] The preprocessing step clips the standard deviation at the 99th percentile to limit cosmic-ray influence. This can bias the reported standard deviations, especially for bands with intrinsically broad distributions. Please specify whether clipping was applied before or after the Lorentzian fitting and discuss its potential effect on the variance metrics that underpin the central comparison.
- [Supporting Information, Equation (10)] The Lorentzian broadening width γ = 5.0 cm−1 used to generate simulated spectra is a free cosmetic parameter. A sentence on how the simulated spectra change with γ, or a justification for the chosen value, would strengthen the DFT-to-experiment comparison.
- [Throughout] There are several typographical errors and style inconsistencies: 'distribtution' should be 'distribution', 'Zeonodo' should be 'Zenodo', 'principle components' should be 'principal components', and 'minization' should be 'minimization'.
- [Table 1] The 2918 rel.cm−1 band is listed as 'This Work' with no simulated mode. Since the assignment to a contaminant is explicitly unconfirmed, consider labeling this band 'unassigned surface species' rather than implying an established origin.
Circularity Check
No significant circularity: the spectral assignments, variance comparison, and PCA are empirical or standard ab-initio procedures, not reductions to inputs.
full rationale
The paper's derivation chain is self-contained. The Raman spectrum is validated by comparing measured bands to DFT-simulated modes, with no frequency fitting; the 3090 rel.cm-1 band is assigned to the out-of-phase CH stretch because a simulated mode appears nearby, which is standard mode-matching rather than circular. The 2918 rel.cm-1 band is explicitly acknowledged as unconfirmed surface contaminants, and the claim about heterogeneity is an empirical comparison of directly measured standard deviations, not a fitted parameter renamed as a prediction. No load-bearing self-citation chain is used: references to prior work are for standard methods or independent assignments. The main limitations—single-crystal representativeness, confounded position sampling, and absence of significance tests—are statistical and external-validity concerns, not circularity. The PCA is applied to the same measured data to reveal features, which is exploratory rather than circular. Overall, no step reduces to its own inputs by construction.
Assumptions & free parameters
free parameters (3)
- Lorentzian broadening width for simulated spectra =
5.0 cm-1
- Standard deviation clipping threshold =
99th percentile
- Number of principal components =
6
assumptions (4)
- domain assumption Harmonic approximation for lattice dynamics and finite-displacement phonons (SI Eq. 1)
- domain assumption PBE-D3 exchange-correlation functional with GTH pseudopotentials as implemented in CP2K
- domain assumption Placzek approximation for Raman intensities (SI Eqs. 4-9)
- domain assumption Measured Raman bands are described by Lorentzian line shapes.
Cite this review
Pith. "Pith review of Mode-selective Raman imaging of metal-organic frameworks reveals surface heterogeneities of single HKUST-1 crystals." pith.science (2026). https://pith.science/paper/3EMX324Q
@misc{pith2026241114644,
author = {Pith},
title = {Pith review of: Mode-selective Raman imaging of metal-organic frameworks reveals surface heterogeneities of single HKUST-1 crystals},
year = {2026},
howpublished = {\url{https://pith.science/paper/3EMX324Q}},
note = {Machine review of arXiv:2411.14644}
}
read the original abstract
Metal organic frameworks (MOFs) are nanoporous materials with high surface-to-volume ratio that have potential applications as gas sorbents. Sample quality is, however, often compromised and it is unclear how defects and surface contaminants affect the spectral properties of single MOF crystals. Raman micro-spectroscopy is a powerful tool for characterizing MOFs, yet spatial spectral heterogeneity distributions of single MOF crystals have not been reported so far. In this work, we use Raman micro-spectroscopy to characterize spatially isolated, single crystals of the MOF species HKUST-1. In a first step, we validate HKUST-1's Raman spectrum based on DFT simulations and we identify a previously unreported vibrational feature. In a second step, we acquire diffraction-limited, mode-selective Raman images of a single HKUST-1 crystals that reveal how the spectral variations are distributed across the crystal surface. In a third step, we statistically analyze the measured spectral peak positions and line widths for quantifying the variability occurring within the same crystal as well as between different crystals taken from the same batch. Finally, we explore how multivariate data analysis can aid feature identification in Raman images of single MOF crystals. For enabling validation and reuse, we have made the spectroscopic data and simulation code publicly available.
Reference graph
Works this paper leans on
- [1]
-
[2]
O.K. Farha, I. Eryazici, N.C. Jeong, B.G. Hauser, C.E. Wilmer, A.A. Sarjeant, R.Q. Snurr, S.T. Nguyen, A.O. Yazaydin, and J.T. Hupp. Metal-organic framework materials with ultra- high surface areas: Is the sky the limit? Journal of the American Chemical Society , 134(36):15016–15021, 2012
work page 2012
-
[3]
H. Furukawa, N. Ko, Y.B. Go, N. Aratani, S.B. Choi, E. Choi, A.O. Yazaydin, R.Q. Snurr, M. O’Keeffe, J. Kim, and O.M. Yaghi. Ultrahigh porosity in metal-organic frameworks. Sci- ence, 329(5990):424–428, 2010
work page 2010
-
[4]
M. Eddaoudi, D.B. Moler, H. Li, B. Chen, T.M. Reineke, M. O’Keeffe, and O.M. Yaghi. Modular chemistry: Secondary building units as a basis for the design of highly porous and robust metal-organic carboxylate frameworks. Accounts of Chemical Research, 34(4):319–330, 2001
work page 2001
-
[5]
F.S. Gentile, M. Pannico, M. Causà, G. Mensitieri, G. Di Palma, G. Scherillo, and P. Musto. Metal defects in hkust-1 mof revealed by vibrational spectroscopy: a combined quantum mechanical and experimental study. Journal of Materials Chemistry A , 8(21):10796–10812, 2020
work page 2020
- [6]
-
[7]
N.R. Dhumal, M.P. Singh, J.A. Anderson, J. Kiefer, and H.J. Kim. Molecular interactions of a cu-based metal-organic framework with a confined imidazolium-based ionic liquid: A combined density functional theory and experimental vibrational spectroscopy study. Journal of Physical Chemistry C , 120(6):3295–3304, 2016
work page 2016
-
[8]
J. Ethiraj, F. Bonino, C. Lamberti, and S. Bordiga. H2s interaction with hkust-1 and zif-8 mofs: A multitechnique study. Microporous and Mesoporous Materials, 207:90–94, 2015
work page 2015
Show all 39 references
-
[9]
Bae, S.H
J. Bae, S.H. Park, D. Moon, and N.C. Jeong. Crystalline hydrogen bonding of water molecules confined in a metal-organic framework. Communications Chemistry, 5(1), 2022
2022
-
[10]
Rivera-Torrente, M
M. Rivera-Torrente, M. Filez, F. Meirer, and B.M. Weckhuysen. Multi-spectroscopic interrogation of the spatial linker distribution in defect-engineered metal–organic frame- work crystals: The [cu3(btc)2x(cydc)x] showcase. Chemistry - A European Journal , 26(16):3614–3625, 2020
2020
-
[11]
Hadjiivanov, D.A
K.I. Hadjiivanov, D.A. Panayotov, M.Y. Mihaylov, E.Z. Ivanova, K.K. Chakarova, S.M. Andonova, and N.L. Drenchev. Power of infrared and raman spectroscopies to characterize metal-organic frameworks and investigate their interaction with guest molecules. Chemical Reviews, 121(3)...
2021
-
[12]
Prestipino, L
C. Prestipino, L. Regli, J.G. Vitillo, F. Bonino, A. Damin, C. Lamberti, A. Zecchina, P.L. Solari, K.O. Kongshaug, and S. Bordiga. Local structure of framework cu(ii) in hkust-1 metallorganic framework: Spectroscopic characterization upon activation and interaction with adsorb...
2006
-
[13]
Kim, W.S
H.K. Kim, W.S. Yun, M.-B. Kim, J.Y. Kim, Y.-S. Bae, J. Lee, and N.C. Jeong. A chemical route to activation of open metal sites in the copper-based metal-organic framework materials hkust-1 and cu-mof-2. Journal of the American Chemical Society , 137(31):10009–10015, 2015
2015
-
[14]
Todaro, A
M. Todaro, A. Alessi, L. Sciortino, S. Agnello, M. Cannas, F.M. Gelardi, and G. Buscarino. Investigation by raman spectroscopy of the decomposition process of hkust-1 upon exposure to air. Journal of Spectroscopy, 2016, 2016
2016
-
[15]
Fuchs, F
A. Fuchs, F. Knechtel, H. Wang, Z. Ji, S. Wuttke, O.M. Yaghi, and E. Ploetz. Water harvesting at the single-crystal level. Journal of the American Chemical Society , 145(26):14324–14334, 2023
2023
-
[16]
Allen, and Antonino Ingargiola
Matthew Newville, Till Stensitzki, Daniel B. Allen, and Antonino Ingargiola. LMFIT: Non-Linear Least-Square Minimization and Curve-Fitting for Python, October 2015
2015
-
[17]
Pedregosa, G
F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, M. Blondel, P. Prettenhofer, R. Weiss, V. Dubourg, J. Vanderplas, A. Passos, D. Cournapeau, M. Brucher, M. Perrot, and E. Duchesnay. Scikit-learn: Machine learning in Python. Journal of Machine Learning...
2011
-
[18]
Generalized gradient approximation made simple
John P Perdew, Kieron Burke, and Matthias Ernzerhof. Generalized gradient approximation made simple. Physical review letters, 77(18):3865, 1996
1996
-
[19]
A consistent and accurate ab initio parametrization of density functional dispersion correction (dft-d) for the 94 elements h-pu
Stefan Grimme, Jens Antony, Stephan Ehrlich, and Helge Krieg. A consistent and accurate ab initio parametrization of density functional dispersion correction (dft-d) for the 94 elements h-pu. The Journal of chemical physics , 132(15), 2010
2010
-
[20]
Effect of the damping function in dispersion corrected density functional theory
Stefan Grimme, Stephan Ehrlich, and Lars Goerigk. Effect of the damping function in dispersion corrected density functional theory. Journal of computational chemistry , 32(7):1456–1465, 2011
2011
-
[21]
Quickstep: Fast and accurate density functional calcula- tions using a mixed gaussian and plane waves approach
Joost VandeVondele, Matthias Krack, Fawzi Mohamed, Michele Parrinello, Thomas Chassaing, and Jürg Hutter. Quickstep: Fast and accurate density functional calcula- tions using a mixed gaussian and plane waves approach. Computer Physics Communications, 167(2):103–128, 2005
2005
-
[22]
cp2k: atomistic simulations of condensed matter systems
Jürg Hutter, Marcella Iannuzzi, Florian Schiffmann, and Joost VandeVondele. cp2k: atomistic simulations of condensed matter systems. Wiley Interdisciplinary Reviews: Computational Molecular Science, 4(1):15–25, 2014
2014
-
[23]
Cp2k: An electronic structure and molecular dynamics software package-quickstep: Efficient and accurate electronic structure calculations
Thomas D Kühne, Marcella Iannuzzi, Mauro Del Ben, Vladimir V Rybkin, Patrick Seewald, Frederick Stein, Teodoro Laino, Rustam Z Khaliullin, Ole Schütt, Florian Schiff- mann, et al. Cp2k: An electronic structure and molecular dynamics software package-quickstep: Efficient and ac...
2020
-
[24]
An efficient orbital transformation method for electronic structure calculations
Joost VandeVondele and Jürg Hutter. An efficient orbital transformation method for electronic structure calculations. The Journal of chemical physics , 118(10):4365– 4369, 2003
2003
-
[25]
Direct energy functional minimization under orthogonality constraints
Valéry Weber, Joost VandeVondele, Jürg Hutter, and Anders Niklasson. Direct energy functional minimization under orthogonality constraints. The Journal of chemical physics, 128(8), 2008
2008
-
[26]
Separable dual-space gaussian pseudopotentials
Stefan Goedecker, Michael Teter, and Jürg Hutter. Separable dual-space gaussian pseudopotentials. Physical Review B, 54(3):1703, 1996
1996
-
[27]
Relativistic separable dual-space gaussian pseudopotentials from h to rn
Christian Hartwigsen, Sephen Gœdecker, and Jürg Hutter. Relativistic separable dual-space gaussian pseudopotentials from h to rn. Physical Review B, 58(7):3641, 1998
1998
-
[28]
Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases
Joost VandeVondele and Jürg Hutter. Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. The Journal of chemical physics , 127(11), 2007
2007
-
[29]
Camp, and David S
Dalar Nazarian, Jeffrey S. Camp, and David S. Sholl. Computation-Ready Experimental Metal-Organic Framework (CoRE MOF) 2014 DDEC Database, August 2020
2014
-
[30]
First-principles phonon calculations with phonopy and phono3py
Atsushi Togo. First-principles phonon calculations with phonopy and phono3py. Journal of the Physical Society of Japan , 92(1):012001, 2023
2023
-
[31]
Implementation strategies in phonopy and phono3py
Atsushi Togo, Laurent Chaput, Terumasa Tadano, and Isao Tanaka. Implementation strategies in phonopy and phono3py. Journal of Physics: Condensed Matter , 2023
2023
-
[32]
Supplementary dataset to mode-selective raman imaging of metal-organic frameworks reveals surface heterogeneities of single hkust-1 crys- tals
Matheus Esteves Ferreira, Mariana Del Grande, Felipe Lopes Oliveira, Rodrigo Neumann Barros Ferreira, Ademir Ferreira da Silva, Pamela Costa Carvalho, Geisa Lima, Ado Jorio, and Mathias Steiner. Supplementary dataset to mode-selective raman imaging of metal-organic frameworks ...
2024
-
[33]
Benefits and complexity of defects in metal-organic frameworks
NS Portillo-Vélez, Juan L Obeso, José Antonio de los Reyes, Ricardo A Peralta, Ilich A Ibarra, and Michael T Huxley. Benefits and complexity of defects in metal-organic frameworks. Communications Materials, 5(1):247, 2024
2024
-
[34]
Defects as color centers: The apparent color of metal–organic frameworks contain- ing cu2+-based paddle-wheel units
Kai Müller, Karin Fink, Ludger Schöttner, Meike Koenig, Lars Heinke, and Christof Wöll. Defects as color centers: The apparent color of metal–organic frameworks contain- ing cu2+-based paddle-wheel units. ACS applied materials & interfaces , 9(42):37463–37467, 2017
2017
-
[35]
Metal–organic frame- works/heterojunction structures for surface-enhanced raman scattering with enhanced sensitivity and tailorability
Wenwen Yuan, Keran Jiao, Hang Yuan, Hongzhao Sun, Eng Gee Lim, Ivona Mitrovic, Sixuan Duan, Shan Cong, Ruiqi Yong, Feifan Li, et al. Metal–organic frame- works/heterojunction structures for surface-enhanced raman scattering with enhanced sensitivity and tailorability. ACS Appl...
2024
-
[36]
H.I. Song, J. Bae, E.J. Lee, K.O. Kirlikovali, O.K. Farha, and N.C. Jeong. Vibrational paddlewheel cu-cu node in metal-organic frameworks: Probe of nonradiative relaxation. Journal of Physical Chemistry C , 124(24):13187–13195, 2020
2020
-
[37]
Karl Pearson. Liii. on lines and planes of closest fit to systems of points in space. The London, Edinburgh, and Dublin philosophical magazine and journal of science , 2(11):559–572, 1901
1901
-
[38]
Analysis of a complex of statistical variables into principal components
Harold Hotelling. Analysis of a complex of statistical variables into principal components. Journal of educational psychology, 24(6):417, 1933
1933
-
[39]
Derek A. Long. The Raman Effect: A Unified Treatment of the Theory of Raman Scattering by Molecules . John Wiley & Sons, 2002. Supporting Information: Mode-selective Raman imaging of metal-organic frameworks reveals surface heterogeneities of single HKUST-1 crystals Matheus Es...
2002
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.