Search NASA⌕ Search

SEARCH · Search NASA

Results for “Screw”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Accelerating Laboratory-Based, All-Solid-State Battery Research and Development: Industrially Relevant Small-Batch Dry-Processing, Small and Low-Cost Test Fixtures, and Short-Tolerant Separators

This work presents improvements to mixing (for dry-processed electrodes), separator robustness, and cell testing jigs for all-solid-state batteries (ASSBs). These developments combine synergistically to enable rapid research and development of ASSB catholytes. A mechanical "kneader" is designed and built which emulates the mixing and fibrillation of poly-tetrafluoroethylene binder that occurs in industrially relevant twin-screw extruders. A modest improvement in the mixing and dispersion catholyte materials is observed using micro-resolution X-ray computed tomography (X-ray CT). Catholytes are paired with an In-metal anode and separated by a novel, dual-layer, polyaramid-fiber-supported, sulfide-solid electrolyte (Li6 PS5 Cl) separator. This separator achieves high short tolerance which enabled a 95% success rate across 20 attempted cells and to date over 100 successful cells have been built. ASSBs were tested in a 2032 coin-cell format. Simultaneous cycling of a large number of cells is further enabled by a compact and low cost (~ 38 USD per jig) pressure application jig of which 100 have been constructed to date. The utility of these advances is demonstrated through a brief study on the effect cathode-side conductive carbon and current collector type has on capacity and rate performance.

25 ENERGY STORAGE↗
DOE OSTI2026

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)↗
DOE OSTI2022

Hot-melt extrudability of amorphous solid dispersions of flubendazole-copovidone: An exploratory study of the effect of drug loading and the balance of adjuvants on extrudability and dissolution

The FDA-approved anthelmintic flubendazole has shown potential to be repositioned to treat cancer and dry macular degeneration; however, its poor water solubility limits its use. Amorphous solid dispersions may overcome this challenge, but the balance of excipients may impact the preparation method and drug release. Here, we evaluate the influence of adjuvants and drug loading on the development of an amorphous solid dispersion of flubendazole-copovidone by hot-melt extrusion. The drug, copovidone, and adjuvants (magnesium stearate and hydroxypropyl cellulose) mixtures were statistically designed, and the process was performed in a twin-screw extruder. The study showed that flubendazole and copovidone mixtures were highly extrudable, except when drug loading was high (>40%). Furthermore, magnesium stearate positively impacted the extrusion and was more effective than hydroxypropyl cellulose. The extruded materials were evaluated by modulated differential scanning calorimetry and X-ray powder diffraction, obtaining positive amorphization and physical stability results. Pair distribution function analysis indicated the presence of drug rich domains with medium-range order structure and no evidence of polymer-drug interaction. All extrudates presented faster dissolution (HCl, pH 1.2) than pure flubendazole, and both adjuvants had a notable influence on the dissolution rate. In conclusion, hot-melt extrusion may be a viable option to obtain stable flubendazole: copovidone amorphous dispersions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

The effect of local chemical ordering on Frank-Read source activation in a refractory multi-principal element alloy

In this work, we investigate the operation of Frank-Read (FR) sources in a refractory multi-principal element alloy (MPEA). Simulations of discrete dislocation motion in MPEAs is enabled by the development of a phase field dislocation dynamics model that treats the atomic-scale fluctuations in lattice energies across the glide plane, present due to local ordering in the chemical composition within the nominally random MPEA atomic structure. We consider, through simulation, a range of length scales over which ordering occurs, varying from short-range lengths, a few times dislocation core width, to long-range lengths, an order of magnitude longer than the core. Characteristic of this body-centered cubic MPEA, the simulations also include screw/edge character dependence in glide resistance, as informed by atomic scale simulation. The critical stresses to activate the source for the same source size are found to be statistically distributed, as a direct consequence of the underlying variation in lattice energy. Furthermore, analysis of critical state for activating edge and screw FR sources in the MPEA reveals that FR source operation occurs via a two-step mechanism, involving athermal kink-pair formation, unlike the conventional FR source operation in a material with no composition fluctuations. This mechanism lowers the average critical stress required to activate the FR source and causes the statistical dispersion in critical stress to depend on the range of composition ordering. More importantly, it leads to a more severe dependence of source strength on FR source length than predicted by line tension alone.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Cross-slip and easy-glide CRSS in titanium: Theoretical predictions and in-situ TEM measurements

