Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boundary layer thickness”

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

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

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

Constraining the Aerosol Effects on Deep Convective Clouds by Considering the Coupling Between Clouds and the Planetary Boundary Layer

Several mechanisms have been proposed for the aerosol invigoration effect. Although their principles are well established, their actual magnitudes and roles in cloud development remain uncertain and debatable. This uncertainty partly stems from observational-based studies, in which it has been challenging to separate the co-variability between aerosols and meteorology. Addressing this problem requires large data samples. To this end, this study employs the Atmospheric Radiation Measurement data set expanding to 16 years (some up to 17 years, compared to 10 years in previous work) in the U.S. Southern Great Plains. It also conducts a more careful and rigorous analysis to isolate the influences of convective available potential energy (CAPE) and synoptic patterns to address a previously raised concern. We incorporated a new key process affecting aerosol-cloud interaction: cloud-surface coupling. The state/degree of the coupling relationship turns out to play an important role in the invigoration effect. Our analysis reinforces earlier findings of a robust positive relationship between cloud thickness and aerosol loading across CAPE percentiles—but only under cloud-surface coupled conditions. The increase in cloud thickness with aerosol loading is most pronounced in coupled clouds with high CAPE and bases below 1 km. Coupled clouds with bases below 1 km thicken between 1 and 4 km, depending on the CAPE percentile. Decoupled clouds show no such systematic changes. Synoptic patterns also lead to different strengths of the invigoration effect. Clean and polluted air masses are predominantly associated with northerly and southerly winds, respectively, with a stronger invigoration effect in cleaner air masses.

Geosciences↗
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Estimates of global recycling coefficients for LTX- β discharges

We report the first observation of global recycling coefficient R near 0.5 in the Lithium Tokamak eXperiment-β (LTX-β), significantly below the minimum R previously reported in other devices. In a series of experiments with varied Li wall conditioning, estimates of the recycling coefficient have been made using a Lyman-α array and DEGAS2 modeling. A progressive reduction in Lyman-α emission with increased lithium and an increase in edge electron temperature are observed. It is also observed that with increasing Li coating thickness, the effective particle confinement time τp* is reduced and approaches TRANSP calculated energy confinement time (τE), with τp* near τE,TRANSP for the lowest recycling coefficients. Edge temperatures approaching core plasma temperatures, first reported in LTX, can now be directly connected to estimates of the recycling coefficient and qualitatively agree with previous UEDGE simulations. The particle flux to the limiting surfaces appears to be significantly reduced in comparison with fluid scrape-off layer (SOL) models, indicating that a large fraction of the SOL ions are mirror trapped. SOL collisionality drops more than an order of magnitude below the banana regime boundary, indicating the importance of kinetic effects. Full-f 1x2v gyrokinetic simulations of SOL field lines with the GKEYLL code indicate that the fraction of ions trapped along field lines increases as collisionality drops, as a result of increased lithium evaporation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Reexamination of hydrodynamic phonon transport in thin graphite

The recent hydrodynamic phonon transport theory for graphitic materials has been supported by the measurements of the second sound at temperatures up to about 100 K. When boundary scattering becomes comparable to momentum-conserving normal phonon scattering processes that are responsible for phonon hydrodynamics, Poiseuille phonon flow phenomena can emerge to give rise to unique size-dependent thermal conductivity in thin graphite. Here, we examine the thickness range for the Poiseuille phonon flow to become observable in thin graphite with the use of both deviational Monte Carlo simulation of the Peierls-Boltzmann transport equation and four-probe thermal transport measurements. As the basal-plane thermal conductivity calculated by prior first-principles theories saturates to the graphite value when the thickness is increased to five graphene layers, the phonon dispersion of graphite is used in the current calculations of thin graphite of micrometer thickness and a 23-layer thick ultrathin graphite (UTG) sample. The calculations show that diffuse surface scattering by surface defects can lead to Poiseuille phonon flow at 50 K in thin graphite with the thickness close to several micrometers but not in the 65 μm thin graphite and 23-layer UTG, where phonon scattering by the top and bottom surfaces become, respectively, much less and more frequent than the normal processes. In addition, the calculation results with the bulk graphite dispersion and diffuse surface scattering show decreased basal-plane thermal conductivity with decreasing thickness, opposite to recent thermocouple measurements of thin graphite samples. In comparison, the calculation results reveal that partially diffuse surface defect scattering can yield the four-probe measurement results of UTG samples, which are prepared here with an improved process to minimize surface contamination.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ferroelectric Al1 − xBxN–GaN heterostructures

