Search for: All records

This investigation presents a generally applicable framework for parameterizing interatomic potentials to accurately capture large deformation pathways. It incorporates a multi-objective genetic algorithm, training and screening property sets, and correlation and principal component analyses. The framework enables iterative definition of properties in the training and screening sets, guided by correlation relationships between properties, aiming to achieve optimal parametrizations for properties of interest. Specifically, the performance of increasingly complex potentials, Buckingham, Stillinger-Weber, Tersoff, and modified reactive empirical bond-order potentials are compared. Using MoSe₂as a case study, we demonstrate good reproducibility of training/screening properties and superior transferability. For MoSe₂, the best performance is achieved using the Tersoff potential, which is ascribed to its apparent higher flexibility embedded in its functional form. These results should facilitate the selection and parametrization of interatomic potentials for exploring mechanical and phononic properties of a large library of two-dimensional and bulk materials.

 

Vacancy engineering of 2H-transition metal dichalcogenides (2H-TMDs) has recently attracted great attention due to its potential to fine-tune the phonon and opto-electric properties of these materials. From a mechanical perspective, this symmetry-breaking process typically reduces the overall crack resistance of the material and adversely affects its reliability. However, vacancies can trigger the formation of heterogeneous phases that synergistically improve fracture properties. In this study, using MoSe2 as an example, we characterize the types and density of vacancies that can emerge under electron irradiation and quantify their effect on fracture. Molecular dynamic (MD) simulations, employing a re-parameterized Tersoff potential capable of accurately capturing bond dissociation and structural phase changes, reveal that isolated transition metal monovacancies or chalcogenide divacancies tend to arrest the crack tip and hence enhance the monolayer toughness. In contrast, isolated chalcogenide monovacancies do not significantly affect toughness. The investigation further reveals that selenium vacancy lines, formed by high electron dose rates, alter the crack propagating direction and lead to multiple crack kinking. Using atomic displacements and virial stresses together with a continuum mapping, displacement, strain, and stress fields are computed to extract mechanistic information, e.g., conditions for crack kinking and size effects in fracture events. The study also reveals the potential of specific defect patterns, “vacancy engineering,” to improve the toughness of 2H-TMDs materials. 

Combining experimental and computational studies of nanocomposite interfaces is highly needed to gain insight into their performance. However, there are very few literature reports, combining well-controlled atomic force microscopy experiments with molecular dynamic simulations, which explore the role of polymer chemistry and assembly on interface adhesion and shear strength. In this work, we investigate graphene oxide (GO)-polymer interfaces prevalent in nanocomposites based on a nacre-like architectures. We examine the interfacial strength resulting from van der Waals and hydrogen bonding interactions by comparing the out-of-plane separation and in-plane shear deformations of GO-polyethylene glycol (PEG) and GO-polyvinyl alcohol (PVA). The investigation reveals an overall better mechanical performance for the anhydrous GO-PVA system in both out-of-plane and in-plane deformation modes, highlighting the benefits of the donor-acceptor hydrogen bond formation present in GO-PVA. Such bond formation results in interchain hydrogen bond networks leading to stronger interfaces. By contrast, PEG, a hydrogen bond acceptor only, relies primarily on van der Waals inter-chain interactions, typically resulting in weaker interactions. The study also predicts that water addition increases the adhesion of GOPEG but decreases the adhesion of GO-PVA, and slightly increases the shear strength in both systems. Furthermore, by comparing simulations and experiments, we show that the CHARMM force field has enough accuracy to capture the effect of polymer content, water distribution, and to provide quantitative guidance for achieving optimum interfacial properties. Therefore, the study demonstrates an effective methodology, in the Materials Genome spirit, toward the design of 2D materials-polymer nanocomposites system for applications demanding mechanical robustness. 

Annihilation of vacancy clusters in monolayer molybdenum diselenide (MoSe2) under electron beam irradiation is reported. In situ high-resolution transmission electron microscopy observation reveals that the annihilation is achieved by diffusion of vacancies to the free edge near the vacancy clusters. Monte Carlo simulations confirm that it is energetically favorable for the vacancies to locate at the free edge. By computing the minimum energy path for the annihilation of one vacancy cluster as a case study, it is further shown that electron beam irradiation and pre-stress in the suspended MoSe2 monolayer are necessary for the vacancies to overcome the energy barriers for diffusion. The findings suggest a new mechanism of vacancy healing in 2D materials and broaden the capability of electron beam for defect engineering of 2D materials, a promising way of tuning their properties for engineering applications.