Chain and ring structures in smectic phases of molecules with transverse dipoles

Comparative Study of the Optical and Dielectric Anisotropy of a Difluoroterphenyl Dimer and Trimer Forming Two Nematic Phases

We present a comparative study of the optical and dielectric anisotropy of a laterally fluorinated liquid crystal dimer and its homologous trimer, both exhibiting two nematic phases. In the high-temperature nematic phase, both oligomers exhibit positive optical anisotropy with similar magnitude, which, however, is lower in comparison with the optical anisotropy of the monomer. In the same temperature range, the dielectric permittivity along and perpendicular to the nematic director, measured on magnetically aligned samples, reveals negative dielectric anisotropy for both oligomers, which saturates as the temperature approaches the N-N phase transition temperature. Comparison of the dielectric anisotropies of the oligomers with the corresponding anisotropy of the monomer indicates a systematic variation of its magnitude with the number of the linked mesogenic units. Results are compared with the corresponding anisotropies of the cyanobiphenyl dimers, the archetypal compounds with two nematic phases, and are discussed in terms of the dipolar structure of the mesogens and the dipolar correlations in their nematic phases.

Pre-Smectic Ordering and the Unwinding Helix in Monte Carlo Simulations of Cholesteric Liquid-Crystals

In this paper, molecular chirality is studied for liquid-crystal fluids represented by hard rods with the addition of an attractive chiral dispersion term. Chiral forces between molecular pairs are assumed to be long-ranged and are described in terms of the pseudotensor of Goossens [W. J. A. Goossens, Mol. Cryst. Liq. Cryst. 1971, 12, 237-244]. Following Varga and Jackson [S. Varga and G. Jackson, Chem. Phys. Lett. 2003, 377, 6-12], this is combined with a hard-spherocylinder core. We investigate the relationship between molecular chirality and the helical pitch of the system, which occurs in the absence of full three-dimensional periodic boundary conditions. The dependence of the wavenumber of this pitch on the thermodynamic variables, temperature, and density is measured. We also explore the use of a novel surface boundary interaction model. As a result of this approach, we are able to lower the temperature of the system without the occurrence of nematic droplets, which would interfere with the formation of a uniaxial pitch. Regarding the theoretical predictions of Wensink and Jackson [H. H. Wensink and G. Jackson, J. Chem. Phys. 2009, 130, 234911], on the one hand, we have qualitative agreement with the observed non-monotonic density dependence of the wavenumber. Initially increasing with density, the wavenumber reaches a maximum, before falling as the density moves toward the point of phase transition from cholesteric to smectic. However, further analysis for shorter rods, in the presence of novel boundary conditions, reveals some disagreement with the theory, at least in this case; the unwinding of the cholesteric helix in the cholesteric phase occurs simultaneously with subtle increases in smectic ordering. These pre-smectic fluctuations have not been accounted for so far in theories on cholesterics but turn out to play a key role in controlling the pitch of cholesteric phases of rod-shaped mesogens with a small to moderate aspect ratio.

Nonconventional Phases of Colloidal Nanorods with a Soft Corona

Using computer simulations, we investigate the phase behavior of hard-core spherocylinders with a length-to-diameter ratio L/σ=5 and coated by a soft deformable corona of length λ/σ=1.35. When quasi-two-dimensional layers are formed in smectic and solid phases at low temperatures, the competition between the two intrinsic length scales of the parallel aligned particles leads to the stabilization of different in-plane lattices of nonconventional symmetry, including low-density hexagonal, square, and high-density hexagonal crystals, as well as an intriguing dodecagonal quasicrystal. Our Letter opens up the opportunity to control the assembly of anisotropic nanoparticles into structures with preengineered symmetry-dependent physical properties.

Machine learning for condensed matter physics

Condensed matter physics (CMP) seeks to understand the microscopic interactions of matter at the quantum and atomistic levels, and describes how these interactions result in both mesoscopic and macroscopic properties. CMP overlaps with many other important branches of science, such as chemistry, materials science, statistical physics, and high-performance computing. With the advancements in modern machine learning (ML) technology, a keen interest in applying these algorithms to further CMP research has created a compelling new area of research at the intersection of both fields. In this review, we aim to explore the main areas within CMP, which have successfully applied ML techniques to further research, such as the description and use of ML schemes for potential energy surfaces, the characterization of topological phases of matter in lattice systems, the prediction of phase transitions in off-lattice and atomistic simulations, the interpretation of ML theories with physics-inspired frameworks and the enhancement of simulation methods with ML algorithms. We also discuss in detail the main challenges and drawbacks of using ML methods on CMP problems, as well as some perspectives for future developments.

