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

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)

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

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)

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

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.

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)

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

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

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

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

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

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

Energetics of Hi-C EUV Brightenings

Authors: Srividya Subramanian, Vinay L. Kashyap, Durgesh Tripathi, Maria S. Madjarska, J. G. Doyle

Abstract: We study the thermal structure and energetics of the point-like EUV brightenings within a system of fan loops observed in the active region \textsl{AR~11520}. These brightenings were simultaneously observed on 2012 July 11 by the HIgh-resolution Coronal (Hi-C) imager and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). We identified 27 brightenings by automatic… ▽ More We study the thermal structure and energetics of the point-like EUV brightenings within a system of fan loops observed in the active region \textsl{AR~11520}. These brightenings were simultaneously observed on 2012 July 11 by the HIgh-resolution Coronal (Hi-C) imager and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). We identified 27 brightenings by automatically determining intensity enhancements in both Hi-C and AIA ~193~Å light curves. The energetics of these brightenings were studied by using the Differential Emission Measure (DEM) diagnostics. The DEM weighted temperatures of these transients are in the range $\log{T} (K) = 6.2 - 6.6$ with radiative energies ${\approx}10^{24-25}$~ergs and densities ${\approx}$ a few times $10^{9}$~cm$^{-3}$. To the best of our knowledge, these are the smallest brightenings in EUV ever detected. We used these results to determine the mechanism of energy loss in these brightenings. Our analysis reveals that the dominant mechanism of energy loss for all the identified brightenings is conduction rather than radiation. △ Less

Submitted 26 March, 2018; originally announced March 2018.

Comments: 17 pages, 8 Figures, 3 tables, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 615, A47 (2018)

Multidimensional Data Driven Classification of Emission-line Galaxies

Authors: Vasileios Stampoulis, David A. van Dyk, Vinay L. Kashyap, Andreas Zezas

Abstract: We propose a new soft clustering scheme for classifying galaxies in different activity classes using simultaneously 4 emission-line ratios; log([NII ]/Ha), log([SII]/Ha), log([OI]/Ha) and log([OIII]/Hb). We fit 20 multivariate Gaussian distributions to the 4-dimensional distribution of these lines obtained from the Sloan Digital Sky Survey (SDSS) in order to capture local structures and subsequent… ▽ More We propose a new soft clustering scheme for classifying galaxies in different activity classes using simultaneously 4 emission-line ratios; log([NII ]/Ha), log([SII]/Ha), log([OI]/Ha) and log([OIII]/Hb). We fit 20 multivariate Gaussian distributions to the 4-dimensional distribution of these lines obtained from the Sloan Digital Sky Survey (SDSS) in order to capture local structures and subsequently group the multivariate Gaussian distributions to represent the complex multi-dimensional structure of the joint distribution of galaxy spectra in the 4 dimensional line ratio space. The main advantages of this method are the use of all four optical-line ratios simultaneously and the adoption of a clustering scheme. This maximises the available information, avoids contradicting classifications, and treats each class as a distribution resulting in soft classification boundaries and providing the probability for an object to belong to each class. We also introduce linear multi-dimensional decision surfaces using support vector machines based on the classification of our soft clustering scheme. This linear multi-dimensional hard clustering technique shows high classification accuracy with respect to our soft-clustering scheme. △ Less

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

Giant coronal loops dominate the quiescent X-ray emission in rapidly rotating M stars

Authors: O. Cohen, R. Yadav, C. Garraffo, S. H. Saar, S. J. Wolk, V. L. Kashyap, J. J. Drake, I. Pillitteri

Abstract: Observations indicate that magnetic fields in rapidly rotating stars are very strong, on both small and large scales. What is the nature of the resulting corona? Here we seek to shed some light on this question. We use the results of an anelastic dynamo simulation of a rapidly rotating fully-convective M-star to drive a physics-based model for the stellar corona. We find that due to the several ki… ▽ More Observations indicate that magnetic fields in rapidly rotating stars are very strong, on both small and large scales. What is the nature of the resulting corona? Here we seek to shed some light on this question. We use the results of an anelastic dynamo simulation of a rapidly rotating fully-convective M-star to drive a physics-based model for the stellar corona. We find that due to the several kilo Gauss large-scale magnetic fields at high latitudes, the corona and its X-ray emission are dominated by star-size large hot loops, while the smaller, underlying colder loops are not visible much in the X-ray. Based on this result we propose that, in rapidly rotating stars, emission from such coronal structures dominates the quiescent, cooler but saturated X-ray emission. △ Less

Submitted 7 November, 2016; originally announced November 2016.

Comments: 9 pages, 4 figures, accepted to ApJ

Optical, UV, and X-Ray Evidence for a 7-Year Stellar Cycle in Proxima Centauri

Authors: B. J. Wargelin, S. H. Saar, G. Pojmański, J. J. Drake, V. L. Kashyap

Abstract: Stars of stellar type later than about M3.5 are believed to be fully convective and therefore unable to support magnetic dynamos like the one that produces the 11-year solar cycle. Because of their intrinsic faintness, very few late M stars have undergone long-term monitoring to test this prediction, which is critical to our understanding of magnetic field generation in such stars. Magnetic activi… ▽ More Stars of stellar type later than about M3.5 are believed to be fully convective and therefore unable to support magnetic dynamos like the one that produces the 11-year solar cycle. Because of their intrinsic faintness, very few late M stars have undergone long-term monitoring to test this prediction, which is critical to our understanding of magnetic field generation in such stars. Magnetic activity is also of interest as the driver of UV and X-ray radiation, as well as energetic particles and stellar winds, that affect the atmospheres of close-in planets that lie within habitable zones, such as the recently discovered Proxima b. We report here on several years of optical, UV, and X-ray observations of Proxima Centauri (GJ 551; dM5.5e): 15 years of ASAS photometry in the V band (1085 nights) and 3 years in the I band (196 nights), 4 years of Swift XRT and UVOT observations (more than 120 exposures), and 9 sets of X-ray observations from other X-ray missions (ASCA, XMM-Newton, and three Chandra instruments) spanning 22 years. We confirm previous reports of an 83-day rotational period and find strong evidence for a 7-year stellar cycle, along with indications of differential rotation at about the solar level. X-ray/UV intensity is anti-correlated with optical V-band brightness for both rotational and cyclical variations. From comparison with other stars observed to have X-ray cycles we deduce a simple empirical relationship between X-ray cyclic modulation and Rossby number, and we also present Swift UV grism spectra covering 2300-6000 Å. △ Less

