Search NASA⌕ Search

SEARCH · Search NASA

Results for “crystal orientation maps”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

72 records · Page 4

Coupling Microstructural Evolution Simulations to Material Property Degradation Predictions for Plasma-Facing Materials

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

Connect microstructure evolution to property degradation with validated simulation

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

X-ray linear dichroic tomography of crystallographic and topological defects

Abstract The functionality of materials is determined by their composition 1–4 and microstructure, that is, the distribution and orientation of crystalline grains, grain boundaries and the defects within them 5,6 . Until now, characterization techniques that map the distribution of grains, their orientation and the presence of defects have been limited to surface investigations, to spatial resolutions of a few hundred nanometres or to systems of thickness around 100 nm, thus requiring destructive sample preparation for measurements and preventing the study of system-representative volumes or the investigation of materials under operational conditions 7–15 . Here we present X-ray linear dichroic orientation tomography (XL-DOT), a quantitative, non-invasive technique that allows for an intragranular and intergranular characterization of extended polycrystalline and non-crystalline 16 materials in three dimensions. We present the detailed characterization of a polycrystalline sample of vanadium pentoxide (V 2 O 5 ), a key catalyst in the production of sulfuric acid 17 . We determine the nanoscale composition, microstructure and crystal orientation throughout the polycrystalline sample with 73 nm spatial resolution. We identify and characterize grains, as well as twist, tilt and twin grain boundaries. We further observe the creation and annihilation of topological defects promoted by the presence of volume crystallographic defects. The non-destructive and spectroscopic nature of our method opens the door to operando combined chemical and microstructural investigations 11,18 of functional materials, including energy, mechanical and quantum materials.

Science & Technology - Other Topics↗

Multibeam Electron Diffraction

One of the primary uses for transmission electron microscopy (TEM) is to measure diffraction pattern images in order to determine a crystal structure and orientation. In nanobeam electron diffraction (NBED), we scan a moderately converged electron probe over the sample to acquire thousands or even millions of sequential diffraction images, a technique that is especially appropriate for polycrystalline samples. However, due to the large Ewald sphere of TEM, excitation of Bragg peaks can be extremely sensitive to sample tilt, varying strongly for even a few degrees of sample tilt for crystalline samples. In this paper, we present multibeam electron diffraction (MBED), where multiple probe-forming apertures are used to create multiple scanning transmission electron microscopy (STEM) probes, all of which interact with the sample simultaneously. We detail designs for MBED experiments, and a method for using a focused ion beam to produce MBED apertures. We show the efficacy of the MBED technique for crystalline orientation mapping using both simulations and proof-of-principle experiments. We also show how the angular information in MBED can be used to perform 3D tomographic reconstruction of samples without needing to tilt or scan the sample multiple times. Finally, we also discuss future opportunities for the MBED method.

47 OTHER INSTRUMENTATION↗

Mapping periplasmic binding protein oligosaccharide recognition with neutron crystallography

Numerous studies have shown how periplasmic binding proteins (PBPs) bind substrates with exquisite specificity, even distinguishing between sugar epimers and anomers, or structurally similar ions. Yet, marked substrate promiscuity is also a feature encoded in some PBPs. Except for three sub-Ångström crystal structures, there are no reports of hydrogen atom positions in the remaining (> 1000) PBP structures. The previous X-ray crystal structure of the maltodextrin periplasmic-binding protein from Thermotoga maritima (tmMBP) complexed with oligosaccharide showed a large network of interconnected water molecules stretching from one end of the substrate binding pocket to the other. These water molecules are positioned to form multiple hydrogen bonds, as well as forming interactions between the protein and substrate. Here we present the neutron crystal structure of tmMBP to a resolution of 2.1 Å. This is the first neutron crystal structure from the PBP superfamily and here we unambiguously identify the nature and orientation of the hydrogen bonding and water-mediated interactions involved in stabilizing a tetrasaccharide in the binding site. More broadly, these results demonstrate the conserved intricate mechanisms that underlie substrate-specificity and affinity in PBPs.

59 BASIC BIOLOGICAL SCIENCES↗

