X射线粉末衍射在保护科学:迈向文物艺术品腐蚀产物的常规晶体结构测定

Modern high resolution X-ray powder diffraction (XRPD) in the laboratory is used as an efficient tool to determine crystal structures of long-known corrosion products on historic objects.

该晶体结构测定和使用实验室的高分辨率X射线粉末衍射（XRPD）历史悠久的艺术对象的腐蚀产物细化过程是通过两个案例呈现的细节。

被调查的第一材料是铜甲酸钠氢氧化物水合物，铜₄娜_{4 O（HCOO）8（OH）2∙4H 2} O（样品1）构成的钠钙玻璃/铜合金复合历史性的对象（ 例如 ，搪瓷）在博物馆藏品，暴露于甲醛和从木制储藏柜，粘合剂等这种退化现象最近被定性为"玻璃诱导金属腐蚀"发射甲酸。

对于第二个案例，thecotrichite，CA _{3（CH 3 COO）3} CL（NO _{3）2∙6H 2} O（示例2），被选中，这是一种风化盐形成在其上存储在木制橱柜和陈列柜瓷砖和石灰石对象针状微晶。在这种情况下，木充当源为乙酸与来自工件或其环境可溶性氯化物和硝酸盐的盐反应。

几何结构的知识有助于保护科学，以便更好地了解生产和衰变反应并允许以混合物的常见的情况充分定量分析。

保护科学应用于文物的保护科学（通常化学）方法。这包括制作器物的调查（"技术艺术史"？这怎么在那个时候做的），它们的衰变途径为前提，以发展适当的保护治疗。通常情况下，这些研究处理像碳酸盐，甲和乙酸盐金属有机盐。其中一些已用合适的化合物（ 例如，醋）被故意制造的，他人与大气（二氧化碳或从室内空气污染的羰基化合物）劣化反应得到^1。作为事实上，许多这些腐蚀材料的晶体结构仍不清楚。这是不幸的事实，因为几何结构的知识有助于保护科学，以便更好地了解生产和衰变反应并允许在混合物的情况下充分定量分析。

根据所关注的材料形成足够的尺寸和质量的单晶的情况下，单晶衍射是选择用于晶体结构的测定方法。如果这些边界条件不满足，粉末衍射是最接近的替代产品。粉末衍射相比单晶衍射作为最大的缺点在于的倒数ð -矢量D *（散射矢量）的取向的信息丢失。换句话说，一个单一的衍射斑点强度被涂污的球体的表面上。这可被认为是三维衍射（=倒数）空间上的投影的粉末图案的一维2θ轴。因此，散射方向不同，但相同或相似的长度载体，系统地重叠或意外使得它很难甚至不可能对这些反射^2个独立^的 （Figure 1）。这也是最主要的原因，粉末衍射，尽管仅仅四年后的第一单晶实验^3,4，主要用于相位识别和量化超过半世纪早期的发明。尽管如此，一个粉末图案的信息内容是巨大的，如从图2中可以很容易地推导出来。真正的挑战，但是，是揭示尽可能多的信息有可能在一个常规的方法。

实现这一目标的一个关键步骤，毫无疑问，来自特维德雨果于1969年⁵想法是谁发明从粉末衍射数据的晶体结构的细化局部优化技术。该方法不修改单强度而是针对日益复杂，因此采取的高峰重叠本质兼顾的典范整个粉末格局。从那个时候起，用粉末衍射技术，科学家们不再局限于数据分析bŸ方法单晶调查开发。 Rietveld方法的发明几年后，为从头计算结构测定粉末衍射法的功率被识别。如今，自然科学和工程用粉末衍射的几乎所有的分支，以确定更多，更复杂的晶体结构，尽管该方法仍不能视为常规。在过去的十年中，在实验室的新一代粉末衍射仪的设计提供高分辨率，高能量和高强度。更高的分辨率立即带来更好的峰的分离，而更高的能量吸收战斗。基于基本的物理参数更好的峰形说明（ 图3），这样做的好处是布拉格反射从而提供更详细的结构调查更准确的强度。拥有现代化的设备和软件，甚至微观结构参数，如域的大小和microstra在例行的粉末衍射数据推断。

