Because of its very large c / a ratio, zinc has proven to be a difficult element to model using semi-empirical classical potentials. It has been shown, in particular, that for the modified embedded atom method (MEAM), a potential cannot simultaneously have an hcp ground state and c / a ratio greater than ideal. As an alloying element, however, useful zinc potentials can be generated by relaxing the condition that hcp be the lowest energy structure. In this paper, we present a MEAM zinc potential, which gives accurate material properties for the pure state, as well as a MEAM ternary potential for the Mg–Al–Zn system which will allow the atomistic modeling of a wide class of alloys containing zinc. The effects of zinc in simple Mg–Zn for this potential is demonstrated and these results verify the accuracy for the new potential in these systems.
To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) -Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.
27 June 2018 in Modelling and Simulation in Materials Science and Engineering
In this paper, we demonstrate that previously observed β to α transitions for titanium interatomic potentials available in the literature arose from a mechanical instability and thus the potentials underestimated the correct thermodynamic phase transition temperature by hundreds of degrees. Using a relative free energy method for the two phases to calculate the true transition temperature, we present a new modified embedded atom method potential for titanium that shows a transition temperature of 1155 ± 2 K in excellent agreement with the experimentally observed transition. This free energy approach avoids the problems of irreversibility which occur when one relies on direct observation of the phase transition in molecular dynamics simulation. Other transformation mechanisms in addition to the mechanical instability are also considered. Finally, the new potential predicts the proper c -axis plastic twinning for titanium under compression making it the only potential that correctly predicts the phase transition temperature and the plastic behavior of α Ti.
Despite their wide spread applications, the mechanical behavior of helically coiled structures has evaded an accurate understanding at any length scale (nano to macro) mainly due to their geometrical complexity. The advent of helically coiled micro/nanoscale structures in nano-robotics, nano-inductors, and impact protection coatings has necessitated the development of new methodologies for determining their shear and tensile properties. Accordingly, we developed a synergistic protocol which (i) integrates analytical, numerical (i.e., finite element using COMSOL[reg]) and experimental (harmonic detection of resonance; HDR) methods to obtain an empirically validated closed form expression for the shear modulus and resonance frequency of a singly clamped helically coiled carbon nanowire (HCNW), and (ii) circumvents the need for solving 12th order differential equations. From the experimental standpoint, a visual detection of resonances (using in situ scanning electron microscopy) combined with HDR revealed intriguing non-planar resonance modes at much lower driving forces relative to those needed for linear carbon nanotube cantilevers. Interestingly, despite the presence of mechanical and geometrical nonlinearities in the HCNW resonance behavior the ratio of the first two transverse modes f2/f1 was found to be similar to the ratio predicted by the Euler-Bernoulli theorem for linear cantilevers.
6 November 2018 in The Journal of Physical Chemistry A