Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic & molecular 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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Predicting and Understanding Non-Covalent Interactions Using Novel Forms of Symmetry-Adapted Perturbation Theory

Although sometimes derided as “weak” interactions, non-covalent forces play a critical role in ligand binding and crystal packing and in determining the conformational landscape of flexible molecules. Symmetry-adapted perturbation theory (SAPT) provides a framework for accurate ab initio calculation of intermolecular interactions and furnishes a natural decomposition of the interaction energy into physically meaningful components: semiclassical electrostatics (rigorously obtained from monomer charge densities), Pauli or steric repulsion, induction (including both polarization and charge transfer), and dispersion. This decomposition helps to foster deeper understanding of non-covalent interactions and can be used to construct transferable, physics-based force fields. Separability of the SAPT interaction energy also provides the flexibility to construct composite methods, a feature that we exploit to improve the description of dispersion interactions. Furthermore, these are challenging to describe accurately because they arise from nonlocal electron correlation effects that appear for the first time at second order in perturbation theory but are not quantitatively described at that level.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Toward Automated Sampling of Polymorph Nucleation and Free Energies with the SGOOP and Metadynamics

Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the “spectral gap optimization of order parameters (SGOOP)” approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. Here, we believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.

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

Probing High-Pressure Structural Evolution in Polyurea with In Situ Energy-Dispersive X-ray Diffraction and Molecular Dynamics Simulations

Polyurea, an elastomer with a phase-segregated microstructure, has been proven as an effective coating in defense applications. To gain a more complete understanding of the high-pressure atomic-level morphology of these phases and to validate molecular dynamics (MD) simulations, multi-angle energy-dispersive X-ray diffraction experiments were performed in situ up to pressures of ~6 GPa at room temperature. Structure factors were obtained and compared to MD simulations with an average error of less than 5% between major peak positions. The first sharp diffraction peak shifted from 4.56 Å to lower d-spacing with pressure, indicating compression between hard segments. This was further supported by the behavior of a peak at ~3.86 Å from the pair distribution function (PDF), suspected to represent π-stacking and separation between soft segments. Compression within the hard segments themselves is minimal as low-r peaks in the PDF are not greatly affected by pressure.

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

Epitaxial (Al x Ga 1-x-y In y ) 2 O 3 Alloys Lattice Matched to Monoclinic Ga 2 O 3 Substrates

We have epitaxially stabilized a series of monoclinic (Al x Ga 1-x-y In y ) 2 O 3 alloys by careful choice of molecular beam epitaxy growth conditions, which balance alloy growth with suboxide desorption. The films are pseudomorphic to (010) β-Ga 2 O 3 substrates at thicknesses up to 150 nm with compositions ranging from (Al 0.01 Ga 0.83 In 0.16 ) 2 O 3 to (Al 0.24 Ga 0.75 In 0.03 ) 2 O 3 . The absorption edge shifts from approximately 4.62-5.14 eV with coincidently increasing Al and decreasing In mole fractions. J-V measurements reveal an increase in resistivity over four orders of magnitude with a maximum value of 4.2 x 10 5 Ω-cm for (Al 0.17 Ga 0.76 In 0.07 ) 2 O 3 , which has nearly identical lattice parameters (both in-plane and out-of-plane) to the underlying β-Ga 2 O 3 . Scanning transmission electron microscopy of this sample reveals a mostly uniform and single crystalline film, though we identify areas of non-uniform In incorporation and some γ-phase inclusions. This work demonstrates the feasibility of thick layers lattice-matched to β-Ga 2 O 3 with increased bandgap compared to phase-separation limited (Al,Ga) 2 O 3 . These alloys can enable higher bandgap epitaxial dielectrics and high sheet charge density transistors by increasing the conduction band offset with respect to β-Ga 2 O 3 .

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

Coupled Cluster Green's function formulations based on the effective Hamiltonians

