Halogen bonding: a theoretical study based on atomic multipoles derived from quantum theory of atoms in molecules

Atomic multipole moments derived from quantum theory of atoms in molecules are used to study halogen bonds in dihalogens (with general formula YX, in which X refers to the halogen directly interacted with the Lewis base) and some molecules containing C–X group. Multipole expansion is used to calculate the electrostatic potential in a vicinity of halogen atom (which is involved in halogen bonding) in terms of atomic monopole, dipole, and quadrupole moments. In all the cases, the zz component of atomic traceless quadrupole moments (where z axis taken along Y–X or C–X bonds) of the halogens plays a stabilizing role in halogen bond formation. The effects of atomic monopole and dipole moments on the formation of a halogen bond in YX molecules depend on Y and X atoms. In Br₂ and Cl₂, the monopole moment of halogens is zero and has no contribution in electrostatic potential and hence in halogen bonding, while in ClBr, FBr, and FCl it is positive and therefore stabilize the halogen bonds. On the other hand, the negative sign of dipole moment of X in all the YX molecules weakens the corresponding halogen bonds. In the C–X-containing molecules, monopole and dipole moments of X atom are negative and consequently destabilize the halogen bonds. So, in these molecules the quadrupole moment of X atom is the only electrostatic term which strengthens the halogen bonds. In addition, we found good linear correlations between halogen bonds strength and electrostatic potentials calculated from multipole expansion.     

A Hirshfeld-I interpretation of the charge distribution, dipole and quadrupole moments of the halogenated acetylenes FCCH, ClCCH, BrCCH, and ICCH

We report the dipole and quadrupole moments of the halogenated acetylenes calculated using large basis sets and the SCF, DFT(B3LYP), and CCSD methods, and we analyze the charge density using the Hirshfeld and Hirshfeld-I techniques. The atomic charges, dipoles, and quadrupoles resulting from the Hirshfeld-I analysis are used to interpret the unusually small molecular dipole moments in the sequence as well as the molecular quadrupole moments. The very small dipoles obtain for two reasons. First, the dipole moment associated with the σ and π electron densities is comparable in magnitude and opposite in direction. Second, the charge and induced dipole contributions for ClCCH, BrCCH, and ICCH have opposite signs further reducing the molecular dipoles. The molecular quadrupole moments are the sum of a charge, atomic dipole, and in situ quadrupole terms, and are dominated by the atomic dipoles and in situ quadrupoles with the charge contributions playing an unexpectedly minor role.     

Molecular dynamics study of polarizable point dipole models for molten sodium iodide

The structure, the ionic transport properties, and the dynamics of long-wavelength charge-density fluctuations, for two polarizable point dipole models of molten NaI, have been studied by molecular dynamics simulations. These models are based on a rigid ion potential to which the induced dipole polarization of the anions is added. The polarization is added in such a way that point dipoles are induced on the anions by both local electric field and short-range damping interactions that oppose the electrically induced dipole moments. The two polarizable ion models differ only in the range of the damping polarization interactions. The influence of the induced anion polarization on the different properties of simulated molten NaI is discussed.     

Electrophoretic motion of a slightly deformed sphere with a nonuniform zeta potential distribution

Electrophoretic motion is analyzed for a rigid, slightly deformed sphere with a nonuniform zeta potential distribution. Hydrodynamics and electrostatics solutions for the deformed sphere with an arbitrary double-layer thickness are determined by using the domain perturbation method. The surface shape and the zeta potential distribution for the deformed sphere are expressed by using the multipole expansion representation. In terms of monopole, dipole, and quadrupole moments of the surface shape and the zeta potential distribution, explicit expressions are obtained for the translational and rotational electrophoretic mobility tensors. The ensemble average for the mobility of the deformed sphere with a uniform orientation distribution is also derived. The utility of the general mobility expression is demonstrated by studying the electrophoretic motion of axisymmetric and ellipsoidal particles. The translational and rotational mobilities of axisymmetric particles are both affected by the monopole, dipole, and quadrupole moments of the zeta potential. For ellipsoidal particles, however, the dipole moment of the zeta potential does not affect the translational mobility, while the rotational mobility depends only on the dipole moment. The mobility of the deformed sphere with either a thick or a thin double layer is also derived.     

