In this work we introduce contact Hamiltonian mechanics, an extension of symplectic Hamiltonian mechanics, and show that it is a natural candidate for a geometric description of non-dissipative and dissipative systems. For this purpose we review in detail the major features of standard symplectic Hamiltonian dynamics and show that all of them can be generalized to the contact case. (Alessandro Bravetti), hans@ciencias.unam.mx (Hans Cruz), diego.tapias@nucleares.unam.mx (Diego Tapias) 4 Conclusions and perspectives 29 Appendix A Invariants for the damped parametric oscillator 32 Appendix B Equivalence between the contact Hamilton-Jacobi equation and the contact Hamiltonian equations 34

We present a deterministic algorithm called contact density dynamics that generates any prescribed target distribution in the physical phase space. Akin to the famous model of Nosé and Hoover, our algorithm is based on a non-Hamiltonian system in an extended phase space. However, the equations of motion in our case follow from contact geometry and we show that in general they have a similar form to those of the so-called density dynamics algorithm. As a prototypical example, we apply our algorithm to produce a Gibbs canonical distribution for a one-dimensional harmonic oscillator.

Standard statistical mechanics of conservative systems relies on the symplectic geometry of the phase space. This is exploited to derive Hamilton's equations, Liouville's theorem and to find the canonical invariant measure. In this work we analyze the statistical mechanics of a class of nonconservative systems stemming from contact geometry. In particular, we find out the generalized Hamilton's equations, Liouville's theorem and the microcanonical and canonical measures invariant under the contact flow. Remarkably, the latter measure has a power law density distribution with respect to the standard contact volume form. Finally, we argue on the several possible applications of our results. * bravetti@correo.nucleares.unam.mx † diego.tapias@nucleares.unam.mx 2

Trap models describe glassy dynamics as a stochastic process on a network of configurations representing local energy minima. We study within this class the paradigmatic Barrat-Mézard model, which has Glauber transition rates. Our focus is on the effects of the network connectivity, where we go beyond the usual mean field (fully connected) approximation and consider sparse networks, specifically random regular graphs (RRG). We obtain the spectral density of relaxation rates of the master operator using the cavity method, revealing very rich behaviour as a function of network connectivity c and temperature T . We trace this back to a crossover from initially entropic barriers, resulting from a paucity of downhill directions, to energy barriers that govern the escape from local minima at long times. The insights gained are used to rationalize the relaxation of the energy after a quench from high T , as well as the corresponding correlation and persistence functions.

Motivated by the recent growing interest about the thermodynamic cost of shortcuts to adiabaticity, we consider the cost of driving a classical system by the so-called counterdiabatic driving (CD). To do so, we proceed in three steps: first we review a general definition recently put forward in the literature for the thermodynamic cost of driving a Hamiltonian system; then we provide a new complementary definition of cost, which is of particular relevance for cases where the average excess work vanishes; finally, we apply our general framework to the case of CD. Interestingly, we find that in such a case our results are the exact classical counterparts of those reported by Funo et al. [Phys. Rev. Lett. 118, 100602 (2017)PRLTAO0031-900710.1103/PhysRevLett.118.100602]. In particular we show that a universal trade-off between speed and cost for CD also exists in the classical case. To illustrate our points we consider the example of a time-dependent harmonic oscillator subject to different strategies of adiabatic control.