Search NASA⌕ Search

SEARCH · Search NASA

Results for “height function”

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 271 records · Page 15

Melt Pool characteristics on surface roughness and printability of 316L stainless steel in laser powder bed fusion

Purpose Surface quality and porosity significantly influence the structural and functional properties of the final product. This study aims to establish and explain the underlying relationships among processing parameters, top surface roughness and porosity level in additively manufactured 316L stainless steel. Design/methodology/approach A systematic variation of printing process parameters was conducted to print cubic samples based on laser power, speed and their combinations of energy density. Melt pool morphologies and dimensions, surface roughness quantified by arithmetic mean height (Sa) and porosity levels were characterized via optical confocal microscopy. Findings The study reveals that the laser power required to achieve optimal top surface quality increases with the volumetric energy density (VED) levels. A smooth top surface (Sa < 15 µm) or a rough surface with humps at high VEDs (VED > 133.3 J/mm 3 ) can serve as indicators for fully dense bulk samples, while rough top surfaces resulting from melt pool discontinuity correlate with high porosity levels. Under insufficient VED, melt pool discontinuity dominates the top surface. At high VEDs, surface quality improves with increased power as mitigation of melt pool discontinuity, followed by the deterioration with hump formation. Originality/value This study reveals and summarizes the formation mechanism of dominant features on top surface features and offers a potential method to predict the porosity by observing the top surface features with consideration of processing conditions.

Engineering↗

Benchmarking and parameter sensitivity of physiological and vegetation dynamics using the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) at Barro Colorado Island, Panama

Plant functional traits determine vegetation responses to environmental variation, but variation in trait values is large, even within a single site. Likewise, uncertainty in how these traits map to Earth system feedbacks is large. We use a vegetation demographic model (VDM), the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), to explore parameter sensitivity of model predictions, and comparison to observations, at a tropical forest site: Barro Colorado Island in Panama. We define a single 12-dimensional distribution of plant trait variation, derived primarily from observations in Panama, and define plant functional types (PFTs) as random draws from this distribution. We compare several model ensembles, where individual ensemble members vary only in the plant traits that define PFTs, and separate ensembles differ from each other based on either model structural assumptions or non-trait, ecosystem-level parameters, which include (a) the number of competing PFTs present in any simulation and (b) parameters that govern disturbance and height-based light competition. While single-PFT simulations are roughly consistent with observations of productivity at Barro Colorado Island, increasing the number of competing PFTs strongly shifts model predictions towards higher productivity and biomass forests. Different ecosystem variables show greater sensitivity than others to the number of competing PFTs, with the predictions that are most dominated by large trees, such as biomass, being the most sensitive. Changing disturbance and height-sorting parameters, i.e., the rules of competitive trait filtering, shifts regimes of dominance or coexistence between early- and late-successional PFTs in the model. Increases to the extent or severity of disturbance, or to the degree of determinism in height-based light competition, all act to shift the community towards early-successional PFTs. In turn, these shifts in competitive outcomes alter predictions of ecosystem states and fluxes, with more early-successional-dominated forests having lower biomass. It is thus crucial to differentiate between plant traits, which are under competitive pressure in VDMs, from those model parameters that are not and to better understand the relationships between these two types of model parameters to quantify sources of uncertainty in VDMs.

54 ENVIRONMENTAL SCIENCES↗

Stream discharge and water quality data for East Fork Poplar Creek beginning 2012

This data note describes 15-min discharge and in situ water quality data at two locations along East Fork Poplar Creek in east Tennessee, USA. Data records include temperature, gauge height, water surface elevation above mean sea level, and volumetric discharge. Water quality measurements include temperature, specific conductance, pH, dissolved oxygen (percent saturation and concentration), turbidity, and less extensive fDOM data at one site. The data records begin in 2012 at one site and 2015 at the second site; monitoring at both sites is ongoing (as of 2021). The goal of this data collection is to improve understanding of watershed functions, hydrologic dynamics, and material flux. The data will contribute to site conceptual and numerical models, exposure and risk evaluation, remediation selection and design, and performance monitoring. The data are publicly available and can be accessed via unique url or DOI.

54 ENVIRONMENTAL SCIENCES↗

Data-Driven Whole-Genome Clustering to Detect Geospatial, Temporal, and Functional Trends in SARS-CoV-2 Evolution

Current methods for defining SARS-CoV-2 lineages ignore the vast majority of the SARS-CoV-2 genome. We develop and apply an exhaustive vector comparison method that directly compares all known SARS-CoV-2 genome sequences to produce novel lineage classifications. We utilize data-driven models that (i) accurately capture the complex interactions across the set of all known SARS-CoV-2 genomes, (ii) scale to leadership-class computing systems, and (iii) enable tracking how such strains evolve geospatially over time. We show that during the height of the original Omicron surge, countries across Europe, Asia, and the Americas had a spatially asynchronous distribution of Omicron sub-strains. Moreover, neighboring countries were often dominated by either different clusters of the same variant or different variants altogether throughout the pandemic. Analyses of this kind may suggest a different pattern of epidemiological risk than was understood from conventional data, as well as produce actionable insights and transform our ability to prepare for and respond to current and future biological threats.

Jacobson, Daniel↗

All-atom molecular dynamics simulations of weak polyionic brushes: influence of charge density on the properties of polyelectrolyte chains, brush-supported counterions, and water molecules

All atom molecular dynamics (MD) simulations of planar Na+-counterion-neutralized polyacrylic acid (PAA) brushes are performed for varying degrees of ionization (and thereby varying charge density) and varying grafting density. Variation in the PE charge density (or degree of ionization) and grafting density leads to massive changes of the properties of the PE molecules (quantified by the changes in the height and the mobility of the PE brushes) as well as the local arrangement and distribution of the brush-supported counterions and water molecules within the brushes. The effect on the counterions is manifested by the corresponding variation of the counterion mobility, counterion concentration, extent of counterion binding to the charged site of the PE brushes, water-in-salt-like structure formation, and counterion-water–oxygen radial distribution function within the PE brushes. On the other hand, the effect on water molecules is manifested by the corresponding variation of water–oxygen–water–oxygen RDF, local water density, water–water and water–PE functional group hydrogen bond networks, static dielectric constant of water molecules, orientational tetrahedral order parameter, and water mobility. Enforcing such varying degree of ionization of weak polyelectrolytes is possible by changing the pH of the surrounding medium. Thus, our results provide insights into the changes in microstructure (at the atomistic level) of weak polyionic brushes at varying pH. Finally, we anticipate that this knowledge will prove to be vital for the efficient design of several nano-scale systems employing PE brushes such as nanomechanical gates, current rectifiers, etc.

42 ENGINEERING↗

The Spatial and Temporal Variability of the Clear Convective Boundary Layer at the ARM SGP Supersite

The convective boundary layer (CBL), also known as the mixing layer, constitutes the critical lower segment of the atmosphere that significantly influences daily human activities. The growing demand for precise weather forecasts is driven by the requirements of agriculture, transportation, and routine societal functions. Here, to enhance understanding of the CBL, this study investigates the spatiotemporal variability in the CBL and its controlling factors using four-year Doppler lidar, surface flux, and profiling measurements at five ARM Southern Great Plains sites within a 100 km radius. This investigation utilizes data collected exclusively under clear-sky conditions or scattered low-cloud conditions. Results reveal significant spatial differences in CBL evolutions. Daily mixing layer heights (MLHs) vary up to 1 km (30% of the mean) in late afternoon. There is a clear east–west contrast: western sites (C1, E32, E37) exhibit higher summer MLH (1.9–2.1 km) and vertical velocity variances (1.0–1.2 m 2 s −2 ) than eastern sites (1.6–1.8 km), reversing in winter. Temporally, the MLH peaks at 70% of the sunrise–sunset interval, the lagging heat flux (HF) peaks at 50%; and the seasonal MLH maxima lag the HF by approximately one month, influenced by nighttime PBL (planetary boundary layer) properties. The HF and lower tropospheric stability are the main factors of influence for the CBL, but site-specific dependencies highlight the critical roles of local factors, underscoring the need for including them in CBL modeling.