Submitted 11 October, 2016; originally announced October 2016.

Comments: 19 pages with 15 figures. Submitted to MNRAS 2016 Jul 1; revised version accepted 2016 Oct 3 and published online 2016 Oct 8. Available at http://dx.doi.org/10.1093/mnras/stw2570

Detecting Relativistic X-ray Jets in High-Redshift Quasars

Authors: Kathryn McKeough, Aneta Siemiginowska, C. C. Cheung, Lukasz Stawarz, Vinay L. Kashyap, Nathan Stein, Vasileios Stampoulis, David A. van Dyk, J. F. C. Wardle, N. P. Lee, D. E. Harris, D. A. Schwartz, Davide Donato, Laura Maraschi, Fabrizio Tavecchio

Abstract: We analyze Chandra X-ray images of a sample of 11 quasars that are known to contain kiloparsec scale radio jets. The sample consists of five high-redshift (z >= 3.6) flat-spectrum radio quasars, and six intermediate redshift (2.1 < z < 2.9) quasars. The dataset includes four sources with integrated steep radio spectra and seven with flat radio spectra. A total of 25 radio jet features are present… ▽ More We analyze Chandra X-ray images of a sample of 11 quasars that are known to contain kiloparsec scale radio jets. The sample consists of five high-redshift (z >= 3.6) flat-spectrum radio quasars, and six intermediate redshift (2.1 < z < 2.9) quasars. The dataset includes four sources with integrated steep radio spectra and seven with flat radio spectra. A total of 25 radio jet features are present in this sample. We apply a Bayesian multi-scale image reconstruction method to detect and measure the X-ray emission from the jets. We compute deviations from a baseline model that does not include the jet, and compare observed X-ray images with those computed with simulated images where no jet features exist. This allows us to compute p-value upper bounds on the significance that an X- ray jet is detected in a pre-determined region of interest. We detected 12 of the features unambiguously, and an additional 6 marginally. We also find residual emission in the cores of 3 quasars and in the background of 1 quasar that suggest the existence of unresolved X-ray jets. The dependence of the X-ray to radio luminosity ratio on redshift is a potential diagnostic of the emission mechanism, since the inverse Compton scattering of cosmic microwave background photons (IC/CMB) is thought to be redshift dependent, whereas in synchrotron models no clear redshift dependence is expected. We find that the high-redshift jets have X-ray to radio flux ratios that are marginally inconsistent with those from lower redshifts, suggesting that either the X-ray emissions is due to the IC/CMB rather than the synchrotron process, or that high redshift jets are qualitatively different. △ Less

Submitted 12 September, 2016; originally announced September 2016.

Comments: 42 pages, 14 figures, submitted to ApJ

Bayesian Estimates of Astronomical Time Delays between Gravitationally Lensed Stochastic Light Curves

Authors: Hyungsuk Tak, Kaisey Mandel, David A. van Dyk, Vinay L. Kashyap, Xiao-Li Meng, Aneta Siemiginowska

Abstract: The gravitational field of a galaxy can act as a lens and deflect the light emitted by a more distant object such as a quasar. Strong gravitational lensing causes multiple images of the same quasar to appear in the sky. Since the light in each gravitationally lensed image traverses a different path length from the quasar to the Earth, fluctuations in the source brightness are observed in the sever… ▽ More The gravitational field of a galaxy can act as a lens and deflect the light emitted by a more distant object such as a quasar. Strong gravitational lensing causes multiple images of the same quasar to appear in the sky. Since the light in each gravitationally lensed image traverses a different path length from the quasar to the Earth, fluctuations in the source brightness are observed in the several images at different times. The time delay between these fluctuations can be used to constrain cosmological parameters and can be inferred from the time series of brightness data or light curves of each image. To estimate the time delay, we construct a model based on a state-space representation for irregularly observed time series generated by a latent continuous-time Ornstein-Uhlenbeck process. We account for microlensing, an additional source of independent long-term extrinsic variability, via a polynomial regression. Our Bayesian strategy adopts a Metropolis-Hastings within Gibbs sampler. We improve the sampler by using an ancillarity-sufficiency interweaving strategy and adaptive Markov chain Monte Carlo. We introduce a profile likelihood of the time delay as an approximation of its marginal posterior distribution. The Bayesian and profile likelihood approaches complement each other, producing almost identical results; the Bayesian method is more principled but the profile likelihood is simpler to implement. We demonstrate our estimation strategy using simulated data of doubly- and quadruply-lensed quasars, and observed data from quasars Q0957+561 and J1029+2623. △ Less

Submitted 30 January, 2017; v1 submitted 2 February, 2016; originally announced February 2016.

Comments: Accepted for publication in the Annals of Applied Statistics

Journal ref: The Annals of Applied Statistics, 2017, 11 (3), 1309-1348

Preprocessing Solar Images while Preserving their Latent Structure

Authors: Nathan M Stein, David A van Dyk, Vinay L Kashyap

Abstract: Telescopes such as the Atmospheric Imaging Assembly aboard the Solar Dynamics Observatory, a NASA satellite, collect massive streams of high resolution images of the Sun through multiple wavelength filters. Reconstructing pixel-by-pixel thermal properties based on these images can be framed as an ill-posed inverse problem with Poisson noise, but this reconstruction is computationally expensive and… ▽ More Telescopes such as the Atmospheric Imaging Assembly aboard the Solar Dynamics Observatory, a NASA satellite, collect massive streams of high resolution images of the Sun through multiple wavelength filters. Reconstructing pixel-by-pixel thermal properties based on these images can be framed as an ill-posed inverse problem with Poisson noise, but this reconstruction is computationally expensive and there is disagreement among researchers about what regularization or prior assumptions are most appropriate. This article presents an image segmentation framework for preprocessing such images in order to reduce the data volume while preserving as much thermal information as possible for later downstream analyses. The resulting segmented images reflect thermal properties but do not depend on solving the ill-posed inverse problem. This allows users to avoid the Poisson inverse problem altogether or to tackle it on each of $\sim$10 segments rather than on each of $\sim$10$^7$ pixels, reducing computing time by a factor of $\sim$10$^6$. We employ a parametric class of dissimilarities that can be expressed as cosine dissimilarity functions or Hellinger distances between nonlinearly transformed vectors of multi-passband observations in each pixel. We develop a decision theoretic framework for choosing the dissimilarity that minimizes the expected loss that arises when estimating identifiable thermal properties based on segmented images rather than on a pixel-by-pixel basis. We also examine the efficacy of different dissimilarities for recovering clusters in the underlying thermal properties. The expected losses are computed under scientifically motivated prior distributions. Two simulation studies guide our choices of dissimilarity function. We illustrate our method by segmenting images of a coronal hole observed on 26 February 2015. △ Less

Submitted 14 December, 2015; originally announced December 2015.

Simulating the sensitivity to stellar point sources of Chandra X-ray observations

Authors: Nicholas J. Wright, Jeremy J. Drake, Mario G. Guarcello, Vinay L. Kashyap, Andreas Zezas

Abstract: The Chandra Cygnus OB2 Legacy Survey is a wide and deep X-ray survey of the nearby and massive Cygnus OB2 association. The survey has detected ~8,000 X-ray sources, the majority of which are pre-main sequence X-ray emitting young stars in the association itself. To facilitate quantitative scientific studies of these sources as well as the underlying OB association it is important to understand the… ▽ More The Chandra Cygnus OB2 Legacy Survey is a wide and deep X-ray survey of the nearby and massive Cygnus OB2 association. The survey has detected ~8,000 X-ray sources, the majority of which are pre-main sequence X-ray emitting young stars in the association itself. To facilitate quantitative scientific studies of these sources as well as the underlying OB association it is important to understand the sensitivity of the observations and the level of completeness the observations have obtained. Here we describe the use of a hierarchical Monte Carlo simulation to achieve this goal by combining the empirical properties of the observations, analytic estimates of the source verification process, and an extensive set of source detection simulations. We find that our survey reaches a 90% completeness level for a pre-main-sequence population at the distance of Cyg OB2 at an X-ray luminosity of 4 x 10^30 ergs/s and a stellar mass of 1.3 Msun for a randomly distributed population. For a spatially clustered population such as Cyg~OB2 the 90% completeness level is reached at 1.1 Msun instead, as the sources are more concentrated in areas of our survey with a high exposure. These simulations can easily be adapted for use with other X-ray observations and surveys, and we provide X-ray detection efficiency curves for a very wide array of source and background properties to allow these simulations to be easily exploited by other users. △ Less

Submitted 12 November, 2015; originally announced November 2015.

Comments: Accepted for publication in ApJS as part of the special issue on the Chandra Cygnus OB2 Legacy Survey. 10 pages, 7 figures, 1 table. A full version of Table 1 is available at http://www.astro.keele.ac.uk/~njw/wright2015b_table1.dat

Detecting Unspecified Structure in Low-Count Images

Authors: Nathan M. Stein, David A. van Dyk, Vinay L. Kashyap, Aneta Siemiginowska

Abstract: Unexpected structure in images of astronomical sources often presents itself upon visual inspection of the image, but such apparent structure may either correspond to true features in the source or be due to noise in the data. This paper presents a method for testing whether inferred structure in an image with Poisson noise represents a significant departure from a baseline (null) model of the ima… ▽ More Unexpected structure in images of astronomical sources often presents itself upon visual inspection of the image, but such apparent structure may either correspond to true features in the source or be due to noise in the data. This paper presents a method for testing whether inferred structure in an image with Poisson noise represents a significant departure from a baseline (null) model of the image. To infer image structure, we conduct a Bayesian analysis of a full model that uses a multiscale component to allow flexible departures from the posited null model. As a test statistic, we use a tail probability of the posterior distribution under the full model. This choice of test statistic allows us to estimate a computationally efficient upper bound on a p-value that enables us to draw strong conclusions even when there are limited computational resources that can be devoted to simulations under the null model. We demonstrate the statistical performance of our method on simulated images. Applying our method to an X-ray image of the quasar 0730+257, we find significant evidence against the null model of a single point source and uniform background, lending support to the claim of an X-ray jet. △ Less

Submitted 15 October, 2015; originally announced October 2015.

Journal ref: The Astrophysical Journal (2015), 813, 66 (15pp)

Detecting Abrupt Changes in the Spectra of High-Energy Astrophysical Sources

Authors: Raymond K. W. Wong, Vinay L. Kashyap, Thomas C. M. Lee, David A. van Dyk

Abstract: Variable-intensity astronomical sources are the result of complex and often extreme physical processes. Abrupt changes in source intensity are typically accompanied by equally sudden spectral shifts, i.e., sudden changes in the wavelength distribution of the emission. This article develops a method for modeling photon counts collected from observation of such sources. We embed change points into a… ▽ More Variable-intensity astronomical sources are the result of complex and often extreme physical processes. Abrupt changes in source intensity are typically accompanied by equally sudden spectral shifts, i.e., sudden changes in the wavelength distribution of the emission. This article develops a method for modeling photon counts collected from observation of such sources. We embed change points into a marked Poisson process, where photon wavelengths are regarded as marks and both the Poisson intensity parameter and the distribution of the marks are allowed to change. To the best of our knowledge this is the first effort to embed change points into a marked Poisson process. Between the change points, the spectrum is modeled non-parametrically using a mixture of a smooth radial basis expansion and a number of local deviations from the smooth term representing spectral emission lines. Because the model is over parameterized we employ an $\ell_1$ penalty. The tuning parameter in the penalty and the number of change points are determined via the minimum description length principle. Our method is validated via a series of simulation studies and its practical utility is illustrated in the analysis of the ultra-fast rotating yellow giant star known as FK Com. △ Less

Submitted 10 December, 2015; v1 submitted 27 August, 2015; originally announced August 2015.

Comments: 30 pages, 6 figures

X-ray Evidence for a Pole-Dominated Corona on AB Dor

Authors: Jeremy J. Drake, Sun Mi Chung, Vinay L. Kashyap, David Garcia-Alvarez