从粉末衍射数据为晶体结构测定所有算法使用单峰强度，整个粉末图案或两者的组合。传统的单晶倒易空间技术往往不能由于可用的观察和结构参数之间的不利比例。这种情况下，与引进的"充电翻转"技术^6（ 图4）和全局优化方法中的直接空间，其中模拟退火技术^7（ 图5）是最突出的代表发展的显着改变。特别是，引入化学知识到结构确定处理使用刚体或关于键长分子化合物的已知的连接和角度强烈减少必要的参数的数量。换言之，在为每一个原子三个位置参数代替，只有外部（和一些内部）度的原子团的自由的需要来确定。它是该还原结构复杂，这使得该粉末方法的真正替代单晶分析。

作者^8,9的两个开创性案例研究证明，它可以解决使用粉末衍射数据的复杂的腐蚀产物复杂的晶体结构。相比于其它方法的晶体学研究的优越性是由以下事实，在这两种情况下，所报告的公式有考虑到解决晶体结构后待校正等等证明。

这两种材料的受调查在博物馆的发生有关，其存储在木制橱柜或暴露于羰基其它形式的污染源。被调查的第一材料是铜甲酸钠氢氧化物九头蛇德，铜₄娜_{4 O（HCOO）8（OH）2∙4H 2} O（样品1），形成于博物馆藏品苏打玻璃/铜合金复合历史性的对象（ 例如，搪瓷），暴露于甲醛和甲酸从木制的储存柜，胶粘剂等 。这种退化现象最近被定性为"玻璃引起的金属腐蚀^"10。对于第二个案例，thecotrichite，CA _{3（CH 3 COO）3} CL（NO _{3）2∙6H 2} O（示例2），被选中。 Thecotrichite是一个经常出现风化形成盐瓷砖和石灰石博物馆的对象，它们存储在橡木柜和陈列柜针状结晶。在这种情况下，木充当源为乙酸与来自工件可溶性氯化物和硝酸盐的盐反应。

在文本中，结构D的各个步骤的以下部分从保护科学应用到腐蚀产物确定过程采用粉末衍射数据进行了详细的介绍。

1.样品制备

1. 材料收集
   1. 使用从不透明的蓝绿色凸圆形设置手术刀和镊子上的一个历史性扣仔细挑选样品1少量的（小于1毫克）的数字显微镜下，属于Rosgartenmuseum康斯坦茨（RMK-1964.79）（收集图6）。
   2. 小心划伤样品2几毫克与釉面瓷砖表面的手术刀，在近代早期，在德国南部生产的约会，与41 X 29×3.5厘米大小，和国家博物馆符腾堡州的集合的一部分（没有的E 3004）（ 图7）。
      注:瓷砖储存从20世纪80年代至2004年的木碗柜。背侧很大程度上受到存在作为相纯品（ 图7b）thecotrichite风化。
2. 样品保持器的制备
   1. 研为s在一个小玛瑙研杵仔细amples。
   2. 分发样品1，它由几个附聚颗粒（<1毫克），两个薄的X射线透明的聚酰亚胺箔片之间，并与直径为8毫米的中心开口它们安装在透射样品架。固定在衍射的θ-圆的发送样本保持器。
   3. 在0.5mm直径的硼硅玻璃毛细管填充样品2（约5毫克）。
      1. 这样做，将粉末少量成使用刮刀毛细管的漏斗。然后将粉末向下到有电振动器的帮助毛细管的尖端。通过将毛细管厚壁玻璃管小巧的粉和人工挖掘它放在桌子上。
      2. 继续，直到达到3厘米左右的填充高度。小心切开玻璃毛细管在使用薄刚玉叶片的大约为4厘米的高度，并用打火机密封的开口端。
      3. 放置少量蜂蜡的黄铜销的开口和熔融使用烙铁它。然后将在熔化的蜡毛细管，并保持直立，直到凝固。安装在测角仪头上的铜针和衍射仪的θ-圈固定测角仪头。
      4. 最后，通过叠加一个十字数码相机投影支持的扳手手动中心通过反复对准四度的测角头自由的（使用XY使用翻译阶段弯曲滑梯和两个译本2转）安装的毛细血管。

