Next Article in Journal
CRISPR/Cas9-Mediated Multiplex Genome Editing of the BnWRKY11 and BnWRKY70 Genes in Brassica napus L.
Previous Article in Journal
Mixed Peptide-Conjugated Chitosan Matrices as Multi-Receptor Targeted Cell-Adhesive Scaffolds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(9), 2715; https://doi.org/10.3390/ijms19092715
Submission received: 17 August 2018 / Revised: 5 September 2018 / Accepted: 6 September 2018 / Published: 11 September 2018
(This article belongs to the Section Materials Science)

Abstract

:
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.

Graphical Abstract

1. Introduction

The different types of thermotropic liquid-crystals can be characterized by varying degrees of molecular orientational and positional order. One of the most common is the nematic liquid-crystal phase. In the nematic phase, there is only orientational order but no positional order between molecules. The molecules of the smectic phase, on the other hand, are positioned such that they form a layered structure, in addition to their orientational order. There are many kinds of smectic phases, which differ in the particular arrangement of molecules within the layers.
Chirality is an important molecular characteristic that, when present, often results in formation of additional classes of liquid-crystal phases with interesting and useful optical properties. Under certain conditions, rod-like chiral molecules or molecular systems with a chiral dopant may arrange themselves into twisted super-structures. An example of this is the chiral-nematic phase, which has a continuously twisting helical structure. This phase is also referred to as the cholesteric phase, since it was first observed for cholesterol derivatives. Analogous to the relationship between nematic and cholesteric phases, the smectic twisted grain boundary (TGB) phase retains strong smectic ordering on a local scale, while exhibiting a larger-scale chiral rotation of the ordering director. In this case, the rotation occurs between adjacent slabs of smectic ordering, which are separated by grain boundaries. These boundaries are regularly spaced and this permits rotation of the director between slabs. Representative diagrams of these phases are shown in Figure 1. Further to these phases, the chiral line liquid phase has been observed from high-resolution calorimetry and X-ray experiments [1,2,3,4]. The chiral line liquid phase exists between the TGB and cholesteric phases. Previous works have reported that it has a short-range smectic layer formation, like those found in the TGB phases [2,4]. Following these previous works, several studies have been performed to investigate the phase behavior of the chiral line liquid, but these are few in number in comparison with other liquid-crystal phases [5,6,7,8,9,10,11]. The recent development of photonic technologies has enabled the creation tunable nanostructures of these phases. A variety of bulk phase properties are thus easily controllable. Examples include optically-tunable helical twisting power [12,13] and reflections [14], and light-driven handedness of helical superstructures inversion [15,16] as well as other-types of structural changes [13,17,18]. These developments will serve to promote the design and application of intelligent advanced functional materials [19].
The observed macroscopic structures of liquid crystals are strongly influenced by the collective behavior of their constituent molecules. Therefore, the study of these fluids at a molecular scale is important to further understanding of the origin of their bulk-scale physical properties. To study physics at a molecular scale, it is useful to carry out Monte Carlo or Molecular Dynamics simulations. These methods allow a better understanding of complex fluids, such as liquid crystals, at a microscopic scale, and reveal how molecular interactions give rise to the bulk scale phenomena commonly observed in experiments. For example, a recent simulation study on triangular prisms demonstrated how a chiral phase may form from molecules with purely steric interactions [20]. This finding partially supports experimental observations of helical filaments [21] and gives weight to the notion of a purely entropic mechanism behind the formation of some chiral nematic phases. Of course, the Gay–Berne model [22] and the hybrid Gay–Berne model [23,24,25,26,27] are useful to study the general behavior of liquid crystals, and are still commonly used. In part, this is because of their computational efficiency, but it is also due to the ease with which they may be tuned to allow focused examination of the effects of varying different molecular characteristics. There are notable examples of similar works for chiral liquid-crystals in the literature: Memmer et al. developed a simple coarse-grained model for chiral liquid-crystals, and demonstrated their phase behavior [28,29,30,31,32,33,34]; Allen et al. investigated the structure of the TGB phase by introducing modified boundary conditions [35]; Berardi et al. studied chiral interactions of a nematic phase with a chiral surface [36]; thermomechanical coupling and heat conduction in a cholesteric phase were studied by Sarman et al. [37]; Melle and Schlotthauer et al. studied nano-confined systems of chiral liquid-crystals [38,39]; and chiral superstructures were investigated using a linear rigid coarse-grained model by Yan et al. [40]. The studies mentioned above employed simple models, but give a general understanding of the characteristics of chiral liquid-crystals. As such, they are useful to predict physical properties of fluids of real molecules which share similar characteristics.
Dipole interactions are also an important characteristic to consider for liquid crystals, since they promote quick response times to the imposition of electric fields, while dipolar liquid-crystal materials are often used for industrial products. For achiral liquid crystals with dipoles, several studies have been reported. For example, we know from a number of Monte Carlo simulations of dipolar hard-spherocylinders that the addition of central dipoles, be they longitudinal or transverse, stabilizes the smectic-A phase at the expense of the nematic phase [41,42,43]. Conversely, terminal longitudinal dipoles increase the stability of the nematic phase while reducing that of the smectic-A phase [42,43,44,45]. An earlier study by Gwóźdź et al. [46] reported that the Gay–Berne system with central transverse dipoles does not yield an appreciable difference in the phase stability of the smectic phase, relative to the non-polar system, but does enhance the ordering within it.
As noted above, many studies have been conducted to promote a better understanding of chiral liquid crystals at the molecular scale. There are, to date, no theoretical or simulation studies that combine both chiral and dipolar electrostatics. In this study, therefore, we carefully investigate the effect of central longitudinal dipole interactions on the structure of chiral liquid crystals using Monte Carlo simulations. The current paper will serve to provide the first such analysis of the interplay between dipolar and chiral electrostatics and their effect on the phase behavior of liquid crystal fluids. The subsequent sections of this paper are arranged as follows: in Section 2, the results of these simulations are presented; the conclusions drawn from these findings are made in Section 3, and the methods used in our simulations are detailed in Section 4.

2. Results and Discussion

2.1. The Effect of Chirality on Non-Polar Liquid Crystals

