On the calculation of first‐order properties: Expectation value versus energy derivative approach