Showing 1–21 of 21 results for author: Mason, J K

Texture tomography, a versatile framework to study crystalline texture in 3D

Authors: M. P. K. Frewein, J. K. Mason, B. Maier, H. Cölfen, M. Burghammer, A. A. Medjahed, M. Allain, T. A. Grünewald

Abstract: The crystallographic texture is a key organization feature of many technical and biological materials. In these materials, especially hierarchically structured ones, the preferential alignment of the nano constituents is heavily influencing the macroscopic behaviour of the material. In order to study local crystallographic texture with both high spatial and angular resolution, we developed Texture… ▽ More The crystallographic texture is a key organization feature of many technical and biological materials. In these materials, especially hierarchically structured ones, the preferential alignment of the nano constituents is heavily influencing the macroscopic behaviour of the material. In order to study local crystallographic texture with both high spatial and angular resolution, we developed Texture tomography (TexTOM). This approach allows to model the diffraction data of polycrystalline materials by using the full reciprocal space of the ensemble of crystals and describe the texture in each voxel via a orientation distribution function. This means, it provides 3D reconstructions of the local texture by measuring the probabilities of all crystal orientations. The TexTOM approach addresses limitations associated with existing models: It correlates the intensities from several Bragg reflections, thus reduces ambiguities resulting from symmetry. Further, it yields quantitative probability distributions of local real space crystal orientations without further assumptions on the sample structure. Finally, its efficient mathematical formulation enables reconstructions faster than the time-scale of the experiment. In this manuscript, we present the mathematical model, the inversion strategy and its current experimental implementation. We show characterizations of simulated data as well as experimental data obtained from a synthetic, inorganic model sample, the silica-witherite biomorph. In conclusion, Tex-TOM provides a versatile framework to reconstruct 3D quantitative texture information for polycrystalline samples. In this way, it opens the door for unprecedented insights into the nanostructural makeup of natural and technical materials. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Simulated surface diffusion in nanoporous gold and its dependence on surface curvature

Authors: Conner Marie Winkeljohn, Sadi Md Shahriar, Erkin Seker, Jeremy K. Mason

Abstract: The morphological evolution of nanoporous gold is generally believed to be governed by surface diffusion. This work specifically explores the dependence of mass transport by surface diffusion on the curvature of a gold surface. The surface diffusivity is estimated by molecular dynamics simulations for a variety of surfaces of constant mean curvature, eliminating any chemical potential gradients an… ▽ More The morphological evolution of nanoporous gold is generally believed to be governed by surface diffusion. This work specifically explores the dependence of mass transport by surface diffusion on the curvature of a gold surface. The surface diffusivity is estimated by molecular dynamics simulations for a variety of surfaces of constant mean curvature, eliminating any chemical potential gradients and allowing the possible dependence of the surface diffusivity on mean curvature to be isolated. The apparent surface diffusivity is found to have an activation energy of ~0.74 eV with a weak dependence on curvature, but is consistent with the values reported in the literature. The apparent concentration of mobile surface atoms is found to be highly variable, having an Arrhenius dependence on temperature with an activation energy that also has a weak curvature dependence. These activation energies depend on curvature in such a way that the rate of mass transport by surface diffusion is nearly independent of curvature, but with a higher activation energy of ~1.01 eV. The curvature dependencies of the apparent surface diffusivity and concentration of mobile surface atoms is believed to be related to the expected lifetime of a mobile surface atom, and has the practical consequence that a simulation study that does not account for this finite lifetime could underestimate the activation energy for mass transport via surface diffusion by ~0.27 eV. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 14 pages, 12 figures

Energetic contributions to deformation twinning in magnesium

Authors: Enver Kapan, Sertan Alkan, C. Can Aydıner, Jeremy K. Mason

