Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic structure”

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 181 records · Page 10

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗

Phonon Screening of Excitons in Atomically Thin Semiconductors

Atomically thin semiconductors, encompassing both 2D materials and quantum wells, exhibit a pronounced enhancement of excitonic effects due to geometric confinement. Consequently, these materials have become foundational platforms for the exploration and utilization of excitons. Recent ab initio studies have demonstrated that phonons can substantially screen electron-hole interactions in bulk semiconductors and strongly modify the properties of excitons. While excitonic properties of atomically thin semiconductors have been the subject of extensive theoretical investigations, the role of phonon screening on excitons in atomically thin structures remains unexplored. In this Letter, we demonstrate via ab initio GW-Bethe-Salpeter equation calculations that phonon screening can have a significant impact on optical excitations in atomically thin semiconductors. We further show that the degree of phonon screening can be tuned by structural engineering. We focus on atomically thin GaN quantum wells embedded in AlN and identify specific phonons in the surrounding material, AlN, that dramatically alter the lowest-lying exciton in monolayer GaN via screening. Our studies provide new intuition beyond standard models into the interplay among structural properties, phonon characteristics, and exciton properties in atomically thin semiconductors, and have implications for future experiments.

2-dimensional systems↗

Strain-driven oxygen vacancy ordering in LaNiO 3 thin films revealed by integrated differential phase contrast imaging in scanning transmission electron microscopy

Rare-earth nickelates, such as LaNiO 3 (LNO), exhibit complex electronic properties, with ordered oxygen vacancies (OOV) influencing conductivity and magnetic behavior. We investigate the structural stability of strain-induced OOV phases in LNO thin films grown on SrTiO 3 substrates and the impact of Ruddlesden–Popper (RP) faults. Using high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) and integrated differential phase contrast (iDPC) STEM imaging, we conducted atomic-scale structural and compositional analyses of OOV. Geometric phase analysis (GPA) was employed to measure the strain in fault-free and RP fault regions, while density functional theory (DFT) calculations explored different OOV arrangements in the LNO phase. Simulated iDPC-STEM imaging of energy-stabilized structures was performed to correlate with experimental results. Here, our findings reveal superstructure modulation in the chemical composition and atomic-scale lattice structure in LNO, primarily due to the formation of the OOV in Ni–O layers of the LaNiO 2.5 phase. The out-of-plane compressive strain of about 2% stabilizes this phase, reducing the strain, diminishing OOV, and transforming them into LNO.

36 MATERIALS SCIENCE↗

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interplay Between Metastability and Mechanically Induced Structural Instability in CsPbBr 3 Photovoltaic Perovskite

An atomic-level understanding of the underlying structural metastability is still absent in halide perovskite photovoltaic systems. Focusing on the model material CsPbBr 3 , the impact of mechanically induced atomic structure alteration is elucidated through structural modeling and X-ray diffraction measurements. Sudden Cs–Br bond breaking drives the system metastability, where the first-order transition arises from cation-halide bond strain relief by severing the corner connectivity of the PbBr 6 units. The pressure–volume and entropy terms govern the Gibbs free-energy landscape in halide perovskites. In conclusion, metastability is revealed as a critical factor limiting the performance of lead halide perovskites under extreme but natural environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Charge-induced atomic strain as a predictor of structural phase transformation in rare-earth intermetallics

