Novel simulation methods for Coulomb and hydrodynamic interactions

This thesis presents new methods to simulate systems with hydrodynamic and electrostatic interactions. Part 1 is devoted to computer simulations of Brownian particles with hydrodynamic interactions. The main influence of the solvent on the dynamics of Brownian particles is that it mediates hydrodynamic interactions. In the method, this is simulated by numerical solution of the Navier--Stokes equation on a lattice. To this end, the Lattice--Boltzmann method is used, namely its D3Q19 version. This model is capable to simulate compressible flow. It gives us the advantage to treat dense systems, in particular away from thermal equilibrium. The Lattice--Boltzmann equation is coupled to the particles via a friction force. In addition to this force, acting on {it point} particles, we construct another coupling force, which comes from the pressure tensor. The coupling is purely local, i.~e. the algorithm scales linearly with the total number of particles. In order to be able to map the physical properties of the Lattice--Boltzmann fluid onto a Molecular Dynamics (MD) fluid, the case of an almost incompressible flow is considered. The Fluctuation--Dissipation theorem for the hybrid coupling is analyzed, and a geometric interpretation of the friction coefficient in terms of a Stokes radius is given. Part 2 is devoted to the simulation of charged particles. We present a novel method for obtaining Coulomb interactions as the potential of mean force between charges which are dynamically coupled to a local electromagnetic field. This algorithm scales linearly, too. We focus on the Molecular Dynamics version of the method and show that it is intimately related to the Car--Parrinello approach, while being equivalent to solving Maxwell's equations with freely adjustable speed of light. The Lagrangian formulation of the coupled particles--fields system is derived. The quasi--Hamiltonian dynamics of the system is studied in great detail. For implementation on the computer, the equations of motion are discretized with respect to both space and time. The discretization of the electromagnetic fields on a lattice, as well as the interpolation of the particle charges on the lattice is given. The algorithm is as local as possible: Only nearest neighbors sites of the lattice are interacting with a charged particle. Unphysical self--energies arise as a result of the lattice interpolation of charges, and are corrected by a subtraction scheme based on the exact lattice Green's function. The method allows easy parallelization using standard domain decomposition. Some benchmarking results of the algorithm are presented and discussed.

Untersuchung des Aggregationsverhaltens amphiphiler Diblockcopolymere in überkritischem Kohlendioxid mittels dynamischer Lichtstreuung

„Untersuchung des Aggregationsverhaltens amphiphiler Diblockcopolymere in überkritischem Kohlendioxid mittels dynamischer Lichtstreuung“ In der vorliegenden Arbeit wurde die Mizellenbildung von Diblockcopolymeren des Typs PS-b-PDMS in überkritischem Kohlendioxid (CO2,SC) mittels dynamischer Lichtstreuung (DLS) charakterisiert. Zu diesem Zweck wurden Mischungen aus den Diblockcopolymeren in CO2,SC mit Styrol als Monomer druckabhängig auf diese Fähigkeit hin untersucht. Eine Mizellenbildung konnte anhand der gemessenen hydrodynamischen Radien Rh gezeigt werden. Um eine Vergleichsmöglichkeit gegenüber den mit Styrol gefüllten Kern-Hüllen-Mizellen zu bekommen, wurde das Diblockcopolymer PS-b-PDMS (9/27) zunächst ohne Styrol auf die Fähigkeit hin untersucht ungefüllte Mizellen zu bilden. Durch Druckvariation konnte ein kritischer Mizellendruck von ca. 46,7 MPa bei einer Temperatur von 338 K im Experiment bestätigt werden, der gefundene Rh liegt bei ca. 34 nm. Dagegen setzt die Aggregation bei einer PS-b-PDMS (9/27)/Styrol/CO2,SC- Mischung bei einem wesentlich niedrigeren Druck ein. Durch Druckvariation zwischen 38 MPa und 45,7 MPa wurde eine Größenänderung der Mizellen beobachtet. Durch zeitabhängige-DLS-Messungen am gleichen System bei einem bestimmten Druck wurde ein langsames Schrumpfen der Mizellen gefunden. Um den Einfluß der Blockgröße der verwendeten Amphiphile auf die Mizellenbildung zu untersuchen wurde das System PS-b-PDMS(6/37)/Styrol/CO2,SC mit Hilfe der DLS im Bereich zwischen 39,4 MPa und 43,1 MPa untersucht. Die Druckänderung zeigte für Rh ein nahezu invariantes Verhalten, daß durch eine verlängerte PDMS-Blocklänge und eine damit verbundene Kompensation der verschiedenen Wechselwirkungskräfte zwischen Mizellenkern, -hülle und CO2,SC erklärt werden kann. Im System PS-b-PDMS(6/16)/Styrol/CO2,SC konnte experimentell mit Hilfe der DLS erst nach einer ver-änderten molaren Zusammensetzung eine Mizellenbildung ab 40 MPa ermöglicht werden. Allerdings ändert sich auch in diesem System der hydrodynamische Radius ebenfalls mit dem Druck. Je nach Druck-, Temperatur- und molarer Zusammensetzung variiert die Tendenz der Systeme, Mizellen zu bilden die eine Emulsion stabilisieren können. Für die in Dispersions-Polymerisationsreaktionen eingesetzten Diblockcopolymere bedeutet dieses Ergebnis differenzierte Applikationsmöglichkeiten. Mit den ermittelten Konzentrationsverhältnissen an Amphiphil und Monomer konnte ein Bereich gefunden werden, in dem die thermodynamischen Bedingungen für die Mizellenbildung einerseits und die Vorraussetzungen für die DLS andererseits gegeben sind.

