Hermes,
With your assistance I have added the Gamow suppression factors; many thanks!
QuoteDisplay More% python scripts/reactions.py --model induced-decay Am --gamow --unstable | sort | grep "\["
e- + 230Am → e- + 4He + 226Np + 7725 keV [22] α, →α, →β+
e- + 231Am → e- + 4He + 227Np + 7455 keV [23] α, →α, →β+
e- + 232Am → e- + 4He + 228Np + 7245 keV [24] α, →α, →β+SF, →ε
e- + 233Am → e- + 4He + 229Np + 7055 keV [25] α, →α
e- + 234Am → e- + 4He + 230Np + 6795 keV [26] α, →β+
e- + 235Am → e- + 4He + 231Np + 6575 keV [27] α, →α, →β+
e- + 236Am → e- + 4He + 232Np + 6255 keV [28] α, →α, →β+
e- + 237Am → e- + 4He + 233Np + 6195 keV [29] α, →α, →β+
e- + 238Am → e- + 4He + 234Np + 6038 keV [29] α, →β+
e- + 239Am → e- + 4He + 235Np + 5922 keV [30] α, →α, →ε
e- + 240Am → e- + 4He + 236Np + 5707 keV [31] α, →α, →β-, →ε
e- + 241Am → e- + 4He + 237Np + 5638 keV [31] trace, α, →SF, →α
e- + 242Am → e- + 4He + 238Np + 5588 keV [32] α, →β-
e- + 243Am → e- + 4He + 239Np + 5439 keV [33] trace, α, →α, →β-
e- + 244Am → e- + 4He + 240Np + 5138 keV [34] α, →β-
e- + 245Am → e- + 4He + 241Np + 5217 keV [34] α, →α, →β-
e- + 246Am → e- + 4He + 242Np + 5150 keV [34] α, →β-
e- + 247Am → e- + 4He + 243Np + 4845 keV [36] α, →β-
e- + 248Am → e- + 4He + 244Np + 4935 keV [36] α, →β-
e- + 249Am → e- + 4He + 245Np + 4725 keV [37] α
I think you must have modeled the screening by adjusting the atomic number of the heavier daughter (the "Z" argument), but it would be good to get confirmation of this.