Numerical modelling of fluid-structure interaction with application in hemodynamics

The thesis deals with numerical algorithms for fluid-structure interaction problems with application in blood flow modelling. It starts with a short introduction on the mathematical description of incompressible viscous flow with non-Newtonian viscosity and a moving linear viscoelastic structure. The mathematical model consists of the generalized Navier-Stokes equation used for the description of fluid flow and the generalized string model for structure movement. The arbitrary Lagrangian-Eulerian approach is used in order to take into account moving computational domain. A part of the thesis is devoted to the discussion on the non-Newtonian behaviour of shear-thinning fluids, which is in our case blood, and derivation of two non-Newtonian models frequently used in the blood flow modelling. Further we give a brief overview on recent fluid-structure interaction schemes with discussion about the difficulties arising in numerical modelling of blood flow. Our main contribution lies in numerical and experimental study of a new loosely-coupled partitioned scheme called the kinematic splitting fluid-structure interaction algorithm. We present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Here, we assume both, the nonlinearity in convective as well is diffusive term. We analyse the convergence of proposed numerical scheme for a simplified fluid model of the Oseen type. Moreover, we present series of experiments including numerical error analysis, comparison of hemodynamic parameters for the Newtonian and non-Newtonian fluids and comparison of several physiologically relevant computational geometries in terms of wall displacement and wall shear stress. Numerical analysis and extensive experimental study for several standard geometries confirm reliability and accuracy of the proposed kinematic splitting scheme in order to approximate fluid-structure interaction problems.

Systematic studies of the interaction between amyloid beta-protein and lipids

The amyloid peptide (Aß), a normal constituent of neuronal and non-neuronal cells, has been shown to be a major component of the extracellular plaque of Alzheimer’s disease (AD). The interaction of Aß peptides with the lipid matrix of neuronal cell membranes plays an important role in the pathogenesis of AD. In this study, we have developed peptide-tethered artificial lipid membranes by the Langmuir-Blodgett and Langmuir-Schaefer methods. Anti-Aß40-mAb labeled with a fluorophore was used to probe the Aß40 binding to the model membrane system. Systematic studies on the antibody or Aß-membrane interactions were carried out in our model systems by Surface Plasmon Field-Enhanced Fluorescence Spectroscopy (SPFS). Aß adsorption is critically determined by the lipid composition of the membranes. Aß specifically binds with membranes of sphingomyelin, and this preferential adsorption was markedly amplified by the addition of sterols (cholesterol or 25-OH-Chol). Fluorescence microscopy indicated that 25-OH-Chol could also form micro-domains with sphingomyelin as cholesterol does at the conditions used for the built-up of the model membranes. Our findings suggest that micro-domains composed of sphingomyelin and the sterols could be the binding sites of Aß and the role of sphingomyelin in AD should receive much more attention. The artificial membranes provide a novel platform for the study on AD, and SPFS is a potential tool for detecting Aß-membrane interaction. Numerous investigations indicate that the ability of Aß to form fibrils is considerably dependent upon the levels of ß-sheet structure adopted by Aß. Membrane-mediated conformational transition of Aß has been demonstrated. In this study, we focus on the interaction of Aß and the membranes composed of POPC/SM/25-OH-Chol (2:1:1). The artificial membrane system was established by the methods as described above. Immunoassy based on a pair of monoclonal antibodies (mAbs) against different epitopes was employed to detect the orientation of the Aß at the model membranes. Kinetics of antibody-Aß binding was determined by surface plasmon field-enhanced fluorescence spectroscopy (SPFS). The attempt has also been made to probe the change in the conformation of Aß using SPFS combined with immunoassay. Melatonin was employed to induce the conformational change of Aß. The orientation and the conformational change of Aß are evaluated by analysing kinetic/affinity parameters. This work provides novel insight into the investigation on the structure of Aß at the membrane surface.

Stress – structure correlations in grafted polymer films

