Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvation models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

17 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Configurations of base-pair complexes in solutions

A theoretical search for the most stable conformations (i.e., stacked or hydrogen bonded) of the base pairs A-U and G-C in water, CCl4, and CHCl3 solutions is presented. The calculations of free energies indicate a significant role of the solvent in determining the conformations of the base-pair complexes. The application of the continuum method yields preferred conformations in good agreement with experiment. Results of the calculations with this method emphasize the importance of both the electrostatic interactions between the two bases in a complex, and the dipolar interaction of the complex with the entire medium. In calculations with the solvation shell method, the last term, i.e., dipolar interaction of the complex with the entire medium, was added. With this modification the prediction of the solvation shell model agrees both with the continuum model and with experiment, i.e., in water the stacked conformation of the bases is preferred.

Egan, J. T.↗
NASA NTRS2004

Theoretical Investigation of OCN(-) Charge Transfer Complexes in Condensed Phase Media: Spectroscopic Properties in Amorphous Ice

Density functional theory (DFT) calculations of cyanate (OCN(-)) charge-transfer complexes were performed to model the "XCN" feature observed in interstellar icy grain mantles. OCN(-) charge-transfer complexes were formed from precursor combinations of HNCO or HOCN with either NH3 or H2O. Three different solvation strategies for realistically modeling the ice matrix environment were explored, including (1) continuum solvation, (2) pure DFT cluster calculations, and (3) an ONIOM DFT/PM3 cluster calculation. The model complexes were evaluated by their ability to reproduce seven spectroscopic measurements associated with XCN: the band origin of the OCN(-) asymmetric stretching mode, shifts in that frequency due to isotopic substitutions of C, N, O, and H, plus two weak features. The continuum solvent field method produced results consistent with some of the experimental data but failed to account for other behavior due to its limited capacity to describe molecular interactions with solvent. DFT cluster calculations successfully reproduced the available spectroscopic measurements very well. In particular, the deuterium shift showed excellent agreement in complexes where OCN(-) was fully solvated. Detailed studies of representative complexes including from two to twelve water molecules allowed the exploration of various possible solvation structures and provided insights into solvation trends. Moreover, complexes arising from cyanic or isocyanic acid in pure water suggested an alternative mechanism for the formation of OCN(-) charge-transfer complexes without the need for a strong base such as NH3 to be present. An extended ONIOM (B3LYP/PM3) cluster calculation was also performed to assess the impact of a more realistic environment on HNCO dissociation in pure water.

Park, Jin-Young↗
NASA NTRS2001

Substituent and Solvent Effects on Excited State Charge Transfer Behavior of Highly Fluorescent Dyes Containing Thiophenylimidazole-Based Aldehydes

The excited state charge transfer for a series of highly fluorescent dyes containing thiophenylimidazole moiety was investigated. These systems follow the Twisted Intramolecular Charge Transfer (TICT) model. Dual fluorescence was observed for each substituted dye. X-ray structures analysis reveals a twisted ground state geometry for the donor substituted aryl on the 4 and 5 position at the imidazole ring. The excited state charge transfer was modeled by a linear solvation energy relationship using Taft's pi and Dimroth's E(sub T)(30) as solvent parameters. There is linear relation between the energy of the fluorescence transition and solvent polarity. The degree of stabilization of the excited state charge transfer was found to be consistent with the intramolecular molecular charge transfer. Excited dipole moment was studied by utilizing the solvatochromic shift method.

Santos, Javier↗
NASA NTRS2018

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗
NASA NTRS2006

Hybrid MD-Nernst Planck Model of Alpha-hemolysin Conductance Properties

Motivated by experiments in which an applied electric field translocates polynucleotides through an alpha-hemolysin protein channel causing ionic current transient blockade, a hybrid simulation model is proposed to predict the conductance properties of the open channel. Time scales corresponding to ion permeation processes are reached using the Poisson-Nemst-Planck (PNP) electro-diffusion model in which both solvent and local ion concentrations are represented as a continuum. The diffusion coefficients of the ions (K(+) and Cl(-)) input in the PNP model are, however, calculated from all-atom molecular dynamics (MD). In the MD simulations, a reduced representation of the channel is used. The channel is solvated in a 1 M KCI solution, and an external electric field is applied. The pore specific diffusion coefficients for both ionic species are reduced 5-7 times in comparison to bulk values. Significant statistical variations (17-45%) of the pore-ions diffusivities are observed. Within the statistics, the ionic diffusivities remain invariable for a range of external applied voltages between 30 and 240mV. In the 2D-PNP calculations, the pore stem is approximated by a smooth cylinder of radius approx. 9A with two constriction blocks where the radius is reduced to approx. 6A. The electrostatic potential includes the contribution from the atomistic charges. The MD-PNP model shows that the atomic charges are responsible for the rectifying behaviour and for the slight anion selectivity of the a-hemolysin pore. Independent of the hierarchy between the anion and cation diffusivities, the anionic contribution to the total ionic current will dominate. The predictions of the MD-PNP model are in good agreement with experimental data and give confidence in the present approach of bridging time scales by combining a microscopic and macroscopic model.

Cozmuta, Ioana↗
NASA NTRS2004

Computational studies of Ras and PI3K

Until recently, experimental techniques in molecular cell biology have been the primary means to investigate biological risk upon space radiation. However, computational modeling provides an alternative theoretical approach, which utilizes various computational tools to simulate proteins, nucleotides, and their interactions. In this study, we are focused on using molecular mechanics (MM) and molecular dynamics (MD) to study the mechanism of protein-protein binding and to estimate the binding free energy between proteins. Ras is a key element in a variety of cell processes, and its activation of phosphoinositide 3-kinase (PI3K) is important for survival of transformed cells. Different computational approaches for this particular study are presented to calculate the solvation energies and binding free energies of H-Ras and PI3K. The goal of this study is to establish computational methods to investigate the roles of different proteins played in the cellular responses to space radiation, including modification of protein function through gene mutation, and to support the studies in molecular cell biology and theoretical kinetics models for our risk assessment project.

Ren, Lei↗
NASA NTRS2011

Ion-Ice Astrochemistry: Barrierless Low-Energy Deposition Pathways to HCOOH, CH3OH, and CO2 on Icy Grain Mantles from Precursor Cations

A new family of very favorable reaction pathways is explored involving the deposition of ions on icy grain mantles with very low energies. Quantum chemical cluster calculations at the MP2/6-31+G** level in 4H2O clusters and at the B3LYP/6-31+G** level in 17H2O clusters indicate that HCO+ and CH3 + are able to react spontaneously with one of the water molecules in the cluster to form protonated formic acid (HCOOH2 +) and protonated methanol (CH3OH2 +), respectively. It is furthermore found that these initial adducts spontaneously transfer their excess protons to the cluster to form neutral formic acid and methanol, plus solvated hydronium, H3O+. In the final case, if a CO molecule is bound to the surface of the cluster, OH+ may react with it to form protonated carbon dioxide (HCO2 +), which then loses its proton to yield CO2 and H3O+. In the present model, all of these processes were found to occur with no barriers. Discussion includes the analogous gas phase processes, which have been considered in previous studies, as well as the competitive abstraction pathway for HCO(+) + H2O.

Woon, David E.↗
NASA NTRS2002

Protein Solvation in Membranes and at Water-Membrane Interfaces

Different salvation properties of water and membranes mediate a host of biologically important processes, such as folding, insertion into a lipid bilayer, associations and functions of membrane proteins. These processes will be discussed in several examples involving synthetic and natural peptides. In particular, a mechanism by which a helical peptide becomes inserted into a model membrane will be described. Further, the molecular mechanism of recognition and association of protein helical segments in membranes will be discussed. These processes are crucial for proper functioning of a cell. A membrane-spanning domain of glycophorin A, which exists as a helical dimer, serves as the model system. For this system, the free energy of dissociation of the helices is being determined for both the wild type and a mutant, in which dimerization is disrupted.

Pohorille, Andrew↗
NASA NTRS1994

Interaction of a Model Peptide with a Water--Bilayer System

We discuss a molecular dynamics study of the alanine dipeptide at the interface between water and a glycerol-1-monooleate (GMO) bilayer. The dipeptide is interfacially active and incorporates into the bilayer without disrupting its structure. The interfacial region that is readily penetrated by the dipeptide spans the entire head group portion of the bilayer. The polar groups of the alanine dipeptide mostly remain well solvated by water and the oxygen atoms of GMO, and conformations of the dipeptide are characterized by (phi, psi) angles typical of alpha-helix and beta-sheet structures. When the molecule is deeper in the bilayer, the C(sub 7eq) state also becomes stable. The barrier to the isomerization reaction at the interface is lower than in bulk phases. After 7 ns of trajectories, the system is not fully equilibrated, due to slow collective motions involving GMO head groups. These result in decreased mobility and lower rates of isomerization of the dipeptide at the interface.

Pohorille, A.↗
NASA NTRS2013

The First X-ray Diffraction Patterns of Clay Minerals from Gale Crater

The Mars Science Laboratory (MSL) Rover, Curiosity spent approx 150 sols at Yellowknife Bay (YKB) studying a section of fluvio-lacustrine sedimentary rocks (with potential indications of volcanic influence), informally known as the Yellowknife Bay formation. YKB lies in a distal region of the Peace Vallis alluvial fan, which extends from the northern rim of Gale Crater toward the dune field at the base of Mt Sharp. Sedimentological and stratigraphic observations are consistent with the Yellowknife Bay formation being part of a distal fan deposit, which could be as young as middle Hesperian to even early Amazonian in age (approx 3.5 to 2.5 Ga). The Yellowknife Bay formation hosts a unit of mudstone called the Sheepbed member. Curiosity obtained powdered rock samples from two drill holes in the Sheepbed Member, named John Klein and Cumberland, and delivered them to instruments in Curiosity. Data from CheMin, a combined X-ray diffraction (XRD)/X-ray fluorescence instrument (XRF), has allowed detailed mineralogical analysis of mudstone powders revealing a clay mineral component of approx 20 wt.% in each sample. The clay minerals are important indicators of paleoenvironmental conditions and sensitive recorders of post-depositional alteration processes. The XRD pattern of John Klein reveals a 021 band consistent with a trioctahedral phyllosilicate. A broad peak at approx 10A with a slight inflexion at approx 12A indicates the presence of 2:1 type clay minerals in the John Klein sample. The trioctahedral nature of the clay minerals, breadth of the basal reflection, and presence of a minor component with larger basal spacing suggests that John Klein contains a trioctahedral smectite (probably saponite), whose interlayer is largely collapsed because of the low-humidity conditions. The XRD patterns show no evidence of corrensite (mixed-layer chlorite/smectite) or chlorite, which are typical diagenetic products of trioctahedral smectites when subjected to burial and heating >60degC in the presence of water. Given estimated geothermal gradients on Mars temperatures <60 degC might still be consistent with (but do not require) moderate burial. However, our ability to identify interstratified minerals is greatly limited by the lack of access to traditional treatments methods used in the lab (e.g., ethylene glycol solvation). Our preferred explanation for the origin of trioctahedral smectites in Sheepbed mudstone is in situ production via reaction of olivine, water and Si-bearing amorphous material, an important mudstone component detected by XRD. Elevated levels of magnetite in the Sheepbed and the trioctahedral monomineralic nature of the clay minerals support this model. These observations, combined with previous studies of olivine stability, support the persistence of circum-neutral hydrous conditions for thousands of years at YKB.

Bristow, Thomas↗
NASA NTRS2014

Paleoenvironmental Implications of Clay Minerals at Yellowknife Bay, Gale Crater, Mars

