Showing 1–37 of 37 results for author: Deka, P

The MeerKAT Absorption Line Survey Data Release 2: Wideband continuum catalogues and a measurement of the cosmic radio dipole

Authors: J. D. Wagenveld, H-R. Klöckner, N. Gupta, S. Sekhar, P. Jagannathan, P. P. Deka, J. Jose, S. A. Balashev, D. Borgaonkar, A. Chatterjee, F. Combes, K. L. Emig, A. N. Gaunekar, M. Hilton, G. I. G. Józsa, D. Y. Klutse, K. Knowles, J. -K. Krogager, E. Momjian, S. Muller, S. P. Sikhosana

Abstract: We present the second data release of the MeerKAT Absorption Line Survey (MALS), consisting of wideband continuum catalogues of 391 pointings observed at L~band. The full wideband catalogue covers 4344 deg$^2$ of sky, reaches a depth of 10 $μ$Jy beam$^{-1}$, and contains 971,980 sources. With its balance between survey depth and sky coverage, MALS DR2 covers five orders of magnitude of flux densit… ▽ More We present the second data release of the MeerKAT Absorption Line Survey (MALS), consisting of wideband continuum catalogues of 391 pointings observed at L~band. The full wideband catalogue covers 4344 deg$^2$ of sky, reaches a depth of 10 $μ$Jy beam$^{-1}$, and contains 971,980 sources. With its balance between survey depth and sky coverage, MALS DR2 covers five orders of magnitude of flux density, presenting a robust view of the extragalactic radio source population down to 200 $μ$Jy. Using this catalogue, we perform a measurement of the cosmic radio dipole, an anisotropy in the number counts of radio sources with respect to the cosmic background, analogous to the dipole found in the cosmic microwave background (CMB). For this measurement, we present the characterisation of completeness and noise properties of the catalogue, and show that a declination-dependent systematic affects the number density of faint sources. In the dipole measurement on the MALS catalogue, we recover reasonable dipole measurements once we model the declination systematic with a linear fit between the size of the major axis of the restoring beam and the amount of sources of each pointing. The final results are consistent with the CMB dipole in terms of direction and amplitude, unlike many recent measurements of the cosmic radio dipole made with other centimetre wavelength catalogues, which generally show a significantly larger amplitude. This result demonstrates the value of dipole measurements with deeper and more sparse radio surveys, as the population of faint sources probed may have had a significant impact on the measured dipole. △ Less

Submitted 29 August, 2024; originally announced August 2024.

Comments: 25 pages, 14 figures. Accepted for publication in Astronomy & Astrophysics. The MALS wideband catalogues and images are publicly available at https://mals.iucaa.in

AttackER: Towards Enhancing Cyber-Attack Attribution with a Named Entity Recognition Dataset

Authors: Pritam Deka, Sampath Rajapaksha, Ruby Rani, Amirah Almutairi, Erisa Karafili

Abstract: Cyber-attack attribution is an important process that allows experts to put in place attacker-oriented countermeasures and legal actions. The analysts mainly perform attribution manually, given the complex nature of this task. AI and, more specifically, Natural Language Processing (NLP) techniques can be leveraged to support cybersecurity analysts during the attribution process. However powerful t… ▽ More Cyber-attack attribution is an important process that allows experts to put in place attacker-oriented countermeasures and legal actions. The analysts mainly perform attribution manually, given the complex nature of this task. AI and, more specifically, Natural Language Processing (NLP) techniques can be leveraged to support cybersecurity analysts during the attribution process. However powerful these techniques are, they need to deal with the lack of datasets in the attack attribution domain. In this work, we will fill this gap and will provide, to the best of our knowledge, the first dataset on cyber-attack attribution. We designed our dataset with the primary goal of extracting attack attribution information from cybersecurity texts, utilizing named entity recognition (NER) methodologies from the field of NLP. Unlike other cybersecurity NER datasets, ours offers a rich set of annotations with contextual details, including some that span phrases and sentences. We conducted extensive experiments and applied NLP techniques to demonstrate the dataset's effectiveness for attack attribution. These experiments highlight the potential of Large Language Models (LLMs) capabilities to improve the NER tasks in cybersecurity datasets for cyber-attack attribution. △ Less

Submitted 9 August, 2024; originally announced August 2024.

Comments: Submitted to WISE 2024

Classification of synchronization in nonlinear systems using ICO learning

Authors: J. P. Deka

Abstract: In this work, we investigate the implications of the differential Hebbian learning rule known as Input-Correlations (ICO) learning in the classification of synchronization in coupled nonlinear oscillator systems. We are investigating the parity-time symmetric coupled Duffing oscillator system with nonlinear dissipation/amplification. In our investigation of the temporal dynamics of this system, it… ▽ More In this work, we investigate the implications of the differential Hebbian learning rule known as Input-Correlations (ICO) learning in the classification of synchronization in coupled nonlinear oscillator systems. We are investigating the parity-time symmetric coupled Duffing oscillator system with nonlinear dissipation/amplification. In our investigation of the temporal dynamics of this system, it is observed that the system exhibits chaotic as well as quasiperiodic dynamics. On further investigation, it is found that the chaotic dynamics is distorted anti-phase synchronized, whereas the quasiperiodic dynamics is desynchronized. So, on the application of the ICO learning in these two parametric regimes, we observe that the weight associated with the stimulus remains constant when the oscillators are anti-phase synchronized, in spite of there being distortion in the synchronization. But when the oscillators exhibit quasiperiodic dynamics, there is erratic evolution of the weight with time. So, from this, it could be ascertained that the ICO learning could be made use of in the classification of synchronization dynamics in nonlinear systems. △ Less

Submitted 7 August, 2024; originally announced August 2024.

Comments: 6 pages, 5 figures

Temporal Dynamics beyond the Exceptional Point in the Ikeda Map with Balanced Gain and Loss

Authors: Jyoti Prasad Deka, Amarendra K. Sarma

Abstract: We investigate the temporal dynamics of the Ikeda Map with Balanced Gain and Loss and in the presence of feedback loops with saturation nonlinearity. From the bifurcation analysis, we find that the temporal evolution of optical power undergoes period quadrupling at the exceptional point (EP) of the system and beyond that, chaotic dynamics emerge in the system and this has been further corroborated… ▽ More We investigate the temporal dynamics of the Ikeda Map with Balanced Gain and Loss and in the presence of feedback loops with saturation nonlinearity. From the bifurcation analysis, we find that the temporal evolution of optical power undergoes period quadrupling at the exceptional point (EP) of the system and beyond that, chaotic dynamics emerge in the system and this has been further corroborated from the Largest Lyapunov Exponent (LLE) of the model. For a closer inspection, we analyzed the parameter basin of the system, which further leads to our inference that the Ikeda Map with Balanced Gain and Loss exhibits the emergence of chaotic dynamics beyond the exceptional point (EP). Furthermore, we find that the temporal dynamics beyond the EP regime leads to the onset of Extreme Events (EE) in this system via attractor merging crisis. △ Less

Submitted 13 May, 2024; originally announced June 2024.

ICO learning as a measure of transient chaos in PT-symmetric Liénard systems

Authors: J. P. Deka, A. Govindarajan, A. K. Sarma

Abstract: In this article, we investigate the implications of the unsupervised learning rule known as Input-Correlations (ICO) learning in the nonlinear dynamics of two linearly coupled PT-symmetric Liénard oscillators. The fixed points of the oscillator have been evaluated analytically and the Jacobian linearization is employed to study their stability. We find that on increasing the amplitude of the exter… ▽ More In this article, we investigate the implications of the unsupervised learning rule known as Input-Correlations (ICO) learning in the nonlinear dynamics of two linearly coupled PT-symmetric Liénard oscillators. The fixed points of the oscillator have been evaluated analytically and the Jacobian linearization is employed to study their stability. We find that on increasing the amplitude of the external periodic drive, the system exhibits period-doubling cascade to chaos within a specific parametric regime wherein we observe emergent chaotic dynamics. We further notice that the system indicates an intermittency route to chaos in the chaotic regime. Finally, in the period-4 regime of our bifurcation analysis, we predict the emergence of transient chaos which eventually settles down to a period-2 oscillator response which has been further validated by both the maximal Finite-Time Lyapunov Exponent (FTLE) using the well-known Gram-Schmidt orthogonalization technique and the Hilbert Transform of the time-series. In the transiently chaotic regime, we deploy the ICO learning to analyze the time-series from which we identify that when the chaotic evolution transforms into periodic dynamics, the synaptic weight associated with the time-series of the loss oscillator exhibits stationary temporal evolution. This signifies that in the periodic regime, there is no overlap between the filtered signals obtained from the time-series of the coupled PT-symmetric oscillators. In addition, the temporal evolution of the weight associated with the stimulus mimics the behaviour of the Hilbert transform of the time-series. △ Less

