Showing 1–50 of 140 results for author: Kashyap, V

International Astrophysical Consortium for High-energy Calibration: Summary of the 15th IACHEC Workshop

Authors: K. K. Madsen, V. Burwitz, K. Forster, C. E. Grant, M. Guainazzi, V. Kashyap, H. L. Marshall, E. D. Miller, L. Natalucci, P. P. Plucinsky, Y. Terada

Abstract: In this report, we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) from the 15th IACHEC Workshop in Pelham, Germany. Sixty scientists directly involved in the calibration of operational and future high-energy missions gathered for 3.5 days to discuss the status of the cross-calibration between the current international complement of X-ray… ▽ More In this report, we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) from the 15th IACHEC Workshop in Pelham, Germany. Sixty scientists directly involved in the calibration of operational and future high-energy missions gathered for 3.5 days to discuss the status of the cross-calibration between the current international complement of X-ray observatories and the possibilities to improve it. This summary consists of reports from the Working Groups with topics ranging across the identification and characterization of standard calibration sources, multi-observatory cross-calibration campaigns, appropriate and new statistical techniques, calibration of instruments and characterization of background, preservation of knowledge, and results for the benefit of the astronomical community. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 10 pages, 1 figure. arXiv admin note: text overlap with arXiv:2111.01613

Light Dark Matter Constraints from SuperCDMS HVeV Detectors Operated Underground with an Anticoincidence Event Selection

Authors: SuperCDMS Collaboration, M. F. Albakry, I. Alkhatib, D. Alonso-González, D. W. P. Amaral, J. Anczarski, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, C. Bathurst, R. Bhattacharyya, A. J. Biffl, P. L. Brink, M. Buchanan, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, J. -H. Chen , et al. (115 additional authors not shown)

Abstract: This article presents constraints on dark-matter-electron interactions obtained from the first underground data-taking campaign with multiple SuperCDMS HVeV detectors operated in the same housing. An exposure of 7.63 g-days is used to set upper limits on the dark-matter-electron scattering cross section for dark matter masses between 0.5 and 1000 MeV/$c^2$, as well as upper limits on dark photon k… ▽ More This article presents constraints on dark-matter-electron interactions obtained from the first underground data-taking campaign with multiple SuperCDMS HVeV detectors operated in the same housing. An exposure of 7.63 g-days is used to set upper limits on the dark-matter-electron scattering cross section for dark matter masses between 0.5 and 1000 MeV/$c^2$, as well as upper limits on dark photon kinetic mixing and axion-like particle axioelectric coupling for masses between 1.2 and 23.3 eV/$c^2$. Compared to an earlier HVeV search, sensitivity was improved as a result of an increased overburden of 225 meters of water equivalent, an anticoincidence event selection, and better pile-up rejection. In the case of dark-matter-electron scattering via a heavy mediator, an improvement by up to a factor of 25 in cross-section sensitivity was achieved. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 7 pages + title and references, 4 figures, and 1 table

Separating States in Astronomical Sources Using Hidden Markov Models: With a Case Study of Flaring and Quiescence on EV Lac

Authors: Robert Zimmerman, David A. van Dyk, Vinay L. Kashyap, Aneta Siemiginowska

Abstract: We present a new method to distinguish between different states (e.g., high and low, quiescent and flaring) in astronomical sources with count data. The method models the underlying physical process as latent variables following a continuous-space Markov chain that determines the expected Poisson counts in observed light curves in multiple passbands. For the underlying state process, we consider s… ▽ More We present a new method to distinguish between different states (e.g., high and low, quiescent and flaring) in astronomical sources with count data. The method models the underlying physical process as latent variables following a continuous-space Markov chain that determines the expected Poisson counts in observed light curves in multiple passbands. For the underlying state process, we consider several autoregressive processes, yielding continuous-space hidden Markov models of varying complexity. Under these models, we can infer the state that the object is in at any given time. The state predictions from these models are then dichotomized with the help of a finite-mixture model to produce state classifications. We apply these techniques to X-ray data from the active dMe flare star EV Lac, splitting the data into quiescent and flaring states. We find that a first-order vector autoregressive process efficiently separates flaring from quiescence: flaring occurs over 30-40% of the observation durations, a well-defined persistent quiescent state can be identified, and the flaring state is characterized by higher temperatures and emission measures. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Comments: 25 pages, 12 figures, 12 tables

Effect of Systematic Uncertainties on Density and Temperature Estimates in Coronae of Capella

Authors: Xixi Yu, Vinay L. Kashyap, Giulio Del Zanna, David A. van Dyk, David C. Stenning, Connor P. Ballance, Harry P. Warren

Abstract: We estimate the coronal density of Capella using the O VII and Fe XVII line systems in the soft X-ray regime that have been observed over the course of the Chandra mission. Our analysis combines measures of error due to uncertainty in the underlying atomic data with statistical errors in the Chandra data to derive meaningful overall uncertainties on the plasma density of the coronae of Capella. We… ▽ More We estimate the coronal density of Capella using the O VII and Fe XVII line systems in the soft X-ray regime that have been observed over the course of the Chandra mission. Our analysis combines measures of error due to uncertainty in the underlying atomic data with statistical errors in the Chandra data to derive meaningful overall uncertainties on the plasma density of the coronae of Capella. We consider two Bayesian frameworks. First, the so-called pragmatic-Bayesian approach considers the atomic data and their uncertainties as fully specified and uncorrectable. The fully-Bayesian approach, on the other hand, allows the observed spectral data to update the atomic data and their uncertainties, thereby reducing the overall errors on the inferred parameters. To incorporate atomic data uncertainties, we obtain a set of atomic data replicates, the distribution of which captures their uncertainty. A principal component analysis of these replicates allows us to represent the atomic uncertainty with a lower-dimensional multivariate Gaussian distribution. A $t$-distribution approximation of the uncertainties of a subset of plasma parameters including a priori temperature information, obtained from the temperature-sensitive-only Fe XVII spectral line analysis, is carried forward into the density- and temperature-sensitive O VII spectral line analysis. Markov Chain Monte Carlo based model fitting is implemented including Multi-step Monte Carlo Gibbs Sampler and Hamiltonian Monte Carlo. Our analysis recovers an isothermally approximated coronal plasma temperature of $\approx$5 MK and a coronal plasma density of $\approx$10$^{10}$ cm$^{-3}$, with uncertainties of 0.1 and 0.2 dex respectively. △ Less

Submitted 18 June, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Joint Deconvolution of Astronomical Images in the Presence of Poisson Noise

Authors: Axel Donath, Aneta Siemiginowska, Vinay L. Kashyap, David A. van Dyk, Douglas Burke

Abstract: We present a new method for joint likelihood deconvolution (Jolideco) of a set of astronomical observations of the same sky region in the presence of Poisson noise. The observations may be obtained from different instruments with different resolution, and different point spread functions. Jolideco reconstructs a single flux image by optimizing the posterior distribution based on the joint Poisson… ▽ More We present a new method for joint likelihood deconvolution (Jolideco) of a set of astronomical observations of the same sky region in the presence of Poisson noise. The observations may be obtained from different instruments with different resolution, and different point spread functions. Jolideco reconstructs a single flux image by optimizing the posterior distribution based on the joint Poisson likelihood of all observations under a patch-based image prior. The patch prior is parameterised via a Gaussian Mixture model which we train on high-signal-to-noise astronomical images, including data from the James Webb Telescope and the GLEAM radio survey. This prior favors correlation structures among the reconstructed pixel intensities that are characteristic of those observed in the training images. It is, however, not informative for the mean or scale of the reconstruction. By applying the method to simulated data we show that the combination of multiple observations and the patch-based prior leads to much improved reconstruction quality in many different source scenarios and signal to noise regimes. We demonstrate that with the patch prior Jolideco yields superior reconstruction quality relative to alternative standard methods such as the Richardson-Lucy method. We illustrate the results of Jolideco applied to example data from the Chandra X-ray Observatory and the Fermi-LAT Gamma-ray Space Telescope. By comparing the measured width of a counts based and the corresponding Jolideco flux profile of an X-ray filament in SNR 1E 0102.2-721} we find the deconvolved width of 0.58+- 0.02 arcsec to be consistent with the theoretical expectation derived from the known width of the PSF. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 25 pages, 11 figures

Spectral fit residuals as an indicator to increase model complexity

Authors: Anshuman Acharya, Vinay L. Kashyap

Abstract: Spectral fitting of X-ray data usually involves minimizing statistics like the chi-square and the Cash statistic. Here we discuss their limitations and introduce two measures based on the cumulative sum (CuSum) of model residuals to evaluate whether model complexity could be increased: the percentage of bins exceeding a nominal threshold in a CuSum array (pct$_{CuSum}$), and the excess area under… ▽ More Spectral fitting of X-ray data usually involves minimizing statistics like the chi-square and the Cash statistic. Here we discuss their limitations and introduce two measures based on the cumulative sum (CuSum) of model residuals to evaluate whether model complexity could be increased: the percentage of bins exceeding a nominal threshold in a CuSum array (pct$_{CuSum}$), and the excess area under the CuSum compared to the nominal (p$_\textit{area}$). We demonstrate their use with an application to a $\textit{Chandra}$ ACIS spectral fit. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Comments: 3 pages, 1 figure, published in the Research Notes of the American Astronomical Society (RNAAS)

Journal ref: Res. Notes AAS 8 1 (2024)

Front-row seat of the recent R Aqr periastron passage: X-ray multi-epoch spectral and spatial analysis

Authors: A. Sacchi, M. Karovska, J. Raymond, V. Kashyap, T. J. Gaetz, W. Hack, J. Kennea, N. Lee, A. J Mioduszewski, M. J Claussen

Abstract: We report on the X-ray spectral and spatial evolution of the Symbiotic star R Aqr. Through a multi-epoch observational campaign performed with Chandra between 2017 and 2022, we study the X-ray emission of this binary system, composed of an evolved red giant star and a white dwarf (WD). This analysis is particularly timely as the WD approached the periastron in late 2018/early 2019, thus mass trans… ▽ More We report on the X-ray spectral and spatial evolution of the Symbiotic star R Aqr. Through a multi-epoch observational campaign performed with Chandra between 2017 and 2022, we study the X-ray emission of this binary system, composed of an evolved red giant star and a white dwarf (WD). This analysis is particularly timely as the WD approached the periastron in late 2018/early 2019, thus mass transfer, jet emission and outburst phenomena are to be expected. Through detailed spectral analysis, we detect a significant rise in the soft X-ray (0.5-2 keV) emission of R Aqr, likely linked to jet emission, followed by a decay towards the previous quiescent state. The hard X-ray emission (5-8 keV), is not immediately affected by the periastron passage; the hard component, after maintaining the same flux level between 2017 and 2021, rapidly decays after 2022. Possible explanations for this are a change in the reflection properties of the medium surrounding the binary, obscuration of the central region by material ejected during the periastron passage, or even the partial/complete destruction of the inner regions of the accretion disc surrounding the WD. In addition to this activity in the central region, extended emission is also detected, likely linked to a hot spot in a pre-outburst-emitted jet, which can be observed moving away from the system's central region. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 14 pages, 9 figures, 3 tables. Accepted for publication in ApJ

