And because the problem is full of coupled calculations, one poorly described defect might throw out your quantitative thermodynamics completely

I think for me, the main issue that makes it hard to implement for automation and reproducibility is handling edge cases, which are common. Even doped is far from fool-proof and requires a lot of on-the-fly application of domain knowledge to get useful results.