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.

Toughened ZrO2 ceramics sintered with a La2O3-rich glass as additive

In this study, the influence of the glass addition and sintering parameters on the densification and mechanical properties of tetragonal zirconia polycrystals (3Y-TZP) ceramics were evaluated. High-purity tetragonal ZrO2 powder and La2O3-rich glass were used as starting powders. Two compositions based on ZrO2 and containing 5wt.% and 10wt.% of La2O3-rich glass were studied in this work. The starting powders were mixed/milled by planetary milling, dried at 90 degrees C for 24 h, sieved through a 60 mesh screen and uniaxially cold pressed under 80 MPa. The samples were sintered in air at 1200 degrees C, 1300 degrees C, 1400 degrees C for 60 min and at 1450 degrees C for 120 min, with heating and cooling rates of 10 degrees C/min. Sintered samples were characterized by relative density, X-ray diffraction (XRD) and scanningelectron microscopy (SEM). Hardness and fracture toughness were obtained by Vickers indentation method. Dense sintered samples were obtained for all conditions. Furthermore, only tetragonal-ZrO2 was identified as crystalline phase in sintered samples, independently of the conditions studied. Samples sintered at 1300 degrees C for 60 min presented the optimal mechanical properties with hardness and fracture toughness values near to 12 GPa and 8.5 MPa m(1/2) respectively. (c) 2007 Elsevier B.V, All rights reserved.

Protein partitioning in poly(ethylene glycol)/sodium polyacrylate aqueous two-phase systems

The partition of hemoglobin, lysozyme and glucose-6-phospate dehydrogenase (G6PDH) in a novel inexpensive aqueous two-phase system (ATPS) composed by poly(ethylene glycol) (PEG) and sodium polyacrylate (NaPA) has been studied. The effect of NaCl and Na2SO4, pH and PEG molecular size on the partitioning has been studied. At high pH (above 9), hemoglobin partitions strongly to the PEG-phase. Although some precipitation of hemoglobin occurs, high recovery values are obtained particularly for lysozyme and G6PDH. The partitioning forces are dominated by the hydrophobic and electrochemical (salt) effects, since the positively charged lysozyme and negatively charged G6PDH partitions to the non-charged PEG and the strongly negatively charged polyacrylate enriched phase, respectively. (c) 2007 Elsevier B.V. All rights reserved.

In situ gelling hexagonal phases for sustained release of an anti-addiction drug

In this study, fluid precursor formulations for subcutaneous injection and in situ formation of hexagonal phase gels upon water absorption were developed as a strategy to sustain the release of naltrexone, a drug used for treatment of drug addiction. Precursor formulations were obtained by combining BRIJ 97 with propylene glycol (PG, 5-70%, w/w). To study the phase behavior of these formulations, water was added at 10-90% (w/w), and the resulting systems were characterized by polarized light microscopy. Two precursor formulations containing BRIJ:PG at 95:5 (w/w, referred to as BRIJ-95) and at 80:20 (w/w, referred to as BRIJ-80) were chosen. Naltrexone was dissolved at 1% or suspended at 5% (w/w). Precursor formulations were transformed into hexagonal phases when water content exceeded 20%. Water uptake followed second-order kinetics, and after 2-4 h all precursor formulations were transformed into hexagonal phases. Drug release was prolonged by the precursor formulations (compared to a drug solution in PBS), and followed pseudo-first order kinetics regardless of naltrexone concentration. The release from BRIJ-80 was significantly higher than that from BRIJ-95 after 48 h. The relative safety of the precursor formulations was assessed in cultured fibroblasts. Even though BRIJ-95 was more cytotoxic than BRIJ-80, both precursor formulations were significantly less cytotoxic than sodium lauryl sulfate (considered moderate-to-severe irritant) at the same concentration (up to 50 mu g/mL). These results suggest the potential of BRIJ-based precursor formulations for sustained naltrexone release. (C) 2011 Elsevier By. All rights reserved.

Solid-state NMR structure determination of melittin in a lipid environment