Morphological and crystallographic orientation of hematite spindles in an applied magnetic field

The magnetic response of spindle-shaped hematite (α-Fe2O3) nanoparticles was investigated by simultaneous small-angle and wide-angle X-ray scattering (SAXS/WAXS) experiments. The field-dependent magnetic and nematic order parameters of the magnetic single-domain nanospindles in a static magnetic field are fully described by SAXS simulations of an oriented ellipsoid with the implemented Langevin function. The experimental scattering intensities of the spindle-like particles can be modeled simply by using the geometrical (length, radius, size distribution) and magnetic parameters (strength of magnetic field, magnetic moment) obtained from isotropic SAXS and macroscopic magnetization measurements, respectively. Whereas SAXS gives information on the morphological particle orientation in the applied field, WAXS texture analysis elucidates the atomic scale orientation of the magnetic easy direction in the hematite crystal structure. Our results strongly suggest the tendency for uniaxial anisotropy but indicate significant thermal fluctuations of the particle moments within the hematite basal plane.

Effect of Central Longitudinal Dipole Interactions on Chiral Liquid-Crystal Phases

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.

Monte Carlo computer simulation of sedimentation of charged hard spherocylinders

In this article we present a NVT Monte Carlo computer simulation study of sedimentation of an electroneutral mixture of oppositely charged hard spherocylinders (CHSC) with aspect ratio L/σ = 5, where L and σ are the length and diameter of the cylinder and hemispherical caps, respectively, for each particle. This system is an extension of the restricted primitive model for spherical particles, where L/σ = 0, and it is assumed that the ions are immersed in an structureless solvent, i.e., a continuum with dielectric constant D. The system consisted of N = 2000 particles and the Wolf method was implemented to handle the coulombic interactions of the inhomogeneous system. Results are presented for different values of the strength ratio between the gravitational and electrostatic interactions, Γ = (mgσ)/(e(2)/Dσ), where m is the mass per particle, e is the electron's charge and g is the gravitational acceleration value. A semi-infinite simulation cell was used with dimensions Lx ≈ Ly and Lz = 5Lx, where Lx, Ly, and Lz are the box dimensions in Cartesian coordinates, and the gravitational force acts along the z-direction. Sedimentation effects were studied by looking at every layer formed by the CHSC along the gravitational field. By increasing Γ, particles tend to get more packed at each layer and to arrange in local domains with an orientational ordering along two perpendicular axis, a feature not observed in the uncharged system with the same hard-body geometry. This type of arrangement, known as tetratic phase, has been observed in two-dimensional systems of hard-rectangles and rounded hard-squares. In this way, the coupling of gravitational and electric interactions in the CHSC system induces the arrangement of particles in layers, with the formation of quasi-two dimensional tetratic phases near the surface.

Ferroelectric order in liquid crystal phases of polar disk-shaped ellipsoids

The demonstration of a spontaneous macroscopic ferroelectric order in liquid phases in the absence of any long range positional order is considered an outstanding problem of both fundamental and technological interest. Recently, we reported that a system of polar achiral disklike ellipsoids can spontaneously exhibit a long searched ferroelectric nematic phase and a ferroelectric columnar phase with strong axial polarization. The major role is played by the dipolar interactions. The model system of interest consists of attractive-repulsive Gay-Berne oblate ellipsoids embedded with two parallel point dipoles positioned symmetrically on the equatorial plane of the ellipsoids. In the present work, we investigate in detail the profound effects of changing the separation between the two symmetrically placed dipoles and the strength of the dipoles upon the existence of different ferroelectric discotic liquid crystal phases via extensive off-lattice N-P-T Monte Carlo simulations. Ferroelectric biaxial phases are exhibited in addition to the uniaxial ferroelectric fluids where the phase biaxiality results from the dipolar interactions. The structures of all the ferroelectric configurations of interest are presented in detail. Simple phase diagrams are determined which include different polar and apolar discotic fluids generated by the system.

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.

Monte Carlo simulations of spontaneous ferroelectric order in discotic liquid crystals

The demonstration of a spontaneous macroscopic ferroelectric order in liquid phases in the absence of any long-range positional order is considered as an outstanding problem of great fundamental and technological interest. We report here off-lattice Monte Carlo simulations of a system of polar achiral disklike ellipsoids which spontaneously exhibit a novel ferroelectric nematic phase which is a liquid in three dimensions, considering attractive-repulsive pair interaction suitable for the anisotropic particles. At lower temperature, the ferroelectric nematic phase condenses to a ferroelectric hexagonal columnar fluid with an axial macroscopic polarization. A spontaneous ferroelectric order of dipolar origin is established here for the first time in columnar liquid crystals. Our study demonstrates that simple dipolar interactions are indeed sufficient to produce a class of novel ferroelectric fluids of essential interest. The present work reveals the structure-property relationship of achieving long searched ferroelectric liquid crystal phases and transitions between them, and we hope these findings will help in future development of technologically important fluid ferroelectric materials.

