Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic spectra”

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.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Scanning tunneling microscopy study of the antiferromagnetic topological insulator MnBi 2 Se 4

Topological insulators with intrinsic magnetic ordering such as MnBi 2 Se 4 have been predicted to host many exotic topological phenomena including dissipationless edge states and magnetoelectric effects. Here we use scanning tunneling microscopy to study the surface of MnBi 2 Se 4 epitaxial thin films. Large scale topographic images show steps between van der Waals layers and screw dislocations. While we find the Se-termination predominant on the surface, areas attributed to Bi-termination are found at the edges of small islands in as-grown films, or in larger areas upon removal of surface capping. Here, we compare these terminations through an analysis of step heights, atomic resolution images and tunneling spectroscopy. Tunneling spectra on the Se-termination reveals semiconducting behavior, with an apparent bandgap of >0.18eV that depends on the tip height. In contrast, tunneling spectra of the Bi-termination appears gapless, which could be an indication of dangling bond states on the surface. In both cases, a hexagonal 3.9Å lattice is observed that is consistent with the expected crystal structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Structure and Desorption Kinetics of Acetonitrile Thin Films on Pt(111) and on Graphene on Pt(111)

Acetonitrile thin films were prepared on Pt(111) and on graphene on Pt(111) using molecular beam techniques. Temperature programed desorption (TPD) experiments of acetonitrile on Pt(111) displayed first-order kinetics for monolayer desorption and zero-order kinetics for multilayer desorption. We used reflection adsorption infrared spectroscopy (RAIRS) to characterize the orientation of the monolayer on Pt(111). Acetonitrile molecules have been shown to adopt a bridging structure where the C=N bond breaks to form two single bonds to underlying platinum atoms. When Tdeposition = 120 K our RAIRS spectra are consistent with this bridging configuration. However, when Tdeposition was lowered to 60 K molecules retained the C=N bond and adopted a vertical orientation with the nitrogen closest to the platinum. On graphene, the TPD experiments show that both the monolayer and multilayer desorb with zero-order kinetics. Our RAIRS spectra suggest that the monolayer consists of acetonitrile molecules oriented with their molecular axes parallel to the graphene substrate. This orientation did not change when Tdeposition was varied from 60 to 120 K. Our results provide a more complete picture of how acetonitrile molecules interact with Pt(111) and with graphene surfaces.

desorption tin films, molecules, thin films, desor↗
DOE OSTI2023

Tracking Active Phase Behavior on Boron Nitride during the Oxidative Dehydrogenation of Propane Using Operando X-ray Raman Spectroscopy

Hexagonal boron nitride (hBN) is a highly selective catalyst for the oxidative dehydrogenation of propane (ODHP) to propylene. Using a variety of ex situ characterization techniques, the activity of the catalyst has been attributed to the formation of an amorphous boron oxyhydroxide surface layer. The ODHP reaction mechanism proceeds via a combination of surface mediated and gas phase propagated radical reactions with the relative importance of both depending on the surface-to-void-volume ratio. Here we demonstrate the unique capability of operando X-ray Raman spectroscopy (XRS) to investigate the oxyfunctionalization of the catalyst under reaction conditions (1 mm outer diameter reactor, 500 to 550 °C, P = 30 kPa C 3 H 8 , 15 kPa O 2 , 56 kPa He). We probe the effect of a water cofeed on the surface of the activated catalyst and find that water removes boron oxyhydroxide from the surface, resulting in a lower reaction rate when the surface reaction dominates and an enhanced reaction rate when the gas phase contribution dominates. Computational description of the surface transformations at an atomic-level combined with high precision XRS spectra simulations with the OCEAN code rationalize the experimental observations. Finally, this work establishes XRS as a powerful technique for the investigation of light element-containing catalysts under working conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Virtual node graph neural network for full phonon prediction

Understanding the structure-property relationship is crucial for designing materials with desired properties. The past few years have witnessed remarkable progress in machine-learning methods for this connection. However, substantial challenges remain, including the generalizability of models and prediction of properties with materials-dependent output dimensions. Here we present the virtual node graph neural network to address the challenges. By developing three virtual node approaches, we achieve Γ-phonon spectra and full phonon dispersion prediction from atomic coordinates. We show that, compared with the machine-learning interatomic potentials, our approach achieves orders-of-magnitude-higher efficiency with comparable to better accuracy. This allows us to generate databases for Γ-phonon containing over 146,000 materials and phonon band structures of zeolites. Additionally, our work provides an avenue for rapid and high-quality prediction of phonon band structures enabling materials design with desired phonon properties. The virtual node method also provides a generic method for machine-learning design with a high level of flexibility. In this study, the authors present a virtual node graph neural network to enable the prediction of material properties with variable output dimensions. This method offers fast and accurate predictions of phonon band structures in complex solids.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Resolving local structural motifs across the phase evolution of zinc titanates with computational x-ray absorption spectroscopy

