Diffusion coefficient increases with density in hard ellipsoid liquid crystals

Increase in rod diffusivity emerges even in Markovian nature

Rod-shaped particles embedded in certain matrices have been reported to exhibit an increase in their center of mass diffusivity upon increasing the matrix density. This increase has been considered to be caused by a kinetic constraint in analogy with tube models. We investigate a mobile rodlike particle in a sea of immobile point obstacles using a kinetic Monte Carlo scheme equipped with a Markovian process, that generates gaslike collision statistics, so that such kinetic constraints do essentially not exist. Even in such a system, provided the particle's aspect ratio exceeds a threshold value of about 24, the unusual increase in the rod diffusivity emerges. This result implies that the kinetic constraint is not a necessary condition for the increase in the diffusivity.

Enhancement in the diffusivity of Brownian spheroids in the presence of spheres

In the present paper, we have extended the simulation technique Brownian cluster dynamics (BCD) to analyze the dynamics of the binary mixture of hard ellipsoids and spheres. The shape dependent diffusional properties have been incorporated into BCD using Perrin's factor and compared with analytical results of a one-component ellipsoidal system. We have investigated pathways to enhance the diffusivity of spheroids in the binary mixture by manipulating the phase behavior of the system through varying the fraction of spheres in the binary mixture. We show that at low volume fraction the spherical particles have a higher diffusion coefficient than the ellipsoids due to the higher friction coefficient. However, at a higher volume fraction, we show that the diffusion coefficient of the ellipsoids increases irrespective of the aspect ratio due to the anisotropic shape.

Demixing of two species via reciprocally concentration-dependent diffusivity

We propose a model for demixing of two species by assuming a density-dependent effective diffusion coefficient of the particles. Both sorts of microswimmers diffuse as active overdamped Brownian particles with a noise intensity that is determined by the surrounding density of the respective other species within a sensing radius r_{s}. A higher concentration of the first (second) sort will enlarge the diffusion and, in consequence, the intensity of the noise experienced by the second (first) sort. Numerical and analytical investigations of steady states of the macroscopic equations prove the demixing of particles due to this reciprocally concentration-dependent diffusivity. An ambiguity of the numerical integration scheme for the purely local model (r_{s}→0) is resolved by considering nonvanishing sensing radii in a nonlocal model with r_{s}>0.

Structure and rheology of soft hybrid systems of magnetic nanoparticles in liquid-crystalline matrices: results from particle-resolved computer simulations

Hybrid mixtures composed of magnetic nanoparticles (MNP) in liquid crystalline (LC) matrices are a fascinating class of soft materials with intriguing physical properties and a wide range of potential applications, e.g., as stimuli-responsive and adaptive materials. Already in the absence of an external stimulus, these systems can display various types of orientationally disordered and ordered phases, which are enriched by self-assembled structures formed by the MNPs. In the presence of external fields, one typically observes highly nonlinear macroscopic behavior. However, an understanding of the structure and dynamics of such systems on the particle level has, so far, remained elusive. In the present paper we review recent computer simulation studies targeting the structure, equilibrium dynamics and rheology of LC-MNP systems, in which the particle sizes of the two components are comparable. As a numerically tractable model system we consider mixtures of soft spherical or elongated particles with a permanent magnetic dipole moment and ellipsoidal non-magnetic particles interacting via a Gay-Berne potential. We address, first, equilibrium aspects such as structural organization and self-assembly (cluster formation) of the MNPs in dependence of the orientational state of the matrix, the role of the size ratio, the impact of an external magnetic field, and the translational and orientational diffusion of the two components. Second, we discuss the non-equilibrium dynamics of LC-MNP mixtures under planar shear flow, considering both, spherical and non-spherical MNPs. Our results contribute to a detailed understanding of these intriguing hybrid materials, and they may serve as a guide for future experiments.

Anisotropic viscoelastic phase separation in polydisperse hard rods leads to nonsticky gelation

