1. A. Woszczyna, W. Czaja, K. Głód, Z. A. Golda, R. A. Kycia, A. Odrzywołek, P. Plaszczyk, L. M. Sokołowski, S. J. Szybka
ccgrg: The symbolic tensor analysis package with tools for general relativity
Wolfram Library Archive, vol. 8848 (2014).
[abstract] [journal]

Abstract:
Riemann and Weyl curvature, covariant derivative, Lie derivative, the first and the second fundamental form on hyper-surfaces, as well as basic notions of relativistic hydrodynamics (expansion, vorticity, shear) are predefined functions of the package. New tensors are easy to define. Instructions, basic examples, and some more advanced examples are attached to the package. Characteristic feature of the ccgrg package is the specific coupling between the functional programming and the Parker-Christensen index convention. This causes that no particular tools to rising/lowering tensor indices neither to the tensor contractions are needed. Tensor formulas are written in the form close to that of classical textbooks in GRG, with the only difference that the summation symbol appears explicitly. Tensors are functions, not matrixes, and their components are evaluated lazily. This means that only these components which are indispensable to realize the final task are computed. The memoization technique prevents repetitive evaluation of the same quantities. This saves both, time and memory.

2. Wojciech Czaja, Zdzislaw A. Golda, Andrzej Woszczyna
Consistent initial data for CMD perturbations
(2011).
[abstract] [preprint]

Abstract:
We investigate the initial condition which are simultaneously consistent with perturbation equations in both the radiation-dominated and the mater-dominated epochs. The exact formula for the spectrum transfer-function is derived.

3. Wojciech Czaja, Zdzislaw A. Golda, Andrzej Woszczyna
Momentum of cosmological acoustic field
(2011).
[abstract] [preprint]

Abstract:
The acoustic spacetime corresponding to perturbed Friedman-Lemaitre-Robertson-Walker universe inherit the space isometries from the original FLRW model, but essentially differs in dynamics. The scale factor manifestly depends on the equation of state of the matter content. Despite the higher complexity of the background evolution the perturbation equation in this space is substantially simpler: the density perturbations obey d'Alembert equation. Canonical formalism reconstructed in the acoustic spacetime enables one to employ the Klein-Gordon scalar product. Consequently, the Fourier decomposition of the perturbation field provide the time-independent Fourier coefficients and the time-independent spectrum. The perturbation spectrum does not depend of the choice of the Cauchy hypersurface from which the data are collected. Noether constants associated with the six-parameter isometry group define the components of the momentum, hyperbolic momentum and angular momentum of sound.

4. Wojciech Czaja, Zdzislaw A. Golda, Andrzej Woszczyna
The acoustic wave equation in the expanding universe. Sachs-Wolfe theorem
The Mathematica Journal, vol. 13 (2011).
[abstract] [journal] [cdf]

Abstract:
This article considers the acoustic field propagating in the radiation-dominated universe of arbitrary space curvature. The field equations are reduced to the d’Alembert equation in an auxiliary static Robertson-Walker spacetime and dispersion relations are discussed. *supported by the grant from The John Templeton Foundation

5. Marek Szydlowski, Adam Krawiec, Wojciech Czaja
Phantom cosmology as a simple model with dynamical complexity
Phys. Rev., vol. E72, p. 036221 (2005).
[abstract] [preprint] [journal]

Abstract:
We study the Friedman-Robertson-Walker model with phantom fields modeled in terms of scalar fields. We apply the Ziglin theory of integrability and find that the flat model is nonintegrable. Then we cannot expect to determine simple analytical solutions of the Einstein equations. We demonstrate that there is only a discrete set of parameters where this model is integrable. For comparison we describe the phantoms fields in terms of the barotropic equation of state. It is shown that in contrast to the phantoms modeled as scalar fields, the dynamics is always integrable and phase portraits are contracted. In this case we find the duality relation.

6. Marek Szydlowski, Wojciech Czaja
Modified Friedmann cosmologies -Theory and observations
Annals Phys. , vol. 320, pp. 261-281 (2005).
[abstract] [preprint] [journal]

Abstract:
This paper shows the results of the investigation of a class of Cardassian scenarios of the evolution of the universe in the formulation of the qualitative theory of dynamical systems. That theory allowed us to analyze all solutions for all possible initial conditions on the phase plane. In the Cardassian models we find that big-rip singularities are present as a typical behavior in the future if $n < 0$. We were also able to find some exact solution for the flat Cardassian models as well as a duality relation. In turn for the statistical analysis of SNIa data, without any priors on the matter content in the model, we obtain that the big-rip scenario is favored. The potential function for the Hamiltonian description of dynamics was reconstructed from the SNIa data (inverse dynamical problem).