Halide anion capture and recognition by a tetrahedral tetraammonium receptor in water: a molecular dynamics investigation

We report molecular dynamics potential of mean force (PMF) simulations on the capture of halide anions X(-) (F(-), Cl(-), Br(-)) by a tetrahedral receptor L(4+) built from four quaternary ammonium sites connected by six (CH(2))(n) chains, leading to the formation of inclusion complexes X(-) subset L(4+). Simulations performed with a reaction field correction of the electrostatics and with PME-Ewald summation gave very similar energy profiles. In aqueous solution, an energy barrier of 12-17 kcal mol(-1) was found for the three anions, mainly due to their dehydration when they enter through the largest triangular face of L(4+). In the inclusion complexes, the anion is anchored near the center of the cavity due to the electrostatic field of the four positively charged ammonium sites, shielded from the surrounding water molecules. It was predicted that L(4+) is selective for Cl(-) over Br(-) which both form stable inclusion complexes, while the F(-) complex should dissociate. The comparison of PMFs in aqueous solution and in the gas phase and the energy component analysis demonstrates the importance of solvent on the nature of these complexes and on the complexation energy profiles. The Cl(-)/Br(-) selectivity obtained from the dissociation pathways in water was in good agreement with the results of free energy perturbation simulations based on the "alchemical route" of a thermodynamic cycle, and consistent with experimental observations.     

Halide anion recognition in water by a hexaprotonated octaaza-cryptand: a molecular dynamics investigation

Based on molecular dynamics simulations, we describe the F- versus Cl- complexation by an hexaprotonated cryptand L6+ in aqueous solution, in order to elucidate their structures, solvation properties and the status of external halide counterions. In water, F- and Cl- simulated inclusive complexes adopt a structure somewhat different from the solid state structure of the F- complex: The anion binding involves two diammonium bridges only, and the accompanying counterions are dissociated from the +5 charged complex. A remarkable result is obtained for the dissociated L6+,3F- ,3Cl- system, where spontaneous complexation of F- (the anion which forms the most stable complex with L6+) takes place during the dynamics. The resulting complex is of facial type; this suggests that the equilibrium involves multiple binding modes and structures in aqueous solution. The question of F-/Cl- binding selectivity is investigated by free energy perturbations simulations which nicely reproduce the spectacular preference for F- over Cl-. Two different methodologies used for the treatment of electrostatics (standard versus Ewald calculations) yield similar conclusions.     

Molecular dynamics analysis of interfacial structures and sum frequency generation spectra of aqueous hydrogen halide solutions

Molecular dynamics simulation for gas/liquid interfaces of aqueous hydrochloric (HCl) and hydroiodic (HI) acid solutions is performed to calculate and analyze their sum frequency generation (SFG) spectra. The present MD simulation supports the strong preference of hydronium ions at the topmost surface layer and a consequent formation of ionic double layers by the hydronium and halide ions near the interface. Accordingly, the orientational order of surface water in the double layers is reversed in the acid solutions from that in the salt (NaCl or NaI). The calculated SFG spectra of the O-H stretching region reproduce the experimental spectra of ssp and sps polarizations well. In the ssp spectra, the strong enhancement in the hydrogen-bonding region for the acid solutions is elucidated by two mechanisms, ordered orientation of water in the double layer and symmetric OH stretching of the surface hydronium ions. In the sps spectra, reversed orientation of surface water is evidenced in the spectral line shapes, which are quite different from those of the salt solutions.     

The diamond/aqueous electrolyte interface: an impedance investigation

We have investigated the electrochemical interface between diamond electrodes and aqueous electrolytes using electrochemical techniques such as cyclic voltammetry and ac impedance spectroscopy. High-quality CVD-grown boron-doped polycrystalline diamond electrodes and IIa single crystalline natural diamond electrodes have been used in this study. In the case of hydrogen-terminated diamond electrodes, the electrochemical interface is dominated by the electrochemical double layer. Frequency-dependent impedance spectroscopy reveals a potential regime in which the contribution of ion adsorption becomes relevant. We have conducted experiments to evaluate the effect of pH and ionic strength on the double layer. Our results suggest that only ions resulting from water auto-dissociation, i.e., hydroxide and hydronium ions, are responsible for ion adsorption and, thus, able to modify the charge at the double layer. In contrast, no effect of the adsorption of several dissolved ions (such as Na+, K+, Cl-) has been observed On the basis of the electrochemical characterization of H-terminated diamond surfaces, we also discuss the phenomenon of the surface conductivity in diamond, as well as the pH sensitivity of the diamond surface. The influence of the O2/OH- and H2/H3O+ redox couples on the origin of the surface conductivity is discussed.     

