Porous Polymeric Bioresorbable Scaffolds for Tissue Engineering

Tissue engineering is a discipline that aims at regenerating damaged biological tissues by using a cell-construct engineered in vitro made of cells grown into a porous 3D scaffold. The role of the scaffold is to guide cell growth and differentiation by acting as a bioresorbable temporary substrate that will be eventually replaced by new tissue produced by cells. As a matter or fact, the obtainment of a successful engineered tissue requires a multidisciplinary approach that must integrate the basic principles of biology, engineering and material science. The present Ph.D. thesis aimed at developing and characterizing innovative polymeric bioresorbable scaffolds made of hydrolysable polyesters. The potentialities of both commercial polyesters (i.e. poly-e-caprolactone, polylactide and some lactide copolymers) and of non-commercial polyesters (i.e. poly-w-pentadecalactone and some of its copolymers) were explored and discussed. Two techniques were employed to fabricate scaffolds: supercritical carbon dioxide (scCO2) foaming and electrospinning (ES). The former is a powerful technology that enables to produce 3D microporous foams by avoiding the use of solvents that can be toxic to mammalian cells. The scCO2 process, which is commonly applied to amorphous polymers, was successfully modified to foam a highly crystalline poly(w-pentadecalactone-co-e-caprolactone) copolymer and the effect of process parameters on scaffold morphology and thermo-mechanical properties was investigated. In the course of the present research activity, sub-micrometric fibrous non-woven meshes were produced using ES technology. Electrospun materials are considered highly promising scaffolds because they resemble the 3D organization of native extra cellular matrix. A careful control of process parameters allowed to fabricate defect-free fibres with diameters ranging from hundreds of nanometers to several microns, having either smooth or porous surface. Moreover, versatility of ES technology enabled to produce electrospun scaffolds from different polyesters as well as “composite” non-woven meshes by concomitantly electrospinning different fibres in terms of both fibre morphology and polymer material. The 3D-architecture of the electrospun scaffolds fabricated in this research was controlled in terms of mutual fibre orientation by properly modifying the instrumental apparatus. This aspect is particularly interesting since the micro/nano-architecture of the scaffold is known to affect cell behaviour. Since last generation scaffolds are expected to induce specific cell response, the present research activity also explored the possibility to produce electrospun scaffolds bioactive towards cells. Bio-functionalized substrates were obtained by loading polymer fibres with growth factors (i.e. biomolecules that elicit specific cell behaviour) and it was demonstrated that, despite the high voltages applied during electrospinning, the growth factor retains its biological activity once released from the fibres upon contact with cell culture medium. A second fuctionalization approach aiming, at a final stage, at controlling cell adhesion on electrospun scaffolds, consisted in covering fibre surface with highly hydrophilic polymer brushes of glycerol monomethacrylate synthesized by Atom Transfer Radical Polymerization. Future investigations are going to exploit the hydroxyl groups of the polymer brushes for functionalizing the fibre surface with desired biomolecules. Electrospun scaffolds were employed in cell culture experiments performed in collaboration with biochemical laboratories aimed at evaluating the biocompatibility of new electrospun polymers and at investigating the effect of fibre orientation on cell behaviour. Moreover, at a preliminary stage, electrospun scaffolds were also cultured with tumour mammalian cells for developing in vitro tumour models aimed at better understanding the role of natural ECM on tumour malignity in vivo.

Discotic materials for organic electronics

Die Kontrolle der ausgeprägten Aggregationsfähigkeit von alkylsubstituierten Hexa-peri-hexabenzocoronenen (HBC) wurde durch die Reduktion der intermolekularen Wechselwirkungen erreicht. Sterisch anspruchsvolle, verzweigte Alkylketten, mit einem Verzweigungspunkt naher des aromatischen Kerns, wurden in die Corona der aromatischen Scheiben eingebracht und verleihen den Derivaten Schmelzbarkeit ohne thermische Zersetzung. Dies erlaubte eine kostengünstige Verarbeitungstechniken direkt aus der Schmelze wie z.B. Zonenschmelzen, um uniaxial organisierte makroskopische Filme zu erhalten. Abhängig von dem sterischen Anspruch, der durch die Seitenkette erzeugt wird, wurden unterschiedliche molekulare Orientierungen auf Oberflächen erhalten, was eine wichtige Voraussetzung ist, um diskotische Materialien in elektronische Bauteile zu implementieren. Eine weitere Voraussetzung sind hohe Ladungsträgerbeweglichkeiten und Ladungsträgerlebenszeiten in den Halbleitermaterialien, die mit time-resolved pulse-radiolysis microwave conductivity (TR-PRMC) und time-of-flight (TOF) auch für die synthetisieren Materialien bestimmt wurden. Die neuen Materialien zeigten bereits in organischen Solarzellen gute Leistungen. Den Einfluss des Perimeters auf die elektronischen Eigenschaften der polyzyklischen aromatischen Kohlenwasserstoffe (PAKs) wurde theoretisch vorhergesagt und in dieser Arbeit durch die Synthese einer homologe Serie von PAKs experimentell bestätigt. Geht man von der „arm-chair“ Peripherie des HBC sukzessive zu einer partiellen „zickzack“ Peripherie, so findet man eine Abhängigkeit der elektronischen Banden von Symmetrie und Größe des aromatischen Systems. Die spontan ausgebildete Überstruktur dieser Derivate zeigte eine Abhängigkeit von Substitutionsmuster und der Natur der Alkylketten. Zusammenfassend wurden neben der Synthese von neuartigen Materialien für den Einsatz in der organischen Elektronik Synthesen entwickelt, die eine vereinfachte Funktionalisierung von ausgedehnten PAKs ermöglicht. Diese Konzepte erlauben eine Justierung der molekularen und supramolekularen Eigenschaften, eines der wichtigsten Voraussetzungen für den Einsatz von Materialien in elektronischen Bauelementen.

Klassische und ab initio Molekulardynamik-Untersuchungen zu Germaniumdioxidschmelzen