Submitted 18 June, 2024; v1 submitted 14 May, 2024; originally announced May 2024.

Comments: 9 pages, 12 figures

Onset of Identical Synchronization in the Spatial Evolution of Optical Power in a Waveguide Coupler

Authors: J. P. Deka

Abstract: In this work, we investigated the spatial evolution of optical power in a closed-form optical waveguide configuration consisting of six passive waveguides and each of the waveguides exhibits equal strength of Kerr nonlinearity. We considered only nearest neighbor interaction between the waveguides. We found that in the case of low Kerr nonlinearity, evolution of optical power shows synchronization… ▽ More In this work, we investigated the spatial evolution of optical power in a closed-form optical waveguide configuration consisting of six passive waveguides and each of the waveguides exhibits equal strength of Kerr nonlinearity. We considered only nearest neighbor interaction between the waveguides. We found that in the case of low Kerr nonlinearity, evolution of optical power shows synchronization behavior. But when we increased the strength of Kerr nonlinearity, we discovered that spatial evolution of optical power in all waveguides shows independent characteristics. On the other hand, we have studied the impact of the coupling constant on the synchronization dynamics of our system. Our findings showed us that strong coupling can strengthen the collective dynamics in the presence of strong Kerr nonlinearity. From our results, we can conclude that Kerr nonlinearity in our system plays the role of disorder parameter that destroys as well as alters the synchronization behavior of evolution of optical power in the waveguides and coupling constant plays the role of an antagonist and restores synchronization in the model. △ Less

Submitted 7 July, 2024; v1 submitted 1 May, 2024; originally announced May 2024.

Comments: 7 pages, 6 figures

Journal ref: Indian Journal of Physics (2024)

Scalar Dark Matter and Stability of Higgs Vacuum within a Minimal SO(10) GUT Model

Authors: Nilavjyoti Hazarika, Paramita Deka, Kalpana Bora

Abstract: In this work, we delve to investigate a feasible range of dark matter (DM) masses within a non-supersymmetric $SO(10)$ Grand Unified Theory (GUT) scalar dark matter model, in freeze-out scenario. This model includes a singlet scalar denoted as S and an inert doublet represented by $φ$. Being part of SO(10), the quantum numbers of DM particles are assigned and hence we know their nature. These fiel… ▽ More In this work, we delve to investigate a feasible range of dark matter (DM) masses within a non-supersymmetric $SO(10)$ Grand Unified Theory (GUT) scalar dark matter model, in freeze-out scenario. This model includes a singlet scalar denoted as S and an inert doublet represented by $φ$. Being part of SO(10), the quantum numbers of DM particles are assigned and hence we know their nature. These fields are odd under a discrete $Z_2$ matter parity $(-1)^{3(B-L)}$. The dark matter with mass 300 $\leq M_{DM} \leq$ 1000 GeV emerges as a mixture of the $Z_2$-odd scalar singlet $S$ and the neutral element of the doublet $φ$, both residing within a \textbf{16}-dim scalar representation of $SO(10)$. In this work we consider a real scalar \textbf{S} belonging to \textbf{16} of SO(10) as DM. We also investigate the one-loop vacuum stability through the solution of Renormalization Group Equations (RGEs) for the model's parameters. Subsequently, we scrutinize the model's predictions within the confines of contemporary theoretical and experimental restrictions. In addition to achieving vacuum stability in the SO(10) framework, the DM mass is found to lie in the previously unaddressed (theoretically) intermediate mass region of 300$\leqslant$ MDM $\leqslant$ 1000 GeV which also adheres to many current phenomenological constraints, like recent direct experimental bounds from XENON1T, indirect detection bounds from Fermi-LAT experiment, Higgs invisible decay and Electroweak Precision Test. The stability of the electroweak vacuum is seen to be present up to the Planck scale. These model predictions possess the potential for future validation through dark matter search experiments, as they are testable in future DM search experiments, along with the added feature that the model is a part of an elegant grand unified theory. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 17 pages, 6 figures, 1 table

MALS discovery of a rare HI 21-cm absorber at $z\sim1.35$: origin of the absorbing gas in powerful AGN

Authors: P. P. Deka, N. Gupta, H. W. Chen, S. D. Johnson, P. Noterdaeme, F. Combes, E. Boettcher, S. A. Balashev, K. L. Emig, G. I. G. Józsa, H. -R. Klöckner, J-. K. Krogager, E. Momjian, P. Petitjean, G. C. Rudie, J. Wagenveld, F. S. Zahedy

Abstract: We report a new, rare detection of HI 21-cm absorption associated with a quasar (only six known at $1<z<2$) here towards J2339-5523 at $z_{em}$ = 1.3531, discovered through the MeerKAT Absorption Line Survey (MALS). The absorption profile is broad ($\sim 400$ km/s), and the peak is redshifted by $\sim 200$ km/s, from $z_{em}$. Interestingly, optical/FUV spectra of the quasar from Magellan-MIKE/HST… ▽ More We report a new, rare detection of HI 21-cm absorption associated with a quasar (only six known at $1<z<2$) here towards J2339-5523 at $z_{em}$ = 1.3531, discovered through the MeerKAT Absorption Line Survey (MALS). The absorption profile is broad ($\sim 400$ km/s), and the peak is redshifted by $\sim 200$ km/s, from $z_{em}$. Interestingly, optical/FUV spectra of the quasar from Magellan-MIKE/HST-COS spectrographs do not show any absorption features associated with the 21-cm absorption. This is despite the coincident presence of the optical quasar and the radio `core' inferred from a flat spectrum component of flux density $\sim 65$ mJy at high frequencies ($>5$ GHz). The simplest explanation would be that no large HI column (N(HI)$>10^{17}$ cm$^{-2}$) is present towards the radio `core' and the optical AGN. Based on the joint optical and radio analysis of a heterogeneous sample of 16 quasars ($z_{median}$ = 0.7) and 15 radio galaxies ($z_{median}$ = 0.3) with HI 21-cm absorption detection and matched in 1.4 GHz luminosity (L$_{\rm 1.4\,GHz}$), a consistent picture emerges where quasars are primarily tracing the gas in the inner circumnuclear disk and cocoon created by the jet-ISM interaction. These exhibit L$_{1.4\,\rm GHz}$ - $ΔV_{\rm null}$ correlation, and frequent mismatch between the radio and optical spectral lines. The radio galaxies show no such correlation and likely trace the gas from the cocoon and the galaxy-wide ISM outside the photoionization cone. The analysis presented here demonstrates the potential of radio spectroscopic observations to reveal the origin of the absorbing gas associated with AGN that may be missed in optical observations. △ Less

Submitted 19 February, 2024; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: 10 pages, 8 figures, accepted for publication in A&A

Anti-Phase Synchronization of Chaos in PT-Symmetric Nonlinear Oscillators

Authors: Jyoti Prasad Deka

Abstract: We investigate the temporal dynamics of the PT-Symmetric nonlinear oscillators in the presence of Duffing nonlinearity for two forms of oscillator configuration. In the former, we consider two oscillator coupled to each other. One oscillator is amplified and the other is attenuated. From the bifurcation analysis, we find that the temporal evolution of oscillators exhibit the transition from quasip… ▽ More We investigate the temporal dynamics of the PT-Symmetric nonlinear oscillators in the presence of Duffing nonlinearity for two forms of oscillator configuration. In the former, we consider two oscillator coupled to each other. One oscillator is amplified and the other is attenuated. From the bifurcation analysis, we find that the temporal evolution of oscillators exhibit the transition from quasiperiodic to chaotic dynamics. This has been corroborated by the maximal Lyapunov exponent of the system. Furthermore, on investigating the correlation of the time-series using the Pearson's correlation coefficient, it is found that the chaotic system is anti-phase synchronized, whereas the quasiperiodic is not synchronized in any form. The parametric regime where this transition has been observed is from the Unbroken PT regime to the Broken PT regime. Similarly, in the latter configuration with two amplified oscillators coupled to two attenuated oscillators, a similar transition has been observed. But in the neighbourhood of the Exceptional Point (EP) of the system, the system is shown to exhibit in-phase synchronized dynamics as is evident from the correlation analysis. △ Less