Spinodal demixing into two phases having very different viscosities leads to viscoelastic networks-i.e., gels-usually as a result of attractive particle interactions. Here, however, we demonstrate demixing in a colloidal system of polydisperse, rod-like clay particles that is driven by particle repulsions instead. One of the phases is a nematic liquid crystal with a highly anisotropic viscosity, allowing flow along the director, but suppressing it in other directions. This phase coexists with a dilute isotropic phase. Real-space analysis and molecular-dynamics simulations both reveal a long-lived network structure that is locally anisotropic, yet macroscopically isotropic. We show that our system exhibits the characteristics of colloidal gelation, leading to nonsticky gels.

Static structure and dynamical behavior of colloidal liquid crystals consisting of hydroxyapatite-based nanorod hybrids

Biominerals such as bones and teeth have elaborate nanostructures composed of aligned anisotropic hydroxyapatite (HAp) nanocrystals, which results in excellent mechanical properties. Construction of such ordered structures of HAp nanocrystals in synthetic materials is challenging. Recently, we reported that HAp-nanorod-based colloidal liquid crystals could be obtained. In the present study, the static structure and dynamics of liquid-crystalline (LC) colloidal dispersions of HAp nanorods are investigated by using small-angle X-ray scattering (SAXS) and X-ray photon correlation spectroscopy (XPCS). The SAXS results reveal that the interparticle distance decreases with increasing HAp concentration, φHAp, and the decrease of the interparticle distance for the short-axis direction is significantly smaller in the LC phase than the interparticle distance in the isotropic phase. In the dynamical studies of the LC phase using XPCS, we observe the diffusive motion of the HAp colloids, with the diffusion coefficient being dependent on the wave number. The diffusive motion slows down with increasing φHAp. We observe anisotropic dynamics after long-term storage (160 days after sealing), whereas only isotropic dynamics are observed in the initial XPCS measurements after short-term storage (14 days after sealing). Moreover, we have found that the dynamics slows down with increasing storage time.

Anomalous transport of magnetic colloids in a liquid crystal-magnetic colloid mixture

We report an extensive molecular dynamics study on the translational dynamics of a hybrid system composed of dipolar soft spheres (DSS), representing ferromagnetic particles, suspended in a liquid crystal (LC) matrix. We observe that the LC matrix strongly modifies the dynamics of the DSS. In the isotropic regime, the DSS show a crossover from subdiffusive to normal diffusive behavior at long times, with an increase of the subdiffusive regime as the dipolar coupling strength is increased. In the nematic regime, the LC matrix, due to the collective reorientation of LC particles, imposes a cylindrical confinement on the DSS chains. This leads to a diffusive dynamics of DSS along the nematic director and a subdiffusive dynamics (with an exponent of ∼0.5) in the perpendicular direction. The confinement provided by the LC matrix is also reflected by the oscillatory behavior of the components of the velocity autocorrelation function of the DSS in the nematic phase.

Investigation of nematic to smectic phase transition and dynamical properties of strongly confined semiflexible polymers using Langevin dynamics

We investigated the nematic to smectic phase transition for strongly confined semiflexible polymer solutions in slit-like confinements using GPU-accelerated Langevin dynamics. We characterized the phase transitions from the nematic to smectic phases for semi-flexible polymer solutions as the polymer density increased. The dependence for the lyotropic nematic to smectic transition can be collapsed by scaling exponents between 0.2 and 0.3. The smectic C phase is found for all the cases with the polymer orientation director tilted with respect to smectic layer lateral alignment. As the chain rigidity increases, the transition density decreases for systems in which the polymer persistence length (P) to slit height (H) ratios are 1.25, 2.5, 3.75, 5 and 25. We also characterized the polymer dynamics for the isotropic-nematic-smectic transitions. The overall polymer diffusivity decreased steadily as the polymer density increased. We observed anomalous polymer diffusion along the nematic director near the isotropic-nematic transition, similar to previously reported behavior for nematic-forming ellipsoids. Polymer diffusivity decreased sharply by two orders of magnitude upon the nematic-smectic transition.

Single molecule translocation in smectics illustrates the challenge for time-mapping in simulations on multiple scales

