Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic and molecular clusters”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Moderate-resolution spectral standards from lambda 5600 to lambda 9000

We present a grid of stellar classification spectra of moderate resolution (R approximately 1500) in the range lambda lambda 5600-9000 A, compiled from high signal-to noise spectra of 275 stars, most in the open clusters Praesepe and M67. The grid covers dwarfs from types B8 through M5, giants from G8 through M7, and subgiants from F5 through K0. We catalog atomic and molecular absorption features useful for stellar classification, and demonstrate the use of luminosity-sensitive features to distinguish between late-type dwarf and giant stars. The entire database is made available in digital format on anonymous ftp and through the World Wide Web.

Allen, Lori E.↗
DOE OSTI2021

Computational study of the electrostatic potential and charges of multivalent ionic liquid molecules

The electrostatic interactions in ionic liquids (ILs) can be difficult to quantify, but this information is very valuable for improving intermolecular potentials for performing molecular-level simulations. Although polarizable forcefields can provide improved accuracy, the computational efficiency of fixed charge molecular models is often preferred for larger systems. Here, we propose a multi-scale screening approach for analyzing the charge scaling behavior and other electrostatic properties of ILs, which is more attainable than previous analysis methods that rely on the IL crystalline structures. Based on isolated molecules and small clusters, we study the electrostatic potential properties and charges of multivalent IL molecules using density-functional theory calculations and ab initio molecular dynamics simulations. The charges of the IL molecules are calculated and compared using several different charge estimation techniques. Although there are some differences among the estimated partial charges on the atoms, the total anion and cation charges are very similar. In conclusion, the charge scaling factors are substantial, and they are only slightly affected by the IL cluster size model.

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

Combined Machine Learning and Molecular Dynamics Reveal Two States of Hydration of a Single Functional Group of Cationic Polymeric Brushes

The state of hydration of a macromolecular system regulates a plethora of different properties of such a system. In this article, we develop a novel machine learning (ML) approach, based on the unsupervised clustering algorithm, for probing the hydration behavior of the {N(CH 3 ) 3 } + functional group of the PMETAC [Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride] polyelectrolyte (PE) brush system. The PE brushes and the brush-supported water molecules and counterions (chloride ions) are first described using all-atom molecular dynamics (MD) simulations. The simulation data is subsequently used in our ML framework to identify that (1) the {N(CH 3 ) 3 } + functional groups of the PMETAC brushes have two distinct hydration states with one state (state 1) being characterized by less structured water molecules and the other state (state 2) being characterized by more structured water molecules and (2) an enhancement in the brush grafting density leads to the progressive dissapparenace of state 2. An increase in the grafting density increases the number of chloride counterions in a given volume around the {N(CH 3 ) 3 } + functional group and increases the number of shared water molecules between the {N(CH 3 ) 3 } + and Cl - . The chloride counterions are associated with a hydration layer with much less structured water molecules. Therefore, with an increase in the grafting density, an increase in the percentage of shared water molecules leads to the prevalence of the hydration state [of the {N(CH 3 ) 3 } + moiety] with less structured water molecules. Finally, we explain how the present findings are commensurate with two key previous related results, namely a significantly large chloride ion mobility inside the PMETAC brush layer and the {N(CH 3 ) 3 } + -Cl - average distance remaining independent of the PMETAC brush grafting density. Furthermore, we anticipate that the combined ML-MD-simulation approach proposed in this study can be adapted to probe other soft matter systems to reveal new insights of the underlying mechanisms of emergent phenomenon.

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

Unconventional low temperature decomposition of a saturated hydrocarbon over atomically-dispersed titanium-aluminum-boron catalyst

