A new numerical method for nonlocal electrostatics in biomolecular simulations
โ Scribed by S. Weggler; V. Rutka; A. Hildebrandt
- Publisher
- Elsevier Science
- Year
- 2010
- Tongue
- English
- Weight
- 911 KB
- Volume
- 229
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
โฆ Synopsis
The electrostatic behavior of biomolecules solved in water can be described by an elliptic system of partial differential equations for the potential. In previous studies, this system has been solved by the Boundary Element Method (BEM).
In this paper, we apply the Explicit Jump Immersed Interface Method (EJIIM) as an alternative method to the BEM. Such a finite difference approach allows for a completely automatized software for analyzing biomolecules in their natural surrounding.
The new method shows excellent agreement with the BEM results and has good convergence properties and runtimes. In addition, in contrast to the BEM, where the fundamental solutions of operators are necessary, the EJIIM approach can be easily extended to more complex, especially nonlinear models.
๐ SIMILAR VOLUMES
We present a new simulation code for electrostatic waves in one dimension which uses the Vlasov equation to integrate the distribution function and Ampรจre's equation to integrate the electric field forward in time. Previous Vlasov codes used the Vlasov and Poisson equations. Using Ampรจre's equation
We describe an algorithm for simulating reactive flows in porous media, in which the pore space is mapped explicitly. Chemical reactions at the solid-fluid boundaries lead to dissolution (or precipitation), which makes it necessary to track the movement of the solid-fluid interface during the course
A new unstructured-grid high-order characteristics upwind ยฎnite-volume algorithm for accurate numerical simulation of incompressible laminar ยฏow and forced convection heat transfer over a composite system containing simultaneously porous-saturated region and ยฏuid region is presented in this paper. I
The flux of an evolving wavepacket is the definite time integral of its probability current density. A new method for calculating the flux, based on a Chebychev polynomial expansion of the quantum evolution operator is presented. The central point of the development is that the time integration of t