Abstract: Modeling deformation twin nucleation in magnesium has proven to be a challenging task. In particular, the absence of a heterogeneous twin nucleation model which provides accurate energetic descriptions for twin-related structures belies a need to more deeply understand twin energetics. To address this problem, molecular dynamics simulations are performed to follow the energetic evolution of… ▽ More Modeling deformation twin nucleation in magnesium has proven to be a challenging task. In particular, the absence of a heterogeneous twin nucleation model which provides accurate energetic descriptions for twin-related structures belies a need to more deeply understand twin energetics. To address this problem, molecular dynamics simulations are performed to follow the energetic evolution of $\{10\overline{1}2\}$ tension twin embryos nucleating from an asymmetrically-tilted grain boundary. The line, surface and volumetric terms associated with twin nucleation are identified. A micromechanical model is proposed where the stress field around the twin nucleus is estimated using the Eshelby formalism, and the contributions of the various twin-related structures to the total energy of the twin are evaluated. The reduction in the grain boundary energy arising from the change in character of the prior grain boundary is found to be able to offset the energy costs of the other interfaces. The defect structures bounding the stacking faults that form inside the twin are also found to possibly have significant energetic contributions. These results suggest that both of these effects could be critical considerations when predicting twin nucleation sites in magnesium. △ Less

Submitted 12 July, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: 14 pages, 14 figures

Dependence of simulated radiation damage on crystal structure and atomic misfit in metals

Authors: Jared C. Stimac, Chad Serrao, Jeremy K. Mason

Abstract: This study investigates radiation damage in three metals in the low temperature and high radiant flux regime using molecular dynamics and a Frenkel pair accumulation method to simulate up to $2.0$ displacements per atom. The metals considered include Fe, equiatomic CrCoNi, and a fictitious metal with identical bulk properties to the CrCoNi composed of a single atom type referred to as an A-atom. C… ▽ More This study investigates radiation damage in three metals in the low temperature and high radiant flux regime using molecular dynamics and a Frenkel pair accumulation method to simulate up to $2.0$ displacements per atom. The metals considered include Fe, equiatomic CrCoNi, and a fictitious metal with identical bulk properties to the CrCoNi composed of a single atom type referred to as an A-atom. CrCoNi is found to sustain higher concentrations of dislocations than either the Fe or A-atom systems and more stacking faults than the A-atom system. The results suggest that the concentration of vacancies and interstitials are substantially higher for the CrCoNi than the A-atom system, perhaps reflecting that the recombination radius is smaller in CrCoNi due to the roughened potential energy landscape. A model that partitions the major contributions from defects to the stored energy is described, and serves to highlight a general need for higher fidelity approaches to point defect identification. △ Less

Submitted 10 January, 2023; originally announced January 2023.

Report number: LLNL-JRNL-839942-DRAFT

Energy storage under high-rate compression of single crystal tantalum

Authors: Jared C. Stimac, Nicolas Bertin, Jeremy K. Mason, Vasily V. Bulatov

Abstract: When a material is plastically deformed the majority of mechanical work is dissipated as heat, and the fraction of plastic work converted into heat is known as the Taylor-Quinney coefficient (TQC). Large-scale molecular dynamics simulations were performed of high strain rate compression of single-crystal tantalum, and the resulting integral and differential TQC values are reported up to true strai… ▽ More When a material is plastically deformed the majority of mechanical work is dissipated as heat, and the fraction of plastic work converted into heat is known as the Taylor-Quinney coefficient (TQC). Large-scale molecular dynamics simulations were performed of high strain rate compression of single-crystal tantalum, and the resulting integral and differential TQC values are reported up to true strains of 1.0. A phenomenological model is proposed for the energy stored in the material as a function of time with an asymptotic limit for this energy defined by the deformation conditions. The model reasonably describes the convergence of TQC values to 1.0 with increasing plastic strain, but does not directly address the physical nature of thermo-mechanical conversion. This is instead developed in a second more detailed model that accurately accounts for energy storage in two distinct contributions, one being the growing dislocation network and the other the point defect debris left behind by the moving dislocations. The contribution of the point defect debris is found to lag behind that of the dislocation network but to be substantial under the high-rate straining conditions considered here. △ Less

