Abstract:We consider the relaxation of an initial non-equilibrium state in a one-dimensional fluid of hard rods. Since it is an interacting integrable system, we expect it to reach the Generalized Gibbs Ensemble (GGE) at long times for generic initial conditions. Here we show that there exist initial conditions for which the system does not reach GGE even at very long times and in the thermodynamic limit. In particular, we consider an initial condition of uniformly distributed hard-rods in a box with the left half having particles with a singular velocity distribution (all moving with unit velocity) and the right half particles in thermal equilibrium. We find that the density profile for the singular component does not spread to the full extent of the box and keeps moving with a fixed effective speed at long times. We show that such density profiles can be well described by the solution of the Euler equations almost everywhere except at the location of the shocks, where we observe slight discrepancies due to dissipation arising from the initial fluctuations of the thermal background. To demonstrate this effect of dissipation analytically, we consider a second initial condition with a single particle at the origin with unit velocity in a thermal background. We find that the probability distribution of the position of the unit velocity quasi-particle has diffusive spreading which can be understood from the solution of the Navier–Stokes (NS) equation of the hard rods. Finally, we consider an initial condition with a spread in velocity distribution for which we show convergence to GGE. Our conclusions are based on molecular dynamics simulations supported by analytical arguments.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the thermalization and hydrodynamic behavior of a one - dimensional hard - rod fluid system under specific non - equilibrium initial conditions. Specifically, the researchers explored the following issues:
1. **Whether the Generalized Gibbs Ensemble (GGE) is reached**: For certain initial conditions, whether the system will evolve over time and finally reach the Generalized Gibbs Ensemble (GGE). In particular, when the initial condition is that the particles in the left half have a single - velocity distribution (all particles move at unit velocity), and the particles in the right half are in thermal equilibrium, what is the behavior of the system.
2. **The influence of dissipative effects**: The researchers also explored the specific influence of dissipative effects on the evolution of the system. For example, in some cases, the density profile will not completely diffuse throughout the box, but will move at a certain effective velocity, and a slight difference due to the initial thermal background fluctuations is observed at the shock position.
### Specific research content
The paper studies these problems through the following initial conditions:
- **Initial condition A**: The particles in the left half have a single velocity \(v_0 = 1\), while the particles in the right half select velocities from the Maxwell distribution.
- **Initial condition B**: Place a particle with a given velocity \(v_0 = 1\) at the origin, and the remaining space is uniformly distributed with background particles having a Maxwell velocity distribution.
- **Initial condition C**: Free expansion from a half - filled box, that is, the left half is uniformly distributed with particles, and the right half is empty.
### Main findings
- **For initial condition A**:
- The Euler equation solution agrees well with the molecular dynamics simulation results in most regions, but there are obvious differences at the shock position, which are attributed to the dissipative effects described by the Navier - Stokes (NS) terms.
- The densities of the two initial components will not become uniform throughout the system in the long term, but will each move at a constant effective velocity \(v_{\text{eff}}\), so the system will not reach the GGE state.
- **For initial condition B**:
- Provides an analytical understanding based on the NS equation solution, indicating that in this case, the system will reach the GGE state in the thermodynamic limit.
- **For initial condition C**:
- The evolution of the density and momentum density profiles is completely described by the Euler GHD, and since the discontinuity of the initial density profile disappears very early, the influence of the NS correction term is too small to be observed.
- After a long time, the system evolves to a state consistent with GGE.
### Conclusion
The research results show that for some non - equilibrium initial conditions, the one - dimensional hard - rod system will not reach the Generalized Gibbs Ensemble (GGE) state, especially under initial condition A. In addition, the dissipative effects significantly affect the evolution behavior of the system at the shock position. These findings contribute to a deeper understanding of the thermalization and hydrodynamic characteristics in interacting integrable systems.