Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ATOMIC MODELS”

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 919 records · Page 51

A model for the spatial and energy distributions of energetic neutral atoms produced within the Saturn/Titan plasma system

A model for the production of energetic neutral atoms (ENAs) in the Saturn/Titan system, that takes into account existing data on the spatial distribution of the background neutral density, the energy and spatial distributions of the energetic ions, and the energy dependence of the charge-exchange cross sections, is presented. Hydrogen production via H(+), H is presented in detail. The extended Saturn/Titan system is 'imaged' in energetic neutral H at 66 R(s) and 25-deg latitude. From this vantage point, the energetic H flux is 6 + or - 1 H/sq cm per sec between 10 and 100 keV. Comparing the predicted H spectrum with that observed by Kirsch et al. (1981), the production of neutral O is estimated, and it is concluded that O(+), H charge exchange in the inner magnetosphere is sufficient to account for all the ENA O production.

Hsieh, K. C.↗

Quantum theory of atoms in molecules/charge-charge flux-dipole flux models for fundamental vibrational intensity changes on H-bond formation of water and hydrogen fluoride

The Quantum Theory of Atoms In Molecules/Charge-Charge Flux-Dipole Flux (QTAIM/CCFDF) model has been used to investigate the electronic structure variations associated with intensity changes on dimerization for the vibrations of the water and hydrogen fluoride dimers as well as in the water-hydrogen fluoride complex. QCISD/cc-pVTZ wave functions applied in the QTAIM/CCFDF model accurately provide the fundamental band intensities of water and its dimer predicting symmetric and antisymmetric stretching intensity increases for the donor unit of 159 and 47 km mol{sup −1} on H-bond formation compared with the experimental values of 141 and 53 km mol{sup −1}. The symmetric stretching of the proton donor water in the dimer has intensity contributions parallel and perpendicular to its C{sub 2v} axis. The largest calculated increase of 107 km mol{sup −1} is perpendicular to this axis and owes to equilibrium atomic charge displacements on vibration. Charge flux decreases occurring parallel and perpendicular to this axis result in 42 and 40 km mol{sup −1} total intensity increases for the symmetric and antisymmetric stretches, respectively. These decreases in charge flux result in intensity enhancements because of the interaction contributions to the intensities between charge flux and the other quantities. Even though dipole flux contributions are much smaller than the charge and charge flux ones in both monomer and dimer water they are important for calculating the total intensity values for their stretching vibrations since the charge-charge flux interaction term cancels the charge and charge flux contributions. The QTAIM/CCFDF hydrogen-bonded stretching intensity strengthening of 321 km mol{sup −1} on HF dimerization and 592 km mol{sup −1} on HF:H{sub 2}O complexation can essentially be explained by charge, charge flux and their interaction cross term. Atomic contributions to the intensities are also calculated. The bridge hydrogen atomic contributions alone explain 145, 237, and 574 km mol{sup −1} of the H-bond stretching intensity enhancements for the water and HF dimers and their heterodimer compared with total increments of 149, 321, and 592 km mol{sup −1}, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Revisiting dynamics near a liquid-liquid phase transition in Si and Ga: The fragile-to-strong transition

Using molecular dynamics simulations we analyze the dynamics of two atomic liquids that display a liquid-liquid phase transition (LLPT): Si described by the Stillinger-Weber potential and Ga as modeled by the modified embedded-atom model. In particular, our objective is to investigate the extent to which the presence of a dip in the self-intermediate scattering function is a manifestation of an excess of vibrational states at low frequencies and may be associated with a fragile-to-strong transition (FTST) across the LLPT, as suggested recently. Our results suggest a somewhat different picture. First, in the case of Ga we observe the appearance of an excess of vibrational states at low frequencies, even in the absence of the appearance of a dip in the self-intermediate scattering function across the LLPT. Second, studying the behavior of the shear viscosities traversing the LLPTs we find that both substances are fragile in character above and below their respective LLPT temperatures. Instead of a FTST in an absolute sense these findings are more in line with a view in which the LLPTs are accompanied by a transition from a more fragile to a less fragile liquid. Furthermore, we do not find this transition to correlate with the presence of a dip in the intermediate scattering function.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic-Structural Synergy for Catalytic CO Oxidation over Palladium-Nickel Nanoalloys