2，数据收集

1. 为_（α1 -辐射与主光束约翰型Ge（111）单色器为了铜K），在室温下使用一个高分辨率粉末衍射仪的样品1（ 图8）和图2（ 图9）的实验室XRPD图谱其中配备有线性位置敏感SI力诚条探测器大约12°2θ的开口。
   1. 测量样品1的范围与传输模式0.015°2θ一个步长20小时，从5-85°2θ（扫描模式:传输;扫描类型:2ThetaOmega;欧米茄方式:移动，PSD模式:移动）。为了达到更好的颗粒统计圈旋转上。
   2. 为期6小时的覆盖在2θ范围5-60°，在德拜 - 谢乐模式0.015°的步长记录的样品2（扫描模式:德拜谢乐;扫描类型:2Theta;欧米茄模式:固定; PSD模式:移动）。为了达到更好的颗粒统计圈旋转上。

3.晶体结构测定和细化

注意:对于样品1和2的晶体结构的测定和细化，一个复杂的计算机程序用于^11。它用图形用户界面或通过基于文本的输入文件运行秒。后者化妆用一个复杂的脚本语言。使用样品1的结构分析的不同阶段的样品输入文件列在表S1，S2，S4-S8。的一般方法是对样品2相同。

1. 峰值搜索
   1. 使用根据低阶（ 图10）的Savitzky-格雷polynoms的第一和第二导数，以制造商的协议通过在布拉格反射的半最大值设定1-1.5倍全宽的数量级上的回旋范围执行自动峰值搜索（0.12°2θ为样品1），调整噪声阈值，以1.5-2倍的估计的标准偏差（1.74样品1）和限制性的粉末型态其显示清楚区分的峰的部分的检索（5-66° 2θ为背景上面的示例1）。
      1. 启动通过双击该图标的程序。点击加载扫描文件 。拉下拉菜单中选择XY数据文件（* .xy）。双击1-GNM-4145-P4-Kapton.xy。
      2. 扩大范围1-GNM-4145-P4-Kapton.xy。点击发射简介 | 负荷排放简介 | CuK1sharp.lam。
      3. 点击按钮自动插入峰 。取消选中删除K-2阿尔法峰 。用滑杆设置峰宽为0.12。设置有滑动条噪声阈值 1.74。按下按钮添加峰 。按关闭 。
   2. 手动使用人的眼睛和正确的整合特性，删去非布拉格峰这显然是由于计数统计和添加峰它们是可辨认但隐藏在其它峰的尾部的自动检测到的峰。建立索引文件，套30至40为反射使用S每个样品TANDARD设置，但允许所有晶系（ 表S1）。
      1. 鼠标放大模式，使用鼠标滚轮滚动。按F3开峰详细信息窗口。通过按鼠标左键设置峰。按F9删除峰。关闭峰详细信息窗口。
      2. 点击高峰期 。马克所有峰通过点击位置的黄色。在黄左击标柱。单击全部复制/选择 | 创建索引范围 。取消范围1-GNM-4145-P4-Kapton.xy。选择范围索引 。选择所有布拉维点阵（点击使用突出显示列表中，并设置了刻度线）。
         
2. 索引
   1. 索引应用"重复使用奇异值分解"算法^12（ 表S2），直到发现样品1原始四方单元电池和用于样品2（参见晶胞参数见表1）一个简单单斜晶格。
      1. 按运行按钮（F6）。按是保持索引的解决方案。按下按钮解决方案 。通过突出显示左点击按钮1第一个解决方案。右键单击选中的（黄色）解决方案。单击全部复制/选择 。取消选择范围索引 。选择范围1-GNM-4145-P4-Kapton.xy。
         注:其他的单元电池示出优点的所有的一个高的数字相比，在表1中列出的具有多个卷但是未观察到的反射的高数量的程序自动建议基于所观察到的反射灭绝的最可能的空间群是P 4。 ₂ / N（86）样品1和P 2 ₁ / 一 对于样品2（ 图11）。
   2. 通过搜索使用晶体A卷^13（ 表S3）的国际表灭绝手动确认这些发现。公式个单位晶胞从平均体积增量的数量估算Z = 8样品1和Z = 4样品2。
3. 全粉谱拟合
   1. 装修根据波利¹⁴为粉模式进行全粉格局。为的峰轮廓的描述中，使用基本参数（FP）方法^15（ 图12）。由高阶正交切比雪夫多项式（一般为8），并描述了在低衍射角散射空气额外的1 / X项模型的背景。洛伦茨-极化因子设置为使用_Cu-Kα1拉迪27.3这是布拉格角为戈（111）单色器通货膨胀表S4包含了所有相关的输入参数。
      1. 按添加HKL阶段 。扩大范围1-GNM-4145-P4-Kapton.xy。展开hkl_Phase。点击索引详细 | 粘贴索引的详细信息 。单击背景 。改变以 8（精炼）。勾选1 / X BKG（精炼）。
      2. 点击仪器 。 主半径（mm）设置为217.5。 二级半径（mm）设置为217.5。勾选点探测器 。勾选接收狭缝宽度 ，0.1（固定）的价值。勾选全部轴向模型 。设置源长度（mm）到6（固定）。设置样品长度（mm）到6（固定）。集RS长度（mm）到6（固定）。
      3. 点击更正 。勾选零误差值为0（精炼）。蜱LP因素 ，价值27.3（修正）。点击杂项 。设置转换次数。步骤 2。勾选启动X，值8。蜱完成 X，值75.点击高峰期 | 删除高峰期 | 是的 。
      4. 点击hkl_Phase | 微观结构 。勾选惊魂尺寸长 ，价值为200（精炼）。勾选惊魂尺寸克 ，价值为200（精炼）。蜱应变L，值0.1（精炼）。蜱应变G，价值0.1（精炼）。按运行按钮（F6）。
   2. 创建布拉格列表峰适于充电翻转^6。点击hkl_Phase。点击电荷翻转输出 。打钩A为CF，价值CF.A 文件 。勾选为CF HKL文件 ，值CF.hkl。按运行按钮（F6）。菜单File下拉| 导出到INP ​​文件 。文件名 ​​，价值Pawley.INP。按保存 。
      注:对于连续结构测定，保持所有的仪器，峰形和晶格参数固定的。
4. 晶体结构测定
   注意:三种方法（以迭代的方式使用）的组合被用来确定试样1和2的晶体结构。
   1. 首先，使用充电翻转⁶通过包含切线式^16（ 图13）所支持的方法找到的最较重原子的位置。删除观察REF与门经一个d-间距小于1.28（样品1）和所计算的范围扩大到0.9的d间距。所有必要的参数列于表S5中 。挑从结果列表中识别的电子密度，正确的错误分配原子类型相似的散射力碳，氧，氮手动的原子。
      1. 拉下文件文件 | 关闭全部 -点击确认是 。拉下启动启动 | 启动内核 。
      2. 拉下启动启动 | 设置INP文件 。选择CF.INP（制备输入文本文件，在表S5中列出）。通过单击打开进行确认。按R键取消（F6）。按STOP按钮（F8）约了。 20000次。确认充电翻转按完成是。
      3. 按临时输出窗口显示在Z矩阵格式按钮来选择原子。按云的选项对话框按 ​​钮。 N设置为挑到45。蜱对称 。 媒体报道按钮。复制临时输出，并将其保存到一个名为FirstGuess.str一个文本文件中。关闭充电翻转图形窗口。
   2. 其次，应用模拟退火^7,17（ 图14）的全局最优化方法，找到所有缺失的非氢原子（主要是甲酸盐/乙酸酯基团的碳原子）的位置。使用该软件提供的自动退火方案。仅须的比例因子和位置和/或所选的原子模拟退火的职业参数。对于样品1，合并1.1埃的半径范围内的钠原子和氧原子以检测特殊的位置（ 表S6）。
      1. 拉下启动启动 | 设为INP文件 。选择SA.INP（制备输入文本文件，在表S6所列）。通过单击打开进行确认。按运行按钮（F6）。几千个周期后按停止按钮（F8）。按确认是以 更新输入文件 。
   3. 第三，关闭模拟退火和注释掉命令Auto_T（0.1）切换到特维德细化⁵模式。修复所有确定位置参数。包括差异傅立叶图（F _OBS -F _计算 ）的计算（ 图15， 表S7）来检查下落不明电子密度。包括在原子列表中还发现了原子和细化原子的位置和占有率。
      1. 拉下启动启动 | 设置INP文件 。选择Fourier_search_for_C.INP（准备输入的文本文件，在表S7中列出）。通过单击打开进行确认。按运行按钮（F6）。检查两个图形输出窗口显示几乎完全晶体结构和差异傅立叶地图。
         注:循环三种结构求解方法迭代地在故障的情况下。如果需要的话，应用抗凸点约束轻原子（碳原子和氧原子）（ 参见表S7）。
5. Rietveld精修
   1. 对于样品1和2的最终晶体结构精炼中，使用的Rietveld全图案精制方法^5（ 图16）。为了避免原子的无意义原子转移在样品2的乙酸和硝酸基，采用所谓的软约束（也称为约束）的基础上理想化的键长和角度。使用标准的软件计算丢失的氢原子的理想化位置是^18（ 表S8）。
   2. 缩小为每个晶体结构的非氢原子的各向同性原子位移参数。在样品二，模型的情况下，引起的微观的布拉格反射通过包括对称的宽度的视各向异性适于二阶球谐。
   3. 最后创建样品1（ 图17）和样品2的晶体结构（ 图18），并且从现在起，可用于充分定量相分析两个晶体学信息文件（CIFS）的突起的图。这种全定量分析的一个例子是使用样品1的晶体结构在图19中给出。

高分辨率XRPD来测定的具有历史意义的对象的两个长知腐蚀产物先前未知的晶体结构。将样品从两个博物馆对象采取的，并仔细研磨它们在传输和毛细管样品架进行密封前（ 图6,7）。用的传输和使用单色X射线德拜-谢乐几何本领域实验室高分辨率粉末衍射仪的状态测量标准进行（ 图8）。

从粉末衍射数据结构测定的标准化过程，使用最近开发高效的算法开发顺序如下:峰值位置（ 图10）测定，索引和空间群确定（ 图11），全粉图案接头（ 图12 ），结构测定（ 图13-15），和Rietveld精修（ 图16）。通过迭代地结合倒数（电荷翻转）（ 图13）和直接空间（模拟退火）（ 图14）与差分傅里叶分析方法（ 图15）执行的两种化合物的晶体结构的测定。

这些化合物的晶体结构的测定（ 图17，18）提高了我们的衰变机制的理解，并允许腐蚀产物充分定量相分析（ 图19）。

倒易空间图1.粉末衍射 。的往复空间区域的插图即在粉末衍射访问测量。较小的圆表示Ewald球。在粉末测量的倒易晶格旋转以采样所有方向。一个等价的操作是在周围的倒易空间的原点所有可能的方向旋转Ewald球。体积扫出（图中的区域）是在实验中可访问的倒易空间的区域^。2

图2.粉模式的信息内容 。背景，峰位，峰强度和峰形的四个主要贡献粉末衍射图案，信息内容的示意图图片^。2 ，请点击此处查看该图的放大版本。

古尔3"SRC ="/文件/ ftp_upload / 54109 / 54109fig3.jpg"/>
图3. 波利契合 。波利全粉模式契合使用基本参数方法在德拜-谢乐几何戈（220）单色莫_-Kα1辐射（λ= 0.7093埃）测量实验室_6个标准的粉末格局。以下四个圈已应用:纯洛伦兹发射剖面，在赤道平面的接收狭缝的一个帽形函数，一轴向卷积服用filament-，样品 - 和接收狭缝的长度和次级索勒狭缝在内，和一个小高斯的贡献相关位置敏感探测器¹⁹ 请点击此处查看该图的放大版本。

