Computer simulation on electron scavenging in irradiated hydrocarbon liquids

A Monte Carlo simulation method is applied to the calculations of the electron scavenging probability in irradiated hydrocarbon liquids. The simulation runs are carried out for different types of ionization clusters that form a high-energy electron track. The total scavenging probability is calculated by averaging the individual scavenging rates over the distribution of all types of structures that form the track. The scavenging probability for the track is lower by 20–30% than the probability calculated on the basis of the independent pair theory. The decrease is caused by the scavenger depletion in the multi-pair ionization clusters.     

Structure and thermodynamics of polymeric liquids: Polymer kirkwood integral equation method

A statistical mechanical theory is employed to predict the structural and thermodynamic properties of polymeric liquids. The theory consists of a coupled set of self-consistent integral equations for the pair correlation functions which relates the polymer structure to molecular interaction potentials. In this article, the essentials of the method are reviewed and the results are shown in comparison with other theories and simulation data. Quantitative agreements are found between the results of the theory and those of the simulations for a wide range of density and temperature.     

Applications of polymer rism theory to blends

Polymer RISM (reference interaction site model) theory is a theory of polymer systems in the liquid phase in which account for chemical realism can be made. Results are reported here of phase diagrams (spinodals) for blend systems calculated by means of this theory, using the mean spherical approximation as a closure. The systems investigated are an isotopic blend, a set of homopolymer/copolymer blends, and a model blend containing specific interactions.     

Computer simulation of reactive liquids in chemical equilibrium

A grand canonical ensemble Monte Carlo method applicable to reactive liquid mixtures in chemical equilibrium is described. The results of preliminary calculations in which the Br₂, Cl₂, BrCl equilibrium system is modeled using the new technique are reported.     

Computer simulation of quenched liquids: the glassy state of water

Monte Carlo simulation of liquid water quenched to low temperatures has been carried out in the NPT ensemble. The quenched system shows a glass-like transition and a distribution of hydrogen-bonded species different from the liquid.     

Computer simulation of ion recombination in irradiated nonpolar liquids

A review of the results of computer simulation of the diffusion controlled recombination of ions is presented. The ions generated in clusters of two and three pairs of oppositely charged ions were considered. The recombination kinetics and the ion escape probability at infinite time with and without external electric field have been computed. These results are compared with the calculations based on the single-pair theory.     

Local structure and thermodynamics of a core-softened potential fluid: theory and simulation.

Phase behavior and structural properties of homogeneous and inhomogeneous core-softened (CS) fluid consisting of particles interacting via the potential, which combines the hard-core repulsion and double attractive well interaction, are investigated. The vapour-liquid coexistence curves and critical points for various interaction ranges of the potential are determined by discrete molecular dynamics simulations to provide guidance for the choice of the bulk density and potential parameters for the study of homogeneous and inhomogeneous structures. Spatial correlations in the homogeneous CS system are studied by the Ornstein-Zernike integral equation in combination with the modified hypernetted chain (MHNC) approximation. The local structure of CS fluid subjected to diverse external fields maintaining the equilibrium with the bulk CS fluid are studied on the basis of a recently proposed third order+second order perturbation density functional approximation (DFA). The accuracy of DFA predictions is tested against the results of a grand canonical ensemble Monte Carlo simulation. Reasonable agreement between the results of both methods proves that the DFA theory applied in this work is a convenient theoretical tool for the investigation of the CS fluid, which is practically applicable for modeling numerous real systems.     

Complexes based on rigid-chain polyelectrolytes: Computer simulation

The structure and stability of complexes formed by oppositely charged rigid-chain macromolecules, as well as their response to variation in external conditions, were studied using the molecular dynamics method. The conditions of conformational transitions depending on the chain length and configuration, temperature, and dielectric permittivity of the medium were considered. It was shown that the chains involved in a complex may take various conformations of the torus, tennis racket, etc., types.     