ARM SGP site↗

Growth of Cd 0.9 Zn 0.1 Te 1-y Se y Single Crystals for Room Temperature Gamma-Ray Detection

Quaternary Cd 0.9 Zn 0.1 Te 1-y Sey (CZTS) single crystals, a novel room-temperature nuclear radiation detector semiconductor material, have been grown using a modified vertical Bridgman method (VBM) and a travelling heater method (THM). The percentage concentration of selenium in the VBM-grown crystal was 3% and that in the THM-grown crystal was 2%. While the THM Frisch collar detector (4.4 × 4.4 × 10.7 mm 3 ) produced a highly resolved pulse height spectra (PHS) with a resolution of ~1.06% for 662-keV gamma rays without any correction, the VBM-grown detector (10 × 10 × 2 mm 3 ) offered a high energy resolution of ~2% after application of a digital biparametric correction. The high-resolution performance of these detectors has been attributed to the addition of Se in the Cd 0.9 Zn 0.1 Te (CZT) matrix. Ab-initio calculations based on density functional theory (DFT) also confirmed that the addition of Se in the CZT matrix helps to reduce the formation of TeCd and the TeZn anti-sites. The VBM-grown crystals were characterized using powder x-ray diffraction (XRD) and energy dispersive x-ray spectroscopy (EDS). While XRD results revealed sharp diffraction peaks confirming the crystalline nature of the grown crystal, the EDS results confirmed the targeted stoichiometry of the elemental composition. The bulk resistivity of the grown crystal was calculated to be ~3 × 10 10 Ω-cm from current-voltage characteristics recorded at room temperature in a planar configuration, ensuring that the grown CZTS crystals have low dark current as required for detector-grade crystals.

42 ENGINEERING↗

Creating microfluidic channels functionalized with micro- and nano-scale features via femtosecond laser surface processing

We demonstrate the ability to create functionalized microfluidic channels using femtosecond laser surface processing (FLSP). FLSP is an emerging advanced manufacturing technology used to modify the surface properties of materials directly and permanently by producing self-organized quasi-periodic micro- and nano-scale surface features along with surface and subsurface chemical and grain structure changes. We demonstrate on Hastelloy X that by controlling the laser fluence and pulse count, the depth of the microchannels and height of the FLSP microstructures within the microchannels can be controlled independently.

Ultrashort pulse laser applications, micro- and na↗

Etched Silicon Planar CRL Optics for the High-Energy X-ray Diffraction Beamlines 11-ID-B and 11-ID-C at the APS

We report the focusing of high-energy X-rays is most commonly achieved by commercially available compound refractive lenses (CRL), which are added together in larger structures. There is a need for high-energy X-ray beamlines to enable focusing in only one dimension down to a few micrometers to study layered materials and thin films. High-energy X-rays above 50 keV require often 100 lenses or more, therefore it is more efficient to etch the complete arrays of lenses into a monolithic substrate like a silicon wafer. Custom lens arrays were fabricated for the high-energy x-ray beamlines 11-ID-B, and 11-ID-C at the Argonne National Laboratory. The fabrication required tuning parameters of the Bosch process to optimize the etching of deep vertical sidewalls at a reasonable speed without overheating. The beamline 11-ID-B is equipped with lens arrays of 85-110 lenses for 59 keV and 190-250 lenses for 87 keV photon energy, focusing the beam vertically below 5.0 µm at the sample positions between 1600 mm to 2100 mm. 11-ID-C has received arrays of 435 lenses that focus 106 keV X-rays to a vertical height of 2.0 µm at the sample position. In both cases, the flux density gain in the focal spot is larger than a factor of 20, relative to using slits to cut the beam size. The lens setup is permanently installed and can be moved in and out of the beam during the experiment without further alignment. The lenses have enabled operando diffraction studies on interfaces in batteries and the development of pair distribution function (PDF) analysis using grazing incidence geometry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Divide and conquer: separating the two probabilities in seismic phase picking

