Source author record

Marcia C. Barbosa

Marcia C. Barbosa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2022arXiv

Density anomaly in a minimum model of water-alcohol mixtures

We study the change in the temperature of maximum density (TMD) of a water-like solvent when small amounts of solute are added to the mixture. The water is modeled as a two length scales potential, which is known to exhibit water-like characteristic anomalies, while the solute is chosen to have an attractive interaction with the solvent. If the solute is a structure maker the TMD increases with the addition of solute, while if it is a structure breaker the TMD decreases. Both in an exact one dimensional model and using molecular dynamics in three dimensions we show that the necessary ingredient to represent structure makers is a strong solute-solvent attraction.

preprint2022arXiv

Water inside charged nanoslits: Structure and Dielectric study with a novel water model FAB/$ε$

In this work, the dielectric behavior of water inside charged nanoslit of graphene is studied to analized the water molecules under electrical confinement; through polarizing the nanoslit of graphene, creating an electric field inside the nanopore. How the water molecules are structured under this type of electrical confinement is studied with two force fields of water, the three-site water models here used are the SCP/$ε$ and the FAB/$ε$, the first is a rigid model that improves the SPC model and the second is a flexible model that improves all the force fields of three sites non-polarizables and flexibles.

preprint2020arXiv

Salt parameterization can drastically affect the results from classical atomistic simulations of water desalination by MoS$_2$ nanopores

Water scarcity is a reality in our world, and scenarios predicted by leading scientists in this area indicate that it will worsen in the next decades. However, new technologies based in low-cost seawater desalination can prevent the worst scenarios, providing fresh water for humanity. With this goal, membranes based in nanoporous materials have been suggested in recent years. One of the materials suggested is MoS$_2$, and classical Molecular Dynamics (MD) simulation is one of the most powerful tools to explore these nanomaterials. However, distinct Force Fields employed in MD simulations are parameterized based on distinct experimental quantities. In this paper, we compare two models of salt that were build based on distinct properties of water-salt mixtures. One model fits the hydration free energy and lattice properties, the second fits the crystal density and the density and the dielectric constant of water and salt mixtures. To compare the models, MD simulations for salty water flow through two nanopores sizes were used -- one pore big enough to accommodate hydrated ions, and one smaller in which the ion has to dehydrate to enter, and two rigid water models from the TIP4P family -- the TIP4P/2005 and TIP4P/$ε$. Our results indicate that the water permeability and salt rejection by the membrane are more influenced by the salt model than by the water model, especially for the narrow pore. In fact, completely distinct mechanisms were observed, and they are related to the characteristics employed in the ion model parameterization. The results show that not only the water model can influence the outcomes, but the ion model plays a crucial role.

preprint2020arXiv

Water diffusion in carbon nanotubes: interplay between confinement, surface deformation and temperature

In this article we investigate through molecular dynamics simulations the diffusion behavior of the TIP4P/2005 water when confined in pristine and deformed carbon nanotubes (armchair and zigzag). To analyze different diffusive mechanisms, the water temperature was varied from $210\leq T\leq 380$~K. The results of our simulations reveal that water present a non-Arrhenius to Arrhenius diffusion crossover. The confinement shifts the diffusion transition to higher temperatures when compared with the bulk system. In addition, for narrower nanotubes, water diffuses in a single line which leads to a mobility independent of the activation energy.

preprint2020arXiv

Water diffusion in rough carbon nanotubes

We use molecular dynamics simulations to study the diffusion of water inside deformed carbon nanotubes with different degrees of deformation at 300 K. We found that the number of hydrogen bonds that water forms depends on nanotube topology, leading to enhancement or suppression of water diffusion. The simulation results reveal that more realistic nanotubes should be considered to understand the confined water diffusion behavior, at least for the narrowest nanotubes, when the interaction between water molecules and carbon atoms is relevant.

preprint2016arXiv

Flow and structure of fluids in functionalized nanopores

We investigate through non-equilibrium Molecular Dynamics simulations the structure and flow of fluids in functionalized nanopores. The nanopores are modeled as cylindrical structures with solvophilic and solvophobic sites. Two fluids are modeled. The first is a standard Lennard Jones fluid. The second one is modeled with a isotropic two-length scale potential, which exhibits in bulk water-like anomalies. Our results indicates distinct dependence of the overall mass flux for each species of fluid with the number of solvophilic sites for different nanotubes radii. Also, the density and fluid structure are dependent from the nanotube radius and the solvophilic properties of the nanotube. This indicates that the presence of a second length scale in the fluid-fluid interaction will lead to distinct behavior. Also, our results shows that chemically functionalized nanotubes with different radius will have distinct nano-fluidic features. Our results are explained in the basis of the characteristic scale fluid properties and the effects of nanoconfinement.

preprint2016arXiv

Potassium bromide, KBr/ε : New Force Field

The correct description of the ionic interaction and stable equilibrium of the simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model of the monovalent ions is very important. The force field needs to reproduce coincident many properties of ions, like their structure, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. Using a similar strategy employed in the parameterization of the NaCl/epsilon , in this paper, we first propose a force field for the Potassium Bromide, the KBr/epsilon. This new model is compared with the experimental values of cristal density and structure for the salt and the density, the viscosity, the dielectric constant and the solubility in the water solution for a range of concentrations. Next, the transferability, of this new model KBr/epsilon and the NaCl/epsilon, is verified by creating the KCl/epsilon and the NaBr/epsilon models. The strategy is to employ the same parameters obtained for the NaCl/epsilon and for the KBr/epsilon force fields. The two new models derived are also compared with the experimental values for the density, the viscosity, the dielectric constant and the solubility in the water solution for a range of concentrations.

preprint2015arXiv

Effects of nanopore and fluid structure on anomalies and phase transitions of confined core-softened fluids

We use Molecular Dynamics simulations to study how the nanopore and the fluid structures affects the dynamic, thermodynamic and structural properties of a confined anomalous fluid. The fluid is modeled using an effective pair potential derived from the ST4 atomistic model for water. This system exhibits density, structural and dynamical anomalies and the vapor-liquid and liquid-liquid critical points similar to the quantities observed in bulk water. The confinement is modeled both by smooth and structured walls. The temperatures of extremum density and diffusion for the confined fluid show a shift to lower values while the pressures move to higher amounts for both smooth and structured confinement. In the case of smooth walls, the critical points and the limit between fluid and amorphous phases show a non-monotonic change in the temperatures and pressures when the the nanopore size is increase. In the case of structured walls the pressures and temperatures of the critical points varies monotonicaly with the porous size. Our results are explained on basis of the competition between the different length scales of the fluid and the wall-fluid interaction.

preprint2015arXiv

Order-Disorder structural transition in a confined fluid

In this paper the amorphous/solid to disorder liquid structural phase transitions of an anomalous confined fluid is analyzed using their local fractal dimension. The model is a system of particles interacting through a two length scales potentials confined by two infinite plates. In the bulk, this fluid exhibit water-like anomalies and under confinement forms layers of particles. The particle distributions of them, present different arrangements related to amorphous/solid phases. Here only the contact layer is analyzed through fractal singularity spectrum. At high densities the structural transition its quantify by the order degree to determine the phases affected by the confinement. This mapping shows that the system as the temperature increased, the fractal dimension decreases, which is consistent with the behavior studying in such systems. This result suggests that under thermodynamic perturbations, an anomalous confined liquid, presents different phase transitions achieving be characterized by its fractality.

preprint2015arXiv

Self-Assembly and Waterlike Anomalies in Janus Nanoparticles

We explore the pressure versus temperature phase diagram of dimeric Janus nanoparticles using Molecular Dynamics simulations. The nanoparticle was modeled as a dumbbells particle, and have one monomer that interacts by a standard Lennard Jones potential and another monomer that is modeled using a two-length scale shoulder potential. Monomeric and dimeric systems modeled by this shoulder potential show waterlike anomalies, and we investigate if a Janus nanoparticle composed by one anomalous monomer will exhibit anomalous behavior and self-assembly structures. The influence of the non-anomalous monomer in the dimeric system properties was explored. We show that the diffusion anomaly is maintained, while the density anomaly can disappear depending on the non-anomalous monomer characteristics. As well, the self-assembled structures are affected. Our results are discussed in the basis of the distinct monomer-monomer interactions and on the two-length scale fluid characteristics.