The Mars Science Laboratory (MSL) Rover, Curiosity spent approx 150 sols at Yellowknife Bay (YKB) studying a section of fluvio-lacustrine sedimentary rocks (with potential indications of volcanic influence), informally known as the Yellowknife Bay formation. YKB lies in a distal region of the Peace Vallis alluvial fan, which extends from the northern rim of Gale Crater toward the dune field at the base of Mt Sharp. Sedimentological and stratigraphic observations are consistent with the Yellowknife Bay formation being part of a distal fan deposit, which could be as young as middle Hesperian to even early Amazonian in age (approx. 3.5 to 2.5 Ga). The Yellowknife Bay formation hosts a unit of mudstone called the Sheepbed member. Curiosity obtained powdered rock samples from two drill holes in the Sheepbed Member, named John Klein and Cumberland, and delivered them to instruments in Curiosity. Data from CheMin, a combined X-ray diffraction (XRD)/X-ray fluorescence instrument (XRF), has allowed detailed mineralogical analysis of mudstone powders revealing a clay mineral component of approx. 20 wt.% in each sample. The clay minerals are important indicators of paleoenvironmental conditions and sensitive recorders of post-depositional alteration processes. The XRD pattern of John Klein reveals a 02l band consistent with a trioctahedral phyllosilicate. A broad peak at approx. 10A with a slight inflexion at approx. 12A indicates the presence of 2:1 type clay minerals in the John Klein sample. The trioctahedral nature of the clay minerals, breadth of the basal reflection, and presence of a minor component with larger basal spacing suggests that John Klein contains a trioctahedral smectite (probably saponite), whose interlayer is largely collapsed because of the low-humidity conditions. The XRD patterns show no evidence of corrensite (mixed-layer chlorite/smectite) or chlorite, which are typical diagenetic products of trioctahedral smectites when subjected to burial and heating >60 C in the presence of water. Given estimated geothermal gradients on Mars temperatures <60 C might still be consistent with (but do not require) moderate burial. However, our ability to identify interstratified minerals is greatly limited by the lack of access to traditional treatments methods used in the lab (e.g., ethylene glycol solvation). Our preferred explanation for the origin of trioctahedral smectites in Sheepbed mudstone is in situ production via reaction of olivine, water and Si-bearing amorphous material, an important mudstone component detected by XRD. Elevated levels of magnetite in the Sheepbed and the trioctahedral monomineralic nature of the clay minerals support this model. These observations, combined with previous studies of olivine stability, support the persistence of circum-neutral hydrous conditions for thousands of years at YKB.

MSL↗
NASA NTRS1997

Adsorption and solvation of ethanol at the water liquid-vapor interface: a molecular dynamics study

The free energy profiles of methanol and ethanol at the water liquid-vapor interface at 310K were calculated using molecular dynamics computer simulations. Both alcohols exhibit a pronounced free energy minimum at the interface and, therefore, have positive adsorption at this interface. The surface excess was computed from the Gibbs adsorption isotherm and was found to be in good agreement with experimental results. Neither compound exhibits a free energy barrier between the bulk and the surface adsorbed state. Scattering calculations of ethanol molecules from a gas phase thermal distribution indicate that the mass accommodation coefficient is 0.98, and the molecules become thermalized within 10 ps of striking the interface. It was determined that the formation of the solvation structure around the ethanol molecule at the interface is not the rate-determining step in its uptake into water droplets. The motion of an ethanol molecule in a water lamella was followed for 30 ns. The time evolution of the probability distribution of finding an ethanol molecule that was initially located at the interface is very well described by the diffusion equation on the free energy surface.

NASA Discipline Exobiology↗
NASA NTRS1996

Mechanism of unassisted ion transport across membrane bilayers

To establish how charged species move from water to the nonpolar membrane interior and to determine the energetic and structural effects accompanying this process, we performed molecular dynamics simulations of the transport of Na+ and Cl- across a lipid bilayer located between two water lamellae. The total length of molecular dynamics trajectories generated for each ion was 10 ns. Our simulations demonstrate that permeation of ions into the membrane is accompanied by the formation of deep, asymmetric thinning defects in the bilayer, whereby polar lipid head groups and water penetrate the nonpolar membrane interior. Once the ion crosses the midplane of the bilayer the deformation "switches sides"; the initial defect slowly relaxes, and a defect forms in the outgoing side of the bilayer. As a result, the ion remains well solvated during the process; the total number of oxygen atoms from water and lipid head groups in the first solvation shell remains constant. A similar membrane deformation is formed when the ion is instantaneously inserted into the interior of the bilayer. The formation of defects considerably lowers the free energy barrier to transfer of the ion across the bilayer and, consequently, increases the permeabilities of the membrane to ions, compared to the rigid, planar structure, by approximately 14 orders of magnitude. Our results have implications for drug delivery using liposomes and peptide insertion into membranes.

