Our collective knowledge of nuclear cross sections is recorded as resonance parameters in nuclear data libraries. To evaluate these parameters, campaigns of measurements are fitted with a parametric model of nuclear cross sections called R-matrix theory. R-matrix theory can parametrize the energy dependence of nuclear cross sections in different ways. Historically, the Wigner-Eisenbud parametrization has been used, because its resonance parameters are real-valued and there is a one-to-one correspondance between each cross-section resonance (called levels) and the Wigner-Eisenbud resonance energies (the poles of the R-matrix). This means that for each level, every nuclear interaction channel has one resonance width and one resonance energy. The drawback of the Wigner-Eisenbud parametrization is that it introduces in each channel an arbitrary boundary parameter upon which all other resonance parameters depend. Thus, if two evaluations of the same experiments are done with different boundary parameters, both will yield the same fit, but a different set of resonance parameters. This is a challenge for nuclear data libraries.To overcome such arbitrarity, Brune proposed an alternative parametrization preserving most benefits of the Wigner-Eisenbud parameters while eliminating the arbitrary boundary parameter. Consequently, the community is considering converting all nuclear data libraries to Brune resonance parameters. The Brune resonance energies are the poles of the Brune alternative level matrix, and Brune proved that, above the channel threshold energy, there is one pole per level (or resonance).In this article, we unveil the existence of more Brune parameters than previously thought. Below the channel threshold, we prove the theoretical existence of two types of additional shadow polesbranch shadow poles and analytic shadow poles -depending on what convention is chosen to continue the R-matrix operators to complex wavenumbers (to do so we establish the first derivations of the Mittag-Leffler expansions of external region R-matrix operators). This entails there are more Brune poles (or Brune resonance energies) than levels. Importantly, we also prove that choosing any subset of Brune poles will yield the same cross sections than using the entire set of Brune poles, as long it is bigger or equal the number of levels (resonances). In practice, this means that shadow Brune poles can safely be discarded from the new nuclear data libraries.Many isotopes in nuclear cross section libraries are evaluated under the Reich-Moore approximation, which introduces complex resonance energies to eliminate certain channels. We generalize Brune's parameterization to encompass both the Reich-Moore approximation and the addional shadow poles. In the process, we show that all Brune parameters values depend on what convention we choose to continue the R-matrix operators to complex wavenumbers. To convert nuclear data libraries to Brune parameters, the nuclear scientists community must thus first decide on such a convention. The authors a...