This study investigates the mechanics of prismatic and first-order pyramidal $\langle$a$\rangle$ slip in titanium (Ti), elucidating the physics of easy-glide and cross-slip through a combination of theory and experiments. Screw-character prismatic (Pr) dislocations in Ti are of particular interest because their complex cores can be stable or unstable, leading to activation by either cross-slip or planar glide. To investigate these mechanisms, site-specific micro tensile samples were prepared using focused ion beam (FIB) milling and mounted on a push-to-pull (PTP) device for in-situ transmission electron microscopy (TEM) tensile testing. The in-situ experiments provide direct observations of the onset of dislocation motion and the precise determination of the critical resolved shear stress (CRSS) for the activated mechanisms, and its evolution with load cycling. A comprehensive theory has been developed to predict the CRSS values for easy glide, cross-slip, and multiplication of dislocations. Predicted critical stresses for pyramidal (π)-to-Pr and reverse cross-slip agree closely with the experimental measurements. The latter cross-slip stress is a factor of two higher than that of unobstructed planar slip. The model accounts for overlapping dislocation cores and employs a Wigner-Seitz based cell to evaluate misfit energies. By combining ab initio density functional theory (DFT) with anisotropic elasticity, the framework identifies minimum energy pathways for dislocation glide, which can be intermittent and zig-zag. A simplified expression utilizing (π) and Pr Schmid factor ratios is proposed for critical stress corresponding to (π)-to-Pr cross-slip transition. The results are strongly dependent on crystal orientation, underscoring non-Schmid behavior. Overall, this study explores key critical stress parameters essential for informing higher-scale simulations of plasticity in Ti.

Critical stress↗
DOE OSTI2023

Investigation of defect formation at the early stage of PVT-grown 4H-SiC crystals

Here, several 4° off-axis 4H-SiC wafers with several hundred microns of initial-stage growth by PVT method are investigated by Synchrotron Monochromatic Beam X-ray Topography (SMBXT). Defect behavior across the seed/newly grown layer interface are demonstrated. Comparison of early stage grown layers to seed sample indicates generation of threading edge dislocation (TED) and threading screw dislocation (TSD)/ threading mixed dislocation (TMD) pairs at the interface while most basal plane dislocations (BPDs) are deflected into TEDs. The (0001) facet of the crystal is already formed at early stage growth at the edge of the wafer and high nitrogen incorporation in facet leads to conditions favorable for nucleation and glide of Shockley/double Shockley faults with layers of 3C-SiC deposited on facet acting as nuclei. Unique-shaped dislocations are observed at early stage growth, which are caused by deflection of TSDs/TMDs and TEDs by macrosteps near the periphery of the sample and the subsequent glide of the a components. The effect of the quality of the seed surface before growth is manifested as randomly oriented arrays of pairs of TEDs and TSDs/TMDs on the as grown surface resulting from residual surface damage from scratches.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stress effects on the energy barrier and mechanisms of cross-slip in FCC nickel

The energy barrier for homogeneous cross-slip of a screw dislocation in face-centered cubic (FCC) nickel is calculated using atomistic simulations as a function of the Escaig stress on the glide plane, the Escaig stress on the cross-slip plane, and the Schmid stress on the cross-slip plane. Two cross-slip mechanisms, Friedel-Escaig (FE) and Fleischer, are examined and their energy barriers are calculated for a large number of stress combinations. For each mechanism, the energy barrier as a function of three stress components can be reduced into a one-dimensional function of an effective stress. The stress domains in which FE and Fleischer mechanisms operate respectively are determined. Additionally, the FE mechanism dominates when the Escaig stress on the glide plane (in the direction that reduces the stacking fault width) is the largest stress component. Increasing the Schmid stress and Escaig stress (in the direction that expands the stacking fault width) on the cross slip plane promotes the Fleischer mechanism. The cross slip energy barrier functions obtained here can be used as input functions for computing cross slip rates in mesoscale dislocation dynamics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase field dislocation dynamics (PFDD) modeling of non-Schmid behavior in BCC metals informed by atomistic simulations

