Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular electrostatic potential”

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

Multipole correction of atomic monopole models of molecular charge distribution. I. Peptides

The defects in atomic monopole models of molecular charge distribution have been analyzed for several model-blocked peptides and compared with accurate quantum chemical values. The results indicate that the angular characteristics of the molecular electrostatic potential around functional groups capable of forming hydrogen bonds can be considerably distorted within various models relying upon isotropic atomic charges only. It is shown that these defects can be corrected by augmenting the atomic point charge models by cumulative atomic multipole moments (CAMMs). Alternatively, sets of off-center atomic point charges could be automatically derived from respective multipoles, providing approximately equivalent corrections. For the first time, correlated atomic multipoles have been calculated for N-acetyl, N'-methylamide-blocked derivatives of glycine, alanine, cysteine, threonine, leucine, lysine, and serine using the MP2 method. The role of the correlation effects in the peptide molecular charge distribution are discussed.

NASA Discipline Exobiology↗

Cumulative atomic multipole moments complement any atomic charge model to obtain more accurate electrostatic properties

The quality of several atomic charge models based on different definitions has been analyzed using cumulative atomic multipole moments (CAMM). This formalism can generate higher atomic moments starting from any atomic charges, while preserving the corresponding molecular moments. The atomic charge contribution to the higher molecular moments, as well as to the electrostatic potentials, has been examined for CO and HCN molecules at several different levels of theory. The results clearly show that the electrostatic potential obtained from CAMM expansion is convergent up to R-5 term for all atomic charge models used. This illustrates that higher atomic moments can be used to supplement any atomic charge model to obtain more accurate description of electrostatic properties.

NASA Discipline Number 52-20↗

Research in Computational Astrobiology

We present results from several projects in the new field of computational astrobiology, which is devoted to advancing our understanding of the origin, evolution and distribution of life in the Universe using theoretical and computational tools. We have developed a procedure for calculating long-range effects in molecular dynamics using a plane wave expansion of the electrostatic potential. This method is expected to be highly efficient for simulating biological systems on massively parallel supercomputers. We have perform genomics analysis on a family of actin binding proteins. We have performed quantum mechanical calculations on carbon nanotubes and nucleic acids, which simulations will allow us to investigate possible sources of organic material on the early earth. Finally, we have developed a model of protobiological chemistry using neural networks.

Chaban, Galina↗

Roto-Translational Collision-Induced Absorption of CO2 for the Atmosphere of Venus at Frequencies from 0 to 250 cm(exp -1), at Temperatures from 200 to 800 K

The collision-induced absorption (CIA)of gaseous CO2 is the primary source of far-infrared opacity of the atmosphere of Venus. At the temperatures and densities of the venusian atmosphere, the absorption is due mainly to binary collisions of CO2 molecules. Using a realistic anisotropic intermolecular potential and assuming the absorbing dipole to be due to the electrostatic induction and a quantum overlap, a series of molecular dynamics simulations were performed for the temperature range 200 to 800 K, and the roto-translational (RT) collision-induced absorption spectra at frequencies from 0 to 250 cm(exp -1) were derived. The absorption coefficient in the submillimeter region, used in constituency retrieval studies, decreases more than 10 times in the temperature range 200 to 800 K. On the other hand, the absorption coefficient at 800 K and at the frequency range above 150 cm(exp -1) was found to be almost 10 times higher than at 200 K. Earlier works relied on experimental RT CIA data at a fixed temperature of 300 K. The new, temperature-dependent absorption bands may, when included in the analysis of the atmospheric radiative transfer of the planet, help explain the observed high far-infrared opacity of the lower layers of the atmosphere. To make the results of the simulations readily available for atmospheric abundance and radiative transfer analysis, an analytic model of the roto-translational collision-induced absorption spectral profile, applicable from 200 to 800 K, is being proposed here. The FORTRAN computer code of this newly developed model is available from the authors on request.