NASA Discipline Exobiology↗
NASA NTRS1990

Cavities in molecular liquids and the theory of hydrophobic solubilities

Thermal configurational data on neat liquids are used to obtain the work of formation of hard spherical cavities of atomic size in six molecular solvents: n-hexane, n-dodecane, n-undecyl alcohol, chloroform, carbon tetrachloride, and water. These results are used to test a recent suggestion that the differences between nonaqueous solvents and liquid water in solvation of inert gases are not principally due to the hydrogen-bonded structure of liquid water but rather to the comparatively small size of the water molecule. The frequencies of occurrence of cavities in liquid water can be meaningfully distinguished from those in the organic solvents. Liquid water has a larger fractional free volume, but that free volume is distributed in smaller packets. With respect to cavity work, water is compared to a solvent of the same molecular density and composed of hard spheres of the same size as the water molecule. That comparison indicates that the hard-sphere liquid finds more ways to configure its free volume in order to accommodate an atomic solute of substantial size and thus, would be more favorable solvent for inert gases. The scaled particle model of inert gas solubility in liquid water predicts cavity works 20% below the numerical data for TIP4P water at 300 K and 1.0 g/cm3 for cavity radii near 2.0 angstroms. It is argued that the sign of this difference is just the sign that ought to be expected and that the magnitude of this difference measures structural differences between water and the directly comparable hard-sphere liquid. In conjunction with previous data, these results indicate that atomic sized cavities should be considered submacroscopic.

NASA Discipline Exobiology↗
NASA NTRS2005

Structural and energetic consequences of mutations in a solvated hydrophobic cavity

The structural and energetic consequences of modifications to the hydrophobic cavity of interleukin 1-beta (IL-1beta) are described. Previous reports demonstrated that the entirely hydrophobic cavity of IL-1beta contains positionally disordered water. To gain a better understanding of the nature of this cavity and the water therein, a number of mutant proteins were constructed by site-directed mutagenesis, designed to result in altered hydrophobicity of the cavity. These mutations involve the replacement of specific phenylalanine residues, which circumscribe the cavity, with tyrosine, tryptophan, leucine and isoleucine. Using differential scanning calorimetry to determine the relative stabilities of the wild-type and mutant proteins, we found all of the mutants to be destabilizing. X-ray crystallography was used to identify the structural consequences of the mutations. No clear correlation between the hydrophobicities of the specific side-chains introduced and the resulting stabilities was found.

Solvents/chemistry↗
NASA NTRS1997

Interactions of anesthetics with the water-hexane interface. A molecular dynamics study

The free energy profiles characterizing the transfer of nine solutes across the liquid-vapor interfaces of water and hexane and across the water-hexane interface were calculated from molecular dynamics simulations. Among the solutes were n-butane and three of its halogenated derivatives, as well as three halogenated cyclobutanes. The two remaining molecules, dichlorodifluoromethane and 1,2-dichloroperfluoroethane, belong to series of halo-substituted methanes and ethanes, described in previous studies (J. Chem. Phys. 1996, 104, 3760; Chem. Phys. 1996, 204, 337). Each series of molecules contains structurally similar compounds that differ greatly in anesthetic potency. The accuracy of the simulations was tested by comparing the calculated and the experimental free energies of solvation of all nine compounds in water and in hexane. In addition. the calculated and the measured surface excess concentrations of n-butane at the water liquid-vapor interface were compared. In all cases, good agreement with experimental results was found. At the water-hexane interface, the free energy profiles for polar molecules exhibited significant interfacial minima, whereas the profiles for nonpolar molecules did not. The existence of these minima was interpreted in terms of a balance between the free energy contribution arising from solute-solvent interactions and the work to form a cavity that accommodates the solute. These two contributions change monotonically, but oppositely, across the interface. The interfacial solubilities of the solutes, obtained from the free energy profiles, correlate very well with their anesthetic potencies. This is the case even when the Meyer-Overton hypothesis, which predicts a correlation between anesthetic potency and solubility in oil, fails.

NASA Center ARC↗