The parallelization, design and scalability of the \sky code to search for periodic gravitational waves from rotating neutron stars is discussed. The code is based on an efficient implementation of the F-statistic using the Fast Fourier Transform algorithm. To perform an analysis of data from the advanced LIGO and Virgo gravitational wave detectors' network, which will start operating in 2015, hundreds of millions of CPU hours will be required - the code utilizing the potential of massively parallel supercomputers is therefore mandatory. We have parallelized the code using the Message Passing Interface standard, implemented a mechanism for combining the searches at different sky-positions and frequency bands into one extremely scalable program. The parallel I/O interface is used to escape bottlenecks, when writing the generated data into file system. This allowed to develop a highly scalable computation code, which would enable the data analysis at large scales on acceptable time scales. Benchmarking of the code on a Cray XE6 system was performed to show efficiency of our parallelization concept and to demonstrate scaling up to 50 thousand cores in parallel.

We discuss an idea for how accreting millisecond pulsars could contribute to the understanding of the QCD phase transition in the high-density nuclear matter equation of state (EoS). It is based on two ingredients, the first one being a ``phase diagram'' of rapidly rotating compact star configurations in the plane of spin frequency and mass, determined with state-of-the-art hybrid equations of state, allowing for a transition to color superconducting quark matter. The second is the study of spin-up and accretion evolution in this phase diagram. We show that the quark matter phase transition leads to a characteristic line in the Omega-M plane, the phase border between neutron stars and hybrid stars with a quark matter core. Along this line a change in the pulsar's moment of inertia entails a waiting point phenomenon in the accreting millisecond X-ray pulsar (AMXP) evolution: most of these objects should therefore be found along the phase border in the Omega-M plane, which may be viewed as the AMXP analog of the main sequence in the Hertzsprung-Russell diagram for normal stars. In order to prove the existence of a high-density phase transition in the cores of compact stars we need population statistics for AMXP's with sufficiently accurate determination of their masses and spin frequencies.

We consider the influence of potential quark matter existing at high densities in neutron star interiors on gravitational waves (GW) emitted in a binary neutron star merger event. Two types of equations of state (EoS) at zero temperatures are used, one describing pure nuclear matter, the other nuclear matter with a phase transition to quark matter at very high densities. Binary equilibrium sequences close to the innermost stable circular orbit (ISCO) are calculated to determine the GW frequencies just before merger. It is found that EoS effects begin to play a role for gravitational masses larger than $M_\infty\simeq1.5M_\odot$. The difference in the gravitational wave frequency at the ISCO grows to up to $\simeq 10%$ for the maximal allowed mass given by the EoSs used. Then, we perform 3D hydrodynamic simulations for each EoS varying the initial mass and determine the characteristic GW frequencies of the merger remnants. The implications of quark matter show up mainly in a different collapse behaviour of the merger remnant. If the collapse does not take place immediately after merger, we find a phase difference between two EoS's in the post-merger GW signal. We also compare the GW frequencies emitted by the merger remnant to values from simulations using a polytropic EoS and find an imprint of the non-constant adiabatic index of our EoSs. All calculations are based on the conformally flat (CF) approximation to general relativity and the GW signal from the merger simulation is extracted up to quadrupole order.

Binary neutron star mergers are expected to be one of the most promising source of gravitational waves (GW) for the network of laser interferometric and bar detectors becoming operational in the next few years. The merger wave signal is expected to be sensitive to the interior structure of the neutron star (NS). The structure of high density phases of matter is under current experimental investigation in heavy-ion collisions. We investigate the dependence of the merger process and its GW signal on the presence of quarks in these phases by performing numerical simulations, where the smoothed particle hydrodynamics (SPH) method and the conformally flat approximation for the 3-geometry in general relativistic gravity are implemented.

For the classification of rotating compact stars with two high density phases a phase diagram in the angular velocity (Omega) - baryon number (N) plane is investigated. The dividing line N_crit(Omega) between configurations with one and two phases is correlated to a local maximum of the moment of inertia and can thus be subject to experimental verification by observation of the rotational behavior of accreting compact stars. Another characteristic line, which also can be measured is the transition line to black holes that of the maximum mass configurations. The positions and the shape of these lines are sensitive to changes in the equation of state (EoS) of stellar matter. A comparison of the regional structure of phase diagrams is performed for polytropic and relativistic mean field type EoS and correlations between the topology of the transition lines and the properties of two-phase EoS are obtained. Different scenarios of compact star evolution are discussed as trajectories in the phase diagram. It is shown that a population gap in the Omega - N plane for accreting compact stars would signal a high density phase transition and could be observed in the distribution of so called Z sources of quasi periodic oscillations in low-mass X-ray binaries.

Neutron star configurations are considered as thermodynamical systems for which a phase diagram in the angular velocity (Omega) - baryon number (N) plane is obtained with a dividing line N_{crit}(Omega) for quark core configurations. Trajectories of neutron star evolution in this diagram are studied for different scenarios defined by the external torque acting on the star due to radiation and/or mass accretion. They show a characteristic change in the rotational kinematics when the star enters the quark core regime. For isolated pulsars the braking index signal for deconfinement has been studied in its dependence on the mass of the star. Model calculations of the spin evolution of accreting low-mass X-ray binaries in the phase diagram have been performed for different values of the initial magnetic field, its decay time as well as initial mass and mass accretion rate. Population clustering of these objects at the line N_{crit}(Omega) in the phase diagram is suggested as an observable signal for the deconfinement phase transition if it exists for spinnning and accreting neutron stars.

It has been suggested that the QPO phenomenon in LMXB's could be explained when the central compact object is a strange star. In this work we investigate within a standard model for disk accretion whether the observed clustering of spin frequencies in a narrow band is in accordance with this hypothesis. We show that frequency clustering occurs for accreting strange stars when typical values of the parameters of magnetic field initial strength and decay time, accretion rate are chosen. In contrast to hybrid star accretion no mass clustering effect is found.

We study the evolution of the rotation frequency for accreting compact stars. The discontinuous change of the moment of inertia of a rapidly rotating star due to the possible quark core appearance entails a characteristic change in the spin evolution. Numerical solutions have been performed using a model equation of state describing the deconfinement phase transition. Trajectories of spin evolution are discussed in the angular velocity - baryon number plane (phase diagram) for different accretion scenarios defined by the initial values of mass and magnetic field of the star, as well as mass accretion rate and magnetic field decay time. We observe a characteristic increase in the waiting time when a configuration enters the quark core regime. Overclustering of the population of Z sources of LMXBs in the phase diagram is suggested as a direct measurement of the critical line for the deconfinement phase transition since it is related to the behaviour of the moment of inertia of the compact star.

Using the formalism of general relativity for axially symmetric gravitational fields and their sources - rotating compact stars - a perturbation theory with respect to angular velocity is developed and physical quantities such as mass, shape, momentum of inertia and total energy of the star are defined. The change of the internal structure of the star due to rotation has been investigated and the different contributions to the moment of inertia have been evaluated separately. Numerical solutions have been performed using a two-flavor model equation of state describing the deconfinement phase transition as constrained by the conservation of total baryon number and electric charge. During the spin down evolution of the rotating neutron star, below critical values of angular velocity a quark matter core can appear which might be detected as a characteristic signal in the pulsar timing. Within the spin-down scenario due to magnetic dipole radiation it is shown that the deviation of the braking index from $n=3$ could signal not only the occurrence but also the size of a quark core in the pulsar. A new scenario is proposed where, due to mass accretion onto the rapidly rotating compact star, a spin-down to spin-up transition might signal a deconfinement transition in its interior.

The occurence of a quark matter core in rotating compact stars has been investigated within general relativity as a function of both the rotational frequency and the total baryon number. We demonstrate that the deviation of the braking index from $n=3$ signals not only the occurence but also the size of a quark matter core in a pulsar. We suggest that in systems with mass accretion onto a rapidly rotating compact star a spin-down to spin-up transition might signal a deconfinement transition in its interior.

We explore the consequences of an equation of state (EOS) obtained in a confining Dyson-Schwinger equation model of QCD for the structure and stability of nonstrange quark stars at finite-T, and compare the results with those obtained using a bag-model EOS. Both models support a temperature profile that varies over the star's volume and the consequences of this are model independent. However, in our model the analogue of the bag pressure is (T,mu)-dependent, which is not the case in the bag model. This is a significant qualitative difference and comparing the results effects a primary goal of elucidating the sensitivity of quark star properties to the form of the EOS.