Thermodynamics of the Polyakov-Nambu-Jona-Lasinio model with nonzero baryon and isospin chemical potentials

We have extended the Polyakov-Nambu-Jona-Lasinio (PNJL) model for two degenerate flavors to include the isospin chemical potential (mu(I)). All the diagonal and mixed derivatives of pressure with respect to the quark number (proportional to baryon number) chemical potential (mu(0)) and isospin chemical potential up to sixth order have been extracted at mu(0)=mu(I)=0. These derivatives give the generalized susceptibilities with respect to quark and isospin numbers. Similar estimates for the flavor diagonal and off-diagonal susceptibilities are also presented. Comparison to lattice QCD (LQCD) data of some of these susceptibilities for which LQCD data are available show similar temperature dependence, though there are some quantitative deviations above the crossover temperature. We have also looked at the effects of instanton induced flavor mixing coming from the U-A(1) chiral symmetry breaking 't Hooft determinantlike term in the NJL part of the model. The diagonal quark number and isospin susceptibilities are completely unaffected. The off-diagonal susceptibilities show significant dependence near the crossover. Finally we present the chemical potential dependence of specific heat and speed of sound within the limits of chemical potentials where neither diquarks nor pions can condense.