We present the first fully relativistic numerical calculations of differentially rotating Strange Quark Stars models for broad ranges of the maximum density and of the degree of differential rotation. Our simulations are performed with the very accurate and stable multi-domain spectral code FlatStar and use the MIT Bag model for describing strange quark matter. Our calculations based on a thorough exploration of the solution space show that the maximum mass of strange stars depends on both the degree of differential rotation and a type of solution, similarly to neutron stars described by a polytropic equations of state. The highest increase of the maximum mass (compared to the value for a non-rotating star) is obtained for models with a low degree of differential rotation. This highest mass is over four times larger than that of the equivalent non rotating configuration. Comparing our results with calculations done for realistic models of neutron stars, we conclude that with the help of differential rotation, strange stars can sustain masses much larger than stars made from nuclear matter for low degree of differential rotation which reinforces the hope of demonstrating, or of ruling out, the existence of strange matter through the observation by gravitational waves, by gamma rays or by neutrinos of the massive material object born from the merger of a compact binary system or during some supernova events.

The dense environment of globular clusters (GCs) can facilitate the formation of binary black holes (BBHs), some of which can merge with gravitational waves (GW) within the age of the Universe. We have performed a survey of Monte-Carlo simulations following the dynamical evolution of GCs with different masses, sizes and binary fractions and explored the impact of the host GC properties on the formation of BBH mergers. We find that the number of BBH mergers from GCs is determined by the GC's initial mass, size and primordial binary fraction. We identify two groups of BBH mergers: a primordial group whose formation does not depend on cluster's dynamics and a dynamical group whose formation is driven by the cluster's dynamical evolution. We show how the BBH origin affects the BBH mergers' main properties such as the chirp mass and merging time distributions. We provide analytic expressions for the dependence of the number of BBH mergers from individual GCs on the main cluster's structural properties and the time evolution of the merger rates of these BBHs. These expressions provide an essential ingredient for a general framework allowing to estimate the merger rate density. Using the relations found in our study, we find a local merger rate density of 0.18$-$1.8 ${\rm Gpc}^{-3}{\rm yr}^{-1}$ for primordial BBH mergers and 0.6$-$18 ${\rm Gpc}^{-3}{\rm yr}^{-1}$ for dynamical BBH mergers, depending on the GC mass and size distributions, initial binary fraction and the number density of GCs in the Universe.

In this first of a series of papers, we utilize results for around two thousand star cluster models simulated using the MOCCA code for star cluster evolution (Survey Database I) to determine the astrophysical properties and local merger rate densities for coalescing binary black holes (BBHs) originating from globular clusters (GCs). We extracted information for all coalescing BBHs that escape the cluster models and subsequently merge within a Hubble time along with BBHs that are retained in our GC models and merge inside the cluster via gravitational wave (GW) emission. By obtaining results from a substantial number of realistic star cluster models that cover different initial parameters, we have an extremely large statistical sample of BBHs with stellar mass and massive stellar BH ($\lesssim 100M_{\odot}$) components that merge within a Hubble time. Using this data, we estimate local merger rate densities for these BBHs originating from GCs to be at least 5.4 ${\rm Gpc}^{-3}\,{\rm yr}^{-1}$