In this work, we demonstrate that the effective Hamiltonians obtained with the downfolding procedure based on double unitary coupled cluster (DUCC) ansatz can be used in the context of Green’s function coupled cluster (GFCC) formalism to calculate spectral functions of molecular systems. This combined approach (DUCC-GFCC) provides a significant reduction of numerical effort and good agreement with the corresponding all-orbital GFCC methods in energy windows that are consistent with the choice of active space. These features are demonstrated on the example of two benchmark systems: H 2 O and N 2 , where DUCC-GFCC calculations were performed for active spaces of various sizes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Generalized representative structures for atomistic systems

A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.

atomic cluster expansion↗
DOE OSTI2023

Non-Empirical and Self-Interaction Corrections for DFTB: Towards Accurate Quantum Simulations for Large Mesoscale Systems (Final Report)

This project was comprised of two complementary (but parallel) thrusts: (1) implementing massively-parallelized computing hardware (with new computational hardware that may replace GPUs) in the density functional tight binding (DFTB) approach and (2) developing new capabilities in DFTB to calculate the electronic structure and dynamics of large chemical systems. While classical molecular dynamics can handle hundreds of thousands of atoms, it cannot provide a first-principles based description of chemical systems at the quantum level. At the other extreme, conventional Kohn-Sham DFT methods can probe the true quantum mechanical nature of chemical systems; however, these methods cannot tackle the large sizes and length scales relevant to dynamics simulations of realistic systems. The DFTB formalism utilized in this project provides a viable approach for probing these large systems at a quantum mechanical level of detail. However, to utilize the DFTB approach for accurate calculations of electronic properties, it is crucial to incorporate quantum-based corrections in DFTB since exchange-correlation effects can still remain very strong in these large systems. At the same time, enhancing the computational efficiency of DFTB is also essential since optimal computational performance is required for addressing the large size scales associated with realistic chemical systems. As such, the new non-empirical corrections and computing hardware enhancements implemented in this project will enable accurate and computationally efficient approaches to directly probe electronic properties in these large, complex systems.

36 MATERIALS SCIENCE↗
NASA NTRS1981

Cluster adsorption on amorphous and crystalline surfaces - A molecular dynamics study of model Pt on Cu and model Pd on Pt

Molecular dynamics was used to study the structure, dispersion and short-time behavior of ten-atom clusters adsorbed onto amorphous and crystalline substrates, in which the cluster atoms differed from the substrate atoms. Two adatom-substrate model systems were chosen; one, in which the interaction energy between adatom pairs was greater than that between substrate pairs, and the other, in which the reverse was true. At relatively low temperature ranges, increased dispersion of cluster atoms occurred: (a) on the amorphous substrate as compared to the FCC(100) surface, (b) with increasing reduced temperature, and (c) with adatom-substrate interaction energy stronger than adatom-adatom interaction. Two-dimensional clusters (rafts) on the FCC(100) surface displayed migration of edge atoms only, indicating a mechanism for the cluster rotation and shape changes found in experimental studies.

Garofalini, S. H.↗
NASA NTRS1985

Spectral identification/elimination of molecular species in spacecraft glow

Computer models of molecular electronic and vibrational emission intensities were developed. Known radiative emission rates (Einstein coefficients) permit the determination of relative excited state densities from spectral intensities. These codes were applied to the published spectra of glow above shuttle surface and to the Spacelab 1 results of Torr and Torr. The theoretical high-resolution spectra were convolved with the appropriate instrumental slit functions to allow accurate comparison with data. The published spacelab spectrum is complex but N2+ Meinel emission can be clearly identified in the ram spectrum. M2 First Positive emission does not correlate well with observed features, nor does the CN Red System. Spectral overlay comparisons are presented. The spectrum of glow above shuttle surfaces, in contrast to the ISO data, is not highly structured. Diatomic molecular emission was matched to the observed spectral shape. Source excitation mechanisms such as (oxygen atom)-(surface species) reaction product chemiluminescence, surface recombination, or resonance fluorescent re-emission will be discussed for each tentative assignment. These assignments are the necessary first analytical step toward mechanism identification. Different glow mechanisms will occur above surfaces under different orbital conditions.