g4.jpg"/>
图4. 充电翻转方案 。在翻转用于从粉末衍射数据结构的决心倒易空间电荷翻转过程的方案和流程图（如小图）。

图5. 模拟退火方案 。在用于从粉末衍射数据结构直接决定一个空间模拟退火过程的流程图^。19

图6.产地样品1。携带历史悠久的艺术对象，铜₄娜_{4 O（HCOO）8（OH）2∙4H 2} O（样品1）。历史扣，属于Rosgartenmuseum康斯坦茨（RMK-1964.79）的集合。TTPS://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/54109/54109fig6large.jpg"目标="_空白">点击此处查看该图的放大版本。

图7.产地的样品2。携带Thecotrichite（样品2）历史悠久的艺术对象。 Thecotrichite从国家博物馆符腾堡州（a）及其布满白色晶体thecotrichite背面（ 二 ）收集琉璃瓦^。9 ，请点击此处查看该图的放大版本。

示出的Cu ₄的Na的散射的X射线强度的样品1。屏幕快照的图8.粉末衍射图案_{4 O（HCOO）8（OH）2∙4H 2} O在室温条件下为衍射角的函数（样品1）。 请点击此处查看该图的放大版本。

显示thecotrichite的散射X射线强度样品2截屏图 9. 粉末衍射图案在环境条件（示例2），作为衍射角的函数。 请点击此处查看该图的放大版本。

图10. 峰值搜索样品1 荣>。屏幕截图显示铜₄的X射线散射强度娜_{4 O（HCOO）8（OH）2∙4H 2} O和使用平滑滤波器Savitzky -格雷的第一和第二衍生物自动寻峰算法的结果。 请点击此处查看该图的放大版本。

图11.索引结果显示的索引和空间群测定结果样品1。屏幕截图铜₄娜_{4 O（HCOO）8（OH）2∙4H 2} O. 请点击此处查看该图的放大版本。

期间结构测定过程样品1。屏幕截图图13.充电翻转铜钠_{4 4 O（HCOO）8（OH）2∙4H 2}在空间群P O 4 ₂ / N使用费与直方图翻转方法匹配。与预先指定的原子类型的晶体结构的一部分，是人准备可见。 请点击此处查看该图的放大版本。

图14. 模拟的过程中为结构测定过程样品1。屏幕截图退火铜钠_{4 4 O（HCOO）8（OH）2∙4H 2}在空间群P O 4 ₂ / N采用模拟退火的全局优化方法。晶体结构的一部分是已经显现。 请点击此处查看该图的放大版本。

图15。 差异傅里叶分析为期间铜₄娜_{4 O（HCOO）8（OH）2∙4H 2}在空间群P O 4 ₂ / N采用了结构测定过程中丢失原子搜索的样品1。屏幕截图差值Fourier方法。晶体结构是和额外的电子密度绘制。 请点击此处查看该图的放大版本。

图16. 特维德适合显示₄铜的特维德情节娜_{4 O（HCOO）8（OH）2∙4H 2} O空间群P 4 ₂ / N 样品1。屏幕截图的 。所观察到的图案（蓝色），最好的Rietveld拟合轮廓（红色）和观测和计算的轮廓（下灰色）显示之间的差异曲线。 请点击此处查看该图的放大版本。

图17. ₄铜的晶体结构的样品1。投影娜_{4 O（HCOO）8（OH）2∙4H 2} O沿晶体c轴的晶体结构 。含有铜和中央分子是钠多面体。 请点击此处查看该图的放大版本。

翁>图18. 晶体结构样品2。thecotrichite的晶体结构的突起，沿着（a）该c轴和沿着将b轴（ 二 ）提出。多面体颜色:Ca1的:品红，钙:青色的Ca3:绿色⁸ 请点击此处查看该图的放大版本。

含样品1 图19. 定量分析 。含Cu ₄娜_{4 O（HCOO）8（OH）2∙4H 2} O作为主相和Cu _2（OH）3腐蚀样品全定量相分析的情节特维德（HCOO）和Cu ₂ O轻微阶段。/54109fig19large.jpg"目标="_空白">点击此处查看该图的放大版本。