A molecular dynamics investigation of the titration of a trivalent aqueous ion

We carried out a series of molecular dynamics simulations of the hydrolysis of a model trivalent metal ion in aqueous solution. We use a dissociative model for water and examine the spontaneous speciation of M³⁺ into M(OH) _n ^(3-n)+ (n =1,4) both in neutral solution and as a function of added protons and hydroxide ions. The species distributions in neutral solution correspond reasonably well with those expected for real trivalent metal ions at neutral pH. However, the change in the species distributions as a function of either added protons or hydroxide ions is much less than expected with very large concentrations of protons or hydroxide ions required to shift the species equilibria in either direction. The influence of added protons and hydroxide ions on the species distributions appears to be proportional to the average charge of the hydrolysis couples, being highest for the 3+/2+ couple and lowest for the 1+/0 and 0/1- couples. Proton exchange rates vary with proton/hydroxide ion concentration giving a minimum at intermediate values ([H+]≈ 0.166) with increasing rates at both lower and higher pH.     

Fluorescent charge-assisted halogen-bonding macrocyclic halo-imidazolium receptors for anion recognition and sensing in aqueous media

The synthesis and anion binding properties of a new family of fluorescent halogen bonding (XB) macrocyclic halo-imidazolium receptors are described. The receptors contain chloro-, bromo-, and iodo-imidazolium motifs incorporated into a cyclic structure using naphthalene spacer groups. The large size of the iodine atom substituents resulted in the isolation of anti and syn conformers of the iodo-imidazoliophane, whereas the chloro- and bromo-imidazoliophane analogues exhibit solution dynamic conformational behavior. The syn iodo-imidazoliophane isomer forms novel dimeric isostructural XB complexes of 2:2 stoichiometry with bromide and iodide anions in the solid state. Solution phase DOSY NMR experiments indicate iodide recognition takes place via cooperative convergent XB-iodide 1:1 stoichiometric binding in aqueous solvent mixtures. (1)H NMR and fluorescence spectroscopic titration experiments with a variety of anions in the competitive CD(3)OD/D(2)O (9:1) aqueous solvent mixture demonstrated the bromo- and syn iodo-imidazoliophane XB receptors to bind selectively iodide and bromide respectively, and sense these halide anions exclusively via a fluorescence response. The protic-, chloro-, and anti iodo-imidazoliophane receptors proved to be ineffectual anion complexants in this aqueous methanolic solvent mixture. Computational DFT and molecular dynamics simulations corroborate the experimental observations that bromo- and syn iodo-imidazoliophane XB receptors form stable cooperative convergent XB associations with bromide and iodide.     

A simple operational relationship between drop-time and interfacial tension at a charged mercury/aqueous solution interface

This paper describes a simple operational relationship between the drop-time of a dropping mercury electrode and the interfacial tension at a charged mercury/aqueous solution interface. An apparatus and technique for measuring drop-times is reported, and examples are given of the fit of experimental drop-times to literature values of interfacial tensions. The operational relationship is independent of temperature in the range 293-313 K.     

A molecular dynamics study of sodium chenodeoxycholate in an aqueous solution

Hydration structure and dynamics of sodium chenodeoxycholate (CDC) in water are studied by a long-time molecular dynamics calculation. Strong hydration shell around the hydrophobic region of this large solute and strong hydrogen bonds of water with both hydroxyl and carboxyl oxygen atoms have been identified. The rotation of CDC around its longitudinal axis is found to be particularly active in comparison with that around other axes of the molecule. The diffusion coefficient of CDC calculated from the slope of the mean-square displacement, 0.95 × 10⁻⁹ m²/s, is only 1/6 of that for water in the solution, 5.4 × 10⁻⁹ m²/s.     

Structure and Energetics of model amphiphilic molecules at the water liquid-vapor interface. A molecular dynamics study