Abstract: A fine analysis of spectral line widths and Doppler shifts employing Fourier transform and cross-correlation techniques has been applied to Chandra HETG spectra obtained in 1999 October of the rapidly rotating young star AB Doradus in order to investigate its coronal topology. The observation lasted 52.3ks, covering 1.2 rotation periods. The X-ray light curve obtained from integrating the disperse… ▽ More A fine analysis of spectral line widths and Doppler shifts employing Fourier transform and cross-correlation techniques has been applied to Chandra HETG spectra obtained in 1999 October of the rapidly rotating young star AB Doradus in order to investigate its coronal topology. The observation lasted 52.3ks, covering 1.2 rotation periods. The X-ray light curve obtained from integrating the dispersed signal revealed a moderate intensity flare midway through the exposure in which the count rate increased sharply by about 50% and subsequently decayed over the next 10ks. We find no significant Doppler shifts in the spectra or modulation of the light curve that could be attributed to rotation of dominant coronal structures at this epoch. Individual spectral line widths are statistically consistent with thermal broadening and formally require no rotational broadening, while the $1σ$ limit to rotational broadening corresponds to a compact corona restricted to latitudes $>57°$. Fourier analysis suggests a small amount of rotational broadening is present consistent with a corona restricted largely to the poles, and excludes models with surface rotational broadening or greater. These results present direct spectroscopic evidence that the dominant coronal activity on rapidly-rotating active stars is associated with the dark polar spots commonly seen in photospheric Doppler images, and support models in which these spots are of mixed magnetic polarity that forms closed loops. △ Less

Submitted 23 January, 2015; originally announced January 2015.

Comments: To appear in the Astrophysical Journal

Optical and infrared counterparts of the X-ray sources detected in the Chandra Cygnus OB2 Legacy Survey

Authors: M. G. Guarcello, J. J. Drake, N. J. Wright, T. Naylor, E. Flaccomio, V. L. Kashyap, D. Garcia-Alvarez

Abstract: The young massive OB association Cygnus OB2, in the Cygnus X complex, is the closest (1400 pc) star forming region to the Sun hosting thousands of young low mass stars and up to 1000 OB stars, among which are some of the most massive stars known in our Galaxy. This region holds great importance for several fields of modern astrophysics, such as the study of the physical properties of massive and y… ▽ More The young massive OB association Cygnus OB2, in the Cygnus X complex, is the closest (1400 pc) star forming region to the Sun hosting thousands of young low mass stars and up to 1000 OB stars, among which are some of the most massive stars known in our Galaxy. This region holds great importance for several fields of modern astrophysics, such as the study of the physical properties of massive and young low-mass stars and the feedback provided by massive stars on star and planet formation process. Cygnus OB2 has been recently observed with Chandra/ACIS-I as part of the 1.08Msec Chandra Cygnus OB2 Legacy Project. This survey detected 7924 X-ray sources in a square degree area centered on Cyg OB2. Since a proper classification and study of the observed X-ray sources also requires the analysis of their optical and infrared counterparts, we combined a large and deep set of optical and infrared catalogs available for this region with our new X-ray catalog. In this paper we describe the matching procedure and present the combined catalog containing 5703 sources. We also briefly discuss the nature of the X-ray sources with optical and infrared counterparts using their position in the color-magnitude and color-color diagrams. △ Less

Submitted 15 January, 2015; originally announced January 2015.

Comments: Matched catalog can be provided by the authors; Accepted for publication by ApJ

Disentangling Overlapping Astronomical Sources using Spatial and Spectral Information

Authors: David E. Jones, Vinay L. Kashyap, David A. van Dyk

Abstract: We present a powerful new algorithm that combines both spatial information (event locations and the point spread function) and spectral information (photon energies) to separate photons from overlapping sources. We use Bayesian statistical methods to simultaneously infer the number of overlapping sources, to probabilistically separate the photons among the sources, and to fit the parameters descri… ▽ More We present a powerful new algorithm that combines both spatial information (event locations and the point spread function) and spectral information (photon energies) to separate photons from overlapping sources. We use Bayesian statistical methods to simultaneously infer the number of overlapping sources, to probabilistically separate the photons among the sources, and to fit the parameters describing the individual sources. Using the Bayesian joint posterior distribution, we are able to coherently quantify the uncertainties associated with all these parameters. The advantages of combining spatial and spectral information are demonstrated through a simulation study. The utility of the approach is then illustrated by analysis of observations of FK Aqr and FL Aqr with the XMM-Newton Observatory and the central region of the Orion Nebula Cluster with the Chandra X-ray Observatory. △ Less

Submitted 21 May, 2015; v1 submitted 26 November, 2014; originally announced November 2014.

Comments: 34 pages, 15 figures, references added

Journal ref: Astrophysical Journal, 2015, 808 (2), 137-160

Determining X-Ray Source Intensity and Confidence Bounds in Crowded Fields

Authors: F. A. Primini, V. L. Kashyap

Abstract: We present a rigorous description of the general problem of aperture photometry in high energy astrophysics photon-count images, in which the statistical noise model is Poisson, not Gaussian. We compute the full posterior probability density function for the expected source intensity for various cases of interest, including the important cases in which both source and background apertures contain… ▽ More We present a rigorous description of the general problem of aperture photometry in high energy astrophysics photon-count images, in which the statistical noise model is Poisson, not Gaussian. We compute the full posterior probability density function for the expected source intensity for various cases of interest, including the important cases in which both source and background apertures contain contributions from the source, and when multiple source apertures partially overlap. A Bayesian approach offers the advantages that it allows one to (a) include explicit prior information on source intensities, (b) propagate posterior distributions as priors for future observations, and (c) use Poisson likelihoods, making the treatment valid in the low counts regime. Elements of this approach have been implemented in the Chandra Source Catalog. △ Less

Submitted 25 November, 2014; v1 submitted 9 October, 2014; originally announced October 2014.

Comments: 32 pages, 12 figures. Fixed typos and incorrect statement at end of 4.2.3. Results unchanged

The Chandra Cygnus OB2 Legacy Survey: Design and X-ray Point Source Catalog

Authors: Nicholas J. Wright, Jeremy J. Drake, Mario G. Guarcello, Tom L. Aldcroft, Vinay L. Kashyap, Francesco Damiani, Joe DePasquale, Antonella Fruscione