Submitted 5 August, 2022; v1 submitted 9 May, 2022; originally announced May 2022.

Comments: 28 pages, 12 figures

Report number: LLNL-JRNL-835041-DRAFT

A geometric conjecture about phase transitions

Authors: O. B. Ericok, J. K. Mason

Abstract: As phenomena that necessarily emerge from the collective behavior of interacting particles, phase transitions continue to be difficult to predict using statistical thermodynamics. A recent proposal called the topological hypothesis suggests that the existence of a phase transition could perhaps be inferred from changes to the topology of the accessible part of the configuration space. This paper i… ▽ More As phenomena that necessarily emerge from the collective behavior of interacting particles, phase transitions continue to be difficult to predict using statistical thermodynamics. A recent proposal called the topological hypothesis suggests that the existence of a phase transition could perhaps be inferred from changes to the topology of the accessible part of the configuration space. This paper instead suggests that such a topological change is often associated with a dramatic change in the configuration space geometry, and that the geometric change is the actual driver of the phase transition. More precisely, a geometric change that brings about a discontinuity in the mixing time required for an initial probability distribution on the configuration space to reach steady-state is conjectured to be related to the onset of a phase transition in the thermodynamic limit. This conjecture is tested by evaluating the diffusion diameter and $ε$-mixing time of the configuration spaces of hard disk and hard sphere systems of increasing size. Explicit geometries are constructed for the configuration spaces of these systems, and numerical evidence suggests that a discontinuity in the $ε$-mixing time coincides with the solid-fluid phase transition in the thermodynamic limit. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 12 pages, 13 figures

Journal ref: Geometric conjecture about phase transitions. Physical Review E, 107(6) (2023)

Foundations of a Finite Non-Equilibrium Statistical Thermodynamics: Extrinsic Quantities

Authors: O. B. Ericok, J. K. Mason

Abstract: Statistical thermodynamics is valuable as a conceptual structure that shapes our thinking about equilibrium thermodynamic states. A cloud of unresolved questions surrounding the foundations of the theory could lead an impartial observer to conclude that statistical thermodynamics is in a state of crisis though. Indeed, the discussion about the microscopic origins of irreversibility has continued i… ▽ More Statistical thermodynamics is valuable as a conceptual structure that shapes our thinking about equilibrium thermodynamic states. A cloud of unresolved questions surrounding the foundations of the theory could lead an impartial observer to conclude that statistical thermodynamics is in a state of crisis though. Indeed, the discussion about the microscopic origins of irreversibility has continued in the scientific community for more than a hundred years. This paper considers these questions while beginning to develop a statistical thermodynamics for finite non-equilibrium systems. Definitions are proposed for all of the extrinsic variables of the fundamental thermodynamic relation that are consistent with existing results in the equilibrium thermodynamic limit. The probability density function on the phase space is interpreted as a subjective uncertainty about the microstate, and the Gibbs entropy formula is modified to allow for entropy creation without introducing additional physics or modifying the phase space dynamics. Resolutions are proposed to the mixing paradox, Gibbs' paradox, Loschmidt's paradox, and Maxwell's demon thought experiment. Finally, the extrinsic variables of the fundamental thermodynamic relation are evaluated as functions of time and space for a diffusing ideal gas, and the initial and final values are shown to coincide with the expected equilibrium values when interpreted in a classical context. △ Less

Submitted 8 January, 2022; originally announced January 2022.

Comments: 18 pages, 10 figures

A constant of motion for ideal grain growth in three dimensions

Authors: Erdem Eren, Jeremy Kyle Mason

