Regarding 1: The expansion coefficients of a multiconfigurational wavefunction depend on the choice of single-particle basis; Ref. 2 shows a particularly stark example in which delocalized vs localized orbitals imply qualitatively different assessments of multiconfigurational character. It has also been shown that most mono-nuclear transition metal coordination complexes are qualitatively well-described by a single determinant plus an appropriate level of dynamic correlation, 3 and that exchange-coupled polynuclear transition metal compounds can often be described by a single spin-adapted configuration state function. Regarding 3: In our view, the characterization of a system as strongly correlated should not depend on the partitioning of a Hilbert space into distinguishable subspaces, as is often required to compute metrics such as entanglement entropy. As in 1, orbital entanglement entropy metrics have been shown to be strongly dependent on the choice of single-particle basis employed. Regarding 4: While the overlap of a given wavefunction with FCI, ⟨Φ|Ψ F CI ⟩, is useful for finite systems such as molecules (as it reports on the percentage of the FCI Hilbert space that is relevant to that wavefunction), the focus on the case where |Φ⟩ is the mean-field state is not relevant in the thermodynamic limit in light of the well-known "orthogonality catastrophe" initially described by Van Vleck: Given N infinitely separated He atoms, and assuming that ⟨Φ i |Ψ F CI ⟩ = 1-ϵ for each of the He atoms, then the overlap for the composite system (1 -ϵ) N will decay exponentially with N , even for weakly correlated systems.