Showing 1–50 of 93 results for author: Bhattacharjee, A

Magnetized compressible turbulence with a fluctuation dynamo and Reynolds numbers over a million

Authors: James R. Beattie, Christoph Federrath, Ralf S. Klessen, Salvatore Cielo, Amitava Bhattacharjee

Abstract: Supersonic magnetohydrodynamic (MHD) turbulence is a ubiquitous state for many astrophysical plasmas. However, even the basic statistics for this type of turbulence remains uncertain. We present results from supersonic MHD turbulence simulations at unparalleled resolutions, with plasma Reynolds numbers of over a million. In the kinetic energy spectrum we find a break between the scales that are do… ▽ More Supersonic magnetohydrodynamic (MHD) turbulence is a ubiquitous state for many astrophysical plasmas. However, even the basic statistics for this type of turbulence remains uncertain. We present results from supersonic MHD turbulence simulations at unparalleled resolutions, with plasma Reynolds numbers of over a million. In the kinetic energy spectrum we find a break between the scales that are dominated by kinetic energy, with spectral index $-2$, and those that become strongly magnetized, with spectral index $-3/2$. By analyzing the Helmholtz decomposed kinetic energy spectrum, we find that the compressible modes are not passively mixed through the cascade of the incompressible modes. At high magnetic Reynolds number, above $10^5$, we find a power law in the magnetic energy spectrum with spectral index $-9/5$. On the strongly magnetized, subsonic scales the plasma tends to self-organize into locally relaxed regions, where there is strong alignment between the current density, magnetic field, velocity field and vorticity field, depleting both the nonlinearities and magnetic terms in the MHD equations, which we attribute to plasma relaxation on scales where the magnetic fluctuations evolve on shorter timescales than the velocity fluctuations. This process constrains the cascade to inhomogenous, volume-poor, fractal surfaces between relaxed regions, which has significant repercussions for understanding the nature of magnetized turbulence in astrophysical plasmas and the saturation of the fluctuation dynamo. △ Less

Submitted 29 May, 2024; v1 submitted 26 May, 2024; originally announced May 2024.

Comments: 22 pages. 10 figures. Currently the world's largest MHD turbulence simulation on a grid of 10,080^3 (Re and Rm over a million). Comments welcome, email james.beattie@princeton.edu

Three-dimensional plasmoid-mediated reconnection and turbulence in Hall magnetohydrodynamics

Authors: Yi-Min Huang, Amitava Bhattacharjee

Abstract: Plasmoid instability accelerates reconnection in collisional plasmas by transforming a laminar reconnection layer into numerous plasmoids connected by secondary current sheets in two dimensions (2D) and by fostering self-generated turbulent reconnection in three dimensions (3D). In large-scale astrophysical and space systems, plasmoid instability likely initiates in the collisional regime but may… ▽ More Plasmoid instability accelerates reconnection in collisional plasmas by transforming a laminar reconnection layer into numerous plasmoids connected by secondary current sheets in two dimensions (2D) and by fostering self-generated turbulent reconnection in three dimensions (3D). In large-scale astrophysical and space systems, plasmoid instability likely initiates in the collisional regime but may transition into the collisionless regime as the fragmentation of the current sheet progresses toward kinetic scales. Hall MHD models are widely regarded as a simplified yet effective representation of the transition from collisional to collisionless reconnection. However, plasmoid instability in 2D Hall MHD simulations often leads to a single-X-line reconnection configuration, which significantly differs from fully kinetic particle-in-cell simulation results. This study shows that single-X-line reconnection is less likely to occur in 3D compared to 2D. Moreover, depending on the Lundquist number and the ratio between the system size and the kinetic scale, Hall MHD can also realize 3D self-generated turbulent reconnection. We analyze the features of the self-generated turbulent state, including the energy power spectra and the scale dependence of turbulent eddy anisotropy. △ Less

Submitted 30 April, 2024; originally announced April 2024.

X-ray imaging and electron temperature evolution in laser-driven magnetic reconnection experiments at the National Ignition Facility

Authors: V. Valenzuela-Villaseca, J. M. Molina, D. B. Schaeffer, S. Malko, J. Griff-McMahon, K. Lezhnin, M. J. Rosenberg, S. X. Hu, D. Kalantar, C. Trosseille, H. -S. Park, B. A. Remington, G. Fiksel, D. Uzdensky, A. Bhattacharjee, W. Fox

Abstract: We present results from X-ray imaging of high-aspect-ratio magnetic reconnection experiments driven at the National Ignition Facility. Two parallel, self-magnetized, elongated laser-driven plumes are produced by tiling 40 laser beams. A magnetic reconnection layer is formed by the collision of the plumes. A gated X-ray framing pinhole camera with micro-channel plate (MCP) detector produces multipl… ▽ More We present results from X-ray imaging of high-aspect-ratio magnetic reconnection experiments driven at the National Ignition Facility. Two parallel, self-magnetized, elongated laser-driven plumes are produced by tiling 40 laser beams. A magnetic reconnection layer is formed by the collision of the plumes. A gated X-ray framing pinhole camera with micro-channel plate (MCP) detector produces multiple images through various filters of the formation and evolution of both the plumes and current sheet. As the diagnostic integrates plasma self-emission along the line of sight, 2-dimensional electron temperature maps $\langle T_e \rangle_Y$ are constructed by taking the ratio of intensity of these images obtained with different filters. The plumes have a characteristic temperature $\langle T_e \rangle_Y = 240 \pm 20$ eV at 2 ns after the initial laser irradiation and exhibit a slow cooling up to 4 ns. The reconnection layer forms at 3 ns with a temperature $\langle T_e \rangle_Y = 280 \pm 50$ eV as the result of the collision of the plumes. The error bars of the plumes and current sheet temperatures separate at $4$ ns, showing the heating of the current sheet from colder inflows. Using a semi-analytical model, we find that the observed heating of the current sheet is consistent with being produced by electron-ion drag, rather than the conversion of magnetic to kinetic energy. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: Submitted to Physics of Plasmas. 19 pages (total), 14 figures, 2 tables

Relaxation of weakly collisional plasma: continuous spectra, Landau eigenmodes, and transition from the collisionless to the fluid limit

Authors: Uddipan Banik, Amitava Bhattacharjee

Abstract: The relaxation of a weakly collisional plasma is described by the Boltzmann-Poisson equations with the Lenard-Bernstein collision operator. We perform a perturbative analysis of these equations, and obtain, for the first time, exact analytic solutions, enabling definitive resolutions to long-standing controversies regarding the impact of weak collisions on continuous spectra and Landau eigenmodes.… ▽ More The relaxation of a weakly collisional plasma is described by the Boltzmann-Poisson equations with the Lenard-Bernstein collision operator. We perform a perturbative analysis of these equations, and obtain, for the first time, exact analytic solutions, enabling definitive resolutions to long-standing controversies regarding the impact of weak collisions on continuous spectra and Landau eigenmodes. Unlike some previous studies, we retain both damping and diffusion terms in the collision operator. We find that the linear response is a temporal convolution of a continuum that depends on the continuous velocities of particles, and discrete normal modes that encapsulate coherent oscillations. The normal modes are exponentially damped over time due to collective effects (Landau damping) as well as collisional dissipation. The continuum is also damped by collisions but somewhat differently. Up to a collision time, which is the inverse of the collision frequency $ν_{\mathrm{c}}$, the continuum decay is driven by velocity diffusion and occurs super-exponentially over a timescale $\sim ν^{-1/3}_{\mathrm{c}}$. After a collision time, however, the continuum decay is driven by the collisional damping of particle velocities and diffusion of their positions, and occurs exponentially over a timescale $\sim ν_{\mathrm{c}}$. This hitherto unknown, slow exponential decay causes perturbations to damp the most on scales comparable to the mean free path, but very slowly on larger scales, which establishes the local thermal equilibrium, the essence of the fluid limit. The long-term decay of the response is driven by the normal modes on scales smaller than the mean free path, but, on larger scales, is governed by the slowly decaying continuum and the least damped normal mode. Our analysis firmly establishes a long-sought connection between the collisionless and fluid limits of weakly collisional plasmas. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: Submitted to Physical Review X; 24 pages, 5 figures

Taking control of compressible modes: bulk viscosity and the turbulent dynamo

Authors: James R. Beattie, Christoph Federrath, Neco Kriel, Justin Kin Jun Hew, Amitava Bhattacharjee

Abstract: Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, and yet, due to Stokes' hypothesis, a means to carefully control the decay of compressible modes in these systems has largely been neglected. This is especially important for small-scale, turbulent dynamo processes, which are known to be sensitive to the effects of compression. To control the viscou… ▽ More Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, and yet, due to Stokes' hypothesis, a means to carefully control the decay of compressible modes in these systems has largely been neglected. This is especially important for small-scale, turbulent dynamo processes, which are known to be sensitive to the effects of compression. To control the viscous properties of the compressible modes, we perform supersonic, visco-resistive dynamo simulations with additional bulk viscosity $ν_{\rm bulk}$, deriving a new $ν_{\rm bulk}$ Reynolds number $\rm{Re}_{\rm bulk}$, and viscous Prandtl number $\rm{P}ν\equiv \rm{Re}_{\rm bulk} / \rm{Re}_{\rm shear}$, where $\rm{Re}_{\rm shear}$ is the shear viscosity Reynolds number. For $10^{-3} \leq \rm{P}ν\leq \infty$, we explore a broad range of statistics critical to the dynamo problem, including the integral and spectral energy ratios, growth rates, and the magnetic $E_{\rm mag}(k)$ and kinetic $E_{\rm kin}(k)$ energy spectrum. We derive a general framework for decomposing $E_{\rm mag}$ growth rates into incompressible and compressible terms via orthogonal tensor decompositions of $\nabla\otimes\mathbf{v}$, where $\mathbf{v}$ is the fluid velocity. We find that compressible modes play a dual role, growing and decaying $E_{\rm mag}$, and that field-line stretching is the main driver of growth, even in supersonic dynamos. In the absence of $ν_{\rm bulk}$, compressible modes pile up on small-scales, creating an apparent spectral bottleneck, which disappears for $\rm{P}ν\approx 1$. As $\rm{P}ν$ decreases, compressible modes are dissipated at increasingly larger scales, in turn suppressing incompressible modes through a coupling between viscosity operators. We emphasise the importance of further understanding the role of $ν_{\rm bulk}$ in compressible astrophysical plasmas. △ Less

Submitted 11 December, 2023; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 28 pages, 18 figures, submitted to MNRAS. Comments welcome. Re-uploaded on 11 Dec with compressed large figures to reduce file size

Exact Calculation of Nonideal Fields Demonstrates Dominance of Injection in Relativistic Reconnection

