2D and 3D numerical modelling of multilayer detachment folding and salt tectonics

Numerical modelling was performed to study the dynamics of multilayer detachment folding and salt tectonics. In the case of multilayer detachment folding, analytically derived diagrams show several folding modes, half of which are applicable to crustal scale folding. 3D numerical simulations are in agreement with 2D predictions, yet fold interactions result in complex fold patterns. Pre-existing salt diapirs change folding patterns as they localize the initial deformation. If diapir spacing is much smaller than the dominant folding wavelength, diapirs appear in fold synclines or limbs.rnNumerical models of 3D down-building diapirism show that sedimentation rate controls whether diapirs will form and influences the overall patterns of diapirism. Numerical codes were used to retrodeform modelled salt diapirs. Reverse modelling can retrieve the initial geometries of a 2D Rayleigh-Taylor instability with non-linear rheologies. Although intermediate geometries of down-built diapirs are retrieved, forward and reverse modelling solutions deviate. rnFinally, the dynamics of fold-and-thrusts belts formed over a tilted viscous detachment is studied and it is demonstrated that mechanical stratigraphy has an impact on the deformation style, switching from thrust- to folding-dominated. The basal angle of the detachment controls the deformation sequence of the fold-and-thrust belt and results are consistent with critical wedge theory.rn

Nucleation studies on graphics processing units

Ein System in einem metastabilen Zustand muss eine bestimmte Barriere in derrnfreien Energie überwinden um einen Tropfen der stabilen Phase zu formen.rnHerkömmliche Untersuchungen nehmen hierbei kugelförmige Tropfen an. Inrnanisotropen Systemen (wie z.B. Kristallen) ist diese Annahme aber nicht ange-rnbracht. Bei tiefen Temperaturen wirkt sich die Anisotropie des Systems starkrnauf die freie Energie ihrer Oberfläche aus. Diese Wirkung wird oberhalb derrnAufrauungstemperatur T R schwächer. Das Ising-Modell ist ein einfaches Mo-rndell, welches eine solche Anisotropie aufweist. Wir führen großangelegte Sim-rnulationen durch, um die Effekte, die mit einer endlichen Simulationsbox ein-rnhergehen, sowie statistische Ungenauigkeiten möglichst klein zu halten. DasrnAusmaß der Simulationen die benötigt werden um sinnvolle Ergebnisse zu pro-rnduzieren, erfordert die Entwicklung eines skalierbaren Simulationsprogrammsrnfür das Ising-Modell, welcher auf verschiedenen parallelen Architekturen (z.B.rnGrafikkarten) verwendet werden kann. Plattformunabhängigkeit wird durch ab-rnstrakte Schnittstellen erreicht, welche plattformspezifische Implementierungs-rndetails verstecken. Wir benutzen eine Systemgeometrie die es erlaubt eine Ober-rnfläche mit einem variablen Winkel zur Kristallebene zu untersuchen. Die Ober-rnfläche ist in Kontakt mit einer harten Wand, wobei der Kontaktwinkel Θ durchrnein Oberflächenfeld eingestellt werden kann. Wir leiten eine Differenzialglei-rnchung ab, welche das Verhalten der freien Energie der Oberfläche in einemrnanisotropen System beschreibt. Kombiniert mit thermodynamischer Integrationrnkann die Gleichung benutzt werden, um die anisotrope Oberflächenspannungrnüber einen großen Winkelbereich zu integrieren. Vergleiche mit früheren Mes-rnsungen in anderen Geometrien und anderen Methoden zeigen hohe Überein-rnstimung und Genauigkeit, welche vor allem durch die im Vergleich zu früherenrnMessungen wesentlich größeren Simulationsdomänen erreicht wird. Die Temper-rnaturabhängigkeit der Oberflächensteifheit κ wird oberhalb von T R durch diernKrümmung der freien Energie der Oberfläche für kleine Winkel gemessen. DiesernMessung lässt sich mit Simulationsergebnissen in der Literatur vergleichen undrnhat bessere Übereinstimmung mit theoretischen Voraussagen über das Skalen-rnverhalten von κ. Darüber hinaus entwickeln wir ein Tieftemperatur-Modell fürrndas Verhalten um Θ = 90 Grad weit unterhalb von T R. Der Winkel bleibt bis zu einemrnkritischen Feld H C quasi null; oberhalb des kritischen Feldes steigt der Winkelrnrapide an. H C wird mit der freien Energie einer Stufe in Verbindung gebracht,rnwas es ermöglicht, das kritische Verhalten dieser Größe zu analysieren. Die harternWand muss in die Analyse einbezogen werden. Durch den Vergleich freier En-rnergien bei geschickt gewählten Systemgrößen ist es möglich, den Beitrag derrnKontaktlinie zur freien Energie in Abhängigkeit von Θ zu messen. Diese Anal-rnyse wird bei verschiedenen Temperaturen durchgeführt. Im letzten Kapitel wirdrneine 2D Fluiddynamik Simulation für Grafikkarten parallelisiert, welche u. a.rnbenutzt werden kann um die Dynamik der Atmosphäre zu simulieren. Wir im-rnplementieren einen parallelen Evolution Galerkin Operator und erreichen

