Search NASA⌕ Search

SEARCH · Search NASA

Results for “developable planes”

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 1,369 records · Page 76

WE-AB-207A-12: HLCC Based Quantitative Evaluation Method of Image Artifact in Dental CBCT

Purpose: Image artifacts are usually evaluated qualitatively via visual observation of the reconstructed images, which is susceptible to subjective factors due to the lack of an objective evaluation criterion. In this work, we propose a Helgason-Ludwig consistency condition (HLCC) based evaluation method to quantify the severity level of different image artifacts in dental CBCT. Methods: Our evaluation method consists of four step: 1) Acquire Cone beam CT(CBCT) projection; 2) Convert 3D CBCT projection to fan-beam projection by extracting its central plane projection; 3) Convert fan-beam projection to parallel-beam projection utilizing sinogram-based rebinning algorithm or detail-based rebinning algorithm; 4) Obtain HLCC profile by integrating parallel-beam projection per view and calculate wave percentage and variance of the HLCC profile, which can be used to describe the severity level of image artifacts. Results: Several sets of dental CBCT projections containing only one type of artifact (i.e. geometry, scatter, beam hardening, lag and noise artifact), were simulated using gDRR, a GPU tool developed for efficient, accurate, and realistic simulation of CBCT Projections. These simulated CBCT projections were used to test our proposed method. HLCC profile wave percentage and variance induced by geometry distortion are about 3∼21 times and 16∼393 times as large as that of the artifact-free projection, respectively. The increase factor of wave percentage and variance are 6 and133 times for beam hardening, 19 and 1184 times for scatter, and 4 and16 times for lag artifacts, respectively. In contrast, for noisy projection the wave percentage, variance and inconsistency level are almost the same with those of the noise-free one. Conclusion: We have proposed a quantitative evaluation method of image artifact based on HLCC theory. According to our simulation results, the severity of different artifact types is found to be in a following order: Scatter>Geometry>Beam hardening>Lag>Noise>Artifact-free in dental CBCT.

60 APPLIED LIFE SCIENCES↗

Beam-, Target- and Double- Spin Asymmetry Measurements in Deeply Virtual π 0 , π + and π - Meson Production

This thesis presents measurements of spin asymmetries in single pion electroproduction from the EG1-DVCS experiment carried out at the Thomas Jefferson National Accelerator Facility in 2009 with use of the CEBAF Large Acceptance Spectrometer. The experiment made use of a 6 GeV longitudinally polarised electron beam and (dynamically) longitudinally polarised targets of 14 NH 3 and 14 ND 3 (giving access to protons and deuterons), providing the means to measure three spin asymmetries: beam-, target-, and double-spin. Deeply Virtual Meson Production for the π 0 , π + and π - channels were studied for the kinematic range of Q 2 > 1 GeV 2 and W > 2 GeV. The collection of measurements includes the analysis of numerous reconstruction topologies, as well as making use of different target materials. These are the world's first charged pion spin asymmetry measurements in this kinematic regime. The sPlot maximum likelihood fit based background subtraction technique was implemented in the data analysis, showing strong agreement with a traditional and commonly used cuts-based approach. This novel technique was thoroughly tested to develop implementation and execution. The results obtained will provide a reference point for future analyses at the upgraded experimental facility and motivate the use of the sPlot technique. These measurements will provide insight into the structure of nucleons through their interpretation in relation to Generalised Parton Distributions which describe substructure correlations between longitudinal momentum and transverse position of partons. In particular, here the spin asymmetries dependence on t , the squared four-momentum transfer to the target nucleon, carries information on the spatial parton distributions in the transverse plane.

Murdoch, Gavin↗

A Variational Nodal Approach to 2D/1D Pin Resolved Neutron Transport: II Spherical Harmonics

