We propose the use of partial dynamical symmetry (PDS) as a selection criterion for higher-order terms in situations when a prescribed symmetry is obeyed by some states and is strongly broken in others. The procedure is demonstrated in a first systematic classification of many-body interactions with SU(3) PDS that can improve the description of deformed nuclei. As an example, the triaxial features of the nucleus 156 Gd are analyzed. Many-body forces play an important role in quantum many-body systems [1]. They appear either at a fundamental level or as effective interactions which arise due to restriction of degrees of freedom and truncation of model spaces. A known example is the structure of light nuclei, where two-nucleon interactions are insufficient to achieve an accurate description and higher-order interactions between the nucleons must be included [2]. Given the difficulty in constraining the nature of such higher-order terms from experiments, one is faced with the problem of their determination. One way, currently the subject of active research [3], is to determine them from chiral effective field theory applied to quantum chromodynamics. This establishes a hierarchy of inter-nucleon interactions according to their order. In light-medium nuclei, these interactions serve as input for ab-initio methods (e.g., the no-core shell model (NCSM) [4]) to generate, by means of similarity transformations, A-body effective Hamiltonians in computational tractable model spaces.The situation is more complex in heavy nuclei, where ab-initio methods are limited by the enormous increase in size of the model spaces required to accommodate correlated collective motion of many nucleons. One possible approach to circumvent this problem, is to augment the NCSM method through a symplectic symmetry-adapted choice of basis [5]. A second approach is to employ energy density functionals and incorporate beyond mean-field effects by mapping to collective Hamiltonians [6], e.g., the interacting boson model (IBM) [7]. In both approaches the Hilbert spaces are based on particular dynamical algebras which lead to a dramatic reduction of the basis dimension. Nevertheless, even with such simplification, the number of possible interactions in the effective Hamiltonians grows rapidly with their order, and a selection criterion is called for. In this Rapid Communication, we suggest a method to select possible higher-order terms which is based on the idea of partial dynamical symmetry (PDS).The concept of PDS [8] is a generalization of that of a dynamical symmetry (DS) [9] where the conditions of the latter (solvability of the complete spectrum, existence of exact quantum numbers for all eigenstates, and predetermined structure of the eigenfunctions) are relaxed and apply to only part of the eigenstates and/or of the quantum numbers. PDSs have been identified in various dynamical systems involving bosons and fermions (for a review, see Ref.[8]). They play a role in diverse phenomena including nuclear and molecular spectroscopy [10][11][12], q...