Showing 1–41 of 41 results for author: Parra, M

A second view on the X-ray polarization of NGC 4151 with IXPE

Authors: V. E. Gianolli, S. Bianchi, E. Kammoun, A. Gnarini, A. Marinucci, F. Ursini, M. Parra, A. Tortosa, A. De Rosa, D. E. Kim, F. Marin, G. Matt, R. Serafinelli, P. Soffitta, D. Tagliacozzo, L. Di Gesu, C. Done, H. L. Marshall, R. Middei, R. Mikusincova, P-O. Petrucci, S. Ravi, J. Svoboda, F. Tombesi

Abstract: We report on the second observing program of the active galactic nucleus NGC 4151 with simultaneous Imaging X-ray Polarimetry Explorer (IXPE; {\sim}750 ks), NuSTAR ({\sim}60 ks), XMM-Newton ({\sim}75 ks), and NICER ({\sim}65 ks) pointings. NGC 4151 is the first Type 1 radio-quiet Seyfert galaxy with constrained polarization properties for the X-ray corona. Despite the lower flux state in which the… ▽ More We report on the second observing program of the active galactic nucleus NGC 4151 with simultaneous Imaging X-ray Polarimetry Explorer (IXPE; {\sim}750 ks), NuSTAR ({\sim}60 ks), XMM-Newton ({\sim}75 ks), and NICER ({\sim}65 ks) pointings. NGC 4151 is the first Type 1 radio-quiet Seyfert galaxy with constrained polarization properties for the X-ray corona. Despite the lower flux state in which the source is re-observed and the resulting higher contribution of the constant reflection component in the IXPE energy band, our results are in agreement with the first detection. From the polarimetric analysis, a polarization degree Π = 4.7 {\pm} 1.3 per cent and angle Ψ = 77° {\pm} 8° east of north (68 per cent confidence level) are derived in the 2.0 - 8.0 keV energy range. Combining the two observations leads to polarization properties that are more constrained than those of the individual detections, showing Π = 4.5 {\pm} 0.9 per cent and Ψ = 81° {\pm} 6° (with detection significance {\sim}4.6σ). The observed polarization angle aligns very well with the radio emission in this source, supporting, together with the significant polarization degree, a slab or wedge geometry for the X-ray corona. However, a switch in the polarization angle at low energies (37° {\pm} 7° in the 2 - 3.5 keV bin) suggests the presence of another component. When it is included in the spectro-polarimetric fit, a high polarization degree disfavors an interpretation in terms of a leakage through the absorbers, rather pointing to scattering from some kind of mirror. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: 10 pages, 8 figures, 3 tables. Submitted to Astronomy & Astrophysics

Graph Permutation Entropy: Extensions to the Continuous Case, A step towards Ordinal Deep Learning, and More

Authors: Om Roy, Avalon Campbell-Cousins, John Stewart Fabila Carrasco, Mario A Parra, Javier Escudero

Abstract: Nonlinear dynamics play an important role in the analysis of signals. A popular, readily interpretable nonlinear measure is Permutation Entropy. It has recently been extended for the analysis of graph signals, thus providing a framework for non-linear analysis of data sampled on irregular domains. Here, we introduce a continuous version of Permutation Entropy, extend it to the graph domain, and de… ▽ More Nonlinear dynamics play an important role in the analysis of signals. A popular, readily interpretable nonlinear measure is Permutation Entropy. It has recently been extended for the analysis of graph signals, thus providing a framework for non-linear analysis of data sampled on irregular domains. Here, we introduce a continuous version of Permutation Entropy, extend it to the graph domain, and develop a ordinal activation function akin to the one of neural networks. This is a step towards Ordinal Deep Learning, a potentially effective and very recently posited concept. We also formally extend ordinal contrasts to the graph domain. Continuous versions of ordinal contrasts of length 3 are also introduced and their advantage is shown in experiments. We also integrate specific contrasts for the analysis of images and show that it generalizes well to the graph domain allowing a representation of images, represented as graph signals, in a plane similar to the entropy-complexity one. Applications to synthetic data, including fractal patterns and popular non-linear maps, and real-life MRI data show the validity of these novel extensions and potential benefits over the state of the art. By extending very recent concepts related to permutation entropy to the graph domain, we expect to accelerate the development of more graph-based entropy methods to enable nonlinear analysis of a broader kind of data and establishing relationships with emerging ideas in data science. △ Less

Submitted 11 July, 2024; v1 submitted 10 July, 2024; originally announced July 2024.

An IXPE-Led X-ray Spectro-Polarimetric Campaign on the Soft State of Cygnus X-1: X-ray Polarimetric Evidence for Strong Gravitational Lensing

Authors: James F. Steiner, Edward Nathan, Kun Hu, Henric Krawczynski, Michal Dovciak, Alexandra Veledina, Fabio Muleri, Jiri Svoboda, Kevin Alabarta, Maxime Parra, Yash Bhargava, Giorgio Matt, Juri Poutanen, Pierre-Olivier Petrucci, Allyn F. Tennant, M. Cristina Baglio, Luca Baldini, Samuel Barnier, Sudip Bhattacharyya, Stefano Bianchi, Maimouna Brigitte, Mauricio Cabezas, Floriane Cangemi, Fiamma Capitanio, Jacob Casey , et al. (112 additional authors not shown)

Abstract: We present the first X-ray spectropolarimetric results for Cygnus X-1 in its soft state from a campaign of five IXPE observations conducted during 2023 May-June. Companion multiwavelength data during the campaign are likewise shown. The 2-8 keV X-rays exhibit a net polarization degree PD=1.99%+/-0.13% (68% confidence). The polarization signal is found to increase with energy across IXPE's 2-8 keV… ▽ More We present the first X-ray spectropolarimetric results for Cygnus X-1 in its soft state from a campaign of five IXPE observations conducted during 2023 May-June. Companion multiwavelength data during the campaign are likewise shown. The 2-8 keV X-rays exhibit a net polarization degree PD=1.99%+/-0.13% (68% confidence). The polarization signal is found to increase with energy across IXPE's 2-8 keV bandpass. The polarized X-rays exhibit an energy-independent polarization angle of PA=-25.7+/-1.8 deg. East of North (68% confidence). This is consistent with being aligned to Cyg X-1's AU-scale compact radio jet and its pc-scale radio lobes. In comparison to earlier hard-state observations, the soft state exhibits a factor of 2 lower polarization degree, but a similar trend with energy and a similar (also energy-independent) position angle. When scaling by the natural unit of the disk temperature, we find the appearance of a consistent trendline in the polarization degree between soft and hard states. Our favored polarimetric model indicates Cyg X-1's spin is likely high (a* above ~0.96). The substantial X-ray polarization in Cyg X-1's soft state is most readily explained as resulting from a large portion of X-rays emitted from the disk returning and reflecting off the disk surface, generating a high polarization degree and a polarization direction parallel to the black hole spin axis and radio jet. In IXPE's bandpass, the polarization signal is dominated by the returning reflection emission. This constitutes polarimetric evidence for strong gravitational lensing of X-rays close to the black hole. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 20 pages, accepted for publication in ApJL

Impact of the disk magnetization on MHD disk wind signature

Authors: Sudeb Ranjan Datta, Susmita Chakravorty, Jonathan Ferreira, Pierre-Olivier Petrucci, Timothy R Kallman, Jonatan Jacquemin-Ide, Nathan Zimniak, Joern Wilms, Stefano Bianchi, Maxime Parra, Maïca Clavel

Abstract: Observation of blue-shifted X-ray absorption lines indicates the presence of wind from the accretion disk in X-ray binaries. Magnetohydrodynamic (MHD) driving is one of the possible wind launching mechanisms. Recent theoretical development makes magnetic accretion-ejection self-similar solutions much more generalized, and wind can be launched even at much lower magnetization compared to equipartit… ▽ More Observation of blue-shifted X-ray absorption lines indicates the presence of wind from the accretion disk in X-ray binaries. Magnetohydrodynamic (MHD) driving is one of the possible wind launching mechanisms. Recent theoretical development makes magnetic accretion-ejection self-similar solutions much more generalized, and wind can be launched even at much lower magnetization compared to equipartition value, which was the only possibility beforehand. Here, we model the transmitted spectra through MHD driven photoionized wind - models which have different values of magnetizations. We investigate the possibility of detecting absorption lines by the upcoming instruments XRISM and Athena. Attempts are made to find the robustness of the method of fitting asymmetric line profiles by multiple Gaussians. We use photoionization code XSTAR to simulate the transmitted model spectra. Fake observed spectra are finally produced by convolving model spectra with instruments' responses. Since the line asymmetries are apparent in the convolved spectra as well, this can be used as an observable diagnostic to fit for, in future XRISM and Athena spectra. We demonstrate some amount of rigor in assessing the equivalent widths of the major absorption lines, including the Fe XXVI Ly$α$ doublets which can be clearly distinguished in the superior quality, future high resolution spectra. Disk magnetization becomes another crucial MHD variable that can significantly alter the absorption line profiles. Low magnetization pure MHD outflow models are dense enough to be observed by the existing or upcoming instruments. Thus these models become simpler alternatives to MHD-thermal models. Fitting with multiple Gaussians is a promising method to handle asymmetric line profiles, as well as the Fe XXVI Ly$α$ doublets. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 17 pages (including 4 pages of Appendix), 13 figures, Accepted for publication in A&A