Understanding the connections between the characteristic dynamical time scales associated with a coarse-grained (CG) and a detailed representation is central to the applicability of the coarse-graining methods to understand molecular processes. The process of coarse graining leads to an accelerated dynamics, owing to the smoothening of the underlying free-energy landscapes. Often a single time-mapping factor is used to relate the time scales associated with the two representations. We critically examine this idea using a model system ideally suited for this purpose. Single molecular transport properties are studied via molecular dynamics simulations of the CG and atomistic representations of a liquid crystalline, azobenzene containing mesogen, simulated in the smectic and the isotropic phases. The out-of-plane dynamics in the smectic phase occurs via molecular hops from one smectic layer to the next. Hopping can occur via two mechanisms, with and without significant reorientation. The out-of-plane transport can be understood as a superposition of two (one associated with each mode of transport) independent continuous time random walks for which a single time-mapping factor would be rather inadequate. A comparison of the free-energy surfaces, relevant to the out-of-plane transport, qualitatively supports the above observations. Thus, this work underlines the need for building CG models that exhibit both structural and dynamical consistency to the underlying atomistic model.

Crowding-facilitated macromolecular transport in attractive micropost arrays

Our study of DNA dynamics in weakly attractive nanofabricated post arrays revealed crowding enhances polymer transport, contrary to hindered transport in repulsive medium. The coupling of DNA diffusion and adsorption to the microposts results in more frequent cross-post hopping and increased long-term diffusivity with increased crowding density. We performed Langevin dynamics simulations and found maximum long-term diffusivity in post arrays with gap sizes comparable to the polymer radius of gyration. We found that macromolecular transport in weakly attractive post arrays is faster than in non-attractive dense medium. Furthermore, we employed hidden Markov analysis to determine the transition of macromolecular adsorption-desorption on posts and hopping between posts. The apparent free energy barriers are comparable to theoretical estimates determined from polymer conformational fluctuations.

The inherent dynamics of isotropic- and nematic-phase liquid crystals

The geodesic (shortest) pathways through the potential energy landscape of a liquid can be thought of as defining what its dynamics would be if thermal noise were removed, revealing what we have called the "inherent dynamics" of the liquid. We show how these inherent paths can be located for a model liquid crystal former, showing, in the process, how the molecular mechanisms of translation and reorientation compare in the isotropic and nematic phases of these systems. These mechanisms turn out to favor the preservation of local orientational order even under macroscopically isotropic conditions (a finding consistent with the experimental observation of pseudonematic domains in these cases), but disfavor the maintenance of macroscopic orientational order, even in the nematic phase. While the most efficient nematic pathways that maintain nematic order are indeed shorter than those that do not, it is apparently difficult for the system to locate these paths, suggesting that molecular motion in liquid-crystal formers is dynamically frustrated, and reinforcing the sense that there are strong analogies between liquid crystals and supercooled liquids.

Abnormal polymer transport in crowded attractive micropost arrays

We investigate polymer diffusion in a quasi-two-dimensional environment decorated with attractive cylindrical posts using Langevin dynamics simulation. We find that the polymer diffusivity has non-monotonic dependence on the post array density. This diffusive behavior strongly depends on the adsorption-desorption transition and the critical adsorption strength ε_c. For ε < ε_c, the polymer undergoes normal diffusion and the diffusivity decreases as the post density increases due to the reduction of the void volume. For ε > ε_c, polymer dynamics is strongly mediated by post adsorption, and we observe a regime where the polymer diffusivity increases as the post density increases. The polymer diffusivity reaches a maximum, which can be attributed to cross-post translation enabled by large polymer conformation fluctuations. We find both cross-post transport and polymer conformation fluctuations strongly depend on the post absorption strength and the chain length.

Diffusion of a nanowire rod through an obstacle field

We report the experimental realization of a rod diffusing in a two-dimensional obstacle field following the single rod dynamics. We use a silver nanowire as our rod and two types of obstacles: repelling light beams and polymer pillars. We study the effect of hydrodynamic interactions on the transport of the rod, comparing both experimental realizations and recent simulations. We propose a framework for analyzing the transport through such systems, and we predict a new superdiffusive regime of rod transport at high obstacle concentration and short times.

Hard-body models of bulk liquid crystals

