Automated DFTB Repulsive Potential Fitting via Differential Evolution
Find Innovative SolutionsGenerate Solutions
Solution Overview
Problem
Existing DFTB methods for fitting repulsive potentials require extensive manual intervention and are inefficient due to the need for multiple bonding forms and parameter compatibility, leading to slow development and inaccuracy in simulating different bonding forms like C—O and C═O interactions.
Innovation Solution
A method using a differential evolution algorithm to automatically fit repulsive potentials by generating saturated molecules, performing high-accuracy energy calculations, and applying spline interpolation with optimized division points and sample points to achieve accurate and efficient fitting of multiple repulsive potential curves.
Engineering Contradictions & Design Principles
Engineering Contradiction Analysis
1Manufacturing precision
If manual spline fitting with large quantity of trained sample points is used, then fitting accuracy can be improved, but the development speed and efficiency deteriorate due to extensive manual intervention
Solution Approach 1:
The system performs self-service by automatically generating training data through computational chemistry calculations and executing the fitting process without manual intervention. The automated workflow includes generating molecular structures, performing energy calculations, and fitting repulsive potentials sequentially, eliminating the need for manual train set construction and parameter adjustment
Solution Approach 2:
The system performs preliminary actions by pre-generating saturated molecules with hydrogens and pre-calculating training data using high-accuracy methods before the actual fitting process. This prepares all necessary data and parameters in advance, enabling the fitting algorithm to proceed efficiently without manual data preparation
2Adaptability or versatility
If large quantity of trained sample points are used to cover different bonding forms, then fitting coverage is improved, but the complexity of train set construction and parameter compatibility deteriorate
Solution Approach 1:
The system achieves universality by using a unified automated fitting framework that handles multiple bonding forms (single bonds, double bonds, triple bonds) and different atomic pairs through the same process. The differential evolution algorithm and spline fitting methodology work consistently across all atom types and bonding scenarios, eliminating the need for separate manual fitting procedures for each case
Solution Approach 2:
The system applies segmentation by dividing the fitting process into distinct automated stages: generating saturated molecules for each atomic pair, performing spacing scan calculations, calculating energy differences, and executing the differential evolution fitting algorithm. This structured segmentation handles complexity through systematic breakdown rather than manual coordination
3Reliability
If repeated iteration of parameter development process is performed to ensure compatibility, then parameter reliability is improved, but the development speed deteriorates
Solution Approach 1:
The system implements feedback through the differential evolution algorithm, which iteratively optimizes repulsive potential parameters by evaluating fitness based on how well the parameters reproduce high-accuracy reference data. The algorithm automatically adjusts parameters across generations, providing continuous feedback on parameter quality and convergence, ensuring compatibility without manual iteration
Solution Approach 2:
The system replaces the mechanical manual iteration process with an automated computational optimization system. Instead of manually adjusting parameters and re-running calculations to check compatibility, the differential evolution algorithm automatically performs the iterative optimization using computational mechanics, significantly reducing the time and human effort required while maintaining or improving parameter reliability
Data Source
AI summary
A method for automatically and efficiently fitting repulsive potentials through DFTB includes: optimizing a molecule containing an atomic pair, and performing spacing scan on the atomic pair according to an acting interval of repulsive potentials; performing high-accuracy energy calculation on a scanning result; performing difference calculation on obtained energy to obtain second-order derivatives of force and energy; saving a structure, the energy, the second-order derivatives of force and energy into a database; acquiring energy without a repulsive potential as well as the second-order derivatives of force and energy to obtain target values of the repulsive potentials; transforming the separated target values of the repulsive potentials into splines; splicing multiple splines on a mean position between equilibrium position points to obtain new train data; constructing a train set by means of necessary samples and choosable samples; and fitting the repulsive potentials through singular value decomposition.