Microscopic theory of orientational order, structure and thermodynamics in strained polymer liquids and networks

A microscopic integral equation theory of the segmental orientational order parameter, structural correlations and thermodynamics of strained polymer solutions, melts and networks has been developed. The nonclassical problem of the consequences of intermolecular excluded volume repulsions and chain connectivity is addressed. The theory makes several novel predictions, including effective power law dependences of the orientational order parameter on monomer concentration and chain degree of polymerization, and strain hardening of the bulk modulus. The predictions of a nearly classical strain dependence, and supralinear scaling with segment concentration, of the strain-induced nematic order parameter is in agreement with nuclear magnetic resonance experiments. The absolute magnitudes of the a priori calculated orientational order parameter agree with simulations and experiments to within a factor of 2. The possible complicating influence of "trapped entanglements" in crosslinked networks is discussed. Extensions of the theory are possible to treat the mechanical response of flexible polymer liquids and rubbers, and the structure, thermodynamics, and mechanical properties of strained liquid crystal forming polymers.     

Solvation thermodynamics: theory and applications

Potential distribution and coupling parameter theories are combined to interrelate previous solvation thermodynamic results and derive several new expressions for the solvent reorganization energy at both constant volume and constant pressure. We further demonstrate that the usual decomposition of the chemical potential into noncompensating energetic and entropic contributions may be extended to obtain a Gaussian fluctuation approximation for the chemical potential plus an exact cumulant expansion for the remainder. These exact expressions are further related to approximate first-order thermodynamic perturbation theory predictions and used to obtain a coupling-parameter integral expression for the sum of all higher-order terms in the perturbation series. The results are compared with the experimental global solvation thermodynamic functions for xenon dissolved in n-hexane and water (under ambient conditions). These comparisons imply that the constant-volume solvent reorganization energy has a magnitude of at most approximately kT in both experimental solutions. The results are used to extract numerical values of the solute-solvent mean interaction energy and associated fluctuation entropy directly from experimental solvation thermodynamic measurements.     

Computer simulation of random polymer networks: Structure and properties

An original rigorous computer model of a crosslinked polymer network simulating vulcanized rubbers has been developed. Owing to the parallel code and to dissipative particle dynamics, a modern simulation technique, it has become possible not only to describe in detail the microstructure of the resulting random network of polymer chains with different initial lengths but also to obtain well-reproducible stress-strain relationships in a wide range of λ values.     

Structure of polydisperse inverse ferrofluids: theory and computer simulation

By using theoretical analysis and molecular dynamics simulations, we investigate the structure of colloidal crystals formed by nonmagnetic microparticles (or magnetic holes) suspended in ferrofluids (called inverse ferrofluids), by taking into account the effect of polydispersity in size of the nonmagnetic microparticles. Such polydispersity often exists in real situations. We obtain an analytical expression for the interaction energy of monodisperse, bidisperse, and polydisperse inverse ferrofluids. Body-centered tetragonal (bct) lattices are shown to possess the lowest energy when compared with other sorts of lattices and thus serve as the ground state of the systems. Also, the effect of microparticle size distributions (namely, polydispersity in size) plays an important role in the formation of various kinds of structural configurations. Thus, it seems possible to fabricate colloidal crystals by choosing appropriate polydispersity in size.     

Fragility and thermodynamics in nonpolymeric glass-forming liquids

For nonpolymeric supercooled liquids, the empirical correlation m = 56Tg DeltaCp(Tg)/DeltaHm provides a reliable means of correlating dynamic and thermodynamic variables. The dynamics are characterized by the fragility or steepness index m and the glass transition temperature Tg, while thermodynamics enter in terms of the heat capacity step DeltaCp at Tg and the melting enthalpy DeltaHm. The combination of the above correlation with the 23 rule for the Tg/Tm ratio yields an expression, m = 40DeltaCp(Tg)/DeltaSm, which was rationalized as the correlation of the thermodynamic and kinetic fragilities. Defining a thermodynamic fragility via DeltaCp(Tg)/DeltaSm also reveals that the slopes in Kauzmann's original DeltaS(T)/DeltaSm versus T/Tm plot reflect the fragility concept [Chem. Rev. 43, 219 (1948)], so long as Tm/Tg = 1.5. For the many liquids whose excess heat capacity is a hyperbolic function of temperature, we deduce that the fragility cannot exceed m = 170, unless the Tg/Tm = 2/3 rule breaks down.     

