Semisoft elasticity and director reorientation in stretched sheets of nematic elastomers

Effects of Operating Mechanical Conditions and Polymer Networks of Nematic Elastomers on Photo-Induced Mechanical Performances

Photoresponsive liquid-crystalline elastomers (LCEs) are promising candidates for light-controlled soft actuators. Photoinduced stress/strain originates from the changes in mechanical properties after light irradiation. However, the correlation between the photoinduced mechanical performance and in-use conditions such as stress/strain states and polymer network properties (such as effective crosslink density and dangling chain density) remains unexplored for practical applications. Here, isometric photo-induced stress or isotonic strain is investigated at different operating strains or stresses, respectively, on LCEs with polymer network variations, produced by different amounts of solvent during polymerization. As the solvent volume increases, the moduli and photoinduced stresses decrease. However, the photo-induced strain, fracture strain, fracture stress, and viscosity increase. The optical response performance initially increases with the operating strain/stress, peaks at a higher actuation strain/stress, and then, decreases depending on the polymer network. The maximum work densities, which also depend on the operating stress, are in the range of ≈200-300 kJm-3. These findings, highlighting the significant variations in the mechanical performance with the operating stress/strain ranges and amount of solvent used in the synthesis, are critical for designing LCE-based mechanical devices.

Surface Instability in a Nematic Elastomer

Liquid crystal elastomers (LCEs) are soft phase-changing solids that exhibit large reversible contractions upon heating, Goldstone-like soft modes, and resultant microstructural instabilities. We heat a planar LCE slab to isotropic, clamp the lower surface, then cool back to nematic. Clamping prevents macroscopic elongation, producing compression and microstructure. We see that the free surface destabilizes, adopting topography with amplitude and wavelength similar to thickness. To understand the instability, we numerically compute the microstructural relaxation of a "nonideal" LCE energy. Linear stability reveals a Biot-like scale-free instability, but with oblique wave vector. However, simulation and experiment show that, unlike classic elastic creasing, instability culminates in a crosshatch without cusps or hysteresis, and is constructed entirely from low-stress soft modes.

Recent Trends in Continuum Modeling of Liquid Crystal Networks: A Mini-Review

This work aims to provide a comprehensive review of the continuum models of the phase behaviors of liquid crystal networks (LCNs), novel materials with various engineering applications thanks to their unique composition of polymer and liquid crystal. Two distinct behaviors are primarily considered: soft elasticity and spontaneous deformation found in the material. First, we revisit these characteristic phase behaviors, followed by an introduction of various constitutive models with diverse techniques and fidelities in describing the phase behaviors. We also present finite element models that predict these behaviors, emphasizing the importance of such models in predicting the material's behavior. By disseminating various models essential to understanding the underlying physics of the behavior, we hope to help researchers and engineers harness the material's full potential. Finally, we discuss future research directions necessary to advance our understanding of LCNs further and enable more sophisticated and precise control of their properties. Overall, this review provides a comprehensive understanding of the state-of-the-art techniques and models used to analyze the behavior of LCNs and their potential for various engineering applications.

Interfacial metric mechanics: stitching patterns of shape change in active sheets

A flat sheet programmed with a planar pattern of spontaneous shape change will morph into a curved surface. Such metric mechanics is seen in growing biological sheets, and may be engineered in actuating soft matter sheets such as phase-changing liquid crystal elastomers (LCEs), swelling gels and inflating baromorphs. Here, we show how to combine multiple patterns in a sheet by stitching regions of different shape changes together piecewise along interfaces. This approach allows simple patterns to be used as building blocks, and enables the design of multi-material or active/passive sheets. We give a general condition for an interface to be geometrically compatible, and explore its consequences for LCE/LCE, gel/gel and active/passive interfaces. In contraction/elongation systems such as LCEs, we find an infinite set of compatible interfaces between any pair of patterns along which the metric is discontinuous, and a finite number across which the metric is continuous. As an example, we find all possible interfaces between pairs of LCE logarithmic spiral patterns. By contrast, in isotropic systems such as swelling gels, only a finite number of continuous interfaces are available, greatly limiting the potential of stitching. In both continuous and discontinuous cases, we find the stitched interfaces generically carry singular Gaussian curvature, leading to intrinsically curved folds in the actuated surface. We give a general expression for the distribution of this curvature, and a more specialized form for interfaces in LCE patterns. The interfaces thus also have rich geometric and mechanical properties in their own right.

Light-Activated Elongation/Shortening and Twisting of a Nematic Elastomer Balloon

Nematic elastomer balloons with inflation-induced axial contraction and shear/torsion effect can be used as actuators for soft robots, artificial muscles, and biomedical instruments. The nematic elastomer can also generate drastic shape changes under illumination, and thus light can be utilized to activate the deformation of nematic elastomer balloons with huge advantages of being accurate, fast, untethered, and environmentally sustainable without chemical byproducts. To explore light-activated deformation behaviors of the balloon, a phenomenological relationship between light intensity and material parameters describing polymer backbone anisotropy is proposed from experiments, and a theoretical model of an optically-responsive nematic elastomer balloon is established based on the nematic elastomer theory. Various light-activated elongation/shortening and twisting behaviors in the cases of free-standing and axial-loading are presented and their mechanisms are elucidated. The light intensity and initial mesogen angle have great influences on the light-activated deformations including the radius, length, shearing angle and mesogen angle. Light can be easily controlled to trigger rich deformation processes, including elongation/shortening and torsion. The results of this paper are expected to promote the understanding of the light-activated deformation behaviors of the nematic elastomer balloon, and the applications in light-activated actuators and machines.

Nematic liquid crystalline elastomers are aeolotropic materials

Continuum models describing ideal nematic solids are widely used in theoretical studies of liquid crystal elastomers. However, experiments on nematic elastomers show a type of anisotropic response that is not predicted by the ideal models. Therefore, their description requires an additional term coupling elastic and nematic responses, to account for aeolotropic effects. In order to better understand the observed elastic response of liquid crystal elastomers, we analyse theoretically and computationally different stretch and shear deformations. We then compare the elastic moduli in the infinitesimal elastic strain limit obtained from the molecular dynamics simulations with the ones derived theoretically, and show that they are better explained by including nematic order effects within the continuum framework.

Probing the in-plane liquid-like behavior of liquid crystal elastomers

When isotropic solids are unequally stretched in two orthogonal directions, the true stress (force per actual cross-sectional area) in the larger strain direction is typically higher than that in the smaller one. We show that thiol-acrylate liquid crystal elastomers with polydomain texture exhibit an unusual tendency: The true stresses in the two directions are always identical and governed only by the area change in the loading plane, independently of the combination of imposed strains in the two directions. This feature proves a previously unidentified state of matter that can vary its shape freely with no extra mechanical energy like liquids when deformed in the plane. The theory and simulation that explain the unique behavior are also provided. The in-plane liquid-like behavior opens doors for manifold applications, including wrinkle-free membranes and adaptable materials.

Effect of stretching angle on the stress plateau behavior of main-chain liquid crystal elastomers

The equilibrium nonlinear stress-stretch relationships for a monodomain main-chain nematic elastomer (MNE) are investigated by varying the angle between the stretching and initial director axes (θ₀). Angle θ₀ has pronounced effects on the ultimate elongation as well as on the width of the low stress plateau regime (Λ_p) during director rotation, whereas θ₀ has no appreciable effect on the plateau stress (σ_p). In the stretching normal to the initial director (θ₀ = 90°), the plateau end exceeds 200% strain. At oblique angles of 90° > θ₀≥ 35°, Λ_p decreases with decreasing θ₀, whereas the definite plateau regime vanishes at θ₀ < 24°. Wide-angle X-ray scattering and polarized optical microscopy measurements reveal that the director rotates uniformly in the biased direction for the MNE of θ₀°≪ 90°, whereas directors rotating clockwise and counterclockwise are coexistent for θ₀ = 90°. Over the entire plateau regime, the MNEs exhibit pure shear deformation characterized by a Poisson's ratio of zero in the direction of the rotation axis. The Λ_p for the corresponding polydomain NE (PNE) undergoing a transition to the monodomain alignment is smaller than that of the MNE of θ₀ = 90°, while the σ_p values for both NEs are almost similar. The semi-soft elasticity concept satisfactorily explains the effects of θ₀ on Λ_p, and the Λ_p value of the PNE, using a single anisotropy parameter which is evaluated from the degree of thermally induced deformation of MNEs.

