Source author record

Rosalba Perna

Rosalba Perna appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

59works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

59 published item(s)

preprint2023arXiv

Close Encounters of Tight Binary Stars with Stellar-mass Black Holes

Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite of 3D hydrodynamics simulations using the moving-mesh code AREPO, we investigate the formation of transient phenomena and their properties in close encounters between an $2M_{\odot}$ or $20M_{\odot}$ equal-mass circular binary star and single $20M_{\odot}$ black hole (BH). Stars can be disrupted by the BH during dynamical interactions, naturally producing electromagnetic transient phenomena. Encounters with impact parameters smaller than the semimajor axis of the initial binary frequently lead to a variety of transients whose electromagnetic signatures are qualitatively different from those of ordinary disruption events involving just two bodies. These include the simultaneous or successive disruptions of both stars and one full disruption of one star accompanied by successive partial disruptions of the other star. On the other hand, when the impact parameter is larger than the semimajor axis of the initial binary, the binary is either simply tidally perturbed or dissociated into bound and unbound single stars ("micro-Hills" mechanism). The dissociation of $20M_{\odot}$ binaries can produce a runaway star and an active BH moving away from one another. Also, the binary dissociation can either produce an interacting binary with the BH, or a non-interacting, hard binary; both could be candidates of BH high- and low-mass X-ray binaries. Hence our simulations especially confirm that strong encounters can lead to the formation of the (generally difficult to form) BH low-mass X-ray binaries.

preprint2022arXiv

Close Encounters of Stars with Stellar-mass Black Hole Binaries

Many astrophysical environments, from star clusters and globular clusters to the disks of Active Galactic Nuclei, are characterized by frequent interactions between stars and the compact objects that they leave behind. Here, using a suite of $3-D$ hydrodynamics simulations, we explore the outcome of close interactions between $1M_{\odot}$ stars and binary black holes (BBHs) in the gravitational wave regime, resulting in a tidal disruption event (TDE) or a pure scattering, focusing on the accretion rates, the back reaction on the BH binary orbital parameters and the increase in the binary BH effective spin. We find that TDEs can make a significant impact on the binary orbit, which is often different from that of pure scattering. Binaries experiencing a prograde (retrograde) TDE tend to be widened (hardened) by up to $\simeq 20\%$. Initially circular binaries become more eccentric by $\lesssim 10\%$ by a prograde or retrograde TDE, whereas the eccentricity of initially eccentric binaries increases (decreases) by a retrograde (prograde) TDE by $\lesssim 5\%$. Overall a single TDE can generally result in changes of the gravitational wave-driven merger time scale by order unity. The accretion rates of both black holes are very highly super-Eddington, showing modulations (preferentially for retrograde TDEs) on a time scale of the orbital period, which can be a characteristic feature of BBH-driven TDEs. Prograde TDEs result in the effective spin parameter $χ$ to vary by $\lesssim 0.02$ while $χ\gtrsim -0.005$ for retrograde TDEs.

preprint2022arXiv

Searching for Ultra-light Bosons and Constraining Black Hole Spin Distributions with Stellar Tidal Disruption Events

Stars that pass close to the supermassive black holes located in the center of galaxies can be violently disrupted by tidal forces, leading to flares that are observed as bright transient events in sky surveys. The rate for these events to occur depends on the black hole spins, which in turn can be affected by ultra-light bosons due to superradiance. We perform a detailed analysis of these effects and show that searches for stellar tidal disruptions have a significant potential to uncover the existence of ultra-light bosons. In particular, we find that upcoming stellar tidal disruption rate measurements by the Vera Rubin Observatory's Legacy Survey of Space and Time can be used to either discover or rule out bosons with masses ranging from $10^{-20}$ to $10^{-18}$ eV. Our analysis also indicates that these measurements may be used to constrain a variety of supermassive black hole spin distributions and determine if close-to maximal spins are preferred.

preprint2022arXiv

The emergence of diffused Gamma-Ray Burst afterglows from the disks of Active Galactic Nuclei

The disks of Active Galactic Nuclei (AGNs) have emerged as rich environments for the production and capture of stars and the compact objects that they leave behind. These stars produce long Gamma-Ray Bursts (LGRBs) at their deaths, while frequent interactions among compact objects form binary neutron stars and neutron star-black hole binaries, leading to short GRBs (SGRBs) upon their merger. Predicting the properties of these transients as they emerge from the dense environments of AGN disks is key to their proper identification and to better constrain the star and compact object population in AGN disks. Some of these transients would appear unusual because they take place in much higher densities than the interstellar medium. Others, which are the subject of this paper, would additionally be modified by radiation diffusion since they are generated within optically thick regions of the accretion disks. Here we compute the GRB afterglow light curves for diffused GRB sources for a representative variety of central black-hole masses and disk locations. We find that the radiation from radio to UV and soft X-rays can be strongly suppressed by synchrotron self-absorption in the dense medium of the AGN disk. In addition, photon diffusion can significantly delay the emergence of the emission peak, turning a beamed, fast transient into a slow, isotropic, and dimmer one. These would appear as broadband-correlated AGN variability with dominance at the higher frequencies. Their properties can constrain both the stellar populations within AGN disks as well as the disk structure.

preprint2021arXiv

Can stellar-mass black hole growth disrupt disks of active galactic nuclei? The role of mechanical feedback

Stellar-mass BHs (sBHs) are predicted to be embedded in active galactic nuclei (AGN) disks due to gravitational drag and in-situ star formation. However, we find that due to a high gas density in an AGN disk environment, compact objects may rapidly grow to intermediate-mass BHs and deplete matter from the AGN disk unless accretion is suppressed by some feedback process(es). These consequences are inconsistent with AGN observations and the dynamics of the Galactic center. Here we consider mechanical feedback mechanisms for the reduction of gas accretion. Rapidly accreting sBHs launch winds and/or jets via the Blandford-Znajek mechanism, which produce high-pressure shocks and cocoons. Such a shock and cocoon can spread laterally in the plane of the disk, eject the outer regions of a circum-sBH disk (CsBD) and puncture a hole in the AGN disk with horizontal size comparable to the disk scale-height. Since the depletion timescale of the bound CsBD is much shorter than the resupply timescale of gas to the sBH, the time-averaged accretion rate onto sBHs is reduced by this process by a factor of $\sim 10$--$100$. This feedback mechanism can therefore help alleviate the sBH over-growth and AGN-disk depletion problems. On the other hand, we find that cocoons of jets can unbind a large fraction of the gas accreting in the disks of less massive SMBHs, which may help explain the dearth of high-Eddington ratio AGNs with SMBH mass $\lesssim10^5{\rm M_\odot}$.

preprint2021arXiv

Host galaxies and electromagnetic counterparts to binary neutron star mergers across the cosmic time: Detectability of GW170817-like events

The detection of electromagnetic radiation (EM) accompanying the gravitational wave (GW) signal from the binary neutron star (BNS) merger GW170817 has revealed that these systems constitute at least a fraction of the progenitors of short gamma-ray bursts (SGRBs). As gravitational wave detectors keep pushing their detection horizons, it is important to assess coupled GW/EM probabilities, and how to maximize observational prospects. Here we perform population synthesis calculations of BNS evolution with the code MOBSE, and seed the binaries in galaxies at three representative redshifts (z=0.01,0.1,1) of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the astrophysical parameters of GRB170817A as a prototype, we numerically compute the broadband lightcurves of jets from BNS mergers, with the afterglow brightness depending on the local medium density at the merger sites. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find that ~70-80% of BNSs from z=0.01 should be detectable in gamma-rays. The afterglow detection probabilities of GW-triggered BNS mergers vary between ~0.3-0.7%, with higher values for jets aligned with the BNS angular momentum, and are comparable across the high and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses.

preprint2021arXiv

Symmetry Breaking in Dynamical Encounters in the Disks of Active Galactic Nuclei

Active galactic nucleus (AGN) disks may be important sites of binary black hole (BBH) mergers. Here we show via numerical experiments with the high-accuracy, high precision code {\tt SpaceHub} that broken symmetry in dynamical encounters in AGN disks can lead to an asymmetry between prograde and retrograde BBH mergers. The direction of the hardening asymmetry depends on the initial binary semi-major axis. An asymmetric distribution of mass-weighted projected spin $χ_{\rm eff}$ should therefore be expected in LIGO-Virgo detections of BBH mergers from AGN disks. This channel further predicts that negative $χ_{\rm eff}$ BBH mergers are most likely for massive binaries.

preprint2020arXiv

Electromagnetic Signatures of Relativistic Explosions in AGN Disks