Germaniumdioxid (GeO2) ist ein Glasbildner, der wie das homologe SiO2 ein ungeordnetes tetraedrisches Netzwerk ausbildet. In dieser Arbeit werden mit Hilfe von Molekulardynamik-Computersimulationen die Struktur und Dynamik von GeO2 in Abhängigkeit von der Temperatur untersucht. Dazu werden sowohl Simulationen mit einem klassischen Paarpotentialmodell von Oeffner und Elliott als auch ab initio-Simulationen gemäß der Car-Parrinello-Molekulardynamik (CPMD), bei der elektronische Freiheitsgrade mittels Dichtefunktionaltheorie beschrieben werden, durchgeführt. In der klassischen Simulation werden dazu ein Temperaturen zwischen 6100 K und 2530 K betrachtet. Darüberhinaus ermöglichen Abkühlläufe auf T=300 K das Studium der Struktur des Glases. Zum Vergleich werden CPMD-Simulationen für kleinere Systeme mit 60 bzw. 120 Teilchen bei den Temperaturen 3760 K und 3000 K durchgeführt. In den klassischen Simulationen kann die im Experiment bis 1700 K nachgewiesene, im Vergleich zu SiO2 starke, Temperaturabhängigkeit der Dichte auch bei höheren Temperaturen beobachtet werden. Gute Übereinstimmungen der Simulationen mit experimentellen Daten zeigen sich bei der Untersuchung verschiedener struktureller Größen, wie z.B. Paarkorrelationsfunktionen, Winkelverteilungen, Koordinationszahlen und Strukturfaktoren. Es können leichte strukturelle Abweichungen der CPMD-Simulationen von den klassischen Simulationen aufgezeigt werden: 1. Die Paarabstände in CPMD sind durchweg etwas kleiner. 2. Es zeigt sich, daß die Bindungen in den ab initio-Simulationen weicher sind, was sich auch in einer etwas stärkeren Temperaturabhängigkeit der strukturellen Größen im Vergleich zu den klassischen Simulationen niederschlägt. 3. Für CPMD kann ein vermehrtes Auftreten von Dreierringstrukturen gezeigt werden. 4. In der CPMD werden temperaturabhängige Defektstrukturen in Form von Sauerstoffpaaren beobachtet, die vor allem bei 3760 K, kaum jedoch bei 3000 K auftreten. Alle strukturellen Unterschiede zwischen klassischer und CPMD-Simulation sind eindeutig nicht auf Finite-Size-Effekte aufgrund der kleinen Systemgrößen in den CPMD-Simulationen zurückzuführen, d.h. sie sind tatsächlich methodisch bedingt. Bei der Dynamik von GeO2 wird in den klassischen Simulationen ebenfalls eine gute Übereinstimmung mit experimentellen Daten beobachtet, was ein Vergleich der Diffusionskonstanten mit Viskositätsmessungen bei hohen Temperaturen belegt. Die Diffusionskonstanten zeigen teilweise ein verschiedenes Verhalten zum homologen SiO2. Sie folgen in GeO2 bei Temperaturen unter 3000 K einem Arrheniusgesetz mit einer deutlich niedrigeren Aktivierungsenergie. Darüberhinaus werden die Möglichkeiten der Parametrisierung eines neuen klassischen Paarpotentials mittels der Kräfte entlang der CPMD-Trajektorien untersucht. Es zeigt sich, daß derartige Parametrisierungen sehr stark von den gewählten Startparametern abhängen. Ferner führen sämtliche an die Schmelze parametrisierten Potentiale zu zu hohen Dichten im Vergleich zum Experiment. Zum einen liegt dies sehr wahrscheinlich daran,daß für das System GeO2 Kraftdaten allein nicht ausreichen, um grundlegende strukturelle Größen, wie z.B. Paarkorrelationen und Winkelverteilungen, der CPMD-Simulationen gut reproduzieren zu können. Zum anderen ist wohl die Beschreibung mittels Paarpotentialen nicht ausreichend und es ist erforderlich, Merkörperwechselwirkungen in Betracht zu ziehen.

Characterization of inclusions and their distribution in natural and artificial samples by synchrotron cryo-micro-tomography (SCXRT)

Zur geometrischen Vermessung und Beschreibung von Einschlüssen in natürlichen sowie im Labor geschaffenen Eispartikeln wurde ein neuartiger Versuchaufbau an der Tomographie-Endstation der Material Science Beam Line an der Swiss Light Source (SLS, Paul Scherrer Institut, Villigen, Schweiz) entwickelt. Dieser besteht aus einer Plexiglas-Tasse und einem doppelwandigen Kaptonfolien-Käfig, der wiederum auf die Düse eines CryojetXL (Oxford Instruments) montiert wurde. Abgesehen von dem hohen Maß an Flexibilit¨at bez¨uglich der Installation erlaubt es dieser Aufbau, die Temperatur des Experiments mit einer Genauigkeit von ± 1 K über einen Bereich von 271 K bis 220 K zu regeln. In den hier beschriebenen Experimenten wurde eine räumliche Auflösung von 1.4 µm erzielt.

New eco-friendly polyesters from renewable resources

Nowadays the development of sustainable polymers, with convenient properties to substitute the traditional petroleum-based materials, is one of the major issues for material science. The utilization of renewable resources as feedstock for biopolyesters is a challenging target.The research work described in the present thesis is strictly connected to these urgent necessities and is focused mainly in finding new biopolymers, in particular biopolyesters, which are obtainable from biomass and characterized by a wide range of properties, in order to potentially substitute polyolefins and aromatic polyesters (for example, poly(ethylene terephthalate))

