SELF-DIFFUSION AND HEAT-FLOW IN ISOTROPIC AND LIQUID-CRYSTAL PHASES OF THE GAY-BERNE FLUID

We use the so-called Evans heat flow algorithm and the color conductivity algorithm to calculate the thermal conductivity, lambda, and self-diffusion, D, tensors of a model liquid crystal forming fluid, the Gay-Berne fluid. We compared the conductivities of the isotropic phase, the nematic liquid crystal phase, and the smectic B phase and found that the transverse components of lambda and D in the nematic phase are about the same as lambda and D for the isotropic phase at the same temperature and slightly lower density. The parallel components of A and D are about two and four times as large as the respective transverse components. In the smectic B phase the parallel thermal conductivity ratio increases to about three, whereas the diffusion coefficient is virtually zero because of the solid like nature of this phase. We cross checked the results by deriving and evaluating the Green-Kubo relations for lambda and D using conventional equilibrium molecular dynamics. We also discovered that in the presence of a temperature gradient, the director of the nematic phase, spontaneously aligns itself into the plane perpendicular to the temperature gradient. This effect which is nonlinear in the temperature gradient, minimizes the entropy production in the heat conducting steady state.