Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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 127 records · Page 7

Ab Initio -Based Bond Order Potential to Investigate Low Thermal Conductivity of Stanene Nanostructures

We introduce a bond order potential (BOP) for stanene based on an ab initio derived training data set. The potential is optimized to accurately describe the energetics, as well as thermal and mechanical properties of a free-standing sheet, and used to study diverse nanostructures of stanene, including tubes and ribbons. As a representative case study, using the potential, we perform molecular dynamics simulations to study stanene’s structure and temperature-dependent thermal conductivity. We find that the structure of stanene is highly rippled, far in excess of other 2-D materials (e.g., graphene), owing to its low in-plane stiffness (stanene: ~ 25 N/m; graphene: ~ 480 N/ m). The extent of stanene’s rippling also shows stronger temperature dependence compared to that in graphene. Furthermore, we find that stanene based nanostructures have significantly lower thermal conductivity compared to graphene based structures owing to their softness (i.e., low phonon group velocities) and high anharmonic response. Our newly developed BOP will facilitate the exploration of stanene based low dimensional heterostructures for thermoelectric and thermal management applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Isolation and electronic structures of derivatized manganocene, ferrocene and cobaltocene anions

The discovery of ferrocene nearly 70 years ago marked the genesis of metallocene chemistry. Although the ferrocenium cation was discovered soon afterwards, a derivatized ferrocenium dication was only isolated in 2016 and the monoanion of ferrocene has only been observed in low-temperature electrochemical studies. Here we report the isolation of a derivatized ferrocene anion in the solid state as part of an isostructural family of 3d metallocenates, which consist of anionic complexes of a metal centre (manganese, iron or cobalt) sandwiched between two bulky Cp ttt ligands (where Cp ttt is {1,2,4-C 5 H 2 t Bu 3 }). These thermally and air-sensitive complexes decompose rapidly above –30 °C; however, we were able to characterize all metallocenates by a wide range of physical techniques and ab initio calculations. These data have allowed us to map the electronic structures of this metallocenate family, including an unexpected high-spin S = 3/2 ground state for the 19e– derivatized ferrocene anion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafast flow chemistry for the acid-catalyzed conversion of fructose

Biomass conversion to chemicals and fuels has been epitomized by batch processing and long reaction times that limit the economic viability of future biorefineries and prevent modular manufacturing near the source-production site that is necessary due to the biomass’ large water content. Herein we introduce a continuous flow microreactor to enable process intensification and one thousand-fold device miniaturization for remote chemical processing. We demonstrate the approach in HCl-catalyzed fructose dehydration in water. We characterize and quantify the mixing using laser induced fluorescence (LIF) of fluorescent dyes and particle image velocimetry (PIV) of fluorescent microspheres. The estimated mixing times ranged from 0.03 to 4.8 s for residence times from 1 to 120 s, respectively. Moreover, a curved channel geometry induces secondary (Dean) vortices that produce a three-fold increase in mixing by diffusion through advection along the cross-sectional direction of the channel, creating dissipative chaotic mixing under laminar flow. We then obtain for the first time isothermal kinetics of fructose dehydration at short contact times and high temperatures and compare to a previously published, hybrid ab initio and data-driven kinetic model. Rigorous model optimization is carried out and an optimal HMF yield of 54% is attained in a single aqueous phase. High temperatures and low pH improve the reactor throughput by increasing the HMF space-time yield. We demonstrate the highest ever productivity (HMF yield per unit time). Importantly, this can be achieved at ~100% fructose conversion, eliminating the need for recycling to improve economics. Furthermore, optimization of the reactor dimensions and the specific energy dissipation rate results in ~60% pressure-drop related energy savings. We describe a methodology to scale-up microreactors to millimeter size and increase productivity 16-fold under optimal conditions. The energy efficiency and cost analysis for HMF production from corn stover of a typical corn farm indicate an optimal strategy for numbering-up millimeter size reactors and, importantly, they underscore the potential feasibility of producing bioproducts in a distributed manner.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Internal conversion and intersystem crossing dynamics based on coupled potential energy surfaces with full geometry-dependent spin–orbit and derivative couplings. Nonadiabatic photodissociation dynamics of NH 3 (A) leading to the NH(X 3 Σ – , a 1 Δ) + H 2 channel