Abstract: Most metallic and ceramic materials are comprised of a space-filling collection of crystalline grains separated by grain boundaries. While this grain structure has been studied for more than a century, there few rigorous results regarding its global properties available in the literature. We present a new, rigorous result for three-dimensional grain structures that relates the integral of the Gaus… ▽ More Most metallic and ceramic materials are comprised of a space-filling collection of crystalline grains separated by grain boundaries. While this grain structure has been studied for more than a century, there few rigorous results regarding its global properties available in the literature. We present a new, rigorous result for three-dimensional grain structures that relates the integral of the Gaussian curvature over the grain boundaries to the numbers of grains and quadruple junctions. The result is numerically verified for a grain structure consisting of periodic truncated octahedra. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: 8 pages, 5 figures

Journal ref: Phys. Rev. B 104, L140103 (2021)

Configuration spaces of hard spheres

Authors: O. B. Ericok, K. Ganesan, J. K. Mason

Abstract: Hard sphere systems are often used to model simple fluids. The configuration spaces of hard spheres in a three-dimensional torus modulo various symmetry groups are comparatively simple, and could provide valuable information about the nature of phase transitions. Specifically, the topological changes in the configuration space as a function of packing fraction have been conjectured to be related t… ▽ More Hard sphere systems are often used to model simple fluids. The configuration spaces of hard spheres in a three-dimensional torus modulo various symmetry groups are comparatively simple, and could provide valuable information about the nature of phase transitions. Specifically, the topological changes in the configuration space as a function of packing fraction have been conjectured to be related to the onset of first-order phase transitions. The critical configurations for one to twelve spheres are sampled using a Morse-theoretic approach, and are available in an online, interactive database. Explicit triangulations are constructed for the configuration spaces of the two sphere system, and their topological and geometric properties are studied. The critical configurations are found to be associated with geometric changes to the configuration space that connect previously distant regions and reduce the configuration space diameter as measured by the commute time and diffusion distances. The number of such critical configurations around the packing fraction of the solid-liquid phase transition increases exponentially with the number of spheres, suggesting that the onset of the first-order phase transition in the thermodynamic limit is associated with a discontinuity in the configuration space diameter. △ Less

Submitted 28 October, 2021; v1 submitted 18 June, 2021; originally announced June 2021.

Comments: 15 pages, 16 figures

Journal ref: Phys. Rev. E 104 (2021) 055304

Topological transitions during grain growth on a finite element mesh

Authors: Erdem Eren, Jeremy K. Mason

Abstract: The topological transitions that occur to the grain boundary network during grain growth in a material with uniform grain boundary energies are believed to be known. The same is not true for more realistic materials, since more general grain boundary energies in principle allow many more viable grain boundary configurations. A simulation of grain growth in such a material therefore requires a proc… ▽ More The topological transitions that occur to the grain boundary network during grain growth in a material with uniform grain boundary energies are believed to be known. The same is not true for more realistic materials, since more general grain boundary energies in principle allow many more viable grain boundary configurations. A simulation of grain growth in such a material therefore requires a procedure to enumerate all possible topological transitions and select the most energetically favorable one. Such a procedure is developed and implemented here for a microstructure represented by a volumetric finite element mesh. As a specific example, all possible transitions for a typical configuration with five grains around a junction point are enumerated, and some exceptional transitions are found to be energetically similar to the conventional ones even for a uniform boundary energy. A general discrete formulation to calculate grain boundary velocities is used to simulate grain growth for an example microstructure. The method is implemented as a C++ library based on SCOREC, an open source massively parallelizable library for finite element simulations with adaptive meshing. △ Less

Submitted 28 January, 2021; originally announced January 2021.

Comments: 19 pages(13 main, 6 supplement), 25 figures (19 main, 6 supplement)

Journal ref: Phys. Rev. Materials 5, 103802 (2021)

Quotient Maps and Configuration Spaces of Hard Disks

Authors: O. B. Ericok, J. K. Mason