The disks of Active Galactic Nuclei (AGNs), traditionally studied as the feeders of the supermassive black holes (SMBHs) at their centers, have recently triggered a lot of interest also as hosts to massive stars and hence their neutron star (NS) and black hole (BH) remnants. Migration traps and gas torques in these disks favor binary formation and enhance the rate of compact object mergers. In these environments both long gamma-ray bursts (GRBs) from the death of massive stars and short GRBs from NS-NS and NS-BH mergers are expected. However, their properties in the unique environments of AGN disks have never been studied. Here we show that GRBs in AGNs can display unique features, owing to the unusual relative position of the shocks that characterize the burst evolution and the Thomson photosphere of the AGN disk. In dense environments, for example, the external shock develops before the internal shocks, leading to prompt emission powered by a relativistic reverse shock. The transient's time evolution is also compressed, yielding afterglow emission that is much brighter and peaks much earlier than for GRBs in the interstellar medium. Additionally, in regions of the disk that are sufficiently dense and extended, the light curves are dominated by diffusion, since the fireball is trapped inside the disc photosphere. These diffusion-dominated transients emerge on timescales of days in disks around SMBHs of $\sim 10^6 M_\odot$ to years for SMBHs of $\sim 10^8 M_\odot$. Finally, a large fraction of events, especially in AGNs with SMBHs $\lesssim 10^7 M_\odot$, display time-variable absorption in the X-ray band.

preprint2020arXiv

Electromagnetic transients and gravitational waves from white dwarf disruptions by stellar black holes in triple systems

Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov-Kozai oscillations imposed by the tidal field of a outer tertiary companion orbiting the inner BH-WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semi-major axes and eccentricities of the triples, and stellar mass ratios. We find rates of WD-TDEs in the range $1.2\times 10^{-3}-1.4$ Gpc$^{-3}$ yr$^{-1}$ for $z\leq 0.1$, rarer than stellar TDEs in triples by a factor of $\sim 3$--$30$. The uncertainty in the TDE rates may be greatly reduced in the future using gravitational wave (GW) observations of Galactic binaries and triples with LISA. WD-TDEs may give rise to high energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient which lasts for only days. WD--BH and WD--NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD-TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies.

preprint2020arXiv

Giant planet swaps during close stellar encounters

The discovery of planetary systems outside of the solar system has challenged some of the tenets of planetary formation. Among the difficult-to-explain observations, are systems with a giant planet orbiting a very-low mass star, such as the recently discovered GJ~3512b planetary system, where a Jupiter-like planet orbits an $M$-star in a tight and eccentric orbit. Systems such as this one are not predicted by the core accretion theory of planet formation. Here we suggest a novel mechanism, in which the giant planet is born around a more typical Sun-like star ($M_{*,1}$), but is subsequently exchanged during a dynamical interaction with a flyby low-mass star ($M_{*,2}$). We perform state-of-the-art $N$-body simulations with $M_{*,1}=1M_\odot$ and $M_{*,2}=0.1M_\odot$ to study the statistical outcomes of this interaction, and show that exchanges result in high eccentricities for the new orbit around the low-mass star, while about half of the outcomes result in tighter orbits than the planet had around its birth star. We numerically compute the cross section for planet exchange, and show that an upper limit for the probability per planetary system to have undergone such an event is $Γ\sim 4.4(M_{\rm c}/100M_\odot)^{-2}(a_{\rm p}/{\rm AU}) (σ/1\,{\rm km}\,{\rm s}^{-1})^{5}$Gyr$^{-1}$, where $a_{\rm p}$ is the planet semi-major axis around the birth star, $σ$ the velocity dispersion of the star cluster, and $M_{\rm c}$ the total mass of the star cluster. Hence these planet exchanges could be relatively common for stars born in open clusters and groups, should already be observed in the exoplanet database, and provide new avenues to create unexpected planetary architectures.

preprint2020arXiv

Hot Jupiter and ultra-cold Saturn formation in dense star clusters

The discovery of high incidence of hot Jupiters in dense clusters challenges the field-based hot Jupiter formation theory. In dense clusters, interactions between planetary systems and flyby stars are relatively common. This has a significant impact on planetary systems, dominating hot Jupiter formation. In this paper, we perform high precision, few-body simulations of stellar flybys and subsequent planet migration in clusters. A large parameter space exploration demonstrates that close flybys that change the architecture of the planetary system can activate high eccentricity migration mechanisms: Lidov-Kozai and planet-planet scattering, leading to high hot Jupiter formation rate in dense clusters. Our simulations predict that many of the hot Jupiters are accompanied by "ultra-cold Saturns", expelled to apastra of thousands of AU. This increase is particularly remarkable for planetary systems originally hosting two giant planets with semi-major axis ratios $\sim$ 4 and the flyby star approaching nearly perpendicular to the planetary orbital plane. The estimated lower limit to the hot Jupiter formation rate of a virialized cluster is $\sim 1.6\times10^{-4}(σ/{\rm 1kms^{-1}})^5({a_{\rm p}}/{\rm 20 AU})({M_{\rm c}}/{\rm 1000M_\odot})^{-2}$Gyr$^{-1}$ per star, where $σ$ is the cluster velocity dispersion, $a_{\rm p}$ is the size of the planetary system and $M_{\rm c}$ is the mass of the cluster. Our simulations yield a hot Jupiter abundance which is $\sim$ 50 times smaller than that observed in the old open cluster M67. We expect that interactions involving binary stars, as well as a third or more giant planets, will close the discrepancy.

preprint2020arXiv

Intrinsic properties of the engine and jet that powered the short gamma-ray burst associated with GW170817

GRB 170817A was a subluminous short gamma-ray burst detected about 1.74 s after the gravitational wave signal GW170817 from a binary neutron star (BNS) merger. It is now understood as an off-axis event powered by the cocoon of a relativistic jet pointing 15 to 30 degrees away from the direction of observation. The cocoon was energized by the interaction of the incipient jet with the non-relativistic baryon wind from the merger remnant, resulting in a structured outflow with a narrow core and broad wings. In this paper, we couple the observational constraints on the structured outflow with a model for the jet-wind interaction to constrain the intrinsic properties with which the jet was launched by the central engine, including its time delay from the merger event. Using wind prescriptions inspired by magnetized BNS merger simulations, we find that the jet was launched within about 0.4 s from the merger, implying that the 1.74 s observed delay was dominated by the fireball propagation up to the photospheric radius. We also constrain, for the first time for any gamma-ray burst, the jet opening angle at injection and set a lower limit to its asymptotic Lorentz factor. These findings suggest an initially Poynting-flux dominated jet, launched via electromagnetic processes. If the jet was powered by an accreting black hole, they also provide a significant constraint on the survival time of the metastable neutron star remnant.

preprint2020arXiv

Planetary Architectures in Interacting Stellar Environments

The discovery of Exoplanetary Systems has challenged some of the theories of planet formation, which assume unperturbed evolution of the host star and its planets. However, in star clusters the interactions with flyby stars and binaries may be relatively common during the lifetime of a planetary system. Here, via high-resolution $N$-body simulations of star-planet systems perturbed by interlopers (stars and binaries), we explore the reconfiguration to the planetary system due to the encounters. In particular, via an exploration focused on the strong scattering regime, we derive the fraction of encounters which result in planet ejections, planet transfers and collisions by the interloper star/binary, as a function of the characteristics of the environment (density, velocity dispersion), and for different masses of the flyby star/binary. We find that binary interlopers can significantly increase the cross section of planet ejections and collisions, while they only slightly change the cross section for planet transfers. Therefore, in environments with high binary fractions, floating planets are expected to be relatively common, while in environments with low binary fractions, where the cross sections of planet ejection and transfer are comparable, the rate of planet exchanges between two stars will be comparable to the rate of production of free-floating planets.

preprint2020arXiv

Radio follow-up of a candidate gamma-ray transient in the sky localization area of GW170608

After the identification of a candidate $γ$-ray transient in the error region of the binary black hole (BBH) merger GW150914 by the \textit{Fermi} satellite, the question of whether BBH mergers can be associated to electromagnetic counterparts remains highly debated. Here, we present radio follow-up observations of GW170608, a BBH merger that occurred during the second observing run (O2) of the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO). Our radio follow up focused on a specific field contained in the GW170608 sky localization area, where a candidate high-energy transient was detected by the \textit{Fermi} Large Area Telescope (LAT). We make use of data collected at 1.4\,GHz with the Karl G. Jansky Very Large Array (VLA), as well as with the VLA Low-band Ionosphere and Transient Experiment (VLITE). Our analysis is sensitive to potential radio afterglows with luminosity densities $L_{\rm 1.4\,GHz}\gtrsim 6\times10^{28}$\,erg\,s$^{-1}$\,Hz$^{-1}$. In the most optimistic theoretical models, $\approx 20\%$ of BBH events occurring in massive hosts could be associated with outflows as radio luminous as this. Although we find no evidence for the presence of a radio counterpart associated with the LAT candidate in the GW170608 error region, our analysis demonstrates the feasibility of future radio follow-up observations of well localized BBHs. Comparing our radio upper-limits with theoretical expectations for the radio afterglows potentially associated with jets launched in BBH mergers, we find that for jets of energy $\approx 10^{49}$\,erg seen on-axis, only jet angles $θ_{jet}\gtrsim 40^{\circ}$ are compatible with the observations.

preprint2016arXiv

Binary Neutron Star Mergers and Short Gamma-Ray Bursts: Effects of Magnetic Field Orientation, Equation of State, and Mass Ratio

