Shengfeng Yang

Nanoscale toughening mechanism of nacre tablet.

Nacre has attracted widespread interest because its unique hierarchical structure, which is assembled by 95 wt% brittle aragonite and 5 wt% soft organic materials, leads to several orders of improvement in fracture toughness. Apart from the well proposed toughening mechanisms such as mineral bridges and tablets interlocks, the organic materials including biopolymers between tablets and proteins exist within a tablet can also potentially improve the toughness. In this work, we employ a novel approach combining steered molecular dynamics (SMD) and classical molecular dynamics (MD) to build a model of mineral-protein composite to mimic nacre tablet. The critical role of protein in improving the fracture toughness of nacre is investigated for the first time. MD simulations of single crystalline aragonite, polycrystalline aragonite and mineral-protein composite under uniaxial tensile loading are performed, and the obtained constitutive responses are compared with experimental measurements of nacre under tension. It is shown that the fracture toughness of mineral-protein composite is significantly larger than that of single crystalline or polycrystalline aragonite. Detailed atomic configuration analyses reveal that the fracture of individual computer model is governed by its unique failure mechanisms. Dislocation motion and phase transformation are observed during the failure of single crystalline aragonite. Polycrystalline aragonite fails by the inter-granular cleavage, as well as phase transformation within grain. It is surprisingly noted that other than the stretching of protein chains on grain boundaries, intra-granular fracture is triggered in mineral-protein composites. Proteins serve as strong glue between the inorganic nanograins. It is believed that the strong electrostatic interaction between protein and aragonite nanograins, combined with the remarkable plastic ductility of protein lead to the intra-granular failure, which consequently enhance the fracture toughness of the whole specimen.

Concurrent atomistic and continuum simulation of bi-crystal strontium titanate with tilt grain boundary.

In this paper, we present the development of a concurrent atomistic–continuum (CAC) methodology for simulation of the grain boundary (GB) structures and their interaction with other defects in ionic materials. Simulation results show that the CAC simulation allows a smooth passage of cracks through the atomistic–continuum interface without the need for additional constitutive rules or special numerical treatment; both the atomic-scale structures and the energies of the four different [001] tilt GBs in bi-crystal strontium titanate obtained by CAC compare well with those obtained by existing experiments and density function theory calculations. Although 98.4% of the degrees of freedom of the simulated atomistic system have been eliminated in a coarsely meshed finite-element region, the CAC results, including the stress–strain responses, the GB–crack interaction mechanisms and the effect of the interaction on the fracture strength, are comparable with that of all-atom molecular dynamics simulation results. In addition, CAC simulation results show that the GB–crack interaction has a significant effect on the fracture behaviour of bi-crystal strontium titanate; not only the misorientation angle but also the atomic-level details of the GB structure influence the effect of the GB on impeding crack propagation.

Concurrent atomistic–continuum simulation of polycrystalline strontium titanate

This paper presents the new development of a concurrent atomistic–continuum (CAC) method in simulation of the dynamic evolution of defects in polycrystalline polyatomic materials. The CAC method is based on a theoretical formulation that extends Kirkwood’s statistical mechanical theory of transport processes to a multiscale description of crystalline materials. It solves for both the deformation of lattice cells and the internal deformation within each lattice cell, making it a suitable method for simulations of polyatomic materials. The simulation results of this work demonstrate that CAC can simulate the nucleation of dislocations and cracks from atomistically resolved grain boundary (GB) regions and the subsequent propagation into coarsely meshed grain interiors in polycrystalline strontium titanate without the need of supplemental constitutive equations or additional numerical treatments. With a significantly reduced computational cost, CAC predicts not only the GB structures, but also the dynamic behaviour of dislocations, cracks and GBs, all of which are comparable with those obtained from atomic-level molecular dynamics simulations. Simulation results also show that dislocations tend to initiate from GBs and triple junctions. The angle between the slip planes and the GB planes plays a key role in determining the GB-dislocation reactions.

Phonon thermal transport through tilt grain boundaries in strontium titanate

In this work, we perform nonequilibrium molecular dynamics simulations to study phonon scattering at two tilt grain boundaries (GBs) in SrTiO3. Mode-wise energy transmission coefficients are obtained based on phonon wave-packet dynamics simulations. The Kapitza conductance is then quantified using a lattice dynamics approach. The obtained results of the Kapitza conductance of both GBs compare well with those obtained by the direct method, except for the temperature dependence. Contrary to common belief, the results of this work show that the optical modes in SrTiO3 contribute significantly to phonon thermal transport, accounting for over 50% of the Kapitza conductance. To understand the effect of the GB structural disorder on phonon transport, we compare the local phonon density of states of the atoms in the GB region with that in the single crystalline grain region. Our results show that the excess vibrational modes introduced by the structural disorder do not have a significant effect on phonon scatteri...