Electric parameters of photosystem II particles — electric light scattering study

Electric light scattering and microelectrophoresis were applied to investigate the electric moments (permanent dipole moment and electric polarizability and electrophoretic mobility of envelope-free chloroplasts and photosystem II (PS II particles. The effect of the removal of the extrinsic polypeptides (18, 24 and 33 kDa) on the electric moments was also studied. A significant difference was observed between the orientation behaviour of chloroplasts and PS II preparations. The data indicate that the permanent and induced dipole moments contribute to the orientation of the PS II particles, whereas chloroplasts possess induced dipole moment only.

NaCl and Tris treatments of PS II preparations influence both the transverse permanent dipole moment and the electric polarizability of PS II particles. The increase in the electrophoretic mobility of PS II particles on removal of the extrinsic proteins corresponds to an increase in the electric polarizability value, demonstrating its interfacial nature.     

Structure and dynamics of the aqueous liquid-vapor interface: a comprehensive particle-based simulation study

This research addresses a comprehensive particle-based simulation study of the structural, dynamic, and electronic properties of the liquid-vapor interface of water utilizing both ab initio (based on density functional theory) and empirical (fixed charge and polarizable) models. Numerous properties such as interfacial width, hydrogen bond populations, dipole moments, and correlation times will be characterized with identical schemes to draw useful conclusions on the strengths and weakness of the proposed models for interfacial water. Our findings indicate that all models considered in this study yield similar results for the radial distribution functions, hydrogen bond populations, and orientational relaxation times. Significant differences in the models appear when examining both the dipole moments and surface relaxation near the aqueous liquid-vapor interface. Here, the ab initio interaction potential predicts a significant decrease in the molecular dipole moment and expansion in the oxygen-oxygen distance as one approaches the interface in accordance with recent experiments. All classical polarizable interaction potentials show a less dramatic drop in the molecular dipole moment, and all empirical interaction potentials studied yield an oxygen-oxygen contraction as the interface is approached.     

A polarizable model for ethylene oxide

A series of interaction models for ethylene oxide are developed for use in molecular simulation of the thermal properties of both the gas and liquid phases. While it is possible to develop nonpolarizable models to accurately generate either the gas or liquid properties separately, it was not possible to do so using a single model for both phases. A polarizable, rigid all-atom model was developed that reproduces the temperature dependence of the second virial coefficient B(T) and the pressure of the liquid at ambient conditions. The model consists of Lennard-Jones and Coulomb interactions between intermolecular atomic sites plus a scalar polarizability located at the midpoint of the line joining the carbon sites. The electrostatic charges and the polarizability are set to match the experimentally determined dipole and quadrupole moments and the molecular polarizability.     

On the performance of molecular polarization methods. I. Water and carbon tetrachloride close to a point charge

The three main methods to implement molecular polarization (point dipoles, fluctuating charges, and shell model) are tested against high level ab initio calculations for a molecule (water, carbon tetrachloride) close to a point charge (at the distance of a lithium or magnesium ion). The goal is to check whether an approximation (linear polarization) strictly valid at large intermolecular distances is sufficiently accurate for liquid state molecular dynamics simulations, where strong polarization effects are to be expected at short separations. The monitored observable is the molecular dipole moment as a function of the charge-molecule distance for selected molecular orientations. Analytic formulas are derived for the components of the molecular polarization tensor, facilitating the optimization of the performance for each polarization method as a function of its underlying parameters. Overall, the methods studied provide a remarkably good representation of the induced dipole, with no divergences appearing even at the shortest distances. For water close to a monovalent point charge the point dipole model, implemented with one or three dipoles, accurately reproduces the water dipole moment at all distances. Deficiencies appear as the molecular polarizability and/or charge increase: basically, the ab initio induced moments grow faster at intermediate distances than the linear increase characteristic of the phenomenological polarization methods, suggesting that nonlinear effects (hyperpolarizability) cannot be neglected in these cases. Regarding the capabilities of each method, the point dipole method is the one that performs best overall, with the shell model achieving acceptable results in most instances. The fluctuating charge method shows some noticeable limitations for implementations of comparable complexity (in terms of the number of sites required).     

Polarizable contributions to the surface tension of liquid water

Surface tension, gamma, strongly affects interfacial properties in fluids. The degree to which polarizability affects gamma in water is thus far not well established. To address this situation, we carry out molecular dynamics simulations to study the interfacial forces acting on a slab of liquid water surrounded by vacuum using the Gaussian charge polarizable (GCP) model at 298.15 K. The GCP model incorporates both a fixed dipole due to Gaussian distributed charges and a polarizable dipole. We find a well-defined bulklike region forms with a width of approximately 31 A. The average density of the bulklike region agrees with the experimental value of 0.997 g/cm3. However, we find that the orientation of the molecules in the bulklike region is strongly influenced by the interfaces, even at a distance five molecular diameters from the interface. Specifically, the orientations of both the permanent and induced dipoles show a preferred orientation parallel to the interface. Near the interface, the preferred orientation of the dipoles becomes more pronounced and the average magnitude of the induced dipoles decreases monotonically. To quantify the degree to which molecular orientation affects gamma, we calculate the contributions to gamma from permanent dipolar interactions, induced dipolar interactions, and dispersion forces. We find that the induced dipole interactions and the permanent dipole interactions, as well as the cross interactions, have positive contributions to gamma, and therefore contribute stability to the interface. The repulsive core interactions result in a negative contribution to gamma, which nearly cancels the positive contributions from the dipoles. The large negative core contributions to gamma are the result of small oxygen-oxygen separation between molecules. These small separations occur due to the strong attractions between hydrogen and oxygen atoms. The final predicted value for gamma (68.65 m/Nm) shows a deviation of approximately 4% of the experimental value of 71.972 m/Nm. The inclusion of polarization is critical for this model to produce an accurate value.     

A theoretical study of the hydration of Rb+ by Monte Carlo simulations with refined ab initio-based model potentials

An infinitely diluted aqueous solution of Rb⁺ was studied using ab initio-based model potentials in classical Monte Carlo simulations to describe its structural and thermodynamic features. An existing flexible and polarizable model [Saint-Martin et al. in J Chem Phys 113(24) 10899, 2000] was used for water–water interactions, and the parameters of the Rb⁺–water potential were fitted to reproduce the polarizability of the cation and a sample of ab initio pair interaction energies. It was necessary to calibrate the basis set to be employed as a reference, which resulted in a new determination of the complete basis set (CBS) limit energy of the optimal Rb⁺–OH₂ configuration. Good agreement was found for the values produced by the model with ab initio calculations of three- and four-body nonadditive contributions to the energy, as well as with ab initio and experimental data for the energies, the enthalpies and the geometric parameters of Rb⁺(H₂O) _n clusters, with n = 1,  2,…, 8. Thus validated, the potential was used for simulations of the aqueous solution with three versions of the MCDHO water model; this allowed to assess the relative importance of including flexibility and polarizability in the molecular model. In agreement with experimental data, the Rb⁺–O radial distribution function (RDF) showed three maxima, and hence three hydration shells. The average coordination number was found to be 6.9, with a broad distribution from 4 to 12. The dipole moment of the water molecules in the first hydration shell was tilted to 55° with respect to the ion’s electric field and had a lower value than the average in bulk water; this latter value was recovered at the second shell. The use of the nonpolarizable version of the MCDHO water model resulted in an enhanced alignment to the ion’s electric field, not only in the first, but also in the second hydration shell. The hydration enthalpy was determined from the numerical simulation, taking into account corrections to the interfacial potential and to the spurious effects due to the periodicity imposed by the Ewald sums; the resulting value lied within the range of the various different experimental data. An analysis of the interaction energies between the ion and the water molecules in the different hydration shells and the bulk showed the same partition of the hydration enthalpy as for K⁺. The reason for this similarity is that at distances longer than 3 Å, the ion–water interaction is dominated by the charge-(enhanced) dipole term. Thus, it was concluded that starting at K⁺, the hydration properties of the heavier alkali metal cations should be very similar.     

Molecular electric properties in electronic excited states: multipole moments and polarizabilities of H2O in the lowest1 B 1 and3 B 1 excited states

Summary The dipole and quadrupole moments and the dipole polarizability tensor components are calculated for the¹ B ₁ and³ B ₁ excited states of the water molecule by using the complete active space (CAS) SCF method and an extended basis set of atomic natural orbitals. The dipole moment in the lowest¹ B ₁ (0.640 a.u.) and³ B ₁ (0.416 a.u.) states is found to be antiparallel to that in the ground electronic state of H₂O. The shape of the quadrupole moment ellipsoid is significantly modified by the electronic excitation to both states investigated in this paper. All components of the excited state dipole polarizability tensor increase by about an order of magnitude compared to their values in the ground electronic state. The present results are used to discuss some aspects of intermolecular interactions involving molecules in their excited electronic states.     

Structure, thermodynamics, and liquid-vapor equilibrium of ethanol from molecular-dynamics simulations using nonadditive interactions

We present a molecular-dynamics simulation study of the bulk and liquid-vapor interfacial properties of ethanol using a polarizable force field based on the fluctuating charge (FQ) formalism, as well as the nonpolarizable CHARMM22 force field. Both models are competitive with respect to the prediction of ambient liquid properties such as liquid density, enthalpy of vaporization, dielectric constant, and self-diffusion constants. The polarizable model predicts an average condensed-phase dipole moment of 2.2 D associated with an induced liquid-phase dipole moment of 0.6 D; though qualitatively in agreement with earlier nonadditive models as well as recent Car-Parinello calculations, the current FQ model underestimates the condensed-phase dipole moment. In terms of liquid structure, both models are in agreement with recent neutron-diffraction results of liquid ethanol structure, although the polarizable model predicts the hydroxyl-hydrogen-hydroxyl-hydrogen structure factor in closer agreement with the experimental data. In terms of interfacial properties, both models predict ambient surface tension to within 4% of the experimental value of 22.8 dyncm, while overestimating the surface excess entropy by almost a factor of 2. Both models display the characteristic preferential orientation of interfacial molecules. The polarizable model allows for a monotonic variation of the average molecular dipole moment from the bulk value to that of the vapor phase. Consequently, there is a dramatic difference in the surface potential predicted by the polarizable and nonpolarizable models. The polarizable model estimates a surface potential of -209+/-3 mV, while the nonpolarizable model yields a value of -944+/-10 mV. Finally, based on the vapor-liquid equilibrium simulation data from several temperatures, we estimate the critical properties of both models. As observed with other FQ models for associating fluids (such as water and methanol), and counter to what one would anticipate by modeling more physically the electrostatic response to local environment, the current FQ model underestimates the critical temperature and overestimates the critical density of ethanol; moreover, the FQ model is, in this respect, equivalent to the underlying fixed-charge model. These results further suggest the need to revisit polarizable models in terms of quantitative vapor-liquid equilibrium prediction.     

Studies on the electrical double layer structure at the water/air interface

Effective dipole moments (calculated from experimental data of surface tension and electric surface potential) of some homologous normal alcohols and carboxylic acid were found to vary linearly with the number of carbon atoms in the hydrocarbon chain. Values of effective dipole moments were used for the determination of the effective dipole moments of water molecules , and the dielectric permittivity of the water subphase (₁), as well as in the vicinity of the hydrophobic part of adsorbed molecule (₂). The latter was found to decrease with the increase of the hydrocarbon chain length. Knowing the effective dipole moment of surface water dipoles, the average orientation angle () of water molecules at the inteface was estimated. The calculated potential drop of water varies within the range –0.038 to –2.38 V for two extreme orientations of water dipoles at the surface.     

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.     

Orientational ordering of a polymer chain in a strong dipole field. Tetrahedral lattice model

Orientational and conformational properties have been investigated of a rotational isomeric model of a polar polymer chain on a tetrahedral lattice in a strong electric or mechanical field of dipole symmetry. Two types of dipole moment distribution along the chain are discussed: (A) constant signs of longitudinal components of the dipole moments, and (B) alternating signs of the longitudinal components of the dipole moments. The second case represents polymer chains such as $ \rlap{--} ({\rm CH}_{\rm 2} \hbox{---} {\rm CR}_{\rm 2} {\rm \rlap{--} )}_n $>/UEQN> when the dipole moments are oriented along the bisector of the CR₂-angle, i.e., normal to the extended trans-chain conformation. It is shown that only a discrete most probable orientation of the lattice relative to the field should be considered, namely that coinciding with one of the symmetry axes of the lattices. The average dipole order parameter and dichroic functions (quadrupole order parameter) are calculated in a strong external dipole field for unit vectors with different orientation relative to the chain backbone. The quadrupole order parameter for different unit vectors is obtained also as a function of chain elongation. The polarizability induced by an additional weak dipole field is calculated as a function of the magnitude of the strong external dipole field. For the model considered here the order parameters are more strongly influenced by the external field than those for the freely jointed chain (FJC) model having the same distribution of dipole moments along the chain. The orientational ordering of the chain in a dipole field is higher than in a quadrupole field of the same magnitude.     

Dipole preserving and polarization consistent charges

A method for estimating dipole preserving and polarization consistent (DPPC) charges is described, which reproduces exactly the molecular dipole moment as well as the local, atomic hybridization dipoles determined from the corresponding wave function and can yield accurate molecular polarization. The method is based on a model described by Thole and van Duijnen and a new feature is introduced to treat molecular polarization. Thus, the DPPC method offers a convenient procedure to describe molecular polarization in applications using semiempirical models and ab initio molecular orbital theory with relatively small basis functions such as 6‐31+G(d,p) or without inclusion of electron correlation; these methods tend to underestimate molecular polarizability. The trends of the DPPC partial atomic charges are found to be in good accord with those of the CM2 model, a class IV charge analysis method that has been used in a variety of applications. The DPPC method is illustrated to mimic the correct molecular polarizability in a water dimer test case and in water‐halide ion complexes using the explicit polarization (X‐Pol) potential with the Austin model 1 Hamiltonian. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Study of multipole contributions to the structure of water around ions in solution using the soft sticky dipole-quadrupole-octupole (SSDQO) model of water

The solvation of ions in the soft sticky dipole-quadrupole-octupole (SSDQO) model for liquid water is presented here. This new potential energy function for liquid water describes water-water interactions by a Lennard-Jones term plus a sticky potential consisting of an approximate moment expansion with point dipole, quadrupole, and octupole moments. The SSDQO potential energy function using the moments from extended simple point charge (SPC/E), TIP3P, or TIP5P reproduces the pair potential energy functions and radial distribution functions of the respective multipoint model but it is much faster than even the three-point models. Here, the solvation of ions in SSDQO water is studied using ion-water potential energy functions consisting of moment expansions up to the charge-quadrupole term, up to the charge-octupole term, and up to an approximate charge-hexadecapole term using the moments of SPC/E water. The radial distributions from Monte Carlo simulations show the best agreement with the results for ions in SPC/E water for the expansion up to the charge-hexadecapole term. Thus, the best results are obtained when the water-water and ion-water potentials are exact up to the 1r(4) term and also contain an approximate 1r(5) term. Overall, the simplicity, efficiency, and accuracy of the SSDQO potential make it potentially very useful for computer simulations of aqueous solvation.     

Water structure in aqueous solutions of alkali halide salts: FTIR spectroscopy of the OD stretching band

Water structure making/breaking studies in solutions of five alkali halide salts (KF, KI, NaI, CsF and CsCl) in 4 wt% D(2)O in H(2)O mixtures have been performed by FTIR analysis of the OD stretching band in the full solubility range. The proposed method gives a microscopic picture of the water structure making/breaking character of the salts in terms of the hydrogen bonding between the water molecules in the solution. With the exception of CsCl, there is a very good correlation of the structure making/breaking character of the salts determined by FTIR analysis, and the viscosity coefficients of the solutions. The results fully support and explain previous studies of bubble attachment to microscopic salt particles of the above salts. The investigations support the primary importance of interfacial water structure in the explanation of the flotation of alkali halide salts in their brines.