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.

Mixed Mode Fracture Behaviour of Piezoelectric Materials

Piezoelectrics present an interactive electromechanical behaviour that, especially in recent years, has generated much interest since it renders these materials adapt for use in a variety of electronic and industrial applications like sensors, actuators, transducers, smart structures. Both mechanical and electric loads are generally applied on these devices and can cause high concentrations of stress, particularly in proximity of defects or inhomogeneities, such as flaws, cavities or included particles. A thorough understanding of their fracture behaviour is crucial in order to improve their performances and avoid unexpected failures. Therefore, a considerable number of research works have addressed this topic in the last decades. Most of the theoretical studies on this subject find their analytical background in the complex variable formulation of plane anisotropic elasticity. This theoretical approach bases its main origins in the pioneering works of Muskelishvili and Lekhnitskii who obtained the solution of the elastic problem in terms of independent analytic functions of complex variables. In the present work, the expressions of stresses and elastic and electric displacements are obtained as functions of complex potentials through an analytical formulation which is the application to the piezoelectric static case of an approach introduced for orthotropic materials to solve elastodynamics problems. This method can be considered an alternative to other formalisms currently used, like the Stroh’s formalism. The equilibrium equations are reduced to a first order system involving a six-dimensional vector field. After that, a similarity transformation is induced to reach three independent Cauchy-Riemann systems, so justifying the introduction of the complex variable notation. Closed form expressions of near tip stress and displacement fields are therefore obtained. In the theoretical study of cracked piezoelectric bodies, the issue of assigning consistent electric boundary conditions on the crack faces is of central importance and has been addressed by many researchers. Three different boundary conditions are commonly accepted in literature: the permeable, the impermeable and the semipermeable (“exact”) crack model. This thesis takes into considerations all the three models, comparing the results obtained and analysing the effects of the boundary condition choice on the solution. The influence of load biaxiality and of the application of a remote electric field has been studied, pointing out that both can affect to a various extent the stress fields and the angle of initial crack extension, especially when non-singular terms are retained in the expressions of the electro-elastic solution. Furthermore, two different fracture criteria are applied to the piezoelectric case, and their outcomes are compared and discussed. The work is organized as follows: Chapter 1 briefly introduces the fundamental concepts of Fracture Mechanics. Chapter 2 describes plane elasticity formalisms for an anisotropic continuum (Eshelby-Read-Shockley and Stroh) and introduces for the simplified orthotropic case the alternative formalism we want to propose. Chapter 3 outlines the Linear Theory of Piezoelectricity, its basic relations and electro-elastic equations. Chapter 4 introduces the proposed method for obtaining the expressions of stresses and elastic and electric displacements, given as functions of complex potentials. The solution is obtained in close form and non-singular terms are retained as well. Chapter 5 presents several numerical applications aimed at estimating the effect of load biaxiality, electric field, considered permittivity of the crack. Through the application of fracture criteria the influence of the above listed conditions on the response of the system and in particular on the direction of crack branching is thoroughly discussed.

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).

Speleogenesis and secondary cave minerals in quartz-sandstone and quartzite environment

The main objective of this research is to improve the comprehension of the processes controlling the formation of caves and karst-like morphologies in quartz-rich lithologies (more than 90% quartz), like quartz-sandstones and metamorphic quartzites. In the scientific community the processes actually most retained to be responsible of these formations are explained in the “Arenisation Theory”. This implies a slow but pervasive dissolution of the quartz grain/mineral boundaries increasing the general porosity until the rock becomes incohesive and can be easily eroded by running waters. The loose sands produced by the weathering processes are then evacuated to the surface through processes of piping due to the infiltration of waters from the fracture network or the bedding planes. To deal with these problems we adopted a multidisciplinary approach through the exploration and the study of several cave systems in different tepuis. The first step was to build a theoretical model of the arenisation process, considering the most recent knowledge about the dissolution kinetics of quartz, the intergranular/grain boundaries diffusion processes, the primary diffusion porosity, in the simplified conditions of an open fracture crossed by a continuous flow of undersatured water. The results of the model were then compared with the world’s widest dataset (more than 150 analyses) of water geochemistry collected till now on the tepui, in superficial and cave settings. All these studies allowed verifying the importance and the effectiveness of the arenisation process that is confirmed to be the main process responsible of the primary formation of these caves and of the karst-like superficial morphologies. The numerical modelling and the field observations allowed evaluating a possible age of the cave systems around 20-30 million of years.

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.

