Helix-coil-transition in solid polypeptides

Abstract Poly(L-glutamic acid) (PLGA) was synthesized by living anionic ring-opening polymerization of the NCA monomer, which was obtained by reacting diphosgene with an amino acid derivative. The chemical structures and thermal properties were characterized by 1H-NMR, 13C-NMR, TGA and DSC. XRD powder patterns found to be amorphous for all polymers obtained. The molecular weights could be determined under severe limitations due to low solubility and high aggregation tendency. The secondary structure of the PLGA films was analyzed in the solid state by IR spectroscopy; the order was determined mainly by XRD. Uniform bulk films (1-5 µm) were produced by drop-casting of PLGA solutions in TFA on silica. The XRD film analysis indicated the absence of a long range order or an orientation even if a helical microstructure was confirmed by IR spectroscopy. The coil solvent TFA delivered constantly a helical or a β-sheet structure in the solid state depending on the water content of the solvent which was observed for the first time to exhibit a high influence on the crystallization process for PLGA. Temperature dependent in-situ IR measurements were examined to analyze if a helix-coil transition occurs, but there could be no solvent system determined, which resulted in a disordered coil structure in the solid state. General parameters like solvent systems, evaporation conditions, concentration, substrates etc. were analyzed. New crystallizations were obtained on silica prepared by drop-casting of solutions of PLGA in DMF, DMA, TMU, NMP, and pyridine/water mixtures, respectively. PSCBC in DMF, CDCl3/TFA-d, and PSBC in CDCl3/TFA-d exhibited the same crystalline diffraction patterns like PLGA. The long range order in the X-ray diffraction pattern is proven by extremely sharp crystalline signals, which are not changing the shape or the position of the peak by increasing the temperature up to 160°C. The substrate seems to play a decisive role because the crystalline structures were not obtainable on glass. The crystal structure consists probably of two different layered structures based on the intensity ratios of the two series of crystalline signals in the X-ray diffraction patterns. The source of the layered structure remains unclear and needs further studies to investigate the spatial arrangement of the chains in more detail. The secondary structure was still not changing upon heating even if a highly crystalline diffraction pattern occurs. Concluding that even the newly investigated crystallization did not show a helix-coil transition in the solid state by annealing, the phenomenon known in solution has to be claimed as unachievable in the solid state based on the results of this work. A remaining open question represents the observation that the same crystalline pattern can be reproducibly prepared with exhibiting two different ordered secondary structures (helix and β-sheet). After the investigation that the evaporation time cannot be decisive for the crystal growth, the choice of a strong hydrogen bonding interrupting solvent is most probably the key to support and induce the crystallization process.

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).

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.