Temperature-dependent phase transitions for non-polar, chiral, liquid crystals were extensively studied by Memmer et al.; however, their simulations were performed under full 3D-periodic boundary conditions that may have induced serious artifacts in the chiral phase due to unrealistic deformation of the equilibrium pitch. Therefore, we performed MC simulations for non-polar liquid crystals at similar conditions with planar walls to investigate the phase behavior carefully. All averaged physical properties for confined systems are calculated excluding molecules near the substrates because those molecules are strongly affected by the repulsive interactions and anchoring effects of the surface, as shown in Figure A1 and Figure A2. The simulation domain is divided into slices of equal volume, with each stacked on top of a neighbouring slice, in the direction of the z-axis. All physical properties are calculated separately for each slice, and averaged over all slices from z / L z = 0.2 to 0.8.
To examine the effect of chirality on phase behavior and system density, the temperature dependence of the averaged local density, as shown in Figure 2, is examined. First order phase transitions can be characterized by a sudden and sharp change in density relative to the change in temperature. The density decreases with temperature and two changes in gradient were observed, which implies the existence of three distinct phases for both achiral and chiral systems ( c = 0.0 and 0.5 ). Examination of the snapshots (Figure 3) and director profiles (Figure 4) reveals their identity. For achiral systems ( c = 0.0 ), smectic phases occur at low temperature, nematic phases at the intermediate temperature range, and isotropic phases at high temperature. The smectic-nematic phase transition occurs at 1.4 T 1.45 , and the nematic-isotropic phase transitions occur at 1.7 T 1.75 . For chiral systems ( c = 0.5 ), isotropic phases are present at high temperatures, while, in the intermediate temperature range, we see cholesteric phases. At the lowest temperatures, we also observe smectic phases. In fact, it is difficult to distinguish TGB phases from smectic phases in our studies due to limitation of systems size and wetting effects by solid substrates. The director correlation profile for z-axis s ( z ) decays as shown in Figure 4 and twisted structures may appear when the system is enough large; however, we here denote these phases as smectic tentatively. The smectic-cholesteric transition occurs at temperatures of 1.35 T 1.4 , and the cholesteric-isotropic phase transition occurs at 1.7 T 1.75 . Interestingly, the density of chiral systems are generally higher than those of the achiral, non-smectic, phases at the same temperatures. This increased density may be caused by the additive chiral potential ( 1 / r ) 7 , which may act to reduce molecule-to-molecule separation distances and potential well depth (see Appendix B). One should note this phenomenon is not found only for chiral Gay–Berne fluids, but also for other additive-type pair potentials, which are often used to model chiral interactions.
The structural properties of our simulated systems were examined by comparison of their orientational ordering. In Figure 5, the temperature dependences of the orientational order parameter for each of the simulated systems are plotted. For the achiral systems ( c = 0.0 ), slight decreases in P 2 were observed at 1.4 T 1.45 , which corresponds to the smectic-nematic phase transition. Significant decreases in P 2 were observed at 1.65 T 1.75 which corresponded to nematic-isotropic phase transitions. For chiral system ( c = 0.5 ), similarly, slight decreases in P 2 were observed at 1.35 T 1.4 corresponding to the smectic-cholesteric phase transition. Significant decreases in P 2 were observed at 1.65 T 1.75 which corresponded to cholesteric-isotropic phase transitions. Figure 5 clearly indicates that local orientational order in the chiral systems is almost equivalent to those found in the bulk of the achiral systems.
In the second stage of our structural analysis, the translational order parameter of the equilibrium systems was calculated using Equation (9). The temperature dependence of translational order parameter τ are plotted in Figure 6. For achiral system ( c = 0.0 ), significant decreases in τ at 1.4 T 1.45 corresponded to smectic-nematic phase transitions. For chiral system ( c = 0.5 ), significant decreases in τ at 1.35 T 1.4 corresponded to smectic-cholesteric phase transitions.
In the third stage of our structural analysis, the bond order parameter of the equilibrium systems was calculated using Equation (10). In Figure 7, the temperature dependences of the six-fold bond order parameter B 6 for each system are plotted. For achiral system ( c = 0.0 ), significant decreases in B 6 at 1.4 T 1.45 corresponded to smectic-nematic phase transitions. For chiral system ( c = 0.5 ), significant decreases in B 6 at 1.35 T 1.4 corresponded to smectic-cholesteric phase transitions. This indicates that the smectic phases obtained from our simulations have a local hexagonal packing structure, which has similarities to that of the smectic-B phase.

2.2. Effect of Dipole Interactions

In this section, the effect of central longitudinal dipole interactions on the structural properties of chiral liquid crystals was studied. At low temperatures, we observed the occurrence of metastable structures. For this reason, our focus here is restricted to a slightly higher temperature range: 1.5 T 2.0 .
We now examine how the phase behavior and system density are effected by the addition of dipoles. As with the non-polar systems above, this is done by inspection of the temperature dependence of the averaged local density, which is displayed in Figure 8. For achiral systems ( c = 0.0 ), two phases were observed in the temperature range; nematic phases occur at 1.5 T 1.65 , and isotropic phases at 1.75 T 2.0 for weak dipoles ( μ = 0.5 and 1.0 ), similar to our results for non-polar systems. For strong dipoles ( μ = 1.5 and 2.0 ), three phases were observed; smectic phases occur at low temperature, nematic phases at the intermediate temperature range, and isotropic phases at high temperature. The results clearly indicate that the smectic-nematic and nematic-isotropic phase transition points are shifted to higher temperatures with increased dipole strength. This tendency is consistent with previous work by Satoh et al. [47] and by Houssa et al. [48]. For chiral systems ( c = 0.5 ), on the other hand, only two phases were observed for the four dipole strengths investigated ( μ = 0.5 , 1.0 , 1.5 and 2.0 ); cholesteric phases occur at low temperature, and isotropic phases at high temperature. This highlights a difference in how the addition of dipole interactions affects achiral and chiral systems. Normally, the central longitudinal dipole interactions enhance the formation of layered structures, as demonstrated by the enhanced stability of the smectic phase with increasing dipole strength in the achiral systems. In the case of the chiral systems, however, the chiral interactions suppress the formation of such layered structures. In Figure 8, we see clearly that the system density increases with dipole strength; however, the effect is not significant for isotropic phase, relative to the ordered phases. In the isotropic phase, molecular behavior is almost entirely random and thus we do not see the creation of that stable dipolar domains which appear in ordered phases. It is for this reason that the density of the isotropic phase is less sensitive to the presence or absence of dipoles.
The structural properties of our simulated systems were examined by comparison of their orientational, translational and bond ordering. In Figure 9, the temperature dependences of the orientational order parameter for dipolar liquid crystals are plotted. For the achiral system ( c = 0.0 ), significant decreases in P 2 were observed, which corresponded to nematic-isotropic phase transitions for all dipole strengths ( μ = 0.5 , 1.0 , 1.5 and 2.0 ). For strong dipole interactions ( μ = 1.5 and 2.0 ), slight decreases in P 2 were also observed, which correspond to smectic-nematic phase transitions. For the chiral system ( c = 0.5 ), significant decreases in P 2 were observed which corresponded to cholesteric-isotropic phase transitions, for all dipole strengths ( μ = 0.5 , 1.0 , 1.5 and 2.0 ).
The temperature dependences of the translational order parameter τ are plotted in Figure 10. For the achiral system ( c = 0.0 ) with strong dipole interactions ( μ = 1.5 and 2.0 ), sharp decreases in τ are observed at low temperature. This corresponds to smectic-nematic phase transitions. On the other hand, the values of τ for chiral systems at low temperature are much smaller than those of achiral systems. Interestingly, τ increases with dipole strength both for achiral and chiral systems. This clearly indicates that dipole interactions promote layered structures; however, the transition to a full smectic phase is inhibited when chiral interactions are present. Moreover, τ for cholesteric phases with μ = 2.0 at 1.5 T 1.6 is appreciably larger than the other cases, a clear sign of possible presmectic ordering. Presmectic ordering [49] in liquid crystals, which may occur at high fluid densities, just below the nematic to smectic phase transition point, have been observed in several experiments [50,51]. The existence of short-range presmectic peaks in structure factor profiles was first predicted theoretically [52], and then subsequently observed in simulation [53]. In Figure 11, the temperature dependences of the six-fold bond order parameter B 6 for dipolar liquid crystals are plotted. In agreement with the study by Houssa et al. [54], we note that smectic-A phases were not obtained in our studies, either with or without a central dipole. For the achiral system ( c = 0.0 ), there are significant decreases in B 6 at low temperature for strong dipole conditions ( μ = 1.5 and 2.0 ), which correspond to smectic-nematic phase transitions. On the other hand, there is no such change for the chiral system ( c = 0.5 ). This also indicates that there is a difference in the structures obtained at low temperature between the achiral and chiral systems. Examination of the snapshots (Figure 12) reveals differences in τ of cholesteric phases. The cholesteric phases for μ = 2.0 at 1.5 T 1.6 appears to have short-range smectic ordering as shown in Figure 12a. These phases appear to exhibit regions of presmetic ordering, and/or the characteristics of the chiral line liquid phase, but further investigation will be required to identify this phase precisely. One should also note that there may be hysteresis for the smectic-nematic and nematic-isotropic phase transitions. According to previous studies [55,56], this hysteresis will be negligibly small, at least for non-polar Gay–Berne systems (see Figures 4 and 5 of [55] and Figure 2 [56]). Even with dipolar interactions, hysteresis at the nematic-isotropic phase transition has not been observed [48]. However, some hysteresis may occur at the smectic-nematic phase transition. A significant hysteresis has been observed for a mixture of Lennard-Jones and Gay–Berne models at strong charges or high density conditions (see Figures 2 and 7 of [57]), thus hysteresis may be increased with these cases. Further investigation will be required to determine the precise locations of the phase transition points between smectic-nematic or smectic-cholesteric phases.