Abstract: Hard disks systems are often considered as prototypes for simple fluids. In a statistical mechanics context, the hard disk configuration space is generally quotiented by the action of various symmetry groups. The changes in the topological and geometric properties of the configuration spaces effected by such quotient maps are studied for small numbers of disks on a square and hexagonal torus. A me… ▽ More Hard disks systems are often considered as prototypes for simple fluids. In a statistical mechanics context, the hard disk configuration space is generally quotiented by the action of various symmetry groups. The changes in the topological and geometric properties of the configuration spaces effected by such quotient maps are studied for small numbers of disks on a square and hexagonal torus. A metric is defined on the configuration space and the various quotient spaces that respects the desired symmetries. This is used to construct explicit triangulations of the configuration spaces as $α$-complexes. Critical points in a configuration space are associated with changes in the topology as a function of disk radius, are conjectured to be related to the configurational entropy of glassy systems, and could reveal the origins of phase transitions in other systems. The number and topological and geometric properties of the critical points are found to depend on the symmetries by which the configuration space is quotiented. △ Less

Submitted 29 May, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Journal ref: Granular Matter 24, 76 (2022)

Classification of atomic environments via the Gromov-Wasserstein distance

Authors: Sakura Kawano, Jeremy K. Mason

Abstract: Interpreting molecular dynamics simulations usually involves automated classification of local atomic environments to identify regions of interest. Existing approaches are generally limited to a small number of reference structures and only include limited information about the local chemical composition. This work proposes to use a variant of the Gromov-Wasserstein (GW) distance to quantify the d… ▽ More Interpreting molecular dynamics simulations usually involves automated classification of local atomic environments to identify regions of interest. Existing approaches are generally limited to a small number of reference structures and only include limited information about the local chemical composition. This work proposes to use a variant of the Gromov-Wasserstein (GW) distance to quantify the difference between a local atomic environment and a set of arbitrary reference environments in a way that is sensitive to atomic displacements, missing atoms, and differences in chemical composition. This involves describing a local atomic environment as a finite metric measure space, which has the additional advantages of not requiring the local environment to be centered on an atom and of not making any assumptions about the material class. Numerical examples illustrate the efficacy and versatility of the algorithm. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Journal ref: Comp. Mater. Sci. 188, 110144 (2021)

Distribution of Topological Types in Grain-Growth Microstructures

Authors: Emanuel A. Lazar, Jeremy K. Mason, Robert D. MacPherson, David J. Srolovitz

Abstract: 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 o… ▽ More 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. △ Less

Submitted 4 July, 2020; originally announced July 2020.

Comments: 6 pages, 5 figures

Journal ref: Phys. Rev. Lett. 125, 015501 (2020)

Basis functions on the grain boundary space: Theory

Authors: Jeremy K. Mason, Srikanth Patala

Abstract: With the increasing availability of experimental and computational data concerning the properties and distribution of grain boundaries in polycrystalline materials, there is a corresponding need to efficiently and systematically express functions on the grain boundary space. A grain boundary can be described by the rotations applied to two grains on either side of a fixed boundary plane, suggestin… ▽ More With the increasing availability of experimental and computational data concerning the properties and distribution of grain boundaries in polycrystalline materials, there is a corresponding need to efficiently and systematically express functions on the grain boundary space. A grain boundary can be described by the rotations applied to two grains on either side of a fixed boundary plane, suggesting that the grain boundary space is related to the space of rotations. This observation is used to construct an orthornormal function basis, allowing effectively arbitrary functions on the grain boundary space to be written as linear combinations of the basis functions. Moreover, a procedure is developed to construct a smaller set of basis functions consistent with the crystallographic point group symmetries, grain exchange symmetry, and the null boundary singularity. Functions with the corresponding symmetries can be efficiently expressed as linear combinations of the symmetrized basis functions. An example is provided that shows the efficacy of the symmetrization procedure. △ Less

Submitted 25 September, 2019; originally announced September 2019.

Continuous and Optimally Complete Description of Chemical Environments Using Spherical Bessel Descriptors

Authors: Emir Kocer, Jeremy K. Mason, Hakan Erturk