We present a descriptor based on charge-induced atomic strain in crystalline lattices for predicting structural phase transformations in rare-earth intermetallic compounds containing lanthanides and transition metals. The charge-induced local atomic strain was obtained from structural optimization of experimentally known crystalline phases using state of the art density-functional theory methods. The predictive power of the descriptor was evaluated on 𝑅⁢𝐸 2 ⁢In (𝑅𝐸 = rare earth) compounds, a class known for diverse phase transformations. We show that incorporating quantum-mechanical effects—such as local charge distribution, bonding, symmetry, and electronic structure—enhances the robustness of the descriptor. To gain further insight, we analyzed phononic and electronic behavior in Y 2 ⁢In and demonstrated that experimental phase transformations are captured only when atomic strain effects are included. The descriptor was further used to predict structural phase changes in (Y⁢b 1–𝑥 ⁢E⁢r 𝑥 ) 2 ⁢In and G⁡d 2 ⁡(I⁢n 1–𝑥⁢ A⁢l 𝑥 ), with predictions confirmed by x-ray powder diffraction. Here, while the current study is focused on lanthanide-based intermetallics, the underlying principles of the descriptor suggest potential applicability to other closely related classes of rare-earth intermetallics.

Density functional theory↗

Pyrite (001) Interface Chemistry is Controlled by a Sulfoxy Termination

Pyrite (FeS 2 ) is the most common sulfide mineral on Earth, forming through inorganic reactions in the crust and oceanic hydrothermal systems and via microbially driven processes in anaerobic sediments. The pyrite–water interface is the site of a wide range of adsorption and reaction processes in Earth systems including oxidation that dramatically affects the geochemistry of surface waters and influences global carbon and oxygen cycles. Mechanistic geochemical models of pyrite interfacial reactivity, however, are limited by the lack of experimentally derived atomistic structures of the reduced and reacting surfaces. Here, in this work, we reveal the atomic-scale structure of the pyrite (001)-water interface that forms at very low oxygen partial pressures, relevant to suboxic environments in Earth. The interface structure and surface speciation were obtained using the crystal truncation rod method supported by ambient-pressure photoelectron spectroscopy and density functional theory calculations. The surface is dominantly composed of disulfide groups bound to a single oxygen atom, forming a sulfoxy group that has no known molecular or bulk mineral analog. This surface is interpreted as the first step in the oxidative dissolution of pyrite. The sulfoxy group is readily protonated through surface acid–base reactions that alter the structure of interfacial water and the free energy of interfacial reactions. Surface iron sites are not oxidized. Surprisingly, this interface can likely develop in equilibrium with bulk pyrite in some reducing and acidic solutions. This termination is therefore likely representative of pyrite surfaces under a vast range of experimental, industrial and Earth conditions.

oxidation↗

Cryo2StructData: A Large Labeled Cryo-EM Density Map Dataset for AI-based Modeling of Protein Structures

The advent of single-particle cryo-electron microscopy (cryo-EM) has brought forth a new era of structural biology, enabling the routine determination of large biological molecules and their complexes at atomic resolution. The high-resolution structures of biological macromolecules and their complexes significantly expedite biomedical research and drug discovery. However, automatically and accurately building atomic models from high-resolution cryo-EM density maps is still time-consuming and challenging when template-based models are unavailable. Artificial intelligence (AI) methods such as deep learning trained on limited amount of labeled cryo-EM density maps generate inaccurate atomic models. To address this issue, we created a dataset called Cryo2StructData consisting of 7,600 preprocessed cryo-EM density maps whose voxels are labelled according to their corresponding known atomic structures for training and testing AI methods to build atomic models from cryo-EM density maps. Cryo2StructData is larger than existing, publicly available datasets for training AI methods to build atomic protein structures from cryo-EM density maps. We trained and tested deep learning models on Cryo2StructData to validate its quality showing that it is ready for being used to train and test AI methods for building atomic models.

59 BASIC BIOLOGICAL SCIENCES↗

Uncovering multiscale structure-property correlations via active learning in scanning tunneling microscopy

