Atomistic investigation on the structure-property relationship during thermal spray nanoparticle impact

During thermal spraying, hot particles impact on a colder substrate. This interaction of crystalline copper nanoparticles and copper substrate is modelled, using MD simulation. The quantitative results of the impacts at different velocities and temperatures are evaluated using a newly defined flattening aspect ratio. This ratio between the maximum diameter after the impact and the height of the splat increases with increasing Reynolds numbers until a critical value is reached. At higher Reynolds numbers the flattening aspect ratio decreases again, as the kinetic energy of the particle leads to increasing substrate temperature and, therefore, decreases the substrate resistance. Thus, the particle penetrates into the substrate and deforms less.

A theoretical assessment of surface defect machining and hot machining of nanocrystalline silicon carbide

In this paper, a newly proposed machining method named “surface defect machining” (SDM) [Wear, 302, 2013 (1124-1135)] was explored for machining of nanocrystalline beta silicon carbide (3C-SiC) at 300K using MD simulation. The results were compared with isothermal high temperature machining at 1200K under the same machining parameters, emulating ductile mode micro laser assisted machining (µ-LAM) and with conventional cutting at 300 K. In the MD simulation, surface defects were generated on the top of the (010) surface of the 3C-SiC work piece prior to cutting, and the workpiece was then cut along the <100> direction using a single point diamond tool at a cutting speed of 10 m/sec. Cutting forces, sub-surface deformation layer depth, temperature in the shear zone, shear plane angle and friction coefficient were used to characterize the response of the workpiece. Simulation results showed that SDM provides a unique advantage of decreased shear plane angle which eases the shearing action. This in turn causes an increased value of average coefficient of friction in contrast to the isothermal cutting (carried at 1200 K) and normal cutting (carried at 300K). The increase of friction coefficient however was found to aid the cutting action of the tool due to an intermittent dropping in the cutting forces, lowering stresses on the cutting tool and reducing operational temperature. Analysis shows that the introduction of surface defects prior to conventional machining can be a viable choice for machining a wide range of ceramics, hard steels and composites compared to hot machining.

Influence of test methodology and probe geometry on nanoscale fatigue mechanisms of diamond-like carbon thin film

The aim of this paper is to investigate the mechanism of nanoscale fatigue using nano-impact and multiple-loading cycle nanoindentation tests, and compare it to previously reported findings of nanoscale fatigue using integrated stiffness and depth sensing approach. Two different film loading mechanism, loading history and indenter shapes are compared to comprehend the influence of test methodology on the nanoscale fatigue failure mechanisms of DLC film. An amorphous 100 nm thick DLC film was deposited on a 500 μm silicon substrate using sputtering of graphite target in pure argon atmosphere. Nano-impact and multiple-load cycle indentations were performed in the load range of 100 μN to 1000 μN and 0.1 mN to 100 mN, respectively. Both test types were conducted using conical and Berkovich indenters. Results indicate that for the case of conical indenter, the combination of nano-impact and multiple-loading cycle nanoindentation tests provide information on the life and failure mechanism of DLC film, which is comparable to the previously reported findings using the integrated stiffness and depth sensing approach. However, the comparison of results is sensitive to the applied load, loading mechanism, test-type and probe geometry. The loading mechanism and load history is therefore critical which also leads to two different definitions of film failure. The choice of exact test methodology, load and probe geometry should therefore be dictated by the in-service tribological conditions, and where necessary both test methodologies can be used to provide better insights of failure mechanism. Molecular dynamics (MD) simulations of the elastic response of nanoindentation is reported, which indicates that the elastic modulus of the film measured using MD simulation was higher than that experimentally measured. This difference is attributed to the factors related to the presence of material defects, crystal structure, residual stress, indenter geometry and loading/unloading rate differences between the MD and experimental results.

Influence of microstructure on the cutting behaviour of silicon