There are two fundamental probabilities in the seismic phase picking process—the probability of the existence of a seismic phase (detection probability) and the probability associated with the phase arrival time estimation (timing probability). The nearly ubiquitous approach in developing deep learning phase picking models is to use a kernel, such as a truncated Gaussian, to mask the labelled phase arrival time and train a segmentation model. Once a model is trained, the times of the peaks in the output are taken as phase arrival times (picks), and the height of the peaks are taken as ‘probability’ of the picks. Here, we show that this ‘probability’ represents neither the detection nor the timing probability because this approach forces the output to follow the shape of the kernel. We introduce an approach using two models to estimate these two distinct probabilities. We use a binary classifier with a calibrated confidence to address the detection probability and a multiclass classifier to obtain a probability mass function to address the timing probability. This new approach can make the deep learning-based phase picking process more interpretable and provide options to logically control seismic monitoring workflows.

58 GEOSCIENCES↗

Gamma-ray spectroscopy using angular distribution of Compton scattering

We present a novel method to perform gamma-ray spectroscopy suitable for extreme pileup situations, wherein traditional pulse-height spectroscopy is infeasible or impractical. One example of such application is spectral characterization of intense laser-based radiation sources that employ laser wakefield acceleration, in which a large number of photons can be produced in a well-collimated beam on a picosecond or shorter timescale. The method relies on the angularly-resolved calorimetric measurement of Compton scattering. The probability for a photon to undergo Compton scattering into a given solid angle is a function of its energy and is described by the Klein–Nishina formula. By first Compton scattering the incident beam and then measuring the total energies deposited in multiple detectors placed at various scattering angles, the incident beam energy spectrum can be reconstructed without resolving the individual events in detectors. Here, we compare the performance of an maximum likelihood expectation maximization algorithm with that of an artificial neural network to measure the energy spectrum of a beam of monoenergetic gamma rays from 137Cs decay. We demonstrate that the method can be used to reconstruct the characteristic inverse Compton scattering spectra.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Morphological traits for allometric scaling of the European Sea Bass Dicentrarchus labrax (Linnaeus, 1758) from Southern Portugal population

Abstract The present study aimed to determine the allometric scaling among a selection of morphological traits in European sea bass ( Dicentrarchus labrax ) to estimate fish body weight. A set of morphological traits (fish body weight, length, height, and width) were directly measured in 146 fish of a recirculating aquaculture system, with body weights ranging from 17.11 to 652.21 g. In addition, a collection of digital imagery of each anesthetized fish from the side and top views were used to estimate other traits (indirect measures). Multiple regression analysis and regression coefficients were calculated using all possible combinations of biometric data (predictors) to estimate fish body weight, applying different numerical fitting models (linear, log‐linear, quadratic, exponential). The results showed that the best combination of traits for estimating fish body weight were fish body width, length and height, collected from direct measure ( R 2 = 0.995), for a log‐linear model fitting, which revealed more accurate determinations than the most commonly used length–weight relationship. Nevertheless, other combinations of morphological traits and fitting models were also found to be suitable in successfully predict fish body weight, with variability ranging between 92.5% and 98.5%. For indirect measures, the best predictor was a combination of traits from top view (width, eye distance and area without fins) fitted with a log‐linear function. These results comprise a relevant baseline in supporting the high potential of noninvasive methods to accurately follow the growth of European sea bass juveniles, recurring to imagery analysis of anesthetized fish. It has major potential applications in feeding consumption trials and fish growth models, as it allows for continuously following up fish growth under different experimental conditions without therein distress derived from manipulation.

Azevedo, Ana↗

Heat transport and temperature boundary-layer profiles in closed turbulent Rayleigh–Bénard convection with slippery conducting surfaces