Atomic arrangements and local sub-structures fundamentally influence emergent material functionalities. These structures are conventionally probed using spatially resolved studies and the property correlations are deciphered by a researcher based on sequential explorations, thereby limiting the efficiency and scope. Here we demonstrate a multi-scale Bayesian deep-learning based framework that automatically correlates material structure with its electronic properties using scanning tunneling microscopy (STM) measurements in real-time. Its predictions are used to autonomously direct exploration toward regions of the sample that optimize a given material property. This method is deployed on a low-temperature ultra-high vacuum STM to understand the structure-property relationship in a europium-based semimetal, EuZn 2 As 2 , a promising candidate relevant to magnetism-driven topological phenomena. The framework employs a sparse-sampling approach to efficiently construct the scalar-property space using minimal measurements, about 1–10% of the data required in standard hyperspectral methods. Moreover, we formulate the problem hierarchically across length scales, implementing autonomous workflow to locate mesoscopic and atomic structures that correspond to a target material property. This framework offers the choice to design scalar-property from the spectroscopic data to steer sample exploration. Our findings reveal correlations of the electronic properties unique to surface terminations, local defect density, and point defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exploration of the Electronic and Catalytic Properties of [Co 5 MS 8 (PEt 3 ) 5 ] 1+ Nanoclusters: A Computational Study

Recent studies have demonstrated the relative stability of undercoordinated hexanuclear cobalt sulfide nanoclusters (NCs) with different charge states. Considering that these small metal NCs have atomically precise structures and high reactivity due to the open shell of the transition metals, and provide selectivity toward ligand loss, they are a vital model for catalysis. In this paper, the electronic structures of these NCs are investigated. These NCs are then used as the reference state to analyze the catalytic properties with respect to hydrogen evolution reaction (HER) and CO 2 reduction (CO 2 R). Further, to understand the effect of heteroatom incorporation, the geometry and reactivity of ten different metal dopants are analyzed. This work shows that the type of metal incorporation greatly affects the electronic structure and formation energies for ligand binding and catalysis. Particularly, the d-orbital occupancy in the cobalt atoms remains largely unchanged, while the heteroatom greatly influences the reactivity of the undercoordinated NCs. Most notably, this work highlights that transition metals in [Co 5 MS 8 (PEt 3 ) 5 ] 1+ NCs would competitively prefer electrochemical adsorption of H over COOH, while the main group metals prefer COOH adsorption.

catalysis↗

Locating the atoms at the hard-soft interface of gold nanoparticles

Surface structure affects the growth, shape and properties of nanoparticles. In wet chemical syntheses, metal additives and surfactants are used to modify surfaces and guide nanocrystal growth. To understand this process, it is critical to understand how the surface structure, and hence its energy, is modified. However, measuring the type and arrangement of atoms at hard-soft interfaces on nanoscale surfaces, especially in the presence of surfactants, is extremely challenging. Here, we determine the atomic structure of the hard-soft interface in a metallic nanoparticle by developing low-dose imaging conditions in four-dimensional scanning transmission electron microscopy that are preferentially sensitive to surface adatoms. By revealing experimentally the copper additives and bromide surfactant counterion at the surface of a gold nanocuboid and quantifying their interatomic distances, our direct, low-dose imaging method provides atomic-level understanding of chemically sophisticated nanomaterial surface structures. These measurements of the atomic structure of the hard-soft interface provide the information necessary to understand and quantify surface chemistries and energies and their pivotal role in nanocrystal growth.

