Showing 1–50 of 170 results for author: Kazanas, D

Dual Role of Accretion Disk Winds as X-ray Obscurers and UV Line Absorbers in AGN

Authors: Keigo Fukumura, Missagh Mehdipour, Ehud Behar, Chris Shrader, Mauro Dadina, Demosthenes Kazanas, Stefano Marchesi, Francesco Tombesi

Abstract: X-ray obscuration of active galactic nuclei (AGNs) is considered in the context of ionized winds of stratified structure launched from accretion disks. We argue that a Compton-thick layer of a large-scale disk wind can obscure continuum X-rays and also lead to broad UV absorption such as in the blue wing of Civ; the former originates from the inner wind while the latter from the outer wind as a du… ▽ More X-ray obscuration of active galactic nuclei (AGNs) is considered in the context of ionized winds of stratified structure launched from accretion disks. We argue that a Compton-thick layer of a large-scale disk wind can obscure continuum X-rays and also lead to broad UV absorption such as in the blue wing of Civ; the former originates from the inner wind while the latter from the outer wind as a dual role. Motivated by a number of observational evidence showing strong AGN obscuration phenomena in Seyfert 1 AGNs such as NGC 5548, we demonstrate in this work, by utilizing a physically-motivated wind model coupled to post-process radiative transfer calculations, that an extended disk wind under certain physical conditions (e.g. morphology and density) could naturally cause a sufficient obscuration qualitatively consistent with UV/X-ray observations. Predicted UV/X-ray correlation is also presented as a consequence of variable spatial size of the wind in this scenario. △ Less

Submitted 17 April, 2024; v1 submitted 11 March, 2024; originally announced March 2024.

Comments: Comments welcome! (10 pages, 4 figures)

Journal ref: revised and under 2nd review for ApJ, 2024

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 43 pages JCAP style (excluding author list and references), 19 figures

Gravitational Waves from the Pulsar Magnetosphere

Authors: Ioannis Contopoulos, Demosthenes Kazanas, Demetrios B. Papadopoulos

Abstract: We investigate the generation of gravitational waves from the rotation of an orthogonal pulsar magnetosphere in flat space time. We calculate the first order metric perturbation due to the rotation of the non-axisymmetric distribution of electromagnetic energy density around the central star. We show that gravitational waves from a strong magnetic field pulsar right after its formation within a di… ▽ More We investigate the generation of gravitational waves from the rotation of an orthogonal pulsar magnetosphere in flat space time. We calculate the first order metric perturbation due to the rotation of the non-axisymmetric distribution of electromagnetic energy density around the central star. We show that gravitational waves from a strong magnetic field pulsar right after its formation within a distance of 1 kpc may be detectable with the new generation of gravitational wave detectors. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 9 pages, 1 figure, accepted for publication in Monthly Notices

Accreting Black Holes Skewing and Bending the Optical Emission from Massive Wolf-Rayet Companions -- A Case Study of IC10 X-1

Authors: Sayantan Bhattacharya, Dimitris M. Christodoulou, Andre-Nicolas Chene, Silas G. T. Laycock, Breanna A. Binder, Demosthenes Kazanas

Abstract: We present a statistical analysis of the He ii 4686 emission line in the spectra of the black hole and Wolf-Rayet (WR) star of the high-mass X-ray binary IC10 X-1. This line is visibly skewed, and the third moment (skewness) varies with the binary's orbital phase. We describe a new method of extracting such weak/faint features lying barely above a noisy continuum. Using the moments of these featur… ▽ More We present a statistical analysis of the He ii 4686 emission line in the spectra of the black hole and Wolf-Rayet (WR) star of the high-mass X-ray binary IC10 X-1. This line is visibly skewed, and the third moment (skewness) varies with the binary's orbital phase. We describe a new method of extracting such weak/faint features lying barely above a noisy continuum. Using the moments of these features, we have been able to decompose these skewed lines into two symmetric Gaussian profiles as a function of the orbital phase. The astrophysical implications of this decomposition are significant due to the complex nature of wind-accretion stream interactions in such binary systems. Previous studies have already shown a 0.25 phase lag in the radial velocity curve of the star and the X-ray eclipse, which indicates that the He ii emitters might be in the stellar wind, hence not tracing the star's orbital motion. Results from this work further suggest the existence of two separate emitting regions, one in the stellar wind in the shadow of the WR star, and another in the accretion stream that impacts the black hole's outer accretion disk; and the observed skewed He ii lines can be reproduced by superposition of the two corresponding time-dependent Gaussian emission profiles. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: 13 pages, 5 figures, accepted for publication in MNRAS Journal

The Gamma-Ray Pulsar Phenomenology in View of 3D Kinetic Global Magnetosphere Models

Authors: Constantinos Kalapotharakos, Zorawar Wadiasingh, Alice K. Harding, Demosthenes Kazanas

Abstract: We develop kinetic plasma models of pulsar magnetospheres with magnetic-field-line-dependent plasma injection that reveal the importance of various magnetosphere regions in regulating the gamma-ray emission. We set different particle injection rates for the so-called open, closed, and separatrix zones. Moderate particle injection rates in open and closed zones ensure a global field structure close… ▽ More We develop kinetic plasma models of pulsar magnetospheres with magnetic-field-line-dependent plasma injection that reveal the importance of various magnetosphere regions in regulating the gamma-ray emission. We set different particle injection rates for the so-called open, closed, and separatrix zones. Moderate particle injection rates in open and closed zones ensure a global field structure close to the force-free one, while the dissipation occurs mainly in and around the equatorial current sheet. The particles injected in the separatrix zone affect the particle populations that enter the equatorial current sheet region and, therefore, the corresponding accelerating electric fields, particle energies, the spectral cutoff energy, and gamma-ray efficiency. The separatrix zone models reproduce the recently discovered fundamental plane of gamma-ray pulsars consistent with curvature radiation emission, the gamma-ray light-curve shapes, and the radio-lag vs. peak-separation correlation reported in the Fermi second pulsar catalog. The model beaming factors indicate that the pulsar total gamma-ray luminosities listed in the Fermi catalogs are overestimations of the actual ones. We find that the radiation reaction limited regime starts ceasing to govern the high-energy emission for spin-down powers less than $10^{34}$ erg/s. Our results also indicate that toward high magnetic inclination angles, the "Y point" around the rotational equator migrates well inside the light cylinder sparking additional peaks in the gamma-ray pulse profiles. We find that an equivalent enhanced particle injection beyond the Y point strengthens these features making the model gamma-ray light curves inconsistent with those observed. △ Less

Submitted 27 July, 2023; v1 submitted 7 March, 2023; originally announced March 2023.

Comments: 28 pages, 22 figures, accepted for publication in ApJ

Poynting-Robertson effect on black-hole-driven winds

Authors: M. Marzi, F. Tombesi, A. Luminari, K. Fukumura, D. Kazanas

Abstract: Layers of ionized plasma, in the form of winds ejected from the accretion disk of Supermassive Black Holes (SMBHs) are frequently observed in Active Galactic Nuclei (AGNs). Winds with a velocity often exceeding $0.1c$ are called Ultra-Fast-Outflows (UFOs) and thanks to their high power they can play a key role in the co-evolution between the SMBH and the host galaxy. In order to construct a proper… ▽ More Layers of ionized plasma, in the form of winds ejected from the accretion disk of Supermassive Black Holes (SMBHs) are frequently observed in Active Galactic Nuclei (AGNs). Winds with a velocity often exceeding $0.1c$ are called Ultra-Fast-Outflows (UFOs) and thanks to their high power they can play a key role in the co-evolution between the SMBH and the host galaxy. In order to construct a proper model of the properties of these winds, it is necessary to consider special relativistic corrections due to their very high velocities. We present a derivation of the Poynting-Robertson effect (P-R effect) and apply it to the description of the dynamics of UFOs. The P-R effect is a special relativistic correction which breaks the isotropy of the radiation emitted by a moving particle funneling the radiation in the direction of motion. As a result of the conservation of the four-momentum, the emitting particles are subjected to a drag force and decelerate. We provide a derivation of the drag force caused by the P-R effect starting from general Lorentz transformations and assuming isotropic emission in the gas reference frame. Then, we derive the equations to easily implement this drag force in future simulations. Finally, we apply them in a toy model in which the gas particles move radially under the influence of the gravitation force, the radiation pressure and the drag due to the P-R effect. P-R effect plays an important role in determining the velocity profile of the wind. For a wind launched from $r_0=10r_s$ (where $r_S$ stands for the Schwarzschild radius), the asymptotic velocity reached by the wind is between $10$% and $24$% smaller than the one it would possess if we neglect the effect. This shows that the P-R effect should be taken into account when studying the dynamics of high-velocity, photoionized outflows in general. △ Less

Submitted 14 February, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: Accepted for publication on Astronomy & Astrophysics. 7 pages, 4 figures

Journal ref: A&A 670, A122 (2023)

Constraining f(R) models with cosmic chronometers and the HII galaxy Hubble diagram

Authors: Joseph Sultana, Manoj K. Yennapureddy, Fulvio Melia, Demosthenes Kazanas

Abstract: We consider several well-known f(R) cosmological models and constrain their parameters, namely the deviation parameter b and the cosmological parameters Ω_m and h. We first obtain analytical approximations for the Hubble rate H(z) and the luminosity distance d_L(z) in terms of these parameters, and then test these against the observational expansion rate derived from cosmic chronometers and the di… ▽ More We consider several well-known f(R) cosmological models and constrain their parameters, namely the deviation parameter b and the cosmological parameters Ω_m and h. We first obtain analytical approximations for the Hubble rate H(z) and the luminosity distance d_L(z) in terms of these parameters, and then test these against the observational expansion rate derived from cosmic chronometers and the distance modulus in the HII galaxy Hubble diagram, obtained in a model-independent way using Gaussian Processes (GP). We first optimize the models based solely on the cosmic chronometers and then repeat this process with a joint analysis using both the cosmic chronometers and HII galaxies. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 16 pages

Tell-Tale Spectral Signatures of MHD-driven Ultra-Fast Outflows in AGNs