Self-Organized Stress Distributions in Polycrystalline Materials [Dissertation]

Understanding stress distributions in solid materials is complicated by the fact that most materials are polycrystalline in nature, with each crystal having an elastically anisotropic reaction to force. This study is to gain a better understanding on how external forces placed on a polycrystal are related to internal reactions within and between the grains. The hypothesis is that the stress distribution in porous to fully dense materials are self-organized based on strong contacts between and within the individual grains created by force chains. Force chains, commonly known in loaded granular materials, and could be the phenomenon that connect micro to macro deformation. Scale bridging measurements conducted through Raman spectroscopy, Atomic Force Microscopy, and Digital Image Correlation will be used to create stress maps, modulus maps, and elastic strain maps across a variety of geological and pharmaceutical polycrystals. When possible, the resulting maps will be compared to current homogenization schemes and a full field models. Finite element modeling will be used to assess whether the patterning seen in the experimental map is a reasonable approximation based on the orientation data of the samples used. A minimum of three publications is projected to be accomplished focusing each on a different method to experimentally test and analyze stress distributions.

36 MATERIALS SCIENCE↗

Machine learning inversion of interatomic force constants from single-crystal inelastic neutron scattering

Atomic vibrations govern many macroscopic properties of materials, but experiments to comprehensively probe them remain challenging. Inelastic neutron scattering (INS) is a powerful technique to map phonon dispersions in crystals, especially when leveraging modern time-of-flight (ToF) spectrometers with large detectors. However, efficiently and robustly extracting interatomic force constants (FCs) parameterizing phonon dynamics from experimental spectra remains a bottleneck due to the complexity and high dimensionality of ToF INS datasets. Here, we present a machine learning approach for the direct inversion of FCs from single-crystal INS measurements. The framework leverages synthetic training data generated using universal machine-learned force fields and an efficient physics-based forward model. We benchmark two neural architectures–one emphasizing structured latent representation learning and the other direct, supervised spectral regression–across simulated datasets for two materials under idealized and noisy conditions. The latent-representation model is subsequently applied to experimental single-crystal INS data on germanium. The model is shown to reproduce FCs derived from both first-principles simulations and from iterative optimization, and furthermore achieves reliable inference even from sparse, single-orientation measurements representing short data acquisitions. Analysis of the learned latent space reveals semantically continuous and physically interpretable encodings that support strong cross-domain generalization. By bridging theoretical and experimental domains, we establish a path toward rapid inversion of experimental spectra and data-driven interpretation of temperature-dependent lattice dynamics.

42 ENGINEERING↗

Analysis of X-ray images and spectra (aXis2000): A toolkit for the analysis of X-ray spectromicroscopy data

Spectromicroscopy refers to analytical methods that combine imaging and spectroscopy to provide detailed, spatially resolved analytical information about a sample, such as the type and quantitative spatial distributions of chemical components, geometric or magnetic alignment information, crystal structure, etc. The analysis of X-ray images and spectra (aXis2000) software described in this work provides a set of routines within a single, integrated, graphical-oriented package to read, display, manipulate and analyze spectromicroscopy data, with particular focus on soft X-ray spectromicroscopy methods such as scanning transmission X-ray microscopy (STXM), X-ray photoemission electron microscopy (XPEEM), scanning photoelectron X-ray microscopy (SPEM) and transmission X-ray microscopy (TXM). Here, this free software is described and compared to other software that can provide similar or complementary capabilities. Examples of spectromicroscopic analyses using advanced features of aXis2000 are provided.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NXREFINE

NXRefine implements a complete workflow for both data acquisition and reduction of single crystal diffuse x-ray scattering collected as a series of area detector frames during continuous sample rotation. The data are stored in HDF5 files, which conform to the NeXus data format standard. These files contain a comprehensive set of metadata, including the incident wavelength, powder calibration scans used to define detector distances, beam centers, and yaw, pitch, and roll corrections. Automated peak searches of the three-dimensional data arrays define a set of Bragg peaks that are used to define an orientation matrix, defined according to Busing and Levy, which is then used to transform the data from angular coordinates to reciprocal space coordinates using the CCTW program. Multiple rotations at different detector translations and offset rotation axes are combined to ensure that there are no gaps in the reconstructed data, and allow the construction of masks to eliminate detector artifacts due to Compton scattering within the sensor layers. Finally, the data can be transformed into 3D-ΔPDF maps using the punch-and-fill method after symmetrization.

OSBORN, RAYMOND↗

Methods for Mapping of Surface Charge Density on Arbitrary Dielectric Objects

A system to map surface charge densities on dielectric objects – e.g. PTFE, PMMA, and PA6 – is described to study maximum surface charge limits discharge and decay characteristics. Two electromechanical movements are used to achieve these goals: a three-axis robotic arm with an electrostatic probe and a three-axis auxiliary movement fitted with a current viewing resistor and laser displacement sensor. Dielectric objects no larger than a cylinder of eight-inch diameter and twelve inches in height are analyzed to create surface point clouds and generate a set of scan points. The arbitrary object is positioned upon a rotational axis with a laser displacement sensor oriented radially towards the rotational axis of the platform to determine the object’s physical limits.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Interface-driven magnetic anisotropy in relaxed La 0.7 Sr 0.3 CrO 3 /La 0.7 Sr 0.3 MnO 3 heterostructures on MgO

We investigate the structural and magnetic properties of La 0.7 Sr 0.3 CrO 3 /La 0.7 Sr 0.3 MnO 3 (LSMO) heterostructures grown on (001)-oriented MgO by molecular beam epitaxy. Due to the large film-substrate lattice mismatch, strain relaxation is found to occur within the first 2–3 unit cells (uc) of the film as evidenced by reflection high energy electron diffraction and high-resolution synchrotron x-ray reciprocal space mapping. Furthermore, we find that the presence of the LSCO spacer and capping layers leads to ferromagnetism in ultra-thin LSMO layers with thicknesses of the order of 2 uc with the magnetic easy axis oriented in the film plane. Net magnetic moments of 1.4 and 2.4 μ B /Mn are measured for (2 uc LSCO/2 uc LSMO) and (2 uc LSCO/4 uc LSMO) superlattices, respectively, by superconducting quantum interference device magnetometry. The effective magnetic anisotropy of the relaxed (2 uc LSCO/4 uc LSMO) heterostructure is found to be an order of magnitude higher than bulk LSMO highlighting the critical role of interfacial interactions in tuning magnetic anisotropy at complex oxide interfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Torsional force microscopy of van der Waals moirés and atomic lattices

In a stack of atomically thin van der Waals layers, introducing interlayer twist creates a moiré superlattice whose period is a function of twist angle. Changes in that twist angle of even hundredths of a degree can dramatically transform the system’s electronic properties. Setting a precise and uniform twist angle for a stack remains difficult; hence, determining that twist angle and mapping its spatial variation is very important. Techniques have emerged to do this by imaging the moiré, but most of these require sophisticated infrastructure, time-consuming sample preparation beyond stack synthesis, or both. In this work, we show that torsional force microscopy (TFM), a scanning probe technique sensitive to dynamic friction, can reveal surface and shallow subsurface structure of van der Waals stacks on multiple length scales: the moirés formed between bi-layers of graphene and between graphene and hexagonal boron nitride (hBN) and also the atomic crystal lattices of graphene and hBN. In TFM, torsional motion of an Atomic Force Microscope (AFM) cantilever is monitored as it is actively driven at a torsional resonance while a feedback loop maintains contact at a set force with the sample surface. TFM works at room temperature in air, with no need for an electrical bias between the tip and the sample, making it applicable to a wide array of samples. It should enable determination of precise structural information including twist angles and strain in moiré superlattices and crystallographic orientation of van der Waals flakes to support predictable moiré heterostructure fabrication.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The role of magnetic anisotropy in the magnetoresistance of Cr 2 O 3 /Al 2 O 3 thin film antiferromagnets