Li, Weilun [Monash University, Melbourne, VIC (Aus↗

Diffraction-quality, ultraflexible protein single crystals engineered with DNA

DNA-functionalized colloidal nanoparticles assemble through flexible, nanoscale DNA hybridization interactions that limit atomic-level structural order. Here, we report a valence-centric strategy that enables DNA-bonded, protein single crystals with unconventional mechanical properties. An octameric enzyme, glutarate L-2-hydroxylase, was site- and number-selectively conjugated with eight self-complementary single-stranded DNA, yielding octavalent molecular bonds. The resulting conjugate assembled into the designed body-centered tetragonal crystals that diffracted to 1.42- to 2.61-angstrom resolution, with contacts mediated by B-form DNA helices spanning 17 to 25 angstroms. Increasing oligonucleotide length induces anisotropic lattice expansion while preserving atomic periodicity, even with partial DNA occupancy. Mechanistic studies suggest that the dynamic motion of unhybridized DNA facilitates crystallization, analogous to fluctuating electron clouds in atomic bonding. Compared with native protein crystals, DNA-hybridized crystals are 23-fold softer. These results challenge the assumption that flexibility is incompatible with structural order and establish a programmable framework for biomolecular crystallization and nanomaterials engineering with atomic precision.

Han, Zhenyu [Department of Chemistry, Northwestern↗

Unrecoverable lattice rotation governs structural degradation of single-crystalline cathodes

Transitioning from polycrystalline to single-crystalline nickel-rich cathodes has garnered considerable attention in both academia and industry, driven by advantages of high tap density and enhanced mechanical properties. However, cathodes with high nickel content (>70%) suffer from substantial capacity degradation, which poses a challenge to their commercial viability. Here, leveraging multiscale spatial resolution diffraction and imaging techniques, we observe that lattice rotations occur universally in single-crystalline cathodes and play a pivotal role in the structure degradation. These lattice rotations prove unrecoverable and govern the accumulation of adverse lattice distortions over repeated cycles, contributing to structural and mechanical degradation and fast capacity fade. These findings bridge the previous knowledge gap that exists in the mechanistic link between fast performance failure and atomic-scale structure degradation.

36 MATERIALS SCIENCE↗

Enhancing Chiroptoelectronic Activity in Chiral 2D Perovskites via Chiral–Achiral Cation Mixing

Rational design of chiral two-dimensional hybrid organic–inorganic perovskites is crucial to achieve chiroptoelecronic, spintronic, and ferroelectric applications. Here, in this study, an efficient way to manipulate the chiroptoelectronic activity of 2D lead iodide perovskites is reported by forming mixed chiral (R- or S-methylbenzylammonium (R-MBA + or S-MBA + )) and achiral (n-butylammonium (nBA + )) cations in the organic layer. The strongest and flipped circular dichroism signals are observed in (R/S-MBA 0.5 nBA 0.5 ) 2 PbI 4 films compared to (R/S-MBA) 2 PbI 4 . Moreover, the (R/S-MBA 0.5 nBA 0.5 ) 2 PbI 4 films exhibit pseudo-symmetric, unchanged circularly polarized photoluminescence peak as temperature increases. First-principles calculations reveal that mixed chiral–achiral cations enhance the asymmetric hydrogen-bonding interaction between the organic and inorganic layers, causing more structural distortion, thus, larger spin-polarized band-splitting than pure chiral cations. Temperature-dependent powder X-ray diffraction and pair distribution function structure studies show the compressed intralayer lattice with enlarged interlayer spacing and increased local ordering. Overall, this work demonstrates a new method to tune chiral and chiroptoelectronic properties and reveals their atomic scale structural origins.

2D perovskites↗

Data for: Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Recently realized high-mobility semiconducting epitaxial graphene on silicon carbide, provided an important step towards integration of the graphene-based system into active components in post-silicon micro- and nano-electronics. However, the exact atomic-scale structure and the complex bonding configurations of the first epitaxial graphene carbon layer remain an open problem. Our recent report has shed new light on understanding this interface, where the external transverse electric field-dependent dynamic switching behavior of the Cbuffer-SiC bonds was observed. Here, using scanning tunneling microscopy and spectroscopy (STM and STS), we present the direct evidence of silicon (Si) vacancies at the interface and provide their distribution at the topmost reconstructed SiC(0001) layer. Experimental STM and density functional theory modeling data were used in the preparation of figures in a published article in the Journal of Physical Chemistry Letters. Files related to the figures and supplementary materials in the article are present in this dataset in .txt format.

Condensed matter imaging↗

Data for: Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Recently realized high-mobility semiconducting epitaxial graphene on silicon carbide, provided an important step towards integration of the graphene-based system into active components in post-silicon micro- and nano-electronics. However, the exact atomic-scale structure and the complex bonding configurations of the first epitaxial graphene carbon layer remain an open problem. Our recent report has shed new light on understanding this interface, where the external transverse electric field-dependent dynamic switching behavior of the Cbuffer-SiC bonds was observed. Here, using scanning tunneling microscopy and spectroscopy (STM and STS), we present the direct evidence of silicon (Si) vacancies at the interface and provide their distribution at the topmost reconstructed SiC(0001) layer. Experimental STM and density functional theory modeling data were used in the preparation of figures in a published article in the Journal of Physical Chemistry Letters. Files related to the figures and supplementary materials in the article are present in this dataset in .txt format.

Condensed matter imaging↗

Evolution of Electronic Properties of Graphene Nanoribbons with Progressive Carving: From Straight to Porous to Chevron Ribbons

Graphene nanoribbons (GNRs) are highly versatile materials due to their unique electronic, magnetic, and optical properties, which can be precisely tuned by controlling their width, edge structure, and topology. Here, we report the on-surface synthesis and characterization of a straight N = 15 armchair GNR with periodic annulene nanopores (15-pGNR). It serves as a structural link between two well-established GNRs: the pristine N = 15 armchair GNR without pores (15-AGNR) and the chevron GNR (cGNR). With the addition of the 15-pGNR reported in this study, these three GNRs form a rare experimentally accessible series of ribbons, in which the evolution of electronic properties can be tracked upon progressive carving of a basic 15-AGNR: first, by creating periodic nanopores to form 15-pGNR and then by extending the pore area and producing meandering cGNR. We have designed a molecular precursor for the 15-pGNR and grown the nanoribbons on single-crystal gold substrates by on-surface synthesis in ultra-high vacuum (UHV) conditions. The atomically precise structure of 15-pGNR was confirmed by scanning tunneling microscopy (STM) and non-contact atomic force microscopy (nc-AFM). The band gap of 15-pGNR was studied by scanning tunneling spectroscopy (STS) and dI/dV mapping, and the occupied electronic levels were investigated by angle-resolved photoemission spectroscopy (ARPES). A theoretical and experimental comparison of 15-pGNRs, 15-AGNRs, and cGNRs demonstrates that the introduction of periodic nanopores into 15-AGNR leads to a more than 2-fold increase in its band gap. In contrast, the band gaps of 15-pGNR and cGNR differ only by about 15%. Such band gap increase can be qualitatively understood to arise from two combined effects, the periodic perforation of the graphene lattice and the confinement effect induced by the GNR width.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Recently realized high-mobility semiconducting epitaxial graphene on silicon carbide (SiC) [Zhao, J. Nature 2024, 625 (7993), 60−65, 10.1038/s41586-023-06811-0] provided an important step toward integration of the graphene-based system into active components in postsilicon micro- and nanoelectronics. However, the exact atomic-scale structure and complex bonding configurations of the first epitaxial graphene carbon layer (C buffer ) remain an open problem. Our recent report [Kolmer, M. Communications Physics 2024, 7 (1), 16, 10.1038/s42005-023-01515-3] has shed new light on understanding this interface, where the external transverse electric field-dependent dynamic switching behavior of the C buffer –SiC bonds was observed. Here, using scanning tunneling microscopy and spectroscopy (STM and STS), we present direct evidence of silicon (Si) vacancies at the interface and provide their distribution at the topmost reconstructed SiC(0001) layer. Bias voltage and epitaxial graphene thickness-dependent characterization of the collective C buffer –SiC interface showed that “Si” vacancy sites beneath C buffer are stable under STM electric fields. Moreover, the vacancies introduce localized electronic states below the Fermi level, thereby enhancing the charge-transfer phenomenon across the interface.

Thupakula, Umamahesh [Ames Laboratory (AMES), Ames↗