Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic gradient structure”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Microsized Pore Structure Determination in EPDM Rubbers Using High-Pressure 129 Xe NMR Techniques

Microsized pore parameters, such as pore size and distance between pores in a series of model EPDM rubbers, were determined in situ under the pressure of 500 psi using 129 Xe nuclear magnetic resonance (NMR) techniques: spin–lattice (T 1 ) and spin–spin (T 2 ) relaxation measurements, pulsed-field gradient (PFG) NMR, and two-dimensional exchange spectroscopy (2D EXSY). The T 1 /T 2 (>>1) ratio for the xenon confined in the pores is larger than that for nonconfined free xenon. This suggests that almost the entire pore surface interacts with xenon atoms like a closed pore. While these pores still connect each other through very narrow diffusion/exchange channels, it is possible to observe the echo decay in PFG-NMR and cross-peaks in 2D EXSY. The results show that both diffusion (D pore ≈ 2.1 × 10 –10 m 2 /s) and exchange (exchange rate, τ exch = a few tens of milliseconds) of xenon between a pore within the material and outer surface are prolonged. The exchange distances (l), which correspond to the xenon gas penetration depth, were estimated to be 70–100 μm based on the measured diffusion coefficients and exchange rate (1/τ exch ). NMR diffraction analysis reveals that pore size (a) and pore distance (b) are on the order of magnitude of micrometers and tens of micrometers, while the diffusion coefficients of xenon gas in the diffusion channels (D eff ) are about 10 –8 m 2 /s. Overall, this study suggests that the pores with a few micrometers connected through very narrow flowing channels with the length of several tens of micrometers are developed 70 to 100 μm below the rubber surface. Furthermore, the overall steady-state diffusion of xenon is slower, approximately 2 orders of magnitudes, than the diffusion in the channel between the pores. Finally, the pore and exchange distances correlated with the composition of rubbers showed that the properties of EPDM rubber as a high-pressure gas barrier could be improved by reducing the size of cracks and the depth of gas penetration by the addition of both carbon black and silica fillers.

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

CASPT2 molecular geometries of Fe( II ) spin-crossover complexes

Using fully internally contracted (FIC)-CASPT2 analytical gradients, geometry optimizations of spin-crossover complexes are reported. This approach is tested on a series of Fe(II) complexes with different sizes, ranging from 13 to 61 atoms. A combination of active space and basis set choices are employed to investigate their role in determining reliable molecular geometries. The reported strategy demonstrates that a wave function-based level of theory can be used to optimize the geometries of metal complexes in reasonable times and enables one to treat the molecular geometry and electronic structure of the complexes using the same level of theory. For a series of smaller Fe(II) SCO complexes, strong field ligands in the LS state result in geometries with the largest differences between DFT and CASPT2; however, good agreement overall is observed between DFT and CASPT2. For the larger complexes, moderate sized basis sets yield geometries that compare well with DFT and available experimental data. Here we recommend using the (10e,12o) active space since convergence to a minimum structure was more efficient than with truncated active spaces despite having similar Fe–ligand bond distances.

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

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

Quantum chemical calculation of the equilibrium structures of small metal atom clusters

The application of ab initio quantum mechanical approaches in the study of metal atom clusters requires simplifying techniques that do not compromise the reliability of the calculations. Various aspects of the implementation of the effective core potential (ECP) technique for the removal of the metal atom core electrons from the calculation were examined. The ECP molecular integral formulae were modified to bring out the shell characteristics as a first step towards fulfilling the increasing need to speed up the computation of the ECP integrals. Work on the relationships among the derivatives of the molecular integrals that extends some of the techniques pioneered by Komornicki for the calculation of the gradients of the electronic energy was completed and a formulation of the ECP approach that quite naturally unifies the various state-of-the-art "shape- and Hamiltonian-consistent" techniques was discovered.

Kahn, L. R.↗
DOE OSTI2020

He ion irradiation response of a gradient T91 steel

Metallic materials with a gradient microstructure usually exhibit excellent mechanical properties. How- ever, the radiation response of gradient structural materials is less well understood. In this work, room- temperature He ion irradiation up to ~4.5 dpa with ~10 at% He injection was performed on a gradient T91 steel processed by surface severe plastic deformation. In comparison to the coarse-grained ferritic T91, the gradient T91 with the nanocrystalline layers shows improved radiation tolerance in terms of less bubble swelling and radiation hardening. Additionally, bubble distribution along grain boundaries depends on misorientation angle suggesting a strong capacity of non-equilibrium grain boundaries in storing He atoms. The present study provides insight into the design of radiation tolerant gradient steels for nuclear industry applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Core–Shell Gold Nanoparticle-Star Copolymer Composites with Gradient Transfer and Transport Properties: Toward Electro-Optical Sensors and Catalysis

