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 739 records · Page 41

Adapting SAFT-γ perturbation theory to site-based molecular dynamics simulation. II. Confined fluids and vapor-liquid interfaces

In this work, a new classical density functional theory is developed for group-contribution equations of state (EOS). Details of implementation are demonstrated for the recently-developed SAFT-γ WCA EOS and selective applications are studied for confined fluids and vapor-liquid interfaces. The acronym WCA (Weeks-Chandler-Andersen) refers to the characterization of the reference part of the third-order thermodynamic perturbation theory applied in formulating the EOS. SAFT-γ refers to the particular form of “statistical associating fluid theory” that is applied to the fused-sphere, heteronuclear, united-atom molecular models of interest. For the monomer term, the modified fundamental measure theory is extended to WCA-spheres. A new chain functional is also introduced for fused and soft heteronuclear chains. The attractive interactions are taken into account by considering the structure of the fluid, thus elevating the theory beyond the mean field approximation. The fluctuations of energy are also included via a non-local third-order perturbation theory. The theory includes resolution of the density profiles of individual groups such as CH{sub 2} and CH{sub 3} and satisfies stoichiometric constraints for the density profiles. New molecular simulations are conducted to demonstrate the accuracy of each Helmholtz free energy contribution in reproducing the microstructure of inhomogeneous systems at the united-atom level of coarse graining. At each stage, comparisons are made to assess where the present theory stands relative to the current state of the art for studying inhomogeneous fluids. Overall, it is shown that the characteristic features of real molecular fluids are captured both qualitatively and quantitatively. For example, the average pore density deviates ∼2% from simulation data for attractive pentadecane in a 2-nm slit pore. Another example is the surface tension of ethane/heptane mixture, which deviates ∼1% from simulation data while the theory reproduces the excess accumulation of ethane at the interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomically thin interlayer phase from first principles enables defect-free incommensurate SnO2/CdTe interface

Advancing optoelectronic and emerging technologies increasingly requires control and design of interfaces between dissimilar materials. However, incommensurate interfaces are notoriously defective and rarely benefit from first-principles predictions, because no explicit atomic-structure models exist. Here, we adopt a bulk crystal structure prediction method to the interface geometry and apply it to SnO2/CdTe heterojunctions without and with the addition of CdCl2, a ubiquitous and beneficial, but abstruse processing step in CdTe photovoltaics. Whereas the direct SnO2/CdTe interface is highly defective, we discover a unique two-dimensional CdCl2 interphase, unrelated to the respective bulk structure. It facilitates a seamless transition from the rutile to zincblende lattices and removes defect-states from the interface bandgap. Implementing the predicted interface electronic structure in device simulations, we demonstrate the theoretical feasibility of bufferless oxide-CdTe heterojunction solar cells approaching the Shockley–Queisser limit. Our results highlight the broader potential of designing atomically thin interlayers to enable defect-free incommensurate interfaces.

14 SOLAR ENERGY↗

Assessing the Role of the Support in Acetic Acid Hydrodeoxygenation Selectivity on Pt/TiO2

To produce drop-in quality biofuels following the thermochemical conversion of biomass via catalytic fast pyrolysis (CFP), hydrodeoxygenation (HDO) reactions can be performed to remove excess oxygen and create a more stable bio-oil product. HDO reactions involve co-feeding the CFP vapor-phase product and H2 gas over bi-functional catalysts, such as noble-metal catalysts supported on reducible metal oxides (e.g., Pt/TiO2). Recent model-compound studies, focusing on important classes of species in the CFP vapor mixture, have sought to determine the role of the various Pt/TiO2 actives sites (i.e., Pt-metal, TiO2-support, Pt-TiO2-interfacial sites) in producing the observed desired deoxygenation and undesired decarboxylation/decarbonylation products. One important class of compounds in the CFP vapor-phase product is carboxylic acids (e.g., acetic acid). Prior experimental work on acetic acid HDO (AA-HDO) found that Pt/C and Pt/TiO2 catalysts favored the formation of undesired C-C and desired C-O bond-dissociation products; however, the fundamental surface chemistry driving this shift in selectivity has not been established. This presentation employs atomic-scale modeling to determine, through comparisons of adsorption and reaction energetics, how Pt-metal, Pt-TiO2-interface, and interfacial-vacancy sites catalyze the competing reaction pathways for AA-HDO. Our analysis indicates that hydroxyl vacancies at the Pt-TiO2 interface are critical for lowering barriers for C-O bond-cleavage steps, such that they become favorable/competitive with respect to C-C bond-dissociation steps.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Theoretical modeling and experimental observations of the atomic layer deposition of SrO using a cyclopentadienyl Sr precursor

