Unconditionally stable time splitting methods for the electrostatic analysis of solvated biomolecules

Simple item page
dc.contributorHadji, Layachi
dc.contributorLeClair, Patrick R.
dc.contributorZhu, Wei
dc.contributor.advisorZhao, Shan
dc.contributor.authorWilson, Leighton Wayne
dc.contributor.otherUniversity of Alabama Tuscaloosa
dc.date.accessioned2017-03-01T17:22:52Z
dc.date.available2017-03-01T17:22:52Z
dc.date.issued2015
dc.descriptionElectronic Thesis or Dissertationen_US
dc.description.abstractThis work introduces novel unconditionally stable operator splitting methods for solving the time dependent nonlinear Poisson-Boltzmann (NPB) equation for the electrostatic analysis of solvated biomolecules. In a pseudo-transient continuation solution of the NPB equation, a long time integration is needed to reach the steady state. This calls for time stepping schemes that are stable and accurate for large time increments. The existing alternating direction implicit (ADI) methods for the NPB equation are known to be conditionally stable, although fully implicit. To overcome this difficulty, we propose several new operator splitting schemes, including the multiplicative locally one-dimensional (LOD) schemes and additive operator splitting (AOS) schemes. The nonlinear term is integrated analytically in these schemes, while standard discretizations with finite differences in space and implicit time integrations are used. The proposed schemes become much more stable than the ADI methods, and some of them are indeed unconditionally stable in dealing with solvated proteins with source singularities and non-smooth solutions. Numerically, the orders of convergence in both space and time are found to be one. Nevertheless, the precision in calculating the electrostatic free energy is low, unless a small time increment is used. Further accuracy improvements are thus considered, through constructing a Richardson extrapolation procedure and a tailored recovery scheme that replaces the fast Fourier transform method by the operator splitting method in the vacuum case. After acceleration, the optimized LOD method can produce a reliable energy estimate by integrating for a small and fixed number of time steps. Since one only needs to solve a tridiagonal linear system in each independent one dimensional process, the overall computation is very efficient. The unconditionally stable LOD method scales linearly with respect to the number of atoms in the protein studies, and is over 20 times faster than the conditionally stable ADI methods. In addition, some preliminary results on increased stability for ADI methods using a regularization scheme are presented.en_US
dc.format.extent63 p.
dc.format.mediumelectronic
dc.format.mimetypeapplication/pdf
dc.identifier.otheru0015_0000001_0001874
dc.identifier.otherWilson_alatus_0004M_12291
dc.identifier.urihttps://ir.ua.edu/handle/123456789/2308
dc.languageEnglish
dc.language.isoen_US
dc.publisherUniversity of Alabama Libraries
dc.relation.hasversionborn digital
dc.relation.ispartofThe University of Alabama Electronic Theses and Dissertations
dc.relation.ispartofThe University of Alabama Libraries Digital Collections
dc.rightsAll rights reserved by the author unless otherwise indicated.en_US
dc.subjectApplied mathematics
dc.subjectBiophysics
dc.subjectBiochemistry
dc.titleUnconditionally stable time splitting methods for the electrostatic analysis of solvated biomoleculesen_US
dc.typethesis
dc.typetext
etdms.degree.departmentUniversity of Alabama. Department of Mathematics
etdms.degree.disciplineMathematics
etdms.degree.grantorThe University of Alabama
etdms.degree.levelmaster's
etdms.degree.nameM.A.
Files
Original bundle
Now showing 1 - 1 of 1
Thumbnail Image
Name:
file_1.pdf
Size:
1.4 MB
Format:
Adobe Portable Document Format
Download
Collections
Theses and Dissertations
Theses and Dissertations - Department of Mathematics