Liquid-vapor interfaces in XY -spin fluids: an inhomogeneous anisotropic integral-equation approach

The BKT transition and its dynamics in a spin fluid

We study the effect of particle mobility on phase transitions in a spin fluid in two dimensions. The presence of a phase transition of the BKT universality class is shown in an off-lattice model of particles with purely repulsive interaction employing computer simulations. A critical spin wave region 0 < T < T_BKT is found with a nonuniversal exponent η(T) that follows the shape suggested by BKT theory, including a critical value consistent with η_BKT = 1/4. One can observe a transition from power-law decay to exponential decay in the static correlation functions at the transition temperature T_BKT, which is supported by finite-size scaling analysis. A critical temperature T_BKT = 0.17(1) is suggested. Investigations into the dynamic aspects of the phase transition are carried out. The short-time behavior of the incoherent spin autocorrelation function agrees with the Nelson-Fisher prediction, whereas the long-time behavior differs from the finite-size scaling known for the static XY model. Analysis of coherent spin wave dynamics shows that the spin wave peak is a propagating mode that can be reasonably well fitted by hydrodynamic theory. The mobility of the particles strongly enhances damping of the spin waves, but the model still lies within the dynamic universality class of the standard XY model.

Ferromagnetism-induced phase separation in a two-dimensional spin fluid

We study the liquid-gas phase separation observed in a system of repulsive particles dressed with ferromagnetically aligning spins, a so-called "spin fluid." Microcanonical ensemble numerical simulations of finite-size systems reveal that magnetization sets in and induces a liquid-gas phase separation between a disordered gas and a ferromagnetic dense phase at low enough energies and large enough densities. The dynamics after a quench into the coexistence region show that the order parameter associated with the liquid-vapor phase separation follows an algebraic law with an unusual exponent, as it is forced to synchronize with the growth of the magnetization: this suggests that for finite size systems the magnetization sets in along a Curie line, which is also the gas-side spinodal line, and that the coexistence region ends at a tricritical point. This picture is confirmed at the mean-field level with different approximation schemes, namely, a Bethe lattice resolution and a virial expansion complemented by the introduction of a self-consistent Weiss-like molecular field. However, a detailed finite-size scaling analysis shows that in two dimensions the ferromagnetic phase escapes the Berezinskii-Kosterlitz-Thouless scenario and that the long-range order is not destroyed by the unbinding of topological defects. The Curie line thus becomes a magnetic crossover in the thermodynamic limit. Finally, the effects of the magnetic interaction range and those of the interaction softness are characterized within a mean-field semianalytical low-density approach.

Phase Behavior of Magnetic Nanocolloids of Different Sizes Suspended in an Apolar Solvent

We employ classical density functional theory (DFT) to investigate the phase behavior and composition of binary mixtures; each compound consists of hard spheres of different sizes with superimposed dispersion attraction. In addition to the dispersion attraction, molecules of one component carry an additional three-dimensional magnetic "spin" where the orientation-dependent spin-spin interaction is accounted for by the Heisenberg model. We are treating the excess free energy using a modified mean-field approximation (second virial coefficient) for the orientation-dependent pair correlation function. Depending on the concentration of the magnetic particles, the strength of the spin-spin coupling, and the size ratio of the particles, the model predicts the formation of ordered (polar) phases in addition to the more conventional gas and isotropic liquid phases. Key features of our model are a particle-size dependent shift of the gas-liquid critical point (critical temperature and density) and a change in the width of the phase diagram. In the near-critical region, the latter can be analyzed quantitativly in terms of an effective critical exponent β_eff that may differ from the classical critical exponent [Formula: see text]; the classical value is attained in the immediate vicinity of the critical point as it must. The deviation between β_eff and β can be linked to nontrivial composition effects along the phase boundaries.

A closure relation to molecular theory of solvation for macromolecules

We propose a closure to the integral equations of molecular theory of solvation, particularly suitable for polar and charged macromolecules in electrolyte solution. This includes such systems as oligomeric polyelectrolytes at a finite concentration in aqueous and various non-aqueous solutions, as well as drug-like compounds in solution. The new closure by Kobryn, Gusarov, and Kovalenko (KGK closure) imposes the mean spherical approximation (MSA) almost everywhere in the solvation shell but levels out the density distribution function to zero (with the continuity at joint boundaries) inside the repulsive core and in the spatial regions of strong density depletion emerging due to molecular associative interactions. Similarly to MSA, the KGK closure reduces the problem to a linear equation for the direct correlation function which is predefined analytically on most of the solvation shells and has to be determined numerically on a relatively small (three-dimensional) domain of strong depletion, typically within the repulsive core. The KGK closure leads to the solvation free energy in the form of the Gaussian fluctuation (GF) functional. We first test the performance of the KGK closure coupled to the reference interaction site model (RISM) integral equations on the examples of Lennard-Jones liquids, polar and nonpolar molecular solvents, including water, and aqueous solutions of simple ions. The solvation structure, solvation chemical potential, and compressibility obtained from RISM with the KGK closure favorably compare to the results of the hypernetted chain (HNC) and Kovalenko-Hirata (KH) closures, including their combination with the GF solvation free energy. We then use the KGK closure coupled to RISM to obtain the solvation structure and thermodynamics of oligomeric polyelectrolytes and drug-like compounds at a finite concentration in electrolyte solution, for which no convergence is obtained with other closures. For comparison, we calculate their solvation structure from molecular dynamics (MD) simulations. We further couple the 3D-RISM integral equation with the 3D-version of the KGK closure, and solve it for molecular mixtures as well as oligomeric polyelectrolytes and drug-like molecules in electrolyte solutions.

Condensation, demixing, and orientational ordering of magnetic colloidal suspensions

In this work we study the phase behavior of magnetic particles suspended in a simple nonmagnetic solvent. Magnetic particles are modelled as spherical particles carrying a three-dimensional, classical Heisenberg spin, whereas solvent molecules are treated as spherically symmetric Lennard-Jones particles. The binary mixture of magnetic particles and solvent is studied within the framework of classical density functional theory (DFT). Within DFT pair correlations are treated at the modified mean-field level at which they are approximated by orientation dependent Mayer f functions. In the absence of an external magnetic field four generic types of phase diagrams are observed depending on the concentration of magnetic particles. In this case we observe liquid-liquid phase coexistence between an orientationally ordered (polarized) and a disordered phase characterized by slightly different concentrations of magnetic particles. Liquid-liquid phase coexistence is suppressed by an external field and vanishes completely if the strength of the field is sufficiently large.