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.

Near-field mediated enhancement effects on plasmonic nanostructures

In dieser Arbeit wird eine detaillierte Untersuchung und Charakterisierung der Zwei-Photonen-induzierten Fluoreszenzverstärkung von organischen Farbstoffen auf plasmonischen Nanostrukturen vorgestellt. Diese Fluoreszenzverstärkung ist insbesondere für hochaufgelöste Fluoreszenzmikroskopie und Einzelmolekülspektroskopie von großer Bedeutung. Durch die Zwei-Photonen-Anregung resultiert eine Begrenzung des Absorptionsprozesses auf das fokale Volumen. In Kombination mit dem elektrischen Nahfeld der Nanostrukturen als Anregungsquelle entsteht eine noch stärkere Verringerung des Anregungsvolumens auf eine Größe unterhalb der Beugungsgrenze. Dies erlaubt die selektive Messung ausgewählter Farbstoffe. Durch die Herstellung der Nanopartikel mittels Kolloidlithografie wird eine definierte, reproduzierbare Geometrie erhalten. Polymermultischichten dienen als Abstandshalter, um die Farbstoffe an einer exakten Distanz zum Metall zu positionieren. Durch die kovalente Anbindung des Farbstoffs an die oberste Schicht wird eine gleichmäßige Verteilung des Farbstoffs in geringer Konzentration erhalten. rnEs wird eine Verstärkung der Fluoreszenz um den Faktor 30 für Farbstoffe auf Goldellipsen detektiert, verglichen mit Farbstoffen außerhalb des Nahfelds. Sichelförmige Nanostrukturen erzeugen eine Verstärkung von 120. Dies belegt, dass das Ausmaß der Fluoreszenzverstärkung entscheidend von der Stärke des elektrischen Nahfelds der Nanostruktur abhängt. Auch das Material der Nanostruktur ist hierbei von Bedeutung. So erzeugen Silberellipsen eine 1,5-fach höhere Fluoreszenzverstärkung als identische Goldellipsen. Distanzabhängige Fluoreszenzmessungen zeigen, dass die Zwei-Photonen-angeregte Fluoreszenzverstärkung an strukturspezifischen Abständen zum Metall maximiert wird. Elliptische Strukturen zeigen ein Maximum bei einem Abstand von 8 nm zum Metall, wohingegen bei sichelförmigen Nanostrukturen die höchste Fluoreszenzintensität bei 12 nm gemessen wird. Bei kleineren Abständen unterliegt der Farbstoff einem starken Löschprozess, sogenanntes Quenching. Dieses konkurriert mit dem Verstärkungsprozess, wodurch es zu einer geringen Nettoverstärkung kommt. Hat die untersuchte Struktur Dimensionen größer als das Auflösungsvermögen des Mikroskops, ist eine direkte Visualisierung des elektrischen Nahfelds der Nanostruktur möglich. rnrnEin weiterer Fokus dieser Arbeit lag auf der Herstellung neuartiger Nanostrukturen durch kolloidlithografische Methoden. Gestapelte Dimere sichelförmiger Nanostrukturen mit exakter vertikaler Ausrichtung und einem Separationsabstand von etwa 10 nm wurden hergestellt. Die räumliche Nähe der beiden Strukturen führt zu einem Kopplungsprozess, der neue optische Resonanzen hervorruft. Diese können als Superpositionen der Plasmonenmoden der einzelnen Sicheln beschrieben werden. Ein Hybridisierungsmodell wird angewandt, um die spektralen Unterschiede zu erklären. Computersimulationen belegen die zugrunde liegende Theorie und erweitern das Modell um experimentell nicht aufgelöste Resonanzen. rnWeiterhin wird ein neuer Herstellungsprozess für sichelförmige Nanostrukturen vorgestellt, der eine präzise Formanpassung ermöglicht. Hierdurch kann die Lage der Plasmonenresonanz exakt justiert werden. Korrelationen der geometrischen Daten mit den Resonanzwellenlängen tragen zum grundlegenden Verständnis der Plasmonenresonanzen bei. Die vorgestellten Resultate wurden mittels Computersimulationen verifiziert. Der Fabrikationsprozess erlaubt die Herstellung von Dimeren sichelförmiger Nanostrukturen in einer Ebene. Durch die räumliche Nähe überlappen die elektrischen Nahfelder, wodurch es zu kopplungs-induzierten Shifts der Plasmonenresonanzen kommt. Der Unterschied zu theoretisch berechneten ungekoppelten Nanosicheln kann auch bei den gegenüberliegenden sichelförmigen Nanostrukturen mit Hilfe des Plasmonenhybridisierungsmodells erklärt werden.