A molecular dynamics study of adsorption of p-n-pentylphenol at infinite dilution at the water liquid-vapor interface is reported. The calculated free energy of adsorption is -8.8 +/- 0.7 kcal/mol, in good agreement with the experimental value of -7.3 kcal/mol. The transition between the interfacial region and the bulk solution is sharp and well-defined by energetic, conformational, and orientational criteria. At the water surface, the phenol head group is mostly immersed in aqueous solvent. The most frequent orientation of the hydrocarbon tail is parallel to the interface, due to dispersion interactions with the water surface. This arrangement of the phenol ring and the alkyl chain requires that the chain exhibits a kink. As the polar head group is being moved into the solvent, the chain length increases and the tail becomes increasingly aligned toward the surface normal, such that the nonpolar part of the molecule exposed to water is minimized. The same effect was achieved when phenol was replaced by a more polar head group, phenolate. This result underscores the difference between hydrophobic hydration at the surface and in the bulk solvent, when nonpolar molecular fragments adopt compact conformations.     

Effect of bromide on the interfacial structure of aqueous tetrabutylammonium iodide: Photoelectron spectroscopy and molecular dynamics simulations

Solvation of surface-active tetrabutylammonium iodide (TBAI) in liquid water and in sodium bromide aqueous solution was investigated by VUV photoelectron spectroscopy and by molecular dynamics simulations. The observed signal intensity changes in the photoemission spectra are consistent with the varying propensities of the different ions for the solution interface. While the cations are surface-bound due to hydrophobic interactions, the anions are driven to the vacuum/solution interface by their large polarizability and size. Iodide is more polarizable, and hence more surface-active than bromide, which explains the relatively small decrease of the iodide photoemission signal when TBAI is dissolved in bromide solution.     

Interior and interfacial aqueous solvation of benzene dicarboxylate dianions and their methylated analogues: A combined molecular dynamics and photoelectron spectroscopy study

Aqueous solvation of benzene dicarboxylate dianions (BCD(2-)) was studied by means of photoelectron spectroscopy and molecular dynamics simulations. Photoelectron spectra of hydrated o- and p-BCD(2-) with up to 25 water molecules were obtained. An even-odd effect was observed for the p-BCD(2-) system as a result of the alternate solvation of the two negative charges. However, the high polarizability of the benzene ring makes the two carboxylate groups interact with each other in p-BCD(2-), suppressing the strength of this even-odd effect compared with the linear dicarboxylate dianions linked by an aliphatic chain. No even-odd effect was observed for the o-BCD(2-) system, because each solvent molecule can interact with the two carboxylate groups at the same time due to their proximity. For large solvated clusters, the spectral features of the solute decreased while the solvent features became dominant, suggesting that both o- and p-BCD(2-) are situated in the center of the solvated clusters. Molecular dynamics simulations with both nonpolarizable and polarizable force fields confirmed that all three isomers (o-, m-, and p-BCD(2-)) solvate in the aqueous bulk. However, upon methylation the hydrophobic forces overwhelm electrostatic interactions and, as a result, the calculations predict that the tetramethyl-o-BCD(2-) is located at the water surface with the carboxylate groups anchored in the liquid and the methylated benzene ring tilted away from the aqueous phase.     

Phase behavior of amphiphiles at liquid crystals/water interface: A coarse-grained molecular dynamics study

We present a coarse-grained molecular dynamics simulation study of phase behavior of amphiphilic monolayers at the liquid crystal （LC）/water interface. The results revealed that LCs at interface can influence the lateral ordering of amphiphiles. Particularly, the amphiphile tails along with perpendicularly penetrated LCs between tails undergo a two-dimension phase transition from liquid-expanded into a liquid-condensed phase as their area density at interface reaches 0.93. While, the liquid-condensed phase of the monolayer never appears at oil/water interface with isotropic shape oil particles. These findings reveal the penetration of anisotropic LC can promote ordered lateral organization of amphiphiles. Moreover, we find the phase transition point is shifted to lower surface coverage of amphiphiles when the LCs have larger affinity to the amphiphile tails.     

Conformational equilibria of terminally blocked single amino acids at the water-hexane interface. A molecular dynamics study