Resolving the local structure motifs that characterize phase evolution as a function of composition is a key challenge in structure characterization of complex materials. Here, in this study, we combine first-principles simulations and x-ray absorption near-edge structures (XANES) analysis to gain insights into the structure evolution revealed by measurements across a combinatorial zinc titanate thin film, which was grown with smoothly varying composition over a wide range of the Ti:Zn ratio. Specifically, we propose a cluster blind-signal-separation (cBSS) method for XANES spectral analysis based on a library of the structures and spectra of representative local motifs. In addition to motifs from zinc titanate crystals, two types of Ti-defect models constructed in this study are key to the understanding of the structure characteristics in the Zn-rich region. The cBSS method makes use of both spectral clustering of the simulated site-XANES spectra library and the BSS procedure to construct high-fidelity spectral basis functions from an experimental spectral sequence. The method provides a rigorous measure of the spectral sensitivity and basis completeness. The results of the XANES analysis are corroborated with other experimental modalities, including x-ray diffraction and spectroscopic ellipsometry, to validate the cBSS method. The calculated motif weights resulting from fitting the XANES spectra with the cBSS basis probe the atomic structure characteristics of both crystalline and amorphous phases as a function of the Ti/Zn composition. The insights of the local structure motif evolution are pivotal to the understanding of the nonmonotonic trend in the optical gap, which may lead to potential applications through tuning the optical properties of zinc titanate. The workflow of the XANES spectral analysis developed in this work can be generalized to construct the structure-property relationship in a broad material space.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The R -process Alliance: A Nearly Complete R -process Abundance Template Derived from Ultraviolet Spectroscopy of the R -process-enhanced Metal-poor Star HD 222925 *

We present a nearly complete rapid neutron-capture process (r-process) chemical inventory of the metal-poor ([Fe/H] = -1.46 ± 0.10) r-process-enhanced ([Eu/Fe] = +1.32 ± 0.08) halo star HD 222925. This abundance set is the most complete for any object beyond the solar system, with a total of 63 metals detected and seven with upper limits. It comprises 42 elements from 31 ≤ Z ≤ 90, including elements rarely detected in r-process-enhanced stars, such as Ga, Ge, As, Se, Cd, In, Sn, Sb, Te, W, Re, Os, Ir, Pt, and Au. We derive these abundances from an analysis of 404 absorption lines in ultraviolet spectra collected using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope and previously analyzed optical spectra. A series of appendices discusses the atomic data and quality of fits for these lines. The r-process elements from Ba to Pb, including all elements at the third r-process peak, exhibit remarkable agreement with the solar r-process residuals, with a standard deviation of the differences of only 0.08 dex (17%). In contrast, deviations among the lighter elements from Ga to Te span nearly 1.4 dex, and they show distinct trends from Ga to Se, Nb through Cd, and In through Te. The r-process contribution to Ga, Ge, and As is small, and Se is the lightest element whose production is dominated by the r-process. The lanthanide fraction, log X La = -1.39 ± 0.09, is typical for r-process-enhanced stars and higher than that of the kilonova from the GW170817 neutron-star merger event. We advocate adopting this pattern as an alternative to the solar r-process-element residuals when confronting future theoretical models of heavy-element nucleosynthesis with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Internal Friction Measures to Study Precipitates Formation in EP and N-Doped Bulk Nb for SRF Applications

Main focus of this study is the investigation of thermodynamics phenomena responsible for the High Field Q Slope (HFQS) in SRF cavities by Internal Friction (IF) measurement. Mechanical spectroscopy is, indeed, a well-established technique to study precipitate formations in BCC materials and several works on the effects of impurities as N and O on the Snoek peak have been published so far and will be taken as reference to explain the mechanisms behind the observed dissipation effects. Internal Friction measurements were performed in Belgium at IMCE on Nb rectangular shape samples with different RRR values prepared at Fermilab by using Electro Polishing (EP), N-doping and heat treatments in order to reproduce the same conditions during the standard treatments applied on bulk Nb SRF cavities. From IF spectra, the H trapping mechanism by interstitial atoms (N and O and/or vacancies, depending on the purity level, RRR) can be easily recognized leading to results that perfectly corroborate previous findings on Q-disease, HFQS and RRR phenomena.

