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.

Polymorphs solvates and co-crystals of molecular materials

The scope of my research project is to produce and characterize new crystalline forms of organic compounds, focusing the attention on co-crystals and then transferring these notions on APIs to produce co-crystals of potential interest in the pharmaceutical field. In the first part of this work co-crystallization experiments were performed using as building blocks the family of aliphatic dicarboxylic acids HOOC-(CH2)n-COOH, with n= 2-8. This class of compounds has always been an object of study because it is characterized by an interesting phenomenon of alternation of melting points: the acids with an even number of carbon atoms show a melting point higher than those with an odd one. The acids were co-crystallized with four dipyridyl molecules (formed by two pyridine rings with a different number of bridging carbon atoms) through the formation of intermolecular interactions N•••(H)O. The bases used were: 4,4’-bipyridine (BPY), 1,2-bis(4-pyridyl)ethane (BPA), 1,2-(di-4-pyridyl)ethylene (BPE) and 1,2-bis(4-pyridyl)propane (BPP). The co-crystals obtained by solution synthesis were characterized by different solid-state techniques to determine the structure and to see how the melting points in co-crystals change. In the second part of this study we tried to obtain new crystal forms of compounds of pharmaceutical interest. The APIs studied are: O-desmethylvenlafaxine, Lidocaine, Nalidixic Acid and Sulfadiazine. Each API was subjected to Polymorph Screening and Salt/Co-crystal Screening experiments to identify new crystal forms characterized by different properties. In a typical Salt/Co-crystal Screening the sample was made to react with a co-former (solid or liquid) through different methods: crystallization by solution, grinding, kneading and solid-gas reactions. The new crystal forms obtained were characterized by different solid state techniques (X-ray single crystal diffraction, X-ray powder diffraction, Differential Scanning Calorimetry, Thermogravimetric Analysis, Evolved gas analysis, FT-IR – ATR, Solid State N.M.R).