Let us take an Hamiltonian $H(\xi)$ which depends on a set of parameters $\xi$, and assume that the matrix elements $h_{ij}(\xi)$ of the Hamiltonian are smooth complex functions of the parameters $\xi$ (i.e., each $h_{ij}(\xi)$ has continuous derivatives of any order in the parameters $\xi$). Assume also that such parameters are real ($h_{ij}(\xi)$ are complex in general) and each of them can vary on the whole real line $\mathbb{R}$, and that the Hamiltonian has a fixed rank $n$.
Under these assumptions, the eigenvalues $E_i(\xi)$ of the Hamiltonian $H(\xi)$ can still have a discontinuous derivative in the presence of level crossings. However, is it always possible to "label" the eigenvalues in such a way that all eigenvalues are smooth in the parameters? (or at least continuous with continuous derivative)
If this is not the case, is it at least the total energy a smooth function of the parameter $\xi$, under certain assumptions? One can define the total energy as $E_T(\xi)=\sum_{E_i<E_F} E_i(\xi)$, i.e., the sum of the energy levels below a certain threshold $E_F$ (e.g., Fermi energy in a fermion system), and assume that $E_i(\xi)\neq E_f$ (if any energy level crosses the energy $E_F$ it is clear that the total energy can have a discontinuous derivative).