The recently introduced coupled cluster (CC) downfolding techniques for reducing the dimensionality of quantum many-body problems recast the CC formalism in the form of the renormalization procedure allowing, for the construction of effective (or downfolded) Hamiltonians in small-dimensionality sub-space, usually identified with the so-called active space, of the entire Hilbert space. The resulting downfolded Hamiltonians integrate out the external (out-of-active-space) Fermionic degrees of freedom from the internal (in-the-active-space) parameters of the wave function, which can be determined as components of the eigenvectors of the downfolded Hamiltonians in the active space. This paper will discuss the extension of non-Hermitian (associated with standard CC formulations) and Hermitian (associated with the unitary CC approaches) downfolding formulations to composite quantum systems commonly encountered in materials science and chemistry. The non-Hermitian formulation can provide a platform for developing local CC approaches, while the Hermitian one can serve as an ideal foundation for developing various quantum computing applications based on the limited quantum resources. We also discuss the algorithm for extracting the semi-analytical form of the inter-electron interactions in the active spaces.

Computational methods are increasingly being incorporated into the exploitation of microstructure–property relationships for microstructure-sensitive design of materials. In the present work, we propose non-intrusive materials informatics methods for the high-throughput exploration and analysis of a synthetic microstructure space using a machine learning-reinforced multi-phase-field modeling scheme. We specifically study the interface energy space as one of the most uncertain inputs in phase-field modeling and its impact on the shape and contact angle of a growing phase during heterogeneous solidification of secondary phase between solid and liquid phases. We evaluate and discuss methods for the study of sensitivity and propagation of uncertainty in these input parameters as reflected on the shape of the Cu6Sn5 intermetallic during growth over the Cu substrate inside the liquid Sn solder due to uncertain interface energies. The sensitivity results rank σSI,σIL, and σIL, respectively, as the most influential parameters on the shape of the intermetallic. Furthermore, we use variational autoencoder, a deep generative neural network method, and label spreading, a semi-supervised machine learning method for establishing correlations between inputs of outputs of the computational model. We clustered the microstructures into three categories (“wetting”, “dewetting”, and “invariant”) using the label spreading method and compared it with the trend observed in the Young-Laplace equation. On the other hand, a structure map in the interface energy space is developed that shows σSI and σSL alter the shape of the intermetallic synchronously where an increase in the latter and decrease in the former changes the shape from dewetting structures to wetting structures. The study shows that the machine learning-reinforced phase-field method is a convenient approach to analyze microstructure design space in the framework of the ICME.

Upon publication of the original manuscript (Prudil and Welland 2017) discrepancies were highlighted in Eq. (50) [equation (1) here].The equation in the original manuscript (Prudil and Welland 2017) erroneously included a factor of 2 on theta and read:$$ \cos \left(2\theta \right)=\frac{\sigma_{\alpha \gamma }}{2{\sigma}_{\alpha \beta }} $$(1)This equation should read as below:$$ \cos \left(\theta \right)=\frac{\sigma_{\alpha \gamma }}{2{\sigma}_{\alpha \beta }} $$(1)This error has since been acknowledged and corrected in this erratum.Prudil AA, Welland MJ. A novel model of third phase inclusions on two phase boundaries. Materials Theory. 2017 Dec 1;1(1):4Download referencesAffiliationsFuel & Fuel Channel Safety Branch, Canadian Nuclear Laboratories, Chalk River, Ontario, CanadaAndrew A. Prudil & Michael J. WellandAuthorsAndrew A. PrudilView author publicationsYou can also search for this author in PubMed Google ScholarMichael J. WellandView author publicationsYou can also search for this author in PubMed Google ScholarCorresponding authorCorrespondence to Andrew A. Prudil.The online version of the original article can be found under doi:10.1186/s41313-017-0003-3 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Reprints and PermissionsCite this articlePrudil, A.A., Welland, M.J. Erratum to: a novel model of third phase inclusions on two phase boundaries. Mater Theory 1, 7 (2017). http://doi.org/10.1186/s41313-017-0007-zDownload citationReceived: 12 July 2017Accepted: 12 July 2017Published: 28 July 2017DOI: http://doi.org/10.1186/s41313-017-0007-z