We use molecular dynamics simulation to study the mechanisms of plasticity during cutting of monocrystalline and polycrystalline silicon. Three scenarios are considered: (i) cutting a single crystal silicon workpiece with a single crystal diamond tool, (ii) cutting a polysilicon workpiece with a single crystal diamond tool, and (iii) cutting a single crystal silicon workpiece with a polycrystalline diamond tool. A long-range analytical bond order potential is used in the simulations, providing a more accurate picture of the atomic-scale mechanisms of brittle fracture, ductile plasticity, and structural changes in silicon. The MD simulation results show a unique phenomenon of brittle cracking typically inclined at an angle of 45° to 55° to the cut surface, leading to the formation of periodic arrays of nanogrooves in monocrystalline silicon, which is a new insight into previously published results. Furthermore, during cutting, silicon is found to undergo solid-state directional amorphisation without prior Si-I to Si-II (beta tin) transformation, which is in direct contrast to many previously published MD studies on this topic. Our simulations also predict that the propensity for amorphisation is significantly higher in single crystal silicon than in polysilicon, signifying that grain boundaries eases the material removal process.

Topical Review: The current understanding on the diamond machining of silicon carbide

The Glenn Research Centre of NASA, USA (www.grc.nasa.gov/WWW/SiC/, silicon carbide electronics) is in pursuit of realizing bulk manufacturing of silicon carbide (SiC), specifically by mechanical means. Single point diamond turning (SPDT) technology which employs diamond (the hardest naturally-occurring material realized to date) as a cutting tool to cut a workpiece is a highly productive manufacturing process. However, machining SiC using SPDT is a complex process and, while several experimental and analytical studies presented to date aid in the understanding of several critical processes of machining SiC, the current knowledge on the ductile behaviour of SiC is still sparse. This is due to a number of simultaneously occurring physical phenomena that may take place on multiple length and time scales. For example, nucleation of dislocation can take place at small inclusions that are of a few atoms in size and once nucleated, the interaction of these nucleations can manifest stresses on the micrometre length scales. The understanding of how stresses manifest during fracture in the brittle range, or dislocations/phase transformations in the ductile range, is crucial in understanding the brittle–ductile transition in SiC. Furthermore, there is a need to incorporate an appropriate simulation-based approach in the manufacturing research on SiC, owing primarily to the number of uncertainties in the experimental research that includes wear of the cutting tool, poor controllability of the nano-regime machining scale (effective thickness of cut), and coolant effects (interfacial phenomena between the tool, workpiece/chip and coolant), etc. In this review, these two problems are combined together to posit an improved understanding on the current theoretical knowledge on the SPDT of SiC obtained from molecular dynamics simulation.

Can a carbon nano-coating resist metallic phase transformation in silicon substrate during nanoimpact?

Nanomechanical response of a silicon specimen coated with a sp3 crystalline carbon coating (1.8 nm thickness) was investigated using MD simulation. A sharp conical rigid tip was impacted at the speed of 50 m/sec up to a depth of ~80% of the coating thickness. Unlike pure silicon specimen, no metallic phase transformation was observed i.e. a thin coating was able to resist Si-I to Si-II metallic phase transformation signifying that the coating could alter the stress distribution and thereby the contact tribology of the substrate. The stress state of the system, radial distribution function and the load-displacement curve were all aligned with above observations

Twinning anisotropy of tantalum during nanoindentation

Unlike other BCC metals, the plastic deformation of nanocrystalline Tantalum during compression is regulated by deformation twinning. Whether or not this twinning exhibits anisotropy was investigated through simulation of displacement-controlled nanoindentation test using molecular dynamics simulation. MD data was found to correlate well with the experimental data in terms of surface topography and hardness measurements. The mechanism of the transport of material was identified due to the formation and motion of prismatic dislocations loops (edge dislocations) belonging to the 1/2<111> type and <100> type Burgers vector family. Further analysis of crystal defects using a fully automated dislocation extraction algorithm (DXA) illuminated formation and migration of twin boundaries on the (110) and (111) orientation but not on the (010) orientation and most importantly after retraction all the dislocations disappeared on the (110) orientation suggesting twinning to dominate dislocation nucleation in driving plasticity in tantalum. A significant finding was that the maximum shear stress (critical Tresca stress) in the deformation zone exceeded the theoretical shear strength of tantalum (Shear modulus/ 2π~10.03 GPa) on the (010) orientation but was lower than it on the (110) and the (111) orientations. In light to this, the conventional lore of assuming the maximum shear stress being 0.465 times the mean contact pressure was found to break down at atomic scale.