We present fully GRMHD simulations of the merger of binary neutron star (BNS) systems. We consider BNSs producing a hypermassive neutron star (HMNS) that collapses to a spinning black hole (BH) surrounded by a magnetized accretion disk in a few tens of ms. We investigate whether such systems may launch relativistic jets and power short gamma-ray bursts. We study the effects of different equations of state (EOSs), different mass ratios, and different magnetic field orientations. For all cases, we present a detailed investigation of the matter dynamics and of the magnetic field evolution, with particular attention to its global structure and possible emission of relativistic jets. The main result of this work is that we found the formation of an organized magnetic field structure. This happens independently of EOS, mass ratio, and initial magnetic field orientation. We also show that those models that produce a longer-lived HMNS lead to a stronger magnetic field before collapse to BH. Such larger fields make it possible, for at least one of our models, to resolve the MRI and hence further amplify the magnetic field. However, by the end of our simulations, we do not observe a magnetically dominated funnel and hence neither a relativistic outflow. With respect to the recent simulations of Ruiz et al 2016, we evolve models with lower and more realistic initial magnetic field strengths and, because of computational reasons, we do not evolve the accretion disk for the long timescales that seem to be required in order to see a relativistic outflow. Since all our models produce a similar ordered magnetic field structure, we expect that the results found in Ruiz et al 2016, where they only considered an equal-mass system with an ideal fluid EOS, should be general and, at least from a qualitative point of view, independent from mass-ratio, magnetic field orientation, and EOS.

preprint2016arXiv

Intermediate-mass black holes from Population III remnants in the first galactic nuclei

We report the formation of intermediate-mass black holes (IMBHs) in suites of numerical $N$-body simulations of Population III remnant black holes (BHs) embedded in gas-rich protogalaxies at redshifts $z\gtrsim10$. We model the effects of gas drag on the BHs' orbits, and allow BHs to grow via gas accretion, including a mode of hyper-Eddington accretion in which photon trapping and rapid gas inflow suppress any negative radiative feedback. Most initial BH configurations lead to the formation of one (but never more than one) IMBH in the center of the protogalaxy, reaching a mass of $10^{3-5}\mathrm{M}_{\odot}$ through hyper-Eddington growth. Our results suggest a viable pathway to forming the earliest massive BHs in the centers of early galaxies. We also find that the nuclear IMBH typically captures a stellar-mass BH companion, making these systems observable in gravitational waves as extreme mass-ratio inspirals (EMRIs) with \textit{eLISA}.

preprint2016arXiv

Short Gamma-Ray Bursts from the Merger of Two Black Holes

Short Gamma-Ray Bursts (GRBs) are explosions of cosmic origin believed to be associated with the merger of two compact objects, either two neutron stars, or a neutron star and a black hole. The presence of at least one neutron star has long been thought to be an essential element of the model: its tidal disruption provides the needed baryonic material whose rapid accretion onto the post-merger black hole powers the burst. The recent tentative detection by the Fermi satellite of a short GRB in association with the gravitational wave signal GW150914 produced by the merger of two black holes has challenged this standard paradigm. Here we show that the evolution of two high-mass, low-metallicity stars with main sequence rotational speeds a few tens of percent of the critical speed eventually undergoing a weak supernova explosion {\em can} produce a short gamma-ray burst. The outer layers of the envelope of the last exploding star remain bound and circularize at large radii. With time, the disk cools and becomes neutral, suppressing the magneto-rotational instability, and hence the viscosity. The disk remains 'long-lived dead' until tidal torques and shocks during the pre-merger phase heat it up and re-ignite accretion, rapidly consuming the disk and powering the short gamma-ray burst.

preprint2015arXiv

Constraining the GRB-magnetar model by means of the Galactic pulsar population

A large fraction of Gamma Ray Bursts (GRBs) displays an X-ray plateau phase within <10^{5} s from the prompt emission, proposed to be powered by the spin-down energy of a rapidly spinning newly born magnetar. In this work we use the properties of the Galactic neutron star population to constrain the GRB-magnetar scenario. We re-analyze the X-ray plateaus of all Swift GRBs with known redshift, between January 2005 and August 2014. From the derived initial magnetic field distribution for the possible magnetars left behind by the GRBs, we study the evolution and properties of a simulated GRB-magnetar population using numerical simulations of magnetic field evolution, coupled with Monte Carlo simulations of Pulsar Population Synthesis in our Galaxy. We find that if the GRB X-ray plateaus are powered by the rotational energy of a newly formed magnetar, the current observational properties of the Galactic magnetar population are not compatible with being formed within the GRB scenario (regardless of the GRB type or rate at z=0). Direct consequences would be that we should allow the existence of magnetars and "super-magnetars" having different progenitors, and that Type Ib/c SNe related to Long GRBs form systematically neutron stars with higher initial magnetic fields. We put an upper limit of <16 "super-magnetars" formed by a GRB in our Galaxy in the past Myr (at 99% c.l.). This limit is somewhat smaller than what roughly expected from Long GRB rates, although the very large uncertainties do not allow us to draw strong conclusion in this respect.

preprint2015arXiv

Gamma-Ray Bursts Trace UV Metrics of Star Formation over 3 < z < 5

We present the first uniform treatment of long duration gamma-ray burst (GRB) host galaxy detections and upper limits over the redshift range 3<z<5, a key epoch for observational and theoretical efforts to understand the processes, environments, and consequences of early cosmic star formation. We contribute deep imaging observations of 13 GRB positions yielding the discovery of eight new host galaxies. We use this dataset in tandem with previously published observations of 31 further GRB positions to estimate or constrain the host galaxy rest-frame ultraviolet (UV; 1600 A) absolute magnitudes M_UV. We then use the combined set of 44 M_UV estimates and limits to construct the M_UV luminosity function (LF) for GRB host galaxies over 3<z<5 and compare it to expectations from Lyman break galaxy (LBG) photometric surveys with the Hubble Space Telescope. Adopting standard prescriptions for the luminosity dependence of galaxy dust obscuration (and hence, total star formation rate), we find that our LF is compatible with LBG observations over a factor of 600x in host luminosity, from M_UV = -22.5 mag to >-15.6 mag, and with extrapolations of the assumed Schechter-type LF well beyond this range. We review proposed astrophysical and observational biases for our sample, and find they are for the most part minimal. We therefore conclude, as the simplest interpretation of our results, that GRBs successfully trace UV metrics of cosmic star formation over the range 3<z<5. Our findings suggest GRBs are providing an accurate picture of star formation processes from z ~3 out to the highest redshifts.

preprint2015arXiv

Modelling of the Surface Emission of the Low-Magnetic Field Magnetar SGR 0418+5729

We perform a detailed modelling of the post-outburst surface emission of the low magnetic field magnetar SGR 0418+5729. The dipolar magnetic field of this source, B=6x10^12 G estimated from its spin-down rate, is in the observed range of magnetic fields for normal pulsars. The source is further characterized by a high pulse fraction and a single-peak profile. Using synthetic temperature distribution profiles, and fully accounting for the general-relativistic effects of light deflection and gravitational redshift, we generate synthetic X-ray spectra and pulse profiles that we fit to the observations. We find that asymmetric and symmetric surface temperature distributions can reproduce equally well the observed pulse profiles and spectra of SGR 0418. Nonetheless, the modelling allows us to place constraints on the system geometry (i.e. the angles $ψ$ and $ξ$ that the rotation axis makes with the line of sight and the dipolar axis, respectively), as well as on the spot size and temperature contrast on the neutron star surface. After performing an analysis iterating between the pulse profile and spectra, as done in similar previous works, we further employed, for the first time in this context, a Markov-Chain Monte-Carlo approach to extract constraints on the model parameters from the pulse profiles and spectra, simultaneously. We find that, to reproduce the observed spectrum and flux modulation: (a) the angles must be restricted to $65°< ψ+ξ< 125°$ or $235°< ψ+ξ<295°$; (b) the temperature contrast between the poles and the equator must be at least a factor of $\sim6$, and (c) the size of the hottest region ranges between 0.2-0.7 km (including uncertainties on the source distance). Last, we interpret our findings within the context of internal and external heating models.

preprint2015arXiv

NuSTAR J095551+6940.8: a highly magnetised neutron star with super-Eddington mass accretion

The identification of the Ultraluminous X-ray source (ULX) X-2 in M82 as an accreting pulsar has shed new light on the nature of a subset of ULXs, while rising new questions on the nature of the super-Eddington accretion. Here, by numerically solving the torque equation of the accreting pulsar within the framework of the magnetically threaded-disk scenario, we show that three classes of solutions, corresponding to different values of the magnetic field, are mathematically allowed. We argue that the highest magnetic field one, corresponding to B $\sim 10^{13}$ G, is favoured based on physical considerations and the observed properties of the source. In particular, that is the only solution which can account for the observed variations in $\dot{P}$ (over four time intervals) without requiring major changes in $\dot{M}$, which would be at odds with the approximately constant X-ray emission of the source during the same time. For this solution, we find that the source can only accomodate a moderate amount of beaming, 0.5 $\lesssim b < 1$. Last, we show that the upper limit on the luminosity, L$_X < 2.5 \times 10^{38}$ erg s$^{-1}$ from archival observations, is consistent with a highly-magnetized neutron star being in the propeller phase at that time.

preprint2015arXiv

Producing Magnetar Magnetic Fields in the Merger of Binary Neutron Stars