First-principle calculations are used to model the adsorption and hydration of strontium bis(cyclopentadienyl) [Sr(Cp){sub 2}] on TiO{sub 2}-terminated strontium titanate, SrTiO{sub 3} (STO), for the deposition of strontium oxide, SrO, by atomic layer deposition (ALD). The Sr(Cp){sub 2} precursor is shown to adsorb on the TiO{sub 2}-terminated surface, with the Sr atom assuming essentially the bulk position in STO. The C–Sr bonds are weaker than in the free molecule, with a Ti atom at the surface bonding to one of the C atoms in the cyclopentadienyl rings. The surface does not need to be hydrogenated for precursor adsorption. The calculations are compared with experimental observations for a related Sr cyclopentadienyl precursor, strontium bis(triisopropylcyclopentadienyl) [Sr({sup i}Pr{sub 3}Cp){sub 2}], adsorbed on TiO{sub 2}-terminated STO. High-resolution x-ray photoelectron spectroscopy and low-energy ion scattering spectroscopy show adsorption of the Sr precursor on the TiO{sub 2}-terminated STO after a single precursor dose. This study suggests that ALD growth from the strontium precursors featuring cyclopentadienyl ligands, such as Sr(Cp){sub 2}, may initiate film growth on non-hydroxylated surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Space Plasma Ion Processing of the Lunar Soil: Modeling of Radiation-Damaged Rim Widths on Lunar Grains

Chemically and microstructurally complex altered rims around grains in the finest size fraction (<20 micron) of the lunar regolith are the result of multi-stage processes involving both solar ion radiation damage and nanoscale deposition of impact or sputter-derived vapors. The formation of the rims is an important part of the space weathering process, and is closely linked to key changes in optical reflectance and other bulk properties of the lunar surface. Recent application of field-emission scanning transmission electron microscope techniques, including energy dispersive X-ray spectral imaging, is making it easier to unravel the "nano-stratigraphy" of grain rims, and to delineate the portions of rims that represent Radiation-Amorphized (RA) host grain from overlying amorphous material that represents vapor/sputter deposits. For the portion of rims formed by host grain amorphization (henceforth called RA rims), we have been investigating the feasibility of using Monte Carlo-type ion-atom collision models, combined with experimental ion irradiation data, to derive predictive numerical models linking the width of RA rims to the grain s integrated solar ion radiation exposure time.

Chamberlin, S.↗

Bond-Dependent Thole Model for Polarizability and Spectroscopy

We present a new model for the calculation of molecular polarizabilities from effec- tive atomic polarizabilities. This model is based on the Thole modified dipole inter- action model for molecular polarizabilities, where the total polarizability is computed as a sum of effective atomic polarizabilities modified by dipole-dipole interactions. We extend this model by making the atomic polarizabilities explicit functions of the in- teratomic distances, scaling them by the radius of the volume an atom occupies in a molecule. We use the SCAN functional to show that this model, denoted TholeL, yields accurate molecular polarizabilities with little dependence on the training set. We also demonstrate that the TholeL model yields accurate polarizabilities for configurations far from the ground state structure for a wide range of molecules. Finally, we show that the TholeL model can be used to generate accurate Raman spectra for water, crystalline urea, and urea in water from ab-initio molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure of periodic crystals and quasicrystals in ultrathin films of Ba-Ti-O

