Search

We introduce a method to visualize dislocations along grain boundaries at the atomic level. It uses an atomic-level Nye tensor, representing the dislocation density. To calculate the Nye tensor at grain boundaries, we extend the Hartley-Mishin strain gradient calculation to the displacement shift complete lattice. We show that the method is effective in visualizing disconnections and the dislocation content of grain boundary phase junctions in body-centered cubic tungsten, as well as face-centered cubic copper. In addition, we use the method to characterize the morphology of a two-dimensional grain boundary phase nucleus in a symmetric tilt grain boundary in tungsten. This method can be applied to both bulk dislocations and grain boundary disconnections, which makes it ideal for studying the interactions and reactions of bulk dislocations with grain boundaries, and grain boundary disconnections.

We present the results of the first Charged-Particle Transport Coefficient Code Comparison Workshop, which was held in Albuquerque, NM October 4-6, 2016. In this first workshop, scientists from eight institutions and four countries gathered to compare calculations of transport coefficients including thermal and electrical conduction, electron-ion coupling, inter-ion diffusion, ion viscosity, and charged particle stopping powers. Here, we give general background on Coulomb coupling and computational expense, review where some transport coefficients appear in hydrodynamic equations, and present the submitted data. Large variations are found when either the relevant Coulomb coupling parameter is large or computational expense causes difficulties. Understanding the general accuracy and uncertainty associated with such transport coefficients is important for quantifying errors in hydrodynamic simulations of inertial confinement fusion and high-energy density experiments., Comment: 45 pages, 17 figures

