Deviation from Corresponding States for a Fluid of Square Well Spherocylinders

Liquid crystal phase behaviour of attractive disc-like particles

We employ a generalized van der Waals-Onsager perturbation theory to construct a free energy functional capable of describing the thermodynamic properties and orientational order of the isotropic and nematic phases of attractive disc particles. The model mesogen is a hard (purely repulsive) cylindrical disc particle decorated with an anisotropic square-well attractive potential placed at the centre of mass. Even for isotropic attractive interactions, the resulting overall inter-particle potential is anisotropic, due to the orientation-dependent excluded volume of the underlying hard core. An algebraic equation of state for attractive disc particles is developed by adopting the Onsager trial function to characterize the orientational order in the nematic phase. The theory is then used to represent the fluid-phase behaviour (vapour-liquid, isotropic-nematic, and nematic-nematic) of the oblate attractive particles for varying values of the molecular aspect ratio and parameters of the attractive potential. When compared to the phase diagram of their athermal analogues, it is seen that the addition of an attractive interaction facilitates the formation of orientationally-ordered phases. Most interestingly, for certain aspect ratios, a coexistence between two anisotropic nematic phases is exhibited by the attractive disc-like fluids.

A new insight on the structural changes of linear quadrupole liquids

Molecular-dynamics simulations for linear quadrupole liquids are presented. The study is carried out for two different molecular lengths at constant density and a number of temperatures and quadrupole moments. All the simulated thermodynamic states correspond to the condensed phases and some of them show typical features of a solid structure. Furthermore, a change on the preferred intermolecular orientation in the liquid phase is observed from a shifted parallel molecular arrangement to a perpendicular orientation as the quadrupole raises. This change depends on the quadrupole moment as well as on the molecular length and is put in relation with the solid structure of different "diatomic" molecules such as nitrogen, ethane, and acetylene. The appearance of a plastic solid phase at low quadrupole moment and density is also justified. A thoroughly discussion about the availability of classical perturbation theories for this kind of systems is presented.

Smectic phase in a system of hard ellipsoids with isotropic attractive interactions

The smectic phase is studied for a thermotropic fluid model consisting of aligned hard ellipsoids with superimposed square-well attractive interactions of variable range. The system is analyzed using a density functional theory in which the hard-core contributions to the free-energy functional are treated within a nonlocal weighted density approximation and the attractive contributions are considered at a mean-field level. In the absence of attractions the model reduces, under appropriate scaling, to a fluid of hard spheres and therefore does not exhibit smectic ordering. It is shown that above a certain value of the square-well range, smectic ordering is stable relative to the nematic state at densities well inside the fluid region. The nematic-smectic-A transition is found to be continuous at high temperatures and first order at low temperatures, these two regimes being separated by a tricritical point at an intermediate temperature. These predictions have been confirmed by computer simulation of the model fluid. The results highlight that smectic ordering can be stabilized by coupling anisotropic short-range repulsions with the isotropic contribution of the soft attractive interactions. By increasing the pressure, the range of stability of the smectic phase is seen to decrease. At sufficiently high pressure, the smectic phase is suppressed, and the solid phase dominates. Our calculations show that smectic ordering is no longer stable if the range of the attractions is made too long ranged.

Solution of the Percus-Yevick equation for square well spherocylinders

The Percus-Yevick equation for square-well spherocylinders has been numerically solved for some selected orientations following a methodology proposed previously for different fluids of elongated molecules. The equation is solved for particles of aspect ratios ranging from L/sigma=0.3 up to L/sigma=5.0, attractive range lambda/sigma=1.5, and packing fractions within eta=0.1-0.3. The resulting pair correlation functions are checked against isothermal-isobaric Monte Carlo simulations and good agreement is found for the short-range structure, at intermolecular distances within one molecular diameter sigma to contact for each of the selected orientations. At larger distances, the integral equation tends to overestimate the pair correlations. The results confirm the prediction of reference-system average Mayer-function perturbation theory for short aspect ratios, reaching the Onsager limit for the greater aspect ratios. Some instabilities of the solution for the longest models and higher densities are tentatively discussed in terms of their possible relation to frustration phenomena found in some polymer and complex systems.