Power Converters and Electric Drives for Smart Grid Applications

The present dissertation aims to explore, theoretically and experimentally, the problems and the potential advantages of different types of power converters for “Smart Grid” applications, with particular emphasis on multi-level architectures, which are attracting a rising interest even for industrial requests. The models of the main multilevel architectures (Diode-Clamped and Cascaded) are shown. The best suited modulation strategies to function as a network interface are identified. In particular, the close correlation between PWM (Pulse Width Modulation) approach and SVM (Space Vector Modulation) approach is highlighted. An innovative multilevel topology called MMC (Modular Multilevel Converter) is investigated, and the single-phase, three-phase and "back to back" configurations are analyzed. Specific control techniques that can manage, in an appropriate way, the charge level of the numerous capacitors and handle the power flow in a flexible way are defined and experimentally validated. Another converter that is attracting interest in “Power Conditioning Systems” field is the “Matrix Converter”. Even in this architecture, the output voltage is multilevel. It offers an high quality input current, a bidirectional power flow and has the possibility to control the input power factor (i.e. possibility to participate to active and reactive power regulations). The implemented control system, that allows fast data acquisition for diagnostic purposes, is described and experimentally verified.

Charge-transport simulations in organic semiconductors

In this thesis we have extended the methods for microscopic charge-transport simulations for organic semiconductors. In these materials the weak intermolecular interactions lead to spatially localized charge carriers, and the charge transport occurs as an activated hopping process between diabatic states. In addition to weak electronic couplings between these states, different electrostatic environments in the organic material lead to a broadening of the density of states for the charge energies which limits carrier mobilities.rnThe contributions to the method development includern(i) the derivation of a bimolecular charge-transfer rate,rn(ii) the efficient evaluation of intermolecular (outer-sphere) reorganization energies,rn(iii) the investigation of effects of conformational disorder on intramolecular reorganization energies or internal site energiesrnand (iv) the inclusion of self-consistent polarization interactions for calculation of charge energies.These methods were applied to study charge transport in amorphous phases of small molecules used in the emission layer of organic light emitting diodes (OLED).rnWhen bulky substituents are attached to an aromatic core in order to adjust energy levels or prevent crystallization, a small amount of delocalization of the frontier orbital to the substituents can increase electronic couplings between neighboring molecules. This leads to improved charge-transfer rates and, hence, larger charge-mobility. We therefore suggest using the mesomeric effect (as opposed to the inductive effect) when attaching substituents to aromatic cores, which is necessary for example in deep blue OLEDs, where the energy levels of a host molecule have to be adjusted to those of the emitter.rnFurthermore, the energy landscape for charges in an amorphous phase cannot be predicted by mesoscopic models because they approximate the realistic morphology by a lattice and represent molecular charge distributions in a multipole expansion. The microscopic approach shows that a polarization-induced stabilization of a molecule in its charged and neutral states can lead to large shifts, broadening, and traps in the distribution of charge energies. These results are especially important for multi-component systems (the emission layer of an OLED or the donor-acceptor interface of an organic solar cell), if the change in polarizability upon charging (or excitation in case of energy transport) is different for the components. Thus, the polarizability change upon charging or excitation should be added to the set of molecular parameters essential for understanding charge and energy transport in organic semiconductors.rnWe also studied charge transport in self-assembled systems, where intermolecular packing motives induced by side chains can increase electronic couplings between molecules. This leads to larger charge mobility, which is essential to improve devices such as organic field effect transistors, where low carrier mobilities limit the switching frequency.rnHowever, it is not sufficient to match the average local molecular order induced by the sidernchains (such as the pitch angle between consecutive molecules in a discotic mesophase) with maxima of the electronic couplings.rnIt is also important to make the corresponding distributions as narrow as possible compared to the window determined by the closest minima of thernelectronic couplings. This is especially important in one-dimensional systems, where charge transport is limited by the smallest electronic couplings.rnThe immediate implication for compound design is that the side chains should assist the self-assemblingrnprocess not only via soft entropic interactions, but also via stronger specific interactions, such as hydrogen bonding.rnrnrnrn

