Room Temperature Ionic Liquid-Lithium Salt Mixtures: Optical Kerr Effect Dynamical Measurements

The addition of lithium salts to ionic liquids causes an increase in viscosity and a decrease in ionic mobility that hinders their possible application as an alternative solvent in lithium ion batteries. Optically heterodyne-detected optical Kerr effect spectroscopy was used to study the change in dynamics, principally orientational relaxation, caused by the addition of lithium bis(trifluoromethylsulfonyl)imide to the ionic liquid 1-buty1-3-methylimidazolium bis(trifluoromethylsulfonyl)imide. Over the time scales studied (1 ps-16 ns) for the pure ionic liquid, two temperature-independent power laws were observed: the intermediate power law (1 ps to similar to 1 ns), followed by the von Schweidler power law. The von Schweidler power law is followed by the final complete exponential relaxation, which is highly sensitive to temperature. The lithium salt concentration, however, was found to affect both power laws, and a discontinuity could be found in the trend observed for the intermediate power law when the concentration (mole fraction) of lithium salt is close to chi(LiTf(2)N) = 0.2. A mode coupling theory (MCT) schematic model was also used to fit the data for both the pure ionic liquid and the different salt concentration mixtures. It was found that dynamics in both types of liquids are described very well by MCT.

Molecular dynamics study of orientational cooperativity in water

Recent experiments on liquid water show collective dipole orientation fluctuations dramatically slower than expected (with relaxation time >tation, the self-dipole randomization time tr, which is an upper limit on ta; we find that tr5ta. Third, to check if there are correlated domains of dipoles in water which have large relaxation times compared to the individual dipoles, we calculate the randomization time tbox of the site-dipole field, the net dipole moment formed by a set of molecules belonging to a box of edge Lbox. We find that the site-dipole randomization time tbox2.5ta for Lbox3 , i.e., it is shorter than the same quantity calculated for the self-dipole. Finally, we find that the orientational correlation length is short even at low T.

Dynamics of orientational fluctuations in a homeotropic Fréedericksz transition under a rotating magnetic field

We study the problem of the Fréedericksz transition under a rotating magnetic field by using a dynamical model which incorporates thermal fluctuations into the whole set of nematodynamic equations. In contrast to other geometries, nonuniform textures in the plane of the sample do not appear favored. The proper consideration of thermal noise enables us to describe the dynamics of orientational fluctuations both below and above the shifted instability.

Transient dynamics of orientational fluctuations in the Fréedericksz transition

Transient dynamics of spatial fluctuations of the director field in the pure twist Fréedericksz transition is studied. A nonlinear calculation is presented. Anomalous transient fluctuations are shown. Different stages of evolution and the domain of validity of linear theories are discussed.

Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

A cellulosic liquid crystal pool for cellulose nanocrystals: structure and molecular dynamics at high shear rates

Cellulose and its derivatives, such as hydroxypropylcellulose (HPC) have been studied for a long time but they are still not well understood particularly in liquid crystalline solutions. These systems can be at the origin of networks with properties similar to liquid crystalline (LC) elastomers. The films produced from LC solutions can be manipulated by the action of moisture allowing for instance the development of a soft motor (Geng et al., 2013) driven by humidity. Cellulose nanocrystals (CNC), which combine cellulose properties with the specific characteristics of nanoscale materials, have been mainly studied for their potential as a reinforcing agent. Suspensions of CNC can also self-order originating a liquid-crystalline chiral nematic phases. Considering the liquid crystalline features that both LC-HPC and CNC can acquire, we prepared LC-HPC/CNC solutions with different CNC contents (1,2 and 5 wt.%). The effect of the CNC into the LC-HPC matrix was determined by coupling rheology and NMR spectroscopy - Rheo-NMR a technique tailored to analyse orientational order in sheared systems. (C) 2015 Elsevier Ltd. All rights reserved.

Solvable dynamics in a system of interacting random tops

A new solvable model of synchronization dynamics is introduced. It consists of a system of long range interacting tops or magnetic moments with random precession frequencies. The model allows for an explicit study of orientational effects in synchronization phenomena as well as nonlinear processes in resonance phenomena in strongly coupled magnetic systems. A stability analysis of the incoherent solution is performed for different types of orientational disorder. A system with orientational disorder always synchronizes in the absence of noise.

Nonequilibrium orientational patterns in two-component Langmuir monolayers

A model of a phase-separating two-component Langmuir monolayer in the presence of a photoinduced reaction interconverting two components is formulated. An interplay between phase separation, orientational ordering, and reaction is found to lead to a variety of nonequilibrium self-organized patterns, both stationary and traveling. Examples of the patterns, observed in numerical simulations, include flowing droplets, traveling stripes, wave sources, and vortex defects.

Effect of the position of the double-bond in monounsaturated phospholipids on the dynamics of model membranes

Therllloelynalllics of lllodel 11lel1ll)rane systeills containing 1110nollnsaturatecl I)lloSI)holil) ids is strongly infllienced l)y the I)osition of the C==C dOlll)le })ond in tIle acyl chain. The telllI)eratllres of both chain-nlelting (TM) and La -+ HI! (TH) I)hase traIlsitions are lowered by IIp to 20°C when C==C is Inoved froln positions 6 or 11 to I)osition 9 in an 18-carl)on chain. This work is an attellll)t to ellicidate the uIlderlying Illoleclilar Illechanisllls reSI)Onsi])le for tllese draillatic tllerillodynaillic changes. Mixtllres of di-18: 1 l)hoSI)hatidylethanolanline with C==C at l)ositioIlS 6, 9, 11 were llsed, witll a sI1lall aI1lOlint of I)erdellterated tetradecanol, known to })e a gooel rel)Orter of the cllain Illoleclilar order. SI)ectral second 11I0I1lents were llsed to Illonitor tIle La -+ HII I)hase transition, which was fOllnd to ])e ])road (2-6°C), with a slight llysteresis on heatiIlg/cooling. The orientational order I)rofiles were nleasllred 1lSiIlg 2H Illiclear Illagnetic resonance and changes in these order I)rofiles between La aIld HII I)hases silow l)oth a local increase in order in the vicinity of the C==C bonds and an o\Terall decrease ill the average orientational order of the chain as a whole. These Sll])tle changes recluire })oth high-fidelity SI)ectrosCol)y and a careflll data analysis that takes into aCCOllnt the effects due to l)artiall1lagnetically-indllced orientational ordering of the l)ilayers. In tIle COIltext of SOllle recently rel)Orted cross-relaxation 11leaSlirenlents in Silllilar l)llOSI)llolil)iels, 0111' reslilts sllggest that large-anll)litllde conforlllational changes in the interior of tIle I110del 111eI11])ranes I)lay a 1110re significant role than I)reviollsly thOllght.

Dynamics of shear-induced orientation transitions in block copolymers

This review discusses liquid crystal phase formation by biopolymers in solution. Lyotropic mesophases have been observed for several classes of biopolymer including DNA, peptides, polymer/peptide conjugates, glycopolymers and proteoglycans. Nematic or chiral nematic (cholesteric) phases are the most commonly observed mesophases, in which the rod-like fibrils have only orientational order. Hexagonal columnar phases are observed for several systems (DNA, PBLG, polymer/peptide hybrids) at higher concentration. Lamellar (smectic) phases are reported less often, although there are examples such as the layer arrangement of amylopectin side chains in starch. Possible explanations for the observed structures are discussed. The biological role of liquid crystal phases for several of these systems is outlined. Commonly, they may serve as a template to align fibrils for defined structural roles when the biopolymer is extruded and dried, for instance in the production of silk by spiders or silkworms, or of chitin in arthropod shells. In other cases, liquid crystal phase formation may occur in vivo simply as a consequence of high concentration, for instance the high packing density of DNA within cell nuclei.

