Anisotropy of single crystal 3C-SiC during nanometric cutting

3C–SiC (the only polytype of SiC that resides in a diamond cubic lattice structure) is a relatively new material that exhibits most of the desirable engineering properties required for advanced electronic applications. The anisotropy exhibited by 3C–SiC during its nanometric cutting is significant, and the potential for its exploitation has yet to be fully investigated. This paper aims to understand the influence of crystal anisotropy of 3C–SiC on its cutting behaviour. A molecular dynamics simulation model was developed to simulate the nanometric cutting of single-crystal 3C–SiC in nine (9) distinct combinations of crystal orientations and cutting directions, i.e. (1?1?1) <-1?1?0>, (1?1?1) <-2?1?1>, (1?1?0) <-1?1?0>, (1?1?0) <0?0?1>, (1?1?0) <1?1?-2>, (0?0?1) <-1?1?0>, (0?0?1) <1?0?0>, (1?1?-2) <1?-1?0> and (1?-2?0) <2?1?0>.

In order to ensure the reliability of the simulation results, two separate simulation trials were carried out with different machining parameters. In the first trial, a cutting tool rake angle of -25°, d/r (uncut chip thickness/cutting edge radius) ratio of 0.57 and cutting velocity of 10 m s-1 were used whereas a second trial was done using a cutting tool rake angle of -30°, d/r ratio of 1 and cutting velocity of 4 m s-1. Both the trials showed similar anisotropic variation.

The simulated orthogonal components of thrust force in 3C–SiC showed a variation of up to 45%, while the resultant cutting forces showed a variation of 37%. This suggests that 3C–SiC is highly anisotropic in its ease of deformation. These results corroborate with the experimentally observed anisotropic variation of 43.6% in Young's modulus of 3C–SiC. The recently developed dislocation extraction algorithm (DXA) [1, 2] was employed to detect the nucleation of dislocations in the MD simulations of varying cutting orientations and cutting directions. Based on the overall analysis, it was found that 3C–SiC offers ease of deformation on either (1?1?1) <-1?1?0>, (1?1?0) <0?0?1>, or (1?0?0) <1?0?0> setups.

Shear instability of nanocrystalline silicon carbide during nanometric cutting

The shear instability of the nanoscrystalline 3C-SiC during nanometric cutting at a cutting speed of 100?m/s has been investigated using molecular dynamics simulation. The deviatoric stress in the cutting zone was found to cause sp3-sp2 disorder resulting in the local formation of SiC-graphene and Herzfeld-Mott transitions of 3C-SiC at much lower transition pressures than that required under pure compression. Besides explaining the ductility of SiC at 1500?K, this is a promising phenomenon in general nanoscale engineering of SiC. It shows that modifying the tetrahedral bonding of 3C-SiC, which would otherwise require sophisticated pressure cells, can be achieved more easily by introducing non-hydrostatic stress conditions.

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.

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.

Universal tight binding model for chemical reactions in solution and at surfaces. III. Stoichiometric and reduced surfaces of titania and the adsorption of water

We demonstrate a model for stoichiometric and reduced titanium dioxide intended for use in molecular dynamics and other atomistic simulations and based in the polarizable ion tight binding theory. This extends the model introduced in two previous papers from molecular and liquid applications into the solid state, thus completing the task of providing a comprehensive and unified scheme for studying chemical reactions, particularly aimed at problems in catalysis and electrochemistry. As before, experimental results are given priority over theoretical ones in selecting targets for model fitting, for which we used crystal parameters and band gaps of titania bulk polymorphs, rutile and anatase. The model is applied to six low index titania surfaces, with and without oxygen vacancies and adsorbed water molecules, both in dissociated and non-dissociated states. Finally, we present the results of molecular dynamics simulation of an anatase cluster with a number of adsorbed water molecules and discuss the role of edge and corner atoms of the cluster. (C) 2014 AIP Publishing LLC.

Cyclic nanoindentation and nano-impact fatigue characterisation of functionally graded TiN/TiNi shape memory film

This paper investigates the mechanism of nanoscale fatigue of functionally graded TiN/TiNi films using nano-impact and multiple-loading-cycle nanoindentation tests. The functionally graded films were deposited on silicon substrate, in which TiNi films maintain shape memory and pseudo elastic behavior, while a modified TiN surface layer provides tribological and anti-corrosion properties. Nanomechanical tests were performed to comprehend the localized film performance and failure modes of the functionally graded film using NanoTestTM equipped with Berkovich and conical indenter between 100 μN to 500 mN loads. The loading mechanism and load history are critical to define film failure modes (i.e. backward depth deviation) including the shape memory effect of the functionally graded layer. The results are sensitive to the applied load, loading type (e.g. semi-static, dynamic) and probe geometry. Based on indentation force-depth profiles, depth-time data and post-test surface observations of films, it is concluded that the shape of the nanoindenter is critical in inducing the localized indentation stress and film failure, including shape recovery at the lower load range. Elastic-plastic finite element (FE) simulation during nanoindentation loading indicated that the location of subsurface maximum stress near the interface influences the backward depth deviation type of film failure. A standalone, molecular dynamics simulation was performed with the help of a long range potential energy function to simulate the tensile test of TiN nanowire with two different aspect ratios to investigate the theory of its failure mechanism.

Cloud benchmarking for performance

