Determination of interfacial co-ion concentration in ionic micelles by chemical trapping: Halide concentration at the interface of sodium dodecyl sulfate micelles

Products from the spontaneous reaction of a long-chain arenediazonium salt, 2,6-dimethyl-4-hexadecylbenzenediazonium tetrafluoroborate(16-ArN2BF4), in aqueous micellar solutions of sodium dodecyl sulfate (SDS)? are used to estimate the local concentration of chloride and bromide ions at the micellar surface. The arenediazonium ion, 16-ArN2+, which is totally bound to the SDS micelle, reacts by rate-determining loss of N-2 to give an aryl cation that traps available nucleophiles, i,e., H2O, Cl-, and Br-, to give stable phenol, 16-ArOH, and halobenzene products, 16-ArCl and 16-ArBr, respectively. Product yields, determined by HPLC, are related to local concentrations using calibration curves obtained from independent standards. The local concentrations determined by this method are consistent with co-ion concentrations calculated, using a cell model, by numerical integration of the Poisson-Boltzmann equation (PBE) taking into account salt-induced micellar growth. The salt dependence of the intel facial concentrations of Cl- and Br- are identical. indicating no specific interactions in the interfacial co-ion compartment. PBE calculations predict that, in micellar SDS, increasing the concentration of a particular halide salt (NaX) at constant concentration of another halide (NaY) should result in an increase in the local concentrations of both co-ions. Using this chemical-trapping method, this prediction was demonstrated experimentally.

Determinação do grau de ionização de aminoácidos polares carregados

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Determinação do gradiente de pH na região interna a duas membranas em equilíbrio com eletrólito

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Amphiphilic planar membranes in ionic equilibrium: a study of pH position-dependent values

The pH values near a planar dissociating membrane are studied under a mean field approximation using the Poisson-Boltzmann equation and its linear form. The equations are solved in planar symmetry with the consideration that the charge density on the dissociating membrane surface results from an equilibrium process with the neighboring electrolyte. Results for the membrane dissociation degree are presented as a function of the electrolyte ionic strength and membrane surface charge density. Our calculations indicate that pH values have an appreciable variation within 2 nm from the membrane. It is shown that the dissociation process is enhanced due to the presence of bivalent ions and that pH values acquire better stability than in an electrolyte containing univalent ions.

Proteins in silico-modeling and sampling

Proteins are linear chain molecules made out of amino acids. Only when they fold to their native states, they become functional. This dissertation aims to model the solvent (environment) effect and to develop & implement enhanced sampling methods that enable a reliable study of the protein folding problem in silico. We have developed an enhanced solvation model based on the solution to the Poisson-Boltzmann equation in order to describe the solvent effect. Following the quantum mechanical Polarizable Continuum Model (PCM), we decomposed net solvation free energy into three physical terms– Polarization, Dispersion and Cavitation. All the terms were implemented, analyzed and parametrized individually to obtain a high level of accuracy. In order to describe the thermodynamics of proteins, their conformational space needs to be sampled thoroughly. Simulations of proteins are hampered by slow relaxation due to their rugged free-energy landscape, with the barriers between minima being higher than the thermal energy at physiological temperatures. In order to overcome this problem a number of approaches have been proposed of which replica exchange method (REM) is the most popular. In this dissertation we describe a new variant of canonical replica exchange method in the context of molecular dynamic simulation. The advantage of this new method is the easily tunable high acceptance rate for the replica exchange. We call our method Microcanonical Replica Exchange Molecular Dynamic (MREMD). We have described the theoretical frame work, comment on its actual implementation, and its application to Trp-cage mini-protein in implicit solvent. We have been able to correctly predict the folding thermodynamics of this protein using our approach.

Dressed polyions, counterion condensation, and adsorption excess in polyelectrolyte solutions.

The phenomenon of Manning-Oosawa counterion condensation is given an explicit statistical mechanical and qualitative basis via a dressed polyelectrolyte formalism in connection with the topology of the electrostatic free-energy surface and is derived explicitly in terms of the adsorption excess of ions about the polyion via the nonlinear Poisson-Boltzmann equation. The approach is closely analogous to the theory of ion binding in micelles. Our results not only elucidate a Poisson-Boltzmann analysis, which shows that a fraction of the counterions lie within a finite volume around the polyion even if the volume of the system tends towards infinity, but also provide a direct link between Manning's theta-the number of condensed counterions for each polyion site-and a statistical thermodynamic quantity, namely, the adsorption excess per monomer.

On the magnitude of the electrostatic contribution to ligand-DNA interactions.

A model based on the nonlinear Poisson-Boltzmann equation is used to study the electrostatic contribution to the binding free energy of a simple intercalating ligand, 3,8-diamino-6-phenylphenanthridine, to DNA. We find that the nonlinear Poisson-Boltzmann model accurately describes both the absolute magnitude of the pKa shift of 3,8-diamino-6-phenylphenanthridine observed upon intercalation and its variation with bulk salt concentration. Since the pKa shift is directly related to the total electrostatic binding free energy of the charged and neutral forms of the ligand, the accuracy of the calculations implies that the electrostatic contributions to binding are accurately predicted as well. Based on our results, we have developed a general physical description of the electrostatic contribution to ligand-DNA binding in which the electrostatic binding free energy is described as a balance between the coulombic attraction of a ligand to DNA and the disruption of solvent upon binding. Long-range coulombic forces associated with highly charged nucleic acids provide a strong driving force for the interaction of cationic ligands with DNA. These favorable electrostatic interactions are, however, largely compensated for by unfavorable changes in the solvation of both the ligand and the DNA upon binding. The formation of a ligand-DNA complex removes both charged and polar groups at the binding interface from pure solvent while it displaces salt from around the nucleic acid. As a result, the total electrostatic binding free energy is quite small. Consequently, nonpolar interactions, such as tight packing and hydrophobic forces, must play a significant role in ligand-DNA stability.

The charge effect on the hindrance factors for diffusion and convection of a solute in pores II

The diffusion and convection of a solute suspended in a fluid across porous membranes are known to be reduced compared to those in a bulk solution, owing to the fluid mechanical interaction between the solute and the pore wall as well as steric restriction. If the solute and the pore wall are electrically charged, the electrostatic interaction between them could affect the hindrance to diffusion and convection. In this study, the transport of charged spherical solutes through charged circular cylindrical pores filled with an electrolyte solution containing small ions was studied numerically by using a fluid mechanical and electrostatic model. Based on a mean field theory, the electrostatic interaction energy between the solute and the pore wall was estimated from the Poisson-Boltzmann equation, and the charge effect on the solute transport was examined for the solute and pore wall of like charge. The results were compared with those obtained from the linearized form of the Poisson-Boltzmann equation, i.e.the Debye-Hückel equation. © 2012 The Japan Society of Fluid Mechanics and IOP Publishing Ltd.

Fluid mechanical and electrostatic study on the osmotic flow through cylindrical pores

When two solutions differing in solute concentration are separated by a porous membrane, the osmotic pressure will generate a net volume flux of the suspending fluid across the membrane; this is termed osmotic flow. We consider the osmotic flow across a membrane with circular cylindrical pores when the solute and the pore walls are electrically charged, and the suspending fluid is an electrolytic solution containing small cations and anions. Under the condition in which the radius of the pores and that of the solute molecules greatly exceed those of the solvent as well as the ions, a fluid mechanical and electrostatic theory is introduced to describe the osmotic flow in the presence of electric charge. The interaction energy, including the electrostatic interaction between the solute and the pore wall, plays a key role in determining the osmotic flow. We examine the electrostatic effect on the osmotic flow and discuss the difference in the interaction energy determined from the nonlinear Poisson-Boltzmann equation and from its linearized equation (the Debye-Hückel equation).

Electrical charge effect on osmotic flow through pores

An electrostatic model for osmotic flow through circular cylindrical pores is developed to describe the reflection coefficient for the membrane transport in the presence of surface charges on the pore wall and the solute. For a spherical solute placed at an arbitrary radial position in the pore, the electrical potential was computed by a spectral element method applied to the Poisson-Boltzmann equation together with the condition of electrical neutrality. The interaction energy between the surface charges was used to estimate the osmotic reflection coefficient. The proposed model predicts that even for a small Debye length compared to the pore radius, the repulsive electrostatic interaction between the surface charges could significantly increase the osmotic flow through the pore.

