The accuracy of the numerical findings of Lukes-Gerakopoulos PRD 86, 044013 (2012), regarding the existence of additional first integrals in the Zipoy-Voorhees space-time, was recently questioned by Maciejewski et al. PRD 88, 064003 (2013). In this comment, it is shown that the discrepancy between the results of Lukes-Gerakopoulos (2012) and Maciejewski et al. (2013) is not due to issues related to numerical accuracy, as claimed in Maciejewski et al (2013), but due to different choice of coordinates used in Maciejewski et al.(2013).

We study the possibility of the vortices formation in axion condensate on the galactic scale. Such vortices can occur as a result of global rotation of the early universe. We study analytical models of vortices and calculate exemplary galaxy rotation curves. Depending on the setup it is possible to obtain a variety of shapes which give a good qualitative agreement with observational results. However, as we show, the extremely low velocity dispersions of the axions velocity are required to form the single vortex on the galactic scales. We find that the required velocity dispersion is of the order of $\sigma \approx 10^{-12} \frac{m}{s}$. This is much smaller that predicted within the present understanding of the axion physics. The vortices in axion condensate can be however formed on the much smaller scales and give seeds to the galaxy formation and to their angular momenta. On the other hand, the vortices can be formed on the galactic scales, but only if the mass of the axion-like particles is of the order of $10^{-30}$eV. In this case, the particle de Broglie wavelength is comparable with the galactic diameter. This condition must be fulfilled in order to keep the coherence of the quantum condensate on the galactic scales.

We discuss a more general class of phantom ($p < -\varrho$) cosmologies with various forms of both phantom ($w < -1$), and standard ($w > -1$) matter. We show that many types of evolution which include both Big-Bang and Big-Rip singularities are admitted and give explicit examples. Among some interesting models, there exist non-singular oscillating (or "bounce") cosmologies, which appear due to a competition between positive and negative pressure of variety of matter content. From the point of view of the current observations the most interesting cosmologies are the ones which start with a Big-Bang and terminate at a Big-Rip. A related consequence of having a possibility of two types of singularities is that there exists an unstable static universe approached by the two asymptotic models - one of them reaches Big-Bang, and another reaches Big-Rip. We also give explicit relations between density parameters $\Omega$ and the dynamical characteristics for these generalized phantom models, including higher-order observational characteristics such as jerk and "kerk". Finally, we discuss the observational quantities such as luminosity distance, angular diameter, and source counts, both in series expansion and explicitly, for phantom models. Our series expansion formulas for the luminosity distance and the apparent magnitude go as far as to the fourth-order in redshift $z$ term, which includes explicitly not only the jerk, but also the "kerk" (or "snap") which may serve as an indicator of the curvature of the universe.

We investigate the possibility of testing cosmological models with interaction between matter and energy sector. We assume the standard FRW model while the so called energy conservation condition is interpreted locally in terms of energy transfer. We analyze two forms of dark energy sectors: the cosmological constant and phantom field. We find a simple exact solution of the models in which energy transfer is described by a Cardassian like term in the relation of $H^{2}(z)$, where $H$ is Hubble's function and $z$ is redshift. The considered models have two additional parameters $(\Omega_{\text{int}},n)$ (apart the parameters of the $\Lambda$CDM model) which can be tested using SNIa data. In the estimation of the model parameters Riess et al.'s sample is used. We also confront the quality of statistical fits for both the $\Lambda$CDM model and the interacting models with the help of the Akaike and Bayesian informative criteria. Our conclusion from standard best fit method is that the interacting models explains the acceleration of the Universe better but they give rise to a universe with high matter density. However, using the tools of information criteria we find that the two new parameters play an insufficient role in improving the fit to SNIa data and the standard $\Lambda$CDM model is still preferred. We conclude that high precision detection of high redshift supernovae could supply data capable of justifying adoption of new parameters.

We discuss a class of phantom ($p < - \varrho$) cosmological models. Except for phantom we admit various forms of standard types of matter and discuss the problem of singularities for these cosmologies. The singularities are different from those of standard matter cosmology since they appear for infinite values of the scale factor. We also find an interesting relation between the phantom models and standard matter models which is like the duality symmetry of string cosmology.