Smectic order parameters via liquid crystal NMR spectroscopy: Application to a partial bilayer smectic A phase

Solute molecules were dissolved in the liquid crystal 4-cyano-4'-n-octyloxybiphenyl (8OCB), known to form a partial bilayer smectic-A phase. Through measurement of solutes' and solvent's orientational order parameters via nuclear magnetic resonance spectroscopy, and their analysis via a statistical thermodynamic density functional theory, values of the solvent's positional order parameters and solutes' positional-orientational distribution functions were obtained. Near to the transition to the nematic phase, the main positional order parameter of the smectic liquid crystal turned out to be comprised in the interval 0.4-0.6, though the quality of the fittings assuming the phase as nematic all across the temperature range investigated was only slightly worse. This may be ascribed to the looseness of the partial bilayer smectic structure. Solutes were found to preferentially lie in those regions where liquid crystal molecule terminal chains are located.

Anomalous columnar order of charged colloidal platelets

Monte Carlo computer simulations are carried out for a model system of like-charged colloidal platelets in the isothermal-isobaric ensemble (NpT). The aim is to elucidate the role of electrostatic interactions on the structure of synthetic clay systems at high particle densities. Short-range repulsions between particles are described by a suitable hard-core model representing a discotic particle. This potential is supplemented with an electrostatic potential based on a Yukawa model for the screened Coulombic potential between infinitely thin disklike macro-ions. The particle aspect-ratio and electrostatic parameters were chosen to mimic an aqueous dispersion of thin, like-charged, rigid colloidal platelets at finite salt concentration. An examination of the fluid phase diagram reveals a marked shift in the isotropic-nematic transition compared to the hard cut-sphere reference system. Several statistical functions, such as the pair correlation function for the center-of-mass coordinates and structure factor, are obtained to characterize the structural organization of the platelets phases. At low salinity and high osmotic pressure we observe anomalous hexagonal columnar structures characterized by interpenetrating columns with a typical intercolumnar distance corresponding to about half of that of a regular columnar phase. Increasing the ionic strength leads to the formation of glassy, disordered structures consisting of compact clusters of platelets stacked into finite-sized columns. These so-called "nematic columnar" structures have been recently observed in systems of charge-stabilized gibbsite platelets. Our findings are corroborated by an analysis of the static structure factor from a simple density functional theory.

Cholesteric order in systems of helical Yukawa rods

We consider the interaction potential between two chiral rod-like colloids which consist of a thin cylindrical backbone decorated with a helical charge distribution on the cylinder surface. For sufficiently slender helical rods a simple scaling expression is derived which relates the chiral 'twisting' potential to the microscopic properties of the particles, such as the internal helical pitch, charge density and electrostatic screening parameter. To predict the behaviour of the macroscopic cholesteric pitch of the fluid bulk phase we invoke a simple second-virial theory generalized to treat anisotropic states with weakly twisted director fields. It is shown that, while particles with weakly coiled helices always form a cholesteric phase whose helical sense is commensurate with that of the internal helix, more strongly coiled rods lead to the formation of a cholesteric state of opposite sense. The correlation between the helical symmetry at the microscopic and macroscopic scale is found to be very sensitive to the pitch of the Yukawa helix. Mixing helical particles of sufficiently disparate length and internal pitch may give rise to a demixing of the uniform cholesteric phase into two fractions with a different macroscopic pitch. Our findings could be relevant to the interpretation of experimental observations in systems of cellulose and chitin microfibres, DNA and fd virus rods.

Crystal structure and microstructure of cholesteryl oleyl carbonate

The crystal structure as well as the microstructure, i.e., size and strain, of crystallites of cholesteryl oleyl carbonate was determined from X-ray powder diffraction data. The X-ray line broadening was analyzed through the refinement of TCH-pseudo-Voigt function parameters (isotropic effects) and the refinement of multipolar functions, i.e., symmetrized cubic harmonics (anisotropic effects). The crystal structure turns out to be primitive monoclinic, space group Pc, type I monolayer having two molecules per unit cell with parameters: a=18.921±0.006Å, b=12.952±0.003Å, c=9.276±0.002Å and β=91.32±0.03°. The average size of a well ground specimen of crystallites was 60nm. The average micro-strain, e.g., 45×10(-4) has been tentatively attributed to fatty chain conformational disorder. The unit cell parameters, including the lamellar thickness, of COC crystal is very closely similar to those of another, structurally similar cholesterol ester, e.g., cholesteryl oleate (CO) crystal, space group P2(1), type II monolayer. Type I monolayer structure has been established for COC on the basis of the intensity calculations of the XRD profiles of both CO and COC. The dipolar and structural disorder in a 4:1 molar, binary mixture of CO and COC can be accommodated in an induced smectic phase with a lamellar thickness, which is nearly equal to that of pure CO or pure COC.