Breakthroughs in Cool Star Physics with the Line Emission Mapper X-ray Probe

Authors: Jeremy J. Drake, Julián Alvarado Gomez, Costanza Argiroffi, Ettore Flaccomio, Cecilia Garraffo, Nicolas Grosso, Nazma Islam, Margarita Karovska, Vinay L. Kashyap, Kristina Monsch, Jan-Uwe Ness, Salvatore Sciortino, Bradford Wargelin

Abstract: We outline some of the highlights of the scientific case for the advancement of stellar high energy physics using the Line Emission Mapper X-ray Probe ({\it LEM}). The key to advancements with LEM lie in its large effective area -- up to 100 times that of the {\it Chandra} MEG -- and 1~eV spectral resolution. The large effective area opens up for the first time the ability to study time-dependent… ▽ More We outline some of the highlights of the scientific case for the advancement of stellar high energy physics using the Line Emission Mapper X-ray Probe ({\it LEM}). The key to advancements with LEM lie in its large effective area -- up to 100 times that of the {\it Chandra} MEG -- and 1~eV spectral resolution. The large effective area opens up for the first time the ability to study time-dependent phenomena on their natural timescales at high resolution, such as flares and coronal mass ejections, and also opens the sky to much fainter targets than available to {\it Chandra} or {\it XMM-Newton}. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: A Line Emission Mapper X-ray Probe White Paper

A Possible Third Body in the X-Ray System GRS 1747-312 and Models with Higher-Order Multiplicity

Authors: Caleb Painter, Rosanne Di Stefano, Vinay L. Kashyap, Roberto Soria, Jose Lopez-Miralles, Ryan Urquhart, James F. Steiner, Sara Motta, Darin Ragozzine, Hideyuki Mori

Abstract: GRS 1747-312 is a bright Low-Mass X-ray Binary in the globular cluster Terzan 6, located at a distance of 9.5 kpc from the Earth. It exhibits regular outbursts approximately every 4.5 months, during which periodic eclipses are known to occur. These eclipses have only been observed in the outburst phase, and are not clearly seen when the source is quiescent. Recent Chandra observations of the sourc… ▽ More GRS 1747-312 is a bright Low-Mass X-ray Binary in the globular cluster Terzan 6, located at a distance of 9.5 kpc from the Earth. It exhibits regular outbursts approximately every 4.5 months, during which periodic eclipses are known to occur. These eclipses have only been observed in the outburst phase, and are not clearly seen when the source is quiescent. Recent Chandra observations of the source were performed in June 2019 and April, June, and August of 2021. Two of these observations captured the source during its outburst, and showed clear flux decreases at the expected time of eclipse. The other two observations occurred when the source was quiescent. We present the discovery of a dip that occurred during the quiescent state. The dip is of longer duration and its time of occurrence does not fit the ephemeris of the shorter eclipses. We study the physical characteristics of the dip and determine that it has all the properties of an eclipse by an object with a well defined surface. We find that there are several possibilities for the nature of the object causing the 5.3 ks eclipse. First, GRS 1747-312 may be an X-ray triple, with an LMXB orbited by an outer third object, which could be an M-dwarf, brown dwarf, or planet. Second, there could be two LMXBs in close proximity to each other, likely bound together. Whatever the true nature of the eclipser, its presence suggests that the GRS 1747-312 system is unique. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 38 pages, 30 figures, submitted to MNRAS

Estimating Ejecta Mass Ratios in Kepler's SNR: Global X-Ray Spectral Analysis Including Suzaku Systematics and Emitting Volume Uncertainties

Authors: Tyler Holland-Ashford, Patrick Slane, Laura A. Lopez, Katie Auchettl, Vinay Kashyap

Abstract: The exact origins of many Type Ia supernovae$\unicode{x2013}$progenitor scenarios and explosive mechanisms$\unicode{x2013}$remain uncertain. In this work, we analyze the global Suzaku X-Ray spectrum of Kepler's supernova remnant in order to constrain mass ratios of various ejecta species synthesized during explosion. Critically, we account for the Suzaku telescope effective area calibration uncert… ▽ More The exact origins of many Type Ia supernovae$\unicode{x2013}$progenitor scenarios and explosive mechanisms$\unicode{x2013}$remain uncertain. In this work, we analyze the global Suzaku X-Ray spectrum of Kepler's supernova remnant in order to constrain mass ratios of various ejecta species synthesized during explosion. Critically, we account for the Suzaku telescope effective area calibration uncertainties of 5$\unicode{x2013}$20% by generating 100 mock effective area curves and using Markov Chain Monte Carlo based spectral fitting to produce 100 sets of best-fit parameter values. Additionally, we characterize the uncertainties from assumptions made about the emitting volumes of each model plasma component: finding that these uncertainties can be the dominant source of error. We then compare our calculated mass ratios to previous observational studies of Kepler's SNR and to the predictions of SN Ia simulations. Our mass ratio estimates require a $\sim$90% attenuated $^{12}$C$+^{16}$O reaction rate and are potentially consistent with both near- and sub-M$_{\rm Ch}$ progenitors, but are inconsistent with the dynamically stable double detonation origin scenario and only marginally consistent with the dynamically unstable dynamically-driven double-degenerate double detonation (D$^6$) scenario. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 25 pages, 12 figures, 4 tables

Classification of Chandra X-ray Sources in Cygnus OB2

Authors: Vinay L. Kashyap, Mario G. Guarcello, Nicholas J. Wright, Jeremy J. Drake, Ettore Flaccomio, Tom L. Aldcroft, Juan F. Albacete Colombo, Kevin Briggs, Francesco Damiani, Janet E. Drew, Eduardo L. Martin, Giusi Micela, Tim Naylor, Salvatore Sciortino

Abstract: We have devised a predominantly Naive Bayes method to classify the optical/IR matches to X-ray sources detected by Chandra in the Cygnus OB2 association into foreground, member, and background objects. We employ a variety of X-ray, optical, and infrared characteristics to construct likelihoods using training sets defined by well-measured sources. Combinations of optical photometry from SDSS (riz)… ▽ More We have devised a predominantly Naive Bayes method to classify the optical/IR matches to X-ray sources detected by Chandra in the Cygnus OB2 association into foreground, member, and background objects. We employ a variety of X-ray, optical, and infrared characteristics to construct likelihoods using training sets defined by well-measured sources. Combinations of optical photometry from SDSS (riz) and IPHAS (riHa), IR magnitudes from UKIDSS and 2MASS (JHK), X-ray quantiles and hardness ratios, and estimates of extinction Av are used to compute the relative probabilities that a given source belongs to one of the classes. We use Principal Component Analysis of photometric magnitude combinations to isolate the best axes for classification. We incorporate measurement errors into the classification. We evaluate the accuracy of the classification by inspection and reclassify a number of sources based on IR magnitudes, presence of disks, and X-ray spectral hardness. We also consider systematic errors due to extinction. We find that about 6100 objects are association members, 1400 are background, and 500 are foreground objects. The overall classification accuracy is 95%. △ Less

Submitted 11 June, 2023; originally announced June 2023.

Comments: 27 pages, 23 figures, 6 tables; accepted for publication in ApJS. Full Table 3 is in Zenodo at https://zenodo.org/record/8025756

MSC Class: 85A35 (Primary); 85A15; 62H30; 62P99 (Secondary)

The extremely X-ray luminous radio-loud quasar CFHQS J142952+544717 at $z=6.18$ under Chandra high-angular resolution lens

Authors: G. Migliori, A. Siemiginowska, M. Sobolewska, C. C. Cheung, Ł. Stawarz, D. Schwartz, B. Snios, A. Saxena, V. Kashyap

Abstract: We present the first X-ray observation at sub-arcsecond resolution of the high-redshift ($z=6.18$) radio-loud quasar CFHQS J142952+544717 (J1429). The ~100 net-count 0.3-7 keV spectrum obtained from $\sim 30$ ksec Chandra exposure is best fit by a single power-law model with a photon index $Γ=2.0\pm0.2$ and no indication of an intrinsic absorber, implying a 3.6-72 keV rest-frame luminosity… ▽ More We present the first X-ray observation at sub-arcsecond resolution of the high-redshift ($z=6.18$) radio-loud quasar CFHQS J142952+544717 (J1429). The ~100 net-count 0.3-7 keV spectrum obtained from $\sim 30$ ksec Chandra exposure is best fit by a single power-law model with a photon index $Γ=2.0\pm0.2$ and no indication of an intrinsic absorber, implying a 3.6-72 keV rest-frame luminosity $L_{\rm X}=(2.3^{+0.6}_{-0.5})\times10^{46}$ erg s$^{-1}$. We identify a second X-ray source at 30 arcsec, distance from J1429 position, with a soft ($Γ\simeq 2.8$) and absorbed (equivalent hydrogen column density $N_{\rm H} <13.4\times 10^{20}$ cm$^{-2}$) spectrum, which likely contaminated J1429 spectra obtained in lower angular resolution observations. Based on the analysis of the Chandra image, the bulk of the X-ray luminosity is produced within the central $\sim 3$ kpc region, either by the disk/corona system, or by a moderately aligned jet. In this context, we discuss the source properties in comparison with samples of low- and high-redshift quasars. We find indication of a possible excess of counts over the expectations for a point-like source in a 0.5 arcsec-1.5 arcsec ($\sim 3-8$ kpc) annular region. The corresponding X-ray luminosity at J1429 redshift is $4\times 10^{45}$ erg s$^{-1}$. If confirmed, this emission could be related to either a large-scale X-ray jet, or a separate X-ray source. △ Less

Submitted 30 June, 2023; v1 submitted 12 May, 2023; originally announced May 2023.

Comments: 9 pages, 4 figures, accepted for publication in MNRAS

First measurement of the nuclear-recoil ionization yield in silicon at 100 eV

Authors: M. F. Albakry, I. Alkhatib, D. Alonso, D. W. P. Amaral, P. An, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, P. S. Barbeau, C. Bathurst, R. Bhattacharyya, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen, N. Chott , et al. (115 additional authors not shown)

Abstract: We measured the nuclear--recoil ionization yield in silicon with a cryogenic phonon-sensitive gram-scale detector. Neutrons from a mono-energetic beam scatter off of the silicon nuclei at angles corresponding to energy depositions from 4\,keV down to 100\,eV, the lowest energy probed so far. The results show no sign of an ionization production threshold above 100\,eV. These results call for furthe… ▽ More We measured the nuclear--recoil ionization yield in silicon with a cryogenic phonon-sensitive gram-scale detector. Neutrons from a mono-energetic beam scatter off of the silicon nuclei at angles corresponding to energy depositions from 4\,keV down to 100\,eV, the lowest energy probed so far. The results show no sign of an ionization production threshold above 100\,eV. These results call for further investigation of the ionization yield theory and a comprehensive determination of the detector response function at energies below the keV scale. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Journal ref: Physical Review Letters 131.9 (2023): 091801

Shocks and Photoionization of the Inner 650 AU Jet of the Interacting Binary Star R Aquarii from Multiwavelength Hubble Space Telescope Observations

Authors: Caroline D. Huang, Margarita Karovska, Warren Hack, John C. Raymond, Rodolfo Montez Jr., Vinay L. Kashyap

Abstract: Astrophysical jets are present in a range of environments, including young stellar objects, X-ray binaries, and active galactic nuclei, but their formation is still not fully understood. As one of the nearest symbiotic binary stars, R Aquarii ($D \sim 220$ pc) offers a unique opportunity to study the inner region within $\sim$ 600 AU of the jet source, which is particularly crucial to our understa… ▽ More Astrophysical jets are present in a range of environments, including young stellar objects, X-ray binaries, and active galactic nuclei, but their formation is still not fully understood. As one of the nearest symbiotic binary stars, R Aquarii ($D \sim 220$ pc) offers a unique opportunity to study the inner region within $\sim$ 600 AU of the jet source, which is particularly crucial to our understanding of non-relativistic jet formation and origin. We present high-angular resolution ultraviolet and optical imaging from the \emph{Hubble} Space Telescope in six emission-line regions of the inner jet. Using these observations to obtain a range of representative line ratios for our system and kinematic data derived from a comparison with previous studies, we model the shocked gas in order to determine the relative roles of shock heating and photoionization in the R Aquarii system. We find that our shock models suggest a nonzero magnetic field is needed to describe the measured line ratios. We also find that the Mg~II$λλ$2795,2802 intensities are overpredicted by our models for most of the jet regions, perhaps because of depletion onto grains or to opacity in these resonance lines. △ Less

Submitted 1 March, 2023; originally announced March 2023.

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

A novel approach to detect line emission under high background in high-resolution X-ray spectra

Authors: Xiangyu Zhang, Sara Algeri, Vinay Kashyap, Margarita Karovska

Abstract: We develop a novel statistical approach to identify emission features or set upper limits in high-resolution spectra in the presence of high background. The method relies on detecting differences from the background using smooth tests and using classical likelihood ratio tests to characterise known shapes like emission lines. We perform signal detection or place upper limits on line fluxes while a… ▽ More We develop a novel statistical approach to identify emission features or set upper limits in high-resolution spectra in the presence of high background. The method relies on detecting differences from the background using smooth tests and using classical likelihood ratio tests to characterise known shapes like emission lines. We perform signal detection or place upper limits on line fluxes while accounting for the problem of multiple comparisons. We illustrate the method by applying it to a Chandra LETGS+HRC-S observation of symbiotic star RT Cru, successfully detecting previously known features like the Fe line emission in the 6-7 keV range and the Iridium-edge due to the mirror coating on Chandra. We search for thermal emission lines from Ne X, Fe XVII, O VIII, and O VII, but do not detect them, and place upper limits on their intensities consistent with a $\approx$1 keV plasma. We serendipitously detect a line at 16.93 $\unicode{x212B}$ that we attribute to photoionisation or a reflection component. △ Less

Submitted 9 March, 2023; v1 submitted 1 February, 2023; originally announced February 2023.

Comments: Accepted by Monthly Notices of the Royal Astronomical Society

Center to limb variation of transition region Doppler shift in active regions

Authors: Abhishek Rajhans, Durgesh Tripathi, Vinay L. Kashyap, James A. Klimchuk, Avyarthana Ghosh