Here we simulate the photodissociation of NH 3 originating from its first excited singlet state S 1 into the NH 2 + H (radical) and NH + H 2 (molecular) channels. The states considered are the ground singlet state S 0 , the first excited singlet state S 1 and the lowest-lying triplet state T 1 , which permit for the first time a uniform treatment of the internal conversion and intersystem crossing. The simulations are based on a diabatic potential energy matrix (DPEM) of S 0 , S 1 coupled by a conical intersection seam, as well as a potential energy surface (PES) for T 1 coupled by spin-orbit coupling (SOC) to the two singlet states. The DPEM and PES are fitted to ab initio electronic structure data (ESD) including energies, energy gradients, and derivative couplings. The DPEM also defines an adiabatic to diabatic state (AtD) transformation, which is used to transform the singular adiabatic SOC into a smooth function of the nuclear coordinates in the diabatic representation, allowing the diabatic SOC to be fit to an analytical functional form. ESD and SOC data obtained from these surfaces can serve as input for either quantum or semi-classical characterization of the nonadiabatic dynamics. Using the SHARC suite of programs, nonadiabatic simulations based on over 40 000 semi-classical trajectories assess the convergence of our results. The production of NH + H 2 is not direct, but is only achieved through a quasi-statistical dissociation mechanism after internal conversion to the ground electronic state. This leads to a much lower yield comparing with the main NH 2 + H channel. The NH(X 3 Σ_) radical produced through the intersystem crossing from S 0 to T 1 is rare (~0.2%) compared to NH(a 1 Δ) due to the process being spin forbidden.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio calculations on the magnetic properties of transition metal complexes

We present a protocol for the ab initio determination of the magnetic properties of mono- and polynuclear transition metal compounds. First, we obtain the low lying electronic states by multireference methods. Then, we include spin-orbit coupling and an external magnetic field for the determination of zero-field splitting and g-tensors. For the polynuclear complexes the magnetic exchange coupling constants are determined by a modified complete active space self consistent field method. Based on the results of the ab initio calculations, magnetic data such as magnetic susceptibility or magnetization are simulated and compared to experimental data. The results obtained for the polynuclear complexes are further analysed by calculations on model complexes where part of the magnetic centers are substituted by diamagnetic ions. The methods are applied to different Co and Ni containing transition metal complexes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Structural and optical properties of the naked and passivated Al{sub 5}Au{sub 5} bimetallic nanoclusters

The structural and optical properties of both the naked and passivated bimetallic Al{sub 5}Au{sub 5} nanoclusters have been analyzed based on data obtained from ab initio density functional theory and quantum molecular dynamics simulations. It has been found that the Al{sub 5}Au{sub 5} nanocluster possesses a hollow shaped minimum energy structure with segregated Al and Au layered domains, the former representing the electrophilic domain and the latter the nucleophilic domain. In particular, it has been shown that alkali metal cations attach in the nucleophilic domain and hop from one Au site to the next one in the picoseconds time scale, while anions are bound tightly to the Al atoms of the electrophilic domain. Simulating annealing studies are very suggestive of the proneness of the nanocluster towards coalescence into large cluster units, when the cluster is left unprotected by appropriate ligands. Further passivation studies with NaF salt suggest, nonetheless, the possibility of the isolation of the Al{sub 5}Au{sub 5} cluster in molten salts or ionic liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design principles for high–pressure force fields: Aqueous TMAO solutions from ambient to kilobar pressures

Accurate force fields are one of the major pillars on which successful molecular dynamics simulations of complex biomolecular processes rest. They have been optimized for ambient conditions, whereas high-pressure simulations become increasingly important in pressure perturbation studies, using pressure as an independent thermodynamic variable. Here, we explore the design of non-polarizable force fields tailored to work well in the realm of kilobar pressures – while avoiding complete reparameterization. Our key is to first compute the pressure-induced electronic and structural response of a solute by combining an integral equation approach to include pressure effects on solvent structure with a quantum-chemical treatment of the solute within the embedded cluster reference interaction site model (EC-RISM) framework. Next, the solute’s response to compression is taken into account by introducing pressure-dependence into selected parameters of a well-established force field. In our proof-of-principle study, the full machinery is applied to N,N,N-trimethylamine-N-oxide (TMAO) in water being a potent osmolyte that counteracts pressure denaturation. EC-RISM theory is shown to describe well the charge redistribution upon compression of TMAO(aq) to 10 kbar, which is then embodied in force field molecular dynamics by pressure-dependent partial charges. The performance of the high pressure force field is assessed by comparing to experimental and ab initio molecular dynamics data. Beyond its broad usefulness for designing non-polarizable force fields for extreme thermodynamic conditions, a good description of the pressure-response of solutions is highly recommended when constructing and validating polarizable force fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A new force field including charge directionality for TMAO in aqueous solution

