Molecular dynamics and strengthening of liquid-crystal polymers

The role of large-scale molecular motion in the self-organization and strengthening of liquid-crystal polymer fibers is discussed. It is shown that, at high temperatures, these objects are oriented liquid-crystal melts in which macromolecules remain extended but execute high-frequency conformational motions without leaving the tube approximately 20 Å in diameter. This large-scale motion is referred to as quasi-segmental motion. During annealing, the chains involved in quasi-segmental motion can accomplish longitudinal displacements (reptate) over considerable distances. It is this reptation that favors spontaneous self-organization and, consequently, strengthening of liquid-crystal polymer fibers upon heat treatment. The role played by the quasi-segmental motion of rigid macromolecules in the strengthening of polymers of different types is compared with the role played by the segmental motion of flexible chains in this process.     

Molecular dynamics analysis of nanoimprinted Cu-Ni alloys

The nanoimprinted behaviors of Cu-Ni alloys were studied through molecular dynamics simulation. The results from the simulations showed that the punching force and the internal energy of the specimen increased rapidly with increasing punch position and both of they at the elevated temperature were lower than that at the low temperature. The punching force of the specimen with a lower Ni content was higher than that with a higher Ni content, while the internal energy of specimen was reverse. In addition, the spring-back phenomenon was more obvious with increasing Ni components of Cu-Ni alloy after the punch was retrieved during the nanoimprinting process and that induced a smaller residual stress within the specimens.     

Multiphase modelling of the growth kinetics of precipitates in Al-Cu alloys during artificial aging

ABSTRACT

In the aluminum industry, the needs of predictability of the kinetics of precipitation during the artificial aging processes increase as the targeted applications require the maximisation of properties at the lowest costs possible. In this regard, kinetics modelling can be helpful to design the heat treatment processes. Despite using many fitting parameters, available models show a lack of fitting with experimental data, especially for the apparent heat capacity measured at high temperatures by a differential scanning calorimeter (DSC). To address this issue, a mixed-mode model was recently developed for isothermal heat treatment, whereas non-isothermal heat treatment must be considered to compare the calculated results with those measured by isochronal heating in a DSC. In this contribution, the model is extended to non-isothermal heat treatments. To this end, the growth kinetics pathway and sequence of precipitation in a binary Al-Cu alloy have been simulated, optimising the pre-exponential factor and the activation energy of the interfacial mobility of the secondary phases. This calibration of the interfacial mobilities allowed a very good reproduction of the evolution of the apparent heat capacity with temperature. The model and calibrated interfacial mobilities were then used to compute the size evolution of θ′ precipitates in an Al-4?wt%Cu. The isothermal growth rates calculated at 4 temperatures were in good agreement with those measured and reported by independent researchers. The good predictability of the model indicates that the assumptions made were suitable and well funded, especially regarding the growth rates of embryos emerging from the subcritical growth regime.     

Accumulation of dislocations and thermal strengthening of alloys having an L12 superstructure

This paper discusses the dislocation-accumulation mechanism in alloys having an L1₂ superstructure, which is associated with the formation of Kira-Wilsdorf barriers and the retardation of superdislocations during plastic deformation. A model of the dislocation-accumulation kinetics during plastic deformation is constructed, on the basis of which a mathematical model is formulated for the thermal and deformation strengthening of single crystals of alloys having the L1₂ superstructure. The results of numerical calculations based on the model are compared with the experimentally observed regularities of the deformation and thermal strengthening of single crystals of Ni₃Ge. Fiz. Tverd. Tela (St. Petersburg) 41, 454–460 (March 1999)     

Molecular dynamics simulation of Cr-precipitate demixing in FeCr alloys

The mechanical properties of FeCr alloys rely heavily on atomic distribution and can be affected by phenomena such as Cr precipitation. While precipitation of FeCr alloys of various Cr concentrations has been studied before, dissolution of already existing Cr precipitates in FeCr alloys has not; this was the focus of this study. Our means of investigation was molecular dynamics computer simulation: we set up a number of configurations of FeCr alloys containing Cr precipitates of various sizes embedded in matrices of either pure Fe or with a 15% random Cr distribution, and examined their behaviour after thermal ageing at temperatures T ranging between 600 and 2000 K. The T range was selected so that it would include the (α+α^′)?α transition in the standard FeCr phase diagram. High-T results provide insight into the mechanisms that govern the dissolution: Cr precipitates dissolve by vacancy exchange, leading to a random distribution of Cr atoms in an Fe matrix, as the short-range order parameter shifts from a positive value (clustering) to zero (random atomic arrangement). Precipitates at low T were found to be stable, as were those at intermediate T (～1000 K), in agreement with previous experimental and simulation studies that challenge the standard phase diagram's reliability.     

Molecular dynamics study for the thermal conductivity of diatomic liquid

