Search NASA⌕ Search

SEARCH · Search NASA

Results for “CHARMM”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Update 1.1 to “pysimm: A python package for simulation of molecular systems”, (PII: S2352711016300395)

This update of the pysimm application programming interface, pysimm 1.1, provides both infrastructural as well as functional updates. Moreover, improvements to the random walk application that allow it to construct polymers with controlled tacticity are highlighted. Additions to the forcefield module include an update to enable working with the family of CHARMM forcefields and automated typing with the CHARMM generalized forcefield (CGenFF). Finally, new detailed examples demonstrating new features are also provided.

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

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials

Abstract The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

Kanhaiya, Krishan (ORCID:0000000236222655)↗
DOE OSTI2017

Glycan Reader is improved to recognize most sugar types and chemical modifications in the Protein Data Bank

Motivation: Glycans play a central role in many essential biological processes. Glycan Reader was originally developed to simplify the reading of Protein Data Bank (PDB) files containing glycans through the automatic detection and annotation of sugars and glycosidic linkages between sugar units and to proteins, all based on atomic coordinates and connectivity information. Carbohydrates can have various chemical modifications at different positions, making their chemical space much diverse. Unfortunately, current PDB files do not provide exact annotations for most carbohydrate derivatives and more than 50% of PDB glycan chains have at least one carbohydrate derivative that could not be correctly recognized by the original Glycan Reader. Results: Glycan Reader has been improved and now identifies most sugar types and chemical modifications (including various glycolipids) in the PDB, and both PDB and PDBx/mmCIF formats are supported. CHARMM-GUI Glycan Reader is updated to generate the simulation system and input of various glycoconjugates with most sugar types and chemical modifications. It also offers a new functionality to edit the glycan structures through addition/deletion/modification of glycosylation types, sugar types, chemical modifications, glycosidic linkages, and anomeric states. The simulation system and input files can be used for CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Glycan Fragment Database in GlycanStructure. Org is also updated to provide an intuitive glycan sequence search tool for complex glycan structures with various chemical modifications in the PDB.

Park, Sang-Jun↗
DOE OSTI2018

Lignin Forcefield Source Datafiles

Three archives are provided here comprising the development of a CHARMM36-compatible lignin forcefield suitable for use in biomolecular simulations. One archive parameters.zip contains the output topology and parameter files in CHARMM format suitable for use in constructing simulation systems. The remaining two hold the logic of the parameterization featuring scripts that do the optimization including the source for the GPU-accelerated objective function evaluation and target data creation scripts.zip as well as minimal outputs needed to recreate the work datafiles.zip. Due to licensing restrictions Gaussian log files that created the target data are not provided here. Instead only binary-formatted parsed target data as well as the input decks used to create the target data are provided in datafiles.zip. Readme files for individual directories created within the directory structure are contained within scripts.zip.

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

Force field for calcium sulfate minerals to predict structural, hydration, and interfacial properties

Calcium sulfates such as anhydrite, hemihydrate, and gypsum find widespread use in building materials, implants, and tissue healing. We introduce a simple and compatible atomistic force field for all calcium sulfate phases that reproduces a wide range of experimental data including lattice parameters, surface, hydration, mechanical, and thermal properties in 1% to 5% accuracy relative to experiments. The performance is several times better than prior force fields and DFT methods, which lead to errors in structures and energies up to 100%. We explain (hkl) cleavage energies, the dynamics of (hkl) water interfaces, and new insights into molecular origins of crystal-facet specific hydration and solubility. Impressive agreement of computed and experimentally measured hydration energies is shown. The models add to the Interface force field (IFF) and are compatible with multiple force fields (CHARMM, AMBER, GROMOS, CVFF, PCFF, OPLS-AA) for property predictions of sulfate-containing materials from atoms to the large nanometer scale.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Redox potential tuning through differential quinone binding in the photosynthetic reaction center of Rhodobacter sphaeroides