Structure and dynamics of surfactant and hydrocarbon aggregates on graphite: a molecular dynamics simulation study

We have examined the structure and dynamics of sodium dodecyl sulfate (SDS) and dodecane (C12) molecular aggregates at varying surface coverages on the basal plane of graphite via classical molecular dynamics simulations. Our results suggest that graphite-hydrocarbon chain interactions favor specific molecular orientations at the single-molecule level via alignment of the tail along the crystallographic directions. This orientational bias is reduced greatly upon increasing the surface coverage for both molecules due to intermolecular interactions, leading to very weak bias at intermediate surface coverages. Interestingly, for complete monolayers, we find a re-emergent orientational bias. Furthermore, by comparing the SDS behavior with C12, we demonstrate that the charged head group plays a key role in the aggregate structures: SDS molecules display a tendency to form linear file-like aggregates while C12 forms tightly bound planar ones. The observed orientational bias for SDS molecules is in agreement with experimental observations of hemimicelle orientation and provides support for the belief that an initial oriented layer governs the orientation of hemimicellar aggregates.     

Scaled particle theory study of the length scale dependence of cavity thermodynamics in different liquids

It has been noted that the work of cavity creation in water exhibits a crossover behavior, in that its cavity size dependence changes from volume dependence for small cavities to area dependence for larger cavities [Lum, K.; Chandler, D.; Weeks, J. D. J. Phys. Chem. B 1999, 103, 4570]. It is shown here that this behavior can be reproduced using the scaled particle theory in a straightforward manner for six different liquids (water, methanol, ethanol, benzene, cyclohexane, and carbon tetrachloride). It has also been suggested that the crossover is due to a change in the physical mechanism of the process, from one entropy-dominated to another enthalpy-dominated. However, the crossover behavior can be produced using the scaled particle theory without invoking any change in any physical mechanism. Also, the crossover occurs at a length scale of the size of the liquid molecules, as has been pointed out by others. This is the length regime where the work of cavity creation bears little relation to the bulk liquid surface tension. In addition, it is pointed out that cavity creation can always be considered as a purely entropy-driven process, which is usually accompanied by another process with compensating enthalpy and entropy changes.     

Computer simulation study of rotational diffusion in polar liquids of different types

Rotational diffusion in liquid acetonitrile, dimethylsulphoxide (DMSO), water, and methanol is studied with molecular dynamics simulations. The effects of hydrogen bonding and local dipole-dipole correlations (Kirkwood g-factor) on the relationship between the single molecule and collective relaxation are examined. The first rank single molecule dipole moment autocorrelation functions (ACFs) are constructed in the molecule-fixed coordinate frame and the principal components of rotation diffusion tensor are reported. Higher rank orientational ACFs are computed. These ACFs, as a rule, are strongly nonexponential (at least not single exponential) at longer times and the decomposition of these functions into a series of single exponentials results in broad distributions of relaxation times, with the broadening being particularly prominent in the case of higher rank ACFs. The rank dependence of characteristic times calculated as weighted averages over the relaxation time distributions does not follow the pattern of small angle (Debye) diffusion model for all liquids studied in this work except methanol. In contradiction, the same rank dependence computed by direct integration of ACFs leads to good agreement with the Debye diffusion model in the case of acetonitrile, DMSO, and water (but not methanol). The linear-angular momentum cross correlation functions are also computed and the effect of rototranslational coupling on reorientaional relaxation at longer times (>1.0?ps) is found to be small.