Molecular Dynamic Simulation of Side-Chain Liquid Crystalline Elastomer Under Load

Summary: We present a molecular dynamic simulation of a side chain liquid crystalline elastomer (LCE) under load. The LCE is composed of a flexible tetrafunctional diamond like network with rod-like mesogens attached to the network. As a precursor of the LC elastomer a flexible polymer network in a low molecular liquid-crystal (LC) solvent was used. The phase behavior of the LCE under uniaxial stretching up to the deformations of λ = 1.5 and 2.0 at different densities was studied. As in the non-stretched case upon density increase an isotropic to nematic phase transition occurs. However, in contrast to thermotropic side chain LC elastomers the stress induced shift transition is not observed. The stretching slightly increases the anisotropy of translational diffusion of mesogens in the nematic state. The stress-strain dependence for LCE both in the isotropic and the nematic states is obtained. Elastic modulus increases at high values of order parameter.     

Rheology of a thermotropic polyester in the nematic and isotropic states

The rheological properties of a thermotropic polyester were determined in the nematic and isotropic states. In the isotropic state, the viscosity is almost constant and the polymer is only slightly elastic. The nematic phase has a lower viscosity than the isotropic, except at low frequencies or shear rates, where the viscosity increases as though the polymer had a yield stress. There is a marked dependence of the rheology on shear history. The effects of shearing can be erased by returning the material first to the isotropic state and then back to the nematic state. The results are discussed with reference to analogous observations in small-molecule liquid crystals and in thermotropic aromatic co-polyesters.     

Orientation and anchoring effects in stretched polymer dispersed nematic liquid crystals

The induction of liquid crystal orientation through mechanical stretching was investigated for polymer dispersed liquid crystals (PDLCs) by means of infrared dichroism. Using a nematic liquid crystal BL006 and polyacrylic acid as the polymer matrix, it was possible to stretch the PDLC films with BL006 in either the isotropic or the nematic phase. After cooling the films under strain to room temperature, the molecular orientation of BL006 was found to be much higher for films that contained isotropic liquid droplets of BL006 at the time of stretching than for films that had nematic droplets. Stretching PDLC films with isotropic droplets results in no molecular orientation, but the orientation is induced during the subsequent cooling when BL006 goes through the isotropic-to-nematic phase transition. Interestingly for PAA/BL006, the nematic director orients along the long axes of the elongated droplets despite liquid crystal anchoring perpendicular to the polymer interface.     

Orientation and anchoring effects in stretched polymer dispersed nematic liquid crystals

The induction of liquid crystal orientation through mechanical stretching was investigated for polymer dispersed liquid crystals (PDLCs) by means of infrared dichroism. Using a nematic liquid crystal BL006 and polyacrylic acid as the polymer matrix, it was possible to stretch the PDLC films with BL006 in either the isotropic or the nematic phase. After cooling the films under strain to room temperature, the molecular orientation of BL006 was found to be much higher for films that contained isotropic liquid droplets of BL006 at the time of stretching than for films that had nematic droplets. Stretching PDLC films with isotropic droplets results in no molecular orientation, but the orientation is induced during the subsequent cooling when BL006 goes through the isotropic-to-nematic phase transition. Interestingly for PAA/BL006, the nematic director orients along the long axes of the elongated droplets despite liquid crystal anchoring perpendicular to the polymer interface.     

Phase equilibria and photopolymerisation-induced phase transitions of mesogenic diacrylate monomer and low molecular mass liquid crystal mixture

Experimental phase diagrams of binary mesogenic mixtures of reactive mesogenic diacrylate (RM257) monomer and low molar mass liquid crystals (E7) were determined by means of differential scanning calorimetry and optical microscopy. The combined free energy densities of Flory–Huggins for liquid–liquid demixing, Maier–Saupe for nematic ordering, and phase field free energy for crystal solidification was proposed to describe the phase diagrams of the starting E7/RM257 mixtures. The phase diagram thus constructed is an ideal mixing type, exhibiting a narrow loop of isotropic + nematic (I + N) coexistence region followed by the crystal + nematic (Cr1 + N) region in descending order of temperature. Of particular interest is the permanent fixation of the mesophase structures upon photopolymerisation of neat RM257 in the corresponding nematic and crystalline phases. Upon photopolymerisation of a low RM257 content mixture in both isotropic and nematic states, the nematic–isotropic transition of E7 was found to persist. The permanent structural anchoring is seen upon photo-curing of the 90/10 RM257/E7 mixture in the crystalline state.     

Rheological behavior of thermotropic flexible main-chain polyesters based on 4,4′-dioxy-2,2′-dimethyl azoxybenzene mesogens. I. Oscillatory shear behavior

