Approximate ab initio calculations and the method of molecular fragments

A two stage approach to performing ab initio calculations on medium and large sized molecules is described. The first step is to perform SCF calculations on small molecules or molecular fragments using the OPIT Program. This employs a small basis set of spherical and p-type Gaussian functions. The Gaussian functions can be identified very closely with atomic cores, bond pairs, lone pairs, etc. The position and exponent of any of the Gaussian functions can be varied by OPIT to produce a small but fully optimised basis set. The second stage is the molecular fragments method. As an example of this, Gaussian exponents and distances are taken from an OPIT calculation on ethylene and used unchanged in a single SCF calculation on benzene. Approximate ab initio calculations of this type give much useful information and are often preferable to semi-empirical approaches, since the nature of the approximations involved is much better defined.

Tile-Based Two-Dimensional Phase Unwrapping for Digital Holography Using a Modular Framework

A variety of physical and biomedical imaging techniques, such as digital holography, interferometric synthetic aperture radar (InSAR), or magnetic resonance imaging (MRI) enable measurement of the phase of a physical quantity additionally to its amplitude. However, the phase can commonly only be measured modulo 2π, as a so called wrapped phase map. Phase unwrapping is the process of obtaining the underlying physical phase map from the wrapped phase. Tile-based phase unwrapping algorithms operate by first tessellating the phase map, then unwrapping individual tiles, and finally merging them to a continuous phase map. They can be implemented computationally efficiently and are robust to noise. However, they are prone to failure in the presence of phase residues or erroneous unwraps of single tiles. We tried to overcome these shortcomings by creating novel tile unwrapping and merging algorithms as well as creating a framework that allows to combine them in modular fashion. To increase the robustness of the tile unwrapping step, we implemented a model-based algorithm that makes efficient use of linear algebra to unwrap individual tiles. Furthermore, we adapted an established pixel-based unwrapping algorithm to create a quality guided tile merger. These original algorithms as well as previously existing ones were implemented in a modular phase unwrapping C++ framework. By examining different combinations of unwrapping and merging algorithms we compared our method to existing approaches. We could show that the appropriate choice of unwrapping and merging algorithms can significantly improve the unwrapped result in the presence of phase residues and noise. Beyond that, our modular framework allows for efficient design and test of new tile-based phase unwrapping algorithms. The software developed in this study is freely available.

Finite difference calculations of permeability in large domains in a wide porosity range.

Determining effective hydraulic, thermal, mechanical and electrical properties of porous materials by means of classical physical experiments is often time-consuming and expensive. Thus, accurate numerical calculations of material properties are of increasing interest in geophysical, manufacturing, bio-mechanical and environmental applications, among other fields. Characteristic material properties (e.g. intrinsic permeability, thermal conductivity and elastic moduli) depend on morphological details on the porescale such as shape and size of pores and pore throats or cracks. To obtain reliable predictions of these properties it is necessary to perform numerical analyses of sufficiently large unit cells. Such representative volume elements require optimized numerical simulation techniques. Current state-of-the-art simulation tools to calculate effective permeabilities of porous materials are based on various methods, e.g. lattice Boltzmann, finite volumes or explicit jump Stokes methods. All approaches still have limitations in the maximum size of the simulation domain. In response to these deficits of the well-established methods we propose an efficient and reliable numerical method which allows to calculate intrinsic permeabilities directly from voxel-based data obtained from 3D imaging techniques like X-ray microtomography. We present a modelling framework based on a parallel finite differences solver, allowing the calculation of large domains with relative low computing requirements (i.e. desktop computers). The presented method is validated in a diverse selection of materials, obtaining accurate results for a large range of porosities, wider than the ranges previously reported. Ongoing work includes the estimation of other effective properties of porous media.

Raman Spectroscopy and DFT calculations for the Electronic Structure Characterization of Conjugated Polymers