Alloys expand or contract as concentrations change, and the resulting relationship between atomic volume and alloy content is an important property of the solid. While a well-known approximation posits that the atomic volume varies linearly with concentration (Zen`s law), the actual variation is more complicated. Here we use an apparent size of the solute (solvent) atom and the elasticity to derive explicit analytical expressions for the atomic volume of binary solid alloys. Two approximations, continuum and terminal, are proposed. Deviations from Zen`s law are studied for 22 binary alloy systems.

We report a computational discovery of novel grain boundary structures and multiple grain boundary phases in elemental bcc tungsten. While grain boundary structures created by the \gamma-surface method as a union of two perfect half crystals have been studied extensively, it is known that the method has limitations and does not always predict the correct ground states. Here, we use a newly developed computational tool, based on evolutionary algorithms, to perform a grand-canonical search of a high-angle symmetric tilt boundary in tungsten, and we find new ground states and multiple phases that cannot be described using the conventional structural unit model. We use MD simulations to demonstrate that the new structures can coexist at finite temperature in a closed system, confirming these are examples of different GB phases. The new ground state is confirmed by first-principles calculations., Comment: 19 pages, 6 figures

Quantization techniques have been applied in many challenging finance applications, including pricing claims with path dependence and early exercise features, stochastic optimal control, filtering problems and efficient calibration of large derivative books. Recursive Marginal Quantization of the Euler scheme has recently been proposed as an efficient numerical method for evaluating functionals of solutions of stochastic differential equations. This method involves recursively quantizing the conditional marginals of the discrete-time Euler approximation of the underlying process. By generalizing this approach, we show that it is possible to perform recursive marginal quantization for two higher-order schemes: the Milstein scheme and a simplified weak order 2.0 scheme. As part of this generalization a simple matrix formulation is presented, allowing efficient implementation. We further extend the applicability of recursive marginal quantization by showing how absorption and reflection at the zero boundary may be incorporated, when this is necessary. To illustrate the improved accuracy of the higher order schemes, various computations are performed using geometric Brownian motion and its generalization, the constant elasticity of variance model. For both processes, we show numerical evidence of improved weak order convergence and we compare the marginal distributions implied by the three schemes to the known analytical distributions. By pricing European, Bermudan and Barrier options, further evidence of improved accuracy of the higher order schemes is demonstrated.

We present a fluid dynamics video showing the results of a 9-billion atom molecular dynamics simulation of complex fluid flow in molten copper and aluminum. Starting with an atomically flat interface, a shear is imposed along the copper-aluminum interface and random atomic fluctuations seed the formation of vortices. These vortices grow due to the Kelvin-Helmholtz instability. The resulting vortical structures are beautifully intricate, decorated with secondary instabilities and complex mixing phenomena. This work performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344., Comment: Description of video submitted to APS DFD Gallery of Fluid Motion 2008

We present results from ab initio calculations of the mechanical properties of the rhombohedral phase (beta) of vanadium metal reported in recent experiments, and other predicted high-pressure phases (gamma and bcc), focusing on properties relevant to dynamic experiments. We find that the volume change associated with these transitions is small: no more than 0.15% (for beta - gamma). Calculations of the single crystal and polycrystal elastic moduli (stress-strain coefficients) reveal a remarkably small discontinuity in the shear modulus and other elastic properties across the phase transitions even at zero temperature where the transitions are first order., Comment: 6 pages, 3 figures

Dissipation in the flexural dynamics of doubly clamped nanomechanical bar resonators is investigated using molecular dynamics simulation. The dependence of the quality factor Q on temperature and the size of the resonator is calculated from direct simulation of the oscillation of a series of Si <001> bars with bare {100} dimerized surfaces. The bar widths range from 3.3 to 8.7nm, all with a fixed length of 22nm. The fundamental mode frequencies range from 40 to 90GHz and Q from 10^2 near 1000K to 10^4 near 50K. The quality factor is shown to be limited by defects in the reconstructed surface., Comment: 4 pages, 5 figures

Void coalescence and interaction in dynamic fracture of ductile metals have been investigated using three-dimensional strain-controlled multi-million atom molecular dynamics simulations of copper. The correlated growth of two voids during the coalescence process leading to fracture is investigated, both in terms of its onset and the ensuing dynamical interactions. Void interactions are quantified through the rate of reduction of the distance between the voids, through the correlated directional growth of the voids, and through correlated shape evolution of the voids. The critical inter-void ligament distance marking the onset of coalescence is shown to be approximately one void radius based on the quantification measurements used, independent of the initial separation distance between the voids and the strain-rate of the expansion of the system. The interaction of the voids is not reflected in the volumetric asymptotic growth rate of the voids, as demonstrated here. Finally, the practice of using a single void and periodic boundary conditions to study coalescence is examined critically and shown to produce results markedly different than the coalescence of a pair of isolated voids., Comment: Accepted for publication in Physical Review B

Molecular dynamics simulations in three-dimensional copper have been performed to quantify the void coalescence process leading to fracture. The correlated growth of the voids during their linking is investigated both in terms of the onset of coalescence and the ensuing dynamical interactions through the rate of reduction of the distance between the voids and the directional growth of the voids. The critical inter-void ligament distance marking the onset of coalescence is shown to be approximately one void radius in both measures., Comment: Accepted for publication in Phys. Rev. Lett

We study the nucleation of GaN islands grown by plasma-assisted molecular-beam epitaxy on AlN(0001) in a Stranski-Krastanov mode. In particular, we assess the variation of their height and density as a function of GaN coverage. We show that the GaN growth passes four stages: initially, the growth is layer-by-layer; subsequently, two-dimensional precursor islands form, which transform into genuine three-dimensional islands. During the latter stage, island height and density increase with GaN coverage until the density saturates. During further GaN growth, the density remains constant and a bimodal height distribution appears. The variation of island height and density as a function of substrate temperature is discussed in the framework of an equilibrium model for Stranski-Krastanov growth., Comment: Submitted to PRB, 10 pages, 15 figures

The effect of stress-triaxiality on growth of a void in a three dimensional single-crystal face-centered-cubic (FCC) lattice has been studied. Molecular dynamics (MD) simulations using an embedded-atom (EAM) potential for copper have been performed at room temperature and using strain controlling with high strain rates ranging from 10^7/sec to 10^10/sec. Strain-rates of these magnitudes can be studied experimentally, e.g. using shock waves induced by laser ablation. Void growth has been simulated in three different conditions, namely uniaxial, biaxial, and triaxial expansion. The response of the system in the three cases have been compared in terms of the void growth rate, the detailed void shape evolution, and the stress-strain behavior including the development of plastic strain. Also macroscopic observables as plastic work and porosity have been computed from the atomistic level. The stress thresholds for void growth are found to be comparable with spall strength values determined by dynamic fracture experiments. The conventional macroscopic assumption that the mean plastic strain results from the growth of the void is validated. The evolution of the system in the uniaxial case is found to exhibit four different regimes: elastic expansion; plastic yielding, when the mean stress is nearly constant, but the stress-triaxiality increases rapidly together with exponential growth of the void; saturation of the stress-triaxiality; and finally the failure., Comment: 35 figures, which are small (and blurry) due to the space limitations; submitted (with original figures) to Physical Review B. Final version

We describe in greater detail the exactly solvable microscopic model we have developed for analyzing the strain-mediated interaction of vacancy lines in a pseudomorphic adsorbate system (Phys. Rev. Lett., to appear). The model is applied to Ga/Si(112) by extracting values for the microscopic parameters from total-energy calculations. The results, which are in good agreement with experimental observations, reveal an unexpectedly complex interplay between compressive and tensile strain within the mixed Ga-Si surface layer., Comment: 4 pages, 4 figures, to appear in Phys. Rev. Lett; this version includes a 4-page appendix giving the full solution of the Frenkel-Kontorova model