Dynamical arrest of electron transfer in liquid crystalline solvents

We argue that electron transfer reactions in slowly relaxing solvents proceed in the nonergodic regime, making the reaction activation barrier strongly dependent on the solvent dynamics. For typical dielectric relaxation times of polar nematics, electron transfer reactions in the subnanosecond time scale fall into nonergodic regime in which nuclear solvation energies entering the activation barrier are significantly lower than their thermodynamic values. The transition from isotropic to nematic phase results in weak discontinuities of the solvation energies at the transition point and the appearance of solvation anisotropy weakening with increasing solute size. The theory is applied to analyze experimental kinetic data for the electron transfer kinetics in the isotropic phase of 5CB liquid crystalline solvent. We predict that the energy gap law of electron transfer reactions in slowly relaxing solvents is characterized by regions of fast change of the rate at points where the reaction switches between the ergodic and nonergodic regimes. The dependence of the rate on the donor-acceptor separation may also be affected in a way of producing low values for the exponential falloff parameter.     

Electronic absorption spectra of cresyl violet acetate in anisotropic and isotropic solvents

The isotropic and anisotropic solvation characteristics of cresyl violet acetate (CVA) were investigated in isotropic liquid solutions and in polar nematic matrices as a function of the solvent type and concentration. The interaction of the ionic dye with the anisotropic surrounding and with that of the isotropic solvents was investigated and compared. The experimental result suggests that the nematic liquid crystalline solvents might create stronger solvation than the isotropic solvents. The spectral shifts were correlated by the solvent permittivity and Kamlet-Taft parameters. The polarized absorption spectra of cresyl violet acetate were measured between 400 and 800 nm and the dichroic ratio R and degree of anisotropy S of this dye in the liquid crystalline host determined.     

Density functional theory of solvation and its relation to implicit solvent models

We describe a density functional theory approach to solvation in molecular solvents. The solvation free energy of a complex solute can be obtained by direct minimization of a density functional, instead of the thermodynamic integration scheme necessary when using atomistic simulations. In the homogeneous reference fluid approximation, the expression of the free-energy functional relies on the knowledge of the direct correlation function of the pure solvent. After discussing general molecular solvents, we present a generic density functional describing a dipolar solvent and we show how it can be reduced to the conventional implicit solvent models when the solvent microscopic structure is neglected. With respect to those models, the functional includes additional effects such as the microscopic structure of the solvent, the dipolar saturation effect, and the nonlocal character of the dielectric constant. We also show how this functional can be minimized numerically on a three-dimensional grid around a solute of complex shape to provide, in a single shot, both the average solvent structure and the absolute solvation free energy.     

Free energy of solvation from molecular dynamics simulation applying Voronoi-Delaunay triangulation to the cavity creation

The free energy of solvation for a large number of representative solutes in various solvents has been calculated from the polarizable continuum model coupled to molecular dynamics computer simulation. A new algorithm based on the Voronoi-Delaunay triangulation of atom-atom contact points between the solute and the solvent molecules is presented for the estimation of the solvent-accessible surface surrounding the solute. The volume of the inscribed cavity is used to rescale the cavitational contribution to the solvation free energy for each atom of the solute atom within scaled particle theory. The computation of the electrostatic free energy of solvation is performed using the Voronoi-Delaunay surface around the solute as the boundary for the polarizable continuum model. Additional short-range contributions to the solvation free energy are included directly from the solute-solvent force field for the van der Waals-type interactions. Calculated solvation free energies for neutral molecules dissolved in benzene, water, CCl4, and octanol are compared with experimental data. We found an excellent correlation between the experimental and computed free energies of solvation for all the solvents. In addition, the employed algorithm for the cavity creation by Voronoi-Delaunay triangulation is compared with the GEPOL algorithm and is shown to predict more accurate free energies of solvation, especially in solvents composed by molecules with nonspherical molecular shapes.     

Transient cavities and the excess chemical potentials of hard-spheroid solutes in dipolar hard-sphere solvents

Monte Carlo computer simulations are used to study transient cavities and the solvation of hard-spheroid solutes in dipolar hard-sphere solvents. The probability distribution of spheroidal cavities in the solvent is shown to be well described by a Gaussian function, and the variations of fit parameters with cavity elongation and solvent properties are analyzed. The excess chemical potentials of hard-spheroid solutes with aspect ratios x in the range of 15< or =x< or =5, and with volumes between 1 and 20 times that of a solvent molecule, are presented. It is shown that for a given molecular volume and solvent dipole moment (or temperature) a spherical solute has the lowest excess chemical potential and hence the highest solubility, while a prolate solute with aspect ratio x should be more soluble than an oblate solute with aspect ratio 1x. For a given solute molecule, the excess chemical potential increases with increasing temperature; this same trend can be observed in hydrophobic solvation. A scaled-particle theory based on the solvent equation of state and a fitted solute-solvent interfacial tension shows excellent agreement with the simulation results over the whole range of solute elongations and volumes considered. An information-theoretic model based on the solvent density and radial distribution function is less successful, being accurate only for small solute volumes and low solvent densities.     

