高聚物亚浓溶液区标度律的蒙特卡洛模拟

利用蒙特卡洛模拟方法对线形链和星形键在亚浓溶液区的标度律作了验证与讨论。我们的计算结果证明线形链渗透压Π与体积分数Φ在亚浓溶液区的标度规律是与德热纳的理论结果相一致的。同时还证明了星形链在亚浓溶液区具有与线形链相同的标度行为。     

Identifying potential binding modes and explaining partitioning behavior using flexible alignments and multidimensional scaling

A method is described for the rapid and automatic analysis of flexible molecular alignments using multidimensional scaling and a normalized scoring scheme. A projection scheme was devised to separate orientational and conformational effects. It is shown that the approach can be utilized for the identification of common binding orientations or to the study of differences in partioning behavior. It is suggested that the method can be employed as a novel approach exploring molecular similarity as a dynamic property, so that it includes aspects of motion (by way of mutual orientations), conformations and molecular properties.     

Importance of the Debye interaction in organic solutions: Henry's law constants for polar liquids in nonpolar solvents and vice versa

Henry's law constants have been determined for -butyrolactone (BL), ethyl acetate (EA), and 2-methyl-3-pentanol (MEP) in mixtures of iso-octane (ISO) and toluene (TOL), for BL, EA, TOL, and ISO in cinnamaldehyde (CIN) and for TOL and ISO in each other and in BL. From these data and published vapor pressures, the activity coefficients at infinite dilution and the standard molar Gibbs free energy of transfer, G ₂ ⁰ of the solutes from dilute solution in ISO to dilute solution in each solvent medium have been calculated. The different behavior patterns of BL and EA are attributed to differences in their abilities to exist in different conformations possessing different dipole moments. For polar solutes, G ₂ ⁰ decreases with increasing polarizability of the solvent and with increasing dipole moment of the solute, suggesting increased contributions from dipole-induced dipole (Debye) interactions. The sigmoidal plot of G ₂ ⁰ against the change in pair potential energy calculated from the classical expressions suggests that G ₂ ⁰ seriously underestimates the strength of the Debye interactions in comparison with the London interactions.     

Dynamic percolation transition induced by phase separation: A Monte Carlo analysis

The percolation transition of geometric clusters in the three-dimensional, simple cubic, nearest neighbor Ising lattice gas model is investigated in the temperature and concentration region inside the coexistence curve. We consider quenching experiments, where the system starts from an initially completely random configuration (corresponding to equilibrium at infinite temperature), letting the system evolve at the considered temperature according to the Kawasaki spinexchange dynamics. Analyzing the distributionn _l(t) of clusters of sizel at timet, we find that after a time of the order of about 100 Monte Carlo steps per site a percolation transition occurs at a concentration distinctly lower than the percolation concentration of the initial random state. This dynamic percolation transition is analyzed with finite-size scaling methods. While at zero temperature, where the system settles down at a frozen-in cluster distribution and further phase separation stops, the critical exponents associated with this percolation transition are consistent with the universality class of random percolation, the critical behavior of the transient time-dependent percolation occurring at nonzero temperature possibly belongs to a different, new universality class.     

Hydrodynamic equations for attractive particle systems on ℤ

Hydrodynamic properties for a class of nondiffusive particle systems are investigated. The method allows one to study local equilibria for a class of asymmetric zero-range processes, and applies as well to other models, such as asymmetric simple exclusion and misanthropes. Attractiveness is an essential ingredient. The hydrodynamic equations present shock wave phenomena. Preservation of local equilibrium is proven to hold away from the shocks. The problem of breakdown of local ergodicity at the shocks, which was investigated by D. Wick in a particular model, remains open in this more general setup.     

Modeling high-pressure densities at wide temperature range with volume scaling: Cyclohexane + n-hexadecane mixtures

High-pressure density data for cyclohexane + n-hexadecane mixtures at a wide temperature range was modeled with several classical equations of state (EOS) and correlative models. A modification for softening the co-volume and another for a volume scaling of the Peng–Robinson EOS (VS-PR) were proposed. The VS-PR model is able to correlate the pure component experimental data employing only five adjustable parameters, with root-mean-square deviation (RMSD) between calculated and experimental densities essentially within the experimental error. This result is superior to widely used approaches, i.e., a six parameter Tait model and six parameter volume translations (temperature and pressure dependent) for Peng–Robinson and Patel–Teja EOS. The VS-PR model also represents well the isobaric thermal expansion and the isothermal compressibility coefficients of the pure cyclohexane, a small naphthenic substance as well as a long chain n-alkane hydrocarbon, n-hexadecane. When modeling the mixture data, the use of VS-PR model of pure components along with the Redlich–Kister expansion, truncated at the first term, the density was correlated within a RMSD only 60% greater than the experimental error. The proposed model is able to accurately represent all the tested mixture data with a relatively small number of parameters.     