Sonochemically-synthesized atomically-dispersed titanium-aluminum-boron nanopowder (TiAlB NP) exhibits a remarkable low-temperature catalytic activation of aliphatic C-H bonds at 750 K followed by C-C bond activation thus emerging as a potent low-cost alternative to expensive platinum group metals. Here, the model saturated hydrocarbon, exo-tetrahydrodicyclopentadiene (C10H16), undergoes catalytic decomposition on TiAlB NPs in a chemical microreactor to produce 1,3-cyclopentadiene (c-C5H6), cyclopentene (c-C5H8), and molecular hydrogen (H2) as detected in situ via isomer-selective, single-photon ionization time-of-flight mass spectrometry. Extensive electronic structure theory calculations on model clusters of the catalyst decode a unique synergy among the atomic constituents of the catalyst and chemical bonding in this stepwise, retro Diels Alder reaction: Ti, although insensitive to C-H activation in its metallic state, initiates the catalysis via chemisorption of the hydrocarbon, adjacent B centers readily abstract hydrogen atoms and store them during the catalytic cycle, while Al stabilizes the catalyst structure yet providing space for critical docking sites for the departing hydrocarbons.

Biswas, Souvick↗
DOE OSTI2023

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Ab initio studies of the impact of the Debye-Waller factor on the structural and dynamical properties of amorphous semiconductors: The case of a-Si

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (a-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in a-Si networks. It has been observed that the MSDs associated with tetrahedrally-bonded Si atoms are considerably lower than their dangling-bond counterparts – originating from isolated and vacancy-induced clustered defects – and those on the surface of microvoids, leading to an asymmetric non-gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in a-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of a-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of a-Si at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simplified tuning of long-range corrected density functionals for use in symmetry-adapted perturbation theory

Long considered a failure, second-order symmetry-adapted perturbation theory (SAPT) based on Kohn–Sham orbitals, or SAPT0(KS), can be resurrected for semiquantitative purposes using long-range corrected density functionals whose asymptotic behavior is adjusted separately for each monomer. Here, as in other contexts, correct asymptotic behavior can be enforced via “optimal tuning” based on the ionization energy theorem of density functional theory, but the tuning procedure is tedious, expensive for large systems, and comes with a troubling dependence on system size. Here, we show that essentially identical results are obtained using a fast, convenient, and automated tuning procedure based on the size of the exchange hole. In conjunction with “extended” (X)SAPT methods that improve the description of dispersion, this procedure achieves benchmark-quality interaction energies, along with the usual SAPT energy decomposition, without the hassle of system-specific tuning.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Fast calculation of diffraction patterns from an ensemble of aligned molecules

We report an algorithm to calculate electron diffraction patterns for molecules with anisotropic angular distribution, which is significantly faster than existing methods. The algorithm uses a transform to convert the molecular orientation distribution, which is a function of three Euler angles, to the atom-pair distribution functions which depend on the polar and azimuthal angles. The diffraction signal can then be calculated from the atom-pair distributions. We demonstrate the computation method numerically by calculating electron diffraction patterns for a symmetric top molecule (trifluoroiodomethane) and an asymmetric top molecule (formaldehyde) and show that it reduces the calculation time by approximately two orders of magnitude compared to the standard brute-force method. Here, the method can also be applied to the calculation of x-ray diffraction patterns.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Materials Data on U(NO7)2 by Materials Project

UNO7NO3(O2)2 crystallizes in the orthorhombic Cmc2_1 space group. The structure is zero-dimensional and consists of eight molecular oxygen molecules, four nitric acid molecules, and four UNO7 clusters. In each UNO7 cluster, U is bonded in a distorted rectangular see-saw-like geometry to six O atoms. There are a spread of U–O bond distances ranging from 1.83–2.45 Å. N is bonded in a trigonal planar geometry to three O atoms. There is one shorter (1.22 Å) and two longer (1.29 Å) N–O bond length. There are five inequivalent O sites. In the first O site, O is bonded in a single-bond geometry to one U atom. In the second O site, O is bonded in a single-bond geometry to one U atom. In the third O site, O is bonded in a distorted single-bond geometry to one U and one N atom. In the fourth O site, O is bonded in a single-bond geometry to one N atom. In the fifth O site, O is bonded in a single-bond geometry to one U atom.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Time correlation functions in classical dense fluids

The velocity autocorrelation functions and memory functions of dense classical fluids may be directly obtained from the static radial distribution function g(r) in an approximate way. Following the Mori projection operator formalism, the memory functions may be related to the fluctuating force correlation. At low densities, these functions may be evaluated by following the trajectories of particle pairs in the interatomic potential. At higher densities, the force correlation functions can be evaluated approximately from particle pair trajectories via the potential of the mean force. The results for argon fluids at various densities and temperatures agree satisfactorily with the molecular dynamics and the Enskog values. The decrease of the diffusion coefficient with density is partly due to the nature of g(r) which reflects the stronger clustering of atoms at higher densities.

Tsang, T.↗
DOE OSTI2023

The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures

The primary focus of GAMESS over the last 5 years has been the development of new high-performance codes that are able to take effective and efficient advantage of the most advanced computer architectures, both CPU and accelerators. These efforts include employing density fitting and fragmentation methods to reduce the high scaling of well-correlated (e.g., coupled-cluster) methods as well as developing novel codes that can take optimal advantage of graphical processing units and other modern accelerators. Because accurate wave functions can be very complex, an important new functionality in GAMESS is the quasi-atomic orbital analysis, an unbiased approach to the understanding of covalent bonds embedded in the wave function. Finally, best practices for the maintenance and distribution of GAMESS are also discussed.

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

A new self-adaptive reconstruction method to identify defects through Wigner–Seitz approach

A new self-adaptive reconstruction method based on local atomic structure at any given molecular dynamics (MD) step has been developed in this article. The method can be used in Wigner–Seitz defect analysis approach to correctly and efficiently explore the information of both point defects and complex defect clusters (e.g. dislocation loops and voids) formed after a displacement cascade where the cascade interacts with grain boundaries and/or dislocations. The algorithm and validation are provided in detail. Results for identification of radiation defects during and after cascades interacting with a dislocation network show that the new method can well recognize all simple and complex defects and defect clusters. Thus, this new method provides a totally new way to explore the density and size of radiation defects at atomic scale after complex MD evolution processes, providing correct information to understand and predict radiation damage in materials through atomic simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Materials Data on MgFe4(S3O23)2 by Materials Project

MgO10(Fe2S3O16)2(O2)2 crystallizes in the triclinic P-1 space group. The structure is one-dimensional and consists of two molecular oxygen molecules; one MgO10 cluster; and two Fe2S3O16 ribbons oriented in the (-1, 1, 0) direction. In the MgO10 cluster, Mg is bonded in an octahedral geometry to six O atoms. There are a spread of Mg–O bond distances ranging from 1.97–2.48 Å. There are five inequivalent O sites. In the first O site, O is bonded in a distorted bent 150 degrees geometry to one Mg and one O atom. The O–O bond length is 1.23 Å. In the second O site, O is bonded in a distorted bent 150 degrees geometry to one Mg and one O atom. The O–O bond length is 1.23 Å. In the third O site, O is bonded in a single-bond geometry to one Mg atom. In the fourth O site, O is bonded in a single-bond geometry to one O atom. In the fifth O site, O is bonded in a single-bond geometry to one O atom. In each Fe2S3O16 ribbon, there are two inequivalent Fe sites. In the first Fe site, Fe is bonded to five O atoms to form FeO5 trigonal bipyramids that share corners with three SO4 tetrahedra. There are a spread of Fe–O bond distances ranging from 1.76–2.01 Å. In the second Fe site, Fe is bonded in a 6-coordinate geometry to six O atoms. There are a spread of Fe–O bond distances ranging from 1.93–2.24 Å. There are three inequivalent S sites. In the first S site, S is bonded to four O atoms to form SO4 tetrahedra that share a cornercorner with one FeO5 trigonal bipyramid. There are a spread of S–O bond distances ranging from 1.46–1.53 Å. In the second S site, S is bonded to four O atoms to form SO4 tetrahedra that share a cornercorner with one FeO5 trigonal bipyramid. There are a spread of S–O bond distances ranging from 1.45–1.53 Å. In the third S site, S is bonded to four O atoms to form SO4 tetrahedra that share a cornercorner with one FeO5 trigonal bipyramid. There are a spread of S–O bond distances ranging from 1.46–1.54 Å. There are sixteen inequivalent O sites. In the first O site, O is bonded in a bent 150 degrees geometry to one Fe and one S atom. In the second O site, O is bonded in a single-bond geometry to one S atom. In the third O site, O is bonded in a distorted bent 120 degrees geometry to one Fe and one S atom. In the fourth O site, O is bonded in a single-bond geometry to one S atom. In the fifth O site, O is bonded in a distorted bent 150 degrees geometry to one Fe and one S atom. In the sixth O site, O is bonded in a single-bond geometry to one S atom. In the seventh O site, O is bonded in a bent 120 degrees geometry to one Fe and one S atom. In the eighth O site, O is bonded in a single-bond geometry to one S atom. In the ninth O site, O is bonded in a single-bond geometry to one S atom. In the tenth O site, O is bonded in a bent 120 degrees geometry to one Fe and one S atom. In the eleventh O site, O is bonded in a single-bond geometry to one S atom. In the twelfth O site, O is bonded in a bent 120 degrees geometry to one Fe and one S atom. In the thirteenth O site, O is bonded in a bent 120 degrees geometry to two Fe atoms. In the fourteenth O site, O is bonded in a single-bond geometry to one Fe atom. In the fifteenth O site, O is bonded in a distorted L-shaped geometry to one Fe and one O atom. The O–O bond length is 1.26 Å. In the sixteenth O site, O is bonded in a 2-coordinate geometry to one Fe and one O atom.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Co-existence of atomically dispersed Ru and Ce 3+ sites is responsible for excellent low temperature N 2 O reduction activity of Ru/CeO 2

Nitrous oxide N 2 O reduction is a big challenge due to high global warming potential of N 2 O. (~300 times higher compared with CO 2 ). The best known catalysts, such as Rh/ceria, require relatively high temperatures for N 2 O decomposition. Herein, we report that Ru/ceria catalysts with low Ru loading of ~0.25 wt% efficiently catalyze low temperature N 2 O reduction by CO starting at 100 °C (full N 2 O conversion below 200 °C) under industrially relevant flow rates and gas concentrations. Further, this remarkable performance stems from maintaining isolated Ru cations even on reduced ceria surface and, simultaneously, the propensity of Ru to affect ceria surface to form labile surface oxygen thereby creating large number of oxygen vacancies (Ce +3 cations) in the presence of CO. In contrast, for Rh/CeO 2 catalysts with equivalent metal loading, the activity is much lower because atomically dispersed Rh sinters into metallic clusters at the onset reaction temperature (~200 °C): these clusters are much less effective than isolated single Ru ions, with lower Ce +3 concentration maintained on reduced Rh/CeO 2 catalyst. Our study highlights the benefits of gaining molecular-level insight into the dynamic nature of catalytically active sites under reaction conditions for preparing catalysts containing low loading of precious metals with unsurpassed low temperature activity.

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

Breaking covalent bonds in the context of the many-body expansion (MBE). I. The purported “first row anomaly” in XH n (X = C, Si, Ge, Sn; n = 1–4)

We present a new, novel implementation of the Many-Body Expansion (MBE) to account for the breaking of covalent bonds, thus extending the range of applications from its previous popular usage in the breaking of hydrogen bonds in clusters to molecules. A central concept of the new implementation is the in situ atomic electronic state of an atom in a molecule that casts the one-body term as the energy required to promote it to that state from its ground state. The rest of the terms correspond to the individual diatomic, triatomic, etc., fragments. Its application to the atomization energies of the XH n series, X = C, Si, Ge, Sn and n = 1–4, suggests that the (negative, stabilizing) 2-B is by far the largest term in the MBE with the higher order terms oscillating between positive and negative values and decreasing dramatically in size with increasing rank of the expansion. The analysis offers an alternative explanation for the purported “first row anomaly” in the incremental H n–1 X–H bond energies seen when these energies are evaluated with respect to the lowest energy among the states of the XH n molecules. Due to the “flipping” of the ground/first excited state between CH 2 ( 3 B 1 ground state, 1 A 1 first excited state) and XH 2 , X = Si, Ge, Sn ( 1 A 1 ground state, 3 B 1 first excited state), the overall picture does not exhibit a “first row anomaly” when the incremental bond energies are evaluated with respect to the molecular states having the same in situ atomic states.

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

In Situ Studies of Copper-Zirconia-Zinc Oxide Model Catalysts for CO 2 Hydrogenation

A combination of in situ X-ray photoelectron spectroscopy (XPS) and infrared reflection absorption spectroscopy (IRAS) was used to investigate the formation of surface intermediates from CO 2 hydrogenation on copper-zirconia-zinc oxide model catalysts under reaction conditions. Copper clusters with different numbers of atoms (n = 1, 4, 13) were deposited onto bare and ZrO 2 -modified ZnO powder supports to systematically examine the effects of Cu cluster size and the synergy between metal and metal-oxide components. Under low pressure CO 2 hydrogenation conditions (CO 2 :H 2 = 1:9, 0.4 mbar, 300–600 K), XPS and IRAS identify the most prominent intermediates as carbonate (CO 3 *), formate (HCOO*) and methoxy (CH 3 O*), which is the final surface-bound intermediate leading to methanol. The temperature profiles are consistent with a mechanism in which CO 2 is adsorbed as carbonate species (HCO 3 *, CO 3 *) followed by hydrogenation reactions to formate (HCOO*) and methoxy (CH 3 O*), but the relative yields strongly depend on surface composition. Specifically, the presence of ZrO 2 promotes CO 2 adsorption and activation and improves the thermal stability of the Cu clusters against loss of surface area. Moreover, the ternary Cu 4 /ZrO 2 /ZnO surface is significantly more active than Cu 4 /ZnO and ZrO 2 /ZnO surfaces for the formation of methoxy (CH 3 O*), indicating that Cu–ZrO 2 interfaces promote the formation of key intermediates leading to methanol. Finally, the yields of intermediates are similar for all Cu cluster sizes (Cu 1 , Cu 4 and Cu 13 ), indicating that the primary role of Cu is to provide H atoms via H 2 dissociation and spillover. In conclusion, these molecular-level insights provide a fundamental understanding of the enhanced efficiency of ternary Cu–ZrO 2 –ZnO catalysts and establish design principles for developing improved catalysts for CO 2 conversion.

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

Perspective on Kramers symmetry breaking and restoration in relativistic electronic structure methods for open-shell systems

Without rigorous symmetry constraints, solutions to approximate electronic structure methods may artificially break symmetry. In the case of the relativistic electronic structure, if time-reversal symmetry is not enforced in calculations of molecules not subject to a magnetic field, it is possible to artificially break Kramers degeneracy in open shell systems. This leads to a description of excited states that may be qualitatively incorrect. Despite this, different electronic structure methods to incorporate correlation and excited states can partially restore Kramers degeneracy from a broken symmetry solution. For single-reference techniques, the inclusion of double and possibly triple excitations in the ground state provides much of the needed correction. Formally, however, this imbalanced treatment of the Kramers-paired spaces is a multi-reference problem, and so methods such as complete-active-space methods perform much better at recovering much of the correct symmetry by state averaging. Using multi-reference configuration interaction, any additional corrections can be obtained as the solution approaches the full configuration interaction limit. A recently proposed “Kramers contamination” value is also used to assess the magnitude of symmetry breaking.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