Hard models for particle interactions have played a crucial role in the understanding of the structure of condensed matter. In particular, they help to explain the formation of oriented phases in liquids made of anisotropic molecules or colloidal particles and continue to be of great interest in the formulation of theories for liquids in bulk, near interfaces and in biophysical environments. Hard models of anisotropic particles give rise to complex phase diagrams, including uniaxial and biaxial nematic phases, discotic phases and spatially ordered phases such as smectic, columnar or crystal. Also, their mixtures exhibit additional interesting behaviours where demixing competes with orientational order. Here we review the different models of hard particles used in the theory of bulk anisotropic liquids, leaving aside interfacial properties and discuss the associated theoretical approaches and computer simulations, focusing on applications in equilibrium situations. The latter include one-component bulk fluids, mixtures and polydisperse fluids, both in two and three dimensions, and emphasis is put on liquid-crystal phase transitions and complex phase behaviour in general.

Two-dimensional system of hard ellipses: a molecular dynamics study

We have simulated the dynamics of a two-dimensional system of hard ellipses by event-oriented molecular dynamics in microcanonical NVE ensemble. Various quantities, namely longitudinal and transverse velocity auto-correlation functions, translational and rotational diffusion mean-squared displacements, pressure, intermediate self-scattering function, radial distribution function, and angular spatial correlation, have been obtained and their dependence on packing fraction is characterized. Despite absence of prominent positional ordering, the orientational degree of freedom behaves nontrivially and exhibits interesting features. Slowing down is observed in the angular part of the motion near isotropic-nematic phase transition. It is shown that above a certain packing fraction the rotational mean-squared displacement exhibits a three-stage temporal regime including a plateau. Comparison to 2D system of hard needles is made and it is shown that from positional viewpoint, the ellipse system is more ordered.

Anisotropic Brownian motion in ordered phases of DNA fragments

Using Fluorescence Recovery After Photobleaching, we investigate the Brownian motion of DNA rod-like fragments in two distinct anisotropic phases with a local nematic symmetry. The height of the measurement volume ensures the averaging of the anisotropy of the in-plane diffusive motion parallel or perpendicular to the local nematic director in aligned domains. Still, as shown in using a model specifically designed to handle such a situation and predicting a non-Gaussian shape for the bleached spot as fluorescence recovery proceeds, the two distinct diffusion coefficients of the DNA particles can be retrieved from data analysis. In the first system investigated (a ternary DNA-lipid lamellar complex), the magnitude and anisotropy of the diffusion coefficient of the DNA fragments confined by the lipid bilayers are obtained for the first time. In the second, binary DNA-solvent system, the magnitude of the diffusion coefficient is found to decrease markedly as DNA concentration is increased from isotropic to cholesteric phase. In addition, the diffusion coefficient anisotropy measured within cholesteric domains in the phase coexistence region increases with concentration, and eventually reaches a high value in the cholesteric phase.

Aspect-ratio-dependent phase transitions and concentration fluctuations in aqueous colloidal dispersions of charged platelike particles

Phase transitions of aqueous colloidal dispersions of charged platelike particles of niobate nanosheets were investigated as a function of the aspect ratio (r(asp)) and particle volume concentration (φ(p)) by means of small-angle neutron scattering and small-angle x-ray scattering. The results elucidated the following three pieces of evidence: (1) the macroscopic phase separation of the dispersions into an isotropic phase and a liquid crystalline (LC) phase under the conditions of (a) varying r(asp) (1.3×10(-4) ≤ r(asp) ≤ 2.5×10(-3)) at a constant φ(p) = 0.01 and (b) varying φ(p) (0.01 ≤ φ(p) ≤ 0.025) at a constant r(asp) = 2.5×10(-3), a mechanism of which is proposed in the text, where r(asp) ≡ d/ ̅L, with d and ̅L being thickness and the average lateral size of the plates, respectively; (2) the r(asp)-induced phase transition of the LC phase from a nematic phase to a highly periodic layered phase, the line shapes of the scattering peaks of which were examined by Caillé's analysis, upon increasing r(asp) under the condition (a); (3) the LC phase having remarkable concentration fluctuations of the particles which are totally unexpected for the conventional lyotropic molecular LC but which are anticipated to be general for the platelike colloidal particles.

Entropy-driven enhanced self-diffusion in confined reentrant supernematics