We propose a new force field for trimethylamine N-oxide (TMAO), which is designed to reproduce the long-lived and highly directional hydrogen bond between the TMAO oxygen (O{sub TMAO}) atom and surrounding water molecules. Based on the data obtained by ab initio molecular dynamics simulations, we introduce three dummy sites around O{sub TMAO} to mimic the O{sub TMAO} lone pairs and we migrate the negative charge on the O{sub TMAO} to the dummy sites. The force field model developed here improves both structural and dynamical properties of aqueous TMAO solutions. Moreover, it reproduces the experimentally observed dependence of viscosity upon increasing TMAO concentration quantitatively. The simple procedure of the force field construction makes it easy to implement in molecular dynamics simulation packages and makes it compatible with the existing biomolecular force fields. This paves the path for further investigation of protein-TMAO interaction in aqueous solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Permutation invariant polynomial neural network based diabatic ansatz for the ( E + A ) × ( e + a ) Jahn–Teller and Pseudo-Jahn–Teller systems

In this work, the permutation invariant polynomial neural network (PIP-NN) approach is employed to construct a quasi-diabatic Hamiltonian for system with non-Abelian symmetries. It provides a flexible and compact NN-based diabatic ansatz from the related approach of Williams, Eisfeld, and co-workers. The example of $H^3_+$ is studied, which is an (E + A) × (e + a) Jahn–Teller and Pseudo-Jahn–Teller system. The PIP-NN diabatic ansatz is based on the symmetric polynomial expansion of Viel and Eisfeld, the coefficients of which are expressed with neural network functions that take permutation-invariant polynomials as input. This PIP-NN-based diabatic ansatz not only preserves the correct symmetry but also provides functional flexibility to accurately reproduce ab initio electronic structure data, thus resulting in excellent fits. Further, the adiabatic energies, energy gradients, and derivative couplings are well reproduced. A good description of the local topology of the conical intersection seam is also achieved. Therefore, this diabatic ansatz completes the PIP-NN based representation of DPEM with correct symmetries and will enable us to diabatize even more complicated systems with complex symmetries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elasticity of two-dimensional ferroelectrics across their paraelectric phase transformation

The mechanical behavior of two-dimensional (2D) materials across 2D phase changes is unknown, and the finite temperature (T) elasticity of paradigmatic SnSe monolayers—ferroelectric 2D materials turning paraelectric as their unit cell turns from a rectangle into a square—is described here in a progressive manner. To begin with, their zero–T total energy landscape gives way to (Boltzmann-like) averages from which the elastic behavior is determined. Furthermore, these estimates are complemented with results from the strain-fluctuation method, which employs the energy landscape or ab initio molecular dynamics data. All approaches capture the coalescence of elastic moduli < C 11 (T) > = < C 22 (T) > due to the structural transformation. The broad evolution and sudden changes of elastic parameters < C 11 (T) >, < C 22 (T) >, and < C 12 (T) > of these atomically thin phase-change membranes establishes a heretofore overlooked connection among 2D materials and soft matter.

2-dimensional systems↗

BC8 Silicon (Si-III) is a Narrow-Gap Semiconductor

Large-volume, phase-pure synthesis of BC8 silicon ($Ia\bar{3}$, cI16) has allowed bulk measurements of optical, electronic, and thermal properties. Unlike previous reports that conclude BC8-Si is semimetallic, we demonstrate that this phase is a direct band gap semiconductor with a very small energy gap and moderate carrier concentration and mobility at room temperature, based on far- and midinfrared optical spectroscopy, temperature-dependent electrical conductivity, Seebeck and heat capacity measurements. Samples exhibit a plasma wavelength near 11 μm, indicating potential for infrared plasmonic applications. Thermal conductivity is reduced by 1–2 orders of magnitude depending on temperature as compared with the diamond cubic (DC-Si) phase. The electronic structure and dielectric properties can be reproduced by first-principles calculations with hybrid functionals after adjusting the level of exact Hartree–Fock (HF) exchange mixing. These results clarify existing limited and controversial experimental data sets and ab initio calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Design of CO 2 Selective Type 3 Porous Liquids Through Porous Host Morphology