Green, B. D.↗
NASA NTRS1997

In Situ identification of mineral resources with an X-ray-optical "Hand-Lens" instrument

The recognition of material resources on a planetary surface requires exploration strategies not dissimilar to those employed by early field geologists who searched for ore deposits primarily from surface clues. In order to determine the location of mineral ores or other materials, it will be necessary to characterize host terranes at regional or subregional scales. This requires geographically broad surveys in which statistically significant numbers of samples are rapidly scanned from a roving platform. To enable broad-scale, yet power-conservative planetary-surface exploration, we are developing an instrument that combines x-ray diffractometry (XRD), x-ray fluorescence spectrometry (XRF), and optical capabilities; the instrument can be deployed at the end of a rover's robotic arm, without the need for sample capture or preparation. The instrument provides XRD data for identification of mineral species and lithological types; diffractometry of minerals is conducted by ascertaining the characteristic lattice parameters or "d-spacings" of mineral compounds. D-spacings of 1.4 to 25 angstroms can be determined to include the large molecular structures of hydrated minerals such as clays. The XRF data will identify elements ranging from carbon (Atomic Number = 6) to elements as heavy as barium (Atomic Number = 56).

Marshall, J.↗
NASA NTRS1999

In Situ Identification of Mineral Resources with an X-Ray-Optical "Hands-Lens" Instrument

The recognition of material resources on a planetary surface requires exploration strategies not dissimilar to those employed by early field geologists who searched for ore deposits primarily from surface clues. In order to determine the location of mineral ores or other materials, it will be necessary to characterize host terranes at regional or subregional scales. This requires geographically broad surveys in which statistically significant numbers of samples are rapidly scanned from a roving platform. To enable broad-scale, yet power-conservative planetary-surface exploration, we are developing an instrument that combines x-ray diffractometry (XRD), x-ray fluorescence spectrometry (XRF), and optical capabilities; the instrument can be deployed at the end of a rover's robotic arm, without the need for sample capture or preparation. The instrument provides XRD data for identification of mineral species and lithological types; diffractometry of minerals is conducted by ascertaining the characteristic lattice parameters or "d-spacings" of mineral compounds. D-spacings of 1.4 to 25 angstroms can be determined to include the large molecular structures of hydrated minerals such as clays. The XRF data will identify elements ranging from carbon (Atomic Number = 6) to elements as heavy as barium (Atomic Number = 56). While a sample is being x-rayed, the instrument simultaneously acquires an optical image of the sample surface at magnifications from lx to at least 50x (200x being feasible, depending on the sample surface). We believe that imaging the sample is extremely important as corroborative sample-identification data (the need for this capability having been illustrated by the experience of the Pathfinder rover). Very few geologists would rely on instrument data for sample identification without having seen the sample. Visual inspection provides critical recognition data such as texture, crystallinity, granularity, porosity, vesicularity, color, lustre, opacity, and so forth. These data can immediately distinguish sedimentary from igneous rocks, for example, and can thus eliminate geochemical or mineral ambiguities arising, say between arkose and granite. It would be important to know if the clay being analyzed was part of a uniform varve deposit laid down in a quiescent lake, or the matrix of a megabreccia diamictite deposited as a catastrophic impact ejecta blanket. The unique design of the instrument, which combines Debye-Scherrer geometry with elements of standard goniometry, negates the need for sample preparation of any kind, and thus negates the need for power-hungry and mechanically-complex sampling systems that would have to chip, crush, sieve, and mount the sample for x-ray analysis. Instead, the instrument is simply rested on the sample surface of interest (like a hand lens); the device can interrogate rough rock surfaces, coarse granular material, or fine rock flour. A breadboard version of the instrument has been deployed from the robotic arm of the Marsokhod rover in field trials at NASA Ames, where large vesicular boulders were x-rayed to demonstrate the functionality of the instrument design, and the ability of such a device to comply with constraints imposed by a roving platform. Currently under development is a flight prototype concept of this instrument that will weigh 0.3 kg, using about 4500 J of energy per sample analysis. It requires about 5 min. for XRD analysis, and about 30 min. for XRF interrogation. Its small mass and rugged design make it ideal for deployment on small rovers of the type currently envisaged for the exploration of Mars (e.g., Sojourner-scale platforms). The design utilizes a monolithic P-N junction photodiode pixel array for XRD, a Si PIN photodiode/avalanche photodiode system for XRF, and an endoscopic imaging camera system unobtrusively embedded between the detectors and the x-ray source (the endoscope with its board-mounted camera can be adapted for IR light in addition to visible wavelenths. A rugged, miniature (7 cu cm) x-ray source for the instrument has already been breadboarded.

