Computer simulations of adsorbed liquid crystal films

Ultraconfined oblate hard particles between hybrid penetrable walls

We have investigated, by Monte Carlo simulation, the orientational structure of very thin films of a discotic liquid crystal (DLC) confined between hybrid walls of controllable penetrability, as a function of wall separation L_{z}. Our purpose was to clarify whether, as predicted by continuum theory, the preferred orientation of the DLC is uniform, changes linearly, or changes discontinuously, when L_{z} and the anchoring strengths at either wall are changed. The model consists of oblate hard Gaussian overlap (HGO) particles: each wall sees a particle as a disk of zero thickness and diameter D less than or equal to that of the actual particle σ_{0}, embedded inside the particle and located halfway along, and perpendicular to, its minor axis. This provides a particle-level mechanism to control the anchoring properties of the walls, from planar (edge-on) for D∼0 to homeotropic (face-on) for D∼σ_{0}, which can be done independently at either wall. As in our earlier work [C. Anquetil-Deck et al., J. Phys. Chem. B 124, 7709 (2020)1520-610610.1021/acs.jpcb.0c05027], which was restricted to L_{z}=6σ_{0}, depending on the values of D_{s}≡D/σ_{0} at the top (D_{s}^{t}) and bottom (D_{s}^{b}) walls, we find domains in (D_{s}^{b},D_{s}^{t}) space in which particle alignment is uniform planar (UP), uniform homeotropic (UH), or varies linearly from planar at one wall to homeotropic at the other (L), but no bistable or tristable regions are identified between these domains. Most importantly, there appears never to occur an abrupt change of the LC orientation when the walls strongly favor different anchorings, in general agreement with the scenario proposed by Velasco and co-workers [D. de las Heras et al., Phys. Rev. E 79, 011712 (2009)1539-375510.1103/PhysRevE.79.011712], but in contrast to the behavior of equivalent calamitic systems [F. Barmes et al., Phys. Rev. E 69, 061705 (2004)1539-375510.1103/PhysRevE.69.061705; Phys. Rev. E, 71, 021705 (2005)1539-375510.1103/PhysRevE.71.021705; C. Anquetil-Deck et al., Phys. Rev. E 86, 041707 (2012)1539-375510.1103/PhysRevE.86.041707]. However, for the thinnest films investigated (L_{z}=2σ_{0}), the system is unable to accommodate a rotation of the preferred particle orientation from one wall to the other and adopts instead a tilted configuration, similar to that reported earlier for Gay-Berne films in symmetric confinement [T. Gruhn et al., Thin Solid Films 330, 46 (1998)0040-609010.1016/S0040-6090(98)00799-8; Mol. Phys. 93, 681 (1998)10.1080/002689798169014] but which, as far as we know, has been missed in most earlier work.

Structures, thermodynamics and dynamics of topological defects in Gay-Berne nematic liquid crystals

Topological defects are a ubiquitous phenomenon across different physical systems. A better understanding of defects can be helpful in elucidating the physical behaviors of many real materials systems. In nematic liquid crystals, defects exhibit unique optical signatures and can segregate impurities, showing their promise as molecular carriers and nano-reactors. Continuum theory and simulations have been successfully applied to link static and dynamical behaviors of topological defects to the material constants of the underlying nematic. However, further evidence and molecular details are still lacking. Here we perform molecular dynamics simulations of Gay-Berne particles, a model nematic, to examine the molecular structures and dynamics of +1/2 defects in a thin-film nematic. Specifically, we measure the bend-to-splay ratio K₃/K₁ using two independent, indirect measurements, showing good agreement. Next, we study the annihilation event of a pair of ±1/2 defects, of which the trajectories are consistent with experiments and hydrodynamic simulations. We further examine the thermodynamics of defect annihilation in an NVE ensemble, leading us to correctly estimate the elastic modulus by using the energy conservation law. Finally, we explore effects of defect annihilation in regions of nonuniform temperature within these coarse-grained molecular models which cannot be analysed by existing continuum level simulations. We find that +1/2 defects tend to move toward hotter areas and their change of speed in a temperature gradient can be quantitatively understood through a term derived from the temperature dependence of the elastic modulus. As such, our work has provided molecular insights into structures and dynamics of topological defects, presented unique and accessible methods to measure elastic constants by inspecting defects, and proposed an alternative control parameter of defects using temperature gradient.

A simple predictor of interface orientation of fluids of disk-like anisotropic particles and its implications for organic semiconductors

From classical molecular dynamics simulations, we identify a simple and general predictor of molecular orientation at solid and vapour interfaces of isotropic fluids of disk-like anisotropic particles based on their shape and interaction anisotropy. For a wide variety of inter-particle interactions, temperatures, and substrate types within the range of typical organic semiconductors and their processing conditions, we find remarkable universal scaling of the orientation at the interface with the free energy calculated from pair interactions between close-packed nearest neighbours and an empirically derived universal relationship between the entropy and the shape anisotropy and bulk volume fraction of the fluid particles. The face-on orientation of fluid particles at the solid interface is generally predicted to be the equilibrium structure, although the alignment can be controlled by tuning the particle shape and substrate type, while changing the strength of fluid-fluid interactions is likely to play a less effective role. At the vapour interface, only the side-on structure is predicted, and conditions for which the face-on structure may be preferred, such as low temperature, low interaction anisotropy, or low shape anisotropy, are likely to result in little orientation preference (due to the low anisotropy) or be associated with a phase transition to an anisotropic bulk phase for systems with interactions in the range of typical organic semiconductors. Based on these results, we propose a set of guidelines for the rational design and processing of organic semiconductors to achieve a target orientation at a solid or vapour interface.

Pressure-tensor method evaluation of the interfacial tension between Gay-Berne isotropic fluid and a smooth repulsive wall

The interfacial properties of a confined thermotropic liquid crystalline material are investigated using a molecular dynamics simulation technique. The pairwise interaction among the soft ellipsoidal particles is modeled by the Gay-Berne (GB) potential. The GB ellipsoids are confined by two soft, smooth, repulsive walls defined by the Weeks-Chandler-Andersen (WCA) potential. The aperiodic confinement due to walls makes the system mechanically anisotropic. Hence using the pressure-tensor method, the interfacial tension of an interface between the bulk isotropic (I) phase and WCA wall at various number densities (ρ) is calculated. From the pressure tensor and orientational order profiles, the arrangement of ellipsoids in the bulk and the vicinity of the wall is determined. The effect of system size and the wall-particle interaction strength (ε_W) on is also analyzed by varying the system size and ε_W.

The interplay of interfaces, supramolecular assembly, and electronics in organic semiconductors

Organic semiconductors, which include a diverse range of carbon-based small molecules and polymers with interesting optoelectronic properties, offer many advantages over conventional inorganic semiconductors such as silicon and are growing in importance in electronic applications. Although these materials are now the basis of a lucrative industry in electronic displays, many promising applications such as photovoltaics remain largely untapped. One major impediment to more rapid development and widespread adoption of organic semiconductor technologies is that device performance is not easily predicted from the chemical structure of the constituent molecules. Fundamentally, this is because organic semiconductor molecules, unlike inorganic materials, interact by weak non-covalent forces, resulting in significant structural disorder that can strongly impact electronic properties. Nevertheless, directional forces between generally anisotropic organic-semiconductor molecules, combined with translational symmetry breaking at interfaces, can be exploited to control supramolecular order and consequent electronic properties in these materials. This review surveys recent advances in understanding of supramolecular assembly at organic-semiconductor interfaces and its impact on device properties in a number of applications, including transistors, light-emitting diodes, and photovoltaics. Recent progress and challenges in computer simulations of supramolecular assembly and orientational anisotropy at these interfaces is also addressed.

Structure and Interfacial Tension of a Hard-Rod Fluid in Planar Confinement

The structural properties and interfacial tension of a fluid of rodlike hard-spherocylinder particles in contact with hard structureless flat walls are studied by means of Monte Carlo simulation. The calculated surface tension between the rod fluid and the substrate is characterized by a nonmonotonic trend as a function of the bulk concentration (density) over the range of isotropic bulk concentrations. As suggested by earlier theoretical studies, a surface-ordering scenario is confirmed by our simulations: the local orientational order close to the wall changes from uniaxial to biaxial nematic when the bulk concentration reaches about 85% of the value at the onset of the isotropic-nematic phase transition. The surface ordering coincides with a wetting transition whereby the hard wall is wetted by a nematic film. Accurate values of the fluid-solid surface tension, the adsorption, and the average particle-wall contact distance are reported (over a broad range of densities into the dense nematic region for the first time), which can serve as a useful benchmark for future theoretical and experimental studies on confined rod fluids. The simulation data are supplemented with predictions from second-virial density functional theory, which are in good qualitative agreement with the simulation results.

Rupture mechanism of liquid crystal thin films realized by large-scale molecular simulations

The ability of liquid crystal (LC) molecules to respond to changes in their environment makes them an interesting candidate for thin film applications, particularly in bio-sensing, bio-mimicking devices, and optics. Yet the understanding of the (in)stability of this family of thin films has been limited by the inherent challenges encountered by experiment and continuum models. Using unprecedented large-scale molecular dynamics (MD) simulations, we address the rupture origin of LC thin films wetting a solid substrate at length scales similar to those in experiment. Our simulations show the key signatures of spinodal instability in isotropic and nematic films on top of thermal nucleation, and importantly, for the first time, evidence of a common rupture mechanism independent of initial thickness and LC orientational ordering. We further demonstrate that the primary driving force for rupture is closely related to the tendency of the LC mesogens to recover their local environment in the bulk state. Our study not only provides new insights into the rupture mechanism of liquid crystal films, but also sets the stage for future investigations of thin film systems using peta-scale molecular dynamics simulations.

Density functional theory of liquid crystals and surface anchoring: hard Gaussian overlap-sphere and hard Gaussian overlap-surface potentials

This article applies the density functional theory to confined liquid crystals, comprised of ellipsoidal shaped particles interacting through the hard Gaussian overlap (HGO) potential. The extended restricted orientation model proposed by Moradi and co-workers [J. Phys.: Condens. Matter 17, 5625 (2005)] is used to study the surface anchoring. The excess free energy is calculated as a functional expansion of density around a reference homogeneous fluid. The pair direct correlation function (DCF) of a homogeneous HGO fluid is approximated, based on the optimized sum of Percus-Yevick and Roth DCF for hard spheres; the anisotropy introduced by means of the closest approach parameter, the expression proposed by Marko [Physica B 392, 242 (2007)] for DCF of HGO, and hard ellipsoids were used. In this study we extend an our previous work [Phys. Rev. E 72, 061706 (2005)] on the anchoring behavior of hard particle liquid crystal model, by studying the effect of changing the particle-substrate contact function instead of hard needle-wall potentials. We use the two particle-surface potentials: the HGO-sphere and the HGO-surface potentials. The average number density and order parameter profiles of a confined HGO fluid are obtained using the two particle-wall potentials. For bulk isotropic liquid, the results are in agreement with the Monte Carlo simulation of Barmes and Cleaver [Phys. Rev. E 71, 021705 (2005)]. Also, for the bulk nematic phase, the theory gives the correct density profile and order parameter between the walls.

Control of the nematic-isotropic phase transition by an electric field

We use a relatively simple continuum model to investigate the effects of dielectric inhomogeneity within confined liquid-crystal cells. Specifically, we consider, in planar, cylindrical, and spherical geometries, the stability of a nematic-isotropic interface subject to an applied voltage when the nematic liquid crystal has a positive dielectric anisotropy. Depending on the magnitude of this voltage, the temperature, and the geometry of the cell, the nematic region may shrink until the material is completely isotropic within the cell, grow until the nematic phase fills the cell, or, in certain geometries, coexist with the isotropic phase. For planar geometry, no coexistence is found, but we are able to give analytical expressions for the critical voltage for an electric-field-induced phase transition as well as the critical wetting layer thickness for arbitrary applied voltage. In cells with cylindrical and spherical geometries, however, locally stable nematic-isotropic coexistence is predicted, the thickness of the nematic region being controllable by alteration of the applied voltage.

Ordering of p-n-alkoxybenzoic acids at phase transition temperatures: a comparative computational analysis

A comparative analysis of molecular ordering of nematogenic p-n-alkoxybenzoic acids has been carried out with respect to translatory and orientational motions for the acids with seven (7OBAC), eight (8OBAC), nine (9OBAC) and 10 (10OBAC) carbon atoms in the alkyl chain. The CNDO/2 method has been used to compute the net atomic charge and dipole moment components at each atomic center. Modified Rayleigh-Schrodinger perturbation theory with multicentered-multipole expansion method has been used to evaluate long-range intermolecular interactions while a '6-exp' potential function has been assumed for short-range interactions. The total interaction-energy values obtained by these computations were used to calculate the probability of each configuration at the phase-transition temperature using the Maxwell-Boltzmann formula. The flexibility of various configurations has been studied in terms of variation of probability due to small departures from the most probable configuration. A comparative picture of molecular parameters like total energy, binding energy and total dipole moment has been given. An attempt has been made to explain the nematogenicity of these acids in terms of their relative order with the molecular parameter introduced in this paper.

A computational analysis of ordering in ABCHN at the nematic–isotropic transition

A computational analysis of ordering in the nematogenic compound 4-alkenyl bicyclohexylnitrile has been carried out based on quantum mechanics and intermolecular forces. The evaluation of atomic charge and dipole moment at each atomic center has been carried out using the complete neglect differential overlap (CNDO/2) method. Modified Rayleigh-Schrodinger perturbation theory along with a multicentered- multipole expansion method has been employed to evaluate long-range intermolecular interactions, while a '6-exp' potential function has been assumed for short-range interactions. The total interaction energy values obtained through these computations were used to calculate the probability of each configuration at room temperature (300 K), the nematic-isotropic transition temperature (364.7 K) and above transition temperature (450 K) using the Maxwell-Boltzmann formula. The various possible configurations during the different modes (i.e., stacking, in-plane and terminal) of interactions have been studied in terms of variation of probability due to small departures from the most probable configurations. An attempt has been made to analyze the characteristic features of liquid crystallinity in terms of their relative order with molecular parameters introduced in this paper.

Using particle shape to induce tilted and bistable liquid crystal anchoring

We use Monte Carlo simulations of hard Gaussian overlap (HGO) particles symmetrically confined in slab geometry to investigate the role of particle-substrate interactions on liquid crystalline anchoring. Despite the restriction here to purely steric interactions and smooth substrates, a range of behaviors are captured, including tilted anchoring and homeotropic-planar bistability. These macroscopic behaviors are all achieved through appropriate tuning of the microscopics of the HGO-substrate interaction, based upon nonadditive descriptions for the HGO-substrate shape parameter.

Surface ordering of diskotic liquid crystals

We present Monte Carlo simulations of diskotic molecules using the Gay-Berne potential in a slab geometry. The disk-wall interaction is described by two different functions according to whether or not the equilibrium distance is dependent on the relative orientation of the disk to the wall. Furthermore, by changing the parameters of these potentials, we model either homeotropic (face-on) or planar (edge-on) anchoring of the disks. We have found that the isotropic-nematic transition does not change in comparison with the bulk situation. The temperature of the nematic-columnar transition, on the contrary, is found to increase for homeotropic anchoring, and decrease for planar anchoring, independently of the details of the potential. We explain the decrease of the transition temperature in the planar anchoring situation as the result of an induced frustration, due to the competition between the two orientations induced independently by the upper and lower walls.

Computer simulation of a liquid-crystal anchoring transition

We present a study of the effects of confinement on a system of hard Gaussian overlap particles interacting with planar substrates through the hard-needle-wall potential. Using geometrical arguments to calculate the molecular volume absorbed at the substrates, we show that both planar and homeotropic arrangements can be obtained using this model. Monte Carlo simulations are then used to perform a systematic study of the model's behavior as a function of the system density and the hard-needle-wall interaction parameter. As well as showing the homeotropic to planar anchoring transition, the anchoring phase diagrams computed from these simulations indicate regions of bistability. This bistable behavior is examined further through the explicit simulation of field-induced two-way switching between the two arrangements.

The role of attractive interactions in rod–sphere mixtures

We present a computer simulation study of binary mixtures of prolate Gay-Berne particles and Lennard-Jones spheres. Results are presented for three such rod-sphere systems which differ from each other only in the interaction between unlike particles. Both the mixing-demixing behavior and the transitions between the isotropic and any liquid crystalline phases are studied for each system, as a function of temperature and concentration ratio. For systems which show macroscopic demixing, the rod-sphere interaction is shown to give direct control over interfacial anchoring properties, giving rise to the possibility of micellar phase formation in the case of homeotropic anchoring. Additionally, it is shown that on incorporating high concentrations of spheres into a system of rods with weak demixing properties, microphase-separated structures can be induced, including bicontinuous and lamellar arrangements.