3. Conclusions

To design highly functional liquid-crystal devices, a clear and systematic understanding of how molecular characteristics influence the structure and phase behavior is necessary. In this study, molecular simulations for model liquid-crystals were performed to investigate the effect of chirality and dipole strength upon the phase behavior. Using a series of analysis methods, accurate phase identities and phase transitions were studied in a systematic manner. Interestingly, the effect of introducing dipole interactions on structure and phase behavior is different for achiral and chiral systems. For achiral systems, dipole interactions promote liquid-crystal ordering, and smectic-nematic and nematic-isotropic phase transition points are shifted to higher temperatures with increased dipole strength. For chiral systems, on the other hand, the transition to a smectic phase does not occur in the temperature range investigated, although the cholesteric-isotropic phase transition is shifted to higher temperatures, as with an achiral system. The chiral interactions inhibit the formation of smectic phases. Moreover, strong dipoles may induce a short-range presmectic ordering in cholesteric phases at low temperature, and may result in a phase that is structurally similar to the chiral line liquid phase. Future work will of course be required to confirm this hypothesis.
For the chiral systems, the focus of our study was restricted to the cholesteric and isotropic phases since many metastable structures were observed at lower temperatures. The sampling of TGB-like phases by standard MC simulation is challenging; however, an extended ensemble approach such as multi-canonical simulation [58] or the generalized replica exchange method [59] may overcome this problem. Furthermore, these methods will also address the issue of hysteresis between smectic and nematic phases. System size is also a factor that limits the scope of our investigations. If found, the TGB phase will have a local smectic structure with a very long pitch distance. Large-scale simulation of such systems may be required to obtain reliable equilibrium data for the TGB phase, which presents a considerable computational burden.

4. Methods

4.1. Fluid–Fluid Interactions

In this work, we make use of a modified version of the Gay–Berne (GB) potential, developed by Memmer et al. [32]. This potential allows for the simulation of chiral liquid-crystal fluids and includes additive chiral and dipole interaction terms. The full potential is
U f f ( u ^ i , u ^ j , r i j ) = U GB ( u ^ i , u ^ j , r i j ) + c U c ( u ^ i , u ^ j , r i j ) + U d ( u ^ i , u ^ j , r i j ) ,
where u ^ i and u ^ j are the orientation vectors of molecules i and j, respectively, r i j is the unit vector of the inter-molecular center-to-center vector r i j = r j r i , and c is the chiral strength parameter. U GB ( u ^ i , u ^ j , r i j ) is the standard GB potential:
U GB ( u ^ i , u ^ j , r i j ) = 4 ϵ ( u ^ i , u ^ j , r ^ i j ) σ 0 r i j σ ( u ^ i , u ^ j , r ^ i j ) + σ 0 12 σ 0 r i j σ ( u ^ i , u ^ j , r ^ i j ) + σ 0 6 .
The chiral term, U c ( u ^ i , u ^ j , r i j ) , is given by
U c ( u ^ i , u ^ j , r i j ) = 4 ϵ ( u ^ i , u ^ j , r ^ i j ) σ 0 r i j σ ( u ^ i , u ^ j , r ^ i j ) + σ 0 7 ( u ^ i × u ^ j ) · r ^ i j ( u ^ i · u ^ j ) .
It should be noted that the parameters σ and ϵ in Equation (3) are identical to those of the standard GB potential.
The dipole term, U d ( u ^ i , u ^ j , r i j ) is given by
U d ( u ^ i , u ^ j , r i j ) = μ i μ j r i j 3 u ^ i · u ^ j 3 ( u ^ i · r ^ i j ) ( u ^ j · r ^ i j ) ,
where μ ^ i and μ ^ j are dipole moments of molecules i and j, namely, μ ^ i = μ i u ^ i and μ ^ j = μ j u ^ j where μ i and μ j are the dipole strengths.

4.2. Fluid–Substrate Interactions

Three-dimensional periodic boundary conditions are often used in such simulations to represent a bulk state, but they sometimes induce serious artifacts in chiral phases due to inappropriate periodic images, unless the box length along the twisted axis is an integer of half the pitch distance. In this work, therefore, we used solid substrates with planar anchoring [38,60,61] along the axis of helical rotation in order to break symmetry and eliminate boundary artifacts in that axis. The fluid-substrate potential energy is expressed as
U f s ( u ^ i , z i ) = U f s ( 1 ) ( u ^ i , z i ) + U f s ( 2 ) ( u ^ i , z i ) = ϵ f s { a 1 σ 0 z i 10 a 2 exp ( η z i ) z i g ( u ^ i ) + a 1 σ 0 L z z i 10 a 2 exp ( η ( L z z i ) ) L z z i g ( u ^ i ) } ,
where ϵ f s is the depth of the attractive well and z i is the perpendicular distance between molecule i and the lower substrate. L z is the cell length of the z-axis, thus L z z i is the perpendicular distance from the upper substrate to the center of molecule i. g ( u ^ i ) is the anchoring function, which depends on molecular orientation. In this study, we employed the following function:
g ( u ^ i ) = ( u ^ i · e ^ x ) 2 + ( u ^ i · e ^ y ) 2 .
Here, e ^ x and e ^ y are the unit vectors along the x-axis and the y-axis of the Cartesian coordinate system, respectively. The parameters a 1 and a 2 determine the location of the minimum of the fluid-substrate potential. They are unique functions of the screening length η 1 and are given by
a 1 = 1 + η σ 9 η σ , a 2 = 10 exp ( η σ ) 9 η σ .

4.3. Orientational Order Parameter

In order to perform orientational analysis for liquid-crystal phases, the second rank orientational order parameter, P 2 , is often used. The order parameter P 2 is the largest eigenvalue of the ordering matrix Q , which takes the form of
Q = 1 2 N i = 1 N 3 u ^ i u ^ i I ,
where N is the number of molecules and I is the identity matrix. If all molecules in the system are oriented along the same axis, P 2 is equal to 1, while P 2 is equal to 0 if molecule oritentations are distributed randomly.

4.4. Translational Order Parameter

To identify the layered structures of smectic phases, the translational order parameter τ is commonly used. In this study, the following equation is used to compute τ :
τ ( d ) = j = 1 N exp ( 2 π i r j , | | / d ) ,
where r j , | | = r j · n ^ is the coordinate parallel to the director n ^ and d is layer spacing. The translational order parameter τ and the layer spacing are given by the maximum value of τ ( d ) .

4.5. Bond Order Parameter

Some liquid-crystal phases have positional order within layered structures, in addition to orientational order. To characterize such positional ordering, it is useful to consider the bond order parameter. This function is defined as
B ν = 1 ν N j = 1 N k = 1 ν exp ( i ν ϕ j k ) ,
where ϕ j k is the angle between bond-linked particles j and k, and a fixed reference frame. ν is the number of nearest-neighbor bonds. In this study, ν = 6 . This value is chosen in order to identify local hexagonal packing, of the type typically found in smectic-B-like phases.

4.6. Simulation Conditions

