Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface energy”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Totally chlorine-free peracetic acid pulping for nanocellulose isolation from hemp and poplar

Nanocellulose is a promising and sustainable feedstock for developing advanced and functional materials. However, the characteristics of nanocellulose, such as crystallinity, surface energy, and aspect ratio, can vary depending on biomass source and pretreatment methods, leading to variable performance of the nanocellu-lose-based materials. In this study, cellulose nanocrystals (CNCs) were isolated from hemp and poplar using totally chlorine free (TCF) peracetic acid and sodium chlorite delignification and bleaching pretreatments to probe the influences of biomass source and treatment methods on the isolation and characteristics of CNCs. Our results showed that hemp and poplar were almost completely delignified by peracetic acid treatment, whereas sodium chlorite treatment left 5%–6% lignin in the pulp. The yields of CNCs from raw hemp and poplar biomass ranged from 9.8% to 21.9% and 10.9% to 28.3%, respectively, depending on the treatment methods. The dimensions of CNCs from TCF-treated biomass generally maintained a larger width and aspect ratio than those from sodium chlorite-treated biomass. The poplar-derived CNCs exhibited slightly higher crystallinity of 53%–58% than hemp-derived CNCs of 49%–54%. The zeta potential of the CNCs, ranging from -20.1 mV to -31.1 mV, ensured a well-dispersed aqueous solution. The surface energy (dispersive energy of 40–80 mJ/m2 and specific energy of 2–10 mJ/m2), water interaction, and thermal stability of the CNCs were comparable, regardless of the biomass source and pretreatment methods. Our finding suggests that the TCF technique with peracetic acid treatment is a promising delignification and bleaching approach to obtain cellulose-rich pulps from herbaceous and hardwood biomass for nanocellulose isolation.

09 BIOMASS FUELS↗
DOE OSTI2024

Photolytic Degradation of Water‐Soluble Organic Carbon in Snowmelts: Changes in Molecular Characteristics, Brown Carbon Chromophores, and Radiative Effects

Abstract Water‐soluble organic carbon (WSOC) deposited in ambient snowpack play key roles in regional carbon cycle and surface energy budget, but the impacts of photo‐induced processes on its optical and chemical properties are poorly understood yet. In this study, melted samples of the seasonal snow collected from northern Xinjiang, northwestern China, were exposed to ultraviolet (UV) radiation to investigate the photolytic transformations of WSOC. Molecular characteristics and chemical composition of WSOC and its brown carbon (BrC) constituents were investigated using high‐performance liquid chromatography interfaced with a photodiode array detector and a high‐resolution mass spectrometer. Upon illumination, formation of nitrogen‐ and sulfur‐containing species with high molecular weight was observed in snow samples influenced by soil‐ and plant‐derived organics. In contrast, the representative sample collected from remote region showed the lowest molecular diversity and photolytic reactivity among all samples, in which no identified BrC chromophores decomposed upon illumination. Approximately 65% of chromophores in urban samples endured UV irradiation. However, most of BrC composed of phenolic/lignin‐derived compounds and flavonoids disappeared in the illuminated samples containing WSOC from soil‐ and plant‐related sources. Effects of the photochemical degradation of WSOC on the potential modulation of snow albedo were estimated. Apparent half‐lives of WSOC estimated as albedo reduction in 300–400 nm indicated 0.1–0.4 atmospheric equivalent days, which are shorter than typical photolysis half‐lives of ambient biomass smoke aerosol. This study provides new insights into the roles of WSOC in snow photochemistry and snow surface energy balance.

Meteorology & Atmospheric Sciences↗
NASA NTRSDate not supplied

The Influence of Sea ice on Arctic Cloud Properties: What Can We Learn by Applying a In Situ Observational Strategy to Satellite Data?

Sea ice is declining because of anthropogenic climate change. This change alters many aspects of the Arctic climate system, including the way that the surface and atmosphere interact. Atmosphere-surface coupling processes represent an important cloud feedback mechanism that can alter the Arctic surface energy budget. For the Arctic, it has been hypothesized that a reduction in sea ice cover could lead to an increase in clouds. If this process were to occur as originally hypothesized, sea ice loss in all seasons would lead to an increase in clouds. Recent observational studies find a cloud response to sea ice loss in non-summer months and no cloud response in summer months. However, previous studies rely on inter-annual variability and reanalysis to control for the influence of meteorology, reducing the confidence in the resulting conclusions. We adopt a phenomenological, event-based approach that does not need to use meteorological reanalysis. The approach analyzes cloud properties derived from CALIPSO-CloudSat over sea ice and adjacent ice-free footprints by compositing individual satellite ground tracks that cross the Arctic sea ice edge. The underlying assumption, which we verify, is that footprints that are close to each other in space and time experience similar large-scale meteorological conditions. Our results show larger cloud fraction and more cloud liquid water over ice-free than over sea ice footprints and provide additional evidence for a seasonal dependence of cloud-sea ice coupling that is in line with previous work. We find a different result where the maximum cloud property differences between sea ice and ice-free ocean occurs in spring, not fall as earlier studies suggest. We argue that these cloud differences between sea ice and ice-free ocean are primarily caused by the influence of the surface type on the thermodynamic stability of the lower troposphere and not principally from an increase in surface evaporation. In addition, we explore the sensitivity of these results to marginal ice zone width, season, Atlantic vs. Pacific sector, and the cloud property dependence on the distance from the sea ice edge. Overall, our results provide further evidence that cloud-sea ice coupling processes are not offsetting the observed surface energy budget perturbation due to sea ice loss in summer; thus, the cloud response to declining sea ice appears to contribute to amplified Arctic warming.

Patrick C Taylor↗
DOE OSTI2023

Observation of Slow Eigen-Zundel Interconversion in H + (H 2 O) 6 Clusters upon Isomer-Selective Vibrational Excitation and Buffer Gas Cooling in a Cryogenic Ion Trap

The formation of isomers when trapping floppy cluster ions in a temperature-controlled ion trap is a generally observed phenomenon. This involves collisional quenching of the ions initially formed at high temperature by buffer gas cooling until their internal energies fall below the barriers in the potential energy surface that separate them. Here we explore the kinetics at play in the case of the two isomers adopted by the H + (H 2 O) 6 cluster ion that differ in the proton accommodation motif. One of these is most like the Eigen cation with a tricoordinated hydronium motif (denoted E), and the other is most like the Zundel ion with the proton equally shared between two water molecules (denoted Z). After initial cooling to about 20 K in the radiofrequency (Paul) trap, the relative populations of these two spectroscopically distinct isomers are abruptly changed through isomer-selective photoexcitation of bands in the OH stretching region with a pulsed (~6 ns) infrared laser while the ions are in the trap. We then monitor the relaxation of the vibrationally excited clusters and reformation of the two cold isomers by recording infrared photodissociation spectra with a second IR laser as a function of delay time from the initial excitation. The latter spectra are obtained after ejecting the trapped ions into a time-of-flight photofragmentation mass spectrometer, thus enabling long (~0.1 s) delay times. Excitation of the Z isomer is observed to display long-lived vibrationally excited states that are collisionally cooled on a ms time scale, some of which quench into the E isomer. These excited E species then display spontaneous interconversion to the Z form on a ~10 ms time scale. Furthermore, these qualitative observations set the stage for a series of experimental measurements that can provide quantitative benchmarks for theoretical simulations of cluster dynamics and the potential energy surfaces that underlie them.

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

Reactive Molecular Dynamics Simulations and Quantum Chemistry Calculations To Investigate Soot-Relevant Reaction Pathways for Hexylamine Isomers

Sooting tendencies of a series of nitrogen-containing hydrocarbons (NHCs) have been recently characterized experimentally using the yield sooting index (YSI) methodology. This work aims to identify soot-relevant reaction pathways for three selected C6H15N amines, namely, dipropylamine (DPA), diisopropylamine (DIPA), and 3,3-dimethylbutylamine (DMBA) using ReaxFF molecular dynamics (MD) simulations and quantum mechanical (QM) calculations and to interpret the experimentally observed trends. ReaxFF MD simulations are performed to determine the important intermediate species and radicals involved in the fuel decomposition and soot formation processes. QM calculations are employed to extensively search for chemical reactions involving these species and radicals based on the ReaxFF MD results and also to quantitatively characterize the potential energy surfaces. Specifically, ReaxFF simulations are carried out in the NVT ensemble at 1400, 1600, and 1800 K, where soot has been identified to form in the YSI experiment. These simulations account for the interactions among test fuel molecules and pre-existing radicals and intermediate species generated from rich methane combustion, using a recently proposed simulation framework. ReaxFF simulations predict that the reactivity of the amines decrease in the order DIPA > DPA > DMBA, independent of temperature. Both QM calculations and ReaxFF simulations predict that C2H4, C3H6, and C4H8 are the main nonaromatic soot precursors formed during the decomposition of DPA, DIPA, and DMBA, respectively, and the associated reaction pathways are identified for each amine. Both theoretical methods predict that sooting tendency increases in the order DPA, DIPA, and DMBA, consistent with the experimentally measured trend in YSI. This work demonstrates that sooting tendencies and soot-relevant reaction pathways of fuels with unknown chemical kinetics can be identified efficiently through combined ReaxFF and QM simulations. Overall, predictions from ReaxFF simulations and QM calculations are consistent, in terms of fuel reactivity, major intermediates, and major nonaromatic soot precursors.

fuel decomposition↗
DOE OSTI2023

Improving snow albedo modeling in the E3SM land model (version 2.0) and assessing its impacts on snow and surface fluxes over the Tibetan Plateau

With the highest albedo of the land surface, snow plays a vital role in Earth's surface energy budget and water cycle. Snow albedo is primarily controlled by snow grain properties (e.g., size and shape) and light-absorbing particles (LAPs) such as black carbon (BC) and dust. The mixing state of LAPs in snow also has impacts on LAP-induced snow albedo reduction and surface radiative forcing (RF). However, most land surface models assume that snow grain shape is spherical and LAPs are externally mixed with the snow grains. This study improves the snow radiative transfer model in the Energy Exascale Earth System Model version 2.0 (E3SM v2.0) Land Model (ELM v2.0) by considering non-spherical snow grain shapes (i.e., spheroid, hexagonal plate, and Koch snowflake) and internal mixing of dust–snow, and it systematically evaluates the impacts on the surface energy budget and water cycle over the Tibetan Plateau (TP). A series of ELM simulations with different treatments of snow grain shape, mixing state of BC–snow and dust–snow, and sub-grid topographic effects (TOP) on solar radiation are performed. Compared with two remote sensing snow products derived from the Moderate Resolution Imaging Spectroradiometer, the control ELM simulation (ELM_Control) with the default configurations of spherical snow grain shape, internal mixing of BC–snow, external mixing of dust–snow, and without TOP as well as the ELM simulation with new model features (ELM_New) can both capture the overall snow distribution reasonably. Additionally, ELM_New overall shows smaller biases in snow cover fraction than ELM_Control in spring when snowmelt is important for water management. The estimated LAP-induced RF in ELM_New ranges from 0 to 19.3 W m -2 with the area-weighted average value of 1.5 W m -2 that is comparable to the reported values in existing studies. The Koch snowflake shape, among other non-spherical shapes, shows the largest difference from the spherical shape in spring when snow processes related to the surface energy budget and water cycle have high importance. The impacts of the mixing state of LAP in snow are smaller than the shape effects and depend on snow grain shape. Compared to external mixing, internal mixing of LAP–snow can lead to larger snow albedo reduction and snowmelt, which further affect the surface energy budget and water cycle. The individual contributions of non-spherical snow shape, mixing state of LAP–snow, and local topography impacts on the snow and surface fluxes have different signs and magnitudes, and their combined effects may be negative or positive due to complex and nonlinear interactions among the factors. Overall, the changes in net solar radiation in spring due to individual and combined effects range from -28.6 to 16.9 W m -2 and -29.7 to 12.2 W m -2 , respectively. This study advances understanding of the role of snow grain shape and mixing state of LAP–snow in land surface processes and offers guidance for improving snow simulations and RF estimates in Earth system models under climate change.

58 GEOSCIENCES↗
DOE OSTI2020

Artificial linear brush abrasion of coatings for photovoltaic module first-surfaces

Natural soiling and the subsequent necessary cleaning of photovoltaic (PV) modules result in abrasion damage to the cover glass. The durability of the front glass has important economic consequences, including determining the use of anti-reflective and/or anti-soiling coatings as well as the method and frequency of operational maintenance (cleaning). Artificial linear brush abrasion using Nylon 6/12 bristles was therefore examined to explore the durability of representative PV first-surfaces, i.e., the surface of a module incident to direct solar radiation. Specimens examined include silane surface functionalized-, roughened (etched)-, porous silica-coated-, fluoropolymer-coated-, and ceramic (TiO 2 or ZrO 2 /SiO 2 /ZrO 2 /SiO 2 )-coated-glass, which are compared to monolithic-poly(methyl methacrylate) and -glass coupons. Characterization methods used in this study include: optical microscopy, ultraviolet–visible–near-infrared (UV-VIS-NIR) spectroscopy, sessile drop goniometry, white-light interferometry, atomic force microscopy (AFM), and depth-profiling X-ray photoelectron spectroscopy (XPS). The corresponding characteristics examined include: surface morphology, transmittance (i.e., optical performance), surface energy (water contact angle), surface roughness, scratch width and depth, and chemical composition, respectively. The study here was performed to determine coating failure modes; identify characterization methods that can detect nascent failures; compare the durability of popular contemporary coating materials; identify their corresponding damage characteristics; and compare slurry and dry-dust abrasion. This study will also aid in developing an abrasion standard for the PV industry.

14 SOLAR ENERGY↗
DOE OSTI2021

Artificial Linear Brush Abrasion of Coatings for Photovoltaic Module First Surfaces

Natural soiling and the subsequent requisite cleaning of photovoltaic (PV) modules result in abrasion damage to the cover glass. The durability of the front glass has important economic consequences, including determining the use of anti-reflective and/or anti-soiling coatings as well as the method and frequency of operational maintenance (cleaning). Artificial linear brush abrasion using Nylon 6/12 bristles was therefore examined to explore the durability of representative PV first-surfaces, i.e., the surface of a module incident to direct solar radiation. Specimens examined include silane surface functionalized-, roughened (etched)-, porous silica-coated-, fluoropolymer-coated-, and ceramic (TiO2 or ZrO2/SiO2/ZrO2/SiO2)-coated-glass, which are compared to monolithic-poly(methyl methacrylate) and -glass coupons. Characterization methods used in this study include: optical microscopy, ultraviolet-visible-near-infrared (UV-VIS-NIR) spectroscopy, sessile drop goniometry, white-light interferometry, atomic force microscopy (AFM), and depth-profiling X-ray photoelectron spectroscopy (XPS). The corresponding characteristics examined include: surface morphology, transmittance (i.e., optical performance), surface energy (water contact angle), surface roughness, scratch width and depth, and chemical composition, respectively. The study here was performed to determine coating failure modes; identify characterization methods that can detect nascent failures; compare the durability of popular contemporary coating materials; identify their corresponding damage characteristics; and compare slurry and dry-dust abrasion. This study will also aid in developing an abrasion standard for the PV industry.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

Phase-field modeling and experiments of dynamic fracture in single crystal quartz

Predicting the onset and characteristics of brittle fracture is important for a wide range of engineering and geological material applications. In this paper, we study important aspects of brittle fracture in α-quartz by phase-field modeling and experiments using a top-down approach. In the modeling framework, the work term in the Griffith energy balance is replaced with internal energy contributions that represent surface energy, thermal energy, and elastic strain energy stored in defects. This allows parametrization of individual energy contributions in terms of internal state variables and keeps track of energy partitioning after the onset of fracture. The path and history dependence of fracture is included in evolution laws for internal state variables, e.g., entropy evolution, while the energy remains a true potential. In the experimental part, dynamic compression experiments coupled with X-ray phase contrast imaging are performed on cube-like samples with a hole. In the top-down analysis, dynamic compression and three point bending experiments from the literature are simulated with the developed phase-field damage model. In conclusion, the fitted model highlights the strain rate, size, and stress state dependence of damage nucleation and evolution in single crystal α-quartz.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

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

Shadow Molecular Dynamics and Atomic Cluster Expansions for Flexible Charge Models

Here, a shadow molecular dynamics scheme for flexible charge models is presented where the shadow Born–Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, is modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born–Oppenheimer molecular dynamics (BOMD). XL-BOMD provides stable dynamics while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO 2 ) and a molecular system of liquid water. The combined ACE+XL-QEq molecular dynamics simulations are stable over a wide range of temperatures both for the oxide and for the molecular systems and provide a precise sampling of the Born–Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

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

Interplay between Facets and Defects during the Dissociative and Molecular Adsorption of Water on Metal Oxide Surfaces

Both surface terminations and defects play a central role in determining how water interacts with metal oxides, thereby setting important properties of the interface that govern reactivity such as the type and distribution of hydroxyl groups. However, the interconnections between facets and defects remain poorly understood, limiting the usefulness of conventional notions such as that hydroxylation is controlled by metal cation exposure at the surface. Here, using hematite (α-Fe 2 O 3 ) as a model system, we show how oxygen vacancies overwhelm surface cation-dependent hydroxylation behavior. Synchrotron-based ambient-pressure X-ray photoelectron spectroscopy was used to monitor the adsorption of molecular water and its dissociation to form hydroxyl groups in situ on (001), (012), or (104) facet-engineered hematite nanoparticles. Supported by density functional theory calculations of the respective surface energies and oxygen vacancy formation energies, the findings show how oxygen vacancies are more prone to form on higher energy facets and induce surface hydroxylation at extremely low relative humidity values of 5 x10 - 5 %. Further, when these vacancies are eliminated, the extent of surface hydroxylation across the facets is as expected from the areal density of exposed iron cations at the surface. These findings help answer fundamental questions about the nature of reducible metal oxide-water interfaces in natural and technological settings and lay the groundwork for rational design of improved oxide-based catalysts.

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

Control of Photoswitching Kinetics with Strong Light–Matter Coupling in a Cavity

Most photochemistry occurs in the regime of weak light–matter coupling, in which a molecule absorbs a photon and then performs photochemistry from its excited state. In the strong coupling regime, enhanced light–matter interactions between an optical field and multiple molecules lead to collective hybrid light–matter states called polaritons. This strong coupling leads to fundamental changes in the nature of the excited states including multi-molecule delocalized excitations, modified potential energy surfaces, and dramatically altered energy levels relative to non-coupled molecules. Here, the effect of strong light–matter coupling on covalent photochemistry has not been well explored. Photoswitches undergo reversible intramolecular photoreactions that can be readily monitored spectroscopically. In this work, we study the effect of strong light–matter coupling on the kinetics of photoswitching within optical cavities. Reproducing prior experiments, photoswitching of spiropyran/merocyanine photoswitches is decelerated in a cavity. Fulgide photoswitches, however, show the opposite effect, with strong coupling accelerating photoswitching. While modified merocyanine switching can be explained by changes in radiative decay rates or the amount of light in the cavity, modified fulgide switching kinetics suggest direct changes to excited-state reaction kinetics.

Cavities↗
DOE OSTI2020

Precise Experimental Test of the Luttinger Theorem and Particle-Hole Symmetry for a Strongly Correlated Fermionic System

