Done. See output-h-chain-peierls.out.