Structure, Dynamics and Reactivity in the Organic Solid State: Anthracene Derivatives and Charge Transfer Crystals

The work presented in this thesis tackles some important points concerning the collective properties of two typical categories of molecular crystals, i.e., anthracene derivatives and charge transfer crystals. Anthracene derivatives have constituted the class of materials from which systematical investigations of crystal-to-crystal photodimerization reactions started, developed and have been the subject of a new awakening in the recent years. In this work some of these compounds, namely, 9-cyanoanthacene, 9-anthacenecarboxylic acid and 9-methylanthracene, have been selected as model systems for a phenomenological approach to some key properties of the solid state, investigated by spectroscopic methods. The present results show that, on the basis of the solid state organization and the chemical nature of each compound, photo-reaction dynamics and kinetics display distinctive behaviors, which allows for a classification of the various processes in topochemical, non topochemical, reversible or topophysical. The second part of the thesis was focused on charge transfer crystals, binary systems formed by stoichiometric combinations of the charge donating perylene (D) and the charge accepting tetracyano-quinodimethane (A), this latter also in its fluorinated derivatives. The work was focused on the growth of single crystals, some of which not yet reported in the literature, by PVT technique. Structural and spectroscopic characterizations have been performed, with the aim of determining the degree of charge transfer between donor and acceptor in the co-crystals. An interesting outcome of the systematic search performed in this work is the definition of the experimental conditions which drive the crystal growth of the binary systems either towards the low (1:1) or the high ratio (3:1 or 3:2) stoichiometries.

Atomistic simulations of liquid crystals in the bulk and at their interfaces

In this thesis, we have dealt with several problems concerning liquid crystals (LC) phases, either in the bulk or at their interfaces, by the use of atomistic molecular dynamics (MD) simulations. We first focused our attention on simulating and characterizing the bulk smectic phase of 4-n-octyl-4'-cyanobiphenyl (8CB), allowing us to investigate the antiparallel molecular arrangement typical of SmAd smectic phases. A second topic of study was the characterization of the 8CB interface with vacuum by simulating freely suspended thin films, which allowed us to determine the influence of the interface on the orientational and positional order. Then we investigated the LC-water and LC-electrolyte water solution interface. This interface has recently found application in the development of sensors for several compounds, including biological molecules, and here we tried to understand the re-orientation mechanism of LC molecules at the interface which is behind the functioning of these sensors. The characterization of this peculiar interface has incidentally led us to develop a polarizable force field for the pentyl-cyanobiphenyl mesogen, whose process of parametrization and validation is reported here in detail. We have shown that this force field is a significant improvement over its previous, static charge non polarizable version in terms of density, orientational order parameter and translational diffusion.

Optimization of molecular and crystalline forms of drugs, agrochemicals, pesticides in relation to activity, bioavailability, patentability and to the fabrication of polymorphs, solvates, co-crystals with green chemistry methods

This doctorate was funded by the Regione Emilia Romagna, within a Spinner PhD project coordinated by the University of Parma, and involving the universities of Bologna, Ferrara and Modena. The aim of the project was: - Production of polymorphs, solvates, hydrates and co-crystals of active pharmaceutical ingredients (APIs) and agrochemicals with green chemistry methods; - Optimization of molecular and crystalline forms of APIs and pesticides in relation to activity, bioavailability and patentability. In the last decades, a growing interest in the solid-state properties of drugs in addition to their solution chemistry has blossomed. The achievement of the desired and/or the more stable polymorph during the production process can be a challenge for the industry. The study of crystalline forms could be a valuable step to produce new polymorphs and/or co-crystals with better physical-chemical properties such as solubility, permeability, thermal stability, habit, bulk density, compressibility, friability, hygroscopicity and dissolution rate in order to have potential industrial applications. Selected APIs (active pharmaceutical ingredients) were studied and their relationship between crystal structure and properties investigated, both in the solid state and in solution. Polymorph screening and synthesis of solvates and molecular/ionic co-crystals were performed according to green chemistry principles. Part of this project was developed in collaboration with chemical/pharmaceutical companies such as BASF (Germany) and UCB (Belgium). We focused on on the optimization of conditions and parameters of crystallization processes (additives, concentration, temperature), and on the synthesis and characterization of ionic co-crystals. Moreover, during a four-months research period in the laboratories of Professor Nair Rodriguez-Hormedo (University of Michigan), the stability in aqueous solution at the equilibrium of ionic co-crystals (ICCs) of the API piracetam was investigated, to understand the relationship between their solid-state and solution properties, in view of future design of new crystalline drugs with predefined solid and solution properties.

New Application of Piezoelectric Ultrasounds in Maxillo-facial bone surgery

Background: Piezoelectric instrumentation seems to offer 3 important advantages for cutting bone structures. Be more precise because it is produced by micro-vibrations from the cutting insert. Be safer because the ultrasonic frequency used does not affect soft tissue. Thirdly, the less invasive cutting action produces minor tissue damage and consequently probably a better healing Aim of the Study: The aim of this study is to evaluate the effectiveness of piezoelectric device capability in maxillo-facial surgery, in order to take advantage of these favourable capacity. Material and Methods: Considering the several potential application of the piezoelectric technology in Orthognathic, Oncologic and Extractive surgery, we would like to design protocols in order to verify how this new device can modify the surgical technique, the surgical time, the patients healing and its quality of life. Results: Due to the precise Piezosurgery cut, we can manage the Cad-Cam-Custom Made plates protocol in Oncologic Surgery and in Orthognatic Surgery increasing our percentage of comparison between the 3D preoperative plan and the surgical execution. We also found a better quality of life impaction in Patient who underwent and extractive surgery Conclusion: Piezosurgery device seems to be a strong surgical aid were safe and precise cut are needed and its capability to reduce the discomfort Patients need to be study in deep also in major surgery like Orthognatic and Oncologic surgery.

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”).

Nanostructured piezoelectric materials for the design and development of self-sensing composite materials and energy harvesting devices

The work activities reported in this PhD thesis regard the functionalization of composite materials and the realization of energy harvesting devices by using nanostructured piezoelectric materials, which can be integrated in the composite without affecting its mechanical properties. The self-sensing composite materials were fabricated by interleaving between the plies of the laminate the piezoelectric elements. The problem of negatively impacting on the mechanical properties of the hosting structure was addressed by shaping the piezoelectric materials in appropriate ways. In the case of polymeric piezoelectric materials, the electrospinning technique allowed to produce highly-porous nanofibrous membranes which can be immerged in the hosting matrix without inducing delamination risk. The flexibility of the polymers was exploited also for the production of flexible tactile sensors. The sensing performances of the specimens were evaluated also in terms of lifetime with fatigue tests. In the case of ceramic piezo-materials, the production and the interleaving of nanometric piezoelectric powder limitedly affected the impact resistance of the laminate, which showed enhanced sensing properties. In addition to this, a model was proposed to predict the piezoelectric response of the self-sensing composite materials as function of the amount of the piezo-phase within the laminate and to adapt its sensing functionalities also for quasi-static loads. Indeed, one final application of the work was to integrate the piezoelectric nanofibers in the sole of a prosthetic foot in order to detect the walking cycle, which has a period in the order of 1 second. In the end, the energy harvesting capabilities of the piezoelectric materials were investigated, with the aim to design wearable devices able to collect energy from the environment and from the body movements. The research activities focused both on the power transfer capability to an external load and the charging of an energy storage unit, like, e.g., a supercapacitor.

