We investigated the nematic to smectic transition undergone by parallel hard spherocylinders in the framework provided by the residual multi-particle entropy (RMPE) formalism. The RMPE is defined as the sum of all contributions to the configurational entropy of the fluid which arise from density correlations involving more than two particles. The vanishing of the RMPE signals the structural changes which take place in the system for increasing pressures. Monte Carlo simulations carried out for parallel hard spherocylinders show that such a one-phase ordering criterion accurately predicts also the nematic-smectic transition threshold notwithstanding the almost continuous character of the transition. A similar quantitative correspondence had been already noted in the case of an isotropic fluid of freely rotating hard spherocylinders undergoing a transition to a nematic, smectic or solid phase. The present analysis confirms the flexibility of the RMPE approach as a practical and reliable tool for detecting the formation of mesophases in model liquid-crystal systems.
Hard spherocylinders (cylinders of length $L$ and diameter $D$ capped at both ends with two hemispheres) provide a suitable model for investigating entropy-driven, mesophase formations in real colloidal fluids that are composed of rigid rodlike molecules. We performed extensive Monte Carlo simulations of this model fluid for elongations in the range $3 leq L/D leq 5$ and up to $L/D = 20$, in order to investigate the relative importance of translational and orientational correlations allowing for the emergence of nematic or smectic order in the framework of the so-called residual multi-particle entropy (RMPE). The vanishing of this quantity, which includes the re-summed contributions of all spatial correlations involving more than two particles, signals the structural changes which take place, at increasing densities, in the isotropic fluid. We found that the ordering thresholds detected through the zero-RMPE condition systematically correlate with the corresponding phase-transition points, whatever the nature of the higher-density phase coexisting with the isotropic fluid.
Using computer simulations we investigate the microscopic structure of the singular director field within a nematic droplet. As a theoretical model for nematic liquid crystals we take hard spherocylinders. To induce an overall topological charge, the particles are either confined to a two-dimensional circular cavity with homeotropic boundary or to the surface of a three-dimensional sphere. Both systems exhibit half-integer topological point defects. The isotropic defect core has a radius of the order of one particle length and is surrounded by free-standing density oscillations. The effective interaction between two defects is investigated. All results should be experimentally observable in thin sheets of colloidal liquid crystals.
A new Monte Carlo approach is proposed to investigate the fluid-solid phase transition of the polydisperse system. By using the extended ensemble, a reversible path was constructed to link the monodisperse and corresponding polydisperse system. Once the fluid-solid coexistence point of the monodisperse system is known, the fluid-solid coexistence point of the polydisperse system can be obtained from the simulation. The validity of the method is checked by the simulation of the fluid-solid phase transition of a size-polydisperse hard sphere colloid. The results are in agreement with the previous studies.
We design generative neural networks that generate Monte Carlo configurations with complete absence of autocorrelation and from which direct measurements of physical observables can be employed, irrespective of the system locating at the classical critical point, fermionic Mott insulator, Dirac semimetal and quantum critical point. We further propose a generic parallel-chain Monte Carlo scheme based on such neural networks, which provides independent samplings and accelerates the Monte Carlo simulations by reducing the thermalization process. We demonstrate the performance of our approach on the two-dimensional Ising and fermion Hubbard models.
Modification of the hard jet substructure in terms of the Soft Drop jet grooming algorithm observables is studied for three different scenarios of jet quenching in a quark-gluon plasma: i) an explicit enhancement of the parton splitting functions, ii) increased soft gluon emissions induced by an in-medium virtuality gain, and iii) energy loss due to a drag force. Despite the fact that first two scenarios both correspond to a radiative energy loss mechanism and lead to similar modifications of parton showers, they are shown to have very different impacts on the momentum balance of hard subjets. Simulations for heavy-ion collisions based on the second scenario are presented and found to be in a good agreement with the experimental data.