The effect of differential mineral shrinkage on crack formation and network geometry.

Rock, concrete, and other engineered materials are often composed of several minerals that change volumetrically in response to variations in the moisture content of the local environment. Such differential shrinkage is caused by varying shrinkage rates between mineral compositions during dehydration. Using both 3D X-ray imaging of geo-architected samples and peridynamic (PD) numerical simulations, we show that the spatial distribution of the clay affects the crack network geometry with distributed clay particles yielding the most complex crack networks and percent damage (99.56%), along with a 60% reduction in material strength. We also demonstrate that crack formation, growth, coalescence, and distribution during dehydration, are controlled by the differential shrinkage rates between a highly shrinkable clay and a homogeneous mortar matrix. Sensitivity tests performed with the PD models show a clay shrinkage parameter of 0.4 yields considerable damage, and reductions in the parameter can result in a significant reduction in fracturing and an increase in material strength. Additionally, isolated clay inclusions induced localized fracturing predominantly due to debonding between the clay and matrix. These insights indicate differential shrinkage is a source of potential failure in natural and engineered barriers used to sequester anthropogenic waste.

Correlating structure and transport behavior in Li+ and O2 containing pyrrolidinium ionic liquids.

Ionic liquids are a unique class of materials with several potential applications in electrochemical energy storage. When used in electrolytes, these highly coordinating solvents can influence device performance through their high viscosities and strong solvation behaviors. In this work, we explore the effects of pyrrolidinium cation structure and Li+ concentration on transport processes in ionic liquid electrolytes. We present correlated experimental measurements and molecular simulations of Li+ mobility and O2 diffusivity, and connect these results to dynamic molecular structural information and device performance. In the context of Li-O2/Li-air battery chemistries, we find that Li+ mobility is largely influenced by Li+-anion coordination, but that both Li+ and O2 diffusion may be affected by variations of the pyrrolidinium cation and Li+ concentration.

Two-dimensional metal-organic frameworks with high thermoelectric efficiency through metal ion selection.

Two-dimensional (2D) materials have attracted much attention due to their novel properties. An exciting new class of 2D materials based on metal-organic frameworks (MOFs) has recently emerged, displaying high electrical conductivity, a rarity among organic nanoporous materials. The emergence of these materials raises intriguing questions about their fundamental electronic, optical, and thermal properties, but few studies exist in this regard. Here we present an atomistic study of the thermoelectric properties of crystalline 2D MOFs X3(HITP)2 with X = Ni, Pd or Pt, and HITP = 2,3,6,7,10,11-hexaiminotriphenylene, using both ab initio transport models and classical molecular dynamics simulations. We find that these materials have a high Seebeck coefficient and low thermal conductivity, making them promising for thermoelectric applications. Furthermore, we explore the dependence of thermoelectric transport properties on the atomic structure by comparing the calculated band structure, band alignment, and electronic density of states of the three 2D MOFs, and find that the thermoelectric transport properties strongly depend on both the interaction between the ligands and the metal ions, and the d orbital splitting of the metal ions induced by the ligands. This demonstrates that selection of the metal ion is a powerful approach to control and enhance the thermoelectric properties. Interestingly we reveal an unexpected effect where, unlike for electrons, the thermal and electrical current may not be equally carried by the holes, leading to a significant deviation from the Wiedemann-Franz law. The results of this work provide fundamental guidance to optimize the existing 2D MOFs, and to design and discover new families of MOF-like materials for thermoelectric applications.

Thermal boundary conductance between Al films and GaN nanowires investigated with molecular dynamics.

GaN nanowires are being pursued for optoelectronic and high-power applications. In either use, increases in operating temperature reduce both performance and reliability making it imperative to minimize thermal resistances. Since interfaces significantly influence the thermal response of nanosystems, the thermal boundary resistance between GaN nanowires and metal contacts has major significance. In response, we have performed systematic molecular dynamics simulations to study the thermal boundary conductance between GaN nanowires and Al films as a function of nanowire dimensions, packing density, and the depth the nanowire is embedded into the metal contact. At low packing densities, the apparent Kapitza conductance between GaN nanowires and an aluminum film is shown to be larger than when contact is made between films of these same materials. This enhancement decreases toward the film-film limit, however, as the packing density increases. For densely packed nanowires, maximizing the Kapitza conductance can be achieved by embedding the nanowires into the films, as the conductance is found to be proportional to the total contact area.

The application of an atomistic J-integral to a ductile crack.

In this work we apply a Lagrangian kernel-based estimator of continuum fields to atomic data to estimate the J-integral for the emission dislocations from a crack tip. Face-centered cubic (fcc) gold and body-centered cubic (bcc) iron modeled with embedded atom method (EAM) potentials are used as example systems. The results of a single crack with a K-loading compare well to an analytical solution from anisotropic linear elastic fracture mechanics. We also discovered that in the post-emission of dislocations from the crack tip there is a loop size-dependent contribution to the J-integral. For a system with a finite width crack loaded in simple tension, the finite size effects for the systems that were feasible to compute prevented precise agreement with theory. However, our results indicate that there is a trend towards convergence.

Adaptive Green-Kubo estimates of transport coefficients from molecular dynamics based on robust error analysis.