Abstract: The Cygnus OB2 association is the largest concentration of young and massive stars within 2 kpc of the Sun, including an estimated 65 O-type stars and hundreds of OB stars. The Chandra Cygnus OB2 Legacy Survey is a large imaging program undertaken with the Advanced CCD Imaging Spectrometer onboard the Chandra X-ray Observatory. The survey has imaged the central 0.5 deg^2 of the Cyg OB2 association… ▽ More The Cygnus OB2 association is the largest concentration of young and massive stars within 2 kpc of the Sun, including an estimated 65 O-type stars and hundreds of OB stars. The Chandra Cygnus OB2 Legacy Survey is a large imaging program undertaken with the Advanced CCD Imaging Spectrometer onboard the Chandra X-ray Observatory. The survey has imaged the central 0.5 deg^2 of the Cyg OB2 association with an effective exposure of 120ks and an outer 0.35 deg^2 area with an exposure of 60ks. Here we describe the survey design and observations, the data reduction and source detection, and present a catalog of 8,000 X-ray point sources. The survey design employs a grid of 36 heavily (~50%) overlapping pointings, a method that overcomes Chandra's low off-axis sensitivity and produces a highly uniform exposure over the inner 0.5 deg^2. The full X-ray catalog is described here and is made available online. △ Less

Submitted 27 August, 2014; originally announced August 2014.

Comments: Accepted for publication in ApJS as part of the special issue on the Chandra Cygnus OB2 Legacy Survey. 15 pages, 8 figures, full catalog available on request

The protoplanetary disks in the nearby massive star forming region Cygnus OB2

Authors: M. G. Guarcello, J. J. Drake, N. J. Wright, J. E. Drew, R. A. Gutermuth, J. L. Hora, T. Naylor, T. Aldcroft, A. Fruscione, D. Garcia-Alvarez, V. L. Kashyap, R. King

Abstract: The formation of stars in massive clusters is one of the main modes of the star formation process. However, the study of massive star forming regions is hampered by their typically large distances to the Sun. One exception to this is the massive star forming region Cygnus OB2 in the Cygnus X region, at the distance of about 1400 pc. Cygnus OB2 hosts very rich populations of massive and low-mass st… ▽ More The formation of stars in massive clusters is one of the main modes of the star formation process. However, the study of massive star forming regions is hampered by their typically large distances to the Sun. One exception to this is the massive star forming region Cygnus OB2 in the Cygnus X region, at the distance of about 1400 pc. Cygnus OB2 hosts very rich populations of massive and low-mass stars, being the best target in our Galaxy to study the formation of stars, circumstellar disks, and planets in presence of massive stars. In this paper we combine a wide and deep set of photometric data, from the r band to 24 micron, in order to select the disk bearing population of stars in Cygnus OB2 and identify the class I, class II, and stars with transition and pre-transition disks. We selected 1843 sources with infrared excesses in an area of 1 degree x 1 degree centered on Cyg OB2 in several evolutionary stages: 8.4% class I, 13.1% flat-spectrum sources, 72.9% class II, 2.3% pre-transition disks, and 3.3% transition disks. The spatial distribution of these sources shows a central cluster surrounded by a annular overdensity and some clumps of recent star formation in the outer region. Several candidate subclusters are identified, both along the overdensity and in the rest of the association. △ Less

Submitted 24 June, 2013; originally announced June 2013.

Comments: Accepted by ApJ

Optical photometric GTC/OSIRIS observations of the young massive association Cygnus OB2

Authors: M. G. Guarcello, N. J. Wright, J. J. Drake, D. Garcia-Alvarez, J. E. Drew, T. Aldcroft, V. L. Kashyap

Abstract: In order to fully understand the gravitational collapse of molecular clouds, the star formation process and the evolution of circumstellar disks, these phenomena must be studied in different Galactic environments with a range of stellar contents and positions in the Galaxy. The young massive association Cygnus OB2, in the Cygnus-X region, is an unique target to study how star formation and the evo… ▽ More In order to fully understand the gravitational collapse of molecular clouds, the star formation process and the evolution of circumstellar disks, these phenomena must be studied in different Galactic environments with a range of stellar contents and positions in the Galaxy. The young massive association Cygnus OB2, in the Cygnus-X region, is an unique target to study how star formation and the evolution of circumstellar disks proceed in the presence of a large number of massive stars. We present a catalog obtained with recent optical observations in r,i,z filters with OSIRIS, mounted on the $10.4\,m$ GTC telescope, which is the deepest optical catalog of Cyg OB2 to date. The catalog consist of 64157 sources down to M=0.15 solar masses at the adopted distance and age of Cyg OB2. A total of 38300 sources have good photometry in all three bands. We combined the optical catalog with existing X-ray data of this region, in order to define the cluster locus in the optical diagrams. The cluster locus in the r-i vs. i-z diagram is compatible with an extinction of the optically selected cluster members in the 2.64<AV<5.57 range. We derive an extinction map of the region, finding a median value of AV=4.33 in the center of the association, decreasing toward the north-west. In the color-magnitude diagrams, the shape of the distribution of main sequence stars is compatible with the presence of an obscuring cloud in the foreground at about 850+/-25 pc from the Sun. △ Less

Submitted 30 August, 2012; originally announced August 2012.

Comments: Accepted for publication ApJS 2012

A Bayesian Analysis of the Correlations Among Sunspot Cycles

Authors: Yaming Yu, David A. van Dyk, Vinay L. Kashyap, C. Alex Young

Abstract: Sunspot numbers form a comprehensive, long-duration proxy of solar activity and have been used numerous times to empirically investigate the properties of the solar cycle. A number of correlations have been discovered over the 24 cycles for which observational records are available. Here we carry out a sophisticated statistical analysis of the sunspot record that reaffirms these correlations, and… ▽ More Sunspot numbers form a comprehensive, long-duration proxy of solar activity and have been used numerous times to empirically investigate the properties of the solar cycle. A number of correlations have been discovered over the 24 cycles for which observational records are available. Here we carry out a sophisticated statistical analysis of the sunspot record that reaffirms these correlations, and sets up an empirical predictive framework for future cycles. An advantage of our approach is that it allows for rigorous assessment of both the statistical significance of various cycle features and the uncertainty associated with predictions. We summarize the data into three sequential relations that estimate the amplitude, duration, and time of rise to maximum for any cycle, given the values from the previous cycle. We find that there is no indication of a persistence in predictive power beyond one cycle, and conclude that the dynamo does not retain memory beyond one cycle. Based on sunspot records up to October 2011, we obtain, for Cycle 24, an estimated maximum smoothed monthly sunspot number of 97 +- 15, to occur in January--February 2014 +- 6 months. △ Less

