SEARCH · Engineering Papers
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.
In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners.
Electrochemical advanced oxidative processes (EAOP) are a promising route to destroy recalcitrant organic contaminants such as per- and polyfluoroalkyl substances (PFAS) in drinking water. Central to EAOP are catalysis-induced reactive free radicals for breaking the carbon fluorine bonds in PFAS. Generating these reactive species electrochemically at electrodes provides an advantage over other oxidation processes that rely on chemicals or other harsh conditions. Herein, we report on the performance of niobium (Nb) doped rutile titanium oxide (TiO 2 ) as a novel EAOP catalytic material, combining theoretical modeling with experimental synthesis and characterization. Calculations based on density functional theory are used to predict the overpotential for oxygen evolution at these candidate electrodes, which must be high in order to oxidize PFAS. The results indicate a non-monotonic trend in which Nb doping below 6.25 at.% is expected to reduce performance relative to TiO 2 , while higher concentrations up to 12.5 at.% lead to increased performance, approaching that of state-of-the-art Magnéli Ti 4 O 7 . TiO 2 samples were synthesized with Nb doping concentration at 10 at.%, heat treated at temperatures from 800 to 1100 °C, and found to exhibit high oxidative stability and high generation of reactive oxygen free radical species. The capability of Nb-doped TiO 2 to destroy two common species of PFAS in challenge water was tested, and moderate reduction by ~ 30% was observed, comparable to that of Ti 4 O 7 using a simple three-electrode configuration. We conclude that Nb-doped TiO 2 is a promising alternative EAOP catalytic material with increased activity towards generating reactive oxygen species and warrants further development for electrochemically destroying PFAS contaminants.
In this work, we present an accurate and efficient finite-difference formulation and parallel implementation of Kohn-Sham Density (Operator) Functional Theory (DFT) for non periodic systems embedded in a bulk environment. Specifically, employing non-local pseudopotentials, local reformulation of electrostatics, and truncation of the spatial Kohn-Sham Hamiltonian, and the Linear Scaling Spectral Quadrature method to solve for the pointwise electronic fields in real-space and the non-local component of the atomic force, we develop a parallel finite difference framework suitable for distributed memory computing architectures to simulate non-periodic systems embedded in a bulk environment. Choosing examples from magnesium-aluminum alloys, we first demonstrate the convergence of energies and forces with respect to spectral quadrature polynomial order, and the width of the spatially truncated Hamiltonian. Next, we demonstrate the parallel scaling of our framework, and show that the computation time and memory scale linearly with respect to the number of atoms. Next, we use the developed framework to simulate isolated point defects and their interactions in magnesium-aluminum alloys. Our findings conclude that the binding energies of divacancies, Al solute-vacancy and two Al solute atoms are anisotropic and are dependent on cell size. Furthermore, the binding is favorable in all three cases.
Here, using real-time quantum dynamics calculations, we perform theoretical investigations of light-induced interactions and electronic excitation transfer in a silver nanoparticle dimer. Real-time time-dependent density functional tight-binding (RT-TDDFTB) calculations provide details of the quantum dynamical processes at an electronic/atomistic level with attosecond resolution. The computational efficiency of RT-TDDFTB allows us to examine electronic dynamics up to picosecond time scales. With time scales varying over six orders of magnitude, we provide insight into interactions between the nanoparticle and laser and between nanoparticles. Our results show that the coupling between nanoparticle monomers is dependent on the separation distance between the nanoparticles in the dimer. As the interparticle distance is varied, the dipole–dipole interactions and electronic excitation transfer mechanisms are markedly different. At large distances (from 50 to 20 Å), the energy transfer from NP1 to NP2 becomes more efficient as the interparticle distance decreases. The total dipole moment of the Ag 14 nanoparticle dimer increases linearly at an interparticle distance of 20 Å and reaches its maximum after 1.2 ps. The electronic excitation transfer is also the most efficient at 20 Å. At short distances, back-transfer effects reduce the ability of the dimer and NP1 to accept energy from the incident electric field. We attribute the distance-dependent features of the nanoparticle dimer to the beating between the laser acting on NP1 and the back transfer from NP2 to NP1.
Here, recent studies show the nonrelativistic antiferromagnetic ordering could generate momentum-dependent spin splitting analogous to the Rashba effect but free from the requirement of relativistic spin-orbit coupling. Whereas the classification of such compounds can be illustrated by different spin-splitting prototypes (SSTs) from symmetry analysis and density-functional-theory calculations, the huge variation in chemical bonding and structures of these diverse compounds possibly clouds the issue of how much of the variation in spin splitting can be traced back to the symmetry-defined characteristics, rather to the underlining chemical and structural diversity. The alternative model Hamiltonian approaches do not confront the issues of chemical and structural complexity but often consider only the magnetic sublattice, dealing with the all-important effects of the nonmagnetic ligands via renormalizing the interactions between the magnetic sites. To this end, we constructed a DFT model Hamiltonian that allows us to study SSTs at constant chemistry while retaining the realistic atomic-scale structure including ligands. This is accomplished by using a single, universal magnetic skeletal lattice (Ni 2+ ions in rocksalt NiO) and designing small displacements of the nonmagnetic (oxygen) sublattice which produce, by design, the different SST magnetic symmetries. We show that (i) even similar crystal structures having very similar band structures can lead to contrasting behavior of spin splitting vs momentum, and (ii) even subtle deformations of the nonmagnetic ligand sublattice could cause a giant spin splitting in AFM-induced SST. This is a paradigm shift relative to the convention of modeling magnets without considering the nonmagnetic ligand that mediates indirect magnetic interaction (e.g., superexchange).
This paper presents a new methodology for adjoint sensitivity analysis, previously developed in general terms by Cacuci, into a form directly applicable 10 meteorological problems. This technique is illustrated by examining the sensitivity of a blocking index in a two-layer isentropic model. The index represents a response function for the sensitivity analysis that unlike previous meteorological applications, is an operator and not a functional, and thus extends the scope of adjoint sensitivity to general operator-type responses depending on time and/or space. The sensitivity of the blocking index to perturbations introduced into the model atmosphere, as well as to model parameters, is discussed. The methodology of generalized adjoint sensitivity analysis described in this paper constitutes a prototype for further applications in the atmospheric and/or oceanic sciences.
In this write-up, we have outlined the general methodology for treating relativity in the FP-LMTO method. The Schrodinger equation is insufficient for describing the physics of heavy elements because it leaves out any relativistic effects. Therefore DFT based on the Schrodinger equation is not sufficient for studying heavy materials. It is, however, possible to incorporate all of the effects of relativity in DFT if starting with a reformulation of the Hohenberg-Kohn theorem that proceeds from the Dirac equation, rather than the Schrodinger equation. The corresponding theorem in relativistic density functional theory (RDFT) shows that the ground-state density of the system is a functional of not only the charge density as in Schrodinger-based DFT, but also of the current density. This leads to a more elaborate procedure in practice when solving the Dirac-Kohn-Sham equations, the relativistic analogs of the Kohn-Sham equations.
A theory for obtaining approximate solutions to nonlinear problems whose exact solutions require the use of large computational procedures is described. The technique represents in some respects a generalization of the method of base and comparison solutions for flows depending on a parameter. For the generalized problem, the input variable is no longer a parameter but a function that is incremented over its entire domain. After performing calculations for a base configuration and a small number of variations of it, solutions for a large class of configurations can be obtained by forming linear combinations of the solution increments. For a restricted class of problems, approximate solutions can be obtained for general variations of a base configuration by using a function-space derivative estimate obtained from a base solution and a single variation.
Dynamics of dust particles in a comet atmosphere was investigated using Monte-Carlo simulation method taking into account radiation and gravitational forces applied to the particles. The model distribution function found out in VEGA experiments was chosen for calculations. Dust particles with masses more then 10(exp -10) g were considered. The delay effect for the big mass particles is shown. In this paper the boundary distribution function for sputtered particles drifting off the surface was determined using Nochilla's model. General dependences for the normal component of the relative mass velocity of sputtered particles are given as the function of particle energy normalized by Tomas-Fermi energy.
The primary goal of this research was to improve monitoring of temperate forest change using remote sensing. In this context, change includes both clearing of forest due to effects such as fire, logging, or land conversion and forest growth and succession. The Landsat 7 ETM+ proved an extremely valuable research tool in this domain. The Landsat 7 program has generated an extremely valuable transformation in the land remote sensing community by making high quality images available for relatively low cost. In addition, the tremendous improvements in the acquisition strategy greatly improved the overall availability of remote sensing images. I believe that from an historical prespective, the Landsat 7 mission will be considered extremely important as the improved image availability will stimulate the use of multitemporal imagery at resolutions useful for local to regional mapping. Also, Landsat 7 has opened the way to global applications of remote sensing at spatial scales where important surface processes and change can be directly monitored. It has been a wonderful experience to have participated on the Landsat 7 Science Team. The research conducted under this project led to contributions in four general domains: I. Improved understanding of the information content of images as a function of spatial resolution; II. Monitoring Forest Change and Succession; III. Development and Integration of Advanced Analysis Methods; and IV. General support of the remote sensing of forests and environmental change. This report is organized according to these topics. This report does not attempt to provide the complete details of the research conducted with support from this grant. That level of detail is provided in the 16 peer reviewed journal articles, 7 book chapters and 5 conference proceedings papers published as part of this grant. This report attempts to explain how the various publications fit together to improve our understanding of how forests are changing and how to monitor forest change with remote sensing. There were no new inventions that resulted from this grant.
We propose a probabilistic model discovery method for identifying ordinary differential equations governing the dynamics of observed multivariate data. Our method is based on the sparse identification of nonlinear dynamics (SINDy) framework, where models are expressed as sparse linear combinations of pre-specified candidate functions. Promoting parsimony through sparsity leads to interpretable models that generalize to unknown data. Instead of targeting point estimates of the SINDy coefficients, we estimate these coefficients via sparse Bayesian inference. The resulting method, uncertainty quantification SINDy (UQ-SINDy), quantifies not only the uncertainty in the values of the SINDy coefficients due to observation errors and limited data, but also the probability of inclusion of each candidate function in the linear combination. UQ-SINDy promotes robustness against observation noise and limited data, interpretability (in terms of model selection and inclusion probabilities) and generalization capacity for out-of-sample forecast. Sparse inference for UQ-SINDy employs Markov chain Monte Carlo, and we explore two sparsifying priors: the spike and slab prior, and the regularized horseshoe prior. UQ-SINDy is shown to discover accurate models in the presence of noise and with orders-of-magnitude less data than current model discovery methods, thus providing a transformative method for real-world applications which have limited data.
Three forms of the high-frequency asymptotic Green's function for Lilley's equation are reviewed and compared to the exact solution over wide range of Strouhal numbers. The asymmetric approximation, which applies to sources away form the jet axis, and the quasi-symmetric approximation, which is arrived at by making a near-axis source assumption, are both obtained for parallel round jets from a formal Fourier-transform solution. The ray-theory solution, which is the only high-frequency approximation that can be applied to more general mean flows, follows from a WKB ansatz and is shown to be closely related to the asymmetric approximation. The comparisons show that the best overall prediction of the exact Green's function is given by the asymmetric approximation which remains accurate down to a Strouhal number of 1/2. The close relationship between the asymmetric and ray-theory approximations suggests that the high-frequency asymptotic Green's function for more general mean flows would be similarly successful.
The general structure of semi-inclusive polarized electron scattering from polarized spin-1/2 targets is developed for use at all energy scales, from modest-energy nuclear physics applications to use in very high energy particle physics. The leptonic and hadronic tensors that enter in the formalism are constructed in a general covariant way in terms of kinematic factors that are frame dependent but model independent and invariant response functions which contain all of the model-dependent dynamics. In the process of developing the general problem the relationships to the conventional responses expressed in terms of the helicity components of the exchanged virtual photon are presented. For semi-inclusive electron scattering with polarized electrons and polarized spin-1/2 targets one finds that 18 invariant response functions are required, each depending on four Lorentz scalar invariants. Additionally it is shown how the semi-inclusive cross sections are related via integrations over the momentum of the selected coincidence particle and sums over open channels.
Reaction networks calculated using density functional theory (DFT) play a crucial role in understanding and predicting heterogeneous catalyst performance. However, results obtained using different density functionals can vary significantly. Accordingly, it is important to assess the accuracy of available functionals in capturing adsorption energetics derived from single-crystal adsorption microcalorimetry experiments. Here, we evaluate how different density functionals describe the molecular and dissociative adsorption of CH 3 I, the adsorption of CH 3 , I, and H, and the dissociative adsorption of CH 4 on Ni(111). First, we identify energetically-preferred adsorbed states for each adsorbate at three surface coverages and calculate the corresponding enthalpies of adsorption at 160 K. Then, we compare the results to experimental measurements published in the literature. For each adsorption process we find at least one functional that is quantitatively accurate within experimental error, but no functional is accurate (within experimental error) for all adsorption processes studied. Even when assuming an additional DFT error of ±20 kJ/mol on top of the experimental error, only the PBE-D3 and RPBE-D3 functionals are accurate for all considered adsorption systems. Overall, these results indicate that quantitative agreement between density functional theory calculations and experimental measurements is both system- and functional-dependent. Importantly, single-crystal adsorption microcalorimetry experiments will continue to play an important role for benchmarking density functionals.
Covalent organic frameworks (COFs) show great potential for many advanced applications on account of their structural uniqueness. To address the synthetic challenges, facile chemical routes to engineer the porosity, crystallinity, and functionality of COFs are highly sought after. Herein, we report a synthetic approach that employs the Cadogan reaction to introduce nitrogen-containing heterocycles as the linkages in the framework. Irreversible indazole and benzimidazolylidene (BIY) linkages are introduced into COFs for the first time via phosphine-induced reductive cyclization of the common imine linkages following either stepwise or one-pot reaction protocols. Further, the successful linkage transformation introduces new functionalities, as demonstrated in the case of BIY-COF, which displays excellent intrinsic proton conductivity without the need of impregnation with external proton transfer reagents. Such a general strategy will open the window to a broader class of functional porous crystalline materials.
As climate change expands the world’s arid and semiarid regions, sustainable systems that integrate food and energy production are becoming increasingly critical. Agrivoltaics—co-locating crops with photovoltaic (PV) panels—offers a dual land-use strategy that mitigates environmental stress by shading crops, conserving soil moisture, and enhancing PV efficiency. While climate-smart crops like the tepary bean ( Phaseolus acutifolius ) are well adapted to heat and drought, little is known about how these crops and their associated soil microbiomes respond to the unique microclimates created by PV shading. This study evaluated tepary bean performance and plant–microbial interactions under PV-shade vs. no shade across three soil amendment treatments at two experimental sites. We assessed plant traits including germination, phenology, biomass, height, as well as yield and bean morphology, alongside shifts in soil microbial composition and functional potential. Plants grown under PV-shade were generally taller, with extended reproductive periods and higher yields: 42% of shaded plants produced beans compared to only 8% under full sun. Shaded plants also produced rounder, higher-quality beans, whereas non-shaded plants yielded flatter, less developed beans. Microbial community composition was more strongly influenced by amendment and site conditions than by shading alone. Key microbial taxa (e.g., Glomeromycetes, Desulfobacterota ) and predicted functions (e.g., denitrification, nitrogen-respiration, sulfate reduction) were associated with differences in plant performance. Finally, combining agrivoltaic systems with targeted soil amendments can enhance crop yield and soil microbial functionality—offering a promising strategy for sustainable agriculture in arid landscapes.
The optimization concepts are described in relation to an overall design process as opposed to a detailed, part-design process where the requirements are firmly stated, the optimization criteria are well established, and a design is known to be feasible. The overall design process starts with the stated requirements. Some of the design criteria are derived directly from the requirements, but others are affected by the design concept. It is these design criteria that define the performance index, or objective function, that is to be minimized within some constraints. In general, there will be multiple objectives, some mutually exclusive, with no clear statement of their relative importance. The optimization loop that is given adjusts the design variables and analyzes the resulting design, in an iterative fashion, until the objective function is minimized within the constraints. This provides a solution, but it is only the beginning. In effect, the problem definition evolves as information is derived from the results. It becomes a learning process as we determine what the physics of the system can deliver in relation to the desirable system characteristics. As with any learning process, an interactive capability is a real attriubute for investigating the many alternatives that will be suggested as learning progresses.
A human mission to Mars and other deep space destinations will require a consistent habitation strategy that applies to disparate environments preferably using multiples of identical hardware. The crew will endure high-g forces during launch, microgravity during transit, and various partial gravities on the destination planet, moon, or asteroid. Habitable volume will likely need to be broken up into smaller modules that do not exceed the capacity of launch vehicles and entry-descent-landing technologies. Though it may be possible to design a different unique, optimized habitat for each environment that the crew may encounter, it is more likely that space agencies and other mission sponsors will only have the resources to develop and produce single multi-use hardware that can be used equally well in all environments, with minimal adaptation by the crew. When sizing habitable volume, mission duration and number of crew become important factors to consider. For example, a pressurized rover cabin may only need to support two crew members for a few days at a time, requiring less volume than a habitat that must support the same number of persons for an entire year. In this study the authors consider constraining factors such as launch vehicle capacity, mission duration, crew size, and Mars lander capacity and study how to divide total pressurized volume, subsystems, equipment, consumables, and supplies into multiple manifests. Preference has been given to solutions that use multiples of identical modules over unique optimized volumes. For example at one extreme, can all the habitat pressure vessels, logistic modules, Mars ascent stage cabins, rover cabins, and airlocks be constructed from multiples of the same small-diameter cylindrical modules? On the other extreme, can logistics, habitation volume, mobility, all be enclosed in a single monolithic volume habitat design? Though there may be advantages and disadvantages to either extreme, the conclusion of this study is that most functions can be distilled into two sizes for hardware: multiples of large diameter modules that generally stay in one place, and multiples of smaller volume cabins that can be moved around or function as cockpits for ascent stages, support vehicles, rovers, airlocks, and logistics delivery.