corrosim — multiscale corrosion-inhibitor report

Substrate Cu(111)  |  Medium 1 M HNO3  |  DFT level B3LYP/6-311++G(d,p) (ddCOSMO:water)  |  Generated 2026-07-18 19:37

The molecules modelled here are screened in silico as representative candidates, not a verified analysis of any real sample. Confirm the actual composition experimentally (e.g. by LC-MS/GC-MS) before drawing firm conclusions.

1. Overview

corrosim pipeline

2. Summary & ranking

Of the 4 molecules screened, 1-phenyl-5-mercaptotetrazole ranks highest on the composite electronic score (+0.88), which averages three independent axes — a small HOMO–LUMO gap (5.11 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 ≈ -10 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).

1-Phenyl-5-mercaptotetrazole5-Aminotetrazole5-PhenyltetrazoleTetrazole
Gap (eV)5.107 5.445.2687.186
ΔN0.0880.1530.108-0.039
Dipole (D)6.8323.1092.7722.946
E_ads (kJ/mol)-10.26 -4.48-9.09-3.76
Cu–O (Å)
Score0.88 0.3180.109-1.308
Ranking basisTop candidate
force-field geometry, neutral1-Phenyl-5-mercaptotetrazole
force-field geometry, pH-weighted1-Phenyl-5-mercaptotetrazole

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

Modelled molecules
Frontier-orbital energies vs the metal work function

3.1 Frontier-orbital isosurfaces (HOMO / LUMO)

tetrazole HOMO
5-aminotetrazole HOMO
5-phenyltetrazole HOMO
1-phenyl-5-mercaptotetrazole HOMO
tetrazole LUMO
5-aminotetrazole LUMO
5-phenyltetrazole LUMO
1-phenyl-5-mercaptotetrazole LUMO
Reactivity descriptors
Protonation effect (DFT-optimised cations)

3.2 Full descriptor table (neutral, aqueous)

DescriptorTetrazole5-Aminotetrazole5-Phenyltetrazole1-Phenyl-5-mercaptotetrazole
FormulaCH2N4CH3N5C7H6N4C7H6N4S
HOMO (eV)-8.755-6.784-6.977-7.026
LUMO (eV)-1.572-1.329-1.683-1.918
Gap ΔE (eV)7.1835.4555.2945.108
η hardness (eV)3.5912.7282.6472.554
σ softness (1/eV)0.2780.3670.3780.392
χ electronegativity (eV)5.1644.0564.334.472
ω electrophilicity (eV)3.7123.0163.5413.915
ΔN-0.0310.1620.1150.092
E_back-donation (eV)-0.898-0.682-0.662-0.638
Dipole (D)2.9463.1092.7726.832
TNC-0.608-1.0-3.031-1.83
E_ads (kJ/mol)-3.76-4.48-9.09-10.26

3.3 Species in the acidic medium (protonated cation)

DescriptorTetrazole+H+5-Aminotetrazole+H+5-Phenyltetrazole+H+1-Phenyl-5-mercaptotetrazole+H+
FormulaCH3N4+CH4N5+C7H7N4+C7H7N4S+
HOMO (eV)-10.733-8.127-7.736-7.704
LUMO (eV)-3.457-3.167-3.301-2.619
Gap ΔE (eV)7.2764.964.4355.085
η hardness (eV)3.6382.482.2172.542
σ softness (1/eV)0.2750.4030.4510.393
χ electronegativity (eV)7.0955.6475.5185.162
ω electrophilicity (eV)6.9186.4296.8675.24
ΔN-0.296-0.143-0.13-0.044
E_back-donation (eV)-0.909-0.62-0.554-0.636
Dipole (D)
TNC-0.299-0.843-3.228-1.487

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

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

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

DescriptorTetrazole5-Aminotetrazole5-Phenyltetrazole1-Phenyl-5-mercaptotetrazole
FormulaCH2N4CH3N5C7H6N4C7H6N4S
HOMO (eV)-8.816-6.825-7.0-7.046
LUMO (eV)-1.63-1.385-1.732-1.939
Gap ΔE (eV)7.1865.445.2685.107
η hardness (eV)3.5932.722.6342.554
σ softness (1/eV)0.2780.3680.380.392
χ electronegativity (eV)5.2234.1054.3664.493
ω electrophilicity (eV)3.8113.1213.6433.955
ΔN-0.0390.1530.1080.088
E_back-donation (eV)-0.898-0.68-0.658-0.638
Dipole (D)2.9463.1092.7726.832
TNC-0.598-0.995-3.037-1.82

3.5 Local reactivity (Fukui)

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

tetrazole — condensed Fukui
5-aminotetrazole — condensed Fukui
5-phenyltetrazole — condensed Fukui
1-phenyl-5-mercaptotetrazole — condensed Fukui

3.6 Electrostatic-potential (ESP) map

tetrazole — ESP map
5-aminotetrazole — ESP map
5-phenyltetrazole — ESP map
1-phenyl-5-mercaptotetrazole — ESP map

4. Monte Carlo adsorption

tetrazole — best pose
5-aminotetrazole — best pose
5-phenyltetrazole — best pose
1-phenyl-5-mercaptotetrazole — best pose
tetrazole — MC annealing
5-aminotetrazole — MC annealing
5-phenyltetrazole — MC annealing
1-phenyl-5-mercaptotetrazole — MC annealing

5. Brownian MD — Cu–donor RDF

tetrazole — Cu–donor RDF
5-aminotetrazole — Cu–donor RDF
5-phenyltetrazole — Cu–donor RDF
1-phenyl-5-mercaptotetrazole — Cu–donor RDF

6. Method & caveats

DFT level: B3LYP/6-311++G(d,p) (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.