Chaos in Nonequilibrium Two-Temperature $(T_x, T_y)$ Nosé-Hoover Cell Models

We revisit a two-temperature Nosé-Hoover wanderer particle embedded in a two-dimensional periodic 2x2 cell with four smooth repulsive corners at $(x,y) = (\pm 1, \pm 1)$ to explore chaos with anisotropic thermostatting. The model employs separate thermostats in the x and y directions, enabling controlled deviations from equilibrium. By integrating the full six-dimensional equations of motion and computing the complete Lyapunov spectrum, we confirm chaos and quantify phase-space contraction with high numerical precision. The total contraction rate, interpreted as entropy production, grows nonlinearly with the thermostat anisotropy and follows a superquadratic power law, $Λ\propto -δ^{2.44}$, deviating from linear-response theory. The approximate Kaplan-Yorke dimension reveals a fractal attractor that concentrates as $|T_x - T_y|$ increases. Momentum statistics show significant non-Gaussian behavior under strong driving. Despite its dissipative nature, the model remains strictly time-reversible, offering a pedagogically rich example of microscopic reversibility coexisting with macroscopic entropy production.