Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MOLECULAR”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Generalized adsorption isotherms for molecular and dissociative adsorption of a polar molecular species on two polar surface geometries: Perovskite (100) (Pm-3m) and fluorite (111) (Fm-3m)

Lattice based kinetic Monte Carlo simulations have been used to determine a functional form for the second order adsorption isotherms on two commonly investigated crystal surfaces: the (111) fluorite surface and the (100) perovskite surface which has the same geometric symmetry as the NaCl (100) surface. The functional form is generalized to be applicable to all values of the equilibrium constant by a shift along the pressure axis. Functions have been determined for estimating the pressure at which a desired coverage would be achieved and, conversely, for estimating the coverage at a certain pressure. The generalized form has been calculated by investigating the surface adsorbate coverage across a range of thermodynamic equilibrium constants that span the range 10{sup −26} to 10{sup 13}. The equations have been shown to be general for any value of the adsorption equilibrium constant.

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

Ab Initio Molecular Dynamics Simulations of Low-Energy Recoil Events in ThO2, CeO2, and ZrO2

Ab initio molecular dynamics simulations of low-energy recoil events in ThO2, CeO2, and ZrO2 have been carried out to determine the threshold displacement energies, resulting defect configurations, dynamics of defect generation, and role of charge transfer during the process. The results reveal that, in most cases, these fluorite structure oxides exhibit a similar response to low-energy recoils. A variety of different defect configurations are created, consisting mainly of vacancies and interstitials. Charge transfer occurs during the dynamic displacement process. Local charge redistribution leads to cation and O vacancies being negatively and positively charged, respectively. Likewise, due to charge redistribution, the cation and O interstitials are less positively and negatively charged, respectively, than the ions on lattice sites in perfect MO2.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2014

Phonon-magnon interactions in body centered cubic iron: A combined molecular and spin dynamics study

Combining an atomistic many-body potential with a classical spin Hamiltonian parameterized by first principles calculations, molecular-spin dynamics computer simulations were performed to investigate phonon-magnon interactions in body centered cubic iron. Results obtained for spin-spin and density-density dynamic structure factors show that noticeable softening and damping of magnon modes occur due to the presence of lattice vibrations. Furthermore, as a result of the phonon-magnon coupling, additional longitudinal spin wave excitations are observed, with the same frequencies as the longitudinal phonon modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Chemical insight from density functional modeling of molecular adsorption: Tracking the bonding and diffusion of anthracene derivatives on Cu(111) with molecular orbitals

We present a method of analyzing the results of density functional modeling of molecular adsorption in terms of an analogue of molecular orbitals. This approach permits intuitive chemical insight into the adsorption process. Applied to a set of anthracene derivates (anthracene, 9,10-anthraquinone, 9,10-dithioanthracene, and 9,10-diselenonanthracene), we follow the electronic states of the molecules that are involved in the bonding process and correlate them to both the molecular adsorption geometry and the species’ diffusive behavior. We additionally provide computational code to easily repeat this analysis on any system.

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

The Role of Pendant Amines in the Breaking and Forming of Molecular Hydrogen Catalyzed by Nickel Complexes

We present the results of a comprehensive theoretical investigation of the role of pendant amine ligands in the oxidation of H2 and formation of H2 by [Ni(P2RN2R’)2]2+ electrocatalysts (where P2RN2R’ is the 1,5-R’-3,7-R derivative of 1,5-diaza-3,7-diphosphacyclooctane with R and R’ are aryl or alkyl groups). We focus our analysis on the thermolchemical part of the catalytic cycle, as it is known to be rate determining for both H2 oxidation and production. We find that the presence of pendant amine functional groups greatly facilitates the heterolytic H2 bond cleavage resulting in a protonated amine functional group and a Ni hydride. Only one single well-positioned pendant amine is required to serve this function. The pendant amine can also effectively shuttle protons to the active site making the redistribution of protons and the H2 evolution a very facile processes. An important requirement for the overall catalytic process is the positioning of at least one amine in close proximity to the metal center. Indeed, only protonation of the pendant amines on the metal center side (endo position) leads to catalytically active intermediates, whereas protonation on the opposite side of the metal center (exo position) leads to a variety of isomers, which are detrimental to catalysis. This material is based upon work supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the US Department of Energy, Office of Science, Office of Basic Energy Sciences under FWP56073.

