I present a set of long-term, direct, relativistic many-body computations of model dense stellar clusters with up-to-date stellar-evolutionary, supernova (SN), and remnant natal-kick models, including pair instability and pulsation pair instability supernova (PSN and PPSN), using an updated version of NBODY7 N-body simulation program. The N-body model also includes stellar evolution-based natal spins of BHs and treatments of binary black hole (BBH) mergers based on numerical relativity. These, for the first time in a direct N-body simulation, allow for second-generation BBH mergers. The set of 65 evolutionary models have initial masses $10^4M_\odot-10^5M_\odot$, sizes 1 pc-3 pc, metallicity $0.0001-0.02$, with the massive stars in primordial binaries and they represent young massive clusters (YMC) and moderately massive open clusters (OC). Such models produce dynamically-paired BBH mergers that agree well with the observed masses, mass ratios, effective spin parameters, and final spins of the LVC O1/O2 merger events, provided BHs are born with low or no spin but spin up after undergoing a BBH merger or matter accretion onto it. In particular, the distinctly higher mass, effective spin parameter, and final spin of GW170729 merger event is naturally reproduced, as also the mass asymmetry of the O3 event GW190412. The computed models produce massive, $\sim100M_\odot$ BBH mergers with primary mass within the 'PSN gap' and also yield mergers involving remnants in the 'mass gap'. They also suggest that YMCs and OCs produce persistent, Local-Universe GW sources detectable by LISA. Such clusters are also capable of producing eccentric LIGO-Virgo mergers.

I study the potential role of young massive (YMCs) and open star clusters (OCs) in assembling stellar-mass binary black holes (BBHs) which would be detectable as persistent gravitational-wave (GW) sources by the forthcoming LISA mission. The energetic dynamical interactions inside star clusters make them factories of assembling BBHs and other types of double-compact binaries that undergo general-relativistic (GR) inspiral and merger. The initial phase of such inspirals would, typically, sweep through the LISA GW band. Here, such LISA sources are studied from a set of evolutionary models of star clusters with masses ranging over $10^4M_\odot-10^5M_\odot$ that represent YMCs and intermediate-aged OCs in metal-rich and metal-poor environments of the Local Universe. These models are evolved with long-term, direct, relativistic many-body computations incorporating state-of-the-art stellar-evolutionary and remnant-formation models. Based on models of Local Universe constructed with such model clusters, it is shown that YMCs and intermediate-aged OCs would yield several 10s to 100s of LISA BBH sources at the current cosmic epoch with GW frequency within $10^{-3}{\rm~Hz} - 10^{-1}{\rm~Hz}$ and signal-to-noise-ratio (S/N) $>5$, assuming a mission lifetime of 5 or 10 years. Such LISA BBHs would have a bimodal distribution in total mass, be generally eccentric ($\lesssim0.7$), and typically have similar component masses although mass-asymmetric systems are possible. Intrinsically, there would be 1000s of present-day, LISA-detectable BBHs from YMCs and OCs. That way, YMCs and OCs would provide a significant and the dominant contribution to the stellar-mass BBH population detectable by LISA. A small fraction, $<5$%, of these BBHs would undergo GR inspiral to make it to LIGO-Virgo GW frequency band and merge, within the mission timespan; $<15$% would do so within twice the timespan.

Owing to their high central densities, star clusters appear to be the ideal environment for the assembly of neutron star-neutron star (NS-NS) and black hole-neutron star (BH-NS) binaries. These binaries are among the most interesting astrophysical objects, being potential sources of several high-energy phenomena, such as gravitational waves (GWs) and gamma-ray bursts. In this paper, we use for the first time high-precision N-body simulations of young massive and open clusters to study the origin and dynamical evolution of NSs, within clusters with different initial masses, metallicities, primordial binary fractions, and prescriptions for the natal kicks imparted to NSs and BHs at birth. We find that the radial profile of NSs is shaped by the BH content of the cluster, which partially quenches the NS segregation due to the BH-burning process. This leaves most of the NSs out of the densest cluster regions, where NS-NS and BH-NS binaries could potentially form. Due to a large velocity kick that they receive at birth, most of the NSs escape the host clusters, with the bulk of their retained population made up of NSs of $\sim 1.3$ M$_\odot$ coming from the electron-capture supernova process. The details of the primordial binary fraction and pairing can smear out this trend. Finally, we find that a subset of our models produce NS-NS mergers, leading to a rate of $\sim 0.01$ Gpc$^{-3}$ yr$^{-1}$ in the local Universe, and compute an upper limit of $\sim 3\times 10^{-3}$ Gpc$^{-3}$ yr$^{-1}$ for the BH-NS merger rate. Our estimates are several orders of magnitude smaller than the current empirical merger rate from LIGO/Virgo, in agreement with the recent rate estimates for old globular clusters.

The limiting mass is a significant characteristic for compact exotic stars including quark stars and can be expressed in terms of fundamental constants and the Bag constant. In the present work, using bag model description, maximum mass of a rotating quark star is found to depend on the rotational frequency apart from other fundamental parameters. The analytical results obtained agrees with the result of several relevant numerical estimates as well as observational evidences.

The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster.

It is speculated that a merger of two massive astrophysical black holes in dense stellar environment may lead to the formation of a massive black hole in the pair-instability mass gap (70/80-125 Msun). Such a merger-formed black hole is expected to typically have a high spin (a=0.7), however, with a broad spectrum of spins being allowed (0<a<1). If such massive black hole acquires another black hole it may lead to another merger detectable by LIGO/Virgo in gravitational waves. We show that it is highly unlikely to form in this way and retain a 100 Msun black hole in a globular cluster if the black hole's spin is low (a<0.3). Massive merger-formed black holes with low spins acquire high recoil speeds (> 200 km/s) from gravitational-wave kick during formation that exceed typical escape speeds from globular clusters (50 km/s). However, a very low-spinning (a=0.1) and massive (100 Msun) black hole could be formed this way and retained in a galactic nuclear star cluster. Even though massive, merger-formed black holes with such low spins acquire high speeds during formation (400 km/s), they may avoid ejection since massive nuclear clusters have high escape velocities (300-500 km/s). This adds to already existing astrophysical scenarios of the formation of massive black holes with low spins. A future detection of a massive black hole in the pair-instability mass gap with low spin is therefore not a proof of the existence of primordial black holes, which are sometimes claimed to have low spins and arbitrarily high masses (not affected by pair-instability).

The Lambda-Cold Dark Matter (LCDM) model agrees with most of the cosmological observations, but has some hindrances from observed data at smaller scales such as galaxies. Recently, Berezhiani and Khoury (2015) proposed a new theory involving interacting superfluid dark matter with three model parameters, which explains galactic dynamics with great accuracy. In the present work, we study the cosmological behaviour of this model in the linear regime of cosmological perturbations. In particular, we compute both analytically and numerically the matter linear growth factor and obtain new bounds for the model parameters which are significantly stronger than previously found. These new constraints come from the fact that structures within the superfluid dark matter framework grow quicker than in LCDM, and quite rapidly when the DM-baryon interactions are strong.

The recent inference of a 70M_sun black hole (BH) in the Galactic, detached binary LB-1 has sparked cross-disciplinary debate since a stellar remnant of such large mass is well above what can be expected from stellar-evolutionary theory, especially in an enriched environment like that of the Milky Way. This study focusses on the possibilities of formation of extraordinarily massive BHs at solar and globular cluster (GC)-like metallicities via evolution of massive stellar binaries. A population-synthesis approach is followed utilizing the recently updated BSE program. BHs in the mass range of 50M_sun-80M_sun could be formed at the solar metallicity only if a large fraction, >70%, of matter is allowed to accrete onto a low-mass BH, in a BH-star merger product (a "black hole Thorne-Zytkow object"; BH-TZO). Their counterparts at GC-like metallicities can reach 100 M_sun. Although post-accretion BHs can, generally, be expected to be of high spin parameter, they can potentially be of low spin in the case of a BH-TZO. This spin aspect remains speculative in this work and deserves detailed hydrodynamic studies.