How can applications be deployed on the cloud to achieve maximum performance? This question has become significant and challenging with the availability of a wide variety of Virtual Machines (VMs) with different performance capabilities in the cloud. The above question is addressed by proposing a six step benchmarking methodology in which a user provides a set of four weights that indicate how important each of the following groups: memory, processor, computation and storage are to the application that needs to be executed on the cloud. The weights along with cloud benchmarking data are used to generate a ranking of VMs that can maximise performance of the application. The rankings are validated through an empirical analysis using two case study applications, the first is a financial risk application and the second is a molecular dynamics simulation, which are both representative of workloads that can benefit from execution on the cloud. Both case studies validate the feasibility of the methodology and highlight that maximum performance can be achieved on the cloud by selecting the top ranked VMs produced by the methodology.

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.

Differential peptide binding to CD40 evokes counteractive responses.

The antigen-presenting cell-expressed CD40 is implied in the regulation of counteractive immune responses such as induction of pro-inflammatory and anti-inflammatory cytokines interleukin (IL)-12 and IL-10, respectively. The mechanism of this duality in CD40 function remains unknown. Here, we investigated whether such duality depends on ligand binding. Based on CD40 binding, we identifed two dodecameric peptides, peptide-7 and peptide-19, from the phage peptide library. Peptide-7 induces IL-10 and increases Leishmania donovani infection in macrophages, whereas peptide-19 induces IL-12 and reduces L. donovani infection. CD40-peptide interaction analyses by surface plasmon resonance and atomic force microscopy suggest that the functional differences are not associated with the studied interaction parameters. The molecular dynamic simulation of the CD40-peptides interaction suggests that these two peptides bind to two different places on CD40. Thus, we suggest for the first time that differential binding of the ligands imparts functional duality to CD40.

Collective dynamic properties in simple crystals: Influence of the structural disorder

Collective dynamic properties in Lennard-Jones crystals are investigated by molecular dynamics simulation. The study is focused on properties such as the dynamic structure factors, the longitudinal and transverse currents and the density of states. The influence on these properties of the structural disorder is analyzed by comparing the results for one-component crystals with those for liquids and supercooled liquids at analogous conditions. The effects of species-disorder on the collective properties of binary crystals are also discussed.

Theory of laser-induced ultrafast structural changes in solids

The present thesis is a contribution to the study of laser-solid interaction. Despite the numerous applications resulting from the recent use of laser technology, there is still a lack of satisfactory answers to theoretical questions regarding the mechanism leading to the structural changes induced by femtosecond lasers in materials. We provide here theoretical approaches for the description of the structural response of different solids (cerium, samarium sulfide, bismuth and germanium) to femtosecond laser excitation. Particular interest is given to the description of the effects of the laser pulse on the electronic systems and changes of the potential energy surface for the ions. Although the general approach of laser-excited solids remains the same, the potential energy surface which drives the structural changes is calculated with different theoretical models for each material. This is due to the difference of the electronic properties of the studied systems. We use the Falicov model combined with an hydrodynamic method to study photoinduced phase changes in cerium. The local density approximation (LDA) together with the Hubbard-type Hamiltonian (LDA+U) in the framework of density functional theory (DFT) is used to describe the structural properties of samarium sulfide. We parametrize the time-dependent potential energy surface (calculated using DFT+ LDA) of bismuth on which we perform quantum dynamical simulations to study the experimentally observed amplitude collapse and revival of coherent $A_{1g}$ phonons. On the basis of a time-dependent potential energy surface calculated from a non-orthogonal tight binding Hamiltonian, we perform molecular dynamics simulation to analyze the time evolution (coherent phonons, ultrafast nonthermal melting) of germanium under laser excitation. The thermodynamic equilibrium properties of germanium are also reported. With the obtained results we are able to give many clarifications and interpretations of experimental results and also make predictions.

Structure and magnetic properties of mono- and bi-dispersed ferrofluids as revealed by simulations

Using a Langevin molecular dynamics simulation, we show that the magnetic properties of a mono- and bi-dispersed ferrofluid system depend on the volume fraction and the dipolar coupling parameter. For the bi-dispersed system, most of the chains are formed by the large particles, but the aggregation behavior of the large particles is hindered by the presence of the small particles, which are predominantly attached to the end of the particle chain. To further elucidate the microscopic fluid structure, anisotropic structure factors are calculated.

The effect of electric field direction on the-shear stress of ER fluids

The effect of the direction of external electric field on the shear stress of an ER fluid has been studied by molecular-dynamics simulation. Due to the formation of inclined chains, the shear stress strongly depends on the direction of the field, and it may be very large under some special field direction. And theoretical model of ideal microstructure of ER fluids has proved this result. Thus the ER effect may be greatly enhanced just by choosing an optimum direction for the field without any additional requirement, suggesting a promising way to the practical application of ER fluids.

Spin current in the Mobius cyclacene belts

We report on density functional theory studies of the electronic structure and magnetic properties of Mobius-[n]cyclacenes. The geometry of Mobius bands presents a modulation of bond lengths that is needed to accommodate the twist. This modulation takes the form of bond alternation defects analogous to those of solitons in polyacetylene. The ground state of all Mobius bands is a triplet, with a spin density distribution that follows the bond length modulation. A molecular dynamics simulation of the Mobius cyclacene at 300 K shows that the twist travels around the belt inducing a magnetic current. (C) 2009 Elsevier B. V. All rights reserved.