Orientational effects of hydrogen bonding in liquid-crystalline solutions containing Schiff bases

Intermolecular hydrogen bonding in binary mixtures containing nematogenic Schiff bases as solvents and proton-donating non-mesomorphic solutes has been considered. Reasons for the anomalous concentration dependences of solute order parameters are discussed. A solution structure model of acetic acid in nematic solvents is proposed; constants of complex-dimer equilibrium and coefficients of the orientational correlation of the non-mesogenic solute are calculated on the basis of this model. Hydrogen bonded complex structure using ¹³C NMR has been studied and stability constants in isotropic solutions in chloroform have been calculated. The influence of the solvent orientational ordering on the complex stability is discussed. Data on the solvation isotopic effects in the solutions investigated, which confirm the adequacy of the model are given.     

Linear solvation energy relationship of the solubility of very polar gases: Sulfur dioxide,hydrogen chloride,hydrogen bromide,and ammonia

A linear free energy relationship was found for the log (mole fraction) of solutes in a wide variety of organic solvents with the solvatochromic parameters and the Hildebrand solubility parameter. The solutes were the highly dipolar gases sulfur dioxide, hydrogen chloride, hydrogen bromide, and ammonia at 25°C and 1 atm. partial pressure of the solute. It was found that correlations were greatly improved if solvatochromic parameters for the solvent as a monomer were used rather than the values for the bulk solvent. In solutions with these very dipolar gases, the mole ratio of solute to solvents approaches unity in many of the solutions, so a molecule of solute is interacting primarily with a particular molecule of the solvent. Therefore, the use of the solvatochromic parameters for the solvent as monomer is physically reasonable.     

Towards a modellisation of the solvation energy in multi-component solvents. The interesting case of a charged solute imbedded in a polymer-containing electrolyte solution

In this paper we propose a mean-field theory to calculate the solvation free energy of a charged solute imbedded in a complex multi-component solvent. We considered a solvent made up of a mixture of small (electrolyte solution) and large (polymer) components. The presence of macromolecules ensures reduced mixing entropy among the different solvent components, an effect due to polymer connectivity. The reduced entropy favours strong preferential distribution of a particular solvent even in the presence of weak preferential solute–solvent interactions. In addition, two energy terms must be considered: (a) the interaction between the solute electrostatic potential and the electrolyte solution and (b) the formation of a polymer–solute interface. Because of the different dielectric permittivity of the solvent components, the electrolyte and polymer distribution functions are strongly coupled: ions, indeed, are more solvated in regions of higher local dielectric permittivity arising from the inhomogeneous mixing of solvent and polymer. We combined together the different energy terms in the framework of the de Gennes free energy functional for polymer solutions along with a generalised Poisson–Boltzmann equation developed for inhomogeneous dielectric media. Moreover, the preferential electrolyte solvation in regions of greater polarity was considered by an extension of the Born equation. Setting the polymer dielectric permittivity smaller than the solvent one and making null the specific polymer–solute interactions, we calculated enhanced electrolyte concentration and reduced polymer concentration near the solute surface on raising the solute surface charge density. The theory shows also the breakdown of the widely used separation between electrostatic and surface tension-dependent contributions to solvation energy when non-ideal mixed solvents are considered. In fact, according to the model, the surface tension of such mixed solvents strongly depends on the solute surface charge density: at high potentials the interfacial tension may increase rather than decrease on raising the polymer volume fraction. The theoretical results have been compared with experimental data on polymer+electrolyte solution surface tension and with solubility data of colloidal particles. The comparison evidences the complex behaviour of multi-component solvents going well beyond the trivial weighted average of the dielectric permittivity and surface tension of the isolated chemical components. Deviations from the simple behaviour predicted by an average picture of multi-component solvents could be understood by developing more sophisticated, but still simple, approaches like that proposed in this paper.Contribution to the Jacopo Tomasi Honorary Issue. This paper is dedicated to Jacopo Tomasi. I learned much of the difficult art of transforming complex problems into simple models after reading his early works on solvation energy.     

Detection of Significant Aprotic Solvent Effects on the Conformational Distribution of Methyl 4‐Nitrophenyl Sulfoxide: From Gas‐Phase Rotational to Liquid‐Crystal NMR Spectroscopy

The conformational equilibrium of methyl 4‐nitrophenyl sulfoxide (MNPSO) was experimentally investigated in the gas phase by using microwave spectroscopy and in isotropic and nematic liquid‐crystal solutions, in which the solvents are nonaqueous and aprotic, by using NMR spectroscopy; moreover, it was theoretically studied in vacuo and in solution at different levels of theory. The overall set of results indicates a significant dependence of the solute conformational distribution on the solvent dielectric permittivity constant: when dissolved in low‐polarity media, the most stable conformation of MNPSO proved to be strongly twisted with respect to that in more polar solvents, in which the conformational distribution maximum essentially coincides with that obtained in the gas phase. We discuss a possible explanation of this behavior, which rests on electrostatic solute–solvent interactions and is supported by calculations of the solute electric dipole moment as a function of the torsional angle. This function shows that the least polar conformation of MNPSO is located at a twist angle close to that of the conformational distribution maximum found in less‐polar solvents. This fact, associated with a relatively flat torsional potential, can justify the stabilization of the twisted conformation by the less‐polar solvents.     