Abstract: Studying Doppler shifts provides deeper insights into the flow of mass and energy in the solar atmosphere. We perform a comprehensive measurement of Doppler shifts in the transition region and its center-to-limb variation (CLV) in the strong field regions ($|\textbf{B}| \geq$ 50 G) of 50 active regions (ARs), using the \ion{Si}{4} 1394~Å line recorded by the Interface Region Imaging Spectrometer(I… ▽ More Studying Doppler shifts provides deeper insights into the flow of mass and energy in the solar atmosphere. We perform a comprehensive measurement of Doppler shifts in the transition region and its center-to-limb variation (CLV) in the strong field regions ($|\textbf{B}| \geq$ 50 G) of 50 active regions (ARs), using the \ion{Si}{4} 1394~Å line recorded by the Interface Region Imaging Spectrometer(IRIS). To locate the ARs and identify strong field regions, we have used the magnetograms obtained by the Helioseismic and Magnetic Imager (HMI). We find that in strong field regions, on average, all the ARs show mean redshifts ranging between 4{--}11~ km/s, which varies with ARs. These flows show a mild CLV, with sizable magnitudes at the limb and substantial scatter at the mid-longitude range. Our observations do not support the idea that redshifts in the lower transition region (T $<\sim$ 0.1 MK) are produced by field-aligned downflows as a result of impulsive heating and warrant alternative interpretation, such as downflow of type-\rm{II} spicules in the presence of a chromospheric wall created by cooler type-\rm{I} spicules. △ Less

Submitted 18 February, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: Accepted for publication in the Astrophysical Journal on 16th January 2023

X-ray Activity Variations and Coronal Abundances of the Star-Planet Interaction candidate HD 179949

Authors: Anshuman Acharya, Vinay L. Kashyap, Steven H. Saar, Kulinder Pal Singh, Manfred Cuntz

Abstract: We carry out detailed spectral and timing analyses of the $Chandra$ X-ray data of HD 179949, a prototypical example of a star with a close-in giant planet with possible star-planet interaction (SPI) effects. We find a low coronal abundance $A({\rm Fe})/A({\rm H}){\approx}0.2$ relative to the solar photospheric baseline of Anders & Grevesse (1989), and significantly lower than the stellar photosphe… ▽ More We carry out detailed spectral and timing analyses of the $Chandra$ X-ray data of HD 179949, a prototypical example of a star with a close-in giant planet with possible star-planet interaction (SPI) effects. We find a low coronal abundance $A({\rm Fe})/A({\rm H}){\approx}0.2$ relative to the solar photospheric baseline of Anders & Grevesse (1989), and significantly lower than the stellar photosphere as well. We further find low abundances of high First Ionization Potential (FIP) elements $A({\rm O})/A({\rm Fe}){\lesssim}1$, $A({\rm Ne})/A({\rm Fe}){\lesssim}0.1$, but with indications of higher abundances of $A({\rm N})/A({\rm Fe}){\gg}1, A({\rm Al})/A({\rm Fe}){\lesssim}10$. We estimate a FIP bias for this star in the range $\approx{-0.3}$ to ${-0.1}$, larger than the ${\lesssim}-$0.5 expected for stars of this type, but similar to stars hosting close-in hot Jupiters. We detect significant intensity variability over time scales ranging from 100 s - 10 ks, and also evidence for spectral variability over time scales of 1-10 ks. We combine the $Chandra$ flux measurements with $Swift$ and XMM-$Newton$ measurements to detect periodicities and determine that the dominant signal is tied to the stellar polar rotational period, consistent with expectations that the corona is rotational-pole dominated. We also find evidence for periodicity at both the planetary orbital frequency and at its beat frequency with the stellar polar rotational period, suggesting the presence of a magnetic connection between the planet and the stellar pole. If these periodicities represent an SPI signal, it is likely driven by a quasi-continuous form of heating (e.g., magnetic field stretching) rather than sporadic, hot, impulsive flare-like reconnections. △ Less

Submitted 14 July, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 26 pages, 12 figures, 11 tables. Accepted for publication in the Astrophysical Journal

Identifying diffuse spatial structures in high-energy photon lists

Authors: Minjie Fan, Jue Wang, Vinay L. Kashyap, Thomas C. M. Lee, David A. van Dyk, Andreas Zezas

Abstract: Data from high-energy observations are usually obtained as lists of photon events. A common analysis task for such data is to identify whether diffuse emission exists, and to estimate its surface brightness, even in the presence of point sources that may be superposed. We have developed a novel non-parametric event list segmentation algorithm to divide up the field of view into distinct emission c… ▽ More Data from high-energy observations are usually obtained as lists of photon events. A common analysis task for such data is to identify whether diffuse emission exists, and to estimate its surface brightness, even in the presence of point sources that may be superposed. We have developed a novel non-parametric event list segmentation algorithm to divide up the field of view into distinct emission components. We use photon location data directly, without binning them into an image. We first construct a graph from the Voronoi tessellation of the observed photon locations and then grow segments using a new adaptation of seeded region growing, that we call Seeded Region Growing on Graph, after which the overall method is named SRGonG. Starting with a set of seed locations, this results in an over-segmented dataset, which SRGonG then coalesces using a greedy algorithm where adjacent segments are merged to minimize a model comparison statistic; we use the Bayesian Information Criterion. Using SRGonG we are able to identify point-like and diffuse extended sources in the data with equal facility. We validate SRGonG using simulations, demonstrating that it is capable of discerning irregularly shaped low surface-brightness emission structures as well as point-like sources with strengths comparable to that seen in typical X-ray data. We demonstrate SRGonG's use on the Chandra data of the Antennae galaxies, and show that it segments the complex structures appropriately. △ Less

Submitted 4 November, 2022; v1 submitted 15 August, 2022; originally announced August 2022.

Effective Field Theory Analysis of CDMSlite Run 2 Data

Authors: SuperCDMS Collaboration, M. F. Albakry, I. Alkhatib, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, C. Bathurst, D. A. Bauer, L. V. S. Bezerra, R. Bhattacharyya, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen, N. Chott , et al. (105 additional authors not shown)

Abstract: CDMSlite Run 2 was a search for weakly interacting massive particles (WIMPs) with a cryogenic 600 g Ge detector operated in a high-voltage mode to optimize sensitivity to WIMPs of relatively low mass from 2 - 20 GeV/$c^2$. In this article, we present an effective field theory (EFT) analysis of the CDMSlite Run 2 data using an extended energy range and a comprehensive treatment of the expected back… ▽ More CDMSlite Run 2 was a search for weakly interacting massive particles (WIMPs) with a cryogenic 600 g Ge detector operated in a high-voltage mode to optimize sensitivity to WIMPs of relatively low mass from 2 - 20 GeV/$c^2$. In this article, we present an effective field theory (EFT) analysis of the CDMSlite Run 2 data using an extended energy range and a comprehensive treatment of the expected background. A binned likelihood Bayesian analysis was performed on the recoil energy data, taking into account the parameters of the EFT interactions and optimizing the data selection with respect to the dominant background components. Energy regions within 5$σ$ of known activation peaks were removed from the analysis. The Bayesian evidences resulting from the different operator hypotheses show that the CDMSlite Run 2 data are consistent with the background-only models and do not allow for a signal interpretation assuming any additional EFT interaction. Consequently, upper limits on the WIMP mass and coupling-coefficient amplitudes and phases are presented for each EFT operator. These limits improve previous CDMSlite Run 2 bounds for WIMP masses above 5 GeV/$c^2$. △ Less

Submitted 23 May, 2022; originally announced May 2022.

Comments: 16 pages, 8 figures

Investigating the sources of low-energy events in a SuperCDMS-HVeV detector

Authors: SuperCDMS Collaboration, M. F. Albakry, I. Alkhatib, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, C. Bathurst, D. A. Bauer, R. Bhattacharyya, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen, N. Chott, J. Cooley , et al. (104 additional authors not shown)

Abstract: Recent experiments searching for sub-GeV/$c^2$ dark matter have observed event excesses close to their respective energy thresholds. Although specific to the individual technologies, the measured excess event rates have been consistently reported at or below event energies of a few-hundred eV, or with charges of a few electron-hole pairs. In the present work, we operated a 1-gram silicon SuperCDMS… ▽ More Recent experiments searching for sub-GeV/$c^2$ dark matter have observed event excesses close to their respective energy thresholds. Although specific to the individual technologies, the measured excess event rates have been consistently reported at or below event energies of a few-hundred eV, or with charges of a few electron-hole pairs. In the present work, we operated a 1-gram silicon SuperCDMS-HVeV detector at three voltages across the crystal (0 V, 60 V and 100 V). The 0 V data show an excess of events in the tens of eV region. Despite this event excess, we demonstrate the ability to set a competitive exclusion limit on the spin-independent dark matter--nucleon elastic scattering cross section for dark matter masses of $\mathcal{O}(100)$ MeV/$c^2$, enabled by operation of the detector at 0 V potential and achievement of a very low $\mathcal{O}(10)$ eV threshold for nuclear recoils. Comparing the data acquired at 0 V, 60 V and 100 V potentials across the crystal, we investigated possible sources of the unexpected events observed at low energy. The data indicate that the dominant contribution to the excess is consistent with a hypothesized luminescence from the printed circuit boards used in the detector holder. △ Less

Submitted 11 October, 2022; v1 submitted 17 April, 2022; originally announced April 2022.

Time domain methods for X-ray and gamma-ray astronomy

Authors: Eric D. Feigelson, Vinay L. Kashyap, Aneta Siemiginowska

Abstract: A variety of statistical methods for understanding variability in the time domain for low count rate X-ray and gamma-ray sources are explored. Variability can be detected using nonparametric (Anderson-Darling and overdispersion tests) and parametric (sequential likelihood-based tests) tools. Once detected, variability can be characterized by nonparametric (autocorrelation function, structure funct… ▽ More A variety of statistical methods for understanding variability in the time domain for low count rate X-ray and gamma-ray sources are explored. Variability can be detected using nonparametric (Anderson-Darling and overdispersion tests) and parametric (sequential likelihood-based tests) tools. Once detected, variability can be characterized by nonparametric (autocorrelation function, structure function,wavelet analysis) and parametric (multiple change point model such as Bayesian Blocks, integer autoregressive models, C-statistic and Poisson regression) methods. New multidimensional variability detection approaches are outlined. Software packages designed for high energy data analysis are deficient but tools are available in the R statistical software environment. Most of the methods presented here are not commonly used in high energy astronomy. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: 26 pages, 9 figures. To appear in Handbook for X-ray and Gamma-Ray Astrophysics, Volume 4: Analysis techniques, Section XVIII: Timing Analysis (Belloni & Bhattacharya, eds., Springer)

A Strategy for Low-Mass Dark Matter Searches with Cryogenic Detectors in the SuperCDMS SNOLAB Facility

Authors: SuperCDMS Collaboration, M. F. Albakry, I. Alkhatib, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, C. Bathurst, D. A. Bauer, R. Bhattacharyya, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeno, Y. -Y. Chang, M. Chaudhuri, R. Chen, N. Chott, J. Cooley , et al. (103 additional authors not shown)

Abstract: The SuperCDMS Collaboration is currently building SuperCDMS SNOLAB, a dark matter search focused on nucleon-coupled dark matter in the 1-5 GeV/c$^2$ mass range. Looking to the future, the Collaboration has developed a set of experience-based upgrade scenarios, as well as novel directions, to extend the search for dark matter using the SuperCDMS technology in the SNOLAB facility. The experienced-ba… ▽ More The SuperCDMS Collaboration is currently building SuperCDMS SNOLAB, a dark matter search focused on nucleon-coupled dark matter in the 1-5 GeV/c$^2$ mass range. Looking to the future, the Collaboration has developed a set of experience-based upgrade scenarios, as well as novel directions, to extend the search for dark matter using the SuperCDMS technology in the SNOLAB facility. The experienced-based scenarios are forecasted to probe many square decades of unexplored dark matter parameter space below 5 GeV/c$^2$, covering over 6 decades in mass: 1-100 eV/c$^2$ for dark photons and axion-like particles, 1-100 MeV/c$^2$ for dark-photon-coupled light dark matter, and 0.05-5 GeV/c$^2$ for nucleon-coupled dark matter. They will reach the neutrino fog in the 0.5-5 GeV/c$^2$ mass range and test a variety of benchmark models and sharp targets. The novel directions involve greater departures from current SuperCDMS technology but promise even greater reach in the long run, and their development must begin now for them to be available in a timely fashion. The experienced-based upgrade scenarios rely mainly on dramatic improvements in detector performance based on demonstrated scaling laws and reasonable extrapolations of current performance. Importantly, these improvements in detector performance obviate significant reductions in background levels beyond current expectations for the SuperCDMS SNOLAB experiment. Given that the dominant limiting backgrounds for SuperCDMS SNOLAB are cosmogenically created radioisotopes in the detectors, likely amenable only to isotopic purification and an underground detector life-cycle from before crystal growth to detector testing, the potential cost and time savings are enormous and the necessary improvements much easier to prototype. △ Less

Submitted 1 April, 2023; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: contribution to Snowmass 2021; v2 updated (assorted corrections and improvements to forecasts) October 2022; v3 updated (corrected SuperCDMS SNOLAB sensitivity curves in upgrade forecast plots in body of text) April 2023

EXCESS workshop: Descriptions of rising low-energy spectra

Authors: P. Adari, A. Aguilar-Arevalo, D. Amidei, G. Angloher, E. Armengaud, C. Augier, L. Balogh, S. Banik, D. Baxter, C. Beaufort, G. Beaulieu, V. Belov, Y. Ben Gal, G. Benato, A. Benoît, A. Bento, L. Bergé, A. Bertolini, R. Bhattacharyya, J. Billard, I. M. Bloch, A. Botti, R. Breier, G. Bres, J-. L. Bret , et al. (281 additional authors not shown)

Abstract: Many low-threshold experiments observe sharply rising event rates of yet unknown origins below a few hundred eV, and larger than expected from known backgrounds. Due to the significant impact of this excess on the dark matter or neutrino sensitivity of these experiments, a collective effort has been started to share the knowledge about the individual observations. For this, the EXCESS Workshop was… ▽ More Many low-threshold experiments observe sharply rising event rates of yet unknown origins below a few hundred eV, and larger than expected from known backgrounds. Due to the significant impact of this excess on the dark matter or neutrino sensitivity of these experiments, a collective effort has been started to share the knowledge about the individual observations. For this, the EXCESS Workshop was initiated. In its first iteration in June 2021, ten rare event search collaborations contributed to this initiative via talks and discussions. The contributing collaborations were CONNIE, CRESST, DAMIC, EDELWEISS, MINER, NEWS-G, NUCLEUS, RICOCHET, SENSEI and SuperCDMS. They presented data about their observed energy spectra and known backgrounds together with details about the respective measurements. In this paper, we summarize the presented information and give a comprehensive overview of the similarities and differences between the distinct measurements. The provided data is furthermore publicly available on the workshop's data repository together with a plotting tool for visualization. △ Less

Submitted 4 March, 2022; v1 submitted 10 February, 2022; originally announced February 2022.

Comments: 44 pages, 20 figures; Editors: A. Fuss, M. Kaznacheeva, F. Reindl, F. Wagner; updated copyright statements and funding information

Journal ref: SciPost Phys. Proc. 9, 001 (2022)

IACHEC 2020/2021 Pandemic Report

Authors: K. K. Madsen, V. Burwitz, K. Forster, C. E. Grant, M. Guainazzi, V. Kashyap, H. L. Marshall, E. D. Miller, L. Natalucci, P. P. Plucinsky, Y. Terada

Abstract: In this report we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) and the work done since the last in-person meeting in Japan (Shonan Village Center), May 2019, through two virtual meetings that were held in November 2020 and May 2021. The on-line only meetings divided the contents of the usual in-person workshop between mission updates an… ▽ More In this report we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) and the work done since the last in-person meeting in Japan (Shonan Village Center), May 2019, through two virtual meetings that were held in November 2020 and May 2021. The on-line only meetings divided the contents of the usual in-person workshop between mission updates and working group updates. The November meeting was dedicated to mission calibration updates and the current status of the cross-calibration between NuSTAR, Swift, and NICER, which frequently join together in observations of bright transients, and a review of the XMM-Newton and Chandra cross-calibration. Results between \nustar\ and \swift\ overall show good agreement, but issues persist in the overlap region 3--5 keV for bright source with large dust scattering halos. The NICER cross-calibration is still progressing and evolving, while for the XMM-Newton and Chandra cross-calibration systematic differences both in the absolute flux and spectral shape determination still exists on different classes of sources. The meeting in May was focused on the Working Group progress and reports summarized here. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: Yearly activity report of the IACHEC, 16 pages, 2 figures

Flows in Enthalpy Based Thermal Evolution of Loops

Authors: Abhishek Rajhans, Durgesh Tripathi, Stephen Bradshaw, Vinay Kashyap, James Klimchuk

Abstract: Plasma-filled loop structures are common in the solar corona. Because detailed modeling of the dynamical evolution of these structures is computationally costly, an efficient method for computing approximate but quick physics-based solutions is to rely on space integrated 0D simulations. The enthalpy-based thermal evolution of loops EBTEL framework is a commonly used method to study the exchange o… ▽ More Plasma-filled loop structures are common in the solar corona. Because detailed modeling of the dynamical evolution of these structures is computationally costly, an efficient method for computing approximate but quick physics-based solutions is to rely on space integrated 0D simulations. The enthalpy-based thermal evolution of loops EBTEL framework is a commonly used method to study the exchange of mass and energy between the corona and transition region. EBTEL solves for density, temperature, and pressure averaged over the coronal part of the loop, the velocity at the coronal base, and the instantaneous differential emission measure distribution in the transition region. The current single-fluid version of the code, EBTEL2, assumes that at all stages the flows are subsonic. However, sometimes the solutions show the presence of supersonic flows during the impulsive phase of heat input. It is thus necessary to account for this effect. Here, we upgrade EBTEL2 to EBTEL3 by including the kinetic energy term in the Navier-Stokes equation. We compare the solutions from EBTEL3 with those obtained using EBTEL2 as well as the state-of-the-art field-aligned hydrodynamics code HYDRAD. We find that the match in pressure between EBTEL3 and HYDRAD is better than that between EBTEL2 and HYDRAD. Additionally, the velocities predicted by EBTEL3 are in close agreement with those obtained with HYDRAD when the flows are subsonic. However, EBTEL3 solutions deviate substantially from HYDRAD's when the latter predicts supersonic flows. Using the mismatches in the solution, we propose a criterion to determine the conditions under which EBTEL can be used to study flows in the system. △ Less