We present a rigorous Green-Kubo methodology for calculating transport coefficients based on on-the-fly estimates of: (a) statistical stationarity of the relevant process, and (b) error in the resulting coefficient. The methodology uses time samples efficiently across an ensemble of parallel replicas to yield accurate estimates, which is particularly useful for estimating the thermal conductivity of semi-conductors near their Debye temperatures where the characteristic decay times of the heat flux correlation functions are large. Employing and extending the error analysis of Zwanzig and Ailawadi [Phys. Rev. 182, 280 (1969)] and Frenkel [in Proceedings of the International School of Physics "Enrico Fermi", Course LXXV (North-Holland Publishing Company, Amsterdam, 1980)] to the integral of correlation, we are able to provide tight theoretical bounds for the error in the estimate of the transport coefficient. To demonstrate the performance of the method, four test cases of increasing computational cost and complexity are presented: the viscosity of Ar and water, and the thermal conductivity of Si and GaN. In addition to producing accurate estimates of the transport coefficients for these materials, this work demonstrates precise agreement of the computed variances in the estimates of the correlation and the transport coefficient with the extended theory based on the assumption that fluctuations follow a Gaussian process. The proposed algorithm in conjunction with the extended theory enables the calculation of transport coefficients with the Green-Kubo method accurately and efficiently.

A Long-Range Electric Field Solver for Molecular Dynamics Based on Atomistic-to-Continuum Modeling.

Understanding charge transport processes at a molecular level is currently hindered by a lack of appropriate models for incorporating nonperiodic, anisotropic electric fields in molecular dynamics (MD) simulations. In this work, we develop a model for including electric fields in MD using an atomistic-to-continuum framework. This framework provides the mathematical and the algorithmic infrastructure to couple finite element (FE) representations of continuous data with atomic data. Our model represents the electric potential on a FE mesh satisfying a Poisson equation with source terms determined by the distribution of the atomic charges. Boundary conditions can be imposed naturally using the FE description of the potential, which then propagate to each atom through modified forces. The method is verified using simulations where analytical solutions are known or comparisons can be made to existing techniques. In addition, a calculation of a salt water solution in a silicon nanochannel is performed to demonstrate the method in a target scientific application in which ions are attracted to charged surfaces in the presence of electric fields and interfering media.

A homogeneous nonequilibrium molecular dynamics method for calculating the heat transport coefficient of mixtures and alloys.

This work generalizes Evans' homogeneous nonequilibrium method for estimating heat transport coefficient to multispecies molecular systems described by general multibody potentials. The proposed method, in addition to being compatible with periodic boundary conditions, is shown to satisfy all the requirements of Evans' original method, namely, adiabatic incompressibility of phase space, equivalence of the dissipative and heat fluxes, and momentum preservation. The difference between the new equations of motion, suitable for mixtures and alloys, and those of Evans' original work are quantified by means of simulations for fluid Ar-Kr and solid GaN test systems.

Generalization of the homogeneous nonequilibrium molecular dynamics method for calculating thermal conductivity to multibody potentials.

This work provides a generalization of Evans' homogeneous nonequilibrium method for estimating thermal conductivity to molecular systems that are described by general multibody potentials. A perturbed form of the usual Nose-Hoover equations of motion is formally constructed and is shown to satisfy the requirements of Evans' original method. These include adiabatic incompressibility of phase space, equivalence of the dissipative and heat fluxes, and momentum preservation.

A homogeneous nonequilibrium molecular dynamics method for calculating thermal conductivity with a three-body potential.

In this work, Evans' homogeneous nonequilibrium molecular dynamics method for estimating thermal conductivity is extended to systems employing three-body potentials. This extension is put on a firm theoretical basis and applied to a silicon lattice modeled by the Stillinger-Weber potential. Two new methods are suggested for estimating the thermal conductivity based on a range of values of the fictitious force. Also, kinetic theory is used to estimate the linear range of the fictitious force necessary to bias the heat flow, thereby potentially reducing the number of simulations needed to estimate thermal conductivity.

Full-field deformation of bovine cornea under constrained inflation conditions.

The viscoelastic response of bovine corneas was characterized using in vitro inflation (bulge) experiments combined with spatially-resolved deformation mapping via digital image correlation. A complex fixture conforming to the limbal annulus was developed to hold the attached sclera rigid while allowing deformation only in the cornea. A statistical set of experiments was performed for a pressure range of 3.6-8 kPa (27-60 mmHg), representing nominal bovine intraocular pressure (IOP) to acute glaucoma conditions. A broader pressure range of 0-32 kPa (0-240 mmHg) was also examined to characterize the nonlinear finite deformation behavior of the tissue. Results showed that for pressures near and above IOP, the majority of the deformation was localized in the limbus and peripheral regions, which left the central cornea largely undeformed. This observation was consistent with the known preferred circumferential alignment of collagen fibrils outside of the central cornea. In general, the inflation experiments observed viscoelastic behavior in the form of rate-dependent hysteresis in the pressure-deformation response of the apex of the cornea, creep in the apex deformation at a constant inflation pressure, and relaxation in the pressure response at a constant inflation volume. The 3.6-8 kPa (27-60 mmHg) pressure range produced small viscoelastic deformations and a nearly linear pressure-deformation response, which suggests that for physiological pressure ranges, the cornea can be approximated as a linear viscoelastic or linear pseudo-elastic material.