Abstract: Recently, machine learning potentials have been advanced as candidates to combine the high-accuracy of quantum mechanical simulations with the speed of classical interatomic potentials. A crucial component of a machine learning potential is the description of local atomic environments by some set of descriptors. These should ideally be continuous throughout the specified local atomic environment,… ▽ More Recently, machine learning potentials have been advanced as candidates to combine the high-accuracy of quantum mechanical simulations with the speed of classical interatomic potentials. A crucial component of a machine learning potential is the description of local atomic environments by some set of descriptors. These should ideally be continuous throughout the specified local atomic environment, twice-differentiable with respect to atomic positions and complete in the sense of containing all possible information about the neighborhood. An updated version of the recently proposed Spherical Bessel descriptors satisfies all three of these properties, and moreover is optimally complete in the sense of encoding all configurational information with the smallest possible number of descriptors. The Smooth Overlap of Atomic Position descriptors that are frequently visited in the literature and the Zernike descriptors that are built upon a similar basis are included into the discussion as being the natural counterparts of the Spherical Bessel descriptors, and shown to be incapable of satisfying the full list of core requirements for an accurate description. Aside being mathematically and physically superior, the Spherical Bessel descriptors have also the advantage of allowing machine learning potentials of comparable accuracy that require roughly an order of magnitude less computation time per evaluation than the Smooth Overlap of Atomic Position descriptors, which appear to be the common choice of descriptors in recent studies. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Comments: 15 pages, 5 figures, under review for Journal of Chemical Physics

A Novel Approach to Describe Chemical Environments in High Dimensional Neural Network Potentials

Authors: Emir Kocer, Jeremy K. Mason, Hakan Erturk

Abstract: A central concern of molecular dynamics simulations are the potential energy surfaces that govern atomic interactions. These hypersurfaces define the potential energy of the system, and have generally been calculated using either predefined analytical formulas (classical) or quantum mechanical simulations (ab initio). The former can accurately reproduce only a selection of material properties, whe… ▽ More A central concern of molecular dynamics simulations are the potential energy surfaces that govern atomic interactions. These hypersurfaces define the potential energy of the system, and have generally been calculated using either predefined analytical formulas (classical) or quantum mechanical simulations (ab initio). The former can accurately reproduce only a selection of material properties, whereas the latter is restricted to short simulation times and small systems. Machine learning potentials have recently emerged as a third approach to model atomic interactions, and are purported to offer the accuracy of ab initio simulations with the speed of classical potentials. However, the performance of machine learning potentials depends crucially on the description of a local atomic environment. A set of invariant, orthogonal and differentiable descriptors for an atomic environment is proposed, implemented in a neural network potential for solid-state silicon, and tested in molecular dynamics simulations. Neural networks using the proposed descriptors are found to outperform ones using the Behler Parinello and SOAP descriptors currently in the literature. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Comments: 23 Pages, 5 figures, 2 tables, journal article

Journal ref: Journal of Chemical Physics, vol.150, 154102, 2019

Improved conditioning of the Floater--Hormann interpolants

Authors: Jeremy K. Mason

Abstract: The Floater--Hormann family of rational interpolants do not have spurious poles or unattainable points, are efficient to calculate, and have arbitrarily high approximation orders. One concern when using them is that the amplification of rounding errors increases with approximation order, and can make balancing the interpolation error and rounding error difficult. This article proposes to modify th… ▽ More The Floater--Hormann family of rational interpolants do not have spurious poles or unattainable points, are efficient to calculate, and have arbitrarily high approximation orders. One concern when using them is that the amplification of rounding errors increases with approximation order, and can make balancing the interpolation error and rounding error difficult. This article proposes to modify the Floater--Hormann interpolants by including additional local polynomial interpolants at the ends of the interval. This appears to improve the conditioning of the interpolants and allow higher approximation orders to be used in practice. △ Less

Submitted 23 June, 2017; originally announced June 2017.

Comments: 13 pages, 4 figures, 1 table

MSC Class: 65D05; 41A05; 41A20

Roundness of grains in cellular microstructures

Authors: Frank H. Lutz, Jeremy K. Mason, Emanuel A. Lazar, Robert D. MacPherson