The magnetic states of antiferromagnetic insulating thin films are a promising medium for information storage, but characterization of these states has proven to be challenging. One approach is via magnetotransport measurements in an adjacent heavy metal layer. To this end, we synthesized and characterized a series of Cr 2 O 3 films and bilayers on Al 2 O 3 substrates with three different orientations: m-plane, a-plane, and c-plane. X-ray diffraction results demonstrated orientation control of the Cr 2 O 3 thin film, with m-plane films displaying a higher degree of mosaic spread than the a- and c-plane films. Reciprocal space maps showed that the films are mostly relaxed, although there was a small and different degree of strain in each orientation. The m-plane films were under 2% compressive strain, the a-plane film was under 0.5% compressive strain, and the c-plane film was completely relaxed to bulk values. To probe the magnetic state of the films, we measured the angular dependent magnetoresistance of Cr 2 O 3 /Pt bilayers for each orientation. We found a nontrivial temperature dependence of the sign of the magnetoresistance, pointing to the complex interplay between the exchange and anisotropy energies that vary with orientation. We propose that strain and mosaic spread may contribute to a difference in magnetic anisotropies among the samples and the resulting temperature dependence of the magnetoresistance. This work demonstrates the importance of considering the competition between antiferromagnetic exchange and magnetic anisotropy when storing information in the spin state of an antiferromagnetic insulator.

42 ENGINEERING↗

Atomic-Scale Scanning of Domain Network in the Ferroelectric HfO 2 Thin Film

Ferroelectric HfO 2 -based thin films have attracted much interest in the utilization of ferroelectricity at the nanoscale for next-generation electronic devices. However, the structural origin and stabilization mechanism of the ferroelectric phase are not understood because the film is typically nanocrystalline with active yet stochastic ferroelectric domains. Here, in this study, electron microscopy is used to map the in-plane domain network structures of epitaxially grown ferroelectric Y:HfO 2 films in atomic resolution. The ferroelectricity is confirmed in free-standing Y:HfO 2 films, allowing for investigating the structural origin for their ferroelectricity by 4D-STEM, high-resolution STEM, and iDPC-STEM. At the grain boundaries of <111>-oriented Pca2 1 orthorhombic grains, a high-symmetry mixed-(R3m, Pnm2 1 ) phase is induced, exhibiting enhanced polarization due to in-plane compressive strain. Nanoscale Pca2 1 orthorhombic grains and their grain boundaries with mixed-(R3m, Pnm2 1 ) phases of higher symmetry cooperatively determine the ferroelectricity of the Y:HfO 2 film. It is also found that such ferroelectric domain networks emerge when the film thickness is beyond a finite value. Furthermore, in-plane mapping of oxygen positions overlaid on ferroelectric domains discloses that polarization is suppressed at vertical domain walls, while it is active when domains are aligned horizontally with subangstrom domain walls. In addition, randomly distributed 180° charged domain walls are confined by spacer layers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing Fundamental Mechanisms of Plastic Deformation with High Energy X-rays (Final Progress Report)