Journal ref: A&A 687, A2 (2024)

Supermassive Black Hole Winds in X-rays -- SUBWAYS. III. A population study on ultra-fast outflows

Authors: V. E. Gianolli, S. Bianchi, P-O Petrucci, M. Brusa, G. Chartas, G. Lanzuisi, G. A. Matzeu, M. Parra, F. Ursini, E. Behar, M. Bischetti, A. Comastri, E. Costantini, G. Cresci, M. Dadina, B. De Marco, A. De Rosa, F. Fiore, M. Gaspari, R. Gilli, M. Giustini, M. Guainazzi, A. R. King, S. Kraemer, G. Kriss , et al. (22 additional authors not shown)

Abstract: The detection of blue-shifted absorption lines likely associated with ionized Iron K-shell transitions in the X-ray spectra of many Active Galactic Nuclei (AGN) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c-0.6c), named Ultra-Fast Outflow (UFO). Within the SUBWAYS project we characterized these winds starting from a sample of 22 radio-quiet qua… ▽ More The detection of blue-shifted absorption lines likely associated with ionized Iron K-shell transitions in the X-ray spectra of many Active Galactic Nuclei (AGN) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c-0.6c), named Ultra-Fast Outflow (UFO). Within the SUBWAYS project we characterized these winds starting from a sample of 22 radio-quiet quasars at 0.1 < z < 0.4, and compared the results with similar studies in the literature on samples of 42 local radio-quiet Seyfert galaxies and 14 high redshift radio-quiet quasars. The scope of our work is a statistical study of UFO parameters and incidence, considering key physical properties of the sources, e.g. supermassive black hole (SMBH) mass, bolometric luminosity, accretion rates and Spectral Energy Distribution, with the aim of gaining new insights into the UFO launching mechanisms. We find indications that highly luminous AGN with steeper X-ray/UV ratio, are more likely to host UFO. The presence of UFO is not significantly related to any other AGN property in our sample. These findings suggest that the UFO phenomenon may be transient. Focusing on AGN with UFO, other important results are: (1) faster UFO have larger ionization parameters and column densities; (2) X-ray radiation plays a more crucial role in driving highly ionized winds compared to UV; (3) the correlation between outflow velocity and luminosity is significantly flatter than what expected for radiatively driven winds; (4) more massive BH experience higher wind mass-losses, suppressing accretion of matter onto the BH; (5) the UFO launching radius is positively correlated with the Eddington ratio. Furthermore, our analysis suggest the involvement of multiple launching mechanisms, including radiation pressure and magneto-hydrodynamic processes, rather than pointing to a single, universally applicable mechanism. △ Less

Submitted 11 April, 2024; v1 submitted 14 March, 2024; originally announced March 2024.

Comments: 52 pages, 154 figures and 7 tables. Accepted in Astronomy & Astrophysics, in press

Journal ref: A&A 687, A235 (2024)

Classical dynamical $r$-matrices for the Chern-Simons formulation of generalised 3d gravity

Authors: Juan Carlos Morales Parra, Bernd Schroers

Abstract: Classical dynamical $r$-matrices arise naturally in the combinatorial description of the phase space of Chern-Simons theories, either through the inclusion of dynamical sources or through a gauge-fixing procedure involving two punctures. Here we consider classical dynamical $r$-matrices for the family of Lie algebras which arise in the Chern-Simons formulation of 3d gravity, for any value of the c… ▽ More Classical dynamical $r$-matrices arise naturally in the combinatorial description of the phase space of Chern-Simons theories, either through the inclusion of dynamical sources or through a gauge-fixing procedure involving two punctures. Here we consider classical dynamical $r$-matrices for the family of Lie algebras which arise in the Chern-Simons formulation of 3d gravity, for any value of the cosmological constant. We derive differential equations for classical dynamical $r$-matrices in this case, and show that they can be viewed as generalised complexifications, in a sense which we define, of the equations governing dynamical $r$-matrices for $\mathfrak{su}(2)$ and $\mathfrak{sl}(2,\mathbb{R})$. We obtain explicit families of solutions and relate them, via Weierstrass factorisation, to solutions found by Feher, Gabor, Marshall, Palla and Pusztai in the context of chiral WZWN models. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 46 pages, no figures

FAST functional connectivity implicates P300 connectivity in working memory deficits in Alzheimer's disease

Authors: Om Roy, Yashar Moshfeghi, Agustin Ibanez, Francisco Lopera, Mario A Parra, Keith M Smith

Abstract: Measuring transient functional connectivity is an important challenge in Electroencephalogram (EEG) research. Here, the rich potential for insightful, discriminative information of brain activity offered by high temporal resolution is confounded by the inherent noise of the medium and the spurious nature of correlations computed over short temporal windows. We propose a novel methodology to overco… ▽ More Measuring transient functional connectivity is an important challenge in Electroencephalogram (EEG) research. Here, the rich potential for insightful, discriminative information of brain activity offered by high temporal resolution is confounded by the inherent noise of the medium and the spurious nature of correlations computed over short temporal windows. We propose a novel methodology to overcome these problems called Filter Average Short-Term (FAST) functional connectivity. First, long-term, stable, functional connectivity is averaged across an entire study cohort for a given pair of Visual Short Term Memory (VSTM) tasks. The resulting average connectivity matrix, containing information on the strongest general connections for the tasks, is used as a filter to analyse the transient high temporal resolution functional connectivity of individual subjects. In simulations, we show that this method accurately discriminates differences in noisy Event-Related Potentials (ERPs) between two conditions where standard connectivity and other comparable methods fail. We then apply this to analyse activity related to visual short-term memory binding deficits in two cohorts of familial and sporadic Alzheimer's disease. Reproducible significant differences were found in the binding task with no significant difference in the shape task in the P300 ERP range. This allows new sensitive measurements of transient functional connectivity, which can be implemented to obtain results of clinical significance. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Tracking the X-ray Polarization of the Black Hole Transient Swift J1727.8-1613 during a State Transition

Authors: Adam Ingram, Niek Bollemeijer, Alexandra Veledina, Michal Dovciak, Juri Poutanen, Elise Egron, Thomas D. Russell, Sergei A. Trushkin, Michela Negro, Ajay Ratheesh, Fiamma Capitanio, Riley Connors, Joseph Neilsen, Alexander Kraus, Maria Noemi Iacolina, Alberto Pellizzoni, Maura Pilia, Francesco Carotenuto, Giorgio Matt, Guglielmo Mastroserio, Philip Kaaret, Stefano Bianchi, Javier A. Garcia, Matteo Bachetti, Kinwah Wu , et al. (98 additional authors not shown)

Abstract: We report on an observational campaign on the bright black hole X-ray binary Swift J1727.8$-$1613 centered around five observations by the Imaging X-ray Polarimetry Explorer (IXPE). These observations track for the first time the evolution of the X-ray polarization of a black hole X-ray binary across a hard to soft state transition. The 2--8 keV polarization degree decreased from $\sim$4\% to… ▽ More We report on an observational campaign on the bright black hole X-ray binary Swift J1727.8$-$1613 centered around five observations by the Imaging X-ray Polarimetry Explorer (IXPE). These observations track for the first time the evolution of the X-ray polarization of a black hole X-ray binary across a hard to soft state transition. The 2--8 keV polarization degree decreased from $\sim$4\% to $\sim$3\% across the five observations, but the polarization angle remained oriented in the North-South direction throughout. Based on observations with the Australia Telescope Compact Array (ATCA), we find that the intrinsic 7.25 GHz radio polarization aligns with the X-ray polarization. Assuming the radio polarization aligns with the jet direction (which can be tested in the future with higher spatial resolution images of the jet), our results imply that the X-ray corona is extended in the disk plane, rather than along the jet axis, for the entire hard intermediate state. This in turn implies that the long ($\gtrsim$10 ms) soft lags that we measure with the Neutron star Interior Composition ExploreR (NICER) are dominated by processes other than pure light-crossing delays. Moreover, we find that the evolution of the soft lag amplitude with spectral state does not follow the trend seen for other sources, implying that Swift J1727.8$-$1613 is a member of a hitherto under-sampled sub-population. △ Less

Submitted 24 April, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: Accepted for publication in ApJ. 21 pages, 8 figures

Entropy of Compositions of Covering Correspondences

Authors: Vanessa Matus De La Parra