Crystal engineering of co-crystals and host-guest systems for environmentally friendly applications: sunscreens stabilization, pollutant sequestering and herbicide reformulation.

Solid state engineered materials have proven to be useful and suitable tools in the quest of new materials. In this thesis different crystalline compounds were synthesized to provide more sustainable products for different applications, as in cosmetics or in agrochemistry, to propose pollutants removal strategy or to obtain materials for electrocatalysis. Therefore, the research projects presented here can be divided into three main topics: (i) sustainable preparation of solid materials of widely used active ingredients aimed at the reduction of their occurrence in the natural environment. The systems studied in this section are cyclodextrins host-guest compounds, obtained via mechanochemical and slurry synthesis. The first chemicals studied are sunscreens inclusion complexes, that proved to have enhanced photostability and desired photoprotection. The same synthetic methods were applied to obtain inclusion complexes of bentazon, a herbicide often found to leach in groundwaters. The resulting products showed to have desired water solubility properties. The same herbicide was also adsorbed on amorphous calcium phosphate nanoparticles, to obtain a biocompatible formulation of this agrochemical. This herbicide could benefit by the adsorption on nanoparticles for what concerns its kinetic release in different media as well as its photostability. (ii) Sustainable synthesis of co-crystals based on polycyclic aromatic hydrocarbons, for the proposal of a sequestering method with a resulting material with enhanced properties. The co-crystallization via mechanochemical means proved that these pollutants can be sequestered via simple solvent-free synthesis and the obtained materials present better photochemical properties when compared to the starting co-formers. (iii) Crystallization from mild solvents of nanosized materials useful for the application in electrocatalysis. The study of compounds based on nickel and cobalt metal ions resulted in the obtainment of 2D and 1D coordination polymers. Moreover, solid solutions were obtained. These crystals showed layered structures and, according to preliminary results, they can be exfoliated.

Design and characterization of polymorphs and co-crystals of organic molecular semiconductors

In the last decades, organic semiconductors have attracted attention due to their possible employment in solution-processed optoelectronic and electronic devices. One of the advantages of solution processing is the possibility to process into flexible substrates at low cost. Organic molecular materials tend to form polymorphs, which can exhibit very different properties. In most cases, the control of the crystal structure is decisive to maximize the performance of the final device. Although organic electronics have progressed a lot, n-type organic semiconductors still lag behind p-type, presenting challenges such as air instability and poor solubility. NDI derivatives are promising candidates for applications in organic electronics due to their characteristics. Recently, the structure-properties relationship and the polymorphism of these molecules have gained attention. In the first part of this thesis, NDI-C6 thermal behavior was extensively explored which revealed two different behaviors depending on the annealing process. This study allowed to define the stability ranking of the NDI-C6 bulk forms and to determine the crystal structure of Form γ at 54°C. Additionally, the polymorphic and thermal behavior of thin films of NDI-C6 was also explored. It was possible to isolate pure Form α, Form β, Form γ and a new metastable Form ε. It was also possible to determine the stability ranking of the phases in thin films. OFETs were fabricated having different polymorphs as active layer, unfortunately the performance was not ideal. During the second part of this thesis, core-chlorinated NDIs with fluoroalkyl chains were studied. Initially, the focus was on the polymorphism of CF3-NDI that revealed a solvate form with a very interesting molecular arrangement suggesting the possibility to form charge transfer co-crystals. In the last part of the thesis, the synthesis and characterization of CT co-crystal with different NDI derivatives, and acceptor and as donor BTBT and ditBu-BTBT were explored.