Body-Centered Cubic (BCC) metals exhibit anomalous mechanical properties such as twinning/anti-twinning and tension/compression asymmetry, which are attributed to asymmetric behavior of screw dislocations. Unlike Face-Centered Cubic (FCC) metals, the Critical Resolved Shear Stress (CRSS) of BCC metals deviates from the Schmid law. Here, we use a mesoscale modeling approach called Phase Field Dislocation Dynamics (PFDD) to understand the critical factors that control this non-Schmid behavior and reproduce atomistic predictions of the CRSS (Peierls stress). All inputs to the PFDD model are obtained from Molecular Statics (MS) simulations. Multiple pathways for modeling the non-Schmid behavior are investigated by incorporating the representative dislocation properties into different energy terms in the PFDD model. One way to understand non-Schmid behavior is to incorporate stress components projected on inclined planes into the external energy term within PFDD. Alternatively, we propose that non-Schmid behavior can also be accounted for by considering the variation of the {110} unstable stacking fault energy and the dislocation core width as a function of the applied tensorial stress field. The CRSS predicted using PFDD modeling is in excellent agreement with MS predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A model of thermal creep and annealing in finite domains based on coupled dislocation climb and vacancy diffusion

Here, we develop a framework to investigate thermal creep and annealing in finite domains, where the climb motion of discrete dislocations is coupled to the diffusion of a continuum vacancy field. The model is first formulated in a continuum finite-deformation setting. All governing equations and boundary conditions are obtained from a unified irreversible thermodynamics principle. The resulting model couples a mechanical boundary value problem (BVP), a vacancy diffusion BVP, and the climb and glide motion of the discrete dislocation network within the crystal. The framework is then linearized for implementation in three-dimensional (3D) discrete dislocation dynamics (DDD) simulations for arbitrary anisotropic crystals. A solution scheme is developed based on the superposition principle, which is imposed weakly on the dislocation network to obtain a Galerkin solution for the nodal climb velocities. The framework includes diffusional (Nabarro–Herring) creep deformation as well as dislocation creep by climb-assisted-glide. The method is applied to simulate the annealing of vacancy loops in Al, with good agreement to experimental measurements by Silcox and Hirsch. We further consider the effects of annealing under stress, and of the proximity of the vacancy loops to loaded and free boundaries Simulations in polycrystalline materials are carried out to highlight the effects of the grain size on dislocation climb and vacancy loop annealing. The method is also applied to estimate the creep rate due to climb-assisted glide of jogged-screw dislocations in γ-TiAl, and results are compared to experiments by Viswanathan et al. Finally, we discuss the effects of uniaxial and hydrostatic stresses on the two diffusive deformation pathways of the material, namely Nabarro–Herring creep and dislocation climb.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic determination of Peierls barriers of dislocation glide in nickel

The Peierls barrier measures the lattice resistance to dislocation glide in crystalline solids. We use the nudged elastic band (NEB) method to calculate the Peierls barriers for screw and edge dislocation glide in a face-centered cubic (FCC) metal of Ni. The minimum energy paths (MEPs) across single or sequential Peierls barriers are determined under shear loading. The NEB results show the decreasing Peierls barrier with increasing shear stress, giving the Peierls stress at which the Peierls barrier vanishes. The effects of boundary condition and system size on Peierls barriers are studied by comparing strain- and stress-controlled NEB results. Furthermore, the free-end NEB methods are applied to determine MEPs with improved computational efficiency. The NEB results are also used to evaluate the energetic driving force of dislocation glide, which is consistent with that determined from the Peach-Koehler force. The accuracy of the present NEB results based on an empirical interatomic potential is assessed by comparison with a machine-learning potential. This work demonstrates the robust and efficient quantification of Peierls barriers to dislocation glide in an FCC metal, and it lays a solid foundation for the atomistic determination of Peierls barriers in compositionally complex alloys with the FCC structure in future studies.

42 ENGINEERING↗
DOE OSTI2023

Phase-field modeling of dislocation–interstitial interactions