Abstract: We use the measure theoretical entropy from Vivas-Sirvent and the topological entropies from Dinh-Sibony and Kelly-Tennant to prove that the equidistribution measures $μ_-$ and $μ_+$ found by the author maximize entropy for the $1$-parameter family $\lbrace\mathcal{F}_a\rbrace_a$ studied by Bullett-Lomonaco. We find as well a measure of maximal entropy for the composition of two deleted covering c… ▽ More We use the measure theoretical entropy from Vivas-Sirvent and the topological entropies from Dinh-Sibony and Kelly-Tennant to prove that the equidistribution measures $μ_-$ and $μ_+$ found by the author maximize entropy for the $1$-parameter family $\lbrace\mathcal{F}_a\rbrace_a$ studied by Bullett-Lomonaco. We find as well a measure of maximal entropy for the composition of two deleted covering correspondences studied by Bullett. △ Less

Submitted 22 October, 2023; originally announced October 2023.

IXPE observation confirms a high spin in the accreting black hole 4U 1957+115

Authors: L. Marra, M. Brigitte, N. Rodriguez Cavero, S. Chun, J. F. Steiner, M. Dovčiak, M. Nowak, S. Bianchi, F. Capitanio, A. Ingram, G. Matt, F. Muleri, J. Podgorný, J. Poutanen, J. Svoboda, R. Taverna, F. Ursini, A. Veledina, A. De Rosa, J. A. Garcia, A. A. Lutovinov, I. A. Mereminskiy, R. Farinelli, S. Gunji, P. Kaaret , et al. (91 additional authors not shown)

Abstract: We present the results of the first X-ray polarimetric observation of the low-mass X-ray binary 4U 1957+115, performed with the Imaging X-ray Polarimetry Explorer in May 2023. The binary system has been in a high-soft spectral state since its discovery and is thought to host a black hole. The $\sim$571 ks observation reveals a linear polarisation degree of $1.9\% \pm 0.6\%$ and a polarisation angl… ▽ More We present the results of the first X-ray polarimetric observation of the low-mass X-ray binary 4U 1957+115, performed with the Imaging X-ray Polarimetry Explorer in May 2023. The binary system has been in a high-soft spectral state since its discovery and is thought to host a black hole. The $\sim$571 ks observation reveals a linear polarisation degree of $1.9\% \pm 0.6\%$ and a polarisation angle of $-41^\circ.8 \pm 7^\circ.9$ in the 2-8 keV energy range. Spectral modelling is consistent with the dominant contribution coming from the standard accretion disc, while polarimetric data suggest a significant role of returning radiation: photons that are bent by strong gravity effects and forced to return to the disc surface, where they can be reflected before eventually reaching the observer. In this setting, we find that models with a black hole spin lower than 0.96 and an inclination lower than $50^\circ$ are disfavoured. △ Less

Submitted 8 February, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 12 pages, 10 figures, 2 tables, accepted for publication in A&A

Discovery of X-ray Polarization from the Black Hole Transient Swift J1727.8-1613

Authors: Alexandra Veledina, Fabio Muleri, Michal Dovciak, Juri Poutanen, Ajay Ratheesh, Fiamma Capitanio, Giorgio Matt, Paolo Soffitta, Allyn F. Tennant, Michela Negro, Philip Kaaret, Enrico Costa, Adam Ingram, Jiri Svoboda, Henric Krawczynski, Stefano Bianchi, James F. Steiner, Javier A. Garcia, Vadim Kravtsov, Anagha P. Nitindala, Melissa Ewing, Guglielmo Mastroserio, Andrea Marinucci, Francesco Ursini, Francesco Tombesi , et al. (91 additional authors not shown)

