Superplasticity in intermetallic NiAl nanowires via atomistic simulations

A novel superplastic deformation in an intermetallic B2-NiAl nanowire of cross-sectional dimensions of similar to 20 angstrom with failure strain as high as similar to 700% at 700 K temperature is reported. The minimum temperature under which the superplasticity has been observed is around 0.36 T-m, which is much lower than 0.5 T-m (T-m = melting temperature i.e. 1911 K for bulk B2-NiAl). Superplasticity is observed due to transformation from crystalline phase to amorphous phase after yielding of the nanowire. (C) 2010 Elsevier B.V. All rights reserved.

Melting of an Anchored Bilayer: Molecular Dynamics Simulations of the Structural Transition in (CnH2n+1NH3)(2)PbI4 (n=12, 14, 16, 18)

The thermally driven Structural phase transition in the organic-inorganic hybrid perovskite (CnH2n+1NH3)(2)PbI4 has been investigated using molecular dynamics (MD) simulations. This system consists of positively charged alkyl-amine chains anchored to a rigid negatively charged PbI4 sheet with the chains organized as bilayers with a herringbone arrangement. Atomistic simulations were performed using ail isothermal-isobaric ensemble over a wide temperature range from 65 to 665 K for different alkyl chain lengths, n = 12, 14, 16, and 18. The simulations are able to reproduce the essential Features of the experimental observations of this system, including the existence of a transition, the linear variation of the transition temperature with alkyl chain length, and the expansion of the bilayer thickness at the transition. By use of the distance fluctuation Criteria, it is Shown that the transition is associated With a Melting of the alkyl chains of the anchored bilayer. Ail analysis of the conformation of the alkyl chains shows increased disorder in the form of gauche defects above due melting transition. Simulations also show that the melting transition is characterized by the complete disappearance of all-trans alkyl chains in the anchored bilayer, in agreement with experimental observations. A conformationally disordered chain has a larger effective cross-sectional area, and above due transition a uniformly tilted arrangement of the anchored chains call no longer be Sustained. At the melt the angular distribution of the orientation of the chains are 110 longer uniform; the chains are splayed allowing for increased space for individual chains of the anchored bilayer. This is reflected in a sharp rise in the ratio of the mean head-to-head to tail-to-tail distance of the chains of the bilayer at the transition resulting in in expansion of the bilayer thickness. The present MD simulations provide a simple explanation as to how changes in conformation of individual alkyl-chains gives rise to the observed increase in the interlayer lattice spacing of (CnH2n+1NH3)(2)PbI4 at the melting transition.

Notch sensitivity in nanoscale metallic glass specimens: Insights from continuum simulations

Recent experiments have shown that nano-sized metallic glass (MG) specimens subjected to tensile loading exhibit increased ductility and work hardening. Failure occurs by necking as opposed to shear banding which is seen in bulk samples. Also, the necking is generally observed at shallow notches present on the specimen surface. In this work, continuum finite element analysis of tensile loading of nano-sized notched MG specimens is conducted using a thermodynamically consistent non-local plasticity model to clearly understand the deformation behavior from a mechanics perspective. It is found that plastic zone size in front of the notch attains a saturation level at the stage when a dominant shear band forms extending across the specimen. This size scales with an intrinsic material length associated with the interaction stress between flow defects. A transition in deformation behavior from quasi-brittle to ductile becomes possible when this critical plastic zone size is larger than the uncracked ligament length. These observations corroborate with atomistic simulations and experimental results. (C) 2015 Elsevier Ltd. All rights reserved.

Stress-induced martensitic phase transformation in Cu-Zr nanowires

A novel stress-induced martensitic phase transformation in an initial < 100 >/{100} B2-CuZr nanowire is reported for the first time in this letter. Such behavior is observed in a nanowire with cross-sectional dimensions of 19.44 x 19.44 angstrom(2) over a temperature range of 100-400 K and at a strain rate of 1 x 10(9) s(-1) using atomistic simulations. Phase transformation from an initial B2 phase to a BCT (Body-Centered-Tetragonal) phase is observed via nucleation and propagation of {100} twinning plane under high strain rate tensile deformation. (C) 2009 Elsevier B.V. All rights reserved.

Self-Organization of n-Alkane Chains in Water: Length Dependent Crossover from Helix and Toroid to Molten Globule

We demonstrate a chain length dependent crossover in the structural properties of linear hydrocarbon (n-alkane) chains using detailed atomistic simulations in explicit water. We identify a number of exotic structures of the polymer chain through energy minimization of representative snapshots collected from molecular dynamics trajectory. While the collapsed state is ring-like (circular) for small chains (CnH2n+2; n <= 20) and spherical for very long ones (n = 100), we find the emergence of ordered helical structures at intermediate lengths (n similar to 40). We find different types of disordered helices and toroid-like structures at n = 60. We also report a sharp transition in the stability of the collapsed state as a function of the chain length through relevant free energy calculations. While the collapsed state is only marginally metastable for C20H42, a clear bistable free energy surface emerges only when the chain is about 30 monomers long. For n = 30, the polymer exhibits an intermittent oscillation between the collapsed and the coil structures, characteristic of two stable states separated by a small barrier.

Structure of polyamidoamide dendrimers up to limiting generations: A mesoscale description

The polyamidoamide (PAMAM) class of dendrimers was one of the first dendrimers synthesized by Tomalia and co-workers at Dow. Since its discovery the PAMAMs have stimulated many discussions on the structure and dynamics of such hyperbranched polymers. Many questions remain open because the huge conformation disorder combined with very similar local symmetries have made it difficult to characterize experimentally at the atomistic level the structure and dynamics of PAMAM dendrimers. The higher generation dendrimers have also been difficult to characterize computationally because of the large size (294852 atoms for generation 11) and the huge number of conformations. To help provide a practical means of atomistic computational studies, we have developed an atomistically informed coarse-grained description for the PAMAM dendrimer. We find that a two-bead per monomer representation retains the accuracy of atomistic simulations for predicting size and conformational complexity, while reducing the degrees of freedom by tenfold. This mesoscale description has allowed us to study the structural properties of PAMAM dendrimer up to generation 11 for time scale of up to several nanoseconds. The gross properties such as the radius of gyration compare very well with those from full atomistic simulation and with available small angle x-ray experiment and small angle neutron scattering data. The radial monomer density shows very similar behavior with those obtained from the fully atomistic simulation. Our approach to deriving the coarse-grain model is general and straightforward to apply to other classes of dendrimers.

Atomic Scale Fluctuations Govern Brittle Fracture and Cavitation Behavior in Metallic Glasses

We perform atomistic simulations on the fracture behavior of two typical metallic glasses, one brittle (FeP) and the other ductile (CuZr), and show that brittle fracture in the FeP glass is governed by an intrinsic cavitation mechanism near crack tips in contrast to extensive shear banding in the ductile CuZr glass. We show that a high degree of atomic scale spatial fluctuations in the local properties is the main reason for the observed cavitation behavior in the brittle metallic glass. Our study corroborates with recent experimental observations of nanoscale cavity nucleation found on the brittle fracture surfaces of metallic glasses and provides important insights into the root cause of the ductile versus brittle behavior in such materials.

Non-monotonic, distance-dependent relaxation of water in reverse micelles: Propagation of surface induced frustration along hydrogen bond networks

Layer-wise, distance-dependent orientational relaxation of water confined in reverse micelles (RM) is studied using theoretical and computational tools. We use both a newly constructed ``spins on a ring'' (SOR) Ising-type model (with Shore-Zwanzig rotational dynamics) and atomistic simulations with explicit water. Our study explores the effect of reverse micelle size and role of intermolecular correlations, compromised by the presence of a highly polar surface, on the distance (from the interface) dependence of water relaxation. The ``spins on a ring'' model can capture some aspects of distance dependence of relaxation, such as acceleration of orientational relaxation at intermediate layers. In atomistic simulations, layer-wise decomposition of hydrogen bond formation pattern clearly reveals that hydrogen bond arrangement of water at a certain distance away from the surface can remain frustrated due to the interaction with the polar surface head groups. This layer-wise analysis also reveals the presence of a non-monotonic slow relaxation component which can be attributed to this frustration effect and which is accentuated in small to intermediate size RMs. For large size RMs, the long time component decreases monotonically from the interface to the interior of the RMs with slowest relaxation observed at the interface. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4732095]

Interaction of nucleic acids with carbon nanotubes and dendrimers

Nucleic acid interaction with nanoscale objects like carbon nanotubes (CNTs) and dendrimers is of fundamental interest because of their potential application in CNT separation, gene therapy and antisense therapy. Combining nucleic acids with CNTs and dendrimers also opens the door towards controllable self-assembly to generate various supra-molecular and nano-structures with desired morphologies. The interaction between these nanoscale objects also serve as a model system for studying DNA compaction, which is a fundamental process in chromatin organization. By using fully atomistic simulations, here we report various aspects of the interactions and binding modes of DNA and small interfering RNA (siRNA) with CNTs, graphene and dendrimers. Our results give a microscopic picture and mechanism of the adsorption of single- and double-strand DNA (ssDNA and dsDNA) on CNT and graphene. The nucleic acid-CNT interaction is dominated by the dispersive van der Waals (vdW) interaction. In contrast, the complexation of DNA (both ssDNA and dsDNA) and siRNA with various generations of poly-amido-amine (PAMAM) dendrimers is governed by electrostatic interactions. Our results reveal that both the DNA and siRNA form stable complex with the PAMAM dendrimer at a physiological pH when the dendrimer is positively charged due to the protonation of the primary amines. The size and binding energy of the complex increase with increase in dendrimer generation. We also give a summary of the current status in these fields and discuss future prospects.

Cavitation in materials with distributed weak zones: Implications on the origin of brittle fracture in metallic glasses

Several experimental studies have shown that fracture surfaces in brittle metallic glasses (MGs) generally exhibit nanoscale corrugations which may be attributed to the nucleation and coalescence of nanovoids during crack propagation. Recent atomistic simulations suggest that this phenomenon is due to large spatial fluctuations in material properties in a brittle MG, which leads to void nucleation in regions of low atomic density and then catastrophic fracture through void coalescence. To explain this behavior, we propose a model of a heterogeneous solid containing a distribution of weak zones to represent a brittle MG. Plane strain continuum finite element analysis of cavitation in such an elastic-plastic solid is performed with the weak zones idealized as periodically distributed regions having lower yield strength than the background material. It is found that the presence of weak zones can significantly reduce the critical hydrostatic stress for the onset of cavitation which is controlled uniquely by the local yield properties of these zones. Also, the presence of weak zones diminishes the sensitivity of the cavitation stress to the volume fraction of a preexisting void. These results provide plausible explanations for the observations reported in recent atomistic simulations of brittle MGs. An analytical solution for a composite, incompressible elastic-plastic solid with a weak inner core is used to investigate the effect of volume fraction and yield strength of the core on the nature of cavitation bifurcation. It is shown that snap-cavitation may occur, giving rise to sudden formation of voids with finite size, which does not happen in a homogeneous plastic solid. (c) 2012 Elsevier Ltd. All rights reserved.

A meshless adaptive multiscale method for fracture