Monte Carlo simulations were performed for both achiral and chiral systems with 3456 liquid-crystal molecules. To represent these liquid-crystal molecules, we use the chiral Gay–Berne potential, with a chirality parameter of c = 0.0 for achiral molecules, and c = 0.5 for chiral molecules. The GB potential usually contains four parameters ( κ , κ , μ , ν ) which determine the anisotropy in the molecular interactions. Molecular shape anisotropy is quantified by κ σ e / σ s , where σ e and σ s are the molecule size parameters for end-to-end and side-by-side configurations, respectively. The amount of anisotropy is determined by μ , ν , and κ ϵ e / ϵ s . Here, ϵ e and ϵ s are the well depths for the end-to-end and side-to-side configurations, respectively. For our study on the effects of interaction-strength anisotropy, the potential parameters ( κ , κ , μ , ν ) = (3.0, 0.2, 1.0, 2.0) were used. These potentials are truncated using the spherical cutoff method, and potentials are shifted to reduce cutoff artifacts. In addition, solid substrates with planar anchoring are introduced in the z-axis, and 2D-periodic boundary conditions are used in the x- and y-axes to mimic the effect of large bulk systems while eliminating periodic boundary artifacts along the axis of helical rotation. In this work, all quantities are described using conventional reduced properties: reduced temperature T = k B T / ϵ 0 ; reduced pressure P = P σ 0 3 / ϵ 0 ; reduced potential energy per particle U = U / ϵ 0 ; reduced distance r = r / σ 0 ; reduced density ρ = N / V σ 0 3 , and reduced dipole strength μ = μ / σ 0 3 ϵ 0 . The cutoff radius for the potential was set at a reduced length of r c = κ + 2.0 = 5.0 . All simulations were performed at a constant reduced pressure tangential to the substrates of P x x = P y y = 3.5 . Several isotherms in the range T = 1.2 2.5 were investigated. In order to prevent unrealistic cell deformation, the cell length in the z direction is fixed, and anisotropic volume changes are attempted by rescaling the cell lengths in the other two cartesian axes. For the fluid–substrate interactions, η 1 = σ 0 and ϵ f s = 10 ϵ 0 were chosen. The dipole strength μ was set to the following five values: 0.0 , 0.5 , 1.0 , 1.5 and 2.0 , in order to examine the relationship between the system’s phase behavior and the dipole moment. For the initial state of the dipolar systems, an isotropic configuration was used in order to reduce the chances of the systems becoming trapped at meta-stable local minima during the equilibration stage. If meta-stable states do occur during this process, those simulations are restarted from the most recent configurations of a neighboring temperature. For each state point, a long MC simulation was run for 4 × 10 6 MC cycles to obtain equilibrium conditions. A further 1 × 10 6 MC cycles were then performed to calculate equilibrium averages.

Author Contributions

T.N. and K.Y. designed the research; T.N. performed simulations and data analysis; T.N. and P.E.B. wrote the paper, and all authors contributed to the final version.

Funding

This research was funded by Grant-in-Aid for the Japan Society for the Promotion of Science (JSPS) Fellows, and a MEXT Grant-in-Aid for the Program for Leading Graduate Schools.

Acknowledgments

We acknowledge the financial support by Grant-in-Aid for the Japan Society for the Promotion of Science (JSPS) Fellows, and a MEXT Grant-in-Aid for the Program for Leading Graduate Schools.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Effect of Confinement by Solid Substrates

Confinement by solid substrates affects local structures, and may change the phase behavior of the system as a whole. In this section, therefore, the simulation results for bulk and confined systems were compared in order to investigate the effect of the substrate. Figure A1 and Figure A2 show local density profiles and local profiles of the orientational order parameter P 2 in the direction of the z-axis, respectively. Figure A1 shows that the density decreases near the substrates, regardless of the phase type, because of the repulsive interaction between fluid and substrate. For the smectic phase, the density in the central portion of the confined system is larger than that of the bulk system, but the density of the central region for nematic and isotropic phases under confinement is almost the same as their corresponding bulk phase densities. Figure A2 shows that orientational order increases in proximity to the substrates for nematic and isotropic phases, but this effect is not observed for the smectic phase. The lower substrate induces slightly more surface ordering than the upper substrate. This is due its directional anchoring effect.
Figure A1. Local density profiles for achiral liquid-crystal systems ( c = 0.0 ) in different phases: smectic ( T = 1.2 ), nematic ( T = 1.5 ), and isotropic ( T = 2.0 ). The continuous lines correspond to bulk systems, while the dashed lines represent confined systems in the presence of solid substrates.
Figure A1. Local density profiles for achiral liquid-crystal systems ( c = 0.0 ) in different phases: smectic ( T = 1.2 ), nematic ( T = 1.5 ), and isotropic ( T = 2.0 ). The continuous lines correspond to bulk systems, while the dashed lines represent confined systems in the presence of solid substrates.
Ijms 19 02715 g0a1
For the next step, the temperature dependence of several properties were compared in order to investigate the effect of substrates on phase behavior. All averaged physical properties for confined systems are calculated excluding molecules near the substrates because they are strongly affected by the repulsive interactions and surface anchoring effects, as shown in Figure A1 and Figure A2. The simulation domain is divided into slices of equal volume, each stacked on top of a neighbouring slice, in the direction of the z-axis. All physical properties are calculated for each slice separately, and then averaged over slices from z / L z = 0.2 to 0.8, which corresponds to the bulk region. Figure A3, Figure A4 and Figure A5 show the temperature versus the averaged local density, the averaged local orietational order parameter P 2 , and the averaged local six-fold bond order parameter B 6 , respectively. In Figure A3, there are two significant decreases in the density, one at 1.4 T 1.45 and the other at 1.7 T 1.75 for both bulk and confined systems, which implies the existence of three distinct phases. For both systems, smectic phases occur at low temperature, nematic phases at the intermediate temperature range, and isotropic phases at high temperature. The density of smectic phases for confined systems are higher than those of the bulk system. On the other hand, the density for nematic and isotropic phases in the central region of the confined systems is almost the same as their corresponding bulk system densities. Figure A4 shows that there are significant decreases in P 2 at 1.4 T 1.45 and 1.7 T 1.75 , which corresponded to smectic-nematic and nematic-isotropic phase transitions, respectively. Figure A5 shows that there is a significant decrease in B 6 at 1.4 T 1.45 , which corresponds to the smectic-nematic phase transition. This clearly indicates that, for both systems, the smectic phase has a local hexagonal packing structure like that found in the smectic-B phase, but it is worth noting that the B 6 values are higher for the confined systems.
Figure A2. Local profiles of orientational order parameter P 2 for achiral liquid-crystal systems ( c = 0.0 ) in different phases: smectic ( T = 1.2 ), nematic ( T = 1.5 ), and isotropic ( T = 2.0 ). The continuous lines correspond to bulk systems, and the dashed lines represent confined systems in the presence of solid substrates.
Figure A2. Local profiles of orientational order parameter P 2 for achiral liquid-crystal systems ( c = 0.0 ) in different phases: smectic ( T = 1.2 ), nematic ( T = 1.5 ), and isotropic ( T = 2.0 ). The continuous lines correspond to bulk systems, and the dashed lines represent confined systems in the presence of solid substrates.
Ijms 19 02715 g0a2
These results clearly indicate that liquid-crystal molecules near substrates are strongly affected by the repulsive interactions and surface anchoring effects. This may lead to the formation of locally ordered structures that are different from those found in the bulk phase. This is especially true for the smectic phase; however, the bulk phase behavior of these systems was not found to change with confinement. It is therefore possible to study the phase behavior of liquid-crystal systems, in the presence of confinement, which can be used to eliminate inappropriate periodic images which would otherwise be problematic in chiral systems.
Figure A3. Temperature dependence of the averaged local density for bulk and confined liquid-crystal systems ( c = 0.0 ). The open squares correspond to bulk systems, and the filled triangles represent confined systems with the solid substrates. For both systems, smectic phases (Sm) occur at low temperature, nematic phases (N) at the intermediate temperature range, and isotropic phases (I) at high temperature.
Figure A3. Temperature dependence of the averaged local density for bulk and confined liquid-crystal systems ( c = 0.0 ). The open squares correspond to bulk systems, and the filled triangles represent confined systems with the solid substrates. For both systems, smectic phases (Sm) occur at low temperature, nematic phases (N) at the intermediate temperature range, and isotropic phases (I) at high temperature.
Ijms 19 02715 g0a3
Figure A4. Temperature dependence of the averaged local orientational order parameter P 2 for bulk and confined liquid-crystal systems ( c = 0.0 ). The open squares correspond to bulk systems, and the filled triangles represent confined systems with the solid substrates. For both systems, smectic phases (Sm) occur at low temperature, nematic phases (N) at the intermediate temperature range, and isotropic phases (I) at high temperature.
Figure A4. Temperature dependence of the averaged local orientational order parameter P 2 for bulk and confined liquid-crystal systems ( c = 0.0 ). The open squares correspond to bulk systems, and the filled triangles represent confined systems with the solid substrates. For both systems, smectic phases (Sm) occur at low temperature, nematic phases (N) at the intermediate temperature range, and isotropic phases (I) at high temperature.
Ijms 19 02715 g0a4
Figure A5. Temperature dependence of the averaged local six-fold bond order parameter B 6 for bulk and confined liquid-crystal systems ( c = 0.0 ). The open squares correspond to bulk systems, and the filled triangles represent confined systems with the solid substrates. For both systems, smectic phases (Sm) occur at low temperature, nematic phases (N) at the intermediate temperature range, and isotropic phases (I) at high temperature.
Figure A5. Temperature dependence of the averaged local six-fold bond order parameter B 6 for bulk and confined liquid-crystal systems ( c = 0.0 ). The open squares correspond to bulk systems, and the filled triangles represent confined systems with the solid substrates. For both systems, smectic phases (Sm) occur at low temperature, nematic phases (N) at the intermediate temperature range, and isotropic phases (I) at high temperature.
Ijms 19 02715 g0a5

