We have developed a formalism to study the secular evolution of a binary system, including the interaction due to the tides that each body raises on the other one. We also consider the influence of the $J_2$-related secular terms on the orbital evolution and the triaxiality-caused torque on the rotational evolution, both associated only to one of the bodies. We apply these set of equations to the study of the orbital and rotational evolution of a binary system composed by a rocky planet and its host star in order to characterize the dynamical evolution, particularly near spin-orbit resonances. We use the equations of motion that give the time evolution of the orbital elements and the spin rates of each body to study the time evolution of the Kepler-21 system as an example of application of the formalism developed. We have obtained a set of equations of motion without singularities for vanishing eccentricities and inclinations. This set gives, on one hand, the time evolution of the orbital elements due to the tidal potentials generated by both partners of the system and the triaxiality of one of them. On the other hand, it gives the time evolution of the stellar spin rate due to the corresponding tidal torque and of the planet's rotation angle due to both the tidal and triaxiality-caused torque. We found that, for the parameters and the initial conditions explored, the tidally and triaxiality induced modifications of the tidal modes is more significative than expected and that the time of tidal synchronization depends strongly on the values of the rheological parameters.

In several previous investigations we presented models of triaxial stellar systems, both cuspy and non cuspy, that were highly stable and harboured large fractions of chaotic orbits. All our models had been obtained through cold collapses of initially spherical $N$--body systems, a method that necessarily results in models with strongly radial velocity distributions. Here we investigate a different method that was reported to yield cuspy triaxial models with virtually no chaos. We show that such result was probably due to the use of an inadequate chaos detection technique and that, in fact, models with significant fractions of chaotic orbits result also from that method. Besides, starting with one of the models from the first paper in this series, we obtained three different models by rendering its velocity distribution much less radially biased (i.e., more isotropic) and by modifying its axial ratios through adiabatic compression. All three models yielded much higher fractions of regular orbits than most of those from our previous work. We conclude that it is possible to obtain stable cuspy triaxial models of stellar systems whose velocity distribution is more isotropic than that of the models obtained from cold collapses. Those models still harbour large fractions of chaotic orbits and, although it is difficult to compare the results from different models, we can tentatively conclude that chaoticity is reduced by velocity isotropy.

In the first paper of this series we used the N--body method to build a dozen cuspy (gamma ~ 1) triaxial models of stellar systems, and we showed that they were highly stable over time intervals of the order of a Hubble time, even though they had very large fractions of chaotic orbits (more than 85 per cent in some cases). The models were grouped in four sets, each one comprising models morphologically resembling E2, E3, E4 and E5 galaxies, respectively. The three models within each set, although different, had the same global properties and were statistically equivalent. In the present paper we use frequency analysis to classify the regular orbits of those models. The bulk of those orbits are short axis tubes (SATs), with a significant fraction of long axis tubes (LATs) in the E2 models that decreases in the E3 and E4 models to become negligibly small in the E5 models. Most of the LATs in the E2 and E3 models are outer LATs, but the situation reverses in the E4 and E5 models where the few LATs are mainly inner LATs. As could be expected for cuspy models, most of the boxes are resonant orbits, i.e., boxlets. Nevertheless, only the (x, y) fishes of models E3 and E4 amount to about 10 per cent of the regular orbits, with most of the fractions of the other boxlets being of the order of 1 per cent or less.