Search NASA⌕ Search

SEARCH · Search NASA

Results for “COMPUTING”

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 OSTI2023

Computing the Relative Affinity of Chlorophylls a and b to Light-Harvesting Complex II

In plants and algae, the primary antenna protein bound to photosystem II is light-harvesting complex II (LHCII), a pigment–protein complex that binds eight chlorophyll (Chl) a molecules and six Chl b molecules. Chl a and Chl b differ only in that Chl a has a methyl group (–CH 3 ) on one of its pyrrole rings, while Chl b has a formyl group (–CHO) at that position. This blue-shifts the Chl b absorbance relative to Chl a . It is not known how the protein selectively binds the right Chl type at each site. Knowing the selection criteria would allow the design of light-harvesting complexes that bind different Chl types, modifying an organism to utilize the light of different wavelengths. The difference in the binding affinity of Chl a and Chl b in pea and spinach LHCII was calculated using multiconformation continuum electrostatics and free energy perturbation. Both methods have identified some Chl sites where the bound Chl type ( a or b ) has a significantly higher affinity, especially when the protein provides a hydrogen bond for the Chl b formyl group. However, the Chl a sites often have little calculated preference for one Chl type, so they are predicted to bind a mixture of Chl a and b . The electron density of the spinach LHCII was reanalyzed, which, however, confirmed that there is negligible Chl b in the Chl a -binding sites. Finally, it is suggested that the protein chooses the correct Chl type during folding, segregating the preferred Chl to the correct binding site.

chemical calculations↗
DOE OSTI2021

Computational Investigation of the Catalytic Hydrodeoxygenation of Propanoic Acid over a Cu(111) Surface

Cu-based alloy catalysts have recently been investigated experimentally for the hydrodeoxygenation (HDO) of biomass-derived organic acids. Here, the HDO of propanoic acid (PAc) has been studied over Cu(111) by mean-field microkinetic modeling based on parameters obtained from first-principles calculations. Models were developed for the gas- and liquid-phase HDO in condensed water and 1,4-dioxane. In agreement with experimental observations, the gas-phase PAc conversion rate is low at 573 K and increases in liquid water by 1 order of magnitude. In all reaction environments, the decarboxylation mechanism is dominant at low hydrogen partial pressures less than 0.1 bar, and the C–COO bond dissociation is the rate-controlling elementary step. This observation contrasts with the rate-controlling step identified over most group VIII metal surfaces, which is the C–OH bond dissociation in the decarbonylation mechanism. At high hydrogen (H 2 ) partial pressures greater than 10 bar, the HDO of PAc produces propionaldehyde that can readsorb and further react through decarbonylation to produce C 2 alkane products, which is conceptually different from the low H 2 partial pressure scenario. At high H 2 partial pressures, the initial hydrogenation at the carbonyl carbon of PAc becomes the rate-controlling elementary step.

Catalysts↗
DOE OSTI2017

Surface Interactions and Confinement of Methane: A High Pressure Magic Angle Spinning NMR and Computational Chemistry Study

Characterization and modeling of the molecular-level behavior of simple hydrocarbon gases, such as methane, in the presence of both nonporous and nano-porous mineral matrices allows for predictive understanding of important processes in engineered and natural systems. In this study, changes in local electromagnetic environments of the carbon atoms in methane under conditions of high pressure (up to 130 bar) and moderate temperature (up to 346 K) were observed with 13C magic-angle spinning (MAS) NMR spectroscopy while the methane gas was mixed with two model solid substrates: a fumed non-porous, 12 nm particle size silica and a mesoporous silica with 200 nm particle size and 4 nm average pore diameter. Examination of the interactions between methane and the silica systems over temperatures and pressures that include the supercritical regime was allowed by a novel high pressure MAS sample containment system, which provided high resolution spectra collected under in situ conditions. For pure methane, no significant thermal effects were found for the observed 13C chemical shifts at all pressures studied here (28.2 bar, 32.6 bar, 56.4 bar, 65.1 bar, 112.7 bar, and 130.3 bar). However, the 13C chemical shifts of resonances arising from confined methane changed slightly with changes in temperature in mixtures with mesoporous silica. The chemical shift values of 13C nuclides in methane change measurably as a function of pressure both in the pure state and in mixtures with both silica matrices, with a more pronounced shift when meso-porous silica is present. Molecular-level simulations utilizing GCMC, MD and DFT confirm qualitatively that the experimentally measured changes are attributed to interactions of methane with the hydroxylated silica surfaces as well as densification of methane within nanopores and on pore surfaces.