Solid-state C-13 NMR spectroscopy was used to investigate the three-dimensional structure of melittin as lyophilized powder and in ditetradecylphosphatidylcholine (DTPC) membranes. The distance between specifically labeled carbons in analogs [1-C-13]Gly3-[2-C-13]Ala4, [1-C-13]Gly3-[2-C-13]Leu6, [1-C-13]Leu13-[2-C-13]Ala15, [2-C-13]Leu13-[1-C-13]Ala15, and [1-C-13]Leu13-[2-C-13]Leu16 was measured by rotational resonance. As expected, the internuclear distances measured in [1-C-13]Gly3-[2-C-13]Ala4 and [1-C-13]Gly3-[2-C-13]Leu6 were consistent with alpha -helical structure in the N-terminus irrespective of environment. The Internuclear distances measured in [1-C-13]Leu13-[2-C-13]Ala15, [2-C-13]Leu13-[1-C-13]Ala15, and [1-C-13]Leu13-[2-C-13]Leu16 revealed, via molecular modeling, some dependence upon environment for conformation in the region of the bend in helical structure induced by Pro14. A slightly larger interhelical angle between the N- and C-terminal helices was indicated for peptide in dry or hydrated gel state DTPC (139 degrees -145 degrees) than in lyophilized powder (121 degrees -139 degrees) or crystals (129 degrees). The angle, however, is not as great as deduced for melittin in aligned bilayers of DTPC in the liquid-crystalline state (similar to 160 degrees) (R. Smith, F. Separovic, T. J. Milne, A. Whittaker, F. M. Bennett, B. A. Cornell, and A. Makriyannis, 1994, J. Mol, Biol 241:456-466). The study illustrates the utility of rotational resonance in determining local structure within peptide-lipid complexes.

Experimental investigation of liquid flow shift due to gas cross flow in non-wetted packed beds