Submitted 31 August, 2023; originally announced October 2023.

Comments: 14 pages, 11 figures

LeXInt: GPU-accelerated Exponential Integrators package

Authors: Pranab J. Deka, Alexander Moriggl, Lukas Einkemmer

Abstract: We present an open-source CUDA-based package that consists of a compilation of exponential integrators where the action of the matrix exponential or the $\varphi_l$ functions on a vector is approximated using the method of polynomial interpolation at Leja points. Using a couple of test examples on an NVIDIA A100 GPU, we show that one can achieve significant speedups using CUDA over the correspondi… ▽ More We present an open-source CUDA-based package that consists of a compilation of exponential integrators where the action of the matrix exponential or the $\varphi_l$ functions on a vector is approximated using the method of polynomial interpolation at Leja points. Using a couple of test examples on an NVIDIA A100 GPU, we show that one can achieve significant speedups using CUDA over the corresponding CPU code. LeXInt, written in a modular format, facilitates easy integration into any existing software package, and can be used for temporal integration of any differential equation. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: Code available at https://github.com/Pranab-JD/LeXInt; comments welcome

The MeerKAT Absorption Line Survey (MALS) data release I: Stokes I image catalogs at 1-1.4 GHz

Authors: P. P. Deka, N. Gupta, P. Jagannathan, S. Sekhar, E. Momjian, S. Bhatnagar, J. Wagenveld, H. -R. Klöckner, J. Jose, S. A. Balashev, F. Combes, M. Hilton, D. Borgaonkar, A. Chatterjee, K. L. Emig, A. N. Gaunekar, G. I. G. Józsa, D. Y. Klutse, K. Knowles, J-. K. Krogager, A. Mohapatra, K. Moodley, Sébastien Muller, P. Noterdaeme, P. Petitjean , et al. (2 additional authors not shown)