This has been suggested that the dwarf satellite galaxies in the Milky Way may contain substantial amount of dark matter in them. These dark matters may undergo self-annihilation to produce $\gamma$-rays. The satellite borne $\gamma$-rays telescope such as Fermi-LAT reported the detection of $\gamma$-rays from around 45 Dwarf Spheroidals (dSphs) of Milky Way. In this work, we consider a particle dark matter model and after studying its phenomenology, we calculate the $\gamma$-ray fluxes of each of these 45 dSphs and compare our results with the upper bounds of mass vs annihilation cross-sections of dark matter provided by the Fermi-LAT collaboration. We calculate these fluxes by considering different dark matter density profiles and make a comparison of the results. We also repeat our analysis with another dark matter candidate namely Kaluza-Klein dark matter inspired by extra-dimensional models. We make a critical comparison between the results obtained for each of these models vis-a-vis the observational results.

We show that a nonsingular bounce, free of ghosts and gradient instabilities, can be realized in the framework of Horndeski or generalized Galileon cosmology. In particular, we first review that the theoretical no-go theorem, which states that the above is impossible, is based on two very strong assumptions, namely that a particular quantity cannot be discontinuous during the bounce, and that there is only one bounce. However, as we show in the present work, the first assumption not only can be violated in a general Horndeski/Galileon scenario, but also it is necessarily violated at the bounce point within the subclass of Horndeski/Galileon gravity in which $K(\phi,X)$ becomes zero at $X=0$. Additionally, concerning the second assumption, which is crucial in improved versions of the theorem which claim that even if a nonlinear free of pathologies can be realized it will lead to pathologies in the infinite past or infinite future, we show that if needed it can be evaded by considering cyclic cosmology, with an infinite sequence of nonsingular bounces free of pathologies, which forbids the universe to reach the "problematic" regime at infinite past or infinite future. Finally, in order to make the analysis more transparent we provide explicit examples where nonsingular bounces without theoretical pathologies can be achieved.

A thin viscous accretion disc around a Kerr black hole, which is warped due to the Lense-Thirring (LT) effect, was shown to cause the spin axis of the black hole to precess and align with the outer disc. We calculate the total LT torque acting on the black hole, and compute the alignment and precession time-scales for both persistent and transient accretors. In our analysis, we consider the contribution of the inner disc, as it can stay misaligned with the black hole spin for a reasonable range of parameter values. We find that the alignment time-scale increases with a decrease in the Kerr parameter below a critical Kerr parameter value, contrary to earlier predictions. Besides, the time-scales are generally longer for transience than the time-scales calculated for persistent accretion. From our analysis of the transient case, we find that the black hole in the low mass X-ray binary (LMXB) 4U 1543-47 could be misaligned, whereas that in the LMXB XTE J1550-564 has aligned itself with the outer disc. The age of the LMXB H 1743-322 is estimated assuming a misaligned disc. We also find that the black hole in a typical Galactic LMXB can take a significantly longer time to align than what was estimated in the past. This may have an important implication on the measurement of black hole spin using the continuum X-ray spectral fitting method.

The inner part of a thin accretion disk around a Kerr black hole can serve as an important tool to study the physics of the strong gravity regime. A tilt in such a disk with respect to the black hole spin axis is particularly useful for this purpose, as such a tilt can have a significant effect on the observed X-ray spectral and timing features via Lense-Thirring precession. However, the inner disk has been predicted to become aligned with the spin direction of the black hole by the well-known Bardeen-Petterson effect. Here we calculate, both analytically and numerically, the radial profile of the thin accretion disk tilt angle in the viscous regime (i.e., $\alpha$ > H/R; $\alpha$ is the Shakura-Sunyaev viscosity parameter, H is the disk thickness and R is the radial distance). We show that the inner disk may not be aligned at all for certain reasonable ranges of parameter values. This makes the inner accretion disk particularly promising to probe the black hole parameters, and the accretion process in the strong gravity region.

Long-term solar magnetic activity reconstructions indicate the solar dynamo operates in two distinct -- grand minimum and regular activity -- modes. By employing bifurcation analysis of a reduced dynamo model we establish this to be a consequence of dynamo hysteresis. We reproduce the observed bimodal distribution of sunspots, but only for subcritical dynamos. A theoretical framework consistent with this finding explains observations of an abrupt midlife transition in stellar activity, characterized by reduced angular momentum loss rates and breakdown of gyrochronology relations. Our study indicates that an evolving dynamo bridges a diversity of phenomena in solar-like stars across their lifetime.

Stellar-remnant black holes (BH) in dense stellar clusters comprise a natural setup to trigger general-relativistic (GR) inspiral and merger of binary black holes (BBH), detectable by the LISA and the LIGO-Virgo, through dynamical encounters inside such environments. In this work, the intricacies of such dynamical interactions are probed utilizing realistic, self-consistent, post-Newtonian, direct N-body evolutionary models of young massive and open stellar clusters. Particularly, the configurations of the compact subsystems, that drive the in-cluster GR BBH coalescences, are tracked on the fly. Such an approach reveals that the GR coalescences within the open clusters take place primarily via chaotic interactions involving triple BH systems. Although less frequently, such mergers are found to happen also in higher-order subsystems such as quadruples and in subsystems involving non-BH members; the mergers can themselves be BH---non-BH, which events would leave electromagnetic signatures. Close, fly-by encounters inside the clusters can also make BBHs and other types of double-compact binaries temporarily post-Newtonian; such binaries would potentially contribute to the GW background for the LISA and the PTA. These calculations, furthermore, suggest that open clusters are potential hosts for not only detached BH---main-sequence binaries, as recently identified in the globular cluster NGC 3201, but also a wide variety of other types of remnant---non-remnant binaries, which are assembled via dynamical interactions inside the clusters and which have the prospects of being discovered in radial-velocity surveys.

Globular clusters (GCs) are the ideal environment for the formation of neutron stars (NSs) and millisecond pulsars (MSPs). NSs origin and evolution provide a useful information on stellar dynamics and evolution in star clusters, and are among the most interesting astrophysical objects, being precursors of several high-energy phenomena such as gravitational waves and gamma-ray bursts. Due to a large velocity kick that they receive at birth, most of the NSs escape the local field, affecting the evolution and dynamics of their parent cluster. In this paper, we study the origin and dynamical evolution of NSs within GCs with different initial masses, metallicities and primordial binary fractions. We find that the radial profile of NSs is shaped by the BH content of the cluster, which partially quenches the NS segregation until most of the BHs are ejected from the system. Independently on the cluster mass and initial configuration, the NSs map the average stellar population, as their average radial distance is $\approx 60-80\%$ of the cluster half-mass radius. Finally, by assuming a recycling fraction of $f_\mathrm{rec}=0.1$ and an average MSP gamma-ray emission of $L_\gamma=2\times 10^{33}$ erg s$^{-1}$, we show that the typical gamma-ray emission from our GCs agrees with observations and supports the MSP origin of the gamma-ray excess signal observed by the Fermi-LAT telescope in the Galactic Centre.

In this work we study certain classes of $F(R,{\cal G})$ gravity which have appealing phenomenological features, with respect to the successful realization of the dark energy and of the inflationary era. Particularly, we discuss the general formalism and we demonstrate how several inflationary and dark energy evolutions can be described in the context of $F(R,{\cal G})$ gravity. Also we propose a unified model, in the context of which the early and late-time dynamics are controlled by the $F(R,{\cal G})$ gravity, thus producing inflation and the dark energy era, while the intermediate era is approximately identical with standard Einstein-Hilbert gravity. Also we calculate the power spectrum of the primordial curvature perturbations corresponding to the unified $F(R,{\cal G})$ gravity model we propose, which as we demonstrate is nearly scale invariant and compatible with the latest observational data constraints.