The merger of binary neutron stars (BNSs) can lead to large amplifications of the magnetic field due to the development of turbulence and instabilities in the fluid, such as the Kelvin-Helmholtz shear instability, which drive small-scale dynamo activity. In order to properly resolve such instabilities and obtain the correct magnetic field amplification, one would need to employ resolutions that are currently unfeasible in global general relativistic magnetohydrodynamic (GRMHD) simulations of BNS mergers. Here, we present a subgrid model that allows global simulations to take into account the small-scale amplification of the magnetic field which is caused by the development of turbulence during BNS mergers. Assuming dynamo saturation, we show that magnetar-level fields ($\sim 10^{16}\,{\rm G}$) can be easily reached, and should therefore be expected from the merger of magnetized BNSs. The total magnetic energy can reach values up to $\sim 10^{51}\,{\rm erg}$ and the post-merger remnant can therefore emit strong electromagnetic signals and possibly produce short gamma-ray bursts.

preprint2015arXiv

The accretion regimes of a highly magnetised NS: the unique case of NuSTAR J095551+6940.8

We analyze archival Chandra HRC observations of the ultra luminous accreting pulsar M82-X2 (NuSTAR J095551+6940.8), and determine an upper limit of $< 1.7\times 10^{38}$~erg/s to its luminosity at an epoch at which it was undetected. Combined with other recent measurements, this confirms that the source X-ray emission has been highly variable during the last 15 years, ranging from a maximum of $10^{40}$ erg/s through intermediate values $\sim$ a few $\times 10^{39}$ erg/s, and down to a minimum that must be below the current detection threshold $\sim (2-3) \times 10^{38}$ erg/s . We interpret these results by means of a magnetically-threaded disk model: when at peak luminosity, the neutron star (NS) is close to spin equilibrium, its inner disk edge r_m ~ 10^8 cm is approximately half the corotation radius r_{co}, and radiation pressure dominates the disk out to r_{tr} ~ 10^9 cm. In the radiation pressure-dominated regime, r_m grows very slowly as the mass inflow rate drops: as a result, r_m < r_{co} remains valid until the mass accretion rate becomes ~ the Eddington accretion rate, allowing a wide range of accretion luminosities to the NS. Once the mass accretion rate is below Eddington, accretion onto the NS is inhibited because r_m > r_{co}, and the source luminosity is expected to drop by a large factor. We conclude that a magnetically threaded, radiation pressure-dominated disk, around a highly magnetized NS (B~10^{13} G) offers the best intepretation for all the currently observed properties of NuSTAR J095551+6940.8. This source offers an unprecedented opportunity to study the disk-magnetosphere interaction in a new regime of supercritical accretion, and across the transition between-radiation pressure and gas-pressure dominance inside the disk.

preprint2015arXiv

The imprint of the cosmic supermassive black hole growth history on the 21 cm background radiation

The redshifted 21 cm transition line of hydrogen tracks the thermal evolution of the neutral intergalactic medium (IGM) at "cosmic dawn," during the emergence of the first luminous astrophysical objects (~100 Myr after the Big Bang) but before these objects ionized the IGM (~400-800 Myr after the Big Bang). Because X-rays, in particular, are likely to be the chief energy courier for heating the IGM, measurements of the 21 cm signature can be used to infer knowledge about the first astrophysical X-ray sources. Using analytic arguments and a numerical population synthesis algorithm, we argue that the progenitors of supermassive black holes (SMBHs) should be the dominant source of hard astrophysical X-rays---and thus the primary driver of IGM heating and the 21 cm signature---at redshifts $z < 20$, if (i) they grow readily from the remnants of Population III stars and (ii) produce X-rays in quantities comparable to what is observed from active galactic nuclei and high-mass X-ray binaries. We show that models satisfying these assumptions dominate over contributions to IGM heating from stellar populations, and cause the 21 cm brightness temperature to rise at $z > 20$. An absence of such a signature in the forthcoming observational data would imply that SMBH formation occurred later (e.g. via so-called direct collapse scenarios), that it was not a common occurrence in early galaxies and protogalaxies, or that it produced far fewer X-rays than empirical trends at lower redshifts, either due to intrinsic dimness (radiative inefficiency) or Compton-thick obscuration close to the source.

preprint2015arXiv

The luminosity and stellar mass functions of GRB host galaxies: Insight into the metallicity bias

Long-Duration Gamma-Ray Bursts (GRBs) are powerful probes of the Universe star formation history, but correlation between the two depends on the highly debated presence/strength of a metallicity bias. To investigate this correlation, we use a phenomenological model that successfully describes star formation rates, luminosities and stellar masses of star forming galaxies, applying it to GRB production. We predict luminosities, stellar masses, and metallicities of host galaxies depending on the metallicity bias. Our best-fitting model includes a moderate metallicity bias, broadly consistent with the large majority of long-duration GRBs in metal-poor environments originating from collapsars (probability ~83%), but with a secondary contribution (~17%) from metal-independent production channels, such as binary evolution. Because of the mass-metallicity relation of galaxies, the maximum likelihood model predicts that the metal-independent channel becomes dominant at z<2, where hosts have higher metallicities and collapsars are suppressed. This possibly explains why some studies find no clear evidence of a metal-bias based on low-z samples. However, while metallicity predictions match observations well at high redshift, there is tension with low redshift observations, since a significant fraction of GRB hosts are predicted to have (near-)solar metallicity. This is in contrast to observations, unless obscured, metal-rich hosts are preferentially missed in current datasets, and suggests that lower efficiencies of the metal-independent GRB channel might be preferred following a comprehensive fit from complete samples. Overall, we are able to establish the presence of a metallicity bias for GRB production, but continued characterization of GRB host galaxies is needed to quantify its strength.

preprint2015arXiv

The Outburst Decay of the Low Magnetic Field Magnetar SWIFT J1822.3-1606: Phase-resolved Analysis and Evidence for a Variable Cyclotron Feature

We study the timing and spectral properties of the low-magnetic field, transient magnetar SWIFT J1822.3-1606 as it approached quiescence. We coherently phase-connect the observations over a time-span of ~500 days since the discovery of SWIFT J1822.3-1606 following the Swift-BAT trigger on 2011 July 14, and carried out a detailed pulse phase spectroscopy along the outburst decay. We follow the spectral evolution of different pulse phase intervals and find a phase and energy-variable spectral feature, which we interpret as proton cyclotron resonant scattering of soft photon from currents circulating in a strong (10^14 G) small-scale component of the magnetic field near the neutron star surface, superimposed to the much weaker (~3 x 10^13 G) magnetic dipole field. We discuss also the implications of the pulse-resolved spectral analysis for the emission regions on the surface of the cooling magnetar.

preprint2014arXiv

Gravitational waves from massive magnetars formed in binary neutron star mergers

Binary neutron star (NS) mergers are among the most promising sources of gravitational waves (GWs), as well as candidate progenitors for short Gamma-Ray Bursts (SGRBs). Depending on the total initial mass of the system, and the NS equation of state (EOS), the post-merger phase can see a prompt collapse to a black hole, or the formation of a supramassive NS, or even a stable NS. In the case of post-merger NS (PMNS) formation, magnetic field amplification during the merger will produce a magnetar with a large induced mass quadrupole moment, and millisecond spin. If the timescale for orthogonalization of the magnetic symmetry axis with the spin axis is sufficiently short the NS will radiate its spin down energy primarily via GWs. Here we study this scenario for various outcomes of NS formation: we generalise the set of equilibrium states for a twisted torus magnetic configuration to include solutions that, at a fixed exterior dipole field, carry a larger magnetic energy reservoir; we hence compute their magnetic ellipticity and the strength of the expected GW signal as a function of the magnitude of the dipole and toroidal field. The relative number of GW detections from PMNSs and from binary NSs is a strong function of the NS equation of state (EOS), being higher (~ 1%) for the stiffest EOSs and negligibly small for the softest ones. For intermediate-stiffness EOSs, such as the n=4/7 polytrope recently used by Giacomazzo \& Perna or the GM1 used by Lasky et al., the relative fraction is ~0.3%; correspondingly we estimate a GW detection rate from stable PMNSs of ~ (0.1-1) yr$^{-1}$ with Advanced detectors, and of ~ (100-1000) yr$^{-1}$ with third generation detectors such as the Einstein Telescope. Measurement of such GW signal would provide strong constraints on the NS EOS and on the nature of the binary progenitors giving rise to SGRBs.

preprint2014arXiv

High Resolution Transmission Spectroscopy as a Diagnostic for Jovian Exoplanet Atmospheres: Constraints from Theoretical Models

