Probing Structural Defects in MOFs Using Water Stability.

Defects in the crystal structures of metal-organic frameworks (MOFs), whether present intrinsically or introduced via so-called defect engineering, can play strong roles in the properties of MOFs for various applications. Unfortunately, direct experimental detection and characterization of defects in MOFs are very challenging. We show that in many cases, the differences between experimentally observed and computationally predicted water stabilities of MOFs can be used to deduce information on the presence of point defects in real materials. Most computational studies of MOFs consider these materials to be defect-free, and in many cases, the resulting structures are predicted to be hydrophobic. Systematic experimental studies, however, have shown that many MOFs are hydrophilic. We show that the existence of chemically plausible point defects can often account for this discrepancy and use this observation in combination with detailed molecular simulations to assess the impact of local defects and flexibility in a variety of MOFs for which defects had not been considered previously.

Many-Particle Li Ion Dynamics in LiMPO4 Olivine Phosphates (M = Mn, Fe).

LiMPO4 (M = Mn, Fe) olivine phosphates are important materials for battery applications due to their stability, safety, and reliable recharge cycle. Despite continuous experimental and computational investigations, several aspects of these materials remain challenging, including conductivity dimensionality and how it maps onto Li pathways. In this work, we use a refined version of our finite temperature molecular dynamics "shooting" approach, originally designed to enhance Li hopping probability. We perform a comparative analysis of ion mobility in both materials, focused on many-particle effects. Therein, we identify main [010] diffusion channels, as well as means of interchannel couplings, in the form of Li lateral [001] hopping, which markedly impact the overall mobility efficiency as measured by self-diffusion coefficients. This clearly supports the need of many-particle approaches for reliable mechanistic investigations and for battery materials benchmarking due to the complex nature of the diffusion and transport mechanisms.

Incorporating Flexibility Effects into Metal-Organic Framework Adsorption Simulations Using Different Models.

High-throughput calculations based on molecular simulations to predict the adsorption of molecules inside metal-organic frameworks (MOFs) have become a useful complement to experimental efforts to identify promising adsorbents for chemical separations and storage. For computational convenience, all existing efforts of this kind have relied on simulations in which the MOF is approximated as rigid. In this paper, we use extensive adsorption-relaxation simulations that fully include MOF flexibility effects to explore the validity of the rigid framework approximation. We also examine the accuracy of several approximate methods to incorporate framework flexibility that are more computationally efficient than adsorption-relaxation calculations. We first benchmark various models of MOF flexibility for four MOFs with well-established CO2 experimental consensus isotherms. We then consider a range of adsorption properties, including Henry's constants, nondilute loadings, and adsorption selectivity, for seven adsorbates in 15 MOFs randomly selected from the CoRE MOF database. Our results indicate that in many MOFs adsorption-relaxation simulations are necessary to make quantitative predictions of adsorption, particularly for adsorption at dilute concentrations, although more standard calculations based on rigid structures can provide useful information. Finally, we investigate whether a correlation exists between the elastic properties of empty MOFs and the importance of including framework flexibility in making accurate predictions of molecular adsorption. Our results did not identify a simple correlation of this type.

Metashooting: a novel tool for free energy reconstruction from polymorphic phase transition mechanisms.

We introduce a novel scheme for the mechanistic investigation of solid-solid phase transitions, which we dub 'metashooting'. Combining transition path sampling, molecular dynamics and metadynamics, this scheme allows for both a complete mechanistic analysis and detailed mapping of the free energy surface. This is illustrated by performing metashooting calculations on the pressure-induced B4/B3 → B1 phase transition in ZnO. The resulting free energy map helps to clarify the role of intermediate configurations along this activated process and the competition between different mechanistic regimes with superior accuracy. We argue that metashooting can be efficiently applied to a broader class of activated processes.

Mechanism of the fcc-to-hcp phase transformation in solid Ar.

We present an atomistic description of the fcc-to-hcp transformation mechanism in solid argon (Ar) obtained from transition path sampling molecular dynamics simulation. The phase transition pathways collected during the sampling for an 8000-particle system reveal three transition types according to the lattice deformation and relaxation details. In all three transition types, we see a critical accumulation of defects and uniform growth of a less ordered transition state, followed by a homogeneous growth of an ordered phase. Stacking disorder is discussed to describe the transition process and the cooperative motions of atoms in {111} planes. We investigate nucleation with a larger system: in a system of 18 000 particles, the collective movements of atoms required for this transition are facilitated by the formation and growth of stacking faults. However, the enthalpy barrier is still far beyond the thermal fluctuation. The high barrier explains previous experimental observations of the inaccessibility of the bulk transition at low pressure and its sluggishness even at extremely high pressure. The transition mechanism in bulk Ar is different from Ar nanoclusters as the orthorhombic intermediate structure proposed for the latter is not observed in any of our simulations.