Advanced aspects of radiation protection in the use of particle accelerators in the medical field

In this work, the well-known MC code FLUKA was used to simulate the GE PETrace cyclotron (16.5 MeV) installed at “S. Orsola-Malpighi” University Hospital (Bologna, IT) and routinely used in the production of positron emitting radionuclides. Simulations yielded estimates of various quantities of interest, including: the effective dose distribution around the equipment; the effective number of neutron produced per incident proton and their spectral distribution; the activation of the structure of the cyclotron and the vault walls; the activation of the ambient air, in particular the production of 41Ar, the assessment of the saturation yield of radionuclides used in nuclear medicine. The simulations were validated against experimental measurements in terms of physical and transport parameters to be used at the energy range of interest in the medical field. The validated model was also extensively used in several practical applications uncluding the direct cyclotron production of non-standard radionuclides such as 99mTc, the production of medical radionuclides at TRIUMF (Vancouver, CA) TR13 cyclotron (13 MeV), the complete design of the new PET facility of “Sacro Cuore – Don Calabria” Hospital (Negrar, IT), including the ACSI TR19 (19 MeV) cyclotron, the dose field around the energy selection system (degrader) of a proton therapy cyclotron, the design of plug-doors for a new cyclotron facility, in which a 70 MeV cyclotron will be installed, and the partial decommissioning of a PET facility, including the replacement of a Scanditronix MC17 cyclotron with a new TR19 cyclotron.

Optimization of purification procedures in organic semiconductor synthesis for improving the performances of organic field-effect transistors

The aim of the research activity focused on the investigation of the correlation between the degree of purity in terms of chemical dopants in organic small molecule semiconductors and their electrical and optoelectronic performances once introduced as active material in devices. The first step of the work was addressed to the study of the electrical performances variation of two commercial organic semiconductors after being processed by means of thermal sublimation process. In particular, the p-type 2,2′′′-Dihexyl-2,2′:5′,2′′:5′′,2′′′-quaterthiophene (DH4T) semiconductor and the n-type 2,2′′′- Perfluoro-Dihexyl-2,2′:5′,2′′:5′′,2′′′-quaterthiophene (DFH4T) semiconductor underwent several sublimation cycles, with consequent improvement of the electrical performances in terms of charge mobility and threshold voltage, highlighting the benefits brought by this treatment to the electric properties of the discussed semiconductors in OFET devices by the removal of residual impurities. The second step consisted in the provision of a metal-free synthesis of DH4T, which was successfully prepared without organometallic reagents or catalysts in collaboration with Dr. Manuela Melucci from ISOF-CNR Institute in Bologna. Indeed the experimental work demonstrated that those compounds are responsible for the electrical degradation by intentionally doping the semiconductor obtained by metal-free method by Tetrakis(triphenylphosphine)palladium(0) (Pd(PPh3)4) and Tributyltin chloride (Bu3SnCl), as well as with an organic impurity, like 5-hexyl-2,2':5',2''-terthiophene (HexT3) at, in different concentrations (1, 5 and 10% w/w). After completing the entire evaluation process loop, from fabricating OFET devices by vacuum sublimation with implemented intentionally-doped batches to the final electrical characterization in inherent-atmosphere conditions, commercial DH4T, metal-free DH4T and the intentionally-doped DH4T were systematically compared. Indeed, the fabrication of OFET based on doped DH4T clearly pointed out that the vacuum sublimation is still an inherent and efficient purification method for crude semiconductors, but also a reliable way to fabricate high performing devices.

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.