In der vorliegenden Arbeit wurden experimentelle Untersuchungen zu gepfropften Polymerfilmen durchgeführt. Dabei wurden endgepfropfte poly-methyl-methacrylate (PMMA) Bürsten hergestellt durch „grafting from“ Methoden und polystyrol (PS)/ poly-vinyl-methyl-ether (PVME) Polymerfilme gepfropft auf UV sensitiven Oberflächen untersucht. Zur Strukturuntersuchung wurden die hergestellten Systeme wurden mit Rasterkraftmikroskopie (engl.: Surface Probe Microscopy, SPM), Röntgen - und Neutronenreflektivitätsmessungen, sowie mit Röntgenstreuung unter streifenden Einfall (engl.: Grazing Incidence Small Angle X-Ray Scattering, GISAXS) untersucht. rnEs wurde gezeigt, dass ein aus der Transmissionsstreuung bekanntes Model auch für auch für die GISAXS Analyse polydisperser Polymerdomänen und Kolloidsysteme verwendet werden kann. Der maximale Fehler durch die gemachten Näherungen wurde auf < 20% abgeschätzt.rnErgebnisse aus der Strukturanalyse wurden mit mechanischen Filmeigenschaften verknüpft. Dazu wurden mechanische Spannungsexperimente durchgeführt. Hierzu wurden die zu untersuchenden Filme selektiv auf einzelne Mikro-Federbalken-Sensoren (engl.: Micro Cantilever Sensor, MCS) der MCS Arrays aufgebracht. Dies wurde durch Maskierungstechniken und Mikro-Kontaktdrucken bewerkstelligt. rnPhasenübergansexperimente der gepfropften PS/PVME Filme haben gezeigt, dass die Möglichkeit einer Polymer/Polymer Phasenseparation stark von Propfpunktdichte der gebundenen Polymerketten mit der Oberfläche abhängt. PS/PVME Filmsysteme mit hohen Pfropfpunktdichten zeigten keinen Phasenübergang. Bei niedrig gepfropften Filmsystemen waren hingegen Polymer/Polymer Phasenseparationen zu beobachten. Es wurde geschlussfolgert, dass die gepfropften Polymersysteme einen hinreichenden Grad an entropischen Freiheitsgraden benötigen um eine Phasenseparation zu zeigen. Mechanische Spannungsexperimente haben dabei das Verstehen der Phasenseparationsmechanismen möglich gemacht.rnAus Quellexperimenten dichtgepfropfter PMMA Bürsten, wurden Lösungsmittel-Polymer Wechselwirkungsparameter (-Parameter) bestimmt. Dabei wurde festgestellt, dass sich die erhaltenen Parameter aufgrund von Filmbenetzung und entropischen Effekten maßgeblich von den errechneten Bulkwerten unterscheiden. Weiterhin wurden nicht reversible Kettenverschlaufungseffekt beobachtet.

Nuclear structure studies in the xenon and radon region and the discovery of a new radon isotope by Penning trap mass spectrometry

Heutzutage gewähren hochpräzise Massenmessungen mit Penning-Fallen tiefe Einblicke in die fundamentalen Eigenschaften der Kernmaterie. Zu diesem Zweck wird die freie Zyklotronfrequenz eines Ions bestimmt, das in einem starken, homogenen Magnetfeld gespeichert ist. Am ISOLTRAP-Massenspektrometer an ISOLDE / CERN können die Massen von kurzlebigen, radioaktiven Nukliden mit Halbwertszeiten bis zu einigen zehn ms mit einer Unsicherheit in der Größenordnung von 10^-8 bestimmt werden. ISOLTRAP besteht aus einem Radiofrequenz-Quadrupol zum akkumulieren der von ISOLDE gelieferten Ionen, sowie zwei Penning-Fallen zum säubern und zur Massenbestimmung der Ionen. Innerhalb dieser Arbeit wurden die Massen von neutronenreichen Xenon- und Radonisotopen (138-146Xe und 223-229Rn) gemessen. Für elf davon wurde zum ersten Mal die Masse direkt bestimmt; 229Rn wurde im Zuge dieses Experimentes sogar erstmalig beobachtet und seine Halbwertszeit konnte zu ungefähr 12 s bestimmt werden. Da die Masse eines Nuklids alle Wechselwirkungen innerhalb des Kerns widerspiegelt, ist sie einzigartig für jedes Nuklid. Eine dieser Wechselwirkungen, die Wechselwirkung zwischen Protonen und Neutronen, führt zum Beispiel zu Deformationen. Das Ziel dieser Arbeit ist eine Verbindung zwischen kollektiven Effekten, wie Deformationen und Doppeldifferenzen von Bindungsenergien, sogenannten deltaVpn-Werten zu finden. Insbesondere in den hier untersuchten Regionen zeigen deltaVpn-Werte ein sehr ungewöhnliches Verhalten, das sich nicht mit einfachen Argumenten deuten lässt. Eine Erklärung könnte das Auftreten von Oktupoldeformationen in diesen Gebieten sein. Nichtsdestotrotz ist eine quantitative Beschreibung von deltaVpn-Werten, die den Effekt von solchen Deformationen berücksichtigt mit modernen Theorien noch nicht möglich.