The thermal conductivity of diatomic liquids was analyzed using a nonequilibrium molecular dynamics (NEMD) method. Five liquids, namely, O₂, CO, CS₂, Cl₂ and Br₂, were assumed. The two-center Lennard-Jones (2CLJ) model was used to express the intermolecular potential acting on liquid molecules. First, the equation of state of each liquid was obtained using MD simulation, and the critical temperature, density and pressure of each liquid were determined. Heat conduction of each liquid at various liquid states [metastable (ρ=1.9ρ_cr), saturated (ρ=2.1ρ_cr), and stable (ρ=2.3ρ_cr)] at T=0.7T_cr was simulated and the thermal conductivity was estimated. These values were compared with experimental results and it was confirmed that the simulated results were consistent with the experimental data within 10%. Obtained thermal conductivities at saturated state were reduced by the critical temperature, density and mass of molecules and these values were compared with each other. It was found that the reduced thermal conductivity increased with the increase in the molecular elongation. Detailed analysis of the molecular contribution to the thermal conductivity revealed that the contribution of the heat flux caused by energy transport and by translational energy transfer to the thermal conductivity is independent of the molecular elongation while the contribution of the heat flux caused by rotational energy transfer to the thermal conductivity increases with the increase in the molecular elongation. Moreover, by comparing the reduced thermal conductivity at various states, it was found that the increase of thermal conductivity with the increase in the density, or pressure, was caused by the increase of the contribution of energy transfer due to molecular interaction.     

锆铌合金的特殊准随机结构模型的分子动力学研究