From lipid bilayers to synaptic vesicles: atomic force microscopy on lipid-based systems

The aim of this thesis was to apply the techniques of the atomic force microscope (AFM) to biological samples, namely lipid-based systems. To this end several systems with biological relevance based on self-assembly, such as a solid-supported membrane (SSM) based sensor for transport proteins, a bilayer of the natural lipid extract from an archaebacterium, and synaptic vesicles, were investigated by the AFM. For the characterization of transport proteins with SSM-sensors proteoliposomes are adsorbed that contain the analyte (transport protein). However the forces governing bilayer-bilayer interactions in solution should be repulsive under physiological conditions. I investigated the nature of the interaction forces with AFM force spectroscopy by mimicking the adsorbing proteoliposome with a cantilever tip, which was functionalized with charged alkane thiols. The nature of the interaction is indeed repulsive, but the lipid layers assemble in stacks on the SSM, which expose their unfavourable edges to the medium. I propose a model by which the proteoliposomes interact with these edges and fuse with the bilayer stacks, so forming a uniform layer on the SSM. Furthermore I characterized freestanding bilayers from a synthetic phospholipid with a phase transition at 41°C and from a natural lipid extract of the archaebacterium Methanococcus jannaschii. The synthetic lipid is in the gel-phase at room temperature and changes to the fluid phase when heated to 50°C. The bilayer of the lipid extract shows no phase transition when heated from room temperature to the growth temperature (~ 50°C) of the archeon. Synaptic vesicles are the containers of neurotransmitter in nerve cells and the synapsins are a family of extrinsic membrane proteins, that are associated with them, and believed to control the synaptic vesicle cycle. I used AFM imaging and force spectroscopy together with dynamic light scattering to investigate the influence of synapsin I on synaptic vesicles. To this end I used native, untreated synaptic vesicles and compared them to synapsin-depleted synaptic vesicles. Synapsin-depleted vesicles were larger in size and showed a higher tendency to aggregate compared to native vesicles, although their mechanical properties were alike. I also measured the aggregation kinetics of synaptic vesicles induced by synapsin I and found that the addition of synapsin I promotes a rapid aggregation of synaptic vesicles. The data indicate that synapsin I affects the stability and the aggregation state of synaptic vesicles, and confirm the physiological role of synapsins in the assembly and regulation of synaptic vesicle pools within nerve cells.

Structure and properties of intermetallic ternary rare earth compounds

The so called material science is an always growing field in modern research. For the development of new materials not only the experimental characterization but also theoretical calculation of the electronic structure plays an important role. A class of compounds that has attracted a great deal of attention in recent years is known as REME compounds. These compounds are often referred to with RE designating rare earth, actinide or an element from group 1 - 4, M representing a late transition metal from groups 8 - 12, and E belonging to groups 13 - 15. There are more than 2000 compounds with 1:1:1 stoichiometry belonging to this class of compounds and they offer a broad variety of different structure types. Although many REME compounds are know to exist, mainly only structure and magnetism has been determined for these compounds. In particular, in the field of electronic and transport properties relatively few efforts have been made. The main focus in this study is on compounds crystallizing in MgAgAs and LiGaGe structure. Both structures can only be found among 18 valence electron compounds. The f electrons are localized and therefor not count as valence electrons. A special focus here was also on the magnetoresistance effects and spintronic properties found among the REME compounds. An examination of the following compounds was made: GdAuE (E = In, Cd, Mg), GdPdSb, GdNiSb, REAuSn (RE = Gd, Er, Tm) and RENiBi (RE = Pr, Sm, Gd - Tm, Lu). The experimental results were compared with theoretic band structure calculations. The first half metallic ferromagnet with LiGaGe structure (GdPdSb) was found. All semiconducting REME compounds with MgAgAs structure show giant magnetoresistance (GMR) at low temperatures. The GMR is related to a metal-insulator transition, and the value of the GMR depends on the value of the spin-orbit coupling. Inhomogeneous DyNiBi samples show a small positive MR at low temperature that depends on the amount of metallic impurities. At higher fields the samples show a negative GMR. Inhomogeneous nonmagnetic LuNiBi samples show no negative GMR, but a large positive MR of 27.5% at room temperature, which is interesting for application.

