Motivated by the observation of comets and asteroids rotating in non-principal axis (NPA) states, we investigate the relaxation of a freely precessing triaxial ellipsoidal rotator towards its lowest-energy spin state. Relaxation of the precession arises from internal dissipative stresses generated by self-gravitation and inertial forces from spin. We develop a general theory to determine the viscoelastic stresses in the rotator, under any linear rheology, for both long-axis (LAM) and short-axis (SAM) modes. By the methods of continuum mechanics, we calculate the power dissipated by the stress field and the viscoelastic material strain which enables us to determine the timescale of the precession dampening. To illustrate how the theory is used, we apply our framework to a triaxial 1I/2017 (`Oumuamua) and 4179 Toutatis under the Maxwell regime. For the former, employing viscoelastic parameters typical of very cold monolithic asteroids renders a dampening timescale longer by a factor of $10^{10}$ and higher than the timescales found in the works relying on the $\,Q$-factor approach, whilst the latter yields a significantly shorter timescale as a consequence of including self-gravitation. We further reduce our triaxial theory to bodies of an oblate geometry and derive a family of relatively simple analytic approximations determining the NPA dampening times for Maxwell rotators, as well as a criterion determining whether self-gravitation is negligible in the relaxation process. Our approximations exhibit a relative error no larger than $0.2\%$, when compared to numerical integration, for close to non-dissipative bodies and $0.002\%$ for highly energy dissipating rotators.

The discovery of Pluto's small moons in the last decade brought attention to the dynamics of the dwarf planet's satellites. With such systems in mind, we study a planar $N$-body system in which all the bodies are point masses, except for a single rigid body. We then present a reduced model consisting of a planar $N$-body problem with the rigid body treated as a 1D continuum (i.e. the body is treated as a rod with an arbitrary mass distribution). Such a model provides a good approximation to highly asymmetric geometries, such as the recently observed interstellar asteroid 'Oumuamua, but is also amenable to analysis. We analytically demonstrate the existence of homoclinic chaos in the case where one of the orbits is nearly circular by way of the Melnikov method, and give numerical evidence for chaos when the orbits are more complicated. We show that the extent of chaos in parameter space is strongly tied to the deviations from a purely circular orbit. These results suggest that chaos is ubiquitous in many-body problems when one or more of the rigid bodies exhibits non-spherical and highly asymmetric geometries. The excitation of chaotic rotations does not appear to require tidal dissipation, obliquity variation, or orbital resonance. Such dynamics give a possible explanation for routes to chaotic dynamics observed in $N$-body systems such as the Pluto system where some of the bodies are highly non-spherical.

We discuss a problem concerning the ultrahigh energy neutrino propagation through the Earth. We present calculation of the neutrino-nucleon cross section at high energies, based on the unifed evolution equation at small x. We also show the solution of the transport equation for different neutrino fluxes originating from active galactic nuclei, gamma ray bursts and top-down model.

In this talk we would like to present recent calculations of the cross sections for the ultrahigh energy neutrinos interacting with nucleons. We briefly present the framework of unifed BFKL/DGLAP equations which resums all the leading log(1/x) effects as well as the subleading terms. The few free parameters which specify the input parton distributions are determined by fitting to HERA deep inelastic data. We then use these parton distributions to calculate the cross section for the neutrino-nucleon interactions at very high energies, up to 10^12 GeV. We do also investigate the attenuation of neutrinos when traversing through the Earth. We use the transport equation which we solve for different fluxes of high energy neutrinos (Active Galactic Nucleai, Gamma Ray Bursts, and top-down models). We study the effects of the regeneration of the neutrino flux via neutral current interaction.

We calculate the cross sections for neutrino interactions with (isoscalar) nuclear targets in the energy domain all the way up to 10^12 GeV. Small x QCD effects are included by using a unified BFKL/DGLAP formalism which embodies non-leading log(1/x) contributions. The few free parameters which specify the input parton distributions are determined by fitting to HERA deep inelastic data. The attenuation of neutrinos transversing the Earth at different nadir angles is calculated for a variety of energy spectra for neutrinos originating from different sources (from Active Galactic Nuclei, Gamma ray bursts, top-down models), as well as for atmospheric neutrinos. For this purpose we solve the transport equation which includes regeneration due to neutral current neutrino interactions, besides attenuation.