Abstract: We report the first detection of the X-ray polarization of the bright transient Swift J1727.8-1613 with the Imaging X-ray Polarimetry Explorer. The observation was performed at the beginning of the 2023 discovery outburst, when the source resided in the bright hard state. We find a time- and energy-averaged polarization degree of 4.1%+/-0.2% and a polarization angle of 2.2+/-1.3 degrees (errors at… ▽ More We report the first detection of the X-ray polarization of the bright transient Swift J1727.8-1613 with the Imaging X-ray Polarimetry Explorer. The observation was performed at the beginning of the 2023 discovery outburst, when the source resided in the bright hard state. We find a time- and energy-averaged polarization degree of 4.1%+/-0.2% and a polarization angle of 2.2+/-1.3 degrees (errors at 68% confidence level; this translates to about 20-sigma significance of the polarization detection). This finding suggests that the hot corona emitting the bulk of the detected X-rays is elongated, rather than spherical. The X-ray polarization angle is consistent with that found in sub-mm wavelengths. Since the sub-mm polarization was found to be aligned with the jet direction in other X-ray binaries, this indicates that the corona is elongated orthogonal to the jet. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 10 pages, 3 figures, submitted

The current state of disk wind observations in BHLMXBs through X-ray absorption lines in the iron band

Authors: M. Parra, P. -O. Petrucci, S. Bianchi, V. E. Gianolli, F. Ursini, G. Ponti

Abstract: The presence of blueshifted absorption lines in the X-ray spectra of Black Hole Low Mass X-ray Binaries is the telltale of massive outflows called winds. These signatures are found almost exclusively in soft states of high-inclined systems, hinting at equatorial ejections originating from the accretion disk and deeply intertwined with the evolution of the outburst patterns displayed by these syste… ▽ More The presence of blueshifted absorption lines in the X-ray spectra of Black Hole Low Mass X-ray Binaries is the telltale of massive outflows called winds. These signatures are found almost exclusively in soft states of high-inclined systems, hinting at equatorial ejections originating from the accretion disk and deeply intertwined with the evolution of the outburst patterns displayed by these systems. In the wake of the launch of the new generation of X-ray spectrometers, studies of wind signatures remain mostly restricted to single sources and outbursts, with some of the recent detections departing from the commonly expected behaviors. We thus give an update to the current state of iron band absorption lines detections, through the analysis of all publicly available XMM-$Newton$-PN and $Chandra$-HETG exposures of known Black Hole Low-Mass X-ray Binary candidates. Our results agree with previous studies, with wind detections exclusively found in dipping, high-inclined sources, and almost exclusively in bright ($L_{X}>0.01L_{Edd}$) soft ($HR<0.8$) states, with blueshift values generally restricted to few 100 km s$^{-1}$. The line parameters indicate similar properties between objects and outbursts of single sources, and despite more than 20 years of data, very few sources have the HID sampling necessary to properly study the evolution of the wind during single outbursts. We provide an online tool with details of the wind signatures and outburst evolution data for all sources in the sample. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: 26 pages, 11 figures | Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 681, A49 (2024)

New Bayesian method for estimation of Value at Risk and Conditional Value at Risk

Authors: Jacinto Martín, M. Isabel Parra, Eva L. Sanjuán, Mario M. Pizarro

Abstract: Value at Risk (VaR) and Conditional Value at Risk (CVaR) have become the most popular measures of market risk in Financial and Insurance fields. However, the estimation of both risk measures is challenging, because it requires the knowledge of the tail of the distribution. Therefore, tools from Extreme Value Theory are usually employed, considering that the tail data follow a Generalized Pareto di… ▽ More Value at Risk (VaR) and Conditional Value at Risk (CVaR) have become the most popular measures of market risk in Financial and Insurance fields. However, the estimation of both risk measures is challenging, because it requires the knowledge of the tail of the distribution. Therefore, tools from Extreme Value Theory are usually employed, considering that the tail data follow a Generalized Pareto distribution (GPD). Using the existing relations from the parameters of the baseline distribution and the limit GPD's parameters, we define highly informative priors that incorporate all the information available for the whole set of observations. We show how to perform Metropolis-Hastings (MH) algorithm to estimate VaR and CVaR employing the highly informative priors, in the case of exponential, stable and Gamma distributions. Afterwards, we perform a thorough simulation study to compare the accuracy and precision provided by three different methods. Finally, data from a real example is analyzed to show the practical application of the methods. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Using R for teaching and research

Authors: M. Isabel Parra, Eva L. Sanjuán, M. Carmen Robustillo, Mario M. Pizarro

Abstract: R is a language and environment for statistical computing and graphics, which provides a wide variety of statistical tools (modeling, statistical testing, time series analysis, classification problems, machine learning, ...), together with amazing graphical techniques and the great advantage that it is highly extensible. Nowadays, there is no doubt that it is the software par excellence in statist… ▽ More R is a language and environment for statistical computing and graphics, which provides a wide variety of statistical tools (modeling, statistical testing, time series analysis, classification problems, machine learning, ...), together with amazing graphical techniques and the great advantage that it is highly extensible. Nowadays, there is no doubt that it is the software par excellence in statistical courses for any level, for theoretical and applied subjects alike. Besides, it has become an almost essential tool for every research work that involves any kind of analysis or data visualization. Furthermore, it is one of the most employed programming languages for general purposes. The goal of this work is helping to share ideas and resources to improve teaching and/or research using the statistical software R. We will cover its benefits, show how to get started and where to locate specific resources, and will make interesting recommendations for using R, according to our experience. For the classroom we will develop a curricular and assessment infrastructure to support both dissemination and evaluation, while for research we will offer a broader approach to quantitative studies that provides an excellent support for work in science and technology. △ Less

Submitted 21 June, 2023; originally announced June 2023.

X-ray Polarization of the Black Hole X-ray Binary 4U 1630-47 Challenges Standard Thin Accretion Disk Scenario

Authors: Ajay Ratheesh, Michal Dovčiak, Henric Krawczynski, Jakub Podgorný, Lorenzo Marra, Alexandra Veledina, Valery Suleimanov, Nicole Rodriguez Cavero, James Steiner, Jiri Svoboda, Andrea Marinucci, Stefano Bianchi, Michela Negro, Giorgio Matt, Francesco Tombesi, Juri Poutanen, Adam Ingram, Roberto Taverna, Andrew West, Vladimir Karas, Francesco Ursini, Paolo Soffitta, Fiamma Capitanio, Domenico Viscolo, Alberto Manfreda , et al. (90 additional authors not shown)

Abstract: Large energy-dependent X-ray polarization degree is detected by the Imaging X-ray Polarimetry Explorer ({IXPE}) in the high-soft emission state of the black hole X-ray binary 4U 1630--47. The highly significant detection (at $\approx50σ$ confidence level) of an unexpectedly high polarization, rising from $\sim6\%$ at $2$ keV to $\sim10\%$ at $8$ keV, cannot be easily reconciled with standard model… ▽ More Large energy-dependent X-ray polarization degree is detected by the Imaging X-ray Polarimetry Explorer ({IXPE}) in the high-soft emission state of the black hole X-ray binary 4U 1630--47. The highly significant detection (at $\approx50σ$ confidence level) of an unexpectedly high polarization, rising from $\sim6\%$ at $2$ keV to $\sim10\%$ at $8$ keV, cannot be easily reconciled with standard models of thin accretion discs. In this work we compare the predictions of different theoretical models with the {IXPE} data and conclude that the observed polarization properties are compatible with a scenario in which matter accretes onto the black hole through a thin disc, covered by a partially-ionized atmosphere flowing away at mildly relativistic velocities. △ Less

Submitted 19 March, 2024; v1 submitted 25 April, 2023; originally announced April 2023.

Comments: Published in ApJ (https://doi.org/10.3847/1538-4357/ad226e)

Journal ref: 2024 ApJ 964 77

SKA Science Data Challenge 2: analysis and results

Authors: P. Hartley, A. Bonaldi, R. Braun, J. N. H. S. Aditya, S. Aicardi, L. Alegre, A. Chakraborty, X. Chen, S. Choudhuri, A. O. Clarke, J. Coles, J. S. Collinson, D. Cornu, L. Darriba, M. Delli Veneri, J. Forbrich, B. Fraga, A. Galan, J. Garrido, F. Gubanov, H. Håkansson, M. J. Hardcastle, C. Heneka, D. Herranz, K. M. Hess , et al. (83 additional authors not shown)

Abstract: The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed t… ▽ More The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: Under review by MNRAS; 28 pages, 16 figures

Absorption lines from magnetically driven winds in X-ray binaries II: high resolution observational signatures expected from future X-ray observatories

Authors: Susmita Chakravorty, Pierre-Olivier Petrucci, Sudeb Ranjan Datta, Jonathan Ferreira, Joern Wilms, Jonatan Jacquemin-Ide, Maica Clavel, Gregoire Marcel, Jerome Rodriguez, Julien Malzac, Renaud Belmont, Stephane Corbel, Mickael Coriat, Gilles Henri, Maxime Parra

Abstract: In our self-similar, analytical, magneto-hydrodynamic (MHD) accretion-ejection solution, the density at the base of the outflow is explicitly dependent on the disk accretion rate - a unique property of this class of solutions. We had earlier found that the ejection index $p >\sim 0.1 (\dot{M}_{acc} \propto r^p ) $ is a key MHD parameter that decides if the flow can cause absorption lines in the hi… ▽ More In our self-similar, analytical, magneto-hydrodynamic (MHD) accretion-ejection solution, the density at the base of the outflow is explicitly dependent on the disk accretion rate - a unique property of this class of solutions. We had earlier found that the ejection index $p >\sim 0.1 (\dot{M}_{acc} \propto r^p ) $ is a key MHD parameter that decides if the flow can cause absorption lines in the high resolution X-ray spectra of black hole binaries. Here we choose 3 dense warm solutions with $p = 0.1, 0.3, 0.45$ and carefully develop a methodology to generate spectra which are convolved with the Athena and XRISM response functions to predict what they will observe seeing through such MHD outflows. In this paper two other external parameters were varied - extent of the disk, $\rm{r_o|_{max}} = 10^5, \, 10^6 \,\, \rm{r_G}$, and the angle of the line of sight, $i \sim 10 - 25^{\circ}$. Resultant absorption lines (H and He-like Fe, Ca, Ar) change in strength and their profiles manifest varying degrees of asymmetry. We checked if a) the lines and ii) the line asymmetries are detected, in our suit of synthetic Athena and XRISM spectra. Our analysis shows that Athena should detect the lines and their asymmetries for a standard 100 ksec observation of a 100 mCrab source - lines with equivalent width as low as a few eV should be detected if the 6-8 keV counts are larger than $10^4 - 10^5$ even for the least favourable simulated cases. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 18 pages, 13 figures in the main body and 3 figures in the appendix. Accepted for publication in MNRAS

Equidistribution for matings of quadratic maps with the Modular group

Authors: Vanessa Matus de la Parra

Abstract: We study the asymptotic behavior of the family of holomorphic correspondences $\lbrace\mathcal{F}_a\rbrace_{a\in\mathcal{K}}$, given by $$\left(\frac{az+1}{z+1}\right)^2+\left(\frac{az+1}{z+1}\right)\left(\frac{aw-1}{w-1}\right)+\left(\frac{aw-1}{w-1}\right)^2=3.$$ It was proven by Bullet and Lomonaco that $\mathcal{F}_a$ is a mating between the modular group $\operatorname{PSL}_2(\mathbb{Z})$ and… ▽ More We study the asymptotic behavior of the family of holomorphic correspondences $\lbrace\mathcal{F}_a\rbrace_{a\in\mathcal{K}}$, given by $$\left(\frac{az+1}{z+1}\right)^2+\left(\frac{az+1}{z+1}\right)\left(\frac{aw-1}{w-1}\right)+\left(\frac{aw-1}{w-1}\right)^2=3.$$ It was proven by Bullet and Lomonaco that $\mathcal{F}_a$ is a mating between the modular group $\operatorname{PSL}_2(\mathbb{Z})$ and a quadratic rational map. We show for every $a\in\mathcal{K}$, the iterated images and preimages under $\mathcal{F}_a$ of nonexceptional points equidistribute, in spite of the fact that $\mathcal{F}_a$ is weakly-modular in the sense of Dinh, Kaufmann and Wu but it is not modular. Furthermore, we prove that periodic points equidistribute as well. △ Less

Submitted 27 December, 2022; v1 submitted 14 February, 2022; originally announced February 2022.

Journal ref: Ergodic Theory and Dynamical Systems, 1-29

MUSE spectroscopy of the ULX NGC 1313 X-1: a shock-ionised bubble, an X-ray photoionised nebula, and two supernova remnants

Authors: Andrés Gúrpide, Maxime Parra, Olivier Godet, Thierry Contini, Jean-François Olive

Abstract: The presence of large ionised gaseous nebulae found around some ultraluminous X-ray sources (ULXs) provides the means to assess the mechanical and radiative feedback of the central source, and hence constrain the efficiency and impact on the surroundings of the super-Eddington regime powering most of these sources. NGC 1313 X--1 is an archetypal ULX which has been reported to be surrounded by abno… ▽ More The presence of large ionised gaseous nebulae found around some ultraluminous X-ray sources (ULXs) provides the means to assess the mechanical and radiative feedback of the central source, and hence constrain the efficiency and impact on the surroundings of the super-Eddington regime powering most of these sources. NGC 1313 X--1 is an archetypal ULX which has been reported to be surrounded by abnormally high [O I]$λ$6300/H$α>$ 0.1 ratios and for which high-resolution spectroscopy in X-rays has hinted the presence of powerful outflows. We report observations taken with the Multi-Unit Spectroscopic Explorer of NGC 1313 X--1 in order to confirm the presence of a nebula inflated by the winds, investigate its main sources of ionisation and estimate the mechanical output of the source. We detect a bubble of 452 $\times$ 266\,pc in size, roughly centred around the ULX, which shows clear evidence for shock ionisation in the outer edges. We estimate shock velocities to be in the $\approx160-180$\,km/s range based on the line ratios. This suggests that an average and continuous outflow power of $\sim(2-4.5) \times 10^{40}$ erg/s over a timescale of $(4.5-7.8) \times10^5$\,yr is required to inflate the bubble. In the interior of the bubble and closer to the ULX we detect an extended ($\sim$140\,pc) X-ray ionised region. Additionally, we detect two supernova remnants coincidentally close to the bubble of which we also report age and explosion energy estimates. The elongated morphology and the kinematics of the bubble strongly suggest that the bubble is being inflated by winds and/or jets emanating from the central source, supporting the presence of winds found through X-ray spectroscopy. The estimated mechanical power is comparable or higher than the X-ray luminosity of the source, providing additional evidence in support of NGC 1313 X--1 harbouring a super-Eddington accretor. △ Less