Submitted 7 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ

Concordance: In-flight Calibration of X-ray Telescopes without Absolute References

Authors: Herman L. Marshall, Yang Chen, Jeremy J. Drake, Matteo Guainazzi, Vinay L. Kashyap, Xiao-Li Meng, Paul P. Plucinsky, Peter Ratzlaff, David A. van Dyk, Xufei Wang

Abstract: We describe a process for cross-calibrating the effective areas of X-ray telescopes that observe common targets. The targets are not assumed to be "standard candles" in the classic sense, in that we assume that the source fluxes have well-defined, but {\it a priori} unknown values. Using a technique developed by Chen et al. (2019, arXiv:1711.09429) that involves a statistical method called {\em sh… ▽ More We describe a process for cross-calibrating the effective areas of X-ray telescopes that observe common targets. The targets are not assumed to be "standard candles" in the classic sense, in that we assume that the source fluxes have well-defined, but {\it a priori} unknown values. Using a technique developed by Chen et al. (2019, arXiv:1711.09429) that involves a statistical method called {\em shrinkage estimation}, we determine effective area correction factors for each instrument that brings estimated fluxes into the best agreement, consistent with prior knowledge of their effective areas. We expand the technique to allow unique priors on systematic uncertainties in effective areas for each X-ray astronomy instrument and to allow correlations between effective areas in different energy bands. We demonstrate the method with several data sets from various X-ray telescopes. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 37 pages, 9 figures, and 10 tables plus an appendix with 7 data tables. Accepted for publication in the Astronomical Journal

Hydrodynamics of small transient brightenings in Solar corona

Authors: Abhishek Rajhans, Durgesh Tripathi, Vinay L. Kashyap

Abstract: Small scale transients occur in the Solar corona at much higher frequencies than flares and play a significant role in coronal dynamics. Here we study three well-identified transients discovered by Hi-C and also detected by the EUV channels of Atmospheric Imaging Assembly (AIA) on board Solar Dynamics Observatory (SDO). We use 0-D enthalpy-based hydrodynamical simulations and produce synthetic lig… ▽ More Small scale transients occur in the Solar corona at much higher frequencies than flares and play a significant role in coronal dynamics. Here we study three well-identified transients discovered by Hi-C and also detected by the EUV channels of Atmospheric Imaging Assembly (AIA) on board Solar Dynamics Observatory (SDO). We use 0-D enthalpy-based hydrodynamical simulations and produce synthetic light curves to compare with AIA observations. We have modeled these transients as loops of ~ 1.0~Mm length depositing energies ~ 10^23 ergs in ~ 50 seconds. The simulated synthetic light curves show reasonable agreement with the observed light curves. During the initial phase, conduction flux from the corona dominates over the radiation, like impulsive flaring events. Our results further show that the time-integrated net enthalpy flux is positive, hence into the corona. The fact that we can model the observed light curves of these transients reasonably well by using the same physics as those for nanoflares, microflares, and large flares, suggests that these transients may have a common origin. △ Less

Submitted 30 September, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: Published in ApJ, 16 pages, 13 figures, and 3 tables

Journal ref: The Astrophysical Journal, 917:29 (13pp), 2021 August 10

eBASCS: Disentangling Overlapping Astronomical Sources II, using Spatial, Spectral, and Temporal Information

Authors: Antoine D. Meyer, David A. van Dyk, Vinay L. Kashyap, Luis F. Campos, David E. Jones, Aneta Siemiginowska, Andreas Zezas

Abstract: The analysis of individual X-ray sources that appear in a crowded field can easily be compromised by the misallocation of recorded events to their originating sources. Even with a small number of sources, that nonetheless have overlapping point spread functions, the allocation of events to sources is a complex task that is subject to uncertainty. We develop a Bayesian method designed to sift high-… ▽ More The analysis of individual X-ray sources that appear in a crowded field can easily be compromised by the misallocation of recorded events to their originating sources. Even with a small number of sources, that nonetheless have overlapping point spread functions, the allocation of events to sources is a complex task that is subject to uncertainty. We develop a Bayesian method designed to sift high-energy photon events from multiple sources with overlapping point spread functions, leveraging the differences in their spatial, spectral, and temporal signatures. The method probabilistically assigns each event to a given source. Such a disentanglement allows more detailed spectral or temporal analysis to focus on the individual component in isolation, free of contamination from other sources or the background. We are also able to compute source parameters of interest like their locations, relative brightness, and background contamination, while accounting for the uncertainty in event assignments. Simulation studies that include event arrival time information demonstrate that the temporal component improves event disambiguation beyond using only spatial and spectral information. The proposed methods correctly allocate up to 65% more events than the corresponding algorithms that ignore event arrival time information. We apply our methods to two stellar X-ray binaries, UV Cet and HBC515 A, observed with Chandra. We demonstrate that our methods are capable of removing the contamination due to a strong flare on UV Cet B in its companion approximately 40 times weaker during that event, and that evidence for spectral variability at timescales of a few ks can be determined in HBC515 Aa and HBC515 Ab. △ Less

Submitted 18 May, 2021; originally announced May 2021.

Change point detection and image segmentation for time series of astrophysical images

Authors: Cong Xu, Hans Moritz Günther, Vinay L. Kashyap, Thomas C. M. Lee, Andreas Zezas

Abstract: Many astrophysical phenomena are time-varying, in the sense that their intensity, energy spectrum, and/or the spatial distribution of the emission suddenly change. This paper develops a method for modeling a time series of images. Under the assumption that the arrival times of the photons follow a Poisson process, the data are binned into 4D grids of voxels (time, energy band, and x-y coordinates)… ▽ More Many astrophysical phenomena are time-varying, in the sense that their intensity, energy spectrum, and/or the spatial distribution of the emission suddenly change. This paper develops a method for modeling a time series of images. Under the assumption that the arrival times of the photons follow a Poisson process, the data are binned into 4D grids of voxels (time, energy band, and x-y coordinates), and viewed as a time series of non-homogeneous Poisson images. The method assumes that at each time point, the corresponding multi-band image stack is an unknown 3D piecewise constant function including Poisson noise. It also assumes that all image stacks between any two adjacent change points (in time domain) share the same unknown piecewise constant function. The proposed method is designed to estimate the number and the locations of all the change points (in time domain), as well as all the unknown piecewise constant functions between any pairs of the change points. The method applies the minimum description length (MDL) principle to perform this task. A practical algorithm is also developed to solve the corresponding complicated optimization problem. Simulation experiments and applications to real datasets show that the proposed method enjoys very promising empirical properties. Applications to two real datasets, the XMM observation of a flaring star and an emerging solar coronal loop, illustrate the usage of the proposed method and the scientific insight gained from it. △ Less

Submitted 26 January, 2021; originally announced January 2021.

Comments: 22 pages, 10 figures

Constraints on Lightly Ionizing Particles from CDMSlite

Authors: SuperCDMS Collaboration, I. Alkhatib, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, D. Barker, C. Bathurst, D. A. Bauer, L. V. S. Bezerra, R. Bhattacharyya, M. A. Bowles, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen , et al. (93 additional authors not shown)

Abstract: The Cryogenic Dark Matter Search low ionization threshold experiment (CDMSlite) achieved efficient detection of very small recoil energies in its germanium target, resulting in sensitivity to Lightly Ionizing Particles (LIPs) in a previously unexplored region of charge, mass, and velocity parameter space. We report first direct-detection limits calculated using the optimum interval method on the v… ▽ More The Cryogenic Dark Matter Search low ionization threshold experiment (CDMSlite) achieved efficient detection of very small recoil energies in its germanium target, resulting in sensitivity to Lightly Ionizing Particles (LIPs) in a previously unexplored region of charge, mass, and velocity parameter space. We report first direct-detection limits calculated using the optimum interval method on the vertical intensity of cosmogenically-produced LIPs with an electric charge smaller than $e/(3\times10^5$), as well as the strongest limits for charge $\leq e/160$, with a minimum vertical intensity of $1.36\times10^{-7}$\,cm$^{-2}$s$^{-1}$sr$^{-1}$ at charge $e/160$. These results apply over a wide range of LIP masses (5\,MeV/$c^2$ to 100\,TeV/$c^2$) and cover a wide range of $βγ$ values (0.1 -- $10^6$), thus excluding non-relativistic LIPs with $βγ$ as small as 0.1 for the first time. △ Less

Submitted 19 February, 2022; v1 submitted 18 November, 2020; originally announced November 2020.

Journal ref: Phys. Rev. Lett. 127, 081802 (2021)

Long-term X-ray Variability of the Symbiotic System RT Cru based on Chandra Spectroscopy

Authors: A. Danehkar, M. Karovska, J. J. Drake, V. L. Kashyap

Abstract: RT Cru belongs to the rare class of hard X-ray emitting symbiotics, whose origin is not yet fully understood. In this work, we have conducted a detailed spectroscopic analysis of X-ray emission from RT Cru based on observations taken by the Chandra Observatory using the Low Energy Transmission Grating (LETG) on the High-Resolution Camera Spectrometer (HRC-S) in 2015 and the High Energy Transmissio… ▽ More RT Cru belongs to the rare class of hard X-ray emitting symbiotics, whose origin is not yet fully understood. In this work, we have conducted a detailed spectroscopic analysis of X-ray emission from RT Cru based on observations taken by the Chandra Observatory using the Low Energy Transmission Grating (LETG) on the High-Resolution Camera Spectrometer (HRC-S) in 2015 and the High Energy Transmission Grating (HETG) on the Advanced CCD Imaging Spectrometer S-array (ACIS-S) in 2005. Our thermal plasma modeling of the time-averaged HRC-S/LETG spectrum suggests a mean temperature of $kT \sim 1.3$ keV, whereas $kT \sim 9.6$ keV according to the time-averaged ACIS-S/HETG. The soft thermal plasma emission component ($\sim1.3$ keV) found in the HRC-S is heavily obscured by dense materials ($> 5 \times 10^{23}$ cm$^{-2}$). The aperiodic variability seen in its light curves could be due to changes in either absorbing material covering the hard X-ray source or intrinsic emission mechanism in the inner layers of the accretion disk. To understand the variability, we extracted the spectra in the "low/hard" and "high/soft" spectral states, which indicated higher plasma temperatures in the low/hard states of both the ACIS-S and HRC-S. The source also has a fluorescent iron emission line at 6.4 keV, likely emitted from reflection off an accretion disk or dense absorber, which was twice as bright in the HRC-S epoch compared to the ACIS-S. The soft thermal component identified in the HRC-S might be an indication of a jet that deserves further evaluations using high-resolution imaging observations. △ Less

Submitted 14 December, 2020; v1 submitted 14 November, 2020; originally announced November 2020.

Comments: 28 pages, 18 figures, 4 tables, published in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.500:4801-4817,2021

Simulation of neutron background for a dark matter search experiment at JUSL

Authors: S. Banik, V. K. S. Kashyap, S. Ghosh, S. Dutta, B. Mohanty, Meghna K. K., P. Bhattacharjee, S. Saha

Abstract: Dark matter search experiments demand low to ultralow radiation background to operate. It is very important to understand the nature of the radiation background including knowledge about the sources contributing to it. Sometimes, evaluation of the background becomes very specific to the site chosen for the experiment, and also to the experimental configuration. A dark matter search experiment is p… ▽ More Dark matter search experiments demand low to ultralow radiation background to operate. It is very important to understand the nature of the radiation background including knowledge about the sources contributing to it. Sometimes, evaluation of the background becomes very specific to the site chosen for the experiment, and also to the experimental configuration. A dark matter search experiment is proposed to be set up at the Jaduguda Underground Science Laboratory (JUSL) in India. The laboratory will be located inside an existing mine with 555 m of vertical rock overburden. Neutrons produced from $(α,n)$ reactions, spontaneous fission of natural radioactive impurities in the rocks, and also from cosmic muon induced reactions are considered as the main background which can affect the sensitivity and outcome of the experiment. In this work, simulations based on GEANT4 are done to understand both the radiogenic neutron background caused by natural radioactivity of the surrounding rock and the cosmogenic neutron background due to interaction of the deeply penetrating cosmic muons with the rock. The muon flux in the cavern is obtained to be $4.49(\pm0.25)\times10^{-7} $cm$^{-2}$s$^{-1}$ and the fluxes of radiogenic and cosmogenic neutrons above an energy threshold of 1 MeV in the cavern are obtained to be $5.75(\pm0.58)\times10^{-6}$cm$^{-2}$s$^{-1}$ and $7.25(\pm0.65)\times 10^{-9}$ cm$^{-2}$s$^{-1}$ respectively. The values obtained are comparable with estimates and measurements done for DAMA, WIPP and dark matter experiments at Boulby mine. The effectiveness of different shielding materials are also investigated to obtain the best possible neutron background reduction for a dark matter search experiment at JUSL. We also estimate the sensitivity of a CsI based detector for Weakly Interacting Massive Particle (WIMP) dark matter search at JUSL considering the estimated neutron background. △ Less

Submitted 10 November, 2020; originally announced November 2020.

M51-ULS-1b: The First Candidate for a Planet in an External Galaxy

Authors: R. Di Stefano, Julia Berndtsson, Ryan Urquhart, Roberto Soria, Vinay L. Kashyap, Theron W. Carmichael, Nia Imara

Abstract: Do external galaxies host planetary systems? Many lines of reasoning suggest that the answer must be 'yes'. In the foreseeable future, however, the question cannot be answered by the methods most successful in our own Galaxy. We report on a different approach which focuses on bright X-ray sources (XRSs). M51-ULS-1b is the first planet candidate to be found because it produces a full, short-lived e… ▽ More Do external galaxies host planetary systems? Many lines of reasoning suggest that the answer must be 'yes'. In the foreseeable future, however, the question cannot be answered by the methods most successful in our own Galaxy. We report on a different approach which focuses on bright X-ray sources (XRSs). M51-ULS-1b is the first planet candidate to be found because it produces a full, short-lived eclipse of a bright XRS. M51-ULS-1b has a most probable radius slightly smaller than Saturn. It orbits one of the brightest XRSs in the external galaxy M51, the Whirlpool Galaxy, located 8.6 Megaparsecs from Earth. It is the first candidate for a planet in an external galaxy. The binary it orbits, M51-ULS-1, is young and massive. One of the binary components is a stellar remnant, either a neutron star (NS) or black hole (BH), and the other is a massive star. X-ray transits can now be used to discover more planets in external galaxies and also planets orbiting XRSs inside the Milky Way. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Comments: 19 pages, 13 figures, 2 tables. Submitted

Light Dark Matter Search with a High-Resolution Athermal Phonon Detector Operated Above Ground

Authors: I. Alkhatib, D. W. P. Amaral, T. Aralis, T. Aramaki, I. J. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, D. Barker, C. Bathurst, D. A. Bauer, L. V. S. Bezerra, R. Bhattacharyya, T. Binder, M. A. Bowles, P. L. Brink, R. Bunker, B. Cabrera, R. Calkins, R. A. Cameron, C. Cartaro, D. G. Cerdeño, Y. -Y. Chang, M. Chaudhuri, R. Chen , et al. (99 additional authors not shown)

Abstract: We present limits on spin-independent dark matter-nucleon interactions using a $10.6$ $\mathrm{g}$ Si athermal phonon detector with a baseline energy resolution of $σ_E=3.86 \pm 0.04$ $(\mathrm{stat.})^{+0.19}_{-0.00}$ $(\mathrm{syst.})$ $\mathrm{eV}$. This exclusion analysis sets the most stringent dark matter-nucleon scattering cross-section limits achieved by a cryogenic detector for dark matte… ▽ More We present limits on spin-independent dark matter-nucleon interactions using a $10.6$ $\mathrm{g}$ Si athermal phonon detector with a baseline energy resolution of $σ_E=3.86 \pm 0.04$ $(\mathrm{stat.})^{+0.19}_{-0.00}$ $(\mathrm{syst.})$ $\mathrm{eV}$. This exclusion analysis sets the most stringent dark matter-nucleon scattering cross-section limits achieved by a cryogenic detector for dark matter particle masses from $93$ to $140$ $\mathrm{MeV}/c^2$, with a raw exposure of $9.9$ $\mathrm{g}\cdot\mathrm{d}$ acquired at an above-ground facility. This work illustrates the scientific potential of detectors with athermal phonon sensors with eV-scale energy resolution for future dark matter searches. △ Less

Submitted 12 October, 2021; v1 submitted 21 July, 2020; originally announced July 2020.

Comments: 7 pages, 4 figures, this version includes ancillary files from official data release

Journal ref: Phys. Rev. Lett. 127, 061801 (2021)

X-ray Properties of Young Radio Quasars at z > 4.5

Authors: Bradford Snios, Aneta Siemiginowska, Małgosia Sobolewska, C. C. Cheung, Vinay Kashyap, Giulia Migliori, Daniel A. Schwartz, Łukasz Stawarz, Diana M. Worrall

Abstract: We present a comprehensive analysis of Chandra X-ray observations of 15 young radio quasars at redshifts $4.5 < z < 5.0$. All sources are detected in the $0.5-7.0$ keV energy band. Emission spectra are extracted, and the average photon index for the sample is measured to be $1.5\pm0.1$. Unabsorbed rest-frame $2-10$ keV luminosities are found to range between $(0.5-23.2) \times 10^{45}$ erg s… ▽ More We present a comprehensive analysis of Chandra X-ray observations of 15 young radio quasars at redshifts $4.5 < z < 5.0$. All sources are detected in the $0.5-7.0$ keV energy band. Emission spectra are extracted, and the average photon index for the sample is measured to be $1.5\pm0.1$. Unabsorbed rest-frame $2-10$ keV luminosities are found to range between $(0.5-23.2) \times 10^{45}$ erg s$^{-1}$. The optical-X-ray power-law spectral index $α_{ox}$ is calculated for each source using optical/UV data available in the literature. The $α_{ox}$-UV relationship is compared with other quasar surveys, and an anticorrelation is observed that agrees with independent estimates. Rest-frame radio and X-ray luminosities are established for the sample, and a correlation between the luminosities is detected. These multiwavelength results reinforce a lack of spectral evolution for quasars over a broad redshift range. We additionally identify three quasars from our multiwavelength analysis that are statistically significant outliers, with one source being a Compton-thick candidate in the early universe, and discuss each in detail. △ Less

Submitted 3 August, 2020; v1 submitted 2 July, 2020; originally announced July 2020.

Comments: Accepted to ApJ, 13 pages, 5 figures, 4 tables

Spatially Resolved Chandra Spectroscopy of the Large Magellanic Cloud Supernova Remnant N132D

Authors: Piyush Sharda, Terrance Gaetz, Vinay Kashyap, Paul Plucinsky

Abstract: We perform detailed spectroscopy of the X-ray brightest supernova remnant (SNR) in the Large Magellanic Cloud (LMC), N132D, using Chandra archival observations. By analyzing the spectra of the entire well-defined rim, we determine the mean abundances for O, Ne, Mg, Si, S and Fe for the local LMC environment. We find evidence of enhanced O on the north-western and S on the north-eastern blast wave.… ▽ More We perform detailed spectroscopy of the X-ray brightest supernova remnant (SNR) in the Large Magellanic Cloud (LMC), N132D, using Chandra archival observations. By analyzing the spectra of the entire well-defined rim, we determine the mean abundances for O, Ne, Mg, Si, S and Fe for the local LMC environment. We find evidence of enhanced O on the north-western and S on the north-eastern blast wave. By analyzing spectra interior to the remnant, we confirm the presence of a Si-rich relatively hot plasma (> 1.5 kev) that is also responsible for the Fe K emission. Chandra images show that the Fe K emission is distributed throughout the interior of the southern half of the remnant but does not extend out to the blast wave. We estimate the progenitor mass to be $15\pm5\,M_{\odot}$ using abundance ratios in different regions that collectively cover a large fraction of the remnant, as well as from the radius of the forward shock compared with models of an explosion in a cavity created by stellar winds. We fit ionizing and recombining plasma models to the Fe K emission and find that the current data cannot distinguish between the two, hence the origin of the high-temperature plasma remains uncertain. Our analysis is consistent with N132D being the result of a core-collapse supernova in a cavity created by its intermediate mass progenitor. △ Less

Submitted 15 April, 2020; originally announced April 2020.

Comments: 32 pages, 28 figures, 10 tables, accepted for ApJ

Summary of the 14th IACHEC Meeting

Authors: K. K. Madsen, Y. Terada, V. Burwitz, G. Belanger, C. E. Grant, M. Guainazzi, V. Kashyap, H. L. Marshall, E. D. Miller, L. Natalucci, P. P. Plucinsky

Abstract: We summarize the 14th meeting of the International Astronomical Consortium for High Energy Calibration (IACHEC) held at \textit{Shonan Village} (Kanagawa, Japan) in May 2019. Sixty scientists directly involved in the calibration of operational and future high-energy missions gathered during 3.5 days to discuss the status of the cross-calibration between the current international complement of X-ra… ▽ More We summarize the 14th meeting of the International Astronomical Consortium for High Energy Calibration (IACHEC) held at \textit{Shonan Village} (Kanagawa, Japan) in May 2019. Sixty scientists directly involved in the calibration of operational and future high-energy missions gathered during 3.5 days to discuss the status of the cross-calibration between the current international complement of X-ray observatories, and the possibilities to improve it. This summary consists of reports from the various WGs with topics ranging from the identification and characterization of standard calibration sources, multi-observatory cross-calibration campaigns, appropriate and new statistical techniques, calibration of instruments and characterization of background, communication and preservation of knowledge, and results for the benefit of the astronomical community. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 15 pies

X-ray Observations of the Peculiar Cepheid V473 Lyr Identify A Low-Mass Companion

Authors: Nancy Remage Evans, Ignazio Pillitteri, Laszlo Molnar, Laszlo Szabados, Emese Plachy, Robert Szabo, Scott Engle, Ed Guinan, Scott Wolk, H. Moritz Guenther, Hilding Neilson, Massimo Marengo, Lynn D. Matthews, Sofia Moschou, Jeremy J. Drake, Vinay Kashyap, Pierre Kervella, Tamas Tordai, Peter Somogyi, Gilbert Burki

Abstract: V473 Lyr is a classical Cepheid which is unique in having substantial amplitude variations with a period of approximately 3.3 years, thought to be similar to the Blazhko variations in RR Lyrae stars. We obtained an {\it XMM-Newton} observation of this star to followup a previous detection in X-rays. Rather than the X-ray burst and rapid decline near maximum radius seen in $δ$ Cephei itself, the X-… ▽ More V473 Lyr is a classical Cepheid which is unique in having substantial amplitude variations with a period of approximately 3.3 years, thought to be similar to the Blazhko variations in RR Lyrae stars. We obtained an {\it XMM-Newton} observation of this star to followup a previous detection in X-rays. Rather than the X-ray burst and rapid decline near maximum radius seen in $δ$ Cephei itself, the X-ray flux in V473 Lyr remained constant for a third of the pulsation cycle covered by the observation. Thus the X-rays are most probably not produced by the changes around the pulsation cycle. The X-ray spectrum is soft (kT = 0.6 keV), with X-ray properties which are consistent with a young low mass companion. Previously there was no evidence of a companion in radial velocities or in {\it Gaia} and {\it Hipparcos} proper motions. While this rules out companions which are very close or very distant, a binary companion at a separation between 30 and 300 AU is possible. This is an example of an X-ray observation revealing evidence of a low mass companion, which is important in completing the mass ratio statistics of binary Cepheids. Furthermore, the detection of a young X-ray bright companion is a further indication that the Cepheid (primary) is a Population I star, even though its pulsation behavior differs from other classical Cepheids. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: Accepted by AJ

Algorithms and Statistical Models for Scientific Discovery in the Petabyte Era

Authors: Brian Nord, Andrew J. Connolly, Jamie Kinney, Jeremy Kubica, Gautaum Narayan, Joshua E. G. Peek, Chad Schafer, Erik J. Tollerud, Camille Avestruz, G. Jogesh Babu, Simon Birrer, Douglas Burke, João Caldeira, Douglas A. Caldwell, Joleen K. Carlberg, Yen-Chi Chen, Chuanfei Dong, Eric D. Feigelson, V. Zach Golkhou, Vinay Kashyap, T. S. Li, Thomas Loredo, Luisa Lucie-Smith, Kaisey S. Mandel, J. R. Martínez-Galarza , et al. (13 additional authors not shown)

Abstract: The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our… ▽ More The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/). △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: arXiv admin note: substantial text overlap with arXiv:1905.05116