confinement↗
DOE OSTI2021

Computational Investigation of the Bonding in [(η 5 –Cp′) 3 (η 1 –Cp′)M] 1– (M = Pu, U, Ce)

Despite the similar ionic radii for Ce 3+ , U 3+ , and Pu 3+ , [(η 5 –Cp′) 3 (η 1 –Cp′)Ce] 1– (Cp′ = C 5 H 4 SiMe 3 , 1-Ce) displays a significantly longer η 1 –Cp′ distance in the solid state compared to the U 3+ and Pu 3+ analogues. To better understand this observation, a theoretical investigation was undertaken to examine the differences in bonding between the actinides 1-Pu and 1-U and how they compare with 1-Ce. The results show that although the bonding is largely ionic and dominated by ligand (2p)–metal (6d/5d) interactions, the polarization of 5f orbitals plays a significant role for 1-Pu and 1-U compared to the 4f-orbitals of 1-Ce. The lack of Ce(4f) interactions is compensated for by increased participation of the Ce(5d) orbitals relative to the actinide 6d orbitals, particularly for the σ-bound η 1 –Cp′ ligand. Furthermore, the use of multiple theoretical approaches including topological, localization, and energy decomposition approaches shows that 1-Pu and 1-U are very similar in covalent character compared to 1-Ce, though the composition and energy of the different interactions suggest that 1-U presents the strongest overall interactions.

Actinides↗
DOE OSTI2020

Computational and Experimental Mechanistic Insights into the Ethanol-to-Butanol Upgrading Reaction over MgO

The mechanism of ethanol upgrading to higher products is still under debate, especially regarding intermediate species and hydrogenation and dehydrogenation steps. In this work, we conducted a combined theoretical and experimental approach to contribute to this discussion. For such, detailed electronic structure density functional theory calculations (aiming at probing density of states, infrared spectra, geometric parameters, charge densities, and reaction energetics) and diffuse reflectance infrared Fourier transform spectroscopy experiments were carried out revealing the relevance of an appropriate combination of reactive surface sites to support the formation of several intermediates that are formed in the C-C coupling over MgO. The roles of Mg and O sites were also studied under an electronic perspective and different geometrical arrangements. We found that a kink configuration was the most adequate for ethanol to 1-butanol upgrading. Our calculations also gave us arguments to propose distinct reaction routes, whose mutual predominance would depend upon reaction temperature. At temperatures up to 573 K, the so-called β-route, which goes through scission of a Cβ-H bond and formation of an oxametallacycle-like intermediate, would dominate the coupling, whereas at higher temperatures, up to 673 K, a more usual Guerbet mechanism, via an aldol coupling step and then consecutive hydrogenations, would be expected. Here, the theoretical conclusions were followed by a careful experimental strategy using sequential experimental planning techniques in order to estimate accurate parameters with the lowest possible experimental load. Information from these different sources were coupled to develop a mathematical model for the rate of the ethanol upgrading reaction, using a Langmuir-Hinshelwood-Hougen-Watson approach. The developed and statistically validated model adequately described the experimental data at 673 K and 1.1 bar total pressure for ethanol partial pressures in the range from 0 to 20 kPa.

09 BIOMASS FUELS↗
DOE OSTI2019

Computational strategies for design and discovery of nanostructured thermoelectrics

