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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Extended Lagrangian Born-Oppenheimer molecular dynamics simulations of the shock-induced chemistry of phenylacetylene

The initial chemical events that occur during the shock compression of liquid phenylacetylene have been investigated using self-consistent tight binding molecular dynamics simulations. The extended Lagrangian Born-Oppenheimer molecular dynamics formalism enabled us to compute microcanonical trajectories with precise conservation of the total energy. Our simulations revealed that the first density-increasing step under shock compression arises from the polymerization of phenylacetylene molecules at the acetylene moiety. The application of electronic structure-based molecular dynamics with long-term conservation of the total energy enabled us to identify electronic signatures of reactivity via monitoring changes in the HOMO-LUMO gap, and to capture directly adiabatic shock heating, transient non-equilibrium states, and changes in temperature arising from exothermic chemistry in classical molecular dynamics trajectories.

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

Nonequilibrium solvent effects in Born-Oppenheimer molecular dynamics for ground and excited electronic states

The effects of solvent on molecular processes such as excited state relaxation and photochemical reaction often occurs in a nonequilibrium regime. Dynamic processes such as these can be simulated using excited state molecular dynamics. In this work, we describe methods of simulating nonequilibrium solvent effects in excited state molecular dynamics using linear-response time-dependent density functional theory and apparent surface charge methods. These developments include a propagation method for solvent degrees of freedom and analytical energy gradients for the calculation of forces. Molecular dynamics of acetaldehyde in water or acetonitrile are demonstrated where the solute-solvent system is out of equilibrium due to photoexcitation and emission.

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

Ab initio molecular dynamics with nuclear quantum effects at classical cost: Ring polymer contraction for density functional theory

Path integral molecular dynamics simulations, combined with an ab initio evaluation of interactions using electronic structure theory, incorporate the quantum mechanical nature of both the electrons and nuclei, which are essential to accurately describe systems containing light nuclei. However, path integral simulations have traditionally required a computational cost around two orders of magnitude greater than treating the nuclei classically, making them prohibitively costly for most applications. Here we show that the cost of path integral simulations can be dramatically reduced by extending our ring polymer contraction approach to ab initio molecular dynamics simulations. By using density functional tight binding as a reference system, we show that our ring polymer contraction scheme gives rapid and systematic convergence to the full path integral density functional theory result. We demonstrate the efficiency of this approach in ab initio simulations of liquid water and the reactive protonated and deprotonated water dimer systems. We find that the vast majority of the nuclear quantum effects are accurately captured using contraction to just the ring polymer centroid, which requires the same number of density functional theory calculations as a classical simulation. Combined with a multiple time step scheme using the same reference system, which allows the time step to be increased, this approach is as fast as a typical classical ab initio molecular dynamics simulation and 35× faster than a full path integral calculation, while still exactly including the quantum sampling of nuclei. This development thus offers a route to routinely include nuclear quantum effects in ab initio molecular dynamics simulations at negligible computational cost.

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

A molecular dynamics study of the effect of thermal boundary conductance on thermal transport of ideal crystal of n-alkanes with different number of carbon atoms

Phase change materials such as n-alkanes that exhibit desirable characteristics such as high latent heat, chemical stability, and negligible supercooling are widely used in thermal energy storage applications. However, n-alkanes have the drawback of low thermal conductivity values. The low thermal conductivity of n-alkanes is linked to formation of randomly oriented nano-domains of molecules in their solid structure that is responsible for excessive phonon scattering at the grain boundaries. Thus, understanding the thermal boundary conductance at the grain boundaries can be crucial for improving the effectiveness of thermal storage systems. The concept of the ideal crystal is proposed in this paper, which describes a simplified model such that all the nano-domains of long-chain n-alkanes are artificially aligned perfectly in one direction. In order to study thermal transport of the ideal crystal of long-chain n-alkanes, four (4) systems (C{sub 20}H{sub 42}, C{sub 24}H{sub 50}, C{sub 26}H{sub 54}, and C{sub 30}H{sub 62}) are investigated by the molecular dynamics simulations. Thermal boundary conductance between the layers of ideal crystals is determined using both non-equilibrium molecular dynamics (NEMD) and equilibrium molecular dynamics (EMD) simulations. Both NEMD and EMD simulations exhibit no significant change in thermal conductance with the molecular length. However, the values obtained from the EMD simulations are less than the values from NEMD simulations with the ratio being nearly three (3) in most cases. This difference is due to the nature of EMD simulations where all the phonons are assumed to be in equilibrium at the interface. Thermal conductivity of the n-alkanes in three structures including liquid, solid, and ideal crystal is investigated utilizing NEMD simulations. Our results exhibit a very slight rise in thermal conductivity values as the number of carbon atoms of the chain increases. The key understanding is that thermal transport can be significantly altered by how the molecules and the nano-domains are oriented in the structure rather than by the length of the n-alkane molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