A series of hybrid core–shell nanomaterials was prepared by solvent phase transfer methods with gold (Au) nanoparticles (NPs) as a core and star copolymers (h-PEI-b-PCL-Cbz) as a shell—consisting of hyperbranched polyethyleneimine (h-PEI) with different poly(ε-caprolactone) (PCL) chain lengths and a carbazole-contained terminating agent (Cbz). By varying the PCL chain length, these core–shell nanomaterials exhibited gradient energy transfer properties, surface hydrophobicity, electropolymerizability, and ionic probe transport and diffusion properties. The systematic copolymer formation was verified with quadruple-detector size exclusion chromatography (SEC) and 1H nuclear magnetic resonance (NMR). High-resolution atomic force microscopy-infrared spectroscopy (AFM-IR) chemical mapping verified nanoparticle-substrate adhesion and film formation. The design of these hybrid NPs and elucidation of gradient properties will be useful for future applications of the core–shell systems in electro-optical sensing and catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nucleon localization function in rotating nuclei

An electron localization function was originally introduced to visualize in positional space bond structures in molecules. It became a useful tool to describe electron configurations in atoms, molecules, and solids. In nuclear physics, a nucleon localization function (NLF) has been used to characterize cluster structures in light nuclei, formation of fragments in fission, and pasta phases appearing in the inner crust of neutron stars. In the case of a deformed rotating nucleus, several NLFs can be introduced, depending on the definition of the spin-quantization axis, direction of the total angular momentum, and self-consistent symmetries of the system. Contributions to the NLF from the current density, spin-current tensor density, and density gradient terms are negligible in the nuclear interior. The oscillating pattern of the simplified NLF can be explained in terms of a constructive interference between kinetic-energy and particle densities. The characteristic nodal pattern seen in the NLF in the direction of major axis of a rotating nucleus comes from single-particle orbits carrying large aligned angular momentum. The variation of the NLF along the minor axis of the nucleus can be traced back to deformation-aligned orbits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Electro-chemo-mechanically Driven Ni Exsolution from (Pr,Ce,Ni)O 2−δ : Controlled Nucleation Density and Enhanced Electrode Kinetics

In situ exsolution of metal nanoparticles is a promising strategy to prepare electrocatalysts with enhanced activity and resistance to agglomeration for efficient chemical transformations and energy conversion. Achieving a high nucleation density of nanoparticles under mild conditions and understanding how to tailor the process is important for performance of these electrodes in electrochemical cells. In this work, we demonstrate facile exsolution of Ni nanoparticles using fluorite-structured (Pr,Ce)O 2−δ as the support oxide, driven by electrochemical potential and aided by the metastability of Ni in the solid solution (elastic driving force). We prepare single-phase oriented thin films of (Pr,Ce,Ni)O 2−δ (NPCO) on (Zr,Y)O 2−δ (YSZ) substrates by pulsed laser deposition. With the aid of a high-throughput electrochemical cell that provides a lateral gradient in Nernst voltage, we apply in situ near-ambient pressure synchrotron X-ray photoelectron spectroscopy and ex situ atomic force microscopy to investigate the impact of electrochemical potential on Ni nucleation density. We find that metallic Ni can be successfully exsolved at 550 °C upon cathodic biasing in 20 mTorr O 2 , and its nucleation density increases with increasing electrochemical driving force/decreasing oxygen chemical potential. We further evaluate the electrochemical performance under highly reducing (fuel electrode) conditions by electrochemical impedance spectroscopy. With the exsolved Ni nanoparticles, the surface exchange coefficient of the NPCO is found to be ∼4× higher than for PCO without exsolution. This work confirms mixed conducting fluorites as beneficial host lattices for facile transition-metal exsolution and suggests the possibility for constructing an all ceria-based electrochemical cell with PCO serving as both the cathode and the anode.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

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

Kohn-Sham Solver (KSSOLV) v2.0

KSSOLV is a MATLAB toolbox for solving Kohn-Sham density functional theory based electronic structure eigenvalue problems. It uses an object oriented features of MATLAB to represent atom, molecules, wavefunctions and Hamiltonians and their operations. It is designed to make it easier for users to prototype and test new algorithms for solving the Kohn-Sham problem. KSSOLV2.0 contains significant improvement over the original KSSOLV described in a paper published in ACM Transaction on Mathematical Software (attached). In addition to performing ground state calculation for small molecules, it can also perform geometry optimization for both molecules and solids. It uses standard pseudopotentials and implements local density approximation, generalized gradient approximation and hybrid functionals. Future releases will also include time-dependent DFT and post DFT calculations such as the GW quasi-particle energy calculation and Bethe-Salpeter equation solver for optical absorption.