08 HYDROGEN↗
DOE OSTI2017

Nonadiabatic Ab Initio Molecular Dynamics with the Floating Occupation Molecular Orbital-Complete Active Space Configuration Interaction Method [Non-Adiabatic Ab Initio Molecular Dynamics with Floating Occupation Molecular Orbitals CASCI Method]

Here, we show that the floating occupation molecular orbital complete active space configuration interaction (FOMO-CASCI) method is a promising alternative to the widely used complete active space self-consistent field (CASSCF) method in direct nonadiabatic dynamics simulations. We have simulated photodynamics of three archetypal molecules in photodynamics: ethylene, methaniminium cation, and malonaldehyde. We compared the time evolution of electronic populations and reaction mechanisms as revealed by the FOMO-CASCI and CASSCF approaches. Generally, the two approaches provide similar results. Some dynamical differences are observed, but these can be traced back to energetically minor differences in the potential energy surfaces. We suggest that the FOMO-CASCI method represents, due to its efficiency and stability, a promising approach for direct ab initio dynamics in the excited state.

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

Preparation of Core-Shell Coordination Molecular Assemblies via the Enrichment of Structure-Directing "Codes" of Bridging Ligands and Metathesis of Metal Units

A series of molybdenum- and copper-based MOPs were synthesized through coordination-driven process of a bridging ligand (3,3'-PDBAD, L-1) and dimetal paddlewheel clusters. Three conformers of the ligand exist with an ideal bridging angle between the two carboxylate groups of 0 degrees (H-2 zeta-L(1)), 120 degrees (H-2 beta-L-1), and of 90 degrees (H-2 beta-L-1), respectively. At ambient or lower temperature, (HL1)-L-2 and Mo-2(OAc)(4) or Cu-2(OAc)(4) were crystallized into a molecular square with ?-L-1 and Mo-2/Cu-2 units. With proper temperature elevation, not only the molecular square with ?-L-1 but also a lantern-shaped cage with a-L-1 formed simultaneously. Similar to how Watson-Crick pairs stabilize the helical structure of duplex DNA, the core-shell molecular assembly possesses favorable H-bonding interaction sites. This is dictated by the ligand conformation in the shell, coding for the formation and providing stabilization of the central lantern shaped core, which was not observed without this complementary interaction. On the basis of the crystallographic implications, a heterobimetallic cage was obtained through a postsynthetic metal ion metathesis, showing different reactivity of coordination bonds in the core and shell. As an innovative synthetic strategy, the site-selective metathesis broadens the structural diversity and properties of coordination assemblies.

Park, J↗
DOE OSTI2013

Extending molecular simulation time scales: Parallel in time integrations for high-level quantum chemistry and complex force representations

Parallel in time simulation algorithms are presented and applied to conventional molecular dynamics (MD) and ab initio molecular dynamics (AIMD) models of realistic complexity. Assuming that a forward time integrator, f (e.g., Verlet algorithm), is available to propagate the system from time t{sub i} (trajectory positions and velocities x{sub i} = (r{sub i}, v{sub i})) to time t{sub i+1} (x{sub i+1}) by x{sub i+1} = f{sub i}(x{sub i}), the dynamics problem spanning an interval from t{sub 0}…t{sub M} can be transformed into a root finding problem, F(X) = [x{sub i} − f(x{sub (i−1})]{sub i} {sub =1,M} = 0, for the trajectory variables. The root finding problem is solved using a variety of root finding techniques, including quasi-Newton and preconditioned quasi-Newton schemes that are all unconditionally convergent. The algorithms are parallelized by assigning a processor to each time-step entry in the columns of F(X). The relation of this approach to other recently proposed parallel in time methods is discussed, and the effectiveness of various approaches to solving the root finding problem is tested. We demonstrate that more efficient dynamical models based on simplified interactions or coarsening time-steps provide preconditioners for the root finding problem. However, for MD and AIMD simulations, such preconditioners are not required to obtain reasonable convergence and their cost must be considered in the performance of the algorithm. The parallel in time algorithms developed are tested by applying them to MD and AIMD simulations of size and complexity similar to those encountered in present day applications. These include a 1000 Si atom MD simulation using Stillinger-Weber potentials, and a HCl + 4H{sub 2}O AIMD simulation at the MP2 level. The maximum speedup ((serial execution time)/(parallel execution time) ) obtained by parallelizing the Stillinger-Weber MD simulation was nearly 3.0. For the AIMD MP2 simulations, the algorithms achieved speedups of up to 14.3. The parallel in time algorithms can be implemented in a distributed computing environment using very slow transmission control protocol/Internet protocol networks. Scripts written in Python that make calls to a precompiled quantum chemistry package (NWChem) are demonstrated to provide an actual speedup of 8.2 for a 2.5 ps AIMD simulation of HCl + 4H{sub 2}O at the MP2/6-31G* level. Implemented in this way these algorithms can be used for long time high-level AIMD simulations at a modest cost using machines connected by very slow networks such as WiFi, or in different time zones connected by the Internet. The algorithms can also be used with programs that are already parallel. Using these algorithms, we are able to reduce the cost of a MP2/6-311++G(2d,2p) simulation that had reached its maximum possible speedup in the parallelization of the electronic structure calculation from 32 s/time step to 6.9 s/time step.

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

Systematic evaluation of structure–property relationships in heteroacene – diketopyrrolopyrrole molecular donors for organic solar cells

Improved understanding of fundamental structure–property relationships, particularly the effects of molecular shape and intermolecular packing on film morphology and active layer charge transport characteristics, enables more rational synthesis of new p-type small molecules. Here we investigate a series of small molecules consisting of an acene-based electron-rich core flanked by one or two electron-deficient diketopyrrolopyrrole (DPP) moieties. Through minor changes in the molecule structures, measurable variations in the crystal structure and sizable differences in macroscopic properties are achieved. The molecular symmetry as well as the conformation of the side chains affects the unit cell packing density and strength of the intermolecular electronic coupling in single crystals of all molecules in this series. Here, the addition of a second DPP unit to the benzodithiophene (BDT) core increases molecular planarity leading to decreased reorganization energy, strong cofacial coupling, and moderate hole mobility (2.7 × 10 –4 cm 2 V –1 s –1 ). Increasing the length of the acene core from benzodithiophene to naphthodithiophene (NDT) results in a further reduction in reorganization energy and formation of smaller crystalline domains (~11 nm) when mixed with PCBM. Decreasing the aspect ratio of the core using a “zig-zag” naphthodithiophene (zNDT) isomer results in the highest hole mobility of 1.3 × 10 –3 cm 2 V –1 s –1 due in part to tight lamellar (d = 13.5 Å) and π–π stacking (d = 3.9 Å). The hole mobility is directly correlated with the short-circuit current (11.7 mA cm –2 ) and solar cell efficiency (4.4%) of the highest performing zNDT:PCBM device. For each of these small molecules the calculated π-coupling constant is correlated with the hole mobility as a function of crystal structure and orientation indicating the importance of designing molecules that create extended crystalline networks with maximal π-orbital overlap.

36 MATERIALS SCIENCE↗
DOE OSTI2018

25 Mg NMR and Computational Modeling Studies of the Solvation Structures and Molecular Dynamics in Magnesium Based Liquid Electrolytes

There is increasing evidence that the solvation structure of the active components in a liquid electrolyte solution strongly impacts the performance in electrochemical applications. In this work, the nanoscale solvation structures and dynamics of Mg(BH 4 ) 2 and Mg(TFSI) 2 dissolved in diglyme (DGM) at various concentrations and ratios of Mg(BH 4 ) 2 /Mg(TFSI) 2 were investigated using a combination of natural abundance 25 Mg NMR, quantum chemistry calculations of 25 Mg NMR chemical shifts, classical molecular dynamics (MD) calculations, and electrochemical performance tests. By mixing two competing Mg salts, we were able to reduce the strong covalent interactions between Mg 2+ and BH 4 – anions. A small increase is observed in the coordination number of Mg-TFSI and a significant increase in the interaction of Mg 2+ ions with glymes. Through a combination of NMR, DFT and MD simulations, various stable species around 1 nm in size were detected in the mixed salt solution, which play key roles in the enhanced electrochemical performance of the mixed electrolyte. It is established that for the neat Mg(TFSI) 2 in DGM electrolyte at dilute concentrations the TFSI - is fully dissociated from Mg 2+ . At higher concentrations, Mg 2+ and TFSI - are only partially dissociated as contact ion pairs are formed. In contrast, at 0.01 M Mg(BH 4 ) 2 (saturated concentration) in DGM, the first solvation shell of a Mg 2+ ion contains two BH 4 - anions and one DGM molecule, while the second solvation shell consists of five to six DGM molecules. An exchange mechanism between the solvation structures in the combined electrolyte containing both Mg(BH 4 ) 2 and Mg(TFSI) 2 in DGM was found to result in the observation of a single 25 Mg NMR peak. This exchange is responsible for an increase in uncoordinated anions, as well as improved stability and ionic conductivity as compared to single anion solution. Solvent molecule rearrangement and direct Mg-ion exchange between the basic solvation structures are hypothesized as likely reasons for the exchange. We elucidate that the solvent rearrangement is energetically much more favorable than direct Mg-ion hopping and is thus suggested as the dominant exchange mechanism.

Magnesium battery, electrolytes, Mg(BH4)2, Mg(TFSI↗
DOE OSTI2017

Two-temperature model in molecular dynamics simulations of cascades in Ni-based alloys

In high-energy irradiation events, energy from the fast moving ion is transferred to the system via nuclear and electronic energy loss mechanisms. The nuclear energy loss results in the creation of point defects and clusters, while the energy transferred to the electrons results in the creation of high electronic temperatures, which can affect the damage evolution. We perform molecular dynamics simulations of 30 keV and 50 keV Ni ion cascades in nickel-based alloys without and with the electronic effects taken into account. We compare the results of classical molecular dynamics (MD) simulations, where the electronic effects are ignored, with results from simulations that include the electronic stopping only, as well as simulations where both the electronic stopping and the electron-phonon coupling are incorporated, as described by the two temperature model (2T-MD). Our results indicate that the 2T-MD leads to a smaller amount of damage, more isolated defects and smaller defect clusters.

Zarkadoula, Eva↗
DOE OSTI2014

Lattice thermal conductivity of UO{sub 2} using ab-initio and classical molecular dynamics

We applied the non-equilibrium ab-initio molecular dynamics and predict the lattice thermal conductivity of the pristine uranium dioxide for up to 2000 K. We also use the equilibrium classical molecular dynamics and heat-current autocorrelation decay theory to decompose the lattice thermal conductivity into acoustic and optical components. The predicted optical phonon transport is temperature independent and small, while the acoustic component follows the Slack relation and is in good agreement with the limited single-crystal experimental results. Considering the phonon grain-boundary and pore scatterings, the effective lattice thermal conductivity is reduced, and we show it is in general agreement with the sintered-powder experimental results. The charge and photon thermal conductivities are also addressed, and we find small roles for electron, surface polaron, and photon in the defect-free structures and for temperatures below 1500 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

A hybrid framework of first principles molecular orbital calculations and a three-dimensional integral equation theory for molecular liquids: Multi-center molecular Ornstein–Zernike self-consistent field approach

In this study, we reported the development of a new quantum mechanics/molecular mechanics (QM/MM)-type framework to describe chemical processes in solution by combining standard molecular-orbital calculations with a three-dimensional formalism of integral equation theory for molecular liquids (multi-center molecular Ornstein–Zernike (MC-MOZ) method). The theoretical procedure is very similar to the 3D-reference interaction site model self-consistent field (RISM-SCF) approach. Since the MC-MOZ method is highly parallelized for computation, the present approach has the potential to be one of the most efficient procedures to treat chemical processes in solution. Benchmark tests to check the validity of this approach were performed for two solute (solute water and formaldehyde) systems and a simple S{sub N}2 reaction (Cl{sup −} + CH{sub 3}Cl → ClCH{sub 3} + Cl{sup −}) in aqueous solution. The results for solute molecular properties and solvation structures obtained by the present approach were in reasonable agreement with those obtained by other hybrid frameworks and experiments. In particular, the results of the proposed approach are in excellent agreements with those of 3D-RISM-SCF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