Direct air capture (DAC) of CO 2 is a negative emission technology under development to limit the impacts of climate change. The dilute concentration of CO 2 in the atmosphere (~400 ppm) requires new materials for carbon capture with increased CO 2 selectivity that is not met with current carbon capture materials. Porous liquids (PLs) are an emerging candidate for carbon capture and consists of a combination of solvents and porous hosts that creates a liquid with permanent porosity. The fundamental mechanisms of carbon capture in a PL are relatively unknown. To uncover these mechanisms, PLs were synthesized consisting of three different zeolitic-imidazolate framework (ZIF-8, ZIF-67, or ZIF-69) porous host in a water/glycol/2-methylimidazole solvent. The most stable composition was based on ZIF-8 and exhibited carbon capture following exposure to CO 2 . Density functional theory identified a three-step carbon capture mechanism based on (i) reaction of OH- with ethylene glycol in the solution followed by (ii) formation of 2-hydroxyethyl carbonate, which (iii) further react with OH- to form a carbonate species. This mechanism was validated with experimental nuclear magnetic resonance spectroscopy (NMR) to identify the dissolved carbonate phases and the decrease in the pH during CO 2 exposure. Deuterated samples of the ZIF-8 PLs were synthesized and analyzed via neutron diffraction at the Spallation Neutron Sources at Oak Ridge National Laboratory. Results identified differences in diffraction for PLs pre- and post-CO 2 exposure that will be combined with ab initio molecular dynamics data of the same PL composition to identify how the presence of a solvent-porous host interfaces results in carbon capture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Reduced-order NLTE Kinetic Model for Radiating Plasmas of Outer Envelopes of Stellar Atmospheres

The present work proposes a self-consistent reduced-order NLTE kinetic model for radiating plasmas found in the outer layers of stellar atmospheres. A detailed collisional-radiative kinetic mechanism is constructed by leveraging the most up-to-date set of ab initio and experimental data available in the literature. This constitutes the starting point for the derivation of a reduced-order model, obtained by lumping the bound energy states into groups. In order to determine the needed thermo-physical group properties, uniform and Maxwell–Boltzmann energy distributions are used to reconstruct the energy population of each group. Finally, the reduced set of governing equations for the material gas and the radiation field is obtained based on the moment method. Applications consider the steady flow across a shock wave in partially ionized hydrogen. The results clearly demonstrate that adopting a Maxwell–Boltzmann grouping allows, on the one hand, for a substantial reduction of the number of unknowns and, on the other, to maintain accuracy for both gas and radiation quantities. Also, it is observed that, when neglecting line radiation, the use of two groups already leads to a very accurate resolution of the photo-ionization precursor, internal relaxation, and radiative cooling regions. The inclusion of line radiation requires adopting just one additional group to account for optically thin losses in the α , β , and γ lines of the Balmer and Paschen series. This trend has been observed for a wide range of shock wave velocities.

79 ASTRONOMY AND ASTROPHYSICS↗

Ab Initio Computation of Charge Densities for Sn and Xe Isotopes

We present the first ab initio calculations for open-shell nuclei past the tin isotopic line, focusing on Xe isotopes as well as doubly magic Sn isotopes. We show that, even for moderately hard interactions, it is possible to obtain meaningful predictions and that the NNLO sat chiral interaction predicts radii and charge density distributions close to the experiment. We then make a new prediction for 100 Sn. This paves the way for ab initio studies of exotic charge density distributions at the limit of the present ab initio mass domain, where experimental data is becoming available. The present study closes the gap between the largest isotopes reachable by ab initio methods and the smallest exotic nuclei accessible to electron scattering experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theoretical level energies and transition data for ion W{sup 28+}