Linear scaling local correlation approach for solving the coupled cluster equations of large systems

A linear scaling local correlation approach is proposed for approximately solving the coupled cluster doubles (CCD) equations of large systems in a basis of orthogonal localized molecular orbitals (LMOs). By restricting double excitations from spatially close occupied LMOs into their associated virtual LMOs, the number of significant excitation amplitudes scales only linearly with molecular size in large molecules. Significant amplitudes are obtained to a very good approximation by solving the CCD equations of various subsystems, each of which is made up of a cluster associated with the orbital indices of a subset of significant amplitudes and the local environmental domain of the cluster. The combined effect of these two approximations leads to a linear scaling algorithm for large systems. By using typical thresholds, which are designed to target an energy accuracy, our numerical calculations for a wide range of molecules using the 6-31G or 6-31G* basis set demonstrate that the present local correlation approach recovers more than 98.5% of the conventional CCD correlation energy.     

Hygrometric Determination of Water Activities and Osmotic and Activity Coefficients of NH4Cl–LiCl–H2O at 25°C

The mixed aqueous electrolyte system of ammonium and lithium chlorides has been studied by the hygrometric method at 25°C. The relative humidities of this system are measured at total molalities from 0.3 to 6 mol-kg^{– 1} for different ionic-strength fractions y of NH₄Cl with y = 0.33, 0.50, and 0.67. The data obtained allow the deduction of new water activities and osmotic coefficients. The experimental results are compared with the predictions of the ECA (extended composed additivity) law proposed in our previous work. The Zdanovskii–Stokes–Robinson (ZSR), the Robinson–Stokes (RS), Reilly–Wood–Robinson (RWR), the Pitzer, and the Lietzke–Stoughton (LS II) models are also compared with our results. Predictions made using these models are, in general, consistent with our results. From these measurements, new Pitzer mixing ionic parameters are determined and used to predict the solute activity coefficients in the mixture for different ionic-strength fractions.     

Mechanism and control of molecular energy flow: a modeling perspective

 Vibrational energy flow in organic molecules occurs by a multiple-time-scale mechanism that can be modeled by a single exponential only in its initial stages. The mechanism is a consequence of the hierarchical structure of the vibrational Hamiltonian, which leads to diffusion of vibrational wavepackets on a manifold with far fewer than the 3N−6 dimensions of the full vibrational state space. The dynamics are controlled by a local density of states, which does not keep increasing with molecular size. In addition, the number of vibrational coordinates severely perturbed during chemical reaction is small, leading to preservation of the hierarchical structure at chemically interesting energies. This regularity opens up the possibility of controlling chemical reactions by controlling the vibrational energy flow. Computationally, laser control of intramolecular vibrational energy redistribution can be modeled by quantum-classical, or by purely quantum-mechanical models of the molecule and control field. Received: 26 July 2002 / Accepted: 30 September 2002 / Published online: 2 December 2002 Electronic Supplementary Material to this paper can be obtained by using the Springer Link server located at http://dx.doi.org/10.1007/s00214-002-0394-2. Acknowledgements. This work was supported by NSF grant CHE 9986670. Correspondence to: M. Gruebele e-mail: gruebele@scs.uiuc.edu     

Gas–liquid partition coefficients and Henry's law constants of methyl mercaptan in aqueous solutions of Fe(II)–CDTA chelate complex

Utilization of ferric chelate complex of trans-1,2-cyclohexanediaminetetraacetic acid (CDTA) for the oxidative scrubbing of H₂S and CH₃SH in Kraft mill streams can be beneficial from the standpoints of iron protection against precipitation and oxygen-mediated regenerative oxidation of the ferrous chelate CDTA. The physical solubility of methyl mercaptan in CDTA–Fe(III) complex cannot be measured directly because of oxidation of the sulfur-bearing gaseous species with the ferric chelate. Therefore, this investigation was carried out to determine the gas–liquid partition coefficients and Henry's law constants of methyl mercaptan in aqueous iron-free CDTA solutions and non-reacting ferrous chelate solutions (CDTA–Fe(II) complex), using the static headspace method with an estimated accuracy of about 2%. Experiments with aqueous solutions of chelate concentrations varying between 38 and 300 mol m⁻³ were carried out at temperatures between 298 and 333 K and atmospheric pressure. It was shown that the methyl mercaptan solubility decreases with increasing temperature for all systems but was not much influenced, in the studied conditions, by the chelate concentration especially at larger temperatures.