Luminescent Ionic Transition-Metal Complexes for Light-Emitting Electrochemical Cells

The European Union set the ambitious target of reducing energy consumption by 20% within 2020. This goal demands a tremendous change in how we generate and consume energy and urgently calls for an aggressive policy on energy efficiency. Since 19% of the European electrical energy is used for lighting, considerable savings can be achieved with the development of novel and more efficient lighting systems. In this thesis, accomplished in the frame of the EU project CELLO, I report some selected goals we achieved attempting to develop highly efficient, flat, low cost and flexible light sources using Light-Emitting Electrochemical Cells (LECs), based on ionic cyclometalated iridium(III) complexes. After an extensive introduction about LECs and solid-state lighting in general, I focus on the research we carried out on cyclometalated iridium(III) complexes displaying deep-blue emission, which has turned out to be a rather challenging task. In order to demonstrate the wide versatility of this class of compounds, I also report a case in which some tailored iridium(III) complexes act as near-infrared (NIR) sources. In fact, standard NIR emitting devices are typically expensive and, also in this case, LECs could serve as low-cost alternatives in fields were NIR luminescence is crucial, such as telecommunications and bioimaging. Since LECs are based on only one active material, in the last chapter I stress the importance of an integrated approach toward the right selection of suitable emitters not only from the photophysical, but also from the point of view of material science. An iridium(III) complex, once in the device, is interacting with ionic liquids, metal cathodes, electric fields, etc. All these interactions should be taken in to account if Europe really wants to implement more efficient lighting paradigms, generating light beyond research labs.

Design, Synthesis and Characterization of N-Containing Organic Compounds

The needed of new intermediates/products for screening in the fields of drug discovery and material science is the driving force behind the development of new methodologies and technologies. Organic scaffolds are privileged targets for this scouting. Among them a priority place must be attributed to those including nitrogen functionalities in their scaffolds. It comes out that new methodologies, allowing the introduction of the nitrogen atom for the synthesis of an established target or for the curiosity driven researches, will always be welcome. The target of this PhD Thesis’ work is framed within this goal. Accordingly, Chapter 1 reports the preparation of new N-Heteroarylmethyl 3-carboxy-5-hydroxy piperidine scaffold, as potential and selective α-glucosidase inhibitors. The proposed reversible uncompetitive mechanism of inhibition makes them attractive as interesting candidate for drug development. Chapter 2 is more environmentally method-driven research. Eco-friendly studies on the synthesis of enantiomerically pure 1,4-dihydropyridines using “solid” ammonia (magnesium nitride) is reported via classical Hantzch method. Chapter 3 and Chapter 4 may be targeted as the core of the Thesis’s research work. Chapter 3 reports the studies addressed to the synthesis of N-containing heterocycles by using N-trialkylsilylimine/hetero-Diels–Alder (HAD) approach. New eco-friendly methodology as MAOS (Microwave Assisted Organic Synthesis) has been used as witness of our interest to a sustainable chemistry. Theoretical calculations were adopted to fully clarify the reaction mechanism. Chapter 4 is dedicated to picture the most recent studies performed on the application of N-Metallo-ketene imines (metallo= Si, Sn, Al), relatively new intermediates which are becoming very popular, in the preparation of highly functionalized N-containing derivatives, accordingly to the Thesis’ target. Derivatives obtained are designed in such a way that they could be of interest in the field of drug and new material chemistry.

High frequency acoustics in colloid-based meso- and nanostructures by spontaneous Brillouin light scattering

