On the incompatibility of strains and its application to mesoscopic studies of plasticity

Structural transitions are invariably affected by lattice distortions. If the body is to remain crack-free, the strain field cannot be arbitrary but has to satisfy the Saint-Venant compatibility constraint. Equivalently, an incompatibility constraint consistent with the actual dislocation network has to be satisfied in media with dislocations. This constraint can be incorporated into strain-based free energy functionals to study the influence of dislocations on phase stability. We provide a systematic analysis of this constraint in three dimensions and show how three incompatibility equations accommodate an arbitrary dislocation density. This approach allows the internal stress field to be calculated for an anisotropic material with spatially inhomogeneous microstructure and distribution of dislocations by minimizing the free energy. This is illustrated by calculating the stress field of an edge dislocation and comparing it with that of an edge dislocation in an infinite isotropic medium. We outline how this procedure can be utilized to study the interaction of plasticity with polarization and magnetization.Comment: 6 pages, 2 figures; will appear in Phys. Rev.

The influence of transition metal solutes on dislocation core structure and values of Peierls stress and barrier in tungsten

Several transition metals were examined to evaluate their potential for improving the ductility of tungsten. The dislocation core structure and Peierls stress and barrier of $1/2$ screw dislocations in binary tungsten-transition metal alloys (W$_{1-x}$TM$_{x}$) were investigated using first principles electronic structure calculations. The periodic quadrupole approach was applied to model the structure of $1/2$ dislocation. Alloying with transition metals was modeled using the virtual crystal approximation and the applicability of this approach was assessed by calculating the equilibrium lattice parameter and elastic constants of the tungsten alloys. Reasonable agreement was obtained with experimental data and with results obtained from the conventional supercell approach. Increasing the concentration of a transition metal from the VIIIA group, i.e. the elements in columns headed by Fe, Co and Ni, leads to reduction of the $C^\prime$ elastic constant and increase of elastic anisotropy A=$C_{44}/C^\prime$. Alloying W with a group VIIIA transition metal changes the structure of the dislocation core from symmetric to asymmetric, similar to results obtained for W$_{1-x}$Re$_{x}$ alloys in the earlier work of Romaner {\it et al} (Phys. Rev. Lett. 104, 195503 (2010))\comments{\cite{WRECORE}}. In addition to a change in the core symmetry, the values of the Peierls stress and barrier are reduced. The latter effect could lead to increased ductility in a tungsten-based alloy\comments{\cite{WRECORE}}. Our results demonstrate that alloying with any of the transition metals from the VIIIA group should have similar effect as alloying with Re.Comment: 12 pages, 8 figures, 3 table

Assessment of interatomic potentials for atomistic analysis of static and dynamic properties of screw dislocations in W

Screw dislocations in bcc metals display non-planar cores at zero temperature which result in high lattice friction and thermally activated strain rate behavior. In bcc W, electronic structure molecular statics calculations reveal a compact, non-degenerate core with an associated Peierls stress between 1.7 and 2.8 GPa. However, a full picture of the dynamic behavior of dislocations can only be gained by using more efficient atomistic simulations based on semiempirical interatomic potentials. In this paper we assess the suitability of five different potentials in terms of static properties relevant to screw dislocations in pure W. As well, we perform molecular dynamics simulations of stress-assisted glide using all five potentials to study the dynamic behavior of screw dislocations under shear stress. Dislocations are seen to display thermally-activated motion in most of the applied stress range, with a gradual transition to a viscous damping regime at high stresses. We find that one potential predicts a core transformation from compact to dissociated at finite temperature that affects the energetics of kink-pair production and impacts the mechanism of motion. We conclude that a modified embedded-atom potential achieves the best compromise in terms of static and dynamic screw dislocation properties, although at an expense of about ten-fold compared to central potentials

Stability estimates for resolvents, eigenvalues and eigenfunctions of elliptic operators on variable domains

We consider general second order uniformly elliptic operators subject to homogeneous boundary conditions on open sets $\phi (\Omega)$ parametrized by Lipschitz homeomorphisms $\phi$ defined on a fixed reference domain $\Omega$. Given two open sets $\phi (\Omega)$, $\tilde \phi (\Omega)$ we estimate the variation of resolvents, eigenvalues and eigenfunctions via the Sobolev norm $\|\tilde \phi -\phi \|_{W^{1,p}(\Omega)}$ for finite values of $p$, under natural summability conditions on eigenfunctions and their gradients. We prove that such conditions are satisfied for a wide class of operators and open sets, including open sets with Lipschitz continuous boundaries. We apply these estimates to control the variation of the eigenvalues and eigenfunctions via the measure of the symmetric difference of the open sets. We also discuss an application to the stability of solutions to the Poisson problem.Comment: 34 pages. Minor changes in the introduction and the refercenes. Published in: Around the research of Vladimir Maz'ya II, pp23--60, Int. Math. Ser. (N.Y.), vol. 12, Springer, New York 201

The mixed problem for the Laplacian in Lipschitz domains