Abstract: The MeerKAT Absorption Line Survey (MALS) has observed 391 telescope pointings at L-band (900 - 1670 MHz) at $δ\lesssim$ $+20°$. We present radio continuum images and a catalog of 495,325 (240,321) radio sources detected at a signal-to-noise ratio (SNR) $>$5 over an area of 2289 deg$^2$ (1132 deg$^2$) at 1006 MHz (1381 MHz). Every MALS pointing contains a central bright radio source (… ▽ More The MeerKAT Absorption Line Survey (MALS) has observed 391 telescope pointings at L-band (900 - 1670 MHz) at $δ\lesssim$ $+20°$. We present radio continuum images and a catalog of 495,325 (240,321) radio sources detected at a signal-to-noise ratio (SNR) $>$5 over an area of 2289 deg$^2$ (1132 deg$^2$) at 1006 MHz (1381 MHz). Every MALS pointing contains a central bright radio source ($S_{1\,\mathrm{GHz}} \gtrsim 0.2$ Jy). The median spatial resolution is $12^{\prime\prime}$ ($8^{\prime\prime}$). The median rms noise away from the pointing center is 25 $μ$Jy beam$^{-1}$ (22 $μ$Jy beam$^{-1}$) and is within $\sim$ 15% of the achievable theoretical sensitivity. The flux density scale ratio and astrometric accuracy deduced from multiply observed sources in MALS are less than 1% (8% scatter) and $1^{\prime\prime}$, respectively. Through comparisons with NVSS and FIRST at 1.4 GHz, we establish the catalog's accuracy in the flux density scale and astrometry to be better than 6% (15% scatter) and $0.8^{\prime\prime}$, respectively. The median flux density offset is higher (9%) for an alternate beam model based on holographic measurements. The MALS radio source counts at 1.4 GHz are in agreement with literature. We estimate spectral indices ($α$) of a subset of 125,621 sources (SNR$>$8), confirm the flattening of spectral indices with decreasing flux density and identify 140 ultra steep-spectrum ($α<-1.3$) sources as prospective high-$z$ radio galaxies ($z>2$). We have identified 1308 variable and 122 transient radio sources comprising primarily of AGN that demonstrate long-term (26 years) variability in their observed flux densities. The MALS catalogs and images are publicly available at https://mals.iucaa.in. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Comments: 64 pages, 25 figures, accepted for publication in the ApJS (full version of the paper with complete tables is available at DR1 release notes)

Sensitivity of octant of $θ_{23}$, CP violation and mass hierarchy in NO$ν$A with multinucleon and detector effects

Authors: Paramita Deka, Kalpana Bora

Abstract: In this work, we investigate how multinucleon enhancement and RPA (Random Phase Approximation) suppression can affect the measurement of three unknown neutrino oscillation parameters - the CP-violating phase $δ_{CP}$, the octant of the atmospheric mixing angle $θ_{23}$, and the determination of the mass hierarchy, in the appearance channel of the NO$ν$A experiment. We include the presence of the d… ▽ More In this work, we investigate how multinucleon enhancement and RPA (Random Phase Approximation) suppression can affect the measurement of three unknown neutrino oscillation parameters - the CP-violating phase $δ_{CP}$, the octant of the atmospheric mixing angle $θ_{23}$, and the determination of the mass hierarchy, in the appearance channel of the NO$ν$A experiment. We include the presence of the detector effect as well in the analysis, which is crucial for capturing realistic experimental scenarios. We also conducted a comparison between the nuclear model Effective Spectral Function (calculated within the RFG model) with and without Transverse Enhancement in terms of sensitivity analysis. It is found that the analysis using our comprehensive model QE(+RPA)+2p2h along with Effective Spectral Function+Transverse Enhancement exhibits significantly enhanced sensitivity compared to the pure QE interaction process, in all the cases.Also, the higher octant of $θ_{23}$, the lower half plane of $δ_{CP}$, and the normal mass hierarchy (HO-LHP-NH) exhibit improved sensitivity, enabling a more precise determination of the corresponding parameters. Furthermore, it is also noted that improving the performance of the detector also improves the results. Thus, including multinucleon effects and improving detector efficiency have the potential to enhance the capabilities of the NO$ν$A (and other long baseline) experiment in conducting precise parameter studies. △ Less

Submitted 21 April, 2024; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 25 pages, 8 figures, 3 tables

Journal ref: JHEP 04 (2024) 095

Efficient numerical methods for Anisotropic Diffusion of Galactic Cosmic Rays

Authors: Pranab J. Deka, Ralf Kissmann, Lukas Einkemmer

Abstract: Anisotropic diffusion is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and the interplay of cosmic rays with the Galactic magnetic field. This diffusion term contributes to the highly stiff nature of the cosmic ray transport equation. To conduct numerical simulations of time-dependent cosmic ray transport, implicit integrators (namely, Crank-Nicolson (CN)) ha… ▽ More Anisotropic diffusion is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and the interplay of cosmic rays with the Galactic magnetic field. This diffusion term contributes to the highly stiff nature of the cosmic ray transport equation. To conduct numerical simulations of time-dependent cosmic ray transport, implicit integrators (namely, Crank-Nicolson (CN)) have been traditionally favoured over the CFL-bound explicit integrators in order to be able to take large step sizes. We propose exponential methods to treat the linear anisotropc diffusion equation in the presence of advection and time-independent and time-dependent sources. These methods allow us to take even larger step sizes that can substantially speed-up the simulations whilst generating highly accurate solutions. In or subsequent work, we will use these exponential solvers in the Picard code to study anisotropic cosmic ray diffusion and we will consider additional physical processes such as continuous momentum losses and reacceleration. △ Less

Submitted 8 August, 2023; v1 submitted 23 July, 2023; originally announced July 2023.

Comments: The 38th International Cosmic Ray Conference (ICRC2023)

Journal ref: PoS (ICRC2023), 444, 040

PageRank Nibble on the sparse directed stochastic block model

Authors: Sayan Banerjee, Prabhanka Deka, Mariana Olvera-Cravioto

Abstract: We present new results on community recovery based on the PageRank Nibble algorithm on a sparse directed stochastic block model (dSBM). Our results are based on a characterization of the local weak limit of the dSBM and the limiting PageRank distribution. This characterization allows us to estimate the probability of misclassification for any given connection kernel and any given number of seeds (… ▽ More We present new results on community recovery based on the PageRank Nibble algorithm on a sparse directed stochastic block model (dSBM). Our results are based on a characterization of the local weak limit of the dSBM and the limiting PageRank distribution. This characterization allows us to estimate the probability of misclassification for any given connection kernel and any given number of seeds (vertices whose community label is known). The fact that PageRank is a local algorithm that can be efficiently computed in both a distributed and asynchronous fashion, makes it an appealing method for identifying members of a given community in very large networks where the identity of some vertices is known. △ Less

Submitted 12 March, 2023; originally announced March 2023.

Comments: PageRank Nibble, directed stochastic block model, local weak convergence, community detection

MSC Class: 0580; 05C81; 05C85; 60J80

The MeerKAT Absorption Line Survey: Homogeneous continuum catalogues towards a measurement of the cosmic radio dipole

Authors: J. D. Wagenveld, H. -R. Klöckner, N. Gupta, P. P. Deka, P. Jagannathan, S. Sekhar, S. A. Balashev, E. Boettcher, F. Combes, K. L. Emig, M. Hilton, G. I. G. Józsa, P. Kamphuis, D. Y. Klutse, K. Knowles, J. -K. Krogager, A. Mohapatra, E. Momjian, K. Moodley, S. Muller, P. Petitjean, P. Salas, S. Sikhosana, R. Srianand

Abstract: The number counts of homogeneous samples of radio sources are a tried and true method of probing the large scale structure of the Universe, as most radio sources outside the galactic plane are at cosmological distances. As such they are expected to trace the cosmic radio dipole, an anisotropy analogous to the dipole seen in the cosmic microwave background (CMB). Results have shown that although th… ▽ More The number counts of homogeneous samples of radio sources are a tried and true method of probing the large scale structure of the Universe, as most radio sources outside the galactic plane are at cosmological distances. As such they are expected to trace the cosmic radio dipole, an anisotropy analogous to the dipole seen in the cosmic microwave background (CMB). Results have shown that although the cosmic radio dipole matches the direction of the CMB dipole, it has a significantly larger amplitude. This result challenges our assumption of the Universe being isotropic, which can have large repercussions for the current cosmological paradigm. Though significant measurements have been made, sensitivity to the radio dipole is generally hampered by systematic effects that can cause large biases in the measurement. Here we assess these systematics with data from the MeerKAT Absorption Line Survey (MALS). We present the analysis of ten MALS pointings, focusing on systematic effects that could lead to an inhomogeneous catalogue. We describe the calibration and creation of full band continuum images and catalogues, producing a combined catalogue containing 16,313 sources and covering 37.5 square degrees of sky down to a sensitivity of 10 $μ$Jy/beam. We measure the completeness, purity, and flux recovery statistics for these catalogues using simulated data. We investigate different source populations in the catalogues by looking at flux densities and spectral indices, and how they might influence source counts. Using the noise characteristics of the pointings, we find global measures that can be used to correct for the incompleteness of the catalogue, producing corrected number counts down to 100 - 200 $μ$Jy. We show that we can homogenise the catalogues and properly account for systematic effects. We determine that we can measure the dipole to $3σ$ significance with 100 MALS pointings. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 30 pages, 25 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 673, A113 (2023)

Local weak limits for collapsed branching processes with random out-degrees

Authors: Sayan Banerjee, Prabhanka Deka, Mariana Olvera-Cravioto

Abstract: We obtain local weak limits in probability for Collapsed Branching Processes (CBP), which are directed random networks obtained by collapsing random-sized families of individuals in a general continuous-time branching process. The local weak limit of a given CBP, as the network grows, is shown to be a related continuous-time branching process stopped at an independent exponential time. This is don… ▽ More We obtain local weak limits in probability for Collapsed Branching Processes (CBP), which are directed random networks obtained by collapsing random-sized families of individuals in a general continuous-time branching process. The local weak limit of a given CBP, as the network grows, is shown to be a related continuous-time branching process stopped at an independent exponential time. This is done through an explicit coupling of the in-components of vertices with the limiting object. We also show that the in-components of a finite collection of uniformly chosen vertices locally weakly converge (in probability) to i.i.d. copies of the above limit, reminiscent of propagation of chaos in interacting particle systems. We obtain as special cases novel descriptions of the local weak limits of directed preferential and uniform attachment models. We also outline some applications of our results for analyzing the limiting in-degree and PageRank distributions. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 23 pages, 3 figures

Discovery of Hydrogen Radio Recombination Lines at z=0.89 towards PKS 1830-211

Authors: Kimberly L. Emig, Neeraj Gupta, Pedro Salas, Sebastien Muller, Sergei A. Balashev, Francoise Combes, Emmanuel Momjian, Yiqing Song, Preshanth Jagannathan, Partha P. Deka, Gyula I. G. Jozsa, Hans-Rainer Klockner, Abhisek Mohapatra, Pasquier Noterdaeme, Patrick Petitjean, Raghunathan Srianand, Jonah D. Wagenveld

Abstract: We report the detection of stimulated hydrogen radio recombination line (RRL) emission from ionized gas in a $z=0.89$ galaxy using 580--1670 MHz observations from the MeerKAT Absorption Line Survey (MALS). The RRL emission originates in a galaxy that intercepts and strongly lenses the radio blazar PKS 1830-211 ($z=2.5$). This is the second detection of RRLs outside of the local universe and the fi… ▽ More We report the detection of stimulated hydrogen radio recombination line (RRL) emission from ionized gas in a $z=0.89$ galaxy using 580--1670 MHz observations from the MeerKAT Absorption Line Survey (MALS). The RRL emission originates in a galaxy that intercepts and strongly lenses the radio blazar PKS 1830-211 ($z=2.5$). This is the second detection of RRLs outside of the local universe and the first clearly associated with hydrogen. We detect effective H144$α$ (and H163$α$) transitions at observed frequencies of 1156 (798) MHz by stacking 17 (27) RRLs with 21$σ$ (14$σ$) significance. The RRL emission contains two main velocity components and is coincident in velocity with HI 21 cm and OH 18 cm absorption. We use the RRL spectral line energy distribution and a Bayesian analysis to constrain the density ($n_e$) and the volume-averaged pathlength ($\ell$) of the ionized gas. We determine $\log( n_e ) = 2.0_{-0.7}^{+1.0}$ cm$^{-3}$ and $\log( \ell ) = -0.7_{-1.1}^{+1.1}$ pc towards the north east (NE) lensed image, likely tracing the diffuse thermal phase of the ionized ISM in a thin disk. Towards the south west (SW) lensed image, we determine $\log( n_e ) = 3.2_{-1.0}^{+0.4}$ cm$^{-3}$ and $\log( \ell ) = -2.7_{-0.2}^{+1.8}$ pc, tracing gas that is more reminiscent of H II regions. We estimate a star formation (surface density) rate of $Σ_{\mathrm{SFR}} \sim 0.6$ M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$ or SFR $\sim 50$ M$_{\odot}$ yr$^{-1}$, consistent with a star-forming main sequence galaxy of $M_{\star} \sim 10^{11}$ M$_{\odot}$. The discovery presented here opens up the possibility of studying ionized gas at high redshifts using RRL observations from current and future (e.g., SKA and ngVLA) radio facilities. △ Less

Submitted 8 January, 2023; originally announced January 2023.

Comments: 17 pages, 10 figures, 2 tables, accepted in ApJ

PKS1413+135: OH and HI at z = 0.247 with MeerKAT

Authors: F. Combes, N. Gupta, S. Muller, S. Balashev, P. Deka, K. Emig, H. -R. Kloeckner, D. Klutse, K. Knowles, A. Mohapatra, E. Momjian, P. Noterdaeme, P. Petitjean, P. Salas, R. Srianand, J. Wagenveld

Abstract: The BL Lac PKS 1413+135 was observed by the Large Survey Project "MeerKAT Absorption Line Survey" (MALS) in the L-band, at 1139 MHz and 1293-1379 MHz, targeting the HI and OH lines in absorption at z = 0.24671. The radio continuum is thought to come from a background object at redshift lower than 0.5, as suggested by the absence of gravitational images. The HI absorption line is detected at high s… ▽ More The BL Lac PKS 1413+135 was observed by the Large Survey Project "MeerKAT Absorption Line Survey" (MALS) in the L-band, at 1139 MHz and 1293-1379 MHz, targeting the HI and OH lines in absorption at z = 0.24671. The radio continuum is thought to come from a background object at redshift lower than 0.5, as suggested by the absence of gravitational images. The HI absorption line is detected at high signal-to-noise, with a narrow central component, and a red wing, confirming previous results. The OH 1720 MHz line is clearly detected in (maser) emission, peaking at a velocity shifted by -10 to -15 km/s with respect to the HI peak. The 1612 MHz line is lost due to radio interferences. The OH 1667 MHz main line is tentatively detected in absorption, but not the 1665 MHz one. Over 30 years, a high variability is observed in optical depths, due to the rapid changes of the line of sight, caused by the superluminal motions of the radio knots. The HI line has varied by 20 per cent in depth, while the OH-1720 MHz depth has varied by a factor 4. The position of the central velocity and the widths also varied. The absorbing galaxy is an early-type spiral (maybe S0) seen edge-on, with a prominent dust lane, covering the whole disk. Given the measured mass concentration, and the radio continuum size at centimeter wavelengths (100 mas corresponding to 400 pc at z = 0.25), the width of absorption lines from the nuclear regions are expected up to 250 km/S. The narrowness of the observed lines (< 15 km/s) suggest that the absorption comes from an outer gas ring, as frequently observed in S0 galaxies. The millimetric lines are even narrower (< 1 km/s), which corresponds to the continuum size restricted to the core. The core source is covered by individual 1 pc molecular clouds, of column density a few 10^22 cm-2, which is compatible with the gas screen detected in X-rays. △ Less

Submitted 2 January, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: 7 pages, 6 figures, accepted in A&A

Journal ref: A&A 671, A43 (2023)

Exponential methods for anisotropic diffusion

Authors: Pranab J. Deka, Lukas Einkemmer, Ralf Kissmann

Abstract: The anisotropic diffusion equation is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and its interplay with the ambient magnetic field. This diffusion term contributes to the highly stiff nature of the CR transport equation. In order to conduct numerical simulations of time-dependent cosmic ray transport, implicit integrators have been traditionally favoured o… ▽ More The anisotropic diffusion equation is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and its interplay with the ambient magnetic field. This diffusion term contributes to the highly stiff nature of the CR transport equation. In order to conduct numerical simulations of time-dependent cosmic ray transport, implicit integrators have been traditionally favoured over the CFL-bound explicit integrators in order to be able to take large step sizes. We propose exponential methods that directly compute the exponential of the matrix to solve the linear anisotropic diffusion equation. These methods allow us to take even larger step sizes; in certain cases, we are able to choose a step size as large as the simulation time, i.e., only one time step. This can substantially speed-up the simulations whilst generating highly accurate solutions (l2 error $\leq 10^{-10}$). Additionally, we test an approach based on extracting a constant diffusion coefficient from the anisotropic diffusion equation, where the constant coefficient term is solved implicitly or exponentially and the remainder is treated using some explicit method. We find that this approach, for homogeneous linear problems, is unable to improve on the exponential-based methods that directly evaluate the matrix exponential. △ Less

Submitted 12 December, 2022; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: in submission

A comparison of Leja- and Krylov-based iterative schemes for Exponential Integrators

Authors: Pranab J. Deka, Mayya Tokman, Lukas Einkemmer

Abstract: Krylov-based algorithms have long been preferred to compute the matrix exponential and exponential-like functions appearing in exponential integrators. Of late, direct polynomial interpolation of the action of these exponential-like functions have been shown to be competitive with the Krylov methods. We analyse the performance of the state-of-the-art Krylov algorithm, KIOPS, and the method of poly… ▽ More Krylov-based algorithms have long been preferred to compute the matrix exponential and exponential-like functions appearing in exponential integrators. Of late, direct polynomial interpolation of the action of these exponential-like functions have been shown to be competitive with the Krylov methods. We analyse the performance of the state-of-the-art Krylov algorithm, KIOPS, and the method of polynomial interpolation at Leja points for a number of exponential integrators for various test problems and with varying amounts of stiffness. Additionally, we investigate the performance of an iterative scheme that combines both the KIOPS and Leja approach, named LeKry, that shows substantial improvements over both the Leja- and Krylov-based methods for certain exponential integrators. Whilst we do manage to single out a favoured iterative scheme for each of the exponential integrators that we consider in this study, we do not find any conclusive evidence for preferring either KIOPS or Leja for different classes of exponential integrators. We are unable to identify a superior exponential integrator, one that performs better than all others, for most, if not all of the problems under consideration. We, however, do find that the performance significantly depends on the interplay between the iterative scheme and the specific exponential integrator under consideration. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: in submission; attached supplemental material

LeXInt: Package for Exponential Integrators employing Leja interpolation

Authors: Pranab J. Deka, Lukas Einkemmer, Mayya Tokman

Abstract: We present a publicly available software for exponential integrators that computes the $\varphi_l(z)$ functions using polynomial interpolation. The interpolation method at Leja points have recently been shown to be competitive with the traditionally-used Krylov subspace method. The developed framework facilitates easy adaptation into any Python software package for time integration. We present a publicly available software for exponential integrators that computes the $\varphi_l(z)$ functions using polynomial interpolation. The interpolation method at Leja points have recently been shown to be competitive with the traditionally-used Krylov subspace method. The developed framework facilitates easy adaptation into any Python software package for time integration. △ Less

Submitted 20 January, 2023; v1 submitted 17 August, 2022; originally announced August 2022.

Comments: Publicly available software available at https://github.com/Pranab-JD/LeXInt, in submission

Journal ref: SoftwareX, 21, 101302 (2023)

Uncertainties in neutrino oscillation parameter sensitivity due to resonance processes at NO$ν$A

Authors: Paramita Deka, Nilavjyoti Hazarika, Kalpana Bora

Abstract: The long baseline (LBL) neutrino experiments use heavy nuclear targets for neutrino scattering in which nuclear effects give rise to complications in measuring the neutrino oscillation parameters up to high precision. These nuclear effects are not yet fully understood and therefore need to be quantified as they contribute to the systematic uncertainties. Precision reconstruction of neutrino energy… ▽ More The long baseline (LBL) neutrino experiments use heavy nuclear targets for neutrino scattering in which nuclear effects give rise to complications in measuring the neutrino oscillation parameters up to high precision. These nuclear effects are not yet fully understood and therefore need to be quantified as they contribute to the systematic uncertainties. Precision reconstruction of neutrino energy is one of the main components in measuring the oscillation parameters, and it is required that the neutrino energy is reconstructed with very high precision. In this work, we investigate the effects of the resonance (RES) interaction process using two models for carbon target for $ν_μ\rightarrowν_μ$ disappearance channel of the NO$ν$A experiment, on neutrino-nucleus scattering cross section, events, and neutrino oscillation parameter sensitivity. We also incorporate the realistic detector specifications of NO$ν$A. To quantify the systematic uncertainties due to RES interactions, we compare the cross-sections, events, and sensitivity analysis for two models of RES processes - Rein-Sehgal and Berger-Sehgal, and comment on which model produces more precision. We observe that RES processes contribute significantly, and should be included carefully in models while extracting neutrino oscillation parameters. △ Less

Submitted 7 January, 2023; v1 submitted 15 August, 2022; originally announced August 2022.

Comments: 11 pages, 8 figures, 3 tables

Journal ref: IJMPA 2022

Impact of Nuclear effects in Energy Reconstruction Methods on Sensitivity of Neutrino Oscillation Parameters at NO$ν$A experiment

Authors: Paramita Deka, Jaydip Singh, Kalpana Bora

Abstract: Long baseline (LBL) neutrino experiments aim to measure the neutrino oscillation parameters to high precision. These experiments use nuclear targets for neutrino scattering and hence are inflicted with complexities of nuclear effects. Nuclear effects and their percolation into sensitivity measurement of neutrino oscillations parameters are not yet fully understood and therefore need to be dealt wi… ▽ More Long baseline (LBL) neutrino experiments aim to measure the neutrino oscillation parameters to high precision. These experiments use nuclear targets for neutrino scattering and hence are inflicted with complexities of nuclear effects. Nuclear effects and their percolation into sensitivity measurement of neutrino oscillations parameters are not yet fully understood and therefore need to be dealt with carefully. In a recent work [1], we reported some results on this for NO$ν$A experiment using the kinematic method of neutrino energy reconstruction, where it was observed that the nuclear effects are important in sensitivity analysis, and inclusion of realistic detector setup specifications increases uncertainty in this analysis as compared to ideal detector case. With this motivation, in this work, we use two methods of neutrino energy reconstruction - kinematic and calorimetric, including the nuclear effects, and study their impact on sensitivity analysis. We consider nuclear interactions such as RPA and 2p2h and compare two energy reconstruction methods with reference to events generation, measurement of neutrino oscillation parameters $Δm_{32}^2$ and $θ_{23}$ for disappearance channel, mass hierarchy sensitivity, and CP-violation sensitivity for appearance channel of the NO$ν$A experiment. It is observed that with an ideal detector setup, the kinematic method shows significant dependence on nuclear effects compared to the calorimetric method. We also investigate the impact of realistic detector setup for NO$ν$A in these two methods (with nuclear effects) and find that the calorimetric method shows more bias (uncertainty increases) in sensitivity contours, as compared to the kinematic method. This is found to be true for both the mass hierarchies and for both neutrino and antineutrino incoming beams. △ Less

Submitted 24 January, 2023; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: 20 pages, 16 figures

Journal ref: Nuclear Physics B, Volume 987, February 2023, 116086

Mapping HI 21-cm in the Klemola 31 group at z = 0.029: emission and absorption towards PKS2020-370

Authors: E. K. Maina, Abhisek Mohapatra, G. I. G. Jozsa, N. Gupta, F. Combes, P. Deka, J. D. Wagenveld, R. Srianand, S. A. Balashev, Hsiao-Wen Chen, J. -K. Krogager, E. Momjian, P. Noterdaeme, P. Petitjean

Abstract: We present MeerKAT Absorption Line Survey (MALS) observations of the HI gas in the Klemola31 galaxy group ($z=0.029$), located along the line of sight to the radio-loud quasar PKS2020-370 ($z=1.048$). Four galaxies of the group are detected in HI emission, and HI absorption is also detected in front of PKS2020-370 in Klemola31A. The emission and absorption are somewhat compensating on the line of… ▽ More We present MeerKAT Absorption Line Survey (MALS) observations of the HI gas in the Klemola31 galaxy group ($z=0.029$), located along the line of sight to the radio-loud quasar PKS2020-370 ($z=1.048$). Four galaxies of the group are detected in HI emission, and HI absorption is also detected in front of PKS2020-370 in Klemola31A. The emission and absorption are somewhat compensating on the line of sight of the quasar, and the derived column density of the absorption appears under-estimated, with respect to the neighbouring emission. A symmetric tilted-ring model of Klemola31A, assuming the absorbing gas in regular rotation in the plane, yields a rather high spin temperature of 530K. An alternative interpretation is that the absorbing gas is extra-planar, which will also account for its non-circular motion. The NaI/CaII ratio also suggests that the absorbing gas is unrelated to cold HI disk. Two of the galaxies in the Klemola group are interacting with a small companion, and reveal typical tidal tails, and velocity perturbations. Only one of the galaxies, ESO400-13, reveals a strong HI deficiency, and a characteristic ram-pressure stripping, with a total asymmetry in the distribution of its gas. Since a small galaxy group as Klemola31 is not expected to host a dense intra-group gas, this galaxy must be crossing the group at a very high velocity, mostly in the sky plane. △ Less

Submitted 26 July, 2022; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: 13 pages, 8 figures, 3 tables, accepted in MNRAS (minor corrections)

Exponential Integrators for Resistive Magnetohydrodynamics: Matrix-free Leja Interpolation and Efficient Adaptive Time Stepping

Authors: Pranab Deka, Lukas Einkemmer

Abstract: We propose a novel algorithm for the temporal integration of the resistive magnetohydrodynamics (MHD) equations. The approach is based on exponential Rosenbrock schemes in combination with Leja interpolation. It naturally preserves Gauss's law for magnetism and is unencumbered by the stability constraints observed for explicit methods. Remarkable progress has been achieved in designing exponential… ▽ More We propose a novel algorithm for the temporal integration of the resistive magnetohydrodynamics (MHD) equations. The approach is based on exponential Rosenbrock schemes in combination with Leja interpolation. It naturally preserves Gauss's law for magnetism and is unencumbered by the stability constraints observed for explicit methods. Remarkable progress has been achieved in designing exponential integrators and computing the required matrix functions efficiently. However, employing them in MHD simulations of realistic physical scenarios requires a matrix-free implementation. We show how an efficient algorithm based on Leja interpolation that only uses the right-hand side of the differential equation (i.e. matrix free), can be constructed. We further demonstrate that it outperforms Krylov-based exponential integrators as well as explicit and implicit methods using test models of magnetic reconnection and the Kelvin--Helmholtz instability. Furthermore, an adaptive step-size strategy that gives excellent and predictable performance, particularly in the lenient- to intermediate-tolerance regime that is often of importance in practical applications, is employed. △ Less

Submitted 2 March, 2022; v1 submitted 31 August, 2021; originally announced August 2021.

Comments: Accepted for publication in ApJS; 20 pages, 13 figures

Journal ref: 2022, ApJS, 259, 57

MALS SALT-NOT survey of MIR-selected powerful radio-bright AGN at 0<z<3.5

Authors: N. Gupta, G. Shukla, R. Srianand, J-. K. Krogager, P. Noterdaeme, A. J. Baker, F. Combes, J. P. U. Fynbo, E. Momjian, M. Hilton, T. Hussain, K. Moodley, P. Petitjean, H. -W. Chen, P. Deka, R. Dutta, J. Jose, G. I. G. Jozsa, C. Kaski, H. -R. Klockner, K. Knowles, S. Sikhosana, J. Wagenveld

Abstract: We present results of an optical spectroscopic survey using SALT and NOT to build a WISE mid-infrared color-based, dust-unbiased sample of powerful radio-bright ($>$200 mJy at 1.4 GHz) AGN for the MeerKAT Absorption Line Survey (MALS). Our sample has 250 AGN (median $z=1.8$) showing emission lines, 26 with no emission lines, and 27 without optical counterparts. Overall, our sample is fainter (… ▽ More We present results of an optical spectroscopic survey using SALT and NOT to build a WISE mid-infrared color-based, dust-unbiased sample of powerful radio-bright ($>$200 mJy at 1.4 GHz) AGN for the MeerKAT Absorption Line Survey (MALS). Our sample has 250 AGN (median $z=1.8$) showing emission lines, 26 with no emission lines, and 27 without optical counterparts. Overall, our sample is fainter ($Δi$=0.6 mag) and redder ($Δ(g-i)$=0.2 mag) than radio-selected quasars, and representative of fainter quasar population detected in optical surveys. About 20% of the sources are narrow line AGN (NLAGN) $-$ 65% of these, at $z < 0.5$ are galaxies without strong nuclear emission, and 10% at $z>1.9$, have emission line ratios similar to radio galaxies. The farthest NLAGN in our sample is M1513$-$2524 ($z_{em}=3.132$), and the largest (size$\sim$330 kpc) is M0909$-$3133 ($z_{em}=0.884$). We discuss in detail 110 AGN at $1.9 < z < 3.5$. Despite representing the radio loudest quasars (median $R$=3685), their Eddington ratios are similar to the SDSS quasars having lower $R$. We detect 4 CIV BALQSOs, all among AGN with least $R$, and highest black hole masses and Eddington ratios. The BAL detection rate ($4^{+3}_{-2}$%) is consistent with that seen in extremely powerful ($L_{1.4GHz}>10^{25}$ WHz$^{-1}$) quasars. Using optical light-curves, radio polarization and $γ$-ray detections, we identify 7 high-probability BL Lacs. We also summarize the full MALS footprint to search for HI 21-cm and OH 18-cm lines at $z<2$. △ Less

Submitted 21 February, 2022; v1 submitted 20 July, 2021; originally announced July 2021.

Comments: 62 pages, 15 figures and 3 tables; accepted in ApJ (updated the redshift of M1312-2026 to z=0.977)

Uncertainties in the oscillation parameters measurement due to multi-nucleon effects at NO$ν$A Experiment

Authors: Paramita Deka, Jaydip Singh, Neelakshi Sarma, Kalpana Bora

Abstract: In this work, we investigate the role of multi-nucleon (MN) effects (mainly 2p-2h and RPA) on the sensitivity measurement of various neutrino oscillation parameters, in the disappearance channel of NO$ν$A (USA) experiment. Short-range correlations and detector effects have also been included in the analysis. We use the kinematical method of reconstruction of the incoming neutrino energy, both at t… ▽ More In this work, we investigate the role of multi-nucleon (MN) effects (mainly 2p-2h and RPA) on the sensitivity measurement of various neutrino oscillation parameters, in the disappearance channel of NO$ν$A (USA) experiment. Short-range correlations and detector effects have also been included in the analysis. We use the kinematical method of reconstruction of the incoming neutrino energy, both at the near and far detectors. The extrapolation technique has been used to estimate oscillated events at the far detector. The latest global best fit values of various light neutrino oscillation parameters have been used in the analysis. We find that MN effects increase uncertainty in the measurement of neutrino oscillation parameters, while lower detector efficiency is reflected in more uncertainty. This study can give useful insight into precision studies at long-baseline neutrino experiments in future measurements. △ Less

Submitted 16 July, 2022; v1 submitted 30 June, 2021; originally announced June 2021.

Comments: 17 pages, 14 figures, 2 tables

Journal ref: Nuclear Physics B, 115903, 2022

Efficient adaptive step size control for exponential integrators

Authors: Pranab Jyoti Deka, Lukas Einkemmer

Abstract: Traditional step size controllers make the tacit assumption that the cost of a time step is independent of the step size. This is reasonable for explicit and implicit integrators that use direct solvers. In the context of exponential integrators, however, an iterative approach, such as Krylov methods or polynomial interpolation, to compute the action of the required matrix functions is usually emp… ▽ More Traditional step size controllers make the tacit assumption that the cost of a time step is independent of the step size. This is reasonable for explicit and implicit integrators that use direct solvers. In the context of exponential integrators, however, an iterative approach, such as Krylov methods or polynomial interpolation, to compute the action of the required matrix functions is usually employed. In this case, the assumption of constant cost is not valid. This is, in particular, a problem for higher-order exponential integrators, which are able to take relatively large time steps based on accuracy considerations. In this paper, we consider an adaptive step size controller for exponential Rosenbrock methods that determines the step size based on the premise of minimizing computational cost. The largest allowed step size, given by accuracy considerations, merely acts as a constraint. We test this approach on a range of nonlinear partial differential equations. Our results show significant improvements (up to a factor of 4 reduction in the computational cost) over the traditional step size controller for a wide range of tolerances. △ Less

Submitted 17 August, 2022; v1 submitted 4 February, 2021; originally announced February 2021.

Comments: 23 pages, 14 figures, matches with the published version

Journal ref: Computers and Mathematics with Applications, 123, 59-74, 2022

A pair of UV nuclei or a compact star forming region near the active nucleus in Mrk~766?

Authors: P. P. Deka, G. C. Dewangan, K. P. Singh, J. Postma

Abstract: We report the discovery of a bright, compact ultraviolet source at a projected separation of 1.1~kpc from the known active galactic nucleus (AGN) in Mrk~766 based on Astrosat/UVIT observations. We perform radial profile analysis and derive the UV flux almost free from the nearby contaminating sources. The new source is about 2.5 and 5.6 times fainter than the AGN in the far and near UV bands. The… ▽ More We report the discovery of a bright, compact ultraviolet source at a projected separation of 1.1~kpc from the known active galactic nucleus (AGN) in Mrk~766 based on Astrosat/UVIT observations. We perform radial profile analysis and derive the UV flux almost free from the nearby contaminating sources. The new source is about 2.5 and 5.6 times fainter than the AGN in the far and near UV bands. The two sources appear as a pair of nuclei in Mrk~766. We investigate the nature of the new source based on the UV flux ratio, X-ray and optical emission. The new source is highly unlikely to be another accreting supermassive black hole in Mrk~766 as it lacks X-ray emission. We find that the UV/Optical flux of the new source measured at four different bands closely follow the shape of the template spectrum of starburst galaxies. This strongly suggests that the new source is a compact star-forming region. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Revisit of non-linear Landau damping for electrostatic instability driven by blazar-induced pair beams

Authors: S. Vafin, P. J. Deka, M. Pohl, A. Bohdan

Abstract: We revisit the effect of non-linear Landau (NL) damping on the electrostatic instability of blazar-induced pair beams, using a realistic pair-beam distribution. We employ a simplified 2D model in ${\bf k}$-space to study the evolution of the electric-field spectrum and to calculate the relaxation time of the beam. We demonstrate that the 2D model is an adequate representation of the 3D physics. We… ▽ More We revisit the effect of non-linear Landau (NL) damping on the electrostatic instability of blazar-induced pair beams, using a realistic pair-beam distribution. We employ a simplified 2D model in ${\bf k}$-space to study the evolution of the electric-field spectrum and to calculate the relaxation time of the beam. We demonstrate that the 2D model is an adequate representation of the 3D physics. We find that non-linear Landau damping, once it operates efficiently, transports essentially the entire wave energy to small wavenumbers where wave driving is weak or absent. The relaxation time also strongly depends on the IGM temperature, $T_\mathrm{IGM}$, and for $T_\mathrm{IGM}\ll10$ eV, and in the absence of any other damping mechanism, the relaxation time of the pair beam is longer than the inverse Compton (IC) scattering time. The weak late-time beam energy losses arise from the accumulation of wave energy at small $k$, that non-linearly drains the wave energy at the resonant $\mathbf{k}$ of the pair-beam instability. Any other dissipation process operating at small $k$ would reduce that wave-energy drain and hence lead to stronger pair-beam energy losses. As an example, collisions reduce the relaxation time by an order of magnitude, although their rate is very small. Other non-linear processes, such as the modulation instability, could provide additional damping of the non-resonant waves and dramatically reduce the relaxation time of the pair beam. An accurate description of the spectral evolution of the electrostatic waves is crucial for calculating the relaxation time of the pair beam. △ Less

Submitted 8 July, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Journal ref: ApJ, 873, 10, 2019

Multifaceted nonlinear dynamics in $\mathcal{PT}$-symmetric coupled Liénard oscillators

Authors: Jyoti Prasad Deka, Arjunan Govindarajan, Manas Kulkarni, Amarendra K. Sarma

Abstract: We propose a generalized parity-time ($\mathcal{PT}$) -symmetric Liénard oscillator with two different orders of nonlinear position-dependent dissipation. We study the stability of the stationary states by using the eigenvalues of Jacobian and evaluate the stability threshold thereafter. In the first order nonlinear damping model, we discover that the temporal evolution of both gain and lossy osci… ▽ More We propose a generalized parity-time ($\mathcal{PT}$) -symmetric Liénard oscillator with two different orders of nonlinear position-dependent dissipation. We study the stability of the stationary states by using the eigenvalues of Jacobian and evaluate the stability threshold thereafter. In the first order nonlinear damping model, we discover that the temporal evolution of both gain and lossy oscillators attains a complete convergence towards the stable stationary state leading to the emergence of oscillation and amplitude deaths. Also, the system displays a remarkable manifestation of transient chaos in the lossy oscillator while the gain counterpart exhibits blow-up dynamics for certain choice of initial conditions and control parameters. Employing an external driving force on the loss oscillator, we find that the blow-up dynamics can be controlled and a pure aperiodic state is achievable. On the other hand, the second order nonlinear damping model yields a completely different dynamics on contrary to the first order where the former reveals a conventional quasi-periodic route to chaos upon decreasing the natural frequency of both gain and loss oscillators. An electronic circuit scheme for the experimental realization of the proposed system has also been put forward. △ Less

Submitted 11 September, 2019; v1 submitted 25 December, 2018; originally announced December 2018.

NAIR: Novel Astronomical Instrumentation through photonic Reformatting

Authors: Robert J. Harris, Jan Tepper, John J. Davenport, Ettore Pedretti, Dionne M. Haynes, Philipp Hottinger, Theodoros Anagnos, Abani Shankar Nayak, Yohana Herrero Alonso, Pranab Jyoti Deka, Stefano Minardi, Andreas Quirrenbach, Lucas Labadie, Roger Haynes

Abstract: The project "Novel Astronomical Instrumentation through photonic Reformatting" is a DFG-funded collaboration to exploit the recognized potential of photonics solutions for a radically new approach to astronomical instrumentation for optical/infrared high precision spectroscopy and high angular resolution imaging. We present a project overview and initial development results from our Adaptive Optic… ▽ More The project "Novel Astronomical Instrumentation through photonic Reformatting" is a DFG-funded collaboration to exploit the recognized potential of photonics solutions for a radically new approach to astronomical instrumentation for optical/infrared high precision spectroscopy and high angular resolution imaging. We present a project overview and initial development results from our Adaptive Optics-photonic test bed, Ultrafast Laser Inscribed waveguides for interferometric beam combination and 3D printing structures for astronomical instrumentation. The project is expected to lead to important technological breakthroughs facilitating uniquely functionality and technical solutions for the next generation of instrumentation. △ Less

Submitted 12 July, 2018; originally announced July 2018.

Comments: 15 pages, 11 figures, SPIE Astronomical Telescopes + Instrumentation

Chaotic Dynamics and Optical Power Saturation in Parity-Time (PT) Symmetric Double Ring Resonator

Authors: Jyoti Prasad Deka, Amarendra K. Sarma

Abstract: We report emergence of saturation and chaotic spiking of optical power in a double ring resonator with balanced loss and gain, obeying the so-called parity-time symmetry. We have modeled the system using a discrete-time iterative equation known as the Ikeda Map. In the linear regime, evolution of optical power in the system shows power saturation behavior below the PT threshold and exponential blo… ▽ More We report emergence of saturation and chaotic spiking of optical power in a double ring resonator with balanced loss and gain, obeying the so-called parity-time symmetry. We have modeled the system using a discrete-time iterative equation known as the Ikeda Map. In the linear regime, evolution of optical power in the system shows power saturation behavior below the PT threshold and exponential blow-up above the PT threshold. We found that in the unbroken PT regime, optical power saturation occurs owing to the existence of stable stationary states, which lies on the surface of 4-dimensional hypersphere. Inclusion of Kerr nonlinearity into our model leads to the emergence of a stable, chaotic and divergent region in the parameter basin for period-1 cycle. A closer inspection into the system shows us that the largest Lyapunov exponent blows up in the divergent region. It is found that the existence of high non-negative largest Lyapunov exponent causes chaotic spiking of optical power in the resonators. △ Less

Submitted 3 January, 2018; v1 submitted 29 December, 2017; originally announced December 2017.

Highly Amplified Light Transmission in Parity-Time Symmetric Multilayered Structure

Authors: Jyoti Prasad Deka, Amarendra K. Sarma

Abstract: We propose a parity-time symmetric dielectric-nanofilm-dielectric multilayered structure that could facilitate highly amplified transmission of optical power in the infrared spectrum. We have theoretically studied our model using the transfer matrix formalism. The reflection and the transmission coefficients of the S-matrix are evaluated. The theoretical results are validated by FDTD numerical sim… ▽ More We propose a parity-time symmetric dielectric-nanofilm-dielectric multilayered structure that could facilitate highly amplified transmission of optical power in the infrared spectrum. We have theoretically studied our model using the transfer matrix formalism. The reflection and the transmission coefficients of the S-matrix are evaluated. The theoretical results are validated by FDTD numerical simulation. We have shown how the thickness of the layers and the gain/loss coefficient of the active layers could generate spectral singularities in the S-matrix and how these singularities could be exploited to achieve amplified transmission of a single wavelength through the structure. △ Less

Submitted 3 January, 2018; v1 submitted 8 June, 2017; originally announced June 2017.

Journal ref: Applied Optics Vol. 57, no. 5, p. 1119, 2018

Perturbative Dynamics of Stationary States in Nonlinear Parity-Time Symmetric Coupler

Authors: Jyoti Prasad Deka, Amarendra K. Sarma

Abstract: We investigate the nonlinear parity-time (PT) symmetric coupler from a dynamical perspective. As opposed to linear PT-coupler where the PT threshold dictates the evolutionary characteristics of optical power in the two waveguides, in a nonlinear coupler, the PT threshold governs the existence of stationary points. We have found that the stability of the ground state undergoes a phase transition wh… ▽ More We investigate the nonlinear parity-time (PT) symmetric coupler from a dynamical perspective. As opposed to linear PT-coupler where the PT threshold dictates the evolutionary characteristics of optical power in the two waveguides, in a nonlinear coupler, the PT threshold governs the existence of stationary points. We have found that the stability of the ground state undergoes a phase transition when the gain/loss coefficient is increased from zero to beyond the PT threshold. Moreover, we found that instabilities in initial conditions can lead to aperiodic oscillations as well as exponential growth and decay of optical power. At the PT threshold, we observed the existence of a stable attractor under the influence of fluctuating gain/loss coefficient. Phase plane analysis has shown us the presence of a toroidal chaotic attractor. The chaotic dynamics can be controlled by a judicious choice of the waveguide parameters. △ Less

Submitted 7 December, 2016; v1 submitted 11 May, 2016; originally announced May 2016.

Controllable chaotic dynamics in a nonlinear fiber ring resonators with balanced gain and loss

Authors: Jyoti P. Deka, Samit Kumar Gupta, Amarendra K. Sarma

Abstract: We show the possibility of controlling the dynamical behavior of a single fiber ring (SFR) resonator system with the fiber being an amplified (gain) channel and the ring being attenuated (loss) nonlinear dielectric medium. Our model is based on the simple alterations in the parity time symmetric synthetic coupler structures proposed recently [A. Regensburger et al., Nature 488, 167 (2012)]. The sy… ▽ More We show the possibility of controlling the dynamical behavior of a single fiber ring (SFR) resonator system with the fiber being an amplified (gain) channel and the ring being attenuated (loss) nonlinear dielectric medium. Our model is based on the simple alterations in the parity time symmetric synthetic coupler structures proposed recently [A. Regensburger et al., Nature 488, 167 (2012)]. The system has been modeled using the transfer matrix formalism. We find that this results in a dynamically controllable algorithm for the chaotic dynamics inherent in the system. We have also shown the dependence of the period doubling point on the input amplitude, emphasizing on the dynamical aspects. Moreover, the fact that the resonator essentially plays the role of a damped harmonic oscillator has been elucidated with the non-zero intensity inside the resonator due to constant influx of input light. △ Less

Submitted 5 September, 2016; v1 submitted 21 January, 2015; originally announced January 2015.

Journal ref: Nonlinear Dynamics 87, 1121 (2017)

Nonlinear Parity-Time (PT) symmetric closed-form optical quadrimer waveguides: Attractor perspective

Authors: Samit Kumar Gupta, Jyoti Prasad Deka, Amarendra K. Sarma

Abstract: We report a study on a closed-form nonlinear parity-time symmetric optical quadrimer waveguides system with a specific coupling scheme. The system yields power saturation behavior in the modes, which may be attributed to the inherent attractor in the system. A detailed analysis has been provided to confirm the attractor aspect of the system. This work also addresses a crucial issue regarding choic… ▽ More We report a study on a closed-form nonlinear parity-time symmetric optical quadrimer waveguides system with a specific coupling scheme. The system yields power saturation behavior in the modes, which may be attributed to the inherent attractor in the system. A detailed analysis has been provided to confirm the attractor aspect of the system. This work also addresses a crucial issue regarding choice of initial conditions while carrying out numerical simulation for such systems. △ Less

Submitted 18 August, 2015; v1 submitted 23 October, 2014; originally announced October 2014.

Journal ref: The European Physical Journal D Vol. 69, p. 199 (2015)