My research spans a range of topics in gravitational physics
and theoretical astrophysics. The ultimate goal of my work
is to understand strong field gravitation and solve
long-standing astrophysical puzzles such as the nature of
the progenitors of short gamma-ray bursts, the origin of
X-shaped radio galaxies, the nature of the equation of state
at and above the nuclear saturation densities, and the way
through which planets may form around isolated pulsars, to
name a few. I am interested in studying compact objects as
multimessenger sources, i.e., as sources of gravitational
wave, electromagnetic and neutrino signals. For this reason,
compact object binaries, such as black hole-black
hole (BHBH), neutron star-neutron star (NSNS), black
hole-neutron star (BHNS) , and white dwarf-neutron star
(WDNS) systems are a major theme of my research.

The era of gravitational wave
astronomy and astrophysics has arrived! On Thursday,
February 11, 2016, the LIGO-VIRGO Scientific
collaboration announced
to the world the first direct detection of a gravitational
wave signal! The detection took place on September 14,
2015 at 09:50:45 UTC! What is more, the signal is very
well matched by the inspiral and merger of two black holes
with masses 36 and 29 times the mass of our Sun! This
detection not only provides proof that binary black holes
exist, but that compact (black) objects exist that can pack
more than \(25M_\odot\) in a tiny volume! These are the
heaviest stellar mass black holes ever observed. After the
news were released, I launched a simulation of non-spinning
BHs to simulate this system by solving the full non-linear
Einstein equations in vacuum. A visualization of the
(apparent) horizons of the two black holes and the cross
polarization of the gravitational waves from this source can
be view below

This detection is only the beginning. More gravitational
waves (GWs) from these (and possibly unknown) sources will
soon be detected by excuisite devices. These include the
ground-based detectors, such as
advanced LIGO,
VIRGO
and the
near-future KAGRA,
all targeting stellar mass compact binaries, Pulsar Timing
Arrays such as
NANOGrav, and the
IPTA, targeting
supermassive black hole binaries, and by future space based
detectors such
as eLISA,
targeting massive black hole binaries and galactic compact
binaries (see Fig. 1 for a more complete
list of gravitational wave sources and detectors). Moreover,
transient electromagnetic signatures arising from such
sources will be detectable by current and future telescopes
such as
Fermi,
PanSTARRS,
JWST, and
LSST. These observations
in conjunction with careful theoretical modeling will
provide unprecedented insights into the workings of the
fabric of spacetime.

Astrophysical general relativity

A common thread uniting the astrophysical problems I am
interested in is the crucial role of relativistic
gravitation and I have particular expertise in dynamical
spacetime scenarios. For example, compact objects in
binaries or isolation are a class of astrophysical systems
for which relativistic gravitation is of paramount
importance.

Mergers of NSNS and BHNS binaries are promising
sources of gravitational waves. GWs from such binaries
can potentially constrain the unknown equation of the state
of the above the nuclear saturation density.

In a recent work of mine I discovered that a one-arm (m=1)
instability can develop in the hypermassive neutron star
that forms following a dynamical capture binary neutron star
merger. The following movie demonstrates the last stages of
the merger and the development of the one-arm instability.

Due to this instability strong GWs are emitted, and if the
remnant is long-lived and the instability persists, such GWs
could be detectable by aLIGO at 10Mpc and by the future
Einstein Telescope at 100Mpc. These GWs could potentially be
powerful discriminators of the equation of state of the
matter at super-nuclear densities. You can read more about
this work in my
article arXiv:1510.03432.