Yang, Chao↗
NASA NTRS1994

A study of X-ray and infrared emissions from dusty nonradiative shock waves

We have constructed models that predict the dynamic evolution and infrared (IR) emission of grains behind nonradiative shock waves. We present a self-consistent treatment of the effect of grain destruction and heating on the ionization structure and X-ray emission of the postshock gas. Incorporating thermal sputtering, collisional heating, and deceleration of grains in the postshock flow, we predict the IR and X-ray fluxes from the dusty plasma as a function of swept-up column density. Heavy elements such as C, O, Mg, S, Si and Fe are initially depleted from the gas phase but are gradually returned as the grains are destroyed. The injected neutral atoms require some time to 'catch up' with the ionization state of the ambient gas. The nonequilibrium ionization state and gradient in elemental abundances in the postshock flow produces characteristic X-ray signatures that can be related to the age of the shock and amount of grain destruction. We study the effects of preshock density and shock velocity on the X-ray and IR emission from the shock. We show that the effects of graindestruction on the X-ray spectra of shock waves are substantial. In particular, temperatures derived from X-ray spectra of middle-aged remnants are likely to be overestimated by approximately 15% if cosmic abundances are assumed. Due to the long timescales for grain destruction in X-ray gases over a wide range of temperatures, we suggest that future X-ray spectra studies of supernova remnants be based on depleted abundances instead of cosmic abundances. Our model predictions agree reasonably well with IRAS and Einstein IPC observations of the Cygnus Loop.

Vancura, Olaf↗
DOE OSTI2026

Uncovering the linear boron environment in Na 3 BP 2 through solid-state 11 B NMR spectroscopy

Boron-based compounds exhibit a wide range of structural diversity, with potential applications spanning organic and inorganic chemistry. Herein, we focus on the characterization of the linear boron-phosphorus unit P═B═P in Na 3 BP 2 using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and density functional theory (DFT) calculations. High-resolution 11 B ssNMR spectra were recorded at two fields, and key parameters such as chemical shift anisotropy (CSA), quadrupolar coupling constants (C Q ), and electric field gradient (EFG) tensors were extracted. The 11 B NMR results revealed a distinct chemical environment for the two-coordinate boron atom, with a CSA span (Ω) of 280 ppm and a C Q of 3.0 MHz. These values were further validated through periodic plane-wave DFT calculations, which showed good agreement with experimental results. The obtained spectral parameters are compared to other linear boron units, such as the BO 2 motif, providing a broader context for understanding boron coordination in inorganic compounds. This work expands the body of NMR knowledge on boron-containing materials, particularly for linear boron motifs. The findings contribute to the growing field of boron chemistry and its potential applications in advanced materials.

