Despite the fundamental nature and practical importance of melting, there is still no generally accepted theory of this ubiquitous phenomenon. Even the earliest simulations of melting of hard discs by Alder and Wainwright indicated the active role of collective atomic motion in melting and here we utilize molecular dynamics simulation to determine whether these correlated motions are similar to those found in recent studies of glass-forming (GF) liquids and other condensed, strongly interacting, particle systems. We indeed find string-like collective atomic motion in our simulations of "superheated" Ni crystals, but other observations indicate significant differences from GF liquids. For example, we observe neither stretched exponential structural relaxation, nor any decoupling phenomenon, while we do find a boson peak, findings that have strong implications for understanding the physical origin of these universal properties of GF liquids. Our simulations also provide a novel view of "homogeneous" melting in which a small concentration of interstitial defects exerts a powerful effect on the crystal stability through their initiation and propagation of collective atomic motion. These relatively rare point defects are found to propagate down the strings like solitons, driving the collective motion. Crystal integrity remains preserved when the permutational atomic motions take the form of ring-like atomic exchanges, but a topological transition occurs at higher temperatures where the rings open to form linear chains similar in geometrical form and length distribution to the strings of GF liquids. The local symmetry breaking effect of the open strings apparently destabilizes the local lattice structure and precipitates crystal melting. The crystal defects are thus not static entities under dynamic conditions, such as elevated temperatures or material loading, but rather are active agents exhibiting a rich nonlinear dynamics that is not addressed in conventional "static" defect melting models.

A novel method to calculate the solid-liquid contact angle is introduced in this study. Using the 3D configuration of a liquid droplet on a solid surface, this method calculates the contact angle along the contact line and provides an angular distribution. Although this method uses the 3D configuration of liquid droplets, it does not require the calculation of the 3D density profile to identify the boundaries of the droplet. This decreases the computational cost of the contact angle calculation greatly. Moreover, no presumption about the shape of the liquid droplet is needed when using the method introduced in this study. Using this method, the relationship between the size and the contact angle of water nano-droplets on a graphite substrate was studied. It is shown that the contact angle generally decreases by increasing the size of the nano-droplet. The microscopic contact angle of 83.0° was obtained for water on graphite which is in a good agreement with previous experimental and numerical studies. Neglecting other nanoscale effects which may influence the contact angle, the line tension of SPC/E (extended simple point charge model) water was calculated to be 3.6×10 N, which is also in good agreement with the previously calculated values.

Recently, ZnS quantum dots have attracted a lot of attention since they can be a suitable alternative for cadmium-based quantum dots, which are known to be highly carcinogenic for living systems. However, the structural stability of nanocrystalline ZnS seems to be a challenging issue since ZnS nanoparticles have the potential to undergo uncontrolled structural change at room temperature. Using the molecular dynamics technique, we have studied the structural evolution of 1 to 5 nm freestanding ZnS nanoparticles with zinc-blende and wurtzite crystal structures. Simulation results revealed that relaxed configurations of ZnS nanoparticles larger than 3 nm consist of three regions: a) a crystalline core, b) a distorted network of 4-coordinated atoms environing the crystalline core, and c) a surface structure made entirely of 3-coordinated atoms. Decreasing the size of ZnS nanoparticle to 2 nm will cause the crystalline core to disappear. Further reducing the size will cause all of the atoms to become 3-coordinated. Dipole moments of zinc-blende and wurtzite nanoparticles are in the same range when the nanoparticles are smaller than 3 nm. Increasing the size makes dipole moments converge to the bulk values. This makes zinc-blende and wurtzite nanoparticles less and more polar, respectively.

Suitable optoelectronic properties and the nontoxic nature of ZnS quantum dots capacitate exciting applications for these nanomaterials especially in the field of biomedical imaging. However, the structural stability of ZnS nanoparticles has been shown to be challenging since they potentially are prone to autonomous structural evolutions in ambient conditions. Thus, it is essential to build an understanding about the structural evolution of ZnS nanoparticles, especially in aqueous environment, before implementing them for in vivo applications. In this study we compared the structure of ZnS nanoparticles relaxed in a vacuum and in water using a classical molecular dynamics method. Structural analyses showed that the previously observed threephase structure of bare nanoparticles is not formed in the hydrated state. The bulk of hydrated nanoparticles has more crystalline structure; however, the dynamic heterogeneity in their surface relaxation makes them more polar compared to bare nanoparticles. This heterogeneity is more severe in hydrated wurtzite nanoparticles, causing them to show larger dipole moments. Analyzing the structure of water in the first hydration shell of the surface atoms shows that water is mainly adsorbed to the nanoparticles' surface through Zn−O interaction, which causes the structure of water in the first hydration shell to be discontinuous.

No abstract