Internal constraints and arrested relaxation in main-chain nematic elastomers

Nematic liquid crystal elastomers (N-LCE) exhibit intriguing mechanical properties, such as reversible actuation and soft elasticity, which manifests as a wide plateau of low nearly-constant stress upon stretching. N-LCE also have a characteristically slow stress relaxation, which sometimes prevents their shape recovery. To understand how the inherent nematic order retards and arrests the equilibration, here we examine hysteretic stress-strain characteristics in a series of specifically designed main-chain N-LCE, investigating both macroscopic mechanical properties and the microscopic nematic director distribution under applied strains. The hysteretic features are attributed to the dynamics of thermodynamically unfavoured hairpins, the sharp folds on anisotropic polymer strands, the creation and transition of which are restricted by the nematic order. These findings provide a new avenue for tuning the hysteretic nature of N-LCE at both macro- and microscopic levels via different designs of polymer networks, toward materials with highly nonlinear mechanical properties and shape-memory applications.

Theory of liquid crystal elastomers and polymer networks : Connection between neoclassical theory and differential geometry

In liquid crystal elastomers and polymer networks, the orientational order of liquid crystals is coupled with elastic distortions of crosslinked polymers. Previous theoretical research has described these materials through two different approaches: a neoclassical theory based on the liquid crystal director and the deformation gradient tensor, and a geometric elasticity theory based on the difference between the actual metric tensor and a reference metric. Here, we connect those two approaches using a formalism based on differential geometry. Through this connection, we determine how both the director and the geometry respond to a change of temperature.

Nonlinear photomechanics of nematic networks: upscaling microscopic behaviour to macroscopic deformation

A liquid crystal network whose chromophores are functionalized by photochromic dye exhibits light-induced mechanical behaviour. As a result, the micro-scaled thermotropic traits of the network and the macroscopic phase behaviour are both influenced as light alternates the shape of the dyes. In this paper, we present an analysis of this photomechanical behaviour based on the proposed multiscale framework, which incorporates the molecular details of microstate evolution into a continuum-based understanding. The effects of trans-to-cis photoisomerization driven by actinic light irradiation are first examined using molecular dynamics simulations, and are compared against the predictions of the classical dilution model; this reveals certain characteristics of mesogenic interaction upon isomerization, followed by changes in the polymeric structure. We then upscale the thermotropic phase-related information with the aid of a nonlinear finite element analysis; macroscopic deflection with respect to the wide ranges of temperature and actinic light intensity are thereby examined, which reveals that the classical model underestimates the true deformation. This work therefore provides measures for analysing photomechanics in general by bridging the gap between the micro- and macro-scales.

On the nonlinear behaviour of nematic single crystal elastomers under biaxial mechanic and electrical force fields

A slab of nematic-side-chain-liquid-single-crystal elastomer (NSCLSCE), with the director initially oriented in the z -direction, is subjected to a pair of equal mechanical loads and electrical force fields in the x , y directions. The electric fields tend to make easier the rotation of the director after the application of the mechanical force field. A nonlinear expression for the free energy density is used to obtain the interval of stretching for which the system becomes unstable. However, the elastic energy of the network is assumed to be linear. The stress-strain curves predicted by the model show an unstable zone between two linearly increasing segments. The possibility of bifurcation phenomena has been examined.

Tricritical behavior of soft nematic elastomers