Surface layer formation in polymer melts and in solutions

In this work the surface layer formation in polymer melts and in polymer solutions have been investigated with the atomic force microscope (AFM). In polymer melts, the formation of an immobile surface layer results in a steric repulsion, which can be measured by the AFM. From former work it is know, that polydimethyl siloxane (PDMS) forms a stable surface layer for molecular weights above 12 kDa. In the present thesis, polyisoprene (PI) was investigated apart from PDMS, by a)measuring the steric surface interactions and b)measuring the surface slip in hydrodynamic experiments. If a polymer flows over a surface, the flow velocity at the surface is larger then zero. If case of a surface layer formation the flow plane changes to the top of the adsorbed layer and the surface slip is reduced to zero. By measuring the surface slip in hydrodynamic experiments, it is therefore possible to determine the presence of a stable surface layer. The results show no stable repulsion for PI and only a small decrease of the surface slip. This indicates that PI does not form a stable surface layer, but is only adsorbed weakly to the surface. Furthermore for 8 kDa PDMS the timescale of the formation of a surface layer was investigated by changing themaximal force the tip applied to the surface. With a repulsive force present, applying a higher force than 15 nN resulted in a destruction of the surface layer, indicated by attractive forces. Reducing the applied force below 15 nN then resulted in an increase of the repulsion to the former state during one minute, thus indicating that a surface layer can be formed within one minute even under the influence of continuous measurements. As a next step, mixtures of two PDMS homopolymers with different chain lengths have been investigated. The aim was to verify theoretical predictions that shorter chains should predominate at the surface due to their smaller loss in conformational entropy. The measurements where done in dependence of the volume fractions of short and long chain PMDS. The results confirmed the short chain dominance for all mixtures with less then 90 vol.% long chain PDMS. Surface layer formation in solution was investigated for superplasticizers which are industrially used as an additive to cement. They change the surface interaction between the cement grains from attractive to repulsive and the freshlymixed cement paste therefore becomes liquid. The aimin this part of the thesis was, to investigate cement particle interactions in a close to real environment. Therefore calcium silicate hydrate phases have been precipitated onto an AFM tip and onto a calcite crystal and the interaction between these surfaces have beenmeasured with and without addition of superplasticizers. The measurements confirmed the change from attraction to repulsion upon addition of superplasticizers. The repulsive steric interaction increased with the length of the sidechain of the superplasticizer, and the dependence of the range of the steric interactions on the sidechain length indicated that the sidechains are in a coiled conformation.

Quantum hydrodynamic models for semiconductors with and without collisions

In dieser Arbeit werden Quantum-Hydrodynamische (QHD) Modelle betrachtet, die ihren Einsatz besonders in der Modellierung von Halbleiterbauteilen finden. Das QHD Modell besteht aus den Erhaltungsgleichungen für die Teilchendichte, das Momentum und die Energiedichte, inklusive der Quanten-Korrekturen durch das Bohmsche Potential. Zu Beginn wird eine Übersicht über die bekannten Ergebnisse der QHD Modelle unter Vernachlässigung von Kollisionseffekten gegeben, die aus ein­em Schrödinger-System für den gemischten-Zustand oder aus der Wigner-Glei­chung hergeleitet werden können. Nach der Reformulierung der eindimensionalen QHD Gleichungen mit linearem Potential als stationäre Schrö­din­ger-Gleichung werden die semianalytischen Fassungen der QHD Gleichungen für die Gleichspannungs-Kurve betrachtet. Weiterhin werden die viskosen Stabilisierungen des QHD Modells be­rück­sich­tigt, sowie die von Gardner vorgeschlagene numerische Viskosität für das {sf upwind} Finite-Differenzen Schema berechnet. Im Weiteren wird das viskose QHD Modell aus der Wigner-Glei­chung mit Fokker-Planck Kollisions-Ope­ra­tor hergeleitet. Dieses Modell enthält die physikalische Viskosität, die durch den Kollision-Operator eingeführt wird. Die Existenz der Lösungen (mit strikt positiver Teilchendichte) für das isotherme, stationäre, eindimensionale, viskose Modell für allgemeine Daten und nichthomogene Randbedingungen wird gezeigt. Die dafür notwendigen Abschätzungen hängen von der Viskosität ab und erlauben daher den Grenzübergang zum nicht-viskosen Fall nicht. Numerische Simulationen der Resonanz-Tunneldiode modelliert mit dem nichtisothermen, stationären, eindimensionalen, viskosen QHD Modell zeigen den Einfluss der Viskosität auf die Lösung. Unter Verwendung des von Degond und Ringhofer entwickelten Quanten-Entropie-Minimierungs-Verfahren werden die allgemeinen QHD-Gleichungen aus der Wigner-Boltzmann-Gleichung mit dem BGK-Kollisions-Operator hergeleitet. Die Herleitung basiert auf der vorsichtige Entwicklung des Quanten-Max­well­ians in Potenzen der skalierten Plankschen Konstante. Das so erhaltene Modell enthält auch vertex-Terme und dispersive Terme für die Ge­schwin­dig­keit. Dadurch bleibt die Gleichspannungs-Kurve für die Re­so­nanz-Tunnel­diode unter Verwendung des allgemeinen QHD Modells in einer Dimension numerisch erhalten. Die Ergebnisse zeigen, dass der dispersive Ge­schwin­dig­keits-Term die Lösung des Systems stabilisiert.

Thermal Casimir effect at fluid interfaces: Fluctuation-induced forces between anisotropic colloids

We study the effective interaction between two ellipsoidal particles at the interface of two fluid phases which are mediated by thermal fluctuations of the interface. Within a coarse-grained picture, the properties of fluid interfaces are very well described by an effective capillary wave Hamiltonian which governs both the equilibrium interface configuration and the thermal fluctuations (capillary waves) around this equilibrium (or mean-field) position. As postulated by the Goldstone theorem the capillary waves are long-range correlated. The interface breaks the continuous translational symmetry of the system, and in the limit of vanishing external fields - like gravity - it has to be accompanied by easily excitable long wavelength (Goldstone) modes – precisely the capillary waves. In this system the restriction of the long-ranged interface fluctuations by particles gives rise to fluctuation-induced forces which are equivalent to interactions of Casimir type and which are anisotropic in the interface plane. Since the position and the orientation of the colloids with respect to the interface normal may also fluctuate, this system is an example for the Casimir effect with fluctuating boundary conditions. In the approach taken here, the Casimir interaction is rewritten as the interaction between fluctuating multipole moments of an auxiliary charge density-like field defined on the area enclosed by the contact lines. These fluctuations are coupled to fluctuations of multipole moments of the contact line position (due to the possible position and orientational fluctuations of the colloids). We obtain explicit expressions for the behavior of the Casimir interaction at large distances for arbitrary ellipsoid aspect ratios. If colloid fluctuations are suppressed, the Casimir interaction at large distances is isotropic, attractive and long ranged (double-logarithmic in the distance). If, however, colloid fluctuations are included, the Casimir interaction at large distances changes to a power law in the inverse distance and becomes anisotropic. The leading power is 4 if only vertical fluctuations of the colloid center are allowed, and it becomes 8 if also orientational fluctuations are included.

Numerical studies of colloidal systems

The interplay of hydrodynamic and electrostatic forces is of great importance for the understanding of colloidal dispersions. Theoretical descriptions are often based on the so called standard electrokinetic model. This Mean Field approach combines the Stokes equation for the hydrodynamic flow field, the Poisson equation for electrostatics and a continuity equation describing the evolution of the ion concentration fields. In the first part of this thesis a new lattice method is presented in order to efficiently solve the set of non-linear equations for a charge-stabilized colloidal dispersion in the presence of an external electric field. Within this framework, the research is mainly focused on the calculation of the electrophoretic mobility. Since this transport coefficient is independent of the electric field only for small driving, the algorithm is based upon a linearization of the governing equations. The zeroth order is the well known Poisson-Boltzmann theory and the first order is a coupled set of linear equations. Furthermore, this set of equations is divided into several subproblems. A specialized solver for each subproblem is developed, and various tests and applications are discussed for every particular method. Finally, all solvers are combined in an iterative procedure and applied to several interesting questions, for example, the effect of the screening mechanism on the electrophoretic mobility or the charge dependence of the field-induced dipole moment and ion clouds surrounding a weakly charged sphere. In the second part a quantitative data analysis method is developed for a new experimental approach, known as "Total Internal Reflection Fluorescence Cross-Correlation Spectroscopy" (TIR-FCCS). The TIR-FCCS setup is an optical method using fluorescent colloidal particles to analyze the flow field close to a solid-fluid interface. The interpretation of the experimental results requires a theoretical model, which is usually the solution of a convection-diffusion equation. Since an analytic solution is not available due to the form of the flow field and the boundary conditions, an alternative numerical approach is presented. It is based on stochastic methods, i. e. a combination of a Brownian Dynamics algorithm and Monte Carlo techniques. Finally, experimental measurements for a hydrophilic surface are analyzed using this new numerical approach.

Rupture forces of split aptamers

Rupture forces of ligand-receptor interactions, such as proteins-proteins, proteins-cells, and cells-tissues, have been successfully measured by atomic force spectroscopy (AFS). For these measurements, the ligands and receptors were chemically modified so that they can be immobilized on the tip and on a substrate, respectively. The ligand interact the receptor when the tip approaches the substrate. This interaction can be studied by measuring rupture force upon retraction. However, this technique is not feasible for measurements involving small molecules, since they form only few H-bonds with their corresponding receptors. Modifying small molecules for immobilization on surfaces may block or change binding sites. Thus, recorded rupture forces might not reflect the full scope of the involved small ligand-receptor interactions.rnIn my thesis, a novel concept that allows measuring the rupture force of small involved ligand-receptor interactions and does not require molecular modification for immobilization was introduced. The rupture force of small ligand-receptor interaction is not directly measured but it can be determined from measurements in the presence and in the absence of the ligand. As a model system, the adenosine mono phosphate (AMP) and the aptamer that binds AMP were selected. The aptamer (receptor) is a single stranded DNA that can partially self-hybridize and form binding pockets for AMP molecules (ligands). The bonds between AMP and aptamer are provided by several H-bonds and pair stacking.rnIn the novel concept, the aptamer was split into two parts (oligo a and oligo b). One part was immobilized on the tip and the other one on the substrate. Approaching the tip to the substrate, oligo a and oligo b partially hybridized and the binding pockets were formed. After adding AMP into the buffer solution, the AMP bound in the pockets and additional H-bonds were formed. Upon retraction of the tip, the rupture force of the AMP-split aptamer complex was measured. In the presence of excess AMP, the rupture force increased by about 10 pN. rnThe dissociation constant of the AMP-split aptamer complex was measured on a single molecular level (~ 4 µM) by varying the AMP concentrations and measuring the rupture force at each concentration. Furthermore, the rupture force was amplified when more pockets were added to the split aptamer. rnIn the absence of AMP, the thermal off-rate was slightly reduced compared to that in the presence of AMP, indicating that the AMP stabilized the aptamer. The rupture forces at different loading rates did not follow the logarithmic fit which was usually used to describe the dependence of rupture forces at different loading rates of oligonucleotides. Two distinguished regimes at low and high loading rates were obtained. The two regimes were explained by a model in which the oligos located at the pockets were stretched at high loading rates. rnThe contribution of a single H-bond formed between the AMP molecule and the split aptamer was measured by reducing the binding groups of the AMP. The rupture forces reduce corresponding to the reduction of the binding groups. The phosphate group played the most important role in the formation of H-bond network between the AMP molecule and the split aptamer. rn

The effect of elastic modulus and equilibrium vapor adsorption on capillary forces

Die Kapillarkraft entsteht durch die Bildung eines Meniskus zwischen zwei Festkörpen. In dieser Doktorarbeit wurden die Auswirkungen von elastischer Verformung und Flϋssigkeitadsorption auf die Kapillarkraft sowohl theoretisch als auch experimentell untersucht. Unter Verwendung eines Rasterkraftmikroskops wurde die Kapillarkraft zwischen eines Siliziumoxid Kolloids von 2 µm Radius und eine weiche Oberfläche wie n.a. Polydimethylsiloxan oder Polyisopren, unter normalen Umgebungsbedingungen sowie in variierende Ethanoldampfdrϋcken gemessen. Diese Ergebnisse wurden mit den Kapillarkräften verglichen, die auf einem harten Substrat (Silizium-Wafer) unter denselben Bedingungen gemessen wurden. Wir beobachteten eine monotone Abnahme der Kapillarkraft mit zunehmendem Ethanoldampfdruck (P) fϋr P/Psat > 0,2, wobei Psat der Sättigungsdampfdruck ist.rnUm die experimentellen Ergebnisse zu erklären, wurde ein zuvor entwickeltes analytisches Modell (Soft Matter 2010, 6, 3930) erweitert, um die Ethanoladsorption zu berϋcksichtigen. Dieses neue analytische Modell zeigte zwei verschiedene Abhängigkeiten der Kapillarkraft von P/Psat auf harten und weichen Oberflächen. Fϋr die harte Oberfläche des Siliziumwafers wird die Abhängigkeit der Kapillarkraft vom Dampfdruck vom Verhältnis der Dicke der adsorbierten Ethanolschicht zum Meniskusradius bestimmt. Auf weichen Polymeroberflächen hingegen hängt die Kapillarkraft von der Oberflächenverformung und des Laplace-Drucks innerhalb des Meniskus ab. Eine Abnahme der Kapillarkraft mit zunehmendem Ethanoldampfdruck hat demnach eine Abnahme des Laplace-Drucks mit zunehmendem Meniskusradius zur folge. rnDie analytischen Berechnungen, fϋr die eine Hertzsche Kontakt-deformation angenommen wurde, wurden mit Finit Element Methode Simulationen verglichen, welche die reale Deformation des elastischen Substrats in der Nähe des Meniskuses explizit berϋcksichtigen. Diese zusätzliche nach oben gerichtete oberflächenverformung im Bereich des Meniskus fϋhrt zu einer weiteren Erhöhung der Kapillarkraft, insbesondere fϋr weiche Oberflächen mit Elastizitätsmodulen < 100 MPa.rn

Development of a system measuring adhesion forces in powder collectives

Fine powders commonly have poor flowability and dispersibility due to interparticle adhesion that leads to formation of agglomerates. Knowing about adhesion in particle collectives is indispensable to gain a deeper fundamental understanding of particle behavior in powders. Especially in pharmaceutical industry a control of adhesion forces in powders is mandatory to improve the performance of inhalation products. Typically the size of inhalable particles is in the range of 1 - 5 µm. In this thesis, a new method was developed to measure adhesion forces of particles as an alternative to the established colloidal probe and centrifuge technique, which are both experimentally demanding, time consuming and of limited practical applicability. The new method is based on detachment of individual particles from a surface due to their inertia. The required acceleration in the order of 500 000 g is provided by a Hopkinson bar shock excitation system and measured via laser vibrometry. Particle detachment events are detected on-line by optical video microscopy. Subsequent automated data evaluation allows obtaining a statistical distribution of particle adhesion forces. To validate the new method, adhesion forces for ensembles of single polystyrene and silica microspheres on a polystyrene coated steel surface were measured under ambient conditions. It was possible to investigate more than 150 individual particles in one experiment and obtain adhesion values of particles in a diameter range of 3 - 13 µm. This enables a statistical evaluation while measuring effort and time are considerably lower compared to the established techniques. Measured adhesion forces of smaller particles agreed well with values from colloidal probe measurements and theoretical predictions. However, for the larger particles a stronger increase of adhesion with diameter was observed. This discrepancy might be induced by surface roughness and heterogeneity that influence small and large particles differently. By measuring adhesion forces of corrugated dextran particles with sizes down to 2 µm it was demonstrated that the Hopkinson bar method can be used to characterize more complex sample systems as well. Thus, the new device will be applicable to study a broad variety of different particle-surface combinations on a routine basis, including strongly cohesive powders like pharmaceutical drugs for inhalation.

Wet bioadhesion in tree frogs

Viele Tiere wie etwa Geckos oder Laubfrösche können mittels ihrer Haftscheiben an Oberflächen kleben. Diese Haftscheiben ermöglichen es den Tieren, sich während ihrerrnFortbewegung an Oberflächen anzuheften und wieder zu lösen unabhängig von denrnvorherrschenden Umweltbedingungen. Frösche besitzen mikro- und nanostrukturierternsowie charakteristisch geformte Haftscheiben an Finger- und Zehenenden. Ihre besonderernevolutionäre Errungenschaft, sich stark und zugleich reversibel in sowohl trockenen alsrnauch feuchten Umgebungen anzuhaften, hat die Wissenschaft zur Nachahmung und Untersuchungrndieser Strukturen inspiriert. Zum besseren Verständnis der Mechanismen vonrnAnhaftung und Loslösung bei Laubfröschen wurden weiche, elastische und mikrostrukturierternOberflächen hergestellt, indem PDMS (Polydimethylsiloxan) auf einer Siliziummaskernmit Hexagonstruktur aufgetragen und vernetzt wurde. Dadurch wurden Anordnungenrnvon hexagonalen Mikrosäulen mit spezifischen geometrischen Eigenschaften undrnunterschiedlichen Kontaktgeometrien (normale, flache Form, T-Form und konkave Formrnder Säulenenden) erhalten. Um den Einfluss der van-der-Waals, hydrodynamischen,rnKapillar-und Adhäsionskräfte zu verstehen, wurden verschiedene experimentelle Ansätzernverfolgt: Die auf eine einzelne Säule wirkenden Adhäsionskräfte wurden mittelsrnRasterkraftmikroskopie gemessen. Dazu wurden speziell hergestellte kolloidale Sensorenrnverwendet. Diese Experimente wurden sowohl mit als auch ohne Flüssigkeitsfilm auf derrnSäule durchgeführt. Die Ergebnisse zeigten den Beitrag von Kapillarkraft und direktenrnKontaktkräften zur Adhäsionskraft bei Vorliegen eines Flüssigkeitsfilms. Die Adhäsionrnfiel umso größer aus, je weniger Flüssigkeit zwischen Sensor und Säule vorhanden war.rnIm Falle einer trockenen Adhäsion zeigte die Säule mit T-Form die höchste Adhäsion. Darndie Haftscheiben der Laubfrösche weich sind, können sie dynamisch ihre Form ändern,rnwas zu einer Änderung der hydrodynamischen Kraft zwischen Scheibe und Oberflächernführt. Der Einfluss der Oberflächenverformbarkeit auf die hydrodynamische Kraft wurderndaher am Modellsystem einer Kugel untersucht, welche sich einer weichen und ebenenrnOberfläche annähert. Dieses System wurde sowohl theoretisch über die Simulation finiterrnElemente als auch experimentell über die Messung mit kolloidalen Sonden untersucht.rnSowohl experimentelle Ergebnisse als auch die Simulationen ergaben eine Abnahme derrnhydrodynamischen Kraft bei Annäherung des kolloidalen Sensors an eine weiche undrnelastische Oberfläche. Beim Entfernen der Sensors von der Oberfläche verstärkte sichrndie hydrodynamische Anziehungskraft. Die Kraft, die zur Trennung eines Partikels von einer Oberfläche in Flüssigkeit notwendig ist, ist für weiche und elastischen Oberflächenrngrößer als für harte Oberflächen. In Bezug zur Bioadhäsion bei Laubfröschen konnternfestgestellt somit festgestellt werden, dass sich der hydrodynamische Anteil zur feuchtenrnBioadhäsion aufgrund der weichen Oberfläche erhöht. Weiterhin wurde der Einflussrndes Aspektverhältnisses der Säulen auf die Reibungskraft mittels eines kolloidalen Sensorsrnuntersucht. Gestreckte Säulen zeigten dabei eine höhere Reibung im Vergleich zu.rnSäulen mit einem gestreckten Hexagon als Querschnitt.

Hydrodynamics and solid dosage form disintegration/dissolution : immediate release tablets and novel in situ polyelectrolyte gastroretentive drug delivery systems

Solid oral dosage form disintegration in the human stomach is a highly complex process dependent on physicochemical properties of the stomach contents as well as on physical variables such as hydrodynamics and mechanical stress. Understanding the role of hydrodynamics and forces in disintegration of oral solid dosage forms can help to improve in vitro disintegration testing and the predictive power of the in vitro test. The aim of this work was to obtain a deep understanding of the influence of changing hydrodynamic conditions on solid oral dosage form performance. Therefore, the hydrodynamic conditions and forces present in the compendial PhEur/USP disintegration test device were characterized using a computational fluid dynamics (CFD) approach. Furthermore, a modified device was developed and the hydrodynamic conditions present were simulated using CFD. This modified device was applied in two case studies comprising immediate release (IR) tablets and gastroretentive drug delivery systems (GRDDS). Due to the description of movement provided in the PhEur, the movement velocity of the basket-rack assembly follows a sinusoidal profile. Therefore, hydrodynamic conditions are changing continually throughout the movement cycle. CFD simulations revealed that the dosage form is exposed to a wide range of fluid velocities and shear forces during the test. The hydrodynamic conditions in the compendial device are highly variable and cannot be controlled. A new, modified disintegration test device based on computerized numerical control (CNC) technique was developed. The modified device can be moved in all three dimensions and radial movement is also possible. Simple and complex moving profiles can be developed and the influence of the hydrodynamic conditions on oral solid dosage form performance can be evaluated. Furthermore, a modified basket was designed that allows two-sided fluid flow. CFD simulations of the hydrodynamics and forces in the modified device revealed significant differences in the fluid flow field and forces when compared to the compendial device. Due to the CNC technique moving velocity and direction are arbitrary and hydrodynamics become controllable. The modified disintegration test device was utilized to examine the influence of moving velocity on disintegration times of IR tablets. Insights into the influence of moving speed, medium viscosity and basket design on disintegration times were obtained. An exponential relationship between moving velocity of the modified basket and disintegration times was established in simulated gastric fluid. The same relationship was found between the disintegration times and the CFD predicted average shear stress on the tablet surface. Furthermore, a GRDDS was developed based on the approach of an in situ polyelectrolyte complex (PEC). Different complexes composed of different grades of chitosan and carrageenan and different ratios of those were investigated for their swelling behavior, mechanical stability, and in vitro drug release. With an optimized formulation the influence of changing hydrodynamic conditions on the swelling behavior and the drug release profile was demonstrated using the modified disintegration test device. Both, swelling behavior and drug release, were largely dependent on the hydrodynamic conditions. Concluding, it has been shown within this thesis that the application of the modified disintegration test device allows for detailed insights into the influence of hydrodynamic conditions on solid oral dosage form disintegration and dissolution. By the application of appropriate test conditions, the predictive power of in vitro disintegration testing can be improved using the modified disintegration test device. Furthermore, CFD has proven a powerful tool to examine the hydrodynamics and forces in the compendial as well as in the modified disintegration test device. rn