Ballistic-diffusive phonon heat transport across grain boundaries

Abstract The propagation of a heat pulse in a single crystal and across grain boundaries (GBs) is simulated using a concurrent atomistic-continuum method furnished with a coherent phonon pulse model. With a heat pulse constructed based on a Bose-Einstein distribution of phonons, this work has reproduced the phenomenon of phonon focusing in single and polycrystalline materials. Simulation results provide visual evidence that the propagation of a heat pulse in crystalline solids with or without GBs is partially ballistic and partially diffusive, i.e., there is a co-existence of ballistic and diffusive thermal transport, with the long-wavelength phonons traveling ballistically while the short-wavelength phonons scatter with each other and travel diffusively. To gain a quantitative understanding of GB thermal resistance, the kinetic energy transmitted across GBs is monitored on the fly and the time-dependent energy transmission for each specimen is measured; the contributions of coherent and incoherent phonon transport to the energy transmission are estimated. Simulation results reveal that the presence of GBs modifies the nature of thermal transport, with the coherent long-wavelength phonons dominating the heat conduction in materials with GBs. In addition, it is found that phonon-GB interactions can result in reconstruction of GBs.

Concurrent Atomistic-Continuum Simulation of Defects in Polyatomic Ionic Materials

This chapter reviews a concurrent atomistic-continuum (CAC) method for studying the dynamic behavior of defects in polyatomic materials. The CAC method combines a unified atomistic and continuum formulation of balance laws and a modified finite element method. In this chapter, we show that in the dynamic simulations of strontium titanate (SrTiO3), the CAC method allows the passages of defects, including dislocations and cracks, from the atomistic to the continuum domain without the need for any special numerical treatment. Then the CAC simulation results of the dynamics of defects in single-crystal, bi-crystal, and polycrystalline SrTiO3 are presented. Simulation results show that CAC successfully reproduces crack propagations, dislocation migrations, and the interplay between cracks, dislocations, and grain boundaries in SrTiO3, with a significant reduction in the degrees of freedom compared with atomistically resolved molecular dynamics (MD) simulation. Most importantly, the essential atomistic features of defects, such as the propagation paths of cracks and dislocations, the dissociation of dislocations, the atomic-scale grain boundary (GB) structures, and atomic details for GB interactions with other defects, are retained in the CAC simulations.

Modeling and analysis of the thermal effects of a circular bimorph piezoelectric actuator.

A theoretical analysis of the thermal effects of a circular bimorph piezoelectric actuator (CBPA) was performed. The circular bimorph structure consists of two flexible piezoelectric ceramic layers and one metallic layer in the middle, and is powered to produce flexural deformation. The CBPA, which may be a good match for large adaptive optics telescopes, has a large stroke and a high resonance frequency. We have derived analytical solutions (both the static solution and the dynamic solution) of the thermal effects of introducing (and increasing the thickness of) a metallic layer into the bimorph. Numerical results are presented to illustrate the dependence of the CBPAs performance upon the physical parameters.

Role of disordered bipolar complexions on the sulfur embrittlement of nickel general grain boundaries

Minor impurities can cause catastrophic fracture of normally ductile metals. Here, a classic example is represented by the sulfur embrittlement of nickel, whose atomic-level mechanism has puzzled researchers for nearly a century. In this study, coupled aberration-corrected electron microscopy and semi-grand-canonical-ensemble atomistic simulation reveal, unexpectedly, the universal formation of amorphous-like and bilayer-like facets at the same general grain boundaries. Challenging the traditional view, the orientation of the lower-Miller-index grain surface, instead of the misorientation, dictates the interfacial structure. We also find partial bipolar structural orders in both amorphous-like and bilayer-like complexions (a.k.a. thermodynamically two-dimensional interfacial phases), which cause brittle intergranular fracture. Such bipolar, yet largely disordered, complexions can exist in and affect the properties of various other materials. Beyond the embrittlement mechanism, this study provides deeper insight to better understand abnormal grain growth in sulfur-doped Ni, and generally enriches our fundamental understanding of performance-limiting and more disordered interfaces.Sulfur at nickel grain boundaries can cause catastrophic failure, but the mechanisms behind that embrittlement remain poorly understood. Here, the authors image and model bipolar sulfur–nickel structures at amorphous-like and bilayer-like facets of general grain boundaries that cause embrittlement.