Effect of mixing, particle interaction, and spatial dimension on the glass transition

In this thesis, the influence of composition changes on the glass transition behavior of binary liquids in two and three spatial dimensions (2D/3D) is studied in the framework of mode-coupling theory (MCT).The well-established MCT equations are generalized to isotropic and homogeneous multicomponent liquids in arbitrary spatial dimensions. Furthermore, a new method is introduced which allows a fast and precise determination of special properties of glass transition lines. The new equations are then applied to the following model systems: binary mixtures of hard disks/spheres in 2D/3D, binary mixtures of dipolar point particles in 2D, and binary mixtures of dipolar hard disks in 2D. Some general features of the glass transition lines are also discussed. The direct comparison of the binary hard disk/sphere models in 2D/3D shows similar qualitative behavior. Particularly, for binary mixtures of hard disks in 2D the same four so-called mixing effects are identified as have been found before by Götze and Voigtmann for binary hard spheres in 3D [Phys. Rev. E 67, 021502 (2003)]. For instance, depending on the size disparity, adding a second component to a one-component liquid may lead to a stabilization of either the liquid or the glassy state. The MCT results for the 2D system are on a qualitative level in agreement with available computer simulation data. Furthermore, the glass transition diagram found for binary hard disks in 2D strongly resembles the corresponding random close packing diagram. Concerning dipolar systems, it is demonstrated that the experimental system of König et al. [Eur. Phys. J. E 18, 287 (2005)] is well described by binary point dipoles in 2D through a comparison between the experimental partial structure factors and those from computer simulations. For such mixtures of point particles it is demonstrated that MCT predicts always a plasticization effect, i.e. a stabilization of the liquid state due to mixing, in contrast to binary hard disks in 2D or binary hard spheres in 3D. It is demonstrated that the predicted plasticization effect is in qualitative agreement with experimental results. Finally, a glass transition diagram for binary mixtures of dipolar hard disks in 2D is calculated. These results demonstrate that at higher packing fractions there is a competition between the mixing effects occurring for binary hard disks in 2D and those for binary point dipoles in 2D.

A soft-tetramer model for diblock copolymer melts: structure formation and geometric characterization

The ability of block copolymers to spontaneously self-assemble into a variety of ordered nano-structures not only makes them a scientifically interesting system for the investigation of order-disorder phase transitions, but also offers a wide range of nano-technological applications. The architecture of a diblock is the most simple among the block copolymer systems, hence it is often used as a model system in both experiment and theory. We introduce a new soft-tetramer model for efficient computer simulations of diblock copolymer melts. The instantaneous non-spherical shape of polymer chains in molten state is incorporated by modeling each of the two blocks as two soft spheres. The interactions between the spheres are modeled in a way that the diblock melt tends to microphase separate with decreasing temperature. Using Monte Carlo simulations, we determine the equilibrium structures at variable values of the two relevant control parameters, the diblock composition and the incompatibility of unlike components. The simplicity of the model allows us to scan the control parameter space in a completeness that has not been reached in previous molecular simulations.The resulting phase diagram shows clear similarities with the phase diagram found in experiments. Moreover, we show that structural details of block copolymer chains can be reproduced by our simple model.We develop a novel method for the identification of the observed diblock copolymer mesophases that formalizes the usual approach of direct visual observation,using the characteristic geometry of the structures. A cluster analysis algorithm is used to determine clusters of each component of the diblock, and the number and shape of the clusters can be used to determine the mesophase.We also employ methods from integral geometry for the identification of mesophases and compare their usefulness to the cluster analysis approach.To probe the properties of our model in confinement, we perform molecular dynamics simulations of atomistic polyethylene melts confined between graphite surfaces. The results from these simulations are used as an input for an iterative coarse-graining procedure that yields a surface interaction potential for the soft-tetramer model. Using the interaction potential derived in that way, we perform an initial study on the behavior of the soft-tetramer model in confinement. Comparing with experimental studies, we find that our model can reflect basic features of confined diblock copolymer melts.

Characterization of structure and dynamics of submicrometric systems via fluorescence correlation spectroscopy

The presented thesis revolves around the study of thermally-responsive PNIPAAm-based hydrogels in water/based environments, as studied by Fluorescence Correlation Spectroscopy (FCS).rnThe goal of the project was the engineering of PNIPAAm gels into biosensors. Specifically, a gamma of such gels were both investigated concerning their dynamics and structure at the nanometer scale, and their performance in retaining bound bodies upon thermal collapse (which PNIPAAm undergoes upon heating above 32 ºC).rnFCS’s requirements, as a technique, match the limitations imposed by the system. Namely, the need to intimately probe a system in a solvent, which was also fragile and easy to alter. FCS, on the other hand, both requires a fluid environment to work, and is based on the observation of diffusion of fluorescents at nanomolar concentrations. FCS was applied to probe the hydrogels on the nanometer size with minimal invasivity.rnVariables in the gels were addressed in the project including crosslinking degree; structural changes during thermal collapse; behavior in different buffers; the possibility of decreasing the degree of inhomogeneity; behavior of differently sized probes; and the effectiveness of antibody functionalization upon thermal collapse.rnThe evidenced results included the heightening of structural inhomogeneities during thermal collapse and under different buffer conditions; the use of annealing to decrease the inhomogeneity degree; the use of differently sized probes to address different length scale of the gel; and the successful functionalization before and after collapse.rnThe thesis also addresses two side projects, also carried forward via FCS. One, diffusion in inverse opals, produced a predictive simulation model for diffusion of bodies in confined systems as dependent on the bodies’ size versus the characteristic sizes of the system. The other was the observation of interaction of bodies of opposite charge in a water solution, resulting in a phenomenological theory and an evaluation method for both the average residence time of the different bodies together, and their attachment likelihood.

Structure and properties of nanostructures from polyelectrolytes and porphyrins: Self-Assembly in aqueous solution

In this work self-assembling model systems in aqueous solution were studied. The systems contained charged polymers, polyelectrolytes, that were combined with oppositely charged counterions to build up supramolecular structures. With imaging, scattering and spectroscopic techniques it was investigated how the structure of building units influences the structure of their assemblies. Polyelectrolytes with different chemical structure, molecular weight and morphology were investigated. In addition to linear polyelectrolytes, semi-flexible cylindrical bottle-brush polymers that possess a defined cross-section and a relatively high persistence along the backbone were studied. The polyelectrolytes were combined with structural organic counterions having charge numbers one to four. Especially the self-assembly of polyelectrolytes with different tetravalent water-soluble porphyrins was studied. Porphyrins have a rigid aromatic structure that has a structural effect on their self-assembly behavior and through which porphyrins are capable of self-aggregation via π-π interaction. The main focus of the thesis is the self-assembly of cylindrical bottle-brush polyelectrolytes with tetravalent porphyrins. It was shown that the addition of porphyrins to oppositely charged brush molecules induces a hierarchical formation of stable nanoscale brush-porphyrin networks. The networks can be disconnected by addition of salt and single porphyrin-decoratedrncylindrical brush polymers are obtained. These two new morphologies, brush-porphyrin networks and porphyrin-decorated brush polymers, may have potential as functional materials with interesting mechanical and optical properties.

Development of a field method to assess root parasites of woody plants - ecological analyses of the grapevine-phylloxera interaction

The development and the growth of plants is strongly affected by the interactions between roots, rootrnassociated organisms and rhizosphere communities. Methods to assess such interactions are hardly torndevelop particularly in perennial and woody plants, due to their complex root system structure and theirrntemporal change in physiology patterns. In this respect, grape root systems are not investigated veryrnwell. The aim of the present work was the development of a method to assess and predict interactionsrnat the root system of rootstocks (Vitis berlandieri x Vitis riparia) in field. To achieve this aim, grapernphylloxera (Daktulosphaira vitifoliae Fitch, Hemiptera, Aphidoidea) was used as a graperoot parasitizingrnmodel.rnTo develop the methodical approach, a longt-term trial (2006-2009) was arranged on a commercial usedrnvineyard in Geisenheim/Rheingau. All 2 to 8 weeks the top most 20 cm of soil under the foliage wallrnwere investigated and root material was extracted (n=8-10). To include temporal, spatial and cultivarrnspecific root system dynamics, the extracted root material was analyzed digitally on the morphologicalrnproperties. The grape phylloxera population was quantified and characterized visually on base of theirrnlarvalstages (oviparous, non oviparous and winged preliminary stages). Infection patches (nodosities)rnwere characterized visually as well, partly supported by digital root color analyses. Due to the knownrneffects of fungal endophytes on the vitality of grape phylloxera infested grapevines, fungal endophytesrnwere isolated from nodosity and root tissue and characterized (morphotypes) afterwards. Further abioticrnand biotic soil conditions of the vineyards were assessed. The temporal, spatial and cultivar specificrnsensitivity of single parameters were analyzed by omnibus tests (ANOVAs) and adjacent post-hoc tests.rnThe relations between different parameters were analyzed by multiple regression models.rnQuantitative parameters to assess the degeneration of nodosity, the development nodosity attachedrnroots and to differentiate between nodosities and other root swellings in field were developed. Significantrndifferences were shown between root dynamic including parameters and root dynamic ignoringrnparameters. Regarding the description of grape phylloxera population and root system dynamic, thernmethod showed a high temporal, spatial and cultivar specific sensitivity. Further, specific differencesrncould be shown in the frequency of endophyte morphotypes between root and nodosity tissue as wellrnas between cultivars. Degeneration of nodosities as well as nodosity occupation rates could be relatedrnto the calculated abundances of grape phylloxera population. Further ecological questions consideringrngrape root development (e.g. relation between moisture and root development) and grape phylloxerarnpopulation development (e.g. relation between temperature and population structure) could be answeredrnfor field conditions.rnGenerally, the presented work provides an approach to evaluate vitality of grape root systems. Thisrnapproach can be useful, considering the development of control strategies against soilborne pests inrnviticulture (e.g. grape phylloxera, Sorospheara viticola, Roesleria subterranea (Weinm.) Redhaed) as well as considering the evaluation of integrated management systems in viticulture.