Simulations of charged droplet collisions in shear flow

Acknowledgments This research has been enabled by the use of computing resources provided by WestGrid, the Shared Hierarchical Academic Research Computing Network (SHARCNET: www.sharcnet.ca), and Compute/Calcul Canada. O.S. thanks NSERC for an Alexander Graham Bell Canada Graduate Scholarship.

Determination of local density of states in amorphous oxide semiconductors with Kelvin Probe Force Microscopy

Amorphous semiconductors are important materials as they can be deposited by physical deposition techniques on large areas and even on plastic substrates. Therefore, they are crucial for transistors in large active matrices for imaging and transparent wearable electronics. The most widely applied candidate for amorphous thin film transistors production is Indium Gallium Zinc Oxide (IGZO). It is attracting much interest because of its optical transparency, facile processing by sputtering deposition and notable improved charge carrier mobility with respect to hydrogenated amorphous silicon a-Si:H. Degradation of the device and long-term performance issues have been observed if IGZO thin film transistors are subjected to electrical stress, leading to a modification of IGZO channel properties and subthreshold slope. Therefore, it is of great interest to have a reliable and precise method to study the conduction band tail, and the density of states in amorphous semiconductors. The aim of this thesis is to develop a local technique using Kelvin Probe Force Microscopy to study the evolution of IGZO DOS properties. The work is divided into three main parts. First, solutions to the non-linear Poisson-Boltzmann equation of a metal-insulator-semiconductor junction describing the charge accumulation and its relation to DOS properties are elaborated. Second macroscopic techniques such as capacitance voltage (CV) measurements and photocurrent spectroscopy are applied to obtain a non-local estimate of band-tail DOS properties in thin film transistor samples. The third part of my my thesis is dedicated to the KPFM measurements. By fitting the data to the developed numerical model, important parameters describing the amorphous conduction band tail are obtained. The results are in excellent agreement with the macroscopic characterizations. KPFM result is comparable also with non-local optoelectronic characterizations, such as photocurrent spectroscopy.

A Direct Numerov Sixth-order Numerical Scheme to Accurately Solve the Unidimensional Poisson Equation with Dirichlet Boundary Conditions

In this article, we present an analytical direct method, based on a Numerov three-point scheme, which is sixth order accurate and has a linear execution time on the grid dimension, to solve the discrete one-dimensional Poisson equation with Dirichlet boundary conditions. Our results should improve numerical codes used mainly in self-consistent calculations in solid state physics.

Mixed initial-boundary value problem in particle modeling of microelectronic devices

The Boltzmann equation in presence of boundary and initial conditions, which describes the general case of carrier transport in microelectronic devices is analysed in terms of Monte Carlo theory. The classical Ensemble Monte Carlo algorithm which has been devised by merely phenomenological considerations of the initial and boundary carrier contributions is now derived in a formal way. The approach allows to suggest a set of event-biasing algorithms for statistical enhancement as an alternative of the population control technique, which is virtually the only algorithm currently used in particle simulators. The scheme of the self-consistent coupling of Boltzmann and Poisson equation is considered for the case of weighted particles. It is shown that particles survive the successive iteration steps.

Modelos cinéticos da equação linearizada de Boltzmann e um problema de transferência de calor em microescala

Neste trabalho, um problema de transferência de calor da dinâmica de gases rarefeitos, causado pela diferença de temperaturas nas superfícies de um canal, é abordado. O problema é formulado através dos modelos cinéticos BGK, S e Gross-Jackson da equação linearizada de Boltzmann e resolvido, de forma unificada, pelo método analítico de ordenadas discretas (método ADO). Resultados numéricos para as perturbações de densidade e temperatura e também para o fluxo de calor são apresentados e comparados, mostrando que não se pode dizer que algum dos três modelos seja uma melhor aproximação da solução aos resultados da equação linearizada de Boltzmann.