Structural Anisotropy in Polar Fluids Subjected to Periodic Boundary Conditions

A heuristic model based on dielectric continuum theory for the long-range solvation free energy of a dipolar system possessing periodic boundary conditions (PBCs) is presented. The predictions of the model are compared to simulation results for Stockmayer fluids simulated using three different cell geometries. The boundary effects induced by the PBCs are shown to lead to anisotropies in the apparent dielectric constant and the long-range solvation free energy of as much as 50%. However, the sum of all of the anisotropic energy contributions yields a value that is very close to the isotropic one derived from dielectric continuum theory, leading to a total system energy close to the dielectric value. It is finally shown that the leading-order contribution to the energetic and structural anisotropy is significantly smaller in the noncubic simulation cell geometries compared to when using a cubic simulation cell.     

Anisotropic solvent model of the lipid bilayer. 1. Parameterization of long-range electrostatics and first solvation shell effects

A new implicit solvation model was developed for calculating free energies of transfer of molecules from water to any solvent with defined bulk properties. The transfer energy was calculated as a sum of the first solvation shell energy and the long-range electrostatic contribution. The first term was proportional to solvent accessible surface area and solvation parameters (σ(i)) for different atom types. The electrostatic term was computed as a product of group dipole moments and dipolar solvation parameter (η) for neutral molecules or using a modified Born equation for ions. The regression coefficients in linear dependencies of solvation parameters σ(i) and η on dielectric constant, solvatochromic polarizability parameter π*, and hydrogen-bonding donor and acceptor capacities of solvents were optimized using 1269 experimental transfer energies from 19 organic solvents to water. The root-mean-square errors for neutral compounds and ions were 0.82 and 1.61 kcal/mol, respectively. Quantification of energy components demonstrates the dominant roles of hydrophobic effect for nonpolar atoms and of hydrogen-bonding for polar atoms. The estimated first solvation shell energy outweighs the long-range electrostatics for most compounds including ions. The simplicity and computational efficiency of the model allows its application for modeling of macromolecules in anisotropic environments, such as biological membranes.     

NMR determination of smectic ordering of probe molecules

The NMR spectra of the three solutes ortho-, meta-, and para-dichlorobenzene in the nematic and smectic A phases of the liquid crystals 8CB and 8OCB are analyzed to yield two orientational order parameters for each solute. Extrapolation of the asymmetry in the energy parameters that describe the orientational ordering in the nematic phase are used to provide estimates of the strength of the nematic potential in the smectic A phase. The experimentally determined asymmetry of the orientational order parameters in the smectic A phase is then used in conjunction with Kobayashi-McMillan theory applied to solutes to give information about the smectic A layering and the nematic/smectic A coupling. In both smectic A solvents, the solute smectic coupling constant, tau, is negative (with the origin fixed at the center of the smectic layer) for all solutes. The signs and relative values of tau indicate that the ortho and para solutes favor the interlayer region while the meta solute is more evenly distributed throughout the layers.     

On the microscopic theory of polar solvation dynamics

A microscopic model of the time-resolved Stokes shift is developed. The model calculates the solvation dynamics by combining the atomic resolution of the solute structure with dipolar dynamics from the polarization structure factors of the homogeneous solvent. Calculations are made for coumarin 153 and quinoxaline optical dyes with atomic geometries and charge distributions taken from quantum calculations. Stokes shift dynamics is calculated and compared to experiment in high-temperature acetonitrile and methanol and in low-temperature 2-methyl-tetrahydrofurane using dielectric relaxation data from experiment.     

The statistical thermodynamics of solvation

Statistical mechanics provides the basis for a rigorous theory of molecules and molecular complexes in chemically inert media. It is found that the free energy of solvation may be calculated, provided the distribution function for the solvent molecules around the solute molecules is known. A relation is deduced between the volume change on dissolving, the free energy of solvation, and the isothermal compressibility of the pure solvent. The general equations are applied to a solution of solid spheres in a structureless medium. The theoretical deductions are found to be in qualitative agreement with the thermodynamic solvation parameters for oxygen and argon in inert solvents.     

Rotational dynamics and conformational kinetics in liquid crystals

Two particular aspects of solute dynamics in ordered media are analysed on the basis of the solution of multivariate diffusion equations: the effects of the solvation dynamics on the rotational motions of dipolar probes in liquid crystal solvents, and the alteration of reaction pathways in isomerization kinetics caused by the solvent order. The introduction of a suitable solvent coordinate allows the interpretation of high frequency contributions in the rotational correlation functions observed by spectroscopic techniques, namely dielectric dispersion, IR and Raman spectroscopy, ESR lineshapes and optical Kerr effect. For molecular systems undergoing conformational changes, a method is offered to evaluate the modification of the torsional barriers resulting from the anisotropic torques modulated by the molecular shape changes along the reaction coordinate.