IDEAS home Printed from https://ideas.repec.org/a/eee/apmaco/v380y2020ics0096300320302368.html
   My bibliography  Save this article

Pseudo-transient ghost fluid methods for the Poisson-Boltzmann equation with a two-component regularization

Author

Listed:
  • Ahmed Ullah, Sheik
  • Zhao, Shan

Abstract

The Poisson Boltzmann equation (PBE) is a well-established implicit solvent continuum model for the electrostatic analysis of solvated biomolecules. The numerical solution of the nonlinear PBE is still a challenge due to its exponential nonlinear term, strong singularity by the source terms, and distinct dielectric regions. In this paper, a new pseudo-transient approach is proposed, which combines an analytical treatment of singular charges in a two-component regularization, with an analytical integration of nonlinear term in pseudo-time solution. To ensure efficiency, both fully implicit alternating direction implicit (ADI) and unconditionally stable locally one-dimensional (LOD) methods have been constructed to decompose three-dimensional linear systems into one-dimensional (1D) ones in each pseudo-time step. Moreover, to accommodate the nonzero function and flux jumps across the dielectric interface, a modified ghost fluid method (GFM) has been introduced as a first order accurate sharp interface method in 1D style, which minimizes the information needed for the molecular surface. The 1D finite-difference matrix generated by the GFM is symmetric and diagonally dominant, so that the stability of ADI and LOD methods is boosted. The proposed pseudo-transient GFM schemes have been numerically validated by calculating solvation free energy, binding energy, and salt effect of various proteins. It has been found that with the augmentation of regularization and GFM interface treatment, the ADI method not only enhances the accuracy dramatically, but also improves the stability significantly. By using a large time increment, an efficient protein simulation can be realized in steady-state solutions. Therefore, the proposed GFM-ADI and GFM-LOD methods provide accurate, stable, and efficient tools for biomolecular simulations.

Suggested Citation

  • Ahmed Ullah, Sheik & Zhao, Shan, 2020. "Pseudo-transient ghost fluid methods for the Poisson-Boltzmann equation with a two-component regularization," Applied Mathematics and Computation, Elsevier, vol. 380(C).
  • Handle: RePEc:eee:apmaco:v:380:y:2020:i:c:s0096300320302368
    DOI: 10.1016/j.amc.2020.125267
    as

    Download full text from publisher

    File URL: http://www.sciencedirect.com/science/article/pii/S0096300320302368
    Download Restriction: Full text for ScienceDirect subscribers only

    File URL: https://libkey.io/10.1016/j.amc.2020.125267?utm_source=ideas
    LibKey link: if access is restricted and if your library uses this service, LibKey will redirect you to where you can use your library subscription to access this item
    ---><---

    As the access to this document is restricted, you may want to search for a different version of it.

    References listed on IDEAS

    as
    1. Li, Chuan & Zhao, Shan, 2017. "A matched Peaceman–Rachford ADI method for solving parabolic interface problems," Applied Mathematics and Computation, Elsevier, vol. 299(C), pages 28-44.
    2. Ewing, Richard E. & Li, Zhilin & Lin, Tao & Lin, Yanping, 1999. "The immersed finite volume element methods for the elliptic interface problems," Mathematics and Computers in Simulation (MATCOM), Elsevier, vol. 50(1), pages 63-76.
    Full references (including those not matched with items on IDEAS)

    Citations

    Citations are extracted by the CitEc Project, subscribe to its RSS feed for this item.
    as


    Cited by:

    1. Shao, Yuanzhen & McGowan, Mark & Wang, Siwen & Alexov, Emil & Zhao, Shan, 2023. "Convergence of a diffuse interface Poisson-Boltzmann (PB) model to the sharp interface PB model: A unified regularization formulation," Applied Mathematics and Computation, Elsevier, vol. 436(C).

    Most related items

    These are the items that most often cite the same works as this one and are cited by the same works as this one.
    1. Coco, A. & Semplice, M. & Serra Capizzano, S., 2020. "A level-set multigrid technique for nonlinear diffusion in the numerical simulation of marble degradation under chemical pollutants," Applied Mathematics and Computation, Elsevier, vol. 386(C).
    2. Mengya Su & Liuqing Xie & Zhiyue Zhang, 2022. "Numerical Analysis of Fourier Finite Volume Element Method for Dirichlet Boundary Optimal Control Problems Governed by Elliptic PDEs on Complex Connected Domains," Mathematics, MDPI, vol. 10(24), pages 1-26, December.
    3. Feng, Qiwei & Han, Bin & Minev, Peter, 2022. "A high order compact finite difference scheme for elliptic interface problems with discontinuous and high-contrast coefficients," Applied Mathematics and Computation, Elsevier, vol. 431(C).
    4. Peng, Jie & Shu, Shi & Yu, HaiYuan & Feng, Chunsheng & Kan, Mingxian & Wang, Ganghua, 2017. "Error estimates on a finite volume method for diffusion problems with interface on rectangular grids," Applied Mathematics and Computation, Elsevier, vol. 311(C), pages 335-352.
    5. Kishore Kumar, N. & Biswas, Pankaj, 2021. "Fully discrete least-squares spectral element method for parabolic interface problems," Mathematics and Computers in Simulation (MATCOM), Elsevier, vol. 181(C), pages 364-379.

    Corrections

    All material on this site has been provided by the respective publishers and authors. You can help correct errors and omissions. When requesting a correction, please mention this item's handle: RePEc:eee:apmaco:v:380:y:2020:i:c:s0096300320302368. See general information about how to correct material in RePEc.

    If you have authored this item and are not yet registered with RePEc, we encourage you to do it here. This allows to link your profile to this item. It also allows you to accept potential citations to this item that we are uncertain about.

    If CitEc recognized a bibliographic reference but did not link an item in RePEc to it, you can help with this form .

    If you know of missing items citing this one, you can help us creating those links by adding the relevant references in the same way as above, for each refering item. If you are a registered author of this item, you may also want to check the "citations" tab in your RePEc Author Service profile, as there may be some citations waiting for confirmation.

    For technical questions regarding this item, or to correct its authors, title, abstract, bibliographic or download information, contact: Catherine Liu (email available below). General contact details of provider: https://www.journals.elsevier.com/applied-mathematics-and-computation .

    Please note that corrections may take a couple of weeks to filter through the various RePEc services.

    IDEAS is a RePEc service. RePEc uses bibliographic data supplied by the respective publishers.