The standard way of modeling plasticity in polycrystals is by using the crystal plasticity model for single crystals in each grain, and imposing suitable traction and slip boundary conditions across grain boundaries. In this fashion, the system is modeled as a collection of boundary-value problems with matching boundary conditions. In this paper, we develop a diffuse-interface crystal plasticity model for polycrystalline materials that results in a single boundary-value problem with a single crystal as the reference configuration. Using a multiplicative decomposition of the deformation gradient into lattice and plastic parts, i.e. F(X,t)=F L(X,t)F P(X,t), an initial stress-free polycrystal is constructed by imposing F L to be a piecewise constant rotation field R 0(X), and F P=R 0(X)T, thereby having F(X,0)=I, and zero elastic strain. This model serves as a precursor to higher order crystal plasticity models with grain boundary energy and evolution.

We present a continuum model of ion-induced surface patterning. The model incorporates the atomic processes of sputtering, re-deposition and surface diffusion, and is shown to display the generic features of the damped Kuramoto-Sivashinsky (KS) equation of non-linear dynamics. Linear and non-linear stability analyses of the evolution equation give estimates of the emerging pattern wavelength and spatial symmetry. The analytical theory is confirmed by numerical simulations of the evolution equation with the Fast Fourier Transform method, where we show the influence of the incident ion angle, flux, and substrate surface temperature. It is shown that large local geometry variations resulting in quadratic non-linearities in the evolution equation dominate pattern selection and stability at long time scales.

Direct Molecular Dynamics (MD) simulations are being increasingly employed to model dislocation-mediated crystal plasticity with atomic resolution. Thanks to the dislocation extraction algorithm (DXA), dislocation lines can be now accurately detected and positioned in space and their Burgers vector unambiguously identified in silico, while the simulation is being performed. However, DXA extracts static snapshots of dislocation configurations that by themselves present no information on dislocation motion. Referred to as a sweep-tracing algorithm (STA), here we introduce a practical computational method to observe dislocation motion and to accurately quantify its important characteristics such as preferential slip planes (slip crystallography). STA reconnects pairs of successive snapshots extracted by DXA and computes elementary slip facets thus precisely tracing the motion of dislocation segments from one snapshot to the next. As a testbed for our new method, we apply STA to the analysis of dislocation motion in large-scale MD simulations of single crystal plasticity in BCC metals. We observe that, when the crystal is subjected to uniaxial deformation along its [001] axis, dislocation slip predominantly occurs on the {112} maximum resolved shear stress plane under tension, while in compression slip is non-crystallographic (pencil) resulting in asymmetric mechanical response. The marked contrast in the observed slip crystallography is attributed to the twinning/anti-twinning asymmetry of shears in the {112} planes relatively favoring dislocation motion in the twinning sense while hindering dislocations from moving in the anti-twinning directions.

Cracks are created by massive breakage of molecular or atomic bonds. The latter, in its turn, leads to the highly localized loss of material, which is the reason why even closed cracks are visible by a naked eye. Thus, fracture can be interpreted as the local material sink. Mass conservation is violated locally in the area of material failure. We consider a theoretical formulation of the coupled mass and momenta balance equations for a description of fracture. Our focus is on brittle fracture and we propose a finite strain hyperelastic thermodynamic framework for the coupled mass-flow-elastic boundary value problem. The attractiveness of the proposed framework as compared to the traditional continuum damage theories is that no internal parameters (like damage variables, phase fields, etc.) are used while the regularization of the failure localization is provided by the physically sound law of mass balance.

A new computationally efficient model of an included phase located at the interface between two other phases is developed by projecting the boundaries of the inclusion onto the boundary between the two other phases. This reduces the 3D problem to one on a 2D surface while still being embedded in 3D space, which significantly reduces computational expense of solving the system. The resulting model is similar to conventional phase-field models. The properties of the solution are examined, compared to classical theory, and the numerical behaviour, including a mesh sensitivity analysis, are discussed. The model accurately captures mesoscale effects, such as the Gibbs-Thompson effect, coarsening, and coalescence. An example application of the model simulating the evolution of grain boundary porosity in nuclear fuel is shown on a representative tetrakaidecahedron-shaped fuel grain.

The purpose of the current work is the formulation of models for conservative and non-conservative dynamics in solid systems with the help of the General Equation for the Non-Equilibrium Reversible-Irreversible Coupling (GENERIC: e.g., Grmela and Öttinger, Phys. Rev. E 56(6), 6620 (1997); Öttinger and Grmela, Phys. Rev. E 56(6), 6633 (1997)). In this context, the resulting models are inherently spatially strongly non-local (i.e., functional) and non-isothermal in character. They are applicable in particular to the modeling of phase transitions as well as mass and heat transport in multiphase, multicomponent solids. In the last part of the work, the strongly non-local model formulation is reduced to weakly non-local form with the help of generalized gradient approximation of the energy and entropy functionals. On this basis, the current model formulation is shown to be consistent with and reduce to a recent non-isothermal generalization (Gladkov et al., J. Non-Equilib. Thermodyn. 41(2), 131 (2016)) of the well-known phase-field models of Cahn and Hilliard (J. Chem. Phys. 28(2), 258 (1958)) for conservative dynamics and of Allen and Cahn (Acta Metall. 27(6), 1085 (1979)) for non-conservative dynamics. Finally, the current approach is applied to derive a non-isothermal generalization of a phase-field crystal model for binary alloys (see, e.g., Elder et al., Phys. Rev. B 75(6), 064107 (2007)).

