Molecular thermodynamics of flexible-molecule fluids

Some of the most important problems encountered -in the statistical-mechanical treatment of dense, complex molecular fluids are indicated. Particular attention is given to the application of the generalized van der Waals concept to flexible-molecule fluids. Perturbed-hard-chain theory and chain-of-rotators theory are reviewed. This is followed by a brief summary of recent experimental work on thermodynamic excess properties of liquid mixtures of type (Y + an n-alkane). Several novel effects connected with local correlation of molecular orientation in the n-alkanes, with medium-induced conformational changes and with preferential orientation of strongly polar molecules are discussed (order/disorder phenomena).     

Chiral interaction and thermodynamics

In a recent series of articles a novel interaction between chiral molecules such as proteins and a polar solvent has been proposed. This interaction is criticized by Tannhauser and Kuper on the basis of thermodynamics, and the purpose of this Reply is to show that thermodynamics per se can raise no objection to this interaction. Special attention is given to time-reversal invariance, which is violated by this interaction. A temporal interplay between relaxation and stochastic times τ_R and τ can cause a non-ergodic behavior of the system.     

Chemical thermodynamics and theoretical models

This paper presents a brief discussion of the relationship between classical thermodynamics and theoretical models of chemical systems. Applications of thermodynamics to real chemical systems and to related theoretical models concerning “complexes in nonelectrolyte solutions” and “substituent and solvent effects in organic chemistry” are described, with particular emphasis on research that my colleagues and I have done in these areas. This paper concludes with a brief autobiographical section on “people and places”.     

Ion-induced nucleation in polar one-component fluids

We present a Ginzburg-Landau theory of ion-induced nucleation in a gas phase of polar one-component fluids, where a liquid droplet grows with an ion at its center. By calculating the density profile around an ion, we show that the solvation free energy is larger in gas than in liquid at the same temperature on the coexistence curve. This difference much reduces the nucleation barrier in a metastable gas.     

Dielectric permittivity profiles of confined polar fluids

The dielectric response of a simple model of a polar fluid near neutral interfaces is examined by a combination of linear response theory and extensive molecular dynamics simulations. Fluctuation expressions for a local permittivity tensor epsilon(r) are derived for planar and spherical geometries, based on the assumption of a purely local relationship between polarization and electric field. While the longitudinal component of epsilon exhibits strong oscillations on the molecular scale near interfaces, the transverse component becomes ill defined and unphysical, indicating nonlocality in the dielectric response. Both components go over to the correct bulk permittivity beyond a few molecular diameters. Upon approaching interfaces from the bulk, the permittivity tends to increase, rather than decrease as commonly assumed, and this behavior is confirmed for a simple model of water near a hydrophobic surface. An unexpected finding of the present analysis is the formation of "electrostatic double layers" signaled by a dramatic overscreening of an externally applied field inside the polar fluid close to an interface. The local electric field is of opposite sign to the external field and of significantly larger amplitude within the first layer of polar molecules.     

Statistical thermodynamics of fluids with both dipole and quadrupole moments

New Gibbs ensemble simulation data for a polar fluid modeled by a square-well potential plus dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole interactions are presented. This simulation data is used in order to assess the applicability of the multipolar square-well perturbation theory [A. L. Benavides, Y. Guevara, and F. del Ri?o, Physica A 202, 420 (1994)] to systems where more than one term in the multipole expansion is relevant. It is found that this theory is able to reproduce qualitatively well the vapor-liquid phase diagram for different multipolar moment strengths, corresponding to typical values of real molecules, except in the critical region. Hence, this theory is used to model the behavior of substances with multiple chemical bonds such as carbon monoxide and nitrous oxide and we found that with a suitable choice of the values of the intermolecular parameters, the vapor-liquid equilibrium of these species is adequately estimated.     

可聚合凝胶因子的合成及其有机凝胶热力学研究

合成了一种可聚合凝胶因子（Gelator)4，4’-二（α－甲基丙烯酰氧基－1， 3－亚乙氧基羰基丙酰氨基）二苯甲烷（BMDM）。BMDM能在二苯醚、甲苯、二甲苯 、氯苯等含苯环的低极性溶剂中形成热可逆的物理凝胶，而在二甲基甲酰胺、乙醇 和二氯甲烷等强极性溶剂中则易溶解，在乙醚、石油醚等非极性溶剂中不可溶解。 FT－IR和DSC研究该凝胶因子在二苯醚中形成的凝胶，发现该凝胶因子是通过氢键 等次价键力相互作用而聚集、自我组装形成凝胶的。利用溶胶－凝胶相转变和DSC 数据，研究了该有机凝胶的热力学参数和性质，BMDM在二苯醚中形成的凝胶聚集体 间的范德华弱相互作用焓为ΔH=0.8kJ·mol^-1，凝胶－溶胶相转变热焓为ΔH=22. 3kJ·mol^-1。     

Soft core thermodynamics from self-consistent hard core fluids

In an effort to generalize the self-consistent Ornstein-Zernike approximation (SCOZA)-an accurate liquid state theory that has been restricted so far to hard core systems-to arbitrary soft core systems we study a combination of SCOZA with a recently developed perturbation theory. The latter was constructed by Ben-Amotz and Stell [J. Phys. Chem. B 108, 6877 (2004)] as a reformulation of the Weeks-Chandler-Andersen [J. Chem. Phys. 54, 5237 (1971)] perturbation theory directly in terms of an arbitrary hard sphere reference system. We investigate the accuracy of the combined approach for the Lennard-Jones fluid in comparison with simulation data and pure perturbation theory predictions and determine the dependence of the thermodynamic properties and the phase behavior on the choice of the effective hard core diameter of the reference system.     

Optical manipulation of charged microparticles in polar fluids

In this study, we report a systematic study of the response of a charged microparticle confined in an optical trap and driven by electric fields. The particle is embedded in a polar fluid, hence, the role of ions and counterions forming a double layer around the electrodes and the particle surface itself has been taken into account. We analyze two different cases: (i) electrodes energized by a step‐wise voltage (DC mode) and (ii) electrodes driven by a sinusoidal voltage (AC mode). The experimental outcomes are analyzed in terms of a model that combines the electric response of the electrolytic cell and the motion of the trapped particle. In particular, for the DC mode we analyze the transient particle motion and correlate it with the electric current flowing in the cell. For the AC mode, the stochastic and deterministic motion of the trapped particle is analyzed either in the frequency domain (power spectral density, PSD) or in the time domain (autocorrelation function). Moreover, we will show how these different approaches (DC and AC modes) allow us, assuming predictable the applied electric field (here generated by plane parallel electrodes), to provide accurate estimation (3%) of the net charge carried by the microparticle. Vice versa, we also demonstrate how, once predetermined the charge, the trapped particle acts as a sensitive probe to reveal locally electric fields generated by arbitrary electrode geometries (in this work, wire‐tip geometry).     

Structure and thermodynamics of hard-core Yukawa fluids: thermodynamic perturbation approaches

The thermodynamic perturbation theories, which are based on the power series of a coupling constant (λ-expansion), have been proposed for studying the structural and thermodynamic properties of a hard-core Yukawa (HCY) fluid: one (A1-approximation) is the perturbation theory based on the hard-sphere repulsion as a reference system. The other (A2-approximation) is the perturbation theory based on the reference system which incorporates both the repulsive and short-range attractive interactions. The first-order mean-spherical approximation (FMSA) provided by Tang and Lu [J. Chem. Phys. 99, 9828 (1993)] has been employed for investigating the thermodynamic properties of a HCY fluid using the alternative method via the direct correlation function. The calculated results show that (i) the A1 and A2 approximations are in excellent agreements with previous computer simulation results in the literature and compare with the semi-empirical works of Shukla including the higher-order free energy terms, (ii) the A1 and A2 approximations are better than the FMSA and the mean-spherical approximation, (iii) the A2-approximation compares with the A1-approximation, even though the perturbation effect of an A2-approximation is much smaller than that of an A1-approximation, and that (iv) the FMSA study is particularly of advantage in providing the structure and thermodynamics in a simple and analytic manner.     

An extended rism equation for molecular polar fluids

The RISM integral equation is extended to molecules with charged sites via a renormalization of the Coulomb potentials and the introduction of appropriate closure relations. For a fluid of diatomics with atomic charges of ±0.2 e the equation yields site-site correlation functions in qualitative agreement with those from computer simulation.     

Picosecond observations of electron solvation in dilute polar fluids

Electron solvation has been studied in dilute polar fluids in order to quantify the role of the fluid in the proposed mechanisms of electron trapping and solvation. In a series of dilute alcohol-alkane systems, the picosecond evolution of the absorption spectrum is shown to be a sensitive function of the local liquid structure and dynamics. A solvation mechanism is outlined which correlates the absorption and mobility data from neat and dilute polar fluids.     

Energy controlled insertion of polar molecules in dense fluids

We present a method to search low energy configurations of polar molecules in the complex potential energy surfaces associated with dense fluids. The search is done in the configurational space of the translational and rotational degrees of freedom of the molecule, combining steepest-descent and Newton-Raphson steps which embed information on the average sizes of the potential energy wells obtained from prior inspection of the liquid structure. We perform a molecular dynamics simulation of a liquid water shell which demonstrates that the method enables fast and energy-controlled water molecule insertion in aqueous environments. The algorithm finds low energy configurations of incoming water molecules around three orders of magnitude faster than direct random insertion. This method represents an important step towards dynamic simulations of open systems and it may also prove useful for energy-biased ensemble average calculations of the chemical potential.     

Concentration fluctuations in binary mixtures of model polar fluids

We present calculations of the mean square concentration fluctuations, S_α(0), for binary mixtures of model polar fluids. The assumed pair interactions are taken to be of the forms of a hard core plus either dipole-dipole, quadrupole-quadrupole or dipole- quadrupole terms. The calculations are carried out within a mean field approximation. We have considered in some detail the interplay between size differences and the difference in the strength and range of the potentials in deciding how S_α(0) deviates from the ideal behaviour.     

A new equation of state for polar and nonpolar pure fluids

Chung, T.H., Khan, M.M., Lee, L.L. and Starling, K.E., 1984. A new equation of state for polar and nonpolar pure fluids. Fluid Phase Equilibria, 17: 351–372A new equation of state based on the concept of perturbation theory and the hard-convex-body equation of state has been developed successfully for nonpolar compounds. The equation can predict the thermodynamic properties (density, enthalpy departure and vapor pressure) of a wide range of pure fluids from small, spherical (argon-like) molecules to large, structurally complex molecules. For nonpolar compounds, the equation employs three parameters: the shape, size and energy parameters. For normal paraffins, the size parameter (hard-core volume) is related to the measurable van der Waals volume given by Bondi. For most other compounds, it is related to the critical volume. The shape-parameter values reflect the structure and degree of acentricity of the compound of interest. The equation has been extended to polar and associating compounds by using the mean-potential model. For polar compounds, a fourth parameter is required. The equation has been tested extensively for polar (dipolar and quadrupolar) and hydrogen-bonding compounds. The applicability of this equation for such a wide variety of substances provides an important first step in the development of a composition-dependent equation of state for mixtures.     

Dispersion corrected interaction of polar and nonpolar fluids confined within carbon nanotubes: Density functional theoretical analysis using Grimme's D3 scheme

