Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSolution Poisson-Boltzmann equation: Application in the Human Neuron Membrane
107
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
With already demonstrated in previous work the equations that describe the space dependence of the electric potential are determined by the solution of the equation of Poisson-Boltzmann. In this work we consider these solutions for the membrane of the human neuron, using a model simplified for this structure considering the distribution of electrolytes in each side of the membrane, as well as the effect of glycocalyx and the lipidic bilayer. It was assumed that on both sides of the membrane the charges are homogeneously distributed and that the potential depends only on coordinate z.
rate research
Read More
We consider the problem of computing approximate solution of Poisson equation in the low-parametric tensor formats. We propose a new algorithm to compute the solution based on the cross approximation algorithm in the frequency space, and it has better complexity with respect to ranks in comparison with standard algorithms, which are based on the exponential sums approximation. To illustrate the effectiveness of our solver, we incorporate into a Uzawa solver for the Stokes problem on semi-staggered grid as a subsolver. The resulting solver outperforms the standard method for $n geq 256$.
In numerical simulations of many charged systems at the micro/nano scale, a common theme is the repeated solution of the Poisson-Boltzmann equation. This task proves challenging, if not entirely infeasible, largely due to the nonlinearity of the equation and the high dimensionality of the physical and parametric domains with the latter emulating the system configuration. In this paper, we for the first time adapt a mathematically rigorous and computationally efficient model order reduction paradigm, the so-called reduced basis method (RBM), to mitigate this challenge. We adopt a finite difference method as the mandatory underlying scheme to produce the {em truth approximations} of the RBM upon which the fast algorithm is built and its performance is measured against. Numerical tests presented in this paper demonstrate the high efficiency and accuracy of the fast algorithm, the reliability of its error estimation, as well as its capability in effectively capturing the boundary layer.
The Chapman-Enskog method of solution of the relativistic Boltzmann equation is generalized in order to admit a time-derivative term associated to a thermodynamic force in its first order solution. Both existence and uniqueness of such a solution are proved based on the standard theory of integral equations. The mathematical implications of the generalization here introduced are thoroughly discussed regarding the nature of heat as chaotic energy transfer in the context of relativity theory.
The Poisson-Boltzmann equation is a widely used model to study the electrostatics in molecular solvation. Its numerical solution using a boundary integral formulation requires a mesh on the molecular surface only, yielding accurate representations of the solute, which is usually a complicated geometry. Here, we utilize adjoint-based analyses to form two goal-oriented error estimates that allows us to determine the contribution of each discretization element (panel) to the numerical error in the solvation free energy. This information is useful to identify high-error panels to then refine them adaptively to find optimal surface meshes. We present results for spheres and real molecular geometries, and see that elements with large error tend to be in regions where there is a high electrostatic potential. We also find that even though both estimates predict different total errors, they have similar performance as part of an adaptive mesh refinement scheme. Our test cases suggest that the adaptive mesh refinement scheme is very effective, as we are able to reduce the error one order of magnitude by increasing the mesh size less than 20%. This result sets the basis towards efficient automatic mesh refinement schemes that produce optimal meshes for solvation energy calculations.
A non-linear Boltzmann equation describing the time evolution of a partonic system in the central rapidity region after a heavy ion collision is solved numerically. A particular model of the collinear logarithmic divergences due to small angle scattering is employed in the numerical solution. The system is followed until it reaches kinetic equilibrium where the equilibration time, temperature and chemical potential are determined for both RHIC and LHC.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
