Differential geometry based solvation model I: Eulerian formulation.

Area of Science:

• Computational chemistry
• Theoretical chemistry
• Applied mathematics

Background:

• Implicit solvation models are crucial for molecular simulations.
• Accurate characterization of solvent-solute boundaries and dielectric functions remains a challenge.
• Existing models often lack stability and differentiability in interface representation.

Purpose of the Study:

• To develop a robust differential geometry-based model for solvation analysis.
• To accurately compute electrostatic potentials and solvent-solute boundaries.
• To improve the accuracy, stability, and efficiency of implicit solvation calculations.

Main Methods:

• Utilizing differential geometry for smooth interface construction.
• Employing geometric measure theory for energy functional formulation.
• Deriving and solving coupled generalized Poisson-Boltzmann (GPBE) and generalized geometric flow equations (GGFE).
• Developing efficient second-order numerical schemes and iterative solvers.

Main Results:

• Obtained accurate electrostatic potentials, solvent-solute boundary profiles, and dielectric functions.
• Demonstrated improved accuracy and stability in implicit solvation calculations.
• Validated the model and numerical methods through extensive experiments on small molecules and proteins.

Conclusions:

• The proposed differential geometry model provides a rigorous and stable framework for solvation analysis.
• The coupled GPBE and GGFE approach enhances the accuracy and robustness of computational solvation.
• The developed numerical methods ensure efficient and reliable computation for complex systems.