The structural and flow characteristics of fluids within carbon nanotube (CNT) is dictated by the interaction of fluid molecules within the nanocavity of CNT. Therefore, in the present study, dispersion corrected density functional theory has been used to investigate the structure and interaction of polar and nonpolar molecules within CNT. The present study shows that there is profound effect on the interaction due to dispersion. The interaction energy of the confined water was found to be reduced with increasing distance of the water molecule from the wall of the CNT. The water is preferentially adsorbed over methane due to stronger interaction with CNT over methane. Further, water is preferentially adsorbed over methanol molecule when interaction is calculated without dispersion but after inclusion of dispersion interaction, the calculated results show that the methanol–CNT interaction is stronger than that of water molecule and hence preferentially adsorbed within the CNT as revealed from MD simulation. The present calculation reveals that that the effect of CNT confinement on the IR spectra of the single file water is quite considerable compared to the IR spectra of tetrahedral bulk water cluster. Therefore, the present results might be useful for the separation of polar molecule from nonpolar molecule during fabrication of CNT‐based filter and purification system.     

Solvation force, structure and thermodynamics of fluids confined in geometrically rough pores

The effect of periodic surface roughness on the behavior of confined soft sphere fluids is investigated using grand canonical Monte Carlo simulations. Rough pores are constructed by taking the prototypical slit-shaped pore and introducing unidirectional sinusoidal undulations on one wall. For the above geometry our study reveals that the solvation force response can be phase shifted in a controlled manner by varying the amplitude of roughness. At a fixed amplitude of roughness, a, the solvation force for pores with structured walls was relatively insensitive to the wavelength of the undulation, lambda for 2.3/=0.5. The predictions of the superposition approximation, where the solvation force response for the rough pores is deduced from the solvation force response of the slit-shaped pores, was in excellent agreement with simulation results for the structured pores and for lambda/sigma(ff)>/=7 in the case of smooth walled pores. Grand potential computations illustrate that interactions between the walls of the pore can alter the pore width corresponding to the thermodynamically stable state, with wall-wall interactions playing an important role at smaller pore widths and higher amplitudes of roughness.     

Computationally useful bridge diagram series for the structure and thermodynamics of Lennard-Jones fluids

The first two orders of bridge diagrams, those with two and three field points, have been calculated exactly for the Lennard-Jones fluid for several isotherms. The method of calculation was one of expansion in Legendre polynomials, and the dependence of the method on the number of polynomials needed for accurate results was investigated. Thermodynamic and structural properties of the Lennard-Jones fluid calculated from integral equation methods with the inclusion of bridge diagrams were found to be systematically improved. Two attempts at predicting the missing bridge diagrams of even higher order were discussed. The first, which uses the functional form of those diagrams that were calculated exactly, showed no significant improvement. The second, a series sum based on the first two orders of calculated diagrams and motivated by the success of a similar heuristic sum in the case of hard spheres, was extremely successful. When the series sum was employed, thermodynamic and structural quantities were improved to the point where the difference between simulation results and integral equation results was of the same order as the error in the simulations themselves. Received: 9 January 1997 / Accepted 28 January 1997     

Development of surface-SFED models for polar solvents

We developed surface grid-based solvation free energy density (Surface-SFED) models for 36 commonly used polar solvents. The parametrization was performed with a large and diverse set of experimental solvation free energies mainly consisting of combinations of polar solvent and multipolar solute. Therefore, the contribution of hydrogen bonds was dominant in the model. In order to increase the accuracy of the model, an elaborate version of a previous hydrogen bond acidity and basicity prediction model was introduced. We present two parametrizations for use with experimentally determined (Surface-SFED/HB(exp)) and empirical (Surface-SFED/HB(cal)) hydrogen bond acidity and basicity values. Our computational results agreed well with experimental results, and inaccuracy of empirical hydrogen bond acidity and basicity values was the main source of error in Surface-SFED/HB(cal). The mean absolute errors of Surface-SFED/HB(exp) and Surface-SFED/HB(cal) were 0.49 and 0.54 kcal/mol, respectively.