Source author record

David J. Srolovitz

David J. Srolovitz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

Disconnection-Mediated Twin/Twin-Junction Migration in FCC metals

We present the results of novel, time-resolved, in situ HRTEM observations, molecular dynamics (MD) simulations, and disconnection theory that elucidate the mechanism by which the motion of grain boundaries (GBs) in polycrystalline materials are coupled through disconnection motion/reactions at/adjacent to GB triple junctions (TJs). We focus on TJs composed of a pair of coherent twin boundaries (CTBs) and a Σ9 GB. As for all GBs, disconnection theory implies that multiple modes/local mechanisms for CTB migration are possible and that the mode selection is affected by the nature of the driving force for migration. While we observe (HRTEM and MD) CTB migration through the motion of pure steps driven by chemical potential jump, other experimental observations (and our simulations) show that stress-driven CTB migration occurs through the motion of disconnections with a non-zero Burgers vector; these are pure-step and twinning-partial CTB migration mechanisms. Our experimental observations and simulations demonstrate that the motion of a GB drags its delimiting TJ and may force the motion of the other GBs meeting at the TJ. Our experiments and simulations focus on two types of TJs composed of a pair of CTBs and a Σ9 GB; a 107° TJ readily migrates while a 70° TJ is immobile (experiment, simulation) in agreement with our disconnection theory even though the intrinsic mobilities of the constituent GBs do not depend on TJ-type. We also demonstrate that disconnections may be formed at TJs (chemical potential jump/stress driven) and at GB/free surface junctions (stress-driven).

preprint2021arXiv

Disconnection-Mediated Migration of Interfaces in Microstructures: I. continuum model

A long-standing goal of materials science is to understand, predict and control the evolution of microstructures in crystalline materials. Most microstructure evolution is controlled by interface motion; hence, the establishment of rigorous interface equations of motion is a universal goal of materials science. We present a new model for the motion of arbitrarily curved interfaces that respects the underlying crystallography of the two phases/domains meeting at the interface and is consistent with microscopic mechanisms of interface motion; i.e., disconnection migration (line defects in the interface with step and dislocation character). We derive the equation of motion for interface migration under the influence of a wide range of driving forces. In Part II of this paper [Salvalaglio, Han and Srolovitz, 2021], we implement the interface model and the equation of motion proposed in this paper in a diffuse interface simulation approach for complex morphology and microstructure evolution.

preprint2021arXiv

Disconnection-Mediated Migration of Interfaces in Microstructures: II. diffuse interface simulations

The motion of interfaces is an essential feature of microstructure evolution in crystalline materials. While atomic-scale descriptions provide mechanistic clarity, continuum descriptions are important for understanding microstructural evolution and upon which microscopic features it depends. We develop a microstructure evolution simulation approach that is linked to the underlying microscopic mechanisms of interface migration. We extend the continuum approach describing the disconnection-mediated motion of interfaces introduced in Part I [Han, Srolovitz and Salvalaglio, 2021] to a diffuse interface, phase-field model suitable for large-scale microstructure evolution. A broad range of numerical simulations showcases the capability of the method and the influence of microscopic interface migration mechanisms on microstructure evolution. These include, in particular, the effects of stress and its coupling to interface migration which arises from disconnections, showing how this leads to important differences from classical microstructure evolution represented by mean curvature flow.

preprint2020arXiv

Distribution of Topological Types in Grain-Growth Microstructures

An open question in studying normal grain growth concerns the asymptotic state to which microstructures converge. In particular, the distribution of grain topologies is unknown. We introduce a thermodynamic-like theory to explain these distributions in two- and three-dimensional systems. In particular, a bending-like energy $E_i$ is associated to each grain topology $t_i$, and the probability of observing that particular topology is proportional to $\frac{1}{s(t_i)}e^{-βE_i}$, where $s(t_i)$ is the order of an associated symmetry group and $β$ is a thermodynamic-like constant. We explain the physical origins of this approach, and provide numerical evidence in support.

preprint2020arXiv

Power-law scaling for solid-state dewetting of thin films: an Onsager variational approach

We examine the kinetics of surface diffusion-controlled, solid-state dewetting by consideration of the retraction of the contact in a semi-infinite solid thin film on a flat rigid substrate. The analysis is performed within the framework of the Onsager variational principle applied to surface diffusion-controlled morphology evolution. Based on this approach, we derive a simple, reduced-order model to quantitatively analyse the power-law scaling of the dewetting process. Using asymptotic analysis and numerical simulations for the reduced-order model, we find that the retraction distance grows as the $2/5$ power of time and the height of the ridge, adjacent to the contact, grows as the $1/5$ power of time for late time. While the asymptotic analysis focuses on late time and a relatively simple geometric model, the Onsager approach is applicable to all times and descriptions of the morphology of arbitrary complexity.

preprint2015arXiv

A Theoretical Examination of Diffusive Molecular Dynamics

Diffusive molecular dynamics is a novel model for materials with atomistic resolution that can reach diffusive time scales. The main ideas of diffusive molecular dynamics are to first minimize an approximate variational Gaussian free energy of the system with respect to the mean atomic coordinates (averaging over many vibrational periods), and to then to perform a diffusive step where atoms and vacancies (or two species in a binary alloy) flow on a diffusive time scale via a master equation. We present a mathematical framework for studying this algorithm based upon relative entropy, or Kullback-Leibler divergence. This adds flexibility in how the algorithm is implemented and interpreted. We then compare our formulation, relying on relative entropy and absolute continuity of measures, to existing formulations. The main difference amongst the equations appears in a model for vacancy diffusion, where additional entropic terms appear in our development.

preprint2015arXiv

A Topological Framework for Local Structure Analysis in Condensed Matter

Physical systems are frequently modeled as sets of points in space, each representing the position of an atom, molecule, or mesoscale particle. As many properties of such systems depend on the underlying ordering of their constituent particles, understanding that structure is a primary objective of condensed matter research. Although perfect crystals are fully described by a set of translation and basis vectors, real-world materials are never perfect, as thermal vibrations and defects introduce significant deviation from ideal order. Meanwhile, liquids and glasses present yet more complexity. A complete understanding of structure thus remains a central, open problem. Here we propose a unified mathematical framework, based on the topology of the Voronoi cell of a particle, for classifying local structure in ordered and disordered systems that is powerful and practical. We explain the underlying reason why this topological description of local structure is better suited for structural analysis than continuous descriptions. We demonstrate the connection of this approach to the behavior of physical systems and explore how crystalline structure is compromised at elevated temperatures. We also illustrate potential applications to identifying defects in plastically deformed polycrystals at high temperatures, automating analysis of complex structures, and characterizing general disordered systems.

preprint2015arXiv

Accurate Measurement of the Cleavage Energy of Graphite

The basal plane cleavage energy (CE) of graphite is a key material parameter for understanding many of the unusual properties of graphite, graphene, and carbon nanotubes. The CE is equal to twice the surface energy and is closely related to the interlayer binding energy and exfoliation energy of graphite. Nonetheless, a wide range of values for these properties have been reported and no consensus has yet emerged as to their magnitude. Here, we report the first direct, accurate experimental measurement of the CE of graphite using a novel method based on the recently discovered self-retraction phenomenon in graphite. The measured value, 0.37 +/- 0.01 J/m2 for the incommensurate state of bicrystal graphite, is nearly invariant with respect to temperature (from 22°C to 198°C) and bicrystal twist angle, and insensitive to impurities (from the atmosphere). The cleavage energy for the ideal ABAB graphite stacking, 0.39 +/- 0.02 J/m2, is calculated based upon a combination of the measured CE and a theoretical calculation. These experimental measurements are ideal for use in evaluating the efficacy of competing theoretical approaches.

preprint2015arXiv

Electron-Beam Driven Relaxation Oscillations in Ferroelectric Nanodisks

Using a combination of computational simulations, atomic-scale resolution imaging and phenomenological modelling, we examine the underlying mechanism for nanodomain restructuring in lead zirconate titanate (PZT) nanodisks driven by electron beams. The observed subhertz nanodomain dynamics are identified with relaxation oscillations where the charging/discharging cycle time is determined by saturation of charge traps and nanodomain wall creep. These results are unusual in that they indicate very slow athermal dynamics in nanoscale systems.

preprint2015arXiv

Geometric and topological properties of the canonical grain growth microstructure

Many physical systems can be modeled as large sets of domains "glued" together along boundaries - biological cells meet along cell membranes, soap bubbles meet along thin films, countries meet along geopolitical boundaries, and metallic crystals meet along grain interfaces. Each class of microstructures results from a complex interplay of initial conditions and particular evolutionary dynamics. The statistical steady-state microstructure resulting from isotropic grain growth of a polycrystalline material is canonical in that it is the simplest example of a cellular microstructure resulting from a gradient flow of a simple energy, directly proportional to the total length or area of all cell boundaries. As many properties of polycrystalline materials depend on their underlying microstructure, a more complete understanding of the grain growth steady-state can provide insight into the physics of a broad range of everyday materials. In this paper we report geometric and topological features of these canonical two- and three-dimensional steady-state microstructures obtained through large, accurate simulations of isotropic grain growth.

preprint2015arXiv