We study the bounce and cyclicity realization in the framework of new gravitational scalar-tensor theories. In these theories the Lagrangian contains the Ricci scalar and its first and second derivatives, in a specific combination that makes them free of ghosts, and transformed into the Einstein frame they are proved to be a subclass of bi-scalar extensions of general relativity. We present analytical expressions for the bounce requirements, and we examine the necessary qualitative behavior of the involved functions that can give rise to a given scale factor. Having in mind these qualitative forms, we reverse the procedure and we construct suitable simple Lagrangian functions that can give rise to a bounce or cyclic scale factor.

We explore a reduced Babcock-Leighton (BL) dynamo model based on delay differential equations using numerical bifurcation analysis. This model reveals hysteresis, seen in the recent mean-field dynamo model and the direct numerical simulations of turbulent dynamos. The BL model with 'magnetic noise' as an additional weak-source of the poloidal field recovers the solar cycle every time from grand minima, which BL source alone cannot do. The noise-incorporated model exhibits a bimodal distribution of toroidal field energy confirming two modes of solar activity. It also shows intermittency and reproduces phase space collapse, an experimental signature of the Maunder Minimum. The occurrence statistics of grand minima in our model agree reasonably well with the observed statistics in the reconstructed sunspot number. Finally, we demonstrate that the level of magnetic noise controls the duration of grand minima and even has a handle over its waiting period, suggesting a triggering effect of grand minima by the noise and thus shutting down the global dynamo. Therefore, we conclude that the 'magnetic noise' due to small-scale turbulent dynamo action (or other sources) plays a vital role even in Babcock-Leighton dynamo models.

Three-dimensional, compressible, magnetohydrodynamic turbulence of an isothermal, self-gravitating fluid is analyzed using two-point statistics in the asymptotic limit of large Reynolds numbers (both kinetic and magnetic). Following an alternative formulation proposed by S. Banerjee and S. Galtier (Phys. Rev. E,93, 033120, 2016) and S. Banerjee and S. Galtier (J. Phys. A, Math. and Theor.,50, 015501, 2017), an exact relation has been derived for the total energy transfer. This approach results in a simpler relation expressed entirely in terms of mixed second-order structure functions. The kinetic, thermodynamic, magnetic and gravitational contributions to the energy transfer rate can be easily separated in the present form. By construction, the new formalism includes such additional effects as global rotation, the Hall term in the induction equation, etc. The analysis shows that solid-body rotation cannot alter the energy flux rate of compressible turbulence. However, the contribution of a uniform background magnetic field to the flux is shown to be non-trivial unlike in the incompressible case. Finally, the compressible, turbulent energy flux rate does not vanish completely due to simple alignments, which leads to a zero turbulent energy flux rate in the incompressible case.

Self-gravitating isothermal supersonic turbulence is analyzed in the asymptotic limit of large Reynolds numbers. Based on the inviscid invariance of total energy, an exact relation is derived for homogeneous, (not necessarily isotropic) turbulence. A modified definition for the two-point energy correlation functions is used to comply with the requirement of detailed energy equipartition in the acoustic limit. In contrast to the previous relations (Galtier and Banerjee, Phys. Rev. Lett., 107, 134501, 2011; Banerjee and Galtier, Phys. Rev. E, 87, 013019, 2013), the current exact relation shows that the pressure dilatation terms plays practically no role in the energy cascade. Both the flux and source terms are written in terms of two-point differences. Sources enter the relation in a form of mixed second-order structure functions. Unlike kinetic and thermodynamic potential energy, gravitational contribution is absent from the flux term. An estimate shows that for the isotropic case, the correlation between density and gravitational acceleration may play an important role in modifying the energy transfer in self-gravitating turbulence. The exact relation is also written in an alternative form in terms of two-point correlation functions, which is then used to describe scale-by-scale energy budget in spectral space.

Modified gravity theories can introduce modifications to the Poisson equation in the Newtonian limit. As a result, we expect to see interesting features of these modifications inside stellar objects. White dwarf stars are one of the most well studied stars in stellar astrophysics. We explore the effect of modified gravity theories inside white dwarfs. We derive the modified stellar structure equations and solve them to study the mass-radius relationships for various modified gravity theories. We also constrain the parameter space of these theories from observations.

The study of stellar-remnant black holes (BH) in dense stellar clusters is now in the spotlight, especially due to their intrinsic ability to form binary black holes (BBH) through dynamical encounters, that potentially coalesce via gravitational-wave (GW) radiation. In this work, which is a continuation of a recent study (Paper I), additional models of compact stellar clusters with initial masses $\lesssim10^5M_\odot$ and also those with small fractions of primordial binaries ($\lesssim10$%) are evolved for long term, applying the direct N-body approach, assuming state-of-the-art stellar-wind and remnant-formation prescriptions. That way, a substantially broader range of computed models than that in Paper I is achieved. As in Paper I, the general-relativistic BBH mergers continue to be mostly mediated by triples that are bound to the clusters rather than happen among the ejected BBHs. In fact, the number of such in situ BBH mergers, per cluster, tend to increase significantly with the introduction of a small population of primordial binaries. Despite the presence of massive primordial binaries, the merging BBHs, especially the in situ ones, are found to be exclusively dynamically assembled and hence would be spin-orbit misaligned. The BBHs typically traverse through both the LISA's and the LIGO's detection bands, being audible to both instruments. The "dynamical heating" of the BHs keeps the Electron-Capture-Supernova (ECS) neutron stars (NS) from effectively mass segregating and participating in exchange interactions; the dynamically-active BHs would also exchange into any NS binary within $\lesssim1$ Gyr. Such young massive and open clusters have the potential to contribute to the dynamical BBH merger detection rate to a similar extent as their more massive globular-cluster counterparts.

We have earlier proposed a fourth order gravity model as a possible explanation for late time cosmic acceleration, and for flattened galaxy rotation curves. The model has a free length parameter whose value depends on the scale of the system under study (e.g. the whole Universe, a galaxy, or a compact binary pulsar). In the present work, we investigate the constraints imposed on the free model parameter by Hulse-Taylor binary pulsar data: periastron advance; and emission of gravitational waves and consequent period decay. It is shown that the model is consistent with these observations, provided the length parameter is bounded from above.

We have earlier proposed a fourth order gravity model as a possible explanation for late time cosmic acceleration, and for flattened galaxy rotation curves. The model has a free length parameter whose value depends on the scale of the system under study (e.g. the whole Universe, a galaxy, or a compact binary pulsar). In the present work, we investigate the constraints imposed on the free model parameter by Hulse-Taylor binary pulsar data: periastron advance; and emission of gravitational waves and consequent period decay. It is shown that the model is consistent with these observations, provided the length parameter is bounded from above.

The observed classicality of primordial perturbations, despite their quantum origin during inflation, calls for a mechanism for quantum-to-classical transition of these initial fluctuations. As literature suggests a number of plausible mechanisms which try to address this issue, it is of importance to seek for concrete observational signatures of these several approaches in order to have a better understanding of the early universe dynamics. Among these several approaches, it is the spontaneous collapse dynamics of Quantum Mechanics which is most viable of leaving discrete observational signatures as collapse mechanism inherently changes the generic Quantum dynamics. We observe in this study that the observables from the scalar sector, i.e. scalar tilt $n_s$, running of scalar tilt $\alpha_s$ and running of running of scalar tilt $\beta_s$, can not potentially distinguish a collapse modified inflationary dynamics in the realm of canonical scalar field and $k-$inflationary scenarios. The only distinguishable imprint of collapse mechanism lies in the observables of tensor sector in the form of modified consistency relation and a blue-tilted tensor spectrum only when the collapse parameter $\delta$ is non-zero and positive.

We investigate the bounce and cyclicity realization in the framework of weakly broken galileon theories. We study bouncing and cyclic solutions at the background level, reconstructing the potential and the galileon functions that can give rise to a given scale factor, and presenting analytical expressions for the bounce requirements. We proceed to a detailed investigation of the perturbations, which after crossing the bouncing point give rise to various observables, such as the scalar and tensor spectral indices and the tensor-to-scalar ratio. Although the scenario at hand shares the disadvantage of all bouncing models, namely that it provides a large tensor-to-scalar ratio, introducing an additional light scalar significantly reduces it through the kinetic amplification of the isocurvature fluctuations.

