Protein-Ion Binding Process on Finite Macromolecular Concentration. A Poisson-Boltzmann and Monte Carlo Study

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

pK determination. A mean field, Poisson-Boltzmann approach

A model describing dissociation of monoprotonic acid and a method for the determination of its pK value are presented. The model is based on a mean field approximation. The Poisson-Boltzmann equation, adopting spherical symmetry, is numerically solved, and the solution of its linearized form is written. By use of the pH values of a dilution experiment of galacturonic acid as the entry data, the proposed method allowed estimation of the value of pK = 3.25 at a temperature of 25 degrees C. Values for the complex dimensions and dissociation degree are calculated using experimental pH values for solution concentration values ranging from 0.1 to 60 mM. The present analysis leads to the conclusion that the Poisson-Boltzmann equation or its linear form is equally suited for the description of such systems.

Podolsky electromagnetism at finite temperature: Implications on the Stefan-Boltzmann law

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Solution of the Poisson-Boltzmann equation for a system with four ionic species

The analytical solution of the Poisson-Boltzmann equation in an electrolyte with four ionic species (2:2:1:1), in the presence of a charged planar membrane or surface is presented. The function describing the mean electrical potential provides a convenient description that helps the understanding of electrical processes of biological interest.

INTEGRATION OF THE NONLINEAR POISSON BOLTZMANN-EQUATION FOR CHARGED VESICLES IN ELECTROLYTIC SOLUTIONS

The Poisson-Boltzmann equation (PBE), with specific ion-surface interactions and a cell model, was used to calculate the electrostatic properties of aqueous solutions containing vesicles of ionic amphiphiles. Vesicles are assumed to be water- and ion-permeable hollow spheres and specific ion adsorption at the surfaces was calculated using a Volmer isotherm. We solved the PBE numerically for a range of amphiphile and salt concentrations (up to 0.1 M) and calculated co-ion and counterion distributions in the inside and outside of vesicles as well as the fields and electrical potentials. The calculations yield results that are consistent with measured values for vesicles of synthetic amphiphiles.

pH at the micellar interface: Synthesis of pH probes derived from salicylic acid, acid-base dissociation in sodium dodecyl sulfate micelles, and Poisson-Boltzmann simulation

The study of the H+ concentration at the micellar interface is a convenient system for modeling the distribution of H+ at interfaces. We have synthesized salicylic acid derivatives to analyze the proton dissociation of both the carboxylic and phenol groups of' the probes, determining spectrophotometrically the apparent pK(a)'s (pK(ap)) in sodium dodecyl Sulfate, SDS, micelles with and without added salt. The synthesized probes were 2-hydroxy-5-(2-trimethylammoniumacetyl)benzoate; 2-hydroxy-5-(2-dimethylhexadecylammoniumacetyl)benzoate- 2-hydroxy-5-(2-dimethylhexadecylammoniumhexanoyl)benzoate-, 2-hydroxy-5-(2-diniethylhexadecylammoniumundecanoyl)betizoate; 2-hydroxy-5-acetylbenzoic acids and 2-hydroxy-5-dodecanoylbenzoic acid. Upon incorporation into SDS micelles the pK(ap)'s of both carboxylic and phenol groups increased by ca. 3 pH units and NaCl addition caused a decrease in the probe-incorporated pKap. The experimental results were fitted with a cell model Poisson-Boltzmann (P-B) equation taking in consideration the effect of salt on the aggregation number of SDS and using the distance of' the dissociating group as a parameter. The conformations of the probes were analyzed theoretically using two dielectric constants, e.g., 2 and 78. Both the P-B analysis and conformation calculations can be interpreted by assuming that the acid groups dissociate very close to, or at, the interface. Our results are consistent with the assumption that the intrinsic pK(a)'s of both carboxylic and phenol groups of the salicylic acid probes used here can be taken as those in water. Using this assumption the micellar and salt effects on the pKap's of the (trialkylammonium)benzoate probes were described accurately using a cell model P-B analysis. (c) 2005 Elsevier B.V. All rights reserved.

Aplicação da equação de Paisson-Boltzmann com condição de contorno para regulação carga-potencial a sistemas polieletrolíticos

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

Model selection for discriminative restricted boltzmann machines through meta-heuristic techniques

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Coupled Boltzmann calculation of mixed axion/neutralino cold dark matter production in the early universe

We calculate the relic abundance of mixed axion/neutralino cold dark matter which arises in R-parity conserving supersymmetric (SUSY) models wherein the strong CP problem is solved by the Peccei-Quinn (PQ) mechanism with a concommitant axion/saxion/axino supermultiplet. By numerically solving the coupled Boltzmann equations, we include the combined effects of 1. thermal axino production with cascade decays to a neutralino LSP, 2. thermal saxion production and production via coherent oscillations along with cascade decays and entropy injection, 3. thermal neutralino production and re-annihilation after both axino and saxion decays, 4. gravitino production and decay and 5. axion production both thermally and via oscillations. For SUSY models with too high a standard neutralino thermal abundance, we find the combined effect of SUSY PQ particles is not enough to lower the neutralino abundance down to its measured value, while at the same time respecting bounds on late-decaying neutral particles from BBN. However, models with a standard neutralino underabundance can now be allowed with either neutralino or axion domination of dark matter, and furthermore, these models can allow the PQ breaking scale f(a) to be pushed up into the 10(14) - 10(15) GeV range, which is where it is typically expected to be in string theory models.

Spectrum reconstruction from a scattering measurement using the adjoint Boltzmann transport equation for photons

Quality control of medical radiological systems is of fundamental importance, and requires efficient methods for accurately determine the X-ray source spectrum. Straightforward measurements of X-ray spectra in standard operating require the limitation of the high photon flux, and therefore the measure has to be performed in a laboratory. However, the optimal quality control requires frequent in situ measurements which can be only performed using a portable system. To reduce the photon flux by 3 magnitude orders an indirect technique based on the scattering of the X-ray source beam by a solid target is used. The measured spectrum presents a lack of information because of transport and detection effects. The solution is then unfolded by solving the matrix equation that represents formally the scattering problem. However, the algebraic system is ill-conditioned and, therefore, it is not possible to obtain a satisfactory solution. Special strategies are necessary to circumvent the ill-conditioning. Numerous attempts have been done to solve this problem by using purely mathematical methods. In this thesis, a more physical point of view is adopted. The proposed method uses both the forward and the adjoint solutions of the Boltzmann transport equation to generate a better conditioned linear algebraic system. The procedure has been tested first on numerical experiments, giving excellent results. Then, the method has been verified with experimental measurements performed at the Operational Unit of Health Physics of the University of Bologna. The reconstructed spectra have been compared with the ones obtained with straightforward measurements, showing very good agreement.

Thermal fluctuations and boundary conditions in the lattice Boltzmann method

The lattice Boltzmann method is a popular approach for simulating hydrodynamic interactions in soft matter and complex fluids. The solvent is represented on a discrete lattice whose nodes are populated by particle distributions that propagate on the discrete links between the nodes and undergo local collisions. On large length and time scales, the microdynamics leads to a hydrodynamic flow field that satisfies the Navier-Stokes equation. In this thesis, several extensions to the lattice Boltzmann method are developed. In complex fluids, for example suspensions, Brownian motion of the solutes is of paramount importance. However, it can not be simulated with the original lattice Boltzmann method because the dynamics is completely deterministic. It is possible, though, to introduce thermal fluctuations in order to reproduce the equations of fluctuating hydrodynamics. In this work, a generalized lattice gas model is used to systematically derive the fluctuating lattice Boltzmann equation from statistical mechanics principles. The stochastic part of the dynamics is interpreted as a Monte Carlo process, which is then required to satisfy the condition of detailed balance. This leads to an expression for the thermal fluctuations which implies that it is essential to thermalize all degrees of freedom of the system, including the kinetic modes. The new formalism guarantees that the fluctuating lattice Boltzmann equation is simultaneously consistent with both fluctuating hydrodynamics and statistical mechanics. This establishes a foundation for future extensions, such as the treatment of multi-phase and thermal flows. An important range of applications for the lattice Boltzmann method is formed by microfluidics. Fostered by the "lab-on-a-chip" paradigm, there is an increasing need for computer simulations which are able to complement the achievements of theory and experiment. Microfluidic systems are characterized by a large surface-to-volume ratio and, therefore, boundary conditions are of special relevance. On the microscale, the standard no-slip boundary condition used in hydrodynamics has to be replaced by a slip boundary condition. In this work, a boundary condition for lattice Boltzmann is constructed that allows the slip length to be tuned by a single model parameter. Furthermore, a conceptually new approach for constructing boundary conditions is explored, where the reduced symmetry at the boundary is explicitly incorporated into the lattice model. The lattice Boltzmann method is systematically extended to the reduced symmetry model. In the case of a Poiseuille flow in a plane channel, it is shown that a special choice of the collision operator is required to reproduce the correct flow profile. This systematic approach sheds light on the consequences of the reduced symmetry at the boundary and leads to a deeper understanding of boundary conditions in the lattice Boltzmann method. This can help to develop improved boundary conditions that lead to more accurate simulation results.

Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

Cosmological predictions for a scalar tensor dark energy model by a dedicated Einstein-Boltzmann code

We have extended the Boltzmann code CLASS and studied a specific scalar tensor dark energy model: Induced Gravity