preprint2015arXiv

Thermodynamic and dynamic anomalous behavior in the TIP4P/ε water model

The model Tip4p/ε for water is tested for the presence of thermodynamic and dy- namic anomalies. Molecular dynamic simulations for this model were performed and we show that for this system the density versus temperature at constant pressure ex- hibits a maximum. In addition we also show that the diffusion coefficient versus density at constant temperature has a maximum and a minimum. The anomalous behavior of the density and of the diffusion coefficient obey the water hierachy. The results for the Tip4p- are consistent with experiments and when compared with the Tip4p-2005 model show similar results a variety of physical properties and better performance for the dielectric constant.

preprint2015arXiv

Universality, Correlations, and Rankings in the Brazilian Universities National Admission Examinations

The scores obtained by students that have performed the ENEM exam, the Brazilian High School National Examination used to admit students at the Brazilian universities, is analyzed. The average high school's scores are compared between different disciplines through the Pearson correlation coefficient. The results show a very large correlation between the performance in the different subjects. Even thought the students' scores in the ENEM due to the standardization form a Gaussian, we show that the high schools' scores form a bimodal distribution that can not be used to evaluate and compare performance over time. We also show that this high schools distribution reflects the correlation between school performance and economic level of the students. The ENEM's scores are compared with a Brazilian non standardized exam, the entrance exam at the Universidade Federal do Rio Grande do Sul. The comparison of the performance of the same individuals in both tests is compared showing that the two tests not only select different abilities but chooses a different set of individuals. Our results indicates that standardized exams might be an interesting tool to compare performance over the years but only of individuals and not of institutions.

preprint2014arXiv

Phase Diagram and Thermodynamic and Dynamic Anomalies in a Pure Repulsive Model

Using Monte Carlo simulations a lattice gas model with only repulsive interactions was checked for the presence of anomalies. We show that this system exhibits the density (temperature of maximum density - TMD) and diffusion anomalies as present in liquid water. These anomalous behavior exist in the region of the chemical potential vs temperature phase diagram where two structured phases are present. A fragile-to-strong dynamic transition is also observed in the vicinity of the TMD line.

preprint2014arXiv

Surface Phase Transition in Anomalous Fluid in Nanoconfinement

We explore by molecular dynamic simulations the thermodynamical behavior of an anomalous fluid confined inside rigid and flexible nanopores. The fluid is modeled by a two length scale potential. In the bulk this system exhibits the density and diffusion anomalous behavior observed in liquid water. We show that the anomalous fluid confined inside rigid and flexible nanopores forms layers. As the volume of the nanopore is decreased the rigid surface exhibits three consecutive first order phase transitions associated with the change in the number of layers. These phase transitions are not present for flexible confinement. Our results indicate that the nature of confinement is relevant for the properties of the confined liquid what suggests that confinement in carbon nanotubes should be quite different from confinement in biological channels.

preprint2013arXiv

A water-like model under confinement for hydrophobic and hydrophilic particle-plate interaction potentials

Molecular dynamic simulations were employed to study a water-like model confined between hydrophobic and hydrophilic plates. The phase behavior of this system is obtained for different distances between the plates and particle-plate potentials. For both hydrophobic and hydrophilic walls there are the formation of layers. Crystallization occurs at lower temperature at the contact layer than at the middle layer. In addition, the melting temperature decreases as the plates become more hydrophobic. Similarly, the temperatures of maximum density and extremum diffusivity decrease with hydrophobicity.

preprint2013arXiv

Phase diagram and water-like anomalies in core-softened shoulder-dumbbells system

Using molecular dynamics we studied the role of the anisotropy on the phase boundary and on the anomalous behavior of 250 dimeric particles interacting by a core-softened potential. This study led us to an unexpected result: the introduction of a rather small anisotropy, quantified by the distance between the particles inside each dimer, leads to the increase of the size of the regions in the pressure-temperature phase diagram of anomalies when compared to the isotropic monomeric case. However, as the anisotropy increases beyond a threshold the anomalous regions shrinks. We found that this behavior can be understood by decoupling the translational and non-translational kinetic energy components into different translational and non-translational temperatures.

preprint2012arXiv

Connection Between Minimum of Solubility and Temperature of Maximum Density in an Associating Lattice Gas Model

In this paper we investigate the solubility of a hard - sphere gas in a solvent modeled as an associating lattice gas (ALG). The solution phase diagram for solute at 5% is compared with the phase diagram of the original solute free model. Model properties are investigated thr ough Monte Carlo simulations and a cluster approximation. The model solubility is computed via simulations and shown to exhibit a minimum as a function of temperature. The line of minimum solubility (TmS) coincides with the line of maximum density (TMD) for different solvent chemical potentials.

preprint2012arXiv

Diffusion Enhancement in Core-softened fluid confined in nanotubes

We study the effect of confinement in the dynamical behavior of a core-softened fluid. The fluid is modeled as a two length scales potential. This potential in the bulk reproduces the anomalous behavior observed in the density and in the diffusion of liquid water. A series of $NpT$ Molecular Dynamics simulations for this two length scales fluid confined in a nanotube were performed. We obtain that the diffusion coefficient increases with the increase of the nanotube radius for wide channels as expected for normal fluids. However, for narrow channels, the confinement shows an enhancement in the diffusion coefficient when the nanotube radius decreases. This behavior, observed for water, is explained in the framework of the two length scales potential.

preprint2012arXiv

Hydration and anomalous solubility of the Bell-Lavis model as solvent

We address the investigation of the solvation properties of the minimal orientational model for water, originally proposed by Bell and Lavis. The model presents two liquid phases separated by a critical line. The difference between the two phases is the presence of structure in the liquid of lower density, described through orientational order of particles. We have considered the effect of small inert solute on the solvent thermodynamic phases. Solute stabilizes the structure of solvent, by the organization of solvent particles around solute particles, at low temperatures. Thus, even at very high densities, the solution presents clusters of structured water particles surrounding solute inert particles, in a region in which pure solvent would be free of structure. Solute intercalates with solvent, a feature which has been suggested by experimental and atomistic simulation data. Examination of solute solubility has yielded a minimum in that property, which may be associated with the minimum found for noble gases. We have obtained a line of minimum solubility (TmS) across the phase diagram, accompanying the line of maximum in density (TMD). This coincidence is easily explained for non-interacting solute and it is in agreement with earlier results in the literature. We give a simple argument which suggests that interacting solute would dislocate TmS to higher temperatures.

preprint2012arXiv

Ion fluxes through nano-pores and transmembrane channels

We introduce an implicit solvent Molecular Dynamics approach for calculating ionic fluxes through narrow nano-pores and transmembrane channels. The method relies on a dual-control- volume grand-canonical molecular dynamics (DCV-GCMD) simulation and the analytical solution for the electrostatic potential inside a cylindrical nano-pore recently obtained by Levin [Europhys. Lett., 76, 163 (2006)]. The theory is used to calculate the ionic fluxes through an artificial trans-membrane c hannel which mimics the antibacterial gramicidin A channel. Both current-voltage and current-concentration relations are calculated under various experimental conditions. We show that our results are comparable to the characteristics associated to the gramicidin A pore, specially the existence of two binding sites inside the pore and the observed saturation in the current-concentration profiles.

preprint2011arXiv

Core-softened Fluids, Water-like Anomalies and the Liquid-Liquid Critical Points