The conformational equilibria of the acetyl and methyl amide terminally blocked L-alanine, L-leucine and L-glutamine amino acids are examined in vacuum, in bulk water, and at the water-hexane interface, using multi-nanosecond molecular dynamics simulations. The two-dimensional probability distribution functions of finding the peptides at different dihedral angles of the backbone, phi and psi, are calculated, and free energy differences between different conformational states are determined. All three peptides are interfacially active, i.e. tend to accumulate at the interface even though they are not amphiphilic. Conformational states stable in both gas phase and water are also stable in the interfacial environment. Their populations, however, cannot be simply predicted from the knowledge of conformational equilibria in the bulk phases, indicating that the interface exerts a unique effect on the peptides. Conformational preferences in the interfacial environment arise from the interplay between electrostatic and hydrophobic effects. As in an aqueous solution, electrostatic solute-solvent interactions lead to the stabilization of more polar peptide conformations. The hydrophobic effect is manifested at the interface by a tendency to segregate polar and nonpolar moieties of the solute into the aqueous and the hexane phases, respectively. For the terminally blocked glutamine, this favors conformations for which such a segregation is compatible with the formation of strong, backbone-side chain intramolecular hydrogen bonds on the hexane side of the interface. The influence of the hydrophobic effect can be also noted in the orientational preferences of the peptides at the interface. The terminally blocked leucine is oriented such that its nonpolar side chain is buried in hexane, whereas the polar side chain of glutamine is immersed in water. The free energies of rotating the peptides along the axis parallel to the interface by more than 90 degrees are substantial. This indicates that peptide folding at interfaces is strong by driven by the tendency to adopt amphiphilic structures.     

Layering at an ionic liquid-vapor interface: a molecular dynamics simulation study of [bmim][PF6

The structure of the planar liquid-vapor interface of a room-temperature ionic liquid, 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6]), is studied using atomistic molecular dynamics simulations. Layering of the ions at the interface is observed as oscillations in the corresponding number density profiles. These oscillations, however, are diminished in amplitude in the electron density profile, due to a near cancellation in the contributions from the anions and the cations. An enhancement by 12% in the electron density at the interface over its value in the bulk liquid is observed, in excellent agreement with X-ray reflectivity experiments. The anions are found to predominantly contribute to this increase in the interfacial electron density. The cations present at the interface are oriented anisotropically. Their butyl chains are observed to be preferentially oriented along the interface normal and to project outside the liquid surface, thus imparting a hydrophobic character. In the densest region of the interface, the imidazolium ring plane is found to lie parallel to the surface normal, in agreement with direct recoil spectroscopy experiments.     

A theoretical investigation on structures of tripodal thiourea derivatives and their anion recognition

The structural geometries of three tripodal thiourea receptors, i.e. 1,3,5-triethyl-2,4,6-tris[(N′-methylthioureido)methyl]benzene (1), tris[N′-methyl-N-(2-aminoethyl)thiourea]methane (2), tris[N′-methyl-N-(2-aminoethyl)thiourea]amine (3), and their complexes with F⁻, Cl⁻, Br⁻, I⁻, NO₃ ⁻, CO₃ ²⁻, SO₄ ²⁻, HSO₄ ⁻, PO₄ ³⁻, HPO₄ ²⁻ and H₂PO₄ ⁻ were obtained using the density functional theory calculations. Electronic and thermodynamic properties of anion binding complexes of the receptors 1, 2 and 3 were investigated. Recognition abilities of all the receptors in terms of selectivity coefficients are reported. Intermolecular interactions in all the studied complexes occurring via multi-point hydrogen bonding were found. The receptors 1, 2 and 3 were found to be excellent selectivity for phosphate ion and their binding free energy for the phosphate ion are −292.57, −291.77 and −295.01 kcal/mol, respectively.     

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.     

Molecular recognition of sulfonatocalixarene with organic cations at the self-assembled interface: a thermodynamic investigation

A microcalorimetric study on molecular recognition of p-sulfonatocalix[4]arene derivatives at selfassembled interface in comparison with in bulk water was performed,inspired by the dramatic change in physicochemical characteristics from bulk water to interface.A total of six cationic molecules were screened as model guests,including ammonium(NH_4~+),guanidinium(Gdm~+).N,N'-dimethyl-1,4-diazabicyclo[2.2.2]octane(DMDABCO~(2+)),tropylium(Tpm~+),N-methyl pyridinium(N-mPY*) and methyl viologen(MV~(2+)).The complexation with NH_4~+.Gdm~+ and DMDABCO2* is pronouncedly enhanced when the recognition process moved from bulk water to interface,whereas the complexation stabilities with Tpm~+,N-mPY~+ and MV2* increase slightly or even decrease to some extent.A more interesting phenomenon arises from the NH_4~+/Gdm~+ pair that the thermodynamic origin at interface differs definitely from each other although with similar association constants.The results were discussed in terms of differential driving forces,electrostatic,hydrogen bond as well as π-stacking interactions,originating from the unique physicochemical features of interfaces,mainly the polarity and dielectric constant.