Existence of solutions and asymtptotic limits of the Euler-Poisson and the quantum drift-diffusion systems

My work concerns two different systems of equations used in the mathematical modeling of semiconductors and plasmas: the Euler-Poisson system and the quantum drift-diffusion system. The first is given by the Euler equations for the conservation of mass and momentum, with a Poisson equation for the electrostatic potential. The second one takes into account the physical effects due to the smallness of the devices (quantum effects). It is a simple extension of the classical drift-diffusion model which consists of two continuity equations for the charge densities, with a Poisson equation for the electrostatic potential. Using an asymptotic expansion method, we study (in the steady-state case for a potential flow) the limit to zero of the three physical parameters which arise in the Euler-Poisson system: the electron mass, the relaxation time and the Debye length. For each limit, we prove the existence and uniqueness of profiles to the asymptotic expansion and some error estimates. For a vanishing electron mass or a vanishing relaxation time, this method gives us a new approach in the convergence of the Euler-Poisson system to the incompressible Euler equations. For a vanishing Debye length (also called quasineutral limit), we obtain a new approach in the existence of solutions when boundary layers can appear (i.e. when no compatibility condition is assumed). Moreover, using an iterative method, and a finite volume scheme or a penalized mixed finite volume scheme, we numerically show the smallness condition on the electron mass needed in the existence of solutions to the system, condition which has already been shown in the literature. In the quantum drift-diffusion model for the transient bipolar case in one-space dimension, we show, by using a time discretization and energy estimates, the existence of solutions (for a general doping profile). We also prove rigorously the quasineutral limit (for a vanishing doping profile). Finally, using a new time discretization and an algorithmic construction of entropies, we prove some regularity properties for the solutions of the equation obtained in the quasineutral limit (for a vanishing pressure). This new regularity permits us to prove the positivity of solutions to this equation for at least times large enough.

Factorization method and inclusions of mixed type in an inverse elliptic boundary value problem

In various imaging problems the task is to use the Cauchy data of the solutions to an elliptic boundary value problem to reconstruct the coefficients of the corresponding partial differential equation. Often the examined object has known background properties but is contaminated by inhomogeneities that cause perturbations of the coefficient functions. The factorization method of Kirsch provides a tool for locating such inclusions. In this paper, the factorization technique is studied in the framework of coercive elliptic partial differential equations of the divergence type: Earlier it has been demonstrated that the factorization algorithm can reconstruct the support of a strictly positive (or negative) definite perturbation of the leading order coefficient, or if that remains unperturbed, the support of a strictly positive (or negative) perturbation of the zeroth order coefficient. In this work we show that these two types of inhomogeneities can, in fact, be located simultaneously. Unlike in the earlier articles on the factorization method, our inclusions may have disconnected complements and we also weaken some other a priori assumptions of the method. Our theoretical findings are complemented by two-dimensional numerical experiments that are presented in the framework of the diffusion approximation of optical tomography.

Structure formation as studied by EPR spectroscopy: From simple solutions to nanoparticles

In this thesis, three nitroxide based ionic systems were used to investigate structure and dynamics of their respective solutions in mixed solvents by means of electron paramagnetic resonance (EPR) and electron nuclear double resonance (ENDOR) spectroscopy at X- and W-band (9.5 and 94.5 GHz, respectively). rnFirst, the solvation of the inorganic radical Fremy’s salt (K2ON(SO3)2) in isotope substituted binary solvent mixtures (methanol/water) was investigated by means of high-field (W-band) pulse ENDOR spectroscopy and molecular dynamics (MD) simulations. From the analysis of orientation-selective 1H and 2H ENDOR spectra the principal components of the hyperfine coupling (hfc) tensor for chemically different protons (alcoholic methyl vs. exchangeable protons) were obtained. The methyl protons of the organic solvent approach with a mean distance of 3.5 Å perpendicular to the approximate plane spanned by ON(S)2 of the probe molecule. Exchangeable protons were found to be distributed isotropically, approaching closest to Fremy’s salt from the hydrogen-bonded network around the sulfonate groups. The distribution of exchangeable and methyl protons as found in MD simulations is in full agreement with the ENDOR results. The solvation was found to be similar for the studied solvent ratios between 1:2.3 and 2.3:1 and dominated by an interplay of H-bond (electrostatic) interactions and steric considerations with the NO group merely involved into H-bonds.rnFurther, the conformation of spin labeled poly(diallyldimethylammonium chloride) (PDADMAC) solutions in aqueous alcohol (methanol, ethanol, n-propanol, ethylene glycol, glycerol) mixtures in dependence of divalent sodium sulfate was investigated with double electron-electron resonance (DEER) spectroscopy. The DEER data was analyzed using the worm-like chain model which suggests that in organic-water solvent mixtures the polymer backbones are preferentially solvated by the organic solvent. We found a less serve impact on conformational changes due to salt than usually predicted in polyelectrolyte theory which stresses the importance of a delicate balance of hydrophobic and electrostatic interactions, in particular in the presence of organic solvents.rnFinally, the structure and dynamics of miniemulsions and polymerdispersions prepared with anionic surfactants, that were partially replaced by a spin labeled fatty acid in presence and absence of a lanthanide beta-diketonate complex was characterized by CW EPR spectroscopy. Such miniemulsions form multilayers with the surfactant head group bound to the lanthanide ion. Beta-diketonates were formerly used as NMR shift reagents and nowadays find application as luminescent materials in OLEDs and LCDs and as contrast agent in MRT. The embedding of the complex into a polymer matrix results in an easy processable material. It was found that the structure formation takes place in miniemulsion and is preserved during polymerization. For surfactants with carboxyl-head group a higher order of the alkyl chains and less lateral diffusion is found than for sulfat-head groups, suggesting a more uniform and stronger coordination to the metal ion. The stability of these bilayers depends on the temperature and the used surfactant which should be considered for the used polymerization temperature if a maximum output of the structured regions is wished.