We study the density-density response function of a collection of charged massive Dirac particles and present analytical expressions for the dynamical polarization function in one, two and three dimensions. The polarization function is then used to find the dispersion of the plasmon modes, and electrostatic screening of Coulomb interactions within the random phase approximation. We find that for massive Dirac systems, the oscillating screened potential (or density) decays as r −2 and r −3in two, and three dimensions respectively, and as r −1 for one dimensional non-interacting systems. However for massless Dirac systems there is no electrostatic screening or Friedel oscillation in one dimension, and the oscillating screened potential decays as r −3 and r −4 , in two and three dimensions respectively. Our analytical results for the polarization function will be useful for exploring the physics of massive and massless Dirac electrons in different experimental systems with varying dimensionality.

We explore the collective density oscillations of a collection of charged massive Dirac particles, in one, two and three dimensions and their one dimensional superlattice. We calculate the long wavelength limit of the dynamical polarization function analytically, and use the random phase approximation to obtain the plasmon dispersion. The density dependence of the long wavelength plasmon frequency in massive Dirac systems is found to be different as compared to systems with parabolic, and gapless Dirac dispersion. We also calculate the long wavelength plasmon dispersion of a 1d metamaterial made from 1d and 2d massive Dirac plasma. Our analytical results will be useful for exploring the use of massive Dirac materials as electrostatically tunable plasmonic metamaterials and can be experimentally verified by infrared spectroscopy as in the case of graphene [Nat. Nanotechnol. 6, 630 (2011)].

We analyze the dynamics of a Luttinger model following a quench in the electron-electron interaction strength, where the change in the interaction strength occurs over a finite time scale τ . We study the Loschmidt echo (the overlap between the initial and final state) as a function of time, both numerically and within a perturbation scheme, treating the change in the interaction strength as a small parameter, for all τ . We derive the corrections appearing in (a) the Loschmidt echo for a finite quench duration τ , (b) the scaling of the echo following a sudden (τ → 0) quench, and (c) the scaling of the echo after an adiabatic (τ → ∞) quench. We study in detail the limiting cases of the echo in the early-time and infinite-time limits, and provide scaling arguments to understand these in a general context. We also show that our perturbative results are in good agreement with the exact numerical ones.

Extended Bose Hubbard models with nearest neighbour interaction describe minimally the effect of long range interaction on ultra cold atoms in deep optical lattices. Rotation of such optical lattices subject such neutral cold atoms to the effect of an artificial magnetic field. The modification of the phase boundaries of the density wave and Mott Insulator phases due to this rotation are shown to be related to the edge spectrum of spinorial and scalar Harper equation. Corresponding profiles of the checkerboard vortex states with sublattice modulated superfluid order parameter near density wave phase boundary are calculated.