Coarse-grained descriptions of dislocation motion in crystalline metals inherently represent a loss of information regarding dislocation-dislocation interactions. In the present work, we consider a coarse-graining framework capable of re-capturing these interactions by means of the dislocation-dislocation correlation functions. The framework depends on a convolution length to define slip-system-specific dislocation densities. Following a statistical definition of this coarse-graining process, we define a spatial correlation function which will allow the arrangement of the discrete line system at two points—and thus the strength of their interactions at short range—to be recaptured into a mean field description of dislocation dynamics. Through a statistical homogeneity argument, we present a method of evaluating this correlation function from discrete dislocation dynamics simulations. Finally, results of this evaluation are shown in the form of the correlation of dislocation densities on the same slip-system. These correlation functions are seen to depend weakly on plastic strain, and in turn, the dislocation density, but are seen to depend strongly on the convolution length. Implications of these correlation functions in regard to continuum dislocation dynamics as well as future directions of investigation are also discussed.

Macroscopic properties of structural materials are strongly dependent on their microstructure. However, the modeling of their evolution is a complex task because of the mechanisms involved such as plasticity, recrystallization, and phase transformations, which are common processes taking place in metallic alloys. This complexity led to a growing interest in atomistic simulations formulated without any auxiliary hypotheses beyond the choice of interatomic potential. In this context, we propose here a model based on an overdamped stochastic evolution of particles interacting through inter-atomic forces. The model settles to the correct thermal equilibrium distribution in canonical and grand-canonical ensembles and is used to study the grain boundary migration. Finally, a comparison of our results with those obtained by molecular dynamics shows that our approach reproduces the complex atomic-scale dynamics of grain boundary migration correctly.

We derive the Green tensor of Mindlin’s anisotropic first strain gradient elasticity. The Green tensor is valid for arbitrary anisotropic materials, with up to 21 elastic constants and 171 gradient elastic constants in the general case of triclinic media. In contrast to its classical counterpart, the Green tensor is non-singular at the origin, and it converges to the classical tensor a few characteristic lengths away from the origin. Therefore, the Green tensor of Mindlin’s first strain gradient elasticity can be regarded as a physical regularization of the classical anisotropic Green tensor. The isotropic Green tensor and other special cases are recovered as particular instances of the general anisotropic result. The Green tensor is implemented numerically and applied to the Kelvin problem with elastic constants determined from interatomic potentials. Results are compared to molecular statics calculations carried out with the same potentials.

Plasticity modelling has long relied on phenomenological models based on ad-hoc assumption of constitutive relations, which are then fitted to limited data. Other work is based on the consideration of physical mechanisms which seek to establish a physical foundation of the observed plastic deformation behavior through identification of isolated defect processes (’mechanisms’) which are observed either experimentally or in simulations and then serve to formulate so-called physically based models. Neither of these approaches is adequate to capture the complexity of plastic deformation which belongs into the realm of emergent collective phenomena, and to understand the complex interplay of multiple deformation pathways which is at the core of modern high performance structural materials. Data based approaches offer alternative pathways towards plasticity modelling whose strengths and limitations we explore here for a simple example, namely the interplay between rate and dislocation density dependent strengthening mechanisms in fcc metals.

In the Portevin–Le Chatelier (PLC) effect sample plastic deformation takes place via localized bands. We present a model to account for band dynamics and the variability the bands exhibit. The approach is tuned to account for strain hardening and the strain-rate dependence for the case of so-called type A (propagating) bands. The main experimental features of the fluctuations are a reduction with strain and increase with the strain rate which is reproduced by a model of plastic deformation with Dynamic Strain Aging, including disorder as a key parameter. Extensions are discussed as are the short-comings in reproducing detailed avalanche statistics.