Predicting the Phase Diagram of Two-Dimensional Colloidal Systems with Long-Range Interactions

The phase diagram of a two-dimensional model system for colloidal particles at the air-water interface was determined using Monte Carlo computer simulations in the isothermic-isobaric ensemble. The micrometer-range binary colloidal interaction has been modeled by hard disklike particles interacting via a secondary minimum followed by a weaker longer-range repulsive maximum, both of the order of kBT. The repulsive part of the potential drives the clustering of particles at low densities and low temperatures. Pinned voids are formed at higher densities and intermediate values of the surface pressure. The analysis of isotherms, translational and orientational correlation functions as well as structure factor gives clear evidence of the presence of a melting first-order transition. However, the melting process can be also followed by a metastable route through a hexatic phase at low surface pressures and low temperatures, before crystalization occurs at higher surface pressure.

Dipolar interactions, molecular flexibility, and flexoelectricity in bent-core liquid crystals

The effects of dipolar interactions and molecular flexibility on the structure and phase behavior of bent-core molecular fluids are studied using Monte Carlo computer simulations. Some calculations of flexoelectric coefficients are also reported. The rigid cores of the model molecules consist of either five or seven soft spheres arranged in a "V" shape with external bend angle gamma. With purely repulsive sphere-sphere interactions and gamma = 0 degrees (linear molecules) the seven-sphere model exhibits isotropic, uniaxial nematic, and untilted and tilted smectic phases. With gamma > or = 20 degrees the untilted smectic phases disappear, while the system with gamma > or = 40 degrees shows a direct tilted smectic-isotropic fluid transition. The addition of electrostatic interactions between transverse dipole moments on the apical spheres is generally seen to reduce the degree of molecular inclination in tilted phases, and destabilizes the nematic and untilted smectic phases of linear molecules. The effects of adding three-segment flexible tails to the ends of five-sphere bent-core molecules are examined using configurational-bias Monte Carlo simulations. Only isotropic and smectic phases are observed. On the one hand, molecular flexibility gives rise to pronounced fluctuations in the smectic-layer structure, bringing the simulated system in better correspondence with real materials; on the other hand, the smectic phase shows almost no tilt. Lastly, the flexoelectric coefficients of various nematic phases--with and without attractive sphere-sphere interactions--are presented. The results are encouraging, but a large computational effort is required to evaluate the appropriate fluctuation relations reliably.

Pattern formation in a monolayer of magnetic spheres

Pattern formation is investigated for a vertically vibrated monolayer of magnetic spheres. The spheres of diameter D encase cylindrical magnetic cores of length l. For large D/l, we find that the particles form a hexagonal-close-packed pattern in which the particles' dipole vectors assume a macroscopic circulating vortical pattern. For smaller D/l, the particles form concentric rings. The static configurational magnetic energy (which depends on D/l) appears to be a determining factor in pattern selection even though the experimental system is driven and dissipative.

Computer simulation of polar bent-core molecules

Results are presented from molecular dynamics simulations in the NPT ensemble of novel bent-core liquid crystal systems. Following on from a previous study of bent-core steric shape, this study examines the effect the addition of a transverse electric dipole has on the phase diagram of a bent-core liquid crystal model. A simple model of the interaction employed a two-site Gay-Berne potential with the sites separated by +/-0.5 reduced units with a central transverse point dipole, for all models investigated. The angle between the sites 180 degrees -gamma was varied in a range gamma=10 degrees to gamma=70 degrees suggested by real molecules. The addition of the dipole to the model tended to stabilize smectic phases and increase the angle of tilted phases. As the angle gamma increased, the transition temperature to the first ordered phase decreased markedly. Smectic A, tilted smectic B, and a spontaneously polarized smectic B phases were observed in the gamma=10 degrees bent-core model. The gamma=20 degrees model showed smectic A and tilted antiferroelectric smectic B phases. The gamma=40 degrees model showed an antiferroelectric phase that exhibited unusual packing behavior. Both the gamma=20 degrees and gamma=40 degrees models demonstrated a significant phase biaxiality in the smectic B phases.