Design methods and geometric modelling of scaffolds for bone tissue engineering

Every year, thousand of surgical treatments are performed in order to fix up or completely substitute, where possible, organs or tissues affected by degenerative diseases. Patients with these kind of illnesses stay long times waiting for a donor that could replace, in a short time, the damaged organ or the tissue. The lack of biological alternates, related to conventional surgical treatments as autografts, allografts, e xenografts, led the researchers belonging to different areas to collaborate to find out innovative solutions. This research brought to a new discipline able to merge molecular biology, biomaterial, engineering, biomechanics and, recently, design and architecture knowledges. This discipline is named Tissue Engineering (TE) and it represents a step forward towards the substitutive or regenerative medicine. One of the major challenge of the TE is to design and develop, using a biomimetic approach, an artificial 3D anatomy scaffold, suitable for cells adhesion that are able to proliferate and differentiate themselves as consequence of the biological and biophysical stimulus offered by the specific tissue to be replaced. Nowadays, powerful instruments allow to perform analysis day by day more accurateand defined on patients that need more precise diagnosis and treatments.Starting from patient specific information provided by TC (Computed Tomography) microCT and MRI(Magnetic Resonance Imaging), an image-based approach can be performed in order to reconstruct the site to be replaced. With the aid of the recent Additive Manufacturing techniques that allow to print tridimensional objects with sub millimetric precision, it is now possible to practice an almost complete control of the parametrical characteristics of the scaffold: this is the way to achieve a correct cellular regeneration. In this work, we focalize the attention on a branch of TE known as Bone TE, whose the bone is main subject. Bone TE combines osteoconductive and morphological aspects of the scaffold, whose main properties are pore diameter, structure porosity and interconnectivity. The realization of the ideal values of these parameters represents the main goal of this work: here we'll a create simple and interactive biomimetic design process based on 3D CAD modeling and generative algorithmsthat provide a way to control the main properties and to create a structure morphologically similar to the cancellous bone. Two different typologies of scaffold will be compared: the first is based on Triply Periodic MinimalSurface (T.P.M.S.) whose basic crystalline geometries are nowadays used for Bone TE scaffolding; the second is based on using Voronoi's diagrams and they are more often used in the design of decorations and jewellery for their capacity to decompose and tasselate a volumetric space using an heterogeneous spatial distribution (often frequent in nature). In this work, we will show how to manipulate the main properties (pore diameter, structure porosity and interconnectivity) of the design TE oriented scaffolding using the implementation of generative algorithms: "bringing back the nature to the nature".

Modellkomplexe zur Sauerstoffaktivierung auf Basis des 1,3,4-Thiadiazols

