The $Z_{bus}$ matrix is a full matrix, whereas $Y_{bus}$ is sparse.

$Z_{bus}$ is calculated by inverting $Y_{bus}$ or through $Z_{bus}$ building algorithms.

For large power systems, $Z_{bus}$ building algorithms are preferred to avoid direct matrix inversion.

Short circuit calculations require the Thevenin impedance at the fault location, which is explicitly available in the diagonal elements of $Z_{bus}$.

Direct calculation of fault currents at any bus.

Accounts for all system branches effectively.

Easy to incorporate changes in network topology.

Computational storage requirement is high ($n²$ elements).

Building $Z_{bus}$ for very large systems is computationally intensive.

Symmetrical fault analysis.

Unsymmetrical fault analysis (using sequence networks).

Relay coordination studies.

$Y_{BUS}$ is typically used for Newton-Raphson or Gauss-Seidel load flow due to its sparsity property, which saves memory.

Option A is incorrect because while $Y_{BUS}$ is essential for power flow, it is not the standard choice for direct fault current calculation as it requires matrix inversion.