Materials that can mold the flow of elastic waves of certain energy in certain directions are called phononic materials. The present thesis deals essentially with such phononic systems, which are structured in the mesoscale (<1 µm), and with their individual components. Such systems show interesting phononic properties in the hypersonic region, i.e., at frequencies in the GHz range. It is shown that colloidal systems are excellent model systems for the realization of such phononic materials. Therefore, different structures and particle architectures are investigated by Brillouin light scattering, the inelastic scattering of light by phonons.rnThe experimental part of this work is divided into three chapters: Chapter 4 is concerned with the localized mechanical waves in the individual spherical colloidal particles, i.e., with their resonance- or eigenvibrations. The investigation of these vibrations with regard to the environment of the particles, their chemical composition, and the influence of temperature on nanoscopically structured colloids allows novel insights into the physical properties of colloids at small length scales. Furthermore, some general questions concerning light scattering on such systems, in dispute so far, are convincingly addressed.rnChapter 5 is a study of the traveling of mechanical waves in colloidal systems, consisting of ordered and disordered colloids in liquid or elastic matrix. Such systems show acoustic band gaps, which can be explained geometrically (Bragg gap) or by the interaction of the acoustic band with the eigenvibrations of the individual spheres (hybridization gap).rnWhile the latter has no analogue in photonics, the presence of strong phonon scatterers, when a large elastic mismatch between the composite components exists, can largely impact phonon propagation in analogy to strong multiple light scattering systems. The former is exemplified in silica based phononic structures that opens the door to new ways of sound propagation manipulation.rnChapter 6 describes the first measurement of the elastic moduli in newly fabricated by physical vapor deposition so-called ‘stable organic glasses’. rnIn brief, this thesis explores novel phenomena in colloid-based hypersonic phononic structures, utilizing a versatile microfabrication technique along with different colloid architectures provided by material science, and applying a non-destructive optical experimental tool to record dispersion diagrams.rn

Fracture Mechanics Investigation of Structures with Defects

Fracture mechanics plays an important role in the material science, structure design and industrial production due to the failure of materials and structures are paid high attention in human activities. This dissertation, concentrates on some of the fractural aspects of shaft and composite which have being increasingly used in modern structures, consists four chapters within two parts. Chapters 1 to 4 are included in part 1. In the first chapter, the basic knowledge about the stress and displacement fields in the vicinity of a crack tip is introduced. A review involves the general methods of calculating stress intensity factors are presented. In Chapter 2, two simple engineering methods for a fast and close approximation of stress intensity factors of cracked or notched beams under tension, bending moment, shear force, as well as torque are presented. New formulae for calculating the stress intensity factors are proposed. One of the methods named Section Method is improved and applied to the three dimensional analysis of cracked circular section for calculating stress intensity factors. The comparisons between the present results and the solutions calculated by ABAQUS for single mode and mixed mode are studied. In chapter 3, fracture criteria for a crack subjected to mixed mode loading of two-dimension and three-dimension are reviewed. The crack extension angle for single mode and mixed mode, and the critical loading domain obtained by SEDF and MTS are compared. The effects of the crack depth and the applied force ratio on the crack propagation angle and the critical loading are investigated. Three different methods calculating the crack initiation angle for three-dimension analysis of various crack depth and crack position are compared. It should be noted that the stress intensity factors used in the criteria are calculated in section 2.1.

Synthese organisch-anorganischer Hybridmaterialien in Mehrkomponenten-Lösungsmittelsystemen

