Complexes of polyelectrolytes with defined charge distance and different dendrimer counterions

Complexes of polyelectrolytes with defined charge distance and different dendrimer counterions Magdalena Chelmecka Max Planck Institute for Polymer Research; Ackermannweg 10; D-55128 Mainz ; Tel.: (+49) 06131- 379 – 226 A study of complexes in solution is of interest to investigate whether the formation of well-defined assemblies like in classical surfactant systems is possible. Aim of this thesis is to investigate the electrostatic self-assembly of linear polycations of varying charge distance with “large” counterions of varying architecture. We especially investigate the morphology of objects formed, but also their stability under salt free condition and after low molecular mass salt addition. As polycations, Poly(dialkylimino)-alkylene salts (Ionenes) I65MeBr and I25MeBr were chosen. Ionenes are synthesized via Menschutkin reaction and characterized by standard methods. Counterions are Polyamidoamine (PAMAM) dendrimers of generations G2.5, G5.5, G7.5 with -COONa surface groups and shape-persistent, Polyphenylene dendrimers of generation G1 with surface -COOH groups. A complex interplay of interactions is expected to direct the self assembly via electrostatic interaction, geometric factors, hydrophobic interaction or hydrogen bonds. Methods used for the investigation of complexes are: UV-spectroscopy, pH-metric techniques, dynamic and static light scattering, small angle neutron scattering,  potential measurements and potentiometric titration. Under certain conditions, (i.e. charge ratio of compounds, charge density of ionene and dendrimer also concentration of sample) polyelectrolyte systems composed of ionenes and dendrimers build complexes in solution. System compounds are typical polyelectrolytes, but structures which they build behave not usual for typical polyelectrolytes. In a one diffusion mode regime aggregates of about 100 nm hydrodynamic radius have been found. Such aggregates are core-shell or anisotropic core shell structures in the case of ionenes/PAMAM dendrimers complexes. These complexes are stable even at high ionic strength. In case of ionenes with poly(phenylene) dendrimers, hard sphere-like objects or spherical objects with hairy-like surface have been found in a one diffusion mode regime. Their stability at high ionic strength is lower. For the ionenes/poly(phenylene) dendrimers systems one transition point has been found from one to two diffusion processes, towards increasing ionene concentration, i.e. for the samples with fixed dendrimer concentration towards increasing ionic strength. For the diffusion profile of ionene/PAMAM dendrimers in most cases two transition regimes are observed. One at very low ionene concentration, the second one at high ionene concentrations, which again means for the samples with fixed dendrimer concentration, also at higher ionic strength. Both two mode regimes are separated by the one mode regime. As was confirmed experimentally, the one diffusion mode regime is caused by the motion of well defined assemblies. The two diffusion mode regimes are caused by the movement of different sized species in solution, large aggregates and middle-size aggregates (oligoaggregates). The location and also the number of transition points in the diffusion profiles is dependent on the ionene to dendrimer charge ratio, charge density of the compounds and concentration. No influence of the molecular mass of the ionene has been found. The aggregates are found to be charged on the surface, however this surface charge does not significantly influence the diffusion properties of the system.

Mixed finite-element methods for elliptic convection-dominated problems arising in semiconductor physics

In this work we develop and analyze an adaptive numerical scheme for simulating a class of macroscopic semiconductor models. At first the numerical modelling of semiconductors is reviewed in order to classify the Energy-Transport models for semiconductors that are later simulated in 2D. In this class of models the flow of charged particles, that are negatively charged electrons and so-called holes, which are quasi-particles of positive charge, as well as their energy distributions are described by a coupled system of nonlinear partial differential equations. A considerable difficulty in simulating these convection-dominated equations is posed by the nonlinear coupling as well as due to the fact that the local phenomena such as "hot electron effects" are only partially assessable through the given data. The primary variables that are used in the simulations are the particle density and the particle energy density. The user of these simulations is mostly interested in the current flow through parts of the domain boundary - the contacts. The numerical method considered here utilizes mixed finite-elements as trial functions for the discrete solution. The continuous discretization of the normal fluxes is the most important property of this discretization from the users perspective. It will be proven that under certain assumptions on the triangulation the particle density remains positive in the iterative solution algorithm. Connected to this result an a priori error estimate for the discrete solution of linear convection-diffusion equations is derived. The local charge transport phenomena will be resolved by an adaptive algorithm, which is based on a posteriori error estimators. At that stage a comparison of different estimations is performed. Additionally a method to effectively estimate the error in local quantities derived from the solution, so-called "functional outputs", is developed by transferring the dual weighted residual method to mixed finite elements. For a model problem we present how this method can deliver promising results even when standard error estimator fail completely to reduce the error in an iterative mesh refinement process.

Synthetical engineering of supramolecular properties of large polycyclic aromatic hydrocarbons

An important property for devices is the charge-carrier mobility values for discotic organic materials like hexa-peri-hexabenzocoronenes. A close relation exists between the degree of their columnar self-arrangement of the molecules and their mobilities. Within this first step an induction of a higher order via hydrogen-bonding was considered, which mainly pointed towards the improvement of the intracolumnar stacking of the materials. For the analytics a broad range of methods was used including differential scanning calorimetry (DSC), wide-angle X-ray diffractometry (WAXS), solid-state NMR spectroscopy and scanning tunneling microscopy (STM). Indeed, a specific influence of the hydrogen-bonds could be identified, although in several cases by the cost of a severe reduction of solubility and processability. This effect was dampened by the addition of a long alkyl chain next to the hydrogen-bond exerting functional group, which resulted in an improved columnar arrangement by retention of processability. In contrast to the before mentioned example of inducing a higher intracolumnar order by hydrogen-bonding, the focus was also be set upon larger aromatic systems. The charge-carrier mobility is also in close relation to the size of the aromatic core and larger π-areas are expected to lead to improved mobilities. For photovoltaic applications a high extinction coefficient over a broad range of the spectrum is favorable, which can also be achieved by enlarging the aromatic core component. In addition the stronger π-interactions between the aromatic core components should yield an improved columnar stability and order. However the strengthening of the π-interactions between the aromatic core components led to a reduction of the solubility and the processability due to the stronger aggregation of the molecules. This required the introduction of efficiently solubilizing features in terms of long alkyl chains in the corona of the aromatic entity, in combination of a distortion of the aromatic core moiety by bulky tert-butyl groups. By this approach not only the processing and cleaning of the materials with standard laboratory techniques became possible, but moreover the first structure-rich UV/vis and a resolved 1H-NMR spectra for an aromatic system two times larger than hexa-peri-hexabenzocoronene were recorded. The bulk properties in an extruded fiber as well as on the surface showed a columnar self-assembly including a phase in which a homeotropic alignment on a substrate was observed, which turns the material into an interesting candidate for future applications in electronic devices.

Nuclear charge radius determination of 7,10 Be and the one-neutron halo nucleus 11 Be

The only nuclear model independent method for the determination of nuclear charge radii of short-lived radioactive isotopes is the measurement of the isotope shift. For light elements (Z < 10) extremely high accuracy in experiment and theory is required and was only reached for He and Li so far. The nuclear charge radii of the lightest elements are of great interest because they have isotopes which exhibit so-called halo nuclei. Those nuclei are characterized by a a very exotic nuclear structure: They have a compact core and an area of less dense nuclear matter that extends far from this core. Examples for halo nuclei are 6^He, 8^He, 11^Li and 11^Be that is investigated in this thesis. Furthermore these isotopes are of interest because up to now only for such systems with a few nucleons the nuclear structure can be calculated ab-initio. In the Institut für Kernchemie at the Johannes Gutenberg-Universität Mainz two approaches with different accuracy were developed. The goal of these approaches was the measurement of the isotope shifts between (7,10,11)^Be^+ and 9^Be^+ in the D1 line. The …first approach is laser spectroscopy on laser cooled Be^+ ions that are trapped in a linear Paul trap. The accessible accuracy should be in the order of some 100 kHz. In this thesis two types of linear Paul traps were developed for this purpose. Moreover, the peripheral experimental setup was simulated and constructed. It allows the efficient deceleration of fast ions with an initial energy of 60 keV down to some eV and an effcient transport into the ion trap. For one of the Paul traps the ion trapping could already be demonstrated, while the optical detection of captured 9^Be^+ ions could not be completed, because the development work was delayed by the second approach. The second approach uses the technique of collinear laser spectroscopy that was already applied in the last 30 years for measuring isotope shifts of plenty of heavier isotopes. For light elements (Z < 10), it was so far not possible to reach the accuracy that is required to extract information about nuclear charge radii. The combination of collinear laser spectroscopy with the most modern methods of frequency metrology …finally permitted the …first-time determination of the nuclear charge radii of (7,10)^Be and the one neutron halo nucleus 11^Be at the COLLAPS experiment at ISOLDE/ CERN. In the course of the work reported in this thesis it was possible to measure the absolute transition frequencies and the isotope shifts in the D1 line for the Be isotopes mentioned above with an accuracy of better than 2 MHz. Combination with the most recent calculations of the mass effect allowed the extraction of the nuclear charge radii of (7,10,11)^Be with an relative accuracy better than 1%. The nuclear charge radius decreases from 7^Be continuously to 10^Be and increases again for 11^Be. This result is compared with predictions of ab-initio nuclear models which reproduce the observed trend. Particularly the "Greens Function Monte Carlo" and the "Fermionic Molecular Dynamic" model show very good agreement.

Simulations of charge transport in organic compounds

To aid the design of organic semiconductors, we study the charge transport properties of organic liquid crystals, i.e. hexabenzocoronene and carbazole macrocycle, and single crystals, i.e. rubrene, indolocarbazole and benzothiophene derivatives (BTBT, BBBT). The aim is to find structure-property relationships linking the chemical structure as well as the morphology with the bulk charge carrier mobility of the compounds. To this end, molecular dynamics (MD) simulations are performed yielding realistic equilibrated morphologies. Partial charges and molecular orbitals are calculated based on single molecules in vacuum using quantum chemical methods. The molecular orbitals are then mapped onto the molecular positions and orientations, which allows calculation of the transfer integrals between nearest neighbors using the molecular orbital overlap method. Thus we obtain realistic transfer integral distributions and their autocorrelations. In case of organic crystals the differences between two descriptions of charge transport, namely semi-classical dynamics (SCD) in the small polaron limit and kinetic Monte Carlo (KMC) based on Marcus rates, are studied. The liquid crystals are investigated solely in the hopping limit. To simulate the charge dynamics using KMC, the centers of mass of the molecules are mapped onto lattice sites and the transfer integrals are used to compute the hopping rates. In the small polaron limit, where the electronic wave function is spread over a limited number of neighboring molecules, the Schroedinger equation is solved numerically using a semi-classical approach. The results are compared for the different compounds and methods and, where available, with experimental data. The carbazole macrocycles form columnar structures arranged on a hexagonal lattice with side chains facing inwards, so columns can closely approach each other allowing inter-columnar and thus three-dimensional transport. When taking only intra-columnar transport into account, the mobility is orders of magnitude lower than in the three-dimensional case. BTBT is a promising material for solution-processed organic field-effect transistors. We are able to show that, on the time-scales of charge transport, static disorder due to slow side chain motions is the main factor determining the mobility. The resulting broad transfer integral distributions modify the connectivity of the system but sufficiently many fast percolation paths remain for the charges. Rubrene, indolocarbazole and BBBT are examples of crystals without significant static disorder. The high mobility of rubrene is explained by two main features: first, the shifted cofacial alignment of its molecules, and second, the high center of mass vibrational frequency. In comparsion to SCD, only KMC based on Marcus rates is capable of describing neighbors with low coupling and of taking static disorder into account three-dimensionally. Thus it is the method of choice for crystalline systems dominated by static disorder. However, it is inappropriate for the case of strong coupling and underestimates the mobility of well-ordered crystals. SCD, despite its one-dimensionality, is valuable for crystals with strong coupling and little disorder. It also allows correct treatment of dynamical effects, such as intermolecular vibrations of the molecules. Rate equations are incapable of this, because simulations are performed on static snapshots. We have thus shown strengths and weaknesses of two state of the art models used to study charge transport in organic compounds, partially developed a program to compute and visualize transfer integral distributions and other charge transport properties, and found structure-mobility relations for several promising organic semiconductors.

Morphology and charge transport in conjugated polymers

To assist rational compound design of organic semiconductors, two problems need to be addressed. First, the material morphology has to be known at an atomistic level. Second, with the morphology at hand, an appropriate charge transport model needs to be developed in order to link charge carrier mobility to structure.rnrnThe former can be addressed by generating atomistic morphologies using molecular dynamics simulations. However, the accessible range of time- and length-scales is limited. To overcome these limitations, systematic coarse-graining methods can be used. In the first part of the thesis, the Versatile Object-oriented Toolkit for Coarse-graining Applications is introduced, which provides a platform for the implementation of coarse-graining methods. Tools to perform Boltzmann inversion, iterative Boltzmann inversion, inverse Monte Carlo, and force-matching are available and have been tested on a set of model systems (water, methanol, propane and a single hexane chain). Advantages and problems of each specific method are discussed.rnrnIn partially disordered systems, the second issue is closely connected to constructing appropriate diabatic states between which charge transfer occurs. In the second part of the thesis, the description initially used for small conjugated molecules is extended to conjugated polymers. Here, charge transport is modeled by introducing conjugated segments on which charge carriers are localized. Inter-chain transport is then treated within a high temperature non-adiabatic Marcus theory while an adiabatic rate expression is used for intra-chain transport. The charge dynamics is simulated using the kinetic Monte Carlo method.rnrnThe entire framework is finally employed to establish a relation between the morphology and the charge mobility of the neutral and doped states of polypyrrole, a conjugated polymer. It is shown that for short oligomers, charge carrier mobility is insensitive to the orientational molecular ordering and is determined by the threshold transfer integral which connects percolating clusters of molecules that form interconnected networks. The value of this transfer integral can be related to the radial distribution function. Hence, charge mobility is mainly determined by the local molecular packing and is independent of the global morphology, at least in such a non-crystalline state of a polymer.

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