Here, we report direct numerical simulations (DNS) of the Nusselt number Nu , the vertical profiles of mean temperature Θ(z) and temperature variance Ω(z) across the thermal boundary layer (BL) in closed turbulent Rayleigh–Bénard convection (RBC) with slippery conducting surfaces ( z is the vertical distance from the bottom surface). The DNS study was conducted in three RBC samples: a three-dimensional cuboid with length L=H and width W=H/4 ( H is the sample height), and two-dimensional rectangles with aspect ratios Γ≡L/H=1 and 10. The slip length b for top and bottom plates varied from 0 to ∞ . The Rayleigh numbers Ra were in the range 10 6 $&LessSlantEqual;$Ra$&LessSlantEqual;$10 10 and the Prandtl number Pr was fixed at 4.3 . As b increases, the normalised Nu/Nu 0 ( Nu 0 is the global heat transport for b =0 ) from the three samples for different Ra and Γ can be well described by the same function Nu/Nu 0 = N 0 tanh( b /λ 0 )+1 , with N 0 =0.8±0.03 . Here λ 0 ≡ L /(2 Nu 0 ) is the thermal boundary layer thickness for b =0 . Considering the BL fluctuations for Pr >1, one can derive solutions of temperature profiles Θ(z) and Ω(z) near the thermal BL for b $&GreaterSlantEqual;$ 0 . When b =0 , the solutions are equivalent to those reported by Shishkina et al. ( Phys. Rev. Lett. , vol. 114, 2015, 114302) and Wang et al. ( Phys. Rev. Fluids , vol. 1, 2016, 082301(R)), respectively, for no-slip plates. For b >0, the derived solutions are in excellent agreement with our DNS data for slippery plates.

42 ENGINEERING↗

Assessing the effect of regularization on the molecular properties predicted by SCAN and self-interaction corrected SCAN meta-GGA

Recent regularization of the SCAN meta-GGA functional (rSCAN) has simplified the numerical complexities of the SCAN functional, alleviating SCAN's stringent demand on the numerical integration grids to some extent. The regularization of rSCAN, however, results in the breaking of some constraints such as the uniform electron gas limit, the slowly varying density limit, and coordinate scaling of the iso-orbital indicator. Here, we assess the effects of regularization on the electronic, structural, vibrational, and magnetic properties of molecules by comparing the SCAN and rSCAN predictions. The properties studied include atomic energies, atomization energies, ionization potentials, electron affinities, barrier heights, infrared intensities, dissociation and reaction energies, spin moments of molecular magnets, and isomer ordering of water clusters. Our results show that rSCAN requires less dense numerical grids and gives very similar results to those of SCAN for all properties examined with the exception of atomization energies, which are worsened in rSCAN. We also examine the performance of self-interaction-corrected (SIC) rSCAN with respect to SIC-SCAN using the Perdew–Zunger (PZ) SIC method. The PZSIC method uses orbital densities to compute one-electron self-interaction errors and places an even more stringent demand on numerical grids. Our results show that SIC-rSCAN gives marginally better performance than SIC-SCAN for almost all properties studied in this work with numerical grids that are on average half or less as dense as that needed for SIC-SCAN.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Urban weather modeling using WRF: linking physical assumptions, code implementation, and observational needs

The Weather Research and Forecasting (WRF) model includes urban schemes that simulate the influence of urban surfaces on the atmosphere using parameterizations for flux, and radiative exchanges. Three core schemes – the Bulk urban parameterization, Single-Layer Urban Canopy Model (SLUCM), and Multi-Layer Urban Canopy Model (MLUCM) – represent increasing levels of complexity. Although the parameterizations within these urban schemes are described in the literature, their specific implementation remains poorly documented, thus slowing down model development efforts. This manuscript provides a roadmap to the three urban schemes in WRF version 4.5.2, presenting equations using the same symbols as in the model code, along with references to code lines, and including graphics and explanations that connect the code to its physical foundations. Our thorough review of the urban parameterizations implemented in WRF version 4.5.2 highlighted a handful of parameters that may introduce discontinuities in simulations: (i) in the SLUCM, a 1 mm h −1 rain rate threshold is employed to switch between two minimum moisture availability parameterizations, thus impacting latent heat flux calculations; (ii) in the SLUCM a threshold is used to partition shortwave radiation into direct and diffuse components; (iii) in all three urban schemes, the bulk Richardson number is employed to select the similarity function, which influences the vertical distribution of heat and momentum. We also identified a highly simplified treatment of the radiative balance on roof surfaces. The implications of these simplifications can be assessed through targeted observations across relevant conditions, including varying precipitation rates, cloud cover, and transitions between stability regimes. Furthermore, the widespread application of the Monin-Obukhov similarity theory in these urban schemes warrants model evaluation under highly stable and unstable conditions and in heterogeneous urban settings with variable land cover and building heights on scales finer than model resolution. To address these challenges, we offer guidance on observational strategies, emphasizing the need for multi-parametric measurements to capture potential compensating biases and multi-height measurements that align with the levels where quantities are diagnostic and prognosed in the model (i.e., the lowest atmospheric level of the WRF model). Finally, our inspection of the code revealed implementation bugs that have now been corrected in WRF versions 4.6.0 and 4.6.1. Sensitivity tests over the Atlanta urban area show that these corrections affect surface temperatures, underscoring the importance of performing rigorous documentation and verification of the implementation of parameterizations in model code.

