We compute the correlation of the net baryon number with the electric charge (χBQ) for an interacting hadron gas using the S-matrix formulation of statistical mechanics. The observable χBQ is particularly sensitive to the details of the pion-nucleon interaction, which are consistently incorporated in the current scheme via the empirical scattering phase shifts. Comparing to the recent lattice QCD studies in the (2 + 1)-flavor system, we find that the natural implementation of interactions and the proper treatment of resonances in the S-matrix approach lead to an improved description of the lattice data over that obtained in the hadron resonance gas model. 25.75.Ld, 12.38.Mh, 24.10.Nz Introduction.-Recent lattice QCD (LQCD) results on the equation of states and the fluctuations of conserved charges provide a very detailed description of the QCD thermal medium [1][2][3][4][5]. In particular the local fluctuations of conserved charges can be probed by appropriate combinations of mixed susceptibilities. An accurate determination of these quantities is also needed to reliably extend the LQCD calculations to finite densities using the Taylor's expansion scheme [6].Confinement dictates that hadrons, instead of quarks and gluons, fill the physical spectrum of QCD, while the spontaneous breaking of chiral symmetry makes pions exceptionally light due to their role as (pseudo-) Goldstone bosons. We thus expect that at low temperatures the partition function can be effectively described by an interacting gas of low-mass hadrons such as pions, kaons, and nucleons.A well-known effective approach which adopts the hadronic degrees of freedom in describing the thermodynamics of strongly interacting matter is the hadron resonance gas (HRG) model. This model assumes that resonance formation dominates the interactions of the confined phase, and as a first approximation, treats the resonances as an ideal gas. The approach gives a satisfactory description of the particle yields measured in heavy ion collisions [7][8][9][10][11][12][13][14], and is capable of providing an overall successful interpretation of LQCD results on bulk properties below the transition temperature [1][2][3][4][5][15][16][17].Nevertheless, the HRG model also makes some simplifying assumptions which are not necessarily consistent with the known hadron physics. Some of the problematic cases include the zero-width treatment of broad resonances [18-20] (e.g. the σ-and κ-meson), and the neglect of non-resonant contributions from both attractive and repulsive channels in computing the thermal observables [21].Very precise information about the hadronic interactions has emerged from the impressive volume of experimental data [22], carefully analyzed by theory such as