Authors: Samuel R. Totorica, Seiji Zenitani, Shuichi Matsukiyo, Mami Machida, Kazuhiro Sekiguchi, Amitava Bhattacharjee

Abstract: Magnetic reconnection is an important source of energetic particles in systems ranging from astrophysics to the laboratory. The large separation of spatiotemporal scales involved makes it critical to determine the minimum physical model containing the necessary physics for modeling particle acceleration. By resolving the energy gain from ideal and nonideal magnetohydrodynamic electric fields self-… ▽ More Magnetic reconnection is an important source of energetic particles in systems ranging from astrophysics to the laboratory. The large separation of spatiotemporal scales involved makes it critical to determine the minimum physical model containing the necessary physics for modeling particle acceleration. By resolving the energy gain from ideal and nonideal magnetohydrodynamic electric fields self-consistently in kinetic particle-in-cell simulations of reconnection, we conclusively show the dominant role of the nonideal field for the early stage of energization known as injection. The importance of the nonideal field increases with magnetization, guide field, and in three-dimensions, indicating its general importance for reconnection in natural astrophysical systems. We obtain the statistical properties of the injection process from the simulations, paving the way for the development of extended MHD models capable of accurately modeling particle acceleration in large-scale systems. The novel analysis method developed in this study can be applied broadly to give new insight into a wide range of processes in plasma physics. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Journal ref: The Astrophysical Journal Letters, Volume 952, Number 1 (2023)

Resonant instabilities mediated by drag and electrostatic interactions in laboratory and astrophysical dusty plasmas

Authors: Ben Y. Israeli, Amitava Bhattacharjee, Hong Qin

Abstract: Dusty plasmas are known to support a diverse range of instabilities, including both generalizations of standard plasma instabilities and ones caused by effects specific to dusty systems. It has been recently demonstrated that a novel broad class of streaming instabilities, termed resonant drag instabilities (RDIs), can be attributed to a particular resonance phenomenon, manifested by defective eig… ▽ More Dusty plasmas are known to support a diverse range of instabilities, including both generalizations of standard plasma instabilities and ones caused by effects specific to dusty systems. It has been recently demonstrated that a novel broad class of streaming instabilities, termed resonant drag instabilities (RDIs), can be attributed to a particular resonance phenomenon, manifested by defective eigenvalues of the linearized dust/fluid system. In this work, it is demonstrated that this resonance phenomenon is not unique to RDIs and can be used as a framework to understand a wider range of instabilities, termed resonant instabilities. Particular attention is given to the filamentary ionization instability seen in laboratory dusty plasmas and to the two-stream instability. It is shown that, due to the commonalities in underlying physics between the dust-ion-acoustic two-stream instability and the acoustic RDI, these instabilities should be relevant in strongly overlapping regimes in astrophysical dusty plasmas. It is proposed that a similar overlap in the experimental accessibility of these modes (and of the filamentary instability) allows for the possibility of experimental investigation in the laboratory of complex and astrophysically relevant instability dynamics. △ Less

Submitted 20 July, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: 21 pages, 15 figures

Reconnection-Driven Energy Cascade in Magnetohydrodynamic Turbulence

Authors: Chuanfei Dong, Liang Wang, Yi-Min Huang, Luca Comisso, Timothy A. Sandstrom, Amitava Bhattacharjee

Abstract: Magnetohydrodynamic turbulence regulates the transfer of energy from large to small scales in many astrophysical systems, including the solar atmosphere. We perform three-dimensional magnetohydrodynamic simulations with unprecedentedly large magnetic Reynolds number to reveal how rapid reconnection of magnetic field lines changes the classical paradigm of the turbulent energy cascade. By breaking… ▽ More Magnetohydrodynamic turbulence regulates the transfer of energy from large to small scales in many astrophysical systems, including the solar atmosphere. We perform three-dimensional magnetohydrodynamic simulations with unprecedentedly large magnetic Reynolds number to reveal how rapid reconnection of magnetic field lines changes the classical paradigm of the turbulent energy cascade. By breaking elongated current sheets into chains of small magnetic flux ropes (or plasmoids), magnetic reconnection leads to a new range of turbulent energy cascade, where the rate of energy transfer is controlled by the growth rate of the plasmoids. As a consequence, the turbulent energy spectra steepen and attain a spectral index of -2.2 that is accompanied by changes in the anisotropy of turbulence eddies. The omnipresence of plasmoids and their consequences on, e.g., solar coronal heating, can be further explored with current and future spacecraft and telescopes. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 32 pages, 8 figures, the world's largest 3D MHD turbulence simulation using a fifth-order scheme

Journal ref: Science Advances 8, eabn7627 (2022)

Radio Dichotomy in Quasars with H$β$ FWHM greater than $15,000$ km\,s$^{-1}$

Authors: A. Chakraborty, A. Bhattacharjee, M. S. Brotherton, R. Chatterjee, S. Chatterjee, M. Gilbert

Abstract: It has been inferred from large unbiased samples that $10\%$-$15\%$ of all quasars are radio-loud (RL). Using the quasar catalog from the Sloan Digital Sky Survey, we show that the radio-loud fraction (RLF) for high broad line (HBL) quasars, containing H$β$ FWHM greater than $15,000$ km s$^{-1}$, is $\sim 57 \%$. While there is no significant difference between the RL and radio-quiet (RQ) populati… ▽ More It has been inferred from large unbiased samples that $10\%$-$15\%$ of all quasars are radio-loud (RL). Using the quasar catalog from the Sloan Digital Sky Survey, we show that the radio-loud fraction (RLF) for high broad line (HBL) quasars, containing H$β$ FWHM greater than $15,000$ km s$^{-1}$, is $\sim 57 \%$. While there is no significant difference between the RL and radio-quiet (RQ) populations in our sample in terms of their black hole mass, Eddington ratio, and covering fraction (CF), optical continuum luminosity of the RL quasars are higher. The similarity in the distribution of their CF indicates that our analysis is unbiased in terms of the viewing angle of the HBL RL and RQ quasars. Hence, we conclude that the accretion disc luminosity of the RL quasars in our HBL sample is higher, which indicates a connection between a brighter disc and a more prominent jet. By comparing them with the non-HBL H$β$ broad emission line quasars, we find that the HBL sources have the lowest Eddington ratios in addition to having a very high RLF. That is consistent with the theories of jet formation, in which jets are launched from low Eddington ratio accreting systems. We find that the [O III] narrow emission line is stronger in the RL compared to RQ quasars in our HBL sample, which is consistent with previous findings in the literature, and may be caused by the interaction of the narrow line gas with the jet. △ Less

Submitted 21 August, 2022; originally announced August 2022.

Comments: Accepted for publication in MNRAS

Do chaotic field lines cause fast reconnection in coronal loops?

Authors: Yi-Min Huang, Amitava Bhattacharjee

Abstract: Over the past decade, Boozer has argued that three-dimensional (3D) magnetic reconnection fundamentally differs from two-dimensional (2D) reconnection due to the fact that the separation between any pair of neighboring field lines almost always increases exponentially over distance in a 3D magnetic field. According to Boozer, this feature makes 3D field-line mapping chaotic and exponentially sensi… ▽ More Over the past decade, Boozer has argued that three-dimensional (3D) magnetic reconnection fundamentally differs from two-dimensional (2D) reconnection due to the fact that the separation between any pair of neighboring field lines almost always increases exponentially over distance in a 3D magnetic field. According to Boozer, this feature makes 3D field-line mapping chaotic and exponentially sensitive to small non-ideal effects; consequently, 3D reconnection can occur without intense current sheets. We test Boozer's theory via ideal and resistive reduced magnetohydrodynamic simulations of the Boozer-Elder coronal loop model driven by sub-Alfvenic footpoint motions [A. H. Boozer and T. Elder, Physics of Plasmas 28, 062303 (2021)]. Our simulation results significantly differ from their predictions. The ideal simulation shows that Boozer and Elder under-predict the intensity of current density due to missing terms in their reduced model equations. Furthermore, resistive simulations of varying Lundquist numbers show that the maximal current density scales linearly rather than logarithmically with the Lundquist number. △ Less

Submitted 25 November, 2022; v1 submitted 14 August, 2022; originally announced August 2022.

Comments: 14 pages, 9 figures. Accepted for publication in Physics of Plasmas

Research Opportunities in Plasma Astrophysics

Authors: Stuart Bale, Amitava Bhattacharjee, Fausto Cattaneo, Jemes Drake, Hantao Ji, Marty Lee, Hui Li, Edison Liang, Marc Pound, Stewart Prager, Eliot Quataert, Bruce Remington, Robert Rosner, Dmitri Ryutov, Edward Thomas Jr, Ellen Zweibel