Marshall, J.↗
DOE OSTI2023

Catching the Killer: Dynamic Disorder Design Rules for Small–Molecule Organic Semiconductors

Small-molecule organic semiconductors are promising materials for applications ranging from solar cells to medical sensors. Their nearly infinite design space means that, in theory, it is possible to tailor the material to the exact specifications of a particular application. In reality; however, design rules to improve mobility (μ) remain elusive because of the complex calculations required to understand its limiters. Transient localization theory posits that charge carriers are slowed down by the collective phonon motions, called dynamic disorder, which localize charge carriers temporarily. Traditionally, a mode-by-mode analysis is performed to try and identify “killer” modes. Herein, the flaws are demonstrated with this analysis and present a streamlined simulation workflow that simplifies simulation of dynamic disorder and enables engineering of new molecular structures based on phonon dynamics. This workflow is combined with a novel visualization technique that enables per-atom insights into how μ is limited. Furthermore, this workflow is applied and analyzed to a series of -acenes and a series of BTBT-based molecules. Then design rules are identified in each series that identify possible ways to improve the μ beyond current experimental limits using phonon engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Single-atom Zr promoter boosts oxygen activation on ceria-supported Pt catalysts

Activation of surface lattice oxygen and chemisorbed oxygen on catalyst surfaces constitutes a pivotal step in heterogeneous oxidative catalysis. Herein, we report a strategy for enhancing oxygen activation by rational design of catalysts with single-atom promoters. Single-site Zr species in CeO 2 (Zr 1 -CeO 2 ) are synthesized using the atom-trapping method. The Zr 1 -CeO 2 -supported Pt catalyst exhibits enhanced catalytic performance over the CeO 2 -supported Pt catalyst in the oxidation of CO, C 3 H 8 , and C 3 H 6 , achieving significantly lower T 50 values (temperature required to reach 50% conversion). This enhanced catalytic activity is attributed to the formation of an asymmetric Zr 1 -O-Pt 1 structure, which favors the activation of the adjacent surface lattice oxygen and chemisorbed molecular oxygen. This work exemplifies that incorporating single-site atoms into oxide support facilitates oxygen activation, providing new insights into the role of atomically dispersed promoters in heterogeneous catalysis.

Science & Technology - Other Topics↗
DOE OSTI2020

Chemical and steric effects in simulating noncontact atomic force microscopy images of organic molecules on a Cu (111) substrate