Gruszka, Marcin↗

The Influence of Chemical Descriptors on the Reactivity of Potential Hypergolic Fuels With Hydrogen Peroxide

The state-of-the-art for storable hypergolic bipropellant systems is monomethylhydrazine (MMH) and mixed oxides of nitrogen (MON). While MMH and MON provide fast and highly reliable hypergolic ignition with excellent propulsive performance, they suffer from safety concerns impacting their production, cost, and supply to various test and launch facilities. In contrast, identifying and demonstrating low-toxicity fuels hypergolic with rocket grade hydrogen peroxide (H2O2) has been challenging with long ignition delay (IDT) and low propulsive performance affecting their deployment in flight systems. This work describes a preliminary observational study to assess the effects of molecular descriptors suspected to predict the reactivity and hypergolicity of compounds within the same chemical class (i.e., thioamide). The approach centered on evaluating functional groups, electrostatic potential, number of carbons, dipole moment, and adiabatic ionization energy. Thiourea is hypergolic with 89.2 wt.% H2O2 (minimum IDT of 21.9 ms) and was selected as a core fuel to study the reactivity of derived fuels with different functionalities. Additionally, the chemical properties were calculated with the Gaussian16 suite. The results showed that group substitution and generally increased number of carbons on thiourea increases IDT. It was observed that a more positively charged sulfur atom in the thioamide does not increase the reactivity towards H2O2. The dipole moment analysis revealed an overall trend of decreasing IDT due to increasing dipole moment. Finally, a lower adiabatic ionization energy related to increased IDTs for most cases. Additional chemical classes comprising a wider range of functionalities and electrophilic attributes are being investigated in on-going work.

Propellant↗

Effect of the Electrostatic Field due to the Surface Monolayer Metal Atoms on the Dissociation of Homonuclear Diatomic Gases in Gas-Metal Surface Reactions

The present paper examines the role of an electrostatic field generated by the outermost monolayer of metal ions prior to gas adsorption in aiding the dissociation of homonuclear diatomic gas molecules. The interaction of five homonuclear diatomic gases, Cl 2 , F 2 , H 2 , N 2 and O 2 , with several pure metals are examined using Coulomb’s law to calculate the attractive and repulsive forces between the electrons or ions in the gas molecule and the free surface metallic electrons or ions assuming a Bohr model. These calculations demonstrate that the total energy of the electrostatic fields from the metals can exceed the molecular binding energies of Cl 2 , F 2 , and O 2 at some distance from the metallic surface thereby suggesting that these gases interact primarily with the metallic surface in their atomic states after molecular bond dissociation prior to reaction. In contrast, H 2 and N 2 gases do not dissociate prior to reaching the metal surface due to the fact that the effective charge of the gas ions is less than the total electron charge at the outermost electronic shell. The present results correlate linearly with the electrochemical series standard reduction potential as well as with the Pauling electronegativity for several metals.

Gas-metal reactions↗

Point charge representation of multicenter multipole moments in calculation of electrostatic properties

Distributed Point Charge Models (PCM) for CO, (H2O)2, and HS-SH molecules have been computed from analytical expressions using multi-center multipole moments. The point charges (set of charges including both atomic and non-atomic positions) exactly reproduce both molecular and segmental multipole moments, thus constituting an accurate representation of the local anisotropy of electrostatic properties. In contrast to other known point charge models, PCM can be used to calculate not only intermolecular, but also intramolecular interactions. Comparison of these results with more accurate calculations demonstrated that PCM can correctly represent both weak and strong (intramolecular) interactions, thus indicating the merit of extending PCM to obtain improved potentials for molecular mechanics and molecular dynamics computational methods.

Non-NASA Center↗

Excess chemical potential of small solutes across water--membrane and water--hexane interfaces

The excess chemical potentials of five small, structurally related solutes, CH4, CH3F, CH2F2, CHF3, and CF4, across the water-glycerol 1-monooleate bilayer and water-hexane interfaces were calculated at 300, 310, and 340 K using the particle insertion method. The excess chemical potentials of nonpolar molecules (CH4 and CF4) decrease monotonically or nearly monotonically from water to a nonpolar phase. In contrast, for molecules that possess permanent dipole moments (CH3F, CH2F, and CHF3), the excess chemical potentials exhibit an interfacial minimum that arises from superposition of two monotonically and oppositely changing contributions: electrostatic and nonelectrostatic. The nonelectrostatic term, dominated by the reversible work of creating a cavity that accommodates the solute, decreases, whereas the electrostatic term increases across the interface from water to the membrane interior. In water, the dependence of this term on the dipole moment is accurately described by second order perturbation theory. To achieve the same accuracy at the interface, third order terms must also be included. In the interfacial region, the molecular structure of the solvent influences both the excess chemical potential and solute orientations. The excess chemical potential across the interface increases with temperature, but this effect is rather small. Our analysis indicates that a broad range of small, moderately polar molecules should be surface active at the water-membrane and water-oil interfaces. The biological and medical significance of this result, especially in relation to the mechanism of anesthetic action, is discussed.

NASA Center ARC↗

Membrane rejection of nitrogen compounds

Rejection characteristics of nitrogen compounds were examined for reverse osmosis, nanofiltration, and low-pressure reverse osmosis membranes. The rejection of nitrogen compounds is explained by integrating experimental results with calculations using the extended Nernst-Planck model coupled with a steric hindrance model. The molecular weight and chemical structure of nitrogen compounds appear to be less important in determining rejection than electrostatic properties. The rejection is greatest when the Donnan potential exceeds 0.05 V or when the ratio of the solute radius to the pore radius is greater than 0.8. The transport of solute in the pore is dominated by diffusion, although convective transport is significant for organic nitrogen compounds. Electromigration contributes negligibly to the overall solute transport in the membrane. Urea, a small organic compound, has lower rejection than ionic compounds such as ammonium, nitrate, and nitrite, indicating the critical role of electrostatic interaction in rejection. This suggests that better treatment efficiency for organic nitrogen compounds can be obtained after ammonification of urea.

NASA Discipline Life Support Systems↗

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗

Room-temperature Broadening and Pressure-shift Coefficients in the nu(exp 2) Band of CH3D-O2: Measurements and Semi-classical Calculations

We report measured Lorentz O2-broadening and O2-induced pressure-shift coefficients of CH3D in the nu(exp 2) fundamental band. Using a multispectrum fitting technique we have analyzed 11 laboratory absorption spectra recorded at 0.011 cm(exp 1) resolution using the McMath-Pierce Fourier transform spectrometer, Kitt Peak, Arizona. Two absorption cells with path lengths of 10.2 and 25 cm were used to record the spectra. The total sample pressures ranged from 0.98 to 339.85 Torr with CH3D volume mixing ratios of 0.012 in oxygen. We report measurements for O2 pressure-broadening coefficients of 320 nu(exp 2) transitions with quantum numbers as high as J0(sup w) = 17 and K = 14, where K(sup w) = K' is equivalent to K (for a parallel band). The measured O2-broadening coefficients range from 0.0153 to 0.0645 cm(exp -1) atm(exp -1) at 296 K. All the measured pressure-shifts are negative. The reported O2-induced pressure-shift coefficients vary from about -0.0017 to -0.0068 cm(exp -1) atm(exp -1). We have examined the dependence of the measured broadening and shift parameters on the J(sup W), and K quantum numbers and also developed empirical expressions to describe the broadening coefficients in terms of m (m = -J(sup W), J(sup W), and J(sup w) + 1 in the QP-, QQ-, and QR-branch, respectively) and K. On average, the empirical expressions reproduce the measured broadening coefficients to within 4.4%. The O2-broadening and pressure shift coefficients were calculated on the basis of a semiclassical model of interacting linear molecules performed by considering in addition to the electrostatic contributions the atom-atom Lennard-Jones potential. The theoretical results of the broadening coefficients are generally larger than the experimental data. Using for the trajectory model an isotropic Lennard-Jones potential derived from molecular parameters instead of the spherical average of the atom-atom model, a better agreement is obtained with these data, especially for |m| <= 12 values (11.3% for the first calculation and 8.1% for the second calculation). The O2-pressure shifts whose vibrational contribution are either derived from parameters fitted in the QQ-branch of selfinduced shifts of CH3D or those obtained from pressure shifts induced by Xe in the nu(sup 3) band of CH3D are in reasonable agreement with the scattered experimental data (17.0% for the first calculation and 18.7% for the second calculation).