The molecular gradient using the divide-expand-consolidate resolution of the identity second-order Møller-Plesset perturbation theory: The DEC-RI-MP2 gradient

We report an implementation of the molecular gradient using the divide-expand-consolidate resolution of the identity second-order Møller-Plesset perturbation theory (DEC-RI-MP2). The new DEC-RI-MP2 gradient method combines the precision control as well as the linear-scaling and massively parallel features of the DEC scheme with efficient evaluations of the gradient contributions using the RI approximation. We further demonstrate that the DEC-RI-MP2 gradient method is capable of calculating molecular gradients for very large molecular systems. A test set of supramolecular complexes containing up to 158 atoms and 1960 contracted basis functions has been employed to demonstrate the general applicability of the DEC-RI-MP2 method and to analyze the errors of the DEC approximation. Moreover, the test set contains molecules of complicated electronic structures and is thus deliberately chosen to stress test the DEC-RI-MP2 gradient implementation. Additionally, as a showcase example the full molecular gradient for insulin (787 atoms and 7604 contracted basis functions) has been evaluated.

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

Experimental distinction of the molecularly induced Balmer emission contribution and its application for inferring molecular divertor density with 2D filtered camera measurements during detachment in JET L-mode plasmas

A previously presented model for generating 2D estimates of the divertor plasma conditions at JET from deuterium Balmer line intensity ratios, obtained from tomographic reconstructions of divertor camera images, was amended to consider also the Balmer emission arising from molecular processes. Utilizing the AMJUEL and H2VIBR atomic and molecular databases of EIRENE enabled also inference of the molecular divertor density from the distinguished molecularly induced emission. Analysis of a JET L-mode density scan suggests the molecularly induced emission accounting for up to 60%–70% and 10%–20% of the Balmer Dα and Dγ intensities, respectively, at the onset of detachment, while electron-ion recombination becomes increasingly dominant with deepening detachment. Similar observations were made by post-processing EDGE2D-EIRENE simulations, which indicated significant roles of molecular D$^{+}_{2}$ ions and vibrational excitation of the D 2 molecules as precursors for the molecularly induced emission. The experimentally inferred molecular density at the outer strike point was found to increase monotonously with decreasing strike point temperature, reaching approximately 30%–50% of the local electron density at n mol, osp = 1–2 X 10 20 m -3 at T e,osp ≈ 0.7 eV. A further steep increase by a factor of 3–5 was observed with decrease of T e,osp to 0.5 eV. The observations are in qualitative and reasonable quantitative agreement with EDGE2D-EIRENE predictions of n mol, osp within the uncertainties of the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Toward Molecular Catalysts by Computer