Submitted 8 February, 2022; v1 submitted 23 January, 2022; originally announced January 2022.

Comments: 22 pages, 20 Figures, 5 Tables. Accepted for publication to A&A, replaced after revision

Journal ref: A&A 666, A100 (2022)

Climate-driven trends in the streamflow records of a reference hydrologic network in Southern Spain

Authors: Patricio Yeste, Javier Dorador, Wenceslao Martín-Rosales, Emilio Molero, María Jesús Esteban Parra

Abstract: Monthly streamflow records from a set of gauging stations, selected to form a reference hydrologic network, are analyzed together with precipitation and temperature data to establish whether the streamflows in the Guadalquivir River Basin have experienced changes during the last half of the XXth century that can be attributed to hydrological forcing. The observed streamflows in the reference netwo… ▽ More Monthly streamflow records from a set of gauging stations, selected to form a reference hydrologic network, are analyzed together with precipitation and temperature data to establish whether the streamflows in the Guadalquivir River Basin have experienced changes during the last half of the XXth century that can be attributed to hydrological forcing. The observed streamflows in the reference network have undergone generalized and significant decreases both at seasonal and annual scales during the study period. Annual rainfall, though, did not experienced statistically significant changes. The observed trends in streamflows can be attributed to either land-use changes, or to the statistically significant changes exhibited both by yearly potential evapotranspiration values and by the seasonal distribution of precipitation. In the attribution work conducted using both data-based and simulation-based methods, the intra-annual redistribution of precipitation is shown to be the main statistically significant climate-driver of streamflow change. The contributions of other non-climate factors, such as the change in land cover, to the reduction in annual streamflows are shown to be minor in comparison. △ Less

Submitted 18 August, 2020; originally announced August 2020.

Comments: Published in Journal of Hydrology

Journal ref: Journal of Hydrology 566 (2018) 55-72

Understanding the HERA Phase I receiver system with simulations and its impact on the detectability of the EoR delay power spectrum

Authors: Nicolas Fagnoni, Eloy de Lera Acedo, David R. DeBoer, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Phil Bull, Jacob Burba, Chris L. Carilli, Carina Cheng, Matt Dexter, Joshua S. Dillon, Aaron Ewall-Wice, Randall Fritz, Steve R. Furlanetto, Kingsley Gale-Sides, Brian Glendenning, Deepthi Gorthi , et al. (45 additional authors not shown)

Abstract: The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Applied to the Hydrogen Epoch of Reionization Array (HERA), this paper combines d… ▽ More The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Applied to the Hydrogen Epoch of Reionization Array (HERA), this paper combines detailed electromagnetic and electrical simulations in order to model the chromatic effects of the instrument, and quantify its frequency and time responses. In particular, the effects of the analogue receiver, transmission cables, and mutual coupling are included. These simulations are able to accurately predict the intensity of the reflections occurring in the 150-m cable which links the antenna to the back-end. They also show that electromagnetic waves can propagate from one dish to another one through large sections of the array due to mutual coupling. The simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h\;\rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity. △ Less

Submitted 25 August, 2020; v1 submitted 6 August, 2019; originally announced August 2019.

Comments: 13 pages, 17 figures - Submitted to MNRAS - 2nd revision

An improved method for the estimation of the Gumbel distribution parameters

Authors: Rubén Gómez González, M. Isabel Parra, Francisco Javier Acero, Jacinto Martín

Abstract: Usual estimation methods for the parameters of extreme values distribution employ only a few values, wasting a lot of information. More precisely, in the case of the Gumbel distribution, only the block maxima values are used. In this work, we propose a method to seize all the available information in order to increase the accuracy of the estimations. This intent can be achieved by taking advantage… ▽ More Usual estimation methods for the parameters of extreme values distribution employ only a few values, wasting a lot of information. More precisely, in the case of the Gumbel distribution, only the block maxima values are used. In this work, we propose a method to seize all the available information in order to increase the accuracy of the estimations. This intent can be achieved by taking advantage of the existing relationship between the parameters of the baseline distribution, which generates data from the full sample space, and the ones for the limit Gumbel distribution. In this way, an informative prior distribution can be obtained. Different statistical tests are used to compare the behaviour of our method with the standard one, showing that the proposed method performs well when dealing with very shortened available data. The empirical effectiveness of the approach is demonstrated through a simulation study and a case study. Reduction in the credible interval width and enhancement in parameter location show that the results with improved prior adapt to very shortened data better than standard method does. △ Less

Submitted 21 February, 2019; originally announced February 2019.

Comments: 21 pages, 10 figures

Wave attenuation and dispersion due to floating ice covers

Authors: L. J. Yiew, S. M. Parra, D. Wang, D. K. K. Sree, A. V. Babanin, A. W. -K. Law

Abstract: Experiments investigating the attenuation and dispersion of surface waves in a variety of ice covers are performed using a refrigerated wave flume. The ice conditions tested in the experiments cover naturally occurring combinations of continuous, fragmented, pancake and grease ice. Attenuation rates are shown to be a function of ice thickness, wave frequency, and the general rigidity of the ice co… ▽ More Experiments investigating the attenuation and dispersion of surface waves in a variety of ice covers are performed using a refrigerated wave flume. The ice conditions tested in the experiments cover naturally occurring combinations of continuous, fragmented, pancake and grease ice. Attenuation rates are shown to be a function of ice thickness, wave frequency, and the general rigidity of the ice cover. Dispersion changes were minor except for large wavelength increases when continuous covers were tested. Results are verified and compared with existing literature to show the extended range of investigation in terms of incident wave frequency and ice conditions. △ Less

Submitted 14 February, 2019; originally announced February 2019.

Scaling notifications beyond alerts: from subtly drawing attention up to forcing the user to take action

Authors: Denys J. C. Matthies, Laura Milena Daza Parra, Bodo Urban

Abstract: New computational devices, in particular wearable devices, offer the unique property of always being available and thus to be able to constantly update the user with information, such as by notifications. While research has been done in sophisticated notifications, devices today mainly stick to a binary level of information, while they are either attention drawing or silent. In this paper, we want… ▽ More New computational devices, in particular wearable devices, offer the unique property of always being available and thus to be able to constantly update the user with information, such as by notifications. While research has been done in sophisticated notifications, devices today mainly stick to a binary level of information, while they are either attention drawing or silent. In this paper, we want to go further and propose scalable notifications, which adjust the intensity reaching from subtle to obtrusive and even going beyond that level, while forcing the user to take action. To illustrate the technical feasibility and validity of this concept, we developed three prototypes providing mechano-pressure, thermal, and electrical feedback and evaluated them in different lab studies. Our first prototype provides subtle poking through to high and frequent pressure on the user's spine, which creates a significantly improved back posture. In a second scenario, the users are enabled to perceive the overuse of a drill by an increased temperature on the palm of a hand until the heat is intolerable and the users are forced to eventually put down the tool. The last project comprises a speed control in a driving simulation, while electric muscle stimulation on the users' legs conveys information on changing the car's speed by a perceived tingling until the system independently forces the foot to move. Although our selected scenarios are long way from being realistic, we see these lab studies as a means to validate our proof-of-concept. In conclusion, all studies' findings support the feasibility of our concept of a scalable notification system, including the system of forced intervention. While we envisage the implementation of our proof-of-concept into future wearables, more realistic application scenarios are worthy of exploration. △ Less

Submitted 6 August, 2018; originally announced August 2018.

ACM Class: H.5.m

Complex tensor factorisation with PARAFAC2 for the estimation of brain connectivity from the EEG

Authors: Loukianos Spyrou, Mario Parra, Javier Escudero

Abstract: Objective: The coupling between neuronal populations and its magnitude have been shown to be informative for various clinical applications. One method to estimate brain connectivity is with electroencephalography (EEG) from which the cross-spectrum between different sensor locations is derived. We wish to test the efficacy of tensor factorisation in the estimation of brain connectivity. Methods: C… ▽ More Objective: The coupling between neuronal populations and its magnitude have been shown to be informative for various clinical applications. One method to estimate brain connectivity is with electroencephalography (EEG) from which the cross-spectrum between different sensor locations is derived. We wish to test the efficacy of tensor factorisation in the estimation of brain connectivity. Methods: Complex tensor factorisation based on PARAFAC2 is used to decompose the EEG into scalp components described by the spatial, spectral, and complex trial profiles. An EEG model in the complex domain was derived that shows the suitability of PARAFAC2. A connectivity metric was also derived on the complex trial profiles of the extracted components. Results: Results on a benchmark EEG dataset confirmed that PARAFAC2 can estimate connectivity better than traditional tensor analysis such as PARAFAC within a range of signal-to-noise ratios. The analysis of EEG from patients with mild cognitive impairment or Alzheimer's disease showed that PARAFAC2 identifies loss of brain connectivity better than traditional approaches and agreeing with prior pathological knowledge. Conclusion: The complex PARAFAC2 algorithm is suitable for EEG connectivity estimation since it allows to extract meaningful coupled sources and provides better estimates than complex PARAFAC. Significance: A new paradigm that employs complex tensor factorisation has demonstrated to be successful in identifying brain connectivity and the location of couples sources for both a benchmark and a real-world EEG dataset. This can enable future applications and has the potential to solve some the issues that deteriorate the performance of traditional connectivity metrics. △ Less

Submitted 2 May, 2017; originally announced May 2017.

ICE: a scalable, low-cost FPGA-based telescope signal processing and networking system

Authors: K. Bandura, A. N. Bender, J. F. Cliche, T. de Haan, M. A. Dobbs, A. J. Gilbert, S. Griffin, G. Hsyu, D. Ittah, J. Mena Parra, J. Montgomery, T. Pinsonneault-Marotte, S. Siegel, G. Smecher, Q. Y. Tang, K. Vanderlinde, N. Whitehorn

Abstract: We present an overview of the 'ICE' hardware and software framework that implements large arrays of interconnected FPGA-based data acquisition, signal processing and networking nodes economically. The system was conceived for application to radio, millimeter and sub-millimeter telescope readout systems that have requirements beyond typical off-the-shelf processing systems, such as careful control… ▽ More We present an overview of the 'ICE' hardware and software framework that implements large arrays of interconnected FPGA-based data acquisition, signal processing and networking nodes economically. The system was conceived for application to radio, millimeter and sub-millimeter telescope readout systems that have requirements beyond typical off-the-shelf processing systems, such as careful control of interference signals produced by the digital electronics, and clocking of all elements in the system from a single precise observatory-derived oscillator. A new generation of telescopes operating at these frequency bands and designed with a vastly increased emphasis on digital signal processing to support their detector multiplexing technology or high-bandwidth correlators---data rates exceeding a terabyte per second---are becoming common. The ICE system is built around a custom FPGA motherboard that makes use of an Xilinx Kintex-7 FPGA and ARM-based co-processor. The system is specialized for specific applications through software, firmware, and custom mezzanine daughter boards that interface to the FPGA through the industry-standard FMC specifications. For high density applications, the motherboards are packaged in 16-slot crates with ICE backplanes that implement a low-cost passive full-mesh network between the motherboards in a crate, allow high bandwidth interconnection between crates, and enable data offload to a computer cluster. A Python-based control software library automatically detects and operates the hardware in the array. Examples of specific telescope applications of the ICE framework are presented, namely the frequency-multiplexed bolometer readout systems used for the SPT and Simons Array and the digitizer, F-engine, and networking engine for the CHIME and HIRAX radio interferometers. △ Less

Submitted 22 August, 2016; originally announced August 2016.

Comments: 20 pages, 8 figures. Submitted to JAI special issue on Digital Signal Processing in Radio Astronomy (2016)

ICE-based Custom Full-Mesh Network for the CHIME High Bandwidth Radio Astronomy Correlator

Authors: Kevin Bandura, Jean-Francois Cliche, Matt Dobbs, Adam Gilbert, David Ittah, Juan Mena Parra, Graeme Smecher