Here, we model the remarkable thin-film Ba-Ti-O structures formed by heat treatment of an initial perovskite BaTiO 3 thin film on a Pt(111) surface. All structures contain a rumpled Ti-O network with all Ti threefold coordinated with O, and with Ba occupying the larger. mainly Ti 7 O 7 , pores. The quasicrystal structue is a simple decoration of three types of tiles: square, triangle and 30° rhombus, with edge lengths 6.85 Å, joined edge-to-edge in a quasicrystalline pattern; observed periodic crystals in ultrathin film Ba-Ti-O are built from these and other tiles. Simulated STM images reproduce the patterns seen experimentally, and identify the bright protrusions as Ba atoms. The models are consistent with all experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Towards Resolving the Atomization Process of an Idealized Fire Sprinkler with VOF Modeling

Fire suppression modeling depends on accurate characterization of the atomization processes of fire sprinklers.Numerical modeling, particularly the volume of fluid (VOF) method, has been applied to understanding the atomization behavior of an idealized sprinkler geometry consisting of a 9.5 mm inner-diameter cylindrical nozzle and a flat, 25.4 mm diameter disk with a liquid flow rate of 0.87 L/s. The simulations have been performed with an OpenFOAM based VOF solver, using the isoAdvector scheme for interfacial reconstruction. The sheet breakup distance and film thickness were calculated and compared with measurements from a previous study. A mesh refinement study identified the sensitivities in the predicted quantities to mesh resolution. This study enables furtherapplication of the model to simulation of the fully atomized spray.

Engineering↗

Fatty acid capped, metal oxo clusters as the smallest conceivable nanocrystal prototypes

Metal oxo clusters of the type M 6 O 4 (OH) 4 (OOCR) 12 (M = Zr or Hf) are valuable building blocks for materials science. Here, we synthesize a series of zirconium and hafnium oxo clusters with ligands that are typically used to stabilize oxide nanocrystals (fatty acids with long and/or branched chains). The fatty acid capped oxo clusters have a high solubility but do not crystallize, precluding traditional purification and single-crystal XRD analysis. We thus develop alternative purification strategies and we use X-ray total scattering and Pair Distribution Function (PDF) analysis as our main method to elucidate the structure of the cluster core. We identify the correct structure from a series of possible clusters (Zr3, Zr4, Zr6, Zr12, Zr10, and Zr26). Excellent refinements are only obtained when the ligands are part of the structure model. Further evidence for the cluster composition is provided by nuclear magnetic resonance (NMR), infrared spectroscopy (FTIR), thermogravimetry analysis (TGA), and mass spectrometry (MS). We find that hydrogen bonded carboxylic acid is an intrinsic part of the oxo cluster. Using our analytical tools, we elucidate the conversion from a Zr6 monomer to a Zr12 dimer (and vice versa), induced by carboxylate ligand exchange. Finally, we compare the catalytic performance of Zr12-oleate clusters with oleate capped, 5.5 nm zirconium oxide nanocrystals in the esterification of oleic acid with ethanol. The oxo clusters present a five times higher reaction rate, due to their higher surface area. Since the oxo clusters are the lower limit of downscaling oxide nanocrystals, we present them as appealing catalytic materials, and as atomically precise model systems. In addition, the lessons learned regarding PDF analysis are applicable to other areas of cluster science as well, from semiconductor and metal clusters, to polyoxometalates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The mechanism of the emergence of distinct overstretched DNA states

Although multiple overstretched DNA states were identified in experiments, the mechanism of the emergence of distinct states is still unclear. Molecular dynamics simulation is an ideal tool to clarify the mechanism, but the force loading rates in stretching achieved by conventional all-atom DNA models are much faster, which essentially affect overstretching states. We employed a modified coarse-grained DNA model with an unprecedented low loading rate in simulations to study the overstretching transitions of end-opened double-stranded DNA. We observed two-strand peeling off for DNA with low stability and the S-DNA with high stability under tension. By introducing a melting-forbidden model which prevents base-pair breaking, we still observed the overstretching transition induced by the formation of S-DNA due to the change of dihedral angle. Hence, we confirmed that the competition between the two strain-softening manners, i.e., base-pair breaking and dihedral angle variation, results in the emergence of distinct overstretched DNA states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calcium kinetics with microgram stable isotope doses and saliva sampling