The mechanical behavior of body-centered cubic (BCC) materials can be dramatically affected by the presence of interstitial solute atoms. Here we present a new phase-field dislocation dynamics formulation to include the diffusion of interstitials. Short-range interactions are accounted for by a concentration-dependent lattice energy, and long-range interactions are accounted for by modifications to the elastic energy. The interstitial diffusion law introduces gradients that require methods for minimizing Gibbs oscillations, which is done via a modified Green’s function. The formulation is general to any solute-solvent system and is applied here to Nb-O as a model system, whose interstitial parameters are obtained from ab initio calculations. The effect of O on the core structures of Nb edge and screw dislocations is calculated. The diffusion of O to form interstitial atmospheres around dislocation cores is simulated, as well as the critical stresses required for dislocations to break away or cross slip from these atmospheres. Future applications of the method to simulate complex interstitial embrittlement mechanisms are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dislocation nano-hydrides in nickel: Nucleation, evolution and effects on dislocation behaviors

Nano-hydrides have been predicted to precipitate at the core of edge dislocations in the Ni-H system, a mechanism that may promote hydrogen embrittlement. However, nano-hydride nucleation, growth, and effects on dislocation behavior have seldom been explored. This work combines molecular dynamics grand canonical Monte Carlo (MD-GCMC) simulations and continuum modeling to uncover a wide range of phenomena linked to dislocation nano-hydrides. Simulations reveal that nano-hydrides can be stabilized at dislocation cores with all character angles, including screw segments, due to the hydrostatic stresses around the cores of the Shockley partials. Nano-hydride nucleation takes place in these regions, and growth is dictated by the character angle θ of the perfect dislocation. The equilibrium stacking fault width d eq varies dynamically to increase the local hydrostatic stress field and facilitate the formation of the nano-hydride, forming a constriction-like feature and leading to three distinct behaviors: d eq decreases for θ > 30°, d eq remains unchanged for θ = 30°, and d eq increases for θ < 30°. Remote hydrostatic and Escaig stresses are also shown to influence the nucleation stage, implying stress concentrations such as those ahead of crack tips may facilitate nano-hydride precipitation. Moreover, we identify a new hydrogen-induced 60° dislocation reaction that emits a Shockley partial on a conjugate plane, with potential implications for twin nucleation. Testable predictions from this study are then used to reinterpret previous results from the literature. In conclusion, these findings provide a comprehensive framework to assess nano-hydride formation and evolution at dislocations in nickel and other face-centered cubic metals, with important implications to hydrogen embrittlement.

Dislocations↗
DOE OSTI2022

Helium interaction with solutes and impurities in neutron-irradiated nanostructured ferritic alloys: A first principles study

Density functional theory calculations are performed to explore the binding between He and alloying solutes, impurities, and transmutation products expected in neutron irradiated nanostructured ferritic alloys (NFAs), here 14YWT is taken as an example. Elements that exhibit significant binding (attraction) with an interstitial He are Y (binding energy = 0.46 eV), Mg (0.32), O (0.33), Ti (0.16), and C (0.15). Those that provide significant binding to a substitutional He are O (1.44), Y (1.24), N (0.73), H (0.56), Mg (0.52), Ti (0.34), Si (0.34), C (0.33), Al (0.32), Ni (0.26), Ta (0.23), and Mn (0.16). The presence of these elements in Fe matrix could reduce the transport of He towards oxide particles, dislocations, and internal boundaries, and could promote He bubble nucleation in the matrix. For convenience, we compile existing binding energy data of He with He n and He n V (He-vacancy) clusters. Dissociation pathway analysis reveals that, in general, the most likely dissociation of a He n V cluster is by a sequential emission of individual He atoms. Furthermore, larger bubbles are more prone to dissociation than smaller ones. In addition, we estimate the binding energy (segregation energy) of He in bulk Y 2 Ti 2 O 7 (YTO) single crystal, YTO/Fe interface, and YTO particle embedded in Fe, with respect to interstitial He in Fe, from existing formation energies of He in these structures. We also compile available data of He binding with Fe self-interstitial atom (SIA), SIA clusters, and edge and screw dislocations. Note that given the absence of DFT data, the binding with SIA clusters and dislocations are gathered from simulations with empirical potentials. Finally, the data presented in this paper is important to inform multiscale simulations of He bubble accumulation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Monotonic testing of single shear-plane CLT-to-CLT joint with hardwood dowels

