Featured Research


2D Relativistic MHD Simulations of Kruskal-Schwarzschild Instability (KSI) in a Strongly Magnetized Striped Wind

Figure: Development of the instability from the linear to the non-linear phase and dripping of hot plasma from the current layer. The unmagnetized hot current layer is surrounded by cold strongly magnetized fluid with anti-aligned magnetic field lines going into and out of the page in, respectively, the top and bottom layers.

We study the linear and non-linear development of the Kruskal-Schwarzchild Instability in a relativisitically expanding striped wind. This instability is the generalization of Rayleigh-Taylor instability in the presence of a magnetic field. It has been suggested to produce a self-sustained acceleration mechanism in strongly magnetized outflows found in active galactic nuclei, gamma-ray bursts, and micro-quasars. The instability leads to magnetic reconnection, but in contrast with
steady-state Sweet-Parker reconnection, the dissipation rate is not limited by the current layer`s small aspect ratio. We performed two-dimensional (2D) relativistic magneto-hydrodynamic (RMHD) simulations featuring two cold and highly magnetized plasma layers with an anti-parallel magnetic field separated by a thin layer of relativistically hot plasma with a local effective gravity induced by the outflow`s acceleration. Our simulations show how the heavier relativistically hot plasma in the reconnecting layer drips out and allows oppositely oriented magnetic field lines to reconnect. The instability`s growth rate in the linear regime matches the predictions of linear stability analysis. We find turbulence rather than an ordered bulk flow near the reconnection region, with turbulent velocities up to 0.1c, largely independent of model parameters. However, the magnetic energy dissipation rate is found to be much slower. This occurs due to the slow evacuation of hot plasma from the current layer, largely because of the Kelvin-Helmholtz instability experienced by the dripping plasma.

See my latest article on the KSI: Gill, R., Granot, J., & Lyubarsky, Y. 2018, MNRAS, 474, 3535 


Gamma-Ray Bursts

Figure: Formation of GRB spectrum using my one-zone time-dependent kinetic code. Here, I start by injecting a quasi-thermal soft seed spectrum and let it evolve in time while accounting for Compton scattering, pair production/annihilation, Coulomb interaction and volumetric heating of pairs.
My recent study with Chris Thompson at CITA and two other works co-authored by me present a complete description of the physics that gives rise to the GRB spectrum. We show that the low energy spectrum forms pre-breakout when the jet, powered by a central black hole in our model, is mildly relativistic and working its way out of the confining medium provided by the Wolf-Rayet progenitor. The inertia in the jet is dominated by the magnetic field and the low-energy spectrum is formed by quasi-thermal Comptonization in a photon-rich electron-positron-pair fireball, with no baryon contamination. The high-energy spectrum forms post-breakout when the jet is relativistic and optically thin. The dissipation of magnetic energy is provided by the interaction of the relativistic magnetofluid with baryons entrained in the jet from the confining medium. In both the optically thick and thin phases, the pairs are heated volumetrically due to the damping of hydromagnetic turbulence in the flow. To simulate the high energy spectrum, I developed a one-zone kinetic code with all radiative processes treated in an exact manner. We have developed a cogent explanation for the prompt emission mechanism in a series of articles:


Magnetar Giant Flares

The soft gamma-ray repeaters (SGRs) dissipate large amounts of magnetic energy in the form of highly luminous gamma-ray bursts called giant flares. How the flares are triggered and what controls their duration and energy is still an open question. The energetics of the bursts, however, do point towards the extremely strong fields that these objects are endowed with. According to the Thompson & Duncan (1995) model of magnetar bursts, the giant flares involve major restructuring of the global magnetic field of the star on timescales longer than the Alfven crossing time of the magnetosphere. This takes place in concert with a magnetic reconnection event in the magnetosphere that ultimately triggers the hyperflare.

Figure: It displays the setup of the different reconnecting current layers. The macroscopic Sweet-Parker layer with length L∼1E5 cm and width δ∼0.01 cm is the largest of the three. This layer is then thinned down vertically as strong magnetic flux is convected into the dissipation region. The Hall reconnection layer, represented by the dark gray region, develops when δ becomes comparable to the ion-inertial length. The system makes a transition from the slow to the impulsive  reconnection and powers the main flare. The tiny region embedded inside the Sweet-Parker layer is  the super-hot turbulent current layer, which aids in creating sufficient anomalous resistivity to facilitate the formation of the  Sweet-Parker layer. The strongly accelerated plasma downstream of the reconnection layer is trapped inside magnetic flux lines and forms a  plasmoid moving at some speed V. This plasmoid is then finally ejected during the initial spike when the external field undergoes a sudden relaxation.


One interesting observation these models do not resolve is the  connection between the precursor burst and the main flare. In two out of the total three giant flares that have been  observed, a precursor burst of very short duration was detected a few seconds before the main flare. In my study, I show that the  precursor bursts are not only causally connected to the main flares but also triggers them. I examine two reconnection models  that explain how the wound up flux tubes inside the NS and the sheared field lines in the magnetosphere can form tangential  discontinuities and dissipate stored magnetic energy. In this work, I introduced a novel idea of Hall reconnection  operating in the inner magnetosphere which ultimately leads to an explosive energy release.

Work on magnetars:


Axions and Axion-like Particles

Unlike weak interactions, that are responsible for the radioactive decay of subatomic particles, strong interactions have thus far been experimentally verified to preserve CP-symmetry. However, this observation is in conflict with QCD which necessarily violates CP. This can be remedied in many ways but the one most favored solution to the Strong CP problem has been the introduction of a new particle into the Standard Model, dubbed the axion. In several extensions of the Standard Model axion-like particles (ALPs) are produced naturally. Such particles remain elusive due to their very small mass and extremely weak coupling to matter and radiation. The limit on its coupling to radiation, derived from astrophysical arguments, in the given mass range is a few orders of magnitude larger than that predicted by theoretical axion models. Several laboratory experiments have been conducted and many are currently underway to exclude as much of the parameter space as possible.

In my work on ALPs, I take a novel approach to constrain the properties of ALPs. I use the polarized light of magnetic white dwarfs (mWDs) to determine the absolute upper limit on the ALPs' coupling to radiation as a function of its mass. My treatment is based on the observation that light is polarized as it traverses the magnetosphere encompassing the tenuous plasma permeated by the ∼ 100 MG magnetic field of some mWDs. Since the photon-axion interaction in a strongly magnetized plasma leads to additional polarization, the observed degree of polarization can be used to constrain axion properties if the plasma contribution is correctly modeled. The derived constraints in my work are hitherto  the best estimates of ALP parameters to come from compact objects.

Work on Axions: