Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Directed evolution”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Enhanced Hydrogen Evolution Catalysis of Pentlandite due to the Increases in Coordination Number and Sulfur Vacancy during Cubic‐Hexagonal Phase Transition

Abstract The search for new phases is an important direction in materials science. The phase transition of sulfides results in significant changes in catalytic performance, such as MoS 2 and WS 2 . Cubic pentlandite [cPn, (Fe, Ni) 9 S 8 ] can be a functional material in batteries, solar cells, and catalytic fields. However, no report about the material properties of other phases of pentlandite exists. In this study, the unit‐cell parameters of a new phase of pentlandite, sulfur‐vacancy enriched hexagonal pentlandite (hPn), and the phase boundary between cPn and hPn are determined for the first time. Compared to cPn, the hPn shows a high coordination number, more sulfur vacancies, and high conductivity, which result in significantly higher hydrogen evolution performance of hPn than that of cPn and make the non‐nano rock catalyst hPn superior to other most known nanosulfide catalysts. The increase of sulfur vacancies during phase transition provides a new approach to designing functional materials.

Chemistry↗
DOE OSTI2021

Analysis of the Reactive CO 2 Surface Flux in Electrocatalytic Aqueous Flow Reactors

We study how the mass transfer of a chemical species in shear flow is suppressed by the production of a second species that partially reacts away the first species. The second species is produced at the surface with the first species as a reactant. Our work is directly motivated by electrochemical CO 2 reduction in aqueous flow reactors, where OH– molecules generated by the CO evolution reaction react away CO 2 molecules, ultimately inhibiting the mass transfer of CO 2 to the cathode surface. We derive a simple approximation of the Sherwood number, a dimensionless measure of the mass flux into the surface, as a function of the Péclet number, surface Damköhler number, and bulk Damköhler number, and we compare our approximation to numerical solutions of the governing equations. We find that in the diffusion-limited regime, the Sherwood number is well described by the classical Graetz-Lévêque result reduced by a reaction factor due to the competing bulk reaction; with the stoichiometry relevant to electrochemical CO 2 reduction, this reaction factor is 1/2. While the model problem we solve provides insight into how OH– production affects CO 2 mass transfer, experimental systems often involve more complex chemistry. We thus also show how a common buffered electrolyte, KHCO 3 , affects our results.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Microstructurally resolved electrochemical evolution of mechanical- and irradiation-induced damage in nuclear alloys

Abstract There is a need for high-throughput, scale-relevant, and direct electrochemical analysis to understand the corrosion behavior and sensitivity of nuclear materials that are exposed to extreme (high pressure, temperature, and radiation exposure) environments. We demonstrate the multi-scale, multi-modal application of scanning electrochemical cell microscopy (SECCM) to electrochemically profile corrosion alterations in nuclear alloys in a microstructurally resolved manner. Particularly, we identify that both mechanically deformed and irradiated microstructures show reduced charge-transfer resistance that leads to accelerated oxidation. We highlight that the effects of mechanical deformation and irradiation are synergistic, and may in fact, superimpose each other, with implications including general-, galvanic-, and/or irradiation-activated stress-corrosion cracking. Taken together, we highlight the ability of non-destructive, electrochemical interrogations to ascertain how microstructural alterations result in changes in the corrosion tendency of a nuclear alloy: knowledge which has implications to rank, qualify and examine alloys for use in nuclear construction applications.

Chen, Xin↗
DOE OSTI2022

The Preprocessing of Galaxies in the Early Stages of Cluster Formation in Abell 1882 at z = 0.139

A rare opportunity to distinguish between internal and environmental effects on galaxy evolution is afforded by "SuperGroups," systems that are rich and massive, but include several comparably rich substructures, surrounded by filaments. We present here a multiwavelength photometric and spectroscopic study of the galaxy population in the SuperGroup Abell 1882 (A1882) at z = 0.139, combining new data from the MMT and Hectospec with archival results from the Galaxy And Mass Assembly survey, the Sloan Digital Sky Survey, the Nasa/IPAC Extragalactic Database, the Gemini Multi-Object Spectrograph, and the Galaxy Evolution Explorer. These provide spectroscopic classifications for 526 member galaxies, across wide ranges of local density and velocity dispersion. We identify three prominent filaments along which galaxies seem to be entering the SuperGroup (mostly in E–W directions). A1882 has a well-populated red sequence, containing most galaxies with stellar mass >10 10.5 M Sun , and a pronounced color–density relation even within its substructures. Thus, galaxy evolution responds to the external environment as strongly in these unrelaxed systems as we find in rich and relaxed clusters. From these data, local density remains the primary factor, with a secondary role for distance from the inferred center of the entire structure's potential well. The effects on star formation, as traced by optical and near-UV colors, depend on galaxy mass. We see changes in lower-mass galaxies (M < 10 10.5 M Sun ) at four times the virial radius of major substructures, while the more massive near-UV Green Valley galaxies show low levels of star formation within two virial radii. The suppression of star formation ("quenching") occurs in the infall regions of these structures even before the galaxies enter the denser group environment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Host Star Metallicity of Directly Imaged Wide-orbit Planets: Implications for Planet Formation

Directly imaged planets (DIPs) are self-luminous companions of pre-main-sequence and young main-sequence stars. They reside in wider orbits (∼tens to thousands of astronomical units) and generally are more massive compared to the close-in (≲10 au) planets. Determining the host star properties of these outstretched planetary systems is important to understand and discern various planet formation and evolution scenarios. We present the stellar parameters and metallicity ([Fe/H]) for a subsample of 18 stars known to host planets discovered by the direct imaging technique. We retrieved the high-resolution spectra for these stars from public archives and used the synthetic spectral fitting technique and Bayesian analysis to determine the stellar properties in a uniform and consistent way. For eight sources, the metallicities are reported for the first time, while the results are consistent with the previous estimates for the other sources. Our analysis shows that metallicities of stars hosting DIPs are close to solar with a mean [Fe/H] = −0.04 ± 0.27 dex. The large scatter in metallicity suggests that a metal-rich environment may not be necessary to form massive planets at large orbital distances. We also find that the planet mass–host star metallicity relation for the directly imaged massive planets in wide orbits is very similar to that found for the well-studied population of short-period (≲1 yr) super-Jupiters and brown dwarfs around main-sequence stars.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reconstruction of Thiospinel to Active Sites and Spin Channels for Water Oxidation

Water electrolysis is a promising technique for carbon neutral hydrogen production. A great challenge remains at developing robust and low-cost anode catalysts. Many pre-catalysts are found to undergo surface reconstruction to give high intrinsic activity in the oxygen evolution reaction (OER). The reconstructed oxyhydroxides on the surface are active species and most of them outperform directly synthesized oxyhydroxides. The reason for the high intrinsic activity remains to be explored. Here, a study is reported to showcase the unique reconstruction behaviors of a pre-catalyst, thiospinel CoFe 2 S 4 , and its reconstruction chemistry for a high OER activity. The reconstruction of CoFe 2 S 4 gives a mixture with both Fe–S component and active oxyhydroxide (Co(Fe)O x H y ) because Co is more inclined to reconstruct as oxyhydroxide, while the Fe is more stable in Fe–S component in a major form of Fe 3 S 4 . The interface spin channel is demonstrated in the reconstructed CoFe 2 S 4 , which optimizes the energetics of OER steps on Co(Fe)O x H y species and facilitates the spin sensitive electron transfer to reduce the kinetic barrier of O–O coupling. The advantage is also demonstrated in a membrane electrode assembly (MEA) electrolyzer. Finally, this work introduces the feasibility of engineering the reconstruction chemistry of the precatalyst for high performance and durable MEA electrolyzers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reversible and Selective CO 2 to HCO 2 − Electrocatalysis near the Thermodynamic Potential

Abstract Reversible catalysis is a hallmark of energy‐efficient chemical transformations, but can only be achieved if the changes in free energy of intermediate steps are minimized and the catalytic cycle is devoid of high transition‐state barriers. Using these criteria, we demonstrate reversible CO 2 /HCO 2 − conversion catalyzed by [Pt(depe) 2 ] 2+ (depe=1,2‐ bis (diethylphosphino)ethane). Direct measurement of the free energies associated with each catalytic step correctly predicts a slight bias towards CO 2 reduction. We demonstrate how the experimentally measured free energy of each step directly contributes to the <50 mV overpotential. We also find that for CO 2 reduction, H 2 evolution is negligible and the Faradaic efficiency for HCO 2 − production is nearly quantitative. A free‐energy analysis reveals H 2 evolution is endergonic, providing a thermodynamic basis for highly selective CO 2 reduction.

Cunningham, Drew W.↗
DOE OSTI2023

Engineering a Self-Grown TiO 2 /Ti-MOF Heterojunction with Selectively Anchored High-Density Pt Single-Atomic Cocatalysts for Efficient Visible-Light-Driven Hydrogen Evolution

A photocatalyst TiO 2 /Ti-BPDC-Pt is developed with a self-grown TiO 2 /Ti-metal–organic framework (MOF) heterojunction, i.e., TiO 2 /Ti-BPDC, and selectively anchored high-density Pt single-atomic cocatalysts on Ti-BPDC for photocatalytic hydrogen evolution. This intimate heterojunction, growing from the surface pyrolytic reconstruction of Ti-BPDC, works in a direct Z-scheme, efficiently separating electrons and holes. Pt is selectively anchored on Ti-BPDC by ligands and is found in the form of single atoms with loading up to 1.8 wt %. The selective location of Pt is the electron-enriched domain of the heterojunction, which further enhances the utilization of the separated electrons. This tailored TiO 2 /Ti-BPDC-Pt shows a significantly enhanced activity of 12.4 mmol g -1 h -1 compared to other TiO 2 - or MOF-based catalysts. In conclusion, the structure-activity relationship further proves the balance of two simultaneously exposed domains of heterojunctions is critical to fulfilling this kind of catalyst.

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

Intrinsic and induced quantum quenches for enhancing qubit-based quantum noise spectroscopy

Quantum sensing protocols that exploit the dephasing of a probe qubit are powerful and ubiquitous methods for interrogating an unknown environment. They have a variety of applications, ranging from noise mitigation in quantum processors, to the study of correlated electron states. Here, we discuss a simple strategy for enhancing these methods, based on the fact that they often give rise to an inadvertent quench of the probed system: there is an effective sudden change in the environmental Hamiltonian at the start of the sensing protocol. These quenches are extremely sensitive to the initial environmental state, and lead to observable changes in the sensor qubit evolution. We show how these new features give access to environmental response properties. This enables methods for direct measurement of bath temperature, and for detecting non-thermal equilibrium states. We also discuss how to deliberately control and modulate this quench physics, which enables reconstruction of the bath spectral function. Extensions to non-Gaussian quantum baths are also discussed, as is the application of our ideas to a range of sensing platforms (e.g., nitrogen-vacancy (NV) centers in diamond, semiconductor quantum dots, and superconducting circuits).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Experimental study on the startup of the annular wick type heat pipe using fiber optical temperature measurement technique

This study used optical fiber-distributed temperature sensors to measure the internal and external temperature distributions of a water-cooled heat pipe. The sensor technology used in this study is fiber optical distributed temperature sensing, a distributed sensing technique based on the naturally occurring Rayleigh backscatter in optical fibers. This measurement technique provides maximum spatial resolution for static and semi-static applications. Using this sensor, the temperature distribution of the heat pipe's internal, external, vapor core, and the wick was measured with a spatial resolution of 0.65 mm, a sampling frequency of 40 Hz, and a temperature resolution of 0.1 °C. Through the measured temperature distribution database, the starting phenomenon, the effective length trend, and the limitation onset were observed. From the results, it is found that a high-temperature peak appears at the evaporator if a high initial power (75 W) is imposed on the heat pipe, even after the heat pipe approaches the normal operating status. The peak is not observed in a slower startup (30 W initial power then slowly increased to 75 W). It is also found that the temperature distributions and effective condenser length of the heat pipe highly depend on the cooling conditions. There are variations in the temperature according to the radial direction of the horizontal heat pipe due to gravity. In conclusion, lead and lag of the temperature evolution were observed at the onset of the operating limitations.

42 ENGINEERING↗
DOE OSTI2022