The dynamical processes involving stellar-remnant black holes (BH) in stellar clusters has always drawn attention due to the BHs' potential in a number of astrophysical phenomena, especially the dynamical formation of binary black holes (BBH), which would potentially coalesce via radiation of gravitational waves (GW). This study presents a preliminary set of evolutionary models of compact stellar clusters with initial masses ranging over $1.0\times10^4M_\odot-5.0\times10^4M_\odot$, and half-mass radius of 2 pc or 1 pc, that is typical for young massive and starburst clusters. They have metallicities between $0.05Z_\odot-Z_\odot$. Including contemporary schemes for stellar wind and remnant-formation, such model clusters are evolved, for the first time, using the state-of-the-art direct N-body evolution program NBODY7, until their dissolution or at least for 10 Gyr. That way, a self-regulatory behaviour in the effects of dynamical interactions among the BHs, especially while heating and expanding the cluster and self-depleting the BHs, is demonstrated. The BBH mergers obtained here show a prominence in triple-mediated mergers while being bound to the clusters, compared to those occurring among the BBHs that are dynamically ejected from the clusters. This is in contrast with earlier N-body computations and also with recent Monte-Carlo method based ones. A broader mass spectrum of BHs and ejection of BBHs generally of wider orbits and in lower numbers, for the cluster masses explored here, might cause this which is yet to be fully understood. Among the BBH coalescences obtained here, there are ones that resemble the detected GW151226, LVT151012, and GW150914 events and also ones which are even more massive. A preliminary estimate suggests few 10s - 100s of BBH coalescences per year, originating due to dynamics in stellar clusters, that can be detected by the LIGO at its design sensitivity.

The residual gas within newly formed star clusters is expelled through stellar feedback on timescales ~ 1 Myr. The subsequent expansion of the cluster results in an unbinding of a fraction of stars before the remaining cluster members can re-virialize and form a surviving cluster. We investigate the bound fraction after gas expulsion as a function of initial cluster mass in stars and gauge the influence of primordial mass segregation, stellar evolution and the tidal field at the solar distance. We also assess the impact of the star-formation efficiency and gas expulsion velocity. We perform N-body simulations using Sverre Aarseth's NBODY7 code, starting with compact clusters in their embedded phase and approximate the gas expulsion by means of an exponentially depleting external gravitational field. We follow the process of re-virialization through detailed monitoring of different Lagrange radii over several Myr, examining initial half-mass radii of 0.1 pc, 0.3 pc and 0.5 pc and initial masses usually ranging from $5\times10^3 M_\odot$ to $5\times10^4 M_\odot$. The adopted star-formation efficiency of 0.33 in the cluster volume results in a distinct sensitivity to the gas expulsion velocity over a wide mass range, while a variation of the star-formation efficiency can make the cluster robust to the rapidly decreasing external potential. We confirm that primordial mass segregation leads to a smaller bound fraction, its influence possibly decreasing with mass. Stellar evolution has a higher impact on lower mass clusters, but heating through dynamical friction could expand the cluster to a similar extent. The examined clusters expand well within their tidal radii and would survive gas expulsion even in a strong tidal field.

Current data (LHC direct searches, Higgs mass, dark matter-related bounds) severely affect the constrained minimal SUSY standard model (CMSSM) with neutralinos as dark matter candidates. But the evidence for neutrino masses coming from oscillations requires extending the SM with at least right-handed neutrinos with a Dirac mass term. In turn, this implies extending the CMSSM with right-handed sneutrino superpartners, a scenario we dub $\tilde\nu$CMSSM. These additional states constitute alternative dark matter candidates of the superWIMP type, produced via the decay of the long-lived next-to-lightest SUSY particle (NLSP). Here we consider the interesting and likely case where the NLSP is a $\tilde{\tau}$: despite the modest extension with respect to the CMSSM this scenario has the distinctive signatures of heavy, stable charged particles. After taking into account the role played by neutrino mass bounds and the specific cosmological bounds from the big bang nucleosynthesis in selecting the viable parameter space, we discuss the excellent discovery prospects for this model at the future runs of the LHC. We show that it is possible to probe $\tilde{\tau}$ masses up to 600 GeV at the 14 TeV LHC with $\mathcal{L} = 1100$ fb$^{-1}$ when one considers a pair production of staus with two or more hard jets through all SUSY processes. We also show the complementary discovery prospects from a direct $\tilde{\tau}$ pair production, as well as at the new experiment MoEDAL.

Among the most striking feature of blue straggler stars (BSS) is the presence of multiple sequences of BSSs in the colour-magnitude diagrams (CMDs) of several globular clusters. It is often envisaged that such a multiple BSS sequence would arise due a recent core collapse of the host cluster, triggering a number of stellar collisions and binary mass transfers simultaneously over a brief episode of time. Here we examine this scenario using direct N-body computations of moderately massive star clusters (of order 10^4 Msun ). As a preliminary attempt, these models are initiated with approx. 8-10 Gyr old stellar population and King profiles of high concentrations, being "tuned" to undergo core collapse quickly. BSSs are indeed found to form in a "burst" at the onset of the core collapse and several of such BS-bursts occur during the post-core-collapse phase. In those models that include a few percent primordial binaries, both collisional and binary BSSs form after the onset of the (near) core-collapse. However, there is as such no clear discrimination between the two types of BSSs in the corresponding computed CMDs. We note that this may be due to the less number of BSSs formed in these less massive models than that in actual globular clusters.

Cosmological hysteresis is a purely thermodynamical phenomenon caused by the gradient in pressure, hence the characteristic equation of state during the expansion and contraction phases of the universe are different, provided that the universe bounces and recollapses. During hysteresis pressure asymmetry is created due to the presence of a single scalar field in the dynamical process. Also such an interesting scenario has vivid implications in cosmology when applied to variants of modified gravity models described within the framework of membrane paradigm. Cyclic universe along with scalar field leads to the increase in the amplitude of the cosmological scale factor at each consecutive cycles of the universe. Detailed analysis shows that the conditions which creates a universe with an ever increasing expansion, depend on the signature of the hysteresis loop integral $\oint pdV$ and on membrane model parameters.

Hysteresis is a phenomenon occurring naturally in several magnetic and electric materials in condensed matter physics. When applied to cosmology, aka cosmological hysteresis, has interesting and vivid implications in the scenario of a cyclic bouncy universe. Most importantly, this physical prescription can be treated as an alternative proposal to inflationary paradigm. Cosmological hysteresis is caused by the asymmetry in the equation of state parameter during expansion and contraction phase of the universe, due to the presence of a single scalar field. This process is purely thermodynamical in nature, results in a non-vanishing hysteresis loop integral $(\oint pdV)$ in cosmology. When applied to variants of modified gravity models -1) Dvali-Gabadadze-Porrati (DGP) brane world gravity, 2) Cosmological constant dominated Einstein gravity, 3) Loop Quantum Gravity (LQG), 4) Einstien-Gauss-Bonnet brane world gravity and 5) Randall Sundrum single brane world gravity (RSII), under certain circumstances, this phenomenon leads to the increase in amplitude of the consecutive cycles and to a universe with older and larger successive cycles, provided we have physical mechanisms to make the universe bounce and turnaround. This inculcates an arrow of time in a dissipationless cosmology. Remarkably, this phenomenon appears to be widespread in several cosmological potentials in variants of modified gravity background, which we explicitly study for- i) Hilltop, ii) Natural and iii) Colemann-Weinberg potentials, in this paper. Semi-analytical analysis of these models, for different potentials with minimum/minima, show that the conditions which creates a universe with an ever increasing expansion, depend on the signature of the hysteresis loop integral $(\oint pdV)$ as well as on the variants of model parameters.

Cosmological hysteresis, has interesting and vivid implications in the scenario of a cyclic bouncy universe. This, purely thermodynamical in nature, is caused by the asymmetry in the equation of state parameter during expansion and contraction phase of the universe, due to the presence of a single scalar field. When applied to variants of modified gravity models this phenomenon leads to the increase in amplitude of the consecutive cycles of the universe, provided we have physical mechanisms to make the universe bounce and turnaround. This also shows that the conditions which creates a universe with an ever increasing expansion, depend on the signature of $\oint pdV$ and on model parameters.

How Very Young Massive star Clusters (VYMCs; also known as "starburst" clusters), which typically are of $\gtrsim 10^4M_\odot$ and are a few Myr old, form out of Giant Molecular Clouds is still largely an open question. Increasingly detailed observations of young star clusters and star-forming molecular clouds and computational studies provide clues about their formation scenarios and the underlying physical processes involved. This chapter is focused on reviewing the decade-long studies that attempt to computationally reproduce the well-observed nearby VYMCs, such as the Orion Nebula Cluster, R136 and NGC 3603 young cluster, thereby shedding light on birth conditions of massive star clusters, in general. On this regard, focus is given on direct N-body modeling of real-sized massive star clusters, with a monolithic structure and undergoing residual gas expulsion, which have consistently reproduced the observed characteristics of several VYMCs and also of young star clusters, in general. The connection of these relatively simplified model calculations with the structural richness of dense molecular clouds and the complexity of hydrodynamic calculations of star cluster formation is presented in detail. Furthermore, the connections of such VYMCs with globular clusters, which are nearly as old as our Universe, is discussed. The chapter is concluded by addressing long-term deeply gas-embedded (at least apparently) and substructured systems like W3 Main. While most of the results are quoted from existing and up-to-date literature, in an integrated fashion, several new insights and discussions are provided.

How starburst clusters form out of molecular clouds is still an open question. In this article, I highlight some of the key constraints in this regard, that one can get from the dynamical evolutionary properties of dense stellar systems. I particularly focus on secular expansion of massive star clusters and hierarchical merging of sub-clusters, and discuss their implications vis-a-vis the observed properties of young massive clusters. The analysis suggests that residual gas expulsion is necessary for shaping these clusters as we see them today, irrespective of their monolithic or hierarchical mode of formation.

The classic question that how young massive star clusters attain their shapes and sizes, as we find them today, remains to be a challenge. Both observational and computational studies of star-forming massive molecular gas clouds infer that massive cluster formation is primarily triggered along the small-scale ($\lesssim0.3$ pc) filamentary substructures within the clouds. The present study is intended to investigate the possible ways in which a filament-like-compact, massive star cluster (effective radius 0.1-0.3 pc) can expand $\gtrsim10$ times, still remaining massive enough ($\gtrsim10^4 M_\odot$), to become a young massive star cluster, as we observe today. To that end, model massive clusters (of initially $10^4 M_\odot-10^5 M_\odot$) are evolved using Sverre Aarseth's state-of-the-art N-body code NBODY7. All the computed clusters expand with time, whose sizes (effective radii) are compared with those observed for young massive clusters, of age $\lesssim100$ Myr, in the Milky Way and other nearby galaxies. It is found that beginning from the above compact sizes, a star cluster cannot expand by its own, i.e., due to two-body relaxation, stellar-evolutionary mass loss, dynamical heating by primordial binaries and stellar-mass black holes, up to the observed sizes of young massive clusters; they always remain much more compact compared to the observed ones. This calls for additional mechanisms that can boost the expansion of a massive cluster after its assembly. Using further N-body calculations, it is shown that a substantial residual gas expulsion, with $\approx30$% star formation efficiency, can indeed swell the newborn embedded cluster adequately. The limitations of the present calculations and their consequences are discussed.

We perform a series of direct $N$-body calculations to investigate the effect of residual gas expulsion from the gas-embedded progenitors of present-day globular clusters (GCs) on the stellar mass function (MF). Our models start either tidally filling or underfilling, and either with or without primordial mass segregation. We cover 100 Myr of the evolution of modeled clusters and show that the expulsion of residual gas from initially mass-segregated clusters leads to a significantly shallower slope of the stellar MF in the low- ($m\leq 0.50 M_\odot$) and intermediate-mass ($\simeq 0.50-0.85 M_\odot$) regime. Therefore, the imprint of residual gas expulsion and primordial mass segregation might be visible in the present-day MF. We find that the strength of the external tidal field, as an essential parameter, influences the degree of flattening, such that a primordially mass-segregated tidally-filling cluster with $r_h/r_t\geq 0.1$ shows a strongly depleted MF in the intermediate stellar mass range. Therefore, the shape of the present-day stellar MF in this mass range probes the birth place of clusters in the Galactic environment. We furthermore find that this flattening agrees with the observed correlation between the concentration of a cluster and its MF slope, as found by de Marchi et al.. We show that if the expansion through the residual gas expulsion in primordial mass segregated clusters is the reason for this correlation then GCs most probably formed in strongly fluctuating local tidal fields in the early proto-Milky Way potential, supporting the recent conclusion by Marks \& Kroupa.

We investigate a fourth order model of gravity, having a free length parameter, and no cosmological constant or dark energy. We consider cosmological evolution of a flat Friedmann universe in this model for the case that the length parameter is of the order of present Hubble radius. By making a suitable choice for the present value of the Hubble parameter, and value of third derivative of the scale factor (the jerk) we find that the model can explain cosmic acceleration to the same degree of accuracy as the standard concordance model. If the free length parameter is assumed to be time-dependent, and of the order of the Hubble parameter of the corresponding epoch, the model can still explain cosmic acceleration, and provides a possible resolution of the cosmic coincidence problem. We work out the effective equation of state, and its time evolution, in our model. We also compare redshift drift in our model, with that in the standard model. The equation of state and the redshift drift serve to discriminate our model from the standard model.

The interstellar medium of galaxies is governed by supersonic turbulence, which likely controls the star formation rate (SFR) and the initial mass function (IMF). Interstellar turbulence is non-universal, with a wide range of Mach numbers, magnetic fields strengths, and driving mechanisms. Although some of these parameters were explored, most previous works assumed that the gas is isothermal. However, we know that cold molecular clouds form out of the warm atomic medium, with the gas passing through chemical and thermodynamic phases that are not isothermal. Here we determine the role of temperature variations by modelling non-isothermal turbulence with a polytropic equation of state (EOS), where pressure and temperature are functions of gas density, P~rho^Gamma, T~rho^(Gamma-1). We use grid resolutions of 2048^3 cells and compare polytropic exponents Gamma=0.7 (soft EOS), Gamma=1 (isothermal EOS), and Gamma=5/3 (stiff EOS). We find a complex network of non-isothermal filaments with more small-scale fragmentation occurring for Gamma<1, while Gamma>1 smoothes out density contrasts. The density probability distribution function (PDF) is significantly affected by temperature variations, with a power-law tail developing at low densities for Gamma>1. In contrast, the PDF becomes closer to a lognormal distribution for Gamma<=1. We derive and test a new density variance - Mach number relation that takes Gamma into account. This new relation is relevant for theoretical models of the SFR and IMF, because it determines the dense gas mass fraction of a cloud, from which stars form. We derive the SFR as a function of Gamma and find that it decreases by a factor of ~5 from Gamma=0.7 to Gamma=5/3.