Abstract The contribution of theoretical calculations and predictions in the development of advanced high-performance thermoelectrics has been increasingly significant and has successfully guided experiments to understand as well as achieve record-breaking results. In this review, recent developments in high-performance nanostructured bulk thermoelectric materials are discussed from the viewpoint of theoretical calculations. An effective emerging strategy for boosting thermoelectric performance involves minimizing electron scattering while maximizing heat-carrying phonon scattering on many length scales. We present several important strategies and key examples that highlight the contributions of first-principles-based calculations in revealing the intricate but tractable relationships for this synergistic optimization of thermoelectric performance. The integrated optimization approach results in a fourfold design strategy for improved materials: (1) a significant reduction of the lattice thermal conductivity through multiscale hierarchical architecturing, (2) a large enhancement of the Seebeck coefficient through intramatrix electronic band convergence engineering, (3) control of the carrier mobility through band alignment between the host and second phases, and (4) design of intrinsically low-thermal-conductivity materials by maximizing vibrational anharmonicity and acoustic-mode Gruneisen parameters. These combined effects serve to enhance the power factor while reducing the lattice thermal conductivity. This review provides an improved understanding of how theory is impacting the current state of this field and helps to guide the future search for high-performance thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2016

High-pressure X-ray diffraction, Raman, and computational studies of MgCl 2 up to 1 Mbar: Extensive pressure stability of the β-MgCl 2 layered structure

We studied magnesium chloride (MgCl 2 ) with the rhombohedral layered CdCl 2 -type structure (α-MgCl 2 ), experimentally, using synchrotron angle-dispersive powder x-ray diffraction and Raman spectroscopy using a diamond-anvil cell up to 100 GPa at room temperature and theoretically using first-principles density functional calculations. Our results reveal a pressure-induced second-order structural phase transition to a hexagonal layered CdI 2 -type structure (β-MgCl 2 ) at 0.7 GPa: the stacking sequence of the Cl anions are altered resulting in a reduction of the c-axis length. Theoretical calculations confirm this phase transition sequence and the calculated transition pressure is in excellent agreement with the experiment. Lattice dynamics calculations also reproduce the experimental Raman spectra measured for the ambient and high-pressure phase. According to our experimental results MgCl 2 remains in a 2D layered phase up to 100 GPa and further, the 6-fold coordination of Mg cations is retained. Theoretical calculations of relative enthalpy suggest that this extensive pressure stability is due to a low enthalpy of the layered structure ruling out kinetic barrier effects. Our observation is unusual, as it contradicts with the general structural behavior of highly compressed AB 2 compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2017

In silico environmental chemical science: properties and processes from statistical and computational modelling

Quantitative structure–activity relationships (QSARs) have long been used in the environmental sciences. More recently, molecular modeling and chemoinformatic methods have become widespread. These methods have the potential to expand and accelerate advances in environmental chemistry because they complement observational and experimental data with “in silico” results and analysis. The opportunities and challenges that arise at the intersection between statistical and theoretical in silico methods are most apparent in the context of properties that determine the environmental fate and effects of chemical contaminants (degradation rate constants, partition coefficients, toxicities, etc.). The main example of this is the calibration of QSARs using descriptor variable data calculated from molecular modeling, which can make QSARs more useful for predicting property data that are unavailable, but also can make them more powerful tools for diagnosis of fate determining pathways and mechanisms. Emerging opportunities for “in silico environmental chemical science” are to move beyond the calculation of specific chemical properties using statistical models and toward more fully in silico models, prediction of transformation pathways and products, incorporation of environmental factors into model predictions, integration of databases and predictive models into more comprehensive and efficient tools for exposure assessment, and extending the applicability of all the above from chemicals to biologicals and materials.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Catalytic decomposition of methane into hydrogen and high-value carbons: combined experimental and DFT computational study

Thermocatalytic decomposition (TCD) of methane can produce hydrogen and valuable nanocarbon co-products with low to near-zero CO2 emission. In this study, a series of Pd promoted Ni catalysts, prepared with various Ni/Pd ratios on a CNT support, were evaluated for methane TCD performance. Characterization and calculations using density functional theory (DFT) were carried out to elucidate the activity–structure relationship and growth mechanism of carbon nanomaterials. It was found that the methane conversion and stability of the catalysts were highly dependent on the Ni/Pd ratio and reaction temperature. DFT calculations revealed that the diffusion of carbon in the metal sublayer required for CNT growth was more favorable in the Ni–Pd alloy lattice suggesting that the buildup of carbon in the metal alloy sublayer facilitated the formation of CNTs and CNFs. A cyclic reaction–regeneration process for self-sustained TCD was experimentally demonstrated. In each cycle, a portion of the separated CNT product was used to re-synthesize the Ni–Pd/CNT catalyst for use in the next reaction cycle. After five cycles of operation, the CH4 conversion, morphology and crystallinity of the carbon product remained unchanged.