An experimental study has been carried out for the gas-liquid two-phase flow in a packed bed simulating conditions of the gas and liquid flows in the lower part of blast furnace. The localised liquid flow phenomenon in presence of gas cross flow, which usually occurs around the cohesive zone and raceway in blast furnace, was investigated in detail. Such liquid flow is characterised in terms of liquid shift distance or liquid shift angle that can effectively be measured by the experiments involved in the current study. It is found that liquid shift angle does not significantly increase or decrease with different packing depth. This finding supports the hypothesis of the force balance model where a vectorial relationship among acting forces, i.e. gas drag force, gravitational force and solid-liquid friction force, and liquid shift angle does exist. Liquid shift angle is inversely proportional to particle size and liquid density, and proportional to square of gas superficial velocity, but is almost independent on liquid flowrate and liquid viscosity. The gas-liquid drag coefficient, an important aspect for quantifying the interaction between gas and liquid flows, was conceptually modified based on the discrete feature of liquid flow through a packed bed and evaluated by the combined theoretical and experimental investigation. Experimental measurements suggest that the gas-liquid drag coefficient is approximately a constant (C-DG(')=5.4+/-1.0) and is independent on liquid properties, gas velocity and packing structure. The result shows a good agreement with previous experimental data and prediction of the existing liquid flow model.

Deuterium NMR Study of Orientational Order in Cellulosic Network Microfibers

Deuterium NMR was used to investigate the orientational order in a composite cellulosic formed by liquid crystalline acetoxypropylcellulose (A PC) and demented nematic 4'-penty1-4-cyanobiphenyl (5CB-4 alpha d(2)) with the per centage of 85% A PC by weight Three forms of the composite including electro spun microfibers, thin film and bulk samples were analyzed The NMR results initially suggest two distinct scenarios, one whet e the 503-alpha d(2), is confined to small droplets with dimensions smaller than the magnetic coherence length and the other where the 503-alpha d(2) molecules arc aligned with the A PC network chains Polarized optical microscopy (POW from thin film samples along with all the NMR results show the presence of 5CB-alpha d(2) droplets in the composite systems with a nematic wetting layer at the APC-5CB-alpha d(2) interface that experiences and order disorder transition driven by the polymer network N-I transition The characterization of the APC network I-N transition shows a pronounced subcritical behavior within a heterogeneity scenario.

Water-Based Cellulose Liquid Crystal System Investigated by Rheo-NMR

Water-based cellulose cholesteric liquid crystalline phases at rest can undergo structural changes induced by shear flow. This reflects on the deuterium spectra recorded when the system is investigated by rheo-nuclear magnetic resonance (rheo-NMR) techniques. In this work, the model system hydroxypropylcellulose (HPC)+water is revisited using rheo-NMR to clarify unsettled points regarding its behavior under shear and in relaxation. The NMR spectra allow the identification of five different stable ordering states, within shear and relaxation, which are well integrated in a mesoscopic picture of the system's structural evolution under shear and relaxation. This picture emerging from the large body of studies available for this system by other experimental techniques, accounts well for the NMR data and is in good agreement with the three distinct regions of steady shear flow recognized for some lyotropic LC polymers. Shear rates in between 0.1 and 1.0 s(-1) where investigated using a Taylor-Couette flow and deuterated water was used as solvent for the deuterium NMR (DNMR) analysis.

A cellulosic liquid crystal pool for cellulose nanocrystals: structure and molecular dynamics at high shear rates

Cellulose and its derivatives, such as hydroxypropylcellulose (HPC) have been studied for a long time but they are still not well understood particularly in liquid crystalline solutions. These systems can be at the origin of networks with properties similar to liquid crystalline (LC) elastomers. The films produced from LC solutions can be manipulated by the action of moisture allowing for instance the development of a soft motor (Geng et al., 2013) driven by humidity. Cellulose nanocrystals (CNC), which combine cellulose properties with the specific characteristics of nanoscale materials, have been mainly studied for their potential as a reinforcing agent. Suspensions of CNC can also self-order originating a liquid-crystalline chiral nematic phases. Considering the liquid crystalline features that both LC-HPC and CNC can acquire, we prepared LC-HPC/CNC solutions with different CNC contents (1,2 and 5 wt.%). The effect of the CNC into the LC-HPC matrix was determined by coupling rheology and NMR spectroscopy - Rheo-NMR a technique tailored to analyse orientational order in sheared systems. (C) 2015 Elsevier Ltd. All rights reserved.

Towards cryogenic liquid –vapor energy storage units for space applications

Dissertation to obtain the Doctoral degree in Physics Engineering

Cellulose micro/nano fibers conformational effects probed by nematic liquid crystal droplets

Probing micro-/nano-sized surface conformations, which are ubiquitous in biological systems, by using liquid crystal droplets, which change their ordering and optical appearance in response to the presence of more than ten times smaller cellulose based micro/nano fibers, might find new uses in a range of biological environments and sensors. Previous studies indicate that electrospun micro/nano cellulosic fibers produced from liquid crystalline solutions could present a twisted form [1]. In this work, we study the structures of nematic liquid crystal droplets threaded by cellulose fibers prepared from liquid crystalline and isotropic solutions as well as droplets pierced by spider-made fibers [2]. Planar anchoring at the fibers and planar and homeotropic at the drop surfaces allowed probing cellulose fibers different helical structures as well as aligned filaments.

40 K Neon liquid energy storage unit

Cryocoolers have been progressively replacing the use of the stored cryogens in cryogenic chains used for detector cooling, thanks to their higher and higher reliability. However, the mechanical vibrations, the electromagnetic interferences and the temperature fluctuations inherent to their functioning could reduce the sensor’s sensitivity. In order to minimize this problem, compact thermal energy storage units (ESU) are studied, devices able to store thermal energy without significant temperature increase. These devices can be used as a temporary cold source making it possible to turn the cryocooler OFF providing a proper environment for the sensor. A heat switch is responsible for the thermal decoupling of the ESU from the cryocooler’s temperature that increases when turned OFF. In this work, several prototypes working around 40 K were designed, built and characterized. They consist in a low temperature cell that contains the liquid neon connected to an expansion volume at room temperature for gas storage during the liquid evaporation phase. To turn this system insensitive to the gravity direction, the liquid is retained in the low temperature cell by capillary effect in a porous material. Thanks to pressure regulation of the liquid neon bath, 900 J were stored at 40K. The higher latent heat of the liquid and the inexistence of triple point transitions at 40 K turn the pressure control during the evaporation a versatile and compact alternative to an ESU working at the triple point transitions. A quite compact second prototype ESU directly connected to the cryocooler cold finger was tested as a temperature stabilizer. This device was able to stabilize the cryocooler temperature ((≈ 40K ±1 K) despite sudden heat bursts corresponding to twice the cooling power of the cryocooler. This thesis describes the construction of these devices as well as the tests performed. It is also shown that the thermal model developed to predict the thermal behaviour of these devices, implemented as a software,describes quite well the experimental results. Solutions to improve these devices are also proposed.

40 K neon liquid energy storage unit

Cryocoolers have been progressively replacing the use of the stored cryogens in cryogenic chains used for detector cooling, thanks to their higher and higher reliability. However, the mechanical vibrations, the electromagnetic interferences and the temperature fluctuations inherent to their functioning could reduce the sensor’s sensitivity. In order to minimize this problem, compact thermal energy storage units (ESU) are studied, devices able to store thermal energy without significant temperature increase. These devices can be used as a temporary cold source making it possible to turn the cryocooler OFF providing a proper environment for the sensor. A heat switch is responsible for the thermal decoupling of the ESU from the cryocooler’s temperature that increases when turned OFF. In this work, several prototypes working around 40 K were designed, built and characterized. They consist in a low temperature cell that contains the liquid neon connected to an expansion volume at room temperature for gas storage during the liquid evaporation phase. To turn this system insensitive to the gravity direction, the liquid is retained in the low temperature cell by capillary effect in a porous material. Thanks to pressure regulation of the liquid neon bath, 900 J were stored at 40K. The higher latent heat of the liquid and the inexistence of triple point transitions at 40 K turn the pressure control during the evaporation a versatile and compact alternative to an ESU working at the triple point transitions. A quite compact second prototype ESU directly connected to the cryocooler cold finger was tested as a temperature stabilizer. This device was able to stabilize the cryocooler temperature ((≈ 40K ±1 K) despite sudden heat bursts corresponding to twice the cooling power of the cryocooler. This thesis describes the construction of these devices as well as the tests performed. It is also shown that the thermal model developed to predict the thermal behaviour of these devices,implemented as a software, describes quite well the experimental results. Solutions to improve these devices are also proposed.

Ion exchange dependent electroactive phase content and electrical properties of poly(vinylidene fluoride)/Na(M)Y composites

Different metal-ion exchanged NaY zeolite, Na(M)Y, were used to prepare poly(vinylidene fluoride) based composites by solvent casting and melting crystallization. The effect of different metal ion-exchanged zeolites on polymer crystallization and electrical properties was reported. Cation-framework interactions and hydration energy of the cations determined that K+ is the most efficient exchanged ion in NaY zeolite, followed by Cs+ and Li+. The electroactive phase crystallization strongly depends on the ions present in the zeolite, leading to variations of the surface energy characteristics of the Na(M)Y zeolites and the polymer chain ability of penetration in the zeolite. Thus, Na(Li)Y and NaY induces the complete electroactive -phase crystallization of the crystalline phase of PVDF, while Na(K)Y only induces it partly and Na(Cs)Y is not able to promote the crystallization of the electroactive phase. Furthermore, different ion size/weigh and different interaction with the zeolite framework results in significant variations in the electrical response of the composite. In this way, iinterfacial polarization effects in the zeolite cavities and zeolite-polymer interface, leads to strong increases of the dielectric constant on the composites with lightest ions weakly bound to the zeolite framework. Polymer composite with Na(Li)Y show the highest dielectric response, followed by NaY and Na(K)Y. Zeolite Na(Cs)Y contribute to a decrease of the dielectric constant of the composite. The results show the relevance of the materials for sensor development.

The Gibbs free energy of formation of a glass alloy

A simple expression for the Gibbs free energy of formation of a pure component or a eutectic alloy glass, relative to the stable crystalline phase (or phases) at the same temperature is deduced by use of thermodynamic arguments. The expression obtained is supposed to apply to both monocomponent and multicomponent liquid alloys that might become glasses from the supercooled liquid state, irrespective of the critical cooling rate needed to avoid crystallization.