铜₄ 表1.晶体和结构细节娜_{4 O（HCOO）8（OH）2∙4H 2} O和Ca _{3（CH 3 COO）3} CL（NO _{3）2∙6H 2} O（thecotrichite）。

补充表格

表S1。峰值搜索₄铜娜₄ O后输入文件（HC_{OO）8（OH）2∙4H 2} O（只有1个峰值显示在峰值列表）。 请点击此处下载此文件。

表S2。对Cu ₄索引娜_{4 O（HCOO）8（OH）2∙4H 2} O.输入文件请点击这里下载此文件。

表三，对从结晶量A.国际四方的表空间群反射条件列表请点击这里下载此文件。

表S4。INP拟合根据铜₄的帕利方法娜₄粉末整体格局UT文件_{O（HCOO）8（OH）2∙4H 2} O（只有少数布拉格反射显示在高峰列表）。 请点击此处下载此文件。

表S5。对Cu ₄充电翻转娜_{4 O（HCOO）8（OH）2∙4H 2} O.输入文件请点击这里下载此文件。

表S6。对Cu ₄的模拟退火娜_{4 O（HCOO）8（OH）}输入文件_{2∙4H 2} O. 请点击这里下载此文件。

表S7的输入文件的差异傅里叶分析铜钠_{4 4 O（HCOO）8（OH）2∙4H 2} O.。 请点击此处下载此文件。

表S8。对Cu ₄的最终Rietveld精修娜_{4 O（HCOO）8（OH）2∙4H 2} O.输入文件请点击这里下载此文件。

XRPD is a suitable technique for conservation research as it is non-destructive, fast and easy-to-use. XRPD data can be used in routine qualitative analysis, owing to the fact that the powder pattern is a fingerprint signature to the corresponding crystal structure. The biggest advantage of XRPD over other analytic techniques is the ability of performing simultaneous qualitative and quantitative analysis of crystalline constituents in mixtures by using the Rietveld refinement method⁵. Moreover, the presence of amorphous content can be detected and its amount estimated. However, this procedure requires knowledge on every crystal structure present in the mixture that is a subject to investigation.

To apply the method of XRPD routinely for structure determination to conservation science, several critical boundary conditions for the laboratory powder diffractometer must be fulfilled: 1.) To avoid preferred orientation in powder samples, transmission or even better Debye-Scherrer geometry must be used. 2.) Laboratory powder diffractometers should be equipped with a primary beam monochromator to ensure strict monochromatization and a position sensitive strip detector for high intensity (= good counting statistics) and high resolution. This particular type of instrument leads to sharp peak profiles which can be adequately described by few fundamental parameters being of great benefit for the separation of overlapping reflections.

Indexing of the powder pattern which is often regarded as the bottleneck in the structure determination process should be done with exhaustive methods like "singular value decomposition", which is also insensitive to small amounts of impurities. Due to the strongly reduced information content of a powder pattern as compared to a single crystal data set, a sophisticated combination of direct and reciprocal space structure determination algorithms is needed for a high success rate. The combination of charge flipping, simulated annealing and difference-Fourier analysis has been proven to be among the most promising approaches. Providing that the material under investigation is reasonably crystalline, crystal structures with 20-25 structural parameters can nowadays been solved almost routinely from powder diffraction data if the procedure described above is used. It can be expected that this limit can be pushed to much more complex crystal structures with the advent of better instrumentation, the use of synchrotron radiation, and even more sophisticated structure determination algorithms.

Even after 250 (!) years of conservation research and 100 years of crystal structure analysis, there are still many crystalline corrosion products on artifacts of unknown exact composition and structure. This is mainly due to the unavailability of naturally or synthetically grown single crystals of suitable size. XRPD data analysis as described here can overcome this restriction since powder samples are amenable to investigation. A quantum leap forward in Conservation Science as well as in other fields!

The authors have nothing to disclose.

作者非常感谢克里斯蒂娜·史蒂芬尼女士进行XRPD测量。玛丽安·舒赫和雷秸卡Kuiter（艺术与设计斯图加特州立学院）被确认为瓷砖（ 图7）的图片。