Growth and other dynamical processes in soft materials can create novel types of mesoscopic defects including discontinuities for the second and higher derivatives of the deformation, and terminating defects for these discontinuities. These higher-order defects move “easily", and can thus confer a great degree of flexibility to the material. We develop a general continuum mechanical framework from which we can derive the dynamics of higher order defects in a thermodynamically consistent manner. We illustrate our framework by obtaining the explicit dynamical equations for the next higher order defects in an elastic body beyond dislocations, phase boundaries, and disclinations, namely, surfaces of inflection and branch lines.

Molecular science is governed by the dynamics of electrons and atomic nuclei, and by their interactions with electromagnetic fields. A faithful physicochemical understanding of these processes is crucial for the design and synthesis of chemicals and materials of value for our society and economy. Although some problems in this field can be adequately addressed by classical mechanics, many demand an explicit quantum mechanical description. Such quantum problems require a representation of wave functions that grows exponentially with system size and therefore should naturally benefit from quantum computation on a number of logical qubits that scales only linearly with system size. In this perspective, we elaborate on the potential benefits of quantum computing in the molecular sciences, i.e., in molecular physics, chemistry, biochemistry, and materials science.

A mesoscale model is introduced to study the dynamics of material defects lying at interface junctions. The proposed framework couples the dynamics of discrete dislocation and disclination lines. Disclinations are expected to be natural defects at interface junctions; their presence serving the purpose of accommodating discontinuities in rotation fields at material interface junctions. Crystallography-based rules are proposed to describe the kinematics of disclination motion. A discrete-continuous couple-stress framework, in which discrete defect lines are introduced as plastic eigenstrains and eigencurvatures, is proposed to explicitly follow the dynamics of interfacial defects. The framework is then applied to study $\left (10\bar {1}2\right)$ twin transverse propagation and thickening in magnesium. Focusing first on the case of a twin domain, It is shown that a disclination based representation of twin domains allows for an appropriate mechanistic description of the kinematics of shear transformations. In what concerns twin thickening, the stability of defects at twin interfaces is further studied. To this end, a 3D crater lying on a twin interface is described as a dipole of disclination loops. Upon self-relaxation, it is found that out of plane motion of disclinations followed by the nucleation of twinning dislocations can be activated; thereby showing that conservative non-planar motion of disclinations can be thermodynamically favorable; mechanism that had been postulated some 50 years ago.

An electrochemical cell contains three open thermodynamic systems that, in dynamic equilibrium, equalize their electrochemical potentials with that of their surrounding by forming an electric-double-layer-capacitor at the interface of the electrolyte with each of the two electrodes. Since the electrode/electrolyte interfaces are heterojunctions, the electrochemical potentials or Fermi levels of the two materials that contact the electrolyte at the two electrodes determine the voltage of a cell. The voltage is the sum of the voltages of the two interfacial electric-double-layer capacitors at the two electrode/electrolyte interfaces. A theoretical analysis of the thermodynamics that gives a quantitative prediction of the observed voltages in an asymmetric cell with an S8 relay at the positive electrode is provided. In addition, new discharge data and an X-ray photoelectron spectroscopy analysis of the lithium plated on the positive electrode of a discharged cell is presented. Ab initio, DFT methods were used to calculate the band structure and surface-state energies of the crystalline S8 solid sulfur relay. The theoretical exposition of the thermodynamics of the operative driving force of the chemical reactions in an electrochemical cell demonstrate that our initial experimental data and conclusions are valid. Other reported observations of lithium plating on the positive electrode, observations that were neither exploited nor their origins specified, are also cited.

The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, fV, was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p. The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of fV vs. p did not change significantly. The percolation threshold pc was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Strain energy decomposition methods in phase field fracture models separate strain energy that contributes to fracture from that which does not. However, various decomposition methods have been proposed in the literature, and it can be difficult to determine an appropriate method for a given problem. The goal of this work is to facilitate the choice of strain decomposition method by assessing the performance of three existing methods (spectral decomposition of the stress or the strain and deviatoric decomposition of the strain) and one new method (deviatoric decomposition of the stress) with several benchmark problems. In each benchmark problem, we compare the performance of the four methods using both qualitative and quantitative metrics. In the first benchmark, we compare the predicted mechanical behavior of cracked material. We then use four quasi-static benchmark cases: a single edge notched tension test, a single edge notched shear test, a three-point bending test, and a L-shaped panel test. Finally, we use two dynamic benchmark cases: a dynamic tensile fracture test and a dynamic shear fracture test. All four methods perform well in tension, the two spectral methods perform better in compression and with mixed mode (though the stress spectral method performs the best), and all the methods show minor issues in at least one of the shear cases. In general, whether the strain or the stress is decomposed does not have a significant impact on the predicted behavior.