54 ENVIRONMENTAL SCIENCES↗

Green buffers near industrial plants, examples in Jefferson County TX

Green buffers between industrial plants and residential areas have multiple purposes in mitigating air pollution and protecting the environment. This study focuses on providing information on examples of green space development by industries near residential areas in Jefferson County, Texas. One represents a full buffer where residences were totally removed, another is a partial buffer where some residences were removed in a neighborhood and the third represents vegetation of a formerly industrial facility. Additional estimates on benefits of buffer spaces with respect to a pollutant plume are also considered. The study used geospatial and property data from 1966 to 2020. The transition from a residential area to a full or partial buffers highlights possible mitigation of air pollution impacts while serving additional functions, such as boosting stormwater infiltration and improving the aesthetic attractiveness of the area. The transformation from oil tanks to forested space may display the potential for environmental remediation, but not all spaces are suitable for fully forested buffers, as safety and other limitations should be considered. With the analysis of the classic Plume Equation insights can be provided regarding the possible impact of green buffers on ground-level pollutant concentrations. The calculations for atmospheric stability classes A-F and for various stack heights provide an understanding of how different conditions may affect the pollutant concentrations. The presence of a green buffer can play a crucial role in controlling ground-level pollution concentrations and reduce a community’s perception of risk.

54 ENVIRONMENTAL SCIENCES↗

Electron transfer calculations between edge sharing octahedra in hematite, goethite, and annite

A key reaction underlying the charge transport in iron containing oxides, clays, micas is the Fe$^{2+}$-Fe$^{3+}$ exchange reaction between edge-sharing iron octahedra. These reactions facilitate conduction in these minerals by the thermally-activated hopping of small polarons across the lattice. Depending on the mineral and local charge state the small polaron can either encase an electron or hole. The probability for conduction of small polarons depends strongly on the height and adiabicity of the reaction barrier, with larger and more diabatic barriers yielding slow conduction associated with either weak coupling or a large prerequisite rearrangement of the lattice during charge transport. To model these reactions, a first principle electron transfer (ET) method was developed to model the small polaron hopping between the edge-sharing octahedra sites in hematite ($e^{-}$ polaron), goethite ($e^{-}$ polaron), and annite ($h^{+}$ polaron) bulk structures. The ET method is based on electronic structure methods (i.e., plane-wave Density Functional Theory) capable of performing calculations with periodic cells and large size systems efficiently while at the same time being accurate enough to be used in the estimation of the electron-transfer coupling matrix element, $V_{AB}$, and the electron transfer transmission factor, $\kappa_{el}$. Additionally, the calculations confirmed the existence of small polarons in all three minerals, and the reactions were predicted to be strongly adiabatic. It was found that transfer of a hole in the octahedral layer of annite had an adiabatic barrier of $0.311$ eV, and the transfer of an extra electron in hematite and goethite had adiabatic barriers of $0.242$ eV and $0.232$ eV respectively. The electronic coupling parameters,$V_{AB}$, were found to be $0.188$ eV, $0.196$ eV, and $0.102$ eV respectively for hematite, goethite, and annite. While similar bonding topologies pertain, the underlying basis for the differences is the subtle differences in local structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