The formation of very young massive clusters or "starburst" clusters is currently one of the most widely debated topic in astronomy. The classical notion dictates that a star cluster is formed in-situ in a dense molecular gas clump followed by a substantial residual gas expulsion. On the other hand, based on the observed morphologies of many young stellar associations, a hierarchical formation scenario is alternatively suggested. A very young (age $\approx$ 1 Myr), massive ($>10^4M_\odot$) star cluster like the Galactic NGC 3603 young cluster (HD 97950) is an appropriate testbed for distinguishing between such "monolithic" and "hierarchical" formation scenarios. A recent study by Banerjee and Kroupa (2014) demonstrates that the monolithic scenario remarkably reproduces the HD 97950 cluster. In the present work, we explore the possibility of the formation of the above cluster via hierarchical assembly of subclusters. These subclusters are initially distributed over a wide range of spatial volumes and have various modes of sub-clustering. Unlike the monolithic model that reproduces HD 97950 very well, the same is found to be prohibitive with hierarchical assembly alone. Only those systems which assemble promptly into a single cluster (in $\lesssim1$ Myr) from a close separation (all within $\lesssim2$ pc) could match the observed density profile of HD 97950 after $\approx70$% gas dispersal. These results therefore suggest that NGC 3603 young cluster has formed essentially monolithically, i.e., either in-situ or via a prompt assembly, followed by a substantial residual gas expulsion. Both scenarios are consistent with the inferred young age and the small age spread of this cluster. Future observations of molecular cloud filaments with ALMA and proper motion measurements of young clusters with Gaia will provide more direct tests of such birth environments. [Abridged]

The present work is connected with the investigation of the origin and properties of compact astrophysical objects endowed with strangeness, with the objective of finding out their relevance in the formation and evolution of the universe. In the first part of the thesis, Chap.~1-3, we discuss a model, proposed by us, to describe the propagation of small lumps of Strange Quark Matter (SQM) or strangelets, through the Terrestrial atmosphere. The theoretical results were found to be well correlated with exotic cosmic ray events characterized by very low charge to mass ratio. In the next part, we have investigated the other end of the mass spectrum of SQM. In Chap 5, we have developed an analytical expression for the Chandrasekhar Limit of Strange Quark Stars. The limit is found to depend on the fundamental constants (including the bag constant). In the last chapter we have endeavored to show that the quark nuggets, surviving the quark-hadron phase transition in the millisecond era of the early Universe can provide the required closure density and can merge to form compact quark matter objects, whose maximum mass would be governed by the formulation laid out in the preceding chapter. We have also found that these Cold Dark Matter objects can explain the recent astronomical observations of MACHOS by gravitational micro-lensing techniques in the Large Magellanic clouds in the Halo of our Galaxy.

We propose a method to estimate the location of the Sun in the disk of the Milky Way using a method based on the Hellinger distance and construct confidence sets on our estimate of the unknown location using a bootstrap based method. Assuming the Galactic disk to be two-dimensional, the sought solar location then reduces to the radial distance separating the Sun from the Galactic center and the angular separation of the Galactic center to Sun line, from a pre-fixed line on the disk. On astronomical scales, the unknown solar location is equivalent to the location of us earthlings who observe the velocities of a sample of stars in the neighborhood of the Sun. This unknown location is estimated by undertaking pairwise comparisons of the estimated density of the observed set of velocities of the sampled stars, with densities estimated using synthetic stellar velocity data sets generated at chosen locations in the Milky Way disk according to four base astrophysical models. The "match" between the pair of estimated densities is parameterized by the affinity measure based on the familiar Hellinger distance. We perform a novel cross-validation procedure to establish a desirable "consistency" property of the proposed method.

Inflationary and hence quantum origin of primordial perturbations is on a firmer ground than ever post the BICEP2 observations of primordial gravitational waves. One crucial ingredient of success of this paradigm rests on explaining the observed classicality of cosmological inhomogeinities despite their quantum origin. Although decoherence provides a partial understanding of this issue, the question of single outcome motivates the analysis of quantum collapse models in cosmological context which generically modify the dynamics of primordial perturbations and hence can leave their imprints on observables. We revisit one such recently proposed working model of classicalization by spontaneous collapse \cite{Das:2013qwa} in the light of BICEP2 observations to look for possible modifications to tensor power spectra and their implications. We show that it can potentially change the consistency relation of single-field models and a precise measurement of $n_T$ and its running could serve as a test of such dynamics in the early universe.

Understanding how distinct, near-spherical gas-free clusters of very young, massive stars shape out of vast, complex clouds of molecular hydrogen is one of the biggest challenges in astrophysics. A popular thought dictates that a single gas cloud fragments into many new-born stars which, in turn, energize and rapidly expel the residual gas to form a gas-free cluster. This study demonstrates that the above classical paradigm remarkably reproduces the well-observed central, young cluster (HD 97950) of the Galactic NGC 3603 star-forming region, in particular, its shape, internal motion and the mass distribution of stars, naturally and consistently follow from a single model calculation. Remarkably, the same parameters (star formation efficiency, gas expulsion time scale and delay) reproduce HD 97950 as were found to reproduce the Orion Nebula Cluster, Pleiades and R136. The present results thereby provide intriguing evidences of formation of star clusters through single-starburst events followed by significant residual gas expulsion.

R144 is a recently confirmed very massive, spectroscopic binary which appears isolated from the core of the massive young star cluster R136. The dynamical ejection hypothesis as an origin for its location is claimed improbable by Sana et al. due to its binary nature and high mass. We demonstrate here by means of direct N-body calculations that a very massive binary system can be readily dynamically ejected from a R136-like cluster, through a close encounter with a very massive system. One out of four N-body cluster models produces a dynamically ejected very massive binary system with a mass comparable to R144. The system has a system mass of $\approx$ 355 Msun and is located at 36.8 pc from the centre of its parent cluster, moving away from the cluster with a velocity of 57 km/s at 2 Myr as a result of a binary-binary interaction. This implies that R144 could have been ejected from R136 through a strong encounter with an other massive binary or single star. In addition, we discuss all massive binaries and single stars which are ejected dynamically from their parent cluster in the N-body models.

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

Based on kinematic data observed for very young, massive clusters that appear to be in dynamical equilibrium, it has recently been argued that such young systems set examples where the early residual gas-expulsion did not happen or had no dynamical effect. The intriguing scenario of a star cluster forming through a single starburst has thereby been challenged. Choosing the case of the R136 cluster of the Large Magellanic Cloud, the most cited one in this context, we perform direct N-body computations that mimic the early evolution of this cluster including the gas-removal phase (on a thermal timescale). Our calculations show that under plausible initial conditions as consistent from observational data, a large fraction (> 60%) of a gas-expelled, expanding R136-like cluster is bound to regain dynamical equilibrium by its current age. Therefore, the recent measurements of velocity dispersion in the inner regions of R136, that indicate that the cluster is in dynamical equilibrium, are consistent with an earlier substantial gas expulsion of R136 followed by a rapid re-virialization (in approx. 1 Myr). Additionally, we find that the less massive Galactic NGC 3603 Young Cluster (NYC), with a substantially longer re-virialization time, is likely to be found deviated from dynamical equilibrium at its present age (approx. 1 Myr). The recently obtained stellar proper motions in the central part of the NYC indeed suggest this and are consistent with the computed models. This work significantly extends previous models of the Orion Nebula Cluster which already demonstrated that the re-virialization time of young post-gas-expulsion clusters decreases with increasing pre-expulsion density.

Compressible isothermal magnetohydrodynamic turbulence is analyzed under the assumption of statistical homogeneity and in the asymptotic limit of large kinetic and magnetic Reynolds numbers. Following Kolmogorov we derive an exact relation for some two-point correlation functions which generalizes the expression recently found for hydrodynamics. We show that the magnetic field brings new source and flux terms into the dynamics which may act on the inertial range similarly as a source or a sink for the mean energy transfer rate. The introduction of a uniform magnetic field simplifies significantly the exact relation for which a simple phenomenology may be given. A prediction for axisymmetric energy spectra is eventually proposed.