This report demonstrates Al0.93B0.7N thin films grown epitaxially on n-type GaN (0002)/c-plane sapphire substrates by reactive magnetron sputtering at 300 °C. At 200 nm film thickness, the Al0.93B0.07N layers exhibit partially relaxed substrate-induced epitaxial strain, a 0.16° wide (0002) rocking curve, in-plane crystallographic registry, and sub-nanometer surface roughness. Electrically, the stack shows robust hysteresis over three frequency decades, a remanent polarization of ∼125 μC/cm2, a strongly frequency dependent coercive field, highly uniform dc leakage currents, and endurance >106 field cycles. This report validates possibilities for ferroelectric nitride integration into conventional III-nitride heterostructures with high crystalline fidelity, high electrical resistivity, and persistent hysteresis. Such materials are synthesizable at thermal budgets and temperatures compatible with back-end-of-the-line boundary conditions.

Physics↗
DOE OSTI2024

Stress localization investigation of additively manufactured GRCop-42 thin-wall structure

A full-field crystal plasticity (CP) framework is presented for the GRCop-42 alloy to study microscopic mechanical behavior and local stress heterogeneities. The microstructures of additively manufactured (AM) materials are often unique relative to conventionally processed materials, and the local thermal histories drive these differences during the build process. These thermal histories depend on the process parameters (laser power, scan speed, and scan strategy) and the part geometry. Prior research has shown that the mechanical properties of thin-walled structures can vary significantly with wall thickness due to changes in the thermal boundary conditions during manufacturing. It is, therefore, desirable to perform CP simulations based on the phenomenological constitutive model to predict the local mechanical responses induced by microstructural heterogeneities. This work generates representative microstructures based on experimentally collected grain information (i.e., texture) for grain scale stress analysis, and the material constitutive parameters are calibrated using the experimental mechanical testing data. Here, we specifically investigated the effect of crystallographic texture and grain morphologies on the size-dependent mechanical properties of AM GRCop-42. The selection of appropriate material properties for implementing an effective free surface boundary condition and the influence of adjacent buffer layers are also discussed. Analysis of local field results reveals a strong correlation between stress localization and the initial grain orientation. However, no significant relationship between the misorientation of the individual adjacent grains and the average misorientation is observed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Interference-Limited Absorption in Dense Molecular Nanolayers Near Reflecting Surfaces

We investigate linear resonant absorption by a dense ensemble of molecules confined to a sub-wavelength layer in two geometries: (i) a free-standing film in a homogeneous space and (ii) the same film placed at a controlled distance from a reflecting surface. In both cases, increasing the effective light–matter coupling (via molecular density/oscillator strength) produces a nonmonotonic response: absorption rises to an optimum and then decreases as the film becomes increasingly radiatively bright and reflective. Finite-difference time-domain simulations and analytical transfer-matrix calculations agree quantitatively and yield compact ridge conditions for the optimum. We interpret the trends using a scattering/port picture: the isolated film is a symmetric two-port system (reflection and transmission), which bounds single-sided resonant absorption to ≤50% in the ultrathin limit (reflecting transition saturation), whereas adding a mirror suppresses transmission and converts the structure into an effectively one-port absorber. In the mirror-backed geometry, interference can cancel reflection, and unity absorption is obtained at critical coupling, when radiative leakage is balanced by intrinsic molecular loss. These results clarify fundamental limits and design rules for collective absorption in dense molecular layers near dielectric or metallic boundaries.

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

