Summary
我们目前的技术来实现低速到脆弱的灰尘聚集在实验室之间的中间速度的碰撞。为了这个目的,两个真空落塔的设置已经开发了允许<0.01和〜10米/秒之间的碰撞速度。在碰撞事件是由高速摄像记录。
Abstract
对于调查的灰尘聚集在早期太阳系演化的目的,我们开发了两个真空下降塔中脆弱的灰尘聚集与规模可达约10厘米,孔隙率高达70%,可以发生碰撞。一个下拉塔主要用于非常低的冲击速度降低至低于0.01米/秒,并使用了双释放机构。冲突由两个高速摄像机,沿玻璃真空集热管中的两个灰尘聚集体的质心的框架落在记录在立体视图。另一自由下落塔使得电磁加速器,其能够平缓加速灰尘聚集到长达5米/秒的用途。在与另一个灰尘聚集体的释放,以自由落体组合,碰撞速度高达约10米/秒就可以实现。在这里,两个固定的高速摄影机记录下碰撞事件。在两个降塔,灰尘聚集在自由落体碰撞时,使它们失重和匹配条件在早期太阳系。
Introduction
人们普遍认为行星的形成始于微观小灰尘颗粒的非引力聚集成较大的灰尘聚集(见审查百隆&亚龙)1。灰尘颗粒的原行星盘内,由于布朗运动,相对漂移运动,和星云的气体湍流碰撞(见审查约翰森等人 )2。如果碰撞速度是足够低的,灰尘颗粒粘在一起,形成更大的团聚体。在过去数年丰富的实验室测量,导致了预测,对灰尘聚集任意群众和碰撞速度3的结果产生灰尘聚集碰撞模型。基本的碰撞结果被粘(一般为小骨料质量和低碰撞速度),弹跳,以及碎片(用于高冲击速度)。然而,这些阶段之间的过渡是不尖锐,还有其他结果,像, 如传质或糜烂。这种模式适用于典型的原行星盘预计在数千年的4厘米大小的灰尘聚集体的增长。厘米大小的灰尘聚集物的存在已被广泛地研究由天文观测,在过去几年,现在可以被认为是关系成立(参见审查TESTI 等人 )5,使我们得出结论,通过该机制原理的第一个宏观体在年轻的行星系统的形式已经确定。
然而,进一步增长到至少千米大小的身体不是那么清楚。对于陆地行星区,两个假设目前讨论(也见约翰森等人在这个问题上,最近的评论2和TESTI 等[5]。):(一)浓度厘米大小的灰尘聚集所, 如流不稳定6及后续gravit憩崩溃7,8和几个“幸运儿”,以更大的尺寸与后续质量吸积的传质过程9,10,11(II)的增长。在这两种模式,厘米大小的灰尘聚集发生相互碰撞的巨大数目在低到中等的速度。目前还不清楚这些碰撞的(除了弹跳)可能的结果是。
为了提高通过Güttler 等 3灰尘聚集碰撞模型,并更详细地碰撞调查宏观灰尘聚集在相关的速度政权中,我们设置了两个下拉塔在我们的实验室,其中个人总-总能碰撞被研究的很详细真空和微重力条件。既降塔具有1.5米自由落体的高度,这限制了观察时间来〜0.5秒。因此,我们观察到了碰撞的高速摄像机,百万像素格式和高达每秒7,500帧。为了获得最大的对比度和高的记录速度,亮场照明选择。照明是由高强度的LED面板从而提供和由漫射屏幕均匀。因此,高速摄像机查看碰撞的灰尘聚集的黑色物体在照亮屏幕的正前方。为了避免闪烁时,LED是直流供电。
为了实现低碰撞速度,这两个聚集的灰尘都放在上面彼此的双重释放机制。 吨的时间释放上部骨料下部1的结果中的ν= GT的相对速度之前,用克 =9.81米/秒2为地球的重力加速度。两个高速摄像机,它看到从两个方向90°的碰撞,通常在发布两粉尘总量之间(通常为T / 2上粒子后)。摄像机中的连续记录模式,这是由照相机的影响终止运行人成沙桶。在这种操作模式下,最大帧频为每秒1000张图像万像素的分辨率。采用这种设置,速度降低到低于0.01米/秒已经实现。由于双释放机构的机械设置的限制,最大相对碰撞速度为〜3米/秒。涉及灰尘聚集多达5厘米大小的碰撞进行了调查这个落塔。对于更高的碰撞速度高达约10米/秒,第二分出塔时,里面装有一个电磁加速器,其能够平稳地加速灰尘聚集高达5米/秒在垂直向上的方向。其它灰尘聚集是由一个双翼陷阱门释放机构保持,并且可以旋转自由在任何给定时间释放到自由下落。在这里,它没有任何意义用自由落体摄像机。我们宁愿用两个固定的高速摄像机,每秒百万像素分辨率高达7,500帧。由于较大diamet这个落塔的呃灰尘聚集到(并且可能上文)10cm大小都可以使用。
Protocol
注意:根据所使用的颗粒,它可以在相应的安全数据表中找到的危害性,口腔保护和安全装置必须由具有防尘工作的人佩戴。此外,还建议使用抽气系统,以保持环境空气无尘。
1,制备厘米大小的灰尘聚集样片
- 由米计算所需材料的量=Φρ0 V,其中m是所要求的质量,Φ是所需的体积填充率(体积填充率= 1 -孔隙度),ρ0是材料的密度,V是体积样品。 77克不规则硅酸盐粉尘(ρ0 = 2.6克/厘米3)是必需的,以实现70%的样品的孔隙率(体积填充率= 0.3)为5cm直径和高度的圆柱状样品,分别为。
注:解放军地面形成网开始微米大小的尘埃颗粒的混凝 - 主要由硅酸盐 - 成厘米大小的多孔体。甲充分研究和合适的实验室模拟材料为SiO 2,其可为不规则形状的粉末具有粒度分布范围从0.5到10μm,以及在单分散的球形颗粒获得更好的可比性,以理论模型的形式(见表1和图1)。
SiO 2的单体粒型 | 生产厂家 | 粒径 | 颗粒形状 | 例如图 |
单分散 | Micromod | 1.52±0.06&#181;米 | 球形 | 图1(左) |
多分散 | Sigma-Aldrich公司 | 0.1 - 10微米 | 不规则 | 图1(右) |
表1特性中的灰尘聚集碰撞实验中使用的SiO 2的粒子。
单分散(左)和多分散(右)二氧化硅用于生产宏观灰尘聚集体颗粒2图1。电子-显微镜图像。 请点击此处查看该图的放大版本。
- 采取与微米尺寸的SiO 2颗粒的容器( 见表1),其内容为0.5mm的筛目尺寸倒入筛子上。筛材料的足够量和填补计算质量到模具中。通过用手推压在活塞直到样品高度达到压缩该材料在模具中( 例如,5厘米)。转身对活塞模具,打开底板,轻轻样品推出。
注意:样品可以以多种形状(球形,圆柱形),大小(1毫米至10厘米)和孔隙率(60〜85%)来制备(参见图3)。样品可以被单独使用在碰撞实验或组合成群集,然后撞击其它聚集体或簇。
图2。的照片。灰尘聚集样品的尺寸和形状的变化的下列样品都显示:灰尘气瓶1厘米,2厘米,和直径为5厘米(后排),粉尘球体用1厘米,直径2厘米(中间行),并2-3毫米大小的Al 2 O 3球(正面)。 请点击此处查看该图的放大版本。
- 表征样品相对于孔隙度和均匀性。如果样品落在外面允许的利润,产生新的样本。
- 为了确定灰尘样品的孔隙率,由一个精确的平衡装置,测量其尺寸和它的质量决定了其体积。
- 使用X-射线断层扫描技术(XRT)12获得的均匀性和所生成的样品的孔隙大小分布的信息。
注:对于5厘米大小的灰尘聚集,我们发现偏离平均体积填充事实R, 即该样品的质量密度的比值和单体的灰尘颗粒,只有约1%的体积的样品体积的范围内,并通过到一个稍大的增加的体积填充率的材料密度朝向外边界12 8%。 图3示出了通过直径为5厘米和5厘米的高度的圆柱形灰尘聚集的切口的XRT重建。我们不使用XRT每个灰尘聚集,但检查随机样本的内部结构和均匀性。
为5厘米的高度和XRT分析后的5厘米直径的圆筒状灰尘聚集样品的内部结构的图3。重构。的灰度级表示的体积填充率,这是吨的质量密度的比值他样品和单体的尘埃微粒的物质的密度。从XRT重建,它是清晰可见的,这种高孔隙度样品用毫米大小的尘埃聚集组装。 请点击此处查看该图的放大版本。
落塔安装2。原理
- 释放机制:
在落塔的上部两个释放机制被附接一个在另一个的上面。他们每个人都拥有一个样本,并释放它变成自由落体。上部的释放和较低的颗粒之间的时间差确定碰撞的相对速度。根据形态和粒子的形状,合适的释放机构被选择。如果一个粒子加速机制被使用时,仅需要一个释放机构。- 颗粒上的一个字符串释放机制(球形样品,UPP呃颗粒):
本次发布机构由线性螺线管磁体和固体金属计数器一块。- 连接被释放到一个字符串的颗粒。通过螺线管磁体和实心金属计数器片之间的夹紧它保持在字符串中的位置。
- 以释放粒子,施加电流到螺线管磁体(参见电影1)。
- 活板门释放机制(球形样品,较低的颗粒):
本次发布机构由旋转螺线管磁铁,把一个粒子支架连接的。- 将微粒成半球形模,这是向下的旋转螺线管,当一个电流施加(见电影1)旋转。
- 这一机制也可用于颗粒簇或聚集团块的释放。在后一种情况下,将两片活板门释放装置上面彼此(见电影2)。
- 剪刀型双释放机制(圆柱试样):
本次发布机制包括两对旋转螺线管磁体的金属棒连接。每个释放机构的两个电磁体被放置,使得两个金属棒平行。- 两个样本放置到所述两个平行的杆。
- 施加电流到两个旋转螺线管的颗粒释放到自由下落。 (见电影3)。
- 双翼陷门释放机制(圆柱试样,在与粒子加速机制组合):
这释放机构包括两个弹簧加载的金属板,它们共同形成一个V形的粒子支架。两个金属板,以代替由一个金属棒,它是连接到一个旋转螺线管磁体保持。- 将圆柱形粉尘采样到封闭的活板门。
- 通过应用解锁陷阱门电流给螺线管磁体。为了避免弹跳回门的,涡流制动器阻止他们(见电影4)。
注意:在颗粒释放到自由下落而不初始速度和旋转是非常重要的。出于这个目的,一些发行机制已经发展(2.1.1 - 2.1.4)。
- 颗粒上的一个字符串释放机制(球形样品,UPP呃颗粒):
- 粒子加速机制:
或者通过一个预加载的弹簧或由电磁驱动的直线加速阶段的颗粒。既促进剂可配用于不同形状的粒子的样品架。 - 电子控制设备:
设置定时器,并释放电子到适当的值,以实现所需的碰撞速度和操作相机在质心的一帧。
注:颗粒释放,粒子加速和照相机释放的时间是由一组电子计时器,其功能在M被解释执行OVIE 5。
3,进行实验
- 低速碰撞(小降塔):
- 负载样品放入剪刀型双释放机制,并接近真空的玻璃管。
- 开始疏散,并设置定时器参数。
- 连接相机至其磁脱扣单元。开始连续摄像头录制。
注意:由于LED的亮场照明的高强度,高速相机的足够短的曝光时间,可以选择这样的曝光过程中颗粒的运动可以忽略不计。最重要的是,在光圈的相机目标必须设置为足够高的值,延长深度焦在落塔的整个直径。 - 当达到所需的真空质量,开关照明,按下启动按钮,并下载图像序列。
- 高速碰撞(大落塔):
- 负载样品到双赢克活板门释放机制和加速器,并关闭真空玻璃管。
- 开始疏散,并设置定时器参数。
- 开始连续摄像头录制。当达到所需的真空品质,开启灯光,并按下启动按钮。下载的图像序列。
注意:由于LED的亮场照明的高强度,高速相机的足够短的曝光时间,可以选择这样的曝光过程中颗粒的运动可以忽略不计。最重要的是,在光圈的相机目标必须设置为足够高的值,延长深度焦在落塔的整个直径。
4例实验
- 谨慎地装载样品到适当的释放机制。
- 低速碰撞(双重释放机制;0.09米/秒):5公分比5厘米,弹跳。
加载样本分为两个剪刀式释放机制。至实现0.09米/秒的碰撞速度,将粒子7毫米并设置释放机制的时间延迟到9毫秒。
注:在此冲击速度,灰尘样品弹开对方碰撞后。图像序列是由一个自由落体高速摄像机(见电影6)抓获。 - 高速碰撞(电磁加速器;7.4米/秒):2厘米比2厘米,碎片。
加载一个样本到双翼陷门释放机制;将其他样品的线性阶段加速器的样品夹持器。
注意:为了实现7.4米/秒的碰撞速度,较低的灰尘聚集体顺利地加速上行用2克,同时上部灰尘聚集降低。在7.4米/秒的相对速度,灰尘样品片段(见电影7)。 - 对大骨料小聚合高速碰撞:0.5厘米比5厘米,传质。
负载日Ë大样本到一个剪刀式释放机制;将小样品在弹簧加速器的样品夹持器。
注:为了达到所需的质量传递的碰撞速度,较低的灰尘聚集顺利加速上行,同时上聚集的灰尘落入。在这个相对速度,更小的样品碎片和传输少量的质量到更大的样本。由于相机落在沿上(更大规模的)粒子,采取高速摄像机的图像给大颗粒的印象,或多或少在休息(见电影8),从外面看这是不正确的落塔。
- 低速碰撞(双重释放机制;0.09米/秒):5公分比5厘米,弹跳。
- 关闭真空玻璃管。
- 小心打开真空阀,泵启动缓慢疏散,并设置定时器参数,为所需的碰撞速度所需的时间差。
- 连接相机至其释放器(如果自由落体摄像机使用)。开始连续摄像头录像并开启照明。
- 当达到所需的真空质量,按下释放按钮启动定时器序列。
- 下载记录的高速摄像机到计算机的图像序列。
5。数据分析
- 选择的背景和对象的灰度值之间的一个合适的阈值的灰度值。通过设置具有高于阈值灰度值的像素为白色(二进制值1),像素较低的灰度值,以黑色创建基于此阈值的二值图像(二进制值0)。
- 求出质量的粒子的在每个图像的中心的位置。一个很好的近似,以确定质量中心为对称的粒子是投影区域的中心。这是从二值化后的图像来计算。
- 使用质量的物体的中心的相对位置,时间信息从相机中的图像来计算相对速度(见电影9)。位置曲线的斜率显示在电影9的右手边。
- 如果发生碰撞反弹的,确定前后联系的相对速度。计算恢复系数, 即速度的碰撞和前后的比率。暗算恢复系数的相对速度。这种分析的一个例子示于图4。
图4实施例的碰撞反弹的分析。恢复系数, 即回弹速度与冲击速度的比值,被绘制为碰撞速度的函数。圆圈显示的数据为2厘米直径球形聚集的灰尘之三13(参见图2),三角形表示碰撞直径为5厘米和5厘米的高度的圆柱形灰尘聚集体之间( 见图2),而在0.3和0.4两个不同的体积填充因子,分别为12。数据显示减少恢复原状随着冲击速度的系数的趋势。 请点击此处查看该图的放大版本。
- 如果一个或两个粒子片段,确定尽可能多地通过测量它们各自的投影面积,并假设适当形状的最大片段的大小。
- 如果出现只有一个粒子的碎片,它通常一定量的质量转移到幸存的颗粒。确定传送质量的量通过测量增生体积,假定适当的形状和孔隙率来量化质量transfe效率较低。
Representative Results
使用在协议中描述的充分表征的灰尘聚集的样品(参见图1-3),在实验室中滴塔之一观察到任何碰撞会产生类似的碰撞在原行星盘的结果,科学的有价值的信息。到目前为止,我们已系统地研究了2厘米碰撞结果大小的球形灰尘聚集体(与填充量为0.5因素)之间的0.008和2.02米/秒13和5厘米大小的柱状灰尘聚集在速度范围内(与体积之间的填充因子在速度范围在0.3和0.5)之间的0.004和2米/秒12。我们发现了灰尘聚集体之间弹跳,作为下面速度〜占主导地位的结果0.4米/秒,这两种类型的灰尘聚集体(参见电影6中的示例)。在图4中,显示这些弹跳碰撞的恢复系数。圆圈表示实验2厘米大小的球形样品13和三角形代表之间5厘米大小的灰尘气缸有两个不同的包装从碰撞结果密度12。虽然个别实验归还系数分散广泛,恢复系数的平均值随碰撞速度减小。
兼具防尘通常聚集在冲击片段速度以上〜1米/秒(见电影7的例子)。为之间〜0.4和〜1米/秒的速度,仅两个碰撞灰尘聚集体的1碎裂的发生原因。在这种情况下,无碎片尘埃聚集由传质13获得质量的百分之几。上述速度极限是不清晰,但表示其中约不同的制度之间的界限在于2,11。对于不同大小的灰尘聚集和中等速度之间的碰撞,影响一般不会导致吨邻的两个灰尘聚集体的较大的碎片。在对面,大机构通过较小的撞击(见电影8)质量的一部分转移提高它们的质量。
为的情况下,其中两个灰尘聚集弹开彼此,从平动能传输碰撞之前(介意的灰尘聚集在碰撞之前不旋转)转换成平动能,回转动能,和其他(耗散)的能量通道(灰尘聚集体如压实)可以被确定。我们发现,对于中央碰撞(在该旋转能量可以忽略不计)耗散的能量的相对量的增加速度强烈地增加,并且是更高的灰尘量较低的填充因子聚集体12。这种行为可以通过分子动力学模拟12为蓝本。
TTPS :/ / www.jove.com/files/ftp_upload/51541/string_trapdoor.MP4“目标=”_blank“>电影1。的粒子上的一个字符串(高速动画(回放慢动作)顶部)和陷阱门释放机制(下)。
电影2 。双活板门释放装置的高速电影(回放慢动作)。这两个样本是直径为2毫米的Al 2 O 3颗粒,这仍然在自由落体由于发布在极低的干扰受限的团块。
电影3 。剪刀型双释放机制的高速电影(回放慢动作)。
电影4,高速双翼陷门释放机制的电影(回放慢动作)。
电影5 。计时器电子开关的上部和下部释放机制以及相机释放到自由下落的动画。
电影6 。两个5厘米大小的灰尘聚集气缸之间的碰撞弹跳的高速电影(回放慢动作)。这两个聚集的灰尘通过剪刀式双释放机制释放,碰撞0.09米/秒的速度。
电影7 。两个2厘米大小的柱状尘埃聚集的高速电影(回放慢动作)collidi纳克在7.4米/秒的相对速度。这两个片段的聚集完全。
电影8 。5毫米大小的灰尘聚集影响有5厘米大小的柱状固体靶的高速电影(回放慢动作)。由于4.3米/秒的冲击速度超过了小灰尘聚集的碎片的速度,这会崩解和转让其质量的一部分的目标,这是在电影中清晰可见。
电影9 。测定的粒子轨迹由一个半自动粒子跟踪算法。这里,示出了22厘米大小的球形聚集体尘埃之间的碰撞。
Discussion
由于高的机械精度,这两种滴塔的故障率是极低的。这是非常重要的,因为样品制备可能需要最多几个小时,这取决于尺寸,形状和所希望的灰尘聚集体的孔隙率。它应该提到的是大的灰尘聚集具有非常高的孔隙度是极其脆弱的,因此,难以处理。它可能发生的提取过模具或转移到落塔在这些灰尘聚集打破。在这些情况下,一个新的样本已被制备。因此,重要的是,小落塔可实现可靠的(和可预测的)碰撞速度下降到0.01米/秒11,13。最低的冲击速度迄今取得为0.004米/秒。这些小的冲击速度只能为自由粒子在微重力环境下即可到达。该实验室落塔是一种廉价和灵活的实现这样一个微重力设施。
Alternati已经的方法来实现低的冲击速度下利用悬浮技术14,15(通过电磁或气动悬浮EG),但通常引起碰撞粒子,其具有要考虑到在碰撞的分析之间的力。此外,磁悬浮往往引起旋转运动14,如果没人要,不允许旋转无碰撞,但是,在另一方面,甚至可能让冲突的旋转颗粒间逼真的模拟。在壳体的气动悬浮的,在碰撞过程中的气垫效应可能导致不希望的条件不匹配的那些原行星盘。然而,悬浮允许无限的观察时间和重复的实验,以便它必须考虑的替代落塔,如果时间限制是必要的。我们所有的努力,到目前为止,已经集中在SiO 2的硅酸盐在陆地行星的形成为r的代表EGION年轻的太阳能系统。由于大多数原行星盘的质量集中超出水冰的冷凝点,有必要也研究聚集体组成的微米尺寸的H 2 O的冰颗粒的碰撞行为。我们目前正在建立一个低温真空落塔用于此目的。但必须指出的是,温度在这样的模拟实验必须低于〜150 K,这就是所谓的“雪线”在原行星盘(温度的“雪线”划分,其中水是在气相内的区域相位并从那里被发现作为固体水冰)的外部区域。我们已经表明,微米尺寸的水冰颗粒的形成是可行的,其骨料可以生产16,让我们乐观地未来1-2年之内对他们的碰撞行为第一批成果。
Materials
Name | Company | Catalog Number | Comments |
Monodisperse SiO2 particles | Micromod | 43-00-153 | Particle diameter 1.52 ± 0.06 µm; particle shape spherical |
Polydisperse SiO2 particles | Sigma-Aldrich | S5631 | Particle diameter 0.1 – 10 µm; particle shape irregular |
References
- Blum, J., Wurm, G. The Growth Mechanisms of Macroscopic Bodies in Protoplanetary Disks. Annual Review of Astronomy and Astrophysics. 46, 21-56 (2008).
- Johansen, A., Blum, J., Tanaka, H., Ormel, C. W., Bizzarro, M., Rickman, H. The Multifaceted Planetesimal Formation Process. Protostars & Planets VI. University of Arizona Space Science Series. , (2014).
- Blum, J., Zsom, A., Ormel, C. W., Dullemond, C. P. The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? I. Mapping the zoo of laboratory collision experiments. Astronomy and Astrophysics. 513, (2010).
- Zsom, A., Ormel, C. W., Güttler, C., Blum, J., Dullemond, C. P. The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? II. Introducing the bouncing barrier. Astronomy and Astrophysics. 513, (2010).
- Testi, L., et al.
Dust Evolution in Protoplanetary Disks. Protostars & Planets VI. University of Arizona Space Science Series. , (2014). - Youdin, A. N., Goodman, J.
Streaming Instabilities in Protoplanetary Disks. The Astrophysical Journal. 620, 459-469 (2005). - Johansen, A., Youdin, A. N., Klahr, H. Particle Clumping and Planetesimal Formation Depend Strongly on Metallicity. The Astrophysical Journal Letters. 704, (2009).
- Bai, X. N., Stone, J. M. Dynamics of Solids in the Midplane of Protoplanetary Disks: Implications for Planetesimal Formation. The Astrophysical Journal. 722, 1437-1459 (2010).
- Windmark, F., Birnstiel, T., Güttler, C., Blum, J., Dullemond, C. P., Henning, T. h Planetesimal formation by sweep-up: How the bouncing barrier can be beneficial to growth. Astronomy and Astrophysics. 540, (2012).
- Windmark, F., Birnstiel, T., Ormel, C. W., Dullemond, C. P. Breaking through: The effects of a velocity distribution on barriers to dust growth. Astronomy and Astrophysics. , 544 (2012).
- Garaud, P., Meru, F. From Dust to Planetesimals: An Improved Model for Collisional Growth in Protoplanetary Disks. The Astrophysical Journal. 764, 146 (2013).
- Schräpler, R., Blum, J., Seizinger, A., Kley, W. The physics of protoplanetesimal dust agglomerates. VII. The low-velocity collision behavior of large dust agglomerates. The Astrophysical Journal. 758, 35 (2012).
- Beitz, E., Güttler, C., Blum, J., Meisner, T., Teiser, J., Wurm, G. Low-velocity collisions of centimeter-sized dust aggregates. The Astrophysical Journal. 736, 34 (2011).
- Beitz, E., Blum, J., Mathieu, R., Pack, A., Hezel, D. C. Experimental investigation of the nebular formation of chondrule rims and the formation of chondrite parent bodies. Geochimica et Cosmochimica Acta. 116, 41-51 (2013).
- Jankowski, T., et al. Crossing barriers in planetesimal formation: The growth of mm-dust aggregates with large constituent grains. Astronomy and Astrophysics. 542, (2012).
- Gundlach, B., Kilias, S., Beitz, E., Blum, J. Micrometer-sized ice particles for planetary-science experiments - I. Preparation, critical rolling friction force, and specific surface energy. Icarus. 214, 717-723 (2011).