MOLECULAR DYNAMICS STUDY OF THE COALESCENCE OF EQUAL AND UNEQUAL SIZED Cu NANOPARTICLES

2009 ◽  
Vol 20 (02) ◽  
pp. 179-196 ◽  
Author(s):  
H. H. KART ◽  
G. WANG ◽  
I. KARAMAN ◽  
T. ÇAĞIN

Molecular dynamics simulations technique is used to study the consolidation of two nanoparticles of Cu element. We have studied sintering processes of two nanoparticles at different temperatures. Two model systems with 4 and 10 nm diameter of particles are selected to study the sintering process of the two nanoparticles. Orientation effects on the physical properties of consolidation of two nanoparticles with respect to each other are investigated. Temperature effects on the consolidation of two nanoparticles are also studied. The order of the values obtained in the simulation for the constant volume heat capacity and latent heat of fusion is good agreement with the bulk results. Moreover, we have investigated the size effects on the consolidation of two different sizes of nanoparticles, that is, one particle of diameter with 10 nm is fixed while the other one is changing from 1 to 10 nm. Melting temperatures of the copper nanoparticles are found to be decreased as the size of the particle decreases. It is found that simulation results are compatible with the other theoretical calculations.

Author(s):  
Toshihiro Kaneko ◽  
Kenji Yasuoka ◽  
Ayori Mitsutake ◽  
Xiao Cheng Zeng

Multicanonical molecular dynamics simulations are applied, for the first time, to study the liquid-solid and solid-solid transitions in Lennard-Jones (LJ) clusters. The transition temperatures are estimated based on the peak position in the heat capacity versus temperature curve. For LJ31, LJ58 and LJ98, our results on the solid-solid transition temperature are in good agreement with previous ones. For LJ309, the predicted liquid-solid transition temperature is also in agreement with previous result.


Author(s):  
Adrian Dominguez-Castro ◽  
Thomas Frauenheim

Theoretical calculations are an effective strategy to comple- ment and understand experimental results in atomistic detail. Ehrenfest molecular dynamics simulations based on the real-time time-dependent density functional tight-binding (RT-TDDFTB) approach...


2009 ◽  
Vol 60-61 ◽  
pp. 315-319 ◽  
Author(s):  
W.W. Zhang ◽  
Qing An Huang ◽  
H. Yu ◽  
L.B. Lu

Molecular dynamics simulations are carried out to characterize the mechanical properties of [001] and [110] oriented silicon nanowires, with the thickness ranging from 1.05nm to 3.24 nm. The nanowires are taken to have ideal surfaces and (2×1) reconstructed surfaces, respectively. A series of simulations for square cross-section Si nanowires have been performed and Young’s modulus is calculated from energy–strain relationship. The results show that the elasticity of Si nanowires is strongly depended on size and surface reconstruction. Furthermore, the physical origin of above results is analyzed, consistent with the bond loss and saturation concept. The results obtained from the molecular dynamics simulations are in good agreement with the values of first-principles. The molecular dynamics simulations combine the accuracy and efficiency.


1988 ◽  
Vol 141 ◽  
Author(s):  
Laurent J. Lewis ◽  
Normand Mousseau ◽  
FranÇois Drolet

AbstractA new algorithm for generating fully-coordinated hydrogenated amorphous silicon models with periodic boundary conditions is presented. The hydrogen is incorporated into an a-Si matrix by a bond-switching process similar to that proposed by Wooten, Winer, and Weaire, making sure that four-fold coordination is preserved and that no rings with less than 5 members are created. After each addition of hydrogen, the structure is fully relaxed. The models so obtained, to be used as input to molecular dynamics simulations, are found to be in good agreement with experiment. A model with 12 at.% H is discussed in detail.


1995 ◽  
Vol 408 ◽  
Author(s):  
P. W.-C. Kung ◽  
J. T. Books ◽  
C. M. Freeman ◽  
S. M. Levine ◽  
B. Vessali ◽  
...  

AbstractWe have used constant pressure molecular dynamics calculations to explore the behavior at various temperatures of two molecular crystals: benzene and a brominated phenyl compound. We observed a melting transition by heating the crystals from a low temperature. In the case of benzene, we performed one heating run of about 1 ns and obtained agreement with the experimental melting point to within some 8%. We have also simulated the melting of a more complex molecular crystal that contains bromine and phenyl groups. We performed four heating runs, with different rates of heating. For total simulation times of about 100, 220, 770, and 1 I50ps, the heating runs predicted melting temperatures that differed from the experimental melting temperature by 53%, 33%, 25%, and 9% respectively.