Predoi-Cross, Adriana↗

Electric Potential and Electric Field Imaging with Dynamic Applications: 2017 Research Award Innovation

The technology and methods for remote quantitative imaging of electrostatic potentials and electrostatic fields in and around objects and in free space is presented. Electric field imaging (EFI) technology may be applied to characterize intrinsic or existing electric potentials and electric fields, or an externally generated electrostatic field may be used for illuminating volumes to be inspected with EFI. The baseline sensor technology (e-Sensor) and its construction, optional electric field generation (quasi-static generator), and current e- Sensor enhancements (ephemeral e-Sensor) are discussed. Critical design elements of current linear and real-time two-dimensional (2D) measurement systems are highlighted, and the development of a three dimensional (3D) EFI system is presented. Demonstrations for structural, electronic, human, and memory applications are shown. Recent work demonstrates that phonons may be used to create and annihilate electric dipoles within structures. Phonon induced dipoles are ephemeral and their polarization, strength, and location may be quantitatively characterized by EFI providing a new subsurface Phonon-EFI imaging technology. Initial results from real-time imaging of combustion and ion flow, and their measurement complications, will be discussed. These new EFI capabilities are demonstrated to characterize electric charge distribution creating a new field of study embracing areas of interest including electrostatic discharge (ESD) mitigation, crime scene forensics, design and materials selection for advanced sensors, combustion science, on-orbit space potential, container inspection, remote characterization of electronic circuits and level of activation, dielectric morphology of structures, tether integrity, organic molecular memory, atmospheric science, and medical diagnostic and treatment efficacy applications such as cardiac polarization wave propagation and electromyography imaging.

Generazio, Ed↗

Electric Potential and Electric Field Imaging with Dynamic Applications & Extensions

The technology and methods for remote quantitative imaging of electrostatic potentials and electrostatic fields in and around objects and in free space is presented. Electric field imaging (EFI) technology may be applied to characterize intrinsic or existing electric potentials and electric fields, or an externally generated electrostatic field made be used for volumes to be inspected with EFI. The baseline sensor technology (e-Sensor) and its construction, optional electric field generation (quasi-static generator), and current e- Sensor enhancements (ephemeral e-Sensor) are discussed. Critical design elements of current linear and real-time two-dimensional (2D) measurement systems are highlighted, and the development of a three dimensional (3D) EFI system is presented. Demonstrations for structural, electronic, human, and memory applications are shown. Recent work demonstrates that phonons may be used to create and annihilate electric dipoles within structures. Phonon induced dipoles are ephemeral and their polarization, strength, and location may be quantitatively characterized by EFI providing a new subsurface Phonon-EFI imaging technology. Results from real-time imaging of combustion and ion flow, and their measurement complications, will be discussed. Extensions to environment, Space and subterranean applications will be presented, and initial results for quantitative characterizing material properties are shown. A wearable EFI system has been developed by using fundamental EFI concepts. These new EFI capabilities are demonstrated to characterize electric charge distribution creating a new field of study embracing areas of interest including electrostatic discharge (ESD) mitigation, manufacturing quality control, crime scene forensics, design and materials selection for advanced sensors, combustion science, on-orbit space potential, container inspection, remote characterization of electronic circuits and level of activation, dielectric morphology of structures, tether integrity, organic molecular memory, atmospheric science, weather prediction, earth quake prediction, and medical diagnostic and treatment efficacy applications such as cardiac polarization wave propagation and electromyography imaging.