Report number: FERMILAB-FN-1093-A-AE-SCD

X-ray Studies of Exoplanets: A 2020 Decadal Survey White Paper

Authors: Scott J. Wolk, Jeremy J. Drake, Graziella Branduardi-Raymont, Katja Poppenhaeger, Vladimir Airapetian, Kevin France, Salvatore Sciortino, Ignazio Pillitteri, Rachel A. Osten, Carey M. Lisse, Vinay Kashyap, Brad Wargelin, Brian Wood, Willaim Dunn, David Principe, Moritz Günther, Damian J. Christian, Julian David Alvarado-Gomez, Chuanfei Dong, Lidia Oskinova, Margarita Karovska, Sofia P. Moschou, Peter K. Williams, Randall Smith, Bradford Snios , et al. (6 additional authors not shown)

Abstract: Over the last two decades, the discovery of exoplanets has fundamentally changed our perception of the universe and humanity's place within it. Recent work indicates that a solar system's X-ray and high energy particle environment is of fundamental importance to the formation and development of the atmospheres of close-in planets such as hot Jupiters, and Earth-like planets around M stars. X-ray i… ▽ More Over the last two decades, the discovery of exoplanets has fundamentally changed our perception of the universe and humanity's place within it. Recent work indicates that a solar system's X-ray and high energy particle environment is of fundamental importance to the formation and development of the atmospheres of close-in planets such as hot Jupiters, and Earth-like planets around M stars. X-ray imaging and spectroscopy provide powerful and unique windows into the high energy flux that an exoplanet experiences, and X-ray photons also serve as proxies for potentially transfigurative coronal mass ejections. Finally, if the host star is a bright enough X-ray source, transit measurements akin to those in the optical and infrared are possible and allow for direct characterization of the upper atmospheres of exoplanets. In this brief white paper, we discuss contributions to the study of exoplanets and their environs which can be made by X-ray data of increasingly high quality that are achievable in the next 10--15 years. △ Less