Overview of the COMPASS results *

COMPASS addressed several physical processes that may explain the behaviour of important phenomena. This paper presents results related to the main fields of COMPASS research obtained in the recent two years, including studies of turbulence, L–H transition, plasma material interaction, runaway electron, and disruption physics: Tomographic reconstruction of the edge/SOL turbulence observed by a fast visible camera allowed to visualize turbulent structures without perturbing the plasma. Dependence of the power threshold on the X-point height was studied and related role of radial electric field in the edge/SOL plasma was identified. The effect of high-field-side error fields on the L–H transition was investigated in order to assess the influence of the central solenoid misalignment and the possibility to compensate these error fields by low-field-side coils. Results of fast measurements of electron temperature during ELMs show the ELM peak values at the divertor are around 80% of the initial temperature at the pedestal. Liquid metals were used for the first time as plasma facing material in ELMy H-mode in the tokamak divertor. Good power handling capability was observed for heat fluxes up to 12 MW m –2 and no direct droplet ejection was observed. Partial detachment regime was achieved by impurity seeding in the divertor. The evolution of the heat flux footprint at the outer target was studied. Runaway electrons were studied using new unique systems—impact calorimetry, carbon pellet injection technique, wide variety of magnetic perturbations. Radial feedback control was imposed on the beam. Forces during plasma disruptions were monitored by a number of new diagnostics for vacuum vessel (VV) motion in order to contribute to the scaling laws of sideways disruption forces for ITER. Current flows towards the divertor tiles, incl. possible short-circuiting through PFCs, were investigated during the VDE experiments. The results support ATEC model and improve understanding of disruption loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Nanoscale dynamics during self-organized ion beam patterning of Si. II. Kr + bombardment

Understanding the self-organized ion beam nanopatterning of elemental semiconductors, particularly silicon, is of intrinsic scientific and technological interest. This is the second component of a two-part coherent x-ray scattering and x-ray photon correlation spectroscopy (XPCS) investigation of the kinetics and fluctuation dynamics of nanoscale ripple development on silicon during 1 keV Ar + (part I) and Kr + bombardment at 65° polar angle. Here it is found that the ion-enhanced viscous flow relaxation is essentially equal for Kr + -induced patterning as previously found for Ar + patterning despite the difference in ion masses. However, the magnitude of the surface curvature-dependent roughening rate in the early-stage kinetics is larger for Kr + than for Ar + , consistent with expectations that the heavier ion gives an increased mass redistributive contribution to the initial surface instability. As with the Ar + case, fluctuation dynamics in the late stage show a peak in correlation times at the length scale corresponding to the dominant structural feature on the surface—the ripples. Finally, it is shown that speckle motion during the surface evolution can be analyzed to determine spatial inhomogeneities in erosion rate and ripple velocity. This allows the direction and speed of ripple motion to be measured in a real time experiment. In the present case, ripple motion is found to be into the projected direction of the ion source, in contrast to expectations from an existing sputter erosion driven model with parameters derived from binary collision approximation simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Wettability and capillary effects: Dynamics of pinch-off in unconstricted straight capillary tubes

We study the interfacial evolution of immiscible two-phase flow within a capillary tube in the partial wetting regime using direct numerical simulation. We investigate the flow patterns resulting from the displacement of a more viscous fluid by a less viscous one under a wide range of wettability conditions. We find that beyond a wettability dependent critical capillary number, a uniform displacement by a less viscous fluid can transition into a growing finger that eventually breaks up into discrete blobs by a series of pinch-off events for both wetting and non-wetting contact angles. This study validates previous experimental observations of pinch-off for wetting contact angles and extends those to non-wetting contact angles. We find that the blob length increases with the capillary number. We observe that the time between consecutive pinch-off events decreases with the capillary number and is greater for more wetting conditions in the displaced phase. Here, we further show that the blob separation distance as a function of the difference between the inlet velocity and the contact line speed collapses into two monotonically decreasing curves for wetting and non-wetting contact angles. For the phase separation in the form of pinch-off, this work provides a quantitative study of the emerging length and time scales and their dependence on the wettability conditions, capillary effects, and viscous forces.

42 ENGINEERING↗
DOE OSTI2022

The 3D Kinematics of the Orion Nebula Cluster: NIRSPEC-AO Radial Velocities of the Core Population

The kinematics and dynamics of stellar and substellar populations within young, still-forming clusters provide valuable information for constraining theories of formation mechanisms. Using Keck II NIRSPEC+AO data, we have measured radial velocities for 56 low-mass sources within 4' of the core of the Orion Nebula Cluster (ONC). We also remeasure radial velocities for 172 sources observed with SDSS/APOGEE. These data are combined with proper motions measured using HST ACS/WFPC2/WFC3IR and Keck II NIRC2, creating a sample of 135 sources with all three velocity components. The velocities measured are consistent with a normal distribution in all three components. We measure intrinsic velocity dispersions of (${\sigma }_{{v}_{\alpha }}$, ${\sigma }_{{v}_{\delta }}$, ${\sigma }_{{v}_{r}}$) = (1.64 ± 0.12, 2.03 ± 0.13, ${2.56}_{-0.17}^{+0.16}$) km s –1 . Our computed intrinsic velocity dispersion profiles are consistent with the dynamical equilibrium models from Da Rio et al. (2014) in the tangential direction but not in the line-of-sight direction, possibly indicating that the core of the ONC is not yet virialized, and may require a nonspherical potential to explain the observed velocity dispersion profiles. We also observe a slight elongation along the north–south direction following the filament, which has been well studied in previous literature, and an elongation in the line-of-sight to tangential velocity direction. These 3D kinematics will help in the development of realistic models of the formation and early evolution of massive clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

High-resolution Chemical Abundances of the Nyx Stream

Abstract Nyx is a nearby, prograde, and high-eccentricity stellar stream physically contained in the thick disk, but its origin is unknown. Nyx could be the remnant of a disrupted dwarf galaxy, in which case the associated dark matter substructure could affect terrestrial dark matter direct-detection experiments. Alternatively, Nyx could be a signature of the Milky Way’s disk formation and evolution. To determine the origin of Nyx, we obtained high-resolution spectroscopy of 34 Nyx stars using Keck/HIRES and Magellan/MIKE. A differential chemical abundance analysis shows that most Nyx stars reside in a metal-rich ([Fe/H] > −1) high- α component that is chemically indistinguishable from the thick disk. This rules out the originally suggested scenario that Nyx is the remnant of a single massive dwarf galaxy merger. However, we also identify 5 substantially more metal-poor stars ([Fe/H] ∼ −2.0) whose chemical abundances are similar to those of the metal-weak thick disk. It remains unclear how stars that are chemically identical to the thick disk can be on such prograde, high-eccentricity orbits. We suggest two most likely scenarios: that Nyx is the result of an early minor dwarf galaxy merger, or that it is a record of the early spin-up of the Milky Way disk—although neither perfectly reproduces the chemodynamic observations. The most likely formation scenarios suggest that future spectroscopic surveys should find Nyx-like structures outside of the solar neighborhood.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The influence of kinematics of deformation on polycrystalline halite dynamic recrystallization: Full-field simulation of simple shear versus pure shear

Rock salt, composed mainly of halite, flows viscoplastically over a wide range of geological conditions, strongly impacting the dynamic evolution of sedimentary basins and orogens. Understanding how dislocation creep, which involves dislocation glide, intracrystalline recovery and dynamic recrystallization, influences the microstructure and rheology of halite under various deformation kinematics and temperatures is crucial for enhancing knowledge of salt flow dynamics. Here, this study employs a full-field numerical simulation method to compare the viscoplastic deformation of polycrystalline halite under simple shear and pure shear conditions up to a natural strain of ε = 1.5 at temperatures ranging from 100 °C to 300 °C. The results are presented in terms of crystallographic preferred orientation (CPO), grain shape preferred orientation (SPO), subgrain boundary direction, grain size and strain rate distribution. The results indicate that the crystallographic anisotropy of individual halite crystals is transferred to the polycrystalline scale, resulting in strain localization, particularly in simple shear simulations. The kinematics of deformation affect the evolution and distribution of high strain-rate bands, determining the direction of intragranular substructures and the morphology of strain-induced grain boundaries, with minimal impact on grain size. The intensity of grain boundary migration increases with temperature, significantly influencing grain morphology and size, thereby obscuring strain localization, while having little effect on CPOs. At low strain (ε < 1.0), CPOs relative to both the maximum shortening direction and the grain SPO are similar regardless of the deformation kinematics. At high strain (ε > 1.0), simple shear CPOs exhibit three stronger {100} maxima with a monoclinic symmetry relative to the grain SPO compared to the six {100} maxima with an orthotropic symmetry relative to the grain SPO generated under pure shear. Therefore, microstructures and CPOs can serve as indicators of the strain path in polycrystalline halite under various conditions, aiding in determining the shear sense and elucidating the deformation kinematics of salt structures.