Context: The shape of the stellar initial mass function (IMF) of a star cluster near its upper mass limit is a focal topic of investigation as it determines the high mass stellar content and hence the dynamics of the cluster at its embedded phase as well as during its young gas-free phase. The massive stellar content of a young cluster, however, can be substantially modified due to the dynamical ejections of the massive stars so that the present-day high-mass stellar mass function (hereafter MF) can be different than that with which the cluster is born. Aims: In the present study, we provide a preliminary estimate of this evolution of the high-mass IMF of a young cluster due to early ejections of massive stars, using the Large Magellanic Cloud massive, young cluster R136 as an example. Methods: To that end, we utilize the results of the state-of-the-art calculations by Banerjee, Kroupa & Oh (2012) comprising direct N-body computations of realistic, binary-rich, mass-segregated models of R136. In particular, these calculations provide the ejection fraction of stars as a function of stellar mass. Results: We find that if the measured IMF of R136 is granted to be canonical, as observations indicate, then the "true" high-mass IMF of R136 at its birth must be at least moderately top-heavy when corrected for the dynamical escape of massive stars. Conclusions: The top-heaviness of the true high-mass IMF over the observationally determined one is a general feature of massive, young clusters where the dynamical ejection of massive stars is efficient. We discuss its implications and possible improvements over our current estimate.

[abridged] Among the most remarkable features of the stellar population of R136, the central, young, massive star cluster in the 30 Doradus complex of the Large Magellanic Cloud, are the single stars whose masses substantially exceed the canonical stellar upper mass limit of 150 M_sun. A recent study by us, viz., that of Banerjee, Kroupa & Oh (2012; Paper I), which involves realistic N-body computations of star clusters mimicking R136, indicates that such "super-canonical" (SC) stars can be formed out of a dense stellar population with a canonical initial mass function (IMF) through dynamically induced mergers of the most massive binaries. Here we study the formation of SC stars in the R136 models of Paper I in detail. To avoid forming extraneous SC stars from initially highly eccentric primordial binaries as in Paper I, we compute additional models with only initially circular primordial binaries. We also take into account the mass-evolution of the SC stars using detailed stellar evolutionary models that incorporate updated treatments of stellar winds. We find that SC stars begin to form via dynamical mergers of massive binaries from approx. 1 Myr cluster age. We obtain SC stars with initial masses up to approx. 250 M_sun from these computations. Multiple SC stars are found to remain bound to the cluster simultaneously within a SC-lifetime. These properties of the dynamically formed SC stars are consistent with those observed in R136. In fact, the stellar evolutionary models of SC stars imply that had they formed primordially along with the rest of the R136 cluster, i.e., violating the canonical upper limit, they would have evolved below the canonical 150 M_sun limit by approx. 3 Myr, the likely age of R136, and would not have been observable as SC stars at the present time in R136. This strongly supports the dynamical formation scenario of the observed SC stars in R136.

We conduct a theoretical study on the ejection of runaway massive stars from R136 --- the central massive, star-burst cluster in the 30 Doradus complex of the Large Magellanic Cloud. Specifically, we investigate the possibility of the very massive star (VMS) VFTS 682 being a runaway member of R136. Recent observations of the above VMS, by virtue of its isolated location and its moderate peculiar motion, have raised the fundamental question whether isolated massive star formation is indeed possible. We perform the first realistic N-body computations of fully mass-segregated R136-type star clusters in which all the massive stars are in primordial binary systems. These calculations confirm that the dynamical ejection of a VMS from a R136-like cluster, with kinematic properties similar to those of VFTS 682, is common. Hence the conjecture of isolated massive star formation is unnecessary to account for this VMS. Our results are also quite consistent with the ejection of 30 Dor 016, another suspected runaway VMS from R136. We further note that during the clusters' evolution, mergers of massive binaries produce a few single stars per cluster with masses significantly exceeding the canonical upper-limit of 150 solar mass. The observations of such single super-canonical stars in R136, therefore, do not imply an IMF with an upper limit greatly exceeding the accepted canonical 150 solar mass limit, as has been suggested recently, and they are consistent with the canonical upper limit.

Among the most explored directions in the study of dense stellar systems is the investigation of the effects of the retention of supernova remnants, especially that of the massive stellar remnant black holes (BHs), in star clusters. By virtue of their eventual high central concentration, these stellar mass BHs potentially invoke a wide variety of physical phenomena, the most important ones being emission of gravitational waves (GWs), formation of X-ray binaries, and modification of the dynamical evolution of the cluster. Here we propose, for the first time, that rapid removal of stars from the outer parts of a cluster by the strong tidal field in the inner region of our Galaxy can unveil its BH sub-cluster, which appears as a star cluster that is gravitationally bound by an invisible mass. We study the formation and properties of such systems through direct N-body computations and estimate that they can be present in significant numbers in the inner region of the Milky Way. We call such objects "dark star clusters" (DSCs) as they appear dimmer than normal star clusters of similar mass and they comprise a predicted, new class of entities. The finding of DSCs will robustly cross-check BH retention; they will not only constrain the uncertain natal kicks of BHs, thereby the widely debated theoretical models of BH formation, but will also pinpoint star clusters as potential sites for GW emission for forthcoming ground-based detectors such as the Advanced LIGO. Finally, we also discuss the relevance of DSCs for the nature of IRS 13E.

Compressible isothermal turbulence is analyzed under the assumption of homogeneity and in the asymptotic limit of a high Reynolds number. An exact relation is derived for some two-point correlation functions which reveals a fundamental difference with the incompressible case. The main difference resides in the presence of a new type of term which acts on the inertial range similarly as a source or a sink for the mean energy transfer rate. When isotropy is assumed, compressible turbulence may be described by the relation, $- {2 \over 3} \epsilon_{\rm{eff}} r = {\cal F}_r(r)$, where ${\cal F}_r$ is the radial component of the two-point correlation functions and $\epsilon_{\rm{eff}}$ is an effective mean total energy injection rate. By dimensional arguments we predict that a spectrum in $k^{-5/3}$ may still be preserved at small scales if the density-weighted fluid velocity, $\rho^{1/3} \uu$, is used.

We investigate the dynamics of stellar-mass black holes (BH) in star clusters focusing on the dynamical formation of BH-BH binaries, which are very important sources of gravitational waves (GW). We examine the properties of these BH-BH binaries through direct N-body computations of Plummer clusters, having initially N(0) <= 10^5 low mass stars and a population of stellar mass BHs, using the state-of-the-art N-body integrator "NBODY6". We find that the stellar mass BHs segregate rapidly into the cluster core to form a central dense sub-cluster of BHs in which BH-BH binaries form via 3-body encounters. While most of the BH binaries finally escape from the cluster by recoils due to super-elastic encounters with the single BHs, we find that for clusters with N(0) >= 5 X 10^4, typically a few of them dynamically harden to the extent that they can merge via GW emission within the cluster. Also, for each of such clusters, there are a few escaped BH binaries that merge within a Hubble time, most of the mergers happening within a few Gyr of cluster evolution. These results imply that the intermediate-aged massive clusters constitute the most important class of star clusters that can produce dynamical BH-BH mergers at the present epoch. The BH-BH merger rates obtained from our computations imply a significant detection rate (~ 30/yr) for the "Advanced LIGO" GW detector that will become operative in the near future. Finally, we briefly discuss our ongoing development on this work incorporating the formation of BHs in star clusters from stellar evolution. In particular, we highlight the effect of stellar metallicity on the BH sub-cluster driven expansion of a star cluster's core.

We study the dynamics of stellar-mass black holes (BH) in star clusters with particular attention to the formation of BH-BH binaries, which are interesting as sources of gravitational waves (GW). We examine the properties of these BH-BH binaries through direct N-body simulations of star clusters using the GPU-enabled NBODY6 code. We perform simulations of N <= 10^5 Plummer clusters of low-mass stars with an initial population of BHs. Additionally, we do several calculations of star clusters confined within a reflective boundary mimicking only the core of a massive cluster. We find that stellar-mass BHs with masses ~ 10 solar mass segregate rapidly into the cluster core and form a sub-cluster of BHs within typically 0.2 - 0.5 pc radius, which is dense enough to form BH-BH binaries through 3-body encounters. While most BH binaries are ejected from the cluster by recoils received during super-elastic encounters with the single BHs, few of them harden sufficiently so that they can merge via GW emission within the cluster. We find that for clusters with $N \ga 5\times 10^4$, typically 1 - 2 BH-BH mergers occur within them during the first ~ 4 Gyr of evolution. Also for each of these clusters, there are a few escaping BH binaries that can merge within a Hubble time, most of the merger times being within a few Gyr. These results indicate that intermediate-age massive clusters constitute the most important class of candidates for producing dynamical BH-BH mergers. Old globular clusters cannot contribute significantly to the present-day BH-BH merger rate since most of the mergers from them would have occurred earlier. In contrast, young massive clusters are too young to produce significant number of BH-BH mergers. Our results imply significant BH-BH merger detection rates for the proposed "Advanced LIGO" GW detector. (Abridged)