Reducing embodied carbon emissions in buildings and striving for carbon neutrality through sustainable design has become a primary goal in the construction industry. The industry is expanding its use of wood-engineered products, such as cross-laminated timber (CLT) to achieve those sustainability goals. Furthermore, replacing steel fasteners with lower embodied carbon alternatives is gaining increased attention in timber engineering. This paper aims to quantify the mechanical behavior of single shear-plane CLT-to-CLT joints using hardwood dowels. The study includes 154 single shear-plane experimental tests with four different CLT species, two hardwood dowel species, and two dowel diameters. Further, the test specimens were subjected to cyclic and monotonic loading until failure, and mechanical properties such as yielding and ultimate strength, serviceability and yielding stiffness, and ductility were quantified. The hardwood dowel CLT-to-CLT joints demonstrated high strength with up to 10.3 kN (2.3 kips) capacity per dowel and lateral stiffness of up to 5.2 kN/mm (30.8 k/in) per dowel. This strength and stiffness are roughly equivalent to or exceed values for mass timber wood screws, making hardwood dowels a good alternative. In addition, the observed mean ductility of joints ranged from 2.4 to 4.5. Finally, analytical equations for the joints' yielding strength, lateral serviceability stiffness, and yielding stiffness were derived with non-linear regression analysis and validated using the experimental test results.

42 ENGINEERING↗
DOE OSTI2024

Mechanical properties and deformation mechanisms of single crystal Mg micropillars subjected to high-strain-rate C-axis compression

Here, the mechanical properties and deformation mechanisms of single crystal magnesium under c-axis quasi-static and high-strain rate compressions are investigated through in situ scanning electron microscope (SEM) experiments and post-mortem transmission electron microscope (TEM) characterization. The findings revealed that ductility and high rates of hardening are preserved for pillars as large as 15 μm. Furthermore, rate effects result in a mild increase in flow stress with plastic deformations controlled primarily by the slip of type dislocations. Importantly and in contrast to other literature reports, plastic deformation occurs in the absence of twining. As the strain increases and plastic deformation exceeds about 4%, crystal rotation activates basal slip, <$\mathrm{a}$> type dislocations, resulting in a more rate independent flow stress. TEM observation on micropillars compressed at a strain rate of 250/s, revealed the activation of {${11}$$\bar{2}$$\bar{2}$} < $\bar{1}$$\bar{1}23$ > slip systems and high mobility of screw dislocations as major contributors to plastic strains in excess of 10% without fracture. These findings are relevant to the design of lightweight materials used in transportation systems, e.g., selection of material grain size. Moreover, the experimental data here reported provides the materials science community with a unique opportunity to validate discrete dislocation dynamics (DDD) formulations employed in multiscale design of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Polyethylene density uncertainty in the simulation of neutron detectors

Simulations with the Monte Carlo N-Particle® (MCNP®) code are frequently used to study neutron coincidence counters for international safeguards. Reducing sources of uncertainty enables more accurate simulations, which reduces the need for expensive or impossible measurements. Accurately specifying the density of polyethylene is crucial for accurate simulations. The density of polyethylene ranges from 0.88 to 0.97 g/cm 3 and high-density polyethylene (HDPE) ranges from about 0.944 to 0.965 g/cm 3 . The density of a specific detector should be measured exactly for the best simulations. Here, the density of the blue/white Active Well Coincidence Counter (AWCC) bulk polyethylene was measured to be 0.9594 +- 0.0013 g/cm 3 at 1-sigma confidence. The density was calculated from repeated measurements of the various polyethylene dimensions for volume, and weighing the polyethylene for the mass. The AWCC neutron detection efficiency was simulated for each polyethylene density. The relative range of efficiencies (3σ) was 3.4% for unspecified polyethylene, 0.9% for HDPE, and 0.3% for the measured density. The impact of modeling small features such as screws was also studied and was found to be a negligible 0.2%. Measuring polyethylene density can reduce its role in uncertainty from one of the largest to one of the smallest contributors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
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↗