Solvation phenomena in dilute multicomponent solutions I. Formal results and molecular outlook

We derive second-order thermodynamically consistent truncated composition expansions for the species residual partial molar properties--including volume, enthalpy, entropy, and Gibbs free energy--of dilute ternary systems aimed at the molecular account of solvation phenomena in compressible media. Then, we provide explicit microscopic interpretation of the expansion coefficients in terms of direct and total correlation function integrals over the microstructure of the corresponding infinite dilution reference system, as well as their pressure and temperature derivatives, allowing for the direct prediction of the species partial molar properties from the knowledge of the effective intermolecular interactions. Finally, we apply these formal results (a) to derive consistent expressions for the corresponding properties of the binary system counterparts, (b) to illustrate how the formal expressions converge, at the zero density limit, to those for multicomponent mixtures of imperfect gases obeying the virial equation of state Z = 1 + BPkT, and (c) to discuss, and highlight with examples from the literature, the thermodynamic inconsistencies encountered in the currently available first-order truncated expansions, by pinpointing the mathematical origin and physical meaning of the inconsistencies that render the first-order truncated expansions invalid.     

Forward--backward semiclassical dynamics with information-guided noise reduction for a molecule in solution

The forward--backward semiclassical dynamics (FBSD) methodology is used to obtain expressions for time correlation functions of a system (atom or molecule) in solution. We use information-guided noise reduction (IGNoR) [Makri, N. Chem. Phys. Lett. 2004, 400, 446] to minimize the statistical error associated with the Monte Carlo integration of oscillatory functions. This is possible by reformulating the correlation function in terms of an oscillatory solvent-dependent contribution whose integral can be obtained analytically and a slowly varying function obtained via a grid-based iterative evaluation of solute properties. Knowledge of the exact integral of the oscillatory function, combined with correlated statistics, leads to partial cancellation of the Monte Carlo error. Application on a one-dimensional solute-solvent model shows a substantial improvement of convergence in the IGNoR-enhanced FBSD correlation function for a fixed number of Monte Carlo samples. The reduction of statistical error achieved by using the IGNoR methodology becomes more significant as the number of solvent particles increases.     

Scattering functions of semidilute solutions of polymers in a good solvent

Expressions for analyzing small-angle scattering data from semidilute solutions of polymers in a good solvent over a broad range of scattering vectors are examined. Three different scattering function expressions are derived from Monte Carlo simulations. The expressions are similar to those of polymer reference interaction site models, with a scattering-vector-dependent direct correlation function. In the most advanced model, the screening of excluded-volume interactions beyond the overlap concentration is taken into account. Two simpler expressions, in which the screening of excluded-volume interactions is not included, are also applied. The three models are tested against small-angle neutron scattering (SANS) experiments on polystyrene in deuterated toluene for a broad range of molar masses and concentrations over a wide range of scattering vectors. For each model, simultaneous fits to all the measured scattering data are performed. The most advanced model excellently reproduces the SANS data over the full range of the parameters. The two simpler models fit the data almost equally well. On the basis of an extensive study, an optimal fitting strategy can be recommended for experimentalists, who want to analyze small-angle scattering data from polymers at any concentration. For data sets that do not contain data on the single-chain scattering function, the simpler model is recommended; it uses a direct correlation function equal to the form factor of an infinitely thin rod, which is independent of the concentration and molar mass. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 3081–3094, 2004     

Pair correlation functions and a free energy functional for the nematic phase

In this paper we have presented the calculation of pair correlation functions in a nematic phase for a model of spherical particles with the long-range anisotropic interaction from the mean spherical approximation (MSA) and the Percus-Yevick (PY) integral equation theories. The results found from the MSA theory have been compared with those found analytically by Holovko and Sokolovska [J. Mol. Liq. 82, 161 (1999)]. A free energy functional which involves both the symmetry conserving and symmetry broken parts of the direct pair correlation function has been used to study the properties of the nematic phase. We have also examined the possibility of constructing a free energy functional with the direct pair correlation function which includes only the principal order parameter of the ordered phase and found that the resulting functional gives results that are in good agreement with the original functional. The isotropic-nematic transition has been located using the grand thermodynamic potential. The PY theory has been found to give a nematic phase with pair correlation function harmonic coefficients having all the desired features. In a nematic phase the harmonic coefficient of the total pair correlation function h(x1,x2) connected with the correlations of the director transverse fluctuations should develop a long-range tail. This feature has been found in both the MSA and PY theories.     

