Reply to "Comment on `Quenches in quantum many-body systems: One-dimensional Bose-Hubbard model reexamined' ''

In his Comment [see preceding Comment, Phys. Rev. A 82, 037601 (2010)] on the paper by Roux [Phys. Rev. A 79, 021608(R) (2009)], Rigol argued that the energy distribution after a quench is not related to standard statistical ensembles and cannot explain thermalization. The latter is proposed to stem from what he calls the eigenstate thermalization hypothesis and which boils down to the fact that simple observables are expected to be smooth functions of the energy. In this Reply, we show that there is no contradiction or confusion between the observations and discussions of Roux and the expected thermalization scenario discussed by Rigol. In addition, we emphasize a few other important aspects, in particular the definition of temperature and the equivalence of ensemble, which are much more difficult to show numerically even though we believe they are essential to the discussion of thermalization. These remarks could be of interest to people interested in the interpretation of the data obtained on finite-size systems.Comment: 3 page

Viscosity calculated in simulations of strongly-coupled dusty plasmas with gas friction

A two-dimensional strongly-coupled dusty plasma is modeled using Langevin and frictionless molecular dynamical simulations. The static viscosity $\eta$ and the wave-number-dependent viscosity $\eta(k)$ are calculated from the microscopic shear in the random motion of particles. A recently developed method of calculating the wave-number-dependent viscosity $\eta(k)$ is validated by comparing the results of $\eta(k)$ from the two simulations. It is also verified that the Green-Kubo relation can still yield an accurate measure of the static viscosity $\eta$ in the presence of a modest level of friction as in dusty plasma experiments.Comment: 6 pages, 3 figures, Physics of Plasmas invited pape

Noninteracting Fermions in infinite dimensions

Usually, we study the statistical behaviours of noninteracting Fermions in finite (mainly two and three) dimensions. For a fixed number of fermions, the average energy per fermion is calculated in two and in three dimensions and it becomes equal to 50 and 60 per cent of the fermi energy respectively. However, in the higher dimensions this percentage increases as the dimensionality increases and in infinite dimensions it becomes 100 per cent. This is an intersting result, at least pedagogically. Which implies all fermions are moving with Fermi momentum. This result is not yet discussed in standard text books of quantum statistics. In this paper, this fact is discussed and explained. I hope, this article will be helpful for graduate students to study the behaviours of free fermions in generalised dimensionality.Comment: To appear in European Journal of Physics (2010

Nuclear condensation and the equation of state of nuclear matter

The isothermal compression of a dilute nucleonic gas invoking cluster degrees of freedom is studied in an equilibrium statistical model; this clusterized system is found to be more stable than the pure nucleonic system. The equation of state (EoS) of this matter shows features qualitatively very similar to the one obtained from pure nucleonic gas. In the isothermal compression process, there is a sudden enhancement of clusterization at a transition density rendering features analogous to the gas-liquid phase transition in normal dilute nucleonic matter. Different observables like the caloric curves, heat capacity, isospin distillation, etc. are studied in both the models. Possible changes in the observables due to recently indicated medium modifications in the symmetry energy are also investigated.Comment: 18 pages and 11 figures. Phys. Rev. C (in press

Limits of the equivalence of time and ensemble averages in shear flows

In equilibrium systems, time and ensemble averages of physical quantities are equivalent due to ergodic exploration of phase space. In driven systems, it is unknown if a similar equivalence of time and ensemble averages exists. We explore effective limits of such convergence in a sheared bubble raft using averages of the bubble velocities. In independent experiments, averaging over time leads to well converged velocity profiles. However, the time-averages from independent experiments result in distinct velocity averages. Ensemble averages are approximated by randomly selecting bubble velocities from independent experiments. Increasingly better approximations of ensemble averages converge toward a unique velocity profile. Therefore, the experiments establish that in practical realizations of non-equilibrium systems, temporal averaging and ensemble averaging can yield convergent (stationary) but distinct distributions.Comment: accepted to PRL - final figure revision

Ultracold heteronuclear molecules and ferroelectric superfluids

We analyze the possibility of a ferroelectric transition in heteronuclear molecules consisting of Bose-Bose, Bose-Fermi or Fermi-Fermi atom pairs. This transition is characterized by the appearance of a spontaneous electric polarization below a critical temperature. We discuss the existence of a ferroelectric Fermi liquid phase for Fermi molecules and the existence of a ferroelectric superfluid phase for Bose molecules characterized by the coexistence of ferroelectric and superfluid orders. Lastly, we propose an experiment to detect ferroelectric correlations through the observation of coherent dipole radiation pulses during time of flight.Comment: 4 pages and 3 figure

Itinerant ferromagnetism in an atomic Fermi gas: Influence of population imbalance

We investigate ferromagnetic ordering in an itinerant ultracold atomic Fermi gas with repulsive interactions and population imbalance. In a spatially uniform system, we show that at zero temperature the transition to the itinerant magnetic phase transforms from first to second order with increasing population imbalance. Drawing on these results, we elucidate the phases present in a trapped geometry, finding three characteristic types of behavior with changing population imbalance. Finally, we outline the potential experimental implications of the findings.Comment: 10 pages, 4 figures, typos added, references adde

Equation of state of an interacting Bose gas at finite temperature: a Path Integral Monte Carlo study

By using exact Path Integral Monte Carlo methods we calculate the equation of state of an interacting Bose gas as a function of temperature both below and above the superfluid transition. The universal character of the equation of state for dilute systems and low temperatures is investigated by modeling the interatomic interactions using different repulsive potentials corresponding to the same s-wave scattering length. The results obtained for the energy and the pressure are compared to the virial expansion for temperatures larger than the critical temperature. At very low temperatures we find agreement with the ground-state energy calculated using the diffusion Monte Carlo method.Comment: 7 pages, 6 figure

Microcanonical Origin of the Maximum Entropy Principle for Open Systems

The canonical ensemble describes an open system in equilibrium with a heat bath of fixed temperature. The probability distribution of such a system, the Boltzmann distribution, is derived from the uniform probability distribution of the closed universe consisting of the open system and the heat bath, by taking the limit where the heat bath is much larger than the system of interest. Alternatively, the Boltzmann distribution can be derived from the Maximum Entropy Principle, where the Gibbs-Shannon entropy is maximized under the constraint that the mean energy of the open system is fixed. To make the connection between these two apparently distinct methods for deriving the Boltzmann distribution, it is first shown that the uniform distribution for a microcanonical distribution is obtained from the Maximum Entropy Principle applied to a closed system. Then I show that the target function in the Maximum Entropy Principle for the open system, is obtained by partial maximization of Gibbs-Shannon entropy of the closed universe over the microstate probability distributions of the heat bath. Thus, microcanonical origin of the Entropy Maximization procedure for an open system, is established in a rigorous manner, showing the equivalence between apparently two distinct approaches for deriving the Boltzmann distribution. By extending the mathematical formalism to dynamical paths, the result may also provide an alternative justification for the principle of path entropy maximization as well.Comment: 12 pages, no figur

Interaction effects and quantum phase transitions in topological insulators

We study strong correlation effects in topological insulators via the Lanczos algorithm, which we utilize to calculate the exact many-particle ground-state wave function and its topological properties. We analyze the simple, noninteracting Haldane model on a honeycomb lattice with known topological properties and demonstrate that these properties are already evident in small clusters. Next, we consider interacting fermions by introducing repulsive nearest-neighbor interactions. A first-order quantum phase transition was discovered at finite interaction strength between the topological band insulator and a topologically trivial Mott insulating phase by use of the fidelity metric and the charge-density-wave structure factor. We construct the phase diagram at $T = 0$ as a function of the interaction strength and the complex phase for the next-nearest-neighbor hoppings. Finally, we consider the Haldane model with interacting hard-core bosons, where no evidence for a topological phase is observed. An important general conclusion of our work is that despite the intrinsic nonlocality of topological phases their key topological properties manifest themselves already in small systems and therefore can be studied numerically via exact diagonalization and observed experimentally, e.g., with trapped ions and cold atoms in optical lattices.Comment: 13 pages, 12 figures. Published versio