We present high resolution transmission spectra of giant planet atmospheres from a coupled 3-D atmospheric dynamics and transmission spectrum model that includes Doppler shifts which arise from winds and planetary motion. We model jovian planets covering more than two orders of magnitude in incident flux, corresponding to planets with 0.9 to 55 day orbital periods around solar-type stars. The results of our 3-D dynamical models reveal certain aspects of high resolution transmission spectra that are not present in simple 1-D models. We find that the hottest planets experience strong substellar to anti-stellar (SSAS) winds, resulting in transmission spectra with net blue shifts of up to 3 km s$^{-1}$, whereas less irradiated planets show almost no net Doppler shifts. Compared to 1-D models, peak line strengths are significantly reduced for the hottest atmospheres owing to Doppler broadening from a combination of rotation (which is faster for close-in planets under the assumption of tidal locking) and atmospheric winds. Finally, high resolution transmission spectra may be useful in studying the atmospheres of exoplanets with optically thick clouds since line cores for very strong transitions should remain optically thick to very high altitude. High resolution transmission spectra are an excellent observational test for the validity of 3-D atmospheric dynamics models, because they provide a direct probe of wind structures and heat circulation. Ground-based exoplanet spectroscopy is currently on the verge of being able to verify some of our modeling predictions, most notably the dependence of SSAS winds on insolation. We caution that interpretation of high resolution transmission spectra based on 1-D atmospheric models may be inadequate, as 3-D atmospheric motions can produce a noticeable effect on the absorption signatures.

preprint2014arXiv

Spectral features in isolated neutron stars induced by inhomogeneous surface temperatures

The thermal X-ray spectra of several isolated neutron stars display deviations from a pure blackbody. The accurate physical interpretation of these spectral features bears profound implications for our understanding of the atmospheric composition, magnetic field strength and topology, and equation of state of dense matter. With specific details varying from source to source, common explanations for the features have ranged from atomic transitions in the magnetized atmospheres or condensed surface, to cyclotron lines generated in a hot ionized layer near the surface. Here we quantitatively evaluate the X-ray spectral distortions induced by inhomogeneous temperature distributions of the neutron star surface. To this aim, we explore several surface temperature distributions, we simulate their corresponding general relativistic X-ray spectra (assuming an isotropic, blackbody emission), and fit the latter with a single blackbody model. We find that, in some cases, the presence of a spurious 'spectral line' is required at a high significance level in order to obtain statistically acceptable fits, with central energy and equivalent width similar to the values typically observed. We also perform a fit to a specific object, RX J0806.4-4123, finding several surface temperature distributions able to model the observed spectrum. The explored effect is unlikely to work in all sources with detected lines, but in some cases it can indeed be responsible for the appearance of such lines. Our results enforce the idea that surface temperature anisotropy can be an important factor that should be considered and explored also in combination with more sophisticated emission models like atmospheres.

preprint2013arXiv

Central compact objects in magnetic lethargy

Central Compact Objects are peculiar young neutron stars, with very low external magnetic fields combined with high fluxes in the X-ray band and surface temperature anisotropies. However, in their crust the magnetic field can be strong, result of its burial during a short post-supernova hypercritical accretion episode. The implications of this latter scenario for the temperature anisotropy, pulsed fraction and luminosity are discussed.

preprint2013arXiv

Constraints on Warm Dark Matter models from high-redshift long gamma-ray bursts

Structures in Warm Dark Matter (WDM) models are exponentially suppressed below a certain scale, characterized by the dark matter particle mass, $m_{\rm x}$. Since structures form hierarchically, the presence of collapsed objects at high-redshifts can set strong lower limits on $m_{\rm x}$. We place robust constraints on $m_{\rm x}$ using recent results from the {\it Swift} database of high-redshift gamma-ray bursts (GRBs). We parameterize the redshift evolution of the ratio between the cosmic GRB rate and star formation rate (SFR) as $\propto (1+z)^α$, thereby allowing astrophysical uncertainties to partially mimic the cosmological suppression of structures in WDM models. Using a maximum likelihood estimator on two different $z>4$ GRB subsamples (including two bursts at $z>8$), we constrain $m_{\rm x} \gtrsim 1.6$-1.8 keV at 95% C.L., when marginalized over a flat prior in $α$. We further estimate that 5 years of a SVOM-like mission would tighten these constraints to $m_{\rm x} \gtrsim 2.3 $ keV. Our results show that GRBs are a powerful probe of high-redshift structures, providing robust and competitive constraints on $m_{\rm x}$.

preprint2013arXiv

Formation of Stable Magnetars from Binary Neutron Star Mergers

By performing fully general relativistic magnetohydrodynamic simulations of binary neutron star mergers, we investigate the possibility that the end result of the merger is a stable magnetar. In particular, we show that, for a binary composed of two equal-mass neutron stars (NSs) of gravitational mass M~1.2 Msun and equation of state similar to Shen et al. at high densities, the merger product is a stable NS. Such NS is found to be differentially rotating and ultraspinning with spin parameter J/M^2~0.86, where J is its total angular momentum, and it is surrounded by a disk of ~0.1 Msun. While in our global simulations the magnetic field is amplified by about two orders of magnitude, local simulations have shown that hydrodynamic instabilities and the onset of the magnetorotational instability could further increase the magnetic field strength up to magnetar levels. This leads to the interesting possibility that, for some NS mergers, a stable and magnetized NS surrounded by an accretion disk could be formed. We discuss the impact of these new results for the emission of electromagnetic counterparts of gravitational wave signals and for the central engine of short gamma-ray bursts.

preprint2013arXiv

On the observed duration distribution of gamma-ray bursts from collapsars

The duration of the prompt emission of long gamma-ray bursts is generally considered to be fairly similar to the duration of the activity of the engine in the center of the progenitor star. Here, we investigate the relation between the duration of the engine activity and that of the observed light curve, using inputs from both numerical simulations and observations. We find that the observed burst duration is a good proxy for the engine duration after the time necessary for the jet to break out the star's surface is subtracted. However, the observed duration is a function of the viewing angle and can be significantly shorter than the duration of the engine activity. We also show that the observed, redshift-corrected burst duration evolves only moderately with redshift for both observations and synthetic light curves. We conclude that the broad distribution of the observed duration of long BATSE gamma-ray bursts is mostly accounted for by an engine lasting ~20 s, the dispersion being due to viewing and redshift effects. Our results do not rule out the existence of engines with very long duration. However, we find that they are constrained to be a small minority of the BATSE detected bursts.

preprint2013arXiv

The fate of fallback matter around newly born compact objects

The presence of fallback disks around young neutron stars has been invoked over the years to explain a large variety of phenomena. Here we perform a numerical investigation of the formation of such disks during a supernova explosion, considering both neutron star (NS) and black hole (BH) remnants. Using the public code MESA, we compute the angular momentum distribution of the pre-supernova material, for stars with initial masses M in the range 13 - 40 Msun, initial surface rotational velocities vsurf between 25% and 75% of the critical velocity, and for metallicities Z of 1%, 10% and 100% of the solar value. These pre SN models are exploded with energies E varying between 10^{50} - 3x10^{52} ergs, and the amount of fallback material is computed. We find that, if magnetic torques play an important role in angular momentum transport, then fallback disks around NSs, even for low-metallicity main sequence stars, are not an outcome of SN explosions. Formation of such disks around young NSs can only happen under the condition of negligible magnetic torques and a fine-tuned explosion energy. For those stars which leave behind BH remnants, disk formation is ubiquitous if magnetic fields do not play a strong role; however, unlike the NS case, even with strong magnetic coupling in the interior, a disk can form in a large region of the {Z,M,vsurf,E} parameter space. Together with the compact, hyperaccreting fallback disks widely discussed in the literature, we identify regions in the above parameter space which lead to extended, long-lived disks around BHs. We find that the physical conditions in these disks may be conducive to planet formation, hence leading to the possible existence of planets orbiting black holes.

preprint2013arXiv

The imprint of the crustal magnetic field on the thermal spectra and pulse profiles of isolated neutron stars

Isolated neutron stars (NSs) show a bewildering variety of astrophysical manifestations, presumably shaped by the magnetic field strength and topology at birth. Here, using state-of-the art calculations of the coupled magnetic and thermal evolution of NSs, we compute the thermal spectra and pulse profiles expected for a variety of initial magnetic field configurations. In particular, we contrast models with purely poloidal magnetic fields to models dominated by a strong internal toroidal component. We find that, while the former displays double peaked profiles and very low pulsed fractions, in the latter, the anisotropy in the surface temperature produced by the toroidal field often results in a single pulse profile, with pulsed fractions that can exceed the 50-60% level even for perfectly isotropic local emission. We further use our theoretical results to generate simulated "observed" spectra, and show that blackbody (BB) fits result in inferred radii that can be significantly smaller than the actual NS radius, even as low as ~ 1-2 km for old NSs with strong internal toroidal fields and a high absorption column density along their line of sight. We compute the size of the inferred BB radius for a few representative magnetic field configurations, NS ages, and magnitudes of the column density. Our theoretical results are of direct relevance to the interpretation of X-ray observations of isolated NSs, as well as to the constraints on the equation of state of dense matter through radius measurements.

preprint2013arXiv

The Transient Gravitational-Wave Sky

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

preprint2013arXiv

Unifying the observational diversity of isolated neutron stars via magneto-thermal evolution models