Spina, Tiziana↗
DOE OSTI2021

Transformation of atomically dispersed platinum in SAPO-37 into platinum clusters: catalyst for ethylene hydrogenation

Atomically dispersed supported platinum catalysts were synthesized by the reaction of Pt(acac) 2 (acac = acetylacetonato) with the silicoaluminophosphate molecular sieve SAPO-37, with infrared spectra showing that the reaction involved SAPO OH groups. Extended X-ray absorption fine structure (EXAFS) spectra show that, after heating in air to 623 K, each platinum atom on average was bonded to approximately two oxygen atoms of the SAPO support, with no evidence of platinum clusters. X-ray absorption near edge spectra (XANES) indicate a platinum formal oxidation state of approximately +2. Here, these supported mononuclear species lacked catalytic activity for ethylene hydrogenation at 1 bar and room temperature, but activity arose during catalysis in a once-through flow reactor/X-ray absorption spectroscopy cell, as the XANES white line intensity decreased, indicating reduction of the platinum, and the EXAFS-determined Pt–Pt coordination number (CN) increased from essentially zero to 2.7 ± 0.5 after 2 h, indicating the formation of platinum clusters of only a few atoms each. The data indicate the almost instantaneous formation of the clusters and identify them as the catalytically active species. Subsequent exposure of the catalyst to ethylene as in situ EXAFS and XANES data were recorded led to a decrease in the Pt–Pt coordination number to 1.6 ± 0.3 and an increase in the white line intensity, indicating partial oxidative fragmentation of the clusters by ethylene. Platinum clusters in SAPO-37 formed in separate experiments by exposure to H 2 prior to catalysis were also catalytically active for ethylene hydrogenation. Here, the data thus show how to dial in the catalytic activity by forming platinum clusters by exposure of mononuclear platinum to H 2 and how to dial back the activity by oxidative fragmentation of the clusters by reaction with ethylene, all under mild conditions; the ethylene : H 2 ratio of the reactant gas determines the average platinum nuclearity and thus the catalytic activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Deducing subnanometer cluster size and shape distributions of heterogeneous supported catalysts

Abstract Infrared (IR) spectra of adsorbate vibrational modes are sensitive to adsorbate/metal interactions, accurate, and easily obtainable in-situ or operando. While they are the gold standards for characterizing single-crystals and large nanoparticles, analogous spectra for highly dispersed heterogeneous catalysts consisting of single-atoms and ultra-small clusters are lacking. Here, we combine data-based approaches with physics-driven surrogate models to generate synthetic IR spectra from first-principles. We bypass the vast combinatorial space of clusters by determining viable, low-energy structures using machine-learned Hamiltonians, genetic algorithm optimization, and grand canonical Monte Carlo calculations. We obtain first-principles vibrations on this tractable ensemble and generate single-cluster primary spectra analogous to pure component gas-phase IR spectra. With such spectra as standards, we predict cluster size distributions from computational and experimental data, demonstrated in the case of CO adsorption on Pd/CeO 2 (111) catalysts, and quantify uncertainty using Bayesian Inference. We discuss extensions for characterizing complex materials towards closing the materials gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

IDB Loader

SAND2024-08428O The IDB Loader is a database designed to hold gamma spectra. IDB is a web-accessible database of reference gamma spectra for measuring uranium (U) and plutonium (Pu) isotopic composition. The database, developed by the International Atomic Energy Agency (IAEA), provides access to well-characterized groups of gamma spectra curated by international experts in gamma spectroscopy. It was created to promote sustainability and maintenance of software used to determine the isotopic abundances of U and Pu. The IDB Loader works by uploading one set of spectral data into an existing IDB. The Data Loader package includes documentation describing the IDB's table, the format of the spectral data, and a description of the loader program. The IDB is currently hosted by the IAEA at: https://nds.iaea.org/idb Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Schwartz, Steven↗
DOE OSTI2021

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Characterization of light-induced potentials in the strong-field dissociation of O 2 +