58 GEOSCIENCES↗
DOE OSTI2022

Cloud adjustments from large-scale smoke–circulation interactions strongly modulate the southeastern Atlantic stratocumulus-to-cumulus transition

Smoke from southern Africa blankets the southeastern Atlantic Ocean from June to October, producing strong and competing aerosol radiative effects. Smoke effects on the transition between overcast stratocumulus and scattered cumulus clouds are investigated along a Lagrangian (air-mass-following) trajectory in regional climate and large eddy simulation models. Results are compared with observations from three recent field campaigns that took place in August 2017: ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES), CLouds and Aerosol Radiative Impacts and Forcing: Year 2017 (CLARIFY), and Layered Atlantic Smoke Interactions with Clouds (LASIC). The case study is set up around the joint ORACLES–CLARIFY flight that took place near Ascension Island on 18 August 2017. Smoke sampled upstream on an ORACLES flight on 15 August 2017 likely entrained into the marine boundary layer later sampled during the joint flight. The case is first simulated with the WRF-CAM5 regional climate model in three distinct setups: (1) FireOn, in which smoke emissions and any resulting smoke–cloud–radiation interactions are included; (2) FireOff, in which no smoke emissions are included; (3) RadOff, in which smoke emissions and their microphysical effects are included but aerosol does not interact directly with radiation. Over the course of the Lagrangian trajectory, differences in free tropospheric thermodynamic properties between FireOn and FireOff are nearly identical to those between FireOn and RadOff, showing that aerosol–radiation interactions are primarily responsible for the free tropospheric effects. These effects are non-intuitive: in addition to the expected heating within the core of the smoke plume, there is also a “banding” effect of cooler temperature (~1–2 K) and greatly enhanced moisture (>2 g kg -1 ) at the plume top. This banding effect is caused by a vertical displacement of the former continental boundary layer in the free troposphere in the FireOn simulation resulting from anomalous diabatic heating due to smoke absorption of sunlight that manifests primarily as a few hundred meters per day reduction in large-scale subsidence over the ocean. A large eddy simulation (LES) is then forced with free tropospheric fields taken from the outputs for the WRF-CAM5 FireOn and FireOff runs. Cases are run by selectively perturbing one variable (e.g., aerosol number concentration, temperature, moisture, vertical velocity) at a time to better understand the contributions from different indirect (microphysical), “large-scale” semi-direct (above-cloud thermodynamic and subsidence changes), and “local” semi-direct (below-cloud smoke absorption) effects. Despite a more than 5-fold increase in cloud droplet number concentration when including smoke aerosol concentrations, minimal differences in cloud fraction evolution are simulated by the LES when comparing the base case with a perturbed aerosol case with identical thermodynamic and dynamic forcings. A factor of 2 decrease in background free tropospheric aerosol concentrations from the FireOff simulation shifts the cloud evolution from a classical entrainment-driven “deepening–warming” transition to trade cumulus to a precipitation-driven “drizzle-depletion” transition to open cells, however. The thermodynamic and dynamic changes caused by the WRF-simulated large-scale adjustments to smoke diabatic heating strongly influence cloud evolution in terms of both the rate of deepening (especially for changes in the inversion temperature jump and in subsidence) and in cloud fraction on the final day of the simulation (especially for the moisture “banding” effect). Such large-scale semi-direct effects would not have been possible to simulate using a small-domain LES model alone.

54 ENVIRONMENTAL SCIENCES↗