Molecular dynamics simulation (MDS) to study nanoscale machining processes

Molecular Dynamics Simulations (MDS) are constantly being used to make important contributions to our fundamental understanding of material behaviour, at the atomic scale, for a variety of thermodynamic processes. This chapter shows that molecular dynamics simulation is a robust numerical analysis tool in addressing a range of complex nanofinishing (machining) problems that are otherwise difficult or impossible to understand using other methods. For example the mechanism of nanometric cutting of silicon carbide is influenced by a number of variables such as machine tool performance, machining conditions, material properties, and cutting tool performance (material microstructure and physical geometry of the contact) and all these variables cannot be monitored online through experimental examination. However, these could suitably be studied using an advanced simulation based approach such as MDS. This chapter details how MD simulation can be used as a research and commercial tool to understand key issues of ultra precision manufacturing research problems and a specific case was addressed by studying diamond machining of silicon carbide. While this is appreciable, there are a lot of challenges and opportunities in this fertile area. For example, the world of MD simulations is dependent on present day computers and the accuracy and reliability of potential energy functions [109]. This presents a limitation: Real-world scale simulation models are yet to be developed. The simulated length and timescales are far shorter than the experimental ones which couples further with the fact that contact loading simulations are typically done in the speed range of a few hundreds of m/sec against the experimental speed of typically about 1 m/sec [17]. Consequently, MD simulations suffer from the spurious effects of high cutting speeds and the accuracy of the simulation results has yet to be fully explored. The development of user-friendly software could help facilitate molecular dynamics as an integral part of computer-aided design and manufacturing to tackle a range of machining problems from all perspectives, including materials science (phase of the material formed due to the sub-surface deformation layer), electronics and optics (properties of the finished machined surface due to the metallurgical transformation in comparison to the bulk material), and mechanical engineering (extent of residual stresses in the machined component) [110]. Overall, this chapter provided key information concerning diamond machining of SiC which is classed as hard, brittle material. From the analysis presented in the earlier sections, MD simulation has helped in understanding the effects of crystal anisotropy in nanometric cutting of 3C-SiC by revealing the atomic-level deformation mechanisms for different crystal orientations and cutting directions. In addition to this, the MD simulation revealed that the material removal mechanism on the (111) surface of 3C-SiC (akin to diamond) is dominated by cleavage. These understandings led to the development of a new approach named the “surface defect machining” method which has the potential to be more effective to implement than ductile mode micro laser assisted machining or conventional nanometric cutting.

Surface energy and surface proton order of the ice Ih basal and prism surfaces

Density-functional theory (DFT) is used to examine the basal and prism surfaces of ice Ih. Similar surface energies are obtained for the two surfaces; however, in each case a strong dependence of the surface energy on surface proton order is identified. This dependence, which can be as much as 50% of the absolute surface energy, is significantly larger than the bulk dependence (< 1%) on proton order, suggesting that the thermodynamic ground state of the ice surface will remain proton ordered well above the bulk order-disorder temperature of about 72 K. On the basal surface this suggestion is supported by Monte Carlo simulations with an empirical potential and solution of a 2D Ising model with nearest neighbor interactions taken from DFT. Order parameters that define the surface energy of each surface in terms of nearest neighbor interactions between dangling OH bonds (those which point out of the surface into vacuum) have been identified and are discussed. Overall, these results suggest that proton order-disorder effects have a profound impact on the stability of ice surfaces and will most likely have an effect on ice surface reactivity as well as ice crystal growth and morphology. S Supplementary data are available from stacks.iop.org/JPhysCM/22/074209/mmedia

SrTiO3(001)(2x1) reconstructions: First-principles calculations of surface energy and atomic structure compared with scanning tunneling microscopy images

(1x1) and (2x1) reconstructions of the (001) SrTiO3 surface were studied using the first-principles full-potential linear muffin-tin orbital method. Surface energies were calculated as a function of TiO2 chemical potential, oxygen partial pressure and temperature. The (1x1) unreconstructed surfaces were found to be energetically stable for many of the conditions considered. Under conditions of very low oxygen partial pressure the (2x1) Ti2O3 reconstruction [Martin R. Castell, Surf. Sci. 505, 1 (2002)] is stable. The question as to why STM images of the (1x1) surfaces have not been obtained was addressed by calculating charge densities for each surface. These suggest that the (2x1) reconstructions would be easier to image than the (1x1) surfaces. The possibility that the presence of oxygen vacancies would destabilise the (1x1) surfaces was also investigated. If the (1x1) surfaces are unstable then there exists the further possibility that the (2x1) DL-TiO2 reconstruction [Natasha Erdman Nature (London) 419, 55 (2002)] is stable in a TiO2-rich environment and for p(O2)>10(-18) atm.