Das Ziel der vorliegenden Arbeit waren die Synthese und Untersuchung von Modellverbindungen zur Sauerstoffaktivierung auf der Basis neuer Ligandensysteme des 1,3,4-Thiadiazols unter Ausarbeitung einer Synthesestrategie zur Derivatisierung der heteroaromatischen 1,3,4-Thiadiazol-Liganden, deren Koordinationsverhalten in Abhängigkeit ihres 2,5-Substitutionsmusters untersucht wurde, sowie die fortführende Bearbeitung bereits bekannter Ligandensysteme zur Erzeugung von homo- und heterovalenten Übergangsmetallkomplexverbindungen.rnDie unter der Verwendung der modifizierten Liganden TPDE, H1TPDP und H1BPMP resultierenden dinuklearen Komplexverbindungen zeigen unterschiedlich starke antiferromagnetische Wechselwirkungen in Abhängigkeit der vorhandenen Brückenliganden. In der Verbindung [Fe6O2(OH)(L´)2(OOCMe3)9(OEt)2] trat eine Fragmentierung des Liganden H1TPDP auf. Das cisoide Ligandensubstitutionsmuster der entstandenen sechskernigen Verbindung ist verantwortlich für die interessanten magnetischen Eigenschaften des Komplexes. rnNeue Perspektiven zur Erzeugung von Modellverbindungen zur Sauerstoffaktivierung wurden mit dem Mono-Chelatliganden H1ETHP und den Bis-Chelatliganden HL2H, H2L2H und H2BATP aufgezeigt. Die Umsetzung von H1ETHP mit verschiedenen Übergangsmetallsalzen resultierte für die Metalle Cr(III), Fe(III), Co(III) und Ni(II) in mononuklearen Verbindungen des Typs [M(ETHP)2]X (X = ClO4, FeCl4, OMe, Cl, Br) sowie in zwei tetranuklearen Verbindungen mit Mn(II) und Cu(II). [Mn4(ETHP)6] besitzt ein propellerförmiges, planares [Mn4O6]2+-System mit einen Spingrundzustand von S = 5. In allen Verbindungen von H1ETHP konnte eine mono-κN-Koordination des 1,3,4-Thiadiazol-Rückgrates über eines seiner beiden endozyklischen Stickstoffdonoratome beobachtet werden. rnAus Umsetzungen der Bis-Chelatliganden wurden fast ausschließlich polynukleare Übergangsmetallkomplexe erhalten. Insbesondere der Ligand H2L2H zeigt eine ausgeprägte Tendenz zur Ausbildung trinuklearer, linearer Komplexe, welche auf Grund ihrer ungeraden Anzahl von Übergangsmetallionen einen Spingrundzustand S ≠ 0 aufweisen.rn Die mit dem Liganden HL2H erhaltenen Verbindungen unterstreichen die hohe Flexibilität dieser Systeme hinsichtlich der Erzeugung polynuklearer und heterovalenter Komplexverbindungen. So konnten in Abhängigkeit vom verwendeten Übergangsmetallsalz trinukleare, pentanukleare, aber auch hepta- und oktanukleare Verbindungen synthetisiert werden. Insbesondere die Komplexe des Mangans und des Cobalts zeigen ein heterovalentes [MnIIMnIII4]- bzw. [CoII2CoIII3]-Motiv, was sich in Spingrundzuständen von S ≠ 0 äußert. Der diamagnetische, achtkernige Fe8-Cluster besitzt eine pseudo C3-symmetrische Anordnung der Metall-Zentren, während für die heptanukleare Cu7-Kette durch ihre stark unterschiedlichen Kupfer-Koordinationsgeometrien interessante magnetische Austauschwechselwirkungen beobachtet werden konnten. Der dreikernige µ3-oxo-verbrückte Komplex des Liganden H2BATP zeigt als interessante strukturelle Eigenschaft ein ein µ3-Verbrückungsmuster des eingesetzten Sulfat-Anions. rnIn allen Komplexen der Bis-Chelatliganden HL2H, H2L2H und H2BATP konnte ein µ2-κN,κN-Koordiantionsmodus des 1,3,4-Thiadiazols und somit eine Abhängigkeit der Verbrückung vom Ligandensubstitutionsmuster beobachtet werden.rn

Field-free control of magnetism in nanostructured materials probed with high resolution X-ray microscopy

Magnetic memories are a backbone of today's digital data storage technology, where the digital information is stored as the magnetic configuration of nanostructured ferromagnetic bits. Currently, the writing of the digital information on the magnetic memory is carried out with the help of magnetic fields. This approach, while viable, is not optimal due to its intrinsically high energy consumption and relatively poor scalability. For this reason, the research for different mechanisms that can be used to manipulate the magnetic configuration of a material is of interest. In this thesis, the control of the magnetization of different nanostructured materials with field-free mechanisms is investigated. The magnetic configuration of these nanostructured materials was imaged directly with high resolution x-ray magnetic microscopy. rnFirst of all, the control of the magnetic configuration of nanostructured ferromagnetic Heusler compounds by fabricating nanostructures with different geometries was analyzed. Here, it was observed that the magnetic configuration of the nanostructured elements is given by the competition of magneto-crystalline and shape anisotropy. By fabricating elements with different geometries, we could alter the point where these two effects equilibrate, allowing for the possibility to tailor the magnetic configuration of these nanostructured elements to the required necessities.rnThen, the control of the magnetic configuration of Ni nanostructures fabricated on top of a piezoelectric material with the magneto-elastic effect (i.e. by applying a piezoelectric strain to the Ni nanostructures) was investigated. Here, the magneto-elastic coupling effect gives rise to an additional anisotropy contribution, proportional to the strain applied to the magnetic material. For this system, a reproducible and reversible control of the magnetic configuration of the nanostructured Ni elements with the application of an electric field across the piezoelectric material was achieved.rnFinally, the control of the magnetic configuration of La0.7Sr0.3MnO3 (LSMO) nanostructures with spin-polarized currents was studied. Here, the spin-transfer torque effect was employed to achieve the displacement of magnetic domain walls in the LSMO nanostructures. A high spin-transfer torque efficiency was observed for LSMO at low temperatures, and a Joule-heating induced hopping of the magnetic domain walls was observed at room temperatures, allowing for the analysis of the energetics of the domain walls in LSMO.rnThe results presented in this thesis give thus an overview on the different field-free approaches that can be used to manipulate and tailor the magnetization configuration of a nanostructured material to the various technological requirements, opening up novel interesting possibilities for these materials.

Novel organic semiconductors and their applications in electronics

The main goals of this thesis were the design, synthesis, and characterization of novel organic semiconductors, together with their applications in electronics, such as OFETs, OPVs, and OLEDs. The results can be summarized as follows:rn1. In chapter II, two novel angular n-type molecules were presented. Their different alkyl chains play a pivotal role in the molecular orientation relative to surface. One molecule with longer branched chains is tilted with respect to the substrate, thereby resulting in poor device performance, while the other adopt an edge-on orientation with an OFET electron mobility of 0.01 cm2 V-1 s-1.rn2. In chapter III, fused bis-benzothiadiazoles with different molecular geometries, namely linear benzoquinone-fused bis(benzothiadiazole) and V-shaped sulfone-fused bis(benzothiadiazole), were shown. This work not only contributes to the diversity of electron acceptors based on bis-benzothiadiazole moieties, but also highlights the important role of molecular shape for the solid-state packing of organic conjugated materials. In chapter IV, we demonstrated the synthesis of layered acceptors via dimerization of thiadiazole end-capped acenes. Interestingly, they feature huge differences in their photophysical properties. One compound showed a new strong emission in the near-infrared region introduced by the aggregation effect. The planosymmetric compound featured intramolecular excimer (IEE) fluorescence in solution. rn3. In chapter V and VI, we have demonstrated the synthesis of novel spiro-bifluorene based asymmetric and symmetric cruciform electron acceptors with dicyanovinylene substitutions. The solar cells based on PTB7:asymmetric acceptor yields the highest PCE of 0.80%. Such results demonstrate for the first time that dicyanovinylene substituted acceptor could be an alternative to fullerene-based acceptors. rn4. In chapter VII, two novel blue-emitting compounds were shown, which consist of dihydroindenofluorenyl units and ladder-type poly-p-phenylene groups, respectively. The two novel cruciform rigid compounds present not only excellent thermal and electrochemical stability but also high PLQYs. Through analysis of their triplet energy levels, both molecules can be served as hosts for other normal fluorescent or phosphorescent materials.rn

Dynamic characterisation of four nine-story large-panel R.C. buildings in Bishkek (Kyrgyzstan): A comparison between experimental ambient vibration analysis and numerical finite element modeling

With the outlook of improving seismic vulnerability assessment for the city of Bishkek (Kyrgyzstan), the global dynamic behaviour of four nine-storey r.c. large-panel buildings in elastic regime is studied. The four buildings were built during the Soviet era within a serial production system. Since they all belong to the same series, they have very similar geometries both in plan and in height. Firstly, ambient vibration measurements are performed in the four buildings. The data analysis composed of discrete Fourier transform, modal analysis (frequency domain decomposition) and deconvolution interferometry, yields the modal characteristics and an estimate of the linear impulse response function for the structures of the four buildings. Then, finite element models are set up for all four buildings and the results of the numerical modal analysis are compared with the experimental ones. The numerical models are finally calibrated considering the first three global modes and their results match the experimental ones with an error of less then 20%.

Conformational transitions of glycine induced by vibrational excitation of the O-H stretch

Vibrational energy flow and conformational transitions following excitation of the OH stretching mode of the most stable conformer of glycine are studied by classical trajectories. "On the fly" simulations with the PM3 semiempirical electronic structure method for the potential surface are used. Initial conditions are selected to correspond to the v = 1 excitation of the OH stretch. The main findings are: (1) An an equilibrium-like ratio is established between the populations of the 3 lowest-lying conformers after about 10 picoseconds. (2) There is a high probability throughout the 150 ps of the simulations for finding the molecule in geometries far from the equilibrium structures of the lowest-energy conformers. (3) Energy from the initial excited OH (v = 1) stretch flows preferentially to 5 other vibrational modes, including the bending motion of the H atom. (4) RRK theory yields conformational transition rates that deviate substantially from the classical trajectory results. Possible implication of these results for vibrational energy flow and conformational transitions in small biological molecules are discussed.

Monte Carlo Dose Calculation on Deforming Anatomy

This article presents the implementation and validation of a dose calculation approach for deforming anatomical objects. Deformation is represented by deformation vector fields leading to deformed voxel grids representing the different deformation scenarios. Particle transport in the resulting deformed voxels is handled through the approximation of voxel surfaces by triangles in the geometry implementation of the Swiss Monte Carlo Plan framework. The focus lies on the validation methodology which uses computational phantoms representing the same physical object through regular and irregular voxel grids. These phantoms are chosen such that the new implementation for a deformed voxel grid can be compared directly with an established dose calculation algorithm for regular grids. Furthermore, separate validation of the aspects voxel geometry and the density changes resulting from deformation is achieved through suitable design of the validation phantom. We show that equivalent results are obtained with the proposed method and that no statistically significant errors are introduced through the implementation for irregular voxel geometries. This enables the use of the presented and validated implementation for further investigations of dose calculation on deforming anatomy.

Further Quantum Mechanical Evidence that Difluorotoluene does not Hydrogen Bond

Calculations were run on the methylated DNA base pairs adenine:thymine and adenine:difluorotoluene to further investigate the hydrogen-bonding properties of difluorotoluene (F). Geometries were optimized using hybrid density functional theory. Single-point calculations at the MP2(full) level were performed to obtain more rigorous energies. The functional counterpoise method was used to correct for the basis set superposition error (BSSE), and the interaction energies were also corrected for fragment relaxation. These corrections brought the B3LYP and MP2 interaction energies into excellent agreement. In the gas phase, the Gibbs free energies calculated at the B3LYP and MP2 levels of theory predict that A and T will spontaneously form an A:T pair while A:F spontaneously dissociates into A and F. Solvation effects on the pairing of the bases were explored using implicit solvent models for water and chloroform. In aqueous solution, both A:T and A:F are predicted to dissociate into their component monomers. Semiempirical calculations were performed on small sections of B-form DNA containing the two pairs, and the results provide support for the concept that base stacking is more important than hydrogen bonding for the stability of the A:F pair within a DNA helix.

Semiempirical Study of the Bergman Reaction: Towards a Computationally Efficient and Accurate Method for Modeling Enediyne Anticancer Antibiotics

Over the past 7 years, the enediyne anticancer antibiotics have been widely studied due to their DNA cleaving ability. The focus of these antibiotics, represented by kedarcidin chromophore, neocarzinostatin chromophore, calicheamicin, esperamicin A, and dynemicin A, is on the enediyne moiety contained within each of these antibiotics. In its inactive form, the moiety is benign to its environment. Upon suitable activation, the system undergoes a Bergman cycloaromatization proceeding through a 1,4-dehydrobenzene diradical intermediate. It is this diradical intermediate that is thought to cleave double-stranded dna through hydrogen atom abstraction. Semiempirical, semiempiricalci, Hartree–Fock ab initio, and mp2 electron correlation methods have been used to investigate the inactive hex-3-ene-1,5-diyne reactant, the 1,4-dehydrobenzene diradical, and a transition state structure of the Bergman reaction. Geometries calculated with different basis sets and by semiempirical methods have been used for single-point calculations using electron correlation methods. These results are compared with the best experimental and theoretical results reported in the literature. Implications of these results for computational studies of the enediyne anticancer antibiotics are discussed.

In search of CS2(H2O)n = 1–4 clusters

Gaussian-3 and MP2/aug-cc-pVnZ methods have been used to calculate geometries and thermochemistry of CS2(H2O)n, where n = 1–4. An extensive molecular dynamics search followed by optimization using these two methods located two dimers, six trimers, six tetramers, and two pentamers. The MP2/aug-cc-pVDZ structure matched best with the experimental result for the CS2(H2O) dimer, showing that diffuse functions are necessary to model the interactions found in this complex. For larger CS2(H2O)n clusters, the MP2/aug-cc-pVDZ minima are significantly different from the MP2(full)/6-31G* structures, revealing that the G3 model chemistry is not suitable for investigation of sulfur containing van der Waals complexes. Based on the MP2/aug-cc-pVTZ free energies, the concentration of saturated water in the atmosphere and the average amount of CS2 in the atmosphere, the concentrations of these clusters are predicted to be on the order of 105CS2(H2O) clusters∙cm−3 and 102 CS2(H2O)2 clusters∙cm−3 at 298.15 K. The MP2/aug-cc-pVDZ scaled harmonic and anharmonic frequencies of the most abundant dimer cluster at 298 K are presented, along with the MP2/aug-cc-pVDZ scaled harmonic frequencies for the CS2(H2O)n structures predicted to be present in a low-temperature molecular beam experiment.

Ortho Effect in the Bergman Cyclization: Electronic and Steric Effects in Hydrogen Abstraction by 1-Substituted Naphthalene 5,8-Diradicals

We present a detailed theoretical study of geometries, electronic structure, and energies of transition states and intermediates completing the full Bergman cycloaromatization pathway of ortho-substituted enediynes with a focus on polar and steric contributions to the kinetics and thermodynamics of hydrogen abstraction. This study provides a rare unambiguous example of remote substitution that affects reactivity of a neutral reactive intermediate through an σ framework.

A Computationally Efficient Procedure for Modeling the First Step in the Alkaline Hydrolysis of Esters

A computationally efficient procedure for modeling the alkaline hydrolysis of esters is proposed based on calculations performed on methyl acetate and methyl benzoate systems. Extensive geometry and energy comparisons were performed on the simple ester methyl acetate. The effectiveness of performing high level single point ab initio energy calculations on the geometries obtained from semiempirical and ab initio methods was determined. The AM1 and PM3 semiempirical methods are evaluated for their ability to model the transition states and intermediates for ester hydrolysis. The Cramer/Truhlar SM3 solvation method was used to determine activation energies. The most computationally efficient way to model the transition states of large esters is to use the PM3 method. The PM3 transition structure can then be used as a template for the design of haptens capable of inducing catalytic antibodies.

CCSD(T), W1, and other model chemistry predictions for gas-phase deprotonation reactions

A series of CCSD(T) single-point calculations on MP4(SDQ) geometries and the W1 model chemistry method have been used to calculate ΔH° and ΔG° values for the deprotonation of 17 gas-phase reactions where the experimental values have reported accuracies within 1 kcal/mol. These values have been compared with previous calculations using the G3 and CBS model chemistries and two DFT methods. The most accurate CCSD(T) method uses the aug-cc-pVQZ basis set. Extrapolation of the aug-cc-pVTZ and aug-cc-pVQZ results yields the most accurate agreement with experiment, with a standard deviation of 0.58 kcal/mol for ΔG° and 0.70 kcal/mol for ΔH°. Standard deviations from experiment for ΔG° and ΔH° for the W1 method are 0.95 and 0.83 kcal/mol, respectively. The G3 and CBS-APNO results are competitive with W1 and are much less expensive. Any of the model chemistry methods or the CCSD(T)/aug-cc-pVQZ method can serve as a valuable check on the accuracy of experimental data reported in the National Institutes of Standards and Technology (NIST) database.