Authors: Keigo Fukumura, Mauro Dadina, Gabriele Matzeu, Francesco Tombesi, Chris Shrader, Demosthenes Kazanas

Abstract: We aim to explore spectral signatures of the predicted multi-ion UFOs in the broadband X-ray spectra of active galactic nuclei (AGNs) by exploiting an accretion disk wind model in the context of a simple magnetohydrodynamic (MHD) framework. We are focused primarily on examining the spectral dependences on a number of key properties; (1) ionizing luminosity ratio $λ_{\rm ion}$, (2) line-of-sight wi… ▽ More We aim to explore spectral signatures of the predicted multi-ion UFOs in the broadband X-ray spectra of active galactic nuclei (AGNs) by exploiting an accretion disk wind model in the context of a simple magnetohydrodynamic (MHD) framework. We are focused primarily on examining the spectral dependences on a number of key properties; (1) ionizing luminosity ratio $λ_{\rm ion}$, (2) line-of-sight wind density slope $p$, (3) optical/UV-to-X-ray strength $α_{\rm OX}$, (4) inclination $θ$, (5) X-ray photon index $Γ$ and (6) wind density factor $f_D$. With an emphasis on radio-quiet Seyferts in sub-Eddington regime, multi-ion UFO spectra are systematically calculated as a function of these parameters to show that MHD-driven UFOs imprint a unique asymmetric absorption line profile with a pronounced blue tail structure on average. Such a characteristic line signature is generic to the simplified MHD disk-wind models presented in this work due to their specific kinematics and density structure. The properties of these absorption line profiles could be utilized as a diagnostics to distinguish between different wind driving mechanisms or even the specific values of a given MHD wind parameters. We also present high fidelity microcalorimeter simulations in anticipation of the upcoming {\it XRISM}/Resolve and {\it Athena}/X-IFU instruments to demonstrate that such a "tell-tale" sign may be immune to a spectral contamination by the presence of additional warm absorbers and partially covering gas. △ Less

Submitted 19 September, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: accepted to ApJ

Stabilizing Spherical Energy Shells with Angular Momentum in Gravitational Backgrounds

Authors: I. Antoniou, D. Kazanas, D. Papadopoulos, L. Perivolaropoulos

Abstract: Spherical energy shells in General Relativity tend to collapse due to gravitational effects and/or due to tension effects. Shell stabilization may be achieved by modifying the gravitational properties of the background spacetime. Thus, gravastars consist of stiff matter shells with an interior deSitter space and an exterior Schwarzshild spacetime whose attractive gravity balances the interior repu… ▽ More Spherical energy shells in General Relativity tend to collapse due to gravitational effects and/or due to tension effects. Shell stabilization may be achieved by modifying the gravitational properties of the background spacetime. Thus, gravastars consist of stiff matter shells with an interior deSitter space and an exterior Schwarzshild spacetime whose attractive gravity balances the interior repulsive gravity of the interior deSitter spacetime leading to a stable stiff matter shell. Similar stabilization effects may be achieved by considering rotating shells. Here we study the stability of slowly rotating fluid shells. We show that the angular velocity of the shell has stabilizing properties analogous to the repulsive deSitter gravity of the interior of a gravastar. We thus use the Israel junction conditions and the fluid equation of state of the rotating shell to construct the dynamical equations that determine the evolution of the rotating shell radius. These dynamical equations depend on the parameters of the background spacetime and on the angular velocity of the shell. Assuming a rotating interior and a Schwarzschild exterior spacetime we show that the angular velocity of the shell has interesting stabilizing properties on the evolution of its radius R. Thus rotating matter (or vacuum) shells can imitate black holes while avoiding the presence of a singularity and without the presence of an interior deSitter space. △ Less

Submitted 29 April, 2022; originally announced April 2022.

Comments: 19 pages, 2 figures. Accepted for publication in International Journal of Modern Physics D (to appear)

The Fundamental Plane Relation for Gamma-Ray Pulsars Implied by 4FGL

Authors: Constantinos Kalapotharakos, Zorawar Wadiasingh, Alice K. Harding, Demosthenes Kazanas

Abstract: We explore the validity of the recently reported fundamental plane (FP) relation of gamma-ray pulsars using 190 pulsars included in the latest 4FGL-DR3 catalog. This sample number is more than twice as large as that of the original study. The FP relation incorporates 4 parameters, i.e., the spin-down power, $\dot{\mathcal{E}}$, the surface magnetic field, $B_{\star}$, the total gamma-ray luminosit… ▽ More We explore the validity of the recently reported fundamental plane (FP) relation of gamma-ray pulsars using 190 pulsars included in the latest 4FGL-DR3 catalog. This sample number is more than twice as large as that of the original study. The FP relation incorporates 4 parameters, i.e., the spin-down power, $\dot{\mathcal{E}}$, the surface magnetic field, $B_{\star}$, the total gamma-ray luminosity, $L_γ$, and a spectral cutoff energy, $ε_{\rm cut}$. The derivation of $ε_{\rm cut}$ is the most intriguing one because $ε_{\rm cut}$ depends on the proper interpretation of the available phase-averaged spectra. We construct synthetic phase-averaged spectra, guided by the few existing phase-resolved ones, to find that the best fit cutoff energy, $ε_{\rm c1}$, corresponding to a purely exponential cutoff (plus a power law) spectral form, is the parameter that optimally probes the maximum cutoff energy of the emission that originates from the core of the dissipative region, i.e., the equatorial current sheet. Computing this parameter for the 190 4FGL pulsars, we find that the resulting FP relation, i.e. the gamma-ray luminosity in terms of the other observables, reads $L_γ=10^{14.3\pm 1.3}(ε_{\rm c1}/{\rm MeV})^{1.39\pm0.17}(B_{\star}/{\rm G})^{0.12\pm 0.03}(\dot{\mathcal{E}}/{\rm erg\;s^{-1}})^{0.39\pm 0.05}{\rm ~erg\;s^{-1}}$; this is in good agreement with both the empirical relation reported by Kalapotharakos et al. (2019) and the theoretically predicted relation for curvature radiation. Finally, we revisit the radiation reaction limited condition, to find it is a sufficient but not necessary condition for the theoretical derivation of the FP relation. However, the assumption of the radiation reaction limited acceleration reveals the underlying accelerating electric field component and its scaling with $\dot{\mathcal{E}}$. △ Less

Submitted 15 June, 2022; v1 submitted 24 March, 2022; originally announced March 2022.

Comments: 20 pages, 14 figures, 3 tables, minor revisions and a typo correction in the most important equation of the paper, accepted for publication in ApJ

Gravitational Waves from GRB Core Spindown

Authors: Ioannis Contopoulos, Achilles Strantzalis, Dimitrios Papadopoulos, Demosthenes Kazanas

Abstract: We investigate long Gamma-Ray Bursts (GRB) which manifest a sharp linear rise followed by an exponential decay in their gamma-ray prompt emission observed with the BAT instrument on board the Swift satellite. We offer a simple electrodynamic model that may account for these particular characteristics. We associate the sharp rise with the winding of the magnetic field by the fast rotating core that… ▽ More We investigate long Gamma-Ray Bursts (GRB) which manifest a sharp linear rise followed by an exponential decay in their gamma-ray prompt emission observed with the BAT instrument on board the Swift satellite. We offer a simple electrodynamic model that may account for these particular characteristics. We associate the sharp rise with the winding of the magnetic field by the fast rotating core that formed in the interior of the stellar precursor. We also associate the subsequent exponential decay with the electromagnetic spin-down of the core following the release of the electromagnetic jet from the stellar interior. Any non-axisymmetric distortion in the rotating core will generate gravitational waves with exponentially decreasing frequency, a so-called "down-chirp". We obtain a detailed estimate of the gravitational wave profile if the distortion of spacetime is due to the winding of a non-axisymmetric component of the magnetic field during that particular phase of the burst. We offer 7 particular time intervals during which one may look into LIGO archival data for the presence of our particular predicted waveforms in order to test our interpretation. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 7 pages, 6 figures, accepted for publication in Monthly Notices

A two-zone emission model for Blazars and the role of Accretion Disk MHD winds

Authors: Stela S. Boula, Apostolos Mastichiadis, Demosthenes Kazanas

Abstract: Blazars are a sub-category of radio-loud active galactic nuclei with relativistic jets pointing towards the observer. They exhibit non-thermal variable emission, which practically extends over the whole electromagnetic spectrum. Despite the plethora of multi-wavelength observations, the origin of the emission in blazar jets remains an open question. In this work, we construct a two-zone leptonic m… ▽ More Blazars are a sub-category of radio-loud active galactic nuclei with relativistic jets pointing towards the observer. They exhibit non-thermal variable emission, which practically extends over the whole electromagnetic spectrum. Despite the plethora of multi-wavelength observations, the origin of the emission in blazar jets remains an open question. In this work, we construct a two-zone leptonic model: particles accelerate in a small region and lose energy through synchrotron radiation and inverse Compton Scattering. Consequently, the relativistic electrons escape to a larger area where the ambient photon field, which is related to Accretion Disk MHD Winds, could play a central role in the gamma-ray emission. This model explains the Blazar Sequence and the broader properties of blazars, as determined by Fermi observations, by varying only one parameter, the mass accretion rate onto the central black hole. Flat Spectrum Radio Quasars have a strong ambient photon field and their gamma-ray emission is dominated by the more extensive zone, while in the case of BL Lac objects, the negligible ambient photons make the smaller, i.e. acceleration, zone dominant. △ Less

Submitted 5 August, 2021; originally announced August 2021.

Comments: 8 pages, 5 figures, 37th International Cosmic Ray Conference (ICRC 2021), 2021 Online, Berlin, Germany

The infrared echo of SN2010jl and its implications for shock breakout characteristics

Authors: Eli Dwek, Arkaprabha Sarangi, Richard G. Arendt, Timothy Kallman, Demos Kazanas, Ori D. Fox

Abstract: SN 2010jl is a Type IIn core collapse supernova whose radiative output is powered by the interaction of the SN shock wave with its surrounding dense circumstellar medium (CSM). After day ~60, its light curve developed a NIR excess emission from dust. This excess could be a thermal IR echo from pre-existing CSM dust, or emission from newly-formed dust either in the cooling postshock region of the C… ▽ More SN 2010jl is a Type IIn core collapse supernova whose radiative output is powered by the interaction of the SN shock wave with its surrounding dense circumstellar medium (CSM). After day ~60, its light curve developed a NIR excess emission from dust. This excess could be a thermal IR echo from pre-existing CSM dust, or emission from newly-formed dust either in the cooling postshock region of the CSM, or in the cooling SN ejecta. Recent analysis has shown that dust formation in the CSM can commence only after day ~380, and has also ruled out newly-formed ejecta dust as the source of the NIR emission. The early (< 380 d) NIR emission can therefore only be attributed to an IR echo. The H-K color temperature of the echo is about 1250 K. The best fitting model requires the presence of about 1.6e-4 Msun of amorphous carbon dust at a distance of 2.2e16 cm from the explosion. The CSM-powered luminosity is preceded by an intense burst of hard radiation generated by the breakout of the SN shock through the stellar surface. The peak burst luminosity seen by the CSM dust is significantly reduced by Thomson scattering in the CSM, but still has the potential of evaporating the dust needed to produce the echo. We show that the survival of the echo-producing dust provides important constraints on the intensity, effective temperature, and duration of the burst. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: Accepted for publication in the ApJ, 20 pages, 1 table, and 17 figures

Gravitational analog of the canonical acoustic black hole in Einstein-scalar-Gauss-Bonnet theory

Authors: Pedro Cañate, Joseph Sultana, Demosthenes Kazanas

Abstract: In this work, in the context of modified gravity, a curved spacetime analogous to the "canonical acoustic black hole" is constructed. The source is a self-interacting scalar field which is non-minimally coupled to gravity through the Gauss-Bonnet invariant. The scalar-Gauss-Bonnet coupling function is characterized by three positive parameters: $σ$ with units of $(length)$, $μ$ with units of… ▽ More In this work, in the context of modified gravity, a curved spacetime analogous to the "canonical acoustic black hole" is constructed. The source is a self-interacting scalar field which is non-minimally coupled to gravity through the Gauss-Bonnet invariant. The scalar-Gauss-Bonnet coupling function is characterized by three positive parameters: $σ$ with units of $(length)$, $μ$ with units of $(length)^{4}$, and a dimensionless parameter $s$, thus defining a three-parameter model for which the line element of canonical acoustic black hole is a solution. The spacetime is equipped with spherical and static symmetry and has a single horizon determined in Schwarzschild coordinates by the region $r=μ^{1/4}$. The solution admits a photon sphere at $r=(3μ)^{1/4}$, and it is shown that in the region $(3μ)^{1/4}\leq r<\infty$ the scalar field satisfies the null, weak, and strong energy conditions. Nonetheless, the model with $s=1$ has major physical relevance since for this case the scalar field is well defined in the entire region $r\geqμ^{1/4}$, while for $s\neq1$ the scalar field blows up on the horizon. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: 12 pages, 7 figures

Modeling Magnetic Disk-Wind State Transitions in Black Hole X-ray Binaries

Authors: Keigo Fukumura, Demosthenes Kazanas, Chris Shrader, Francesco Tombesi, Constantinos Kalapotharakos, Ehud Behar

Abstract: We analyze three prototypical black hole (BH) X-ray binaries (XRBs), \4u1630, \gro1655\ and \h1743, in an effort to systematically understand the intrinsic state transition of the observed accretion-disk winds between \windon\ and \windoff\ states by utilizing state-of-the-art {\it Chandra}/HETGS archival data from multi-epoch observations. We apply our magnetically-driven wind models in the conte… ▽ More We analyze three prototypical black hole (BH) X-ray binaries (XRBs), \4u1630, \gro1655\ and \h1743, in an effort to systematically understand the intrinsic state transition of the observed accretion-disk winds between \windon\ and \windoff\ states by utilizing state-of-the-art {\it Chandra}/HETGS archival data from multi-epoch observations. We apply our magnetically-driven wind models in the context of magnetohydrodynamic (MHD) calculations to constrain their (1) global density slope ($p$), (2) their density ($n_{17}$) at the foot point of the innermost launching radius and (3) the abundances of heavier elements ($A_{\rm Fe,S,Si}$). Incorporating the MHD winds into {\tt xstar} photoionization calculations in a self-consistent manner, we create a library of synthetic absorption spectra given the observed X-ray continua. Our analysis clearly indicates a characteristic bi-modal transition of multi-ion X-ray winds; i.e. the wind density gradient is found to steepen (from $p \sim 1.2-1.4$ to $\sim 1.4-1.5$) while its density normalization declines as the source transitions from \windon\ to \windoff\ state. The model implies that the ionized wind {\it remains physically present} even in \windoff\ state, despite its absent appearance in the observed spectra. Super-solar abundances for heavier elements are also favored. Our global multi-ion wind models, taking into account soft X-ray ions as well as Fe K absorbers, show that the internal wind condition plays an important role in wind transitions besides photoionization changes. % Simulated {\it XRISM}/Resolve and {\it Athena}/X-IFU spectra are presented to demonstrate a high fidelity of the multi-ion wind model for better understanding of these powerful ionized winds in the coming decades. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: 47 pages, 19 figures, accepted to ApJ

On the observed clustering of major bodies in solar and extrasolar subsystems

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: Major (exo)planetary and satellite bodies seem to concentrate at intermediate areas of the radial distributions of all the objects present in each (sub)system. We prove rigorously that the secular evolution of (exo)planets and satellites necessarily results in the observed intermediate accumulation of the massive objects in all such subsystems. We quantify a "middle" as the mean of mean motions (o… ▽ More Major (exo)planetary and satellite bodies seem to concentrate at intermediate areas of the radial distributions of all the objects present in each (sub)system. We prove rigorously that the secular evolution of (exo)planets and satellites necessarily results in the observed intermediate accumulation of the massive objects in all such subsystems. We quantify a "middle" as the mean of mean motions (orbital angular velocities) of three or more massive objects involved. Orbital evolution is expected to be halted or severely diminished when the survivors settle near mean-motion resonances and substantial angular-momentum transfer between bodies ceases to occur (gravitational Landau damping). The dynamics is opposite in direction to what has been theorized for viscous and magnetized accretion disks in which gas spreads out and away from either side of any conceivable intermediate area. The results are bound to change the way we think about planet and moon formation and evolution. △ Less

Submitted 20 January, 2021; v1 submitted 6 January, 2021; originally announced January 2021.

Comments: Submitted to RAA

A variable magnetic disc wind in the black hole X-ray binary GRS 1915+105?

Authors: Ajay Ratheesh, Francesco Tombesi, Kiego Fukumura, Paolo Soffitta, Enrico Costa, Demosthenes Kazanas

Abstract: GRS 1915+105 being one of the brightest transient black hole binary (BHB) in the X-rays, offers a unique test-bed for the study of the connection between accretion and ejection mechanisms in BHBs. In particular, this source can be used to study the accretion disc wind and how it depends on the state changes in BHBs. Our aim is to investigate the origin and geometry of the accretion disc wind in GR… ▽ More GRS 1915+105 being one of the brightest transient black hole binary (BHB) in the X-rays, offers a unique test-bed for the study of the connection between accretion and ejection mechanisms in BHBs. In particular, this source can be used to study the accretion disc wind and how it depends on the state changes in BHBs. Our aim is to investigate the origin and geometry of the accretion disc wind in GRS 1915+105. We analysed the spectra of GRS 1915+105 in the soft $φ$ and hard $χ$ classes, using the high resolution spectroscopy offered by Chandra HETGS. In the soft state, we find a series of wind absorption lines that follow a non linear dependence of velocity width, velocity shift and equivalent width with respect to ionisation, indicating a multiple component or stratified outflow. In the hard state we find only a faint Fe XXVI absorption line. We model the absorption lines in both the states using a dedicated MHD wind model to investigate a magnetic origin of the wind and to probe the cause of variability in the observed lines flux between the two states. The MHD disc wind model provides a good fit for both states, indicating the possibility of a magnetic origin of the wind. The multiple ionisation components of the wind are well characterised as a stratification of the same magnetic outflow. We find that the observed variability in the lines flux between soft and hard states cannot be explained by photo-ionisation alone but it is most likely due to a large (three orders of magnitude) increase in the wind density. We find the mass outflow rate of the wind to be comparable to the accretion rate, suggesting a intimate link between accretion and ejection processes that lead to state changes in BHBs. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: 11 pages, 9 figures. Shortened abstract. Accepted for publication in A&A

Journal ref: A&A 646, A154 (2021)

The Multipolar Magnetic Field of Millisecond Pulsar PSR J0030+0451

Authors: Constantinos Kalapotharakos, Zorawar Wadiasingh, Alice K. Harding, Demosthenes Kazanas

Abstract: Modeling of the NICER X-ray waveform of the pulsar PSR J0030+0451, aimed to constrain the neutron star mass and radius, has inferred surface hot-spots (the magnetic polar caps) that imply significantly non-dipolar magnetic fields. To this end, we investigate magnetic field configurations that comprise offset dipole plus quadrupole components using static vacuum field and force-free global magnetos… ▽ More Modeling of the NICER X-ray waveform of the pulsar PSR J0030+0451, aimed to constrain the neutron star mass and radius, has inferred surface hot-spots (the magnetic polar caps) that imply significantly non-dipolar magnetic fields. To this end, we investigate magnetic field configurations that comprise offset dipole plus quadrupole components using static vacuum field and force-free global magnetosphere models. Taking into account the compactness and observer angle values provided by Miller et al. (2019) and Riley et al. (2019), we compute geodesics from the observer plane to the polar caps to compute the resulting X-ray light curve. We explore, through Markov chain Monte Carlo techniques, the detailed magnetic field configurations that can reproduce the observed X-ray light curve and have discovered degeneracies, i.e., diverse field configurations, which can provide sufficient descriptions to the NICER X-ray waveforms. Having obtained the force-free field structures, we then compute the corresponding synchronous gamma-ray light curves following Kalapotharakos et al. (2014) these we compare to those obtained by Fermi-LAT, to provide models consistent with both the X-ray and the gamma-ray data, thereby restricting further the multipole field parameters. An essential aspect of this approach is the proper computation of the relative phase between the synchronous X- and gamma-ray light curves. We conclude with a discussion of the broader implications of our study. △ Less

Submitted 3 December, 2020; v1 submitted 17 September, 2020; originally announced September 2020.

Comments: 21 pages, 11 figures, accepted for publication in ApJ (minor update)

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

Authors: The Cherenkov Telescope Array Consortium, :, A. Acharyya, R. Adam, C. Adams, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, L. Amati, G. Ambrosi, E. O. Angüner, L. A. Antonelli, C. Aramo, A. Araudo, T. Armstrong, F. Arqueros, K. Asano, Y. Ascasíbar, M. Ashley, C. Balazs, O. Ballester , et al. (427 additional authors not shown)

Abstract: We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models giv… ▽ More We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]" △ Less