Observations of magnetars and some of the high magnetic field pulsars have shown that their thermal luminosity is systematically higher than that of classical radio-pulsars, thus confirming the idea that magnetic fields are involved in their X-ray emission. Here we present the results of 2D simulations of the fully-coupled evolution of temperature and magnetic field in neutron stars, including the state-of-the-art kinetic coefficients and, for the first time, the important effect of the Hall term. After gathering and thoroughly re-analysing in a consistent way all the best available data on isolated, thermally emitting neutron stars, we compare our theoretical models to a data sample of 40 sources. We find that our evolutionary models can explain the phenomenological diversity of magnetars, high-B radio-pulsars, and isolated nearby neutron stars by only varying their initial magnetic field, mass and envelope composition. Nearly all sources appear to follow the expectations of the standard theoretical models. Finally, we discuss the expected outburst rates and the evolutionary links between different classes. Our results constitute a major step towards the grand unification of the isolated neutron star zoo.

preprint2012arXiv

Compact Binary Progenitors of Short Gamma-Ray Bursts

In recent years, detailed observations and accurate numerical simulations have provided support to the idea that mergers of compact binaries containing either two neutron stars (NSs) or an NS and a black hole (BH) may constitute the central engine of short gamma-ray bursts (SGRBs). The merger of such compact binaries is expected to lead to the production of a spinning BH surrounded by an accreting torus. Several mechanisms can extract energy from this system and power the SGRBs. Here we connect observations and numerical simulations of compact binary mergers, and use the current sample of SGRBs with measured energies to constrain the mass of their powering tori. By comparing the masses of the tori with the results of fully general-relativistic simulations, we are able to infer the properties of the binary progenitors which yield SGRBs. By assuming a constant efficiency in converting torus mass into jet energy, epsilon_{jet}=10%, we find that most of the tori have masses smaller than 0.01M_{sun}, favoring "high-mass" binary NSs mergers, i.e., binaries with total masses >~1.5 the maximum mass of an isolated NS. This has important consequences for the gravitational-wave signals that may be detected in association with SGRBs, since "high-mass" systems do not form a long-lived hypermassive NS after the merger. While NS-BH systems cannot be excluded to be the engine of at least some of the SGRBs, the BH would need to have an initial spin of ~0.9, or higher.

preprint2012arXiv

Gamma-Ray-Burst Host Galaxy Surveys at Redshift z>4: Probes of Star Formation Rate and Cosmic Reionization

Measuring the star formation rate (SFR) at high redshift is crucial for understanding cosmic reionization and galaxy formation. Two common complementary approaches are Lyman-Break-Galaxy (LBG) surveys for large samples and Gamma-Ray-Burst (GRB) observations for sensitivity to SFR in small galaxies. The z>4 GRB-inferred SFR is higher than the LBG rate, but this difference is difficult to understand, as both methods rely on several modeling assumptions. Using a physically motivated galaxy luminosity function model, with star formation in dark-matter halos with virial temperature Tvir>2e4 K (M_DM>2e8 M_sun), we show that GRB and LBG-derived SFRs are consistent if GRBs extend to faint galaxies (M_AB<-11). To test star formation below the detection limit L_lim~0.05L^*_{z=3} of LBG surveys, we propose to measure the fraction f_det(L>L_lim,z) of GRB hosts with L>L_lim. This fraction quantifies the missing star formation fraction in LBG surveys, constraining the mass-suppression scale for galaxy formation, with weak dependence on modeling assumptions. Because f_det(L>L_lim,z) corresponds to the ratio of star formation rates derived from LBG and GRB surveys, if these estimators are unbiased, measuring f_det(L>L_lim,z) also constrains the redshift evolution of the GRB production rate per unit mass of star formation. Our analysis predicts significant success for GRB host detections at z~5 with f_det(L>L_lim,z)~0.4, but rarer detections at z>6. By analyzing the upper limits on host-galaxy luminosities of six z>5 GRBs from literature data, we infer that galaxies with M_AB>-15 were present at z>5 at 95% confidence, demonstrating the key role played by very faint galaxies during reionization.

preprint2012arXiv

General Relativistic Simulations of Accretion Induced Collapse of Neutron Stars to Black Holes

Neutron stars (NSs) in the astrophysical universe are often surrounded by accretion disks. Accretion of matter onto an NS may increase its mass above the maximum value allowed by its equation of state, inducing its collapse to a black hole (BH). Here we study this process for the first time, in three-dimensions, and in full general relativity. By considering three initial NS configurations, each with and without a surrounding disk (of mass ~7% M_{NS}), we investigate the effect of the accretion disk on the dynamics of the collapse and its imprint on both the gravitational wave (GW) and electromagnetic (EM) signals that can be emitted by these sources. We show in particular that, even if the GW signal is similar for the accretion induced collapse (AIC) and the collapse of an NS in vacuum (and detectable only for Galactic sources), the EM counterpart could allow us to discriminate between these two types of events. In fact, our simulations show that, while the collapse of an NS in vacuum leaves no appreciable baryonic matter outside the event horizon, an AIC is followed by a phase of rapid accretion of the surviving disk onto the newly formed BH. The post-collapse accretion rates, on the order of ~10^{-2} M_{sun} s^{-1}, make these events tantalizing candidates as engines of short gamma-ray bursts.

preprint2012arXiv

Signatures of photon-axion conversion in the thermal spectra and polarization of neutron stars

Conversion of photons into axions under the presence of a strong magnetic field can dim the radiation from magnetized astrophysical objects. Here we perform a detailed calculation aimed at quantifying the signatures of photon-axion conversion in the spectra, light curves, and polarization of neutron stars (NSs). We take into account the energy and angle-dependence of the conversion probability and the surface thermal emission from NSs. The latter is computed from magnetized atmosphere models that include the effect of photon polarization mode conversion due to vacuum polarization. The resulting spectral models, inclusive of the general-relativistic effects of gravitational redshift and light deflection, allow us to make realistic predictions for the effects of photon to axion conversion on observed NS spectra, light curves, and polarization signals. We identify unique signatures of the conversion, such as an increase of the effective area of a hot spot as it rotates away from the observer line of sight. For a star emitting from the entire surface, the conversion produces apparent radii that are either larger or smaller (depending on axion mass and coupling strength) than the limits set by NS equations of state. For an emission region that is observed phase-on, photon-axion conversion results in an inversion of the plane of polarization with respect to the no-conversion case. While the quantitative details of the features that we identify depend on NS properties (magnetic field strength, temperature) and axion parameters, the spectral and polarization signatures induced by photon-axion conversion are distinctive enough to make NSs very interesting and promising probes of axion physics.

preprint2012arXiv

The Effects of Irradiation on Hot Jovian Atmospheres: Heat Redistribution and Energy Dissipation

Hot Jupiters, due to the proximity to their parent stars, are subjected to a strong irradiating flux which governs their radiative and dynamical properties. We compute a suite of 3D circulation models with dual-band radiative transfer, exploring a relevant range of irradiation temperatures, both with and without temperature inversions. We find that, for irradiation temperatures T \lesssim 2000 K, heat redistribution is very efficient, producing comparable day- and night-side fluxes. For Tirr \approx 2200-2400 K, the redistribution starts to break down, resulting in a high day-night flux contrast. Our simulations indicate that the efficiency of redistribution is primarily governed by the ratio of advective to radiative timescales. Models with temperature inversions display a higher day-night contrast due to the deposition of starlight at higher altitudes, but we find this opacity-driven effect to be secondary compared to the effects of irradiation. The hotspot offset from the substellar point is large when insolation is weak and redistribution is efficient, and decreases as redistribution breaks down. The atmospheric flow can be potentially subjected to the Kelvin-Helmholtz instability (as indicated by the Richardson number) only in the uppermost layers, with a depth that penetrates down to pressures of a few millibars at most. Shocks penetrate deeper, down to several bars in the hottest model. Ohmic dissipation generally occurs down to deeper levels than shock dissipation (to tens of bars), but the penetration depth varies with the atmospheric opacity. The total dissipated Ohmic power increases steeply with the strength of the irradiating flux and the dissipation depth recedes into the atmosphere, favoring radius inflation in the most irradiated objects.

preprint2012arXiv

The influence of fallback discs on the spectral and timing properties of neutron stars

Fallback discs around neutron stars (NSs) are believed to be an expected outcome of supernova explosions. Here we investigate the consequences of such a common outcome for the timing and spectral properties of the associated NS population, using Monte Carlo population synthesis models. We find that the long-term torque exerted by the fallback disc can substantially influence the late-time period distribution, but with quantitative differences which depend on whether the initial spin distribution is dominated by slow or fast pulsars. For the latter, a single-peaked initial spin distribution becomes bimodal at later times. Timing ages tend to underestimate the real age of older pulsars, and overestimate the age of younger ones. Braking indices cluster in the range 1.5 <~ n <~ 3 for slow-born pulsars, and -0.5 <~ n <~ 5 for fast-born pulsars, with the younger objects found predominantly below n <~ 3. Large values of n, while not common, are possible, and associated with torque transitions in the NS+disc system. The 0.1-10 keV thermal luminosity of the NS+disc system is found to be generally dominated by the disc emission at early times, t <~ 10^3 yr, but this declines faster than the thermal surface emission of the NS. Depending on the initial parameters, there can be occasional periods in which some NSs switch from the propeller to the accretion phase, increasing their luminosity up to the Eddington limit for ~ 10^3-10^4 years.

preprint2012arXiv

Time Dependent Radiative Transfer for Multi-Level Atoms using Accelerated Lambda Iteration

