Liquid crystal materials offer many opportunities for applications such as data storage , image processing , or optical modulators  calling for media that are sensitive to external stimuli. Today, they are considered to have great potential for the realization of so called "smart materials" through micro- to nano-scale patterning techniques. One interesting material in the flourishing group of smart materials is the Liquid Crystal Elastomer (LCE) which combines the properties of liquid crystals and polymers. An elastomer is formed by a weakly crosslinked network of polymers, which confers it a high elasticity. The polymer network has maximum entropy in his undistorted state, and as a result, it resists deformation . As for liquid crystals, these are materials which present phases intermediate between crystalline and isotropic called mesophases, the molecules responsible for this property are named mesogens.
In this chapter, we will present our results on the optical microstructuring of a LCE and the monitoring of its elastic properties. The liquid crystal based elastomer can be considered as an artificial muscle material. The designation of artificial muscle is attributed to soft actuators which can have muscle like behavior. A muscle is in reality an energy converter which converts chemical energy to mechanical motion. The actuators are converters which accept different types of energy and deliver a mechanical quantity like a displacement, a tension,
The LCE used in this study is a nematic thermotropic (sensitive to thermal stimuli) elastomer. A nematic elastomer can be described as crosslinked polymer chains with incorporated mesogens corresponding to rigid rod-like units which can order nematically. The average shape of the material is coupled to the molecular orientation order : the nematic-isotropic phase transition results in the change of shape of the macroscopic sample. When the elastomer is judiciously prepared to obtain a monodomain sample , meaning that all the mesogens are oriented uniformly, the sample is elongated in the nematic phase where all the mesogens are parallel to each other, and has a contracted shape in the isotropic phase where the nematic order is lost.
First, we will show that advantage can be taken of the crosslinking properties of the material to inscribe refractive index modulations by optical means as in a common photoresist . We present two approaches based respectively on one- and two-photon photopolymerization, that can be implemented using an optical microscopy setup. Such a method provides an accurate way to locally modify the properties of the material. Compared to its one-photon relative, the two-photon excitation mode will yield finer patterns because of its higher intrinsic spatial resolution. This technique relies on the fact that the simultaneous absorption of two low energy photons is equivalent to the absorption of a high energy photon with an energy equal to the sum of these two photons energy. The two-photon absorption (TPA) process is a quadratic phenomenon: the probability of the simultaneous absorption of two photon is not proportional to the photon flux as in the case for the linear absorption, but to the square of the flux. TPA is then significant only for high incident fluences. Thus, the use of ultra-short pulse lasers combined with microscopy techniques allows light-matter interaction to be confined within a micrometric volume around the focal point.
After the presentation of the microstructuring process, we will describe how the contractile properties of the material can be measured by detecting the thermally induced step size changes of an inscribed grating. Finally, these results will be interpreted with the help of Molecular Dynamics (MD) simulations. MD is a powerful tool to elucidate the structure and the behavior of the molecules. It has been used to simulate the behavior of polymers  and has proved its efficiency to describe the liquid crystalline phases , thus it also allows the description of elastomer behavior . MD consists of a computational method which calculates time dependent atomic motions by applying the laws of classical mechanics. We have used a combination of Lennard-Jones and Gay-Berne potentials to represent the anisotropic mesogens and the calculations have been carried out in the Parrinello-Rahman-Nosé-Hoover ensemble.
2.1. Photopolymerizable liquid crystals
LCE are materials combining the elastic properties of the elastomers and the anisotropic properties of liquid crystals. They are formed by reticulated liquid crystalline polymers, which, in general, at low reticulation degrees, conserve the mesomorphic properties of the polymer before reticulation as well as the nature of the phases and the transition temperatures. These elastomers, when subjected to external stimuli like temperature or electromagnetic fields, can mimic the actions of muscles (contraction/extension). They are, thereby, considered as artificial muscle materials. Different types of actuators have been developed using polymers as parent material, like for example polymer gels   , conductive polymers   , carbon nanotubes     or dielectric elastomers . Soft and highly flexible, the polymers are well adapted to be used as artificial muscle materials due to their high strength, deformation capabilities, and their capacity to keep intact their properties after several operating cycles. Recently, LCE took over the field of artificial muscles. They are robust and don’t need any solvent to operate (unlike the gels or conductive polymers). The concept of LCE was proposed by P. G. de Gennes in 1975 . In his studies on reticulated liquid crystal polymers, he mentioned the possibility of the material deformation without constraint. In 1981 Frinkelmann
Indeed, according to the insertion topology of the mesogen (liquid crystal phase-forming unit) in the elastomer, there exist two main families of LCEs: side-chain and main-chain elastomers. In the case of main-chain elastomers, the mesogens are directly integrated in the polymer chain, while in the case of side-chain elastomers, they are laterally attached to the polymer chain, either orthogonal to the chain, attached by one end (end on), or parallel to the chain, attached by the side (side on), as shown in figure 1.
Since the monomer used in this study is thermotropic, we will study the contraction properties of the resulting elastomer as a function of temperature. The studied elastomer is a nematic elastomer. The nematic phase is favored by the side-chain conformation. In nematic crystals, the mesogens have a random position but an orientational order: they are, on average, parallel to each other. This behavior is encouraged by the elongated form of these units.
For a nematic crystal, the amount of order can be characterized by the order parameter defined as the average value of a second-order Legendre polynomial :
where corresponds to the angle between the long axis of the molecule and the nematic director. For a completely ordered phase =1, while =0 corresponds to a completely disordered phase (isotropic). The mesogens interact
In the case of the elastomer, the mesomorphic properties (as the nature of phases of the initial monomers) are conserved. The mesogens attached to the sides of the main chain are free to move. During the nematic-isotropic phase transition the mesogens are disoriented and the main chain drag behind them. These microscopic form changes are transferred to the macroscopic sample. The destruction of the nematic order thus generates a contraction of the system along the mesogens alignment direction and an extension in the orthogonal directions as schematized in figure 3.
2.2. Elastomer thin film preparation
The precursor of the LCE used in this study is a mixture of three compounds: a liquid crystalline acrylate monomer, a crosslinker and an ultraviolet (UV) photoinitiator. The monomer whose structure is represented in figure 4 is the 4’-acryloyloxybutyl 2,5-(4’-butyloxybenzoyloxy)benzoate . The part containing the aromatic chains is the rigid part of the monomer and the attached carbon chains correspond to soft parts. The photoinitiator is the 2-benzyl-2-(dimethylamine)-4’ morpholinobutyrophenone (Irgacure 369) added at a concentration of 1 mol %, the crosslinker is the 1,6-hexanediol diacrylate at a concentration of 10 mol %. The elastomer is obtained by the photopolymerization of the monomers and the crosslinking of the polymeric chains.
Let us recall that in the case of a thermotropic nematic elastomer, the destruction of the nematic order by an increase in temperature will cause the contraction of the material along the nematic director. For the contraction to be maximum, the initial order parameter of the sample should be high. The samples were prepared in several micrometer thick glass cells filled by capillarity with the material in its isotropic phase (above 81.5 C, the nematic-isotropic phase transition temperature). The cell surfaces were coated with rubbed poly-vinyl alcohol to insure the alignment of the liquid crystal moieties by guiding them during the filling process and allowing the production of a fully monodomain aligned sample. Thus the contraction can be unidirectional and along the original mesogens’ alignment. To produce the elastomer the filled cell was slowly cooled down (1 C/min) to the nematic phase (63 C) conserving the alignment of the liquid crystal moieties. The elastomer film was then obtained by radical photopolymerization and photocrosslinking under UV light, fixing in this way the nematic alignment of the mesogens. The resulting monodomain elastomer film is named "Liquid Single Crystal Elastomer" (LSCE). The radical photopolymerization starts with the excitation of the photoinitiator Irgacure 369 which is a benzoyl type chromophore with a large absorbance and a high initiation efficiency in the range of 300-400 nm. The initiator generates a radical (a chemical species possessing one or several unpaired electrons in the outer shell) which can react with the monomer and form the macroradical which is also able to react with other monomers. The propagation of the polymerization proceeds by successive addition of monomers to this macro-radical. The monomers are bound to each other by the CH=CH groups located at the extremities of the molecules. It is the breaking of the double bond that allows the bonding between monomers. The polymerization takes place in such a way as to produce a form of "sheathing" of mesogen around the principal chain, thus avoiding steric hindrance (see figure 5). The crosslinker present in the mixture helps to create a crosslinked network and the elastomer so obtained can be extracted from the cell.
3. Experimental setup
Two different light sources were used to photo-structure the elastomers. We have used 365 nm UV light from an Argon Ion (Ar) laser to polymerize the mixture by one-photon photopolymerization or 780 nm infrared (IR) femtosecond laser pulse for two-photon photopolymerization. The advantage of two-photon absorption is the possibility to spatially confine the process by limiting the light-matter interaction to volumes which can be smaller than 1 . The probability of TPA is very low, its cross section is about 31 orders of magnitude lower than a one-photon absorption cross section. To be efficient TPA therefore necessitates a high density of photons. Thus, we have to spatially and temporally confine the photons. This can be achieved by combining the use of high numerical aperture objectives and ultra-short pulsed lasers such as titanium-sapphire lasers emitting in the near-IR region. When approaching the focal point of the objective, the photon density will sharply increase. The low probability TPA process will thus be automatically confined around this focal point in a volume called "voxel". If the sample is transparent to the IR radiation of the laser, but absorptive for the UV light, only the confined volume around the focal point will be excited, while the rest of the sample will remain transparent to the incident radiation. This microscopy technique is named "two-photon microscopy". Figure 6 presents a comparison between the excitation volumes in the case of conventional microscopy and two-photon microscopy.
From the Rayleigh criterion, in optimal conditions, the radial resolution of a conventional microscope is given by the radius of the Airy disc:
In the case of two-photon microscopy, Webb
The excitation is given by:
The experimental setup used for the microstructurating and the realization of the elastomers is based on a confocal microscope. The cell containing the elastomer precursor mixture is placed on a heating plate mounted on a motorized stage that can execute computer-controlled 3D translations along the X, Y, and Z axes. The observation is done by reflection with a CCD camera. The sample is illuminated in transmission with a white light lamp; a filter is used to stop any UV photons. The sample is also placed between polarizers allowing the monitoring of isotropic and nematic zones through the polarization of transmitted light induced by the mesogen alignment. The photo-patterning of the desired structures is realized in the nematic phase by moving the sample under the focus of the objective with the translation stage.
4. Experimental results
4.1. Photostructuration of the material
We have first studied the creation of a patterned elastomer by photopolymerization initiated by linear or two-photon absorption. Different sets of experiments have been performed depending on the sample alignment state (using treated or untreated glass cells) and the excitation source (UV or IR light). The first set of experiments was made using UV excitation (Argon laser, = 365 nm, power = 8 , objective N.A = 0.45) and a non-aligned nematic monomer sample. In a first step the photopolymerizable mixture was introduced in an untreated cell. Then the cell was heated beyond the isotropic-nematic phase transition temperature (> 81.5 C) for the mixture initially in the form of a powder to liquefy and penetrate by capillarity. The mixture was then cooled down until the nematic phase could form (63 C). A shape representing the capital letter E has been patterned by photopolymerization with the focused UV beam using the translation stage with a displacement velocity of 200 /s. The resulting polymerized and crosslinked part draws the letter "E" as can be seen in figure 8(a). The second set of experiments was performed under the same conditions as described above, except that an aligned nematic monomer sample was used figure 3(b). This polymerization was also successful. The polymerized letter "E" presents thinner (30 ) and smoother lines than those obtained for the non-aligned sample (50-60 ). The patterned lines are undistorted in the aligned sample. Therefore, using aligned nematic samples results in considerably improved micro-pattern formation. In the third set of experiments performed again on aligned samples, the IR source (Ti:Sa laser, pulses duration = 100 fs, repetition rate = 80 MHz, = 780 nm, power = 500 mW, objective N.A = 0.45) was used instead of the Argon laser. The results presented in figure 3(c) show that such two-photon excitation processes produce a smooth pattern with even thinner lines measuring 5 to 10 across. The high resolution achieved in this case is indeed better than UV illumination.
4.2. Obtaining the thermo-active elastomer
In order to stabilize the patterns, a post-photopolymerization of the whole sample can be performed using UV light of weak intensity. Under these conditions, the pattern written by photopolymerization using a high intensity beam is preserved as it can be seen on figure 9. The resulting rubber-like film is then removed from the cell and placed on a heating plate to observe the shape changes as a function of temperature. Figure 10 represents the heating of an elastomer with an inscribed grating of 170 step size. The elastomer lies on a film of glycerol oil to minimize the constraints during the contraction. By heating this sample up from room temperature to 120 C, a clear deformation of the film is observed in the direction of the mesogen alignment.
An elastomer with a circular pattern was also produced to highlight the uniaxial contraction and the reversible aspect of the phenomenon as can be seen in figure 11.
5.1. Contraction properties
Direct observation of the micropatterned elastomer’s contraction can be readily achieved by optical microscopy. Figure 12 shows the temperature dependent deformation of the concentric pattern inscribed in a 13 thick film of elastomer. The value of the width perpendicular to the nematic director (denoted a) shows a small variation (less than 2% of extension) whereas the perpendicular direction (denoted b) experiences about 17% of contraction. This demonstrates that the shape change consists mostly in a contraction along the nematic director and proves that a well aligned monodomain sample was achieved. This contraction factor varies with the sample thickness. About 50% contraction is observed in 30 thick samples but as the thickness decreases, the contraction drops to about 30% for 20 thick samples and to about 17 % for 13 thick samples. This behavior can be attributed to mechanical stress. Indeed, the size of the artificial muscle with respect to the surface of adhesion is larger for a thick film than for a thin film. The contraction forces thus become prevalent (as the samples become thicker) over the constrains resulting from the friction of the film on its substrate. It is noteworthy that no dependency of the rate of contraction on the pattern structure has been observed, suggesting that any diffraction element can be inscribed without altering the physical properties of the material. The contraction of the elastomer starts around 90 C. It stops when the phase becomes isotropic, towards 100 C. Both of these temperatures are higher than the phase transition temperatures of the initial liquid crystal which is 81.5 C. The mesogens being free to move within the elastomer, the thermal characteristics of the mesophases should be preserved. This is, however, only partially true since the mesogens are not completely decoupled from the main chain. However it is reasonable to assume that the longer the spacer, the higher the degree of this decoupling. This aspect may thus bear on the phase transition temperature but there is no doubt that the primary cause behind the observed contraction temperatures lies with the fact that the response of the material to the drop in nematic order is not immediate. Contraction will only occur once the degree of disorientation of the mesogens reaches levels that are enough to overcome the frictional forces between the elastomer and its substrate. In addition, with the monitoring technique that is being put to use, the contraction becomes only apparent beyond a certain level. All this may be the origin of the gap between the phase transition temperature and the temperature at which contraction is being observed.
5.2. Diffraction properties
The contraction of the material can be monitored indirectly by the observation of the diffraction pattern proceeding from a grating inscribed in the elastomer. Temperature induced contraction of the material will result in a change in the grating period which will modify the diffraction pattern. We present the study of the diffraction pattern obtained for a 13 thick sample with a 30 period linear grating using a linearly polarized He-Ne laser beam at 633nm.
Before testing the properties of the obtained gratings, we have examined the surface of the elastomer with a profilometer. In some cases when the sample is removed from the cell, its surface exhibits ridges corresponding to the grating. Such observations are commonly made in photopolymerization experiments . The surface corrugation comes from the relaxation of the mechanical stresses introduced by the cross-linking process and which have been maintained by the glass plates. These ondulations are not observed in two-photon photopolymerized samples where the patterns are much thinner and inscribed inside the volume of the elastomer. In the cases where inscribed gratings led to surface ridges, we have checked that the addition of a compensating liquid did not modify the diffraction properties. Thus we may pretend that we are dealing with refractive index gratings. Index inhomogeneity stems from spatial variations in material density between parts that have been strongly polimerized and the remainder of the sample which has only been subjected to a light post-polymerization. The diffraction regime is indicated by the Klein-Cook parameter , where L is the grating thickness, is the vacuum wavelength of light, is the grating step size, and is the mean refractive index. High values of Q (Q>10) correspond to the Brag regime while low values (Q<1) correspond to the Raman-Nath regime. One usually ascribes the Bragg regime, with essentialy only one diffracted wave being produced, to gratings that are thick with respect to step size and the Raman-Nath regime, with several diffracted waves being produced, to thin gratings. In the case of a grating with a step size of 30 and a thickness of 13 yielding the diffraction pattern shown in figure 14, the Q parameter is about 0.4, placing the diffraction in the Raman-Nath regime as demonstrated by the multiple order diffraction pattern. Figure 15 shows the variation of the angle of diffraction for the +1 order as a function of temperature.
It clearly demonstrates the decrease of the grating step with the temperature rise stemming from the unidirectional contraction accompanying the nematic/isotropic transition. In figure 15, the black dots correspond to measurements of the diffraction angle, while the red dots represent values calculated following the diffraction grating formula from the data shown in figure 13 for the 13 thick sample. The agreement between both sets of data shows the consistency between the grating period and the diffraction angle measurements (discrepancies between data points may come from the fact that they have been gathered over several heating cycles). In conclusion, the contraction induced by the temperature is easily monitored by the widening of the diffraction figure. This allows a feedback on the deformation of the material through direct monitoring of the resulting diffraction.
Since we work with a cross-linked polymer, the contraction/extension process will affect all three directions in space. The initial homogeneous nematic order makes the film behave as a biaxial material with an optical axis oriented along the nematic director. Thus, the material may change the polarization state of the light upon diffraction. To illustrate the birefringence properties of the gratings, we have measured the diffraction efficiency dependence on the orientation of the incident linear polarization. The resulting angular distribution for the polarization of the diffracted light is presented in figure 16b. Here, the simplest configuration is adopted: the incident beam is normal to the surface and its polarization is linear. As one can see in figure 16a, the diffraction efficiency is higher when the direction of the polarization of the incident beam and the grating vector are parallel. The diffracted light intensity has been fitted by :
where the angle corresponds to the rotation of the polarisation and is the light intensity of the first diffracted order. Thus we can determine the parallel and the perpendicular diffraction efficiencies which amount respectively to and .
In the case of figure 16b the incident beam polarization is at 45 of the grating vector. This time a light intensity dependence on the rotation of an analyzer placed after the sample is observed. Equation (6) represents the theoretical function which will fit the intensity of light measured after the analyser.
5.3. Molecular dynamics simulations
Molecular dynamics (MD) simulations are a powerful tool for understanding the properties of a sample in terms of molecular collective phenomena. We used classical MD method to simulate the contraction of the elastomer. First, the monomer molecules composed of 95 atoms was fully optimized by the density functional theory (DFT). The chosen functional was B3LYP, a hybrid functional obtained by linear combination of exchange-correlation energy functionals (Local-density approximations LDA, Generalized gradient approximations GGA) and the Hatree-Fock exchange . The orbital basis set was the contracted Gaussian 6-311** set. All calculations were carried out with NWChem , executed on the IBM Power 6 at the "Institut du Développement et des Ressources en Informatique Scientifique" of Orsay. The obtained minimum energy configuration is represented in figure 17.
This model presents a dipole localized on the mesogen corresponding to the aromatic core. This will allow dipole-dipole interactions between the monomers. Then to simulate the contraction of the elastomer resulting from the collective behavior of the monomers the interaction between pairs of anisotropic rigid mesogenic sites was modeled using Gay-Berne (GB) potential. The Gay-Berne potential corresponds to a modified form of Lennard-Jones potential which can be easily adjusted to modify the shape of the studied system . In that sense it has proven its efficiency to model mesogenic systems   .
The force field modeling intra and inter-molecular contribution to the interaction energy is given by the following expression:
where , , , , are respectively the angles number of the dihedral angles of the numbers, and the number of GB and LJ sites numbers. The variations in dihedral angles , which represent torsions, are modeled by a truncated Fourier series ( coefficients). The alkyne groups of the chain are condensed into single atom sites and represented as spheres. All the chain segment parameters are modeled by resorting to universal force field . , et are respectively Lennard-Jones, Gay-Berne and the mixte Lennard-Jones/Gay-Berne potentials. Bond lengths between the mesogen and the flexible chains were set from
For a separation distance of between two particles, the standard LJ potential is written as:
where corresponds to the energy well depth, is the contact parameter, it corresponds to the separation distance at which the inter-particle potential is zero.represents the orientation dependent interaction energy for two GB particles:
where and are unitary vectors along the molecular axis of GB particles, and correspond respectively to the energy well depth and the contact parameter depending on the relative orientation of the particles. These terms are defined as  :
with and corresponding to the contact parameter along the GB particle length and width and and corresponding to the energy well depth for two GB particles with an "end to end" and a "side by side". is the distance between the two particles i and j, corresponds to the unitary vector along , and are unitary vectors along the molecular axis of the particles i and j. and are adjustable parameters which were set to their canonical values of 1.0 and 2.0, respectively. The ratio was =3.46 and =21.2 meV.corresponds to the potential representing the mixte interaction between LJ and GB sites :
We have considered the molecular dynamics in a modified isothermal-isobaric (NPT) ensemble to allow an anisotropic deformation of the simulation cell. The ensemble used is the Parrinello-Rahman-Nosé-Hoover ensemble based on the Hamiltonian   . Thanks to the introduction of a thermostat and a barostat, this ensemble allows to control the temperature and the constraints acting on the simulation cell. It allows the study of anisotropic deformation of a system as a function of temperature. It is in that sense well adapted to the study of systems containing liquid crystal molecules. The Hamiltonian is given by:
where is the reduced position of the particle in the simulation cell
The forces acting on the system are obtained by deriving the field U. The system’s deformations are obtained by integrating the equations of motion using Beeman’s algorithme :
An autocoherent process is necessary to calculate the conjugated moments. The time derivatives of the moments are calculated with the Parrinello-Rahman-Nosé-Hoover Hamiltonian.
Calculations have been carried out for 100 molecules and for simulation times up to 1 fs. We considered that the molecules are bonded by the CH=CH groups located at their extremity as represented in figure 19.
The bonding of the molecules (the polymerization) takes place in such a way to produce a form of sheathing around the principal chain. Crosslinking between principal chains was modeled by resorting to alkyne chains.
To measure the contraction of the simulated system we have monitored the obtained dimensions of the simulation cell vectors and . In figure 20 and 21 we have represented the results of the simulation and the measured contraction. We can notice a difference between the curvature of the two curves.
Figure 22 represents the behavior of the order parameter obtained by polarized Fourier transform spectroscopy (FTIR) in the case of the same elastomer . If we compare the simulation curve with figure 22 we can see that they have the same apparence. This similarity was to be expected, in fact it is the decrease of the order parameter which generates the contraction of the material. But differences between the behavior of the order parameter and the contraction can be identified: the contraction is not immediate, it starts above a given temperature which is higher than the phase transition temperature. In reality a sufficient contraction force should be achieved before observing the phenomena and these contraction forces should also exceed the friction forces with the glass plate. The fact that the simulation curve looks like the order parameter curve more than the contraction curve encourages us to consider external mechanical forces acting on the elastomer in the future.
We have demonstrated that one photon "UV" photopolymerization as well as two-photon "IR" photopolymerization can be used to microstructure artificial muscle materials made of nematic liquid crystalline elastomers without losing the contraction/extension properties. We have shown that the use of two-photon absorption allows to achieve 65 % greater spatial resolutions. A major advantage of the TPA consists in creating shape-changing volume objects, a property particularly interesting for the domain of microfluidics. We have established the possibility of generating a grating design in the sample which can be used as a step changing grating when subject to a temperature increase. We have shown that the contraction induced by the temperature is easily monitored by the widening of the diffraction figure. We can then consider the use of this kind of grating for temperature adjusted feedback devices. In addition, the birefringence properties of the gratings can open the path for polarization dependent diffractive elements.
Molecular dynamics simulations have been used to describe the contraction of the elastomer as a function of temperature. The proposed model allows the simulation of anisotropic molecules made from a combination of Gay-Berne potentials, representing the rigid mesogenic parts, and Lennard-Jones sites representing the alkyne groups of the flexible chains. Even if experimental and numerical contractions behave somewhat differently as functions of the temperature, a fact that may be attributed in part to the neglect of external mechanical forces, a rather good agreement between the behavior of the simulation curve and the order parameter has been achieved and a reasonable agreement with the contraction curve has been obtained.
Ikeda T. Tsutsumi O. 1995 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based Elastomers 268 1873 1875
Lagerwall S. T. 1999 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersWiley-VCH.
Drzaic P. S. 1999Liquid Crystal Dispersions.World Scientific.
Warner M. Terentjev E. M. 2003 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersOxford University Press, USA.
de Gennes P. G. 1997A semi-fast artificial muscle. CR Acad. Sci. Ser. II B 324 343 348
Yu Y. Ikeda T. 2006Two-Photon Excitation in CaF2Eu2+. Phys. Rev. Lett. 45 5416 5418
Naciri J. Srinivasan A. Jeon H. Nikolov N. Keller P. Ratna B. R. 2003Nematic elastomer fiber actuator. Macromolecules 36 8499 8505
Ikeda Y. Yu T. Mamiya J. 2007 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersAngaw. Chem. Int. Ed, 46: 506.
De Gennes,P.G( 1975Réflexionssurun type de polymèresnématiques. ACR Acad. Sci. Paris, Ser. B, 46 101 103
Kupfer J. Finkelmann H. 1991Nematic liquid single crystal elastomers. Makromol Chem., Rapid Commun. 12 717 726
Lessard R. A. Gurusamy M. 1995Photoreactive Polymers in Advance Applications. Chapman and Hall, New York.
Mc Cammon J. A. Gelin B. R. Karplus M. 1977Dynamics of folded proteins. Nature. 267 585 590
Allen M. P. Warren M. A. Wilson M. R. Sauron A. Smith W. 1996 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersJ. Chem. Phys. 105 2850 2858
Darinskii A. A. Zarembo A. Balabaev N. K. 2007 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based Elastomers 252 101 109
Osada Y. Okuzaki H. Hori H. 1992 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based Elastomers 355 242 244
Osada Y. De Rossi D. E. (200 2000 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersSpringer Berlin.
Osada Y. Khokhlov A. R. (200 2002 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersMarcel Dekker, Inc, 270 Madison Avenue, New York, NY 10016, USA, 1 381
Osada Y. Khokhlov A. R. (199 1996Conducting polymer artificial muscles. Synth. Met. 78 339 354
and Fadeev, A.G. and Qi, B. and Smela, E. and Mattes, B.R. and Ding, J. and Spinks, G.M. and Mazurkiewicz, J. and Zhou, D. andWallace, G.G. and others ( Lu W. Fadeev A. G. Qi B. Smela E. Mattes B. R. Ding J. Spinks G. M. Mazurkiewicz J. Zhou D. and Wallace. G. 2002Use of ionic liquids for pi-conjugated polymer electrochemical devices. Science 297: 983.
Pyo M. Bohn C. C. Smela E. Reynolds J. R. Brennan A. B. (200 2003 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersChem. Mater. 15 916 922
and Cui, C. and Zakhidov, A.A. and Iqbal, Z. and Barisci, J.N. and Spinks, G.M. and Wallace, G.G. and Mazzoldi, A. and De Rossi, D. and Rinzler, A.G. and others( Baughman R. H. Cui C. Zakhidov A. A. Iqbal Z. Barisci J. N. Spinks G. M. Wallace G. G. Mazzoldi A. De Rossi D. Rinzler A. 1999Carbon nanotube actuators. Science 284: 1340.
Kim P. Lieber C. M. (199 1999 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based Elastomers 2148 EOF
Zhang Y. Iijima S. 1999 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersPhys. Rev. Lett. 82 3472 3475
Spinks G. M. Wallace G. G. Fifield L. S. Dalton L. R. Mazzoldi A. De Rossi D. Khayrullin I. I. Baughman R. H. 2002Pneumatic carbon nanotube actuators. Adv. Mater. 14 1728 1732
Pelrine K. R. Kornbluh R. Pei Q. Joseph J. (200 2000High-speed electrically actuated elastomers with strain greater than 100%. Science 287: 836.
Finkelmann H. Kock H. J. Rehage G. (198 1981Investigations on liquid crystalline polysiloxanes 3. Liquid crystalline elastomers-a new type of liquid crystalline material. Makromol. Chem., Rapid Commun. 2 317 322
Bergmann G. H. F. Finkelmann H. Percec V. Zhao M. 1997Liquid-crystalline main-chain elastomers. Macromol. Rapid Commun. 18 353 360
Zannoni C. 1979The molecular physics of liquid crystals. Rijeka: G. R. Luckhurst and G. W. Gray.
Thomsen D. L. Keller P. Naciri J. Pink R. Jeon H. Shenoy D. Ratna B. R. (200 2001Liquid crystal elastomers with mechanical properties of a muscle. Macromolecules 34 5868 5875
Zipfel W. R. Williams R. M. Webb W. W. (200 2003 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based Elastomers 21 1369 1377
Zipfel Naydenova. I. Mihaylova E. Martin S. Toal V. 2005Holographic patterning of acrylamide-based photopolymer surface. Opt. Express 13 4878 4889
( Becke A. D. (199 1993)Density-functional thermochemistry. III. The role of exact exchange. Chem. Phys. 98 5648 5652.
Valiev M. Bylaska E. J. Govind N. Kowalski K. Straatsma T. P. van Dam H. J. J. Wang D. Nieplocha J. Apra E. Windus T. L. de Jong W. A. 2010 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersComput. Phys. Commun. 181: 1477.
Gay J. G. Berne B. J. 1981Modification of the overlap potential to mimic a linear site-site potential. J. Chem. Phys. 74: 3316.
Miguel E. Rull L. F. Chalam M. K. Gubbins K. E. Van Swol F. 1991Location of the isotropic-nematic transition in the Gay-Berne model. Mol. Phys. 72 593 605
Luckhurst G. R. Stephens R. A. Phippen R. W. 1990 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersXIX. Mesophases formed by the Gay-Berne model mesogen. Liq. Cryst. 8: 451.
Smith( Allen M. P. Warren M. A. Wilson M. R. Sauraon A. w. 1996 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based Elastomersem. Phys. 105 2850 2858
Rappe, AK and Casewit, CJ and Colwell, KS and Goddard Iii, WA and Skiff, WM( 1999) UFF, UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations. JOSA 114 10024 10035.
Smith( Allen M. P. Warren M. A. Wilson M. R. Sauraon A. w. 1985 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersPhys. 55: 549.
Cleaver D. J. Care C. M. Allen M. P. Neal M. P. 1996 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersPhys. Rev. E 54 559 567
Parrinello M. Rahman A. 1981 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersJ. Appl. Phys. 52: 7182 EOF 7190 EOF
Nose S. 1984A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 81: 511.
Hoover W. G. 1985 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersRev. A 31 1695 1697
Beeman D. 1976 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersJ. Comput. Phys. 20: 130 EOF 139 EOF
Li M. H. Keller P. Li B. and Wang. X. Brunet M. 2003 Monitoring the Contractile Properties of Optically Patterned Liquid Crystal Based ElastomersAdv. Mater. E 15 5569 572