-
Probing the connection between IceCube neutrinos and MOJAVE AGN
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
S. K. Agarwalla,
J. A. Aguilar,
M. Ahlers,
J. M. Alameddine,
N. M. Amin,
K. Andeen,
C. Argüelles,
Y. Ashida,
S. Athanasiadou,
L. Ausborm,
S. N. Axani,
X. Bai,
A. Balagopal V.,
M. Baricevic,
S. W. Barwick,
S. Bash,
V. Basu,
R. Bay,
J. J. Beatty,
J. Becker Tjus,
J. Beise,
C. Bellenghi
, et al. (399 additional authors not shown)
Abstract:
Active Galactic Nuclei (AGN) are prime candidate sources of the high-energy, astrophysical neutrinos detected by IceCube. This is demonstrated by the real-time multi-messenger detection of the blazar TXS 0506+056 and the recent evidence of neutrino emission from NGC 1068 from a separate time-averaged study. However, the production mechanism of the astrophysical neutrinos in AGN is not well establi…
▽ More
Active Galactic Nuclei (AGN) are prime candidate sources of the high-energy, astrophysical neutrinos detected by IceCube. This is demonstrated by the real-time multi-messenger detection of the blazar TXS 0506+056 and the recent evidence of neutrino emission from NGC 1068 from a separate time-averaged study. However, the production mechanism of the astrophysical neutrinos in AGN is not well established which can be resolved via correlation studies with photon observations. For neutrinos produced due to photohadronic interactions in AGN, in addition to a correlation of neutrinos with high-energy photons, there would also be a correlation of neutrinos with photons emitted at radio wavelengths. In this work, we perform an in-depth stacking study of the correlation between 15 GHz radio observations of AGN reported in the MOJAVE XV catalog, and ten years of neutrino data from IceCube. We also use a time-dependent approach which improves the statistical power of the stacking analysis. No significant correlation was found for both analyses and upper limits are reported. When compared to the IceCube diffuse flux, at 100 TeV and for a spectral index of 2.5, the upper limits derived are $\sim3\%$ and $\sim9\%$ for the time-averaged and time-dependent case, respectively.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Time Non-locality in Dark Matter and LSS
Authors:
Arhum Ansari,
Arka Banerjee,
Sachin Jain,
Shaunak Padhyegurjar
Abstract:
We explore the intriguing phenomenon of time non-locality in the evolution of dark matter and Large Scale Structure (LSS). Recently in\,\cite{Donath:2023sav}, it was shown that time non-locality emerges in bias tracer fluctuations, which are $SO(3)$ scalars in real space, at fifth order in the perturbation expansion in dark matter overdensity. We demonstrate that by breaking the symmetry down to…
▽ More
We explore the intriguing phenomenon of time non-locality in the evolution of dark matter and Large Scale Structure (LSS). Recently in\,\cite{Donath:2023sav}, it was shown that time non-locality emerges in bias tracer fluctuations, which are $SO(3)$ scalars in real space, at fifth order in the perturbation expansion in dark matter overdensity. We demonstrate that by breaking the symmetry down to $SO(2)$, which is the case whenever line-of-sight effects become important, such as for flux fluctuations in the Lyman $α$ forest, the temporal non-locality appears at the third order in expansion. Additionally, within the framework of EFTofLSS, we demonstrate that time non-locality manifests in the effective stress tensor of dark matter, which is a second rank tensor under $SO(3)$ transformations, again at the third order in dark matter overdensity. Furthermore, we highlight the effectiveness of the standard $Π$ basis\,\cite{Mirbabayi:2014zca} in handling time non-local operators.
△ Less
Submitted 24 June, 2024;
originally announced June 2024.
-
Search for neutrino emission from hard X-ray AGN with IceCube
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
S. K. Agarwalla,
J. A. Aguilar,
M. Ahlers,
J. M. Alameddine,
N. M. Amin,
K. Andeen,
C. Argüelles,
Y. Ashida,
S. Athanasiadou,
L. Ausborm,
S. N. Axani,
X. Bai,
A. Balagopal V.,
M. Baricevic,
S. W. Barwick,
S. Bash,
V. Basu,
R. Bay,
J. J. Beatty,
J. Becker Tjus,
J. Beise,
C. Bellenghi
, et al. (401 additional authors not shown)
Abstract:
Active Galactic Nuclei (AGN) are promising candidate sources of high-energy astrophysical neutrinos since they provide environments rich in matter and photon targets where cosmic ray interactions may lead to the production of gamma rays and neutrinos. We searched for high-energy neutrino emission from AGN using the $\textit{Swift}$-BAT Spectroscopic Survey (BASS) catalog of hard X-ray sources and…
▽ More
Active Galactic Nuclei (AGN) are promising candidate sources of high-energy astrophysical neutrinos since they provide environments rich in matter and photon targets where cosmic ray interactions may lead to the production of gamma rays and neutrinos. We searched for high-energy neutrino emission from AGN using the $\textit{Swift}$-BAT Spectroscopic Survey (BASS) catalog of hard X-ray sources and 12 years of IceCube muon track data. First, upon performing a stacked search, no significant emission was found. Second, we searched for neutrinos from a list of 43 candidate sources and found an excess from the direction of two sources, Seyfert galaxies NGC 1068 and NGC 4151. We observed NGC 1068 at flux $φ_{ν_μ+\barν_μ}$ = $4.02_{-1.52}^{+1.58} \times 10^{-11}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ normalized at 1 TeV, with power-law spectral index, $γ$ = 3.10$^{+0.26}_{-0.22}$, consistent with previous IceCube results. The observation of a neutrino excess from the direction of NGC 4151 is at a post-trial significance of 2.9$σ$. If interpreted as an astrophysical signal, the excess observed from NGC 4151 corresponds to a flux $φ_{ν_μ+\barν_μ}$ = $1.51_{-0.81}^{+0.99} \times 10^{-11}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ normalized at 1 TeV and $γ$ = 2.83$^{+0.35}_{-0.28}$.
△ Less
Submitted 12 June, 2024; v1 submitted 10 June, 2024;
originally announced June 2024.
-
Spectroscopic Investigation of Nebular Gas (SING): Instrument Design, Assembly and Calibration
Authors:
Bharat Chandra P,
Binukumar G. Nair,
Shubham Jankiram Ghatul,
Shubhangi Jain,
S. Sriram,
Mahesh Babu S.,
Rekhesh Mohan,
Margarita Safonova,
Jayant Murthy,
Mikhail Sachkov
Abstract:
The Spectroscopic Investigation of Nebular Gas (SING) is a near-ultraviolet (NUV) low-resolution spectrograph payload designed to operate in the NUV range, 1400 $\unicode{x212B}$ -- 2700 $\unicode{x212B}$, from a stable space platform. SING telescope has a primary aperture of 298 mm, feeding the light to the long-slit UV spectrograph. SING has a field of view (FOV) of 1$^{\circ}$, achieving a spat…
▽ More
The Spectroscopic Investigation of Nebular Gas (SING) is a near-ultraviolet (NUV) low-resolution spectrograph payload designed to operate in the NUV range, 1400 $\unicode{x212B}$ -- 2700 $\unicode{x212B}$, from a stable space platform. SING telescope has a primary aperture of 298 mm, feeding the light to the long-slit UV spectrograph. SING has a field of view (FOV) of 1$^{\circ}$, achieving a spatial resolution of 1.33 arc minute and spectral resolution of 3.7 $\unicode{x212B}$ ($R\sim600$) at the central wavelength. SING employs a micro-channel plate (MCP) with a CMOS readout-based photon-counting detector. The instrument is designed to observe diffuse sources such as nebulae, supernova remnants, and the interstellar medium (ISM) to understand their chemistry. SING was selected by the United Nations Office for Outer Space Affairs to be hosted on the Chinese Space Station. The instrument will undergo qualification tests as per the launch requirements. In this paper, we describe the hardware design, optomechanical assembly, and calibration of the instrument.
△ Less
Submitted 26 April, 2024;
originally announced April 2024.
-
Inflationary non-Gaussianities in alpha vacua and consistency with conformal symmetries
Authors:
Arhum Ansari,
Pinak Banerjee,
Prateksh Dhivakar,
Sachin Jain,
Nilay Kundu
Abstract:
We study the conformal invariance of inflationary non-Gaussianities associated with scalar fluctuations in a non-Bunch-Davies initial state, known as the $α$-vacuum, in single-field slow-roll inflation. The $α$-vacuum is a one-parameter family of states, including the Bunch-Davies one, that preserves the conformal symmetry of inflationary dynamics in a nearly de-Sitter space-time. Working within t…
▽ More
We study the conformal invariance of inflationary non-Gaussianities associated with scalar fluctuations in a non-Bunch-Davies initial state, known as the $α$-vacuum, in single-field slow-roll inflation. The $α$-vacuum is a one-parameter family of states, including the Bunch-Davies one, that preserves the conformal symmetry of inflationary dynamics in a nearly de-Sitter space-time. Working within the leading slow-roll approximation, we compute the four-point scalar correlator (the trispectrum) in $α$-vacuum using the in-in formalism. We check that the conformal Ward identities are met between the three and four-point scalar $α$-vacua correlators. Surprisingly, this contrasts the previously reported negative result of the Ward identities being violated between the two and the three-point correlators. We have also extended the wave-functional method, previously used for correlators with Bunch-Davies initial condition, to compute the three and four-point scalar correlators in $α$-vacua. The results obtained from the wave-function method match the corresponding in-in results, adding further justification to our check of Ward identities with $α$-vacua correlators.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
The McNish and Lincoln Solar Activity Predictions: The Method and its Performance
Authors:
Frédéric Clette,
Shantanu Jain,
Tatiana Podladchikova
Abstract:
The McNish and Lincoln (ML) method, introduced in 1949, was one of the first attempts to produce mid-term forecasts of solar activity, up to 12 months ahead. However, it has been poorly described and evaluated in the past literature, in particular its actual operational implementation by NOAA.
Here, we reconstruct the exact formulation of the method, as it was applied since the early 1970s, and…
▽ More
The McNish and Lincoln (ML) method, introduced in 1949, was one of the first attempts to produce mid-term forecasts of solar activity, up to 12 months ahead. However, it has been poorly described and evaluated in the past literature, in particular its actual operational implementation by NOAA.
Here, we reconstruct the exact formulation of the method, as it was applied since the early 1970s, and we provide a full mathematical derivation of the prediction errors. For bench-marking the method, we also produce monthly predictions over the past 190 years, from 1833 (Cycle 8) to 2023 (Cycle 25), and develop statistics of the differences between the predictions and the observed 13-month smoothed sunspot number (SSN) time series, according to the phase in the solar cycle.
Our analysis shows that the ML method is heavily constrained because it is primarily based on the mean of all past cycles, which imposes a fixed amplitude and length and suffers from a temporal smearing that grows towards the end of the solar cycle. We find that predictions are completely unreliable in the first 12 months of the cycle, and over the last two years preceding the ending minimum (around 130 months), and beyond this minimum. By contrast, in the course of the cycle (months 18 to 65), ML predictions prove to be reliable over a time range of up to 50 months (4.2 years), thus much longer than the 12-month conventional range used so far. However, we find that predictions then suffer from systematic under-(over-)estimates for cycles that have a higher (lower) amplitude than the base mean cycle.
Overall, we conclude that although the ML method provides valid prediction errors, it suffers from strong limitations, with very little room for improvement, as it indifferently merges all past cycles into a single fixed statistics.
△ Less
Submitted 14 February, 2024;
originally announced February 2024.
-
Development Of Raspberry Pi-based Processing Unit for UV Photon-Counting Detectors
Authors:
Bharat Chandra P,
Binukumar G.,
Shubham Janakiram,
Mahesh Babu,
Shubhangi Jain,
Richa Rai,
Rekhesh Mohan,
Margarita Safonova,
Jayant Murthy
Abstract:
In ultraviolet (UV) astronomical observations, photons from the sources are very few compared to the visible or infrared (IR) wavelength ranges. Detectors operating in the UV usually employ a photon-counting mode of operation. These detectors usually have an image intensifier sensitive to UV photons and a readout mechanism that employs photon counting. The development of readouts for these detecto…
▽ More
In ultraviolet (UV) astronomical observations, photons from the sources are very few compared to the visible or infrared (IR) wavelength ranges. Detectors operating in the UV usually employ a photon-counting mode of operation. These detectors usually have an image intensifier sensitive to UV photons and a readout mechanism that employs photon counting. The development of readouts for these detectors is resource-intensive and expensive. In this paper, we describe the development of a low-cost UV photon-counting detector processing unit that employs a Raspberry Pi with its in built readout to perform the photon-counting operation. Our system can operate in both 3x3 and 5x5 window modes at 30 frames per sec (fps), where 5x5 window mode also enables the provision of detection of double events. The system can be built quickly from readily available custom-off-the-shelf (COTS) components and is thus used in inexpensive CubeSats or small satellite missions. This low-cost solution promises to broaden access to UV observations, advancing research possibilities in space-based astronomy.
△ Less
Submitted 5 January, 2024; v1 submitted 2 January, 2024;
originally announced January 2024.
-
Coronal dimmings as indicators of early CME propagation direction
Authors:
Shantanu Jain,
Tatiana Podladchikova,
Galina Chikunova,
Karin Dissauer,
Astrid M. Veronig
Abstract:
Coronal mass ejections (CMEs) are solar eruptions of plasma and magnetic fields that significantly impact Space Weather, causing disruptions in technological systems and potential damage to power grids when directed towards Earth. Traditional coronagraphs along the Sun-Earth line struggle to precisely track the early evolution of Earth-directed CMEs. Coronal dimmings, localized reductions in extre…
▽ More
Coronal mass ejections (CMEs) are solar eruptions of plasma and magnetic fields that significantly impact Space Weather, causing disruptions in technological systems and potential damage to power grids when directed towards Earth. Traditional coronagraphs along the Sun-Earth line struggle to precisely track the early evolution of Earth-directed CMEs. Coronal dimmings, localized reductions in extreme-ultraviolet (EUV) and soft X-ray emissions, are key indicators of CMEs in the low corona, resulting from mass loss and expansion during the eruption. This study introduces a novel method, DIRECD (Dimming InfeRred Estimate of CME Direction), to estimate the early propagation direction of CMEs based on the expansion of coronal dimmings. The approach involves 3D simulations of CMEs using a geometric cone model, exploring parameters like width, height, source location, and deflection from the radial direction. The dominant direction of dimming evolution is then determined, and an inverse problem is solved to reconstruct an ensemble of CME cones at various heights, widths, and deflections. By comparing the CME orthogonal projections onto the solar sphere with the dimming geometry, the 3D CME direction is derived.Validated through case studies on October 1, 2011, and September 6, 2011, the DIRECD method reveals the early propagation directions of CMEs. The CME on October 1, 2011, predominantly expands towards the South-East, while the CME on September 6, 2011, inclines towards the North-West. These findings align with previous studies using multi-viewpoint coronagraphic observations. The study demonstrates the utility of coronal dimming information for early CME direction estimation, providing valuable data for space weather forecasting and mitigating potential adverse impacts on Earth before observation in coronographs' field-of-view.
△ Less
Submitted 11 December, 2023; v1 submitted 23 November, 2023;
originally announced November 2023.
-
The motivation for flexible star-formation histories from spatially resolved scales within galaxies
Authors:
Shweta Jain,
Sandro Tacchella,
Moein Mosleh
Abstract:
The estimation of galaxy stellar masses depends on the assumed prior of the star-formation history (SFH) and spatial scale of the analysis (spatially resolved versus integrated scales). In this paper, we connect the prescription of the SFH in the Spectral Energy Distribution (SED) fitting to spatially resolved scales ($\sim\mathrm{kpc}$) to shed light on the systematics involved when estimating st…
▽ More
The estimation of galaxy stellar masses depends on the assumed prior of the star-formation history (SFH) and spatial scale of the analysis (spatially resolved versus integrated scales). In this paper, we connect the prescription of the SFH in the Spectral Energy Distribution (SED) fitting to spatially resolved scales ($\sim\mathrm{kpc}$) to shed light on the systematics involved when estimating stellar masses. Specifically, we fit the integrated photometry of $\sim970$ massive (log (M$_{\star}$/M$_{\odot}) = 9.8-11.6$), intermediate redshift ($z=0.5-2.0$) galaxies with $\texttt{Prospector}$, assuming both exponentially declining tau model and flexible SFHs. We complement these fits with the results of spatially resolved SFH estimates obtained by pixel-by-pixel SED fitting, which assume tau models for individual pixels. These spatially resolved SFHs show a large diversity in shapes, which can largely be accounted for by the flexible SFHs with $\texttt{Prospector}$. The differences in the stellar masses from those two approaches are overall in good agreement (average difference of $\sim 0.07$ dex). Contrarily, the simpler tau model SFHs typically miss the oldest episode of star formation, leading to an underestimation of the stellar mass by $\sim 0.3$ dex. We further compare the derived global specific star-formation rate (sSFR), the mass-weighted stellar age (t$_{50}$), and the star-formation timescale ($τ_{\mathrm{SF}}$) obtained from the different SFH approaches. We conclude that the spatially resolved scales within galaxies motivate a flexible SFH on global scales to account for the diversity of SFHs and counteract the effects of outshining of older stellar populations by younger ones.
△ Less
Submitted 27 October, 2023;
originally announced October 2023.
-
The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
Authors:
The Astropy Collaboration,
Adrian M. Price-Whelan,
Pey Lian Lim,
Nicholas Earl,
Nathaniel Starkman,
Larry Bradley,
David L. Shupe,
Aarya A. Patil,
Lia Corrales,
C. E. Brasseur,
Maximilian Nöthe,
Axel Donath,
Erik Tollerud,
Brett M. Morris,
Adam Ginsburg,
Eero Vaher,
Benjamin A. Weaver,
James Tocknell,
William Jamieson,
Marten H. van Kerkwijk,
Thomas P. Robitaille,
Bruce Merry,
Matteo Bachetti,
H. Moritz Günther,
Thomas L. Aldcroft
, et al. (111 additional authors not shown)
Abstract:
The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as…
▽ More
The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project.
△ Less
Submitted 28 June, 2022;
originally announced June 2022.
-
Maximal growth rate of the ascending phase of a sunspot cycle for predicting its amplitude
Authors:
Tatiana Podladchikova,
Shantanu Jain,
Astrid M. Veronig,
Olga Sutyrina,
Mateja Dumbovic,
Frederic Clette,
Werner Poetzi
Abstract:
Forecasting the solar cycle amplitude is important for a better understanding of the solar dynamo as well as for many space weather applications. We demonstrated a steady relationship between the maximal growth rate of sunspot activity in the ascending phase of a cycle and the subsequent cycle amplitude on the basis of four data sets of solar activity indices: total sunspot numbers, hemispheric su…
▽ More
Forecasting the solar cycle amplitude is important for a better understanding of the solar dynamo as well as for many space weather applications. We demonstrated a steady relationship between the maximal growth rate of sunspot activity in the ascending phase of a cycle and the subsequent cycle amplitude on the basis of four data sets of solar activity indices: total sunspot numbers, hemispheric sunspot numbers from the new catalogue from 1874 onwards, total sunspot areas, and hemispheric sunspot areas. For all the data sets, a linear regression based on the maximal growth rate precursor shows a significant correlation. Validation of predictions for cycles 1-24 shows high correlations between the true and predicted cycle amplitudes reaching r = 0.93 for the total sunspot numbers. The lead time of the predictions varies from 2 to 49 months, with a mean value of 21 months. Furthermore, we demonstrated that the sum of maximal growth rate indicators determined separately for the north and the south hemispheric sunspot numbers provides more accurate predictions than that using total sunspot numbers. The advantages reach 27% and 11% on average in terms of rms and correlation coefficient, respectively. The superior performance is also confirmed with hemispheric sunspot areas with respect to total sunspot areas. The maximal growth rate of sunspot activity in the ascending phase of a solar cycle serves as a reliable precursor of the subsequent cycle amplitude. Furthermore, our findings provide a strong foundation for supporting regular monitoring, recording, and predictions of solar activity with hemispheric sunspot data, which capture the asymmetric behaviour of the solar activity and solar magnetic field and enhance solar cycle prediction methods.
△ Less
Submitted 25 June, 2022;
originally announced June 2022.
-
Effect of damped oscillations in the inflationary potential
Authors:
Akhil Antony,
Shweta Jain
Abstract:
We investigate the effect of damped oscillations on a nearly flat inflationary potential and the features they produce in the power-spectrum and bi-spectrum. We compare the model with the Planck data using Plik unbinned and CamSpec clean likelihood and we are able to obtain noticeable improvement in fit compared to the power-law $Λ$CDM model. We are able to identify three plausible candidates each…
▽ More
We investigate the effect of damped oscillations on a nearly flat inflationary potential and the features they produce in the power-spectrum and bi-spectrum. We compare the model with the Planck data using Plik unbinned and CamSpec clean likelihood and we are able to obtain noticeable improvement in fit compared to the power-law $Λ$CDM model. We are able to identify three plausible candidates each for the two likelihoods used. We find that the best-fit to Plik and CamSpec likelihoods match closely to each other. The improvement comes from various possible outliers at the intermediate to small scales. We also compute the bi-spectrum for the best-fits. At all limits, the amplitude of bi-spectrum, $f_{NL}$ is oscillatory in nature and its peak value is determined by the amplitude and frequency of the oscillations in the potential, as expected. We find that the bi-spectrum consistency relation strictly holds at all scales in all the best-fit candidates.
△ Less
Submitted 17 March, 2023; v1 submitted 13 October, 2021;
originally announced October 2021.
-
Hemispheric sunspot numbers 1874--2020
Authors:
Astrid M. Veronig,
Shantanu Jain,
Tatiana Podladchikova,
Werner Poetzi,
Frederic Clette
Abstract:
We create a continuous series of daily and monthly hemispheric sunspot numbers (HSNs) from 1874 to 2020, which will be continuously expanded in the future with the HSNs provided by SILSO. Based on the available daily measurements of hemispheric sunspot areas from 1874 to 2016 from Greenwich Royal Observatory and NOAA, we derive the relative fractions of the northern and southern activity. These fr…
▽ More
We create a continuous series of daily and monthly hemispheric sunspot numbers (HSNs) from 1874 to 2020, which will be continuously expanded in the future with the HSNs provided by SILSO. Based on the available daily measurements of hemispheric sunspot areas from 1874 to 2016 from Greenwich Royal Observatory and NOAA, we derive the relative fractions of the northern and southern activity. These fractions are applied to the international sunspot number (ISN) to derive the HSNs. This method and obtained data are validated against published HSNs for the period 1945--2020. We provide a continuous data series and catalogue of daily, monthly mean, and 13-month smoothed monthly mean HSNs for the time range 1874--2020 that are consistent with the newly calibrated ISN. Validation of the reconstructed HSNs against the direct data available since 1945 reveals a high level of consistency, with a correlation of r=0.94 (0.97) for the daily (monthly) data. The cumulative hemispheric asymmetries for cycles 12-24 give a mean value of 16%, with no obvious pattern in north-south predominance over the cycle evolution. The strongest asymmetry occurs for cycle no. 19, in which the northern hemisphere shows a cumulated predominance of 42%. The phase shift between the peaks of solar activity in the two hemispheres may be up to 28 months, with a mean absolute value of 16.4 months. The phase shifts reveal an overall asymmetry of the northern hemisphere reaching its cycle maximum earlier (in 10 out of 13 cases). Relating the ISN and HSN peak growth rates during the cycle rise phase with the cycle amplitude reveals higher correlations when considering the two hemispheres individually, with r = 0.9. Our findings demonstrate that empirical solar cycle prediction methods can be improved by investigating the solar cycle dynamics in terms of the hemispheric sunspot numbers.
△ Less
Submitted 1 July, 2021;
originally announced July 2021.
-
Forbidden atomic oxygen emissions in the Martian dayside upper atmosphere
Authors:
Susarla Raghuram,
Sonal Kumar Jain,
Anil Bhardwaj
Abstract:
Recently, Nadir and Occultation for Mars Discovery (NOMAD) ultraviolet and visible spectrometer instrument on board the European Space Agency's ExoMars Trace Gas Orbiter (TGO) simultaneously measured the limb emission intensities for both [OI] 2972 and 5577 Å (green) emissions in the dayside of Martian upper atmosphere. We aim to explore the photochemistry of all these forbidden atomic oxygen emis…
▽ More
Recently, Nadir and Occultation for Mars Discovery (NOMAD) ultraviolet and visible spectrometer instrument on board the European Space Agency's ExoMars Trace Gas Orbiter (TGO) simultaneously measured the limb emission intensities for both [OI] 2972 and 5577 Å (green) emissions in the dayside of Martian upper atmosphere. We aim to explore the photochemistry of all these forbidden atomic oxygen emissions ([OI] 2972, 5577, 6300, 6464 Å) in the Martian daylight upper atmosphere and suitable conditions for the simultaneous detection of these emissions lines in the dayside visible spectra. A photochemical model is developed to study the production and loss processes of O(1S) and O(1D) by incorporating various chemical reactions of different O-bearing species in the upper atmosphere of Mars. By reducing Fox (2004) modelled neutral density profiles by a factor of 2, the calculated limb intensity profiles for [OI] 5577 and 2972 Å emissions are found to be consistent with the NOMAD-TGO observations. In this case, at altitudes below 120 km, our modelled limb intensity for [OI] 6300 Å emission is smaller by a factor 2 to 5 compared to that of NOMAD-TGO observation for [OI] 2972 Å emission, and above this distance it is comparable with the upper limit of the observation. We studied various parameters which can influence the limb intensities of these atomic oxygen forbidden emission lines. Our calculated limb intensity for [OI] 6300 Å emission, when the Mars is at near perihelion and for solar maximum condition, suggests that all these forbidden emissions should be observable in the NOMAD-TGO visible spectra taken on the dayside of Martian upper atmosphere. More simultaneous observations of forbidden atomic oxygen emission lines will help to understand the photochemical processes of oxygen-bearing species in the dayside Martian upper atmosphere.
△ Less
Submitted 10 December, 2020;
originally announced December 2020.
-
EinsteinPy: A Community Python Package for General Relativity
Authors:
Shreyas Bapat,
Ritwik Saha,
Bhavya Bhatt,
Shilpi Jain,
Akshita Jain,
Sofía Ortín Vela,
Priyanshu Khandelwal,
Jyotirmaya Shivottam,
Jialin Ma,
Gim Seng Ng,
Pratyush Kerhalkar,
Hrishikesh Sudam Sarode,
Rishi Sharma,
Manvi Gupta,
Divya Gupta,
Tushar Tyagi,
Tanmay Rustagi,
Varun Singh,
Saurabh Bansal,
Naman Tayal,
Abhijeet Manhas,
Raphael Reyna,
Gaurav Kumar,
Govind Dixit,
Ratin Kumar
, et al. (24 additional authors not shown)
Abstract:
This paper presents EinsteinPy (version 0.3), a community-developed Python package for gravitational and relativistic astrophysics. Python is a free, easy to use a high-level programming language which has seen a huge expansion in the number of its users and developers in recent years. Specifically, a lot of recent studies show that the use of Python in Astrophysics and general physics has increas…
▽ More
This paper presents EinsteinPy (version 0.3), a community-developed Python package for gravitational and relativistic astrophysics. Python is a free, easy to use a high-level programming language which has seen a huge expansion in the number of its users and developers in recent years. Specifically, a lot of recent studies show that the use of Python in Astrophysics and general physics has increased exponentially. We aim to provide a very high level of abstraction, an easy to use interface and pleasing user experience. EinsteinPy is developed keeping in mind the state of a theoretical gravitational physicist with little or no background in computer programming and trying to work in the field of numerical relativity or trying to use simulations in their research. Currently, EinsteinPy supports simulation of time-like and null geodesics and calculates trajectories in different background geometries some of which are Schwarzschild, Kerr, and KerrNewmann along with coordinate inter-conversion pipeline. It has a partially developed pipeline for plotting and visualization with dependencies on libraries like Plotly, matplotlib, etc. One of the unique features of EinsteinPy is a sufficiently developed symbolic tensor manipulation utilities which are a great tool in itself for teaching yourself tensor algebra which for many beginner students can be overwhelmingly tricky. EinsteinPy also provides few utility functions for hypersurface embedding of Schwarzschild spacetime which further will be extended to model gravitational lensing simulation.
△ Less
Submitted 22 May, 2020;
originally announced May 2020.
-
Neutrino spin-flavor oscillations in solar environment
Authors:
Sandeep Joshi,
Sudhir R. Jain
Abstract:
We study the phenomenon of neutrino spin-flavor oscillations due to solar magnetic fields. This allows us to examine how significantly the electron neutrinos produced in the solar interior undergo a resonant spin-flavor conversion. We construct analytical models for the solar magnetic field in all the three regions of the Sun. Neutrino spin-flavor oscillations in these magnetic fields are examined…
▽ More
We study the phenomenon of neutrino spin-flavor oscillations due to solar magnetic fields. This allows us to examine how significantly the electron neutrinos produced in the solar interior undergo a resonant spin-flavor conversion. We construct analytical models for the solar magnetic field in all the three regions of the Sun. Neutrino spin-flavor oscillations in these magnetic fields are examined by studying the level crossing phenomenon and comparing the two cases of zero and non-zero vacuum mixing respectively. Results from the Borexino experiment are used to place an upper limit on the magnetic field in the solar core. Related phenomena such as effects of matter on neutrino spin transitions and differences between Dirac and Majorana transitions in the solar magnetic fields are also discussed.
△ Less
Submitted 18 August, 2020; v1 submitted 22 June, 2019;
originally announced June 2019.
-
Noncyclic geometric phases and helicity transitions for neutrino oscillations in a magnetic field
Authors:
Sandeep Joshi,
Sudhir R. Jain
Abstract:
We show that neutrino spin and spin-flavor transitions involve nonvanishing geometric phases. The geometric character of neutrino spin rotation is explored by studying the neutrino spin trajectory in the projective Hilbert space representation and its relation to the geometric phase. Analytical expressions are derived for noncyclic geometric phases. Several calculations are performed for different…
▽ More
We show that neutrino spin and spin-flavor transitions involve nonvanishing geometric phases. The geometric character of neutrino spin rotation is explored by studying the neutrino spin trajectory in the projective Hilbert space representation and its relation to the geometric phase. Analytical expressions are derived for noncyclic geometric phases. Several calculations are performed for different cases of rotating and nonrotating magnetic fields in the context of solar neutrinos and neutrinos produced inside neutron stars. Also the effects of adiabaticity, critical magnetic fields and cross boundary effects in the case of neutrinos emanating out of neutron stars are examined.
△ Less
Submitted 3 November, 2017; v1 submitted 15 March, 2017;
originally announced March 2017.
-
Production of N2 Vegard-Kaplan and Lyman-Birge-Hopfield emissions on Pluto
Authors:
Sonal Kumar Jain,
Anil Bhardwaj
Abstract:
We have developed a model to calculate the emission intensities of various vibrational transitions of N$_2$ triplet band and Lyman-Birge-Hopfield (LBH) band emissions in the dayglow of Pluto for solar minimum, moderate, and maximum conditions. The calculated overhead intensities of Vegard-Kaplan ($ A^3Σ_u^+ - X^1Σ^+_g $), First Positive ($ B^3Π_g - A^3Σ^+_u $), Second Positive (…
▽ More
We have developed a model to calculate the emission intensities of various vibrational transitions of N$_2$ triplet band and Lyman-Birge-Hopfield (LBH) band emissions in the dayglow of Pluto for solar minimum, moderate, and maximum conditions. The calculated overhead intensities of Vegard-Kaplan ($ A^3Σ_u^+ - X^1Σ^+_g $), First Positive ($ B^3Π_g - A^3Σ^+_u $), Second Positive ($ C^3Π_u - B^3Π_g $), Wu-Benesch ($W^3Δ_u - B^3Π_g$), Reverse First Positive, and LBH ($a^1Π_g$-- $X^1Σ^+_g$) bands of N$_2$ are 17 (74), 14.8 (64), 2.4 (10.8), 2.9 (12.7), 2.9 (12.5), and 2.3 (10) R, respectively, for solar minimum (maximum) condition. We have predicted the overhead and limb intensities of VK (150-190 nm) and LBH (120-190 nm) bands of N$_2$ on Pluto for the New Horizons (NH) flyby condition that can be observed by Alice: the ultraviolet imaging spectrograph also know as P-Alice. The predicted limb intensities of VK and LBH bands peak at radial distance of $ \sim$2000 km with the value of about 5 (13) and 9.5 (22) R for solar zenith angle 60$ ^\circ $ (0$ ^\circ $), respectively. We have also calculated overhead and limp intensities of few prominent transition of CO Fourth Positive bands for NH flyby condition.
△ Less
Submitted 26 December, 2014;
originally announced December 2014.
-
CO Cameron band and CO2+ UV doublet emissions in the dayglow of Venus: Role of CO in the Cameron band production
Authors:
Anil Bhardwaj,
Sonal Kumar Jain
Abstract:
Present study deals with the model calculations of CO Cameron band and CO2+ ultraviolet doublet emissions in the dayglow of Venus. The overhead and limb intensities of CO Cameron band and CO2+ UV doublet emissions are calculated for low, moderate, and high solar activity conditions. Using updated cross sections, the impact of dierent e-CO cross section for Cameron band production is estimated. The…
▽ More
Present study deals with the model calculations of CO Cameron band and CO2+ ultraviolet doublet emissions in the dayglow of Venus. The overhead and limb intensities of CO Cameron band and CO2+ UV doublet emissions are calculated for low, moderate, and high solar activity conditions. Using updated cross sections, the impact of dierent e-CO cross section for Cameron band production is estimated. The electron impact on CO is the major source mechanism of Cameron band, followed by electron and photon impact dissociation of CO2. The overhead intensities of CO Cameron band and CO2+ UV doublet emissions are about a factor of 2 higher in solar maximum than those in solar minimum condition. The effect of solar EUV flux models on the emission intensity is ~30-40% in solar minimum condition and ~2-10% in solar maximum condition. At the altitude of emission peak (135 km), the model predicted limb intensity of CO Cameron band and CO2+ UV doublet emissions in moderate (F10.7 = 130) solar activity condition is about 2400 and 300 kR, respectively, which is in agreement with the very recently published SPICAV/Venus Express observation. The model limb intensity profiles of CO Cameron band and CO2+ UV doublet are compared with SPICAV observation. We also calculated intensities of N2 Vegard-Kaplan UV bands and OI 2972 Å emissions during moderate and high solar activity conditions.
△ Less
Submitted 27 May, 2013;
originally announced May 2013.
-
Production of N2 Vegard-Kaplan and other triplet band emissions in the dayglow of Titan
Authors:
Anil Bhardwaj,
Sonal Kumar Jain
Abstract:
Recently the Cassini Ultraviolet Imaging Spectrograph has revealed the presence of N2 Vegard-Kaplan band emissions in Titan's dayglow limb observation. We present model calculations for the production of various N2 triplet states in the upper atmosphere of Titan. The Analytical Yield Spectra technique is used to calculate steady state photoelectron fluxes in Titan's atmosphere, which are in agreem…
▽ More
Recently the Cassini Ultraviolet Imaging Spectrograph has revealed the presence of N2 Vegard-Kaplan band emissions in Titan's dayglow limb observation. We present model calculations for the production of various N2 triplet states in the upper atmosphere of Titan. The Analytical Yield Spectra technique is used to calculate steady state photoelectron fluxes in Titan's atmosphere, which are in agreement with those observed by the Cassini's CAPS instrument. Considering direct electron impact excitation, inter-state cascading, and quenching effects, the population of different levels of N2 triplet states are calculated under statistical equilibrium. Densities of all vibrational levels of each triplet state and volume production rates for various triplet states are calculated in the model. Vertically integrated overhead intensities for the same date and lighting conditions as the reported by UVIS observations for N2 VK, 1P, 2P, Wu-Benesch, and Reverse First Positive bands of N2 are found to be 132, 114, 19, 22, and 22 R, respectively. Overhead intensities are calculated for each vibrational transition of all the triplet band emissions of N2, which span a wider spectrum of wavelengths from ultraviolet to infrared. The calculated limb intensities of total and prominent transitions of VK band are presented. The model limb intensity of VK emission within the 150-190 nm wavelength region is in good agreement with the Cassini UVIS observed limb profile. An assessment of the impact of solar EUV flux on the N2 triplet band emission intensity has been made by using three different solar flux models, viz., Solar EUV Experiment, SOLAR2000 model of Tobiska (2004), and HEUVAC model of Richards et al, (2006). The calculated N2 VK band intensity at the peak of limb intensity due to S2K and HEUVAC solar flux models is a factor of 1.2 and 0.9, respectively, of that obtained using SEE solar EUV flux.
△ Less
Submitted 10 February, 2012;
originally announced February 2012.
-
Impact of solar EUV flux on CO Cameron band and CO2+ UV doublet emissions in the dayglow of Mars
Authors:
Sonal Kumar Jain,
Anil Bhardwaj
Abstract:
This study is aimed at making a calculation about the impact of the two most commonly used solar EUV flux models -- SOLAR2000 (S2K) of \cite{Tobiska04} and EUVAC model of \cite{Richards94} -- on photoelectron fluxes, volume emission rates, ion densities and CO Cameron and CO$_2^+$ UV doublet band dayglow emissions on Mars in three solar activity conditions: minimum, moderate, and maximum. Calculat…
▽ More
This study is aimed at making a calculation about the impact of the two most commonly used solar EUV flux models -- SOLAR2000 (S2K) of \cite{Tobiska04} and EUVAC model of \cite{Richards94} -- on photoelectron fluxes, volume emission rates, ion densities and CO Cameron and CO$_2^+$ UV doublet band dayglow emissions on Mars in three solar activity conditions: minimum, moderate, and maximum. Calculated limb intensities profiles are compared with SPICAM/Mars Express and Mariner observations. Analytical yield spectrum (AYS) approach has been used to calculate photoelectron fluxes in Martian upper atmosphere. Densities of prominent ions and CO molecule in excited triplet a$^3Π$ state are calculated using major ion-neutral reactions. Volume emission rates of CO Cameron and CO$_2^+$ UV doublet bands have been calculated for dif{}ferent observations (Viking condition, Mariner and Mars Express SPICAM observations) on Mars. For the low solar activity condition, dayglow intensities calculated using the S2K model are $\sim$40% higher than those calculated using the EUVAC model. During high solar activity, due to the higher EUV fluxes at wavelengths below 250 Å in the EUVAC model, intensities calculated using EUVAC model are slightly higher ($\sim$20%) than those calculated using S2K model. Irrespective of the solar activity condition, production of Cameron band due to photodissociative excitation of CO$_2$ is around 50% higher when S2K model is used. Altitude of peak limb brightness of CO Cameron and CO$_2^+$ UV doublet band is found to be independent of solar EUV flux models. Calculated limb intensities of CO Cameron and CO$_2^+$ UV doublet bands are on an average a factor of $\sim$2 and $\sim$1.5, respectively, higher than the SPICAM Mars Express observation, while they are consistent with the Mariner observations.
△ Less
Submitted 18 October, 2011;
originally announced October 2011.
-
Calculations of N2 triplet states vibrational populations and band emissions in Venusian dayglow
Authors:
Anil Bhardwaj,
Sonal Kumar Jain
Abstract:
A model for N2 triplet states band emissions in the Venusian dayglow has been developed for low and high solar activity conditions. Steady state photoelectron fluxes and volume excitation rates for N2 triplet states have been calculated using the Analytical Yield Spectra (AYS) technique. Model calculated photoelectron flux is in good agreement with Pioneer Venus Orbiter-observed electron flux. Sin…
▽ More
A model for N2 triplet states band emissions in the Venusian dayglow has been developed for low and high solar activity conditions. Steady state photoelectron fluxes and volume excitation rates for N2 triplet states have been calculated using the Analytical Yield Spectra (AYS) technique. Model calculated photoelectron flux is in good agreement with Pioneer Venus Orbiter-observed electron flux. Since inter-state cascading is important for the triplet states of N2, populations of different levels of N2 triplet states are calculated under statistical equilibrium considering direct electron impact excitation, and cascading and quenching effects. Densities of all vibrational levels of each triplet state are calculated in the model. Height-integrated overhead intensities of N2 triplet band emissions are calculated, the values for Vegard-Kaplan (A^3Sigma_u^+ - X^1Pi_g^+), First Positive (B^3Pi_g - A^3Sigma_u^+), Second Positive (C^3Pi_u - B^3Pi_g), and Wu-Benesch (W^3Delta_u - B^3Pi_g) bands of N2, are 1.9 (3.2), 3 (6), 0.4 (0.8), and 0.5 (1.1) kR, respectively, for solar minimum (maximum) conditions. The intensities of the three strong Vegard-Kaplan bands (0, 5), (0, 6), and (0, 7) are 94 (160), 120 (204), and 114 (194) R, respectively, for solar minimum (maximum) conditions. Limb profiles are calculated for VK (0, 4), (0, 5), (0, 6) and (0, 7) bands. The calculated intensities on Venus are about a factor 10 higher than those on Mars. The present study provides a motivation for a search of N2 triplet band emissions in the dayglow of Venus.
△ Less
Submitted 3 August, 2011;
originally announced August 2011.
-
Model calculation of N2 Vegard-Kaplan band emissions in Martian dayglow
Authors:
Sonal Kumar Jain,
Anil Bhardwaj
Abstract:
A model for N2 Vegard-Kaplan (VK) band (A^3Sigma_u^+ - X^1Sigma_g^+) emissions in Martian dayglow has been developed to explain the recent observations made by the Spectroscopy for Investigation of Characteristics of the Atmosphere of Mars (SPICAM) ultraviolet spectrograph aboard Mars Express. Steady state photoelectron fluxes and volume excitation rates have been calculated using the analytical y…
▽ More
A model for N2 Vegard-Kaplan (VK) band (A^3Sigma_u^+ - X^1Sigma_g^+) emissions in Martian dayglow has been developed to explain the recent observations made by the Spectroscopy for Investigation of Characteristics of the Atmosphere of Mars (SPICAM) ultraviolet spectrograph aboard Mars Express. Steady state photoelectron fluxes and volume excitation rates have been calculated using the analytical yield spectra technique. Since interstate cascading is important for triplet states of N2, the population of any given level of N2 triplet states is calculated under statistical equilibrium considering direct excitation, cascading, and quenching effects. Relative population of all vibrational levels of each triplet state is calculated in the model. Line of sight intensities and height-integrated overhead intensities have been calculated for VK, first positive (B^3Pi_g - A^3Sigma_u^+), second positive (C^3Pi_u - B^3Pi_g), and Wu-Benesch (W^3Delta_u - B^3Pi_g) bands of N2. A reduction in the N2 density by a factor of 3 in the Mars thermospheric general circulation model is required to obtain agreement between calculated limb profiles of VK (0-6) and SPICAM observation. Calculations are carried out to asses the impact of model parameters, namely, electron impact cross sections, solar EUV flux, and model atmosphere, on the emission intensities. Constraining the N2/CO2 ratio by SPICAM observations, we suggest the N2/CO2 ratios to be in the range 1.1-1.4% at 120 km, 1.8-3.2% at 140 km, and 4-7% at 170 km. During high solar activity the overhead intensity of N2 VK band emissions would be ~2.5 times higher than that during low solar activity.
△ Less
Submitted 2 August, 2011;
originally announced August 2011.
-
Monte Carlo model of electron energy degradation in a CO2 atmosphere
Authors:
Anil Bhardwaj,
Sonal Kumar Jain
Abstract:
A Monte Carlo model has been developed to study the degradation of <1000 eV electrons in an atmosphere of CO2, which is one of the most abundant species in Mars' and Venus' atmospheres. The e-CO2 cross sections are presented in an assembled set along with their analytical representations. Monte Carlo simulations are carried out at several energies to calculate the "yield spectra", which embodied a…
▽ More
A Monte Carlo model has been developed to study the degradation of <1000 eV electrons in an atmosphere of CO2, which is one of the most abundant species in Mars' and Venus' atmospheres. The e-CO2 cross sections are presented in an assembled set along with their analytical representations. Monte Carlo simulations are carried out at several energies to calculate the "yield spectra", which embodied all the information related to electron degradation process and can be used to calculate "yield" (or population) for any inelastic process. The numerical yield spectra have been fitted analytically resulting in an analytical yield spectra (AYS). We have calculated the mean energy per ion pair and efficiencies for various inelastic processes, including the double and dissociative double ionization of \car\ and negative ion formation. The energy distribution of the secondary electrons produced per incident electron is also presented at few incident energies. The mean energy per ion pair for CO2 is 37.5 (35.8) eV at 200 (1000) eV, compared to experimental value 32.7 eV at high energies. Ionization is the dominant loss process at energies above 50 eV with contribution of ~50%. Among the excitation processes, 13.6 eV and 12.4 eV states are the dominant loss processes consuming ~28% energy above 200 eV. Around and below ionization threshold, 13.6 eV, 12.4 eV, and 11.1 eV, followed by 8.6 eV and 9.3 eV excitation states are important loss processes, while below 10 eV vibrational excitation dominates.
△ Less
Submitted 19 November, 2010;
originally announced November 2010.