In this work, we investigate theoretically the imprints of light-induced potentials (LIPs) on the dissociation dynamics of O 2 + molecular ions, as observable in angle-resolved fragment kinetic-energy-release (KER) spectra. Following the vibrational and rotational dynamics of the initial pump-laser-excited cationic nuclear wave packet, while accounting for the dipole coupling between the O 2 + ( a 4 Π u ) and O 2 + ( f 4 Π g ) electronic states in 800-nm 40-fs probe-laser pulses with peak intensities between 10 13 and 10 14 W / cm 2 , we calculate angle-resolved KER spectra which reveal characteristic energy- and angle-dependent fringe structures. These fringes shift downward in energy as the molecular alignment angle θ relative to the probe-pulse polarization direction increases from 0 to π / 2 . The angle-dependent shifts in the KER fringes increase for larger probe-pulse peak intensities and follow the angle and light-wave-intensity dependence of the vibrational spectrum in the associated Floquet bond-hardening well, which is a manifestation of transient O 2 + nuclear-probability trapping in the LIP during dissociation. By examining the rovibrational dynamics of the dissociating molecular cation near the light-induced conical intersection (LICI) in the cationic LIP surface at θ = π / 2 , we identify related angle-dependent structures in the KER spectra, suggesting a means for assessing the significance of LICIs in molecular dissociation pathways.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

NEXAFS spectroscopy of alkylated benzothienobenzothiophene thin films at the carbon and sulfur K-edges

Alkylated benzothienobenzothiophenes are an important class of organic semiconductors that exhibit high performance in solution-processed organic field-effect transistors. In this work, we study the near-edge x-ray absorption fine-structure (NEXAFS) spectra of 2,7-didecyl[1]benzothieno[3,2-b][1]benzothiophene (C 10 -BTBT) at both the carbon and sulfur K-edges. Angle-resolved experiments of thin films are performed to characterize the dichroism associated with molecular orientation. First-principles calculations using the density functional theory-based many-body x-ray absorption spectroscopy (MBXAS) method are also performed to correlate the peaks observed and their dichroism with transitions to specific antibonding molecular orbitals. Interestingly, the dichroism of the dominant, lowest energy peak is opposite at the carbon and sulfur K-edges. While the low-energy peak at the carbon K-edge is assigned to carbon 1s → π* transitions with transition dipole moment (TDM) perpendicular to the planar BTBT core, the dominant low energy peak at the sulfur K-edge is assigned to sulfur 1s → σ* transitions with TDM oriented along the long axis of the BTBT core. These differences at the sulfur and carbon K-edges are understood through the MBAXS simulations that find a reordering of the energy of the lowest energy π* and σ* transitions at the sulfur K-edge due to the strong localization of the σ* orbital over the sulfur atom. This work highlights differences in the NEXAFS spectra of organic semiconductors at carbon and sulfur K-edges and provides new insights into peak assignment and x-ray dichroism relevant for studying the molecular orientation of organic semiconductor films.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The microwave spectra of the conformers of $\mathcal{n}$-butyl nitrate

We report the microwave spectrum of n-butyl nitrate was recorded in the 5 to 20 GHz frequency range using broadband chirp and narrowband pulse excitation molecular jet Fourier transform microwave spectrometers. A quantum chemistry structural analysis yielded thirteen stable conformers. Among them, the five most energetically stable conformers were observed in the experimental spectra. The most stable conformer features a butyl chain with an anti-gauche-anti conformation (AGA) where the γ-carbon atom is about 64° out of the nitrate plane. For this conformer, spectra of all 13 C and 15 N minor isotopologues could be measured. The conformer with a straight butyl chain (AAA), and three other conformers (GAA, GGA, and AGG) were also observed. Accurate rotational constants, centrifugal distortion constants, and 14 N nuclear quadrupole coupling constants could be deduced and compared to the theoretical values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Lyα Observations of Comet C/2013 A1 (Siding Spring) Using MAVEN IUVS Echelle

The close approach of comet C/2013 A1 (Siding Spring) to Mars in 2014 October provided a unique opportunity to observe a dynamically new Oort cloud comet with potential for interaction with a planet’s atmosphere. The water-originating hydrogen coma of the comet extended to over 20 million km from the nucleus. Determining the properties of this coma contributes to characterizing the comet’s water content and production rate. The present study analyzes a unique data set of high spectral resolution UV observations of comet C/2013 A1 Siding Spring measured by the Mars Atmosphere and Volatile Evolution spacecraft. The Siding Spring observations capture Lyα emissions from the Martian corona, the interplanetary medium, as well as the cometary H and D reservoirs. The isolated cometary spectra are analyzed to reveal a velocity distribution of H atoms that are consistent with model estimates of H{sub 2}O photodissociated H emissions and of OH photodissociated H emissions, Doppler shifted from the main comet H emission line center by 18 km s{sup −1} and 8 km s{sup −1}, respectively. The variations in comet H brightness with distance from the nucleus are used to constrain cometary water production to a rate of 0.5 × 10{sup 28} molecules s{sup −1} at a time when Siding Spring was at 1.5 au, pre-perihelion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Enhanced cutoff energies for direct and rescattered strong-field photoelectron emission of plasmonic nanoparticles