We present a general formalism for computing self-consistent, numerical solutions to the time-dependent radiative transfer equation in low velocity, multi-level ions undergoing radiative interactions. Recent studies of time-dependent radiative transfer have focused on radiation hydrodynamic and magnetohydrodynamic effects without lines, or have solved time-independent equations for the radiation field simultaneously with time-dependent equations for the state of the medium. In this paper, we provide a fully time-dependent numerical solution to the radiative transfer and atomic rate equations for a medium irradiated by an external source of photons. We use Accelerated Lambda Iteration to achieve convergence of the radiation field and atomic states. We perform calculations for a three-level atomic model that illustrates important time-dependent effects. We demonstrate that our method provides an efficient, accurate solution to the time-dependent radiative transfer problem. Finally, we characterize astrophysical scenarios in which we expect our solutions to be important.

preprint2012arXiv

X-ray emission from high-redshift miniquasars: self-regulating the population of massive black holes through global warming

Observations of high-redshift quasars at z>6 imply that supermassive black holes (SMBHs) with masses over a billion solar masses were in place less than 1 Gyr after the Big Bang. If these SMBHs assembled from "seed" BHs left behind by the first stars, then they must have accreted gas at close to the Eddington limit during a large fraction (>50%) of the time. A generic problem with this scenario, however, is that the mass density in million-solar-mass SMBHs at z=6 already exceeds the locally observed SMBH mass density by several orders of magnitude; in order to avoid this overproduction, BH seed formation and growth must become significantly less efficient in less massive protogalaxies, while proceeding uninterrupted in the most massive galaxies that formed first. Using Monte-Carlo realizations of the merger and growth history of BHs, we show that X-rays from the earliest accreting BHs can provide such a feedback mechanism. Our calculations paint a self-consistent picture of black-hole-made climate change, in which the first miniquasars---among them the ancestors of the z>6 quasar SMBHs---globally warm the IGM and suppress the formation and growth of subsequent generations of BHs. We present two specific models with global miniquasar feedback that provide excellent agreement with recent estimates of the z=6 SMBH mass function. For each of these models, we estimate the rate of BH mergers at z>6 that could be detected by the proposed gravitational-wave observatory eLISA/NGO.

preprint2011arXiv

A unified model of the magnetar and radio pulsar bursting phenomenology

Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) are young neutron stars (NSs) characterized by high X-ray quiescent luminosities, outbursts, and, in the case of SGRs, sporadic giant flares. They are believed to be powered by ultra-strong magnetic fields (hence dubbed magnetars). The diversity of their observed behaviours is however not understood, and made even more puzzling by the discovery of magnetar-like bursts from "low-field" pulsars. Here we perform long-term 2D simulations that follow the evolution of magnetic stresses in the crust; these, together with recent calculations of the breaking stress of the neutron star crust, allow us to establish when starquakes occur. For the first time, we provide a quantitative estimate of the burst energetics, event rate, and location on the neutron star surface, which bear a direct relevance for the interpretation of the overall magnetar phenomenology. Typically, an "SGR-like" object tends to be more active than an "AXP-like" object or a "high-$B$ radio pulsar", but there is no fundamental separation among what constitutes the apparent different classes. Among the key elements that create the variety of observed phenomena, age is more important than a small variation in magnetic field strength. We find that outbursts can also be produced in old, lower-field pulsars (B ~ a few x 10^{12} G), but those events are much less frequent than in young, high-field magnetars.

preprint2011arXiv

Emission geometry, radiation pattern, and magnetic topology of the magnetar XTE J1810-197 in its quiescent state

The return to the quiescent state of the Anomalous X-ray pulsar XTE J1810-197 following its 2003 outburst represents a unique opportunity to probe the surface emission properties of a magnetar. The quiescent emission of XTE J1810-197 is composed of two thermal components, one arising from the whole star surface, and the other from a small warm spot on it. By modeling the magnitude and shape of the pulse profile in narrow spectral bands, we have been able to constrain the physical characteristics and geometrical parameters of the system: the two angles that the line of sight and the spin axis make with respect to the warm spot axis (ψ and ξ respectively), the angular size of the spot, and the overall surface temperature distribution. Our modeling accounts for the general relativistic effects of gravitational redshift and light bending near the stellar surface, and allows for local anisotropic emission. We found that the surface temperature distribution on the neutron star is consistent with the expectations of a dipole magnetic field configuration; the local radiation requires a pencil-beamed emission pattern, suggesting the presence of a magnetized atmosphere. For a typical value of the radius, R=13 km, the viewing parameters (symmetric for an interchange between ψ and ξ), range from ψ=ξ=38 deg to (ψ,ξ)=(52 deg, 29 deg). These angles are consistent with those obtained by modeling the AXP in outburst, with uncertainty contours reduced by a factor of 2.5.

preprint2011arXiv

Magnetars vs. high magnetic field pulsars: a theoretical interpretation of the apparent dichotomy

Highly magnetized neutron stars (NSs) are characterized by a bewildering range of astrophysical manifestations. Here, building on our simulations of the evolution of magnetic stresses in the NS crust and its ensuing fractures (Perna & Pons 2011), we explore in detail, for the middle-age and old NSs, the dependence of starquake frequency and energetics on the relative strength of the poloidal (B_p) and toroidal (B_tor) components. We find that, for B_p >~10^{14}G, since a strong crustal toroidal field B_tor B_p is quickly formed on a Hall timescale, the initial toroidal field needs to be B_tor >> B_p to have a clear influence on the outbursting behaviour. For initial fields B_p <~ 10^{14}G, it is very unlikely that a middle-age (t~10^5 years) NS shows any bursting activity. This study allows us to solve the apparent puzzle of how NSs with similar dipolar magnetic fields can behave in a remarkably different way: an outbursting 'magnetar' with a high X-ray luminosity, or a quiet, low-luminosity, "high-$B$" radio pulsar. As an example, we consider the specific cases of the magnetar 1E2259+586 and the radio pulsar PSRJ1814-1744, which at present have a similar dipolar field ~6x10^{13}G. We determine for each object an initial magnetic field configuration that reproduces the observed timing parameters at their current age. The same two configurations also account for the differences in quiescent X-ray luminosity and for the 'magnetar/outbursting' behaviour of 1E2259+586 but not of PSRJ1814-1744. We further use the theoretically predicted surface temperature distribution to compute the light-curve for these objects. In the case of 1E2259+586, for which data are available, our predicted temperature distribution gives rise to a pulse profile whose double-peaked nature and modulation level is consistent with the observations.

preprint2010arXiv

Flare-less long Gamma-ray Bursts and the properties of their massive star progenitors

While there is mounting evidence that long Gamma-Ray Bursts (GRBs) are associated with the collapse of massive stars, the detailed structure of their pre-supernova stage is still debatable. Particularly uncertain is the degree of mixing among shells of different composition, and hence the role of magnetic torques and convection in transporting angular momentum. Here we show that early-time afterglow observations with the Swift satellite place constraints on the allowed GRB pre-supernova models. In particular, they argue against pre-supernova models in which different elemental shells are unmixed. These types of models would produce energy injections into the GRB engine on timescales between several hundreds of seconds to a few hours. Flaring activity has {\em not} been observed in a large fraction of well-monitored long GRBs. Therefore, if the progenitors of long GRBs have common properties, then the lack of flares indicates that the massive stars which produce GRBs are mostly well mixed, as expected in low-metallicity, rapidly rotating massive stars.

preprint2010arXiv

Instabilities in the Gamma Ray Burst central engine. What makes the jet variable?

Both types of long and short gamma ray bursts involve a stage of a hyper-Eddington accretion of hot and dense plasma torus onto a newly born black hole. The prompt gamma ray emission originates in jets at some distance from this 'central engine' and in most events is rapidly variable, having a form of spikes and subpulses. This indicates at the variable nature of the engine itself, for which a plausible mechanism is an internal instability in the accreting flow. We solve numerically the structure and evolution of the neutrino-cooled torus. We take into account the detailed treatment of the microphysics in the nuclear equation of state that includes the neutrino trapping effect. The models are calculated for both Schwarzschild and Kerr black holes. We find that for sufficiently large accretion rates (> 10 Msun/s for non-rotating black hole, and >1 Msun/s for rotating black hole, depending on its spin), the inner regions of the disk become opaque, while the helium nuclei are being photodissociated. The sudden change of pressure in this region leads to the development of a viscous and thermal instability, and the neutrino pressure acts similarly to the radiation pressure in sub-Eddington disks. In the case of rapidly rotating black holes, the instability is enhanced and appears for much lower accretion rates. We also find the important and possibly further destabilizing role of the energy transfer from the rotating black hole to the torus via the magnetic coupling.

preprint2010arXiv

Magnetic Drag on Hot Jupiter Atmospheric Winds