Appendix B. Chiral Gay–Berne potential

Energy contour maps are plotted in Figure A6 to help explain the reasons behind the density increases for chiral systems. The additive-type chiral potential has potential minima at a molecule-to-molecule separation distance of r = 0 , thus the separation distance and potential well depth of chiral Gay–Berne model decreases, relative to the achiral model. For c = 0.5 , the separation distance is reduced from r = 1.12 to r = 1.1 , and the well depth from 2.777289 to 3.080405 .
Figure A6. Energy contour maps: (upper) standard Gay–Berne potential with potential parameters ( κ , κ , μ , ν ) = (3.0, 0.2, 1.0, 2.0); (below) chiral Gay–Berne potential with potential parameters ( κ , κ , μ , ν ) = (3.0, 0.2, 1.0, 2.0) and c = 0.5 .
Figure A6. Energy contour maps: (upper) standard Gay–Berne potential with potential parameters ( κ , κ , μ , ν ) = (3.0, 0.2, 1.0, 2.0); (below) chiral Gay–Berne potential with potential parameters ( κ , κ , μ , ν ) = (3.0, 0.2, 1.0, 2.0) and c = 0.5 .
Ijms 19 02715 g0a6

References

  1. Kamien, R.D.; Lubensky, T. Twisted line liquids. J. Phys. I 1993, 3, 2131–2138. [Google Scholar] [CrossRef] [Green Version]
  2. Navailles, L.; Pansu, B.; Gorre-Talini, L.; Nguyen, H. Structural study of a commensurate TGB A phase and of a presumed chiral line liquid phase. Phys. Rev. Lett. 1998, 81, 4168–4171. [Google Scholar] [CrossRef]
  3. Chan, T.; Garland, C.W.; Nguyen, H. Calorimetric study of chiral liquid crystals with a twist-grain-boundary phase. Phys. Rev. E 1995, 52, 5000–5003. [Google Scholar] [CrossRef]
  4. Ybert, C.; Navailles, L.; Pansu, B.; Rieutord, F.; Nguyen, H.; Barois, P. Structural study of the liquid-crystal analog of the vortex liquid phase in type-II superconductors. Europhys. Lett. 2003, 63, 840–845. [Google Scholar] [CrossRef]
  5. Iannacchione, G.; Garland, C.; Mieczkowski, J.; Gorecka, E. Phase behavior in the reentrant-nematic region of chiral frustrated smectic liquid crystals. Phys. Rev. E 1998, 58, 595–601. [Google Scholar] [CrossRef]
  6. Anakkar, A.; Isaert, N.; Ismaili, M.; Buisine, J.M.; Nguyen, H. Experimental observation of multicritical behaviors in pressure-temperature diagrams of pure compounds with twist grain boundary phases. Phys. Rev. E 1999, 60, 620–627. [Google Scholar] [CrossRef]
  7. Ismaïli, M.; Anakkar, A.; Joly, G.; Isaert, N.; Nguyen, H. Multicritical behaviors and an induced twist grain boundary phase in a binary liquid crystalline mixture. Phys. Rev. E 2000, 61, 519–525. [Google Scholar] [CrossRef]
  8. Adorjan, A.; Stojadinovic, S.; Sukhomlinova, L.; Twieg, R.; Sprunt, S. Light scattering study of a twist grain boundary liquid crystal. Phys. Rev. Lett. 2003, 90, 035503. [Google Scholar] [CrossRef] [PubMed]
  9. Yamamoto, J.; Nishiyama, I.; Inoue, M.; Yokoyama, H. Optical isotropy and iridescence in a smectic ‘blue phase’. Nature 2005, 437, 525–528. [Google Scholar] [CrossRef] [PubMed]
  10. Cordoyiannis, G.; Jampani, V.S.R.; Kralj, S.; Dhara, S.; Tzitzios, V.; Basina, G.; Nounesis, G.; Kutnjak, Z.; Tripathi, C.S.P.; Losada-Pérez, P.; et al. Different modulated structures of topological defects stabilized by adaptive targeting nanoparticles. Soft Matter 2013, 9, 3956–3964. [Google Scholar] [CrossRef]
  11. Trček, M.; Cordoyiannis, G.; Tzitzios, V.; Kralj, S.; Nounesis, G.; Lelidis, I.; Kutnjak, Z. Nanoparticle-induced twist-grain boundary phase. Phys. Rev. E 2014, 90, 032501. [Google Scholar] [CrossRef] [PubMed]
  12. Sato, H.; Yamagishi, A. Application of Δ- and Λ-isomerism of Octahedral Metal Complexes for Inducing Chiral Nematic Phases. Int. J. Mol. Sci. 2009, 10, 4559–4574. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, L.; Wang, M.; Wang, L.; Yang, D.k.; Yu, H.; Yang, H. Polymeric infrared reflective thin films with ultra-broad bandwidth. Liq. Cryst. 2016, 43, 750–757. [Google Scholar] [CrossRef]
  14. Wang, L.; Dong, H.; Li, Y.; Xue, C.; Sun, L.D.; Yan, C.H.; Li, Q. Reversible near-infrared light directed reflection in a self-organized helical superstructure loaded with upconversion nanoparticles. J. Am. Chem. Soc. 2014, 136, 4480–4483. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, L.; Dong, H.; Li, Y.; Liu, R.; Wang, Y.F.; Bisoyi, H.K.; Sun, L.D.; Yan, C.H.; Li, Q. Luminescence-driven reversible handedness inversion of self-organized helical superstructures enabled by a novel near-infrared light nanotransducer. Adv. Mater. 2015, 27, 2065–2069. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, L.; Chen, D.; Gutierrez-Cuevas, K.G.; Bisoyi, H.K.; Fan, J.; Zola, R.S.; Li, G.; Urbas, A.M.; Bunning, T.J.; Weitz, D.A.; et al. Optically reconfigurable chiral microspheres of self-organized helical superstructures with handedness inversion. Mater. Horiz. 2017, 4, 1190–1195. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, L.; Bisoyi, H.K.; Zheng, Z.; Gutierrez-Cuevas, K.G.; Singh, G.; Kumar, S.; Bunning, T.J.; Li, Q. Stimuli-directed self-organized chiral superstructures for adaptive windows enabled by mesogen-functionalized graphene. Mater. Today 2017, 20, 230–237. [Google Scholar] [CrossRef]
  18. Wang, L.; Li, Q. Photochromism into nanosystems: Towards lighting up the future nanoworld. Chem. Soc. Rev. 2018, 47, 1044–1097. [Google Scholar] [CrossRef] [PubMed]
  19. Wang, L.; Li, Q. Stimuli-Directing Self-Organized 3D Liquid-Crystalline Nanostructures: From Materials Design to Photonic Applications. Adv. Func. Mater. 2016, 26, 10–28. [Google Scholar] [CrossRef]
  20. Dussi, S.; Dijkstra, M. Entropy-driven formation of chiral nematic phases by computer simulations. Nat. Commun. 2016, 7, 11175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Barry, E.; Hensel, Z.; Dogic, Z.; Shribak, M.; Oldenbourg, R. Entropy-Driven Formation of a Chiral Liquid- Crystalline Phase of Helical Filaments. Phys. Rev. Lett. 2006, 96, 018305. [Google Scholar] [CrossRef] [PubMed]
  22. Gay, J.G.; Berne, B.J. Modification of the overlap potential to mimic a linear site-site potential. J. Chem. Phys. 1981, 74, 3316–3319. [Google Scholar] [CrossRef]
  23. Wilson, M.R. Molecular dynamics simulations of flexible liquid crystal molecules using a Gay–Berne/Lennard-Jones model. J. Chem. Phys. 1997, 107, 8654–8663. [Google Scholar] [CrossRef]
  24. McBride, C.; Wilson, M.R. Molecular dynamics simulations of a flexible liquid crystal. Mol. Phys. 1999, 97, 511–522. [Google Scholar] [CrossRef]
  25. Fukunaga, H.; Takimoto, J.; Doi, M. Molecular dynamics simulation study on the phase behavior of the Gay–Berne model with a terminal dipole and a flexible tail. J. Chem. Phys. 2004, 120, 7792–7800. [Google Scholar] [CrossRef] [PubMed]
  26. Prampolini, G. Parametrization and validation of coarse grained force-fields derived from ab initio calculations. J. Chem. Theory Comput. 2006, 2, 556–567. [Google Scholar] [CrossRef] [PubMed]
  27. Tagashira, K.; Takahashi, K.Z.; Fukuda, J.I.; Aoyagi, T. Development of Coarse-Grained Liquid-Crystal Polymer Model with Efficient Electrostatic Interaction: Toward Molecular Dynamics Simulations of Electroactive Materials. Materials 2018, 11, 83. [Google Scholar] [CrossRef] [PubMed]
  28. Memmer, R.; Kuball, H.G.; Schönhofer, A. Computer simulation of chiral liquid crystal phases. I. The polymorphism of the chiral Gay–Berne fluid. Liq. Cryst. 1993, 15, 345–360. [Google Scholar] [CrossRef]
  29. Memmer, R.; Kuball, H.G.; Schönhofer, A. Computer simulation of chiral liquid crystal phases. IV. Intermolecular chirality transfer to rotamers in a cholesteric phase. Liq. Cryst. 1995, 19, 749–753. [Google Scholar] [CrossRef]
  30. Memmer, R. Computer simulation of chiral liquid crystal phases III. A cholesteric phase formed by chiral Gay–Berne atropisomers. Mol. Phys. 1996, 89, 1633–1649. [Google Scholar]
  31. Memmer, R. Computer simulation of chiral liquid crystal phases—VII. The Chiral Gay–Berne Discogen. Ber. Bunsen-Ges. Phys. Chem. 1998, 102, 1002–1010. [Google Scholar] [CrossRef]
  32. Memmer, R. Determination of equilibrium pitch of cholesteric phases by isobaric-isothermal Monte Carlo simulation. J. Chem. Phys. 2001, 114, 8210–8222. [Google Scholar] [CrossRef]
  33. Memmer, R. Liquid crystal phases of achiral banana-shaped molecules: A computer simulation study. Liq. Cryst. 2002, 29, 483–496. [Google Scholar] [CrossRef]
  34. Memmer, R.; Fliegans, O. Monte Carlo simulation of twisted nematic and supertwisted nematic liquid crystal cells. Phys. Chem. Chem. Phys. 2003, 5, 558–566. [Google Scholar] [CrossRef]
  35. Allen, M.P.; Warren, M.A.; Wilson, M.R. Molecular-dynamics simulation of the smectic-A* twist grain-boundary phase. Phys. Rev. E 1998, 57, 5585–5596. [Google Scholar] [CrossRef]
  36. Berardi, R.; Kuball, H.G.; Memmer, R.; Zannoni, C. Chiral induction in nematics a computer simulation study. J. Chem. Soc. Faraday Trans. 1998, 94, 1229–1234. [Google Scholar] [CrossRef]
  37. Sarman, S.; Laaksonen, A. Thermomechanical coupling, heat conduction and director rotation in cholesteric liquid crystals studied by molecular dynamics simulation. Phys. Chem. Chem. Phys. 2013, 15, 3442–3453. [Google Scholar] [CrossRef] [PubMed]
  38. Melle, M.; Theile, M.; Hall, C.K.; Schoen, M. Nanoconfinement-induced structures in chiral liquid crystals. Int. J. Mol. Sci. 2013, 14, 17584–17607. [Google Scholar] [CrossRef] [PubMed]
  39. Schlotthauer, S.; Skutnik, R.A.; Stieger, T.; Schoen, M. Defect topologies in chiral liquid crystals confined to mesoscopic channels. J. Chem. Phys. 2015, 142, 194704. [Google Scholar] [CrossRef] [PubMed]
  40. Yan, F.; Hixson, C.A.; Earl, D.J. Computer simulations of linear rigid particles that form chiral superstructures and tilted smectic phases. Soft Matter 2009, 5, 4477–4483. [Google Scholar] [CrossRef]
  41. Gil-Villegas, A.; McGrother, S.C.; Jackson, G. Chain and ring structures in smectic phases of molecules with transverse dipoles. Chem. Phys. Lett. 1997, 269, 441–447. [Google Scholar] [CrossRef]
  42. Gil-Villegas, A.; McGrother, S.C.; Jackson, G. Reaction-field and Ewald summation methods in Monte Carlo simulations of dipolar liquid crystals. Mol. Phys. 1997, 92, 723–734. [Google Scholar] [CrossRef]
  43. Gil-Villegas, A.; Jackson, G.; McGrother, S.C. Computer simulation of dipolar liquid crystals. J. Mol. Liq. 1998, 76, 171–181. [Google Scholar] [CrossRef]
  44. McGrother, S.C.; Williamson, D.C.; Jackson, G. A re-examination of the phase diagram of hard spherocylinders. J. Chem. Phys. 1996, 104, 6755–6771. [Google Scholar] [CrossRef]
  45. Satoh, K.; Mita, S.; Kondo, S. Monte Carlo simulations using the dipolar Gay–Berne model: Effect of terminal dipole moment on mesophase formation. Chem. Phys. Lett. 1996, 255, 99–104. [Google Scholar] [CrossRef]
  46. Gwóźdź, E.; Bródka, A.; Pasterny, K. Molecular dynamics simulation study of a liquid crystal model: Gay–Berne particles with transverse dipole moments. Chem. Phys. Lett. 1997, 267, 557–562. [Google Scholar] [CrossRef]
  47. Satoh, K.; Mita, S.; Kondo, S. Monte Carlo simulations on mesophase formation using dipolar Gay–Berne model. Liq. Cryst. 1996, 20, 757–763. [Google Scholar] [CrossRef]
  48. Houssa, M.; Rull, L.F.; McGrother, S.C. Effect of dipolar interactions on the phase behavior of the Gay–Berne liquid crystal model. J. Chem. Phys. 1998, 109, 9529–9542. [Google Scholar] [CrossRef]
  49. De Gennes, P. Interactions between solid surfaces in a presmectic fluid. Langmuir 1990, 6, 1448–1450. [Google Scholar] [CrossRef]
  50. Dogic, Z.; Fraden, S. Cholesteric Phase in Virus Suspensions. Langmuir 2000, 16, 7820–7824. [Google Scholar] [CrossRef] [Green Version]
  51. Kralj, S.; Cordoyiannis, G.; Zidanšek, A.; Lahajnar, G.; Amenitsch, H.; Žumer, S.; Kutnjak, Z. Presmectic wetting and supercritical-like phase behavior of octylcyanobiphenyl liquid crystal confined to controlled-pore glass matrices. J. Chem. Phys. 2007, 127, 154905. [Google Scholar] [CrossRef] [PubMed]
  52. van der Schoot, P. Self-assembly of globular particles in a nematic dispersion of colloidal rods. J. Chem. Phys. 2002, 117, 3537–3540. [Google Scholar] [CrossRef]
  53. Savenko, S.V.; Dijkstra, M. Accuracy of measuring the nematic order from intensity scatter: A simulation study. Phys. Rev. E 2004, 70, 011705. [Google Scholar] [CrossRef] [PubMed]
  54. Houssa, M.; Rull, L.F.; Romero-Enrique, J.M. Bilayered smectic phase polymorphism in the dipolar Gay–Berne liquid crystal model. J. Chem. Phys. 2009, 130, 154504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Bates, M.; Luckhurst, G. Computer simulation studies of anisotropic systems. XXX. The phase behavior and structure of a Gay–Berne mesogen. J. Chem. Phys. 1999, 110, 7087–7108. [Google Scholar] [CrossRef]
  56. de Miguel, E.; Martin del Rio, E.; Blas, F.J. Stability of smectic phases in the Gay–Berne model. J. Chem. Phys. 2004, 121, 11183–11194. [Google Scholar] [CrossRef] [PubMed]
  57. Saielli, G.; Margola, T.; Satoh, K. Tuning Coulombic interactions to stabilize nematic and smectic ionic liquid crystal phases in mixtures of charged soft ellipsoids and spheres. Soft Matter 2017, 13, 5204–5213. [Google Scholar] [CrossRef] [PubMed]
  58. Berg, B.A.; Neuhaus, T. Multicanonical ensemble: A new approach to simulate first-order phase transitions. Phys. Rev. Lett. 1992, 68, 9. [Google Scholar] [CrossRef] [PubMed]
  59. Kim, J.; Keyes, T.; Straub, J.E. Generalized replica exchange method. J. Chem. Phys. 2010, 132, 224107. [Google Scholar] [CrossRef] [PubMed]
  60. Mazza, M.G.; Schoen, M. Structure and dynamics of reentrant nematics: Any open questions after almost 40 years? Int. J. Mol. Sci. 2011, 12, 5352–5372. [Google Scholar] [CrossRef] [PubMed]
  61. Stieger, T.; Mazza, M.G.; Schoen, M. Diffusivity Maximum in a Reentrant Nematic Phase. Int. J. Mol. Sci. 2012, 13, 7854–7871. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic view of representative liquid-crystal phases. (a) smectic phase; (b) nematic phase; (c) cholesteric phase and (d) twisted grain boundary (TGB) phase.
