Neste trabalho estudamos a dinâmica de discos torcidos finos e espessos para compreender melhor a propagação da deformação nestes discos. No caso dos discos finos, estudamos a física do efeito Bardeen-Petterson e aplicamos este modelo para explicar o jato em escalas de parsec e kilo-parsec da galáxia NGC 1275. Encotramos que o efeito Bardeen-Petterson reproduziu muito bem a forma do jato e com isto derivamos os parâmetros do disco como raio, valores das viscosidades azimutal e vertical, lei de potência da densidade superficial e spin do buraco negro. Para uma melhor compreensão da física destes discos, realizamos simulações GRMHD de discos moderadamente finos tanto planos como inclinados para estudar a evolução do ângulo de inclinação entre os momentos angular do buraco negro e do disco de acresção assim como o ângulo de torção que está associado com a precessão do disco. Encontramos que quando o disco de acresção e o buraco negro rotacionam no mesmo sentido, o ângulo de inclinação entre os momentos angular apresentou um comportamento oscilatório na parte interna do disco e permaneceu constante na parte externa em acordo com as previsões teóricas. Já quando o buraco negro rotacina no sentido oposto ao disco de acresção, encontramos pela primeira vez numa simulação GRMHD evidências de alinhamento, ocorrendo um alinhamento de 10\\% do angulo entre os momentos angulares do disco e buraco negro. Além disso, comprovamos pela primeira vez numa simulação GRMHD a não isotropia do stress. Utilizando um modelo semi-analítico, comparamos os resultados de nossas simulações com este modelo, utilizando os dados da simulações de disco plano como entrada e obitivemos os mesmos comportamentos das simulações tanto no caso prógrado quanto no caso retrógrado mostrando que o alinhamento é devido ao regime onda. / In this work we studied the dynamics of twisted thin and thick disks to better understand how the warp propagates in these discs. In the case of thin discs, we studied the physics of the Bardeen-Petterson effect and we applied this model to explain the shape of the jet in both parsec and kilo-parsec scales of the galaxy NGC 1275. We found that the Bardeen-Petterson effect could explain very well the shape of the jet and with that we derived the disc parameters such as its radius, the values of the kinematic azimutal and vertical viscosities, the power-law of the surface density and the spin of the black hole. To better understand the physics of such discs, we have performed GRMHD simulations of moderatelly thin tilted disks to study the evolution of the tilt angle between the angular momentum of the accretion disk and black hole and also the twist angle which is associated with the precession of the disc. We found that when the accretion disc and the black hole are rotating in the same direction, the tilt angle showed an oscillatory behavior in the inner parts of the disk while in the outer parts it remained constant in agreement with the theorical modelos. However, when both rotate in the opposite direction, we found for the very first time in a GRMHD simulation, evidences of alignment of 10\\% of the tilt angle. Besides that, we prove for the first time in a GRMHD simulation that the stress is far from being isotropic. Using a semi-analitic model, we compared the results of our simulations with this model, using the datas of the untilted simulations as inputs and we found the same behaviors found in the simulations even in prograde case as in the retrograde case showing that the alignment is due to bending waves.

Neste trabalho estudamos a dinâmica de discos torcidos finos e espessos para compreender melhor a propagação da deformação nestes discos. No caso dos discos finos, estudamos a física do efeito Bardeen-Petterson e aplicamos este modelo para explicar o jato em escalas de parsec e kilo-parsec da galáxia NGC 1275. Encotramos que o efeito Bardeen-Petterson reproduziu muito bem a forma do jato e com isto derivamos os parâmetros do disco como raio, valores das viscosidades azimutal e vertical, lei de potência da densidade superficial e spin do buraco negro. Para uma melhor compreensão da física destes discos, realizamos simulações GRMHD de discos moderadamente finos tanto planos como inclinados para estudar a evolução do ângulo de inclinação entre os momentos angular do buraco negro e do disco de acresção assim como o ângulo de torção que está associado com a precessão do disco. Encontramos que quando o disco de acresção e o buraco negro rotacionam no mesmo sentido, o ângulo de inclinação entre os momentos angular apresentou um comportamento oscilatório na parte interna do disco e permaneceu constante na parte externa em acordo com as previsões teóricas. Já quando o buraco negro rotacina no sentido oposto ao disco de acresção, encontramos pela primeira vez numa simulação GRMHD evidências de alinhamento, ocorrendo um alinhamento de 10\\% do angulo entre os momentos angulares do disco e buraco negro. Além disso, comprovamos pela primeira vez numa simulação GRMHD a não isotropia do stress. Utilizando um modelo semi-analítico, comparamos os resultados de nossas simulações com este modelo, utilizando os dados da simulações de disco plano como entrada e obitivemos os mesmos comportamentos das simulações tanto no caso prógrado quanto no caso retrógrado mostrando que o alinhamento é devido ao regime onda. / In this work we studied the dynamics of twisted thin and thick disks to better understand how the warp propagates in these discs. In the case of thin discs, we studied the physics of the Bardeen-Petterson effect and we applied this model to explain the shape of the jet in both parsec and kilo-parsec scales of the galaxy NGC 1275. We found that the Bardeen-Petterson effect could explain very well the shape of the jet and with that we derived the disc parameters such as its radius, the values of the kinematic azimutal and vertical viscosities, the power-law of the surface density and the spin of the black hole. To better understand the physics of such discs, we have performed GRMHD simulations of moderatelly thin tilted disks to study the evolution of the tilt angle between the angular momentum of the accretion disk and black hole and also the twist angle which is associated with the precession of the disc. We found that when the accretion disc and the black hole are rotating in the same direction, the tilt angle showed an oscillatory behavior in the inner parts of the disk while in the outer parts it remained constant in agreement with the theorical modelos. However, when both rotate in the opposite direction, we found for the very first time in a GRMHD simulation, evidences of alignment of 10\\% of the tilt angle. Besides that, we prove for the first time in a GRMHD simulation that the stress is far from being isotropic. Using a semi-analitic model, we compared the results of our simulations with this model, using the datas of the untilted simulations as inputs and we found the same behaviors found in the simulations even in prograde case as in the retrograde case showing that the alignment is due to bending waves.

We show our work to refine the process of evolutions in general relativistic magnetohydrodynamics. We investigate several areas in order to improve the overall accuracy of our results. We test several versions of conversion methodologies between different sets of variables. We compare both single equation and two equations solvers to do the conversion. We find no significant improvement for multiple equation conversion solvers when compared to single equation solvers. We also investigate the construction of initial data and the conversion of coordinate systems between initial data code and evolution code. In addition to the conversion work, we have improved some methodologies to ensure data integrity when moving data from the initial data code to the evolution code. Additionally we add into the system of MHD equations a new field to help control the no monopole constraint. We perform a characteristic decomposition of the system of equations in order to derive the associated boundary condition for this new field. Finally, we implement a WENO (weighted non-oscillatory) system. This is done so we can evolve and track shocks that are generated during an evolution of our GRMHD equations.

This thesis is concerned with several problems in general relativity and low energy string theory that are pertinent to the time evolution of the gravitational field. We present a formulation of the Einstein field equations in terms of variational techniques borrowed from geometric analysis. These equations yield the evolution equations for the Cauchy problems of both general relativity and low energy string theory. We then proceed to investigate the evolutionary linear stability of Schwarzschild-like solutions in higher dimensional relativity called black strings. These objects are determined to be linearly unstable. This motivates a further stability analysis of the charged p-brane solutions of low energy string theory. We show that one can eliminate linear instabilities in p-branes for sufficiently large values of charge. We also consider the characteristic problem of general relativistic magnetohydrodynamics (GRMHD). We compute the eigenvalues and eigenvectors of GRMHD and establish degeneracy conditions. Finally, we consider the initial value problem for axisymmetric GRMHD. We formulate the general Einstein and MHD equations under the assumption of a stationary axisymmetric spacetime without assuming the circularity condition.