锆合金(如:锆铌(Zr-Nb)合金)的辐照损伤问题是裂变堆结构材料和燃料棒包壳材料设计的关键,而深入理解辐照损伤的物理机制,往往需借助于原子尺度的计算模拟,如:分子动力学和第一性原理等.针对随机置换固溶体合金的模拟,首先需构建能反映合金元素随机分布特征的大尺寸超胞,然而第一性原理计算量大,不宜使用过大(如 ≥200原子...     

Molecular dynamics simulation of thermal conductivity of silicone rubber

Silicone rubber is widely used as a kind of thermal interface material(TIM)in electronic devices.However few studies have been carried out on the thermal conductivity mechanism of silicone rubber.This paper investigates the thermal conductivity mechanism by non-equilibrium molecular dynamics(NEMD)in three aspects:chain length,morphology,and temperature.It is found that the effect of chain length on thermal conductivity varies with morphologies.In crystalline state where the chains are aligned,the thermal conductivity increases apparently with the length of the silicone-oxygen chain,the thermal conductivity of 79 nm-long crystalline silicone rubber could reach 1.49 W/(m·K).The thermal conductivity of amorphous silicone rubber is less affected by the chain length.The temperature dependence of thermal conductivity of silicone rubbers with different morphologies is trivial.The phonon density of states(DOS)is calculated and analyzed.The results indicate that crystalline silicone rubber with aligned orientation has more low frequency phonons,longer phonon MFP,and shorter conducting path,which contribute to a larger thermal conductivity.     

硅功能化石墨烯热导率的分子动力学模拟

采用Tersoff势函数与Lennard-Jones势函数,结合速度形式的Verlet算法和Fourier定律,对单层和两层硅功能化石墨烯沿长度方向的导热性能进行了正向非平衡态分子动力学模拟.通过模拟发现,硅原子的加入改变了石墨烯声子的模式、平均自由程和移动速度,使得单层硅功能化石墨烯模型的热导率随着硅原子数目的增加而急剧地减小.在300 K至1000 K温度变化范围内,单层硅功能化石墨烯的热导率呈下降趋势,具有明显的温度效应.对双层硅功能化石墨烯而言,少量的硅原子嵌入,起到了提高热导率的作用,但当硅原子数目达到一定数量后,材料的导热性能下降.     

Dislocation dynamics and thermal hardening of alloys with L12 superlattice

We studied the effect of superdislocation velocities on thermal hardening of alloys with the L1₂ superlattice. We estimate the activation energy for self-stopping of dislocations due to changes in the elastic field of moving dislocations as compared with immobile dislocations. We calculate the velocity of a dislocation loop element, taking into account interactions with dislocations. We show that, because of deformation stress, superdislocations can move at nearly the speed of sound. This, in turn, can help reduce the activation energy for self-stopping processes of superdislocations by a factor of 1.3 to 5.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 2, pp. 57–61, February, 1996.     

Molecular dynamics studies of thermal conductivity time correlation functions

The Green–Kubo time correlation function for the thermal conductivity in liquid argon is studied for a thermodynamic state close to the triple point by standard molecular dynamics simulations. The collective heat flux vector has been separated into contributions originated at the kinetic energy, the intermolecular potential and the pair virial function. Furthermore, the Green–Kubo time correlation functions have been broken down into partial n-body terms (n=1,2,3,4). The most important contribution to the thermal conductivity is represented by the auto correlation of the virial term. In contrast to other collective phenomena described by time correlation functions involving n-body terms, the partial Green–Kubo time correlation functions for the thermal conductivity are not affected by exponential long-time tails.     

Molecular dynamics studies on the thermal conductivity of single-walled carbon nanotubes

We have studied the thermal conductivity of single-walled carbon nanotubes (SWCNTs) using the NEMD method. The results indicate that the thermal conductivity values are not profoundly influenced by the specific simulation-technique used in the MD simulations. Some possible reasons, which could be responsible for the discrepancy on thermal conductivity values of SWCNTs in the literatures, are discussed.      

Long-term low-temperature ageing of amorphous alloys

Long term stability of the coercivity and⁵⁷Fe Mössbauer parameters (magnitude and average orientation of hyperfine induction) at ca 150 °C was tested at three TM-M type amorphous alloys (Fe₈₃B₁₇, Fe₄₀Ni₄₀B₂₀ and Co₆₆Fe₅Cr₇Si₈B₁₄) up to 5000 hrs. Correlations between the Mössbauer and magnetic quantities were found and their possible causes are discussed.     

CO_2水合物导热特性的分子动力学模拟

使用分子动力学模拟方法在NVT系综下对结构完整CO_2水合物以及结构缺陷CO_2水合物进行了导热模拟计算.对于结构完整的CO_2水合物,在200-230 K温度区间内,体系导热系数由0.4684 W·m~(-1)·K~(-1)变化到0.4836 W·m~(-1)·K~(-1),温度相关性较弱;而在230-280 K温度区间内,体系导热系数由0.4836 W·m~(-1)·K~(-1)变化到0.7494 W·m~(-1)·K~(-1),温度相关性变强;另外,通过计算功率图谱发现主体分子对水合物体系的导热贡献更大.对于结构缺陷CO_2水合物,发现晶穴占有率和笼形结构缺陷对体系导热均有一定影响,空笼晶胞导热系数约为完整晶胞导热系数的86.67%,体系的导热能力主要取决于主体结构的性质.     

Pt-Au核-壳结构纳米粒子热稳定性的分子动力学研究

采用分子动力学方法结合嵌入原子势, 对Pt-Au核-壳纳米粒子的热稳定性进行了研究. 计算结果表明: Pt-Au纳米粒子的熔点明显高于Au纳米粒子而低于Pt纳米粒子. 通过计算Lindemann指数发现: 壳层中的Au首先熔化, 然后逐渐向内部扩展, 最终导致核中的Pt完全熔化; 熔化所经历的温度区间明显宽于单质纳米粒子, 而且该熔化过程呈现典型的两阶段熔化特征; 在两次熔化之间, 存在着固(核)液(壳)共存的结构. 关键词： 纳米粒子 熔化 分子动力学     

Effect of thermal annealing on nanoimprinted Cu-Ni alloys using molecular dynamics simulation

The mechanical behaviors of nanoimprinted Cu-Ni alloys before and after annealing were studied using molecular dynamics simulation with a tight-binding potential. The results showed that when the punch is advancing, the punching force obtained from the simulation with a tight-binding potential is lower than with the Morse potential. During and after withdrawing the punch from the specimen, the adhesive phenomena are observed and the large residual stress in the Cu-Ni alloys is induced. During the annealing process, the internal energy of Cu-Ni alloys decreased with increasing the temperature and the component of Cu. In addition, comparing the maximum residual stress in the Cu-Ni alloys with and without annealing treatment, the stress is significantly released after annealing, especially in the higher component of Ni.     

Study of the phenomenon of recovery in ageing alloys

A structural study of the process of recovery in the ageing alloys Al-20wt.% Zn, and Al-40 wt.% Zn was carried out. It was shown that in recovery, diffusion processes take place in the alloy, which are associated with the growth of regions having the structure of '-phase, their dimensions being above the critical size, and with the dissolving of regions whose dimensions are below the critical size. It was established that as a result of cyclic recovery, localized plastic deformation regions are formed in the Al-20 wt.% Zn alloy.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, Vol. 16, No. 9, pp. 41–45, September, 1973.     

用分子动力学模拟甲烷水合物热激法分解

用分子动力学模拟方法研究甲烷水合物热激法分解，系统地研究注入340 K液态水的结构Ⅰ型甲烷水合物的分解机理.模拟显示水合物表层水分子与高温液态水分子接触获得热能，分子运动激烈，摆脱水分子间的氢键束缚，笼状结构被破坏.甲烷分子获得热能从笼中挣脱，向外体系扩散.热能通过分子碰撞从外层传递给内层水分子，水合物逐层分解.对比注入277K液态水体系模拟结果，得出热激法促进水合物分解. 关键词： 甲烷水合物 分子动力学模拟 热激法     

Cu50Ni50合金快速凝固的分子动力学模拟

采用分子动力学模拟的方法研究了Cu₅₀Ni₅₀合金在不同冷却速度下的凝固过程,利用均方位移、径向分布函数和结构可视化等方法分析其微观结构.并对凝固模型进行拉伸模拟,通过应力应变曲线和直观结构变化分析其性能.研究表明：冷却速度对Cu₅₀Ni₅₀合金凝固形成的结构有较大影响,随着冷却速度的升高,凝固形成的结构中晶体含量减少,在较低的冷却速度下,如冷却1×10¹²K/s时,Cu₅₀Ni₅₀合金凝固形成晶体结构;在较高的冷却速度下,如1×10¹⁴K/s时,Cu₅₀Ni₅₀合金凝固形成非晶体结构,且非晶Cu₅₀Ni₅₀合金的抗拉性能要优于晶体Cu₅₀Ni₅₀合金.