Studies of calcium kinetics require administration of tracer doses of calcium and subsequent repeated sampling of biological fluids. This study was designed to develop techniques that would allow estimation of calcium kinetics by using small (micrograms) doses of isotopes instead of the more common large (mg) doses to minimize tracer perturbation of the system and reduce cost, and to explore the use of saliva sampling as an alternative to blood sampling. Subjects received an oral dose (133 micrograms) of 43Ca and an i.v. dose (7.7 micrograms) of 46Ca. Isotopic enrichment in blood, urine, saliva and feces was well above thermal ionization mass spectrometry measurement precision up to 170 h after dosing. Fractional calcium absorptions determined from isotopic ratios in blood, urine and saliva were similar. Compartmental modeling revealed that kinetic parameters determined from serum or saliva data were similar, decreasing the necessity for blood samples. It is concluded from these results that calcium kinetics can be assessed with micrograms doses of stable isotopes, thereby reducing tracer costs and with saliva samples, thereby reducing the amount of blood needed.

Clinical Trial↗

One-dimensional time-dependent fluid model of a very high density low-pressure inductively coupled plasma

A time-dependent two-fluid model has been developed to understand axial variations in the plasma parameters in a very high density (peak n e~ > 5x10 19 m –3 ) argon inductively coupled discharge in a long 1.1 cm radius tube. The model equations are written in 1D, with radial losses to the tube walls accounted for by the inclusion of effective particle and energy sink terms. The ambipolar diffusion equation and electron energy equation are solved to find the electron density n e (z,t) and temperature T e (z,t), and the populations of the neutral argon 4s metastable, 4s resonant, and 4p excited state manifolds are calculated in order to determine the stepwise ionization rate and calculate radiative energy losses. The model has been validated through comparisons with Langmuir probe ion saturation current measurements; close agreement between the simulated and measured axial plasma density profiles and the initial density rise rate at each location was obtained at p Ar = 30-60 mTorr. Lastly, we present detailed results from calculations at 60 mTorr, including the time-dependent electron temperature, excited state populations, and energy budget within and downstream of the radiofrequency (RF) antenna.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ATOMIC DATA OF Cu I FOR THE INVESTIGATION OF ELEMENT ABUNDANCE

A complete set of atomic data of Cu I, including the energy levels, oscillator strengths, and photoionization cross sections, is theoretically studied to investigate element abundance including nonlocal thermodynamic equilibrium (NLTE) effects. The calculations are carried out by using the R-matrix method in the LS-coupling scheme. Twenty terms of Cu II are utilized as target states, and extensive configuration interactions are included to properly delineate the quantum states of Cu II and Cu I. One hundred thirteen bound states and 1699 oscillator strengths for E1 transitions between these states are obtained. Photoionization cross sections for all bound states are calculated in a photon energy range covering 1.28 Ry from the threshold of the respective state. Resonances shown in the photoionization cross sections are identified, and some strong resonances are expected to play an important role in NLTE modeling. The atomic data in this work represent the first complete data set for copper abundance studies. Our results are compared with the experimental and other theoretical data wherever available.

74 ATOMIC AND MOLECULAR PHYSICS↗

Time-dependent restricted-active-space self-consistent-field singles method for many-electron dynamics

The time-dependent restricted-active-space self-consistent-field singles (TD-RASSCF-S) method is presented for investigating TD many-electron dynamics in atoms and molecules. Adopting the SCF notion from the muticonfigurational TD Hartree-Fock (MCTDHF) method and the RAS scheme (single-orbital excitation concept) from the TD configuration-interaction singles (TDCIS) method, the TD-RASSCF-S method can be regarded as a hybrid of them. We prove that, for closed-shell N{sub e}-electron systems, the TD-RASSCF-S wave function can be fully converged using only N{sub e}/2 + 1 ⩽ M ⩽ N{sub e} spatial orbitals. Importantly, based on the TD variational principle, the converged TD-RASSCF-S wave function with M = N{sub e} is more accurate than the TDCIS wave function. The accuracy of the TD-RASSCF-S approach over the TDCIS is illustrated by the calculation of high-order harmonic generation spectra for one-dimensional models of atomic helium, beryllium, and carbon in an intense laser pulse. The electronic dynamics during the process is investigated by analyzing the behavior of electron density and orbitals. The TD-RASSCF-S method is accurate, numerically tractable, and applicable for large systems beyond the capability of the MCTDHF method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Graph theory approach to determine configurations of multidentate and high coverage adsorbates for heterogeneous catalysis