A fundamental concept in physics is the Fermi surface, the constant-energy surface in momentum space encompassing all the occupied quantum states at absolute zero temperature. In 1960, Luttinger postulated that the area enclosed by the Fermi surface should remain unaffected even when electron-electron interaction is turned on, so long as the interaction does not cause a phase transition. Understanding what determines the Fermi surface size is a crucial and yet unsolved problem in strongly interacting systems such as high-$T_c$ superconductors. Here we present a precise test of the Luttinger theorem for a two-dimensional Fermi liquid system where the exotic quasiparticles themselves emerge from the strong interaction, namely, for the Fermi sea of composite fermions (CFs). Via direct, geometric resonance measurements of the CFs’ Fermi wave vector down to very low electron densities, we show that the Luttinger theorem is obeyed over a significant range of interaction strengths, in the sense that the Fermi sea area is determined by the density of the minority carriers in the lowest Landau level. Our data in this work also address the ongoing debates on whether or not CFs obey particle-hole symmetry, and if they are Dirac particles. We find that particle-hole symmetry is obeyed, but the measured Fermi sea area differs quantitatively from that predicted by the Dirac model for CFs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Insights on sources and formation mechanisms of liquid-bearing clouds over MOSAiC examined from a Lagrangian framework

Understanding Arctic stratiform liquid-bearing cloud life cycles and properly representing these life cycles in models is crucial for evaluations of cloud feedbacks as well as the faithfulness of climate projections for this rapidly warming region. Examination of cloud life cycles typically requires analyses of cloud evolution and origins on short time scales, on the order of hours to several days. Measurements from the recent Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition provide a unique view of the current state of the central Arctic over an annual cycle. Here, we use the MOSAiC radiosonde measurements to detect liquid-bearing cloud layers over full atmospheric columns and to examine the cloud-generating air masses’ properties. We perform 5-day (120 h) back-trajectory calculations for every detected cloud and cluster them using a unique set of variables extracted from these trajectories informed by ERA5 reanalysis data. This clustering method enables us to separate between the air mass source regions such as ice-covered Arctic and midlatitude open water. We find that moisture intrusions into the central Arctic typically result in multilayer liquid-bearing cloud structures and that more than half of multilayer profiles include overlying liquid-bearing clouds originating in different types of air masses. Finally, we conclude that Arctic cloud formation via prolonged radiative cooling of elevated stable subsaturated air masses circulating over the Arctic can occur frequently (up to 20% of detected clouds in the sounding data set) and may lead to a significant impact of ensuing clouds on the surface energy budget, including net surface warming in some cases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The LAKE model input dataset for three Arctic lakes

This dataset contains meteorological data collected for three Arctic lakes and compiled to satisfy input requirements of the LAKE 2.0 model. The dataset was generated to act as a benchmarking dataset for future model-data inter-comparisons. The LAKE 2.0 model simulates temperatures within the water later and the sedimentary layer of a lake. The LAKE2.0. is an open-source code and available to download via this weblike http://tesla.parallel.ru/Viktor/LAKE/-/wikis/LAKE-model (last visit July 14, 2021). The meteorological data are required to simulate the surface energy balance at the surface of a lake. This dataset includes a compilation of the meteorological data pulled from multiple data streams, including National Oceanic and Atmospheric Administration (NOAA) climate data, Circumarctic Lakes Observation Network (CALON) data, and the United States Geological Survey (USGS) data. The data were compiled for three Arctic lakes: FoxDen (66.55877, -164.45670), Atqasuk (70.452497, -156.951984), and Toolik (68.63150, -149.60740). Each meteorological data is in comma-delimited format (file extension ‘.dat’) and includes eight columns: Temperature [K], Pressure [Pa], longwave downward radiation [W/m2], shortwave downward radiation [W/m2], “U” wind speed [m/s], ”V” wind speed [m/s], humidity [kg/kg], precipitation [m/s]. In addition to the meteorological data file, we included setup and driver files. The Toolik lake is the deepest out of three lakes and has inflowing and outflowing groundwater data. InflowOutflowREADME.txt has more information about inflow and outflow flies. The other two lakes are much shallower and modeled as a closed system (i.e. no water inflow or outflow).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