We propose a lattice statistical model to investigate the phase diagrams and the soft responses of nematic liquid-crystal elastomers. Using suitably scaled infinite-range interactions, we obtain exact self-consistent equations for the tensor components of the nematic order parameter in terms of temperature, the distortion and stress tensors, and the initial nematic order. These equations are amenable to simple numerical calculations, which are used to characterize the low-temperature soft regime. We find a peculiar phase diagram, in terms of temperature and the diagonal component of the distortion tensor along the stretching direction, with first- and second-order transitions to the soft phase, and the prediction of tricritical points. This behavior is not qualitatively changed if we use different values of the initial nematic order parameter.

Phase separation and disorder in doped nematic elastomers

We formulate and analyse a model describing the combined effect of mechanical deformation, dynamics of the nematic order parameter, and concentration inhomogeneities in an elastomeric mixture of a mesogenic and an isotropic component. The uniform nematic state may exhibit a long-wave instability corresponding to nematic-isotropic demixing. Numerical simulations starting from either a perfectly ordered nematic state or a quenched isotropic state show that coupling between the mesogen concentration and the nematic order parameter influences the shape and orientation of the domains formed during the demixing process.

Numerical study of stretched smectic-A elastomer sheets

We present a numerical study of stretching monodomain smectic-A elastomer sheets, computed using the finite element method. When stretched parallel to their smectic layer normal the smectic layers are unstable to a transition to a buckled state. We model macroscopic deformations by replacing the microscopic energy with a coarse grained effective free energy that accounts for the fine-scale layer buckling. We augment this model with a term to describe the energy of deforming buckled layers, which is necessary to reproduce the experimentally observed Poisson ratios postbuckling. We examine the spatial distribution of the microstructure phases for various stretching angles relative to the layer normal and for different length-to-width aspect ratios. When stretching parallel to the layer normal the majority of the sample forms a bidirectionally buckled microstructure, except at the clamps where a unidirectionally buckled microstructure is predicted. When stretching at small inclinations to the layer normal the phase of the sample is sensitive to the aspect ratio of the sample, with the bidirectionally buckled phase persistent to large angles only for small aspect ratios. We relate these theoretical results to experiments on smectic-A elastomers.

Molecular simulations elucidate electric field actuation in swollen liquid crystal elastomers

Swollen elastomer liquid crystals undergo significant deformations by application of an electric field perpendicular to their alignment axis, as shown in experiments by Urayama et al. [Urayama K, Honda S, Takigawa T (2006) Macromolecules 39:1943-1949]. Here we clarify this surprising effect at the molecular level using large-scale Monte Carlo simulations of an off-lattice model based on a soft Gay-Berne potential. We provide the internal change of molecular organization, as well as the key observables during the actuation cycle.

Negative Poisson's ratio and semisoft elasticity of smectic-C liquid-crystal elastomers

Models of smectic-C liquid-crystal elastomers predict that they can display soft elasticity, in which the shape of the elastomer changes at no energy cost. The amplitude of the soft mode and the accompanying shears are dependent on the orientation of the layer normal and the director with respect to the stretch axis. We demonstrate that in some geometries the director is forced to rotate perpendicular to the stretch axis, causing lateral expansion of the sample-a negative Poisson's ratio. Current models do not include the effect of imperfections that must be present in the physical sample. We investigate the effect of a simple model of these imperfections on the soft modes in monodomain smectic-C elastomers in a variety of geometries. When stretching parallel to the layer normal (with imposed strain) the elastomer has a negative stiffness once the director starts to rotate. We show that this is a result of the negative Poisson's ratio in this geometry through a simple scalar model.

Mechanical properties of monodomain nematic side-chain liquid-crystalline elastomers with homeotropic and in-plane orientation of the director