Structure formation of amphiphilic molecules at the air,water interface and after film transfer

One of the basic concepts of molecular self-assembly is that the morphology of the aggregate is directly related to the structure and interaction of the aggregating molecules. This is not only true for the aggregation in bulk solution, but also for the formation of Langmuir films at the air/water interface. Thus, molecules at the interface do not necessarily form flat monomolecular films but can also aggregate into multilayers or surface micelles. In this context, various novel synthetic molecules were investigated in terms of their morphology at the air/water interface and in transferred films. rnFirst, the self-assembly of semifluorinated alkanes and their molecular orientation at the air/water interface and in transferred films was studied employing scanning force microscopy (SFM) and Kelvin potential force microscopy. Here it was found, that the investigated semifluorinated alkanes aggregate to form circular surface micelles with a diameter of 30 nm, which are constituted of smaller muffin-shaped subunits with a diameter of 10 nm. A further result is that the introduction of an aromatic core into the molecular structure leads to the formation of elongated surface micelles and thus implements a directionality to the self-assembly. rnSecond, the self-assembly of two different amphiphilic hybrid materials containing a short single stranded desoxyribonucleic acid (DNA) sequence was investigated at the air/water interface. The first molecule was a single stranded DNA (11mer) molecule with two hydrophobically modified 5-(dodec-1-ynyl)uracil nucleobases at the terminal 5'-end of the oligonucleotide sequence. Isotherm measurements revealed the formation of semi-stable films at the air/water interface. SFM imaging of films transferred via Langmuir-Blodgett technique supported this finding and indicated mono-, bi- and multilayer formation, according to the surface pressure applied upon transfer. Within these films, the hydrophilic DNA sequence was oriented towards air covering 95% of the substrate.rnSimilar results were obtained with a second type of amphiphile, a DNA block copolymer. Furthermore, the potential to perform molecular recognition experiments at the air/water interface with these DNA hybrid materials was evaluated.rnThird, polyglycerol ester molecules (PGE), which are known to form very stable foams, were studies. Aim was to elucidate the molecular structure of PGE molecules at the air/water interface in order to comprehend the foam stabilization mechanism. Several model systems mimicking the air/water interface of a PGE foam and methods for a noninvasive transfer were tested and characterized by SFM. It could be shown, that PGE stabilizes the air/water interface of a foam bubble by formation of multiple surfactant layers. Additionally, a new transfer technique, the bubble film transfer was established and characterized by high speed camera imaging.The results demonstrate the diversity of structures, which can be formed by amphiphilic molecules at the air/water interface and after film transfer, as well as the impact of the chemical structure on the aggregate morphology.

Nucleon electromagnetic structure studies in the spacelike and timelike regions

The thesis investigates the nucleon structure probed by the electromagnetic interaction. One of the most basic observables, reflecting the electromagnetic structure of the nucleon, are the form factors, which have been studied by means of elastic electron-proton scattering with ever increasing precision for several decades. In the timelike region, corresponding with the proton-antiproton annihilation into a electron-positron pair, the present experimental information is much less accurate. However, in the near future high-precision form factor measurements are planned. About 50 years after the first pioneering measurements of the electromagnetic form factors, polarization experiments stirred up the field since the results were found to be in striking contradiction to the findings of previous form factor investigations from unpolarized measurements. Triggered by the conflicting results, a whole new field studying the influence of two-photon exchange corrections to elastic electron-proton scattering emerged, which appeared as the most likely explanation of the discrepancy. The main part of this thesis deals with theoretical studies of two-photon exchange, which is investigated particularly with regard to form factor measurements in the spacelike as well as in the timelike region. An extraction of the two-photon amplitudes in the spacelike region through a combined analysis using the results of unpolarized cross section measurements and polarization experiments is presented. Furthermore, predictions of the two-photon exchange effects on the e+p/e-p cross section ratio are given for several new experiments, which are currently ongoing. The two-photon exchange corrections are also investigated in the timelike region in the process pbar{p} -> e+ e- by means of two factorization approaches. These corrections are found to be smaller than those obtained for the spacelike scattering process. The influence of the two-photon exchange corrections on cross section measurements as well as asymmetries, which allow a direct access of the two-photon exchange contribution, is discussed. Furthermore, one of the factorization approaches is applied for investigating the two-boson exchange effects in parity-violating electron-proton scattering. In the last part of the underlying work, the process pbar{p} -> pi0 e+e- is analyzed with the aim of determining the form factors in the so-called unphysical, timelike region below the two-nucleon production threshold. For this purpose, a phenomenological model is used, which provides a good description of the available data of the real photoproduction process pbar{p} -> pi0 gamma.

Organic lasing: correlation between molecular structure, optical and optoelectronic properties

In this thesis mainly two alternating indenofluorene-phenanthrene copolymers were investigated with a variety of spectroscopic and optoelectronic experiments. The different experimental techniques allowed to retrieve deeper insights into their unique optical as well as optoelectronic properties. The motivation of the research presented in this work was to correlate their photophysical properties with respect to their application in electrically pumped lasing. This thesis begins with the description of optical properties studied by classical absorption and emission spectroscopy and successively describes an overall picture regarding their excited state dynamics occurring after photoexcitation studied by time-resolved spectroscopy. The different spectroscopic methods do not only allow to elucidate the different optical transitions occurring in this class of materials, but also contribute to a better understanding of exciton dynamics and exciton interaction with respect to the molecular structure as well as aggregation and photooxidation of the polymers. Furthermore, the stimulated emission properties were analyzed by amplified spontaneous emission (ASE) experiments. Especially one of the investigated materials, called BLUE-1, showed outstanding optical properties including a high optical gain, a low threshold for ASE and low optical losses. Apart from the optical experiments, the charge carrier mobility was measured with the time-of-flight technique and a comparably high hole mobility on the order of 1 x 10-² cm²/(Vs) was determined for BLUE-1 which makes this material promising for organic lasing. The impact of the high charge carrier mobility in this material class was further analyzed in different optoelectronic devices such as organic LEDs (OLEDs) and organic solar cells.

Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.

Numerical analysis of the relation between interactions and structure in a molecular fluid

Coarse graining is a popular technique used in physics to speed up the computer simulation of molecular fluids. An essential part of this technique is a method that solves the inverse problem of determining the interaction potential or its parameters from the given structural data. Due to discrepancies between model and reality, the potential is not unique, such that stability of such method and its convergence to a meaningful solution are issues.rnrnIn this work, we investigate empirically whether coarse graining can be improved by applying the theory of inverse problems from applied mathematics. In particular, we use the singular value analysis to reveal the weak interaction parameters, that have a negligible influence on the structure of the fluid and which cause non-uniqueness of the solution. Further, we apply a regularizing Levenberg-Marquardt method, which is stable against the mentioned discrepancies. Then, we compare it to the existing physical methods - the Iterative Boltzmann Inversion and the Inverse Monte Carlo method, which are fast and well adapted to the problem, but sometimes have convergence problems.rnrnFrom analysis of the Iterative Boltzmann Inversion, we elaborate a meaningful approximation of the structure and use it to derive a modification of the Levenberg-Marquardt method. We engage the latter for reconstruction of the interaction parameters from experimental data for liquid argon and nitrogen. We show that the modified method is stable, convergent and fast. Further, the singular value analysis of the structure and its approximation allows to determine the crucial interaction parameters, that is, to simplify the modeling of interactions. Therefore, our results build a rigorous bridge between the inverse problem from physics and the powerful solution tools from mathematics. rn