A white dwarf is most likely the end stage of a low mass star like our Sun, which results when the parent star consumes all the hydrogen in its core, thus bringing fusion to a halt. It is a dense and compact object, where the inward gravitational pull is balanced by the outward pressure arising due to the motion of its constituent degenerate electrons. The theory of non-magnetized and non-rotating white dwarfs was formulated extensively by S. Chandrasekhar in the 1930s, who also proposed a maximum possible mass for this objects, known as the Chandrasekhar limit (Chandrasekhar 1935)1. White dwarfs are believed to be the progenitors of extremely bright explosions called type Ia supernovae (SNeIa). SNeIa are extremely important and popular astronomical events, which are hypothesized to be triggered in white dwarfs having mass close to the famous Chandrasekhar limit ∼ 1.44M⊙. The characteristic nature of the variation of luminosity with time of SNeIa is believed to be powered by the decay of 56Ni to 56Co and, finally, to 56Fe. This feature, along with the consistent mass of the exploding white dwarf, is deeply linked with their utilization as “standard candles” for cosmic distance measurement. In fact, SNeIa measurements were instrumental in establishing the accelerated nature of the current expansion of the universe (Perlmutter et al. 1999). However, several recently observed peculiar SNeIa do not conform to this traditional explanation. Some of these SNeIa are highly over-luminous, e.g. SN 2003fg, SN 2006gz, SN 2007if, SN 2009dc (Howell et al. 2006; Scalzo et al. 2010), and some others are highly under-luminous, e.g. SN 1991bg, SN 1997cn, SN 1998de, SN 1999by, SN 2005bl (Filippenko et al. 1992; Taubenberger et al. 2008). The luminosity of the former group of SNeIa implies a huge Ni-mass (often itself super-Chandrasekhar), invoking highly super-Chandrasekhar white dwarfs, having mass 2.1 − 2.8M⊙, as their most plausible progenitors (Howell et al. 2006; Scalzo et al. 2010). On the other hand, the latter group produces as low as ∼ 0.1M⊙ of Ni (Stritzinger et al. 2006), which rather seem to favor sub-Chandrasekhar explosion scenarios. In this thesis, as the title suggests, we have endeavored to establish the existence of exotic, super- and sub-Chandrasekhar limiting mass white dwarfs, in order to explain the aforementioned peculiar SNeIa. This is an extremely important puzzle to solve in order to comprehensively understand the phenomena of SNeIa, which in turn is essential for the correct interpretation of the evolutionary history of the universe. Effects of magnetic field: White dwarfs have been observed to be magnetized, having surface fields as high as 105 − 109 G (Vanlandingham et al. 2005). The interior field of a white dwarf cannot be probed directly but it is quite likely that it is several orders of magnitude higher than the surface field. The theory of weakly magnetized white dwarfs has been investigated by a few authors, however, their properties do not starkly contrast with that of the non-magnetized cases (Ostriker & Hartwick 1968). In our venture to find a fundamental basis behind the formation of super-Chandrasekhar white dwarfs, we have explored in this thesis the impact of stronger magnetic fields on the properties of white dwarfs, which has so far been overlooked. We have progressed from a simplistic to a more rigorous, self-consistent model, by adding complexities step by step, as follows: • spherically symmetric Newtonian model with constant (central) magnetic field • spherically symmetric general relativistic model with varying magnetic field • model with self-consistent departure from spherical symmetry by general relativis-tic magnetohydrodynamic (GRMHD) numerical modeling. We have started by exploiting the quantum mechanical effect of Landau quanti-zation due to a maximum allowed equipartition central field greater than a critical value Bc = 4.414 × 1013 G. To begin with, we have carried out the calculations in a Newtonian framework assuming spherically symmetric white dwarfs. The primary ef-fect of Landau quantization is to stiffen the equation of state (EoS) of the underlying electron degenerate matter in the high density regime, and, hence, yield significantly super-Chandrasekhar white dwarfs having mass much & 2M⊙ (Das & Mukhopadhyay 2012a,b). Consequently, we have proposed a new mass limit for magnetized white dwarfs which may establish the aforementioned peculiar, over-luminous SNeIa as new standard candles (Das & Mukhopadhyay 2013a,b). We have furthermore predicted possible evo-lutionary scenarios by which super-Chandrasekhar white dwarfs could form by accretion on to a commonly observed magnetized white dwarf, by invoking the phenomenon of flux freezing, subsequently ending in over-luminous, super-Chandrasekhar SNeIa (Das et al. 2013). Before moving on to a more complex model, we have justified the assumptions in our simplistic model, in the light of various related physics issues (Das & Mukhopad-hyay 2014b), and have also clarified, and, hence, removed some serious misconceptions regarding our work (Das & Mukhopadhyay 2015c). Next, we have considered a more self-consistent general relativistic framework. We have obtained stable solutions of magnetostatic equilibrium models for white dwarfs pertaining to various magnetic field profiles, however, still in spherical symmetry. We have showed that in this framework, a maximum stable mass as high as ∼ 3.3M⊙ can be realized (Das & Mukhopadhyay 2014a). However, it is likely that the anisotropic effect due to a strong magnetic field may cause a deformation in the spherical structure of the white dwarfs. Hence, in order to most self-consistently take into account this departure from spherical symmetry, we have constructed equilibrium models of strongly magnetized, static, white dwarfs in a general relativistic framework, first time in the literature to the best of our knowledge. In order to achieve this, we have modified the GRMHD code XNS (Pili et al. 2014), to apply it in the context of white dwarfs. Interestingly, we have found that signifi-cantly super-Chandrasekhar white dwarfs, in the range ∼ 1.7 − 3.4M⊙, are obtained for many possible field configurations, namely, poloidal, toroidal and mixed (Das & Mukhopadhyay 2015a). Furthermore, due to the inclusion of deformation caused by a strong magnetic field, super-Chandrasekhar white dwarfs are obtained for relatively lower central magnetic field strengths (∼ 1014 G) compared to that in the simplistic model — as correctly speculated in our first work of this series (Das & Mukhopadhyay 2012a). We have also found that although the characteristic deformation induced by a purely toroidal field is prolate, the overall shape remains quasi-spherical — justifying our earlier spherically symmetric assumption while constructing at least some models of strongly magnetized white dwarfs (Das & Mukhopadhyay 2014a). Indeed more accurate and extensive numerical analysis seems to have validated our analytical findings. Thus, very interestingly, our investigation has established that magnetized white dwarfs can indeed have mass that significantly exceeds the Chandrasekhar limit, irre-spective of the origin of the underlying magnetic effect — a discovery which is not only of theoretical importance, but also has a direct astrophysical implication in explaining the progenitors of the peculiar, over-luminous, super-Chandrasekhar SNeIa. Effects of modified Einstein’s gravity: A large array of models has been required to explain the peculiar, over- and under- luminous SNeIa. However, it is unlikely that nature would seek mutually antagonistic scenarios to exhibit sub-classes of apparently the same phenomena, i.e., triggering of thermonuclear explosions in white dwarfs. Hence, driven by the aim to establish a unification theory of SNeIa, we have invoked in the last part of this thesis a modification to Einstein’s theory of general relativity in white dwarfs. The validity of general relativity has been tested mainly in the weak field regime, for example, through laboratory experiments and solar system tests. However, the question remains, whether general relativity requires modification in the strong gravity regime, such as, the expanding universe, the region close to a black hole and neutron star. For instance, there is evidence from observational cosmology that the universe has undergone two epochs of cosmic acceleration, the theory behind which is not yet well understood. The period of acceleration in the early universe is known as inflation, while the current accelerated expansion is often explained by invoking a mysterious dark energy. An alternative approach to explain the mysteries of inflation and dark energy is to modify the underlying gravitational theory itself, as it conveniently avoids involving any exotic form of matter. Several modified gravity theories have been proposed which are extensions of Einstein’s theory of general relativity. A popular class of such theories is known as f (R) gravity (e.g. see de Felice & Tsujikawa 2010), where the Lagrangian density f of the gravitational field is an arbitrary function of the Ricci scalar R. In the context of astrophysical compact objects, so far, modified gravity theories have been applied only to neutron stars, which are much more compact than white dwarfs, in order to test the validity of such theories in the strong field regime (e.g. Cooney et al. 2010; Arapoˇglu et al. 2011). Moreover, a general relativistic correction itself does not seem to modify the properties of a white dwarf appreciably when compared to Newtonian calculations. Our venture of exploring modified gravity in white dwarfs in this thesis, is a first in the literature to the best of our knowledge. We have exploited the advantage that white dwarfs have over neutron stars, i.e., their EoS is well established. Hence, any change in the properties of white dwarfs can be solely attributed to the modification of the underlying gravity, unlike in neutron stars, where similar effects could be produced by invoking a different EoS. We have explored a popular, yet simple, model of f (R) gravity, known as the Starobinsky model (Starobinsky 1980) or R−squared model, which was originally pro-posed to explain inflation. Based on this model, we have first shown that modified gravity reproduces those results which are already explained in the paradigm of general relativity (and Newtonian framework), namely, low density white dwarfs in this context. This is a very important test of the modified gravity model and is furthermore necessary to constrain the underlying model parameter. Next, depending on the magnitude and sign of a single model parameter, we have not only obtained both highly super-Chandrasekhar and highly sub-Chandrasekhar limiting mass white dwarfs, but we have also established them as progenitors of the peculiar, over- and under-luminous SNeIa, respectively (Das & Mukhopadhyay 2015b). Thus, an effectively single underlying the-ory unifies the two apparently disjoint sub-classes of SNeIa, which have so far hugely puzzled astronomers. To summarize, in the first part of the thesis, we have established the enormous significance of magnetic fields in white dwarfs in revealing the existence of significantly super-Chandrasekhar white dwarfs. These super-Chandrasekhar white dwarfs could be ideal progenitors of the peculiar, over-luminous SNeIa, which can, hence, be used as new standard candles of cosmic distance measurements. In the latter part of the thesis, we have established the importance of a modified theory of Einstein’s gravity in revealing both highly super- and highly sub-Chandrasekhar limiting mass white dwarfs. We have furthermore demonstrated how such a theory can serve as a missing link between the peculiar, super- and sub-Chandrasekhar SNeIa. Thus, the significance of the current thesis lies in the fact that it not only questions the uniqueness of the Chandrasekhar mass-limit for white dwarfs, but it also argues for the need of a modified theory of Einstein’s gravity to explain astrophysical observations.