Submitted 8 August, 2012; originally announced August 2012.

Comments: Accepted for publication in Solar Physics

The Role of White Dwarfs in Cataclysmic Variable Spin-down

Authors: O. Cohen, J. J. Drake, V. L. Kashyap

Abstract: We study the effect of a white dwarf on the spin-down of a cataclysmic variable system using a three-dimensional magnetohydrodynamic numerical model. The model includes the stellar corona, the stellar wind, and the WD mass and magnetic field. The existence of the WD modifies the system spin-down by physically blocking the stellar wind, restructuring the wind, channeling the wind towards the WD sur… ▽ More We study the effect of a white dwarf on the spin-down of a cataclysmic variable system using a three-dimensional magnetohydrodynamic numerical model. The model includes the stellar corona, the stellar wind, and the WD mass and magnetic field. The existence of the WD modifies the system spin-down by physically blocking the stellar wind, restructuring the wind, channeling the wind towards the WD surface, and by modifying the shape and size of the Alfvén surface. The combination of these processes differs among a set of simple test cases, and the resulting angular momentum loss rates vary by factors of 2-3, and by factors of two relative to a test model with a single M dwarf. While the model employs some simplifications, the results suggest angular momentum loss schemes currently employed in cataclysmic variable studies do not require drastic revision. Insights are also gained on wind accretion. We find that efficient accretion switches on quite rapidly with decreasing orbital separation. Accretion rates depend on magnetic field alignment and should be modulated by magnetic cycles on the M dwarf. For particular values of white dwarf magnetic field strength, an efficient syphoning of coronal plasma from the inward facing M dwarf hemisphere occurs. Wind accretion rates are expected to vary by factors of 10 or more between fairly similar close binaries, depending on magnetic field strengths and orbital separation. △ Less

Submitted 4 January, 2012; originally announced January 2012.

Comments: 7 pages, 3 figures, 1 table, accepted to ApJ Letters

Starspots on the fastest rotators in the Beta Pic moving group

Authors: D. Garcia-Alvarez, A. F. Lanza, S. Messina, J. J. Drake, F. van Wyk, R. R. Shobbrook, C. J. Butler, D. Kilkenny, J. G. Doyle, V. L. Kashyap

Abstract: Aims: We carried out high-resolution spectroscopy and BV(I)_C photometric monitoring of the two fastest late-type rotators in the nearby Beta Pictoris moving group, HD199143 (F7V) and CD-641208 (K7V). The motivation for this work is to investigate the rotation periods and photospheric spot patterns of these very young stars, with a longer term view to probing the evolution of rotation and magnetic… ▽ More Aims: We carried out high-resolution spectroscopy and BV(I)_C photometric monitoring of the two fastest late-type rotators in the nearby Beta Pictoris moving group, HD199143 (F7V) and CD-641208 (K7V). The motivation for this work is to investigate the rotation periods and photospheric spot patterns of these very young stars, with a longer term view to probing the evolution of rotation and magnetic activity during the early phases of main-sequence evolution. We also aim to derive information on key physical parameters, such as rotational velocity and rotation period. Methods: We applied maximum entropy (ME) and Tikhonov regularizing (TR) criteria to derive the surface spot map distributions of the optical modulation observed in HD199143 (F7 V) and CD-641208 (K7V). We also used cross-correlation techniques to determine stellar parameters such as radial velocities and rotational velocities. Lomb-Scargle periodograms were used to obtain the rotational periods from differential magnitude time series. Results: We find periods and inclinations of 0.356 days and 21.5deg for HD199143, and 0.355 days and 50.1deg for CD-641208. The spot maps of HD199143 obtained from the ME and TR methods are very similar, although the latter gives a smoother distribution of the filling factor. Maps obtained at two different epochs three weeks apart show a remarkable increase in spot coverage amounting to ~7% of the surface of the photosphere over a time period of only ~20 days. The spot maps of CD-641208 from the two methods show good longitudinal agreement, whereas the latitude range of the spots is extended to cover the whole visible hemisphere in the TR map. The distributions obtained from the first light curve of HD199143 show the presence of an extended and asymmetric active longitude with the maximum filling factor at longitude ~325degree. △ Less

Submitted 11 August, 2011; v1 submitted 28 July, 2011; originally announced July 2011.

Comments: Accepted by A&A. 13 pages, 13 figures (4 online included), 5 Tables

Statistical Characterization of the Chandra Source Catalog

Authors: Francis A. Primini, John C. Houck, John E. Davis, Michael A. Nowak, Ian N. Evans, Kenny J. Glotfelty, Craig S. Anderson, Nina R. Bonaventura, Judy C. Chen, Stephen M. Doe, Janet D. Evans, Giuseppina Fabbiano, Elizabeth C. Galle, Danny G. Gibbs II, John D. Grier, Roger M. Hain, Diane M. Hall, Peter N. Harbo, Xiangqun, He, Margarita Karovska, Vinay L. Kashyap, Jennifer Lauer, Michael L. McCollough, Jonathan C. McDowell , et al. (14 additional authors not shown)

Abstract: The first release of the Chandra Source Catalog (CSC) contains ~95,000 X-ray sources in a total area of ~0.75% of the entire sky, using data from ~3,900 separate ACIS observations of a multitude of different types of X-ray sources. In order to maximize the scientific benefit of such a large, heterogeneous data-set, careful characterization of the statistical properties of the catalog, i.e., comple… ▽ More The first release of the Chandra Source Catalog (CSC) contains ~95,000 X-ray sources in a total area of ~0.75% of the entire sky, using data from ~3,900 separate ACIS observations of a multitude of different types of X-ray sources. In order to maximize the scientific benefit of such a large, heterogeneous data-set, careful characterization of the statistical properties of the catalog, i.e., completeness, sensitivity, false source rate, and accuracy of source properties, is required. Characterization efforts of other, large Chandra catalogs, such as the ChaMP Point Source Catalog (Kim et al. 2007) or the 2 Mega-second Deep Field Surveys (Alexander et al. 2003), while informative, cannot serve this purpose, since the CSC analysis procedures are significantly different and the range of allowable data is much less restrictive. We describe here the characterization process for the CSC. This process includes both a comparison of real CSC results with those of other, deeper Chandra catalogs of the same targets and extensive simulations of blank-sky and point source populations. △ Less

Submitted 5 May, 2011; v1 submitted 3 May, 2011; originally announced May 2011.

Comments: To be published in the Astrophysical Journal Supplement Series (Fig. 52 replaced with a version which astro-ph can convert to PDF without issues.)

Accounting for Calibration Uncertainties in X-ray Analysis: Effective Areas in Spectral Fitting

Authors: Hyunsook Lee, Vinay L. Kashyap, David A. van Dyk, Alanna Connors, Jeremy J. Drake, Rima Izem, Xiao-Li Meng, Shandong Min, Taeyoung Park, Pete Ratzlaff, Aneta Siemiginowska, Andreas Zezas

Abstract: While considerable advance has been made to account for statistical uncertainties in astronomical analyses, systematic instrumental uncertainties have been generally ignored. This can be crucial to a proper interpretation of analysis results because instrumental calibration uncertainty is a form of systematic uncertainty. Ignoring it can underestimate error bars and introduce bias into the fitted… ▽ More While considerable advance has been made to account for statistical uncertainties in astronomical analyses, systematic instrumental uncertainties have been generally ignored. This can be crucial to a proper interpretation of analysis results because instrumental calibration uncertainty is a form of systematic uncertainty. Ignoring it can underestimate error bars and introduce bias into the fitted values of model parameters. Accounting for such uncertainties currently requires extensive case-specific simulations if using existing analysis packages. Here we present general statistical methods that incorporate calibration uncertainties into spectral analysis of high-energy data. We first present a method based on multiple imputation that can be applied with any fitting method, but is necessarily approximate. We then describe a more exact Bayesian approach that works in conjunction with a Markov chain Monte Carlo based fitting. We explore methods for improving computational efficiency, and in particular detail a method of summarizing calibration uncertainties with a principal component analysis of samples of plausible calibration files. This method is implemented using recently codified Chandra effective area uncertainties for low-resolution spectral analysis and is verified using both simulated and actual Chandra data. Our procedure for incorporating effective area uncertainty is easily generalized to other types of calibration uncertainties. △ Less

Submitted 22 February, 2011; originally announced February 2011.

Comments: 61 pages double spaced, 8 figures, accepted for publication in ApJ

Journal ref: The Astrophysical Journal (2011), 731, 126

The Dynamics of Stellar Coronae Harboring Hot-jupiters II. A Space Weather Event on A Hot-jupiter

Authors: O. Cohen, V. L. Kashyap, J. J. Drake, I. V. Sokolov, T. I. Gombosi

Abstract: We carry out a numerical simulation depicting the effects of a Coronal Mass Ejection (CME) event on a close-in giant planet in an extrasolar system. We drive the CME in a similar manner as in simulations of space weather events on Earth. The simulation includes the planetary orbital motion, which leads to the forming of a comet-like planetary magnetotail which is oriented almost perpendicular to t… ▽ More We carry out a numerical simulation depicting the effects of a Coronal Mass Ejection (CME) event on a close-in giant planet in an extrasolar system. We drive the CME in a similar manner as in simulations of space weather events on Earth. The simulation includes the planetary orbital motion, which leads to the forming of a comet-like planetary magnetotail which is oriented almost perpendicular to the direction of propagation of the CME. The combination of this feature and the fact that the CME does not expand much by the time it reaches the planet leads to a unique CME-magnetosphere interaction, where the CME itself is highly affected by the presence of the planetary magnetosphere. We find that the planet is well-shielded from CME penetration, even for a relatively weak internal magnetic field. The planetary angular momentum loss associated with such an event is negligible compared to the total planetary angular momentum. We also find that the energy which is deposited in the magnetosphere is much higher than in the case of the Earth, and our simulation suggests there is a large-scale change in the orientation of the magnetosphere-ionosphere current system during the CME event. △ Less

Submitted 7 July, 2011; v1 submitted 20 February, 2011; originally announced February 2011.

Comments: 16 pages, 10 figures, accepted to ApJ

The Dynamics of Stellar Coronae Harboring Hot-jupiters I. A Time-dependent MHD Simulation of the Interplanetary Environment in the HD 189733 Planetary System

Authors: O. Cohen, V. L. Kashyap, J. J. Drake, I. V. Sokolov, C. Garraffo, T. I. Gombosi

Abstract: We carry out the first time-dependent numerical MagnetoHydroDynamic modeling of an extrasolar planetary system to study the interaction of the stellar magnetic field and wind with the planetary magnetosphere and outflow. We base our model on the parameters of the HD 189733 system, which harbors a close-in giant planet. Our simulation reveals a highly structured stellar corona characterized by sect… ▽ More We carry out the first time-dependent numerical MagnetoHydroDynamic modeling of an extrasolar planetary system to study the interaction of the stellar magnetic field and wind with the planetary magnetosphere and outflow. We base our model on the parameters of the HD 189733 system, which harbors a close-in giant planet. Our simulation reveals a highly structured stellar corona characterized by sectors with different plasma properties. The star-planet interaction varies in magnitude and complexity, depending on the planetary phase, planetary magnetic field strength, and the relative orientation of the stellar and planetary fields. It also reveals a long, comet-like tail which is a result of the wrapping of the planetary magnetospheric tail by its fast orbital motion. A reconnection event occurs at a specific orbital phase, causing mass loss from the planetary magnetosphere that can generate a hot spot on the stellar surface. The simulation also shows that the system has sufficient energy to produce hot-spots observed in Ca II lines in giant planet hosting stars. However, the short duration of the reconnection event suggests that such SPI cannot be observed persistently. △ Less

Submitted 25 January, 2011; originally announced January 2011.

Comments: Submitted to ApJ, 14 pages, 9 figures

The Impact of Hot Jupiters on the Spin-down of their Host Stars

Authors: O. Cohen, J. J Drake, V. L. Kashyap, I. V. Sokolov, T. I. Gombosi

Abstract: We present a numerical Magnetohydrodynamic (MHD) study of the dependence of stellar mass and angular momentum- loss rates on the orbital distance to close-in giant planets. We find that the mass loss rate drops by a factor of $\approx$1.5-2, while the angular momentum loss rate drops by a factor of $\approx$4 as the distance decreases past the Alfvén surface. This reduction in angular momentum los… ▽ More We present a numerical Magnetohydrodynamic (MHD) study of the dependence of stellar mass and angular momentum- loss rates on the orbital distance to close-in giant planets. We find that the mass loss rate drops by a factor of $\approx$1.5-2, while the angular momentum loss rate drops by a factor of $\approx$4 as the distance decreases past the Alfvén surface. This reduction in angular momentum loss is due to the interaction between the stellar and planetary Alfvén surfaces, which modifies the global structure of the stellar corona and stellar wind on the hemisphere facing the planet, as well as the opposite hemisphere. The simulation also shows that the magnitude of change in angular momentum loss rate depends mostly on the strength of the planetary magnetic field and not on its polarity. The interaction however, begins at greater separation if the overall field topology of the star and the planet are of anti-aligned. Our results are consistent with evidence for excess angular momentum in stars harboring close-in giant planets, and show that the reduction in wind-driven angular momentum loss can compete with, and perhaps dominate, spin-up due to tidal interaction. △ Less

Submitted 29 September, 2010; originally announced September 2010.

Comments: 6 pages, 3 figures

Close Binaries with Infrared Excess: Destroyers of Worlds?

Authors: M. Matranga, J. J. Drake, V. L. Kashyap, M. Marengo, M. J. Kuchner

Abstract: We present the results of a Spitzer photometric investigation into the IR excesses of close binary systems. In a sample of 10 objects, excesses in IRAC and MIPS24 bands implying the presence of warm dust are found for 3. For 2 objects we do not find excesses reported in earlier IRAS studies. We discuss the results in the context of the scenario suggested by Rhee and co-workers, in which warm dust… ▽ More We present the results of a Spitzer photometric investigation into the IR excesses of close binary systems. In a sample of 10 objects, excesses in IRAC and MIPS24 bands implying the presence of warm dust are found for 3. For 2 objects we do not find excesses reported in earlier IRAS studies. We discuss the results in the context of the scenario suggested by Rhee and co-workers, in which warm dust is continuously created by destructive collisions between planetary bodies. A simple numerical model for the steady-state distribution of dust in one IR excess system shows a central clearing of radius 0.22 AU caused by dynamical perturbations from the binary star. This is consistent with the size of the central clearing derived from the Spitzer spectral energy distribution. We conclude that close binaries could be efficient "destroyers of worlds", lead to destabilize the orbits of their planetary progeny by magnetically-driven angular momentum loss and secular shrinkage of the binary separation. △ Less

Submitted 28 July, 2010; originally announced July 2010.

Comments: Accepted by The Astrophysical Journal

The Coronal Structure of AB Doradus

Authors: O. Cohen, J. J. Drake, V. L. Kashyap, G. A. J. Hussain, T. I. Gombosi

Abstract: We perform a numerical simulation of the corona of the young, rapidly rotating K0 dwarf AB Doradus using a global MHD model. The model is driven by a surface map of the radial magnetic field constructed using Zeeman-Doppler Imaging. We find that the global structure of the stellar corona is dominated by strong azimuthal tangling of the magnetic field due to the rapid rotation. The MHD solution ena… ▽ More We perform a numerical simulation of the corona of the young, rapidly rotating K0 dwarf AB Doradus using a global MHD model. The model is driven by a surface map of the radial magnetic field constructed using Zeeman-Doppler Imaging. We find that the global structure of the stellar corona is dominated by strong azimuthal tangling of the magnetic field due to the rapid rotation. The MHD solution enables us to calculate realistic Alfvén surfaces and we can therefore estimate the stellar mass loss rate and angular momentum loss rate without making undue theoretical simplifications. We consider three cases, parametrized by the base density of the corona, that span the range of possible solutions for the system. We find that overall, the mass and angular-momentum loss rates are higher than in the solar case; the mass loss rates are 10 to 500 times higher, and the angular momentum loss rate can be up to $3\times{10}^4$ higher than present day solar values. Our simulations show that this model can be use to constrain the wide parameter space of stellar systems. It also shows that an MHD approach can provide more information about the physical system over the commonly used potential field extrapolation. △ Less

Submitted 25 July, 2010; originally announced July 2010.

Comments: 13 pages, 7 figures

On Computing Upper Limits to Source Intensities

Authors: Vinay L. Kashyap, David A. van Dyk, Alanna Connors, Peter Freeman, Aneta Siemiginowska, Jin Xu, Andreas Zezas

Abstract: A common problem in astrophysics is determining how bright a source could be and still not be detected. Despite the simplicity with which the problem can be stated, the solution involves complex statistical issues that require careful analysis. In contrast to the confidence bound, this concept has never been formally analyzed, leading to a great variety of often ad hoc solutions. Here we formulate… ▽ More A common problem in astrophysics is determining how bright a source could be and still not be detected. Despite the simplicity with which the problem can be stated, the solution involves complex statistical issues that require careful analysis. In contrast to the confidence bound, this concept has never been formally analyzed, leading to a great variety of often ad hoc solutions. Here we formulate and describe the problem in a self-consistent manner. Detection significance is usually defined by the acceptable proportion of false positives (the TypeI error), and we invoke the complementary concept of false negatives (the TypeII error), based on the statistical power of a test, to compute an upper limit to the detectable source intensity. To determine the minimum intensity that a source must have for it to be detected, we first define a detection threshold, and then compute the probabilities of detecting sources of various intensities at the given threshold. The intensity that corresponds to the specified TypeII error probability defines that minimum intensity, and is identified as the upper limit. Thus, an upper limit is a characteristic of the detection procedure rather than the strength of any particular source and should not be confused with confidence intervals or other estimates of source intensity. This is particularly important given the large number of catalogs that are being generated from increasingly sensitive surveys. We discuss the differences between these upper limits and confidence bounds. Both measures are useful quantities that should be reported in order to extract the most science from catalogs, though they answer different statistical questions: an upper bound describes an inference range on the source intensity, while an upper limit calibrates the detection process. We provide a recipe for computing upper limits that applies to all detection algorithms. △ Less

Submitted 22 June, 2010; originally announced June 2010.

Comments: 30 pages, 12 figures, accepted in ApJ

Journal ref: The Astrophysical Journal (2010), 719, 900