We present the first study of the shear mechanical properties of monodomain nematic side-chain liquid-crystal elastomers (SCLCEs) prepared by cross-linking with UV irradiation a nematic side-chain liquid-crystal polymer oriented with an electric or a magnetic field. Their elastic behavior was studied in the dry, swollen and stretched states, in order to check the various theoretical descriptions of these systems. The shear measurements taken on the dry samples show that the shear anisotropy is much smaller than that of the usual twice cross-linked samples oriented by a mechanical stretching of the network formed after the first cross-linking step, demonstrating that the elasticity of the networks strongly depends on the preparation procedure used. The shear experiments performed on the swollen state of these two different types of elastomers reveal that the elasticity of the network is Gaussian for the elastomers oriented with the electric or the magnetic field, and non-Gaussian for the elastomers oriented with the usual stretching procedure. The analysis of the stress-strain curves of both types of elastomers with the neoclassical model based on Gaussian rubber elasticity confirms the Gaussian and non-Gaussian nature of their elasticity. The shear experiments performed as a function of the elongation of the homeotropically oriented elastomer when the shear is applied in a direction parallel to the elongation, do not show the decrease of the associated shear modulus, which is theoretically expected when the strain approaches the threshold value marking the beginning of the elastic plateau. However, the observation of this effect could be prevented by possible small misalignments of the director, as suggested by a calculation presented in one of the theories describing this effect.

Modeling elastic instabilities in nematic elastomers

Liquid crystal elastomers are cross-linked polymer networks covalently bonded with liquid crystal mesogens. In the nematic phase, due to strong coupling between mechanical strain and orientational order, these materials display strain-induced instabilities associated with formation and evolution of orientational domains. Using a three-dimensional finite element elastodynamics simulation, we investigate one such instability, the onset of stripe formation in a monodomain film stretched along an axis perpendicular to the nematic director. In our simulation, we observe the formation of striped domains with alternating director rotation. This model allows us to explore the fundamental physics governing dynamic mechanical response of nematic elastomers and also provides a potentially useful computational tool for engineering device applications.

Elastic response and Ward identities in stressed nematic elastomers

Nematic elastomers exhibit a rich elastic response to external stresses. Of particular interest is the semisoft response of elastomers with an anisotropy direction (z) frozen in by a double cross-linking process. This response is characterized by a stress-strain curve for stresses along x perpendicular to z that rises initially, exhibits a nearly flat plateau between two critical values of strain, and then rises again. This paper explores elastic response in semisoft elastomers as a function of externally applied strain. It derives general Ward identities for elastic moduli and shows that the elastic modulus measuring response to xz shears vanishes at the boundaries of the semisoft plateau whereas moduli measuring response to shears perpendicular to the xz plane do not. It then calculates all relevant moduli in a simple model of elastomers and verifies the general Ward-identity predictions.

Observation of a soft mode of elastic instability in liquid crystal elastomers

In monodomain liquid crystal elastomers a symmetry-breaking locked-in anisotropy causes a semisoft elastic response characterized by a plateau in the stress-strain curve. We show by dynamic light scattering performed as a function of deformation that the relaxation rate of the nematic director fluctuations decreases with strain to a very small value at the onset of the soft elastic response, revealing the existence of a dynamic soft mode. The results are in complete agreement with the theory of semisoft elasticity and allow us to determine all the constants of the model.

Mechanical switching of ferroelectric rubber

At the A to C transition, smectic elastomers have recently been observed to undergo approximately 35% spontaneous shear strains. We first explicitly describe how strains of up to twice this value could be mechanically or electrically induced in Sm-C elastomers by rotation of the director on a cone around the layer normal with an elastic cost dependent on constraints. Second, for typical sample geometries, we give the various microstructures in Sm-C akin to those seen in nematic elastomers under distortions with constraints. It is possible to give explicit results for the nature of the textures. Chiral Sm-C elastomers are ferroelectric. We calculate how the polarization could be mechanically reversed by large, hard, or soft strains of the rubber depending upon sample geometry.

Strain-induced polarization in non-ideal chiral nematic elastomers

