Development of empirical potentials for amorphous silica

Development of empirical potentials for amorphous silica Amorphous silica (SiO2) is of great importance in geoscience and mineralogy as well as a raw material in glass industry. Its structure is characterized as a disordered continuous network of SiO4 tetrahedra. Many efforts have been undertaken to understand the microscopic properties of silica by classical molecular dynamics (MD) simulations. In this method the interatomic interactions are modeled by an effective potential that does not take explicitely into account the electronic degrees of freedom. In this work, we propose a new methodology to parameterize such a potential for silica using ab initio simulations, namely Car-Parrinello (CP) method [Phys. Rev. Lett. 55, 2471 (1985)]. The new potential proposed is compared to the BKS potential [Phys. Rev. Lett. 64, 1955 (1990)] that is considered as the benchmark potential for silica. First, CP simulations have been performed on a liquid silica sample at 3600 K. The structural features so obtained have been compared to the ones predicted by the classical BKS potential. Regarding the bond lengths the BKS tends to underestimate the Si-O bond whereas the Si-Si bond is overestimated. The inter-tetrahedral angular distribution functions are also not well described by the BKS potential. The corresponding mean value of theSiOSi angle is found to be ≃ 147◦, while the CP yields to aSiOSi angle centered around 135◦. Our aim is to fit a classical Born-Mayer/Coulomb pair potential using ab initio calculations. To this end, we use the force-matching method proposed by Ercolessi and Adams [Europhys. Lett. 26, 583 (1994)]. The CP configurations and their corresponding interatomic forces have been considered for a least square fitting procedure. The classical MD simulations with the resulting potential have lead to a structure that is very different from the CP one. Therefore, a different fitting criterion based on the CP partial pair correlation functions was applied. Using this approach the resulting potential shows a better agreement with the CP data than the BKS ones: pair correlation functions, angular distribution functions, structure factors, density of states and pressure/density were improved. At low temperature, the diffusion coefficients appear to be three times higher than those predicted by the BKS model, however showing a similar temperature dependence. Calculations have also been carried out on crystalline samples in order to check the transferability of the potential. The equilibrium geometry as well as the elastic constants of α-quartz at 0 K are well described by our new potential although the crystalline phases have not been considered for the parameterization. We have developed a new potential for silica which represents an improvement over the pair potentials class proposed so far. Furthermore, the fitting methodology that has been developed in this work can be applied to other network forming systems such as germania as well as mixtures of SiO2 with other oxides (e.g. Al2O3, K2O, Na2O).

Spin-crossover in metallomesogens of iron(II)

Covalent grafting mesogenic groups to the coordination cores of the parent mononuclear low-spin and spin-crossover compounds afforded metallomesogenic complexes of iron(II). In comparison with the parent complexes the spin-crossover properties of the alkylated derivatives are substantially modified. The type of the modification was found to be dependent on the properties of the parent system and the nature of the used anion, however, the general tendency is the destabilization of the low-spin state at the favor of spin-crossover or high-spin behavior below 400 K. The structural insight revealed the micro-segregated layered organization. The effect of the alkylation of the parent compounds consists first of all in the change of the lattice to a two-dimensional lamellar one retaining significant intermolecular contacts only within the ionic bilayers. The comprehensive analysis of the structural and thermodynamic data in the homologous series pointed at the mechanism of the interplay between the structural modification on melting and the induced anomalous change of the magnetic properties. A family of one-dimensional spin-crossover polymers was synthesized and characterized using a series of spectroscopic methods, X-ray powder diffraction, magnetic susceptibility measurements and differential scanning calorimetry. The copper analogue of was also synthesized and its crystal structure solved. In comparison with the mononuclear systems, the polymeric mesogens of iron(II) are less sensitive to the glass transition, which was attributed to the moderate concomitant variation of the structure. Nevertheless, the observed increase of the magnetic hysteresis with lengthening of the alkyl substituents was ascribed to the interplay of the structural reorganization of the coordination core due to spin-crossover with the structural delay in the spatial reorganization of the mesogenic substituents. The classification of mononuclear and polymeric metallomesogens according to the interactions between the structural- and the spin-transition and analysis of the data on the reported spin-crossover metallomesogens led to the separation of three types, namely: Type i: systems with coupling between the electronic structure of the iron(II) ions and the mesomorphic behavior of the substance; Type ii: systems where both transitions coexist in the same temperature region but are not coupled due to competition with the dehydration or due to negligible structural transformation; Type iii: systems where both transitions occur in different temperature regions and therefore are uncoupled. Fine-tuning, in particular regarding the temperature at which the spin-transition occurs with hysteresis properties responsible for the memory effect, are still a major challenge towards practical implementation of spin-crossover materials. A possible answer to the problem could be materials in which the spin-crossover transition is coupled with another transition easily controllable by external stimuli. In the present thesis we have shown the viability of the approach realized in the mesogenic systems with coupled phase- and spin-transitions.

Charge transport in amorphous organic semiconductors

Organic semiconductors with the unique combination of electronic and mechanical properties may offer cost-effective ways of realizing many electronic applications, e.g. large-area flexible displays, printed integrated circuits and plastic solar cells. In order to facilitate the rational compound design of organic semiconductors, it is essential to understand relevant physical properties e.g. charge transport. This, however, is not straightforward, since physical models operating on different time and length scales need to be combined. First, the material morphology has to be known at an atomistic scale. For this atomistic molecular dynamics simulations can be employed, provided that an atomistic force field is available. Otherwise it has to be developed based on the existing force fields and first principle calculations. However, atomistic simulations are typically limited to the nanometer length- and nanosecond time-scales. To overcome these limitations, systematic coarse-graining techniques can be used. In the first part of this thesis, it is demonstrated how a force field can be parameterized for a typical organic molecule. Then different coarse-graining approaches are introduced together with the analysis of their advantages and problems. When atomistic morphology is available, charge transport can be studied by combining the high-temperature Marcus theory with kinetic Monte Carlo simulations. The approach is applied to the hole transport in amorphous films of tris(8-hydroxyquinoline)aluminium (Alq3). First the influence of the force field parameters and the corresponding morphological changes on charge transport is studied. It is shown that the energetic disorder plays an important role for amorphous Alq3, defining charge carrier dynamics. Its spatial correlations govern the Poole-Frenkel behavior of the charge carrier mobility. It is found that hole transport is dispersive for system sizes accessible to simulations, meaning that calculated mobilities depend strongly on the system size. A method for extrapolating calculated mobilities to the infinite system size is proposed, allowing direct comparison of simulation results and time-of-flight experiments. The extracted value of the nondispersive hole mobility and its electric field dependence for amorphous Alq3 agree well with the experimental results.

Nanoparticular iron complex drugs for parenteral administration - physicochemical characterization, biological distribution and pharmacological safety

Iron deficiency is the most common deficiency disease worldwide with many patients who require intravenous iron. Within the last years new kind of parenteral iron complexes as well as generic preparations entered the market. There is a high demand for methods clarifying benefit to risk profiles of old and new iron complexes. It is also necessary to disclose interchangeability between originator and intended copies to avoid severe anaphylactic and anaphylactoid side reaction and assure equivalence of therapeutic effect.rnrnThe investigations presented in this work include physicochemical characterization of nine different parenteral iron containing non-biological complex drugs. rnWe developed an in-vitro assay, which allows the quantification of labile iron in the different complexes and thus it is a useful tool to estimate the pharmaclogical safety regarding iron related adverse drug events. This assay additionally allowed the estimation of complex stability by evaluation of degradation kinetics at the applied conditions.rnrnAn in-ovo study was performed to additionally compare different complexes in respect to body distribution. This in combination with complex stability information allowed the risk estimation of potential local acute and chronic reactions to iron overload.rnrnInformation obtained by the combination of the methods within this work are helpful to estimate the safety and efficacy profile of different iron containing non-biological complex drugs. rnrnPhysicochemical differences between the complexes were demonstrated in respect to size of the inorganic fraction, size and size distribution of the complete particles, structure of the inorganic iron fraction, morphology of the complexes and charge of the complexes. And furthermore significant differences in the biological behavior of different complexes were demonstrated. rnrnThe combination of complex stability and biodistribution as well as the combination of structure, size and stability represent helpful tools for the physicochemical characterization of iron containing non-biological complex drugs and for the estimation of pharmacological safety. This work thus represents an up to date summary of some relevant methods for the characterization of intravenous iron complex drugs in respect to pharmaceutical quality, pharmacological safety and aspects of efficacy. rnrnProspectively, it is worthwhile that the methods within this work will contribute to the development and/or characterization of iron containing nanoparticular formulations with beneficial efficacy and safety profiles.rn