Abstract: New generation radio interferometers encode signals from thousands of antenna feeds across large bandwidth. Channelizing and correlating this data requires networking capabilities that can handle unprecedented data rates with reasonable cost. The Canadian Hydrogen Intensity Mapping Experiment (CHIME) correlator processes 8-bits from N=2048 digitizer inputs across 400~MHz of bandwidth. Measured in… ▽ More New generation radio interferometers encode signals from thousands of antenna feeds across large bandwidth. Channelizing and correlating this data requires networking capabilities that can handle unprecedented data rates with reasonable cost. The Canadian Hydrogen Intensity Mapping Experiment (CHIME) correlator processes 8-bits from N=2048 digitizer inputs across 400~MHz of bandwidth. Measured in $N^2~\times $ bandwidth, it is the largest radio correlator that has been built. Its digital back-end must exchange and reorganize the 6.6~terabit/s produced by its 128 digitizing and channelizing nodes, and feed it to the 256-node spatial correlator in a way that each node obtains data from all digitizer inputs but across a small fraction of the bandwidth (i.e. `corner-turn'). In order to maximize performance and reliability of the corner-turn system while minimizing cost, a custom networking solution has been implemented. The system makes use of Field Programmable Gate Array (FPGA) transceivers to implement direct, passive, full-mesh, high speed serial connections between sixteen circuit boards in a crate, to exchange data between crates, and to offload the data to a cluster of 256 graphics processing unit (GPU) nodes using standard 10~Gbit/s Ethernet links. The GPU nodes complete the corner-turn by combining data from all crates and then computing visibilities. Eye diagrams and frame error counters confirm error-free operation of the corner-turn network in both the currently operating CHIME Pathfinder telescope (a prototype for the full CHIME telescope) and a representative fraction of the full CHIME hardware providing an end-to-end system validation. An analysis of an equivalent corner-turn system built with Ethernet switches instead of custom passive data links is provided. △ Less

Submitted 15 August, 2016; originally announced August 2016.

HIRAX: A Probe of Dark Energy and Radio Transients

Authors: L. B. Newburgh, K. Bandura, M. A. Bucher, T. -C. Chang, H. C. Chiang, J. F. Cliche, R. Dave, M. Dobbs, C. Clarkson, K. M. Ganga, T. Gogo, A. Gumba, N. Gupta, M. Hilton, B. Johnstone, A. Karastergiou, M. Kunz, D. Lokhorst, R. Maartens, S. Macpherson, M. Mdlalose, K. Moodley, L. Ngwenya, J. M. Parra, J. Peterson , et al. (11 additional authors not shown)

Abstract: The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a new 400-800MHz radio interferometer under development for deployment in South Africa. HIRAX will comprise 1024 six meter parabolic dishes on a compact grid and will map most of the southern sky over the course of four years. HIRAX has two primary science goals: to constrain Dark Energy and measure structure at high redshift, and… ▽ More The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a new 400-800MHz radio interferometer under development for deployment in South Africa. HIRAX will comprise 1024 six meter parabolic dishes on a compact grid and will map most of the southern sky over the course of four years. HIRAX has two primary science goals: to constrain Dark Energy and measure structure at high redshift, and to study radio transients and pulsars. HIRAX will observe unresolved sources of neutral hydrogen via their redshifted 21-cm emission line (`hydrogen intensity mapping'). The resulting maps of large-scale structure at redshifts 0.8-2.5 will be used to measure Baryon Acoustic Oscillations (BAO). HIRAX will improve upon current BAO measurements from galaxy surveys by observing a larger cosmological volume (larger in both survey area and redshift range) and by measuring BAO at higher redshift when the expansion of the universe transitioned to Dark Energy domination. HIRAX will complement CHIME, a hydrogen intensity mapping experiment in the Northern Hemisphere, by completing the sky coverage in the same redshift range. HIRAX's location in the Southern Hemisphere also allows a variety of cross-correlation measurements with large-scale structure surveys at many wavelengths. Daily maps of a few thousand square degrees of the Southern Hemisphere, encompassing much of the Milky Way galaxy, will also open new opportunities for discovering and monitoring radio transients. The HIRAX correlator will have the ability to rapidly and eXperimentciently detect transient events. This new data will shed light on the poorly understood nature of fast radio bursts (FRBs), enable pulsar monitoring to enhance long-wavelength gravitational wave searches, and provide a rich data set for new radio transient phenomena searches. This paper discusses the HIRAX instrument, science goals, and current status. △ Less

Submitted 7 July, 2016; originally announced July 2016.

Comments: 11 pages, 5 figures

Journal ref: Proc. SPIE 9906, Ground-based and Airborne Telescopes VI (2016)

Holographic Beam Mapping of the CHIME Pathfinder Array

Authors: Philippe Berger, Laura B. Newburgh, Mandana Amiri, Kevin Bandura, Jean-Francois Cliche, Liam Connor, Meiling Deng, Nolan Denman, Matt Dobbs, Mateus Fandino, Adam J. Gilbert, Deborah Good, Mark Halpern, David Hanna, Adam D. Hincks, Gary Hinshaw, Carolin Hofer, Andre M. Johnson, Tom L. Landecker, Kiyoshi W. Masui, Juan Mena Parra, Niels Oppermann, Ue-Li Pen, Jeffrey B. Peterson, Andre Recnik , et al. (10 additional authors not shown)

Abstract: The Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder radio telescope is currently surveying the northern hemisphere between 400 and 800 MHz. By mapping the large scale structure of neutral hydrogen through its redshifted 21 cm line emission between $z \sim 0.8-2.5$ CHIME will contribute to our understanding of Dark Energy. Bright astrophysical foregrounds must be separated from th… ▽ More The Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder radio telescope is currently surveying the northern hemisphere between 400 and 800 MHz. By mapping the large scale structure of neutral hydrogen through its redshifted 21 cm line emission between $z \sim 0.8-2.5$ CHIME will contribute to our understanding of Dark Energy. Bright astrophysical foregrounds must be separated from the neutral hydrogen signal, a task which requires precise characterization of the polarized telescope beams. Using the DRAO John A. Galt 26 m telescope, we have developed a holography instrument and technique for mapping the CHIME Pathfinder beams. We report the status of the instrument and initial results of this effort. △ Less

Submitted 5 July, 2016; originally announced July 2016.

Comments: 16 pages, 9 figures. Submitted to Proc. SPIE, Astronomical Telescopes + Instrumentation (2016)

Report number: Proc. SPIE 9906, Ground-based and Airborne Telescopes VI, 99060D (August 18, 2016)

Locating Temporal Functional Dynamics of Visual Short-Term Memory Binding using Graph Modular Dirichlet Energy

Authors: Keith Smith, Benjamin Ricaud, Nauman Shahid, Stephen Rhodes, John M. Starr, Agustin Ibanez, Mario A. Parra, Javier Escudero, Pierre Vandergheynst

Abstract: Visual short-term memory binding tasks are a promising early marker for Alzheimer's disease (AD). To uncover functional deficits of AD in these tasks it is meaningful to first study unimpaired brain function. Electroencephalogram recordings were obtained from encoding and maintenance periods of tasks performed by healthy young volunteers. We probe the task's transient physiological underpinnings b… ▽ More Visual short-term memory binding tasks are a promising early marker for Alzheimer's disease (AD). To uncover functional deficits of AD in these tasks it is meaningful to first study unimpaired brain function. Electroencephalogram recordings were obtained from encoding and maintenance periods of tasks performed by healthy young volunteers. We probe the task's transient physiological underpinnings by contrasting shape only (Shape) and shape-colour binding (Bind) conditions, displayed in the left and right sides of the screen, separately. Particularly, we introduce and implement a novel technique named Modular Dirichlet Energy (MDE) which allows robust and flexible analysis of the functional network with unprecedented temporal precision. We find that connectivity in the Bind condition is less integrated with the global network than in the Shape condition in occipital and frontal modules during the encoding period of the right screen condition. Using MDE we are able to discern driving effects in the occipital module between 100-140ms, coinciding with the P100 visually evoked potential, followed by a driving effect in the frontal module between 140-180ms, suggesting that the differences found constitute an information processing difference between these modules. This provides temporally precise information over a heterogeneous population in promising tasks for the detection of AD. △ Less

Submitted 8 September, 2016; v1 submitted 8 June, 2016; originally announced June 2016.

Journal ref: K. Smith et al., Scientific Reports, 7: 42013, 2017

Comparison of Network Analysis Approaches on EEG Connectivity in Beta during Visual Short-Term Memory Binding Tasks

Authors: Keith Smith, Hamed Azami, Mario A. Parra, Javier Escudero, John M. Starr

Abstract: We analyse the electroencephalogram signals in the beta band of working memory representation recorded from young healthy volunteers performing several different Visual Short-Term Memory (VSTM) tasks which have proven useful in the assessment of clinical and preclinical Alzheimer's disease. We compare network analysis using Maximum Spanning Trees (MSTs) with network analysis obtained using 20% and… ▽ More We analyse the electroencephalogram signals in the beta band of working memory representation recorded from young healthy volunteers performing several different Visual Short-Term Memory (VSTM) tasks which have proven useful in the assessment of clinical and preclinical Alzheimer's disease. We compare network analysis using Maximum Spanning Trees (MSTs) with network analysis obtained using 20% and 25% connection thresholds on the VSTM data. MSTs are a promising method of network analysis negating the more classical use of thresholds which are so far chosen arbitrarily. However, we find that the threshold analyses outperforms MSTs for detection of functional network differences. Particularly, MSTs fail to find any significant differences. Further, the thresholds detect significant differences between shape and shape-colour binding tasks when these are tested in the left side of the display screen, but no such differences are detected when these tasks are tested for in the right side of the display screen. This provides evidence that contralateral activity is a significant factor in sensitivity for detection of cognitive task differences. △ Less

Submitted 8 April, 2016; originally announced April 2016.

Comments: Proceeds from the IEEE Engineering in Medicine and Biology Conference 2015

Cluster-Span Threshold: An unbiased threshold for binarising weighted complete networks in functional connectivity analysis

Authors: Keith Smith, Hamed Azami, Mario A. Parra, John M. Starr, Javier Escudero

Abstract: We propose a new unbiased threshold for network analysis named the Cluster-Span Threshold (CST). This is based on the clustering coefficient, C, following logic that a balance of `clustering' to `spanning' triples results in a useful topology for network analysis and that the product of complementing properties has a unique value only when perfectly balanced. We threshold networks by fixing C at t… ▽ More We propose a new unbiased threshold for network analysis named the Cluster-Span Threshold (CST). This is based on the clustering coefficient, C, following logic that a balance of `clustering' to `spanning' triples results in a useful topology for network analysis and that the product of complementing properties has a unique value only when perfectly balanced. We threshold networks by fixing C at this balanced value, rather than fixing connection density at an arbitrary value, as has been the trend. We compare results from an electroencephalogram data set of volunteers performing visual short term memory tasks of the CST alongside other thresholds, including maximum spanning trees. We find that the CST holds as a sensitive threshold for distinguishing differences in the functional connectivity between tasks. This provides a sensitive and objective method for setting a threshold on weighted complete networks which may prove influential on the future of functional connectivity research. △ Less

Submitted 8 April, 2016; originally announced April 2016.

Comments: conference paper from the IEEE Engineering in Medicine and Biology Conference 2015

A GPU-based Correlator X-engine Implemented on the CHIME Pathfinder

Authors: Nolan Denman, Mandana Amiri, Kevin Bandura, Jean-François Cliche, Liam Connor, Matt Dobbs, Mateus Fandino, Mark Halpern, Adam Hincks, Gary Hinshaw, Carolin Höfer, Peter Klages, Kiyoshi Masui, Juan Mena Parra, Laura Newburgh, Andre Recnik, J. Richard Shaw, Kris Sigurdson, Kendrick Smith, Keith Vanderlinde

Abstract: We present the design and implementation of a custom GPU-based compute cluster that provides the correlation X-engine of the CHIME Pathfinder radio telescope. It is among the largest such systems in operation, correlating 32,896 baselines (256 inputs) over 400MHz of radio bandwidth. Making heavy use of consumer-grade parts and a custom software stack, the system was developed at a small fraction o… ▽ More We present the design and implementation of a custom GPU-based compute cluster that provides the correlation X-engine of the CHIME Pathfinder radio telescope. It is among the largest such systems in operation, correlating 32,896 baselines (256 inputs) over 400MHz of radio bandwidth. Making heavy use of consumer-grade parts and a custom software stack, the system was developed at a small fraction of the cost of comparable installations. Unlike existing GPU backends, this system is built around OpenCL kernels running on consumer-level AMD GPUs, taking advantage of low-cost hardware and leveraging packed integer operations to double algorithmic efficiency. The system achieves the required 105TOPS in a 10kW power envelope, making it among the most power-efficient X-engines in use today. △ Less

Submitted 11 June, 2015; v1 submitted 20 March, 2015; originally announced March 2015.

Comments: 6 pages, 5 figures. Accepted by IEEE ASAP 2015

A compression scheme for radio data in high performance computing

Authors: Kiyoshi Masui, Mandana Amiri, Liam Connor, Meiling Deng, Mateus Fandino, Carolin Höfer, Mark Halpern, David Hanna, Adam D. Hincks, Gary Hinshaw, Juan Mena Parra, Laura B. Newburgh, J. Richard Shaw, Keith Vanderlinde

Abstract: We present a procedure for efficiently compressing astronomical radio data for high performance applications. Integrated, post-correlation data are first passed through a nearly lossless rounding step which compares the precision of the data to a generalized and calibration-independent form of the radiometer equation. This allows the precision of the data to be reduced in a way that has an insigni… ▽ More We present a procedure for efficiently compressing astronomical radio data for high performance applications. Integrated, post-correlation data are first passed through a nearly lossless rounding step which compares the precision of the data to a generalized and calibration-independent form of the radiometer equation. This allows the precision of the data to be reduced in a way that has an insignificant impact on the data. The newly developed Bitshuffle lossless compression algorithm is subsequently applied. When the algorithm is used in conjunction with the HDF5 library and data format, data produced by the CHIME Pathfinder telescope is compressed to 28% of its original size and decompression throughputs in excess of 1 GB/s are obtained on a single core. △ Less

Submitted 16 September, 2015; v1 submitted 2 March, 2015; originally announced March 2015.

Comments: 12 pages, 3 figures, 2 tables. As published in Astronomy and Computing special issue on "The future of astronomical data formats"

Calibrating CHIME, A New Radio Interferometer to Probe Dark Energy

Authors: Laura B. Newburgh, Graeme E. Addison, Mandana Amiri, Kevin Bandura, J. Richard Bond, Liam Connor, Jean-François Cliche, Greg Davis, Meiling Deng, Nolan Denman, Matt Dobbs, Mateus Fandino, Heather Fong, Kenneth Gibbs, Adam Gilbert, Elizabeth Griffin, Mark Halpern, David Hanna, Adam D. Hincks, Gary Hinshaw, Carolin Höfer, Peter Klages, Tom Landecker, Kiyoshi Masui, Juan Mena Parra , et al. (10 additional authors not shown)

Abstract: The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark… ▽ More The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark energy. We have deployed a pathfinder version of CHIME that will yield constraints on the BAO power spectrum and provide a test-bed for our calibration scheme. I will discuss the CHIME calibration requirements and describe instrumentation we are developing to meet these requirements. △ Less

Submitted 9 June, 2014; originally announced June 2014.

Journal ref: Proc. SPIE 9145, Ground-based and Airborne Telescopes V (2014)

The Jacobian cocycle and equidistribution towards the Green current

Authors: Manuel Rodrigo Parra

Abstract: Given a holomorphic self-map of complex projective space of degree larger than one, we prove that there exists a finite collection of totally invariant algebraic sets with the following property: given any positive closed (1,1)-current of mass 1 with no mass on any element of this family, the sequence of normalized pull-backs of the current converges to the Green current. Under suitable geometric… ▽ More Given a holomorphic self-map of complex projective space of degree larger than one, we prove that there exists a finite collection of totally invariant algebraic sets with the following property: given any positive closed (1,1)-current of mass 1 with no mass on any element of this family, the sequence of normalized pull-backs of the current converges to the Green current. Under suitable geometric conditions on the collection of totally invariant algebraic sets, we prove a sharper equidistribution result. △ Less

Submitted 23 March, 2011; originally announced March 2011.

Lelong numbers on projective varieties

Authors: Manuel Rodrigo Parra

Abstract: Given a positive closed (1,1)-current $T$ defined on the regular locus of a projective variety $X$ with bounded mass near the singular part of $X$ and $Y$ an irreducible algebraic subset of $X$, we present uniform estimates for the locus inside $Y$ where the Lelong numbers of $T$ are larger than the generic Lelong number of $T$ along $Y$. Given a positive closed (1,1)-current $T$ defined on the regular locus of a projective variety $X$ with bounded mass near the singular part of $X$ and $Y$ an irreducible algebraic subset of $X$, we present uniform estimates for the locus inside $Y$ where the Lelong numbers of $T$ are larger than the generic Lelong number of $T$ along $Y$. △ Less

Submitted 23 November, 2010; originally announced November 2010.

Bouchaud walks with variable drift

Authors: Manuel Cabezas Parra

Abstract: In this paper we study a sequence of Bouchaud trap models on $\mathbb{Z}$ with drift. We analyze the possible scaling limits for a sequence of walks, where we make the drift decay to 0 as we rescale the walks. Depending on the speed of the decay of the drift we obtain three different scaling limits. If the drift decays slowly as we rescale the walks we obtain the inverse of an α… ▽ More In this paper we study a sequence of Bouchaud trap models on $\mathbb{Z}$ with drift. We analyze the possible scaling limits for a sequence of walks, where we make the drift decay to 0 as we rescale the walks. Depending on the speed of the decay of the drift we obtain three different scaling limits. If the drift decays slowly as we rescale the walks we obtain the inverse of an α$-stable subordinator as scaling limit. If the drift decays quickly as we rescale the walks, we obtain the F.I.N. diffusion as scaling limit. There is a critical speed of decay separating these two main regimes, where a new process appears as scaling limit. This critical speed is related to the index $α$ of the inhomogeneity of the environment. △ Less

Submitted 26 November, 2010; v1 submitted 25 March, 2010; originally announced March 2010.

MSC Class: 60K37; 60G52; 60F17

Positivity and conservation of superenergy tensors

Authors: Jose M. Pozo, Josep M. Parra

Abstract: Two essential properties of energy-momentum tensors T_{μν} are their positivity and conservation. This is mathematically formalized by, respectively, an energy condition, as the dominant energy condition, and the vanishing of their divergence \nabla^μT_{μν}=0. The classical Bel and Bel-Robinson superenergy tensors, generated from the Riemann and Weyl tensors, respectively, are rank-4 tensors. Bu… ▽ More Two essential properties of energy-momentum tensors T_{μν} are their positivity and conservation. This is mathematically formalized by, respectively, an energy condition, as the dominant energy condition, and the vanishing of their divergence \nabla^μT_{μν}=0. The classical Bel and Bel-Robinson superenergy tensors, generated from the Riemann and Weyl tensors, respectively, are rank-4 tensors. But they share these two properties with energy momentum tensors: the Dominant Property (DP) and the divergence-free property in the absence of sources (vacuum). Senovilla defined a universal algebraic construction which generates a basic superenergy tensor T{A} from any arbitrary tensor A. In this construction the seed tensor A is structured as an r-fold multivector, which can always be done. The most important feature of the basic superenergy tensors is that they satisfy automatically the DP, independently of the generating tensor A. In a previous paper we presented a more compact definition of T{A} using the r-fold Clifford algebra. This form for the superenergy tensors allowed to obtain an easy proof of the DP valid for any dimension. In this paper we include this proof. We explain which new elements appear when we consider the tensor T{A} generated by a non-degree-defined r-fold multivector A and how orthogonal Lorentz transformations and bilinear observables of spinor fields are included as particular cases of superenergy tensors. We find some sufficient conditions for the seed tensor A, which guarantee that the generated tensor T{A} is divergence-free. These sufficient conditions are satisfied by some physical fields, which are presented as examples. △ Less

Submitted 25 February, 2002; v1 submitted 9 October, 2001; originally announced October 2001.

Comments: 19 pages, no figures. Language and minor changes. Published version

Journal ref: Class.Quant.Grav. 19 (2002) 967-984

Signature Change and Clifford Algebras

Authors: David Miralles, Josep Manel Parra, Jayme Vaz Jr

Abstract: Given the real Clifford algebra of a quadratic space with a given signature, we define a new product in this structure such that it simulates the Clifford product of a quadratic space with another signature different from the original one. Among the possible applications of this new product, we use it in order to write the minkowskian Dirac equation over the euclidean spacetime and to define a n… ▽ More Given the real Clifford algebra of a quadratic space with a given signature, we define a new product in this structure such that it simulates the Clifford product of a quadratic space with another signature different from the original one. Among the possible applications of this new product, we use it in order to write the minkowskian Dirac equation over the euclidean spacetime and to define a new duality operation in terms of which one can find self-dual and anti-self-dual solutions of gauge fields over Minkowski spacetime analogous to the ones over Euclidean spacetime and without needing to complexify the original real algebra. △ Less

Submitted 29 March, 2000; originally announced March 2000.

Comments: 11 pages, Latex file, no figures. To appear in International Journal of Theoretical Physics

MSC Class: 15A66 (Primary); 15A63 (Secondary)

Journal ref: Int.J.Theor.Phys. 40 (2001) 229-242

Clifford Algebra Approach to Superenergy Tensors

Authors: Jose M. Pozo, Josep M. Parra

Abstract: Senovilla has recently defined an algebraic construction of a superenergy tensor T{A} from any arbitrary tensor A, by structuring it as an r-fold form. This superenergy tensor satisfies automatically the dominant superenergy property. We present a more compact definition using the r-direct product Clifford algebra r-Cl(p,q). This form for the superenergy tensors allows to obtain an easy proof of… ▽ More Senovilla has recently defined an algebraic construction of a superenergy tensor T{A} from any arbitrary tensor A, by structuring it as an r-fold form. This superenergy tensor satisfies automatically the dominant superenergy property. We present a more compact definition using the r-direct product Clifford algebra r-Cl(p,q). This form for the superenergy tensors allows to obtain an easy proof of the dominant superenergy property valid for any dimension. △ Less

Submitted 10 November, 1999; originally announced November 1999.

Comments: 6 pages, LaTeX, uses ere99.sty, contribution to proceedings of ERE99, Bilbao, Sep 7th-10th