The paper presents a multiscale method for crack propagation. The coarse region is modelled by the differential reproducing kernel particle method. Fracture in the coarse scale region is modelled with the Phantom node method. A molecular statics approach is employed in the fine scale where crack propagation is modelled naturally by breaking of bonds. The triangular lattice corresponds to the lattice structure of the (111) plane of an FCC crystal in the fine scale region. The Lennard-Jones potential is used to model the atom-atom interactions. The coupling between the coarse scale and fine scale is realized through ghost atoms. The ghost atom positions are interpolated from the coarse scale solution and enforced as boundary conditions on the fine scale. The fine scale region is adaptively refined and coarsened as the crack propagates. The centro symmetry parameter is used to detect the crack tip location. The method is implemented in two dimensions. The results are compared to pure atomistic simulations and show excellent agreement. (C) 2014 Elsevier B. V. All rights reserved.

Cavitation in brittle metallic glasses - Effects of stress state and distributed weak zones

Experimental studies and atomistic simulations have shown that brittle metallic glasses fail by a cavitation mechanism whose origin has been traced to the presence of intrinsic atomic density fluctuations which give rise to weak zones with reduced yield strength. It has been shown recently through continuum analysis that the presence of these zones can lower the cavitation stress considerably under equibiaxial loading. The objective of the present work is to study the effect of the applied stress state on the cavitation behavior of such a heterogeneous plastic solid with distributed weak zones. To this end, 2D plane strain finite element simulations are performed by subjecting a unit cell containing a weak zone to different (biaxiality) stress ratios. The volume fraction and yield strength of the weak zone are varied over a wide range. The results show that unlike in a homogeneous plastic solid, the cavitation stress of the heterogeneous aggregate does not reduce appreciably as the stress ratio decreases from unity when the yield strength of the weak zone is low. It is found that a non-dimensional parameter characterizing the stress state prevailing in the weak zone and its yield properties uniquely control the cavitation stress. The nature of cavitation bifurcation may change from unstable bifurcation to the left at sufficiently low stress ratio to one involving snap cavitation at high stress ratio. (C) 2014 Elsevier Ltd. All rights reserved.

Prospects of zero Schottky barrier height in a graphene-inserted MoS2-metal interface

A low Schottky barrier height (SBH) at source/drain contact is essential for achieving high drive current in atomic layer MoS(2-)channel-based field effect transistors. Approaches such as choosing metals with appropriate work functions and chemical doping are employed previously to improve the carrier injection from the contact electrodes to the channel and to mitigate the SBH between the MoS2 and metal. Recent experiments demonstrate significant SBH reduction when graphene layer is inserted between metal slab (Ti and Ni) and MoS2. However, the physical or chemical origin of this phenomenon is not yet clearly understood. In this work, density functional theory simulations are performed, employing pseudopotentials with very high basis sets to get insights of the charge transfer between metal and monolayer MoS2 through the inserted graphene layer. Our atomistic simulations on 16 different interfaces involving five different metals (Ti, Ag, Ru, Au, and Pt) reveal that (i) such a decrease in SBH is not consistent among various metals, rather an increase in SBH is observed in case of Au and Pt; (ii) unlike MoS2-metal interface, the projected dispersion of MoS2 remains preserved in any MoS2-graphene- metal system with shift in the bands on the energy axis. (iii) A proper choice of metal (e.g., Ru) may exhibit ohmic nature in a graphene-inserted MoS2-metal contact. These understandings would provide a direction in developing high-performance transistors involving heteroatomic layers as contact electrodes. (c) 2016 AIP Publishing LLC.

Understanding the Growth of Interfacial Reaction Product Layers between Dissimilar Materials

When one starts to analyze the evolution of the interfacial reaction product layers between dissimilar materials it is often found out that as the number of interacting species grows, the complexity of the analysis increases extremely rapidly. It may even appear that the task is just too difficult to be completed. In this article we present the thermodynamic-kinetic method, which can be used to rationalize the evolution of interfacial reaction layers and bring back the physics to the analyses. The method is conceptually very simple. It combines energetics-what can happen-with kinetics-how fast things take place. Yet the method is flexible enough that it can utilize quantitative and qualitative data starting from the atomistic simulations up to the experiments carried out with bulk materials. Several examples about how to utilize this method in material scientific problems are given.