Hierarchical thermoelectrics: crystal grain boundaries as scalable phonon scatterers.

Thermoelectric materials are strategically valuable for sustainable development, as they allow for the generation of electrical energy from wasted heat. In recent years several strategies have demonstrated some efficiency in improving thermoelectric properties. Dopants affect carrier concentration, while thermal conductivity can be influenced by alloying and nanostructuring. Features at the nanoscale positively contribute to scattering phonons, however those with long mean free paths remain difficult to alter. Here we use the concept of hierarchical nano-grains to demonstrate thermal conductivity reduction in rocksalt lead chalcogenides. We demonstrate that grains can be obtained by taking advantage of the reconstructions along the phase transition path that connects the rocksalt structure to its high-pressure form. Since grain features naturally change as a function of size, they impact thermal conductivity over different length scales. To understand this effect we use a combination of advanced molecular dynamics techniques to engineer grains and to evaluate thermal conductivity in PbSe. By affecting grain morphologies only, i.e. at constant chemistry, two distinct effects emerge: the lattice thermal conductivity is significantly lowered with respect to the perfect crystal, and its temperature dependence is markedly suppressed. This is due to an increased scattering of low-frequency phonons by grain boundaries over different size scales. Along this line we propose a viable process to produce hierarchical thermoelectric materials by applying pressure via a mechanical load or a shockwave as a novel paradigm for material design.

Unexpected stable stoichiometries of sodium chlorides.

Sodium chloride (NaCl), or rocksalt, is well characterized at ambient pressure. As a result of the large electronegativity difference between Na and Cl atoms, it has highly ionic chemical bonding (with 1:1 stoichiometry dictated by charge balance) and B1-type crystal structure. By combining theoretical predictions and diamond anvil cell experiments, we found that new materials with different stoichiometries emerge at high pressures. Compounds such as Na3Cl, Na2Cl, Na3Cl2, NaCl3, and NaCl7 are theoretically stable and have unusual bonding and electronic properties. To test this prediction, we synthesized cubic and orthorhombic NaCl3 and two-dimensional metallic tetragonal Na3Cl. These experiments establish that compounds violating chemical intuition can be thermodynamically stable even in simple systems at nonambient conditions.

Understanding the nature of "superhard graphite".

Numerous experiments showed that on cold compression graphite transforms into a new superhard and transparent allotrope. Several structures with different topologies have been proposed for this phase. While experimental data are compatible with most of these models, the only way to solve this puzzle is to find which structure is kinetically easiest to form. Using state-of-the-art molecular-dynamics transition path sampling simulations, we investigate kinetic pathways of the pressure-induced transformation of graphite to various superhard candidate structures. Unlike hitherto applied methods for elucidating nature of superhard graphite, transition path sampling realistically models nucleation events necessary for physically meaningful transformation kinetics. We demonstrate that nucleation mechanism and kinetics lead to M-carbon as the final product. W-carbon, initially competitor to M-carbon, is ruled out by phase growth. Bct-C4 structure is not expected to be produced by cold compression due to less probable nucleation and higher barrier of formation.

Polarized cluster dynamics at the paraelectric to ferroelectric phase transition in BaTiO3.

The mechanism of the prototype ferroelectric phase transition in BaTiO(3) is a matter of intense debate and to a large extent still wrapped in mystery. Precursor phenomena in the form of polarized clusters in paraelectric BaTiO(3) are by now believed to represent a key step into the ferroelectric phenomenon. The determination of a slower dynamics of cluster polarization flipping along with a faster order-disorder Ti hopping mode among <111> off-center sites suggests coexistence, instead of mutual exclusion, of displacive and order-disorder types, initially proposed as distinct models. However, no clear picture of the transition state has been proposed so far, which is able to provide insight into the coexistence of the paraelectric and ferroelectric phenomena. Here, by means of a dedicated molecular dynamics approach, we provide a detailed atomistic picture of intermediate regions along the transition. Therein, different time and length scales coexist as they characterize different portions of the same material. From an imbalance of dynamically and more statically polarized clusters in this highly inhomogeneous intermediate, a symmetry breaking step naturally results. Further, we find that ferroelectric nanodomains may host antiferroelectric defects, which appear as an intrinsic feature of the growing BaTiO(3) ferroelectric material.