The output energy_density causes a segfault when computing the xc contribution. This output is also not tested in the testsuite at the moment. A test should be added when this bug is fixed.