Rational design of molecular catalysts requires a systematic approach to designing ligands with specific functionality and precisely tailored electronic and steric properties. It then becomes possible to devise computer protocols to predict accurately the required properties and ultimately to design catalysts by computer. In this account we first review how thermodynamic properties such as oxidation-reduction potentials (E0), acidities (pKa), and hydride donor abilities (ΔGH-) form the basis for a systematic design of molecular catalysts for reactions that are critical for a secure energy future (hydrogen evolution and oxidation, oxygen and nitrogen reduction, and carbon dioxide reduction). We highlight how density functional theory allows us to determine and predict these properties within “chemical” accuracy (~ 0.06 eV for redox potentials, ~ 1 pKa unit for pKa values, and ~ 1.5 kcal/mol for hydricities). These quantities determine free energy maps and profiles associated with catalytic cycles, i.e. the relative energies of intermediates, and help us distinguish between desirable and high-energy pathways and mechanisms. Good catalysts have flat profiles that avoid high activation barriers due to low and high energy intermediates. We illustrate how the criterion of a flat energy profile lends itself to the prediction of design points by computer for optimum catalysts. This research was carried out in the Center for Molecular Electro-catalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences. Pacific Northwest National Laboratory (PNNL) is operated for the DOE by Battelle.

acidity↗
DOE OSTI2019

Efficiency of low-molecular-weight and high-molecular-weight alkoxyamines in the synthesis of polystyrene

The influence of the structure of high-molecular-weight alkoxyamines (obtained in situ from C-phenyl-N-tert-butylnitrone and 2-(benzylideneamino)-2-methyl-1-phenylpropanol-1-N-oxide) and their low-molecular-weight analogs on styrene polymerization in the temperature range 70—110 ° was analyzed. The controlled synthesis of polystyrene in the presence of low-molecular-weight alkoxyamines was revealed to occur most efficiently at 110 °. The use of high-molecular-weight alkoxyamines makes it possible to perform the controlled process under lower-temperature conditions (70—90 °), when the constant polymerization rate and linear increase in the molecular weight with conversion are observed, and the synthesis of polymers with fairly low polydispersity (M{sub w}/M{sub n} ~ 1.26) is also possible. The quantum chemical calculations showed that the energy of the C-ON bond in alkoxyamines is determined by the length of the chains at the nitroxide fragment rather than the structure of the substituents. It is found that the equilibrium dissociation constants of high-molecular-weight alkoxyamines are caused by both the structure of the nitroxide fragment and molecular weights.

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

Origin of molecular conformational stability: Perspectives from molecular orbital interactions and density functional reactivity theory

To have a quantitative understanding about the origin of conformation stability for molecular systems is still an unaccomplished task. Frontier orbital interactions from molecular orbital theory and energy partition schemes from density functional reactivity theory are the two approaches available in the literature that can be used for this purpose. In this work, we compare the performance of these approaches for a total of 48 simple molecules. We also conduct studies to flexibly bend bond angles for water, carbon dioxide, borane, and ammonia molecules to obtain energy profiles for these systems over a wide range of conformations. We find that results from molecular orbital interactions using frontier occupied orbitals such as the highest occupied molecular orbital and its neighbors are only qualitatively, at most semi-qualitatively, trustworthy. To obtain quantitative insights into relative stability of different conformations, the energy partition approach from density functional reactivity theory is much more reliable. We also find that the electrostatic interaction is the dominant descriptor for conformational stability, and steric and quantum effects are smaller in contribution but their contributions are indispensable. Stable molecular conformations prefer to have a strong electrostatic interaction, small molecular size, and large exchange-correlation effect. This work should shed new light towards establishing a general theoretical framework for molecular stability.

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

Recent Developments in the General Atomic and Molecular Electronic Structure System

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods like the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite like GAMESS is maintained and developed is considered. Future developments are briefly summarized.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Surface Immobilization of Molecular Electrocatalysts for Energy Conversion

Electrocatalysts are critically important for a secure energy future, as they facilitate the conversion between electrical energy and chemical energy. Molecular catalysts offer precise control of their structure, and the ability to modify the substituents to understand structure-reactivity relationships that are more difficult to achieve with heterogeneous catalysts. Molecular electrocatalysts can be immobilized on surfaces by covalent bonds or through non-covalent interactions. Advantages of surface immobilization include the need for less catalyst, avoidance of bimolecular decomposition pathways, and easier determination of catalyst lifetime. Copper-catalyzed click reactions are often used to form covalent bonds to surfaces, and pi-pi stacking of pyrene substituents appended to the ligand of a molecular complex is a frequently used method to achieve non-covalent surface immobilization. This mini-review highlights surface confinement of molecular electrocatalysts for reduction of O2, oxidation of H2O, production of H2, and reduction of CO2.

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

