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.

Transport Properties and Novel Sensing Applications of Organic Semiconducting Crystals

The present thesis is focused on the study of Organic Semiconducting Single Crystals (OSSCs) and crystalline thin films. In particular solution-grown OSSC, e.g. 4-hdroxycyanobenzene (4HCB) have been characterized in view of their applications as novel sensors of X-rays, gamma-rays, alpha particles radiations and chemical sensors. In the field of ionizing radiation detection, organic semiconductors have been proposed so far mainly as indirect detectors, i.e. as scintillators or as photodiodes. I first study the performance of 4HCB single crystals as direct X-ray detector i.e. the direct photon conversion into an electrical signal, assessing that they can operate at room temperature and in atmosphere, showing a stable and linear response with increasing dose rate. A dedicated study of the collecting electrodes geometry, crystal thickness and interaction volume allowed us to maximize the charge collection efficiency and sensitivity, thus assessing how OSSCs perform at low operating voltages and offer a great potential in the development of novel ionizing radiation sensors. To better understand the processes generating the observed X-ray signal, a comparative study is presented on OSSCs based on several small-molecules: 1,5-dinitronaphthalene (DNN), 1,8-naphthaleneimide (NTI), Rubrene and TIPS-pentacene. In addition, the proof of principle of gamma-rays and alpha particles has been assessed for 4HCB single crystals. I have also carried out an investigation of the electrical response of OSSCs exposed to vapour of volatile molecules, polar and non-polar. The last chapter deals with rubrene, the highest performing molecular crystals for electronic applications. We present an investigation on high quality, millimeter-sized, crystalline thin films (10 – 100 nm thick) realized by exploiting organic molecular beam epitaxy on water-soluble substrates. Space-Charge-Limited Current (SCLC) and photocurrent spectroscopy measurements have been carried out. A thin film transistor was fabricated onto a Cytop® dielectric layer. The FET mobility exceeding 2 cm2/Vs, definitely assess the quality of RUB films.

Computing the structural and vibrational properties of polymorphic organic molecular crystals through van der Waals corrected density functional theory and the electronic properties of organic thin films through microelectrostatic calculations.

The present Thesis reports on the various research projects to which I have contributed during my PhD period, working with several research groups, and whose results have been communicated in a number of scientific publications. The main focus of my research activity was to learn, test, exploit and extend the recently developed vdW-DFT (van der Waals corrected Density Functional Theory) methods for computing the structural, vibrational and electronic properties of ordered molecular crystals from first principles. A secondary, and more recent, research activity has been the analysis with microelectrostatic methods of Molecular Dynamics (MD) simulations of disordered molecular systems. While only very unreliable methods based on empirical models were practically usable until a few years ago, accurate calculations of the crystal energy are now possible, thanks to very fast modern computers and to the excellent performance of the best vdW-DFT methods. Accurate energies are particularly important for describing organic molecular solids, since they often exhibit several alternative crystal structures (polymorphs), with very different packing arrangements but very small energy differences. Standard DFT methods do not describe the long-range electron correlations which give rise to the vdW interactions. Although weak, these interactions are extremely sensitive to the packing arrangement, and neglecting them used to be a problem. The calculations of reliable crystal structures and vibrational frequencies has been made possible only recently, thanks to development of some good representations of the vdW contribution to the energy (known as “vdW corrections”).