Alloying palladium (Pd) with other transition metals at the nanoscale has become an important pathway for preparation of low-cost, highly-active and stable catalysts. However the lack of understanding of how the alloying phase state, chemical composition and atomic-scale structure of the alloys at the nanoscale influence their catalytic activity impedes the rational design of Pd-nanoalloy catalysts. This work addresses this challenge by a novel approach to investigating the catalytic oxidation of carbon monoxide (CO) over palladium-nickel (PdNi) nanoalloys with well-defined bimetallic composition, which reveals a remarkable a maximal catalytic activity at Pd:Ni ratio of ~50:50. Key to understanding the structural-catalytic synergy is the use of high-energy synchrotron X-ray diffraction coupled to atomic pair distribution function (HE-XRD/PDF) analysis to probe the atomic structure of PdNi nanoalloys under controlled thermochemical treatments and CO reaction conditions. Three-dimensional (3D) models of the atomic structure of the nanoalloy particles were generated by reverse Monte Carlo simulations (RMC) guided by the experimental HE-XRD/PDF data. Structural details of the PdNi nanoalloys were extracted from the respective 3D models and compared with the measured catalytic properties. The comparison revealed a strong correlation between the phase state, chemical composition and atomic-scale structure of PdNi nanoalloys and their catalytic activity for CO oxidation. This correlation is further substantiated by analyzing the first atomic neighbor distances and coordination numbers inside the nanoalloy particles and at their surfaces. These findings have provided new insights into the structural synergy of nanoalloy catalysts by controlling the phase state, composition and atomic structure, complementing findings of traditional density functional theory studies.

Environmental Molecular Sciences Laboratory↗

Atomic oxygen interaction with solar array blankets at protective coating defect sites

Atomic oxygen in the low-Earth-orbital environment oxidizes SiOx protected polyimide Kapton solar array blankets at sites which are not protected such as pin windows or scratches in the protective coatings. The magnitude and shape of the atomic oxygen undercutting which occurs at these sites is dependent upon the exposure environment details such as arrival direction and reaction probability. The geometry of atomic oxygen undercutting at defect sites exposed to atomic oxygen in plasma asher was used to develop a Monte Carlo model to simulate atomic oxygen erosion processes at defect sites in protected Kapton. Comparisons of Monte Carlo predictions and experimental results are presented for plasma asher atomic oxygen exposures for large and small defects as well as for protective coatings on one or both sides of Kapton. The model is used to predict in-space exposure results at defect sites for both directed and sweeping atomic oxygen exposure. A comparison of surface textures predicted by the Monte Carlo model and those experimentally observed from both directed space ram and laboratory plasma asher atomic oxygen exposure indicate substantial agreement.

Banks, Bruce A.↗

Solution-based thermodynamic modeling of the Ni-Al-Mo system using first-principles calculations

A solution-based thermodynamic description of the ternary Ni–Al–Mo system is developed here, incorporating first-principles calculations and reported modeling of the binary Ni–Al, Ni–Mo and Al–Mo systems. To search for the configurations with the lowest energies of the N phase, the Alloy Theoretic Automated Toolkit (ATAT) was employed and combined with VASP. The liquid, bcc and γ-fcc phases are modeled as random atomic solutions, and the γ'-Ni3Al phase is modeled by describing the ordering within the fcc structure using two sublattices, summarized as (Al,Mo,Ni)0.75(Al,Mo,Ni)0.25. Thus, γ-fcc and γ'-Ni3Al are modeled with a single Gibbs free energy function with appropriate treatment of the chemical ordering contribution. In addition, notable improvements are the following: first, the ternary effects of Mo and Al in the B2-NiAl and D0a-Ni3Mo phases, respectively, are considered; second, the N-NiAl8Mo3 phase is described as a solid solution using a three-sublattice model; third, the X-Ni14Al75Mo11 phase is treated as a stoichiometric compound. Model parameters are evaluated using first-principles calculations of zero-Kelvin formation enthalpies and reported experimental data. In comparison with the enthalpies of formation for the compounds ψ-AlMo, θ-Al8Mo3 and B2-NiAl, the first-principles results indicate that the N-NiAl8Mo3 phase, which is stable at high temperatures, decomposes into other phases at low temperature. Resulting phase equilibria are summarized in the form of isothermal sections and liquidus projections. To clearly identify the relationship between the γ-fcc and γ'-Ni3Al phases in the ternary Ni–Al–Mo system, the specific γ-fcc and γ'-Ni3Al phase fields are plotted in x(Al)–x(Mo)–T space for a temperature range 1200–1800 K.

