in silico
Phase 1 : In Silico Feasibility & Transferability Screening
- Objectif
- To falsify the core hypothesis at minimal cost by testing the two most critical assumptions: (1) the fragmentation error is ≤1.5 kcal/mol for representative transition states and minima, and (2) a preliminary MLFF trained on small fragments can predict DFT forces of a full system with MAE ≤ 1.5 kcal/mol/Å.
- Coût estimé
- €1,500 - €3,000 (mainly for HPC compute time)
- Durée estimée
- 5-7 weeks
- Critères de réussite
- Fragmentation Error (Prediction 3) · ≤ 1.5 kcal/mol in at least 8 out of 9 test systems (both minima and TS) · (Single-point DFT energy difference: E(full_system) - Σ E(fragments) + E(overlap_correction))
- MLFF Force MAE on Full System · ≤ 1.5 kcal/mol/Å (initial target, to be refined) · (RMSE between MLFF-predicted forces and DFT forces on the full-system AIMD trajectory)
- On continue si
- Both success criteria are met: fragmentation error is ≤1.5 kcal/mol in ≥8/9 systems AND the MLFF force MAE on the full system is ≤1.5 kcal/mol/Å.
- On arrête si
- Fragmentation error is >1.5 kcal/mol in ≥2 systems OR the MLFF force MAE is >2.5 kcal/mol/Å, indicating a fundamental breakdown of the locality assumption.
- On réoriente si
- If the fragmentation error is acceptable but the MLFF force MAE is high (1.5-2.5 kcal/mol/Å), we will pivot to testing larger fragments (e.g., 40-60 atoms) or a more sophisticated embedding scheme (e.g., electrostatic embedding via point charges) in Phase 1b before proceeding.
- Risques
- The fragmentation error for transition states is significantly larger than for minima, potentially exceeding the 1.5 kcal/mol threshold.Probabilité : mediumParade : Use a larger fragment size for TS structures. If this fails, consider a “fragmentation correction” term trained on a small set of full-system TS calculations.
- The MLFF trained on fragments fails to generalize to the full system due to long-range electrostatic or dispersion interactions not captured in the fragments.Probabilité : mediumParade : Include a fixed electrostatic embedding (e.g., from a background charge distribution) in the fragment calculations. If this is insufficient, the “pivot” strategy will be activated.
- AIMD sampling of the full system is too short to provide a reliable force MAE estimate.Probabilité : lowParade : Extend the AIMD trajectory to 5 ps and use block averaging to estimate the error. We will also compare the MLFF’s potential energy surface against a few static DFT single points.