A computational study of ionic liquids used in the production of fuels and biofuels

For the past decades it has been a worldwide concern to reduce the emission of harmful gases released during the combustion of fossil fuels. This goal has been addressed through the reduction of sulfur-containing compounds, and the replacement of fossil fuels by biofuels, such as bioethanol, produced in large scale from biomass. For this purpose, a new class of solvents, the Ionic Liquids (ILs), has been applied, aiming at developing new processes and replacing common organic solvents in the current processes. ILs can be composed by a large number of different combinations of cations and anions, which confer unique but desired properties to ILs. The ability of fine-tuning the properties of ILs to meet the requirements of a specific application range by mixing different cations and anions arises as the most relevant aspect for rendering ILs so attractive to researchers. Nonetheless, due to the huge number of possible combinations between the ions it is required the use of cheap predictive approaches for anticipating how they will act in a given situation. Molecular dynamics (MD) simulation is a statistical mechanics computational approach, based on Newton’s equations of motion, which can be used to study macroscopic systems at the atomic level, through the prediction of their properties, and other structural information. In the case of ILs, MD simulations have been extensively applied. The slow dynamics associated to ILs constitutes a challenge for their correct description that requires improvements and developments of existent force fields, as well as larger computational efforts (longer times of simulation). The present document reports studies based on MD simulations devoted to disclose the mechanisms of interaction established by ILs in systems representative of fuel and biofuels streams, and at biomass pre-treatment process. Hence, MD simulations were used to evaluate different systems composed of ILs and thiophene, benzene, water, ethanol and also glucose molecules. For the latter molecules, it was carried out a study aiming to ascertain the performance of a recently proposed force field (GROMOS 56ACARBO) to reproduce the dynamic behavior of such molecules in aqueous solution. The results here reported reveal that the interactions established by ILs are dependent on the individual characteristics of each IL. Generally, the polar character of ILs is deterministic in their propensity to interact with the other molecules. Although it is unquestionable the advantage of using MD simulations, it is necessary to recognize the need for improvements and developments of force fields, not only for a successful description of ILs, but also for other relevant compounds such as the carbohydrates.

High-affinity small molecule-phospholipid complex formation: Binding of siramesine to phosphatidic acid

Siramesine (SRM) is a sigma-2 receptor agonist which has been recently shown to inhibit growth of cancer cells. Fluorescence spectroscopy experiments revealed two distinct binding sites for this drug in phospholipid membranes. More specifically, acidic phospholipids retain siramesine on the bilayer surface due to a high-affinity interaction, reaching saturation at an apparent 1:1 drug-acidic phospholipid stoichiometry, where after the drug penetrates into the hydrocarbon core of the membrane. This behavior was confirmed using Langmuir films. Of the anionic phospholipids, the highest affinity, comparable to the affinities for the binding of small molecule ligands to proteins, was measured for phosphatidic acid (PA, mole fraction Of X-PA = 0.2 in phosphatidylcholine vesicles), yielding a molecular partition coefficient of 240 +/- 80 x 10(6). An MD simulation on the siramesine:PA interaction was in agreement with the above data. Taking into account the key role of PA as a signaling molecule promoting cell growth our results suggest a new paradigm for the development of anticancer drugs, viz. design of small molecules specifically scavenging phospholipids involved in the signaling cascades controlling cell behavior.

Effects of quenching rate on amorphous structures of Cu46Zr54 metallic glass