Application of the zero-range potential model to positron annihilation on molecules

In this paper we use a zero-range potential (ZRP) method to model positron interaction with molecules. This allows us to investigate the e?ect of molecular vibrations on positron–molecule annihilation using the van der Waals dimer Kr2 as an example. We also use the ZRP to explore positron binding to polyatomics and examine the dependence of the binding energy on the size of the molecule for alkanes. We ?nd that a second bound state appears for a molecule with ten carbons, similar to recent experimental evidence for such a state emerging in alkanes with twelve carbons.

Exchange-correlation energy and potential as approximate functionals of occupied and virtual Kohn-Sham orbitals: Application to dissociating H-2

The standard local density approximation and generalized gradient approximations fail to properly describe the dissociation of an electron pair bond, yielding large errors (on the order of 50 kcal/mol) at long bond distances. To remedy this failure, a self-consistent Kohn-Sham (KS) method is proposed with the exchange-correlation (xc) energy and potential depending on both occupied and virtual KS orbitals. The xc energy functional of Buijse and Baerends [Mol. Phys. 100, 401 (2002); Phys. Rev. Lett. 87, 133004 (2001)] is employed, which, based on an ansatz for the xc-hole amplitude, is able to reproduce the important dynamical and nondynamical effects of Coulomb correlation through the efficient use of virtual orbitals. Self-consistent calculations require the corresponding xc potential to be obtained, to which end the optimized effective potential (OEP) method is used within the common energy denominator approximation for the static orbital Green's function. The problem of the asymptotic divergence of the xc potential of the OEP when a finite number of virtual orbitals is used is addressed. The self-consistent calculations reproduce very well the entire H-2 potential curve, describing correctly the gradual buildup of strong left-right correlation in stretched H-2. (C) 2003 American Institute of Physics.

Exchange potential from the common energy denominator approximation for the Kohn-Sham Green's function: Application to (hyper)polarizabilities of molecular chains

An approximate Kohn-Sham (KS) exchange potential v(xsigma)(CEDA) is developed, based on the common energy denominator approximation (CEDA) for the static orbital Green's function, which preserves the essential structure of the density response function. v(xsigma)(CEDA) is an explicit functional of the occupied KS orbitals, which has the Slater v(Ssigma) and response v(respsigma)(CEDA) potentials as its components. The latter exhibits the characteristic step structure with "diagonal" contributions from the orbital densities \psi(isigma)\(2), as well as "off-diagonal" ones from the occupied-occupied orbital products psi(isigma)psi(j(not equal1)sigma). Comparison of the results of atomic and molecular ground-state CEDA calculations with those of the Krieger-Li-Iafrate (KLI), exact exchange (EXX), and Hartree-Fock (HF) methods show, that both KLI and CEDA potentials can be considered as very good analytical "closure approximations" to the exact KS exchange potential. The total CEDA and KLI energies nearly coincide with the EXX ones and the corresponding orbital energies epsilon(isigma) are rather close to each other for the light atoms and small molecules considered. The CEDA, KLI, EXX-epsilon(isigma) values provide the qualitatively correct order of ionizations and they give an estimate of VIPs comparable to that of the HF Koopmans' theorem. However, the additional off-diagonal orbital structure of v(xsigma)(CEDA) appears to be essential for the calculated response properties of molecular chains. KLI already considerably improves the calculated (hyper)polarizabilities of the prototype hydrogen chains H-n over local density approximation (LDA) and standard generalized gradient approximations (GGAs), while the CEDA results are definitely an improvement over the KLI ones. The reasons of this success are the specific orbital structures of the CEDA and KLI response potentials, which produce in an external field an ultranonlocal field-counteracting exchange potential. (C) 2002 American Institute of Physics.

Electrochemical versus gas-phase oxidation of ru single-crystal surfaces