Anisotropic model for granulated sea ice dynamics

A continuum model describing sea ice as a layer of granulated thick ice, consisting of many rigid, brittle floes, intersected by long and narrow regions of thinner ice, known as leads, is developed. We consider the evolution of mesoscale leads, formed under extension, whose lengths span many floes, so that the surrounding ice is treated as a granular plastic. The leads are sufficiently small with respect to basin scales of sea ice deformation that they may be modelled using a continuum approach. The model includes evolution equations for the orientational distribution of leads, their thickness and width expressed through second-rank tensors and terms requiring closures. The closing assumptions are constructed for the case of negligibly small lead ice thickness and the canonical deformation types of pure and simple shear, pure divergence and pure convergence. We present a new continuum-scale sea ice rheology that depends upon the isotropic, material rheology of sea ice, the orientational distribution of lead properties and the thick ice thickness. A new model of lead and thick ice interaction is presented that successfully describes a number of effects: (i) because of its brittle nature, thick ice does not thin under extension and (ii) the consideration of the thick sea ice as a granular material determines finite lead opening under pure shear, when granular dilation is unimportant.

Temperature-dependent orientational ordering on a spherical surface modeled with a lattice spin model

We study the orientational ordering on the surface of a sphere using Monte Carlo and Brownian dynamics simulations of rods interacting with an anisotropic potential. We restrict the orientations to the local tangent plane of the spherical surface and fix the position of each rod to be at a discrete point on the spherical surface. On the surface of a sphere, orientational ordering cannot be perfectly nematic due to the inevitable presence of defects. We find that the ground state of four +1/2 point defects is stable across a broad range of temperatures. We investigate the transition from disordered to ordered phase by decreasing the temperature and find a very smooth transition. We use fluctuations of the local directors to estimate the Frank elastic constant on the surface of a sphere and compare it to the planar case. We observe subdiffusive behavior in the mean square displacement of the defect cores and estimate their diffusion constants.

Computer simulation of ordering and dynamics in liquid crystals in the bulk and close to the surface

The aim of this PhD thesis is to investigate the orientational and dynamical properties of liquid crystalline systems, at molecular level and using atomistic computer simulations, to reach a better understanding of material behavior from a microscopic point view. In perspective this should allow to clarify the relation between the micro and macroscopic properties with the objective of predicting or confirming experimental results on these systems. In this context, we developed four different lines of work in the thesis. The first one concerns the orientational order and alignment mechanism of rigid solutes of small dimensions dissolved in a nematic phase formed by the 4-pentyl,4 cyanobiphenyl (5CB) nematic liquid crystal. The orientational distribution of solutes have been obtained with Molecular Dynamics Simulation (MD) and have been compared with experimental data reported in literature. we have also verified the agreement between order parameters and dipolar coupling values measured in NMR experiments. The MD determined effective orientational potentials have been compared with the predictions of Maier­Saupe and Surface tensor models. The second line concerns the development of a correct parametrization able to reproduce the phase transition properties of a prototype of the oligothiophene semiconductor family: sexithiophene (T6). T6 forms two crystalline polymorphs largely studied, and possesses liquid crystalline phases still not well characterized, From simulations we detected a phase transition from crystal to liquid crystal at about 580 K, in agreement with available experiments, and in particular we found two LC phases, smectic and nematic. The crystal­smectic transition is associated to a relevant density variation and to strong conformational changes of T6, namely the molecules in the liquid crystal phase easily assume a bent shape, deviating from the planar structure typical of the crystal. The third line explores a new approach for calculating the viscosity in a nematic through a virtual exper- iment resembling the classical falling sphere experiment. The falling sphere is replaced by an hydrogenated silicon nanoparticle of spherical shape suspended in 5CB, and gravity effects are replaced by a constant force applied to the nanoparticle in a selected direction. Once the nanoparticle reaches a constant velocity, the viscosity of the medium can be evaluated using Stokes' law. With this method we successfully reproduced experimental viscosities and viscosity anisotropy for the solvent 5CB. The last line deals with the study of order induction on nematic molecules by an hydrogenated silicon surface. Gaining predicting power for the anchoring behavior of liquid crystals at surfaces will be a very desirable capability, as many properties related to devices depend on molecular organization close to surfaces. Here we studied, by means of atomistic MD simulations, the flat interface between an hydrogenated (001) silicon surface in contact with a sample of 5CB molecules. We found a planar anchoring of the first layers of 5CB where surface interactions are dominating with respect to the mesogen intermolecular interactions. We also analyzed the interface 5CB­vacuum, finding a homeotropic orientation of the nematic at this interface.

