摘要
月面物质中的水可以促进月球地质和环境的演变,并为月球资源利用提供所需的条件。由于月球遥感探测手段存在分辨率低的缺点,使获得水冰的直接证据或确定水冰的赋存形态存在困难。激光拉曼光谱具有无需光照、无需样品预处理且可无损地获得月岩和矿物以及水体的指纹信息的优势,能为月面物质的类型、分布范围甚至含量提供直接有利的信息。本文利用拉曼光谱法对典型月岩/矿物以及吸附水、冰、结晶水和羟基结构水等形态的含水特征进行检测,并对水体含量进行定量分析。首先,利用532 nm激光显微拉曼对月面典型矿物和模拟月壤中多种形态的水进行含水信号识别与特征分析。其次,检验并分析了模拟月壤中吸附水、结晶水和羟基结构水的检测限,及其含量与信号强度间存在的规律。然后,利用岭回归、偏最小二乘回归、线性回归对模拟月壤中三种形态水的不同含量进行定量分析。实验结果表明,能明确识别月壤模拟物中四种形态水的特征谱峰,其中模拟月壤组分和水体的峰位分布区域分别位于100~1 700 c
月球表面的物质成分探测可以帮助人们全面了解月球的化学成分、矿物组成和地质结构,推演月球及太阳系的起源和演化过程,为未来月球人类活动提供必要条件。水/冰作为溶
针对月面是否存在水冰的问题,早在1961年Watson等首次提出月球的永久阴影区可能存在水冰的假
受限于航天器的遥感探测的探测方式、尺度、分辨率等因素,很难提供低储量水冰的直接有利证据;特别是对于水冰数量和浓度较低的环境,这就需要采用质谱、红外光谱和激光光谱等原位检测技术。基于质谱-光谱检测的原位探测手段,能够对物质成分进行更精细的探测、识别和定量,获得更直接有利的证据;目前也将被用于四个探测任务,关注月球水冰的原位探测。包括NASA的Artemis
本文旨在测试拉曼光谱对月壤、典型矿物及混合物中吸附水、结晶水、结构水和冰的光谱特征识别与检测能力。首先,利用532 nm激光显微拉曼光谱对含水月壤模拟物进行检测分析,获得模拟月球材料中不同质量比和不同赋存形态水的原位拉曼光谱特征及检测限,探索拉曼光谱对水冰的检测能力与检测参数的关系。然后,通过研究不同相态和元素键的羟基振动模式建立水含量的定量关系。最后,尝试建立线性、偏最小二乘、支持向量等回归模型对模拟月壤中不同含水浓度进行定量分析。
拉曼光谱是利用激发光照射到样品时,分子与激发光子碰撞后产生的一种非弹性散射,散射光频率较激发光源频率发生变化,即拉曼位移。每种分子或官能团受到激发光源照射产生独特的振动、转动模式,可以通过拉曼谱峰的位置、强度、半高宽、峰形等信息来反映所探测物质的分子结构、组分、性质乃至含量等重要信息。由于拉曼光谱对分子结构敏感,被称为指纹光谱;又因其具有高灵敏度、无损、无需样品制备等优点,成为样品的结构、成分和性质信息探测的有效手段。实际中通常使用激光与样品相互作用,产生拉曼散射,进而使用收集光路及分光器将散射出的光子分解成不同波长的光,光子探测器测量不同波长频率的光的强度,得到拉曼光谱图。根据探测距离与激发收集光路的不同,通常分为显微拉曼和远程拉曼等两种形式。显微拉曼提供了高空间分辨率,可以对样品的局部区域进行高灵敏度的拉曼光谱分析,远程拉曼技术可以实现对远距离或难以接触的样品进行拉曼光谱分析。
由于拉曼谱峰信息与物质组分及结构间存在直接的关联,对于吸附水、结晶水、羟基化合水等物质,主要的拉曼活性官能团是自由羟基以及参与分子间氢键网络的羟基的伸缩振动模式,其振动频率通常在2 900~3 700 c
3 900 c
拉曼光谱是散射光谱,峰强与拉曼散射截面强相关。拉曼散射截面是描述样品在拉曼散射过程中与入射光相互作用的强度的物理量。与溶液或气体相比,面向模拟月球物质含水固体样品的拉曼散射截面通常较小,因为固体样品中的分子受到晶格约束和相互作用的影响,使得分子的振动模式更为复杂。固体样品的拉曼散射截面受到晶体结构、粒径、晶格振动模式、晶体缺陷或杂质、晶体取向、表面光学吸收特性、表面形貌等多种因素的影响。为了降低上述因素对固体样品的拉曼光谱的干扰,提高拉曼信号的检测灵敏度和准确性,可以将含水粉末颗粒样本进行拉曼特征谱峰的解译与识别,然后控制部分变量的含水模拟月壤均匀样本压片,对吸附水、结晶水、羟基水定性识别和水体含量的规律分析,具体是根据不同的分析对象分别观察峰形、提取峰面积(像元去基线数据累加值)、降噪后提取峰高等方法。最终以具有代表性的含水波段的拉曼特征峰来建立回归模型,进而对水体的含量定量反演。不同于尖锐且不重叠的矿物拉曼谱峰,水体拉曼谱峰的特征数量多且可能存在混叠,由此常规基于最小二乘的线性回归(linear regression, LR)可能出现过拟合,而岭回归(ridge regression, RR)和偏最小二乘回归(partial least squares regression, PLSR)可以降低过拟合。其中,岭回归模型是在线性回归模型的损失函数中加入 L2 正则化,加入参数的范数的惩罚项,新的损失函数如
, | (1) |
式中,为岭回归模型的损失函数,为样本数量,为模型对第i个样本的预测值,为第i个样本的真实标签,为正则化参数,用于控制正则化项的影响程度,为模型的真实权重,为L2范数,表示权重参数的平方和。该损失函数由两部分组成:残差平方和和正则化项。残差平方和用于衡量模型的拟合程度,正则化项用于控制模型的复杂度。通过调节这个超参数,可以平衡模型的拟合效果和正则化效果,从而得到更好的模型性能。偏最小二乘回归通过在自变量和因变量之间找到最大协方差的主成分来建立模型,从而减少了特征之间的相关性,降低了模型的复杂度。这种方法可以帮助提高模型的泛化能力,降低过拟合的风险。
对定量精度的评估,可采用K折交叉验证,其原理如
, | (2) |
式中,表示第k次交叉验证中的训练集,由除第k个子集外的所有其他子集组成;表示第k次交叉验证中的测试集,即第k个子集。
最后,使用预测含量和真实含量间的偏差的均方根误差(RMSE)和回归模型的决定系数(R²)来定量评估,这里将5次验证结果取平均作为最终性能指标,如
, | (3) |
, | (4) |
式中,为第i折交叉验证中的样本数量,为第i折中第j个样本的真实标签,为第i折中第j个样本的预测值,为第i折中所有样本标签的均值。
实验中,选用斜长石、辉石、钛铁矿、橄榄石、玄武岩玻璃等五种月面的典型矿物混合作为模拟月壤,分别加入吸附水、结晶水、羟基化合水来对不同形态的水的拉曼光谱特征进行测试;然后利用不同含量的水体来对拉曼-含水量的检测能力进行测试。近年来月球遥感与原位探测结果表明,月面主要的矿物可分为硅酸盐矿物、氧化物矿物、自然金属、硫化物矿物和磷酸盐矿
为了模拟物的有效粒径与实际月壤粒径相近,块状样品经过粉碎和研磨,做成平均粒径50~100 μm的粉
测试目的 | 物质 | 样品 状态 | 含水 形态 | 含水质量比 / wt% |
---|---|---|---|---|
1235 | 模拟月壤+吸附水 | 压片 | 吸附水,冰 | 注射器实时注水 |
1235 | 模拟月壤+七水硫酸镁 | 压片 | 结晶水 | 1/2/3/6/9/12/15 |
123 | 模拟月壤+高岭土 | 压片 | 羟基水 | 1/2/3/6/9/12/15 |
4 | 辉石 | 压片 | 吸附水,冰 | 注射器实时注水 |
4 | 橄榄石 | 压片 | 吸附水,冰 | 注射器实时注水 |
1 | 七水硫酸镁 | 粉末 | 结晶水 | 51 |
1 | 十八水硫酸铝 | 粉末 | 结晶水 | 48 |
1 | 七水硫酸亚铁 | 粉末 | 结晶水 | 45 |
1 | 二水硫酸钙 | 粉末 | 结晶水 | 21 |
1 | 半水硫酸钙 | 粉末 | 结晶水 | 6 |
1 | 高岭土 | 粉末 | 羟基水 | 26 |
1 | 氢氧化钙 | 粉末 | 羟基水 | 46 |
1 | 氢氧化铝 | 粉末 | 羟基水 | 65 |
1 | 羟基磷灰石 | 粉末 | 羟基水 | 3.4 |
测试目的列表明了该样品所参与的实验流程及目的,1为不同形态、结构的含水信号识别及特征分析,2为确定含水量检测限,3为探究同种含水矿物中,谱峰强度与含水量的关系,4为探究不同含水矿物中,谱峰强度与含水量的关系,5为使用RR、PLSR、LR进行含水量反演
这里使用搭载有徕卡2 700 m金相显微镜的renishaw invia显微拉曼光谱系统,样品采用后向散射几何结构的532 nm倍频Nd:YAG连续激光器进行分析,激光器输出功率0~50 mW,使用数值孔径为0.4的50倍物镜将激光聚焦到样品,光斑理论直径为1.3 μm,实际光斑随样品粗糙度变化直径约3~10 μm,单次测量摄谱时间为30 s,光谱分辨率1 c
实验首先测得模拟月壤样品中各矿物组分的拉曼特征谱峰,其中模拟月壤压片样品(含水3 wt%)的矿物特征谱峰如