Overall, substantial progress was made on developing BCDI techniques in the direction of imaging defects in multi-crystal samples. Section 2.1 describes a new method to reconstruct the 3D atomic displacement field within a grain from multiple coherently resolved diffraction peaks. The key development is the ability to recover a best fit to the grain shape at the same time as the internal displacement (strain) variation with no prior knowledge of the microstructure. This is complemented by the development of a method, Sec. 2.2, to separate out sets of Laue reflections that belong to a single grain (and thus be able to index its orientation) in the commonly encountered situation in which the beam illuminates several grains in a given location. This is essential for enabling collection of coherent diffraction data (at APS 34-ID-C) from neighboring grains. All experiments depend on the quality of samples used. Although good progress has been made with thin film samples, we are working with a broader range of geometries. A new sample design uses a high-Z metal (e.g., Pt) deposited in a sub-micron trench in a low-Z metal (e.g., Al) so that the grain size is prevented from exceeding the current limit for BCDI. We also report in section 3.1 on successfully applying NF-HEDM to a fine-grain Mg sample that would normally be considered out of scope for the method. This represents progress towards being able to map out polycrystalline samples in 3D before performing BCDI on individual grains. We further report in section 3.2 on preliminary NF-HEDM characterization of a sample of UO2 that has sub-micron grains; the data imply the capability of mapping fine-grain polycrystals at 1-ID-E, again for the purposes of BCDI in the bulk. Lastly, our plan is to load grains in polycrystalline samples by inducing thermal stress. We also report the experimental measurements of the Ti- 6Al-4V alloy obtained at the European Synchrotron Radiation Facility (ESRF), which had been indexed and mapped via HEDM at APS in section 9. As discussed below, however, technique development pre-occupied the team throughout the project which meant that we did not make any significant progress towards answering the original hypotheses as stated in the proposal. In summary, we assert that we have made good progress in methods and samples and future work taking advantage of this progress is likely to obtain and analyze BCDI data from multiple adjacent grains in 3D polycrystals. Obtaining such data depends on re-gaining access to the Advanced Photon Source, which is anticipated in late summer or the Fall. In all cases, work has involved close collaboration with faculty, staff, postdocs and students from Brigham Young University (BYU), the Los Alamos National Laboratory and APS beamlines 1-ID and 34-ID.

36 MATERIALS SCIENCE↗

In-situ fabrication and transport properties of (111) Y 2 Ir 2 O 7 epitaxial thin film

We report on the successful fabrication of (111) oriented thin films of pyrochlore iridate Y 2 Ir 2 O 7 by means of solid phase epitaxy. Essentially, we developed an entirely in situ annealing protocol, which is superior to the conventional ex-situ routine that requires multi-hour annealing to stabilize the proper pyrochlore structure. The morphological quality of the films has been confirmed by x-ray diffraction, reflectivity, and reciprocal space mapping, demonstrating their high crystallinity with a pure pyrochlore phase and an expected epitaxial relation to the substrate. Strikingly, below the magnetic phase transition, the temperature dependence of resistivity shows a power-law behavior with the power exponent of similar to 3/2, characteristic of a Weyl semimetal in the presence of impurities and electron-electron correlations. Finally, below 10K, the magneto-transport measurements reveal the emergence of a non-saturated negative magnetoresistance up to 9T with a small "valley" around zero field. These findings are further interpreted in relation to the weak anti-localization effect and the non-collinear antiferromagnetic ordering on the Ir sublattice.

36 MATERIALS SCIENCE↗

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗

Combining Theory and Experiment to Map the Atomic-Level Structure–Energy Pathways of Adsorbate-Mediated Phase Changes in a Cooperatively Flexible Metal–Organic Framework

An important subclass of metal–organic frameworks (MOFs) exhibits cooperative flexibility, wherein individual crystallites undergo global structural phase changes in response to external stimuli. Where cooperative flexibility results in reversible changes between crystalline states of distinct accessible porosity, these frameworks can exhibit rare yet desirable behaviors that cannot be explained by local dynamics alone. Yet, the chemical and structural origins of cooperative flexibility and how frameworks undergo these reversible phase changes at the atomic level remain poorly understood. Deliberate design for specific applications is therefore exceedingly difficult, and there is great impetus to develop a fundamental understanding of this phenomenon. Here, an effective and widely accessible computational approach is developed, which is designed to provide microscopic resolution via direct comparison to experimental data along the desorption-guided pathway. The strategy is applied to explain the desorption-induced phase change in an experimentally well-characterized framework, CdIF-13 (sod-Cd(benzimidazolate)2), where experiment alone was unable to resolve the atomistically detailed phase change landscape. Our findings reveal that the cooperative phase change pathways are adsorbate dependent with thermodynamics of intermediate structural states dictated by a nuanced interplay of ligand orientation, skeletal symmetry, and modes of surface adsorption. The results reveal that this isotropically flexible framework is “chaperoned” through a complex energy landscape by specific adsorbates, revealed by the reported computational approach with atomic-level insight and validated by experimentally determined structures. Thus, this work facilitates both understanding and future design of flexible materials for applications in gas storage, transport, delivery, and separation technologies.

03 NATURAL GAS↗