Ubiquinone forms an integral part of the electron transport chain in cellular respiration and photosynthesis across a vast number of organisms. Prior experimental results have shown that the photosynthetic reaction center (RC) from Rhodobacter sphaeroides is only fully functional with a limited set of methoxy-bearing quinones, suggesting that specific interactions with this substituent are required to drive electron transport and the formation of quinol. The nature of these interactions has yet to be determined. Through parameterization of a CHARMM-compatible quinone force field and subsequent molecular dynamics simulations of the quinone-bound RC, in this paper we have investigated and characterized the interactions of the protein with the quinones in the Q A and Q B sites using both equilibrium simulation and thermodynamic integration. In particular, we identify a specific interaction between the 2-methoxy group of ubiquinone in the Q B site and the amide nitrogen of GlyL225 that we implicate in locking the orientation of the 2-methoxy group, thereby tuning the redox potential difference between the quinones occupying the Q A and Q B sites. Finally, disruption of this interaction leads to weaker binding in a ubiquinone analogue that lacks a 2-methoxy group, a finding supported by reverse electron transfer electron paramagnetic resonance experiments of the Q A –Q B – biradical and competitive binding assays.

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

Nanopolysaccharide Builder: A User-Friendly Tool for Atomistic Models of Polysaccharide-Based Nanostructures

Here, we introduce Nanopolysaccharide Builder (NPB), a user-friendly software tool designed to construct polysaccharide nanostructures─mainly those based on cellulose, chitin, and chitosan─using experimental data or user-defined parameters. NPB enables the generation of cellulose and chitin allomorphs with customizable biochemical topologies and also facilitates the construction of large bundles that replicate nanostructures found in biological support systems, including plant cell walls and arthropod cuticles. The software outputs atomic Cartesian coordinates in Protein Data Bank (PDB) format and also provides atom connectivity files in PSF and PARM formats, ensuring seamless integration with major molecular dynamics (MD) engines such as NAMD, CHARMM, GROMACS, AMBER, OpenMM, and LAMMPS. Built on an interactive visualization framework, NPB features a graphical user interface (GUI) and supports both macOS and Linux operating systems. By enabling detailed atomic-scale studies of polysaccharide evolution in extracellular matrices and cell walls of algae, bacteria, fungi, and plants, NPB is poised to advance AI-guided research in sustainable chemical development and biomass utilization.

Wan, Zhangmin [Univ. of British Columbia, Vancouve↗
DOE OSTI2020

Molecular Structure of Sphingomyelin in Fluid Phase Bilayers Determined by the Joint Analysis of Small-Angle Neutron and X-ray Scattering Data

In this work, we have determined the fluid bilayer structure of palmitoyl sphingomyelin (PSM) and stearoyl sphingomyelin (SSM) by simultaneously analyzing small-angle neutron and X-ray scattering data. Using a newly developed scattering density profile (SDP) model for sphingomyelin lipids, we report structural parameters including the area per lipid, total bilayer thickness, and hydrocarbon thickness, in addition to lipid volumes determined by densitometry. Unconstrained all-atom simulations of PSM bilayers at 55 °C using the C36 CHARMM force field produced a lipid area of 56 Å 2 , a value that is 10% lower than the one determined experimentally by SDP analysis (61.9 Å2). Furthermore, scattering form factors calculated from the unconstrained simulations were in poor agreement with experimental form factors, even though segmental order parameter (S CD ) profiles calculated from the simulations were in relatively good agreement with S CD profiles obtained from NMR experiments. Conversely, constrained area simulations at 61.9 Å 2 resulted in good agreement between the simulation and experimental scattering form factors, but not with S CD profiles from NMR. We discuss possible reasons for the discrepancies between these two types of data that are frequently used as validation metrics for molecular dynamics force fields.

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

Characterization of Aβ Monomers through the Convergence of Ensemble Properties among Simulations with Multiple Force Fields

We found that amyloid β (Aβ) monomers represent a base state in the pathways of aggregation that result in the fibrils and oligomers implicated in the pathogenesis of Alzheimer’s disease (AD). The structural properties of these intrinsically disordered peptides remain unclear despite extensive experimental and computational investigations. Further, there are mutations within Aβ that change the way the peptide aggregates and are known to cause familial AD (FAD). Here, we analyze the ensembles of different isoforms (Aβ42 and Aβ40) and mutants (E22Δ, D23N, E22K, E22G, and A2T in Aβ40) of Aβ generated with all-atom replica exchange molecular dynamics (REMD) simulations on the μs/replica time scale. These were run using three different force field/water model combinations: OPLS-AA/L and TIP3P (“OPLS”), AMBER99sb-ILDN and TIP4P-Ew (“ILDN”), as well as CHARMM22* and TIP3SP (“CHARMM”). Despite fundamental changes in simulation parameters, we find that the resulting ensembles demonstrate a strong convergence in structural properties. In particular, antiparallel contacts between L17–A21 and A30–L34 are prevalent in ensembles of Aβ40, directly forming β sheets in the OPLS and ILDN combinations. A21–A30 commonly forms an interceding region that rarely interacts with the rest of the peptide. Further, Aβ42 contributes new β hairpin motifs involving V40–I41 in both OPLS and ILDN. However, the structural flexibility of the central region and the electrostatic interactions that characterize it are notably different between the different conditions. Further, for OPLS, each of the FAD mutations disrupts central bend character and increases the polymorphism of antiparallel contacts across the central region. However, the studied mutations in the ILDN set primarily encourage more global contacts involving the N-terminus and the central region, and promote the formation of new β topologies that may seed different aggregates involved in disease phenotypes. Furthermore, these differences aside, the large degree of agreement between simulation sets across multiple force fields provides a generalizable characterization of Aβ that is also consistent with experimental data and models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

All-Atom Molecular Dynamics Study of Water–Dodecane Interface in the Presence of Octanol

In this work, classical atomistic molecular dynamics (MD) simulations were performed to investigate the structural properties of octanol at water/dodecane interface. Multiple combinations of simulation conditions and force fields were employed in order to :obtain high fidelity interfacial structure. The analysis of noncovalent interactions shows that the convergence of interfacial octanol distribution requires a simulation time scale larger than 100 ns. Finite size effect diminishes with increasing dimension of simulation box both laterally and vertically. In the largest simulated systems with vanishing finite size effect, the additive CHARMM force field shows similar interfacial aggregation behavior with the DRUDE polarizable model. Self-assembly of octanol molecules at water/oil interface is observed, expected to reduce;the free energy barrier between water phase and oil phase. This work provides an MD simulation benchmark for structural study of complex liquid/liquid interface with cosurfactant reducing interfacial tension between two immiscible phases.

Qiao, Baofu↗
DOE OSTI2018

Understanding Chemical Bonding in Alloys and the Representation in Atomistic Simulations

Alloys are widely used in catalysts and structural materials. The nature of chemical bonding and the origin of alloy formation energies, defect energies, and interfacial properties have not been well understood to date but are critical to material performance. In this contribution, we explain the polar nature of chemical bonding and an implementation in classical and reactive atomistic simulations to understand such properties more quantitatively. Electronegativity differences between metal atoms lead to polar bonding, and exothermic alloy formation energies are related to charge transfer between the different elements. These differences can be quantified by atomic charges using pairwise charge increments, determined by matching the computed alloy formation energy to experimentally measured alloy formation energies using pair potentials for the pure metals. The polar character of alloys is comparable to organic molecules and partially ionic minerals, for example, AlNi and AlNi3 alloys assume significant atomic charges of ±0.40e and +0.60e/–0.20e, respectively. The subsequent analysis of defect sites and defect energies using force-field-based calculations shows excellent agreement with calculations using density functional theory and embedded atom models (EAM). The formation of vacancy and antisite defects is characterized by a redistribution of charge in the first shell of neighbor atoms in the classical models whereby electroneutrality is maintained and charge increments correlate with differences in electronegativity. The proposed atomic charges represent internal dipole and multipole moments, consistent with existing definitions for organic and inorganic compounds and with the extended Born model (Heinz, H.; Suter, U. W. J. Phys. Chem. B2004,108 (47), 18341–18352). The method can be applied to any alloy and has a reproducibility of ±10%. In contrast, quantum mechanical charge schemes remain associated with deviations exceeding ±100%. The atomic charges for alloys provide a simple initial measure for the internal electronic structure, surface adsorption of molecules, and reactivity in catalysis and corrosion. The models are compatible with the Interface force field (IFF), CHARMM, AMBER, OPLS-AA, PCFF, CVFF, and GROMOS for reliable atomistic simulations of alloys and their interfaces with minerals and electrolytes from the nanometer scale to the micrometer scale.

Chemistry↗
DOE OSTI2026

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
DOE OSTI2023

Performance efficient macromolecular mechanics via sub-nanometer shape based coarse graining

Dimensionality reduction via coarse grain modeling is a valuable tool in biomolecular research. For large assemblies, ultra coarse models are often knowledge-based, relying on a priori information to parameterize models thus hindering general predictive capability. Here, we present substantial advances to the shape based coarse graining (SBCG) method, which we refer to as SBCG2. SBCG2 utilizes a revitalized formulation of the topology representing network which makes high-granularity modeling possible, preserving atomistic details that maintain assembly characteristics. Further, we present a method of granularity selection based on charge density Fourier Shell Correlation and have additionally developed a refinement method to optimize, adjust and validate high-granularity models. We demonstrate our approach with the conical HIV-1 capsid and heteromultimeric cofilin-2 bound actin filaments. Our approach is available in the Visual Molecular Dynamics (VMD) software suite, and employs a CHARMM-compatible Hamiltonian that enables high-performance simulation in the GPU-resident NAMD3 molecular dynamics engine.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Potential of mean force analysis of the self-association of leucine-rich transmembrane α-helices: Difference between atomistic and coarse-grained simulations

Interaction of transmembrane (TM) proteins is important in many biological processes. Large-scale computational studies using coarse-grained (CG) simulations are becoming popular. However, most CG model parameters have not fully been calibrated with respect to lateral interactions of TM peptide segments. Here, we compare the potential of mean forces (PMFs) of dimerization of TM helices obtained using a MARTINI CG model and an atomistic (AT) Berger lipids-OPLS/AA model (AT{sup OPLS}). For helical, tryptophan-flanked, leucine-rich peptides (WL15 and WALP15) embedded in a parallel configuration in an octane slab, the AT{sup OPLS} PMF profiles showed a shallow minimum (with a depth of approximately 3 kJ/mol; i.e., a weak tendency to dimerize). A similar analysis using the CHARMM36 all-atom model (AT{sup CHARMM}) showed comparable results. In contrast, the CG analysis generally showed steep PMF curves with depths of approximately 16–22 kJ/mol, suggesting a stronger tendency to dimerize compared to the AT model. This CG > AT discrepancy in the propensity for dimerization was also seen for dilauroylphosphatidylcholine (DLPC)-embedded peptides. For a WL15 (and WALP15)/DLPC bilayer system, AT{sup OPLS} PMF showed a repulsive mean force for a wide range of interhelical distances, in contrast to the attractive forces observed in the octane system. The change from the octane slab to the DLPC bilayer also mitigated the dimerization propensity in the CG system. The dimerization energies of CG (AALALAA){sub 3} peptides in DLPC and dioleoylphosphatidylcholine bilayers were in good agreement with previous experimental data. The lipid headgroup, but not the length of the lipid tails, was a key causative factor contributing to the differences between octane and DLPC. Furthermore, the CG model, but not the AT model, showed high sensitivity to changes in amino acid residues located near the lipid-water interface and hydrophobic mismatch between the peptides and membrane. These findings may help interpret CG and AT simulation results on membrane proteins.

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

Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD{sub 3}CN treated with a parallel multi-state EVB model

We describe a parallelized linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM and TINKER. Forces are obtained using the Hellmann-Feynman relationship, giving continuous gradients, and good energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to explicitly correlated coupled cluster theory, we built a 64-state MS-EVB model designed to study the F + CD{sub 3}CN → DF + CD{sub 2}CN reaction in CD{sub 3}CN solvent (recently reported in Dunning et al. [Science 347(6221), 530 (2015)]). This approach allows us to build a reactive potential energy surface whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We ran molecular dynamics simulations to examine a range of observables which follow in the wake of the reactive event: energy deposition in the nascent reaction products, vibrational relaxation rates of excited DF in CD{sub 3}CN solvent, equilibrium power spectra of DF in CD{sub 3}CN, and time dependent spectral shifts associated with relaxation of the nascent DF. Many of our results are in good agreement with time-resolved experimental observations, providing evidence for the accuracy of our MS-EVB framework in treating both the solute and solute/solvent interactions. The simulations provide additional insight into the dynamics at sub-picosecond time scales that are difficult to resolve experimentally. In particular, the simulations show that (immediately following deuterium abstraction) the nascent DF finds itself in a non-equilibrium regime in two different respects: (1) it is highly vibrationally excited, with ∼23 kcal mol{sup −1} localized in the stretch and (2) its post-reaction solvation environment, in which it is not yet hydrogen-bonded to CD{sub 3}CN solvent molecules, is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of the post-reaction solvation environment results in a red shift. These two competing effects mean that the post-reaction relaxation profile is distinct from what is observed when Franck-Condon vibrational excitation of DF occurs within a microsolvation environment initially at equilibrium. Our conclusions, along with the theoretical and parallel software framework presented in this paper, should be more broadly applicable to a range of complex reactive systems.

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

Polarizable charge equilibration model for predicting accurate electrostatic interactions in molecules and solids

Electrostatic interactions play a critical role in determining the properties, structures, and dynamics of chemical, biochemical, and material systems. These interactions are described well at the level of quantum mechanics (QM) but not so well for the various models used in force field simulations of these systems. In this paper, we propose and validate a new general methodology, denoted PQEq, to predict rapidly and dynamically the atomic charges and polarization underlying the electrostatic interactions. Here the polarization is described using an atomic sized Gaussian shaped electron density that can polarize away from the core in response to internal and external electric fields, while at the same time adjusting the charge on each core (described as a Gaussian function) so as to achieve a constant chemical potential across all atoms of the system. The parameters for PQEq are derived from experimental atomic properties of all elements up to Nobelium (atomic no. = 102). Additionally, we validate PQEq by comparing to QM interaction energy as probe dipoles are brought along various directions up to 30 molecules containing H, C, N, O, F, Si, P, S, and Cl atoms. We find that PQEq predicts interaction energies in excellent agreement with QM, much better than other common charge models such as obtained from QM using Mulliken or ESP charges and those from standard force fields (OPLS and AMBER). Since PQEq increases the accuracy of electrostatic interactions and the response to external electric fields, we expect that PQEq will be useful for a large range of applications including ligand docking to proteins, catalytic reactions, electrocatalysis, ferroelectrics, and growth of ceramics and films, where it could be incorporated into standard force fields as OPLS, AMBER, CHARMM, Dreiding, ReaxFF, and UFF.

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

Electrostatic and induction effects in the solubility of water in alkanes

Experiments show that at 298 K and 1 atm pressure, the transfer free energy, μ ex , of water from its vapor to liquid normal alkanes C n H 2n+2 (n=5…12) is negative. Earlier it was found that with the united-atom TraPPE model for alkanes and the SPC/E model for water, one had to artificially enhance the attractive alkane-water cross interaction to capture this behavior. Here we revisit the calculation of μ ex using the polarizable AMOEBA and the non-polarizable Charmm General (CGenFF) forcefields. We test both the AMOEBA03 and AMOEBA14 water models; the former has been validated with the AMOEBA alkane model while the latter is a revision of AMOEBA03 to better describe liquid water. We calculate μ ex using the test particle method. With CGenFF, μ ex is positive and the error relative to experiments is about 1.5 k B T. With AMOEBA, μ ex is negative and deviations relative to experiments are between 0.25 k B T (AMOEBA14) and 0.5 k B T (AMOEBA03). Quantum chemical calculations in a continuum solvent suggest that zero point effects may account for some of the deviation. In conclusion, forcefield limitations notwithstanding, electrostatic and induction effects, commonly ignored in consideration of water-alkane interactions, appear to be decisive in the solubility of water in alkanes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