Abstract: Major scientific questions and research opportunities are described on 10 unprioritized plasma astrophysics topics: (1) magnetic reconnection, (2) collisionless shocks and particle acceleration, (3) waves and turbulence, (4) magnetic dynamos, (5) interface and shear instabilities, (6) angular momentum transport, (7) dusty plasmas, (8) radiative hydrodynamics, (9) relativistic, pair-dominated and s… ▽ More Major scientific questions and research opportunities are described on 10 unprioritized plasma astrophysics topics: (1) magnetic reconnection, (2) collisionless shocks and particle acceleration, (3) waves and turbulence, (4) magnetic dynamos, (5) interface and shear instabilities, (6) angular momentum transport, (7) dusty plasmas, (8) radiative hydrodynamics, (9) relativistic, pair-dominated and strongly magnetized plasmas, (10) jets and outflows. Note that this is a conference report from a Workshop on Opportunities in Plasma Astrophysics (WOPA, https://w3.pppl.gov/conferences/2010/WOPA/) in January 2010, that attracted broad representation from the community and was supported by the U.S. Department of Energy, National Aeronautics and Space Administration, National Science Foundation, American Physical Society's Topical Group for Plasma Astrophysics and Division of Plasma Physics, and Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas. Although there has been much planning and many developments in both science and infrastructure since the report was written, most of the motivation, priorities, problems and technical challenges discussed therein remain unaddressed and are relevant at the time of posting. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: The WOPA committee members are alphabetically listed as authors above. Each chapter was authored by a working group with their full membership listed in Appendix D. Conference participants are listed in Appendix B. Editor is Patti Wieser. 130 pages, 10 chapters, 4 appendices, 26 figures

On large-scale dynamos with stable stratification and the application to stellar radiative zones

Authors: Valentin Skoutnev, Jonathan Squire, Amitava Bhattacharjee

Abstract: Our understanding of large-scale magnetic fields in stellar radiative zones remains fragmented and incomplete. Such magnetic fields, which must be produced by some form of dynamo mechanism, are thought to dominate angular-momentum transport, making them crucial to stellar evolution. A major difficulty is the effect of stable stratification, which generally suppresses dynamo action. We explore the… ▽ More Our understanding of large-scale magnetic fields in stellar radiative zones remains fragmented and incomplete. Such magnetic fields, which must be produced by some form of dynamo mechanism, are thought to dominate angular-momentum transport, making them crucial to stellar evolution. A major difficulty is the effect of stable stratification, which generally suppresses dynamo action. We explore the effects of stable stratification on mean-field dynamo theory with a particular focus on a non-helical large-scale dynamo (LSD) mechanism known as the magnetic shear-current effect. We find that the mechanism is robust to increasing stable stratification as long as the original requirements for its operation are met: a source of shear and non-helical magnetic fluctuations (e.g. from a small-scale dynamo). Both are plausibly sourced in the presence of differential rotation. Our idealized direct numerical simulations, supported by mean-field theory, demonstrate the generation of near equipartition large-scale toroidal fields. Additionally, a scan over magnetic Reynolds number shows no change in the growth or saturation of the LSD, providing good numerical evidence of a dynamo mechanism resilient to catastrophic quenching, which has been an issue for helical dynamos. These properties -- the absence of catastrophic quenching and robustness to stable stratification -- make the mechanism a plausible candidate for generating in-situ large-scale magnetic fields in stellar radiative zones. △ Less

Submitted 16 September, 2022; v1 submitted 3 March, 2022; originally announced March 2022.

Transonic accretion and winds around Pseudo-Kerr black holes and comparison with general relativistic solutions

Authors: Abhrajit Bhattacharjee, Sandip K. Chakrabarti, Dipak Debnath

Abstract: Spectral and timing properties of accretion flows on a black hole depend on their density and temperature distributions, which in turn come from the underlying dynamics. Thus, an accurate description of the flow which includes hydrodynamics and radiative transfer is a must to interpret the observational results. In the case of non-rotating black holes, Pseudo-Newtonian description of surrounding s… ▽ More Spectral and timing properties of accretion flows on a black hole depend on their density and temperature distributions, which in turn come from the underlying dynamics. Thus, an accurate description of the flow which includes hydrodynamics and radiative transfer is a must to interpret the observational results. In the case of non-rotating black holes, Pseudo-Newtonian description of surrounding space-time enables one to make a significant progress in predicting spectral and timing properties. This formalism is lacking for spinning black holes. In this paper, we show that there exists an exact form of 'natural' potential derivable from the general relativistic (GR) radial momentum equation. Use of this potential in an otherwise Newtonian set of equations allows to describe transonic flows very accurately as is evidenced by comparing with solutions obtained from the full GR framework. We study the properties of the critical points and the centrifugal pressure supported shocks in the parameter space spanned by the specific energy and the angular momentum, and compare with the results of GR hydrodynamics. We show that this potential can safely be used for the entire range of Kerr parameter $-1<a<1$ for modeling of observational results around spinning black holes. We assume the flow to be inviscid. Thus, it is non-dissipative with constant energy and angular momentum. These assumptions are valid very close to the black hole as the infall timescale is much shorter as compared to the viscous timescale. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: 14 pages, 7 figures, accepted for publication in RAA

Properties of Very Broad Line MgII Radio-Loud and Radio-Quiet Quasars

Authors: Avinanda Chakraborty, Anirban Bhattacharjee, Suchetana Chatterjee

Abstract: We perform an analysis of the properties of radio-loud (RL) and radio-quiet (RQ) quasars with MgII broad emission line (i-band magnitude $\leq 19.1$ and z $\leq 1.9$), selected from the parent sample of SDSS DR7 catalogue. For sources with full-width half maxima (FWHM) greater than \mbox{15,000 km s$^{-1}$} (very broad line sample; VBL) we find the radio loud fraction (RLF) to be about 40\%. To fu… ▽ More We perform an analysis of the properties of radio-loud (RL) and radio-quiet (RQ) quasars with MgII broad emission line (i-band magnitude $\leq 19.1$ and z $\leq 1.9$), selected from the parent sample of SDSS DR7 catalogue. For sources with full-width half maxima (FWHM) greater than \mbox{15,000 km s$^{-1}$} (very broad line sample; VBL) we find the radio loud fraction (RLF) to be about 40\%. To further investigate this result we compare the bolometric luminosity, optical continuum luminosity, black hole (BH) mass and Eddington ratios of our VBL sample of RL and RQ quasars. Our analysis shows that in our VBL sample space, RL quasars have higher luminosities and BH mass than RQ quasars. The similarity in the distribution of their covering fraction (CF) shows that there is no difference in dust distribution between VBL RL and RQ quasars and hence dust is not affecting our results. We also find that there is no correlation of RL quasar properties with optical continuum luminosity and BH mass. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 6 pages, 4 figures, Accepted to the Special Issue "A New Window on the Radio Emission from Galaxies, Galaxy Clusters and Cosmic Web: Current Status and Perspectives" of MDPI Galaxies

Journal ref: Galaxies 2021, 9, 74

Landau Modes are Eigenmodes of Stellar Systems in the Limit of Zero Collisions

Authors: C. S. Ng, A. Bhattacharjee

Abstract: We consider the spectrum of eigenmodes in a stellar system dominated by gravitational forces in the limit of zero collisions. We show analytically and numerically using the Lenard-Bernstein collision operator that the Landau modes, which are not true eigenmodes in a strictly collisionless system (except for the Jeans unstable mode), become part of the true eigenmode spectrum in the limit of zero c… ▽ More We consider the spectrum of eigenmodes in a stellar system dominated by gravitational forces in the limit of zero collisions. We show analytically and numerically using the Lenard-Bernstein collision operator that the Landau modes, which are not true eigenmodes in a strictly collisionless system (except for the Jeans unstable mode), become part of the true eigenmode spectrum in the limit of zero collisions. Under these conditions, the continuous spectrum of true eigenmodes in the collisionless system, also known as the Case-van Kampen modes, is eliminated. Furthermore, since the background distribution function in a weakly collisional system can exhibit significant deviations from a Maxwellian distribution function over long times, we show that the spectrum of Landau modes can change drastically even in the presence of slight deviations from a Maxwellian, primarily through the appearance of weakly damped modes that may be otherwise heavily damped for a Maxwellian distribution. Our results provide important insights for developing statistical theories to describe thermal fluctuations in a stellar system, which are currently a subject of great interest for N-body simulations as well as observations of gravitational systems. △ Less

Submitted 19 October, 2021; v1 submitted 16 September, 2021; originally announced September 2021.

Comments: Accepted for publication in Astrophysical Journal

Journal ref: The Astrophysical Journal, 923:271 (2021)

Scaling of Small-Scale Dynamo Properties in the Rayleigh-Taylor Instability

Authors: Valentin Skoutnev, Elias R. Most, Amitava Bhattacharjee, Alexander A. Philippov

Abstract: We derive scaling relations based on free-fall and isotropy assumptions for the kinematic small-scale dynamo growth rate and amplification factor over the course of the mixing, saturation, and decay phases of the Rayleigh-Taylor instability (RTI) in a fully-ionized plasma. The scaling relations are tested using sets of three dimensional, visco-resistive MHD simulations of the RTI and found to hold… ▽ More We derive scaling relations based on free-fall and isotropy assumptions for the kinematic small-scale dynamo growth rate and amplification factor over the course of the mixing, saturation, and decay phases of the Rayleigh-Taylor instability (RTI) in a fully-ionized plasma. The scaling relations are tested using sets of three dimensional, visco-resistive MHD simulations of the RTI and found to hold in the saturation phase, but exhibit discrepancies during the mixing and decays phases, suggesting a need to relax either the free-fall or isotropy assumptions. Application of the scaling relations allows for quantitative prediction of the net amplification of magnetic energy in the kinematic dynamo phase and therefore a determination of whether the magnetic energy either remains sub-equipartition at all velocity scales or reaches equipartition with at least some scales of the turbulent kinetic energy in laboratory and astrophysical scenarios. As an example, we consider the dynamo in RTI-unstable regions of the outer envelope of a binary neutron star merger and predict that the kinematic regime of the small-scale dynamo ends on the time scale of nanoseconds and then reaches saturation on a timescale of microseconds, which are both fast compared to the millisecond relaxation time of the post-merger. △ Less

Submitted 3 November, 2021; v1 submitted 8 June, 2021; originally announced June 2021.

Journal ref: The Astrophysical Journal (2021), 921(1), 75

Weak Alfvénic turbulence in relativistic plasmas. Part 1. Dynamical equations and basic dynamics of interacting resonant triads

Authors: J. M. TenBarge, B. Ripperda, A. Chernoglazov, A. Bhattacharjee, J. F. Mahlmann, E. R. Most, J. Juno, Y. Yuan, A. A. Philippov

Abstract: Alfvén wave collisions are the primary building blocks of the non-relativistic turbulence that permeates the heliosphere and low-to-moderate energy astrophysical systems. However, many astrophysical systems such as gamma-ray bursts, pulsar and magnetar magnetospheres, and active galactic nuclei have relativistic flows or energy densities. To better understand these high energy systems, we derive r… ▽ More Alfvén wave collisions are the primary building blocks of the non-relativistic turbulence that permeates the heliosphere and low-to-moderate energy astrophysical systems. However, many astrophysical systems such as gamma-ray bursts, pulsar and magnetar magnetospheres, and active galactic nuclei have relativistic flows or energy densities. To better understand these high energy systems, we derive reduced relativistic MHD equations and employ them to examine weak Alfvénic turbulence, dominated by three-wave interactions, in reduced relativistic magnetohydrodynamics, including the force-free, infinitely magnetized limit. We compare both numerical and analytical solutions to demonstrate that many of the findings from non-relativistic weak turbulence are retained in the relativistic system. But, an important distinction in the relativistic limit is the inapplicability of a formally incompressible limit, i.e, there exists finite coupling to the compressible fast mode regardless of the strength of the magnetic field. Since fast modes can propagate across field lines, this mechanism provides a route for energy to escape strongly magnetized systems, e.g., magnetar magnetospheres. However, we find that the fast-Alfvén coupling is diminished in the limit of oblique propagation. △ Less

Submitted 23 February, 2022; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: Accepted

Weak Alfvénic turbulence in relativistic plasmas II: Current sheets and dissipation

Authors: B. Ripperda, J. F. Mahlmann, A. Chernoglazov, J. M. TenBarge, E. R. Most, J. Juno, Y. Yuan, A. A. Philippov, A. Bhattacharjee

Abstract: Alfvén waves as excited in black hole accretion disks and neutron star magnetospheres are the building blocks of turbulence in relativistic, magnetized plasmas. A large reservoir of magnetic energy is available in these systems, such that the plasma can be heated significantly even in the weak turbulence regime. We perform high-resolution three-dimensional simulations of counter-propagating Alfvén… ▽ More Alfvén waves as excited in black hole accretion disks and neutron star magnetospheres are the building blocks of turbulence in relativistic, magnetized plasmas. A large reservoir of magnetic energy is available in these systems, such that the plasma can be heated significantly even in the weak turbulence regime. We perform high-resolution three-dimensional simulations of counter-propagating Alfvén waves, showing that an $E_{B_{\perp}}(k_{\perp}) \propto k_{\perp}^{-2}$ energy spectrum develops as a result of the weak turbulence cascade in relativistic magnetohydrodynamics and its infinitely magnetized (force-free) limit. The plasma turbulence ubiquitously generates current sheets, which act as locations where magnetic energy dissipates. We show that current sheets form as a natural result of nonlinear interactions between counter-propagating Alfvén waves. These current sheets form due to the compression of elongated eddies, driven by the shear induced by growing higher order modes, and undergo a thinning process until they break-up into small-scale turbulent structures. We explore the formation of {current sheets} both in overlapping waves and in localized wave packet collisions. The relativistic interaction of localized Alfvén waves induces both Alfvén waves and fast waves and efficiently mediates the conversion and dissipation of electromagnetic energy in astrophysical systems. Plasma energization through reconnection in current sheets emerging during the interaction of Alfvén waves can potentially explain X-ray emission in black hole accretion coronae and neutron star magnetospheres. △ Less

Submitted 25 October, 2021; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: Accepted in JPP, accompanied by Paper I: Weak Alfvénic turbulence in relativistic plasmas I: dynamical equations and basic dynamics of interacting resonant triads

Similarities and Differences in Accretion Flow Properties between GRS 1915+105 and IGR J17091-3624: a Case Study

Authors: Anuvab Banerjee, Ayan Bhattacharjee, Dipak Debnath, Sandip K. Chakrabarti

Abstract: We perform a comparative spectro-temporal analysis on the variability classes of GRS 1915+105 and IGR J17091-3624 to draw inferences regarding the underlying accretion flow mechanism. The $ν$, as well as C2 class Rossi X-Ray Timing Explorer observation, have been considered for analysis. We investigate the intensity variation of the source in different energy domains that correspond to different c… ▽ More We perform a comparative spectro-temporal analysis on the variability classes of GRS 1915+105 and IGR J17091-3624 to draw inferences regarding the underlying accretion flow mechanism. The $ν$, as well as C2 class Rossi X-Ray Timing Explorer observation, have been considered for analysis. We investigate the intensity variation of the source in different energy domains that correspond to different components of the accretion flow and infer the relative dominance of these flow components during the dip/flare events. We correlate the dependence of the dynamic photon index ($Θ$) with intensities in different energy bands and comment on the transition of the source to hard/soft phases during soft dips/flares. We also report the presence of sharp QPOs at \sim7.1 Hz corresponding to both softer and harder domain in the case of $ν$ variability class of GRS 1915+105 and discuss the possible accretion flow configuration it suggests. Sharp QPO around \sim20 mHz is observed in $ν$ and C2 classes of IGR J17091-3624 in low and mid energy band (2.0-6.0 keV and 6.0-15.0 keV), but remains undetected in high energy (15.0-60.0 keV). The 2.5-25.0 keV background-subtracted spectra have also been fitted with TCAF along with a Compton reflection component. A plausible accretion flow mechanism in order to explain the observed variability has been proposed. △ Less

Submitted 19 January, 2022; v1 submitted 28 January, 2021; originally announced January 2021.

Comments: 23 Pages, 17 Figures, 3 Tables, Accepted for publication in Advances in Space Research on Jan 11, 2022

Two-Component Advective Flows around Neutron Stars

Authors: Ayan Bhattacharjee, Sandip K. Chakrabarti

Abstract: The fundamental difference between accretion around black holes and neutron stars is the inner boundary condition, which affects the behavior of matter very close to the compact objects. This leads to formation of additional shocks and boundary layers for neutron stars. Previous studies on the formation of such boundary layers focused on Keplerian flows that reached the surface of the star, either… ▽ More The fundamental difference between accretion around black holes and neutron stars is the inner boundary condition, which affects the behavior of matter very close to the compact objects. This leads to formation of additional shocks and boundary layers for neutron stars. Previous studies on the formation of such boundary layers focused on Keplerian flows that reached the surface of the star, either directly or through the formation of a transition layer. However, behavior of sub-Keplerian matter near the surface of a neutron star has not been studied in detail. Here, we study the effect of viscosity, in presence of cooling, on the sub-Keplerian flows around neutron stars, using Smoothed Particle Hydrodynamics. Our time-dependent study shows that multiple shocks, transition and boundary layers form in such type of accretion, when viscosity is significant, and one or more layers are absent when the viscosity is moderate. These flows are particularly of interest for the wind dominated systems such as Cir X-1. We also report the formation of a generalized flow configuration, Two-Component Advective Flow, for the first time. △ Less

Submitted 5 January, 2021; originally announced January 2021.

Comments: 4 Pages, 5 Figures, 1 Table. Submitted to Journal of Korean Astronomical Society. arXiv admin note: substantial text overlap with arXiv:1906.10630

A Unified Spectral Model for Accreting Neutron Stars

Authors: Ayan Bhattacharjee, Sandip K. Chakrabarti

Abstract: Recent studies show that accretion flows around weakly magnetic Neutron Stars undergo multiple shocks, before reaching the surface of the star, which contribute to the spectral and timing variabilities observed in the X-Rays. Here, we report, for the first time, the spectral properties of a unified model of shocked accretion flows around Neutron Stars, based on the Two-Component Advective Flow par… ▽ More Recent studies show that accretion flows around weakly magnetic Neutron Stars undergo multiple shocks, before reaching the surface of the star, which contribute to the spectral and timing variabilities observed in the X-Rays. Here, we report, for the first time, the spectral properties of a unified model of shocked accretion flows around Neutron Stars, based on the Two-Component Advective Flow paradigm. We compare our theoretical results with the X-Ray spectral features of Z and Atoll sources, across different states. We also fit RXTE/PCA spectra of Sco X-1 and 4U 1705-44 to show the potential application of this new model. △ Less

Submitted 28 December, 2020; originally announced December 2020.

Comments: 4 pages, 5 figures, 3 tables, submitted to JKAS

BH Mass, Jet and Accretion Disk Connection: An Analysis of Radio-loud and Radio-quiet Quasars

Authors: Avinanda Chakraborty, Anirban Bhattacharjee

Abstract: Surveys have shown radio-loud (RL) quasars constitute $10 \%$-$15 \%$ of the total quasar population. However, it is unknown if the radio-loud fraction (RL quasars/Total quasars) remains consistent among different parameter spaces. This study shows that radio-loud fraction increases for increasing full width half maximum (FWHM) velocity of the H$β$ and MgII broad emission line. Our data has been o… ▽ More Surveys have shown radio-loud (RL) quasars constitute $10 \%$-$15 \%$ of the total quasar population. However, it is unknown if the radio-loud fraction (RL quasars/Total quasars) remains consistent among different parameter spaces. This study shows that radio-loud fraction increases for increasing full width half maximum (FWHM) velocity of the H$β$ and MgII broad emission line. Our data has been obtained from Shen et al. (2011) catalogue. To investigate the reason, in this preliminary study we analyse various properties like bolometric luminosity, optical continuum luminosity, black hole (BH) mass and accretion rate of RL quasars and RQ quasars sample which have FWHM greater than 15000km/s (High broad line). From the distributions we can conclude for all the properties in high broad line, RL quasars are having higher values than RQ quasars. We have predicted RL quasars are intrinsically brighter than RQ quasars and also predicted BH mass-jet connection and accretion disk-jet connection from our results but to conclude anything more analysis is needed. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: Accepted in Astronomische Nachrichten (AN) Journal. IWARA 2020 conference proceeding

Observations and modeling of the onset of fast reconnection in the solar transition region

Authors: L. -J. Guo, B. De Pontieu, Y. -M. Huang, H. Peter, A. Bhattacharjee

Abstract: Magnetic reconnection is a fundamental plasma process that plays a critical role not only in energy release in the solar atmosphere, but also in fusion, astrophysical, and other space plasma environments. One of the challenges in explaining solar observations in which reconnection is thought to play a critical role is to account for the transition of the dynamics from a slow quasi-continuous phase… ▽ More Magnetic reconnection is a fundamental plasma process that plays a critical role not only in energy release in the solar atmosphere, but also in fusion, astrophysical, and other space plasma environments. One of the challenges in explaining solar observations in which reconnection is thought to play a critical role is to account for the transition of the dynamics from a slow quasi-continuous phase to a fast and impulsive energetic burst of much shorter duration. Despite the theoretical progress in identifying mechanisms that might lead to rapid onset, a lack of observations of this transition has left models poorly constrained. High-resolution spectroscopic observations from NASA's Interface Region Imaging Spectrograph (IRIS) now reveal tell-tale signatures of the abrupt transition of reconnection from a slow phase to a fast, impulsive phase during UV bursts or explosive events in the Sun's atmosphere. Our observations are consistent with numerical simulations of the plasmoid instability, and provide evidence for the onset of fast reconnection mediated by plasmoids and new opportunities for remote-sensing diagnostics of reconnection mechanisms on the Sun. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Comments: 9 pages, 5 figures, accepted for publication in Astrophysical Journal

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Authors: H. Ji, J. Karpen, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, A. Bhattacharjee, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, B. Chen, L. -J. Chen, Y. Chen, A. Chien, L. Comisso, D. Craig, J. Dahlin, W. Daughton, E. DeLuca , et al. (83 additional authors not shown)

Abstract: Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagen… ▽ More Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 4 pages (including a title page), white paper submitted to Helio2050 workshop at https://www.hou.usra.edu/meetings/helio2050/. arXiv admin note: substantial text overlap with arXiv:2004.00079

Small-Scale Dynamo in Stably Stratified Turbulence

Authors: Valentin Skoutnev, Jonathan Squire, Amitava Bhattacharjee

Abstract: We present numerical investigations into three principal properties of the small-scale dynamo in stably stratified turbulence: the onset criterion, the growth rate, and the nature of the magnetic field anisotropy in the kinematic regime. The results suggest that all three dynamo properties are controlled by the scale separation between the Ozmidov scale and the viscous or resistive scale. In addit… ▽ More We present numerical investigations into three principal properties of the small-scale dynamo in stably stratified turbulence: the onset criterion, the growth rate, and the nature of the magnetic field anisotropy in the kinematic regime. The results suggest that all three dynamo properties are controlled by the scale separation between the Ozmidov scale and the viscous or resistive scale. In addition to the critical magnetic Reynolds number, this allows for the definition of critical buoyancy and magnetic buoyancy Reynolds numbers for stratified small-scale dynamo onset in the high and low magnetic Prandtl number regimes, respectively. The presence of a small-scale dynamo in stellar radiative zones could affect dynamics through resulting Maxwell stresses and/or influence large-scale dynamo mechanisms in regions of differential rotation. Taking the solar radiative zone as a representative example and applying the onset criterion, we find that the stratification is strong enough to make the small-scale dynamo marginally active in the stably stratified turbulence of the solar tachocline. △ Less

Submitted 3 February, 2021; v1 submitted 3 August, 2020; originally announced August 2020.

Journal ref: ApJ 906 61 (2021)

Accretion Flow Properties of GRS 1915+105 During Its $θ$ Class Using AstroSat Data

Authors: Anuvab Banerjee, Ayan Bhattacharjee, Debjit Chatterjee, Dipak Debnath, Sandip Kumar Chakrabarti, Tilak Katoch, H. M. Antia

Abstract: The Galactic microquasar GRS 1915+105 shows rich variability which is categorized into different classes. In this paper, we report temporal and spectral analysis of GRS 1915+105 to study the properties of the accretion flow when the light curve is showing $θ$ class variability. For this purpose, we use the LAXPC data from the Target of Opportunity observations of India's first multi-wavelength ast… ▽ More The Galactic microquasar GRS 1915+105 shows rich variability which is categorized into different classes. In this paper, we report temporal and spectral analysis of GRS 1915+105 to study the properties of the accretion flow when the light curve is showing $θ$ class variability. For this purpose, we use the LAXPC data from the Target of Opportunity observations of India's first multi-wavelength astronomy satellite AstroSat. The $θ$ class is marked by the recurrent appearance of U-shaped regions in the light curve, where the photon count rate first decreases rapidly and then increases slowly. For our analysis, we use U-shaped regions of first two orbits (02345 and 02346) on 2016 March 04. In both of the cases, the dynamic Power Density Spectra (PDS) showed significant power at around $4-5$ Hz, suggesting the presence of a low-frequency Quasi-Periodic Oscillation (QPO) around that frequency interval. The QPO frequency is found to increase with time when the enhancement of the energy flux also takes place. From the evolution of the spectra, we determine the evolution of the accretion flow parameters in both of these observations. Fitting the spectra with the transonic flow solution based Two-Component Advective Flow (TCAF) model in $4-25$ keV energy band shows that the Keplerian disk accretion rate rises with the rise in the radiation intensity, while the location of the centrifugal pressure driven shock front decreases. In both these data, a gradual increment of power-law photon index with intensity is observed suggesting the progressive softening of the source. △ Less

Submitted 5 July, 2021; v1 submitted 10 July, 2020; originally announced July 2020.

Comments: 16 pages, 6 figures, 2 tables. Accepted in the Astrophysical Journal on May 5, 2021

AstroSat Observation of Non-Resonant Type-C QPOs in MAXI J1535-571

Authors: D. Chatterjee, D. Debnath, A. Jana, J. -R. Shang, S. K. Chakrabarti, H. -K. Chang, A. Banerjee, A. Bhattacharjee, K. Chatterjee, R. Bhowmick, S. K. Nath

Abstract: Galactic transient black hole candidate (BHC) MAXI J1535-571 was discovered on 2017 September 02 simultaneously by {\it MAXI}/GSC and {\it Swift}/BAT instruments. It has also been observed by India's first multi-wavelength astronomy-mission satellite {\it AstroSat}, during the rising phase of its 2017-18 outburst. We make both the spectral and the temporal analysis of the source during 2017 Septem… ▽ More Galactic transient black hole candidate (BHC) MAXI J1535-571 was discovered on 2017 September 02 simultaneously by {\it MAXI}/GSC and {\it Swift}/BAT instruments. It has also been observed by India's first multi-wavelength astronomy-mission satellite {\it AstroSat}, during the rising phase of its 2017-18 outburst. We make both the spectral and the temporal analysis of the source during 2017 September 12-17 using data of {\it AstroSat}'s Large Area X-ray Proportional Counter (LAXPC) in the energy range of $3-40$~keV to infer the accretion flow properties of the source. Spectral analysis is done with the physical two-component advective flow (TCAF) solution-based {\it fits} file. From the nature of the variation of the TCAF model fitted physical flow parameters, we conclude and confirm that the source was in the intermediate spectral state during our analysis period. We observe sharp type-C quasi-periodic oscillations (QPOs) in the frequency range of $\sim 1.75-2.81$~Hz. For a better understanding of the nature and evolution of these type-C QPOs, a dynamic study of the power density spectra is done. We also investigate the origin of these QPOs from the shock oscillation model. We find that non-satisfaction of Rankine-Hugoniot conditions for non-dissipative shocks and not their resonance oscillations is the cause of the observed type-C QPOs. △ Less

Submitted 11 August, 2021; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: 15 pages, 6 figures, 3 tables (accepted for publication in Ap&SS)

Spectroscopic properties of radio-loud and radio-quiet quasars

Authors: Avinanda Chakraborty, Anirban Bhattacharjee

Abstract: Surveys have shown radio-loud (RL) quasars constitute 10%-15% of the total quasar population and rest are radio-quiet (RQ). However, it is unknown if this radio-loud fraction (RLF) remains consistent among different parameter spaces. This study shows that RLF increases for increasing full width half maximum (FWHM) velocity of the H$β$ broad emission line (z < 0.75). To analyse the reason, we compa… ▽ More Surveys have shown radio-loud (RL) quasars constitute 10%-15% of the total quasar population and rest are radio-quiet (RQ). However, it is unknown if this radio-loud fraction (RLF) remains consistent among different parameter spaces. This study shows that RLF increases for increasing full width half maximum (FWHM) velocity of the H$β$ broad emission line (z < 0.75). To analyse the reason, we compared bolometric luminosity of RL and RQ quasars sample which have FWHM of H$β$ broad emission line greater than 15000km/s (High Broad Line or HBL) with which have FWHM of H$β$ emission line less than 2500km/s (Low Broad Line or LBL). From the distributions we can conclude for the HBL, RQ and RL quasars are peaking separately and RL quasars are having higher values whereas for the LBL the peaks are almost indistinguishable. We predicted selection effects could be the possible reason but to conclude anything more analysis is needed. Then we compared our result with Wills & Browne (1986) and have shown that some objects from our sample do not follow the pattern of the logR vs FWHM plot where R is the ratio of 5 GHz radio core flux density with the extended radio lobe flux density. △ Less

Submitted 7 April, 2020; originally announced April 2020.

Comments: Proceedings paper of the IAU symposium "Nuclear Activity in Galaxies Across Cosmic Time" (Ethiopia) accepted to be published under the Cambridge University Press, eds. M. Povic, P. Marziani, J. Masegosa, H. Netzer, S. H. Negu, and S. B. Tessema

Journal ref: Proc. IAU 15 (2019) 72-76

Scalings Pertaining to Current Sheet Disruption Mediated by the Plasmoid Instability

Authors: Yi-Min Huang, Luca Comisso, Amitava Bhattacharjee

Abstract: Analytic scaling relations are derived for a phenomenological model of the plasmoid instability in an evolving current sheet, including the effects of reconnection outflow. Two scenarios are considered, where the plasmoid instability can be triggered either by an injected initial perturbation or by the natural noise of the system (here referred to as the system noise). The two scenarios lead to di… ▽ More Analytic scaling relations are derived for a phenomenological model of the plasmoid instability in an evolving current sheet, including the effects of reconnection outflow. Two scenarios are considered, where the plasmoid instability can be triggered either by an injected initial perturbation or by the natural noise of the system (here referred to as the system noise). The two scenarios lead to different scaling relations because the initial noise decays when the linear growth of the plasmoid instability is not sufficiently fast to overcome the advection loss caused by the reconnection outflow, whereas the system noise represents the lowest level of fluctuations in the system. The leading order approximation for the current sheet width at disruption takes the form of a power-law multiplied by a logarithmic factor, and from that, the scaling relations for the wavenumber and the linear growth rate of the dominant mode are obtained. When the effects of the outflow are neglected, the scaling relations agree, up to the leading order approximation, with previously derived scaling relations based on a principle of least time. The analytic scaling relations are validated with numerical solutions of the model. △ Less

Submitted 16 September, 2019; v1 submitted 6 September, 2019; originally announced September 2019.

Comments: Accepted for publication in Physics of Plasmas

Timing Properties of Shocked Accretion Flows around Neutron Stars -- II. Viscous Disks and Boundary Layers

Authors: Ayan Bhattacharjee, Sandip K. Chakrabarti

Abstract: We use Smoothed Particle Hydrodynamics to study viscous accretion flows around a weakly magnetic neutron star. We show the formation of multiple ``boundary" layers in presence of both cooling and viscosity. We find that with the introduction of a small viscosity in a sub-Keplerian flow, much like the wind accretion in HMXBs such as Cir X-1, only a single Normal Boundary Layer (NBOL) forms to adjus… ▽ More We use Smoothed Particle Hydrodynamics to study viscous accretion flows around a weakly magnetic neutron star. We show the formation of multiple ``boundary" layers in presence of both cooling and viscosity. We find that with the introduction of a small viscosity in a sub-Keplerian flow, much like the wind accretion in HMXBs such as Cir X-1, only a single Normal Boundary Layer (NBOL) forms to adjust the rotational velocity component. With the increase of viscosity, the region extends radially and beyond some critical value, a RAdiative KEplerian Disk/layer (RAKED) forms between the sub-Keplerian flow and the NBOL. When viscosity is increased further only NBOL and RAKED remain. In all such cases, the CENtrifugal pressure dominated BOundary Layer (CENBOL) is formed, away from the star, as in the case of black holes. This is the first self-consistent study where such a transition from sub-Keplerian flows has been reported for neutron stars. We also identify the connection between accretion and ejection of matter, following the Two-Component Advective Flow for black holes, for neutron stars. The results are crucial in the understanding of the formation of disks, boundary layers and outflows in wind dominated neutron star systems. △ Less

Submitted 25 June, 2019; originally announced June 2019.

Comments: 5 Pages, 1 Figure. Submitted in The Astrophysical Journal on June 22, 2019

Spectral Analysis of $χ$ Class Data of GRS 1915+105 Using TCAF Solution

Authors: Anuvab Banerjee, Ayan Bhattacharjee, Dipak Debnath, Sandip K. Chakrabarti

Abstract: The class variable source GRS 1915+105 exhibits a wide range of time variabilities in time scales of a few seconds to a few days. Depending on the count rates in different energy bands and the nature of the conventional color-color diagram, the variabilities were classified into sixteen classes, which were later sequenced in ascending order of Comptonization Efficiency (CE). The CE is the ratio of… ▽ More The class variable source GRS 1915+105 exhibits a wide range of time variabilities in time scales of a few seconds to a few days. Depending on the count rates in different energy bands and the nature of the conventional color-color diagram, the variabilities were classified into sixteen classes, which were later sequenced in ascending order of Comptonization Efficiency (CE). The CE is the ratio of the phenomenological power-law and disk blackbody model fitted spectral photon fluxes. For a better understanding of the flow dynamics in one of the variability classes, namely $χ$, we fitted spectra with the physical two-component advective flow (TCAF) model. In $χ$ class, X-ray flux is steady with no significant variation, however, various $χ$ subclasses are observed at different X-ray flux and variations of count rates across different $χ$ subclasses must be linked to the variation of flow parameters such as the accretion rates, be it the Keplerian disc rate and/or the low angular momentum halo rate. We find that in the $χ_{2,4}$ classes, which are reportedly devoid of significant outflows, the spectra could be fitted well using TCAF solution alone. In the $χ_{1,3}$ class, which are always linked with outflows, a cutoff power-law model is needed in addition to the TCAF solution. At the same time, the normalization required by this model along with the variation of photon index and exponential roll-off factor provides us with the information on the relative dominance of the outflow in the latter two classes. TCAF fit also provides us with the size and location of the Compton cloud along with its optical depth. Thus by fitting with TCAF, a physical understanding of the flow geometry in different $χ$ classes of GRS 1915+105 has been obtained. △ Less

Submitted 13 June, 2020; v1 submitted 4 May, 2019; originally announced May 2019.

Comments: 13 pages, 6 figures, 3 tables, Accepted for publication in Research in Astronomy and Astrophysics on June 4, 2020 (Original version submitted on February 20, 2020)

Constraining Mass of Cygnus X-1 from Analysis of the Hard State Spectral Data using TCAF Solution

Authors: Indrani Banerjee, Ayan Bhattacharjee, Anuvab Banerjee, Dipak Debnath, Sandip K. Chakrabarti

Abstract: The galactic black hole candidate Cygnus X-1, one of the brightest sources in the sky, is the first ever black hole candidate to be discovered. Despite being a very well-studied object due to its persistent brightness in X-rays, there has been much difficulty in determining its mass since its discovery. In spite of continuous efforts since the early nineteen seventies, there is yet no concensus on… ▽ More The galactic black hole candidate Cygnus X-1, one of the brightest sources in the sky, is the first ever black hole candidate to be discovered. Despite being a very well-studied object due to its persistent brightness in X-rays, there has been much difficulty in determining its mass since its discovery. In spite of continuous efforts since the early nineteen seventies, there is yet no concensus on its mass for nearly four decades. The uncertainties in the mass measurements are due to the high degree of error involved in the measurement of its distance. In 2011, Orosz et al. constrained the mass of this object to be M = 14.8 +/- 1.0 M_Sun using dynamical methods. In this paper, we obtained the mass of Cygnus X-1, using a completely independent method, namely, carrying out the spectral analysis using Two Component Advective Flow (TCAF) solution based fits file and the archival data of RXTE PCA instrument. Our result does not require the distance of the source or the information about the companion. Each spectral fit with the TCAF gives one best fitted mass. Averaging fitted masses of Cygnus X-1 over a span of five months of observation during its persistent hard phase, mass of the source comes out to be M_avg = 14.20 +/- 0.36 M_Sun, which is consistent with the dynamically estimated mass. △ Less

Submitted 25 April, 2019; originally announced April 2019.

Comments: 4 pages, 2 figures, 1 table

Global Ten-Moment Multifluid Simulations of the Solar Wind Interaction with Mercury: From the Planetary Conducting Core to the Dynamic Magnetosphere

Authors: Chuanfei Dong, Liang Wang, Ammar Hakim, Amitava Bhattacharjee, James A. Slavin, Gina A. DiBraccio, Kai Germaschewski

Abstract: For the first time, we explore the tightly coupled interior-magnetosphere system of Mercury by employing a three-dimensional ten-moment multifluid model. This novel fluid model incorporates the non-ideal effects including the Hall effect, inertia, and tensorial pressures that are critical for collisionless magnetic reconnection; therefore, it is particularly well suited for investigating… ▽ More For the first time, we explore the tightly coupled interior-magnetosphere system of Mercury by employing a three-dimensional ten-moment multifluid model. This novel fluid model incorporates the non-ideal effects including the Hall effect, inertia, and tensorial pressures that are critical for collisionless magnetic reconnection; therefore, it is particularly well suited for investigating $collisionless$ magnetic reconnection in Mercury's magnetotail and at the planet's magnetopause. The model is able to reproduce the observed magnetic field vectors, field-aligned currents, and cross-tail current sheet asymmetry (beyond the MHD approach) and the simulation results are in good agreement with spacecraft observations. We also study the magnetospheric response of Mercury to a hypothetical extreme event with an enhanced solar wind dynamic pressure, which demonstrates the significance of induction effects resulting from the electromagnetically-coupled interior. More interestingly, plasmoids (or flux ropes) are formed in Mercury's magnetotail during the event, indicating the highly dynamic nature of Mercury's magnetosphere. △ Less

Submitted 21 September, 2019; v1 submitted 4 April, 2019; originally announced April 2019.

Comments: Geophysical Research Letters, in press, 17 pages, 4 (fancy) figures

Journal ref: Geophysical Research Letters, 46, 11584-11596, 2019

An extended MHD study of the 16 October 2015 MMS diffusion region crossing

Authors: J. M. TenBarge, J. Ng, J. Juno, L. Wang, A. H. Hakim, A. Bhattacharjee

Abstract: The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a… ▽ More The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multi-fluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. We focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: 17 pages, 13 figures, submitted to JGR - Space Physics

Timing Properties of Shocked Accretion Flows around Neutron Stars in Presence of Cooling

Authors: Ayan Bhattacharjee, Sandip K. Chakrabarti

Abstract: We carry out the first robust numerical simulation of accretion flows on a weakly magnetized neutron star using Smoothed Particle Hydrodynamics (SPH). We follow the Two-Component Advective Flow (TCAF) paradigm for black holes, and focus only on the advective component for the case of a neutron star. This low viscosity sub-Keplerian flow will create a normal boundary layer (or, NBOL) right on the s… ▽ More We carry out the first robust numerical simulation of accretion flows on a weakly magnetized neutron star using Smoothed Particle Hydrodynamics (SPH). We follow the Two-Component Advective Flow (TCAF) paradigm for black holes, and focus only on the advective component for the case of a neutron star. This low viscosity sub-Keplerian flow will create a normal boundary layer (or, NBOL) right on the star surface in addition to the centrifugal pressure supported boundary layer (or, CENBOL) present in a black hole accretion. These density jumps could give rise to standing or oscillating shock fronts. During a hard spectral state, the incoming flow has a negligible viscosity causing more sub-Keplerian component as compared to the Keplerian disc component. We show that our simulation of flows with a cooling and a negligible viscosity produces precisely two shocks and strong supersonic winds from these boundary layers. We find that the specific angular momentum of matter dictates the locations and the nature of oscillations of these shocks. For low angular momentum flows, the radial oscillation appears to be preferred. For flows with higher angular momentum, the vertical oscillation appears to become dominant. In all the cases, asymmetries w.r.t. the Z=0 plane are seen and instabilities set in due to the interaction of inflow and outgoing strong winds. Our results capture both the low and high-frequency quasi-periodic oscillations without invoking magnetic fields or any precession mechanism. Most importantly, these solutions directly corroborate observed features of wind dominated high-mass X-ray binaries, such as Cir X-1. △ Less

Submitted 14 February, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: 13 pages (two-column), 4 figures (16 files), 1 table. Accepted for publication in the Astrophysical Journal on February 13, 2019

Generalized Flows around Neutron Stars

Authors: Ayan Bhattacharjee

Abstract: In this chapter, we present a brief and non-exhaustive review of the developments of theoretical models for accretion flows around neutron stars. A somewhat chronological summary of crucial observations and modelling of timing and spectral properties are given in sections 2 and 3. In section 4, we argue why and how the Two-Component Advective Flow (TCAF) solution can be applied to the cases of neu… ▽ More In this chapter, we present a brief and non-exhaustive review of the developments of theoretical models for accretion flows around neutron stars. A somewhat chronological summary of crucial observations and modelling of timing and spectral properties are given in sections 2 and 3. In section 4, we argue why and how the Two-Component Advective Flow (TCAF) solution can be applied to the cases of neutron stars when suitable modifications are made for the NSs. We showcase some of our findings from Monte Carlo and Smoothed Particle Hydrodynamic simulations which further strengthens the points raised in section 4. In summary, we remark on the possibility of future works using TCAF for both weakly magnetic and magnetic Neutron Stars. △ Less

Submitted 4 January, 2019; originally announced January 2019.

Comments: 15 pages, 7 figures. arXiv admin note: text overlap with arXiv:1901.00842

Journal ref: Exploring the Universe: From Near Space to Extra-Galactic, Astrophysics and Space Science Proceedings, Volume 53. ISBN 978-3-319-94606-1. Springer International Publishing AG, part of Springer Nature, 2018, p. 93

Monte Carlo Simulations of Thermal Comptonization Process in a Two Component Advective Flow around a Neutron Star

Authors: Ayan Bhattacharjee, Sandip K. Chakrabarti

Abstract: We explore spectral properties of a Two-Component Advective Flow (TCAF) around a neutron star. We compute the effects of thermal Comptonization of soft photons emitted from a Keplerian disc and the boundary layer of the neutron star by the post-shock region of a sub-Keplerian flow, formed due to the centrifugal barrier. The shock location $X_s$ is also the inner edge of the Keplerian disc. We comp… ▽ More We explore spectral properties of a Two-Component Advective Flow (TCAF) around a neutron star. We compute the effects of thermal Comptonization of soft photons emitted from a Keplerian disc and the boundary layer of the neutron star by the post-shock region of a sub-Keplerian flow, formed due to the centrifugal barrier. The shock location $X_s$ is also the inner edge of the Keplerian disc. We compute a series of realistic spectra assuming a set of electron temperatures of the post-shock region $T_{CE}$, the temperature of the normal boundary layer (NBOL) $T_{NS}$ of the neutron star and the shock location $X_s$. These parameters depend on the disc and halo accretion rates ($\dot{m_d}$ and $\dot{m_h}$, respectively) which control the resultant spectra. We find that the spectrum becomes harder when $\dot{m}_h$ is increased. The spectrum is controlled strongly by $T_{NS}$ due to its proximity to the Comptonizing cloud since photons emitted from the NBOL cool down the post-shock region very effectively. We also show the evidence of spectral hardening as the inclination angle of the disc is increased. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: 28 pages, 14 figures, 1 table

Journal ref: MNRAS 472 (2017) 1361-1371

The 2004 Outburst of BHC H1743-322: Analysis of spectral and timing properties using the TCAF Solution

Authors: Ayan Bhattacharjee, Indrani Banerjee, Anuvab Banerjee, Dipak Debnath, Sandip K Chakrabarti

Abstract: The black hole transient H1743-322 exhibited several outbursts with temporal and spectral variability since its discovery in 1977. These outbursts occur at a quasi-regular recurrence period of around $0.5-2$ years, since its rediscovery in March 2003. We investigate accretion flow dynamics around the Low Mass X-ray Binary H1743-322 during its 2004 outburst using the RXTE/PCA archival data. We use… ▽ More The black hole transient H1743-322 exhibited several outbursts with temporal and spectral variability since its discovery in 1977. These outbursts occur at a quasi-regular recurrence period of around $0.5-2$ years, since its rediscovery in March 2003. We investigate accretion flow dynamics around the Low Mass X-ray Binary H1743-322 during its 2004 outburst using the RXTE/PCA archival data. We use Two Component Advective Flow (TCAF) solution to analyse the spectral data. From the fits with TCAF solution, we obtain day to day variation of physical accretion rates of Keplerian and sub-Keplerian components, size of the Compton cloud and its other properties. Analysis of the spectral properties of the 2004 outburst by keeping fitted normalization to be in a narrow range and its timing properties in terms of the presence and absence of QPOs, enable us to constrain the mass of the black hole in a range of $10.31 M_{\odot} - 14.07 M_{\odot}$ which is consistent with other estimates reported in the literature. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: 24 pages, 7 figures, 2 tables (Published in MNRAS)

Journal ref: MNRAS 466 (2017) 1372-1381

X-ray Surface Brightness Profiles of Optically Selected Active Galactic Nuclei: Comparison with X-ray AGN

Authors: Sagnick Mukherjee, Anirban Bhattacharjee, Suchetana Chatterjee, Jeffrey A. Newman, Renbin Yan

Abstract: We use data from the All Wavelength Extended Groth Strip International Survey to construct stacked X-ray maps of optically bright active galaxies (AGN) and an associated control sample of galaxies at high redshift (z less than 0.6). From our analysis of the surface brightness profiles obtained from these X-ray maps, we find evidence of feedback from the active nuclei. We find that excluding galaxi… ▽ More We use data from the All Wavelength Extended Groth Strip International Survey to construct stacked X-ray maps of optically bright active galaxies (AGN) and an associated control sample of galaxies at high redshift (z less than 0.6). From our analysis of the surface brightness profiles obtained from these X-ray maps, we find evidence of feedback from the active nuclei. We find that excluding galaxies and AGN, residing in group environments, from our samples enhances the significance of our detection. Our results support the tentative findings of Chatterjee et al. who use X-ray selected AGN for their analysis. We discuss the implications of these results in the context of quantifying AGN feedback and show that the current method can be used to extract X-ray source population in high redshift galaxies. △ Less

Submitted 7 February, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: The Astrophysical Journal, Volume 872, Number 1, Page- 35

Velocity-resolved reverberation mapping of five bright Seyfert 1 galaxies

Authors: G. De Rosa, M. M. Fausnaugh, C. J. Grier, B. M. Peterson, K. D. Denney, Keith Horne, M. C. Bentz, S. Ciroi, E. Dalla Bonta`, M. D. Joner, S. Kaspi, C. S. Kochanek, R. W. Pogge, S. G. Sergeev, M. Vestergaard, S. M. Adams, J. Antognini, C. Araya Salvo, E. Armstrong, J. Bae, A. J. Barth, T. G. Beatty, A. Bhattacharjee, G. A. Borman, T. A. Boroson , et al. (77 additional authors not shown)

Abstract: We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line reg… ▽ More We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007. △ Less

Submitted 3 August, 2018; v1 submitted 12 July, 2018; originally announced July 2018.

Comments: Accepted for publication on ApJ; 32 pages, 16 figures, 10 tables

Role of the Plasmoid Instability in Magnetohydrodynamic Turbulence

Authors: Chuanfei Dong, Liang Wang, Yi-Min Huang, Luca Comisso, Amitava Bhattacharjee

Abstract: The plasmoid instability in evolving current sheets has been widely studied due to its effects on the disruption of current sheets, the formation of plasmoids, and the resultant fast magnetic reconnection. In this Letter, we study the role of the plasmoid instability in two-dimensional magnetohydrodynamic (MHD) turbulence by means of high-resolution direct numerical simulations. At sufficiently la… ▽ More The plasmoid instability in evolving current sheets has been widely studied due to its effects on the disruption of current sheets, the formation of plasmoids, and the resultant fast magnetic reconnection. In this Letter, we study the role of the plasmoid instability in two-dimensional magnetohydrodynamic (MHD) turbulence by means of high-resolution direct numerical simulations. At sufficiently large magnetic Reynolds number ($R_m=10^6$), the combined effects of dynamic alignment and turbulent intermittency lead to a copious formation of plasmoids in a multitude of intense current sheets. The disruption of current sheet structures facilitates the energy cascade towards small scales, leading to the breaking and steepening of the energy spectrum. In the plasmoid-mediated regime, the energy spectrum displays a scaling that is close to the spectral index $-2.2$ as proposed by recent analytic theories. We also demonstrate that the scale-dependent dynamic alignment exists in 2D MHD turbulence and the corresponding slope of the alignment angle is close to 0.25. △ Less

Submitted 25 September, 2018; v1 submitted 19 April, 2018; originally announced April 2018.

Comments: 6 pages, 5 figures, accepted in Physical Review Letters

Journal ref: Phys. Rev. Lett. 121, 165101 (2018)

Magnetohydrodynamic Turbulence in the Plasmoid-Mediated Regime

Authors: L. Comisso, Y. -M. Huang, M. Lingam, E. Hirvijoki, A. Bhattacharjee

Abstract: Magnetohydrodynamic turbulence and magnetic reconnection are ubiquitous in astrophysical environments. In most situations, these processes do not occur in isolation, but interact with each other. This renders a comprehensive theory of these processes highly challenging. Here, we propose a theory of magnetohydrodynamic turbulence driven at large scale that self-consistently accounts for the mutual… ▽ More Magnetohydrodynamic turbulence and magnetic reconnection are ubiquitous in astrophysical environments. In most situations, these processes do not occur in isolation, but interact with each other. This renders a comprehensive theory of these processes highly challenging. Here, we propose a theory of magnetohydrodynamic turbulence driven at large scale that self-consistently accounts for the mutual interplay with magnetic reconnection occurring at smaller scales. Magnetic reconnection produces plasmoids that grow from turbulence-generated noise and eventually disrupt the sheet-like structures in which they are born. The disruption of these structures leads to a modification of the turbulent energy cascade, which, in turn, exerts a feedback effect on the plasmoid formation via the turbulence-generated noise. The energy spectrum in this plasmoid-mediated range steepens relative to the standard inertial range and does not follow a simple power law. As a result of the complex interplay between turbulence and reconnection, we also find that the length scale which marks the beginning of the plasmoid-mediated range and the dissipation length scale do not obey true power laws. The transitional magnetic Reynolds number above which the plasmoid formation becomes statistically significant enough to affect the turbulent cascade is fairly modest, implying that plasmoids are expected to modify the turbulent path to dissipation in many astrophysical systems. △ Less

Submitted 26 March, 2018; v1 submitted 6 February, 2018; originally announced February 2018.

Comments: Added an important reference: Carbone, Veltri, Mangeney, "Coherent structure formation and magnetic field line reconnection in magnetohydrodynamic turbulence", Phys. Fluids A (1990)

Journal ref: Astrophys. J. 854, 103 (2018)

Continuum Reverberation Mapping of the Accretion Disks in Two Seyfert 1 Galaxies

Authors: M. M. Fausnaugh, D. A. Starkey, Keith Horne, C. S. Kochanek, B. M. Peterson, M. C. Bentz, K. D. Denney, C. J. Grier, D. Grupe, R. W. Pogge, G. DeRosa, S. M. Adams, A. J. Barth, Thomas G. Beatty, A. Bhattacharjee, G. A. Borman, T. A. Boroson, M. C. Bottorff, Jacob E. Brown, Jonathan S. Brown, M. S. Brotherton, C. T. Coker, S. M. Crawford, K. V. Croxall, Sarah Eftekharzadeh , et al. (47 additional authors not shown)

Abstract: We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC\,2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find… ▽ More We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC\,2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of $τ\propto λ^{4/3}$. However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC\,2617. These differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC\,2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC\,2617 precedes the UV/optical variability, the long 2.6 day lag is problematic for coronal reprocessing models. △ Less

Submitted 29 January, 2018; originally announced January 2018.

Comments: Accepted to ApJ, please send comments to faus@mit.edu. 24 pages, 8 figures

A model of solar equilibrium: the hydrodynamic limit

Authors: Lee Gunderson, Amitava Bhattacharjee

Abstract: Helioseismology has revealed the internal density and rotation profiles of the Sun. Yet, knowledge of its magnetic fields and meridional circulation is confined much closer to the surface, and latitudinal entropy gradients are below detectable limits. While numerical simulations can offer insight into the interior dynamics and help identify which ingredients are necessary to reproduce particular o… ▽ More Helioseismology has revealed the internal density and rotation profiles of the Sun. Yet, knowledge of its magnetic fields and meridional circulation is confined much closer to the surface, and latitudinal entropy gradients are below detectable limits. While numerical simulations can offer insight into the interior dynamics and help identify which ingredients are necessary to reproduce particular observations, some features of the Sun can be understood analytically from an equilibrium perspective. Examples of such features include: the 1D density profile arising from steady-state energy transport from the core to the surface, and the tilting of isorotation contours in the convection zone (CZ) due to baroclinic forcing. To help identify which features can be explained by equilibrium, we propose analyzing stationary axisymmetric ideal magnetohydrodynamic flows in the solar regime. By prescribing an appropriate entropy profile at the surface, we recover a rotation profile that reasonably matches observations in the bulk of the CZ. Additionally, by including the effects of poloidal flow, we reproduce a feature that is reminiscent of the near surface shear layer. However, no tachocline-like feature is seen in hydrodynamic equilibrium, suggesting the importance of either dynamics or magnetic fields in its description. △ Less

Submitted 27 January, 2019; v1 submitted 21 December, 2017; originally announced December 2017.

The dehydration of water worlds via atmospheric losses

Authors: Chuanfei Dong, Zhenguang Huang, Manasvi Lingam, Gábor Tóth, Tamas Gombosi, Amitava Bhattacharjee

Abstract: We present a three-species multi-fluid MHD model (H$^+$, H$_2$O$^+$ and e$^-$), endowed with the requisite atmospheric chemistry, that is capable of accurately quantifying the magnitude of water ion losses from exoplanets. We apply this model to a water world with Earth-like parameters orbiting a Sun-like star for three cases: (i) current normal solar wind conditions, (ii) ancient normal solar win… ▽ More We present a three-species multi-fluid MHD model (H$^+$, H$_2$O$^+$ and e$^-$), endowed with the requisite atmospheric chemistry, that is capable of accurately quantifying the magnitude of water ion losses from exoplanets. We apply this model to a water world with Earth-like parameters orbiting a Sun-like star for three cases: (i) current normal solar wind conditions, (ii) ancient normal solar wind conditions, and (iii) one extreme "Carrington-type" space weather event. We demonstrate that the ion escape rate for (ii), with a value of 6.0$\times$10$^{26}$ s$^{-1}$, is about an order of magnitude higher than the corresponding value of 6.7$\times$10$^{25}$ s$^{-1}$ for (i). Studies of ion losses induced by space weather events, where the ion escape rates can reach $\sim$ 10$^{28}$ s$^{-1}$, are crucial for understanding how an active, early solar-type star (e.g., with frequent coronal mass ejections) could have accelerated the depletion of the exoplanet's atmosphere. We briefly explore the ramifications arising from the loss of water ions, especially for planets orbiting M-dwarfs where such effects are likely to be significant. △ Less

Submitted 4 September, 2017; originally announced September 2017.

Comments: ApJ Letters, in press, 8 pages, 2 figures, 3 tables

Journal ref: The Astrophysical Journal Letters, Vol. 847, L4 (2017)

Plasmoid Instability in Evolving Current Sheets and Onset of Fast Reconnection

Authors: Yi-Min Huang, Luca Comisso, A. Bhattacharjee

Abstract: The scaling of plasmoid instability maximum linear growth rate with respect to Lundquist number $S$ in a Sweet-Parker current sheet, $γ_{max}\sim S^{1/4}$, indicates that at high $S$, the current sheet will break apart before it approaches the Sweet-Parker width. Therefore, a proper description for the onset of the plasmoid instability must incorporate the evolving process of the current sheet. We… ▽ More The scaling of plasmoid instability maximum linear growth rate with respect to Lundquist number $S$ in a Sweet-Parker current sheet, $γ_{max}\sim S^{1/4}$, indicates that at high $S$, the current sheet will break apart before it approaches the Sweet-Parker width. Therefore, a proper description for the onset of the plasmoid instability must incorporate the evolving process of the current sheet. We carry out a series of two-dimensional simulations and develop diagnostics to separate fluctuations from an evolving background. It is found that the fluctuation amplitude starts to grow only when the linear growth rate is sufficiently large ($γ_{max}τ_{A}>O(1)$) to overcome convective losses. The linear growth rate continues to rise until the sizes of plasmoids become comparable to the inner layer width of the tearing mode. At this point the current sheet is disrupted and the instability enters the early nonlinear regime. The growth rate suddenly decreases, but the fluctuation amplitude continues to grow until it reaches nonlinear saturation. We identify important time scales of the instability development, as well as scalings for the linear growth rate, current sheet width, and dominant wavenumber at current sheet disruption. These scalings depend on not only the Lundquist number, but also the initial noise amplitude. A phenomenological model that reproduces scalings from simulation results is proposed. The model incorporates the effect of reconnection outflow, which is crucial for yielding a critical Lundquist number $S_{c}$ below which disruption does not occur. The critical Lundquist number $S_{c}$ is not a constant value but has a weak dependence on the noise amplitude. △ Less

Submitted 3 October, 2017; v1 submitted 6 July, 2017; originally announced July 2017.

Journal ref: Astrophysical Journal, 849, 75 (2017)

Plasmoid Instability in Forming Current Sheets

Authors: Luca Comisso, Manasvi Lingam, Yi-Min Huang, Amitava Bhattacharjee

Abstract: The plasmoid instability has revolutionized our understanding of magnetic reconnection in astrophysical environments. By preventing the formation of highly elongated reconnection layers, it is crucial in enabling the rapid energy conversion rates that are characteristic of many astrophysical phenomena. Most of the previous studies have focused on Sweet-Parker current sheets, which, however, are un… ▽ More The plasmoid instability has revolutionized our understanding of magnetic reconnection in astrophysical environments. By preventing the formation of highly elongated reconnection layers, it is crucial in enabling the rapid energy conversion rates that are characteristic of many astrophysical phenomena. Most of the previous studies have focused on Sweet-Parker current sheets, which, however, are unattainable in typical astrophysical systems. Here, we derive a general set of scaling laws for the plasmoid instability in resistive and visco-resistive current sheets that evolve over time. Our method relies on a principle of least time that enables us to determine the properties of the reconnecting current sheet (aspect ratio and elapsed time) and the plasmoid instability (growth rate, wavenumber, inner layer width) at the end of the linear phase. After this phase the reconnecting current sheet is disrupted and fast reconnection can occur. The scaling laws of the plasmoid instability are \emph{not} simple power laws, and depend on the Lundquist number ($S$), the magnetic Prandtl number ($P_m$), the noise of the system ($ψ_0$), the characteristic rate of current sheet evolution ($1/τ$), as well as the thinning process. We also demonstrate that previous scalings are inapplicable to the vast majority of the astrophysical systems. We explore the implications of the new scaling relations in astrophysical systems such as the solar corona and the interstellar medium. In both these systems, we show that our scaling laws yield values for the growth rate, wavenumber, and aspect ratio that are much smaller than the Sweet-Parker based scalings. △ Less

Submitted 6 July, 2017; originally announced July 2017.

Comments: 20 pages; Comments Welcome

Journal ref: Astrophys. J. 850, 142 (2017)

Constructing current singularity in a 3D line-tied plasma

Authors: Yao Zhou, Yi-Min Huang, Hong Qin, A. Bhattacharjee

Abstract: We revisit Parker's conjecture of current singularity formation in 3D line-tied plasmas using a recently developed numerical method, variational integration for ideal magnetohydrodynamics in Lagrangian labeling. With the frozen-in equation built-in, the method is free of artificial reconnection, and hence it is arguably an optimal tool for studying current singularity formation. Using this method,… ▽ More We revisit Parker's conjecture of current singularity formation in 3D line-tied plasmas using a recently developed numerical method, variational integration for ideal magnetohydrodynamics in Lagrangian labeling. With the frozen-in equation built-in, the method is free of artificial reconnection, and hence it is arguably an optimal tool for studying current singularity formation. Using this method, the formation of current singularity has previously been confirmed in the Hahm--Kulsrud--Taylor problem in 2D. In this paper, we extend this problem to 3D line-tied geometry. The linear solution, which is singular in 2D, is found to be smooth for arbitrary system length. However, with finite amplitude, the linear solution can become pathological when the system is sufficiently long. The nonlinear solutions turn out to be smooth for short systems. Nonetheless, the scaling of peak current density versus system length suggests that the nonlinear solution may become singular at finite length. With the results in hand, we can neither confirm nor rule out this possibility conclusively, since we cannot obtain solutions with system length near the extrapolated critical value. △ Less

Submitted 1 January, 2018; v1 submitted 9 May, 2017; originally announced May 2017.

Journal ref: Astrophys. J. 852, 3 (2018)

Nonlinear Resistivity for Magnetohydrodynamical Models

Authors: Manasvi Lingam, Eero Hirvijoki, David Pfefferlé, Luca Comisso, Amitava Bhattacharjee

Abstract: A new formulation of the plasma resistivity that arises from the collisional momentum-transfer rate between electrons and ions is presented. The resistivity computed herein is shown to depend not only on the temperature and density but also on all other polynomial velocity-space moments of the distribution function, such as the pressure tensor and heat flux vector. The exact expression for the col… ▽ More A new formulation of the plasma resistivity that arises from the collisional momentum-transfer rate between electrons and ions is presented. The resistivity computed herein is shown to depend not only on the temperature and density but also on all other polynomial velocity-space moments of the distribution function, such as the pressure tensor and heat flux vector. The exact expression for the collisional momentum-transfer rate is determined, and is used to formulate the nonlinear anisotropic resistivity. The new formalism recovers the Spitzer resistivity, as well as the concept of thermal force if the heat flux is assumed to be proportional to a temperature gradient. Furthermore, if the pressure tensor is related to viscous stress, the latter enters the expression for the resistivity. The relative importance of the nonlinear term(s) with respect to the well-established electron inertia and Hall terms is also examined. The subtle implications of the nonlinear resistivity, and its dependence on the fluid variables, are discussed in the context of magnetized plasma environments and phenomena such as magnetic reconnection. △ Less

Submitted 15 March, 2017; v1 submitted 17 October, 2016; originally announced October 2016.

Comments: 9 pages; 0 figures

Journal ref: Phys. Plasmas 24, 042120 (2017)

Reverberation Mapping of Optical Emission Lines in Five Active Galaxies

Authors: M. M. Fausnaugh, C. J. Grier, M. C. Bentz, K. D. Denney, G. De Rosa, B. M. Peterson, C. S. Kochanek, R. W. Pogge, S. M. Adams, A. J. Barth, Thomas G. Beatty, A. Bhattacharjee, G. A. Borman, T. A. Boroson, M. C. Bottorff, Jacob E. Brown, Jonathan S. Brown, M. S. Brotherton, C. T. Coker, S. M. Crawford, K. V. Croxall, Sarah Eftekharzadeh, Michael Eracleous, M. D. Joner, C. B. Henderson , et al. (46 additional authors not shown)

Abstract: We present the first results from an optical reverberation mapping campaign executed in 2014, targeting the active galactic nuclei (AGN) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a "changing look" AGN and a broad-line radio galaxy. Based on continuum-H$β$ lags, we measure black hole masses for all five targe… ▽ More We present the first results from an optical reverberation mapping campaign executed in 2014, targeting the active galactic nuclei (AGN) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a "changing look" AGN and a broad-line radio galaxy. Based on continuum-H$β$ lags, we measure black hole masses for all five targets. We also obtain H$γ$ and He{\sc ii}\,$λ4686$ lags for all objects except 3C 382. The He{\sc ii}\,$λ4686$ lags indicate radial stratification of the BLR, and the masses derived from different emission lines are in general agreement. The relative responsivities of these lines are also in qualitative agreement with photoionization models. These spectra have extremely high signal-to-noise ratios (100--300 per pixel) and there are excellent prospects for obtaining velocity-resolved reverberation signatures. △ Less

Submitted 18 December, 2017; v1 submitted 30 September, 2016; originally announced October 2016.

Comments: 34 pages, 12 figures, published in ApJ. For a video summarizing the main results, see https://www.youtube.com/watch?v=KaC-jPsIY0Q