图1 模拟月壤中典型矿物的特征谱峰
Fig. 1 Raman characteristic bands of minerals in lunar soil simulants
实验测得含水模拟月壤压片样品和含水化合物粉末样品中,几种状态水的拉曼特征谱峰如

图2 高波数区域(2 600-3 900 c
Fig. 2 Raman spectral features of water in high wavenumber region (2 600-3 900 c
综上所述,拉曼振动光谱对氢键网络和水体的局域结构特征敏感,用拉曼光谱探测含水样品,不仅可以判断水的有无,还能准确识别出水的结晶水、羟基水、纯水冰赋存形态。
如
, | (5) |
式中,S为信号强度,N为背景噪声强度,即信号强度达到背景噪声强度的2倍。对于含七水硫酸镁结晶水样品的拉曼谱峰(

图3 模拟月壤中不同含量的吸附水、结晶水、羟基水、冰的拉曼光谱特征谱峰:(a)液态吸附水;(b)七水硫酸镁结晶水;(c)高岭土羟基化合水;(d)冰
Fig. 3 Raman spectral features of adsorbed water, crystalline water, hydroxyl water, and ice in Lunar soil simulant: (a) liquid adsorbed water;(b) magnesium sulfate heptahydrate crystalline water;(c) Kaolinite hydroxyl water; (d) ice

图4 模拟月壤中三种水体的特征谱峰的面积与含水量间的关系:(a)谱峰面积与水体含量的关系;(b)谱峰面积峰归一化后与水体含量的关系;(c)水体检测下限的估计
Fig. 4 Relationship between peak area and water content of three water forms in lunar soil simulants: (a) relationship between peak area and water content; (b) relationship between peak area peak normalized and water content; (c) estimation of water detection limit
综上所述,拉曼光谱法对模拟月壤基体中的吸附水、结晶水、羟基化合水能准确识别,且谱峰强度与含量呈线性关系,在满足采样点具有粒径、矿物基体成分等先验知识的情况下,有能力对月壤含水量进行半定量线性回归。对于吸附水、结晶水(MgSO4·7H2O)和羟基化合水(Al2Si2O5(OH)4)的检测限分别约为1.3 wt%、0.8 wt%和0.3 wt%,对于自由羟基含量更高的物质,检测限强于0.3 wt%。
在前述结论的基础上,进一步研究不同含水基体下,含水量与谱峰强度的关系。如图

图5 不同基体(橄榄石和辉石)内水的拉曼谱峰与含水量间的变化与比较:(a)橄榄石及液态水拉曼谱峰; (b)辉石及液态水拉曼谱峰; (c)两个基体内含水拉曼谱峰的比较; (d)谱峰强度与含水量间的关系
Fig. 5 Changes and comparison of Raman peaks with water content in different matrices (olivine and pyroxene): (a) Raman peaks of olivine and liquid water; (b) Raman peaks of pyroxene and liquid water; (c) comparison of water-containing Raman peaks in two matrices; (d) relationship between peak intensity and water content
另外,对于含吸附水样品中,水体拉曼谱峰的强度通常会随着含水质量比的增加而上升,水体的拉曼特征信号与荧光基线增强,达到过饱和后水体层荧光下降,如

图6 辉石中吸附水的拉曼谱峰与挥发时间的变化规律
Fig. 6 Changes in Raman spectra of pyroxene with adsorbed water over evaporation time
在实验中,根据2.2节中模拟月壤含水样品中吸附水质量比与高波数区拉曼光谱峰面积之间的线性关系,并考虑到结晶水(七水硫酸镁)与相同浓度的吸附水峰强相近。我们将模拟月壤中含吸附水与结晶水的高波数区域数据对含水量进行了定量分析。在微米级探测尺度下,我们制备的含吸附水样品局部含水量均匀,因此在2 600~3 900 c

图7 RR、PLSR和LR三种回归模型对模拟月壤中含水量的反演:(a) RR; (b) PLSR; (c) LR; (d) RR中最优正则化参数确定; (e) PLSR中最优主成分个数确定
Fig. 7 Inversion of water content in lunar soil simulanta by three regression models: (a) RR; (b) PLSR; (c) LR; (d) determination of optimal regularization parameter in RR; (e) determination of optimal number of principal components in PLSR
RR、PLSR和LR对模拟月壤中梯度含量水体的定量精度如
回归模型 | 训练集均方差 | 测试集均方差 | 测试集相关系数 |
---|---|---|---|
RR | 1.7811 | 1.3505 | 0.92279 |
PLSR | 2.3362 | 1.4059 | 0.91962 |
LR | 2E-18 | 3.0996 | 0.82130 |
为检验拉曼光谱对月面典型矿物中不同形态水的检测性能与探测能力,本文利用532 nm激光显微拉曼分别对模拟月壤、普通辉石、镁橄榄石中吸附水、冰、结晶水和羟基结构水的特征检测、信号识别;检验了模拟月壤中四种状态水的检测限,及其含量与信号强度间存在的规律;利用岭回归、偏最小二乘回归和线性回归对模拟月壤中不同浓度含水量进行定量分析。结果表明:(1)拉曼振动光谱对氢键网络敏感,对水体的局域结构特征敏感,除判断模拟月壤相关组分中水的有无,还识别水的结晶水、羟基水、纯水冰赋存形态。(2)显微拉曼光谱对模拟月壤中纯水、结晶水合物(MgSO4·7H2O)和羟基水合物(Al2Si2O5(OH)4)的检测限分别约为1.3 wt%、0.8 wt%和0.3 wt%;对含有大量氢键网络H2O的物质检测限约为1 wt%,对无氢键羟基占比高、羟基伸缩振动模式单一的物质中羟基的检测限强于0.3 wt%;对过饱和以及覆盖水状态下样品,拉曼光谱能分辨出水冰层-冰壤胶合体的存在形式,并能反演出水冰的厚度。(3)模拟月壤中从干燥到近饱和含水状态(本实验压片样品为0~15 wt%)的水特征峰强与含水量近似呈现线性关系;不同矿物的水峰峰强-含水量曲线斜率有差异,橄榄石和辉石的水体相对拉曼散射效率为1.9:1;使用岭回归和偏最小二乘回归可对模拟月壤中吸附水、结晶水的含水量进行定量反演,其回归决定系数
拉曼光谱可以对月面典型矿物及模拟月壤中含水的形态和含量进行探测,本研究结果能够丰富CE7拉曼光谱仪地面参考数据,为CE7拉曼光谱仪的物质成分及水冰探测提供支撑。
References
Peale S J, Schubert G, Lingenfelter R E. Distribution of sinuous rilles and water on the Moon [J]. Nature, 1968, 220(5173): 1222-1225. [百度学术]
Lammer H, Lichtenegger H I M, Kolb C, et al. Loss of water from Mars: Implications for the oxidation of the soil [J]. Icarus, 2003, 165(1): 9-25. [百度学术]
Reiss P, Warren T, Sefton‐Nash E, et al. Dynamics of subsurface migration of water on the Moon [J]. Journal of Geophysical Research: Planets, 2021, 126(5): e2020JE006742. [百度学术]
Lin Y, Tronche E J, Steenstra E S, et al. Evidence for an early wet Moon from experimental crystallization of the lunar magma ocean [J]. Nature Geoscience, 2017, 10(1): 14-18. [百度学术]
Starukhina L V, Shkuratov Y G. The lunar poles: Water ice or chemically trapped hydrogen? [J]. Icarus, 2000, 147(2): 585-587. [百度学术]
Mukherjee N R, Siscoe G L. Possible sources of water on the Moon [J]. Journal of Geophysical Research, 1973, 78(11): 1741-1752. [百度学术]
Saal A E, Hauri E H, Cascio M L, et al. Volatile content of lunar volcanic glasses and the presence of water in the Moon’s interior [J]. Nature, 2008, 454(7201): 192-195. [百度学术]
Hauri E H, Saal A E, Rutherford M J, et al. Water in the Moon's interior: Truth and consequences [J]. Earth and Planetary Science Letters, 2015, 409: 252-264. [百度学术]
Watson K, Murray B C, Brown H. The behavior of volatiles on the lunar surface [J]. Journal of Geophysical Research, 1961, 66(9): 3033-3045. [百度学术]
Alexander C M O D, Bowden R, Fogel M L, et al. The provenances of asteroids, and their contributions to the volatile inventories of the terrestrial planets [J]. Science, 2012, 337(6095): 721-723. [百度学术]
Andreasen R, Simmons S T, Righter M, et al. Lutetium-hafnium and samarium-neodymium systematics of Apollo 17 Sample 78236: Age and the importance of thermal neutron fluence on the lutetium-hafnium system [C]. 44th Lunar and Planetary Science Conference, 2013, 2887. [百度学术]
Barnes J J, Franchi I A, Anand M, et al. Accurate and precise measurements of the D/H ratio and hydroxyl content in lunar apatites using NanoSIMS [J]. Chemical Geology, 2013, 337: 48-55. [百度学术]
Barnes J J, Tartèse R, Anand M, et al. Investigating the H 2 O content and H isotopic composition of the primitive lunar magma ocean (LMO) cumulates [C]. European Planetary Science Congress, 2013, 105. [百度学术]
Bottke W F, Walker R J, Day J M D, et al. Stochastic late accretion to Earth, the Moon, and Mars [J]. Science, 2010, 330(6010): 1527-1530. [百度学术]
Biswas J, Sheridan S, Pitcher C, et al. Searching for potential ice-rich mining sites on the Moon with the Lunar Volatiles Scout [J]. Planetary and Space Science, 2020, 181: 104826. [百度学术]
Sowers G F. A cislunar transportation system fueled by lunar resources [J]. Space Policy, 2016, 37: 103-109. [百度学术]
Auer B M, Skinner J L. IR and Raman spectra of liquid water: Theory and interpretation [J]. The Journal of Chemical Physics, 2008, 128(22):224511. [百度学术]
Nozette S, Lichtenberg C L, Spudis P, et al. The Clementine bistatic radar experiment [J]. Science, 1996, 274(5292): 1495-1498. [百度学术]
Feldman W C, Maurice S, Binder A B, et al. Fluxes of fast and epithermal neutrons from Lunar Prospector: Evidence for water ice at the lunar poles [J]. Science, 1998, 281(5382): 1496-1500. [百度学术]
Pieters C M, Goswami J N, Clark R N, et al. Character and spatial distribution of OH/H2O on the surface of the Moon seen by M3 on Chandrayaan-1 [J]. Science, 2009, 326(5952): 568-572. [百度学术]
Li S, Lucey P G, Milliken R E, et al. Direct evidence of surface exposed water ice in the lunar polar regions [J]. Proceedings of the National Academy of Sciences, 2018, 115(36): 8907-8912. [百度学术]
Keller J W, Petro N E, Vondrak R R. The lunar reconnaissance orbiter mission–six years of science and exploration at the moon [J]. Icarus, 2016, 273: 2-24. [百度学术]
Gladstone G R, Hurley D M, Retherford K D, et al. LRO-LAMP observations of the LCROSS impact plume [J]. Science, 2010, 330(6003): 472-476. [百度学术]
Colaprete A, Schultz P, Heldmann J, et al. Detection of water in the LCROSS ejecta plume [J]. Science, 2010, 330(6003): 463-468. [百度学术]
Creech S, Guidi J, Elburn D. Artemis: an overview of NASA's activities to return humans to the Moon [C]. 2022 IEEE Aerospace Conference (aero), 2022: 1-7. [百度学术]
Djachkova M V, Mitrofanov I G, Sanin A B, et al. Selecting a landing site for the luna 27 spacecraft [J]. Solar System Research, 2022, 56(3): 145-154. [百度学术]
Kanu N J, Gupta E, Verma G C. An insight into India's moon mission–Chandrayan-3: The first nation to land on the southernmost polar region of the Moon [J]. Planetary and Space Science, 2024: 105864. [百度学术]
Wang C, Jia Y, Xue C, et al. Scientific objectives and payload configuration of the Chang'E-7 mission [J]. National Science Review, 2024, 11(2): nwad329. [百度学术]
Wang A, Jolliff B L, Haskin L A. Raman spectroscopy as a method for mineral identification on lunar robotic exploration missions [J]. Journal of Geophysical Research: Planets, 1995, 100(E10): 21189-21199. [百度学术]
Sharma S K, Angel S M, Ghosh M, et al. Remote pulsed laser Raman spectroscopy system for mineral analysis on planetary surfaces to 66 meters [J]. Applied Spectroscopy, 2002, 56(6): 699-705. [百度学术]
Qi X, Ling Z, Liu P, et al. Quantitative mineralogy of planetary silicate ternary mixtures using Raman spectroscopy [J]. Earth and Space Science, 2023, 10(5): e2023EA002825. [百度学术]
Huang Z. Lunar mineral distribution [J]. Encyclopedia of Lunar Science, 2017: 1-4. [百度学术]
Tompkins S, Pieters C M. Mineralogy of the lunar crust: Results from Clementine [J]. Meteoritics & Planetary Science, 1999, 34(1): 25-41. [百度学术]
Slyuta E N. Physical and mechanical properties of the lunar soil (a review) [J]. Solar System Research, 2014, 48: 330-353. [百度学术]
Zhang P, Dai W, Niu R, et al. Overview of the lunar in situ resource utilization techniques for future lunar Missions [J]. Space: Science & Technology, 2023, 3: 0037. [百度学术]
Sun Q. The Raman OH stretching bands of liquid water [J]. Vibrational Spectroscopy, 2009, 51(2): 213-217. [百度学术]
Kawamoto T, Ochiai S, Kagi H. Changes in the structure of water deduced from the pressure dependence of the Raman OH frequency [J]. The Journal of Chemical Physics, 2004, 120(13): 5867-5870. [百度学术]
Okada T, Komatsu K, Kawamoto T, et al. Pressure response of Raman spectra of water and its implication to the change in hydrogen bond interaction [J]. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2005, 61(10): 2423-2427. [百度学术]
Ponterio R C, Pochylski M, Aliotta F, et al. Raman scattering measurements on a floating water bridge [J]. Journal of Physics D: Applied Physics, 2010, 43(17): 175405. [百度学术]
Suzuki H, Matsuzaki Y, Muraoka A, et al. Raman spectroscopy of optically levitated supercooled water droplet [J]. The Journal of Chemical Physics, 2012, 136(23):234508. [百度学术]
Sun Q. The Raman OH stretching bands of liquid water [J]. Vibrational Spectroscopy, 2009, 51(2): 213-217. [百度学术]
Sun Q. The single donator-single acceptor hydrogen bonding structure in water probed by Raman spectroscopy [J]. The Journal of Chemical Physics, 2010, 132(5):054507. [百度学术]
Sivakumar T C, Schuh D, Sceats M G, et al. The 2500–4000 c
Li R H, Jiang Z P. Effects of anions on OH stretching Raman spectra of water [J]. Acta Physico-Chimica Sinica, 2007, 23(1): 103-106. [百度学术]
Larouche P, Max J J, Chapados C. Isotope effects in liquid water by infrared spectroscopy. II. Factor analysis of the temperature effect on H2O and D2O [J]. The Journal of Chemical Physics, 2008, 129(6):064503. [百度学术]
Lin K, Zhou X, Liu S, et al. Identification of free OH and its implication on structural changes of liquid water [J]. Chinese Journal of Chemical Physics, 2013, 26(2): 121-127. [百度学术]
Anand M, Tartèse R, Barnes J J. Understanding the origin and evolution of water in the Moon through lunar sample studies [J]. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2014, 372(2024): 20130254. [百度学术]