Statistical Topology of Perturbed Two-Dimensional Lattices

The Voronoi cell of any atom in a lattice is identical. If atoms are perturbed from their lattice coordinates, then the topologies of the Voronoi cells of the atoms will change. We consider the distribution of Voronoi cell topologies in two-dimensional perturbed systems. These systems can be thought of as simple models of finite-temperature crystals. We give analytical results for the distribution of Voronoi topologies of points in two-dimensional Bravais lattices under infinitesimal perturbations and present a discussion with numerical results for finite perturbations.

preprint2014arXiv

Statistical Topology of Three-Dimensional Poisson-Voronoi Cells and Cell Boundary Networks

Voronoi tessellations of Poisson point processes are widely used for modeling many types of physical and biological systems. In this paper, we analyze simulated Poisson-Voronoi structures containing a total of 250,000,000 cells to provide topological and geometrical statistics of this important class of networks. We also report correlations between some of these topological and geometrical measures. Using these results, we are able to corroborate several conjectures regarding the properties of three-dimensional Poisson-Voronoi networks and refute others. In many cases, we provide accurate fits to these data to aid further analysis. We also demonstrate that topological measures represent powerful tools for describing cellular networks and for distinguishing among different types of networks.

preprint2013arXiv

Depletion layer-induced size effects in ferroelectric thin films: A Ginzburg-Landau model study

A Ginzburg-Landau model is used to demonstrate how depletion layers give rise to thickness-dependent ferroelectric properties in thin films. It is shown that free charge layers at the film-electrode interface can result in an internal electric field in the bulk of the film even when no external voltage is applied. At high values of the donor dopant density and small thicknesses, this internal electric field can be strong enough to lead to the formation of a domain pattern. This causes a drop in the remnant polarization; a direct demonstration of the important role free charge plays in thin ferroelectric films.

preprint2013arXiv

Flaw-driven Failure in Nanostructures

Understanding failure in nanomaterials is critical for the design of reliable structural materials and small-scale devices that have components or microstructural elements at the nanometer length scale. No consensus exists on the effect of flaws on fracture in bulk nanostructured materials or in nanostructures. Proposed theories include nanoscale flaw tolerance and maintaining macroscopic fracture relationships at the nanoscale with virtually no experimental support. We explore fracture mechanisms in nanomaterials via nanomechanical experiments on nanostructures with pre-fabricated surface flaws in combination with molecular dynamics simulations. Nanocrystalline Pt cylinders with diameters of 120 nm with intentionally introduced surface notches were created using a template-assisted electroplating method and tested in uniaxial tension in in-situ SEM. Experiments demonstrate that 8 out of 12 samples failed at the notches and that tensile failure strengths were 1.8 GPa regardless of whether failure occurred at or away from the flaw. These findings suggest that failure location was sensitive to the presence of flaws, while strength was flaw-insensitive. Molecular dynamics simulations support these observations and show that incipient plastic deformation commences via nucleation and motion of dislocations in concert with grain boundary sliding. We postulate that such local plasticity reduces stress concentration ahead of the flaw to levels comparable with the strengths of intrinsic microstructural features like grain boundary triple junctions, a phenomenon unique to nano-scale solids that contain an internal microstructural energy landscape. This mechanism causes failure to occur at the weakest link, be it an internal inhomogeneity or a surface feature with a high local stress.

preprint2013arXiv

Manipulating Ferroelectric Domains in Nanostructures Under Electron Beams

Freestanding BaTiO3 nanodots exhibit domain structures characterized by distinct quadrants of ferroelastic 90° domains in transmission electron microscopy (TEM) observations. These differ significantly from flux-closure domain patterns in the same systems imaged by piezoresponse force microscopy. Based upon a series of phase field simulations of BaTiO3 nanodots, we suggest that the TEM patterns result from a radial electric field arising from electron beam charging of the nanodot. For sufficiently large charging, this converts flux-closure domain patterns to quadrant patterns with radial net polarizations. Not only does this explain the puzzling patterns that have been observed in TEM studies of ferroelectric nanodots, but also suggests how to manipulate ferroelectric domain patterns via electron beams.

preprint2012arXiv

Continuum Elasticity Theory Approach for Spontaneous Bending and Twisting of Ribbons Induced by Mechanical Anisotropy