We continue exploration of the Boltzmann scheme started in Banerjee and Ghosh (2007, henceforth Paper I) for studying the evolution of compact-binary populations of globular clusters, introducing in this paper our method of handling the stochasticity inherent in dynamical processes of binary formation, destruction and hardening in globular clusters. We describe these stochastic processes as "Wiener processes", whereupon the Boltzmann equation becomes a stochastic partial differential equation, the solution of which requires the use of "Ito calculus" (this use being the first, to our knowledge, in this subject), in addition to ordinary calculus. We focus on the evolution of (a) the number of X-ray binaries $N_{XB}$ in globular clusters, and (b) the orbital-period distribution of these binaries. We show that, although the details of the fluctuations in the above quantities differ from one "realization" to another of the stochastic processes, the general trends follow those found in the continuous-limit study of Paper I, and the average result over many such realizations is close to the continuous-limit result. We investigate the dependence of $N_{XB}$ found by these calculations on two essential globular-cluster parameters, namely, the star-star and star-binary encounter-rate parameters $\Gamma$ and $\gamma$, for which we had coined the name Verbunt parameters in Paper I. We compare our computed results with those from CHANDRA observations of Galactic globular clusters, showing that the expected scalings of $N_{XB}$ with the Verbunt parameters are in good agreement with the observed ones. We indicate what additional features can be incorporated into the scheme in future, and how more elaborate problems can be tackled.

We explore a Boltzmann scheme for studying the evolution of compact binary populations of globular clusters. We include processes of compact-binary formation by tidal capture and exchange encounters, binary destruction by dissociation and other mechanisms, and binary hardening by encounters, gravitational radiation and magnetic braking, as also the orbital evolution during mass transfer, following Roche lobe contact. For the encounter processes which are stochastic in nature, we study the probabilistic, continuous limit in this introductory work, deferring the specific handling of the stochastic terms to the next step. We focus on the evolution of (a) the number of X-ray sources N_{XB} in globular clusters, and (b) the orbital-period distribution of the X-ray binaries, as a result of the above processes. We investigate the dependence of N_{XB} on two essential cluster properties, namely, the star-star and star-binary encounter-rate parameters 'Gamma' and 'gamma', which we call Verbunt parameters. We compare our model results with observation, showing that the model values of N_{XB} and their expected scaling with the Verbunt parameters are in good agreement with results from recent X-ray observations of Galactic globular clusters, encouraging us to build more detailed models.

We consider essential mechanisms for orbit-shrinkage or "hardening" of compact binaries in globular clusters to the point of Roche-lobe contact and X-ray emission phase, focussing on the process of collisional hardening due to encounters between binaries and single stars in the cluster core. The interplay between this kind of hardening and that due to emission of gravitational radiation produces a characteristic scaling of the orbit-shrinkage time with the single-star binary encounter rate $\gamma$ in the cluster which we introduce, clarify, and explore. We investigate possible effects of this scaling on populations of X-ray binaries in globular clusters within the framework of a simple "toy" scheme for describing the evolution of pre-X-ray binaries in globular clusters. We find the expected qualitative trends sufficiently supported by data on X-ray binaries in galactic globular clusters to encourage us toward a more quantitative study.

Recent astrophysical observations indicate that the universe is composed of a large amount of dark energy (DE) responsible for an accelerated expansion of the universe, along with a sizeable amount of cold dark matter (CDM), responsible for structure formation. At present, the explanations for the origin or the nature of both CDM and DE seem to require ideas beyond the standard model of elementary particle interactions. Here, for the first time, we show that CDM and DE can arise entirely from the standard principles of strong interaction physics and quantum entanglement. Quantitative agreement with the present data obtains without the need for any adjustable parameters.

It is now believed that the universe is composed of a small amount of the normal luminous matter, a substantial amount of matter (Cold Dark Matter: CDM) which is non-luminous and a large amount of smooth energy (Dark Energy: DE). Both CDM and DE seem to require ideas beyond the standard model of particle interactions. In this work, we argue that CDM and DE can arise entirely from the standard principles of strong interaction physics out of the same mechanism.

We propose that the cold dark matter (CDM) is composed entirely of quark matter, arising from a cosmic quark-hadron transition. We show that compact gravitational objects, with masses around 0.5 (M_{\odot}), could have evolved out of the such CDM.

The existence of compact gravitational lenses, with masses around 0.5 (M_{\odot}), has been reported in the halo of the Milky Way. The nature of these dark lenses is as yet obscure, particularly because these objects have masses well above the threshold for nuclear fusion. In this work, we show that they find a natural explanation as being the evolutionary product of the metastable false vacuum domains (the so-called strange quark nuggets) formed in a first order cosmic quark-hadron transition.

We propose that the cold dark matter (CDM) is composed entirely of quark matter, arising from a cosmic quark-hadron transition. We denote this phase as "quasibaryonic", distinct from the usual baryons. We show that compact gravitational lenses, with masses around 0.5 (M_{\odot}), could have evolved out of the quasibaryonic CDM.

There have been several reports of exotic nuclear fragments, with highly unusual charge to mass ratio, in cosmic ray experiments. Although there exist experimental uncertainties which make them, at best, only candidate "exotic" events, it is important to understand what they could be, if they are eventually confirmed. Among other possible explanations, some authors have interpreted them to be lumps of strange quark matter (strangelets). A major problem with such an interpretation is that to reach the earth's surface, they must possess an unusually high penetrability through the terrestrial atmosphere. We show that a recently proposed mechanism for the propagation of strangelets through the earth's atmosphere, together with a proper account of charge capture and ionisation loss, would solve this problem. We also argue that this could lead to viable strategies for definitive detection of strange quark matter in cosmic ray flux using a ground based large area array of passive detectors.

The mechanism for the propagation of strangelets with low baryon number through the atmosphere of the Earth has been explored. It has been shown that under suitable initial conditions, such strangelets may indeed reach depths near mountain altitudes with mass numbers and charges close to the observed values in cosmic ray experiments.

One of the abiding mysteries in the so-called standard cosmological model is the nature of the dark matter. It is universally accepted that there is an abundance of matter in the universe which is non-luminous, due to their very weak interaction, if at all, with the other forms of matter, excepting of course the gravitational attraction. Speculations as to the nature of dark matter are numerous, often bordering on exotics, and searches for such exotic matter is a very active field of astroparticle physics at the dawn of the new century. Nevertheless, in recent years, there has been experimental evidence for at least one form of dark matter - the massive compact halo objects detected through gravitational microlensing effects proposed by Paczynski some years ago. To date, no clear consensus as to what these objects, referred to in the literature as well as in the following by the acronym MACHO, are made of; for a brief discussion of some of the suggestions, see below. In this work, we show that they find a natural explanation as leftover relics from the putative first order cosmic quark - hadron phase transition that is predicted by the standard model of particle interactions to have occurred during the microsecond epoch of the early universe.

The Chandrasekhar limit for quark stars is evaluated from simple energy balance relations, as proposed by Landau for white dwarfs or neutron stars. It has been found that the limit for quark stars depends on, in addition to the fundamental constants, the Bag constant.

The theory of inertial induction, which has been quite successful in explaining a number of hitherto unexplained astrophysical phenomena, predicts a redshift in light grazing massive objects. The total solar eclipse on 11 August 1999 offers an opportunity of testing the theory. In this paper we predict the additional redshift of a few stars during the eclipse and urge professional astronomers to do the observation.