We present a molecular dynamics study of reentrant nematic phases using the Gay-Berne-Kihara model of a liquid crystal in nanoconfinement. At densities above those characteristic of smectic A phases, reentrant nematic phases form that are characterized by a large value of the nematic order parameter S≃1. Along the nematic director these "supernematic" phases exhibit a remarkably high self-diffusivity, which exceeds that for ordinary, lower-density nematic phases by an order of magnitude. Enhancement of self-diffusivity is attributed to a decrease of rotational configurational entropy in confinement. Recent developments in the pulsed field gradient NMR technique are shown to provide favorable conditions for an experimental confirmation of our simulations.

The role played by self-orientational properties in nematics of colloids with molecules axially symmetric

The self-orientational structure factor as well as the short-time self-orientational diffusion coefficient is computed for colloids composed by nonspherical molecules. To compute the short-time dynamics the hydrodynamic interactions are not taken into account. The hard molecules with at least one symmetry axis considered are: rods, spherocylinders, and tetragonal parallelepipeds. Because both orientational properties in study are written in terms of the second and fourth order parameters, these automatically hold the features of the order parameters. That is, they present a discontinuity for first order transitions, determining in this way the spinodal line. In order to analyze the nematic phase only, we choose the appropriate values for the representative quantities that characterize the molecules. Different formalisms are used to compute the structural properties: de Gennes-Landau approach, Smoluchowski equation and computer simulations. Some of the necessary inputs are taken from literature. Our results show that the self-orientational properties play an important role in the characterization and the localization of axially symmetric phases. While the self-structure decreases throughout the nematics, the short-time self-diffusion does not decrease but rather increases. We study the evolution of the second and fourth order parameters; we find different responses for axial and biaxial nematics, predicting the possibility of a biaxial nematics in tetragonal parallelepiped molecules. By considering the second order in the axial-biaxial phase transition, with the support of the self-orientational structure factor, we are able to propose the density at which this occurs. The short-time dynamics is able to predict a different value in the axial and the biaxial phases. Because the different behavior of the fourth order parameter, the diffusion coefficient is lower for a biaxial phase than for an axial one. Therefore the self-structure factor is able to localize continuous phase transitions involving axially symmetric phases and the short-time self-orientational diffusion is able to distinguish the ordered phase by considering the degree of alignment, that is, axial or biaxial.

Mechanism of diffusion in the smectic- A phase of wormlike rods studied by computer simulation

Molecular dynamics computer simulations have been carried out in the smectic- A phase of stiff wormlike rods. The analysis of the long trajectories generated has allowed for a detailed insight into that diffusion mechanism which is operative in the above-mentioned liquid-crystalline phase, as recently visualized in a system of colloidal virus rods. Fast particles, i.e., those able to move abruptly out from one into an adjacent layer, have been identified. Their properties, such as the velocity autocorrelation function and the orientational distribution function, have been determined and compared to the corresponding quantities valid for a generic rod.

Interfacial anisotropy in the transport of liquid crystals confined between flat, structureless walls: a molecular dynamics simulation approach

Molecular dynamics simulations of uniaxial Gay-Berne ellipsoids as prolate liquid crystal molecules confined between two flat, structureless walls have been carried out in order to investigate anisotropy in their dynamic properties. Several physical quantities are profiled as a function of distance from a wall. The walls stimulate ellipsoids into different behaviors from those of the bulk system. The profiles of self-diffusion coefficients, which are distinguished in each direction of a director-based coordinate system, show that the ellipsoids are more diffusive parallel to the walls and less diffusive perpendicular to the walls with decreasing distance from the walls. According to the self-rotation coefficient and rotational viscosity profiles, ellipsoids are easy to rotate parallel to the walls and hard to rotate in the plane perpendicular to the walls. The analyses of velocity autocorrelation functions, angular velocity autocorrelation functions, director angular velocity autocorrelation functions, and their spectra are useful for the investigation of anisotropy near the walls. We conclude that the flat, structureless wall not only prevents ellipsoids from diffusing and rotating in the plane perpendicular to the walls, but also stimulates them to diffuse and rotate in the plane parallel to the walls.

Orientational relaxation in a discotic liquid crystal