Theoretical direct correlation function for two-dimensional fluids of monodisperse hard spheres

The direct correlation function plays an important role in describing the effects of the structure of particle systems with respect to light diffraction, x-ray diffraction as well as transmission and transmission fluctuations of radiation through a dense suspension. In this paper, the direct correlation function for a monolayer of monodisperse hard spheres or disks is derived theoretically. Based on the approximation of Baus and Colot [Phys. Rev. A 36, 3912 (1987)] and the equation of state for a fluid of hard disks by Santos et al. [J. Chem. Phys. 103, 4622 (1995)], we propose a new direct correlation function, which compares well to the approximate analytical expressions and gives a good prediction of the structure factor in a wide range of monolayer density or suspension concentration. The resulting radial distribution function also agrees well with Monte Carlo computer simulation data. The corresponding contact values of the radial distribution function compare well with the results of analytic approximations, numerical solutions, and computer simulations. Our proposed direct correlation function is applied to the transmission fluctuation spectrometric study. Experimental results show good agreement with the theory.     

Probability Properties of Multi－contact in Protein Molecules

The compact conformations of polymers are important because the native conformations of all bio-polymers with certain function are highly compact. The properties of mutil-contact bio-polymer chains were studied by Gaussian statistics of the random-flight chain. The theoretical expressions(were given, also), the calculations of probability distributions and correlation functions for different topologic cases were derived and made respectively. Comparison between single, double and triple contacts was also made. By means of setting the parameters, the results of the current calculations of the multiple contacts are just the same as those calculated by single, double or tripe contacts separately. It is a useful method to investigate native conformations of biopolymers. The probabilities of multi contacts and correlation functions between chain‘s contacts were calculated for the Gaussian chains. Because the bond probability distributions are Gaussian‘s distributions, the probability distributions of the separations of various points along the chains are always consecutive. All the contacts may break up into several groups, and each group consists of many contacts. Here we investigated the probability distribution from one group to three groups of contacts.     

The sound velocities in dense fluids from distribution functions

The sound velocities and adiabatic compressibilities in dense fluids have been evaluated using three known analytical expressions for radial distribution functions (RDFs). Using such approach not only tests the power of distribution functions theory in predicting the sound velocities and adiabatic compressibilities, but also specifies better expressions in determining these properties. To calculate these quantities, the variation of RDF with density and temperature is required. Therefore, we should have analytical expressions which explicitly present RDF as a function of temperature, density and interparticle distance. It is shown that if an expression is used which properly presents RDFs as a function of interparticle distance, density and temperature, it is possible to calculate sound velocities and adiabatic compressibilities from distribution function theory.     

The changes in free energies and entropies from analytical radial distribution functions

The changes in Helmholtz free energies and entropies in dense fluids have been evaluated using three known analytical expressions for radial distribution functions (RDFs) of Lennard–Jones (L-J) fluid. This method provides a simpler and a more expeditious way for the calculation of free energy and entropy in L-J dense fluids through statistical mechanics. Previously, integral equations or perturbation theories were used for this purpose. Such approach not only tests the power of analytical distribution functions in predicting the changes in Helmholtz free energies and entropies, but also specifies better expressions in determining these properties. The results are compared with experimental data and an accurate analytic equation of state for the L-J fluid. It is shown if an expression properly presents RDFs as a function of interparticle distance, density and temperature, it is possible to calculate the changes in Helmholtz free energies and entropies from analytical distribution functions.     

Analytic gradients for density cumulant functional theory: The DCFT-06 model

Density cumulant functional theory (DCFT) is one of a number of nascent electron correlation methods that are derived from reduced density matrices and cumulants thereof, instead of the wavefunction. Deriving properties from the density cumulant naturally yields methods that are size extensive and size consistent. In this work, we derive expressions for the analytic gradient, with respect to an external perturbation, for the DCFT-06 variant of density cumulant functional theory. Despite the fact that the DCFT-06 energy functional is stationary with respect to the density cumulant, the analytic gradients of the energy require the solution of perturbation-independent equations for both orbital and cumulant response. These two sets of linear response equations are coupled in nature and are solved iteratively with the solution of orbital and cumulant response equations each macroiteration, exhibiting rapid convergence. The gradients are implemented and benchmarked against coupled cluster theory with single and double excitations (CCSD) and CCSD with perturbative triple excitations [CCSD(T)], as well as accurate empirically corrected experimental data, for a test set comprising 15 small molecules. For most of the test cases, results from DCFT-06 are closer to CCSD(T) and empirical data than those from CCSD. Although the total energy and analytic gradient have the same asymptotic scaling, the present experience shows that the computational cost of the gradient is significantly lower.     

