Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical calculations, Polarization”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Martian terrains

Terrain studies of candidate landing sites for a future rover/sample-return mission to Mars are being conducted to evaluate the geologic and trafficability aspects of each site. An optimum site should have geologic units of widely diverse ages and chemical compositions occurring in close enough proximity and in smooth enough terrain so that a roving vehicle of limited traverse ability (+ or - 100 km) could collect representative samples. In FY 1986, geologic maps were compiled at 1:500,000 and 1:2 million scales of the Mangala Valles, Kasei Valles, Chasma Boreale (north polar), and Planum Australe (south polar) areas, and a study was begun of the topography and surface roughness characteristics of the Mangala Valles site. Geologic mapping has been greatly facilitated by specially enhanced, high-resolution Viking photographs, which clarify stratigraphic relations of units unrecognized earlier. Photoclinometric profiles of topographic features provide width and depth measurements of four classes of channels, the thickness of some volcanic units, and the throw on some faults. Estimates of the surface roughness of units are calculated using a newly developed USGS computer program and using measurements derived from Earth-based radar.

Masursky, Harold↗
DOE OSTI2021

Widespread Negative Longitudinal Piezoelectric Responses in Ferroelectric Crystals with Layered Structures

In this study, we investigate the underlying mechanisms of the universal negative piezoelectricity in low-dimensional layered materials by carrying out first-principles calculations. Two-dimensional layered ferroelectric CuInP 2 S 6 is analyzed in detail as a typical example, but the theory can be applied to any other low-dimensional layered piezoelectrics. Consistent with the theory proposed in [Phys. Rev. Lett. 119, 207601 (2017), the anomalous negative piezoelectricity in CuInP 2 S 6 also results from its negative clamped-ion term, which cannot be compensated by the positive internal-strain part. Here, we focus on a more general rule by proposing that having a negative clamped-ion term should be universal among piezoelectric materials, which is attributed to the “lag of Wannier center” effect. The internal-strain term, which is the change in polarization due to structural relaxation in response to strain, is mostly determined by the spatial structure and chemical bonding of the material. In a low-dimensional layered piezoelectric material such as CuInP 2 S 6 , the internal-strain term is approximately zero. This is because the internal structure of the molecular layers, which are bonded by the weak van der Waals interaction, responds little to the strain. As a result, the magnitude of the dipole, which depends strongly on the dimension and structure of the molecular layer, also has a small response with respect to strain. An equation bridging the internal strain responses in low-dimensional and three-dimensional piezoelectrics is also derived to analytically express this point. This work aims to deepen our understanding about this anomalous piezoelectric effect, especially in low-dimensional layered materials, and provide strategies for discovering materials with novel electromechanical properties.

2-dimensional systems↗
NASA NTRS2002

Gas-Phase Reaction Pathways and Rate Coefficients for the Dichlorosilane-Hydrogen and Trichlorosilane-Hydrogen Systems

As part of NASA Ames Research Center's Integrated Process Team on Device/Process Modeling and Nanotechnology our goal is to create/contribute to a gas-phase chemical database for use in modeling microelectronics devices. In particular, we use ab initio methods to determine chemical reaction pathways and to evaluate reaction rate coefficients. Our initial studies concern reactions involved in the dichlorosilane-hydrogen (SiCl2H2--H2) and trichlorosilane-hydrogen (SiCl2H-H2) systems. Reactant, saddle point (transition state), and product geometries and their vibrational harmonic frequencies are determined using the complete-active-space self-consistent-field (CASSCF) electronic structure method with the correlation consistent polarized valence double-zeta basis set (cc-pVDZ). Reaction pathways are constructed by following the imaginary frequency mode of the saddle point to both the reactant and product. Accurate energetics are determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations (CCSD(T)) extrapolated to the complete basis set limit. Using the data from the electronic structure calculations, reaction rate coefficients are obtained using conventional and variational transition state and RRKM theories.

Dateo, Christopher E.↗
NASA NTRS2001

Rotating reverse osmosis: a dynamic model for flux and rejection

Reverse osmosis (RO) is a compact process for the removal of ionic and organic pollutants from contaminated water. However, flux decline and rejection deterioration due to concentration polarization and membrane fouling hinders the application of RO technology. In this study, a rotating cylindrical RO membrane is theoretically investigated as a novel method to reduce polarization and fouling. A dynamic model based on RO membrane transport incorporating concentration polarization is used to predict the performance of rotating RO system. Operating parameters such as rotational speed and transmembrane pressure play an important role in determining the flux and rejection in rotating RO. For a given geometry, a rotational speed sufficient to generate Taylor vortices in the annulus is essential to maintain high flux as well as high rejection. The flux and rejection were calculated for wide range of operating pressures and rotational speeds. c 2001 Elsevier Science B.V. All rights reserved.

Non-NASA Center↗
NASA NTRS1999

Computational Study of Field Initiated Surface Reactions for Synthesis of Diamond and Silicon

This project involves using quantum chemistry to simulate surface chemical reactions in the presence of an electric field for nanofabrication of diamond and silicon. A field delivered by a scanning tunneling microscope (STM) to a nanometer scale region of a surface affects chemical reaction potential energy surfaces (PES) to direct atomic scale surface modification to fabricate sub-nanometer structures. Our original hypothesis is that the applied voltage polarizes the charge distribution of the valence electrons and that these distorted molecular orbitals can be manipulated with the STM so as to change the relative stabilities of the electronic configurations over the reaction coordinates and thus the topology of the PES and reaction kinetics. Our objective is to investigate the effect of applied bias on surface reactions and the extent to which STM delivered fields can be used to direct surface chemical reactions on an atomic scale on diamond and silicon. To analyze the fundamentals of field induced chemistry and to investigate the application of this technique for the fabrication of nanostructures, we have employed methods capable of accurately describing molecular electronic structure. The methods we employ are density functional theory (DFT) quantum chemical (QC) methods. To determine the effect of applied bias on surface reactions we have calculated the QC PESs in various applied external fields for various reaction steps for depositing or etching diamond and silicon. We have chosen reactions which are thought to play a role in etching and the chemical vapor deposition growth of Si and diamond. The PESs of the elementary reaction steps involved are then calculated under the applied fields, which we vary in magnitude and configuration. We pay special attention to the change in the reaction barriers, and transition state locations, and search for low energy reaction channels which were inaccessible without the applied bias.

Musgrave, Charles Bruce↗
DOE OSTI2020

Enhanced Menshutkin S N 2 Reactivity in Mesoporous Silica: The Influence of Surface Catalysis and Confinement

A significant enhancement in the Menshutkin S N 2 reaction between 1-methylimidazole (MeIm) and methyl thiocyanate (MeSCN) is observed when the reaction is confined in the nanoscale silica pores of MCM41 and SBA15. The experiments in the silica pores are conducted without surrounding bulk reaction mixture. The influences of temperature, pore radius, and surface chemistry on the kinetics of the confined reaction are analyzed with time-dependent infrared spectroscopy, molecular dynamics simulations, and ab initio calculations. The rate constant of the pseudo-first order reaction increases with decreasing pore size, and the activation energy is found to decrease by 5.6 kJ/mol in the smallest pore studied (2.8 nm) relative to the bulk reaction. The rate constant dependence on pore size is accurately described by a two-state model in which molecules within the 4.6 Å interfacial layer experience a 2.4-fold rate constant increase relative to those reacting at the bulk rate further away from the interface. The removal of polar silanol groups from the silica surface via passivation with trimethylsilyl chloride results in bulk-like kinetics despite a reduction in the pore diameter, demonstrating the role of silanols as catalytic sites. Here, electronic structure calculations of the energy profile on a model silica surface confirm that silanol groups, particularly those of the vicinal type, can reduce the activation energy and reaction endothermicity through the donation of hydrogen bonds to the reactant, transition state, and product complexes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electronic origin of half-metal to semiconductor transition and colossal magnetoresistance in spinel HgCr 2 Se 4

Half metals are ferromagnets hosting spin-polarized conducting carriers and are crucial for spintronics applications. The chromium spinel HgCr 2 Se 4 represents a unique type of half-metal, which features a half-metal to semiconductor transition (HMST) and exhibits colossal magnetoresistance (CMR) across the ferromagnetic-paramagnetic (FM-PM) transition. Here, using angle-resolved photoemission spectroscopy, we find that the Fermi surface of n-type HgCr 2 Se 4 (n–HgCr 2 Se 4 ) consists of a single electron pocket which moves above the Fermi level (E F ) upon the FM-PM transition, leading to the HMST. Such a Lifshitz transition manifests a giant band splitting which originates from the exchange interaction unveiled with a specific chemical nonstoichiometry. The exchange band splitting and the chemical nonstoichiometry are two key ingredients to the HMST and CMR, consistent with our ab initio calculation. Our findings provide spectroscopic evidence of the electronic origin of the anomalous properties of HgCr 2 Se 4 , which address the unique phase transition in half-metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Metastable piezoelectric group-IV monochalcogenide monolayers with a buckled honeycomb structure

Multiple two-dimensional materials are being naïvely termed stable on the grounds of displaying phonon dispersions with no negative frequencies and of not collapsing on molecular dynamics calculations at fixed volume. But, if these phases do not possess the smallest possible structural energy, how does one understand and establish their actual meta stability? To answer this question, twelve two-dimensional group-IV monochalcogenide monolayers (SiS, SiSe, SiTe, GeS, GeSe, GeTe, SnS, SnSe, SnTe, PbS, PbSe, and PbTe) with a buckled honeycomb atomistic structure—belonging to symmetry group P3m1—displaying an out-of-plane intrinsic electric polarization are shown to be metastable by three independent methods. First, we uncover a coordination-preserving structural transformation from the low-buckled honeycomb structure onto the lower-energy Pnm2 1 (or Pmmn for PbS, PbSe, and PbTe) phase to estimate energy barriers E B that must be overcome during such structural transformation. Using the curvature of the local minima and E B as inputs to Kramers escape formula, large escape times are found, implying the structural metastability of the buckled honeycomb phase (with the exception of PbS and PbSe, these phases display escape times ranging from 700 years to multiple times the age of the universe and can be considered “stable” for practical purposes in that relative sense). The second demonstration is provided by phonon dispersion relations that include the effect of long-range Coulomb forces and display no negative vibrational modes. The third and final demonstration of structural metastability is furnished by room-temperature ab initio molecular dynamics for selected compounds. Here, the magnitude of the electronic band gap evolves with chemical composition. Different from other binary two-dimensional compounds such as transition metal dichalcogenide monolayers and hexagonal boron nitride monolayers which only develop an in-plane piezoelectric response, the twelve group-IV monochalcogenide monolayers with a buckled honeycomb structure also display out-of-plane piezoelectric properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Linear Combination of Atomic Dipoles to Calculate the Bond and Molecular Dipole Moments of Molecules and Molecular Liquids

Here we report a linear combination of atomic dipole (LCAD) method for calculating the bond dipole moments of molecules. We show that the LCAD method reproduces the known molecular dipole moments of small to large molecules with a small error with respect to experimental and benchmark ab initio calculations, and molecular dipole distributions of bulk water that agree with maximally localized Wannier functions. The bond dipole moments derived from LCAD are also chemically interpretable in terms of the trend in bond ionicity in going from neutral to charged molecules. Moreover, the LCAD method accurately captures the influence of electric fields, supported by the correct trend in the change of the dipole moment under a uniform external electric field. The better grounding of bond dipole calculations indicates that it should also serve as a useful approach to bond dipole-field models used in catalysis or to reconstruct the small dipole of a H-terminated graphene flake.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗
NASA NTRS1989

Photochemistry of the atmosphere of Uranus

A detailed one-dimensional photochemical model incorporating vertical transport by eddy and molecular diffusion is used to study the importance of chemical and physical processes that control the hydrocarbon abundances in the stratosphere of Uranus. The study includes updated photochemical reaction rates and photolysis branching ratios; explicit calculation of the condensation rate of condensing species; and appropriate boundary conditions at the tropopause. It is found that the basic photochemical process generally produces low net conversion efficiencies of methane to higher hydrocarbons. However, even with these low efficiencies, the total condensation rate, which represents the net rate of conversion of methane to hydrocarbons, is roughly 2 x 10 to the -16th g/cm sec, in agreement with the upper range of the low condensation rate for the production of the stratospheric aerosol layer as determined by Pollack et al. (1987). Furthermore, a range of vertical transport rates is required to separately match the observational constraints on hydrocarbon abundances. The solar reflection data refer to the sunlit polar region of Uranus and require an eddy diffusion coefficient of K approx. = 50 sq cm/sec for globally averaged insolation or K approx = 1500 (Nh/N) exp 1/2 for overhead sun. The Voyager 2 Utraviolet Spectrometer solar occultation data refer to the equatorial regions, thus implying a large latitudinal gradient in stratospheric hydrocarbon abundance.

Summers, Michael E.↗
DOE OSTI2021

Chemically and electrically tunable spin polarization in ferroelectric Cd-based hybrid organic-inorganic perovskites

Density functional theory computations are used to predict the electronic structure of Cd-based hybrid organic-inorganic high-T C ferroelectric perovskites with TMCM-CdCl 3 being one representative. Here we report Rashba-Dresselhaus spin splitting in the valence band of these nonmagnetic compounds. Interestingly, we find in computations that the splitting is not necessarily sensitive to the polarization of the material but to the organic molecule itself which opens a way to its chemical tunability through the choice of the molecule. Further chemical tunability of spin splitting is shown to be possible through a substitution of Cl in the CdCl 3 chains as the valence band was found to originate from Cl-Cl weekly bonding orbitals. For example, the substitution of Cl with Br in TMCM-CdCl 3 resulted in a ten times increase of spin splitting. Furthermore, the spin polarization in these materials gives origin to persistent spin textures which are coupled to the polarization direction, and, therefore, can be controlled by the electric field. This is promising for spintronics applications.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Ir 6 In 32 S 21 , a polar, metal-rich semiconducting subchalcogenide

Subchalcogenides are uncommon, and their chemical bonding results from an interplay between metal–metal and metal–chalcogenide interactions. Herein, we present Ir 6 In 32 S 21 , a novel semiconducting subchalcogenide compound that crystallizes in a new structure type in the polar P31m space group, with unit cell parameters a = 13.9378(12) Å, c = 8.2316(8) Å, α = β = 90°, γ = 120°. The compound has a large band gap of 1.48(2) eV, and photoemission and Kelvin probe measurements corroborate this semiconducting behavior with a valence band maximum (VBM) of –4.95(5) eV, conduction band minimum of –3.47(5) eV, and a photoresponse shift of the Fermi level by ~0.2 eV in the presence of white light. X-ray absorption spectroscopy shows absorption edges for In and Ir do not indicate clear oxidation states, suggesting that the numerous coordination environments of Ir 6 In 32 S 21 make such assignments ambiguous. Electronic structure calculations confirm the semiconducting character with a nearly direct band gap, and electron localization function (ELF) analysis suggests that the origin of the gap is the result of electron transfer from the In atoms to the S 3p and Ir 5d orbitals. DFT calculations indicate that the average hole effective masses near the VBM (1.19m e ) are substantially smaller than the average electron masses near the CBM (2.51m e ), an unusual feature for most semiconductors. The crystal and electronic structure of Ir 6 In 32 S 21 , along with spectroscopic data, suggest that it is neither a true intermetallic nor a classical semiconductor, but somewhere in between those two extremes.

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

Photoinduced charge transfer in Zn(II) and Au(III)-ligated symmetric and asymmetric bacteriochlorin dyads: A computational study

The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso–meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4- phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi’s golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.

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

Valence Electron Distributions from Sub-Angstrom Convergent Beam Electron Diffraction

Modern aberration-corrected scanning transmission electron microscopes can acquire four-dimensional data sets (“4D STEM”) by recording convergent beam electron diffraction (CBED) patterns, using precisely positioned, sub-angstrom probes. Here, we demonstrate that these patterns can probe the site symmetry, atomic displacements, and valence electron distributions at individual atomic columns. To this end, 4D STEM CBED patterns were acquired from SrTiO 3 single crystals and compared with patterns calculated using scattering potentials derived from density functional theory. Here, we show that an aspherical valence electron charge build-up at the oxygen sites causes intensity asymmetries in the low-angle scattering portion of the patterns. Using strained SrTiO 3 films containing subtle polar displacements within nanometer-sized domains, it is shown that the high-angle scattering portion in each pattern is sensitive to atomic displacements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural Distortions and Uniaxial Negative Thermal Expansion in the Polar Dion–Jacobson Oxide RbNdTa 2 O 7

We provide deeper insight into the crystal structures, sequential structural phase transitions (I2cm → Cmce → I4/mcm → P4/mmm), thermal expansion, and electronic properties of the n = 2 Dion–Jacobson polar oxide RbNdTa 2 O 7 , through X-ray powder diffraction, neutron powder diffraction, Raman studies, and density functional theory calculations. We observed a uniaxial negative thermal expansion (NTE) across the first-order transition, I2cm → Cmce, where the unit cell contracts along the c-axis, which is driven by a contraction of the NdTa 2 O 6 layer. Here, this NTE occurs within the temperature range of the first-order phase transition and contrasts with the corkscrew mechanism typically observed in Ruddlesden–Popper phases. In RbNdTa 2 O 7 , the I2cm (hybrid improper ferroelectric) → Cmce (antipolar) transition involves crucial changes in the bond lengths of Nd and Ta polyhedra, coupled with polar to antipolar displacement of the Nd ions, leading to a net contraction in the NdTa 2 O 6 layer along the c-axis, while preserving the overall octahedral tilting magnitude. This transition highlights the intricate interplay between the Nd and Ta coordination and the associated TaO 6 distortions. Temperature-dependent Raman spectra analysis further confirms the first-order structural transition and associated NTE, providing evidence for increased bond stiffness across this transition. Additionally, using neutron powder diffraction, we have determined that the transition I4/mcm → P4/mmm occurs at approximately 1150 K. Finally, we have calculated from DFT + U, the partial density of states, the energy bandgaps, and effective masses of the charge carriers of the polar ground structure.

Chemical structure↗
DOE OSTI2021

Helical versus collinear antiferromagnetic order tuned by magnetic anisotropy in polar and chiral (Ni,Mn) 3 TeO 6

In this work, polar and chiral Ni 3 TeO 6 was found to exhibit a colossal magnetoelectric (ME) effect associated with collinear antiferromagnetic order. We have investigated the evolution of its magnetic state with the substitution of the Ni 2+ spins by isotropic Mn 2+ spins, namely (Ni, Mn) 3 TeO 6 (NMTO). The ground state of NMTO maintains an Ising nature, but a new phase with XY-type magnetic anisotropy is discovered at an intermediate temperature range. Neutron powder diffraction experiments reveal that (1) Mn ions tend to occupy a specific transition-metal site, consistent with the first-principle calculation, and (2) the intermediate phase is an incommensurate (IC) helical magnetic state propagating along the $\textit{c}$ axis with $\vec q_{\text{IC}}$ = (0, 0, 1.5 ± δ). We also found that the balance between these two magnetic states can be readily manipulated with magnetic fields, which results in significant ME effects. This selective site occupancy of NMTO allows for a unique competition between collinear and helical magnetic structures with minimal chemical disorder. Thus, NMTO serves as a model system to study site-specific chemical control of noncollinear magnetism and ME coupling.

36 MATERIALS SCIENCE↗
NASA NTRS1996

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↗