Abstract: Many physical systems are composed of polyhedral cells of varying sizes and shapes. These structures are simple in the sense that no more than three faces meet at an edge and no more than four edges meet at a vertex. This means that individual cells can usually be considered as simple, three-dimensional polyhedra. This paper is concerned with determining the distribution of combinatorial types of… ▽ More Many physical systems are composed of polyhedral cells of varying sizes and shapes. These structures are simple in the sense that no more than three faces meet at an edge and no more than four edges meet at a vertex. This means that individual cells can usually be considered as simple, three-dimensional polyhedra. This paper is concerned with determining the distribution of combinatorial types of such polyhedral cells. We introduce the terms \emph{fundamental} and \emph{vertex-truncated} types and apply these concepts to the grain growth microstructure as a testing ground. For these microstructures we demonstrate that most grains are of particular fundamental types, whereas the frequency of vertex-truncated types decreases exponentially with the number of truncations. This can be explained by the evolutionary process through which grain growth structures are formed, and in which energetically unfavorable surfaces are quickly eliminated. Furthermore, we observe that these grain types are `round' in a combinatorial sense: there are no `short' separating cycles that partition the polyhedra into two parts of similar sizes. A particular microstructure derived from the Poisson--Voronoi initial condition is identified as containing an unusually large proportion of round grains. This Round microstructure has an average of $14.036$ faces per grain, and is conjectured to be more resistant to topological change than the steady-state grain growth microstructure. △ Less

Submitted 8 July, 2017; v1 submitted 30 September, 2016; originally announced September 2016.

Comments: 26 pages, 33 figures, 11 tables

Journal ref: Phys. Rev. E 96, 023001 (2017)

Geometric and topological properties of the canonical grain growth microstructure

Authors: Jeremy K. Mason, Emanuel A. Lazar, Robert D. MacPherson, David J. Srolovitz

Abstract: 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. Th… ▽ More 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. △ Less

Submitted 13 July, 2015; originally announced July 2015.

Comments: 19 pages, 21 figures

Journal ref: Phys. Rev. E 92, 063308 (2015)

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

Authors: Emanuel A. Lazar, Jeremy K. Mason, Robert D. MacPherson, David J. Srolovitz

Abstract: 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 measure… ▽ More 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. △ Less

Submitted 8 January, 2014; originally announced January 2014.

Comments: 15 pages, 19 figures, Supplemental Material included

Journal ref: Phys. Rev. E 88, 063309 (2013)

Complete topology of cells, grains, and bubbles in three-dimensional microstructures

Authors: Emanuel A. Lazar, Jeremy K. Mason, Robert D. MacPherson, David J. Srolovitz

Abstract: We introduce a general, efficient method to completely describe the topology of individual grains, bubbles, and cells in three-dimensional polycrystals, foams, and other multicellular microstructures. This approach is applied to a pair of three-dimensional microstructures that are often regarded as close analogues in the literature: one resulting from normal grain growth (mean curvature flow) and… ▽ More We introduce a general, efficient method to completely describe the topology of individual grains, bubbles, and cells in three-dimensional polycrystals, foams, and other multicellular microstructures. This approach is applied to a pair of three-dimensional microstructures that are often regarded as close analogues in the literature: one resulting from normal grain growth (mean curvature flow) and another resulting from a random Poisson-Voronoi tessellation of space. Grain growth strongly favors particular grain topologies, compared with the Poisson-Voronoi model. Moreover, the frequencies of highly symmetric grains are orders of magnitude higher in the the grain growth microstructure than they are in the Poisson-Voronoi one. Grain topology statistics provide a strong, robust differentiator of different cellular microstructures and provide hints to the processes that drive different classes of microstructure evolution. △ Less

Submitted 20 July, 2012; originally announced July 2012.

Comments: 5 pages, 6 figures, 5 supplementary pages

Journal ref: Phys. Rev. Lett. 109, 095505 (2012)