接枝共聚物结构的理论分析

<正> 在非交联型聚合物中,接枝共聚物的分子结构最为复杂,它既取决于整个共聚物的分子量分布和平均分子量,又与骨架和接枝链的链长分布和平均链长紧密相关。此外,接枝产物的接枝数分布、平均接枝数、接枝点间链段的平均分子量及其分子量分布等都是描述这类共聚物的重要参数。因此,接枝共聚物分子结构的精确表征,是一个尚待解决的研究     

Configurational statistics of macromolecules in solution from correlations of liquid molecules

We address the relevant quest for a simple formalism describing the microstructure of liquid solutions of polymer chains. On the basis of a recent relativistic-type picture of self-diffusion in (simple) liquids named Brownian relativity (BWR), a covariant van Hove's distribution function in a Vineyard-like convolution approximation is proposed to relate the statistical features of liquid and chain molecules forming a dilute polymer solution. It provides an extension of the Gaussian statistics of ideal chains to correlated systems, allowing an analysis of macromolecular configurations in solution by the only statistical properties of the liquid units (and vice versa). However, the mathematical solution to this issue is not straightforward because, when the liquid and polymer van Hove's functions are equated, an inverse problem takes place. It presents some conceptual analogies with a scattering experiment in which the correlation of the liquid molecules acts as the radiation source and the macromolecule as the scatterer. After inverting the equation by a theorem coming from the Tikhonov's approach, it turns out that the probability distribution function of a real polymer can be expressed from a static Ornstein-Uhlenbeck process, modified by correlations. This result is used to show that the probability distribution of a true self-avoiding walk polymer (TSWP) can be modeled as a universal Percus-Yevick hard-sphere solution for the total correlation function of the liquid units. This method suits in particular the configurational analysis of single macromolecules. The analytical study of arbitrary many-polymer systems may require further mathematical investigation.     

A Single and a Dual-Fractal Analysis of Analyte-Receptor Binding Kinetics for Surface Plasmon Resonance Biosensor Applications

The diffusion-limited binding kinetics of analyte in solution to either a receptor immobilized on a surface or to a receptorless surface is analyzed within a fractal framework for a surface plasmon resonance biosensor. The data is adequately described by a single- or a dual-fractal analysis. Initially, the data was modeled by a single-fractal analysis. If an inadequate fit was obtained then a dual-fractal analysis was utilized. The regression analysis provided by Sigmaplot (32) was used to determine if a single fractal analysis is sufficient or if a dual-fractal analysis is required. In general, it is of interest to note that the binding rate coefficient and the fractal dimension exhibit changes in the same direction (except for a single example) for the analyte-receptor systems analyzed. Binding rate coefficient expressions as a function of the fractal dimension developed for the analyte-receptor binding systems indicate, in general, the high sensitivity of the binding rate coefficient on the fractal dimension when both a single- and a dual-fractal analysis is used. For example, for a single-fractal analysis and for the binding of human endothelin-1 (ET-1) antibody in solution to ET-115-21.BSA immobilized on a surface plasmon resonance (SPR) surface (33), the order of dependence of the binding rate coefficient, k, on the fractal dimension, Df, is 6.4405. Similarly, for a dual-fractal analysis and for the binding of 10(-6) to 10(-4) M bSA in solution to a receptorless surface (direct binding to SPR surface) (41) the order of dependence of k1 and k2 on Df1 and Df2 were -2.356 and 6.241, respectively. Binding rate coefficient expressions are also developed as a function of the analyte concentration in solution. The binding rate coefficient expressions developed as a function of the fractal dimension(s) are of particular value since they provide a means to better control SPR biosensor performance by linking it to the degree of heterogeneity that exists on the SPR biosensor surface. Copyright 1999 Academic Press.     

Study of the pair correlations between p-nitroaniline molecules in solution by depolarized hyper-Rayleigh scattering

The concentration dependence of the hyper-Rayleigh scattering depolarization ratios of p-nitroaniline in solution was obtained and the results were compared with theory. It was found that the experimental data can be theoretically accounted for by using a pair distribution function that includes only direct correlation, with the molecules interact through a dipolar hard-sphere potential. The results show that short-range dipole-dipole interactions are responsible for the correlation between pairs of p-nitroaniline molecules in solution.     