Submitted 8 April, 2019; originally announced April 2019.

Comments: 7 pages with 2 figures. A 2020 Decadal Survey White Paper

High-Energy Photon and Particle Effects onExoplanet Atmospheres and Habitability

Authors: Jeremy J. Drake, Julián D. Alvarado-Gómez, Vladimir Airapetian, P. Wilson Cauley, Costanza Argiroffi, Matthew K. Browning, Damian J. Christian, Ofer Cohen, Lia Corrales, William Danchi, Miguel de Val-Borro, Chuanfei Dong, William Forman, Kevin France, Elena Gallo, Katherine Garcia-Sage, Cecilia Garraffo, Dawn M. Gelino, Guillaume Gronoff, H. Moritz Günther, Graham M. Harper, Raphaëlle D. Haywood, Margarita Karovska, Vinay Kashyap, Joel Kastner , et al. (19 additional authors not shown)

Abstract: It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current inst… ▽ More It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections. △ Less

Submitted 28 March, 2019; originally announced March 2019.

Comments: Astro2020 Decadal Survey Science White Paper

Astro2020 Science White Paper: The Next Decade of Astroinformatics and Astrostatistics

Authors: A. Siemiginowska, G. Eadie, I. Czekala, E. Feigelson, E. B. Ford, V. Kashyap, M. Kuhn, T. Loredo, M. Ntampaka, A. Stevens, A. Avelino, K. Borne, T. Budavari, B. Burkhart, J. Cisewski-Kehe, F. Civano, I. Chilingarian, D. A. van Dyk, G. Fabbiano, D. P. Finkbeiner, D. Foreman-Mackey, P. Freeman, A. Fruscione, A. A. Goodman, M. Graham , et al. (27 additional authors not shown)

Abstract: Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New met… ▽ More Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 Decadal Survey call for science white papers

Increasing the Discovery Space in Astrophysics - A Collation of Six Submitted White Papers

Authors: G. Fabbiano, M. Elvis, A. Accomazzi, G. B. Berriman, N. Brickhouse, S. Bose, D. Carrera, I. Chilingarian, F. Civano, B. Czerny, R. D'Abrusco, B. Diemer, J. Drake, R. Emami Meibody, J. R. Farah, G. G. Fazio, E. Feigelson, F. Fornasini, Jay Gallagher, J. Grindlay, L. Hernquist, D. J. James, M. Karovska, V. Kashyap, D. -W. Kim , et al. (24 additional authors not shown)

Abstract: We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the kn… ▽ More We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey. △ Less

Submitted 18 March, 2019; v1 submitted 15 March, 2019; originally announced March 2019.

White paper for Chandra cool attitude targets (CAT): Stellar activity with TESS and Chandra

Authors: Hans Moritz Günther, David A. Principe, Carl Melis, Kristina Monsch, P. Christian Schneider, S. Czesla, Nicholas J. Wright, Vinay L. Kashyap, J. H. M. M. Schmitt, E. R. Newton, Jeremy J. Drake, D. P. Huenemoerder

Abstract: All cool stars show magnetic activity, and X-ray emission is the hallmark of this activity. Gaining an understanding of activity aids us in answering fundamental questions about stellar astrophysics and in determining the impact of activity on the exoplanets that orbit these stars. Stellar activity is driven by magnetic fields, which are ultimately powered by convection and stellar rotation. Howev… ▽ More All cool stars show magnetic activity, and X-ray emission is the hallmark of this activity. Gaining an understanding of activity aids us in answering fundamental questions about stellar astrophysics and in determining the impact of activity on the exoplanets that orbit these stars. Stellar activity is driven by magnetic fields, which are ultimately powered by convection and stellar rotation. However, the resulting dynamo properties heavily depend on the stellar interior structure and are far from being understood. X-ray radiation can evaporate exoplanet atmospheres and damage organic materials on the planetary surface, reducing the probability that life can form or be sustained, but also provides an important source of energy for prebiotic chemical reactions. Over the next two years, the TESS mission will deliver a catalog of the closest exoplanets, along with rotation periods and activity diagnostics for millions of stars, whether or not they have a planet. We propose to include all cool stars that are TESS targets and bright enough for Chandra observations, as determined by their detection and flux in the ROSAT all-sky survey (RASS), to the list of Chandra Cool Attitude Targets (CATs). For each target, the signal will be sufficient to fit the coronal plasma with at least two temperature components, and compare abundances of groups of elements with low, medium, or high first ionization potential. Similar to the known relation between X-ray luminosity $L_X$ and rotation period, we can correlate stellar properties with coronal temperatures and abundances to constrain models for stellar activity, coronal heating, and stellar dynamos. Detailed X-ray characterization for even a subset of planet-hosting systems would dramatically advance our knowledge of what impact these emissions have on orbiting planets. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Comments: White paper for Chandra cool attitude targets (4 pages)

Pointing Chandra Toward the Extreme Ultraviolet Fluxes of Very Low-Mass Stars

Authors: Jeremy J. Drake, Vinay L. Kashyap, Bradford J. Wargelin, Scott J. Wolk