Graphene-polyelectrolyte multilayer membranes with tunable structure and internal charge

One great advantage of graphene-polyelectrolyte multilayer (GPM) membranes is their tunable structure and internal charge for improved separation performance. In this study, we synthesized GO-dominant GPM membrane with internal negatively-charged domains, polyethyleneimine (PEI)-dominant GPM membrane with internal positively-charged domains and charge-balanced dense/loose GPM membranes by simply adjusting the ionic strength and pH of the GO and PEI solutions used in layer-by-layer membrane synthesis. A combined system of quartz crystal microbalance with dissipation (QCM-D) and ellipsometry was used to analyze the mass deposition, film thickness, and layer density of the GPM membranes. The performance of the GPM membranes were compared in terms of both permeability and selectivity to determine the optimal membrane structure and synthesis strategy. One effective strategy to improve the GPM membrane permeability-selectivity tradeoff is to assemble charge-balanced dense membranes under weak electrostatic interactions. This balanced membrane exhibits the highest MgCl 2 selectivity (~86%). Another effective strategy for improved cation removal is to create PEI-dominant membranes that provide internal positively-charged barrier to enhance cation selectivity without sacrificing water permeability. Finally, these findings shine lights on the development of a systematic approach to push the boundary of permeability-selectivity tradeoff for GPM membranes.

42 ENGINEERING↗
DOE OSTI2026

A large interlaboratory electron diffraction study of monolayer graphene

Standardisation of data collection and analysis is essential to enable commercialisation of 2D materials in a wide range of technologies. Selected area electron diffraction (SAED) in the transmission electron microscope (TEM) is one of the key methods for distinguishing monolayer from bilayer and few-layer graphene by comparing the 1st and 2nd order diffraction spot intensities. Yet there are many factors that can affect the reliability of data collection and interpretation, causing the measurement of monolayer samples to deviate from the literature boundary condition of $I_{\{\bar{2}110\}}$$/$$I_{\{1\bar{1}00\}}$ < 1 for monolayer graphene (1LG). Here we present the results of a large interlaboratory SAED comparison study, where 15 international laboratories measured and analysed nominally identical samples of chemical vapour deposited graphene. Large variations were observed in the measured ratios of diffraction spot intensities, with the largest variance associated with poor quality SAED data resulting from inadequate specimen handling and storage. To inform the reliable determination of monolayer thickness from SAED patterns we provide a description of best practice for specimen handling, TEM operation, data collection and analysis. This work was undertaken within VAMAS Technical Working Area 41: Graphene and related 2D materials—Project 9, the results of which have been directly incorporated into ISO/TS 21356–2 for the characterisation of graphene sheets. We find that when this methodology is followed, 1LG can be distinguished from bilayer or thicker material with high confidence where analysis of a single SAED pattern gives $I_{\{\bar{2}110\}}$$/$$I_{\{1\bar{1}00\}}$ < 1.2, even in the absence of precise specimen tilting.

2D materials↗
DOE OSTI2025

Design and Analysis of the Open-Surface Slow Li Flow Divertor and Comparison to the Fast Li Flow Divertor

In the ongoing U.S. project, “Liquid Metal Plasma Facing Components,” sponsored by the U.S. Department of Energy, efforts have been taken to develop two open-surface divertor designs for the Fusion Nuclear Science Facility using liquid lithium (Li) as a heat and particle flux removal media. The main focus of this study is the design and analysis of a slow (~1 mm/s) and thin (<1 mm) open-surface Li flow divertor with a Li-cooled substrate, which is then compared with an earlier design of a fast (up to 10 m/s) and thick (~0.5 cm) Li flow divertor with the substrate cooled with helium. The slow Li flow divertor design is based on the original LiWall concept developed at the Princeton Plasma Physics Laboratory. Such a thin and slow Li layer can remove the particle flux by reducing the recycling flux, while the heat flux is removed mainly through the heat sink located beneath. In the present study, the heat sink is provided through a Li cooling flow inside the substrate of reduced activation ferritic/martensitic steel. By performing a multiphysics analysis with COMSOL that included liquid-metal magnetohydrodynamics (MHD), heat transfer, and structural mechanics, the impact of various factors on the divertor heat removal capability, such as Li flow velocity, MHD effects, and inlet velocity boundary condition, were examined. Based on comparisons of the two divertor designs, it was shown that the fast-flow divertor significantly outperformed the slow-flow design, whose heat removal capability was limited to ~1 to 2 MW/m2.