Hot Jupiters, with atmospheric temperatures T ~ 1000 K, have residual thermal ionization levels sufficient for the interaction of the ions with the planetary magnetic field to result in a sizable magnetic drag on the (neutral) atmospheric winds. We evaluate the magnitude of magnetic drag in a representative three-dimensional atmospheric model of the hot Jupiter HD 209458b and find that it is a plausible mechanism to limit wind speeds in this class of atmospheres. Magnetic drag has a strong geometrical dependence, both meridionally and from the day to the night side (in the upper atmosphere), which could have interesting consequences for the atmospheric flow pattern. By extension, close-in eccentric planets with transiently heated atmospheres will experience time-variable levels of magnetic drag. A robust treatment of magnetic drag in circulation models for hot atmospheres may require iterated solutions to the magnetic induction and Saha equations as the hydrodynamical flow is evolved.

preprint2010arXiv

Ohmic Dissipation in the Atmospheres of Hot Jupiters

Hot Jupiter atmospheres exhibit fast, weakly-ionized winds. The interaction of these winds with the planetary magnetic field generates drag on the winds and leads to ohmic dissipation of the induced electric currents. We study the magnitude of ohmic dissipation in representative, three-dimensional atmospheric circulation models of the hot Jupiter HD 209458b. We find that ohmic dissipation can reach or exceed 1% of the stellar insolation power in the deepest atmospheric layers, in models with and without dragged winds. Such power, dissipated in the deep atmosphere, appears sufficient to slow down planetary contraction and explain the typically inflated radii of hot Jupiters. This atmospheric scenario does not require a top insulating layer or radial currents that penetrate deep in the planetary interior. Circulation in the deepest atmospheric layers may actually be driven by spatially non-uniform ohmic dissipation. A consistent treatment of magnetic drag and ohmic dissipation is required to further elucidate the consequences of magnetic effects for the atmospheres and the contracting interiors of hot Jupiters.

preprint2010arXiv

The rest-frame ultraviolet spectra of GRBs from massive rapidly-rotating stellar progenitors

The properties of a massive star prior to its final explosion are imprinted in the circumstellar medium (CSM) created by its wind and termination shock. We perform a detailed, comprehensive calculation of the time-variable and angle-dependent transmission spectra of an average-luminosity Gamma-Ray Burst (GRB) which explodes in the CSM structure produced by the collapse of a 20 Msun, rapidly rotating, Z=0.001 progenitor star. We study both the case in which metals are initially in the gaseous phase, as well as the situation in which they are heavily depleted into dust. We find that high-velocity lines from low-ionization states of silicon, carbon, and iron are initially present in the spectrum only if the metals are heavily depleted into dust prior to the GRB explosion. However, such lines disappear on timescales of a fraction of a second for a burst observed on-axis, and of a few seconds for a burst seen at high-latitude, making their observation virtually impossible. Rest-frame lines produced in the termination shock are instead clearly visible in all conditions. We conclude that time-resolved, early-time spectroscopy is not a promising way in which the properties of the GRB progenitor wind can be routinely studied. Previous detections of high velocity features in GRB UV spectra must have been due either due to a superposition of a physically unrelated absorber or to a progenitor star with very unusual properties.

preprint2009arXiv

Gravitational Lensing of Anisotropic Sources

In strong gravitational lensing, the multiple images we see correspond to light rays that leave the source in slightly different directions. If the source emission is anisotropic, the images may differ from conventional lensing predictions (which assume isotropy). To identify scales on which source anisotropy may be important, we study the angle delta between the light rays emerging from the source, for different lensing configurations. If the lens has a power law profile M propto R^gamma, the angle delta initially increases with lens redshift and then either diverges (for a steep profile gamma<1), remains constant (for an isothermal profile gamma=1), or vanishes (for a shallow profile gamma>1) as zl approaches zs. The scaling with lens mass is roughly delta propto M^(1/(2-gamma)). The results for an NFW profile are qualitatively similar to those for a shallow power law, with delta peaking at about half the redshift of the source (not half the distance). In practice, beaming could modify the statistics of beamed sources lensed by massive clusters: for an opening angle theta_jet, there is a probability as high as P ~ 0.02-0.07 (theta_jet/{0.5} deg)^-1 that one of the lensed images may be missed (for 2 < zs < 6). Differential absorption within Active Galactic Nuclei could modify the flux ratios of AGNs lensed by clusters; a sample of AGNs lensed by clusters could provide further constraints on the sizes of absorbing regions. Source anisotropy is not likely to be a significant effect in galaxy-scale strong lensing.

preprint2009arXiv

Soft X-ray Polarization in Thermal Magnetar Emission

Emission spectra from magnetars in the soft X-ray band likely contain a thermal component emerging directly from the neutron star surface. However, the lack of observed absorption-like features in quiescent spectra makes it difficult to directly constrain physical properties of the atmosphere. We argue that future X-ray polarization measurements represent a promising technique for directly constraining the magnetar magnetic field strength and geometry. We construct models of the observed polarization signal from a finite surface hotspot, using the latest NS atmosphere models for magnetic fields B = 4 x 10^13--5 x 10^14 G. Our calculations are strongly dependent on the NS magnetic field strength and geometry, and are more weakly dependent on the NS equation of state and atmosphere composition. We discuss how the complementary dependencies of phase-resolved spectroscopy and polarimetry might resolve degeneracies that currently hamper the determination of magnetar physical parameters using thermal models.

preprint2009arXiv

The first outburst of the new magnetar candidate SGR 0501+4516

We report here on the outburst onset and evolution of the new Soft Gamma Repeater SGR 0501+4516. We monitored the new SGR with XMM-Newton starting on 2008 August 23, one day after the source became burst-active, and continuing with 4 more observations, with the last one on 2008 September 30. Combining the data with the Swift-XRT and Suzaku data, we modelled the outburst decay over 160 days, and we found that the source flux decreased exponentially with a timescale of t_c=23.8 days. In the first XMM-Newton observation a large number of short X-ray bursts were observed, the rate of which decayed drastically in the following observations. We found large changes in the spectral and timing behavior of the source during the outburst, with softening emission as the flux decayed, and the non-thermal soft X-ray spectral component fading faster than the thermal one. Almost simultaneously to our XMM-Newton observations (on 2008 August 29 and September 2), we observed the source in the hard X-ray range with INTEGRAL, which clearly detected the source up to ~100keV in the first pointing, while giving only upper limits during the second pointing, discovering a variable hard X-ray component fading in less than 10 days after the bursting activation. We performed a phase-coherent X-ray timing analysis over about 160 days starting with the burst activation and found evidence of a strong second derivative period component (\ddot{P} = -1.6(4)x10^{-19} s/s^{-2}). Thanks to the phase-connection, we were able to study the the phase-resolved spectral evolution of SGR 0501+4516 in great detail. We also report on the ROSAT quiescent source data, taken back in 1992 when the source exhibits a flux ~80 times lower than that measured during the outburst, and a rather soft, thermal spectrum.

preprint2007arXiv

Dust Echoes from the Ambient Medium of Gamma-Ray Bursts

Long gamma-ray bursts (GRBs) are likely associated with the collapse of massive stars, which produce dust and are born in dusty environments. Absorption and scattering of ultraviolet/X-ray photons from the prompt, optical flash and afterglow emission of the GRB produce dust echoes. We perform time-dependent calculations of these echoes, accounting for the evolution of the dust grain distribution due to selective grain destruction by the GRB radiation, and for off-axis beaming. We explore cloud configurations of differing density and size -- the echo light curve and spectrum depend on the cloud radius, with larger clouds peaking at longer wavelengths. For a region ~ 3 pc in size with hydrogen density ~ 1000 per cubic centimeter, the echo spectrum peaks at ~ 3.6 microns and ~ 8.8 eV for thermal and scattered components, respectively. Dust echoes should be detectable with the Very Large Telescope up to z ~ 0.1, IRAC onboard the Spitzer Space Telescope up to z ~ 0.2, and NICMOS onboard the Hubble Space Telescope up to z ~ 0.3. Furthermore, the shape of the echo light curve allows one to infer: the jet opening angle; the inclination of the jet axis with respect to the line of sight; the size of the dust-emitting region. For sources with symmetric, bipolar jets, dust echoes exhibit two bumps in the light curve, making them easily distinguishable from the rebrightening due to an underlying supernova.

preprint2002arXiv

Neutrino trapping and accretion models for Gamma-Ray Bursts

Many models of Gamma Ray Bursts invoke a central engine consisting of a black hole of a few solar masses accreting matter from a disk at a rate of a fraction to a few solar masses per second. Popham et al. and Narayan et al. have shown that, for Mdot >~ 0.1 Msun/s, accretion proceeds via neutrino cooling and neutrinos can carry away a significant amount of energy from the inner regions of the disks. We improve on these calculations by including a simple prescription for neutrino transfer and neutrino opacities in such regions. We find that the flows become optically thick to neutrinos inside a radius R~6-40R_s for Mdot in the range of 0.1 -10 Msun/s, where R_s is the black hole Schwarzchild radius. Most of the neutrino emission comes from outside this region and, the neutrino luminosity stays roughly constant at a value L_ν ~ 10^{53} erg/s. We show that, for Mdot > 1 Msun/s, neutrinos are sufficiently trapped that energy advection becomes the dominant cooling mechanism in the flow. These results imply that neutrino annihilation in hyperaccreting black holes is an inefficient mechanism for liberating large amounts of energy. Extraction of rotational energy by magnetic processes remains the most viable mechanism.