We investigate orientational relaxation of a model discotic liquid crystal, consisting of disclike molecules, by molecular dynamics simulations along two isobars starting from the high temperature isotropic phase. The two isobars have been so chosen that (a) the phase sequence isotropic- (I-) nematic- (N-) columnar (C) appears upon cooling along one of them and (b) the sequence isotropic- (I-) columnar- (C) along the other. While the orientational relaxation in the isotropic phase near the I-N phase transition in system (a) shows a power law decay at short to intermediate times, such power law relaxation is not observed in the isotropic phase near the I-C phase boundary in system (b). In order to understand this difference (the existence or the absence of the power law decay), we calculated the growth of the orientational pair distribution functions (OPDFs) near the I-N phase boundary and also near the I-C phase boundary. We find that the OPDF shows a marked growth in long range correlation as the I-N phase boundary is approached in the I-N-C system (a), but such a growth is absent in the I-C system, which appears to be consistent with the result that I-N phase transition in the former is weakly first order while the I-C phase transition in the latter is not weak. As the system settles into the nematic phase, the decay of the single-particle second-rank orientational time correlation function follows a pattern that is similar to what is observed with calamitic liquid crystals and supercooled molecular liquids.

Analysis of the transverse and the longitudinal pseudodiffusion of CO2 in sub- and supercritical states: a molecular-dynamics analysis

We have performed molecular-dynamics simulations of CO(2) system along the gas-liquid coexistence curve and on the isochore 94.22 cm(3) mol(-1) (which corresponds to the critical isochore). The calculation has been carried out in order to analyze the diffusion of CO(2) and particularly to figure out how the diffusion coefficient may be decomposed along the molecular axes. This makes it possible to analyze the anisotropy of the diffusion along these axes and to shed light on the microscopic changes which accompany such behavior. This anisotropy is traced back to the effect of the translation-rotation coupling (TRC) along the molecular axes. Along the liquid-gas coexistence curve, the pseudolongitudinal diffusion is found to be more rapid than the transverse one. The opposite trend is found along the isochore 94.22 cm(3) mol(-1). The role of the local structure was explored by calculating intermediate scattering function and the autocorrelation functions for the forces acting along the molecular axes. It is shown that the strength of the TRC effect is correlated to the difference between the relaxation times of the local structure, that of the reorientation along the molecular axes, and that of the translational motion. The analysis of the correlation time and the average mean square force along the longitudinal and transverse directions confirms the anisotropy of the local environment that determines the translational dynamics of a molecule.

Analysis of the Effect of Translation−Rotation Coupling on Diffusion along the Molecular Axes

The analysis of the microscopic diffusion processes of CO(2) and H(2)O in the coordinate system defined by the molecular orientation allows a new criterion to be introduced which provides information on the short and long time behavior of the rotation-translation coupling as well as to quantify the strength of this coupling. We discuss the general conditions under which this affects the translation diffusion "seen" by the molecule along its molecular axes. The results show that the translation-rotation coupling is correlated to the local environment in shaping the longitudinal and transversal translation dynamics of a molecule at a microscopic level.

Diffusion and viscosity of a calamitic liquid crystal model studied by computer simulation

We report a molecular dynamics simulation study on an ensemble of rod-like particles, each composed of nine soft spheres held rigidly along a line. We have calculated translational mean square displacements and velocity autocorrelation functions in the fluid phases exhibited by the model, i.e., smectic A, nematic and isotropic. These quantities have then been used to compute diffusion coefficients. In addition, we have calculated viscosities in the nematic and isotropic phases. Despite its crude nature, the model is capable of providing a faithful reproduction of many features of the transport behavior observed in real liquid-crystalline materials. The simulation results have been compared with the predictions of the modified affine transformation theory, finding only a fair agreement.

Anisotropic diffusion in nematic liquid crystals and in ferrofluids

A unified, mean-field kinetic theory approach to the anisotropic translational diffusion observed in liquid crystals and in ferrofluids is proposed. In the dilute regime, unified expressions for the parallel as well as for the perpendicular diffusion coefficient in terms of orientational order parameters are found that apply for liquid crystals as well as for ferrofluids. This result explains the common origin of the anisotropic diffusion found in liquid crystals and in ferrofluids. Differences between the two liquids appear in the semi-dilute regime, where the diffusion coefficients depend on the specific interaction potentials. Explicit expressions for the diffusion coefficients are worked out also in this regime within a mean-field approximation. Comparisons with previous theoretical and experimental results are performed, showing satisfactory agreement to the present results.

