Understanding phase transformations in 2D materials can unlock unprecedented developments in nanotechnology, since their unique properties can be dramatically modified by external fields that control the phase change. Here, experiments and simulations are used to investigate the mechanical properties of a 2D diamond boron nitride (BN) phase induced by applying local pressure on atomically thin h-BN on a SiO2 substrate, at room temperature, and without chemical functionalization. Molecular dynamics (MD) simulations show a metastable local rearrangement of the h-BN atoms into diamond crystal clusters when increasing the indentation pressure. Raman spectroscopy experiments confirm the presence of a pressure-induced cubic BN phase, and its metastability upon release of pressure. Å-indentation experiments and simulations show that at pressures of 2–4 GPa, the indentation stiffness of monolayer h-BN on SiO2 is the same of bare SiO2, whereas for two- and three-layer-thick h-BN on SiO2 the stiffness increases of up to 50% compared to bare SiO2, and then it decreases when increasing the number of layers. Up to 4 GPa, the reduced strain in the layers closer to the substrate decreases the probability of the sp2-to-sp3 phase transition, explaining the lower stiffness observed in thicker h-BN.
The accumulation of point defects and defect clusters in materials, as seen in irradiated metals for example, can lead to the formation and growth of voids. Void nucleation is derived from the condensation of supersaturated vacancies and depends strongly on the stress state. It is usually assumed that such stress states can be produced by microstructural defects such dislocations, grain boundaries or triple junctions, however, much less attention has been brought to the formation of voids near microcracks. In this paper, we investigate the coupling between point-defect diffusion/recombination and concentrated stress fields near mode-I crack tips via a spatially-resolved rate theory approach. A modified chemical potential enables point-defect diffusion to be partially driven by the mechanical fields in the vicinity of the crack tip. Simulations are carried out for microcracks using the Griffith model with increasing stress intensity factor K1. Our results show that below a threshold for the stress intensity factor, the microcrack acts purely as a microstructural sink, absorbing point defects. Above this threshold, vacancies accumulate at the crack tip. These results suggest that, even in the absence of plastic deformation, voids can form in the vicinity of a microcrack for a given load when the crack’s characteristic length is above a critical length. While in ductile metals, irradiation damage generally causes hardening and corresponding quasi-brittle cleavage, our results show that irradiation conditions can favor void formation near microstructural stressors such as crack tips leading to lower resistance to crack propagation as predicted by traditional failure analysis.
This report details the current benchmark results to verify, validate and demonstrate the capabilities of the in-house multi-physics phase-field modeling framework Mesoscale Multiphysics Phase Field Simulator (MEMPHIS) developed at the Center for Integrated Nanotechnologies (CINT). MEMPHIS is a general phase-field capability to model various nanoscience and materials science phenomena related to microstructure evolution. MEMPHIS has been benchmarked against a suite of reported classical phase-field benchmark problems to verify and validate the correctness, accuracy and precision of the models and numerical methods currently implemented into the code.
Critical components, such as detonators, in Sandia's stockpile contain heterogeneous materials whose performance and reliability depend on accurate, predictive models of coupled, complex phenomena to predict their synthesis, processing, and operation. Ongoing research in energetic materials has shown that microstructural properties, such as density, pore-size, morphology, and specific surface area are correlated to their initiation threshold and detonation behavior. However, experiments to study these specific characteristics of energetic materials are challenging and time consuming. Therefore, in this work, we turn to mesoscale modeling methods that may be capable of reproducing some observed phenomena to refine and predict outcomes beforehand. Even so, we have no physics-based modeling capability to predict how the microstructure of an energetic material will evolve over near- and long-term time scales. Thus, the goal of this work is to (i) identify any knowledge gaps in how the underlying microstructure forms and evolves during the synthesis process, and (ii) develop and test a mesoscale phase-field model for vapor deposition to capture critical mechanisms of microstructure formation, evolution, and variability in vapor-deposited energetic materials, such as processing conditions, material properties, and substrate interactions. Based on this work, the phase-field method is shown to be a valuable tool for developing the necessary models containing coupled, complex phenomena to investigate and understand the synthesis and processing of energetic materials.
While lattice metamaterials can achieve exceptional energy absorption by tailoring periodically distributed heterogeneous unit cells, relatively little focus has been placed on engineering heterogeneity above the unit-cell level. In this work, the energy-absorption performance of lattice metamaterials with a heterogeneous spatial layout of different unit cell architectures was studied. Such multi-morphology lattices can harness the distinct mechanical properties of different unit cells while being composed out of a single base material. A rational design approach was developed to explore the design space of these lattices, inspiring a non-intuitive design which was evaluated alongside designs based on mixture rules. Fabrication was carried out using two different base materials: 316L stainless steel and Vero White photopolymer. Results show that multi-morphology lattices can be used to achieve higher specific energy absorption than homogeneous lattice metamaterials. Additionally, it is shown that a rational design approach can inspire multi-morphology lattices which exceed rule-of-mixtures expectations.
The self-interstitial atom (SIA) is one of two fundamental point defects in bulk Si. Isolated Si SIAs are extremely difficult to observe experimentally. Even at very low temperatures, they anneal before typical experiments can be performed. Given the challenges associated with experimental characterization, accurate theoretical calculations provide valuable information necessary to elucidate the properties of these defects. Previous studies have applied Kohn-Sham density functional theory (DFT) to the Si SIA, using either the local density approximation or the generalized gradient approximation to the exchange-correlation (XC) energy. The consensus of these studies indicates that a Si SIA may exist in five charge states ranging from -2 to +2 with the defect structure being dependent on the charge state. This study aims to re-examine the existence of these charge states in light of recently derived "approximate bounds"on the defect levels obtained from finite-size supercell calculations and new DFT calculations using both semi-local and hybrid XC approximations. We conclude that only the neutral and +2 charge states are directly supported by DFT as localized charge states of the Si SIA. Within the current accuracy of DFT, our results indicate that the +1 charge state likely consists of an electron in a conduction-band-like state that is coulombically bound to a +2 SIA. Furthermore, the -1 and -2 charge states likely consist of a neutral SIA with one and two additional electrons in the conduction band, respectively.
The interaction of energetic ions with the electronic and ionic system of target materials is an interesting but challenging multiscale problem, and understanding of the early stages after impact of heavy, initially charged ions is particularly poor. At the same time, energy deposition during these early stages determines later formation of damage cascades. We address the multiscale character by combining real-time time-dependent density functional theory for electron dynamics with molecular dynamics simulations of damage cascades. Our first-principles simulations prove that core electrons affect electronic stopping and have an unexpected influence on the charge state of the projectile. We show that this effect is absent for light projectiles, but dominates the stopping physics for heavy projectiles. By parametrizing an inelastic energy loss friction term in the molecular dynamics simulations using our first-principles results, we also show a qualitative influence of electronic stopping physics on radiation-damage cascades.
This report details the current benchmark results to verify, validate and demonstrate the capabilities of the in-house multi-physics phase-field modeling framework Mesoscale Multiphysics Phase Field Simulator (MEMPHIS) developed at the Center for Integrated Nanotechnologies (CINT). MEMPHIS is a general phase-field capability to model various nanoscience and materials science phenomena related to microstructure evolution. MEMPHIS has been benchmarked against a suite of reported 'classical' phase-field benchmark problems to verify and validate the correctness, accuracy and precision of the models and numerical methods currently implemented into the code.
Atomistic modeling of radiation damage through displacement cascades is deceptively non-trivial. Due to the high energy and stochastic nature of atomic collisions, individual primary knock-on atom (PKA) cascade simulations are computationally expensive and ill-suited for length and dose upscaling. Here, we propose a reduced-order atomistic cascade model capable of predicting and replicating radiation events in metals across a wide range of recoil energies. Our methodology approximates cascade and displacement damage production by modeling the cascade as a core-shell atomic structure composed of two damage production estimators, namely an athermal recombination corrected displacements per atom (arc-dpa) in the shell and a replacements per atom (rpa) representing atomic mixing in the core. These estimators are calibrated from explicit PKA simulations and a standard displacement damage model that incorporates cascade defect production efficiency and mixing effects. We illustrate the predictability and accuracy of our reduced-order atomistic cascade method for the cases of copper and niobium by comparing its results with those from full PKA simulations in terms of defect production as well as the resulting cascade evolution and structure. We provide examples for simulating high energy cascade fragmentation and large dose ion-bombardment to demonstrate its possible applicability. Finally, we discuss the various practical considerations and challenges associated with this methodology especially when simulating subcascade formation and dose effects.
Since the landmark development of the Scherrer method a century ago, multiple generations of width methods for X-ray diffraction originated to non-invasively and rapidly characterize the property-controlling sizes of nanoparticles, nanowires, and nanocrystalline materials. However, the predictive power of this approach suffers from inconsistencies among numerous methods and from misinterpretations of the results. Therefore, we systematically evaluated twenty-two width methods on a representative nanomaterial subjected to thermal and mechanical loads. To bypass experimental complications and enable a 1:1 comparison between ground truths and the results of width methods, we produced virtual X-ray diffractograms from atomistic simulations. These simulations realistically captured the trends that we observed in experimental synchrotron diffraction. To comprehensively survey the width methods and to guide future investigations, we introduced a consistent, descriptive nomenclature. Alarmingly, our results demonstrated that popular width methods, especially the Williamson-Hall methods, can produce dramatically incorrect trends. We also showed that the simple Scherrer methods and the rare Energy methods can well characterize unloaded and loaded states, respectively. Overall, this work improved the utility of X-ray diffraction in experimentally evaluating a variety of nanomaterials by guiding the selection and interpretation of width methods.
The predictions of scaling laws for the structure and properties of defect clusters are generally limited to small defect clusters in their ground-state configurations. We investigated the size and geometrical configuration dependence of nano-sized defect clusters in niobium (Nb) using molecular dynamics. We studied the structure and stability of large clusters of size up to fifty defects for vacancies and one hundred defects for interstitials, as well as the role of helium and metastable configurations on the stability of these clusters. We compared three different interatomic potentials in order to determine the relative stability of these clusters as a function of their size and geometrical configurations. Additionally, we conducted a statistical analysis to predict the formation and binding energies of interstitial clusters as a function of both their size and configuration. We find that the size dependence of vacancy and interstitial clusters can be approximated by functional forms that account for bulk and surface effects as well as some considerations of elastic interactions. We also find that helium and metastable configurations can make vacancy and interstitial clusters thermally stable depending on the configuration. Our parameterized functional forms for the formation and binding energies are valid for a very broad range of defect sizes and configurations making it possible to be used directly in a coarse-grained modeling strategy such as Monte Carlo, cluster dynamics or dislocation dynamics which look at defect accumulation and evolution in microstructures.
The embrittling or strengthening effect of solute atoms at grain boundaries (GBs), commonly known as the embrittling potency, is an essential thermodynamic property for characterizing the effects of solute segregation on GB fracture. One of the more technologically relevant material systems related to embrittlement is the Ni-S system where S has a deleterious effect on fracture behavior in polycrystalline Ni. In this work, we develop a Ni-S embedded-atom method (EAM) interatomic potential that accounts for the embrittling behavior of S at Ni GBs. Results using this new interatomic potential are then compared to previous density functional theory studies and a reactive force-field potential via a layer-by-layer segregation analysis. Our potential shows strong agreement with existing literature and performs well in predicting properties that are not included in the fitting database. Finally, we calculate embrittling potencies and segregation energies for six [100] symmetric-tilt GBs using the new EAM potential. We observe that embrittling potency is dependent on GB structure, indicating that specific GBs are more susceptible to sulfur-induced embrittlement.