Abstract Heterogeneous catalysts constitute a crucial component of many industrial processes, and to gain an understanding of the atomic-scale features of such catalysts, ab initio density functional theory is widely employed. Recently, growing computational power has permitted the extension of such studies to complex reaction networks involving either high adsorbate coverages or multidentate adsorbates, which bind to the surface through multiple atoms. Describing all possible adsorbate configurations for such systems, however, is often not possible based on chemical intuition alone. To systematically treat such complexities, we present a generalized Python-based graph theory approach to convert atomic scale models into undirected graph representations. These representations, when combined with workflows such as evolutionary algorithms, can systematically generate high coverage adsorbate models and classify unique minimum energy multidentate adsorbate configurations for surfaces of low symmetry, including multi-elemental alloy surfaces, steps, and kinks. Two case studies are presented which demonstrate these capabilities; first, an analysis of a coverage-dependent phase diagram of absorbate NO on the Pt 3 Sn(111) terrace surface, and second, an investigation of adsorption energies, together with identifying unique minimum energy configurations, for the reaction intermediate propyne (CHCCH 3 *) adsorbed on a PdIn(021) step surface. The evolutionary algorithm approach reproduces high coverage configurations of NO on Pt 3 Sn(111) using only 15% of the number of simulations required for a brute force approach. Furthermore, the screening of potentially hundreds of multidentate adsorbates is shown to be possible without human intervention. The strategy presented is quite general and can be applied to a spectrum of complex atomic systems.

36 MATERIALS SCIENCE↗

DECT evaluation of noncalcified coronary artery plaque

Purpose: Composition of the coronary artery plaque is known to have critical role in heart attack. While calcified plaque can easily be diagnosed by conventional CT, it fails to distinguish between fibrous and lipid rich plaques. In the present paper, the authors discuss the experimental techniques and obtain a numerical algorithm by which the electron density (ρ{sub e}) and the effective atomic number (Z{sub eff}) can be obtained from the dual energy computed tomography (DECT) data. The idea is to use this inversion method to characterize and distinguish between the lipid and fibrous coronary artery plaques. Methods: For the purpose of calibration of the CT machine, the authors prepare aqueous samples whose calculated values of (ρ{sub e}, Z{sub eff}) lie in the range of (2.65 × 10{sup 23} ≤ ρ{sub e} ≤ 3.64 × 10{sup 23}/cm{sup 3}) and (6.80 ≤ Z{sub eff} ≤ 8.90). The authors fill the phantom with these known samples and experimentally determine HU(V{sub 1}) and HU(V{sub 2}), with V{sub 1},V{sub 2} = 100 and 140 kVp, for the same pixels and thus determine the coefficients of inversion that allow us to determine (ρ{sub e}, Z{sub eff}) from the DECT data. The HU(100) and HU(140) for the coronary artery plaque are obtained by filling the channel of the coronary artery with a viscous solution of methyl cellulose in water, containing 2% contrast. These (ρ{sub e}, Z{sub eff}) values of the coronary artery plaque are used for their characterization on the basis of theoretical models of atomic compositions of the plaque materials. These results are compared with histopathological report. Results: The authors find that the calibration gives ρ{sub e} with an accuracy of ±3.5% while Z{sub eff} is found within ±1% of the actual value, the confidence being 95%. The HU(100) and HU(140) are found to be considerably different for the same plaque at the same position and there is a linear trend between these two HU values. It is noted that pure lipid type plaques are practically nonexistent, and microcalcification, as observed in histopathology, has to be taken into account to explain the nature of the observed (ρ{sub e}, Z{sub eff}) data. This also enables us to judge the composition of the plaque in terms of basic model which considers the plaque to be composed of fibres, lipids, and microcalcification. Conclusions: This simple and reliable method has the potential as an effective modality to investigate the composition of noncalcified coronary artery plaques and thus help in their characterization. In this inversion method, (ρ{sub e}, Z{sub eff}) of the scanned sample can be found by eliminating the effects of the CT machine and also by ensuring that the determination of the two unknowns (ρ{sub e}, Z{sub eff}) does not interfere with each other and the nature of the plaque can be identified in terms of a three component model.

60 APPLIED LIFE SCIENCES↗

Permeation Properties of Disordered Metal-Organic Framework Membranes Made by Vapor Phase Ligand Treatment

Zeolitic imidazolate frameworks (ZIFs) are nanoporous molecular sieves offering exciting opportunities for membrane- and adsorption-based gas separations. The Tsapatsis group developed an all-vapor-phase ligand induced permselectivation (LIPS) method for the fabrication of ZIF nanocomposite membranes. The LIPS method consists of a combination of atomic layer deposition of a dense oxide inside the mesopores of a porous support followed by transformation of the dense oxide deposit to nanoporous ZIF by exposure to sublimated vapors of an imidazolate ligand. It enables the formation of thin nanocomposite films consisting of ZIF deposited inside pores with diameters smaller than 10 nm. Vapor-phase conversion of the impermeable ZnO to ZIF increases flux of certain gases more than others. Vapor phase linker/ligand treatment (VPLT) can further modify the properties of LIPS and other membranes. The combination of LIPS with VPLT provides an all-vapor methodology for making and modifying ZIF membranes with the ability to control the ZIF deposit composition and the level of confinement in a variety of mesoporous supports. The propylene/propane separation performance of LIPS/VPLT ZIF membranes confined (entirely or partially) inside mesopores is remarkably superior to that of most conventional ZIF membranes consisting of unconfined (deposited on the external surface of porous supports) polycrystalline films. The accomplishments of this award are: (i) understanding of the LIPS process using mathematical models of atomic layer deposition (ALD), (ii) development of surface modification methods (vapor, e-beam, plasma, and X-ray treatments) to tune permeation properties, (iii) microstructural control of ZIF-8 films using ALD (support modification, direct ALD, and by use of in situ permeability monitoring), (iv) molecular simulations, and (v) thin film characterization and adsorption measurements using spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomic protein structure refinement using all-atom graph representations and SE(3)-equivariant graph transformer

Abstract Motivation The state-of-art protein structure prediction methods such as AlphaFold are being widely used to predict structures of uncharacterized proteins in biomedical research. There is a significant need to further improve the quality and nativeness of the predicted structures to enhance their usability. In this work, we develop ATOMRefine, a deep learning-based, end-to-end, all-atom protein structural model refinement method. It uses a SE(3)-equivariant graph transformer network to directly refine protein atomic coordinates in a predicted tertiary structure represented as a molecular graph. Results The method is first trained and tested on the structural models in AlphaFoldDB whose experimental structures are known, and then blindly tested on 69 CASP14 regular targets and 7 CASP14 refinement targets. ATOMRefine improves the quality of both backbone atoms and all-atom conformation of the initial structural models generated by AlphaFold. It also performs better than two state-of-the-art refinement methods in multiple evaluation metrics including an all-atom model quality score—the MolProbity score based on the analysis of all-atom contacts, bond length, atom clashes, torsion angles, and side-chain rotamers. As ATOMRefine can refine a protein structure quickly, it provides a viable, fast solution for improving protein geometry and fixing structural errors of predicted structures through direct coordinate refinement. Availability and implementation The source code of ATOMRefine is available in the GitHub repository (https://github.com/BioinfoMachineLearning/ATOMRefine). All the required data for training and testing are available at https://doi.org/10.5281/zenodo.6944368.

59 BASIC BIOLOGICAL SCIENCES↗