On modified simple reacting spheres kinetic model for chemically reactive gases

Versão dos autores para esta publicação.

Effect of mixing, particle interaction, and spatial dimension on the glass transition

In this thesis, the influence of composition changes on the glass transition behavior of binary liquids in two and three spatial dimensions (2D/3D) is studied in the framework of mode-coupling theory (MCT).The well-established MCT equations are generalized to isotropic and homogeneous multicomponent liquids in arbitrary spatial dimensions. Furthermore, a new method is introduced which allows a fast and precise determination of special properties of glass transition lines. The new equations are then applied to the following model systems: binary mixtures of hard disks/spheres in 2D/3D, binary mixtures of dipolar point particles in 2D, and binary mixtures of dipolar hard disks in 2D. Some general features of the glass transition lines are also discussed. The direct comparison of the binary hard disk/sphere models in 2D/3D shows similar qualitative behavior. Particularly, for binary mixtures of hard disks in 2D the same four so-called mixing effects are identified as have been found before by Götze and Voigtmann for binary hard spheres in 3D [Phys. Rev. E 67, 021502 (2003)]. For instance, depending on the size disparity, adding a second component to a one-component liquid may lead to a stabilization of either the liquid or the glassy state. The MCT results for the 2D system are on a qualitative level in agreement with available computer simulation data. Furthermore, the glass transition diagram found for binary hard disks in 2D strongly resembles the corresponding random close packing diagram. Concerning dipolar systems, it is demonstrated that the experimental system of König et al. [Eur. Phys. J. E 18, 287 (2005)] is well described by binary point dipoles in 2D through a comparison between the experimental partial structure factors and those from computer simulations. For such mixtures of point particles it is demonstrated that MCT predicts always a plasticization effect, i.e. a stabilization of the liquid state due to mixing, in contrast to binary hard disks in 2D or binary hard spheres in 3D. It is demonstrated that the predicted plasticization effect is in qualitative agreement with experimental results. Finally, a glass transition diagram for binary mixtures of dipolar hard disks in 2D is calculated. These results demonstrate that at higher packing fractions there is a competition between the mixing effects occurring for binary hard disks in 2D and those for binary point dipoles in 2D.

Close packed transitions in slit-shaped pores: Density functional theory study of methane adsorption capacity in carbon

An important feature of improving lattice gas models and classical isotherms is the incorporation of a pore size dependent capacity, which has hitherto been overlooked. In this paper, we develop a model for predicting the temperature dependent variation in capacity with pore size. The model is based on the analysis of a lattice gas model using a density functional theory approach at the close packed limit. Fluid-fluid and solid-fluid interactions are modeled by the Lennard-Jones 12-6 potential and Steele's 10-4-3, potential respectively. The capacity of methane in a slit-shaped carbon pore is calculated from the characteristic parameters of the unit cell, which are extracted by minimizing the grand potential of the unit cell. The capacities predicted by the proposed model are in good agreement with those obtained from grand canonical Monte Carlo simulation, for pores that can accommodate up to three adsorbed layers. Single particle and pair distributions exhibit characteristic features that correspond to the sequence of buckling and rhombic transitions that occur as the slit pore width is increased. The model provides a useful tool to model continuous variation in the microstructure of an adsorbed phase, namely buckling and rhombic transitions, with increasing pore width. (C) 2002 American Institute of Physics.

Modelo de escoamento multifásico para estudo da interação onda/sedimento

Modelos de escoamento multifásico são amplamente usados em diversas áreas de pesquisa ambiental, como leitos fluidizados, dispersão de gás em líquidos e vários outros processos que englobam mais de uma propriedade físico-química do meio. Dessa forma, um modelo multifásico foi desenvolvido e adaptado para o estudo do transporte de sedimentos de fundo devido à ação de ondas de gravidade. Neste trabalho, foi elaborado o acoplamento multifásico de um modelo euleriano não-linear de ondas do tipo Boussinesq, baseado na formulação numérica encontrada em Wei et al. (1995), com um modelo lagrangiano de partículas, fundamentado pelo princípio Newtoniano do movimento com o esquema de colisões do tipo esferas rígidas. O modelo de ondas foi testado quanto à sua fonte geradora, representada por uma função gaussiana, pá-pistão e pá-batedor, e quanto à sua interação com a profundidade, através da não-linearidade e de propriedades dispersivas. Nos testes realizados da fonte geradora, foi observado que a fonte gaussiana, conforme Wei et al. (1999), apresentou melhor consistência e estabilidade na geração das ondas, quando comparada à teoria linear para um kh   . A não-linearidade do modelo de ondas de 2ª ordem para a dispersão apresentou resultados satisfatórios quando confrontados com o experimento de ondas sobre um obstáculo trapezoidal, onde a deformação da onda sobre a estrutura submersa está em concordância com os dados experimentais encontrados na literatura. A partir daí, o modelo granular também foi testado em dois experimentos. O primeiro simula uma quebra de barragem em um tanque contendo água e o segundo, a quebra de barragem é simulada com um obstáculo rígido adicionado ao centro do tanque. Nesses experimentos, o algoritmo de colisão foi eficaz no tratamento da interação entre partícula-partícula e partícula-parede, permitindo a evidência de processos físicos que são complicados de serem simulados por modelos de malhas regulares. Para o acoplamento do modelo de ondas e de sedimentos, o algoritmo foi testado com base de dados da literatura quanto à morfologia do leito. Os resultados foram confrontados com dados analíticos e de modelos numéricos, e se mostraram satisfatórios com relação aos pontos de erosão, de sedimentação e na alteração da forma da barra arenosa

Liquid-vapor interfaces of patchy colloids

We investigate the liquid-vapor interface of a model of patchy colloids. This model consists of hard spheres decorated with short-ranged attractive sites ("patches") of different types on their surfaces. We focus on a one-component fluid with two patches of type A and nine patches of type B (2A9B colloids), which has been found to exhibit reentrant liquid-vapor coexistence curves and very low-density liquid phases. We have used the density-functional theory form of Wertheim's first-order perturbation theory of association, as implemented by Yu and Wu [J. Chem. Phys. 116, 7094 (2002)], to calculate the surface tension, and the density and degree of association profiles, at the liquid-vapor interface of our model. In reentrant systems, where AB bonds dominate, an unusual thickening of the interface is observed at low temperatures. Furthermore, the surface tension versus temperature curve reaches a maximum, in agreement with Bernardino and Telo da Gama's mesoscopic Landau-Safran theory [Phys. Rev. Lett. 109, 116103 (2012)]. If BB attractions are also present, competition between AB and BB bonds gradually restores the monotonic temperature dependence of the surface tension. Lastly, the interface is "hairy," i.e., it contains a region where the average chain length is close to that in the bulk liquid, but where the density is that of the vapor. Sufficiently strong BB attractions remove these features, and the system reverts to the behavior seen in atomic fluids.

Viscosity measurements of compressed liquid dipropyl and dibutyl adipates

The paper reports viscosity measurements of compressed liquid dipropyl (DPA) and dibutyl (DBA) adipates obtained with two vibrating wire sensors developed in our group. The vibrating wire instruments were operated in the forced oscillation, or steady-state mode. The viscosity measurements of DPA were carried out in a range of pressures up to 18. MPa and temperatures from (303 to 333). K, and DBA up to 65. MPa and temperature from (303 to 373). K, covering a total range of viscosities from (1.3 to 8.3). mPa. s. The required density data of the liquid samples were obtained in our laboratory using an Anton Paar vibrating tube densimeter and were reported in a previous paper. The viscosity results were correlated with density, using a modified hard-spheres scheme. The root mean square deviation of the data from the correlation is less than (0.21 and 0.32)% and the maximum absolute relative deviations are within (0.43 and 0.81)%, for DPA and DBA respectively. No data for the viscosity of both adipates could be found in the literature. Independent viscosity measurements were also performed, at atmospheric pressure, using an Ubbelohde capillary in order to compare with the vibrating wire results. The expanded uncertainty of these results is estimated as ±1.5% at a 95% confidence level. The two data sets agree within the uncertainty of both methods. © 2015 Published by Elsevier B.V.

A simple analysis of thermodynamic properties for classical plasmas: I. theory

By eliminating the short range negative divergence of the Debye–Hückel pair distribution function, but retaining the exponential charge screening known to operate at large interparticle separation, the thermodynamic properties of one-component plasmas of point ions or charged hard spheres can be well represented even in the strong coupling regime. Predicted electrostatic free energies agree within 5% of simulation data for typical Coulomb interactions up to a factor of 10 times the average kinetic energy. Here, this idea is extended to the general case of a uniform ionic mixture, comprising an arbitrary number of components, embedded in a rigid neutralizing background. The new theory is implemented in two ways: (i) by an unambiguous iterative algorithm that requires numerical methods and breaks the symmetry of cross correlation functions; and (ii) by invoking generalized matrix inverses that maintain symmetry and yield completely analytic solutions, but which are not uniquely determined. The extreme computational simplicity of the theory is attractive when considering applications to complex inhomogeneous fluids of charged particles.

Complex formation in solutions of oppositely charged polyelectrolytes at different polyion compositions and salt content

The formation of complexes in solutions containing positively charged polyions (polycations) and a variable amount of negatively charged polyions (polyanions) has been investigated by Monte Carlo simulations. The polyions were described as flexible chains of charged hard spheres interacting through a screened Coulomb potential. The systems were analyzed in terms of cluster compositions, structure factors, and radial distribution functions. At 50% charge equivalence or less, complexes involving two polycations and one polyanion were frequent, while closer to charge equivalence, larger clusters were formed. Small and neutral complexes dominated the solution at charge equivalence in a monodisperse system, while larger clusters again dominated the solution when the polyions were made polydisperse. The cluster composition and solution structure were also examined as functions of added salt by varying the electrostatic screening length. The observed formation of clusters could be rationalized by a few simple rules.

A Monte Carlo study of solutions of oppositely charged polyelectrolytes

The formation of complexes appearing in solutions containing oppositely charged polyelectrolytes has been investigated by Monte Carlo simulations using two different models. The polyions are described as flexible chains of 20 connected charged hard spheres immersed in a homogenous dielectric background representing water. The small ions are either explicitly included or their effect described by using a screened Coulomb potential. The simulated solutions contained 10 positively charged polyions with 0, 2, or 5 negatively charged polyions and the respective counterions. Two different linear charge densities were considered, and structure factors, radial distribution functions, and polyion extensions were determined. A redistribution of positively charged polyions involving strong complexes formed between the oppositely charged polyions appeared as the number of negatively charged polyions was increased. The nature of the complexes was found to depend on the linear charge density of the chains. The simplified model involving the screened Coulomb potential gave qualitatively similar results as the model with explicit small ions. Finally, owing to the complex formation, the sampling in configurational space is nontrivial, and the efficiency of different trial moves was examined.

Structure of redispersible SnO2 nanoparticles

Suspensions of undoped SnO2 nanoparticles and containing Eu3+ ions were prepared by a sol-gel procedure. Using the classical synthesis method ( precipitation), the particles tend to grow by a coarsening process in order to minimize the surface free energy. This effect can strongly be reduced by the addition of an amide and surfactant during the synthesis, which decreases the surface free energy of the colloidal particles. These additives promote the formation of powders composed of very small primary particles formed by a crystallite of 10 Angstrom, and exhibit good redispersion properties. The local and long order structures of the redispersible powder were studied by X-rays absorption spectroscopy at Sn L-I edge and X-rays diffraction, respectively. The structure of the colloidal aggregates in suspension was investigated by small angle X-rays scattering (SAXS). SAXS results indicate the sol are composed by a polidisperse system of hard spheres resulting of agglomeration of the primary particles and their size increasing by agglomeration for progressively higher Eu3+ content.

Estudo por simulação computacional da interação entre polianfóteros fracos e macroíons cilíndricos

Pós-graduação em Biofísica Molecular - IBILCE

Properties of a new poly-ether-glycol copolymer.

Triblock copolymers are made of monomer segments, being the central part usually hydrophobic and the outer parts hydrophilic. By varying sizes, molecular weights and monomer types of the segments one obtains different final molecules, with different physico-chemical properties, which are directly related to the performance of the final product. Looking for new products to be used, among other possibilities, in biological applications, a new polymer (Figure 1) was synthesized by the Dow Chemical and studied by Size Exclusion Chromatography, Fourier Transformed Infrared Spectrometry, Small-angle X-ray Scattering (SAXS) and its cloud point was determined by measuring light transmittance. The studies showed low molecular polydispersivety, but different polarities in the macromolecules fractions. Due to the low solubility of Diol in water, a mixture of water/butyl diglycol was used as solvent. An extensive analysis by SAXS was performed for concentrations from 50 wt% to 80 wt% of Diol in solution. Small concentrations showed very low signal to noise ratio, making it impossible to be analysed. The scattering intensity including the form factor of polydisperse non-homogeneous spheres, and the structure factor of interacting hard spheres was fitted to the curves. As the polymer concentration is high, the fitting of form factors of direct and reverse micelles were compared. The results for direct micelles were better up to 80 wt%, whereas at 90 wt% and 95 wt% the curves were better fitted by reverse micelles. It might seem odd that direct micelles are present up to such high concentrations, but it might have been caused by the presence of butyl diglycol, which increases the solubility of Diol in water. The inner and outer radius of the micelles, electron density distribution, and interaction radius of the micelles were obtained. The polydispersivety increases with Diol concentration. Besides, the interaction radius increases with solvent concentration, even when reversed micelles are present. In the last case, accompanied by an increase of inner radius (water content), as there are fewer Diol molecules to involve the water nuclei, which become larger, further apart, and in less number.

Physikalische Eigenschaften kolloidaler Festkörper

In dieser Arbeit wird das Phasenverhalten fluid-kristallin und kristallin-amorph, die elastischen Eigenschaften, das Nukleationsverhalten und das diffusive Verhalten ladungsstabilisierter Kolloide aus sphärischen Polystyrol- und Polytetrafluorethylenpartikeln in wässerigen Dispersionsmitteln bei sehr geringem Fremdionengehalt systematisch untersucht. Die dazugehörigen Messungen werden an einer neuartigen selbstkonstruierten Kombinationslichtstreuapparatur durchgeführt, die die Meßmethoden der dynamischen Lichtstreuung, statischen Lichtstreuung und Torsionsresonanzspektroskopie in sich vereint. Die drei Meßmethoden sind optimal auf die Untersuchung kolloidaler Festkörper abgestimmt. Das elastische Verhalten der Festkörper kann sehr gut durch die Elastizitätstheorie atomarer Kristallsysteme beschrieben werden, wenn ein Debye-Hückel-Potential im Sinne des Poisson-Boltzmann-Cell Modells als Wechselwirkungspotential verwendet wird. Die ermittelten Phasengrenzen fluid-kristallin stehen erstmalig in guter Übereinstimmung mit Ergebnissen aus molekulardynamischen Simulationen, wenn die in der Torsionsresonanzspektroskopie bestimmte Wechselwirkungsenergie zu Grunde gelegt wird. Neben der Gleichgewichtsstruktur sind Aussagen zur Verfestigungskinetik möglich. Das gefundene Nukleationserhalten kann gut durch die klassische Nukleationstheorie beschrieben werden, wenn bei niedriger Unterkühlung der Schmelze ein Untergrund heterogener Keimung berücksichtigt wird. PTFE-Partikel zeigen auch bei hohen Konzentrationen nur geringfügige Mehrfachstreuung. Durch ihren Einsatz ist erstmals eine systematische Untersuchung des Glasübergangs in hochgeladenen ladungsstabilisierten Systemen möglich. Ladungsstabilisierte Kolloide unterscheiden sich vor allem durch ihre extreme Kristallisationstendenz von früher untersuchten Hartkugelsystemen. Bei hohen Partikelkonzentrationen (Volumenbrüche größer 10 Prozent) kann ein glasartiger Festkörper identifiziert werden, dessen physikalisches Verhalten die Existenz eines Bernalglases nahe legt. Der Glasübergang ist im Vergleich mit den in anderen kolloidalen Systemen und atomaren Systemen beobachteten Übergängen von sehr unterschiedlichem Charakter. Im verwendeten PTFE-System ist auf Grund der langreichweitigen stark repulsiven Wechselwirkung kein direkter Zugang des Glaszustandes aus der übersättigten Schmelze möglich. Der amorphe Festkörper entsteht hier aus einer nanokristallinen Phase. Die Keimrate steigt im zugänglichen Meßbereich annähernd exponentiell mit der Partikelanzahldichte, so daß man feststellen kann, daß der Glaszustand nicht durch Unterdrückung der Nukleation, sondern durch eine Forcierung derselben erreicht wird.

Realisierung einer 3D-Kreuzkorrelationsanlage zur Untersuchung von Struktur und Dynamik hochkonzentrierter Kolloide

Realisierung einer 3D-Kreuzkorrelationsanlage zur Untersuchung von Struktur und Dynamik hochkonzentrierter Kolloide Im Rahmen dieser Arbeit wird eine neuartige 3D-Kreuzkorrelationsanlage zur mehrfachstreufreien Untersuchung des diffusiven Verhaltens hochkonzentrierter kolloidaler Suspensionen vorgestellt. Hierzu werden zwei Lichtstreuexperimente gleichzeitig am gleichen Streuvolumen und mit dem gleichen Streuvektor durchgeführt. Aus der so gewonnenen Kreuzkorrelationsfunktion kann das dynamische Verhalten der Kolloide bestimmt werden. Für die Diffusion der Partikel spielen neben der direkten Wechselwirkung elektroviskoser Effekt und die hydrodynamische Wechselwirkung eine entscheidende Rolle. Insbesondere bei hohen Konzentrationen kann keiner der drei Effekte vernachlässigt werden. Die zu messenden Unterschiede in den Diffusionskoeffizienten sind sehr klein. Daher wurde der experimentelle Aufbau detailliert charakterisiert. Hierbei konnten theoretische Überlegungen hinsichtlich des Nachpulsens und der Totzeit der verwendeten Si-Avalanche-Photodioden überprüft werden. Der Kurzzeitselbstdiffusionskoeffizient hochkonzentrierter geladener kolloidaler Suspensionen wurde gemessen. Um die Daten bei hohen Konzentrationen korrekt zu normieren, wurde der elektroviskose Effekt bei geringen Konzentrationen ausführlich untersucht. Hierbei zeigte sich, dass der elektroviskose Einzelteilcheneffekt zu einer monotonen Abnahme des Diffusionskoeffizienten bei abnehmender Ionenstärke führt. Anhand der volumenbruchabhängigen Daten des Kurzzeitselbstdiffusionskoeffizienten konnte zum ersten Mal gezeigt werden, dass die hydrodynamische Wechselwirkung einen geringeren Einfluss auf die Diffusion hat, falls das direkte Wechselwirkungspotential ein Coulomb-Potential anstelle eines Harte-Kugel-Potentials ist.

Untersuchungen zur Verfestigungskinetik in Suspensionen kolloidaler Partikel

Kolloidale Suspensionen, bei denen man die kolloidalen Teilchen als "Makroatome" in einem Kontinuum aus Lösungsmittelmolekülen auffaßt, stellen ein geeignetes Modellsystem zur Untersuchung von Verfestigungsvorgängen dar. Auf Grund der typischen beteiligten Längen- und Zeitskalen können Phasenübergänge bequem mit optischen Verfahren studiert werden. In der vorliegenden Arbeit wurde die Kinetik der Kristallisation in drei kolloidalen Systemen unterschiedlicher Teilchen-Teilchen-Wechselwirkung mit Lichtstreu- und mikroskopischen Methoden untersucht. Zur Untersuchung von Suspensionen aus sterisch stabilisierten PMMA-Teilchen, die in guter Näherung wie harte Kugeln wechselwirken, wurde ein neuartiges Laserlichtstreuexperiment aufgebaut, das die gleichzeitige Detektion von Bragg- und Kleinwinkelstreuung an einer Probe erlaubt. Damit konnte der zeitliche Verlauf der Kristallisation verfolgt sowie u.a. Nukleationsraten und erstmals auch Wachstumsgeschwindigkeiten bestimmt werden; diese wurden mit klassischer Nukleationstheorie sowie Wilson-Frenkel-Wachstum verglichen. In beiden Fällen konnte sehr gute Übereinstimmung mit der Theorie festgestellt werden. In Systemen geladener Partikel wurden mit Bragg-Mikroskopie die Wachstumsgeschwindigkeiten heterogener, an der Wand der Probenzelle aufwachsender Kristalle untersucht. Die Anpassung eines Wilson-Frenkel-Wachstumsgesetzes gelingt auch hier, wenn man die dazu eingeführte reskalierte Energiedichte auf den Schmelzpunkt bezieht. Geeignete Reskalierung der Daten erlaubt den Vergleich mit den Hartkugelsystemen. Zum ersten Mal wurde die Kristallisationskinetik in zwei verschiedenen kolloidalen binären Mischungen bestimmt und ausgewertet: In Beimischungen einer nichtkristallisierenden Teilchensorte zu einer kristallisierenden Suspension konnten die Daten mit einem modifizierten Wilson-Frenkel-Gesetz beschrieben werden, während in Mischungen aus zwei kristallisierenden Partikelsystemen eine unerwartet hohe Abnahme der Wachstumsgeschwindigkeiten beobachtet wurde. Kolloidale Suspensionen hartkugelähnlicher Mikrogel-Partikel konnten mit Hilfe des Lichtstreuaufbaues ebenfalls zum ersten Mal untersucht werden. Es wurde eine ähnliche Kristallisationskinetik wie in den PMMA-Systemen gefunden, jedoch auch einige wichtige Unterschiede, die insbesondere den Streumechanismus im Kleinwinkelbereich betrafen. Hier wurden verschiedene Interpretationsvorschläge diskutiert.