Quantum corrected Langevin dynamics for adsorbates on metal surfaces interacting with hot electrons

We investigate the importance of including quantized initial conditions in Langevin dynamics for adsorbates interacting with a thermal reservoir of electrons. For quadratic potentials the time evolution is exactly described by a classical Langevin equation and it is shown how to rigorously obtain quantum mechanical probabilities from the classical phase space distributions resulting from the dynamics. At short time scales, classical and quasiclassical initial conditions lead to wrong results and only correctly quantized initial conditions give a close agreement with an inherently quantum mechanical master equation approach. With CO on Cu(100) as an example, we demonstrate the effect for a system with ab initio frictional tensor and potential energy surfaces and show that quantizing the initial conditions can have a large impact on both the desorption probability and the distribution of molecular vibrational states

Dislocation nucleation and vacancy formation during high-speed deformation of fcc metals

Recently, a dislocation free deformation mechanism was proposed by Kiritani et al., based on a series of experiments where thin foils of fcc metals were deformed at very high strain rates. In the experimental study, they observed a large density of stacking fault tetrahedra, but very low dislocation densities in the foils after deformation. This was interpreted as evidence for a new dislocation-free deformation mechanism, resulting in a very high vacancy production rate. In this paper we investigate this proposition using large-scale computer simulations of bulk and thin films of copper. To favour such a dislocation-free deformation mechanism, we have made dislocation nucleation very difficult by not introducing any potential dislocation sources in the initial configuration. Nevertheless, we observe the nucleation of dislocation loops, and the deformation is carried by dislocations. The dislocations are nucleated as single Shockley partials. The large stresses required before dislocations are nucleated result in a very high dislocation density, and therefore in many inelastic interactions between the dislocations. These interactions create vacancies, and a very large vacancy concentration is quickly reached.Comment: LaTeX2e, 8 pages, PostScript figures included. Minor modifications only. Final version, to appear in Philos. Mag. Let

Modelling of dislocation generation and interaction during high-speed deformation of metals

Recent experiments by Kiritani et al. have revealed a surprisingly high rate of vacancy production during high-speed deformation of thin foils of fcc metals. Virtually no dislocations are seen after the deformation. This is interpreted as evidence for a dislocation-free deformation mechanism at very high strain rates. We have used molecular-dynamics simulations to investigate high-speed deformation of copper crystals. Even though no pre-existing dislocation sources are present in the initial system, dislocations are quickly nucleated and a very high dislocation density is reached during the deformation. Due to the high density of dislocations, many inelastic interactions occur between dislocations, resulting in the generation of vacancies. After the deformation, a very high density of vacancies is observed, in agreement with the experimental observations. The processes responsible for the generation of vacancies are investigated. The main process is found to be incomplete annihilation of segments of edge dislocations on adjacent slip planes. The dislocations are also seen to be participating in complicated dislocation reactions, where sessile dislocation segments are constantly formed and destroyed.Comment: 8 pages, LaTeX2e + PS figures. Presented at the Third Workshop on High-speed Plastic Deformation, Hiroshima, August 200

Calculation of Elastic Green's Functions for Lattices with Cavities

In this Brief Report, we present an algorithm for calculating the elastic Lattice Greens Function of a regular lattice, in which defects are created by removing lattice points. The method is computationally efficient, since the required matrix operations are on matrices that scale with the size of the defect subspace, and not with the size of the full lattice. This method allows the treatment of force fields with multi-atom interactions.Comment: 3 pages. RevTeX, using epsfig.sty. One figur

Simulations of super-structure domain walls in two dimensional assemblies of magnetic nanoparticles

We simulate the formation of domain walls in two-dimensional assemblies of magnetic nanoparticles. Particle parameters are chosen to match recent electron holography and Lorentz microscopy studies of almost monodisperse cobalt nanoparticles assembled into regular, elongated lattices. As the particles are small enough to consist of a single magnetic domain each, their magnetic interactions can be described by a spin model in which each particle is assigned a macroscopic "superspin." Thus, the magnetic behaviour of these lattices may be compared to magnetic crystals with nanoparticle superspins taking the role of the atomic spins. The coupling is, however, different. The superspins interact only by dipolar interactions as exchange coupling between individual nanoparticles may be neglected due to interparticle spacing. We observe that it is energetically favorable to introduce domain walls oriented along the long dimension of nanoparticle assemblies rather than along the short dimension. This is unlike what is typically observed in continuous magnetic materials, where the exchange interaction introduces an energetic cost proportional to the area of the domain walls. Structural disorder, which will always be present in realistic assemblies, pins longitudinal domain walls when the external field is reversed, and makes a gradual reversal of the magnetization by migration of longitudinal domain walls possible, in agreement with previous experimental results. (C) 2015 AIP Publishing LLC

Effects of crack tip geometry on dislocation emission and cleavage: A possible path to enhanced ductility

We present a systematic study of the effect of crack blunting on subsequent crack propagation and dislocation emission. We show that the stress intensity factor required to propagate the crack is increased as the crack is blunted by up to thirteen atomic layers, but only by a relatively modest amount for a crack with a sharp 60$^\circ$ corner. The effect of the blunting is far less than would be expected from a smoothly blunted crack; the sharp corners preserve the stress concentration, reducing the effect of the blunting. However, for some material parameters blunting changes the preferred deformation mode from brittle cleavage to dislocation emission. In such materials, the absorption of preexisting dislocations by the crack tip can cause the crack tip to be locally arrested, causing a significant increase in the microscopic toughness of the crack tip. Continuum plasticity models have shown that even a moderate increase in the microscopic toughness can lead to an increase in the macroscopic fracture toughness of the material by several orders of magnitude. We thus propose an atomic-scale mechanism at the crack tip, that ultimately may lead to a high fracture toughness in some materials where a sharp crack would seem to be able to propagate in a brittle manner. Results for blunt cracks loaded in mode II are also presented.Comment: 12 pages, REVTeX using epsfig.sty. 13 PostScript figures. Final version to appear in Phys. Rev. B. Main changes: Discussion slightly shortened, one figure remove

Matching conditions in the quasicontinuum method: Removal of the error introduced at the interface between the coarse-grained and fully atomistic region

金沢大学大学院自然科学研究科機能開発システム金沢大学工学部The quasicontinuum method is a way of reducing the number of degrees of freedom in an atomistic simulation by removing the majority of the atoms in regions of slowly varying strain fields. Due to the different ways the energy of the atoms is calculated in the coarse-grained regions and the regions where all the atoms are present, unphysical forces called "ghost forces" arise at the interfaces. Corrections may be used to almost remove the ghost forces, but the correction forces are nonconservative, ruining energy conservation in dynamic simulations. We show that it is possible to formulate the quasicontinuum method without these problems by introducing a buffer layer between the two regions of space. The method is applicable to short-ranged potentials in the face-centered cubic, body-centered cubic, and hexagonal close-packed crystal structures