Die Funktionalisierung anorganischer Nanopartikel stellt einen Schlüsselschritt in der Herstellung von Nanokompositen dar. Nanokomposite erzielen ein wachsendes Interesse im Bereich der Polymer- und der Materialwissenschaften, da die Kombination mehrerer Materialien mit unterschiedlichen Eigenschaften, wie etwa die Kombination anorganischer Nanopartikel mit Polymeren, große Synergieeffekte erhoffen lässt.rnrnDer Einbau anorganischer Nanopartikel in polymere Matrixmaterialien zur Verbesserung oder Einführung mechanischer, optischer oder magnetischer Eigenschaften von Polymeren bedarf allerdings der Modifizierung der Oberfläche des anorganischen Materials, um die für die positiven Synergieeffekte essentielle Kompatibilität zwischen Füllstoff und Matrix zu erreichen.rnrnEine Vielzahl anorganischer Partikel ist bereits als wässrige Dispersion erhältlich (SiO2, Al2O3, CeO2, ZrO2, ...). Mehrkomponenten- Lösungsmittelsysteme ermöglichen den Transfer dieser Partikel in eine unpolare Umgebung und gleichzeitig deren Funktionalisierung mit amphiphilen Copolymeren. Aufgrund der reversiblen Schaltbarkeit dieser Lösungsmittelsysteme zwischen einem einphasigen und zweiphasigen Zustand werden die zu Beginn in zwei nichtmischbaren Phasen vorliegenden Reaktionspartner durch Übergang in einen einphasigen Zustand unter homogenen Bedingungen in Kontakt gebracht und durch eine erneute Phasentrennung isoliert.rnEin weiterer Vorteil dieser Lösungsmittelsysteme ist deren Tolerierung funktioneller Gruppen in den verwendeten amphiphilen Copolymeren, welche nicht in Wechselwirkung mit der Partikeloberfläche stehen. Beispielsweise können Amine in den amphiphilen Copolymeren für die Wechselwirkung der funktionalisierten Partikel mit einer Polyurethanmatrix dienen, Alkine können mittels einer 1,3-dipolaren Cycloaddition umgesetzt werden oder aber perfluorierten Seitenketten in den Seitenketten der amphiphilen Copolymere die Kompatibilisierung der funktionalisierten Partikel mit einem perfluorierten Polymer gewährleisten.

Thermodynamic concepts in adaptive resolution simulations

This thesis work is devoted to the conceptual and technical development of the Adaptive Resolution Scheme (AdResS), a molecular dynamics method that allows the simulation of a system with different levels of resolution simultaneously. The simulation domain is divided into high and low resolution zones and a transition region that links them, through which molecules can freely diffuse.rnThe first issue of this work regards the thermodynamic consistency of the method, which is tested and verified in a model liquid of tetrahedral molecules. The results allow the introduction of the concept of the Thermodynamic Force, an external field able to correct spurious density fluctuations present in the transition region in usual AdResS simulations.rnThe AdResS is also applied to a system where two different representations with the same degree of resolution are confronted. This simple test extends the method from an Adaptive Resolution Scheme to an Adaptive Representation Scheme, providing a way of coupling different force fields based on thermodynamic consistency arguments. The Thermodynamic Force is successfully applied to the example described in this work as well.rnAn alternative approach of deducing the Thermodynamic Force from pressure consistency considerations allows the interpretation of AdResS as a first step towards a molecular dynamics simulation in the Grand Canonical ensemble. Additionally, such a definition leads to a practical way of determining the Thermodynamic Force, tested in the well studied tetrahedral liquid. The effects of AdResS and this correction on the atomistic domain are analyzed by inspecting the local distribution of velocities, radial distribution functions, pressure and particle number fluctuation. Their comparison with analogous results coming from purely atomistic simulations shows good agreement, which is greatly improved under the effect of the external field.rnA further step in the development of AdResS, necessary for several applications in biophysics and material science, consists of its application to multicomponent systems. To this aim, the high-resolution representation of a model binary mixture is confronted with its coarse-grained representation systematically parametrized. The Thermodynamic Force, whose development requires a more delicate treatment, also gives satisfactory results.rnFinally, AdResS is tested in systems including two-body bonded forces, through the simulation of a model polymer allowed to adaptively change its representation. It is shown that the distribution functions that characterize the polymer structure are in practice not affected by the change of resolution.rnThe technical details of the implementation of AdResS in the ESPResSo package conclude this thesis work.