Closed-form expressions for correlated density matrices: application to dispersive interactions and example of He2

Empirically correlated density matrices of N-electron systems are investigated. Closed-form expressions are derived for the one- and two-electron reduced density matrices from a pairwise correlated wave function. Approximate expressions are then proposed which reflect dispersive interactions between closed-shell centrosymmetric subsystems. Said expressions clearly illustrate the consequences of second-order correlation effects on the reduced density matrices. Application is made to a simple example: the He(2) system. Reduced density matrices are explicitly calculated, correct to second order in correlation, and compared with approximations of independent electrons and independent electron pairs. The models proposed allow for variational calculations of interaction energies and equilibrium distance as well as a clear interpretation of dispersive effects on electron distributions. Both exchange and second order correlation effects are shown to play a critical role on the quality of the results.     

Analytic solutions for Baxter's model of sticky hard sphere fluids within closures different from the Percus-Yevick approximation

We discuss structural and thermodynamical properties of Baxter's adhesive hard sphere model within a class of closures which includes the Percus-Yevick (PY) one. The common feature of all these closures is to have a direct correlation function vanishing beyond a certain range, each closure being identified by a different approximation within the original square-well region. This allows a common analytical solution of the Ornstein-Zernike integral equation, with the cavity function playing a privileged role. A careful analytical treatment of the equation of state is reported. Numerical comparison with Monte Carlo simulations shows that the PY approximation lies between simpler closures, which may yield less accurate predictions but are easily extensible to multicomponent fluids, and more sophisticate closures which give more precise predictions but can hardly be extended to mixtures. In regimes typical for colloidal and protein solutions, however, it is found that the perturbative closures, even when limited to first order, produce satisfactory results.     

Energy and capacity time correlation functions for investigation of vibrational energy relaxation

Vibrational energy relaxation of a diatomic solute in a liquid solvent is investigated by means of the generalized Langevin equation. The vibrational energy, velocity and capacity time correlation functions (TCFs) are considered. It is shown that the detailed structure of the energy TCF contains an initial fast (subpicosecond) decay segment that is followed by weak oscillations on the background of an exponential relaxation curve. The direct method for evaluating the relaxation rate constant from equilibrium molecular dynamics simulations of a flexible solute is proposed and implemented. The closed form expressions for the memory function and for the relaxation rate constant in terms of quantities accessible from the simulations are derived. The simulation results for rigid and flexible solutes are compared and analyzed.     

三元多嵌段共聚物分子链结构的统计表征

本文根据三元多嵌段共聚反应的一般模型,运用母函数方法,严格推导出了型三元多嵌段共聚物的分子量分布及平均分子量的解析表达式.对几种具有特殊分布的预聚体进行了讨论.     

Transport properties of long straight nano-channels in electrolyte solutions: a systematic approach

The principle of local thermodynamic equilibrium is systematically employed for obtaining various transport properties of long straight nano-channels. The concept of virtual solution is used to describe situations of non-negligible overlap of diffuse parts of electric double layers (EDLs) in nano-channels. Generic expressions for a variety of transport properties of long straight nano-channels are obtained in terms of quasi-equilibrium distribution coefficients of ions and functionals of quasi-equilibrium distribution of electrostatic potential. Further, the Poisson-Boltzmann approach is used to specify these expressions for long straight slit-like nano-channels. In the approximation of non-overlapped diffuse parts of double electric layers in nano-channels, simple analytical expressions are obtained for the apparent electrophoretic mobilities of (trace) analytes of arbitrary charge as well as for the salt reflection coefficient (osmotic pressure), salt diffusion permeability and electro-viscosity (electrokinetic energy conversion). The approximate solutions are compared with the results of rigorous solution of non-linearized Poisson-Boltzmann equation, and the accuracy of approximation is shown to be typically excellent when the nano-channel half-height exceeds ca.3 Debye screening lengths. Due to non-negligible electrostatic adsorption of ions by nano-channels, the apparent electrophoretic mobilities of counter-ionic analytes in nano-channels are smaller than in micro-channels whereas those of co-ionic analytes are larger. This dependence on the charge is useful for the separation of analytes of close electrophoretic mobilities. The osmotic pressure is shown to be positive, negative or pass through maxima as a function of applied salt-concentration difference within a fairly narrow range of ratios of nano-channel height to the Debye screening length. The diffusion permeability of charged nano-channels to single salts is demonstrated (for the first time) to be typically larger than that of neutral nano-channels of the same dimensions due to electrical facilitation of salt diffusion.