36 MATERIALS SCIENCE↗

Local packing modulates diversity of iron pathways and cooperative behavior in eukaryotic and prokaryotic ferritins

Ferritin-like molecules show a remarkable combination of the evolutionary conserved activity of iron uptake and release that engage different pores in the conserved ferritin shell. It was hypothesized that pore selection and iron traffic depend on dynamic allostery with no conformational changes in the backbone. In this study, we detect the allosteric networks in Pseudomonas aeruginosa bacterioferritin (BfrB), bacterial ferritin (FtnA), and bullfrog M and L ferritins (Ftns) by a network-weaving algorithm (NWA) that passes threads of an allosteric network through highly correlated residues using hierarchical clustering. The residue-residue correlations are calculated in the packing-on elastic network model that introduces atom packing into the common packing-off model. Applying NWA revealed that each of the molecules has an extended allosteric network mostly buried inside the ferritin shell. The structure of the networks is consistent with experimental observations of iron transport: The allosteric networks in BfrB and FtnA connect the ferroxidase center with the 4-fold pores and B-pores, leaving the 3-fold pores unengaged. In contrast, the allosteric network directly links the 3-fold pores with the 4-fold pores in M and L Ftns. The majority of the network residues are either on the inner surface or buried inside the subunit fold or at the subunit interfaces. We hypothesize that the ferritin structures evolved in a way to limit the influence of functionally unrelated events in the cytoplasm on the allosteric network to maintain stability of the translocation mechanisms. We showed that the residue-residue correlations and the resultant long-range cooperativity depend on the ferritin shell packing, which, in turn, depends on protein sequence composition. Switching from the packing-on to the packing-off model reduces correlations by 35%–38% so that no allosteric network can be found. The influence of the side-chain packing on the allosteric networks explains the diversity in mechanisms of iron traffic suggested by experimental approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Escape and ionization of atomic oxygen from Io

Model calculations of the neutral O cloud of Io, based on atom-electron-impact excitation and ionization processes in the plasma torus, are presented. The model is quantitatively anchored to the 6300-A O I emission intensity observations of Brown (1981), and uses a plasma temperature and density structure based on Voyager 1 and 2 and EUV data. Parameters predicted by the model include satellite emission flux = 1.5 x 10 to the 9th/sq cm sec, ion-loading rate = 6.2 x 10 to the 26th ions/sec, O mass-loading rate = 16.6 kg/sec, and O ion-energy input rate = 2.7 x 10 to the 10th W. The spatial morphologies of the parameters are shown. Rough estimates obtained by accounting for a neutral S cloud and plasma-torus charge-exchange reactions include O source flux = 1.2 x 10 to the 10th/sq cm sec, ion-loading rate = 4.0 x 10 to the 27th ions/sec, ion-diffusive-loss time = 200 days, plasma mass-loading rate = 150 kg/sec, satellite mass-loss rate = 270 kg/sec, and maximum ion-energy input = 4 x 10 to the 11th W.

Smyth, W. H.↗

UNCERTAINTIES IN ATOMIC DATA AND THEIR PROPAGATION THROUGH SPECTRAL MODELS. I

We present a method for computing uncertainties in spectral models, i.e., level populations, line emissivities, and emission line ratios, based upon the propagation of uncertainties originating from atomic data. We provide analytic expressions, in the form of linear sets of algebraic equations, for the coupled uncertainties among all levels. These equations can be solved efficiently for any set of physical conditions and uncertainties in the atomic data. We illustrate our method applied to spectral models of O III and Fe II and discuss the impact of the uncertainties on atomic systems under different physical conditions. As to intrinsic uncertainties in theoretical atomic data, we propose that these uncertainties can be estimated from the dispersion in the results from various independent calculations. This technique provides excellent results for the uncertainties in A-values of forbidden transitions in [Fe II].

79 ASTRONOMY AND ASTROPHYSICS↗

IC w22_phadiagractox Highlight: Composition-dependent melting temperature of a compound [Poster]

Figure 1 displays, ambient melting points of the Au-Cu system as a function of Cu atomic percent: present model vs. experimental data from four different sources. Figure 2 displays, ambient melting points of the U-O system as a function of O atomic percent: present model vs. both experiment and quantum molecular dynamics (QMD) simulations on the stoichiometric U-O compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffraction Studies of the Atomic Vibrations of Bulk and Surface Atoms in the Reciprocal and Real Spaces of Nanocrystalline SiC

To describe and evaluate the vibrational properties of nanoparticles it is necessary to distinguish between the surface and the core of the particles. Theoretical calculations show that vibrational density of states of the inner atoms of nanograins is similar to bulk material but shifted to higher energies which can be explained by the fact that the gain core is stressed (hardened) due to the presence of internal pressure. Theoretical calculations also show that there is a difference between vibrational properties of a crystal lattice of the grain interior in isolated particles and in a dense (sintered) nanocrystalline material. This is probably due to a coupling of the modes inside the grains via the grain boundaries in dense nanocrystalline bodies. We examined strains present in the surface shell based on examination of diamond and Sic nanocrystals in reciprocal (Bragg-type scattering) and real (PDF analysis) space analysis of neutron diffraction data. Recently we examined the atomic thermal motions in nanocrystalline Sic based on the assumption of a simple Einstein model for uncorrelated atomic notions. According to this model, the Bragg intensity is attenuated as a function of scattering angle by the Debye-Waller factor. Based on this assumption overall temperature factors were determined from the Wilson plots.

Stelmakh, S.↗

Polarizabilities and van der Waals C{sub 6} coefficients of fullerenes from an atomistic electrodynamics model: Anomalous scaling with number of carbon atoms

The van der Waals C{sub 6} coefficients of fullerenes are shown to exhibit an anomalous dependence on the number of carbon atoms N such that C{sub 6} ∝ N{sup 2.2} as predicted using state-of-the-art quantum mechanical calculations based on fullerenes with small sizes, and N{sup 2.75} as predicted using a classical-metallic spherical-shell approximation of the fullerenes. We use an atomistic electrodynamics model where each carbon atom is described by a polarizable object to extend the quantum mechanical calculations to larger fullerenes. The parameters of this model are optimized to describe accurately the static and complex polarizabilities of the fullerenes by fitting against accurate ab initio calculations. This model shows that C{sub 6} ∝ N{sup 2.8}, which is supportive of the classical-metallic spherical-shell approximation. Additionally, we show that the anomalous dependence of the polarizability on N is attributed to the electric charge term, while the dipole–dipole term scales almost linearly with the number of carbon atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ion-beam-induced atomic mixing

Calculations based on the diffusion model are presented of atomic mixing by ion bombardment. This mixing is assumed to have its basis, as does sputtering, in the collision cascades generated by the primary beam. Sharp interfaces within a target are seen to be smoothed by ion bombardment. Mixing may place fundamental limits on the resolution of ion microprobes.

Haff, P. K.↗

Effects of hydrogen atoms on surface conductivity of diamond film

To investigate the effects of surface chemisorbed hydrogen atoms and hydrogen atoms in the subsurface region of diamond on surface conductivity, models of hydrogen atoms chemisorbed on diamond with (100) orientation and various concentrations of hydrogen atoms in the subsurface layer of the diamond were built. By using the first-principles method based on density functional theory, the equilibrium geometries and densities of states of the models were studied. The results showed that the surface chemisorbed hydrogen alone could not induce high surface conductivity. In addition, isolated hydrogen atoms in the subsurface layer of the diamond prefer to exist at the bond centre site of the C-C bond. However, such a structure would induce deep localized states, which could not improve the surface conductivity. When the hydrogen concentration increases, the C-H-C-H structure and C-3H{sub bc}-C structure in the subsurface region are more stable than other configurations. The former is not beneficial to the increase of the surface conductivity. However, the latter would induce strong surface states near the Fermi level, which would give rise to high surface conductivity. Thus, a high concentration of subsurface hydrogen atoms in diamond would make significant contributions to surface conductivity.