Symmetry arguments are advanced that, although ideal chiral nematic elastomers cannot show strain-induced electrical polarization, non-ideal ones can. Phenomenological arguments are then presented, which predict a simple and universal form for the direction and strain dependence of the polarization. A microscopic minimal model is also developed, which predicts the same form. Finally, an example of a polarization–strain curve is calculated for a typical experimental geometry. In this geometry, the polarization is exactly zero at both small and large strains, but pronounced for a large set of intermediate strains corresponding to the strains that cause incremental rotation of the nematic director.

Semisoft elastic response of nematic elastomers to complex deformations

We consider a relaxed semisoft elastomer with its director oriented along the z axis that is first subjected to a large stretch in the x direction then to a slight x-z shear. We give a general argument that in any theory including director rotation, at the onset and end of the director rotation induced by these large stretches, there will be kinks in the stress-large strain curve (forming a stress-strain plateau) and zeros in the x-z shear modulus (C5) associated with small shears imposed on top of the stretches. We then find the analytical forms of the C5 -strain curves for a particular model of semisoftness (arising from compositional fluctuations) and show that it, together with the known stress-strain curve, provides the basis for a strong test of this theory. Finally, we consider the scope for other semisoft models and show that the compositional fluctuations model in fact yielded a generic form, that is, it is the most general quadratic free energy that does not explicitly include a final state direction other than the director. By introducing such additional directions, a large range of alternative models could be developed.

Constitutive model for plasticity in an amorphous polycarbonate

A constitutive model for describing the mechanical response of an amorphous glassy polycarbonate is proposed. The model is based on an isotropic elastic phase surrounded by an SO(3) continuum of plastic phases onto which the elastic phase can collapse under strain. An approximate relaxed energy is developed for this model on the basis of physical considerations and extensive numerical testing, and it is shown that it corresponds to an ideal elastic-plastic behavior. Kinetic effects are introduced as rate-independent viscoplasticity, and a comparison with experimental data is presented, showing that the proposed model is able to capture the main features of the plastic behavior of amophous glassy polycarbonate.

Semisoft nematic elastomers and nematics in crossed electric and magnetic fields

Nematic elastomers with a locked-in anisotropy direction exhibit semisoft elastic response characterized by a plateau in the stress-strain curve in which stress does not change with strain. We calculate the global phase diagram for a minimal model, which is equivalent to one describing a nematic in crossed electric and magnetic fields, and show that semisoft behavior is associated with a broken symmetry biaxial phase and that it persists well into the supercritical regime. We also consider generalizations beyond the minimal model and find similar results.

Strain-induced reorientation and mobility in nematic liquid-crystalline elastomers as studied by time-resolved FTIR spectroscopy

Polarized Fourier transform infrared (FTIR) spectroscopy is employed to study the segmental orientation and mobility of liquid-crystalline elastomers (LCEs) with a monodomain structure in response to external mechanical fields parallel and perpendicular to the initial nematic director. The mean orientation and the molecular order parameter of the different molecular moieties referring to the mesogen, the spacer and the network are analyzed in detail. Parallel stretch leaves the mean orientation of the different molecular moieties and its molecular order parameter nearly uninfluenced. Perpendicular stretch results in a threshold-like dependence: for elongation ratios lambda < or = lambda(c) = 1.3 (10 mol% crosslinker density), respectively lambda < or = lambda(c) = 1.6 (5 mol% crosslinker density) no change of the mean orientation and the molecular order parameters is observed, while for lambda > or = lambda(c) all molecular units reorient and their molecular order parameters are strongly decreased. The present studies give no indications that the reorientational dynamics of the network and the mesogens differ as long as the elongation ratio is smaller than lambda(c).

Evidence of supercritical behavior in liquid single crystal elastomers

Temperature profiles of the first and the second moment of the nematic order parameter distribution function, as determined from the deuteron nuclear magnetic resonance line shapes, as well as heat capacity response, provide support for the supercritical scenario of the nematic-paranematic phase transition in liquid single crystal elastomers. The relative strength of the locked-in internal mechanical field with respect to the critical field can be decreased by swelling the elastomer samples with low molecular mass nematogen. By increasing the concentration of the dopant, critical and below-critical behavior is promoted.