The rheological properties of fractionated thermotropic nematic polyesters 4,4′-dioxy-2,2′-dimethyl azoxybenzenenonanediyl (n = 7) (AZA-9) and dodecanediyl (n = 10) (DDA-9) have been determined under oscillatory shear in the nematic and isotropic phase. A markedly more pronounced shear thinning characterizes the polydisperse DDA-9 polymer as compared to the fractionated polymer in the nematic state. This difference becomes less pronounced in the isotropic phase. In all cases, the elastic compliance J′ in the isotropic state is much lower than in the nematic state. The values of relaxation time λ, although lower in the isotropic state, are not very different in the nematic state. The loss tangent tan δ is significantly higher in the isotropic phase for lower molecular weights. For molecular weights of 13,000 and above, this difference becomes smaller and shifts to lower angular frequencies. For AZA-9 and DDA-9, η→₀ the limit of the dynamic viscosity can be represented by η^* = M with a value of α similar for both phases and approximately equal to 4.     

Sol/gel and isotropic/nematic transitions in aqueous suspensions of natural nontronite clay. Influence of particle anisotropy. 1. Features of the i/n transition

The phase behavior of a natural nontronite clay was studied for size-selected particles by combining osmotic pressure measurements, visual observations under polarized light, and rheological experiments. In parallel, the positional and orientational correlations of the particles were analyzed by small-angle X-ray scattering. Aqueous suspensions of nontronite exhibit a true isotropic/nematic (I/N) transition that occurs before the sol/gel transition, for ionic strengths below 10(-3) M/L. In this region of the phase diagrams, the system appears to be purely repulsive. The I/N transition shifts toward lower volume fractions for increasing particle anisotropy, and its position in the phase diagram agrees well with the theoretical predictions for platelets. SAXS measurements reveal the presence of characteristic interparticular distances in the isotropic, nematic, and gel phases. The swelling law (separation distance vs swelling law) exhibits two regimes. For high volume fractions, the swelling law is one-dimensional as in layered systems and reveals the presence of isolated platelets. At lower volume fraction, distances scale as phi(-1/3), indicating isotropic volumic swelling. Finally, the experimental osmotic pressure curves can be satisfactorily reproduced by considering the interparticle distances between two charged planes whose effective charge is around 10% of the structural charge.     

Swelling behaviour of poly(butadiene) gels in liquid crystal solvents

The swelling behaviour of poly(butadiene) gels in four different nematogenic liquid crystalline (LC) solvents has been investigated as a function of temperature (T). Microscopy with crossed polarizers reveals that the nematic to isotropic phase transition temperature of the LC solvents inside the gels (T_NI ^g) is slightly lower than that of the surrounding pure LC solvents (T_NI ^o), but the degrees of depression in T_NI ^g in each system are comparable regardless of the considerable differences in the degrees of equilibrium swelling (Q) at T_NI ^g between the various systems. In general, Q in the isotropic phase is larger than that in the nematic phase, but a unique swelling behaviour of the gel is found in the vicinity of T_NI due to the phase transition of the LC solvents. Q remains constant in the temperature range of T_NI ^g ≤ T ≤ T_NI ^o in which the phases of the LCs outside or inside the domain of the gels are different, namely, nematic and isotropic phase, respectively. In addition, a finite abrupt (discontinuous-like) change in Q is observed at around T_NI. The gels swollen in the LCs, having an ability to interact with the crosslinking points via hydrogen bonding, show a significant thermal hysteresis for the temperature dependence of Q in the vicinity of T_NI, while no discernible thermal hysteresis is observed for the gels in the LCs incapable of forming hydrogen bonds.     

Density functional theory study of the nematic-isotropic transition in an hybrid cell

We have employed the density functional theory formalism to investigate the nematic-isotropic capillary transitions of a nematogen confined by walls that favor antagonist orientations to the liquid crystal molecules (hybrid cell). We analyze the behavior of the capillary transition as a function of the fluid-substrate interactions and the pore width. In addition to the usual capillary transition between isotropiclike to nematiclike states, we find that this transition can be suppressed when one substrate is wet by the isotropic phase and the other by the nematic phase. Under this condition the system presents interfacelike states which allow us to continuously transform the nematiclike phase to the isotropiclike phase without undergoing a sharp phase transition. Two different mechanisms for the disappearance of the capillary transition are identified. When the director of the nematiclike state is homogeneously planar-anchored with respect to the substrates, the capillary transition ends up in a critical point. This scenario is analogous to the observed in Ising models when confined in slit pores with opposing surface fields which have critical wetting transitions. When the nematiclike state has a linearly distorted director field, the capillary transition continuously transforms in a transition between two nematiclike states.     

Structure formation of liquid-crystalline isotropic AB block copolymers in nematic solvents

AB block copolymers with an isotropic and a nematic block are synthesized by living anionic copolymerization of isotropic and mesogenic methacrylates. The complete phase diagrams of the corresponding homopolymers with the nematic solvent have been established. In dilute solution (0.2 to 5% by weight) in a low-molar-mass nematic solvent, the copolymers exhibit two different types of structure formation. Polymers with high isotropic block fraction form irreversible networks as a results of phase separation followed by vitrification of the polymer-rich phase. Polymers with high nematic block fraction on the other hand form thermoreversible threadlike structures perfectly ordered parallel to the director of the nematic solvent. These structures cannot be explained by a macroscopic phase separation process but more likely by a formation of micellar-like aggregates.     

The electro-optical response of nematic emulsions

Switchable nematic emulsions are micron-sized droplets of nematic liquid crystal, floating in isotropic fluid matrices. Such droplets can be switched from an opaque (off) to a transparent (on) state by application of very low electric fields. It is known that the electro-optical properties of liquid crystal dispersions are affected by several parameters, including the liquid crystal loading. The electro-optical response of nematic emulsions has been investigated as a function of liquid crystal weight percentage. Almost transparent films with a reduced contrast ratio are obtained with lower liquid crystal contents. A macroscopic phase separation is observed when liquid crystal content exceeds 45 wt %. On the contrary, large contrast ratios and very low switching fields can be obtained if liquid crystal ranges from 25 to 35 wt %. Consequently, nematic emulsions prepared in this liquid crystal range can be used as promising systems for electro-optical applications. In addition to technological developments, these results can help computational and basic studies of phase separation in novel multiphase liquid crystalline materials.     

Conformational analysis of the structure of <Emphasis Type="Italic">p</Emphasis>-butoxybenzylidene-<Emphasis Type="Italic">p</Emphasis>′-propionyloxyaniline at temperatures of crystal → nematic and nematic → isotropic liquid phase transitions

The molecular structure and conformational properties of the p-butoxybenzylidene-p′-propionyloxyphenylaniline molecule in the crystalline state and at temperatures of the crystal → nematic and nematic → isotropic liquid phase transitions were examined by AM1 calculations. It was found that the nematic → isotropic liquid phase transition is accompanied by a change in the molecular conformation.     

A density-functional theory study of the confined soft ellipsoid fluid