Generazio, Ed↗

Modeling the effects of hydrogen and dose rate sensitivity in CMOS and bipolar technologies

Low dose rate (LDR) and high dose rate (HDR) experiments on field-oxide-field-effect-transistors (FOXFETs) and gated lateral pnp (GLPNP) bipolar transistors indicate that there is a dose rate enhancement factor (EF) associated with radiation-induced degradation. This EF is also affect by the presence of hydrogen in the oxide. In this work, we developed a one dimensional (1-D) numerical calculations code to investigate the key mechanisms that describe the dose rate sensitivity and the effect of hydrogen on dose rate effects. We used a finite-difference methodology for the numerical calculations that allow for computing solutions for the densities of the mobile species as well as for the electrostatic potential at nodes contained within a mesh superimposed on the solution domain. Results from calculations of damage EF indicate that oxide thickness, distribution of hole traps and hole capture cross-section affect dose rate sensitivity. In addition, calculations show that molecular hydrogen cracking at positively charged defects may be a key reaction relating hydrogen and dose rate response. Comparison to experimental data on bipolar and CMOS devices results in good agreement with the dose rate calculations of interface trap buildup

Barnaby, Hugh↗

Surface potential of the water liquid-vapor interface

An analysis of an extended molecular dynamics calculation of the surface potential (SP) of the water liquid-vapor interface is presented. The SP predicted by the TIP4P model is -(130 + or - 50) mV. This value is of reasonable magnitude but of opposite sign to the expectations based on laboratory experiments. The electrostatic potential shows a nonmonotonic variation with depth into the liquid.

Wilson, Michael A.↗

A molecular dynamics study of chloride binding by the cryptand SC24

The capture of chloride from water by the tetraprotonated form of the spherical macrotricyclic molecule SC24 was studied using molecular dynamics simulation methods. This model ionophore represents a broad class of molecules which remove ions from water. Two binding sites for the chloride were found, one inside and one outside the ligand. These sites are separated by a potential energy barrier of approximately 20 kcal mol-1. The major contribution to this barrier comes from dehydration of the chloride. The large, unfavorable dehydration effect is compensated for by an increase in electrostatic attraction between the oppositely charged chloride and cryptand, and by energetically favorable rearrangements of water structure. Additional assistance in crossing the barrier and completing the dehydration of the ion is provided by the shift of three positively charged hydrogen atoms of the cryptand towards the chloride. This structural rigidity is partially responsible for its selectivity.

NASA Discipline Exobiology↗

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗

Influence of Molecular Simulation Model Accuracy on the Interfacial Properties of an Ionic Liquid: Overview of Recommended Practices

Increasing the energy storage capability of ionic liquid supercapacitors will require better understanding of ion-electrode interactions. We have probed the influence of these interactions on the structure and differential capacitance of of an ionic liquid ([EMIM][BF4]) at an ideal graphite interface as a function of model accuracy. Of note, differential capacitance is determined through newly derived and validated fluctuation formulas. In terms of model accuracy, we test electrostatic techniques, electrode charging techniques, and electrolyte interatomic potentials. For electrostatic summations, we employ high cost, high fidelity techniques as well as less expensive, approximate techniques for summation in slab geometry. For electrode charging, uniform, constant-charge and environmentally responsive, constant-potential conditions are employed. For the ionic liquid, constant charge and atomically polarizable models are employed. We comment on the role of model accuracy on the structure and energetics of the electric double layer as well as on the magnitude and shape of differential capacitance.

simulations↗