Hydraulic evaluation of the gyro electric wave energy converter

Motivation Thanks for a scholarship offered by ALma Mater Studiorum I could stay in Denmark for six months during which I could do physical tests on the device Gyro PTO at the Departmet of Civil Engineering of Aalborg University. Aim The goal of my thesis is an hydraulic evaluation of the device: Gyro PTO, a gyroscopic device for conversion of mechanical energy in ocean surface waves to electrical energy. The principle of the system is the application of the gyroscopic moment of flywheels equipped on a swing float excited by waves. The laboratory activities were carried out by: Morten Kramer, Jan Olsen, Irene Guaraldi, Morten Thøtt, Nikolaj Holk. The main purpose of the tests was to investigate the power absorption performance in irregular waves, but testing also included performance measures in regular waves and simple tests to get knowledge about characteristics of the device, which could facilitate the possibility of performing numerical simulations and optimizations. Methodology To generate the waves and measure the performance of the device a workstation was created in the laboratory. The workstation consist of four computers in each of wich there was a different program. Programs have been used : Awasys6, LabView, Wave lab, Motive optitrack, Matlab, Autocad Main Results Thanks to the obtained data with the tank testing was possible to make the process of wave analisys. We obtained significant wave height and period through a script Matlab and then the values of power produced, and energy efficiency of the device for two types of waves: regular and irregular. We also got results as: physical size, weight, inertia moments, hydrostatics, eigen periods, mooring stiffness, friction, hydrodynamic coefficients etc. We obtained significant parameters related to the prototype in the laboratory after which we scale up the results obtained for two future applications: one in Nissun Brending and in the North Sea. Conclusions The main conclusion on the testing is that more focus should be put into ensuring a stable and positive power output in a variety of wave conditions. In the irregular waves the power production was negative and therefore it does not make sense to scale up the results directly. The average measured capture width in the regular waves was 0.21 m. As the device width is 0.63 m this corresponds to a capture width ratio of: 0.21/0.63 * 100 = 33 %. Let’s assume that it is possible to get the device to produce as well in irregular waves under any wave conditions, and lets further assume that the yearly absorbed energy can be converted into electricity at a PTO-efficiency of 90 %. Under all those assumptions the results in table are found, i.e. a Nissum Bredning would produce 0.87 MWh/year and a North Sea device 85 MWh/year.

Galactic fountain flows in high-resolution hydrodynamical simulations of galaxy disks

Feedback from the most massive components of a young stellar cluster deeply affects the surrounding ISM driving an expanding over-pressured hot gas cavity in it. In spiral galaxies these structures may have sufficient energy to break the disk and eject large amount of material into the halo. The cycling of this gas, which eventually will fall back onto the disk, is known as galactic fountains. We aim at better understanding the dynamics of such fountain flow in a Galactic context, frame the problem in a more dynamic environment possibly learning about its connection and regulation to the local driving mechanism and understand its role as a metal diffusion channel. The interaction of the fountain with a hot corona is hereby analyzed, trying to understand the properties and evolution of the extraplanar material. We perform high resolution hydrodynamical simulations with the moving-mesh code AREPO to model the multi-phase ISM of a Milky Way type galaxy. A non-equilibrium chemical network is included to self consistently follow the evolution of the main coolants of the ISM. Spiral arm perturbations in the potential are considered so that large molecular gas structures are able to dynamically form here, self shielded from the interstellar radiation field. We model the effect of SN feedback from a new-born stellar cluster inside such a giant molecular cloud, as the driving force of the fountain. Passive Lagrangian tracer particles are used in conjunction to the SN energy deposition to model and study diffusion of freshly synthesized metals. We find that both interactions with hot coronal gas and local ISM properties and motions are equally important in shaping the fountain. We notice a bimodal morphology where most of the ejected gas is in a cold $10^4$ K clumpy state while the majority of the affected volume is occupied by a hot diffuse medium. While only about 20\% of the produced metals stay local, most of them quickly diffuse through this hot regime to great scales.