Jiang, Yuchen [ORNL] (ORCID:0000000283280284)↗
DOE OSTI2020

The Role of Secondary Particle Structures in Surface Phase Transitions of Ni-Rich Cathodes

Nanometer-thick reconstruction layers on layered cathode surfaces have been widely observed on both pristine and cycled materials. However, the mechanisms of reconstruction and the role that these structures play in electrochemical performance are not fully understood. From a crystallographic perspective, such surface reconstruction layers result from cation site-mixing between Li and transition-metal ions, but it remains far from clear as to what the critical factors are that control such cation mixing for various cathode chemistries. In this paper, we report observations on surface transformations during cycling of a Ni-rich cathode-oxide that are governed by direct contact with the liquid electrolyte. Specifically, within a secondary particle, the three-dimensional hierarchical aggregate of primary particles leads to the formation of grain boundaries of different characteristics, including open grain boundaries that allow the permeation of liquid electrolyte, and closed grain boundaries that exclude liquid electrolyte penetration. Upon battery cycling, surfaces in direct contact with the liquid electrolyte showed a surface reconstruction layer while the surfaces that did not have contact with the electrolyte showed no such reconstructions. In addition, the critical role of oxygen in reconstruction of grain boundaries is demonstrated, and the unexpected oxygen depletion at twin boundaries result in severe phase transitions. The present work provides further insights into phase transitions at solid-solid and solid-liquid interfaces within secondary cathode particles and suggests strategies for mitigating surface/interfacial degradation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Suppression of electroconvective and morphological instabilities by an imposed cross flow of the electrolyte

Electroconvective (EC) instability and its influence on surface morphological perturbations (Morph) are important in many applications, including electrodialysis, batteries, and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the EC and Morph instabilities using two approaches. In the bulk analysis, we derive the asymptotic solutions for small and large wave numbers by neglecting the space charge layer and imposing a second kind electroosmosis slip velocity boundary condition on the electroneutral bulk region. In the full analysis, the instability of the entire region of the liquid electrolyte is analyzed using the ultraspherical spectral method. Both studies show that the flow significantly affects the EC instability. The imposed flow distorts the concentration field, causing a sheltering effect which hinders the ion transport from low- to high-concentration regions, and therefore suppresses the EC instability below a certain wave number. In combination with the viscous stabilization of the high wave number modes by the space charge layer, a sufficiently strong imposed flow can fully suppress the EC instability. The increment in the critical voltage for the instability onset is roughly proportional to the square root of the product of imposed velocity and double layer thickness. The imposed flow has a smaller effect on the Morph instability, except that it may remove the morphological modes resulting from the EC instability and thereby change the wave number of the most unstable mode.

42 ENGINEERING↗
DOE OSTI2021

Exploring Sources of Surface Bias in HRRR Using New York State Mesonet