Approcci innovativi per future simulazioni HPC verso l'exascale computing con applicazioni alla fisica dei materiali

In questo lavoro di tesi sono state evidenziate alcune problematiche relative alle macchine exascale (sistemi che sviluppano un exaflops di Potenza di calcolo) e all'evoluzione dei software che saranno eseguiti su questi sistemi, prendendo in esame principalmente la necessità del loro sviluppo, in quanto indispensabili per lo studio di problemi scientifici e tecnologici di più grandi dimensioni, con particolare attenzione alla Material Science, che è uno dei campi che ha avuto maggiori sviluppi grazie all'utilizzo di supercomputer, ed ad uno dei codici HPC più utilizzati in questo contesto: Quantum ESPRESSO. Dal punto di vista del software sono state presentate le prime misure di efficienza energetica su architettura ibrida grazie al prototipo di cluster EURORA sul software Quantum ESPRESSO. Queste misure sono le prime ad essere state pubblicate nel contesto software per la Material Science e serviranno come baseline per future ottimizzazioni basate sull'efficienza energetica. Nelle macchine exascale infatti uno dei requisiti per l'accesso sarà la capacità di essere energeticamente efficiente, così come oggi è un requisito la scalabilità del codice. Un altro aspetto molto importante, riguardante le macchine exascale, è la riduzione del numero di comunicazioni che riduce il costo energetico dell'algoritmo parallelo, poiché in questi nuovi sistemi costerà di più, da un punto di vista energetico, spostare i dati che calcolarli. Per tale motivo in questo lavoro sono state esposte una strategia, e la relativa implementazione, per aumentare la località dei dati in uno degli algoritmi più dispendiosi, dal punto di vista computazionale, in Quantum ESPRESSO: Fast Fourier Transform (FFT). Per portare i software attuali su una macchina exascale bisogna iniziare a testare la robustezza di tali software e i loro workflow su test case che stressino al massimo le macchine attualmente a disposizione. In questa tesi per testare il flusso di lavoro di Quantum ESPRESSO e WanT, un software per calcolo di trasporto, è stato caratterizzato un sistema scientificamente rilevante costituito da un cristallo di PDI - FCN2 che viene utilizzato per la costruzione di transistor organici OFET. Infine è stato simulato un dispositivo ideale costituito da due elettrodi in oro con al centro una singola molecola organica.

Prediction of new crystal structures under extreme conditions

One of the most important challenges in chemistry and material science is the connection between the contents of a compound and its chemical and physical properties. In solids, these are greatly influenced by the crystal structure.rnrnThe prediction of hitherto unknown crystal structures with regard to external conditions like pressure and temperature is therefore one of the most important goals to achieve in theoretical chemistry. The stable structure of a compound is the global minimum of the potential energy surface, which is the high dimensional representation of the enthalpy of the investigated system with respect to its structural parameters. The fact that the complexity of the problem grows exponentially with the system size is the reason why it can only be solved via heuristic strategies.rnrnImprovements to the artificial bee colony method, where the local exploration of the potential energy surface is done by a high number of independent walkers, are developed and implemented. This results in an improved communication scheme between these walkers. This directs the search towards the most promising areas of the potential energy surface.rnrnThe minima hopping method uses short molecular dynamics simulations at elevated temperatures to direct the structure search from one local minimum of the potential energy surface to the next. A modification, where the local information around each minimum is extracted and used in an optimization of the search direction, is developed and implemented. Our method uses this local information to increase the probability of finding new, lower local minima. This leads to an enhanced performance in the global optimization algorithm.rnrnHydrogen is a highly relevant system, due to the possibility of finding a metallic phase and even superconductor with a high critical temperature. An application of a structure prediction method on SiH12 finds stable crystal structures in this material. Additionally, it becomes metallic at relatively low pressures.