Molecular dynamics simulation study of sintering mechanism and thermal conductivity of nano-Ag particles
-
摘要:
针对纳米颗粒的烧结机理不明确及烧结过后热导率显著下降等问题进行探究. 利用分子动力学仿真的方法，通过热导率的计算验证了势函数选取的正确性，模拟了银纳米颗粒的烧结过程，观察了烧结过程中微观结构的变化，得到了烧结颈生长长度与烧结时间拟合关系，并计算了烧结颈处声子热导率变化情况，最后通过透射电子显微镜观察到试验烧结纳米银颗粒具备与仿真结果相同的五重孪晶结构. 结果表明，五重孪晶为纳米颗粒烧结后的稳定结构，纳米颗粒烧结过程中烧结颈生长长度与时间成幂函数关系，烧结颈处微观结构使纳米银的热导率降低了4.42%.
Abstract:Nanosilver paste can achieve low-temperature sintering and high-temperature service, due to its special size effect and high specific surface area, which has a broad application prospect in the third-generation semiconductor packaging. However, the sintering mechanism of nanoparticles has not been clear and the thermal conductivity has decreased significantly after sintering, which needs to be further investigated. In this paper, the sintering process of silver nanoparticles was simulated by molecular dynamics simulation. The correctness of the potential function was verified by the calculation of thermal conductivity, firstly. And then the microstructural changes of silver nanoparticles were observed, the relationship between the sintering neck growth length and the sintering time was obtained, and the change of the phonon thermal conductivity at the sintering neck was also calculated. Finally, the experimental silver nanoparticles with the same characteristics as the simulation results were observed by the transmission electron microscope. The results showed that the fivefold twins were the stable structure of nanoparticles after sintering, the sintering neck length growth during nanoparticle sintering was a power function of time, and the microstructure at the sintering neck reduced the thermal conductivity of nanosilver by 4.42%.
-
Keywords:
- MD simulation /
- nano-particle /
- thermal conductivity
-
0. 序言
近年来随着新能源汽车、5G通讯和云计算等领域的发展，对计算芯片提出了更高计算能力的要求，功率半导体器件不同于其它类型半导体器件，能够耐受高电压、大电流，是实现电路转换与信号传输的核心，对电能的高效处理、传输等起到了关键的作用. 与此同时，高功率15的运行带来的导电导热需求，以及高温的使用环境，也对封装材料有着更高的需求[1-2]. 传统高温钎料，如金锡共晶合金(Au-Sn)、锡银共晶合金(Sn-Ag)等在烧结过程中会出现多种金属间化合物(intermetallic compounds，IMCs)，使焊点可靠性下降，更重要的是上述钎料的热导率较低，限制了大功率器件的高温性能. 银是一种具有良好性能的材料，具有较高的导电导热性，但因其熔点过高，难以成为连接材料[3].
纳米银浆料是指尺寸大小在10 ~ 100 nm之间的银颗粒与多种有机物混合而成的一种连接材料[4]，具有较高的比表面积，因此具有表面效应、体积效应、量子尺寸效应、宏观量子隧道效应和介电限域效应[5]，由于极细的晶粒和大量处于晶界和晶粒缺陷中心的原子，纳米材料的物化性能与微米多晶材料有着巨大的差异，具有奇特的电学、热学和化学等多方面的性能，烧结后的浆料随着纳米银颗粒尺寸增大而熔点增高[6]，能够在高温环境下进行使用，实现低温烧结高温服役的特殊烧结与应用形式，并且烧结后的浆料呈现优于其它高温钎料的导电导热性[7]. 纳米银浆料互连材料在低温互连、多级互连和大功率芯片互连等方向受到愈来愈多的关注，然而银的导热性能良好，块体银的导热率达429 W/(K·m)，烧结后的纳米浆料在导热性能上却有显著降低.
针对于纳米银浆料的制备以及烧结后的导热性能变化已有许多人员进行了系统的研究. 刘彦立[8]通过仿真对比烧结后的纳米颗粒导热性能与纳米颗粒的孔隙率和孔隙形状等的关系，发现孔隙率是影响纳米导热性的主要因素；Wang等人通过不同烧结温度[9]、颗粒粒径[10]与烧结时间[11]的纳米浆料对比发现，通过不同尺寸的颗粒混合，提高烧结温度、增加烧结时间等方法能够有效降低孔隙率，从而提高纳米浆料的物理性能；Bai等人[12]通过对纳米浆料进行280 ℃烧结，热导率仅为238 W/(K·m)，而银材料的热导率为429 W/(K·m)，烧结过后的导热性能远没有达到预期程度；通过Li等人[10]的研究，采用双峰颗粒提高纳米浆料的致密度，也仅能使烧结后的纳米浆料热导率达到块体银的65%，造成烧结后的纳米浆料热导率远低于银材料的原因至今不够明确. 烧结过后颗粒烧结颈处的复杂结构可能对声子的传热有一定阻碍作用，从而降低了整体材料的热导率，但烧结过程中其微观结构尺寸较小，结构变化复杂，常规试验手段难以观察. 针对于纳米尺度原子的运动情况，已有学者[13-15]通过分子动力学仿真手段对其进行一定探索研究，通过仿真模拟指导实际的试验. 通过分子动力学仿真方法模拟纳米颗粒烧结过程中的结构变化，对烧结后的微观复杂结构计算其热导率，探究烧结颈处的缺陷结构是否对材料宏观性能造成一定影响.
1. 仿真方法
有大量学者针对纳米颗粒的熔化进行了分子动力学模拟，均取得了较为良好的成绩[13-15]. 分子动力学基于牛顿第二定律进行计算，将原子假设为实体小球，忽略原子的电子作用，考虑在相互作用下的运动情况，最终统计每一个原子的位移和速度从而计算其性能. 相比于第一性原理，分子动力学也是微观尺度的一种仿真方式，但是忽略了原子中电子的作用，对于一些只需考虑原子运动的仿真节约了大量算力，从而能够计算更大尺度上的仿真过程. 纳米银浆料的烧结过程是在高温的作用下原子扩散的行为，通过模拟原子的运动就可以较好地模拟其烧结过程，因此选择分子动力学进行仿真.
为了探究纳米颗粒烧结过程中的连接行为，建立两个直径为5 nm的纳米银颗粒的模型，为降低周期性边界条件中周围空间内颗粒对模拟的影响，设置颗粒到模拟环境边界的距离为20 nm，如图1所示，时间步长设置为0.001 ps，选择EAM势函数[16]，周期性边界，在NPT系综下进行1 000 K烧结.
在金属材料中，传热作用主要依靠电子和声子的共同作用，其表达式[17]为
$$ \kappa = {\kappa _{\rm{e}}} + {\kappa _{\rm{p}}} $$ (1) 式中：κ为材料的热导率；κe为电子振动传热产生的热导率；κp为声子振动传热产生的热导率.根据Wiedeman-Franz定律，即
$$ \frac{{{\kappa _{\rm{e}}}}}{\sigma } = {{L}}T $$ (2) 式中：σ为材料的电导率；L为洛伦兹常数，取值2.45 × 10−8 W·Ω/K2；T为材料的温度. 声子对导热的贡献主要受晶格结构影响，目前还没有确定的公式能够进行计算，通过非平衡分子动力学方法计算纳米颗粒声子热导率的部分. 该方法通过给材料两端施加恒定热源与冷端，制造温度梯度，模拟得到材料内部在传热过程中的温度分布，根据傅里叶定律计算，即
$$ {\kappa _{{{\rm{p}}}} } = - \frac{Q}{{{\rm{d}}T/{\rm{d}}x}} $$ (3) 式中：dT/dx为传热方向上的温度梯度；Q为热流密度.
通过硝酸银与硼氢化钠溶液氧化还原反应合成纳米银颗粒，并在管式炉中进行烧结，将烧结过后的颗粒分别在泰思肯离子双束电镜TESCAN AMBER和透射电镜FEI Talos F200X下进行电子显微镜 (scanning electron microscope, SEM) 和透射电子显微镜 (transmission electron microscope，TEM)扫描，观察其微观结构.
2. 结果与讨论
2.1 势函数的验证
分子动力学模拟微观与介观尺度的原子运动情况，往往难以使用试验的手段对仿真的正确性进行验证，而仿真的结果差异主要来源于势函数的选择. 在仿真开始之前，首先对势函数的正确性进行验证，根据主要纳米颗粒的烧结规律以及微观结构的导热率问题，在验证过程中主要针对分子动力学能否较好的模拟出材料在导热性能上的尺寸效应，并进行验证.
建立横截面积为5 nm × 5 nm，导热长度Lz为10 ~ 35 nm，根据傅里叶定律与Wiedeman-Franz定律分别计算材料常温下在x方向上的热导率，并得到图2所示的长度与热导率关系曲线. 可以看到随长度增长，银材料的热导率不断增大，但增长趋势不断减缓，这一现象与文献[18]计算Al/Ni材料热导率以及文献[19]计算Si/Ge超晶格界面热阻过程中观察到的现象一致，表现了材料在导热性能上的尺寸效应.
![]()
图 2 银材料热导率与导热长度关系Figure 2. Relationship between thermal conductivity of silver materials and heat transfer distance根据材料的有效声子平均自由程和材料传热长度的关系以及量子力学中材料的热导率近似公式，可以推导出如下关系[20]，即
$$ \frac{1}{{\kappa _{\text{p}}^{}}} = \frac{3}{{{{cv}}}}\left( {\frac{1}{{{{l}_\infty }}} + \frac{1}{{{L_{\text{z}}}}}} \right) $$ (4) 式中：v是声子速度； l∞为无限大系统内声子平均自由程， 为声子传热方向上长度. 式(4)表现了声子热导率的尺寸效应，将声子热导率的倒数与传热长度的倒数进行线性拟合，得到如图3所示的关系，通过线性拟合，表达式为
![]()
图 3 热导率倒数与导热长度倒数关系Figure 3. Reciprocal of thermal conductivity in relation to the reciprocal of the heat transfer distance$$ \frac{1}{\kappa } = 0.002\;32 + \frac{{0.001\;11}}{{{{{L}}_{\text{z}}}}} $$ (5) 通过计算可以得到银材料在无限长传热方向上的热导率为431.03 W/(K·m)，其与文献[21]中银材料热导率的试验值429 W/(K·m)的误差仅有0.5%.
使用嵌入原子势(embedded atom method, EAM)函数计算后的热导率与文献[20]的尺寸效应现象一致，计算得到的宏观热导率与试验值接近，证明所选择的势函数能够较为良好的完成仿真，并对研究纳米颗粒烧结后热导率变化情况具有一定的启示性作用.
2.2 晶格结构变化
将建模好颗粒在1 000 K下模拟烧结1 000 ps后进行降温，并在2 000 ps时观察降温后的稳定结构. 图4为烧结过程中0，25，250，500和1 000 ps时的颗粒形貌和位移云图. 从图中可以看出，随着烧结的进行，两个颗粒的原子都进行了相互扩散，颗粒形态由原来的双球状逐渐转变为椭球棒状.
![]()
图 4 不同烧结时间的颗粒形貌与位移云图Figure 4. Image of particle morphology and displacement at different sintering times. (a) 0 ps; (b) 25 ps; (c) 250 ps; (d) 500 ps; (e) 1 000 ps图4原子模拟过程中运动的距离，颜色越偏向于红色，表明原子距离开始位置越远. 可以看出，在烧结过程中，颗粒表面的原子颜色更偏向于红色，说明在颗粒表面的原子进行了充分的扩散，都运动到了距初始位点较远的位置，而内部的原子位移距离很小，填充两个颗粒之间缝隙的原子主要来自于两边颗粒的表层原子.
进一步分析观察在烧结过程中纳米颗粒的晶格结构变化，将后处理文件在Ovito软件中进行晶格识别，根据每一帧图像中原子间配位数判断其满足的晶格结构，绿色代表规则排列的面心立方结构，红色表示密排六方结构，白色为未标定结构. 根据与试验结果比对分析，仿真软件识别出的片层状红色颗粒并不是密排六方结构，而是在烧结颈形成后在两侧形成的晶界，因其特有的配位结构被软件误认为密排六方结构. 白色为未标定的晶格结构，在模拟过程中存在两种情况，一种是在升温过程中由于扩散脱离晶格位点的原子，而另一种是在仿真开始，处于颗粒表面的原子因为没有满足面心立方的配位数，而被判定为未标定的晶格结构，最后为了观察到颗粒内部的晶格结构，对纳米颗粒进行剖面处理，如图5所示.
![]()
图 5 不同烧结时间的纳米颗粒晶格结构图Figure 5. Lattice structure of nanoparticles at different sintering times. (a) 0 ps; (b) 25 ps; (c) 250 ps; (d) 500 ps; (e) 2 000 ps从图5中可以看出，在加热过程中，颗粒表面白色原子层变厚，原子吸收热能后转变为原子振动的动能，随机振动使原子失去了晶体的规则排列，表面的原子首先脱离了晶格位置，向表面能更小的位置进行扩散，填充到两个纳米颗粒中间的位置上，而颗粒内部仍保持较为完整的晶格结构，随着烧结过程的进行，在烧结颈处两边出现了晶界，纳米颗粒从原本的单一晶粒演变成多晶粒结构, 最后在冷却后形成面心立方的多晶结构，不同晶粒之间取向不同，晶界之间存在一定的夹角，而这一现象也在试验观察中得到了验证. 结合纳米颗粒烧结过程中原子的位移云图和晶格结构变化图，可以总结出形成烧结颈结构过程中主要的原子扩散方式是表面扩散.
图6是合成的纳米银颗粒分别在没有烧结前，200和250 ℃无压烧结1 h的SEM形貌. 从图中可以看出进行烧结后整体颗粒粒径变大，颗粒之间形成了明显的烧结颈结构，且随着烧结温度的提高，烧结颈的长度有明显的增长. 试验结果表明，在实际的烧结过程中存在两种颗粒连接的情况，对于不同尺寸的颗粒，大颗粒通过奥斯瓦尔德熟化现象融合小颗粒，颗粒本身粒径增长[22]；对于尺寸相近的颗粒，其生长符合仿真情况，将产生烧结颈结构，通过一定时间的烧结，烧结颈的长度趋向于略小于颗粒直径.
![]()
图 6 纳米颗粒烧结SEM形貌Figure 6. SEM image of sintered nanoparticles. (a) no sintering; (b) sintered at 200 ℃; (c) sintered at 250 ℃图7是250 ℃无压烧结1 h选取几个尺寸相近的颗粒进行TEM观察的形貌，从图中可以看出有多个颗粒都进行了良好的冶金结合，其中还出现了颗粒堆叠，形成了莫尔条纹. 观察红色虚线内的两个颗粒，通过表面的原子扩散填充到颗粒相接触的地方形成烧结颈，填充后的原子呈规则排列，面间距为0.233 nm，通过计算为银的面心立方(111)晶面的面间距. 黄色箭头指向的多晶结构与仿真结果相对应，晶界之间的夹角一一对应，通过查阅文献[23]，得出是银在纳米尺寸下的五重孪晶结构稳定状态，而烧结后的纳米颗粒内部同样出现了这一稳定结构.
![]()
图 7 纳米银颗粒烧结结构TEM形貌Figure 7. TEM diagram of the sintered structure of nano-silver particles图8为颗粒模拟烧结过程中烧结颈生长长度与时间关系，可以看到在烧结刚开始的阶段，烧结颈迅速增长，经过100 ps的烧结后，烧结颈已经具备一定长度，后续仍在增长，但增长的速度相较之前有所减少. 可以看出烧结颈的生长主要是在开始的阶段，原子迅速扩散，而当烧结颈结构生长到一定程度后，增速变缓，继续烧结将通过原子的滑移调整颗粒内部晶粒取向，最终形成多晶结构.
烧结颈生长规律这一现象也与粉末冶金中总结的烧结颈生长规律相似[15]，即
$$ {\left( {\frac{X}{{{R}}}} \right)^{\text{m}}}{\text{ = }}\frac{{{H}}}{{{{{R}}^{{n}}}}}t $$ (6) 式中：X表示烧结颈长度(nm)；m是与扩散机制有关的常数；R为颗粒半径(nm)；H和n都是与扩散有关的常数. 烧结颈长度和时间成幂函数关系，拟合得出的公式为
$$ X^{8.076}=1.926 t$$ (7) 拟合得到m的数值为8.076，这与传统的各种扩散机制的数值均不同，说明在纳米尺寸下的扩散机制与宏观的扩散行为有所不同，具体的扩散机制还有待进一步探究.
2.3 晶格结构变化对导热的影响
为探究烧结时间对纳米颗粒的热导率的影响，分别选取烧结前、烧结中和烧结完成的纳米颗粒模型，针对于两个球心的截面处施加平行于烧结颈的平面恒定热源与冷端，平衡导热后输出其烧结颈处的温度梯度分布，结合热导率计算公式，计算得到纳米颗粒烧结结构的热导率.
图9为烧结0，25和2 000 ps颗粒内部温度分布及晶格结构剖面. 从图中可以看出，在烧结前，两个颗粒烧结颈位置处具有温度突变，证明没有形成烧结结构的两个颗粒之间具有较大热阻，而结构图中也可以看到两个颗粒连接处较窄小，并没有形成完整的晶格结构，对热量的传递有很大阻碍作用. 而在烧结过程中，两边仍有温度平台，此时具备一定的烧结颈结构，烧结颈处的原子对导热具有一定贡献，但是仍有很多没有规则排列的白色原子，其对热量的传递仍有阻碍作用；当完全烧结后，热量传递过程中形成温度梯度，证明此时热量能够稳定地从热源传递到冷端，两个颗粒之间热阻显著降低. 烧结完全后计算得到纳米颗粒烧结颈处的热导率为410.73 W/(K·m)，与相同尺寸下规则排列的银材料的热导率仍有一定差距，其可能的原因是在烧结颈处形成了不同的晶格结构，阻碍了声子的热扩散运动，从而在一定程度上降低了热导率. 仿真中烧结后的两个纳米银颗粒烧结颈处的热导率相较于块体银下降了4.42%，然而文献[10]表征试验中纳米银浆料烧结后热导率出现大幅下降，甚至下降到了块体银热导率的65%，相较之下，烧结颈处的微观结构对热导率的阻碍作用并不明显.
![]()
图 9 不同烧结时间温度分布及晶格结构剖面图Figure 9. Temperature distribution and lattice structure profile at different sintering times. (a) 0 ps; (b) 25 ps; (c) 2 000 ps3. 结论
(1) 颗粒在烧结过程中表面的原子会扩散到颗粒连接的地方，形成面心立方结构的烧结颈，烧结后的浆料内部会倾向于形成五重孪晶的稳定多晶结构.
(2) 烧结颈的生长主要在初期阶段，烧结颈生长长度与时间关系满足幂函数关系，拟合公式X 8.076=1.926t，符合烧结颈生长经验公式的规律.
(3) 通过仿真3个时间节点的温度分布，发现烧结颈处的微观结构的确一定程度阻碍了热量的传递，热导率下降到410.73 W/(K·m)，相较于纯银块体下降了4.42%. 此外还发现通过完整的烧结可以降低纳米银颗粒之间的热阻.
-
![]()
图 2 银材料热导率与导热长度关系
Figure 2. Relationship between thermal conductivity of silver materials and heat transfer distance
![]()
图 3 热导率倒数与导热长度倒数关系
Figure 3. Reciprocal of thermal conductivity in relation to the reciprocal of the heat transfer distance
![]()
图 4 不同烧结时间的颗粒形貌与位移云图
Figure 4. Image of particle morphology and displacement at different sintering times. (a) 0 ps; (b) 25 ps; (c) 250 ps; (d) 500 ps; (e) 1 000 ps
![]()
图 5 不同烧结时间的纳米颗粒晶格结构图
Figure 5. Lattice structure of nanoparticles at different sintering times. (a) 0 ps; (b) 25 ps; (c) 250 ps; (d) 500 ps; (e) 2 000 ps
![]()
图 6 纳米颗粒烧结SEM形貌
Figure 6. SEM image of sintered nanoparticles. (a) no sintering; (b) sintered at 200 ℃; (c) sintered at 250 ℃
![]()
图 7 纳米银颗粒烧结结构TEM形貌
Figure 7. TEM diagram of the sintered structure of nano-silver particles
-
[1] Chen Hongtao, Hu Tianqi, Li Mingyu. Cu@Sn core–shell structure powder preform for high-temperature applications based on transient liquid phase bonding[J]. IEEE Transactions on Power Electronics, 2017, 32(1): 441 − 451. doi: 10.1109/TPEL.2016.2535365
[2] 杨东升, 张贺, 田艳红, 等. 电子封装微纳连接技术即失效行为研究进展[J]. 焊接学报, 2022, 43(11): 126 − 136. Yang Dongsheng, Zhang He, Tian Yanhong, et al. Progress in the study of micro-nano connection technology, i. e. failure behaviour, in electronic packaging[J]. Transactions of the China Welding Institution, 2022, 43(11): 126 − 136.
[3] Zhong Ying, An Rong, Wang Chunqing. Low temperature sintering Cu6Sn5 nanoparticles for superplastic and super-uniform high temperature circuit interconnections[J]. Small, 2015, 11(33): 4097 − 4103. doi: 10.1002/smll.201500896
[4] Liu Hao, Ma Rui, Wang Jianqiang, et al. Application of silver nanoparticles in electrically conductive adhesives with silver micro flakes[J]. China Welding, 2022, 31(1): 23 − 28.
[5] 李红, 袁俊丽, 粟卓新. 纳米连接过程的分子动力学模拟研究进展[J]. 中国机械工程, 2019, 30(4): 486 − 493. Li Hong, Yuan Junli, Li Zhuoxin. Advances in molecular dynamics simulation of nano-connectivity processes[J]. Chinese Journal of Mechanical Engineering, 2019, 30(4): 486 − 493.
[6] 裴淳. 一种基于空间点过程的银膏焊层等效弹性模量的仿真预测方法[J]. 电子元件与材料, 2020, 39(6): 80 − 83. doi: 10.14106/j.cnki.1001-2028.2020.06.014 Pei Chun. A spatial point process based simulation prediction method for equivalent elastic modulus of silver paste solder layers[J]. Electronic Components & Materials, 2020, 39(6): 80 − 83. doi: 10.14106/j.cnki.1001-2028.2020.06.014
[7] Carr J, Milhert X, Gadaud P, et al. Quantitative characterization of porosity and determination of elastic modules for sintered micro-silver joints[J]. Journal of Materials Processing Technology, 2015, 225: 19 − 23. doi: 10.1016/j.jmatprotec.2015.03.037
[8] 刘彦立. 基于孔洞特征的芯片互联烧结组织性能数值模拟[D]. 哈尔滨: 哈尔滨理工大学, 2021. Liu Yanli . Numerical simulation of chip interconnect sintering tissue properties based on hole characteristics[D]. Harbin: Harbin University of Science and Technology, 2021
[9] Wang M, Mei Y, Li X, et al. Die-attach on nickel substrate by pressureless sintering a trimodal silver paste[J]. Materials Letters, 2019, 253: 131 − 135. doi: 10.1016/j.matlet.2019.06.041
[10] Li M, Xiao Y, Zhang Z, et al. Bimodal sintered silver nanoparticle paste with ultrahigh thermal conductivity and shear strength for high temperature thermal interface material applications[J]. ACS Applied Materials & Interfaces, 2015, 7(17): 9157 − 9168.
[11] Ide E, Angata S, Hirose A, et al. Metal–metal bonding process using Ag metallo-organic nanoparticles[J]. Acta Materialia, 2005, 53(8): 2385 − 2393. doi: 10.1016/j.actamat.2005.01.047
[12] Bai J G, Zhiye Z Z, Calata J N, et al. Low-temperature sintered nanoscale silver as a novel semiconductor device-metallized substrate interconnect material[J]. IEEE Transactions on Components and Packaging Technologies, 2006, 29(3): 589 − 593. doi: 10.1109/TCAPT.2005.853167
[13] 贺晓斌. 金属纳米线连接分子动力学模拟[D]. 哈尔滨: 哈尔滨工业大学, 2013. He Xiaobin . Molecular dynamics simulation of metal nanowire connections[D]. Harbin: Harbin Institute of Technology, 2013.
[14] 吴茂, 常玲玲, 崔亚男. 纳米金颗粒融化与烧结过程的分子动力学模拟[J]. 粉末冶金材料科学与工程, 2013, 18(6): 775 − 782. Wu Mao, Chang Lingling, Cui Yanan. Molecular dynamics simulation of the melting and sintering process of gold nanoparticles[J]. Powder Metallurgy Materials Science and Engineering, 2013, 18(6): 775 − 782.
[15] Suzuki T, Terasaki T, Kawana Y, et al. Effect of manufacturing process on micro-deformation behavion of sintered-silver dia-attach material[J]. IEEE Transactions on Device and Material Reliability, 2016, 16: 588 − 596. doi: 10.1109/TDMR.2016.2614510
[16] Jiang Shan, Zhang Yuwen, Hao Peng. Molecular dynamics study of neck growth in laser sintering of hollow silver nanoparticles with different heating rates[J]. Journal of Physics D:Applied Physics, 2013, 46: 335302. doi: 10.1088/0022-3727/46/33/335302
[17] 王鑫, 鲁丹, 申胜平. Ni/Al 层合结构热传导性能的非平衡分子动力学研究[J]. 中国科学: 物理学 力学 天文学, 2014, 44: 506–513. Wang Xin , Lu Dan, Shen Shengping . Non-equilibrium molecular dynamics study of the thermal conductivity of Ni/Al laminated structures[J] Science China: Physics, Mechanics & Astronomy, 2014, 44: 506-513.
[18] Lange A P, Samanta A, Majidi H, et al. Dislocation mediated alignment during metal nanoparticle Coalescence[J]. Acta Materialia, 2016(120): 364 − 378.
[19] Balasubramanian G, Puri I K. Heat conduction across a solid-solid interface: Understanding nanoscale interfacial effects on thermal resistance[J]. Applied Physics Letters, 2011, 99(1): 013116. doi: 10.1063/1.3607477
[20] Gao Y F, Meng Q Y. Molecular dynamics simulation on thermal conductivity of one dimenison nanomaterials[J]. Acta Metallurgica Sinica, 2010, 46(10): 1244 − 1249.
[21] Sun Zhijian, Li Jiaxiong, Yu Michael. A review of the thermal conductivity of silver-epoxy nanocomposites as encapsulation material for packaging application[J]. Chemical Engineering Journal, 2022, 10(15): 137319.
[22] Fang Z Z, Wang H. Densification and grain growth during sintering of nanosized particles[J]. International Materials Reviews, 2013, 53(6): 326 − 352.
[23] 王尚. 银纳米线导电薄膜制备与电镀修饰及性能[D]. 哈尔滨: 哈尔滨工业大学, 2019. Wang Shang. Preparation and electroplating modification of silver nanowire conductive thin films and properties[D]. Harbin: Harbin Institute of Technology, 2019.
-
期刊类型引用(0)
其他类型引用(1)
下载:
计量
- 文章访问数: 449
- HTML全文浏览量: 180
- PDF下载量: 157
- 被引次数: 1