In the last three decades, there has been a broad academic and industrial interest in conjugated polymers as semiconducting materials for organic electronics. Their applications in polymer light-emitting diodes (PLEDs), polymer solar cells (PSCs), and organic field-effect transistors (OFETs) offer opportunities for the resolution of energy issues as well as the development of display and information technologies1. Conjugated polymers provide several advantages including low cost, light weight, good flexibility, as well as solubility which make them readily processed and easily printed, removing the conventional photolithography for patterning2. A large library of polymer semiconductors have been synthesized and investigated with different building blocks, such as acenes or thiophene and derivatives, which have been employed to design new materials according to individual demands for specific applications. To design ideal conjugated polymers for specific applications, some general principles should be taken into account, including (i) side chains (ii) molecular weights, (iii) band gap and HOMO and LUMO energy levels, and (iv) suited morphology.3-6 The aim of this study is to elucidate the impact that substitution exerts on the molecular and electronic structure of π-conjugated polymers with outstanding performances in organic electronic devices. Different configurations of the π-conjugated backbones are analyzed: (i) donor-acceptor configuration, (ii) 1D lineal or 2D branched conjugated backbones, and (iii) encapsulated polymers (see Figure 1). Our combined vibrational spectroscopy and DFT study shows that small changes in the substitution pattern and in the molecular configuration have a strong impact on the electronic characteristics of these polymers. We hope this study can advance useful structure-property relationships of conjugated polymers and guide the design of new materials for organic electronic applications.

Thermodynamic calculations for the salt crystallisation damage in porous built heritage using PHREEQC

This work considers the crystallisation mechanisms of the most common and aggressive salts that generate stress in porous building stones as a result of changing ambient conditions. These mechanisms include the salt crystallisation that result from decreasing relative humidity and changes in temperature and, in hydrated salts, the dissolution of the lower hydrated form and the subsequent precipitation of the hydrated salt. We propose a new methodology for thermodynamic calculations using PHREEQC that includes these crystallisation mechanisms. This approach permits the calculation of the equilibrium relative humidity and the parameterization of the critical relative humidity and crystallisation pressures for the dissolution–precipitation transitions. The influence of other salts on the effectives of salt crystallisation and chemical weathering is also assessed. We review the sodium and magnesium sulphate and sodium chloride systems, in both single and multicomponent solutions, and they are compared to the sodium carbonate and calcium carbonate systems. The variation of crystallisation pressure, the formation of new minerals and the chemical dissolution by the presence of other salts is also evaluated. Results for hydrated salt systems show that high crystallisation pressures are possible as lower hydrated salts dissolve and more hydrated salts precipitate. High stresses may be also produced by decreasing temperature, although it requires that porous materials are wet for long periods of time. The presence of other salts changes the temperature and relative humidity of salt transitions that generates stress rather than reducing the pressure of crystallisation, if any salt has previously precipitated. Several practical conclusions derive from proposed methodology and provide conservators and architects with information on the potential weathering activity of soluble salts. Furthermore, the model calculations might be coupled with projections of future climate to give as improved understanding of the likely changes in the frequency of phase transitions in salts within porous stone.

Designing intrinsically photostable low band gap polymers:a smart tool combining EPR spectroscopy and DFT calculations

A rapid and efficient method to identify the weak points of the complex chemical structure of low band gap (LBG) polymers, designed for efficient solar cells, when submitted to light exposure is reported. This tool combines Electron Paramagnetic Resonance (EPR) using the 'spin trapping method' coupled with density functional theory modelling (DFT). First, the nature of the short life-time radicals formed during the early-stages of photo-degradation processes are determined by a spin-trapping technique. Two kinds of short life-time radical (R and R′O) are formed after 'short-duration' illumination in an inert atmosphere and in ambient air, respectively. Second, simulation allows the identification of the chemical structures of these radicals revealing the most probable photochemical process, namely homolytical scission between the Si atom of the conjugated skeleton and its pendent side-chains. Finally, DFT calculations confirm the homolytical cleavage observed by EPR, as well as the presence of a group that is highly susceptible to photooxidative attack. Therefore, the synergetic coupling of a spin trapping method with DFT calculations is shown to be a rapid and efficient method for providing unprecedented information on photochemical mechanisms. This approach will allow the design of LBG polymers without the need to trial the material within actual solar cell devices, an often long and costly screening procedure.

First-Principles Calculations of the Electronic and Structural Properties of GaSb

In this paper, we carried out first-principles calculations in order to investigate the structural and electronic properties of the binary compound gallium antimonide (GaSb). This theoretical study was carried out using the Density Functional Theory within the plane-wave pseudopotential method. The effects ofexchange and correlation (XC) were treated using the functional Local Density Approximation (LDA), generalized gradient approximation (GGA): Perdew–Burke–Ernzerhof (PBE), Perdew-Burke-Ernzerhof revised for solids (PBEsol), Perdew-Wang91 (PW91), revised Perdew–Burke–Ernzerhof (rPBE), Armiento–Mattson 2005 (AM05) and meta-generalized gradient approximation (meta-GGA): Tao–Perdew– Staroverov–Scuseria (TPSS) and revised Tao–Perdew–Staroverov–Scuseria (RTPSS) and modified Becke-Johnson (MBJ). We calculated the densities of state (DOS) and band structure with different XC potentials identified and compared them with the theoretical and experimental results reported in the literature. It was discovered that functional: LDA, PBEsol, AM05 and RTPSS provide the best results to calculate the lattice parameters (a) and bulk modulus (B0); while for the cohesive energy (Ecoh), functional: AM05, RTPSS and PW91 are closer to the values obtained experimentally. The MBJ, Rtpss and AM05 values found for the band gap energy is slightly underestimated with those values reported experimentally.

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

Sustainability and environmental impact of the Italian ceramic tile industry

The increasing consumption rates among citizens and the uncontrolled exploitation of natural resources have made environmental pollution and management of waste the main problems facing humanity in its upcoming future. Together with generation of energy and transport, industrial production certainly plays a key role in the genesis of these problems. It is for this reason that the concepts of environmental, social and economic sustainability have emerged over the years as the cornerstones for future development. In light of this, the most forward-looking industries have begun to study their impact on environment and society in order to improve their performances and, at the same time, to anticipate the increasingly rigorous environmental regulations. In this work, various performance indicators related to the Italian ceramic tile sector will be presented and discussed. In particular, the emission factor of characteristic pollutants will be reported on a period of up to fifteen years while data regarding waste management, concentration of pollutants and emission legal limits for the last decade will be here disclosed as a result of a vast analysis on recorded data. The collected information describes the present level of performance of the ceramic tile manufacturing industries in Italy and shows how recycling is now a consolidated reality and how some pollutants, such as particulate matter, fluorine and lead are actually disappearing from production processes and how others, such as volatile organic compounds, are increasing instead. Moreover, the adoption of alternative raw materials for the production of ceramic tiles is discussed and the implementation of the recycling of various waste is addressed at experimental or industrial scale. Finally, the development of a new ceramic engobe with high content of waste glass (20%) is presented as an experimental example of reutilization of resources in the ceramic tile industry.

Geopolymers for tile industry

This project aims to develop energy-saving methods for manufacturing tiles by using geopolymerization as an innovative processing technique. In fact, classical ceramic tiles are usually produced by sintering at temperatures ranging between 1100-1250 °C depending on natural raw materials selected for the mixture. Geopolymers, instead, are amorphous aluminosilicate inorganic polymers obtained by mixing alkali-silicate solutions with reactive aluminosilicate precursors, and geopolymerization may occur at room temperature or at temperature lower than 150 °C, and several types of waste may be used as geopolymer precursors. The possibility to obtain geopolymer tiles with properties comparable with those of ceramic tiles is extremely challenging, but geopolymerization would represent a huge energy-saving in the consolidation process as well as in natural raw materials exploitation.

Theoretical models and ab initio calculations of phonons

Questa tesi intende approfondire da un punto di vista, sia teorico sia computazionale, le proprietà fondamentali dei fononi. A tal fine, sono presentati i modelli quantistici di Einstein e di Debye che permettono la derivazione analitica degli osservabili macroscopici principali di un solido, come l’energia media e la capacità termica. Ciò è possibile tramite una trattazione meccano-statistica basata sull’approssimazione armonica dei modi normali di vibrazione degli ioni reticolari. Quindi, all’inizio si mostrano brevemente i risultati principali riguardanti l’oscillatore armonico quantistico. Successivamente, si approfondiscono i temi della dispersione fononica e della densità degli stati vibrazionali per reticoli cristallini 1D e 3D. Si ottiene che la prima non può essere considerata lineare se non nel limite di alte lunghezze d’onda, e che la seconda può presentare punti di singolarità correlati alla forma della relazione di dispersione. Infine, sono state svolte alcune analisi computazionali ab initio relative alla dispersione fononica, la densità degli stati vibrazionali e la frequenza di Debye del Carbonio (diamante) tramite i programmi VASP e Phonopy, confrontando i risultati con dati sperimentali presenti in letteratura.