The electrochemical uptake of oxygen on a Ru(0001) electrode was investigated by electron diffraction, Auger spectroscopy, and cyclic voltammetry. An ordered (2 × 2)-O overlayer forms at a potential close to the hydrogen region. At +0.42 and +1.12 V vs Ag/AgCl, a (3 × 1) phase and a (1 × 1)-O phase, respectively, emerge. When the Ru electrode potential is maintained at +1.12 V for 2 min, RuO2 grows epitaxially with its (100) plane parallel to the Ru(0001) surface. In contrast to the RuO domains, the non-oxidized regions of the Ru electrode surface are flat. If, however, the electrode potential is increased to +1.98 V for 2 min, the remaining non-oxidized Ru area also becomes rough. These findings are compared with O overlayers and oxides on the Ru(0001) and Ru(101¯1) surfaces created by exposure to gaseous O under UHV conditions. On the other hand, gas-phase oxidation of the Ru(101¯0) surface leads to the formation of RuO with a (100) orientation. It is concluded that the difference in surface energy between RuO(110) and RuO(100) is quite small. RuO again grows epitaxially on Ru(0001), but with the (110) face oriented parallel to the Ru(0001) surface. The electrochemical oxidation of the Ru(0001) electrode surface proceeds via a 3-dimensional growth mechanism with a mean cluster size of 1.6 nm, whereas under UHV conditions, a 2-dimensional oxide film (1-2 nm thick) is epitaxially formed with an average domain size of 20 µm. © 2000 American Chemical Society.

On the potential of significance-driven execution for energy-aware HPC

Dynamic Voltage and Frequency Scaling (DVFS) exhibits fundamental limitations as a method to reduce energy consumption in computing systems. In the HPC domain, where performance is of highest priority and codes are heavily optimized to minimize idle time, DVFS has limited opportunity to achieve substantial energy savings. This paper explores if operating processors Near the transistor Threshold Volt- age (NTV) is a better alternative to DVFS for break- ing the power wall in HPC. NTV presents challenges, since it compromises both performance and reliability to reduce power consumption. We present a first of its kind study of a significance-driven execution paradigm that selectively uses NTV and algorithmic error tolerance to reduce energy consumption in performance- constrained HPC environments. Using an iterative algorithm as a use case, we present an adaptive execution scheme that switches between near-threshold execution on many cores and above-threshold execution on one core, as the computational significance of iterations in the algorithm evolves over time. Using this scheme on state-of-the-art hardware, we demonstrate energy savings ranging between 35% to 67%, while compromising neither correctness nor performance.

Oxygen vacancy formation in CeO2 and Ce1-xZrxO2 solid solutions:electron localization, electrostatic potential and structural relaxation

Ceria (CeO2) and ceria-based composite materials, especially Ce1-xZrxO2 solid solutions, possess a wide range of applications in many important catalytic processes, such as three-way catalysts, owing to their excellent oxygen storage capacity (OSC) through the oxygen vacancy formation and refilling. Much of this activity has focused on the understanding of the electronic and structural properties of defective CeO2 with and without doping, and comprehending the determining factor for oxygen vacancy formation and the rule to tune the formation energy by doping has constituted a central issue in material chemistry related to ceria. However, the calculation on electronic structures and the corresponding relaxation patterns in defective CeO2-x oxides remains at present a challenge in the DFT framework. A pragmatic approach based on density functional theory with the inclusion of on-site Coulomb correction, i.e. the so-called DFT + U technique, has been extensively applied in the majority of recent theoretical investigations. Firstly, we review briefly the latest electronic structure calculations of defective CeO2(111), focusing on the phenomenon of multiple configurations of the localized 4f electrons, as well as the discussions of its formation mechanism and the catalytic role in activating the O-2 molecule. Secondly, aiming at shedding light on the doping effect on tuning the oxygen vacancy formation in ceria-based solid solutions, we summarize the recent theoretical results of Ce1-xZrxO2 solid solutions in terms of the effect of dopant concentrations and crystal phases. A general model on O vacancy formation is also discussed; it consists of electrostatic and structural relaxation terms, and the vital role of the later is emphasized. Particularly, we discuss the crucial role of the localized structural relaxation patterns in determining the superb oxygen storage capacity in kappa-phase Ce1-xZr1-xO2. Thirdly, we briefly discuss some interesting findings for the oxygen vacancy formation in pure ceria nanoparticles (NPs) uncovered by DFT calculations and compare those with the bulk or extended surfaces of ceria as well as different particle sizes, emphasizing the role of the electrostatic field in determining the O vacancy formation.

High-energy monoenergetic proton beams from two stage acceleration with a slow laser pulse

We present a new regime to generate high-energy quasimonoenergetic proton beams in a "slow-pulse" regime, where the laser group velocity vg<c is reduced by an extended near-critical density plasma. In this regime, for properly matched laser intensity and group velocity, ions initially accelerated by the light sail (LS) mode can be further trapped and reflected by the snowplough potential generated by the laser in the near-critical density plasma. These two acceleration stages are connected by the onset of Rayleigh-Taylor-like (RT) instability. The usual ion energy spectrum broadening by RT instability is controlled and high quality proton beams can be generated. It is shown by multidimensional particle-in-cell simulation that quasimonoenergetic proton beams with energy up to hundreds of MeV can be generated at laser intensities of 1021W/cm².