Submitted 30 January, 2021; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 68 pages (including references) and 26 figures; text identical to the version published in JCAP

Journal ref: JCAP01(2021)057

On the importance of special relativistic effects in modelling ultra-fast outflows

Authors: A. Luminari, F. Tombesi, E. Piconcelli, F. Nicastro, K. Fukumura, D. Kazanas, F. Fiore, L. Zappacosta

Abstract: Outflows are observed in a variety of astrophysical sources. Remarkably, ultra-fast ($v\geq 0.1c$), outflows in the UV and X-ray bands are often seen in AGNs. Depending on their energy and mass outflow rate, respectively $\dot{E}_{out}, \dot{M}_{out}$, such outflows may play a key role in regulating the AGN-host galaxy co-evolution process through cosmic time. It is therefore crucial to provide ac… ▽ More Outflows are observed in a variety of astrophysical sources. Remarkably, ultra-fast ($v\geq 0.1c$), outflows in the UV and X-ray bands are often seen in AGNs. Depending on their energy and mass outflow rate, respectively $\dot{E}_{out}, \dot{M}_{out}$, such outflows may play a key role in regulating the AGN-host galaxy co-evolution process through cosmic time. It is therefore crucial to provide accurate estimates of the wind properties. Here, we concentrate on special relativistic effects concerning the interaction of light with matter moving at relativistic speed relatively to the source of radiation. Our aim is to assess the impact of these effects on the observed properties of the outflows and implement a relativistic correction in the existing spectral modelling routines. We define a simple procedure to incorporate relativistic effects in radiative transfer codes. Following this procedure, we run a series of simulations to explore the impact of these effects on the simulated spectra, for different $v$ and column densities of the outflow. The observed optical depth is usually considered a proxy for the wind $N_H$, independently on its velocity. However, our simulations show that the observed optical depth of an outflow with a given column density $N_H$ decreases rapidly as the velocity of the wind approaches relativistic values. This, in turn, implies that when estimating $N_H$ from the optical depth, it is necessary to include a velocity-dependent correction, already for moderate velocities (e.g. $v \geq 0.05c$). This correction linearly propagates to the derived $\dot{M}_{out}, \dot{E}_{out}$. As an example of these effects, we calculate the relativistically corrected values of $\dot{M}_{out}$ and $\dot{E}_{out}$ for a sample of $\sim 30$ Ultra-Fast Outflows taken from the literature, and find correction factors of $20-120 \%$ within the observed range of outflowing velocities. △ Less

Submitted 1 December, 2019; originally announced December 2019.

Comments: 6 pages, 5 figures. Accepted for publication in A&A

Journal ref: A&A 633, A55 (2020)

Generation of massive stellar black holes by rapid gas accretion in primordial dense clusters

Authors: Zacharias Roupas, Demosthenes Kazanas

Abstract: Supernova theory suggests that black holes of a stellar origin cannot attain masses in the range of 50-135 solar masses in isolation. We argue here that this mass gap is filled in by black holes that grow by gas accretion in dense stellar clusters, such as protoglobular clusters. The accretion proceeds rapidly, during the first 10 megayears of the cluster life, before the remnant gas is depleted.… ▽ More Supernova theory suggests that black holes of a stellar origin cannot attain masses in the range of 50-135 solar masses in isolation. We argue here that this mass gap is filled in by black holes that grow by gas accretion in dense stellar clusters, such as protoglobular clusters. The accretion proceeds rapidly, during the first 10 megayears of the cluster life, before the remnant gas is depleted. We predict that binaries of black holes within the mass gap can be observed by LIGO. △ Less

Submitted 23 December, 2019; v1 submitted 10 November, 2019; originally announced November 2019.

Comments: Minor typos corrected in final version

Journal ref: A&A 632, L8 (2019)

The Fast Radio Burst Luminosity Function and Death Line in the Low-Twist Magnetar Model

Authors: Zorawar Wadiasingh, Paz Beniamini, Andrey Timokhin, Matthew G. Baring, Alexander J. van der Horst, Alice K. Harding, Demosthenes Kazanas

Abstract: We explore the burst energy distribution of fast radio bursts (FRBs) in the low-twist magnetar model of Wadiasingh and Timokhin (2019). Motivated by the power-law fluence distributions of FRB 121102, we propose an elementary model for the FRB luminosity function of individual repeaters with an inversion protocol which directly relates the power-law distribution index of magnetar short burst fluenc… ▽ More We explore the burst energy distribution of fast radio bursts (FRBs) in the low-twist magnetar model of Wadiasingh and Timokhin (2019). Motivated by the power-law fluence distributions of FRB 121102, we propose an elementary model for the FRB luminosity function of individual repeaters with an inversion protocol which directly relates the power-law distribution index of magnetar short burst fluences to that for FRBs. The protocol indicates the FRB energy scales virtually linearly with crust/field dislocation amplitude, if magnetar short bursts prevail in the magnetoelastic regime. Charge starvation in the magnetosphere during bursts (required in WT19) for individual repeaters implies the predicted burst fluence distribution is narrow, $\lesssim 3$ decades for yielding strains and oscillation frequencies feasible in magnetar crusts. Requiring magnetic confinement and charge starvation, we obtain a death line for FRBs which segregates magnetars from the normal pulsar population, suggesting only the former will host recurrent FRBs. We convolve the burst energy distribution for individual magnetars to define the distribution of luminosities in evolved magnetar populations. The broken power-law luminosity function's low energy character depends on the population model, while the high energy index traces that of individual repeaters. Independent of the evolved population, the broken power-law isotropic-equivalent energy/luminosity function peaks at $\sim10^{37}-10^{40}$ erg with a low-energy cutoff at $\sim 10^{37}$ erg. Lastly, we consider the local fluence distribution of FRBs, and find that it can constrain the subset of FRB-producing magnetar progenitors. Our model suggests that improvements in sensitivity may reveal flattening of the global FRB fluence distribution and saturation in FRB rates. △ Less

Submitted 17 January, 2020; v1 submitted 15 October, 2019; originally announced October 2019.

Comments: 23 pages, 3 figures, accepted to The Astrophysical Journal

Dust formation in AGN winds

Authors: Arkaprabha Sarangi, Eli Dwek, Demos Kazanas

Abstract: Infrared observations of active galactic nucleus (AGN) reveal emission from the putative dusty circumnuclear 'torus' invoked by AGN unification, that is heated up by radiation from the central accreting black hole (BH). The strong 9.7 and 18 micron silicate features observed in the AGN spectra both in emission and absorption, further indicate the presence of such dusty environments. We present det… ▽ More Infrared observations of active galactic nucleus (AGN) reveal emission from the putative dusty circumnuclear 'torus' invoked by AGN unification, that is heated up by radiation from the central accreting black hole (BH). The strong 9.7 and 18 micron silicate features observed in the AGN spectra both in emission and absorption, further indicate the presence of such dusty environments. We present detailed calculations of the chemistry of silicate dust formation in AGN accretion disk winds. The winds considered herein are magnetohydrodynamic (MHD) winds driven off the entire accretion disk domain that extends from the BH vicinity to the radius of BH influence, of order of 1 to 100 pc depending on the mass of the resident BH. Our results indicate that these winds provide conditions conducive to the formation of significant amounts of dust, especially for objects accreting close to their Eddington limit, making AGN a significant source of dust in the universe, especially for luminous quasars. Our models justify the importance of a r to the power -1 density law in the winds for efficient formation and survival of dust grains. The dust production rate scales linearly with the mass of the central BH and varies as a power law of index between 2 to 2.5 with the dimensionless mass accretion rate. The resultant distribution of the dense dusty gas resembles a toroidal shape, with high column density and optical depths along the equatorial viewing angles, in agreement with the AGN unification picture. △ Less

Submitted 23 September, 2019; originally announced September 2019.

Comments: 22 pages, 17 figures (some have subfigures), Accepted for publication in The Astrophysical Journal

Meta-analysis of Electron Cyclotron Resonance Absorption Features Detected in High-Mass X-ray Binaries

Authors: Dimitris M. Christodoulou, Silas G. T. Laycock, Demosthenes Kazanas

Abstract: Using recent compilations of detailed X-ray observations and spectral models of exceptional quality, we record the electron cyclotron resonance absorption (ECRA) features that have been detected in 45 pulsating high-mass X-ray binaries (HMXBs) and ultraluminous X-ray (ULX) sources harboring neutron stars, although seven of these detections are still questionable and another 21 are single and/or no… ▽ More Using recent compilations of detailed X-ray observations and spectral models of exceptional quality, we record the electron cyclotron resonance absorption (ECRA) features that have been detected in 45 pulsating high-mass X-ray binaries (HMXBs) and ultraluminous X-ray (ULX) sources harboring neutron stars, although seven of these detections are still questionable and another 21 are single and/or not independently confirmed. From the comprehensive catalogs of Jaisawal \& Naik and Staubert et al. and from several additional recent observations, we produce two lists of HMXB ECRA sources: a list of 17 sources in which multiple ECRA lines or single very low-energy lines are seen, in which we can reasonably assume that the lowest energy reveals the fundamental cyclotron level for each source; and a `contaminated' list of 38 sources including the 21 detections of single ECRA lines that may (not) be higher-level harmonics. Both lists confirm a previous result that we have obtained independently by modeling the propeller lines of Magellanic HMXB pulsars: the surface dipolar magnetic fields $B_*$ of HMXB neutron stars are segregated around five distinct values with $B_* = 0.28\pm 0.08, 0.55\pm 0.11, 1.3\pm 0.37, 3.0\pm 0.68$, and $7.9\pm 3.1$, in units of TG. An explanation of this phenomenon is currently lacking. We have found no correlation between these $B_*$ values and the corresponding observed spin periods, spin period derivatives, orbital periods, maximum X-ray luminosities, neutron star masses, or companion star masses. △ Less

Submitted 16 May, 2019; v1 submitted 13 May, 2019; originally announced May 2019.

Comments: Submitted to RAA. This version incorporates the referee's comments and matches the published version

A Fundamental Plane for Gamma-Ray Pulsars

Authors: Constantinos Kalapotharakos, Alice K. Harding, Demosthenes Kazanas, Zorawar Wadiasingh

Abstract: We show that the $γ$-ray pulsar observables, i.e., their total $γ$-ray luminosity, $L_γ$, spectral cut-off energy, $ε_{\rm cut}$, stellar surface magnetic field, $B_{\star}$, and spin-down power $\dot{\mathcal{E}}$, obey a relation of the form $L_γ=f(ε_{\rm cut},B_{\star},\dot{\mathcal{E}})$, which represents a 3D plane in their 4D log-space. Fitting the data of 88 pulsars of the second Fermi puls… ▽ More We show that the $γ$-ray pulsar observables, i.e., their total $γ$-ray luminosity, $L_γ$, spectral cut-off energy, $ε_{\rm cut}$, stellar surface magnetic field, $B_{\star}$, and spin-down power $\dot{\mathcal{E}}$, obey a relation of the form $L_γ=f(ε_{\rm cut},B_{\star},\dot{\mathcal{E}})$, which represents a 3D plane in their 4D log-space. Fitting the data of 88 pulsars of the second Fermi pulsar catalog, we show this relation to be $L_γ\propto ε_{\rm cut}^{1.18\pm 0.24}B_{\star}^{0.17\pm 0.05}\dot{\mathcal{E}}^{0.41\pm 0.08}$, a pulsar fundamental plane (FP). We show that the observed FP is remarkably close to the theoretical relation $L_γ\propto ε_{\rm cut}^{4/3}B_{\star}^{1/6}\dot{\mathcal{E}}^{5/12}$ obtained assuming that the pulsar $γ$-ray emission is due to curvature radiation by particles accelerated at the pulsar equatorial current sheet just outside the light cylinder. Interestingly, the FP seems incompatible with emission by synchrotron radiation. The corresponding scatter about the FP is $\sim 0.35$dex and can only partly be explained by the observational errors while the rest is probably due to the variation of the inclination and observer angles. We predict also that $ε_{\rm cut}\propto \dot{\mathcal{E}}^{7/16}$ toward low $\dot{\mathcal{E}}$ for both young and millisecond pulsars implying that the observed death-line of $γ$-ray pulsars is due to $ε_{\rm cut}$ dropping below the Fermi-band. Our results provide a comprehensive interpretation of the observations of $γ$-ray pulsars, setting requirement for successful theoretical modeling. △ Less

Submitted 24 August, 2019; v1 submitted 3 April, 2019; originally announced April 2019.

Comments: 8 pages, 6 figures, accepted for publication in ApJL

Astro2020 Science White Paper: Do Supermassive Black Hole Winds Impact Galaxy Evolution?

Authors: F. Tombesi, M. Cappi, F. Carrera, G. Chartas, K. Fukumura, M. Guainazzi, D. Kazanas, G. Kriss, D. Proga, T. J. Turner, Y. Ueda, S. Veilleux, M. Brusa, M. Gaspari

Abstract: Powerful winds driven by supermassive black holes (SMBHs) are likely the main mechanism through which SMBHs regulate their own growth and influence the host galaxy evolution. However, their origin and their capability to impact the large-scale environment are still highly debated. Fundamental results will come from high-energy and spatial resolution X-ray observatories. Powerful winds driven by supermassive black holes (SMBHs) are likely the main mechanism through which SMBHs regulate their own growth and influence the host galaxy evolution. However, their origin and their capability to impact the large-scale environment are still highly debated. Fundamental results will come from high-energy and spatial resolution X-ray observatories. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Science White Paper Submitted for the Astro2020 Decadal Survey on Astronomy and Astrophysics

Magnetars as Astrophysical Laboratories of Extreme Quantum Electrodynamics: The Case for a Compton Telescope

Authors: Zorawar Wadiasingh, George Younes, Matthew G. Baring, Alice K. Harding, Peter L. Gonthier, Kun Hu, Alexander van der Horst, Silvia Zane, Chryssa Kouveliotou, Andrei M. Beloborodov, Chanda Prescod-Weinstein, Tanmoy Chattopadhyay, Sunil Chandra, Constantinos Kalapotharakos, Kyle Parfrey, Harsha Blumer, Demos Kazanas

Abstract: A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This p… ▽ More A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: 11 pages, 4 figures, Astro2020 Science White Paper submitted to the National Academies of Sciences

High-Energy Polarimetry - a new window to probe extreme physics in AGN jets

Authors: B. Rani, H. Zhang, S. D. Hunter, F. Kislat, M. Böttcher, J. E. McEnery, D. J. Thompson, D. Giannios, F. Guo, H. Li, M. Baring, I. Agudo, S. Buson, M. Petropoulou, V. Pavlidou, E. Angelakis, I. Myserlis, Z. Wadiasingh, R. M. Curado da Silva, P. Kilian, S. Guiriec, V. V. Bozhilov, J. Hodgson, S. Antón, D. Kazanas , et al. (9 additional authors not shown)

Abstract: The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is m… ▽ More The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: submitted to Astro2020 (Astronomy and Astrophysics Decadal Survey)

Electromagnetic probes of primordial black holes as dark matter

Authors: Y. Ali-Haimoud, S. Clesse, J. Garcia-Bellido, A. Kashlinsky, L. Wyrzykowski, A. Achucarro, L. Amendola, J. Annis, A. Arbey, R. G. Arendt, F. Atrio-Barandela, N. Bellomo, K. Belotsky, J-L. Bernal, S. Bird, V. Bozza, C. Byrnes, S. Calchi Novati, F. Calore, B. J. Carr, J. Chluba, I. Cholis, A. Cieplak, P. Cole, I. Dalianis , et al. (69 additional authors not shown)

Abstract: The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence wit… ▽ More The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing. △ Less

Submitted 12 March, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: Science whitepaper submitted to the Astro2020 Decadal Survey

Theoretical model of HD 163296 presently forming in-situ planets and comparison with the models of AS 209, HL Tau, and TW Hya

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: We fit an isothermal oscillatory density model to the disk of HD 163296 in which planets have presumably already formed and they are orbiting at least within the four observed dark gaps. This 156 AU large axisymmetric disk shows various physical properties comparable to those of AS 209, HL Tau, and TW Hya that we have modeled previously; but it compares best to AS 209. The disks of HD 163296 and A… ▽ More We fit an isothermal oscillatory density model to the disk of HD 163296 in which planets have presumably already formed and they are orbiting at least within the four observed dark gaps. This 156 AU large axisymmetric disk shows various physical properties comparable to those of AS 209, HL Tau, and TW Hya that we have modeled previously; but it compares best to AS 209. The disks of HD 163296 and AS 209 are comparable in size and they share similar values of the power-law index $k\approx 0$ (a radial density profile of the form $ρ(R)\propto R^{-1}$), the rotational parameter $β_0$ (to within a factor of 3); a relatively small inner core radius (although this parameter for HD 163296 is exceptionally small, $R_1\simeq 0.15$ AU, presumably due to unresolved planets in the inner 50 AU); the scale length $R_0$ and the Jeans gravitational frequency $Ω_J$ (to within factors of 1.4); the equation of state ($c_0^2/ρ_0$) and the central density $ρ_0$ (to within factors of 2); and the core angular velocity $Ω_0$ (to within a factor of 4.5). In the end, we compare all six nebular disks that we have modeled so far. △ Less

Submitted 3 March, 2019; originally announced March 2019.

Comments: Part 10 and the final part (HD 163296)

Theoretical model of the outer disk of TW Hya presently forming in-situ planets and comparison with models of AS 209 and HL Tau

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: We fit an isothermal oscillatory density model to the outer disk of TW Hya in which planets have presumably already formed and they are orbiting within four observed dark gaps. At first sight, this 52 AU small disk does not appear to be similar to our solar nebula; it shows several physical properties comparable to those in HL Tau (size $R_{\rm max}=102$ AU) and very few similarities to AS 209 (… ▽ More We fit an isothermal oscillatory density model to the outer disk of TW Hya in which planets have presumably already formed and they are orbiting within four observed dark gaps. At first sight, this 52 AU small disk does not appear to be similar to our solar nebula; it shows several physical properties comparable to those in HL Tau (size $R_{\rm max}=102$ AU) and very few similarities to AS 209 ($R_{\rm max}=144$ AU). We find a power-law density profile with index $k=-0.2$ (radial densities $ρ(R) \propto R^{-1.2}$) and centrifugal support against self-gravity so small that it virtually guarantees dynamical stability for millions of years of evolution to come. Compared to HL Tau, the scale length $R_0$ and the core size $R_1$ of TW Hya are smaller only by factors of $\sim$2, reflecting the disk's half size. On the opposite end, the Jeans frequency $Ω_J$ and the angular velocity $Ω_0$ of the smaller core of TW Hya are larger only by factors of $\sim$2. The only striking difference is that the central density ($ρ_0$) of TW Hya is 5.7 times larger than that of HL Tau, which is understood because the core of TW Hya is only half the size ($R_1$) of HL Tau and about twice as heavy ($Ω_J$). In the end, we compare the protostellar disks that we have modeled so far. △ Less

Submitted 3 March, 2019; v1 submitted 9 February, 2019; originally announced February 2019.

Comments: Updated version, Part 9 (TW Hya). arXiv admin note: text overlap with arXiv:1901.10642

Preliminary model of the outer disk of RU Lup presently showing only four dark gaps

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: Recent ALMA/DSHARP observations reveal that the disk of RU Lup shows presently only four dark gaps at large radii in which planets are presumably already orbiting. This gap arrangement is quite similar to the current structure of the outer planets in our solar system. Although there may be more undetected planets forming in the inner disk of RU Lup, the information obtained from the current observ… ▽ More Recent ALMA/DSHARP observations reveal that the disk of RU Lup shows presently only four dark gaps at large radii in which planets are presumably already orbiting. This gap arrangement is quite similar to the current structure of the outer planets in our solar system. Although there may be more undetected planets forming in the inner disk of RU Lup, the information obtained from the current observations prompted us to fit preliminary models with an isothermal oscillatory density profile to the outer disk of RU Lup. The best-fit model confirms that several physical properties of the disk of RU Lup are comparable to those of our solar nebula to within factors of 1.0-2.3; but there are some differences as well: central density ($ρ_0$), centrifugal support ($β_0$), and core angular velocity ($Ω_0$) in RU Lup are smaller by factors of 5.5, 8.0, and 18, respectively, whereas the inner core radius ($R_1$) of RU Lup is 25 times larger. The smaller values found for RU Lup are linked to the (apparently gap-empty) inner 10 AU of its disk, a region that is presently inadequately resolved. △ Less

Submitted 3 March, 2019; v1 submitted 31 January, 2019; originally announced February 2019.

Comments: Updated version, Part 8 (RU Lup). arXiv admin note: text overlap with arXiv:1901.10642

Theoretical models of the protostellar disks of AS 209 and HL Tau presently forming in-situ planets

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: We fit an isothermal oscillatory density model to two ALMA/DSHARP-observed disks, AS 209 and HL Tau, in which planets have presumably already formed and they are orbiting within the observed seven dark gaps in each system. These large disks are roughly similar to our solar nebula, albeit they exhibit milder radial density profiles and they enjoy lower centrifugal support. We find power-law density… ▽ More We fit an isothermal oscillatory density model to two ALMA/DSHARP-observed disks, AS 209 and HL Tau, in which planets have presumably already formed and they are orbiting within the observed seven dark gaps in each system. These large disks are roughly similar to our solar nebula, albeit they exhibit milder radial density profiles and they enjoy lower centrifugal support. We find power-law density profiles with index $k=0.0$ (radial densities $ρ(R) \propto R^{-1}$) and centrifugal support against self-gravity so small that it guarantees dynamical stability for millions of years of evolution. The scale lengths of the models differ only by a factor of 1.9, but the inner cores of the disks are very different: HL Tau's core is 8.0 times larger and 3.6 times denser than the core of AS 209. This results in four dark gaps having formed within the core of HL Tau, whereas no dark gap is found in the core of AS 209. On the other hand, the Jeans frequencies and the angular velocities of the cores are comparable to within factors of 1.9 and 1.6, respectively. △ Less

Submitted 3 March, 2019; v1 submitted 29 January, 2019; originally announced January 2019.

Comments: Updated version, Part 7 (AS 209 and HL Tau)

Models of Saturn's protoplanetary disk forming in-situ its regular satellites and innermost rings before the planet is formed

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: We fit an isothermal oscillatory density model of Saturn's protoplanetary disk to the present-day major satellites and innermost rings D/C and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Saturnian nebula. This disk does not look like the Jovian and Uranian disks that we modeled previously. It… ▽ More We fit an isothermal oscillatory density model of Saturn's protoplanetary disk to the present-day major satellites and innermost rings D/C and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Saturnian nebula. This disk does not look like the Jovian and Uranian disks that we modeled previously. Its power-law index is extremely steep ($k=-4.5$) and its radial extent is very narrow ($ΔR\lesssim 0.9$ Gm), its rotation parameter that measures centrifugal support against self-gravity is somewhat larger ($β_0=0.0431$), as is its radial scale length (395 km); but, as was expected, the size of the Saturnian disk, $R_{\rm max}=3.6$ Gm, takes just an intermediate value. On the other hand, the central density of the compact Saturnian core and its angular velocity are both comparable to that of Jupiter's core (density of $\approx 0.3$~g~cm$^{-3}$ in both cases, and rotation period of 5.0 d versus 6.8 d); and significantly less than the corresponding parameters of Uranus' core. As with the other primordial nebulae, this rotation is sufficiently slow to guarantee the disk's long-term stability against self-gravity induced instabilities for millions of years of evolution. △ Less

Submitted 3 March, 2019; v1 submitted 22 January, 2019; originally announced January 2019.

Comments: Updated version, Part 5 (Saturn). arXiv admin note: substantial text overlap with arXiv:1901.06448; text overlap with arXiv:1901.05131

Model of Neptune's protoplanetary disk forming in-situ its surviving regular satellites after Triton's capture and comparison of the protoplanetary disks of the four gaseous giants

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: We fit an isothermal oscillatory density model of Neptune's protoplanetary disk to the surviving regular satellites and its innermost ring and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Neptunian nebula. Neptune's regular moons suffered from the retrograde capture of Triton that disrupted th… ▽ More We fit an isothermal oscillatory density model of Neptune's protoplanetary disk to the surviving regular satellites and its innermost ring and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Neptunian nebula. Neptune's regular moons suffered from the retrograde capture of Triton that disrupted the system. Some moons may have been ejected, while others may have survived inside their potential minima. For this reason, the Neptunian nebula does not look like any of the nebulae that we modeled previously. In particular, there must be two density maxima deep inside the core of the nebula where no moons or rings are found nowadays. Even with this strong assumption, the recent discovery of the minor moon N XIV complicates further the modeling effort. With some additional assumptions, the Neptunian nebula still shares many similarities with the Uranian nebula, as was expected from the relative proximity and similar physical conditions of the two systems. For Neptune's primordial disk, we find a steep power-law index ($k=-3.0$), needed to accommodate the arrangement of the outer moons Larissa, N XIV, and Proteus. The rotation parameter that measures centrifugal support against self-gravity is quite small ($β_0=0.00808$), as is its radial scale length (13.6 km). The extent of the disk ($R_{\rm max}=0.12$ Gm) is a lot smaller than that of Uranus ($R_{\rm max}=0.60$ Gm) and Triton appears to be responsible for the truncation of the disk. The central density of the compact Neptunian core and its angular velocity are higher than but comparable to those of Uranus' core. In the end, we compare the models of the protoplanetary disks of the four gaseous giants. △ Less

Submitted 3 March, 2019; v1 submitted 25 January, 2019; originally announced January 2019.

Comments: Updated version, Part 6 (Neptune). arXiv admin note: text overlap with arXiv:1901.06448

Models of a protoplanetary disk forming in-situ the major Uranian satellites before the planet is formed

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: We fit an isothermal oscillatory density model of Uranus' protoplanetary disk to the present-day major satellites and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Uranian nebula. This disk does not at all look like the Jovian disk that we modeled previously. Its rotation parameter that measure… ▽ More We fit an isothermal oscillatory density model of Uranus' protoplanetary disk to the present-day major satellites and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Uranian nebula. This disk does not at all look like the Jovian disk that we modeled previously. Its rotation parameter that measures centrifugal support against self-gravity is a lot smaller ($β_0=0.00507$), as is the radial scale length (only 27.6 km) and the size of the disk (only 0.60 Gm). On the other hand, the central density of the compact Uranian core is higher by a factor of 180 and its core's angular velocity is about 2.3 times that of Jupiter's core (a rotation period of 3.0 d as opposed to 6.8 d). Yet, the rotation of the disk is sufficiently slow to guarantee its long-term stability against self-gravity induced instabilities for millions of years. △ Less

Submitted 3 March, 2019; v1 submitted 18 January, 2019; originally announced January 2019.

Comments: Updated version, Part 4 (Uranus). arXiv admin note: text overlap with arXiv:1901.05131

Models of a protoplanetary disk forming in-situ the Galilean and smaller nearby satellites before Jupiter is formed

Authors: D. M. Christodoulou, D. Kazanas

Abstract: We fit an isothermal oscillatory density model of Jupiter's protoplanetary disk to the present-day Galilean and other nearby satellites and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Jovian nebula. Although the radial density profile of Jupiter's disk was similar to that of the solar nebula,… ▽ More We fit an isothermal oscillatory density model of Jupiter's protoplanetary disk to the present-day Galilean and other nearby satellites and we determine the radial scale length of the disk, the equation of state and the central density of the primordial gas, and the rotational state of the Jovian nebula. Although the radial density profile of Jupiter's disk was similar to that of the solar nebula, its rotational support against self-gravity was very low, a property that also guaranteed its long-term stability against self-gravity induced instabilities for millions of years. △ Less

Submitted 3 March, 2019; v1 submitted 15 January, 2019; originally announced January 2019.

Comments: Updated version, Part 3 (Jupiter)

Conundrums and constraints concerning the formation of our solar system -- An alternative view

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: We have proposed an alternative model for the formation of our solar system that does not predict any mean-motion resonant interactions, planetary migrations, or self-gravitating instabilities in the very early isothermal solar nebula and before the protosun has formed. Within this context of nonviolent protoplanetary evolution over more than 10 million years, we examine some conundrums and constr… ▽ More We have proposed an alternative model for the formation of our solar system that does not predict any mean-motion resonant interactions, planetary migrations, or self-gravitating instabilities in the very early isothermal solar nebula and before the protosun has formed. Within this context of nonviolent protoplanetary evolution over more than 10 million years, we examine some conundrums and constraints that have been discovered from studies of small bodies in the present-day solar system (Jupiter and Neptune's Trojans and their differences from Kuiper belt objects, the irregular satellites of gaseous giants, the stability of the main asteroid belt, and the Late Heavy Bombardment). These issues that have caused substantial difficulties to models of violent formation do not appear to be problematic for the alternative model, and the reason is the complete lack of violent events during the evolution of protoplanets. △ Less

Submitted 3 March, 2019; v1 submitted 8 January, 2019; originally announced January 2019.

Comments: Updated version, Part 2 (Small Bodies in the Solar System)

On the formation of our solar system and many other protoplanetary systems observed by ALMA and SPHERE

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: In view of the many recent observations conducted by ALMA and SPHERE, it is becoming clear that protoplanetary disks form planets in narrow annular gaps at various distances from the central protostars before these protostars are actually fully formed and the gaseous disks have concluded their accretion/dispersal processes. This is in marked contrast to the many multi-planet exoplanetary systems t… ▽ More In view of the many recent observations conducted by ALMA and SPHERE, it is becoming clear that protoplanetary disks form planets in narrow annular gaps at various distances from the central protostars before these protostars are actually fully formed and the gaseous disks have concluded their accretion/dispersal processes. This is in marked contrast to the many multi-planet exoplanetary systems that do not conform to this pristine picture. This major discrepancy calls for an explanation. We provide such an explanation in this work, based on analytical solutions of the cylindrical isothermal Lane-Emden equation with rotation which do not depend on boundary conditions. These ``intrinsic'' solutions of the differential equation attract the solutions of the Cauchy problem and force them to oscillate permanently. The oscillations create density maxima in which dust and planetesimals are trapped and they can form protoplanetary cores during the very early isothermal evolution of such protoplanetary nebulae. We apply this model to our solar nebula that formed in-situ a minimum of eleven protoplanetary cores that have grown to planets which have survived undisturbed to the present day. We are also in the process of applying the same model to the ALMA/DSHARP disks. △ Less

Submitted 3 March, 2019; v1 submitted 8 January, 2019; originally announced January 2019.

Comments: Updated version, Part 1 (Solar Nebula). arXiv admin note: text overlap with arXiv:0706.3205

Gauss's Law and the Source for Poisson's Equation in Modified Gravity with Varying G

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: We have recently shown that the baryonic Tully-Fisher and Faber-Jackson relations imply that the gravitational "constant" $G$ in the force law varies with acceleration $a$ as $G\propto 1/a$ and vice versa. These results prompt us to reconsider every facet of Newtonian dynamics. Here we show that the integral form of Gauss's law in spherical symmetry remains valid in $G(a)$ gravity, but the differe… ▽ More We have recently shown that the baryonic Tully-Fisher and Faber-Jackson relations imply that the gravitational "constant" $G$ in the force law varies with acceleration $a$ as $G\propto 1/a$ and vice versa. These results prompt us to reconsider every facet of Newtonian dynamics. Here we show that the integral form of Gauss's law in spherical symmetry remains valid in $G(a)$ gravity, but the differential form depends on the precise distribution of $G(a)M(r)$, where $r$ is the distance from the origin and $M(r)$ is the mass distribution. We derive the differential form of Gauss's law in spherical symmetry, thus the source for Poisson's equation as well. Modified Newtonian dynamics (MOND) and weak-field Weyl gravity are asymptotic limits of $G(a)$ gravity at low and high accelerations, respectively. In these limits, we derive telling approximations to the source in spherical symmetry. It turns out that the source has a strong dependence on surface density $M/r^2$ everywhere in $a$-space except in the deep Newton-Weyl regime of very high accelerations. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: MNRAS, submitted

Gravitational Potential and Nonrelativistic Lagrangian in Modified Gravity with Varying G

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: We have recently shown that the baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations imply that the gravitational "constant" $G$ in the force law vary with acceleration $a$ as $1/a$. Here we derive the converse from first principles. First we obtain the gravitational potential for all accelerations and we formulate the Lagrangian for the central-force problem. Then action minimization imp… ▽ More We have recently shown that the baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations imply that the gravitational "constant" $G$ in the force law vary with acceleration $a$ as $1/a$. Here we derive the converse from first principles. First we obtain the gravitational potential for all accelerations and we formulate the Lagrangian for the central-force problem. Then action minimization implies the BTF/BFJ relations in the deep MOND limit as well as weak-field Weyl gravity in the Newtonian limit. The results show how we can properly formulate a nonrelativistic conformal theory of modified dynamics that reduces to MOND in its low-acceleration limit and to Weyl gravity in the opposite limit. An unavoidable conclusion is that $a_0$, the transitional acceleration in modified dynamics, does not have a cosmological origin and it may not even be constant among galaxies and galaxy clusters. △ Less

Submitted 21 November, 2018; originally announced November 2018.

Comments: To appear in MNRASL

Accretion Disk MHD Winds and Blazar Classification

Authors: Stella Boula, Demosthenes Kazanas, Apostolos Mastichiadis

Abstract: The Fermi Gamma-Ray Space Telescope observations of blazars show a strong correlation between the spectral index of their gamma-ray spectra and their synchrotron peak frequency $ν_{\rm{pk}}^{\rm{syn}}$; additionally, the rate of Compton Dominance of these sources also seems to be a function of $ν_{\rm{pk}}^{\rm{syn}}$. In this work, we adopt the assumption that the nonthermal emission of blazars i… ▽ More The Fermi Gamma-Ray Space Telescope observations of blazars show a strong correlation between the spectral index of their gamma-ray spectra and their synchrotron peak frequency $ν_{\rm{pk}}^{\rm{syn}}$; additionally, the rate of Compton Dominance of these sources also seems to be a function of $ν_{\rm{pk}}^{\rm{syn}}$. In this work, we adopt the assumption that the nonthermal emission of blazars is primarily due to radiation by a population of Fermi-accelerated electrons in a relativistic outflow (jet) along the symmetry axis of the blazar's accretion disk. Furthermore, we assume that the Compton component is related to an external photon field of photons, which are scattered from particles of the magnetohydrodynamic (MHD) wind emanating from the accretion disk. Our results reproduce well the aforementioned basic observational trends of blazar classification by varying just one parameter, namely the mass accretion rate onto the central black hole. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 5 pages, 3 figures, submitted to MNRAS Letters

Binary black hole growth by gas accretion in stellar clusters

Authors: Zacharias Roupas, Demosthenes Kazanas

Abstract: We show that binaries of stellar-mass black holes formed inside a young protoglobular cluster, can grow rapidly inside the cluster's core by accretion of the intracluster gas, before the gas may be depleted from the core. A black hole with mass of the order of eight solar masses can grow to values of the order of thirty five solar masses in accordance with recent gravitational waves signals observ… ▽ More We show that binaries of stellar-mass black holes formed inside a young protoglobular cluster, can grow rapidly inside the cluster's core by accretion of the intracluster gas, before the gas may be depleted from the core. A black hole with mass of the order of eight solar masses can grow to values of the order of thirty five solar masses in accordance with recent gravitational waves signals observed by LIGO. Due to the black hole mass increase, a binary may also harden. The growth of binary black holes in a dense protoglobular cluster through mass accretion indicates a potentially important formation and hardening channel. △ Less

Submitted 4 January, 2019; v1 submitted 11 September, 2018; originally announced September 2018.

Journal ref: Astronomy & Astrophysics, 621, L1 (2019)

Variable Nature of Magnetically-Driven Ultra-Fast Outflows

Authors: Keigo Fukumura, Demosthenes Kazanas, Chris Shrader, Ehud Behar, Francesco Tombesi, Ioannis Contopoulos

Abstract: Among a number of active galactic nuclei (AGNs) that drive ionized outflows in X-rays, a low-redshift (z = 0.184) quasar, PDS 456, is long known to exhibit one of the exemplary ultra-fast outflows (UFOs). However, the physical process of acceleration mechanism is yet to be definitively constrained. In this work, we model the variations of the Fe K UFO properties in PDS 456 over many epochs in X-ra… ▽ More Among a number of active galactic nuclei (AGNs) that drive ionized outflows in X-rays, a low-redshift (z = 0.184) quasar, PDS 456, is long known to exhibit one of the exemplary ultra-fast outflows (UFOs). However, the physical process of acceleration mechanism is yet to be definitively constrained. In this work, we model the variations of the Fe K UFO properties in PDS 456 over many epochs in X-ray observations in the context of magnetohydrodynamic (MHD) accretion-disk winds employed in our earlier studies of similar X-ray absorbers. We applied the model to the 2013/2014 XMM-Newton/NuSTAR spectra to determine the UFO's condition; namely, velocity, ionization parameter, column density and equivalent width (EW). Under some provisions on the dependence of X-ray luminosity on the accretion rate applicable to near-Eddington state, our photoionization calculations, coupled to a 2.5-dimensional MHD-driven wind model, can further reproduce the observed correlations of the UFO velocity and the anti-correlation of its EW with X-ray strength of PDS 456. This work demonstrates that UFOs, even without radiative pressure, can be driven as an extreme case purely by magnetic interaction while also producing the observed spectrum and correlations. △ Less

Submitted 27 August, 2018; originally announced August 2018.

Comments: accepted to ApJ Letter; 16 pages, 3 figures

Interposing a Varying Gravitational Constant Between Modified Newtonian Dynamics and Weak Weyl Gravity

Authors: Dimitris M. Christodoulou, Demosthenes Kazanas

Abstract: The Newtonian gravitational constant $G$ obeys the dimensional relation $[G] [M] [a] = [v]^4$, where $M$, $a$, and $v$ denote mass, acceleration, and speed, respectively. Since the baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations are observed facts, this relation implies that $G\, a = {\rm constant}$. This result cannot be obtained in Newtonian dynamics which cannot explain the origin… ▽ More The Newtonian gravitational constant $G$ obeys the dimensional relation $[G] [M] [a] = [v]^4$, where $M$, $a$, and $v$ denote mass, acceleration, and speed, respectively. Since the baryonic Tully-Fisher (BTF) and Faber-Jackson (BFJ) relations are observed facts, this relation implies that $G\, a = {\rm constant}$. This result cannot be obtained in Newtonian dynamics which cannot explain the origin of the BTF and BFJ relations. An alternative, modified Newtonian dynamics (MOND) assumes that $G=G_0$ is constant in space and derives naturally a characteristic constant acceleration $a=a_0$, as well as the BTF and BFJ relations. This is overkill and it comes with a penalty: MOND cannot explain the origin of $a_0$. A solid physical resolution of this issue is that $G \propto a^{-1}$, which implies that in lower-acceleration environments the gravitational force is boosted relative to its Newtonian value because $G$ increases. This eliminates all problems related to MOND's empirical cutoff $a_0$ and yields a quantitative method for mapping the detailed variations of $G(a)$ across each individual galaxy as well as on larger and smaller scales. On the opposite end, the large accelerations produced by $G(a)$ appear to be linked to the weak-field limit of the fourth-order theory of conformal Weyl gravity. △ Less

Submitted 25 June, 2018; originally announced June 2018.

Comments: An original letter to appear in MNRAS

Not an Oxymoron: Some X-ray Binary Pulsars with Enormous Spinup Rates Reveal Weak Magnetic Fields

Authors: D. M. Christodoulou, S. G. T. Laycock, D. Kazanas

Abstract: Three high-mass X-ray binaries have been discovered recently exhibiting enormous spinup rates. Conventional accretion theory predicts extremely high surface dipolar magnetic fields that we believe are unphysical. Instead, we propose quite the opposite scenario: some of these pulsars exhibit weak magnetic fields, so much so that their magnetospheres are crushed by the weight of inflowing matter. Th… ▽ More Three high-mass X-ray binaries have been discovered recently exhibiting enormous spinup rates. Conventional accretion theory predicts extremely high surface dipolar magnetic fields that we believe are unphysical. Instead, we propose quite the opposite scenario: some of these pulsars exhibit weak magnetic fields, so much so that their magnetospheres are crushed by the weight of inflowing matter. The enormous spinup rate is achieved before inflowing matter reaches the pulsar's surface as the penetrating inner disk transfers its excess angular momentum to the receding magnetosphere which, in turn, applies a powerful spinup torque to the pulsar. This mechanism also works in reverse: it spins a pulsar down when the magnetosphere expands beyond corotation and finds itself rotating faster than the accretion disk which then exerts a powerful retarding torque to the magnetic field and to the pulsar itself. The above scenaria cannot be accommodated within the context of neutron-star accretion processes occurring near spin equilibrium, thus they constitute a step toward a new theory of extreme (far from equilibrium) accretion phenomena. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: Submitted to MNRAS

Magnetized Disk-Winds in NGC 3783

Authors: Keigo Fukumura, Demosthenes Kazanas, Chris Shrader, Ehud Behar, Francesco Tombesi, Ioannis Contopoulos

Abstract: We analyze a 900-ks stacked Chandra/HETG spectrum of NGC 3783 in the context of magnetically-driven accretion-disk wind models in an effort to provide tight constraints on the global conditions of the underlying absorbers. Motivated by the earlier measurements of its absorption measure distribution (AMD) indicating X-ray-absorbing ionic columns that decrease slowly with decreasing ionization param… ▽ More We analyze a 900-ks stacked Chandra/HETG spectrum of NGC 3783 in the context of magnetically-driven accretion-disk wind models in an effort to provide tight constraints on the global conditions of the underlying absorbers. Motivated by the earlier measurements of its absorption measure distribution (AMD) indicating X-ray-absorbing ionic columns that decrease slowly with decreasing ionization parameter, we employ 2D magnetohydrodynamic (MHD) disk-wind models to describe the global outflow. We compute its photoionization structure along with the wind kinematic properties allowing us to further calculate in a self-consistent fashion the shapes of the major X-ray absorption lines. With the wind radial density profile determined by the AMD, the profiles of the ensemble of the observed absorption features are determined by the two global parameters of the MHD wind; i.e. disk inclination θ_obs and wind density normalization n_o. Considering the most significant absorption features in the (~1.8A-20A) range, we show that the MHD-wind is best described by n(r)~6.9e11(r/ro)^-1.15 [cm^-3] and θ_obs=44deg. We argue that winds launched by X-ray heating, radiation pressure or even MHD winds but with steeper radial density profiles are strongly disfavored by data. Considering the properties of Fe K band absorption features (i.e. Fe xxv and Fe xxvi), while typically prominent in the AGN X-ray spectra, they appear to be weak in NGC 3783. For the specific parameters of our model obtained by fitting the AMD and the rest of absorption features, these features are found to be weak in agreement with observation. △ Less

Submitted 21 December, 2017; originally announced December 2017.

Comments: 37 pages, 7 figs, 3 tables; accepted to ApJ

A Striking Confluence Between Theory and Observations of High-Mass X-ray Binary Pulsars

Authors: D. M. Christodoulou, S. G. T. Laycock, D. Kazanas

Abstract: We analyse the most powerful X-ray outbursts from neutron stars in ten Magellanic high-mass X-ray binaries and three pulsating ultraluminous X-ray sources. Most of the outbursts rise to $L_{max}$ which is about the level of the Eddington luminosity, while the rest and more powerful outbursts also appear to recognize that limit when their emissions are assumed to be anisotropic and beamed toward ou… ▽ More We analyse the most powerful X-ray outbursts from neutron stars in ten Magellanic high-mass X-ray binaries and three pulsating ultraluminous X-ray sources. Most of the outbursts rise to $L_{max}$ which is about the level of the Eddington luminosity, while the rest and more powerful outbursts also appear to recognize that limit when their emissions are assumed to be anisotropic and beamed toward our direction. We use the measurements of pulsar spin periods $P_S$ and their derivatives $\dot{P_S}$ to calculate the X-ray luminosities $L_p$ in their faintest accreting ("propeller") states. In four cases with unknown $\dot{P_S}$, we use the lowest observed X-ray luminosities, which only adds to the heterogeneity of the sample. Then we calculate the ratios $L_p/L_{max}$ and we obtain an outstanding confluence of theory and observations from which we conclude that work done on both fronts is accurate and the results are trustworthy: sources known to reside on the lowest Magellanic propeller line are all located on/near that line, whereas other sources jump higher and reach higher-lying propeller lines. These jumps can be interpreted in only one way, higher-lying pulsars have stronger surface magnetic fields in agreement with empirical results in which $\dot{P_S}$ and $L_p$ values were not used. △ Less

Submitted 9 January, 2018; v1 submitted 12 December, 2017; originally announced December 2017.

Comments: Added LMC X-4 and commented on the cyclotron absorption line of SMC X-2. 4 pages, 1 figure, 2 tables, submitted to MNRASL

Electron positron pair flow and current composition in the pulsar magnetosphere

Authors: Gabriele Brambilla, Constantinos Kalapotharakos, Andrey Timokhin, Alice Harding, Demosthenes Kazanas

Abstract: We performed ab-initio Particle-In-Cell (PIC) simulations of a pulsar magnetosphere with electron-positron plasma produced only in the regions close to the neutron star surface. We study how the magnetosphere transitions from the vacuum to a nearly force-free configuration. We compare the resulting force-free like configuration with ones obtained in a PIC simulation where particles are injected ev… ▽ More We performed ab-initio Particle-In-Cell (PIC) simulations of a pulsar magnetosphere with electron-positron plasma produced only in the regions close to the neutron star surface. We study how the magnetosphere transitions from the vacuum to a nearly force-free configuration. We compare the resulting force-free like configuration with ones obtained in a PIC simulation where particles are injected everywhere as well as with macroscopic force-free simulations. We found that although both PIC solutions have similar structure of electromagnetic fields and current density distributions, they have different particle density distribution. In fact in the injection from the surface solution, electrons and positrons counterstream only along parts of the return current regions and most of the particles leave the magnetosphere without returning to the star. We also found that pair production in the outer magnetosphere is not critical for filling the whole magnetosphere with plasma. We study how the current density distribution supporting the global electromagnetic configuration is formed by analyzing particle trajectories. We found that electrons precipitate to the return current layer inside the light cylinder and positrons precipitate to the current sheet outside the light cylinder by crossing magnetic field lines contributing to the charge density distribution required by the global electrodynamics. Moreover, there is a population of electrons trapped in the region close to the Y-point. On the other hand the most energetic positrons are accelerated close to the Y-point. These processes can have observational signatures that, with further modeling efforts, would help to distinguish this particular magnetosphere configuration from others. △ Less

Submitted 3 March, 2018; v1 submitted 10 October, 2017; originally announced October 2017.

Comments: 17 pages, 15 figures

Journal ref: Gabriele Brambilla et al 2018 ApJ 858 81

3D Kinetic Pulsar Magnetosphere Models: Connecting to Gamma-Ray Observations

Authors: Constantinos Kalapotharakos, Gabriele Brambilla, Andrey Timokhin, Alice K. Harding, Demosthenes Kazanas

Abstract: We present 3D global kinetic pulsar magnetosphere models, where the charged particle trajectories and the corresponding electromagnetic fields are treated self-consistently. For our study, we have developed a cartesian 3D relativistic particle-in-cell code that incorporates the radiation reaction forces. We describe our code and discuss the related technical issues, treatments, and assumptions. In… ▽ More We present 3D global kinetic pulsar magnetosphere models, where the charged particle trajectories and the corresponding electromagnetic fields are treated self-consistently. For our study, we have developed a cartesian 3D relativistic particle-in-cell code that incorporates the radiation reaction forces. We describe our code and discuss the related technical issues, treatments, and assumptions. Injecting particles up to large distances in the magnetosphere, we apply arbitrarily low to high particle injection rates and get an entire spectrum of solutions from close to the Vacuum-Retarded-Dipole to close to the Force-Free solution, respectively. For high particle injection rates (close to FF solutions) significant accelerating electric field components are confined only near the equatorial current sheet outside the light-cylinder. A judicious interpretation of our models allows the calculation of the particle emission and consequently the derivation of the corresponding realistic high-energy sky-maps and spectra. Using model parameters that cover the entire range of spin-down powers of Fermi young and millisecond pulsars, we compare the corresponding model $γ$-ray light-curves, cutoff energies, and total $γ$-ray luminosities with those observed by Fermi to discover a dependence of the particle injection-rate, $\mathcal{F}$, on the spin-down power, $\dot{\mathcal{E}}$, indicating an increase of $\mathcal{F}$ with $\dot{\mathcal{E}}$. Our models guided by Fermi observations provide field-structures and particle distributions that are not only consistent with each other but also able to reproduce a broad range of the observed $γ$-ray phenomenology of both young and millisecond pulsars. △ Less

Submitted 15 March, 2018; v1 submitted 9 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ (24 pages, 22 figures). This version has some corrections and additions