Skip to main content

molecular dynamics

Deformation of Top-Down and Bottom-Up Silver Nanowires

Submitted by amleach on

I wanted to share some our work on the deformation behavior of metal nanowires that was recently published in Advanced Functional Materials. In this work, we considered the tensile deformation of three experimentally observed silver nanowire geometries, including five-fold twinned, pentagonal nanowires. The manuscript abstract and urls to videos of the tensile deformation of the three nanowire geometries are below. A copy of the manuscript is attached.

Tension-Compression Asymmetry in Homogeneous Dislocation Nucleation

Submitted by Mark Tschopp on

Abstract. This letter addresses the dependence of homogeneous dislocation nucleation on the crystallographic orientation of pure copper under uniaxial tension and compression.  Molecular dynamics simulation results with an embedded-atom method potential show that the stress required for homogeneous dislocation nucleation is highly dependent on the crystallographic orientation and the uniaxial loading conditions; certain orientations require a higher stress in compression (e.g., <110> and <111>) and other orientations require a higher stress in tension (<100>).  Furthermore, the resolved shear stress in the slip direction is unable to completely capture the dependence of homogeneous dislocation nucleation on crystal orientation and uniaxial loading conditions.

Atomistic simulations for the evolution of a U-shaped dislocation in fcc Al

Submitted by Xiaoyan Li on

We show, through MD simulations, a new evolution pattern of the U-shaped dislocation in fcc Al that would enrich the FR mechanism. Direct atomistic investigation indicates that a U-shaped dislocation may behave in different manners when it emits the first dislocation loop by bowing out of an extended dislocation. One manner is that the glissile dislocation segment always bows in the original glide plane, as the conventional FR mechanism. Another is that non-coplanar composite dislocations appear owing to conservative motion of polar dislocation segments, and then bow out along each slip plane, creating a closed helical loop. The motion of these segments involves a cross-slip mechanism by which a dislocation with screw component moves from one slip plane into another. Ultimately, such non-coplanar evolution results in the formation of a FR source.

How can we obtain more information from protein structure?

Submitted by Cheng Zhu on

We know - or believe - protein function is determined by structure. Crystallographic and NMR studies can provide protein structures with atomic-level details at equilibrium. MD simulations can follow protein conformational changes in time with fs temporal resolution in the absence or presence of a bias mechanism, e.g., applied force, used to induce such changes.

Why is molecular mechanics simulation at 0K useful?

Submitted by Xi Chen on
Although it is more realistic to study the mechanical properties of nanostructures such as the carbon nanotubes (CNTs) at room temperature, atomistic simulations at finite temperature (such as molecular dynamics, MD) may cause the following problems: (1) Due to the limitation of the time scale achievable in MD (typically at the nanosecond scale), the loading rate in MD simulation at any finite temperature is not realistic. Very often, the loading rate used in MD simulations may well exceed 10m/s at 300K and thus many orders of magnitude higher than the real loading rate used in experiments. (2) A great advantage of simulation is to be able to turn on and turn off certain features and explore their effects, which is otherwise impossible in experiments. For example, the buckling behavior of CNTs is very sensitive to geometrical perturbations, which is prominent at room temperature and such perturbations causes severe uncertainties and makes it difficult to explore the intrinsic buckling behaviors. Therefore, by removing the temperature effect, we could better evaluate other key factors affecting the intrinsic buckling behavior, such as tube chirality, radius, and length, which could be otherwise covered by the thermal fluctuation effect. (3) Due to both time and length scale limitations, the MD simulations of large system are not yet possible, and thus the effective continuum models must be developed which need to be calibrated by atomistic simulations. At present, the temperature factor is still absent in most continuum models. Therefore, atomistic simulations at 0K or near 0K may provide a useful benchmark for the development of parallel continuum models, focusing on the most intrinsic and basic mechanical properties of nanostructures. Based on the above analysis, atomistic simulations at 0K by using the molecular mechanics (MM) method are still very useful, especially to us as mechanicians.

use NMA to get the elastic properties of loop

Submitted by Xi Chen on
(originally written by Yuye Tang
A key procedure of the molecular-dynamics decorated finite element method (MDeFEM) is to determine the effective properties of components of a macromolecule. Here I illustrate how could one use the NMA computed from MD to estimate the elastic properties of loops in mechanosensitive channels, which is related with my research.