Molecular dynamics simulations are used to examine the relationship between water-like anomalies and the liquid-liquid critical point in a family of model fluids with multi-Gaussian, core-softened pair interactions. The core-softened pair interactions have two length scales, such that the longer length scale associated with a shallow, attractive well is kept constant while the shorter length scale associated with the repulsive shoulder is varied from an inflexion point to a minimum of progressively increasing depth. As the shoulder well depth increases, the pressure required to form the high density liquid decreases and the temperature up to which the high-density liquid is stable increases, resulting in the shift of the liquid-liquid critical point to much lower pressures and higher temperatures. The pair correlation entropy is computed to show that the excess entropy anomaly diminishes when the shoulder well depth increases. Excess entropy scaling of diffusivity in this class of fluids is demonstrated, showing that decreasing strength of the excess entropy anomaly with increasing shoulder depth results in the progressive loss of water-like thermodynamic, structural and transport anomalies. Instantaneous normal mode analysis was used to index the overall curvature distribution of the fluid and the fraction of imaginary frequency modes was shown to correlate well with the anomalous behaviour of the diffusivity and the pair correlation entropy. The results suggest in the case of core-softened potentials, in addition to the presence of two length scales, energetic and entropic effects associated with local minima and curvatures of the pair interaction play an important role in determining the presence of water-like anomalies and the liquid-liquid phase transition.

preprint2011arXiv

Thermodynamic, Dynamic, Structural and Excess Entropy Anomalies for core-softened potentials

Using molecular dynamic simulations we study three families of continuous core-softened potentials consisting of two length scales: a shoulder scale and an attractive scale. All the families have the same slope between the two length scales but exhibit different potential energy gap between them. For each family three shoulder depths are analyzed. We show that all these systems exhibit a liquid-liquid phase transition between a high density liquid phase and a low density liquid phase ending at a critical point. The critical temperature is the same for all cases suggesting that the critical temperature is only dependent on the slope between the two scales. The critical pressure decreases with the decrease of the potential energy gap between the two scales suggesting that the pressure is responsible for forming the high density liquid. We also show, using the radial distribution function and the excess entropy analysis, that the density, the diffusion and the structural anomalies are present if particles move from the attractive scale to the shoulder scale with the increase of the temperature indicating that the anomalous behavior depends only in what happens up to the second coordination shell.

preprint2010arXiv

Diffusion anomaly and dynamic transitions in the Bell-Lavis water model

In this paper we investigate the dynamic properties of the minimal Bell-Lavis (BL) water model and their relation to the thermodynamic anomalies. The Bell-Lavis model is defined on a triangular lattice in which water molecules are represented by particles with three symmetric bonding arms interacting through van der Waals and hydrogen bonds. We have studied the model diffusivity in different regions of the phase diagram through Monte Carlo simulations. Our results show that the model displays a region of anomalous diffusion which lies inside the region of anomalous density, englobed by the line of temperatures of maximum density (TMD). Further, we have found that the diffusivity undergoes a dynamic transition which may be classified as fragile-to-strong transition at the critical line only at low pressures. At higher densities, no dynamic transition is seen on crossing the critical line. Thus evidence from this study is that relation of dynamic transitions to criticality may be discarded.

preprint2010arXiv

Entropy, diffusivity and the energy landscape of a water-like fluid

Molecular dynamics simulations and instantaneous normal mode (INM) analysis of a fluid with core-softened pair interactions and water-like liquid-state anomalies are performed to obtain an understanding of the relationship between thermodynamics, transport properties and the poten- tial energy landscape. Rosenfeld-scaling of diffusivities with the thermodynamic excess and pair correlation entropy is demonstrated for this model. The INM spectra are shown to carry infor- mation about the dynamical consequences of the interplay between length scales characteristic of anomalous fluids, such as bimodality of the real and imaginary branches of the frequency distribu- tion. The INM spectral information is used to partition the liquid entropy into two contributions associated with the real and imaginary frequency modes; only the entropy contribution from the imaginary branch captures the non-monotonic behaviour of the excess entropy and diffusivity in the anomalous regime of the fluid.

preprint2009arXiv

Dynamic Transitions in a Two Dimensional Associating Lattice Gas Model

Using Monte Carlo simulations we investigate some new aspects of the phase diagram and the behavior of the diffusion coefficient in an associating lattice gas (ALG) model on different regions of the phase diagram. The ALG model combines a two dimensional lattice gas where particles interact through a soft core potential and orientational degrees of freedom. The competition between soft core potential and directional attractive forces results in a high density liquid phase, a low density liquid phase, and a gas phase. Besides anomalies in the behavior of the density with the temperature at constant pressure and of the diffusion coefficient with density at constant temperature are also found. The two liquid phases are separated by a coexistence line that ends in a bicritical point. The low density liquid phase is separated from the gas phase by a coexistence line that ends in tricritical point. The bicritical and tricritical points are linked by a critical $λ$-line. The high density liquid phase and the fluid phases are separated by a second $τ$ critical line. We then investigate how the diffusion coefficient behaves on different regions of the chemical potential-temperature phase diagram. We find that diffusivity undergoes two types of dynamic transitions: a fragile-to-strong trans ition when the critical $λ$-line is crossed by decreasing the temperature at a constant chemical potential; and a strong-to-strong transition when the $τ$-critical line is crossed by decreasing the temperature at a constant chemical potential.

preprint2009arXiv

Liquid crystal phase and waterlike anomalies in a core-softened shoulder-dumbbells system

Using molecular dynamics we investigate the thermodynamics, dynamics and structure of 250 diatomic molecules interacting by a core-softened potential. This system exhibits thermodynamics, dynamics and structural anomalies: a maximum in density-temperature plane at constante pressure and maximum and minimum points in the diffusivity and translational order parameter against density at constant temperature. Starting with very dense systems and decreasing density the mobility at low temperatures first increases, reach a maximum, then decreases, reach a minimum and finally increases. In the pressure-temperature phase diagram the line of maximum translational order parameter is located outside the line of diffusivity extrema that is enclosing the temperature of maximum density line. We compare our results with the monomeric system showing that the anisotropy due to the dumbbell leads to a much larger solid phase and to the appearance of a liquid crystal phase. the double ranged thermodynamic and dynamic anomalies.

preprint2009arXiv

Liquid polymorphism, order-disorder transitions and anomalous behavior: a Monte Carlo study of the Bell-Lavis model for water

The Bell-Lavis model for liquid water is investigated through numerical simulations. The lattice-gas model on a triangular lattice presents orientational states and is known to present a highly bonded low density phase and a loosely bonded high density phase. We show that the model liquid-liquid transition is continuous, in contradiction with mean-field results on the Husimi cactus and from the cluster variational method. We define an order parameter which allows interpretation of the transition as an order-disorder transition of the bond-network. Our results indicate that the order-disorder transition is in the Ising universality class. Previous proposal of an Ehrenfest second order transition is discarded. A detailed investigation of anomalous properties has also been undertaken. The line of density maxima in the HDL phase is stabilized by fluctuations, absent in the mean-field solution.

preprint2009arXiv

Thermodynamic, Dynamic and Structural Anomalies for Shoulder-like potentials

Using molecular dynamic simulations we study a family of continuous core-softened potentials consisting of a hard core, a shoulder at closest distances and an attractive well at further distance. The repulsive shoulder and the well distances represent two length scales. We show that if the first scale, the shoulder, is repulsive or has a small well, the potential has a region in the pressure-temperature phase diagram with density, diffusion and structural anomalies. However, if the closest scale becomes a deep attractive well the regions in the pressure-temperature phase diagram where the three anomalies are present shrink and disappear. This result enables us to predict by the shape of the core-softened potential if anomalies would or would not be present.

preprint2006arXiv

Liquid polymorphism and density anomaly in a three-dimensional associating lattice gas

We investigate the phase diagram of a three-dimensional associating gas $(ALG)$ model. This model combines orientational ice-like interactions and ``van der Waals'' that might be repulsive, representing, in this case, a penalty for distortion of hydrogen bonds. These interactions can be interpreted as two competing distances making the connection between this model and continuous isotropic soft-core potentials. We present Monte Carlo studies of the $ALG$ model showing the presence of two liquid phase, two critical points and A density anomaly.

preprint2004arXiv

Density anomaly in a competing interactions lattice gas model

We study a very simple model of a short-range attraction and an outer shell repulsion as a test system for demixing phase transition and density anomaly. The phase-diagram is obtained by applying mean field analysis and Monte Carlo simulations to a two dimensional lattice gas with nearest-neighbors attraction and next-nearest-neighbors repulsion (the outer shell). Two liquid phases and density anomaly are found. The coexistence line between these two liquid phases meets a critical line between the fluid and the low density liquid at a tricritical point. The line of maximum density emerges in the vicinity of the tricritical point, close to the demixing transition.

preprint1999arXiv

Density-functional theory for attraction between like-charged plates

We study the interactions between two negatively charged macroscopic surfaces confining positive counterions. A density-functional approach is introduced which, besides the usual mean-field interactions, takes into account the correlations in the positions of counterions. The excess free energy is derived in the framework of the Debye-H{ü}ckel theory of the one-component plasma, with the homogeneous density replaced by a weighted density. The minimization of the total free energy yields the density profile of the microions. The pressure is calculated and compared with the simulations and the results derived from integral equations theories. We find that the interaction between the two plates becomes attractive when their separation distance is sufficiently small and the surface charge density is larger than a threshold value.

preprint1999arXiv

The one-component plasma: a conceptual approach

The one-component plasma (OCP) represents the simplest statistical mechanical model of a Coulomb system. For this reason, it has been extensively studied over the last forty years. The advent of the integral equations has resulted in a dramatic improvement in our ability to carry out numerical calculations, but came at the expense of a physical insight gained in a simpler analytic theory. In this paper we present an extension of the Debye-Hueckel (DH) theory to the OCP. The theory allows for analytic calculations of all the thermodynamic functions, as well as the structure factor. The theory explicitly satisfies the Stillinger-Lovett and, for small couplings, the compressibility sum rules, implying its internal self consistency.

preprint1998arXiv

Debye-Hueckel-Bjerrum theory for charged colloids

We formulate an extension of the Debye-Hueckel-Bjerrum theory [M. E. Fisher and Y. Levin, Phys. Rev. Lett. 71, 3826 (1993)] to the fluid state of a highly asymmetric charged colloid. Allowing for the formation of clusters consisting of one polyion and n condensed counterions, the total Helmholtz free energy of the colloidal suspension is constructed. The thermodynamic properties, such as the cluster-density distribution and the pressure, are obtained by the minimization of the free energy under the constraints of fixed number of polyions and counterions. In agreement with the current experimental and Monte Carlo results, no evidence of any phase transition is encountered.

preprint1998arXiv

Polyelectrolyte Solutions with Multivalent Salts

We investigate the thermodynamic properties of a polyelectrolyte solution in a presence of {\it multivalent} salts. The polyions are modeled as rigid cylinders with the charge distributed uniformly along the major axis. The solution, besides the polyions, contain monovalent and divalent counterions as well as monovalent coions. The strong electrostatic attraction existing between the polyions and the counterions results in formation of clusters consisting of one polyion and a number of associated monovalent and divalent counterions. The theory presented in the paper allows us to explicitly construct the Helmholtz free energy of a polyelectrolyte solution. The characteristic cluster size, as well as any other thermodynamic property can then be determined by an appropriate operation on the free energy.

preprint1997arXiv

Phase-diagram for Irregular and Non-Symmetric Cross-linked Polymer Blends

We consider here a blend made of two types of polymers, $A$ and $B$, of different chemical nature. At high temperature the homogeneous mixture is cross-linked. As the temperature is lowered, the two species try to segregate but are kept together by the cross-links. We show that for inhomogeneous, non-regular and non-permanent cross-links, there is a complete segregation at low temperatures if the system is just weakly cross-linked and partial segregation, otherwise. We also demonstrate that there is no phase transition between the homogeneous phase and the microphase for non-symmetric systems. Our analysis is checked with the experiment.

preprint1997arXiv

Sine-Gordon mean field theory of a Coulomb Gas

Sine-Gordon field theory is used to investigate the phase diagram of a neutral Coulomb gas. A variational mean field free energy is constructed and the corresponding phase diagrams in two (2d) and three dimensions (3d) are obtained. When analyzed in terms of chemical potential, the Sine-Gordon theory predicts the phase diagram topologically identical with the Monte Carlo simulations and a recently developed Debye-Hückel-Bjerrum (DHBj) theory. In 2d we find that the infinite order Kosterlitz-Thouless line terminates in a tricritical point, after which the metal-insulator transition becomes first order. However, when the transformation from chemical potential to the density is made the whole of the insulating phase is mapped onto zero density.