Metallic glass shows some superior properties different from crystalline, but the nature of amorphous structure and structural change during glass transition have not been completely understood yet. Molecular dynamics simulation provides intuitive insight into the microstructure and properties at atomistic level. Before probing into the microstructures of metallic glass with molecular dynamics (MD) simulation, it is important to obtain amorphous state first. In the current work, we reproduce the process of manufacturing metallic glass in laboratory including the melting, equilibrating and quenching procedure with molecular dynamics simulations. The structure changing at melting point and glass transition temperature are investigated with the different cooling processing. The partial radial distribution function (PRDF) is applied as a criterion to judge the final amorphous state obtained considering the quenching at different cooling rates and the effects of cooling rate on the formation of amorphous structures are further discussed.

Effects of quenching rate on amorphous structures of Cu46Zr54 metallic glass

Metallic glass shows some superior properties different from crystalline, but the nature of amorphous structure and structural change during glass transition have not been completely understood yet. Molecular dynamics simulation provides intuitive insight into the microstructure and properties at atomistic level. Before probing into the microstructures of metallic glass with molecular dynamics (MD) simulation, it is important to obtain amorphous state first. In the current work, we reproduce the process of manufacturing metallic glass in laboratory including the melting, equilibrating and quenching procedure with molecular dynamics simulations. The structure changing at melting point and glass transition temperature are investigated with the different cooling processing. The partial radial distribution function (PRDF) is applied as a criterion to judge the final amorphous state obtained considering the quenching at different cooling rates and the effects of cooling rate on the formation of amorphous structures are further discussed.

Characterization of structure and energy of TixCy cluster at early formation stage in iron matrix by molecular dynamics

The structure, energy and bonding property of TixCy clusters formed in iron matrix were studied through molecular dynamics (MD) simulation method. The initial clusters with 1D-linear, 2D-ring, and 3D-tetrahedral structures were determined and their stability was calculated. The effect of temperature on the stability of the clusters was also discussed. In addition, the dissociation path of TiC clusters in iron matrix and the corresponding energy variation were analyzed. © 2014 Elsevier B.V. All rights reserved.

How C-terminal carboxyamidation alters the biological activity of peptides from the venom of the eumenine solitary wasp

Inflammatory peptides display different types of post-transcriptional modifications, such as C-terminal amidation, that alter their biological activity. Here we describe the structural and molecular dynamics features of the mast cell degranulating peptide, eumenine mastoparan-AF (EMP-AF-NH2), found in the venom of the solitary wasp, and of its carboxyl-free C-terminal form (EMP-AF-COO-) characterized by a reduced activity. Circular dichroism indicates that both peptides switch from a random coil conformation in water to a helical structure in TFE and SDS micelles. NMR data, in 30% TFE, reveal that the two peptides fold into an alpha-helix spanning most of their length, while they differ in terms of molecular rigidity. To understand the origins of the conformational flexibility observed in the case of EMP-AF-COO-, a 5 ns MD simulation was carried out for each peptide, in an explicit water/TFE environment. The results show that the two peptides differ in an H-bond between Leu14 NH2 and the backbone carbonyl of Ile11. The loss of that H-bond in EMP-AF-COO- leads to a significant modification of its structural dynamics. In fact, as evidenced by essential dynamics analysis, while EMP-AF-NH2 exists mainly as a rigid structure, EMP-AF-COO- presents two helical stretches that fluctuate in some sort of independent fashion. We conclude that the diverse biological activity of the two peptides is not simply due to the reduction of the net positive charge, as generally suggested, but also to a structural perturbation of the amphipathic alpha-helix that affects their ability to perturb the cell membrane.

Why silver deposition is so fast: Solving the enigma of metal deposition

A perfect match: Silver deposition is one of the fastest electrochemical reactions, even though the Ag+ ion loses more than 5 eV solvation energy in the process. This phenomenon, an example of the enigma of metal deposition, was investigated by a combination of MD simulations, DFT, and specially developed theory. At the surface, the Ag+ ion experiences a strong interaction with the sp band of silver, which catalyzes the reaction. Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Interaction of cyclic and linear labaditin peptides with anionic and zwitterionic micelles

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)