The ab initio quasirelativistic approach was used to derive transition data for Pd-like tungsten W{sup 28+} ion. This approach was developed specifically for the calculation of spectral parameters for highly charged ions. The relativistic effects were taken into account in the Breit–Pauli approximation. The configuration interaction method was applied to include electron correlation effects. The ground configuration 4p{sup 6}4d{sup 10} and the lowest excited 4p{sup 6}4d{sup 9}4f were considered. The calculations for W{sup 28+} were performed with different sets of interacting configurations in order to evaluate reliability and accuracy of produced results. Level energies, radiative lifetimes τ, Landé g-factors were determined for the excited configuration 4p{sup 6}4d{sup 9}4f of W{sup 28+} levels. The radiative transition wavelengths λ and emission transition probabilities A for the electric and magnetic transitions of various multipole order (up to E5 and M4) among the levels of these two configurations were produced. The uncertainties of computed spectroscopic parameters were evaluated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ti-6Al-4V to over 1.2 TPa: Shock Hugoniot experiments, ab initio calculations, and a broad-range multiphase equation of state

Titanium alloys are used in a large array of applications. In this work we focus our attention on the most used alloy, Ti-6Al-4V (Ti64), which has excellent mechanical and biocompatibility properties with applications in aerospace, defense, biomedical, and other fields. Here we present high-fidelity experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for Ti64 to more than threefold compression, up to over 1.2 TPa. We use the data to validate our ab initio molecular dynamics simulations and to develop a highly reliable, multiphase equation of state (EOS) for Ti64, spanning a broad range of temperature and pressures. The first-principles simulations show very good agreement with Z data and with previous three-stage gas gun data from Sandia’s STAR facility. The resulting principal Hugoniot and the broad-range EOS and phase diagram up to 10 TPa and 10 5 K are suitable for use in shock experiments and in hydrodynamic simulations. The high-precision experimental results and high-fidelity simulations demonstrate that the Hugoniot of the Ti64 alloy is stiffer than that of pure Ti and reveal that Ti64 melts on the Hugoniot at a significantly lower pressure and temperature than previously modeled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the accuracy of the MB-pol many-body potential for water: Interaction energies, vibrational frequencies, and classical thermodynamic and dynamical properties from clusters to liquid water and ice

The MB-pol many-body potential has recently emerged as an accurate molecular model for water simulations from the gas to the condensed phase. In this study, the accuracy of MB-pol is systematically assessed across the three phases of water through extensive comparisons with experimental data and high-level ab initio calculations. Individual many-body contributions to the interaction energies as well as vibrational spectra of water clusters calculated with MB-pol are in excellent agreement with reference data obtained at the coupled cluster level. Several structural, thermodynamic, and dynamical properties of the liquid phase at atmospheric pressure are investigated through classical molecular dynamics simulations as a function of temperature. The structural properties of the liquid phase are in nearly quantitative agreement with X-ray diffraction data available over the temperature range from 268 to 368 K. The analysis of other thermodynamic and dynamical quantities emphasizes the importance of explicitly including nuclear quantum effects in the simulations, especially at low temperature, for a physically correct description of the properties of liquid water. Furthermore, both densities and lattice energies of several ice phases are also correctly reproduced by MB-pol. Following a recent study of DFT models for water, a score is assigned to each computed property, which demonstrates the high and, in many respects, unprecedented accuracy of MB-pol in representing all three phases of water. Published by AIP Publishing.

Reddy, Sandeep K. (ORCID:0000000214586853)↗

The comparison of semiempirical and ab initio molecular modeling methods in activity and property evaluation of selected antimicrobial sulfonamides

A group of sulfonamides exhibiting antimicrobial activity was analyzed with the application of quantitative structure–activity/property relationships method. The purpose of this study was to show the common and differentiating characteristics of the analyzed chemical structures alike physicochemically as well as pharmacologically based on the quantum-chemical calculations both in vacuo and in the aquatic environment together with their microbiological activity and chromatographic retention data. The semiempirical and ab initio level of in silico molecular modeling was performed for calculations of statistically significant molecular descriptors to compare obtained results. The relationship between the structure and biological activity and physicochemical parameters data was able to class and describe analyzed molecules and the applied chemometric approaches (principal component analysis, factor analysis, and multiple regression analysis) revealed the influential features of the tested structures responsible for the antimicrobial activity.

60 APPLIED LIFE SCIENCES↗