corrosim — multiscale corrosion-inhibitor report

Substrate Cu(111)  |  Medium 1 M HCl  |  DFT level B3LYP/def2-SVP (ddCOSMO:water)  |  Generated 2026-07-18 19:37

1. Overview

2. Summary & ranking

Of the 5 molecules screened, pyrazolylnucleoside methoxy ranks highest on the composite electronic score (+1.66), which averages three independent axes — a small HOMO–LUMO gap (5.46 eV), the Lukovits ΔN and the molecular dipole — the strongest all-round electron-donating profile of the set. It adsorbs flat on Cu in the physisorption regime (E_ads ≈ -9 kJ/mol). The lead holds across all 2 descriptor bases evaluated (geometry and protonation). This is a computational screening prediction requiring electrochemical confirmation (see the Method section and caveats).

Lead robust: the same candidate tops all 2 descriptor bases (geometry and protonation).

Scored on three independent axes of the canonical basis (force-field geometry, pH-weighted; z-scored; higher = stronger): the HOMO–LUMO gap, the Lukovits ΔN (electron donation to the metal) and the dipole (a weak tie-breaker — its direction vs efficiency is disputed). The geometry / protonation panels below are a sensitivity check, not separate rankings. E_ads and the Cu–O distance validate the lead, they do not enter the score. Full method: docs/pipeline.md.

3. DFT electronic descriptors

3.1 Frontier-orbital isosurfaces (HOMO / LUMO)

3.2 Full descriptor table (neutral, aqueous)

3.3 Species in the acidic medium (protonated cation)

Protonated +1 cation descriptors in 1 M HCl; a component of the pH-weighted canonical basis (see the Summary), shown here on its own.

3.4 Speciation in 1 M HCl (pH ≈ -0.0)

0% neutral / 100% protonated at this pH — the protonated form dominates. Population-weighted descriptors (the speciation axis of the ranking ensemble; the headline ranks on the canonical basis):

3.5 Local reactivity (Fukui)

Strongest electron-donating oxygens (highest f⁻) per molecule:

• Pyrazolylnucleoside methyl: O3 (f⁻=0.039), O22 (f⁻=0.002), O0 (f⁻=0.001)
• Pyrazolylnucleoside methoxy: O17 (f⁻=0.188), O3 (f⁻=0.009), O23 (f⁻=0.005)
• Pyrazolylnucleoside fluoro: O3 (f⁻=0.001), O22 (f⁻=0.000), O0 (f⁻=0.000)
• Pyrazolylnucleoside chloro: O3 (f⁻=0.011), O22 (f⁻=0.003), O0 (f⁻=0.001)
• Pyrazolylnucleoside bromo: O3 (f⁻=0.009), O22 (f⁻=0.003), O0 (f⁻=0.001)

3.6 Electrostatic-potential (ESP) map

4. Monte Carlo adsorption

5. Brownian MD — Cu–donor RDF

6. Method & caveats

DFT level: B3LYP/def2-SVP (ddCOSMO:water). Global descriptors come from the frontier-orbital energies via Koopmans' theorem; ΔN uses the metal work function with η(metal) = 0. The Monte-Carlo and Brownian-MD stages use a classical van-der-Waals adsorption model — a physics-based screening surrogate, not a substitute for periodic DFT or for electrochemical validation. The composite ranking is a heuristic that orders candidates; it does not prove inhibition. Full methodology: docs/pipeline.md; validation record: docs/validation.md.