Helical ribbons arise in many biological and engineered systems, often driven by anisotropic surface stress, residual strain, and geometric or elastic mismatch between layers of a laminated composite. A full mathematical analysis is developed to analytically predict the equilibrium deformed helical shape of an initially flat, straight ribbon, with prescribed magnitudes and orientations of the principal curvatures when subjected to arbitrary surface stress and/or internal residual strain distribution. The helix angle, radius, axis and chirality of the deformed helical ribbons are predicted with a comprehensive, three-dimensional model that incorporates elasticity, differential geometry, and variational principles. In general, the mechanical anisotropy (e.g., in surface/external stress, residual strain or elastic modulus) will lead to spontaneous, three-dimensional helical deformations. Ring shapes are formed when the principle axes of deformation coincide with the geometric axes of the ribbon. The transition from cylindrical helical ribbons to purely twisted ribbons, or tubular structures is controlled by tuning relevant geometric parameters. Analytic, closed-form predictions of the ribbon shapes are validated with simple, table-top experiments. This theoretical approach represent a tool to inform the design of materials and systems in order to achieve desired helical geometries on demand.

preprint2012arXiv

Nonlinear geometric effects in mechanical bistable morphing structures

Bistable structures associated with non-linear deformation behavior, exemplified by the Venus flytrap and slap bracelet, can switch between different functional shapes upon actuation. Despite numerous efforts in modeling such large deformation behavior of shells, the roles of mechanical and nonlinear geometric effects on bistability remain elusive. We demonstrate, through both theoretical analysis and table-top experiments, that two dimensionless parameters control bistability. Our work classifies the conditions for bistability, and extends the large deformation theory of plates and shells.

preprint2011arXiv

First-principles study of graphene edge properties and flake shapes

We use density functional theory to determine the equilibrium shape of graphene flakes, through the calculation of the edge orientation dependence of the edge energy and edge stress of graphene nanoribbons. The edge energy is a nearly linear function of edge orientation angle; increasing from the armchair orientation to the zigzag orientation. Reconstruction of the zigzag edge lowers its energy to less than that of the armchair edge. The edge stress for all edge orientations is compressive, however, reconstruction of the zigzag edge reduces this edge stress to near zero. Hydrogen adsorption is favorable for all edge orientations; dramatically lowering all edge energies and all edge stresses. It also removes the reconstruction of the zigzag edge. Using the new edge energy data, we determine the equilibrium shape of a graphene sheet (with unreconstructed edges) to be hexagonal with straight armchair edges in the presence and absence of hydrogen. However, zigzag edge reconstruction produces graphene flakes with a six-fold symmetry, but with rounded edges. This shape is dominated by near zigzag edges. The compressive edge stresses will lead to edge buckling (out-of-the-plane of the graphene sheet) for all edge orientations, in the absence of hydrogen. Exposing the graphene flake to hydrogen dramatically decreases the buckling amplitude.

preprint2011arXiv

Molecular dynamics simulation of Ga penetration along Al grain boundaries under a constant strain rate condition

While diverse fracture characteristics have been observed in liquid metal embrittlement (LME) depending on the solid-liquid metal pairs, the penetration of nanometer-thick liquid metal films along the grain boundary has been identified as one of the key mechanisms for embrittlement in many classical LME systems, such as Al-Ga, Cu-Bi and Ni-Bi. For example, liquid Ga quickly penetrates deep into grain boundaries in Al, leading to intergranular fracture under very small stresses. We report on a series of molecular dynamics simulations of liquid Ga in contact with an Al bicrystal under a constant strain rate. We identify the grain boundary dislocations that are nucleated at the grain boundary groove tip and climb down along the grain boundary during Ga penetration and characterize their atomic structures based on topological method.

preprint2008arXiv

Grain Boundary Motion Exhibits the Dynamics of Glass-Forming Liquids

Polycrystalline materials can be viewed as composites of crystalline particles or grains separated from one another by thin amorphous grain boundary (GB) regions. While GB have been exhaustively investigated at low temperatures, where these regions resolve into complex ordered structures accessible to measurement, much less is known about them at higher temperatures where the GB can exhibit significant mobility, structural disorder, and where experimental characterization methods are limited. The time and spatial scales accessible to molecular dynamics (MD) simulation make this method appropriate for investigating both the dynamical and structural properties of grain boundaries at elevated temperatures. In the present study, we use MD simulations to determine how the GB dynamics changes with temperature and applied stress. It has long been hypothesized that GB have features in common with glass-forming liquids based on the processing characteristics of polycrystalline materials. We find remarkable support for this suggestion, as evidenced by string-like collective motion, transient caging of atom motion, and non-Arrhenius (Vogel-Fulcher) temperature dependence of GB mobility. Evidently, the frustration caused by the inability of atoms in GB region to simultaneously order with respect to competing grains is responsible for this striking similarity. The paradigm that grains in a polycrystalline material are encapsulated by a frustrated fluid provides a powerful conceptual model of polycrystalline materials, pointing the way to improved control over their material properties.