The variational nodal method (VNM) as implemented in the VARIANT code has found wide application to fast reactor problems where fuel assembly homogenization does not result in unacceptable errors. The objective of this research is to generalize VNM for application to PWR problems, where homogenization even at the pin-cell level often results in unacceptable errors. To accomplish this, we are developing the PANX code (Purdue - Argonne - Northwestern - Xi'an) to implement a 2D/1D VNM transport method without spatial homogenization employing spherical harmonics and diffusion theory, respectively, in the radial and axial directions. In the preceding paper, we employed diffusion theory in conjunction with finite element methods to implement heterogeneous VNM nodes, in which the problem of fuel-coolant smearing at axial interfaces, found in existing 2D/1D approximations, is eliminated. In this work, we generalize the resulting VNM nodes to couple spherical harmonics in the x-y plane with axial diffusion theory. We then apply the method to the three-dimensional C5G7 benchmark problem. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Surface Reordering During Layer‐by‐Layer Growth on SrTiO 3

With the development of layer‐by‐layer growth techniques such as molecular beam epitaxy (MBE), it is now possible to construct materials with atomic‐level precision, with the layer sequence primarily determined by the growth recipe. However, materials can still restructure in the high‐temperature growth environment with surface thermodynamics strongly influencing growth behavior. Here, we demonstrate that growth on ${\rm SrTiO}_{3}$ (001), the principal platform for oxide electronics, does not occur in a layer‐by‐layer fashion but follows a more complex process in which a ${\rm TiO}_{2}$ plane continually diffuses toward the growth surface. Employing in situ synchrotron X‐ray scattering combined with ab initio thermodynamic calculations, we discover the existence of a stable ${\rm TiO}_{2}$ double‐layer structure on the pristine substrate and the occurrence of dynamic layer rearrangement during homoepitaxial growth by oxide MBE. Our findings suggest that the current methods used to precisely control ${\rm SrTiO}_{3}$ surfaces are limited, as well as our ability to dictate the composition of ultrathin films, resulting in important ramifications regarding the surface reactivity of perovskite materials grown on ${\rm SrTiO}_{3}$.

Coherent Bragg Rod Analysis↗

Phase-field modeling of the interactions between an edge dislocation and an array of obstacles

Obstacles, such as voids and precipitates, are prevalent in crystalline materials. They strengthen crystals by serving as barriers to dislocation glide. Here in this work, we develop a phase-field dislocation dynamics (PFDD) technique for investigating the interactions between dislocations and second-phase obstacles, which can be either voids or precipitates. The PFDD technique is constructed to account for elastic heterogeneity, elastic anisotropy, dissociation of the dislocation, and dislocation transmission across bicrystalline interfaces. Within the framework, we present a model for “pseudo-voids”, which are voids shearable by dislocations, in contrast to unphysical, unshearable voids in conventional phase-field dislocation formulations. We employ the PFDD technique to investigate the in-plane interactions between an edge dislocation and an array of nano-scale obstacles with different spacings. In this application, the interactions take place in glide planes of either a face-centered cubic (FCC) Cu or a body-centered cubic (BCC) Nb matrix, while the precipitates have a Cu 1-x Nb x composition, with x varying from 0.1 to 0.9. Our atomistic simulations find that the alloy precipitates can have an FCC, an amorphous, or a BCC phase, depending on the compositional ratio between Cu and Nb, i.e., value of x. Among all types of obstacles, the critical stresses for dislocation bypass are the highest for unshearable amorphous precipitates, followed by shearable crystalline precipitates, and then the pseudo-voids.

42 ENGINEERING↗

Automated Crystal Orientation Mapping in py4DSTEM using Sparse Correlation Matching

Crystalline materials used in technological applications are often complex assemblies composed of multiple phases and differently oriented grains. Robust identification of the phases and orientation relationships from these samples is crucial, but the information extracted from the diffraction condition probed by an electron beam is often incomplete. We have developed an automated crystal orientation mapping (ACOM) procedure which uses a converged electron probe to collect diffraction patterns from multiple locations across a complex sample. We provide an algorithm to determine the orientation of each diffraction pattern based on a fast sparse correlation method. We demonstrate the speed and accuracy of our method by indexing diffraction patterns generated using both kinematical and dynamical simulations. We have also measured orientation maps from an experimental dataset consisting of a complex polycrystalline twisted helical AuAgPd nanowire. From these maps we identify twin planes between adjacent grains, which may be responsible for the twisted helical structure. All of our methods are made freely available as open source code, including tutorials which can be easily adapted to perform ACOM measurements on diffraction pattern datasets.

36 MATERIALS SCIENCE↗

Effects of Particle Mixing and Gravitational Settling on Charge Transport in Carbon Flow-Electrode Cells

Flowable carbon slurries are actively studied and under development for charge transport in various electrochemical systems including flow capacitors, capacitive deionization cells, semi-solid flow batteries, and lithium extraction. However, much less is known about in operando slurry flow dynamics and their corresponding effect on charge transport. We performed an experimental study of mixing and settling dynamics of slurry electrodes within an electrochemical flow cell during continuous operations. The electrochemical cell consisted of two horizontal co-flowing channels, separated by a cation-exchange membrane (CEM). Here we used high-speed optical imaging of planes parallel to gravity and simultaneous electrochemical measurements. At low flow rates, dense yet dynamic particle beds formed on the bottom electrode in each channel, which unexpectedly yielded the highest currents. This approach enables the operation of the flow cell at low system-average particle concentrations while leveraging gravity-driven particle settling to locally enhance carbon concentrations precisely at the current collector sites. Conversely, high flow rates were characterized by thin particle beds and well-mixed particle flows. In the latter case, the electrodes in closest proximity (located on either side of the CEM) achieved a current higher than the other electrode pairs. The observations have implications for slurry control and electrode designs in electrochemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Particle-Based Conditional Sampling Scheme for the Simulation of Transport in Fractured Rock With Diffusion Into Stagnant Water and Rock Matrix

In situ experiments and field-scale characterization studies have pointed out that, in fractured crystalline media, groundwater flow is highly channelized. This implies that, at the scale of a single fracture, only part of the fracture surface area is in contact with flowing water, while the rest of in-plane water is essentially stagnant and can be accessed by solutes via molecular diffusion. Despite their importance for contaminant retention, to date, there are no numerical or analytical approaches that could be used to assess the implication of stagnant water zones on solute transport in realistic large-scale Discrete Fracture Network-based models. In this study, we present an efficient and flexible algorithm for the simulation of transport in fractured rock with diffusion into stagnant water and rock matrix. The algorithm is a generalization of a previously developed numerical framework for time domain particle tracking in sparsely fractured rock. The key of the generalization is that total time in fracture (τ f ) is first evaluated using a Monte Carlo sampling and then a second sampling is performed conditioned on τ f . The algorithm has been successfully validated against existing independent solutions and the implication of diffusion into stagnant water and secondary diffusion into the matrix has been assessed for a realistic modeling scenario. The results show that, due to diffusion into stagnant water, contaminants are more strongly retarded. This increased retention is more significant for sorbing species, as a larger number of sorption sites is accessible. A high sensitivity to the flowing channel/stagnant water zone geometry has also been observed.

54 ENVIRONMENTAL SCIENCES↗

Switching the spin cycloid in BiFeO3 with an electric field

Abstract Bismuth ferrite (BiFeO 3 ) is a multiferroic material that exhibits both ferroelectricity and canted antiferromagnetism at room temperature, making it a unique candidate in the development of electric-field controllable magnetic devices. The magnetic moments in BiFeO 3 are arranged into a spin cycloid, resulting in unique magnetic properties which are tied to the ferroelectric order. Previous understanding of this coupling has relied on average, mesoscale measurements. Using nitrogen vacancy-based diamond magnetometry, we observe the magnetic spin cycloid structure of BiFeO 3 in real space. This structure is magnetoelectrically coupled through symmetry to the ferroelectric polarization and this relationship is maintained through electric field switching. Through a combination of in-plane and out-of-plane electrical switching, coupled with ab initio studies, we have discovered that the epitaxy from the substrate imposes a magnetoelastic anisotropy on the spin cycloid, which establishes preferred cycloid propagation directions. The energy landscape of the cycloid is shaped by both the ferroelectric degree of freedom and strain-induced anisotropy, restricting the spin spiral propagation vector to changes to specific switching events.

36 MATERIALS SCIENCE↗

The next-generation magnetic recoil spectrometer (MRSnext) on OMEGA and NIF for diagnosing ion temperature, yield, areal density, and alpha heating

The next-generation magnetic recoil spectrometer (MRSnext) is being designed to replace the current MRS at the National Ignition Facility and OMEGA for measurements of the neutron spectrum from an inertial confinement fusion implosion. The MRSnext will provide a far-superior performance and faster data turnaround than the current MRS systems, i.e., a 2× and 6× improvement in energy resolution at the NIF and OMEGA, respectively, and 20× improvement in data turnaround time. The substantially improved performance of the MRSnext is enabled by using electromagnets that provide a short focal plane (12–16 cm) and unprecedented flexibility for a wide range of applications. In addition to being able to measure neutron yield, apparent ion temperature, areal density, and plasma-flow velocity over a wide range of yields, the NIF MRSnext will be able to directly, uniquely assess the alpha heating of the fuel ions through measurements of the alpha knock-on tail in the neutron spectrum. The goal is to implement a radiation-hard electronic detection system capable of providing rapid data acquisition and analysis. The development of the MRSnext will also set the foundation for the more advanced, time-resolving MRSt and serve as a testbed for its implementation on the NIF.

47 OTHER INSTRUMENTATION↗

THE THOMSON SURFACE. II. POLARIZATION

The solar corona and heliosphere are visible via sunlight that is Thomson-scattered off free electrons, yielding a radiance against the celestial sphere. In this second part of a three-article series, we discuss linear polarization of this scattered light parallel and perpendicular to the plane of scatter in the context of heliospheric imaging far from the Sun. The difference between these two radiances (pB) varies quite differently with scattering angle, compared to the sum that would be detected in unpolarized light (B). The difference between these two quantities has long been used in a coronagraphic context for background subtraction and to extract some three-dimensional information about the corona; we explore how these effects differ in the wider-field heliospheric imaging case where small-angle approximations do not apply. We develop an appropriately simplified theory of polarized Thomson scattering in the heliosphere, discuss signal-to-noise considerations, invert the scattering equations analytically to solve the three-dimensional object location problem for small objects, discuss exploiting polarization for background subtraction, and generate simple forward models of several classes of heliospheric feature. We conclude that pB measurements of heliospheric material are much more localized to the Thomson surface than are B measurements, that the ratio pB/B can be used to track solar wind features in three dimensions for scientific and space weather applications better in the heliosphere than corona, and that,as an independent measurement of background signal, pB measurements may be used to reduce the effect of background radiances including the stably polarized zodiacal light.

79 ASTRONOMY AND ASTROPHYSICS↗

Diverted Negative Triangularity plasmas on DIII–D: The benefit of high confinement without the liability of an edge pedestal

Diverted discharges at negative triangularity on the DIII-D tokamak sustain normalized confinement and pressure levels typical of standard H-mode scenarios (H 98y,2 ≅1, β N ≅3) without developing an edge pressure pedestal, despite the auxiliary power far exceeding the L → H power threshold expected from conventional scaling laws. The power degradation of confinement is substantially weaker than the ITER-89P scaling, resulting in a confinement factor that improves with increasing auxiliary power. Furthermore, the absence of the edge pedestal is beneficial in several aspects, such as eliminating the need for active mitigation or suppression of edge localized modes, low impurity retention and a reconstructed scrape-off layer heat flux width at the mid-plane that exceeds the ITPA multi-machine scaling law by up to 50%. Together with technological advantages granted by placing the divertor at larger radii, plasmas at Negative Triangularity without an edge pedestal feature both core confinement and power handling characteristics that are potentially suitable for operation in future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental and Monte Carlo studies of fluence corrections for graphite calorimetry in low- and high-energy clinical proton beams

Purpose: The aim of this study was to determine fluence corrections necessary to convert absorbed dose to graphite, measured by graphite calorimetry, to absorbed dose to water. Fluence corrections were obtained from experiments and Monte Carlo simulations in low- and high-energy proton beams. Methods: Fluence corrections were calculated to account for the difference in fluence between water and graphite at equivalent depths. Measurements were performed with narrow proton beams. Plane-parallel-plate ionization chambers with a large collecting area compared to the beam diameter were used to intercept the whole beam. High- and low-energy proton beams were provided by a scanning and double scattering delivery system, respectively. A mathematical formalism was established to relate fluence corrections derived from Monte Carlo simulations, using the FLUKA code [A. Ferrari et al., “FLUKA: A multi-particle transport code,” in CERN 2005-10, INFN/TC 05/11, SLAC-R-773 (2005) and T. T. Böhlen et al., “The FLUKA Code: Developments and challenges for high energy and medical applications,” Nucl. Data Sheets 120, 211–214 (2014)], to partial fluence corrections measured experimentally. Results: A good agreement was found between the partial fluence corrections derived by Monte Carlo simulations and those determined experimentally. For a high-energy beam of 180 MeV, the fluence corrections from Monte Carlo simulations were found to increase from 0.99 to 1.04 with depth. In the case of a low-energy beam of 60 MeV, the magnitude of fluence corrections was approximately 0.99 at all depths when calculated in the sensitive area of the chamber used in the experiments. Fluence correction calculations were also performed for a larger area and found to increase from 0.99 at the surface to 1.01 at greater depths. Conclusions: Fluence corrections obtained experimentally are partial fluence corrections because they account for differences in the primary and part of the secondary particle fluence. A correction factor, F(d), has been established to relate fluence corrections defined theoretically to partial fluence corrections derived experimentally. The findings presented here are also relevant to water and tissue-equivalent-plastic materials given their carbon content.

60 APPLIED LIFE SCIENCES↗

Improved Regional Moment Tensor Inversion for Moderately Large Earthquakes in the Western United States Using a 3D Earth Model Based on Full Waveform Tomography

The nature of seismic sources for moderately large (moment magnitude, M w 5.0–6.5) events are commonly characterized by their moment tensor (MT) solutions and obtained by inversion of regional distance (200–1600 km) long‐period (20–50 s) waveforms. Regional MT estimates are often calculated from average plane‐layered, one‐dimensional (1D) velocity models. However, 1D model calculations can produce misfits in the arrival times and waveform shapes that introduce errors, particularly at longer distances or for shorter periods, which are necessary for analyzing lower magnitude events. Approximate Earth models (e.g., 1D) representing broad areas may be inadequate, particularly in the crust and uppermost mantle of tectonically complex regions. In this study, we show how a three‐dimensional (3D) Earth model obtained from full waveform inversion tomography can improve waveform fits and decrease phase errors. We developed a platform and workflow to perform routine 3D MT inversions and inverted MTs for 25 earthquakes in the western United States and seven nuclear explosions using an average 1D and a recent 3D Earth model, WUS256 (Rodgers et al., 2022). Using the 3D model improves waveform fits (variance reduction and phase time shifts) compared with the 1D model, and the 3D MT solutions are stable across large distances. This study shows that 3D models obtained from full waveform tomography can improve MTs and source characterization especially at far regional distances (>800 km).

Geosciences↗

The Measurement of Electrical and Thermal Conductivities in Warm Dense Matter

This document serves as the Final Technical Report for grant DE-SC0012333, The Measurement of Electrical and Thermal Conductivities in Warm Dense Matter . This work proposed to make the first systematic measurements of the reflectivity of warm dense matter (WDM) at multiple wavelengths by controlling and measuring the experimental factors necessary to use the reflectivity measurements to determine the electrical and thermal conductivity as a function of temperature and, to a lesser degree, density. A new technique was developed to measure the real and imaginary components of the electrical conductivity by making only scalar measurements of the relative reflectivity different angles. A series of preliminary experiments were designed to use the second, third, and fourth harmonics of the 800 nm driving laser. Reflectivity at each harmonic for two different angles would provide the scalar information needed to extract the real and imaginary parts of the index of refraction. Experiments were made using continuous wave lasers to determine the limits of accuracy on the measured reflectance for surrogate plane mirrors. It was determined, that while the technique was sound, the practical application of it was limited for the conditions on Scarlet. Due to significant challenges with the Scarlet Laser Facility at the time, progress was delayed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Blunt Impact Brain Injury using Cellular Injury Criterion

The Advanced Combat Helmet (\ACH") military specification (\mil-spec") requires a helmeted magnesium (\Mg") Department of Transportation (\DOT") headform be dropped vertically, with an impact speed of 3.1 m/s (10 ft/s), onto a steel hemispherical target. The pass/fail criteria are based on translational acceleration (150 G) alone, absent of any rotational component. Without a rotational component, the specification's injury risk application is limited to skull fracture and peripheral hematomas (subdural, subarachnoid), since this translational acceleration injury risk assessment is based on the Wayne State Tolerance Curve (\WSTC"). To provide a more comprehensive view of injury for the entire brain, an alternative approach is needed. To meet this need, we worked with a larger group called PANTHER, a collaboration between national laboratories, industry, and academia. Collaborations specific to research and results presented here come from efforts led by Mr. Ron Szalkowski and Mr. Sushant Malave, Ms. Alice Fawzi, and Dr. Christian Franck. We have developed a prototypical injury risk criterion based on the neuronal response to abrupt changes in general motion (translation, rotation, or both). The cellular-based mild traumatic brain injury (\cbmTBI") criterion utilizes both the strain and strain rate of brain tissue to account for the stretch and rate of stretch that occurs throughout the brain as a result of blunt impact to the head. We conducted physical experiments of an ACH-fitted magnesium headform, which produced repeatable headform peak accelerations. Then, we developed a simulation of the experiment, and validated the simulation output with the experimental data. We then substituted the magnesium headform with a human headform, consisting of skin, muscle, bone, gray matter, white matter, cerebral-spinal fluid, membranes, vasculature, intravertebral discs, airway and sinus. We quantified brain injury risk using the cbmTBI criterion, using the current mil-spec test and a modified test. The modified mil-spec test used an inclined anvil target that was located posterior to the crown of the helmet in the axial plane. While the current mil-spec test produced brain deformation from head translation alone, the modified test produced brain deformation from head translation and rotation, which is closer to most real world and combat theater impacts (e.g., such as occur in tertiary blast exposure). Compared to the current mil-spec test, the modified test produced elevated strains in the human digital twin. These data, mapped to the cbmTBI criterion, suggest increased injury risk for blunt impacts that cause rotation and translation, rather than just translation alone. Moreover, these data may lead to a rotational performance metric, which is rooted in the biology and pathology of the brain's response to impact and blast, and which should be used to improve next-generation helmet designs.

42 ENGINEERING↗

Growth Oscillatory Zoning in Erythrite, Ideally Co3(AsO4)2·8H2O: Structural Variations in Vivianite-Group Minerals

The crystal structure of an oscillatory zoned erythrite sample from Aghbar mine, Bou Azzer, Morocco, was refined using synchrotron high-resolution powder X-ray diffraction (HRPXRD) data, Rietveld refinement, space group C2/m, and Z = 2. The crystal contains two sets of oscillatory zones that appear to have developed during epitaxial growth. The unit-cell parameters obtained are a = 10.24799(3) Å, b = 13.42490(7) Å, c = 4.755885(8) Å, β = 105.1116(3)°, and V = 631.680(4) Å3. The empirical formula for erythrite, obtained with electron-probe micro-analysis (EPMA), is [Co2.78Zn0.11Ni0.07Fe0.04]Σ3.00(AsO4)2·8H2O. Erythrite belongs to the vivianite-type structure that contains M1O2(H2O)4 octahedra and M22O6(H2O)4 octahedral dimers that are linked by TO4 (T5+ = As or P) tetrahedra to form complex layers parallel to the (010) plane. These layers are connected by hydrogen bonds. The average [6] = 2.122(1) Å and average [6] = 2.088(1) Å. With space group C2/m, there are two solid solutions: M3(AsO4)2·8H2O and M3(PO4)2·8H2O where M2+ = Mg, Fe, Co, Ni, or Zn. In these As- and P-series, using data from this study and from the literature, we find that their structural parameters evolve linearly with V and in a nearly parallel manner despite of the large difference in size between P5+ (0.170 Å) and As5+ (0.355 Å) cations. Average [4], [6], and [6] distances increase linearly with V. The average distance is affected by M atoms, whereas the average distance is unaffected because it contains shorter and stronger P–O bonds. Although As- and P-series occur naturally, there is no structural reason why similar V-series vivianite-group minerals do not occur naturally or cannot be synthesized.

36 MATERIALS SCIENCE↗

Enhanced Spontaneous Polarization in Ultrathin SnTe Films with Layered Antipolar Structure

2D SnTe films with a thickness of as little as 2 atomic layers (ALs) have recently been shown to be ferroelectric with in–plane polarization. Remarkably, they exhibit transition temperatures (T c ) much higher than that of bulk SnTe. Here, combining molecular beam epitaxy, variable temperature scanning tunneling microscopy, and ab initio calculations, the underlying mechanism of the Tc enhancement is unveiled, which relies on the formation of γ–SnTe, a van der Waals orthorhombic phase with antipolar inter–layer coupling in few–AL thick SnTe films. In this phase, 4n – 2 AL (n = 1, 2, 3…) thick films are found to possess finite in–plane polarization (space group Pmn2 1 ), while 4n AL thick films have zero total polarization (space group Pnma). Above 8 AL, the γ–SnTe phase becomes metastable, and can convert irreversibly to the bulk rock salt phase as the temperature is increased. This finding unambiguously bridges experiments on ultrathin SnTe films with predictions of robust ferroelectricity in GeS–type monochalcogenide monolayers. Here, the observed high transition temperature, together with the strong spin–orbit coupling and van der Waals structure, underlines the potential of atomically thin γ–SnTe films for the development of novel spontaneous polarization–based devices.

2d ferroelectrics↗