A system of soft ellipsoid molecules confined between two planar walls is studied using classical density-functional theory. Both the isotropic and nematic phases are considered. The excess free energy is evaluated using two different Ans?tze and the intermolecular interaction is incorporated using two different direct correlation functions (DCF's). The first is a numerical DCF obtained from simulations of bulk soft ellipsoid fluids and the second is taken from the Parsons-Lee theory. In both the isotropic and nematic phases the numerical DCF gives density and order parameter profiles in reasonable agreement with simulation. The Parsons-Lee DCF also gives reasonable agreement in the isotropic phase but poor agreement in the nematic phase.     

Structural and dynamic properties of a new type of discotic nematic compounds

Thermodynamic, structural and dynamical properties of a new type of discotic compounds, a hydrocarbon without any heteroatoms, displaying a nematic discotic phase have been investigated by means of X-ray diffraction, electro-optical relaxation, and calorimetric studies. Of particular interest are the strength of the first order nematic—isotropic phase transition and the nature of the orientational fluctuations in the isotropic phase. The short range positional order was found to be biaxial in both the isotropic and the nematic phase. The isotropic phase displays strong pretransitional effects originating from orientational fluctuations in the neighbourhood of the transition to the nematic phase. The character of these pretransitional effects differs from that found for calamitic systems in that the number of correlated molecules g₂ is extremely large, of the order of 600 at the clearing temperature and the electro-optical relaxation time is very large, caused by the large value of g₂.     

Phenomenological description of the crystal-liquid crystal phase diagram

A new Landau-type phenomenological free energy function to describe the phase diagram of uniaxial anisotropic materials is proposed. The resulting phase diagram includes crystalline, smectic A, nematic and isotropic states, as well as a theoretically possible plastic solid state. Depending on the parameters of the free energy, smectic A phase may or may not be present for a given material. The obtained topologies of the temperature-pressure phase diagram qualitatively agree with results from computer simulations of model liquid crystalline materials.     

Modeling a liquid crystal dynamics by atomistic simulation with an ab initio derived force field

Atomistic molecular dynamics (MD) simulations of 4-n-pentyl 4'-cyano-biphenyl (5CB) have been performed, adopting a specific ab initio derived force field. Two state points in the nematic phase and three in the isotropic phase, as determined in a previous work, have been considered. At each state point, at least 10 ns have been produced, allowing us to accurately calculate single-molecule properties. In the isotropic phase, the values of the translational diffusion coefficient, and even more so the activation energy for the process, agree well with experimental data. Qualitatively, also the dynamic anisotropy of the nematic phase is correctly accounted for. Rotational diffusion coefficients, which describe spinning and tumbling motions, fall well within the range of experimental values. The reorientational dynamics of our model 5CB covers diverse time regimes. The longest one is strongly temperature dependent and characterized by a relaxation time in accord with experimental dielectric relaxation data. Shear viscosity and Landau-de Gennes relaxation times, typically collective variables, reproduce the experimental results very well in the isotropic phase. In the nematic phase, despite a large statistical uncertainty due to the extremely slow relaxation of the correlation functions involved, our simulation yields the correct relative order of the three experimental Miesowicz viscosities.     

Phase Behavior of Poly(n-butyl acrylate) and Poly(2-ethylhexyl acrylate) in Nematic Liquid Crystal E7

Summary: This paper reports on the phase behavior of photochemically crosslinked poly(n-butyl acrylate) (PABu) and poly(2-ethylhexyl acrylate) (PEHA) networks in a nematic liquid crystal (LC) solvent. The swelling properties of these networks were studied as a function of temperature using a low molecular weight LC (LMWLC), which is an eutectic mixture of cyanoparaphenylenes named E7. To obtain different polymer network densities, the ratio of the reactive monomers ABu (or EHA) to the crosslinking agent hexanedioldiacrylate (HDDA) was varied prior to polymerization/crosslinking reactions. Immersion of these networks in an excess of LC solvent allowed for the measurement of size increase by polarizing optical microscopy in terms of temperature. Diameter ratios were calculated considering swollen to dry network states of the samples. The uptake of LMWLC inside the network preferentially takes place around the nematic to isotropic transition temperature of the solvent. Phase diagrams in the concentration-temperature framework were given and discussed as a function of crosslinking degree of the polymer network and temperature and phase behavior of the solvent.     

Corresponding states in the smectic-A phase

In the last few years, experimental data of many different sources has been gathered to show that order parameter curves of essentially all nematic liquid crystals (NLCs) compounds can be reduced to a single and universal curve of the corresponding states; experimental data of an impressive number of different nematic compounds has been linearly rescaled in such a way that all them coalesced in a single curve. The surprising novelty observed on NLCs is that this universal behaviour is not restricted to the neighbourhoods of a critical region, as is usually expected from critical phenomena, but it covers the entire domain of the nematic phase, from the nematic isotropic temperature to the nematic crystalline temperature. The aim of this paper is to study the statistical mechanics fundamentals of this behaviour and show that it is a consequence of the universality of the angular part of the potential of the interaction between the particles composing the liquid crystal (LC). We will illustrate this behaviour presenting evidences that it can be found on the smectic-A phase.     

Interplay between shear flow and elastic deformations in liquid crystals

We study shear flow in liquid crystal cells with elastic deformations using a lattice Boltzmann scheme that solves the full, three-dimensional Beris-Edwards equations of hydrodynamics. We consider first twisted and hybrid aligned nematic cells, in which the deformation is imposed by conflicting anchoring at the boundaries. We find that backflow renders the velocity profile non Newtonian, and that the director profile divides into two regions characterized by different director orientations. We next consider a cholesteric liquid crystal, in which a twist deformation is naturally present. We confirm the presence of secondary flow for small shear rates, and are able to follow the dynamical pathway of shear-induced unwinding, for higher shear rates. Finally, we analyze how the coupling between shear and elastic deformation can affect shear banding in an initially isotropic phase. We find that for a nematic liquid crystal, elastic distortions may cause an asymmetry in the dynamics of band formation, whereas for a cholesteric, shear can induce twist in an initially isotropic sample.     

Interactions between spherical colloids mediated by a liquid crystal: a molecular simulation and mesoscale study

Monte Carlo simulations and dynamic field theory (DyFT) are used to study the interactions between dilute spherical particles, dispersed in nematic and isotropic phases of a liquid crystal. A recently developed simulation method (expanded ensemble density of states) was used to determine the potential of mean force (PMF) between the two spheres as a function of their separation and size. The PMF was also calculated by a dynamic field theory that describes the evolution of the local tensor order parameter. Both methods reveal an overall attraction between the colloids in the nematic phase; in the isotropic phase, the overall attraction between the colloids is much weaker, whereas the repulsion at short range is stronger. In addition, both methods predict a new topology of the disclination lines, which arises when the particles approach each other. The theory is found to describe the results of simulations remarkably well, down to length scales comparable to the size of the molecules. At separations corresponding to the width of individual molecular layers on the particles' surface, the two methods yield different defect structures. We attribute this difference to the neglect of density inhomogeneities in the DyFT. We also investigate the effects of the size of spherical colloids on their interactions.