Figure 1. Schematic view of representative liquid-crystal phases. (a) smectic phase; (b) nematic phase; (c) cholesteric phase and (d) twisted grain boundary (TGB) phase.
Ijms 19 02715 g001
Figure 2. Temperature dependence of the averaged local density for achiral and chiral liquid-crystal systems without dipoles. For achiral systems ( c = 0.0 ), smectic phases occur at lower temperature, nematic phases at the intermediate temperature range, and isotropic phases at higher temperature. For chiral systems ( c = 0.5 ), similarly, smectic phases occur at lower temperature, cholesteric phases at the intermediate temperature range, and isotropic phases at higher temperature. The system densities of chiral phases are generally higher than those of achiral phases.
Figure 2. Temperature dependence of the averaged local density for achiral and chiral liquid-crystal systems without dipoles. For achiral systems ( c = 0.0 ), smectic phases occur at lower temperature, nematic phases at the intermediate temperature range, and isotropic phases at higher temperature. For chiral systems ( c = 0.5 ), similarly, smectic phases occur at lower temperature, cholesteric phases at the intermediate temperature range, and isotropic phases at higher temperature. The system densities of chiral phases are generally higher than those of achiral phases.
Ijms 19 02715 g002
Figure 3. Representative snapshots of the non-polar achiral system ( c = 0.0 ): (a) smectic phase, at lower temperature ( T = 1.2 ); (b) nematic phase at the intermediate temperature ( T = 1.5 ); (c) isotropic phase at higher temperature ( T = 2.0 ), and the non-polar chiral system ( c = 0.5 ): (d) smectic phase, at lower temperature ( T = 1.2 ); (e) cholesteric phase at the intermediate temperature ( T = 1.5 ); and (f) isotropic phase at higher temperature ( T = 2.0 ).
Figure 3. Representative snapshots of the non-polar achiral system ( c = 0.0 ): (a) smectic phase, at lower temperature ( T = 1.2 ); (b) nematic phase at the intermediate temperature ( T = 1.5 ); (c) isotropic phase at higher temperature ( T = 2.0 ), and the non-polar chiral system ( c = 0.5 ): (d) smectic phase, at lower temperature ( T = 1.2 ); (e) cholesteric phase at the intermediate temperature ( T = 1.5 ); and (f) isotropic phase at higher temperature ( T = 2.0 ).
Ijms 19 02715 g003
Figure 4. Local profiles of the director correlation function in the z-axis s ( z ) = n 0 · n z for achiral and chiral liquid-crystal systems without dipoles, where n z is the local orientation director along the z-axis. For achiral systems ( c = 0.0 ), strong correlations of director are observed for smectic (Sm) and nematic (N) phases; they are ordered in one direction. For chiral systems ( c = 0.5 ), director correlation changes near the solid substrates for smectic (Sm) because of the wetting effect. The director correlation of cholesteric phase ( N ) decays because of the helical structure. For the isotropic phases (I) in the achiral and chiral systems, random director correlations are a characteristic of this phase.
Figure 4. Local profiles of the director correlation function in the z-axis s ( z ) = n 0 · n z for achiral and chiral liquid-crystal systems without dipoles, where n z is the local orientation director along the z-axis. For achiral systems ( c = 0.0 ), strong correlations of director are observed for smectic (Sm) and nematic (N) phases; they are ordered in one direction. For chiral systems ( c = 0.5 ), director correlation changes near the solid substrates for smectic (Sm) because of the wetting effect. The director correlation of cholesteric phase ( N ) decays because of the helical structure. For the isotropic phases (I) in the achiral and chiral systems, random director correlations are a characteristic of this phase.
Ijms 19 02715 g004
Figure 5. Temperature dependence of the averaged local orientational order parameter P 2 for achiral and chiral liquid-crystal systems without dipoles. For achiral systems ( c = 0.0 ), smectic phases occur at a lower temperature, nematic phases at the intermediate temperature range, and isotropic phases at higher temperature. For chiral systems ( c = 0.5 ), similarly, smectic phases occur at a lower temperature, cholesteric phases at the intermediate temperature range, and isotropic phases at higher temperature. The local orientational order of chiral systems is almost equivalent to those of the achiral systems.
Figure 5. Temperature dependence of the averaged local orientational order parameter P 2 for achiral and chiral liquid-crystal systems without dipoles. For achiral systems ( c = 0.0 ), smectic phases occur at a lower temperature, nematic phases at the intermediate temperature range, and isotropic phases at higher temperature. For chiral systems ( c = 0.5 ), similarly, smectic phases occur at a lower temperature, cholesteric phases at the intermediate temperature range, and isotropic phases at higher temperature. The local orientational order of chiral systems is almost equivalent to those of the achiral systems.
Ijms 19 02715 g005
Figure 6. Temperature dependence of the averaged local translational order parameter τ for achiral and chiral liquid-crystal systems without dipoles. For achiral systems ( c = 0.0 ), smectic phases occur at a lower temperature, nematic phases at the intermediate temperature range, and isotropic phases at a higher temperature. For chiral systems ( c = 0.5 ), similarly, smectic phases occur at lower temperature range, cholesteric phases at the intermediate temperature range, and isotropic phases at higher temperature. This figure reveals clear evidence of layered structures found in smectic phases.
Figure 6. Temperature dependence of the averaged local translational order parameter τ for achiral and chiral liquid-crystal systems without dipoles. For achiral systems ( c = 0.0 ), smectic phases occur at a lower temperature, nematic phases at the intermediate temperature range, and isotropic phases at a higher temperature. For chiral systems ( c = 0.5 ), similarly, smectic phases occur at lower temperature range, cholesteric phases at the intermediate temperature range, and isotropic phases at higher temperature. This figure reveals clear evidence of layered structures found in smectic phases.
Ijms 19 02715 g006
Figure 7. Temperature dependence of the averaged local six-fold bond order parameter B 6 for achiral and chiral liquid-crystal systems without dipoles. For achiral systems ( c = 0.0 ), smectic phases occur at lower temperature, nematic phases at the intermediate temperature range, and isotropic phases at higher temperature. For chiral systems ( c = 0.5 ), similarly, smectic phases occur at lower temperature, cholesteric phases at the intermediate temperature range, and isotropic phases at higher temperature. The figure clearly indicates that the smectic phases obtained from our simulations have a local hexagonal packing structure, which has similarities to that of the smectic-B phase.
Figure 7. Temperature dependence of the averaged local six-fold bond order parameter B 6 for achiral and chiral liquid-crystal systems without dipoles. For achiral systems ( c = 0.0 ), smectic phases occur at lower temperature, nematic phases at the intermediate temperature range, and isotropic phases at higher temperature. For chiral systems ( c = 0.5 ), similarly, smectic phases occur at lower temperature, cholesteric phases at the intermediate temperature range, and isotropic phases at higher temperature. The figure clearly indicates that the smectic phases obtained from our simulations have a local hexagonal packing structure, which has similarities to that of the smectic-B phase.
Ijms 19 02715 g007
Figure 8. Temperature dependence of the averaged local density for (a) achiral ( c = 0.0 ) and (b) chiral ( c = 0.5 ) liquid-crystal systems with dipole interaction. For achiral systems, nematic and isotropic phases occur for μ = 0.5 and 1.0 , and smectic phases occur in addition to these two phases for μ = 1.5 and 2.0 . For chiral systems, on the other hand, cholesteric and isotropic phases occur for all dipole conditions. For both achiral and chiral systems, density increases with dipole strength.
Figure 8. Temperature dependence of the averaged local density for (a) achiral ( c = 0.0 ) and (b) chiral ( c = 0.5 ) liquid-crystal systems with dipole interaction. For achiral systems, nematic and isotropic phases occur for μ = 0.5 and 1.0 , and smectic phases occur in addition to these two phases for μ = 1.5 and 2.0 . For chiral systems, on the other hand, cholesteric and isotropic phases occur for all dipole conditions. For both achiral and chiral systems, density increases with dipole strength.
Ijms 19 02715 g008
Figure 9. Temperature dependence of the averaged local orientational order parameter P 2 for (a) achiral ( c = 0.0 ) and (b) chiral ( c = 0.5 ) liquid-crystal systems with dipole interactions. Slight decreases in P 2 are observed, which correspond to smectic-nematic phase transitions of achiral system, and significant decreases in P 2 are observed, which corresponds to nematic-isotropic transition of an achiral system or cholesteric-isotropic phase transition of a chiral system.
Figure 9. Temperature dependence of the averaged local orientational order parameter P 2 for (a) achiral ( c = 0.0 ) and (b) chiral ( c = 0.5 ) liquid-crystal systems with dipole interactions. Slight decreases in P 2 are observed, which correspond to smectic-nematic phase transitions of achiral system, and significant decreases in P 2 are observed, which corresponds to nematic-isotropic transition of an achiral system or cholesteric-isotropic phase transition of a chiral system.
Ijms 19 02715 g009
Figure 10. Temperature dependence of the averaged translational order parameter τ for (a) achiral ( c = 0.0 ) and (b) chiral ( c = 0.5 ) liquid-crystal systems with dipole interactions. Significant decreases in τ are observed, which correspond to smectic-nematic phase transitions of achiral systems. The figure clearly indicates that τ increases with dipole strength, especially for nematic and cholesteric at low temperature. τ for cholesteric phases with μ = 2.0 at 1.5 T 1.6 is appreciably larger than others.
Figure 10. Temperature dependence of the averaged translational order parameter τ for (a) achiral ( c = 0.0 ) and (b) chiral ( c = 0.5 ) liquid-crystal systems with dipole interactions. Significant decreases in τ are observed, which correspond to smectic-nematic phase transitions of achiral systems. The figure clearly indicates that τ increases with dipole strength, especially for nematic and cholesteric at low temperature. τ for cholesteric phases with μ = 2.0 at 1.5 T 1.6 is appreciably larger than others.
Ijms 19 02715 g010
Figure 11. Temperature dependence of the averaged local six-fold bond order parameter B 6 for (a) achiral ( c = 0.0 ) and (b) chiral ( c = 0.5 ) liquid-crystal systems with dipole interactions. Significant decreases in B 6 are observed, which correspond to smectic-nematic phase transitions of achiral systems.
Figure 11. Temperature dependence of the averaged local six-fold bond order parameter B 6 for (a) achiral ( c = 0.0 ) and (b) chiral ( c = 0.5 ) liquid-crystal systems with dipole interactions. Significant decreases in B 6 are observed, which correspond to smectic-nematic phase transitions of achiral systems.
Ijms 19 02715 g011
Figure 12. Representative snapshots of the chiral ( c = 0.5 ) with μ = 2.0 : (a) cholesteric phase with short-range smectic layer at T = 1.5 ; (b) cholesteric phase at T = 1.7 , and (c) isotropic phase at T = 2.0 . The yellow head shows dipole direction.
Figure 12. Representative snapshots of the chiral ( c = 0.5 ) with μ = 2.0 : (a) cholesteric phase with short-range smectic layer at T = 1.5 ; (b) cholesteric phase at T = 1.7 , and (c) isotropic phase at T = 2.0 . The yellow head shows dipole direction.
Ijms 19 02715 g012

Share and Cite

MDPI and ACS Style

Nozawa, T.; Brumby, P.E.; Yasuoka, K. Effect of Central Longitudinal Dipole Interactions on Chiral Liquid-Crystal Phases. Int. J. Mol. Sci. 2018, 19, 2715. https://doi.org/10.3390/ijms19092715

AMA Style

Nozawa T, Brumby PE, Yasuoka K. Effect of Central Longitudinal Dipole Interactions on Chiral Liquid-Crystal Phases. International Journal of Molecular Sciences. 2018; 19(9):2715. https://doi.org/10.3390/ijms19092715

Chicago/Turabian Style

Nozawa, Takuma, Paul E. Brumby, and Kenji Yasuoka. 2018. "Effect of Central Longitudinal Dipole Interactions on Chiral Liquid-Crystal Phases" International Journal of Molecular Sciences 19, no. 9: 2715. https://doi.org/10.3390/ijms19092715

APA Style

Nozawa, T., Brumby, P. E., & Yasuoka, K. (2018). Effect of Central Longitudinal Dipole Interactions on Chiral Liquid-Crystal Phases. International Journal of Molecular Sciences, 19(9), 2715. https://doi.org/10.3390/ijms19092715

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop