Skip to content

Newton's Method Converges to Minimum Instead of Barrier in dreid-kernel Buckingham::precompute #20

@TKanX

Description

@TKanX

Description:

dreid-kernel ≤0.4.1's Buckingham::precompute finds the repulsive-wall barrier r_max via Newton's method seeded at r_min. The energy surface is monotonically increasing on (0, r_min), so the iteration converges toward r_min rather than the maximum, yielding r_max ≈ r_min and 2E_max ≈ −2D₀. forge::paramgen::generate_vdw_potentials passes this corrupted tuple verbatim into VdwPairPotential::Buckingham; the downstream kernel then triggers energy reflection at equilibrium distance, rejecting valid geometries as close contacts and producing incorrect VdW energies. Fixed in dreid-kernel v0.4.2.

Tasks:

  • Bump dreid-kernel to 0.4.2

Metadata

Metadata

Assignees

Labels

Type

Projects

Status

Done

Milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions