A periodic reaction field (PRF) method is a technique to estimate long-range interactions. The method has the potential to effectively reduce the computational cost while maintaining adequate accuracy. We performed molecular dynamics (MD) simulations of a model liquid-crystal system to assess the accuracy of some variations of the PRF method in low-charge-density systems. All the methods had adequate accuracy compared with the results of the particle mesh Ewald (PME) method, except for a few simulation conditions. Furthermore, in all of the simulation conditions, one of the PRF methods had the same accuracy as the PME method.

In future large-scale molecular dynamics (MD) simulations that will use parallel computing, the isotropic periodic sum (IPS) method is expected to effectively reduce the cost of interaction calculations while maintaining adequate accuracy. To assess the accuracy of this method in estimating low-charge-density polymer systems, we performed atomistic MD simulations of the bulk state of liquid crystal systems based on 4-pentyl-4 0 -cyanobiphenyl (5CB). In conditions of 270 K # T # 320 K and a normal pressure, the temperature dependence of the density, potential energy and order parameter was estimated using the IPS and Ewald sum method. The results of the IPS method and Ewald sum were consistent within the range of error. In conditions close to the phase transition point, however, the averaged values of potential energy and order parameter had a small difference. We concluded that the fundamental physical properties for the bulk state of 5CB systems are determined reasonably by using the IPS method, at least in conditions that are not close to the phase transition point.

The nucleation process of anisotropic particles often differs from that of their spherically symmetric counterparts. Despite a large body of work on the structure of droplets of anisotropic particles, their formation process remains poorly understood. In this study, homogeneous nucleation of uniaxial anisotropic particles was studied. Through structural analysis of cluster development and the formation free energy during the nucleation stage, it was revealed that the nucleation of uniaxial particles begins from highly ordered states. There is, however, a marked decrease in orientational order within the cluster before critical nucleus size is attained. Further investigation on variations in the molecular interactions demonstrates how droplet elongation and the direction of the nematic ordering director relative to the axis of elongation can both be controlled according to the nature of the molecular anisotropy.

Truncation is still chosen for many long-range intermolecular interaction calculations to efficiently compute free-boundary systems, macromolecular systems and net-charge molecular systems, for example. Advanced truncation methods have been developed for long-range intermolecular interactions. Every truncation method can be implemented as one of two basic cut-off schemes, namely either an atom-based or a group-based cut-off scheme. The former computes interactions of “atoms” inside the cut-off radius, whereas the latter computes interactions of “molecules” inside the cut-off radius. In this work, the effect of group-based cut-off is investigated for isotropic periodic sum (IPS) techniques, which are promising cut-off treatments to attain advanced accuracy for many types of molecular system. The effect of group-based cut-off is clearly different from that of atom-based cut-off, and severe artefacts are observed in some cases. However, no severe discrepancy from the Ewald sum is observed with the extended IPS techniques.

Monte Carlo simulations of chiral liquid-crystals, represented by a simple coarse-grained chiral Gay–Berne model, were performed to investigate the effect of central longitudinal dipole interactions on phase behavior. A systematic analysis of the structural properties and phase behavior of both achiral and chiral systems, with dipole interactions, reveals differing effects; strong dipole interactions enhance the formation of layered structures; however, chiral interactions may prevent the formation of such phases under certain conditions. We also observed a short-ranged smectic structure within the cholesteric phases with strong dipole interactions. This constitutes possible evidence of presmectic ordering and/or the existence of chiral line liquid phases, which have previously been observed in X-ray experiments to occur between the smectic twisted grain boundary and cholesteric phases. These results provide a systematic understanding of how the phase behavior of chiral liquid-crystals changes when alterations are made to the strength of dipole interactions.