36 MATERIALS SCIENCE↗

A multiscale approach to model hydrogen bonding: The case of polyamide

We present a simple multiscale model for polymer chains in which it is possible to selectively remove degrees of freedom. The model integrates all-atom and coarse-grained potentials in a simple and systematic way and allows a fast sampling of the complex conformational energy surface typical of polymers whilst maintaining a realistic description of selected atomistic interactions. In particular, we show that it is possible to simultaneously reproduce the structure of highly directional non-bonded interactions such as hydrogen bonds and efficiently explore the large number of conformations accessible to the polymer chain. We apply the method to a melt of polyamide removing from the model only the degrees of freedom associated to the aliphatic segments and keeping at atomistic resolution the amide groups involved in the formation of the hydrogen bonds. The results show that the multiscale model produces structural properties that are comparable with the fully atomistic model despite being five times faster to simulate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

LINE-INTERLOCKING EFFECTS ON POLARIZATION IN SPECTRAL LINES BY RAYLEIGH AND RAMAN SCATTERING

The polarized spectrum of the Sun and stars is formed from the scattering of anisotropic radiation on atoms. Interpretation of this spectrum requires the solution of polarized line transfer in multilevel atomic systems. While sophisticated quantum theories of polarized line formation in multilevel atomic systems exist, they are limited by the approximation of complete frequency redistribution in scattering. The partial frequency redistribution (PRD) in line scattering is a necessary component in modeling the polarized spectra of strong lines. The polarized PRD line scattering theories developed so far confine themselves to a two-level or a two-term atom model. In this paper, we present a heuristic approach to the problem of polarized line formation in multilevel atoms taking into account the effects of PRD and a weak magnetic field. Starting from the unpolarized PRD multilevel atom approach of Hubeny et al., we incorporate the polarization state of the radiation field. However, the lower level polarization is neglected. Two iterative methods of solving the polarized PRD line transfer in multilevel atoms are also presented. Taking the example of a five-level Ca II atom model, we present illustrative results for an isothermal one-dimensional model atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Moments of dipole oscillator-strength distribution for the helium sequence

The moments S(mu) for mu at least -6 but no more than 2 and L(mu) for mu = 0, 1, and 2 are calculated for the helium sequence for atomic numbers (Z) up to 30 under a screened hydrogenic model. The model describes the atom by single-particle hydrogenic wave functions and treats the initial and the final state as characterized by two different effective charge parameters Zi and Zf, respectively. The differential oscillator strength of the screened hydrogenic model is asymptotically expanded. Assuming the value of 287.6 for the coefficient of the term epsilon to the -7/2 for helium atoms, the parameter Zf is determined for the helium sequence.

Khan, F.↗

Suppression of fine-structure splitting and oscillator strength of sodium D-line in a Debye plasma

We investigate theoretically the influence of static plasma screening on relativistic spin-orbit interaction-induced fine-structure splitting of the D-line doublet arising from the transitions 3p{sub 1/2}–3s{sub 1/2} and 3p{sub 3/2}–3s{sub 1/2} of the valence electron of a sodium atom embedded in a model plasma environment. The many-electron atomic problem is formulated first as an effective one-electron problem in which the interaction between the optically active valence electron and the atomic ion core is represented by an accurate parametric model potential including core-polarization correction, and then the plasma effect on the atomic system is simulated by the Debye-screening model for the valence-core interaction. It is observed that the magnitude of spin-orbit energy shift reduces for both the upper component 3p{sub 3/2} and the lower component 3p{sub 1/2} with increasing plasma screening strength, thereby reducing the spin-orbit energy separation between these two components as the screening becomes stronger. As a consequence, the magnitude of fine-structure splitting between the D{sub 1} and D{sub 2} line energies of sodium drops significantly with stronger plasma screening. The optical (absorption) oscillator strength for 3s → 3p transition is seen to reduce with stronger screening and this leads to a screening-induced gradual suppression of the 3p → 3s spontaneous decay rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