Porter, Andrew P. [Ames Laboratory (AMES), Ames, I↗
DOE OSTI2021

Structural transitions of 4:1 methanol–ethanol mixture and silicone oil under high pressure

A 4:1 (volume ratio) methanol–ethanol (ME) mixture and silicone oil are two of the most widely used liquid pressure-transmitting media (PTM) in high-pressure studies. Their hydrostatic limits have been extensively studied using various methods; however, the evolution of the atomic structures associated with their emerging nonhydrostaticity remains unclear. Here, we monitor their structures as functions of pressure up to ∼30 GPa at room temperature using in situ high-pressure synchrotron x-ray diffraction (XRD), optical micro-Raman spectroscopy, and ruby fluorescence spectroscopy in a diamond anvil cell. No crystallization is observed for either PTM. The pressure dependence of the principal diffraction peak position and width indicates the existence of a glass transition in the 4:1 ME mixture at ∼12 GPa and in the silicone oil at ∼3 GPa, beyond which a pressure gradient emerges and grows quickly with pressure. There may be another liquid-to-liquid transition in the 4:1 ME mixture at ∼5 GPa and two more glass-to-glass transitions in the silicone oil at ∼10 GPa and ∼16 GPa. By contrast, Raman signals only show peak weakening and broadening for typical structural disordering, and Raman spectroscopy seems to be less sensitive than XRD in catching these structural transitions related to hydrostaticity variations in both PTM. These results uncover rich pressure-induced transitions in the two PTM and clarify their effects on hydrostaticity with direct structural evidence. The high-pressure XRD and Raman data on the two PTM obtained in this work could also be helpful in distinguishing between signals from samples and those from PTM in future high-pressure experiments.

Chen, Xiehang (ORCID:0000000168541506)↗
NASA NTRS1982

Radiation-induced defects in SiO2 as determined with XPS

An analytical method based on oxide/substrate intensity ratios is used to deduce the product of the atomic number density and electron mean free path as a function of depth, for the cases of device-quality gate oxides grown on Si substrates and irradiated with zero to 20 eV electrons during in situ X-ray photoelectron spectroscopy. These structures had been thinned stepwise to 25-60 A by means of a wet chemical depth-profiling procedure. Si(3+) species are formed in the Si/SiO2 interface region, and their observation during their relaxation and annihilation is correlated with strained Si-O-Si bonds. The observation of bond of cleavage and bond strain gradients in these samples is used to extend silica devitrification models for the explanation of fixed oxide charge generation and interface states.

Grunthaner, F. J.↗
DOE OSTI2020

PyFLOSIC: Python-based Fermi–Löwdin orbital self-interaction correction

We present pyflosic, an open-source, general-purpose python implementation of the Fermi–Löwdin orbital self-interaction correction (FLO-SIC), which is based on the python simulation of chemistry framework (pyscf) electronic structure and quantum chemistry code. Thanks to pyscf, pyflosic can be used with any kind of Gaussian-type basis set, various kinds of radial and angular quadrature grids, and all exchange-correlation functionals within the local density approximation, generalized-gradient approximation (GGA), and meta-GGA provided in the libxc and xcfun libraries. A central aspect of FLO-SIC is the Fermi-orbital descriptors, which are used to estimate the self-interaction correction. Importantly, they can be initialized automatically within pyflosic; they can also be optimized within pyflosic with an interface to the atomic simulation environment, a python library that provides a variety of powerful gradient-based algorithms for geometry optimization. Although pyflosic has already facilitated applications of FLO-SIC to chemical studies, it offers an excellent starting point for further developments in FLO-SIC approaches, thanks to its use of a high-level programming language and pronounced modularity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Atomic oxygen in the Martian thermosphere

Modern models of thermospheric composition and temperature and of excitation and radiative transfer processes are used to simulate the O I 130-nm emission from Mars measured by the Mariner 9 ultraviolet spectrometer. This paper uses the Mars thermospheric general circulation model calculations (MTGCM) of Bougher et al. (1988) and the Monte Carlo partial frequency redistribution multiple scattering code of Meier and Lee (1982). It is found that the decline in atomic oxygen through the daylight hours predicted by the MTGCM cannot be reconciled with the excess afternoon brightness seen in the data. Oxygen concentrations inferred from the data show a positive gradient through the day, in agreement with the original analysis by Strickland et al. (1973). In addition, the data suggest that the oxygen abundance increases toward high southerly latitudes, in contrast with the MTGCM prediction of high values in the Northern Hemisphere. It appears that solar forcing alone cannot account for the observed characteristics of the Martian thermosphere and that wave and tidal effects may profoundly affect the structure, winds, and composition.

Stewart, A. I. F.↗
NASA NTRS1991

Initial drop size and velocity distributions for airblast coaxial atomizers

Phase Doppler measurements were used to determine initial drop size and velocity distributions after a complete disintegration of coaxial liquid jets. The Sauter mean diameter (SMD) distribution was found to be strongly affected by the structure and behavior of the preceding liquid intact jet. The axial measurement stations were determined from the photographs of the coaxial liquid jet at very short distances (1-2 mm) downstream of the observed break-up locations. Minimum droplet mean velocities were found at the center, and maximum velocities were near the spray boundary. Size-velocity correlations show that the velocity of larger drops did not change with drop size. Drop rms velocity distributions have double peaks whose radial positions coincide with the maximum mean velocity gradients.

Eroglu, H.↗
NASA NTRS1987

Synthesis, physical and chemical properties, and potential applications of graphite fluoride fibers

Graphite fluoride fibers can be produced by fluorinating pristine or intercalated graphite fibers. The higher the degree of graphitization of the fibers, the higher the temperature needed to reach the same degree of fluorination. Pitched based fibers were fluorinated to flourine-to-carbon atom rations between 0 and 1. The graphite fluoride fibers with a fluorine-to-carbon atom ration near 1 have extensive visible structural damage. On the other hand, fluorination of fibers pretreated with bromine or fluorine and bromine result in fibers with a fluorine-to-carbon atom ratio nearly equal to 0.5 with no visible structural damage. The electrical resistivity of the fibers is dependent upon the fluorine to carbon atom ratio and ranged from .01 to 10 to the 11th ohm/cm. The thermal conductivity of these fibers ranged from 5 to 73 W/m-k, which is much larger than the thermal conductivity of glass, which is the regular filler in epoxy composites. If graphite fluoride fibers are used as a filler in epoxy or PTFE, the resulting composite may be a high thermal conductivity material with an electrical resistivity in either the insulator or semiconductor range. The electrically insulating product may provide heat transfer with lower temperature gradients than many current electrical insulators. Potential applications are presented.

Hung, Ching-Cheh↗