The Enskog-type theory of the velocity autocorrelations in the two-dimensional nematic of hard needles

As it was shown from molecular dynamics of two-dimensional hard needles, the uniaxial velocity autocorrelation function (VACF) of this system exhibits a two time scale character. This corresponds to the symmetry of the particles. In this paper we provide a theory of the Enskog type that corroborates the idea that the VACF can be successfully described as a sum of two single decays. From the comparison between molecular dynamics and theoretical results, we show that the Enskog approach is a satisfactory kinetic theory, which functions as well in the nematic as in the isotropic phase. Different properties of VACFs have been investigated, in particular, the dependence on the orientational order, temperature, and particles' inertness.

Influence of shape and energy anisotropies on the phase diagram of discotic molecules

We present Monte Carlo simulations of discotic molecules using the Gay-Berne potential with shape (kappa) and energy (kappa(')) anisotropies. Following the previous work of Bates and Luckhurst [J. Chem. Phys. 104, 6696 (1996)] at kappa=0.345, kappa(')=0.2 when we determine the sequence of different phases at the same reduced pressure P(*)=50, we find an additional phase at low temperatures corresponding to an orthorhombic crystalline phase and we characterize it. Keeping the shape anisotropy fixed at kappa=0.2, we determine the evolution of the phase diagram with varying energy anisotropy. At high kappa('), low anisotropy, the system is not able to build columns while at low kappa('), the system exhibits both orthorhombic crystal as well as hexagonal liquid crystal phases over a wide range of pressures and temperatures. The domain of stability of the nematic phase is found to systematically shift towards higher pressures as kappa(') decreases.

Anisotropic self-diffusion in thermotropic liquid crystals studied by 1H and 2H pulse-field-gradient spin-echo NMR

The molecular self-diffusion coefficients in nematic and smectic-A thermotropic liquid crystals are measured using stimulated-echo-type 2H and 1H pulse-field-gradient spin-echo nuclear magnetic resonance (PGSE NMR) combined with multiple-pulse dipolar decoupling and slice selection. The temperature dependence of the principal components of the diffusion tensor in the nematic phase follows a simple Arrhenius relationship except in the region of nematic-isotropic phase transition where it reflects, merely, the decrease of the molecular orientational order. The average of the principal diffusion coefficients in the isotropic-nematic phase transition region is close to the diffusion coefficient in the isotropic phase. At the nematic-smectic-A phase transition the diffusion coefficients change continuously. The results in nematic phase are best described in terms of the affine transformation model for diffusion in nematics formed by hard ellipsoids. In the smectic-A phase the data are interpreted using a modified model for diffusion in presence of a periodic potential along the director.

Transport tensors in perfectly aligned low-density fluids: self-diffusion and thermal conductivity

The modified Taxman equation for the kinetic theory of low-density fluids composed of rigid aspherical molecules possessing internal degrees of freedom is generalized to obtain the transport tensors in a fluid of aligned molecules. The theory takes care of the shape of the particles exactly but the solution has been obtained only for the case of perfectly aligned hard spheroids within the framework of the first Sonine polynomial approximation. The expressions for the thermal-conductivity components have been obtained for the first time whereas the self-diffusion components obtained here turn out to be exactly the same as those derived by Kumar and Masters [Mol. Phys. 81, 491 (1994)] through the solution of the Lorentz-Boltzmann equation. All our expressions yield correct results in the hard-sphere limit.

Simulations using hard particles

A short review is given of recent progress in the computer simulation of liquid crystal phases using hard particles. Emphasis is placed on the richness of phase behaviour that may result from the effects of molecular size and shape alone, and on the role of simulations in testing modern theories of liquid crystal phase transitions, structure and dynamics. Two specific examples are treated in detail: the simulation of twisted nematic liquid crystals, allowing a direct calculation of the twist elastic constant and the helical twisting power of chiral dopant molecules; and the recent quantitative explanation of diffusive behaviour in isotropic and nematic liquids using kinetic theory.