A multigrid method is presented for the numerical solution of the linearized Poisson-Boltzmann equation arising in molecular biophysics. The equation is discretized with the finite volume method, and the numerical solution of the discrete equations is accomplished with multiple grid techniques originally developed for two-dimensional interface problems occurring in reactor physics. A detailed analysis of the resulting method is presented for several computer architectures, including comparisons to diagonally scaled CG, ICCG, vectorized ICCG and MICCG, and to SOR provided with an optimal relaxation parameter. Our results indicate that the multigrid method is superior to the preconditioned CG methods and SOR, and that the advantage of multigrid grows with the problem size.
We present a robust and efficient numerical method for solution of the nonlinear Poisson-Boltzmann equation arising in molecular biophysics. The equation is discretized with the box method, and solution of the discrete equations is accomplished with a global inexact-Newton method, combined with linear multilevel techniques we have described in a paper appearing previously in this journal. A detailed analysis of the resulting method is presented, with comparisons to other methods that have been proposed in the literature, including the classical nonlinear multigrid method, the nonlinear conjugate gradient method, and nonlinear relaxation methods such as successive over-relaxation. Both theoretical and numerical evidence suggests that this method will converge in the case of molecules for which many of the existing methods will not. In addition, for problems which the other methods are able to solve, numerical experiments show that the new method is substantially more efficient, and the superiority of this method grows with the problem size. The method is easy to implement once a linear multilevel solver is available, and can also easily be used in conjunction with linear methods other than multigrid.
The nonlinear Poisson-Boltzmann equation (NPBE) provides a continuum description of the electrostatic field in an ionic medium around a macromolecule. Here, a novel approach to the solution of the full NPBE is developed. This robust and efficient algorithm combines multilevel techniques with a damped inexact Newton's method. The CPU time required for solution of the full NPBE, which is less than that for standard single-grid approaches in solving the corresponding linearized equation, is proportional to the number of unknowns enabling applications to very large macromolecular systems. Convergence of the method is demonstrated for a variety of protein systems. Comparison of the solutions to the linearized Poisson-Boltzmann equation shows that the damping of the electrostatic field around the charge is increased and that the potential scales logarithmically with charge. The inclusion of the full nonlinearity thus reduces the impact of highly charged residues on protein surfaces and provides a more realistic representation of electrostatic effects. This is demonstrated through calculation of potential around the active site regions of the 1,266-residue tryptophan synthase dimer and in the computation of rate constants from Brownian dynamics calculations in the superoxide dismutase-superoxide and antibody-antigen systems.
Binary black hole interactions provide potentially the strongest source of gravitational radiation for detectors currently under development. We present some results from the Binary Black Hole Grand Challenge Alliance three-dimensional Cauchy evolution module. These constitute essential steps towards modeling such interactions and predicting gravitational radiation waveforms. We report on single black hole evolutions and the first successful demonstration of a black hole moving freely through a three-dimensional computational grid via a Cauchy evolution: a hole moving ∼ 6M at 0.1c during a total evolution of duration ∼ 60M . The accurate computational modeling of black-hole interactions is essential to the confident detection of astrophysical gravitational radiation by future space-based detectors such as LISA and by the LIGO/VIRGO/GEO complex of ground-based detectors currently under construction. The sensitivity of these detectors will be significantly enhanced if accurate computer simulations of black-hole mergers can produce predictions of radiation waveforms [1]. The Binary Black Hole Grand Challenge Alliance [2] was funded in September 1993 to develop the computational infrastructure necessary accurately to simulate the coalescence of black-hole binaries. The primary objective of the resulting code will be the production of waveforms from binary black hole mergers. In this Letter we report on an important step towards achieving such simulations.A key difficulty in evolving black-hole spacetimes is handling the curvature singularity contained within each hole. The only viable means of accomplishing this over time scales required for binary coalescence appears to be black-hole excision: exclude all or part of the black-hole interior (and the singularity) from the computational domain and evolve only the exterior region [3,4
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.