Numerical simulations of the effects of residual stresses in adhesively bonded composite specimen

In this work the problem of performing a numerical simulation of quasi-static crack propagation within an adhesive layer of a bonded joint under Mode I loading affected by stress field changes due to thermal-chemical shrinkage induced by cure process is addressed. Secondly, a parametric study on fracture critical energy, cohesive strength and Young's modulus is performed. Finally, a particular case of adhesive layer stiffening is simulated in order to verify qualitatively the major effect.

From the chiral magnetic wave to the charge dependence of elliptic flow

The quark-gluon plasma formed in heavy ion collisions contains charged chiral fermions evolving in an external magnetic field. At finite density of electric charge or baryon number (resulting either from nuclear stopping or from fluctuations), the triangle anomaly induces in the plasma the Chiral Magnetic Wave (CMW). The CMW first induces a separation of the right and left chiral charges along the magnetic field; the resulting dipolar axial charge density in turn induces the oppositely directed vector charge currents leading to an electric quadrupole moment of the quark-gluon plasma. Boosted by the strong collective flow, the electric quadrupole moment translates into the charge dependence of the elliptic flow coefficients, so that $v_2(\pi^+) < v_2(\pi^-)$ (at positive net charge). Using the latest quantitative simulations of the produced magnetic field and solving the CMW equation, we make further quantitative estimates of the produced $v_2$ splitting and its centrality dependence. We compare the results with the available experimental data.

Magnetic field exposure and neurodegenerative diseases--recent epidemiological studies

OBJECTIVES: To analyse the results of recent studies not yet included in a 2003 report of the International Commission on Non-Ionizing Radiation Protection (ICNIRP) on occupational exposure to low-frequency electromagnetic fields as potential risk factor for neurodegenerative diseases. METHODS: A literature search was conducted in the online databases of PubMed, ISI Web of Knowledge, DIMDI and COCHRANE, as well as in specialised databases and journals. Eight studies published between January 2000 and July 2005 were included in the review. RESULTS: The findings of these studies contribute to the evidence of an association between occupational magnetic field exposure and the risk of dementia. Regarding amyotrophic lateral sclerosis, the recent results confirm earlier observations of an association with electric and electronic work and welding. Its relationship with magnetic field exposure remains unsolved. There are only few findings pointing towards an association between magnetic field exposure and Parkinson's disease. CONCLUSIONS: The epidemiological evidence for an association between occupational exposure to low-frequency electromagnetic fields and the risk of dementia has increased during the last five years. The impact of potential confounders should be evaluated in further studies.

Collapsed cone convolution and analytical anisotropic algorithm dose calculations compared to VMC++ Monte Carlo simulations in clinical cases

The purpose of this work was to study and quantify the differences in dose distributions computed with some of the newest dose calculation algorithms available in commercial planning systems. The study was done for clinical cases originally calculated with pencil beam convolution (PBC) where large density inhomogeneities were present. Three other dose algorithms were used: a pencil beam like algorithm, the anisotropic analytic algorithm (AAA), a convolution superposition algorithm, collapsed cone convolution (CCC), and a Monte Carlo program, voxel Monte Carlo (VMC++). The dose calculation algorithms were compared under static field irradiations at 6 MV and 15 MV using multileaf collimators and hard wedges where necessary. Five clinical cases were studied: three lung and two breast cases. We found that, in terms of accuracy, the CCC algorithm performed better overall than AAA compared to VMC++, but AAA remains an attractive option for routine use in the clinic due to its short computation times. Dose differences between the different algorithms and VMC++ for the median value of the planning target volume (PTV) were typically 0.4% (range: 0.0 to 1.4%) in the lung and -1.3% (range: -2.1 to -0.6%) in the breast for the few cases we analysed. As expected, PTV coverage and dose homogeneity turned out to be more critical in the lung than in the breast cases with respect to the accuracy of the dose calculation. This was observed in the dose volume histograms obtained from the Monte Carlo simulations.