Abstract The efficient generation, accurate detection, and detailed physical tracking of energetic electrons are of applied interest for high harmonics generation, electron-impact spectroscopy, and femtosecond time-resolved scanning tunneling microscopy. We here investigate the generation of photoelectrons (PEs) by exposing plasmonic nanostructures to intense laser pulses in the infrared (IR) spectral regime and analyze the sensitivity of PE spectra to competing elementary interactions for direct and rescattered photoemission pathways. Specifically, we measured and numerically simulated emitted PE momentum distributions from prototypical spherical gold nanoparticles (NPs) with diameters between 5 and 70 nm generated by short laser pulses with peak intensities of 8.0 × 10 12 and 1.2 × 10 13 W/cm 2 , demonstrating the shaping of PE spectra by the Coulomb repulsion between PEs, accumulating residual charges on the NP, and induced plasmonic electric fields. Compared to well-understood rescattering PE cutoff energies for strong-field photoemission from gaseous atomic targets (10× the ponderomotive energy), our measured and simulated PE spectra reveal a dramatic cutoff-energy increase of two orders of magnitude with a significantly higher contribution from direct photoemission. Our findings indicate that direct PEs reach up to 93 % of the rescattered electron cutoff energy, in contrast to 20 % for gaseous atoms, suggesting a novel scheme for the development of compact tunable tabletop electron sources.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Ab initio molecular dynamics simulation of Nd3+ incorporation in calcite

Trivalent rare-earth elements (REEs) readily incorporate in Ca-bearing minerals such as calcite (CaCO3), a phenomenon with important implications for the use of REEs as biogeochemical tracers as well as for environmental remediation. Despite substantial efforts, determining the incorporation modes of trivalent REEs in calcite has remained challenging due to their presence in multiple coordination environments. This is particularly problematic for the interpretation of extended X-ray absorption fine structure (EXAFS) spectroscopy. In this work, ab initio molecular dynamics (AIMD) simulations of Nd3+ incorporated in calcite were performed to resolve its incorporation modes. Nd3+ and Ca2+ have very similar ionic radii but the additional positive charge requires a charge compensation scheme (CCS). Five CCSs were considered in AIMD simulations, and the resulting atomic trajectories were then used to compute EXAFS spectra for direct fitting to a published experimental spectrum. The results of the fits indicated that Nd3+ incorporates at the Ca2+ site in calcite in a mixture of six- and seven-fold coordination environments and through three main incorporation modes: associated with a protonated nearest-neighbor Ca2+ vacancy, incorporated as NdOH2+, and charge compensated by a remote species (on the spatial scale probed by EXAFS). The incorporation of Nd3+ in calcite via multiple modes was consistent with previous work, but the AIMD-EXAFS approach allowed for resolving the nature, proportion, and structure of the different coordination environments.

Kerisit, Sebastien N.↗
DOE OSTI2021

Engineering Thermal Transport across Layered Graphene–MoS 2 Superlattices

Layering two-dimensional van der Waals materials provides a high degree of control over atomic placement, which could enable tailoring of vibrational spectra and heat flow at the sub-nanometer scale. In this work, using spatially resolved ultrafast thermoreflectance and spectroscopy, we uncover the design rules governing cross-plane heat transport in superlattices assembled from monolayers of graphene (G) and MoS 2 (M). Using a combinatorial experimental approach, we probe nine different stacking sequences, G, GG, MG, GGG, GMG, GGMG, GMGG, GMMG, and GMGMG, and identify the effects of vibrational mismatch, interlayer adhesion, and junction asymmetry on thermal transport. Pure G sequences display evidence of quasi-ballistic transport, whereas adding even a single M layer strongly disrupts heat conduction. The experimental data are described well by molecular dynamics simulations, which include thermal expansion, accounting for the effect of finite temperature on the interlayer spacing. The simulations show that an increase of ~2.4% in the layer separation of GMGMG, relative to its value at 300 K, can lead to a doubling of the thermal resistance. Using these design rules, we experimentally demonstrate a five-layer GMGMG superlattice “thermal metamaterial” with an ultralow effective cross-plane thermal conductivity comparable to that of air.

2D materials↗