Noncontact mode of atomic force microscopy (nc-AFM) employing a CO-functionalized tip is a very powerful tool for studying molecular structures. However, interpreting nc-AFM images for nonplanar molecules can sometimes be problematic. To illustrate and resolve the nature of such problematic systems, we employ real-space pseudopotentials constructed within density functional theory to simulate nc-AFM images. In this work, we focus on several representative nonplanar organic molecules (pentacene, naphthanthrone, olympicene, and 6-phenylhexa-1,3,5-triynybenzene (PHTB)) on a typical substrate: the Cu (111) surface. This substrate results in significant distortions in the molecular geometries of pentacene and naphthanthrone. Including these distortions in simulated nc-AFM imaging notably improves the agreement between the simulated and measured images. In naphthanthrone, the relatively large interaction between the O atom and the Cu substrate offers a straightforward explanation for the absence of the C = O bond in the measured image. Nonplanar features such as tilting or twisting are also apparent in olympicene and PHTB. A “triangular” bright feature associated with the –CH 2 group in olympicene appears in simulated and measured nc-AFM images. This feature is directly related to the tilting angle of the molecule with respect to the Cu substrate. The “defective” benzene ring feature and the faint ellipsoidal C ≡ C feature in PHTB can be ascribed to its twisted nature. The ability to simulate such subatomic images in nc-AFM reflects the accuracy and efficiency of calculating quantum forces in real space.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of the bonding in tetrahedrane and phosphorus-substituted tetrahedranes

The bonding structures of tetrahedrane, phosphatetrahedrane, diphosphatetrahedrane and triphosphatetrahedrane are studied by employing an intrinsic quasi-atomic orbital analysis. Ethane, cyclopropane and tetrahedral P 4 are employed as reference systems. The orbital analysis is paired with the computation of strain energies via isodesmic reactions. The results show that the increase in geometric strain upon transition from ethane to cyclopropane to tetrahedrane weakens the CC bonds, despite leading to shorter C–C interatomic distances. With the increase in strain, the orbitals centered on C and involved in the bonding of the cage structure are observed to have elevated p-character, and the orbital structure of C deviates from sp 3 hybridization. The systematic substitution of CH groups by P atoms in the cage structure of tetrahedrane leads to stronger CC bonds, larger angles in the cage structures of the resulting phosphatetrahedranes, lower p-character in the orbitals involved in the bonding of the cages, and lower strain energies. It is found that P is more amenable to strained molecular arrangements than is C, and that the propensity of a given atom to hybridize s and p functions, or the lack thereof, has implications in the stability of molecules with strained geometries. The combination of the calculations presented here with the existing literature provides insight into the apparent propensity of tetrahedrane and P 4 to ‘break’ their tetrahedral structures. Trends in the bonding interactions, such as bond strengths, s- and p-orbital characters and charge transfer are identified and related to the strain energy observed in each of the analyzed systems.

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

Efficient analysis of small-angle scattering curves for large biomolecular assemblies using Monte Carlo methods

Structure elucidation from small-angle scattering curves of large biomolecular assemblies is notoriously challenging. This is because the simulation of high-resolution features in the structure of large macromolecular assemblies, such as de novo protein assemblies, is computationally demanding when it needs to cover a broad range of length scales. Conventional methods, such as the numerical approximation to the Debye equation or the use of spherical harmonics, do not scale well as the size of the assembly increases, which limits their application to small structures (e.g. individual proteins). This work explores the effectiveness of a Monte Carlo method to simulate and fit scattering curves for large biomolecular assemblies spanning over ranges covering atomic and molecular detail (e.g. spacing and orientation of proteins in an assembly) as well as large-scale (hundreds of nanometres) features. Owing to its speed and scalability, it can be combined with a fitting algorithm to extract structural features from experimental small-angle scattering curves in biomolecular assemblies that are otherwise intractable for interpretation. This work first demonstrates the effectiveness of the tool using experimental small-angle X-ray scattering (SAXS) data from tile-like proteins that assemble into 1D tube-like macromolecular structures. Here, the diameter distribution of tubes is extracted from SAXS fits, and this is quantitatively compared with distributions from electron microscopy. SAXS data are also obtained from 2D sheet-like protein assemblies, and the proposed method is used to quantify structural features such as the separation distance between protein building blocks and the flexing of the sheet. An open-source implementation of the methodology is provided for use in a broad range of biological systems involving multi-scale scattering analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