Crystal phases and glassy dynamics in monodisperse hard ellipsoids

We have performed Monte Carlo and molecular dynamics simulations of suspensions of monodisperse, hard ellipsoids of revolution. Hard-particle models play a key role in statistical mechanics. They are conceptually and computationally simple, and they offer insight into systems in which particle shape is important, including atomic, molecular, colloidal, and granular systems. In the high density phase diagram of prolate hard ellipsoids we have found a new crystal, which is more stable than the stretched FCC structure proposed previously . The new phase, SM2, has a simple monoclinic unit cell containing a basis of two ellipsoids with unequal orientations. The angle of inclination is very soft for length-to-width (aspect) ratio l/w=3, while the other angles are not. A symmetric state of the unit cell exists, related to the densest-known packings of ellipsoids; it is not always the stable one. Our results remove the stretched FCC structure for aspect ratio l/w=3 from the phase diagram of hard, uni-axial ellipsoids. We provide evidence that this holds between aspect ratios 3 and 6, and possibly beyond. Finally, ellipsoids in SM2 at l/w=1.55 exhibit end-over-end flipping, warranting studies of the cross-over to where this dynamics is not possible. Secondly, we studied the dynamics of nearly spherical ellipsoids. In equilibrium, they show a first-order transition from an isotropic phase to a rotator phase, where positions are crystalline but orientations are free. When over-compressing the isotropic phase into the rotator regime, we observed super-Arrhenius slowing down of diffusion and relaxation, and signatures of the cage effect. These features of glassy dynamics are sufficiently strong that asymptotic scaling laws of the Mode-Coupling Theory of the glass transition (MCT) could be tested, and were found to apply. We found strong coupling of positional and orientational degrees of freedom, leading to a common value for the MCT glass-transition volume fraction. Flipping modes were not slowed down significantly. We demonstrated that the results are independent of simulation method, as predicted by MCT. Further, we determined that even intra-cage motion is cooperative. We confirmed the presence of dynamical heterogeneities associated with the cage effect. The transit between cages was seen to occur on short time scales, compared to the time spent in cages; but the transit was shown not to involve displacements distinguishable in character from intra-cage motion. The presence of glassy dynamics was predicted by molecular MCT (MMCT). However, as MMCT disregards crystallization, a test by simulation was required. Glassy dynamics is unusual in monodisperse systems. Crystallization typically intervenes unless polydispersity, network-forming bonds or other asymmetries are introduced. We argue that particle anisometry acts as a sufficient source of disorder to prevent crystallization. This sheds new light on the question of which ingredients are required for glass formation.

Complexity of classical dynamics of molecular systems. II. Finite statistical complexity of a water-Na+ system

The computational mechanics approach has been applied to the orientational behavior of water molecules in a molecular dynamics simulated water–Na + system. The distinctively different statistical complexity of water molecules in the bulk and in the first solvation shell of the ion is demonstrated. It is shown that the molecules undergo more complex orientational motion when surrounded by other water molecules compared to those constrained by the electric field of the ion. However the spatial coordinates of the oxygen atom shows the opposite complexity behavior in that complexity is higher for the solvation shell molecules. New information about the dynamics of water molecules in the solvation shell is provided that is additional to that given by traditional methods of analysis.