Skip to content

Conversation

@WHUweiqingzhou
Copy link
Collaborator

@WHUweiqingzhou WHUweiqingzhou commented Nov 1, 2023

To improve the behavior of Kerker preconditioner in ABACUS

Linked Issue

Fix #3094

What's changed?

Current implementation is:
double filter_g = std::max(gg / (gg + gg0), 0.1)
change into
double filter_g = std::max(gg / (gg + gg0), 0.1/this->mixing_beta)

This new implementation can systematically improve the convergence for many case, see detail here.

Owing to this new implementation, we have chance to update the optimal mixing setting for different cases. As a benchmark, I update the optimal mixing setting for daily ABACUS-Test, and get a much better convergence, see parameters here.

I suggest this PR is accepted with update of mixing setting for daily ABACUS-Test.

@WHUweiqingzhou WHUweiqingzhou changed the title improved implementation of kerker preconditioner Feature: improved implementation of kerker preconditioner Nov 1, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Feature: change criterion in kerker dielectric matrix calculations for charge mixing

2 participants