We consider the mixed boundary value problem or Zaremba's problem for the Laplacian in a bounded Lipschitz domain in R^n. We specify Dirichlet data on part of the boundary and Neumann data on the remainder of the boundary. We assume that the boundary between the sets where we specify Dirichlet and Neumann data is a Lipschitz surface. We require that the Neumann data is in L^p and the Dirichlet data is in the Sobolev space of functions having one derivative in L^p for some p near 1. Under these conditions, there is a unique solution to the mixed problem with the non-tangential maximal function of the gradient of the solution in L^p of the boundary. We also obtain results with data from Hardy spaces when p=1.Comment: Version 5 includes a correction to one step of the main proof. Since the paper appeared long ago, this submission includes the complete paper, followed by a short section that gives the correction to one step in the proo

Modelling avalanches in martensites

Solids subject to continuous changes of temperature or mechanical load often exhibit discontinuous avalanche-like responses. For instance, avalanche dynamics have been observed during plastic deformation, fracture, domain switching in ferroic materials or martensitic transformations. The statistical analysis of avalanches reveals a very complex scenario with a distinctive lack of characteristic scales. Much effort has been devoted in the last decades to understand the origin and ubiquity of scale-free behaviour in solids and many other systems. This chapter reviews some efforts to understand the characteristics of avalanches in martensites through mathematical modelling.Comment: Chapter in the book "Avalanches in Functional Materials and Geophysics", edited by E. K. H. Salje, A. Saxena, and A. Planes. The final publication is available at Springer via http://dx.doi.org/10.1007/978-3-319-45612-6_

Bond-order potential for simulations of extended defects in tungsten

We present a bond-order potential (BOP) for the bcc transition metal tungsten. The bond-order potentials are a real-space semiempirical scheme for the description of interatomic interactions based on the tight-binding approximation. In the hierarchy of atomic-scale-modeling methods the BOPs thus provide a direct bridge between electronic-structure and atomistic techniques. Two variants of the BOP were constructed and extensively tested against accurate first-principles methods in order to assess the potentials\u27 reliability and applicability. A comparison of the BOP with a central-force potential is used to demonstrate that a correct description of directional mixed covalent and metallic bonds is crucial for a successful and fully transferable model. The potentials are applied in studies of low-index surfaces, symmetrical tilt grain boundaries, and dislocations

Efficient CO2-Reducing Activity of NAD-Dependent Formate Dehydrogenase from Thiobacillus sp KNK65MA for Formate Production from CO2 Gas

NAD-dependent formate dehydrogenase (FDH) from Candida boidinii (CbFDH) has been widely used in various CO2 reduction systems but its practical applications are often impeded due to low CO2-reducing activity. In this study, we demonstrated superior CO2-reducing properties of FDH from Thiobacillus sp. KNK65MA (TsFDH) for production of formate from CO2 gas. To discover more efficient CO2-reducing FDHs than a reference enzyme e. CbFDH, five FDHs were selected with biochemical properties and then, their CO2-reducing activities were evaluated. All FDHs including CbFDH showed better CO2-reducing activities at acidic pHs than at neutral pHs and four FDHs were more active than CbFDH in the CO2 reduction reaction. In particular, the FDH from Thiobacillus sp. KNK65IVIA (TsFDH) exhibited the highest CO2-reducing activity and had a dramatic preference for the reduction reaction, i.e., a 84.2-fold higher ratio of CO2 reduction to formate oxidation in catalytic efficiency (k(cat)/K-B) compared to CbFDH. Formate was produced from CO2 gas using TsFDH and CbFDH, and TsFDH showed a 5.8-fold higher formate production rate than CbFDH. A sequence and structural comparison showed that FDHs with relatively high CO2-reducing activities had elongated N- and C-terminal loops. The experimental results demonstrate that TsFDH can be an alternative to CbFDH as a biocatalyst in CO2 reduction systemsope

Bilevel Parameter Learning for Higher-Order Total Variation Regularisation Models.

We consider a bilevel optimisation approach for parameter learning in higher-order total variation image reconstruction models. Apart from the least squares cost functional, naturally used in bilevel learning, we propose and analyse an alternative cost based on a Huber-regularised TV seminorm. Differentiability properties of the solution operator are verified and a first-order optimality system is derived. Based on the adjoint information, a combined quasi-Newton/semismooth Newton algorithm is proposed for the numerical solution of the bilevel problems. Numerical experiments are carried out to show the suitability of our approach and the improved performance of the new cost functional. Thanks to the bilevel optimisation framework, also a detailed comparison between TGV 2 and ICTV is carried out, showing the advantages and shortcomings of both regularisers, depending on the structure of the processed images and their noise level.King Abdullah University of Science and Technology (KAUST) (Grant ID: KUKI1-007-43), Engineering and Physical Sciences Research Council (Grant IDs: Nr. EP/J009539/1 “Sparse & Higher-order Image Restoration” and Nr. EP/M00483X/1 “Efficient computational tools for inverse imaging problems”), Escuela Politécnica Nacional de Quito (Grant ID: PIS 12-14, MATHAmSud project SOCDE “Sparse Optimal Control of Differential Equations”), Leverhulme Trust (project on “Breaking the non-convexity barrier”), SENESCYT (Ecuadorian Ministry of Higher Education, Science, Technology and Innovation) (Prometeo Fellowship)This is the final version of the article. It first appeared from Springer via http://dx.doi.org/10.1007/s10851-016-0662-