NSNS and BHNS binaries as true multimessenger systems, not
only emit GWs, but can also power a host of electromagnetic
signatures both prior to and following merger that will
accompany the GWs. These binaries may ultimately also solve
long standing astrophysical puzzles such as the nature of
the progenitors of short gamma-ray bursts (sGRB), and the
origin of r-process elements in the Universe. I am
interested in solving these puzzles and have spent
considerable effort trying to prove from a theoretical point
of view that BHNS and NSNS are viable progenitors of sGRB
engines.
The following movie shows the generation of a collimated,
mildly relativistic, Poynting-dominated outflow - an
incipient jet - following the merger of a black hole with a
magnetized neutron star. Relativistic gravitation is used
for all these simulations coupled to the equations of ideal
magnetohydrodynamics. The black hole to neutron star mass
ratio is 3:1 and the black hole initial dimensionless spin
parameter is 0.75. You can read more about this work in my
Letter ApJ,
806, L14 (2015).

A visualization of the accompanying gravitational waves
generated during the inspiral and merger of the previous
BHNS system is shown in the following movie

Apart from these quasicircular mergers, eccentric
mergers of these compact binaries can take place in dense
stellar systems, such as globular clusters. However, neutron
stars in globular clusters tend to be rapidly spinning as
they are typically observed as millisecond pulsars. This is
because globular cluster favor the formation of low-mass
X-ray binaries which are believed to be the progenitors of
millisecond pulsars. In these mergers the spin of the
neutron star can no longer be ignored. In my Letter
ApJ.,
807, L3 (2015), the effects of the neutron star spin
in black hole-neutron star eccentric mergers were studied,
finding that even moderately large neutron star spins can
have a large impact on the amount of neutron-rich matter
that escapes to infinity following tidal disruption of the
neutron star by the black hole, and which will shine
electromagnetically due to the decay of r-process
elements.

Supermassive BHBH binaries emit not only copious
amounts of detectable GWs, but also shine
electromagnetically when surrounded by accretion disks (as
is likely to be the case at the centers of distant active
galactic nuclei which formed after galactic mergers) due to
the binary-disk interaction. There is a world wide effort
tagetted at observing such "binary black hole AGNs", and
some exciting candidates already have been found. Such
systems can drive spectacular galactic jets as is
demonstrated in my articles
Phys.
Rev. D 89 , 064060 (2014), and
Phys.
Rev. D 90, 104030 (2014). A movie showing the
generation of such jets can be viewed below.

A movie of the previous system in equatorial and meridional
2D slices is shown below. It becomes clear that following
merger the jet accelerates and become more strongly
magnetized. The front that propagates changes the features
of the jet from its foot to its head and this phenomenon may
help observers identify past black hole - black hole mergers
based on jet morphology alone.

Moreover, BHBH-disk systems can potentially explain the
origin
of X-shaped
radio galaxies because of spin-flip of the single black
hole that forms following merger (work in progress), and
their gravitational waves can even be used as standard
sirens to determine the Hubble constant to very high
precision, thus constraining cosmological dark energy
models. A visualization of the gravitational waves from
these sources (in the limit where the self-gravity of the
disk can be ignored) can viewed in the following movie

Mergers of WDNS binaries are sources of
low-frequency gravitational waves. The ultimate fate of
massive WDNS mergers is currently unknown. These could lead
either to the formation a Thorn-Zytkow-like object
surrounded by an extended accretion disk or collapse to form
a BH or even explode because of rapid nuclear burning. If
the remnant neither collapses nor explodes, and the disk
survives and fragments, such mergers may offer a plausible
route to planet formation around pulsars. Due to the vast
range of lengthscales and timescales involved in this
problem, simulating these systems in full general relativity
while resolving at the same time both the neutron star
(which is necessary to ascertain if the remnant will
collapse to a black hole) and the white dwarf, is
prohibitive with current resources. However, with a careful
modification of the equation of state to scale down the
white dwarf branch while leaving the neutron star branch
intact and maintaining all inequalities between timescales
and lengthscales unchanged, one can study these mergers with
current resources. Calling pseudo-white dwarfs these scaled
down white dwarf models, and constructing a sequence of
decreasing white dwarf compaction (\(M/R\)), one can obtain
a good answer as to what happens to the neutron star core
following merger. I have devised such an equation of state
in my
article Phys. Rev. D83,
064002 (2011). A movie showing the merger of a neutron
star with a pseudo-white dwarf using hydrodynamic
simulations in full general relativity can viewed below. So
far, these simulations indicate that such mergers do not
lead to the prompt formation of a black hole and that they
are supported against gravitational collapse primarily due
to the centrifugal support from differential rotation and
not due to thermal pressure.