Molecular coupling competing with defects within insulator of the magnetic tunnel junction-based molecular spintronics devices

Nearly 70 years old dream of incorporating molecule as the device element is still challenged by competing defects in almost every experimentally tested molecular device approach. This paper focuses on the magnetic tunnel junction (MTJ) based molecular spintronics device (MTJMSD) method. An MTJMSD utilizes a tunnel barrier to ensure a robust and mass-producible physical gap between two ferromagnetic electrodes. MTJMSD approach may benefit from MTJ's industrial practices; however, the MTJMSD approach still needs to overcome additional challenges arising from the inclusion of magnetic molecules in conjunction with competing defects. Molecular device channels are covalently bonded between two ferromagnets across the insulating barrier. An insulating barrier may possess a variety of potential defects arising during the fabrication or operational phase. This paper describes an experimental and theoretical study of molecular coupling between ferromagnets in the presence of the competing coupling via an insulating tunnel barrier. We discuss the experimental observations of hillocks and pinhole-type defects producing inter-layer coupling that compete with molecular device elements. We performed theoretical simulations to encompass a wide range of competition between molecules and defects. Monte Carlo Simulation (MCS) was used for investigating the defect-induced inter-layer coupling on MTJMSD. Our research may help understand and design molecular spintronics devices utilizing various insulating spacers such as aluminum oxide (AlOx) and magnesium oxide (MgO) on a wide range of metal electrodes. This paper intends to provide practical insights for researchers intending to investigate the molecular device properties via the MTJMSD approach and do not have a background in magnetic tunnel junction fabrication.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Molecular View of Protein Crystal Growth: Molecular Interactions, Surface Reconstruction and Growth Mechanism

Studies of the growth and molecular packing of tetragonal lysozyme crystals suggest that there is an underlying molecular growth mechanism, in addition to the classical one involving screw dislocation/2D) nucleation growth. These crystals are constructed by strongly bonded molecular chains forming helices about the 43 axes. The helices are connected to each other by weaker bonds. Crystal growth proceeds by the formation of these 4(sub 3) helices, which would explain some unexpected observations by earlier investigators, such as bimolecular growth steps on the (110) face. Another consequence of these molecular considerations is that only one of two possible packing arrangements could occur on the crystal faces and that their growth unit was at least a tetramer corresponding to the 4(sub 3) helix. Two new high resolution atomic force microscopy (AFM) techniques were developed to directly confirm these predictions on tetragonal lysozyme crystals. Most earlier investigations of protein crystal growth with AFM were in the low resolution mode which is adequate to investigate the classical growth mechanisms, but cannot resolve molecular features and mechanisms. Employing the first of the newly developed techniques, high resolution AFM images of the (110) face were compared with the theoretically constructed images for the two possible packing arrangements on this face. The prediction that the molecular packing arrangement of these faces corresponded to that for complete 4(sub 3) helices was confirmed in this manner. This investigation also showed the occurrence of surface reconstruction on protein crystals. The molecules on the surface of the (110) face were found to pack closer along the 4(sub 3) axes than those in the interior. The second new AFM technique was used to follow the growth process by measuring the dimensions of individual growth units on the (110) face. Linescans across a growth step, performed near the saturation limit of the crystals, allowed the growth unit dimensions to be measured as they were being incorporated into the crystal. These measurements showed that growth on the (110) face proceeded by the formation of new 4(sub 3) helices from the addition of many growth units of at least tetramer size. This is the first time that a time-resolved AFM technique has been used to observe individual molecular events during crystal growth.

Nadarajah, Arunan↗
DOE OSTI2025

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

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

Atomic and molecular adsorption on Au(111)

Periodic self-consistent density functional theory (DFT-GGA) calculations were used to study the adsorption of several atomic species, molecular species and molecular fragments on the Au(111) surface with a coverage of 1/4 monolayer (ML). Binding geometries, binding energies, and diffusion barriers were calculated for 27 species. Furthermore, we calculated the surface deformation energy associated with the binding events. The binding strength for all the analyzed species can be ordered as follows: NH 3 < NO < CO < CH 3 < HCO < NH 2 < COOH < OH < HCOO < CNH 2 < H < N < NH < NOH < COH < Cl,< HCO 3 < CH 2 < CN b HNO < O < F < S < C < CH. Although the atomic species preferred to bind at the three-fold fcc site, no tendency was observed in site preference for the molecular species and fragments. The intramolecular and adsorbate-surface vibrational frequencies were calculated for all the adsorbates on their most energetically stable adsorption site. Most of the theoretical binding energies and frequencies agreed with experimental values reported in the literature. In general, the values obtained with the PW91 functional are more accurate than RPBE in reproducing these experimental binding energies. The energies of the adsorbed species were used to calculate the thermochemical potential energy surfaces for decomposition of CO, NO, N 2 , NH 3 and CH 4 , oxidation of CO, and hydrogenation of CO, CO 2 and NO, giving insight into the thermochemistry of these reactions on gold nanoparticles. These potential energy surfaces demonstrated that: the decomposition of species is not energetically favorable on Au(111); the desorption of NH 3 , NO and CO are more favorable than their decomposition; the oxidation of CO and hydrogenation of CO and NO on Au(111) to form HCO and HNO, respectively, are also thermodynamically favorable.

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

Combined 3D-QSAR, molecular docking and molecular dynamics study on thyroid hormone activity of hydroxylated polybrominated diphenyl ethers to thyroid receptors β

Several recent reports suggested that hydroxylated polybrominated diphenyl ethers (HO-PBDEs) may disturb thyroid hormone homeostasis. To illuminate the structural features for thyroid hormone activity of HO-PBDEs and the binding mode between HO-PBDEs and thyroid hormone receptor (TR), the hormone activity of a series of HO-PBDEs to thyroid receptors β was studied based on the combination of 3D-QSAR, molecular docking, and molecular dynamics (MD) methods. The ligand- and receptor-based 3D-QSAR models were obtained using Comparative Molecular Similarity Index Analysis (CoMSIA) method. The optimum CoMSIA model with region focusing yielded satisfactory statistical results: leave-one-out cross-validation correlation coefficient (q{sup 2}) was 0.571 and non-cross-validation correlation coefficient (r{sup 2}) was 0.951. Furthermore, the results of internal validation such as bootstrapping, leave-many-out cross-validation, and progressive scrambling as well as external validation indicated the rationality and good predictive ability of the best model. In addition, molecular docking elucidated the conformations of compounds and key amino acid residues at the docking pocket, MD simulation further determined the binding process and validated the rationality of docking results. -- Highlights: ► The thyroid hormone activities of HO-PBDEs were studied by 3D-QSAR. ► The binding modes between HO-PBDEs and TRβ were explored. ► 3D-QSAR, molecular docking, and molecular dynamics (MD) methods were performed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Equation of state of dense plasmas: Orbital-free molecular dynamics as the limit of quantum molecular dynamics for high-Z elements

The applicability of quantum molecular dynamics to the calculation of the equation of state of a dense plasma is limited at high temperature by computational cost. Orbital-free molecular dynamics, based on a semiclassical approximation and possibly on a gradient correction, is a simulation method available at high temperature. For a high-Z element such as lutetium, we examine how orbital-free molecular dynamics applied to the equation of state of a dense plasma can be regarded as the limit of quantum molecular dynamics at high temperature. For the normal mass density and twice the normal mass density, we show that the pressures calculated with the quantum approach converge monotonically towards those calculated with the orbital-free approach; we observe a faster convergence when the orbital-free approach includes the gradient correction. We propose a method to obtain an equation of state reproducing quantum molecular dynamics results up to high temperatures where this approach cannot be directly implemented. With the results already obtained for low-Z plasmas, the present study opens the way for reproducing the quantum molecular dynamics pressure for all elements up to high temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