/// <summary>
        /// Builds Hamiltonian from Broombridge if data is available.
        /// </summary>
        internal static OrbitalIntegralHamiltonian ToOrbitalIntegralHamiltonian(ProblemDescription broombridge)
        {
            // Add the identity terms
            var identityterm = broombridge.CoulombRepulsion.Value + broombridge.EnergyOffset.Value;
            var hamiltonian  = V0_1.ToOrbitalIntegralHamiltonian(broombridge.Hamiltonian);

            hamiltonian.Add(new OrbitalIntegral(), identityterm);
            return(hamiltonian);
        }