The initial data problem is an important problem one
needs to solve to prepare the initial conditions for all
these compact object (binary merger) simulations. This
involves solving the elliptic Einstein constraints - a
highly non-trivial task, especially when (quasi)equilibrium
configurations are the desired solutions. In my
articles Phys. Rev. D83,
064002 (2011)
and Phys. Rev. D
84, 104032 (2011), I developed such an elliptic solver
for corotating binaries. More recently, in my Letter
ApJ.,
807, L3 (2015) a method was developed for generating
eccentric black hole-neutron star initial data treating
rapidly spinning neutron stars. I am interested in the
development of more sophisticated/efficient methods for
generating generic initial data for relativistic
computations.

Non-astrophysical relativistic gravitation

In addition to the above topics, I have recently grown a
keen interest in the stability of black objects (black
holes, black strings, black rings) in higher than 4
dimensions. Moreover, I have always been fascinated by the
mathematical properties of the partial differential
equations describing classical theories of gravitation, such
as their classification (hyperbolic, parabolic, elliptic or
mixed type), and the well-posedness of the initial value
problem. While we know that general relativity admits a
well-posed Cauchy problem, this is currently unknown for
other competing theories of gravitation. Work of mine
that falls in this category can be found in my articles
Phys. Rev. D
75, 024026 (2007)
and Phys. Rev. D
78, 024002 (2008). I have also done some mathematical
work on f(R) theories of gravity in my
article Class. Quantum
Grav. 28 085006 (2011).
Finally, formulating new tests capable of constraining these
alternative gravity theories using the playground provided
by compact objects is among my main interests.

Methods

I put considerable effort in designing computer models for
compact objects in binaries or isolation, and spend time on
smashing compact objects in the computer. The ultimate goal
is to theoretically predict gravitational wave,
electromagnetic, and neutrino signatures from these systems,
so that ultimately we may better understand strong field
gravitation and compact object physics.

To achieve this goal, I develop and use both
(semi)analytic methods and state-of-the-art codes at
supercomputer centers, to study the equations of
relativistic gravitation. These codes are capable of solving
the Einstein equations of general relativity coupled to
Maxwell's equations of electromagnetism, the equations of
relativistic (magneto)hydrodynamics for the motion of the
matter, and those governing radiation transport.

An important theoretical product of my research was the
development of a robust electromagnetic gauge condition,
which I dubbed "Generalized Lorenz gauge". This new gauge
allows for long-term and stable numerical integration of the
equations of general relativistic ideal magnetohydrodynamics
(GRMHD) using a vector potential formulation with mesh refinement
techniques. You can read more about this work in my article
Phys. Rev. D
85, 024013 (2012), where I carried out a mathematical
analysis of electromagnetic gauge conditions and simulations
were performed to show their effects, and in my Letter
Phys. Rev. Lett. 109,
221102 (2012), where I proposed the new Generalized
Lorenz gauge condition.

The Generalized Lorenz gauge is now an essential ingredient
for the Illinois relativity group GRMHD code, of which I am
a main co-developer and co-maintainer, and it has already
been adopted by other groups. Recently, my collaborators and
I (the lion's share credit for re-writing this code goes to
Zach Etienne) embedded this code in the publicly
available Einstein
Toolkit, and this work has been published
in Class. Quantum
Grav. 32 (2015) 175009. More about this code can be
found on
the IllinoisGRMHD
website. I am also interested in the development of highly
accurate computational methods for general relativistic
magnetohydrodynamics on structured and unstructured grids.

Collaborators

I have had the fortune to work with wonderful people and here is an incomplete
list of my collaborators.