Abstract: The X-ray and EUV emission of stars plays a key role in the loss and evolution of the atmospheres of their planets. The coronae of dwarf stars later than M6 appear to behave differently to those of earlier spectral types and are more X-ray dim and radio bright. Too faint to have been observed by the Extreme Ultraviolet Explorer, their EUV behavior is currently highly uncertain. We have devised a m… ▽ More The X-ray and EUV emission of stars plays a key role in the loss and evolution of the atmospheres of their planets. The coronae of dwarf stars later than M6 appear to behave differently to those of earlier spectral types and are more X-ray dim and radio bright. Too faint to have been observed by the Extreme Ultraviolet Explorer, their EUV behavior is currently highly uncertain. We have devised a method to use the Chandra X-ray Observatory High Resolution Camera to provide a measure of EUV emission in the 50-170 Å range and have applied it to the M6.5 dwarf LHS 248 in a pilot 10 ks exposure. Analysis with model spectra using simple, idealised coronal emission measure distributions inspired by an analysis of Chandra HETG spectra of the M5.5 dwarf Proxima Cen and results from the literature, finds greatest consistency with a very shallow emission measure distribution slope, $DEM \propto T^{3/2}$ or shallower, in the range $\log T=5.5$-$6.5$. Within $2σ$ confidence, a much wider range of slopes can be accommodated. Model spectra constrained by this method can provide accurate (within a factor of 2-4) synthesis and extrapolation of EUV spectra for wavelengths $<400$-500 Å. At longer wavelengths models are uncertain by an order of magnitude or more, and depend on the details of the emission measure distribution at temperatures $\log T < 5.5$. The method is sensitive to possible incompleteness of plasma radiative loss models in the 30 170 Å range for which re-examination would be warranted. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: Submitted to ApJ

Summary of the 13th IACHEC Meeting

Authors: K. K. Madsen, L. Natalucci, G. Belanger, C. E. Grant, M. Guainazzi, V. Kashyap, H. L. Marshall, E. D. Miller, J. Nevalainen, P. P. Plucinsky, Y. Terada

Abstract: We summarize the outcome of the 13th meeting of the International Astronomical Consortium for High Energy Calibration (IACHEC), held at Tenuta dei Ciclamini (Avigliano Umbro, Italy) in April 2018. Fifty-one scientists directly involved in the calibration of operational and future high-energy missions gathered during 3.5 days to discuss the current status of the X-ray payload inter-calibration and… ▽ More We summarize the outcome of the 13th meeting of the International Astronomical Consortium for High Energy Calibration (IACHEC), held at Tenuta dei Ciclamini (Avigliano Umbro, Italy) in April 2018. Fifty-one scientists directly involved in the calibration of operational and future high-energy missions gathered during 3.5 days to discuss the current status of the X-ray payload inter-calibration and possible approaches to improve it. This summary consists of reports from the various working groups with topics ranging from the identification and characterization of standard calibration sources, multi-observatory cross-calibration campaigns, appropriate and new statistical techniques, calibration of instruments and characterization of background, and communication and preservation of knowledge and results for the benefit of the astronomical community. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: 12 pages

X-ray spectral characterization of the young Cygnus OB2 population

Authors: E. Flaccomio, J. F. Albacete-Colombo, J. J. Drake, M. G. Guarcello, V. Kashyap, N. J. Wright, K. Briggs, B. Ercolano, M. Mccollogh, S. Sciortino

Abstract: We analyze the X-ray spectra of the $\sim$8000 sources detected in the Cygnus OB2 Chandra Legacy Survey (Drake et al., this issue), with the goals of characterizing the coronal plasma of the young low-mass stars in the region and estimating their intrinsic X-ray luminosities. We adopt two different strategies for X-ray sources for which more or less than 20 photons were detected. For the brighter… ▽ More We analyze the X-ray spectra of the $\sim$8000 sources detected in the Cygnus OB2 Chandra Legacy Survey (Drake et al., this issue), with the goals of characterizing the coronal plasma of the young low-mass stars in the region and estimating their intrinsic X-ray luminosities. We adopt two different strategies for X-ray sources for which more or less than 20 photons were detected. For the brighter sample we fit the spectra with absorbed isothermal models. In order to limit uncertainties, for most of the fainter Cygnus OB2 members in this sample, we constrain the spectral parameters to characteristic ranges defined from the brightest stars. For X-ray sources with $<$20 net photons we adopt a conversion factor from detected photon flux to intrinsic flux. This was defined, building on the results for the previous sample, as a function of the 20% quantile of the detected photon energy distributions, which we prove to also correlate well with extinction. We then use the X-ray extinction from the spectral fits to constrain the ratio between optical and X-ray extinction toward Cygnus OB2, finding it consistent with standard "Galactic"' values, when properly accounting for systematics. Finally we exploit the large number of sources to constrain the average coronal abundances of several elements, through two different ensemble analyses of the X-ray spectra of low-mass Cygnus OB2 members. We find the pattern of abundances to be largely consistent with that derived for the young stellar coronae in the Orion Nebula Cluster. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Comments: 28 pages, 16 figures, accepted for publication in ApJS as part of the special issue on the Chandra Cygnus OB2 Legacy Survey

Incorporating Uncertainties in Atomic Data Into the Analysis of Solar and Stellar Observations: A Case Study in Fe XIII

Authors: Xixi Yu, Giulio Del Zanna, David C. Stenning, Jessi Cisewski-Kehe, Vinay L. Kashyap, Nathan Stein, David A. van Dyk, Harry P. Warren, Mark A. Weber

Abstract: Information about the physical properties of astrophysical objects cannot be measured directly but is inferred by interpreting spectroscopic observations in the context of atomic physics calculations. Ratios of emission lines, for example, can be used to infer the electron density of the emitting plasma. Similarly, the relative intensities of emission lines formed over a wide range of temperatures… ▽ More Information about the physical properties of astrophysical objects cannot be measured directly but is inferred by interpreting spectroscopic observations in the context of atomic physics calculations. Ratios of emission lines, for example, can be used to infer the electron density of the emitting plasma. Similarly, the relative intensities of emission lines formed over a wide range of temperatures yield information on the temperature structure. A critical component of this analysis is understanding how uncertainties in the underlying atomic physics propagates to the uncertainties in the inferred plasma parameters. At present, however, atomic physics databases do not include uncertainties on the atomic parameters and there is no established methodology for using them even if they did. In this paper we develop simple models for the uncertainties in the collision strengths and decay rates for Fe XIII and apply them to the interpretation of density sensitive lines observed with the EUV Imagining spectrometer (EIS) on Hinode. We incorporate these uncertainties in a Bayesian framework. We consider both a pragmatic Bayesian method where the atomic physics information is unaffected by the observed data, and a fully Bayesian method where the data can be used to probe the physics. The former generally increases the uncertainty in the inferred density by about a factor of 5 compared with models that incorporate only statistical uncertainties. The latter reduces the uncertainties on the inferred densities, but identifies areas of possible systematic problems with either the atomic physics or the observed intensities. △ Less

Submitted 17 September, 2018; originally announced September 2018.

Comments: in press at ApJ

Diffuse X-ray emission in the Cygnus OB2 association

Authors: J. F. Albacete Colombo, J. J. Drake, E. Flaccomio, N. J. Wright, V. Kashyap, M. G. Guarcello, K. Briggs, J. E. Drew, D. M. Fenech, G. Micela, M. McCollough, R. K. Prinja, N. Schneider, S. Sciortino, J. S. Vink

Abstract: We present a large-scale study of diffuse X-ray emission in the nearby massive stellar association Cygnus OB2 as part of the Chandra Cygnus OB2 Legacy Program. We used 40 Chandra X-ray ACIS-I observations covering $\sim$1.0 deg$^2$. After removing 7924 point-like sources detected in our survey, background-corrected X-ray emission, the adaptive smoothing reveals large-scale diffuse X-ray emission.… ▽ More We present a large-scale study of diffuse X-ray emission in the nearby massive stellar association Cygnus OB2 as part of the Chandra Cygnus OB2 Legacy Program. We used 40 Chandra X-ray ACIS-I observations covering $\sim$1.0 deg$^2$. After removing 7924 point-like sources detected in our survey, background-corrected X-ray emission, the adaptive smoothing reveals large-scale diffuse X-ray emission. Diffuse emission was detected in the sub-bands Soft [0.5 : 1.2] and Medium [1.2 : 2.5], and marginally in the Hard [2.5 : 7.0] keV band. From X-ray spectral analysis of stacked spectra we compute a total [0.5 : 7.0 keV] diffuse X-ray luminosity of L$_{\rm x}^{\rm diff}\approx$4.2$\times$10$^{\rm 34}$ erg s$^{-1}$, characterized with plasma temperature components at kT$\approx$ 0.11, 0.40 and 1.18 keV, respectively. The HI absorption column density corresponding to these temperatures has a distribution consistent with N$_{\rm H}$ = 0.43, 0.80 and 1.39 $\times$10$^{22}$ cm$^{-2}$. The extended medium band energy emission likely arises from O-type stellar winds thermalized by wind-wind collisions in the most populated regions of the association, while the soft band emission probably arises from less energetic termination shocks against the surrounding Interstellar-Medium. Super-soft and Soft diffuse emission appears more widely dispersed and intense than the medium band emission. The diffuse X-ray emission is generally spatially coincident with low-extinction regions that we attribute to the ubiquitous influence of powerful stellar winds from massive stars and their interaction with the local Interstellar-Medium. Diffuse X-ray emission is volume-filling, rather than edge-brightened, oppositely to other star-forming regions. We reveal the first observational evidence of X-ray haloes around some evolved massive stars. △ Less

Submitted 4 June, 2018; originally announced June 2018.

Comments: Accepted for publication on ApJSS, Cygnus OB2 Special Issue (2018). 20 pages, 19 figures

The Revolution Revolution: magnetic morphology driven spin-down

Authors: C. Garraffo, J. J. Drake, A. Dotter, J. Choi, D. J. Burke, S. P. Moschou, J. D. Alvarado-Gomez, V. L. Kashyap, O. Cohen

Abstract: Observations of young open clusters show a bimodal distribution of rotation periods that has been difficult to explain with existing stellar spin-down models. Detailed MHD stellar wind simulations have demonstrated that surface magnetic field morphology has a strong influence on wind-driven angular momentum loss. Observations suggest that faster rotating stars store a larger fraction of their magn… ▽ More Observations of young open clusters show a bimodal distribution of rotation periods that has been difficult to explain with existing stellar spin-down models. Detailed MHD stellar wind simulations have demonstrated that surface magnetic field morphology has a strong influence on wind-driven angular momentum loss. Observations suggest that faster rotating stars store a larger fraction of their magnetic flux in higher-order multipolar components of the magnetic field. In this work, we present a new model for stellar spin-down that, for the first time, accounts for the stellar surface magnetic field configuration. We show how a magnetic complexity that evolves from complex toward simple configurations as a star spins down can explain the salient features of stellar rotation evolution, including the bimodal distribution of both slow and fast rotators seen in young open clusters. △ Less

Submitted 15 June, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 8 pages, 4 figures, to appear in ApJ