In recent years, there has been increasing demand for applications of short-term forecasting of renewable energy potential and assessments of the likelihood of extreme weather events using the High-Resolution Rapid Refresh (HRRR) model. Examining the biases in the newest version of HRRR is necessary to promote further model development. Using data from one of the most comprehensive and dense monitoring networks, New York State Mesonet (NYSM), we evaluate the HRRR version 3 meteorological fields for an entire year. In this work, the land-atmosphere-cloud coupling system is evaluated as an integrated whole. We investigate the physical processes influencing the soil hydrological balance and the thermodynamic interactions, from surface fluxes up to the level of boundary layer convection from both temporal (seasonal and diurnal) and spatial perspectives. Results show that the model 2 m temperature and humidity biases are seasonally dependent, with warm and dry bias present during the warm season, and an extreme nocturnal cold bias in winter. The summer warm bias includes both a land-surface-induced bias and a cloud-induced bias. Inaccurate representation of energy partition and soil hydrological process across different land use types as well as a hydrological bias in describing spring snowmelt are identified as the main source of the land-surface-induced bias. A feedback loop linking cloud presence, flux changes, and temperature contributes to the cloud-induced bias. The positive solar radiation bias increases from clear sky to overcast sky conditions. The most significant bias occurs during overcast and thick cloud conditions associated with frontal passage and thunderstorms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Modelling dislocation-graphene interactions in a BCC Fe matrix by molecular dynamics simulations and gradient plasticity theory

We report graphene nanosheets (GNS) can enhance the strength and ductility of metal-based composites as they can obstruct the propagation of dislocations. The present article employs Molecular Dynamics (MD) simulations to investigate dislocation-GNS interaction mechanisms and possible influencing factors, including the number of GNS layers, the thickness of the metallic amorphous layer and the C - C bond strength. The results indicated that the shear strength of the metal/GNS interface and the bending stiffness of GNS determined the ability of GNS to block dislocation transmission. A physically based phenomenological parameter that can capture such dislocation-GNS interactions is the mechanical interface energy that has been put forth within gradient plasticity. By fitting the theoretical expressions to the simulation data, it was possible to obtain estimates for the mechanical interface energy for the GNS. It was found that increasing the GNS layers and adding an amorphous layer resulted in a strengthening in the stress–strain response and increased the value of this interfacial parameter. This indicates that the mechanical interfacial energy can be a unified measure for capturing and tuning the strength of various interfaces such as grain boundaries, GNS, amorphous-crystalline interface and bimetal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Surface oxidation and crystallization of FeNi-Based soft magnetic nanocrystalline and amorphous nanocomposite alloys

In this work, we present characterization of surface oxidation and crystallization behavior of Metal Amorphous Nanocomposite (MANC) (Fe 70 Ni 30 ) 80 Nb 4 Si 2 B 14 alloys via complementary x-ray diffraction, electron microscopy, and x-ray photoelectron spectroscopy techniques. An adherent amorphous oxide of roughly 15–20 nm in thickness is observed on the wheel side of the ribbon, enriched in Fe relative to Ni and with significant content of Si and Nb as would be expected based upon a comparison of the Gibbs free energy of formation of the stable oxide phases. The amorphous character of the oxide is likely to result from the significant Si-content coupled with the amorphous and nanocrystalline character of the underlying metallic layer. Higher density of FeNi-enriched nanocrystals is also observed within ~100 nm of the wheel-side surface, presumably due to reduction in local glass former concentration. A rough, crystalline copper oxide layer is also present and forms through contact between the copper wheel and the solidifying melt during the rapid solidification process. On the air side of the ribbon, a crystalline oxide layer of ~15–20 nm is formed as well as intergranular oxidation that penetrates along the grain boundaries between highly textured surface crystallites of ~50–100 nm in size. The crystalline nature of the oxide on these surfaces is asserted to result from the lack of significant Si-content and the pre-existing surface crystallites that form during the rapid solidification process which promote a crystalline oxide surface layer. Furthermore, these results shed new insights into the details of surface oxidation and can help to inform processing and alloy composition adjustments which seek to tailor surface properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Direct Imaging of the Structural and Morphological Evolution of Epitaxial LiCoO 2 Films during Charge and Overcharge