2013 ◽  
Vol 9 ◽  
pp. 118-134 ◽  
Author(s):  
Jutta Erika Helga Köhler ◽  
Nicole Grczelschak-Mick

Four highly ordered hydrogen-bonded models of β-cyclodextrin (β-CD) and its inclusion complex with benzene were investigated by three different theoretical methods: classical quantum mechanics (QM) on AM1 and on the BP/TZVP-DISP3 level of approximation, and thirdly by classical molecular dynamics simulations (MD) at different temperatures (120 K and 273 to 300 K). The hydrogen bonds at the larger O2/O3 rim of empty β-CDs prefer the right-hand orientation, e.g., O3-H…O2-H in the same glucose unit and bifurcated towards …O4 and O3 of the next glucose unit on the right side. On AM1 level the complex energy was −2.75 kcal mol−1 when the benzene molecule was located parallel inside the β-CD cavity and −2.46 kcal mol−1 when it was positioned vertically. The AM1 HOMO/LUMO gap of the empty β-CD with about 12 eV is lowered to about 10 eV in the complex, in agreement with data from the literature. AM1 IR spectra displayed a splitting of the O–H frequencies of cyclodextrin upon complex formation. At the BP/TZVP-DISP3 level the parallel and vertical positions from the starting structures converged to a structure where benzene assumes a more oblique position (−20.16 kcal mol−1 and −20.22 kcal mol−1, resp.) as was reported in the literature. The character of the COSMO-RS σ-surface of β-CD was much more hydrophobic on its O6 rim than on its O2/O3 side when all hydrogen bonds were arranged in a concerted mode. This static QM picture of the β-CD/benzene complex at 0 K was extended by MD simulations. At 120 K benzene was mobile but always stayed inside the cavity of β-CD. The trajectories at 273, 280, 290 and 300 K certainly no longer displayed the highly ordered hydrogen bonds of β-CD and benzene occupied many different positions inside the cavity, before it left the β-CD finally at its O2/O3 side.


2021 ◽  
Author(s):  
Xue-Qi Lv ◽  
Xiong-Ying Li

Abstract The melting at the magnesium/aluminum (Mg/Al) interface is an essential step during the fabrications of Mg-Al structural materials and biomaterials. We carried out molecular dynamics simulations on the melting at the Mg/Al interface in a Mg-Al-Mg nanolayer via analyzing the changes of average atomic potential energy, Lindemann index, heat capacity, atomic density distribution and radial distribution function with temperature. The melting temperatures (T m) of the nanolayer and the slabs near the interface are significantly sensitive to the heating rate (v h) over the range of v h≤4.0 K/ps. The distance (d) range in which the interface affects the melting of the slabs is predicted to be (-98.2, 89.9) Å at v h→0, if the interface is put at d=0 and Mg (Al) is located at the left (right) side of the interface. The (T m) of the Mg (Al) slab just near the interface (e.g., d=4.0 Å) is predicted to be 926.8 K (926.6 K) at v h→0, with 36.9 K (37.1 K) below 963.7 K for the nanolayer. These results highlight the importance of regional research on the melting at an interface in the nanolayers consisting of two different metals.


MRS Advances ◽  
2020 ◽  
Vol 5 (27-28) ◽  
pp. 1475-1482
Author(s):  
Marco A. Maria ◽  
Alexandre F. Fonseca

ABSTRACTSarin gas is a dangerous chemical warfare agent (CWA). It is a nerve agent capable of bringing a person to death in about 15 minutes. A lethal concentration of sarin molecules in air is about 30 mg/m3. Experimental research on this gas requires very careful safety protocols for handling and storage. Therefore, theoretical and computational studies on sarin gas are very welcome and might provide important safe guides towards the management of this lethal substance. In this work, we investigated the interactions between sarin, air and nanoporous graphene, using tools of classical molecular dynamics simulations. Aiming to cast some light in the possible sarin selective filtration by graphene, we designed a bipartite simulation box with a porous graphene nanosheet placed at the middle. Sarin and air molecules were initially placed only on one side of the box so as to create an initial pressure towards the passage of both to the other side. The box dimensions were chosen so that the hole in the graphene was the only possible way through which sarin and air molecules can get to the other side of the box. The number of molecules that passed through the hole in graphene was monitored during 10 ns of simulation and the results for different temperatures were compared. The results show that, as far as the size of the holes are small, van der Waals forces between graphene and the molecules play a significant role on keeping sarin near graphene, at room temperature.


Sign in / Sign up

Export Citation Format

Share Document