Wang, I-Wen↗
DOE OSTI2015

Communication: Automatic code generation enables nuclear gradient computations for fully internally contracted multireference theory

Analytical nuclear gradients for fully internally contracted complete active space second-order perturbation theory (CASPT2) are reported. This implementation has been realized by an automated code generator that can handle spin-free formulas for the CASPT2 energy and its derivatives with respect to variations of molecular orbitals and reference coefficients. The underlying complete active space self-consistent field and the so-called Z-vector equations are solved using density fitting. The implementation has been applied to the vertical and adiabatic ionization potentials of the porphin molecule to illustrate its capability.

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

Computationally efficient dielectric calculations of molecular crystals

The microscopic dielectric response is a key quantity for electronic materials such as organic semiconductors. Calculations of this response for molecular crystals are currently either expensive or rely on extreme simplifications such as multipole expansions which lack microscopic detail. We present an alternate approach using a microscopic analogue of the Clausius-Mossotti equation, which constructs the dielectric response of a crystal from an eigenvalue decomposition of the dielectric response of individual molecules. This method can potentially be used to examine the effects of defects, disorder, and surfaces on the dielectric properties of molecular solids.

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

Computational extended magneto-hydrodynamical study of shock structure generated by flows past an obstacle

The magnetized shock problem is studied in the context where supersonic plasma flows past a solid obstacle. This problem exhibits interesting and important phenomena such as a bow shock, magnetotail formation, reconnection, and plasmoid formation. This study is carried out using a discontinuous Galerkin method to solve an extended magneto-hydrodynamic model (XMHD). The main goals of this paper are to present a reasonably complete picture of the properties of this interaction using the MHD model and then to compare the results to the XMHD model. The inflow parameters, such as the magnetosonic Mach number M{sub f} and the ratio of thermal pressure to magnetic pressure β, can significantly affect the physical structures of the flow-obstacle interaction. The Hall effect can also significantly influence the results in the regime in which the ion inertial length is numerically resolved. Most of the results presented are for the two-dimensional case; however, two three-dimensional simulations are presented to make a connection to the important case in which the solar wind interacts with a solid body and to explore the possibility of performing scaled laboratory experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Density-functional theory computer simulations of CZTS{sub 0.25}Se{sub 0.75} alloy phase diagrams

Density-functional theory simulations of CZTS, CZTSe, and CZTS{sub 0.25}Se{sub 0.75} photovoltaic compounds have been performed to investigate the stability of the CZTS{sub 0.25}Se{sub 0.75} alloy vs. decomposition into CZTS, CZTSe, and other secondary compounds. The Gibbs energy for vibrational contributions was estimated by calculating phonon spectra and thermodynamic properties at finite temperatures. It was demonstrated that the CZTS{sub 0.25}Se{sub 0.75} alloy is stabilized not by enthalpy of formation but primarily by the mixing contributions to the Gibbs energy. The Gibbs energy gains/losses for several decomposition reactions were calculated as a function of temperature with/without intermixing and vibration contributions to the Gibbs energy. A set of phase diagrams was built in the multidimensional space of chemical potentials at 300 K and 900 K temperatures to demonstrate alloy stability and boundary compounds at various chemical conditions. It demonstrated for CZTS{sub 0.25}Se{sub 0.75} that the chemical potentials for stability differ between typical processing temperature (∼900 K) and operating temperature (300 K). This implies that as cooling progresses, the flux/concentration of S should be increased in MBE growth to maintain the CZTS{sub 0.25}Se{sub 0.75} in a thermodynamically stable state to minimize phase decomposition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