The capacity decay of layered cathodes in high-voltage applications underscores the need to utilize accurate and precise techniques to understand the underlying mechanisms. Here, we use well-defined epitaxial LiCoO 2 (LCO) films on SrRuO 3 /SrTiO 3 (SRO/STO) with controlled orientations and defect structures along with in situ electrochemical atomic force microscopy to probe the structural and morphological evolutions during the charge and overcharge processes. We quantitatively show the morphological changes in both the reversible delithiation regime and the irreversible over-delithiation regime and correlate the overall electrochemical behaviors to atomic scale defect evolutions in the films. Further, we also observe a significantly lower charging capacity for LCO/SRO/STO(111) compared to that of LCO/SRO/STO(001) films of the same thickness, which is ascribed to different types of atomic scale defects formed during the film growth process. Our high- resolution scanning transmission electron microscopy (STEM) and electron energy loss spectroscopy (EELS) studies reveal that the antiphase boundaries in LCO/SRO/STO(111) act as viable channels for Li migration but are more susceptible to irreversible phase transitions, which then block subsequent Li diffusion. The failure mechanisms developed here may provide insight into the design of future cathode materials.

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

Whole-lithosphere shear during oblique rifting

Processes controlling the formation of continental whole-lithosphere shear zones are debated, but their existence requires that the lithosphere is mechanically coupled from base to top. We document the formation of a dextral, whole-lithosphere shear zone in the Death Valley region (DVR), southwest United States. Dextral deflections of depth gradients in the lithosphere-asthenosphere boundary and Moho are stacked vertically, defining a 20–50-km-wide, lower lithospheric shear zone with ~60 km of shear. These deflections underlie an upper-crustal fault zone that accrued ~60 km of dextral slip since ca. 8–7 Ma, when we infer that whole-lithosphere shear began. This dextral offset is less than net dextral offset on the upper-crustal fault zone (~90 km, ca. 13–0 Ma) and total upper-crustal extension (~250 km, ca. 16–0 Ma). We show that, before ca. 8–7 Ma, weak middle crust decoupled upper-crustal deformation from deformation in the lower crust and mantle lithosphere. Between 16 and 7 Ma, detachment slip thinned, uplifted, cooled, and thus strengthened the middle crust, which is exposed in metamorphic core complexes collocated with the whole-lithosphere shear zone. Midcrustal strengthening coupled the layered lithosphere vertically and therefore enabled whole-lithosphere dextral shear. Where thick crust exists (as in pre–16 Ma DVR), midcrustal strengthening is probably a necessary condition for whole-lithosphere shear.

58 GEOSCIENCES↗
DOE OSTI2023

Experimental and Theoretical Studies of the Surface Oxidation Process of Rare–Earth Tritellurides

Recent studies have established Van der Waals (vdW) layered and 2D rare-earth tritellurides (RTe 3 ) as superconductors and near room-temperature charge density wave (CDW) materials. Their environmental stability raises natural concern owing to aging/stability effects observed in other tellurium-based layered crystals. Here, the results establish the stability and environmental aging characteristics of these RTe 3 systems involving a variety of metals such as La, Nd, Sm, Gd, Dy, and Ho. The atomic force microscopy (AFM) and scanning electron microscopy (SEM) results show that all the RTe 3 sheets oxidize to form thin TeO x layers that are primarily confined to the surface, edges, and grain boundaries. Time-resolved in situ Raman spectroscopy measurements are used to understand the kinetics of the oxidization process for different lanthanide metal cations and establish their relative stability/resilience to oxidization. Overall results indicate that the vdW layers show higher air stability as the 4f electron number decreases going from Ho to La, resulting in the most stable LaTe 3 compared to the least stable HoTe 3 . Comprehensive quantum mechanical simulations reveal that environmental degradation originates from a strong oxidizing reaction with O 2 molecules, while humidity (H 2 O) plays a negligible role unless Te vacancies are present. Moreover, the simulations explain the effects of 4f electrons on the work function and Te vacancies formation, which directly impact the aging characteristics of RTe 3 layers. Interestingly, optical and electrical measurements show that the CDW response is still observed in aged RTe 3 layers owing to the presence of underlying pristine/nonoxidized RTe 3 layers, except CDW transition temperatures increase due to the thickness effect. Overall results offer the first in-depth environmental aging studies on these materials, which can be applied to engineer and design their chemical stability, surface properties, and overall CDW characteristics.

36 MATERIALS SCIENCE↗