Showing 1–41 of 41 results for author: Kolb, C

How Inverse Conditional Flows Can Serve as a Substitute for Distributional Regression

Authors: Lucas Kook, Chris Kolb, Philipp Schiele, Daniel Dold, Marcel Arpogaus, Cornelius Fritz, Philipp F. Baumann, Philipp Kopper, Tobias Pielok, Emilio Dorigatti, David Rügamer

Abstract: Neural network representations of simple models, such as linear regression, are being studied increasingly to better understand the underlying principles of deep learning algorithms. However, neural representations of distributional regression models, such as the Cox model, have received little attention so far. We close this gap by proposing a framework for distributional regression using inverse… ▽ More Neural network representations of simple models, such as linear regression, are being studied increasingly to better understand the underlying principles of deep learning algorithms. However, neural representations of distributional regression models, such as the Cox model, have received little attention so far. We close this gap by proposing a framework for distributional regression using inverse flow transformations (DRIFT), which includes neural representations of the aforementioned models. We empirically demonstrate that the neural representations of models in DRIFT can serve as a substitute for their classical statistical counterparts in several applications involving continuous, ordered, time-series, and survival outcomes. We confirm that models in DRIFT empirically match the performance of several statistical methods in terms of estimation of partial effects, prediction, and aleatoric uncertainty quantification. DRIFT covers both interpretable statistical models and flexible neural networks opening up new avenues in both statistical modeling and deep learning. △ Less

Submitted 10 July, 2024; v1 submitted 8 May, 2024; originally announced May 2024.

Comments: Accepted at UAI 2024 https://www.auai.org/uai2024/accepted_papers

Generalizing Orthogonalization for Models with Non-Linearities

Authors: David Rügamer, Chris Kolb, Tobias Weber, Lucas Kook, Thomas Nagler

Abstract: The complexity of black-box algorithms can lead to various challenges, including the introduction of biases. These biases present immediate risks in the algorithms' application. It was, for instance, shown that neural networks can deduce racial information solely from a patient's X-ray scan, a task beyond the capability of medical experts. If this fact is not known to the medical expert, automatic… ▽ More The complexity of black-box algorithms can lead to various challenges, including the introduction of biases. These biases present immediate risks in the algorithms' application. It was, for instance, shown that neural networks can deduce racial information solely from a patient's X-ray scan, a task beyond the capability of medical experts. If this fact is not known to the medical expert, automatic decision-making based on this algorithm could lead to prescribing a treatment (purely) based on racial information. While current methodologies allow for the "orthogonalization" or "normalization" of neural networks with respect to such information, existing approaches are grounded in linear models. Our paper advances the discourse by introducing corrections for non-linearities such as ReLU activations. Our approach also encompasses scalar and tensor-valued predictions, facilitating its integration into neural network architectures. Through extensive experiments, we validate our method's effectiveness in safeguarding sensitive data in generalized linear models, normalizing convolutional neural networks for metadata, and rectifying pre-existing embeddings for undesired attributes. △ Less

Submitted 2 June, 2024; v1 submitted 3 May, 2024; originally announced May 2024.

Input Data Adaptive Learning (IDAL) for Sub-acute Ischemic Stroke Lesion Segmentation

Authors: Michael Götz, Christian Weber, Christoph Kolb, Klaus Maier-Hein

Abstract: In machine learning larger databases are usually associated with higher classification accuracy due to better generalization. This generalization may lead to non-optimal classifiers in some medical applications with highly variable expressions of pathologies. This paper presents a method for learning from a large training base by adaptively selecting optimal training samples for given input data.… ▽ More In machine learning larger databases are usually associated with higher classification accuracy due to better generalization. This generalization may lead to non-optimal classifiers in some medical applications with highly variable expressions of pathologies. This paper presents a method for learning from a large training base by adaptively selecting optimal training samples for given input data. In this way heterogeneous databases are supported two-fold. First, by being able to deal with sparsely annotated data allows a quick inclusion of new data set and second, by training an input-dependent classifier. The proposed approach is evaluated using the SISS challenge. The proposed algorithm leads to a significant improvement of the classification accuracy. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Journal ref: Brainlesion: Glioma, Multiple Sclerosis, Stroke and Traumatic Brain Injuries. BrainLes 2015

Smoothing the Edges: Smooth Optimization for Sparse Regularization using Hadamard Overparametrization

Authors: Chris Kolb, Christian L. Müller, Bernd Bischl, David Rügamer

Abstract: We present a framework for smooth optimization of explicitly regularized objectives for (structured) sparsity. These non-smooth and possibly non-convex problems typically rely on solvers tailored to specific models and regularizers. In contrast, our method enables fully differentiable and approximation-free optimization and is thus compatible with the ubiquitous gradient descent paradigm in deep l… ▽ More We present a framework for smooth optimization of explicitly regularized objectives for (structured) sparsity. These non-smooth and possibly non-convex problems typically rely on solvers tailored to specific models and regularizers. In contrast, our method enables fully differentiable and approximation-free optimization and is thus compatible with the ubiquitous gradient descent paradigm in deep learning. The proposed optimization transfer comprises an overparameterization of selected parameters and a change of penalties. In the overparametrized problem, smooth surrogate regularization induces non-smooth, sparse regularization in the base parametrization. We prove that the surrogate objective is equivalent in the sense that it not only has identical global minima but also matching local minima, thereby avoiding the introduction of spurious solutions. Additionally, our theory establishes results of independent interest regarding matching local minima for arbitrary, potentially unregularized, objectives. We comprehensively review sparsity-inducing parametrizations across different fields that are covered by our general theory, extend their scope, and propose improvements in several aspects. Numerical experiments further demonstrate the correctness and effectiveness of our approach on several sparse learning problems ranging from high-dimensional regression to sparse neural network training. △ Less

Submitted 26 April, 2024; v1 submitted 7 July, 2023; originally announced July 2023.

The Effect of Aspect Changes on Near-Earth Asteroid Phase Curves

Authors: Samuel L. Jackson, Benjamin Rozitis, Lord R. Dover, Simon F. Green, Ulrich C. Kolb, Allison E. Andrews, Stephen C. Lowry

Abstract: Phase curves of asteroids are typically considered to depend solely on the scattering properties of airless particulate surfaces and the size of the object being studied. In this study, we demonstrate the additional dependence of phase curves on object shape, rotation pole orientation, and viewing geometry over an apparition. Variations in the phase curve of near-Earth asteroid (159402) 1999 AP10… ▽ More Phase curves of asteroids are typically considered to depend solely on the scattering properties of airless particulate surfaces and the size of the object being studied. In this study, we demonstrate the additional dependence of phase curves on object shape, rotation pole orientation, and viewing geometry over an apparition. Variations in the phase curve of near-Earth asteroid (159402) 1999 AP10 over its apparition from July 2020 - January 2021 are verified to be due to aspect changes over the apparition. This is achieved through shape modelling of the asteroid and simulation of the phase curve over the apparition. We present simulations of asteroid phase curves over a range of geometries to understand the potential magnitude of this aspect effect, and under which circumstances it can begin to dominate in the phase curves. This dependence on aspect may introduce significant additional uncertainty in the properties derived from phase curve data. We provide and demonstrate software code to estimate the aspect-related uncertainty in near-Earth asteroid phase curves through simulation and model fitting of a randomly generated sample of ellipsoidal asteroid models over the observed viewing geometry. We demonstrate how ignoring this effect may lead to misleading interpretations of the data and underestimation of uncertainties in further studies, such as those in the infrared that use phase curve derived parameters when fitting physical properties of an asteroid. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 14 pages, 11 figures, 5 tables

VeSPA: The SuperWASP Variable Star Photometry Archive

Authors: Adam McMaster, Andrew J. Norton, Hugh J. Dickinson, Heidi B. Thiemann, Ulrich C. Kolb

Abstract: We present the first results from the SuperWASP Variable Stars (SVS) citizen science project. The photometry archive of the Wide Angle Search for Planets has previously been searched for periodic variations and the results of this search formed the basis of the SVS project on the Zooniverse. The SVS project asks volunteers to visually inspect light curve plots and categorize each one according to… ▽ More We present the first results from the SuperWASP Variable Stars (SVS) citizen science project. The photometry archive of the Wide Angle Search for Planets has previously been searched for periodic variations and the results of this search formed the basis of the SVS project on the Zooniverse. The SVS project asks volunteers to visually inspect light curve plots and categorize each one according to a broad classification scheme. Results from the first two years of SVS have now been published online as the SuperWASP Variable Star Photometry Archive (VeSPA). The archive can be browsed online, downloaded in full, or queried, filtered, and sorted to export a refined set of results. An interactive light curve viewer also allows any light curve to be folded at a user-defined period. Analysis of citizen science results and development of VeSPA features are both ongoing. Updated results will be published every six months. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Beep-And-Sleep: Message and Energy Efficient Set Cover

Authors: Thorsten Götte, Christina Kolb, Christian Scheideler, Julian Werthmann

Abstract: We observe message-efficient distributed algorithms for the Set Cover problem. Given a ground set $U$ of $n$ elements and $m$ subsets of $U$, we aim to find the minimal number of these subsets that contain all elements. In the default distributed setup of this problem, each set has a bidirected communication link with each element it contains. Our first result is a $\tilde{O}(\log^2(Δ))$-time and… ▽ More We observe message-efficient distributed algorithms for the Set Cover problem. Given a ground set $U$ of $n$ elements and $m$ subsets of $U$, we aim to find the minimal number of these subsets that contain all elements. In the default distributed setup of this problem, each set has a bidirected communication link with each element it contains. Our first result is a $\tilde{O}(\log^2(Δ))$-time and $O(\sqrt{Δ)}(n+m))$-message algorithm with expected approximation ration of $O(\log(Δ))$ in the $KT_0$ model. The value $Δ$ denotes the maximal cardinality of each subset. Our algorithm is \emph{almost} optimal with regard to time and message complexity. Further, we present Set Cover algorithm in the Beeping model that only relies on carrier-sensing and can trade runtime for approximation ratio similar to the celebrated algorithm by Kuhn and Wattenhofer [PODC '03]. △ Less

Submitted 30 July, 2021; originally announced July 2021.

Asteroid Photometry with PIRATE: Optimizations and Techniques for Small Aperture Telescopes

Authors: Samuel L. Jackson, Ulrich C. Kolb, Simon F. Green

Abstract: Small aperture telescopes provide the opportunity to conduct high frequency, targeted observations of near-Earth Asteroids that are not feasible with larger facilities due to highly competitive time allocation requirements. Observations of asteroids with these types of facilities often focus on rotational brightness variations rather than longer-term phase angle dependent variations (phase curves)… ▽ More Small aperture telescopes provide the opportunity to conduct high frequency, targeted observations of near-Earth Asteroids that are not feasible with larger facilities due to highly competitive time allocation requirements. Observations of asteroids with these types of facilities often focus on rotational brightness variations rather than longer-term phase angle dependent variations (phase curves) due to the difficulty of achieving high precision photometric calibration. We have developed an automated asteroid light curve extraction and calibration pipeline for images of moving objects from the 0.43 m Physics Innovations Robotic Telescope Explorer (PIRATE). This allows for the frequency and quality of observations required to construct asteroid phase curves. Optimisations in standard data reduction procedures are identified that may allow for similar small aperture facilities, constructed from commercially available/off-the-shelf components, to improve image and subsequent data quality. A demonstration of the hardware and software capabilities is expressed through observation statistics from a 10 month observing campaign, and through the photometric characterisation of near-Earth Asteroids 8014 (1990 MF) and 19764 (2000 NF5). △ Less

Submitted 9 July, 2021; originally announced July 2021.

Comments: 41 pages, 15 figures, 2 tables, 1 appendix. Accepted for publication in PASP

Journal ref: PASP, 133, 075003 (2021)

Model-based Joint Analysis of Safety and Security: Survey and Identification of Gaps

Authors: Stefano M. Nicoletti, Marijn Peppelman, Christina Kolb, Mariëlle Stoelinga

Abstract: We survey the state-of-the-art on model-based formalisms for safety and security joint analysis, where safety refers to the absence of unintended failures, and security to absence of malicious attacks. We conduct a thorough literature review and - as a result - we consider fourteen model-based formalisms and compare them with respect to several criteria: (1) Modelling capabilities and Expressivene… ▽ More We survey the state-of-the-art on model-based formalisms for safety and security joint analysis, where safety refers to the absence of unintended failures, and security to absence of malicious attacks. We conduct a thorough literature review and - as a result - we consider fourteen model-based formalisms and compare them with respect to several criteria: (1) Modelling capabilities and Expressiveness: which phenomena can be expressed in these formalisms? To which extent can they capture safety-security interactions? (2) Analytical capabilities: which analysis types are supported? (3) Practical applicability: to what extent have the formalisms been used to analyze small or larger case studies? Furthermore, (1) we present more precise definitions for safety-security dependencies in tree-like formalisms; (2) we showcase the potential of each formalism by modelling the same toy example from the literature and (3) we present our findings and reflect on possible ways to narrow highlighted gaps. In summary, our key findings are the following: (1) the majority of approaches combine tree-like formal models; (2) the exact nature of safety-security interaction is still ill-understood and (3) diverse formalisms can capture different interactions; (4) analyzed formalisms merge modelling constructs from existing safety- and security-specific formalisms, without introducing ad hoc constructs to model safety-security interactions, or (5) metrics to analyze trade offs. Moreover, (6) large case studies representing safety-security interactions are still missing. △ Less

Submitted 23 October, 2023; v1 submitted 11 June, 2021; originally announced June 2021.

deepregression: a Flexible Neural Network Framework for Semi-Structured Deep Distributional Regression

Authors: David Rügamer, Chris Kolb, Cornelius Fritz, Florian Pfisterer, Philipp Kopper, Bernd Bischl, Ruolin Shen, Christina Bukas, Lisa Barros de Andrade e Sousa, Dominik Thalmeier, Philipp Baumann, Lucas Kook, Nadja Klein, Christian L. Müller

Abstract: In this paper we describe the implementation of semi-structured deep distributional regression, a flexible framework to learn conditional distributions based on the combination of additive regression models and deep networks. Our implementation encompasses (1) a modular neural network building system based on the deep learning library \pkg{TensorFlow} for the fusion of various statistical and deep… ▽ More In this paper we describe the implementation of semi-structured deep distributional regression, a flexible framework to learn conditional distributions based on the combination of additive regression models and deep networks. Our implementation encompasses (1) a modular neural network building system based on the deep learning library \pkg{TensorFlow} for the fusion of various statistical and deep learning approaches, (2) an orthogonalization cell to allow for an interpretable combination of different subnetworks, as well as (3) pre-processing steps necessary to set up such models. The software package allows to define models in a user-friendly manner via a formula interface that is inspired by classical statistical model frameworks such as \pkg{mgcv}. The packages' modular design and functionality provides a unique resource for both scalable estimation of complex statistical models and the combination of approaches from deep learning and statistics. This allows for state-of-the-art predictive performance while simultaneously retaining the indispensable interpretability of classical statistical models. △ Less

Submitted 10 March, 2022; v1 submitted 6 April, 2021; originally announced April 2021.

SuperWASP Variable Stars: Classifying Light Curves Using Citizen Science

Authors: Heidi B. Thiemann, Andrew J. Norton, Hugh J. Dickinson, Adam McMaster, Ulrich C. Kolb

Abstract: We present the first analysis of results from the SuperWASP Variable Stars Zooniverse project, which is aiming to classify 1.6 million phase-folded light curves of candidate stellar variables observed by the SuperWASP all sky survey with periods detected in the SuperWASP periodicity catalogue. The resultant data set currently contains $>$1 million classifications corresponding to $>$500,000 object… ▽ More We present the first analysis of results from the SuperWASP Variable Stars Zooniverse project, which is aiming to classify 1.6 million phase-folded light curves of candidate stellar variables observed by the SuperWASP all sky survey with periods detected in the SuperWASP periodicity catalogue. The resultant data set currently contains $>$1 million classifications corresponding to $>$500,000 object-period combinations, provided by citizen scientist volunteers. Volunteer-classified light curves have $\sim$89 per cent accuracy for detached and semi-detached eclipsing binaries, but only $\sim$9 per cent accuracy for rotationally modulated variables, based on known objects. We demonstrate that this Zooniverse project will be valuable for both population studies of individual variable types and the identification of stellar variables for follow up. We present preliminary findings on various unique and extreme variables in this analysis, including long period contact binaries and binaries near the short-period cutoff, and we identify 301 previously unknown binaries and pulsators. We are now in the process of developing a web portal to enable other researchers to access the outputs of the SuperWASP Variable Stars project. △ Less

Submitted 15 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS, 13 pages, 16 figures, 4 tables

WASP-186 and WASP-187: two hot Jupiters discovered by SuperWASP and SOPHIE with additional observations by TESS

Authors: N. Schanche, G. Hébrard, A. Collier Cameron, S. Dalal, B. Smalley, T. G. Wilson, I. Boisse, F. Bouchy, D. J. A Brown, O. Demangeon, C. A. Haswell, C. Hellier, U. C. Kolb, T. Lopez, P. F. L. Maxted, D. L. Pollacco, R. G. West, P. J. Wheatley

Abstract: We present the discovery of two new hot Jupiters identified from the WASP survey, WASP-186b and WASP-187b (TOI-1494.01 and TOI-1493.01). Their planetary nature was established from SOPHIE spectroscopic observations, and additional photometry was obtained from TESS. Stellar parameters for the host stars are derived from spectral line, IRFM, and isochrone placement analyses. These parameters are com… ▽ More We present the discovery of two new hot Jupiters identified from the WASP survey, WASP-186b and WASP-187b (TOI-1494.01 and TOI-1493.01). Their planetary nature was established from SOPHIE spectroscopic observations, and additional photometry was obtained from TESS. Stellar parameters for the host stars are derived from spectral line, IRFM, and isochrone placement analyses. These parameters are combined with the photometric and radial velocity data in an MCMC method to determine the planetary properties. WASP-186b is a massive Jupiter (4.22 +/- 0.18 M_J, 1.11 +/-0.03 R_J) orbiting a mid-F star on a 5.03 day eccentric (e=0.327 +/- 0.008) orbit. WASP-187b is a low density (0.80 +/- 0.09 M_J, 1.64 +/- 0.05 R_J) planet in a 5.15 day circular orbit around a slightly evolved early F-type star. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: 7 pages, 2 appendices

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 499, Issue 1, pp.428-440 (2020)

Investigation of the rotation-activity relation for a sample of SuperWASP and ASAS-SN field stars

Authors: Heidi B. Thiemann, Andrew J. Norton, Ulrich C. Kolb

Abstract: It is well established that late-type main-sequence stars display a relationship between X-ray activity and the Rossby number, $Ro$, the ratio of rotation period to the convective turnover time. This manifests itself as a saturated regime (where X-ray activity is constant) and an unsaturated regime (where X-ray activity anti-correlates with the Rossby number). However, this relationship breaks dow… ▽ More It is well established that late-type main-sequence stars display a relationship between X-ray activity and the Rossby number, $Ro$, the ratio of rotation period to the convective turnover time. This manifests itself as a saturated regime (where X-ray activity is constant) and an unsaturated regime (where X-ray activity anti-correlates with the Rossby number). However, this relationship breaks down for the fastest rotators. We cross-correlated SuperWASP visually classified photometric light curves and ASAS-SN automatically classified photometric light curves with XMM-Newton X-ray observations to identify 3,178 stars displaying a photometrically defined rotational modulation in their light curve and corresponding X-ray observations. We fitted a power-law to characterise the rotation-activity relation of 900 main-sequence stars. We identified that automatically classified rotationally modulated light curves are not as reliable as visually classified light curves for this work. We found a power-law index in the unsaturated regime of G- to M-type stars of $β=-1.84\pm0.18$ for the SuperWASP catalogue, in line with the canonical value of $β=-2$. We find evidence of supersaturation in the fastest rotating K-type stars, with a power-law index of $β_{s}=1.42\pm0.26$. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: Accepted for publication in PASA. 14 pages, 12 figures, 3 tables

Monitoring of transiting exoplanets and their host stars with small aperture telescopes

Authors: M. A. Salisbury, U. C. Kolb, A. J. Norton, C. A. Haswell

Abstract: Exoplanet research is now target rich with a wide diversity of systems making it difficult for high demand observatories to undertake follow up observations over extended periods of time. We investigate the effectiveness of using 0.4m-class telescopes for monitoring transiting hot Jupiters and their host stars. We consider two representative case studies: WASP-52b with 13 new transits, and HAT-P-2… ▽ More Exoplanet research is now target rich with a wide diversity of systems making it difficult for high demand observatories to undertake follow up observations over extended periods of time. We investigate the effectiveness of using 0.4m-class telescopes for monitoring transiting hot Jupiters and their host stars. We consider two representative case studies: WASP-52b with 13 new transits, and HAT-P-23b with 17 new transits and concurrent photometric monitoring covering 78 days. We present updated system parameters and combine our new transit times with previously published results to calculate new ephemerides for both systems. Our analysis of transit mid-times for WASP-52b results in a slight preference for a quadratic ephemeris ($Δχ_ν^2 = 0.07$, $ΔBIC = 1.53$ over a linear ephemeris. We discuss the reality of this quadratic ephemeris indicating a period change of $ δP / δt = -38.6\pm4{\rm ms~yr^{-1}}$ and consider possible causes. WASP-52 is known to be an active star with previous publications reporting many spot crossing events, however no such events are seen in our new photometry. Our analysis shows that WASP-52 is still active and that the latitude of the spots has likely migrated away from the transit chord. We confirm the inflated nature and circular orbit for HAT-P-23b. Our monitoring of HAT-P-23 reveals a periodicity of 7.015 days with an amplitude of 0.011 mag which we interpret as the rotation period of HAT-P-23. The photometric and transit timing precision achieved in the case studies shows that this class of telescope is capable of precise characterisation and long-term monitoring of transiting hot Jupiters in support of dedicated ongoing and future ground and space based observations. △ Less

Submitted 5 September, 2020; v1 submitted 19 August, 2020; originally announced August 2020.

Comments: 15 pages, 9 figures, accepted for publication in New Astronomy Updated to correct typographical error on line 1 of introduction

Discovery and characterization of the exoplanets WASP-148b and c. A transiting system with two interacting giant planets

Authors: G. Hebrard, R. F. Diaz, A. C. M. Correia, A. Collier Cameron, J. Laskar, D. Pollacco, J. -M. Almenara, D. R. Anderson, S. C. C. Barros, I. Boisse, A. S. Bonomo, F. Bouchy, G. Boue, P. Boumis, D. J. A. Brown, S. Dalal, M. Deleuil, O. Demangeon, A. P. Doyle, C. A. Haswell, C. Hellier, H. Osborn, F. Kiefer, U. C. Kolb, K. Lam , et al. (17 additional authors not shown)

Abstract: We present the discovery and characterization of WASP-148, a new extrasolar system that includes at least two giant planets. The host star is a slowly rotating inactive late-G dwarf with a V=12 magnitude. The planet WASP-148b is a hot Jupiter of 0.72 R_Jup and 0.29 M_Jup that transits its host with an orbital period of 8.80 days. We found the planetary candidate with the SuperWASP photometric surv… ▽ More We present the discovery and characterization of WASP-148, a new extrasolar system that includes at least two giant planets. The host star is a slowly rotating inactive late-G dwarf with a V=12 magnitude. The planet WASP-148b is a hot Jupiter of 0.72 R_Jup and 0.29 M_Jup that transits its host with an orbital period of 8.80 days. We found the planetary candidate with the SuperWASP photometric survey, then characterized it with the SOPHIE spectrograph. Our radial velocity measurements subsequently revealed a second planet in the system, WASP-148c, with an orbital period of 34.5 days and a minimum mass of 0.40 M_Jup. No transits of this outer planet were detected. The orbits of both planets are eccentric and fall near the 4:1 mean-motion resonances. This configuration is stable on long timescales, but induces dynamical interactions so that the orbits differ slightly from purely Keplerian orbits. In particular, WASP-148b shows transit-timing variations of typically 15 minutes, making it the first interacting system with transit-timing variations that is detected on ground-based light curves. We establish that the mutual inclination of the orbital plane of the two planets cannot be higher than 35 degrees, and the true mass of WASP-148c is below 0.60 M_Jup. We present photometric and spectroscopic observations of this system that cover a time span of ten years. We also provide their Keplerian and Newtonian analyses; these analyses should be significantly improved through future TESS~observations. △ Less

Submitted 24 June, 2020; v1 submitted 30 April, 2020; originally announced April 2020.

Comments: 16 pages, 10 figures, 7 tables, A&A in press

Journal ref: A&A 640, A32 (2020)

Semi-Structured Distributional Regression -- Extending Structured Additive Models by Arbitrary Deep Neural Networks and Data Modalities

Authors: David Rügamer, Chris Kolb, Nadja Klein

Abstract: Combining additive models and neural networks allows to broaden the scope of statistical regression and extend deep learning-based approaches by interpretable structured additive predictors at the same time. Existing attempts uniting the two modeling approaches are, however, limited to very specific combinations and, more importantly, involve an identifiability issue. As a consequence, interpretab… ▽ More Combining additive models and neural networks allows to broaden the scope of statistical regression and extend deep learning-based approaches by interpretable structured additive predictors at the same time. Existing attempts uniting the two modeling approaches are, however, limited to very specific combinations and, more importantly, involve an identifiability issue. As a consequence, interpretability and stable estimation are typically lost. We propose a general framework to combine structured regression models and deep neural networks into a unifying network architecture. To overcome the inherent identifiability issues between different model parts, we construct an orthogonalization cell that projects the deep neural network into the orthogonal complement of the statistical model predictor. This enables proper estimation of structured model parts and thereby interpretability. We demonstrate the framework's efficacy in numerical experiments and illustrate its special merits in benchmarks and real-world applications. △ Less

Submitted 9 July, 2022; v1 submitted 13 February, 2020; originally announced February 2020.

A Bounding Box Overlay for Competitive Routing in Hybrid Communication Networks

Authors: Jannik Castenow, Christina Kolb, Christian Scheideler

Abstract: In this work, we present a new approach for competitive geometric routing in wireless ad hoc networks. In general, it is well-known that any online routing strategy performs very poor in the worst case. The main difficulty are uncovered regions within the wireless ad hoc network, which we denote as radio holes. Complex shapes of radio holes, for example zig-zag-shapes, make local geometric routing… ▽ More In this work, we present a new approach for competitive geometric routing in wireless ad hoc networks. In general, it is well-known that any online routing strategy performs very poor in the worst case. The main difficulty are uncovered regions within the wireless ad hoc network, which we denote as radio holes. Complex shapes of radio holes, for example zig-zag-shapes, make local geometric routing even more difficult, i.e., forwarded messages in direction to the destination might get stuck in a dead end or are routed along very long detours, when there is no knowledge about the ad hoc network. To obtain knowledge about the position and shape of radio holes, we make use of a hybrid network approach. This approach assumes that we can not just make use of the ad hoc network but also of some cellular infrastructure, which is used to gather knowledge about the underlying ad hoc network. Communication via the cellular infrastructure incurs costs as cell phone providers are involved. Therefore, we use the cellular infrastructure only to compute routing paths in the ad hoc network. The actual data transmission takes place in the ad hoc network. In order to find good routing paths we aim at computing an abstraction of the ad hoc network in which radio holes are abstracted by bounding boxes. The advantage of bounding boxes as hole abstraction is that we only have to consider a constant number of nodes per hole. We prove that bounding boxes are a suitable hole abstraction that allows us to find $c$-competitive paths in the ad hoc network in case of non-intersecting bounding boxes. In case of intersecting bounding boxes, we show via simulations that our routing strategy significantly outperforms the so far best online routing strategies for wireless ad hoc networks. Finally, we also present a routing strategy that is $c$-competitive in case of pairwise intersecting bounding boxes. △ Less

Submitted 10 April, 2019; v1 submitted 12 October, 2018; originally announced October 2018.

Self-stabilizing Overlays for high-dimensional Monotonic Searchability

Authors: Michael Feldmann, Christina Kolb, Christian Scheideler

Abstract: We extend the concept of monotonic searchability for self-stabilizing systems from one to multiple dimensions. A system is self-stabilizing if it can recover to a legitimate state from any initial illegal state. These kind of systems are most often used in distributed applications. Monotonic searchability provides guarantees when searching for nodes while the recovery process is going on. More pre… ▽ More We extend the concept of monotonic searchability for self-stabilizing systems from one to multiple dimensions. A system is self-stabilizing if it can recover to a legitimate state from any initial illegal state. These kind of systems are most often used in distributed applications. Monotonic searchability provides guarantees when searching for nodes while the recovery process is going on. More precisely, if a search request started at some node $u$ succeeds in reaching its destination $v$, then all future search requests from $u$ to $v$ succeed as well. Although there already exists a self-stabilizing protocol for a two-dimensional topology and an universal approach for monotonic searchability, it is not clear how both of these concepts fit together effectively. The latter concept even comes with some restrictive assumptions on messages, which is not the case for our protocol. We propose a simple novel protocol for a self-stabilizing two-dimensional quadtree that satisfies monotonic searchability. Our protocol can easily be extended to higher dimensions and offers routing in $\mathcal O(\log n)$ hops for any search request. △ Less

Submitted 30 August, 2018; originally announced August 2018.

Investigating the Structure of Vela X

Authors: P. Slane, I. Lovchinsky, C. Kolb, S. L. Snowden, T. Temim, J. Blondin, F. Bocchino, M. Miceli, R. A. Chevalier, J. P. Hughes, D. J. Patnaude, T. Gaetz

Abstract: Vela X is the prototypical example of a pulsar wind nebula whose morphology and detailed structure have been affected by the interaction with the reverse shock of its host supernova remnant. The resulting complex of filamentary structure and mixed-in ejecta embedded in a nebula that is offset from the pulsar provides the best example we have of this middle-age state that characterizes a significan… ▽ More Vela X is the prototypical example of a pulsar wind nebula whose morphology and detailed structure have been affected by the interaction with the reverse shock of its host supernova remnant. The resulting complex of filamentary structure and mixed-in ejecta embedded in a nebula that is offset from the pulsar provides the best example we have of this middle-age state that characterizes a significant fraction of composite SNRs, and perhaps all of the large-diameter PWNe seen as TeV sources. Here we report on an XMM-Newton Large Project study of Vela X, supplemented by additional Chandra observations. Through broad spectral modeling as well as detailed spectral investigations of discrete emission regions, we confirm previous studies that report evidence for ejecta material within Vela X, and show that equivalent width variations of O VII and O VIII are consistent with temperature maps within the PWN that show low-temperature regions where the projected SNR emission appears to dominate emission from the ejecta. We identify spectral variations in the nonthermal emission, with hard emission being concentrated near the pulsar. We carry out investigations of the Vela X "cocoon" structure and, with hydrodynamical simulations, show that its overall properties are consistent with structures formed in the late-phase evolution of a composite SNR expanding into a surrounding medium with a density gradient, with ejecta material being swept beyond the pulsar and compressed into an elongated structure in the direction opposite the high external density. △ Less

Submitted 11 August, 2018; originally announced August 2018.

Comments: 13 pages, 13 figures. Accepted for publication in the Astrophysical Journal

Proper Motion of the High-Velocity Pulsar in SNR MSH 15-56

Authors: Tea Temim, Patrick Slane, Paul P. Plucinsky, Joseph Gelfand, Daniel Castro, Christopher Kolb

Abstract: We present a measurement of the proper motion of the presumed pulsar in the evolved composite supernova remnant (SNR) MSH 15-56 whose pulsar wind nebula (PWN) has been disrupted by the supernova (SN) reverse shock. Using Chandra X-ray observations acquired over a baseline of 15 years, we measure a pulsar velocity of 720 (+290/-215) km/s and a direction of motion of 14 +/- 22 degrees west of south.… ▽ More We present a measurement of the proper motion of the presumed pulsar in the evolved composite supernova remnant (SNR) MSH 15-56 whose pulsar wind nebula (PWN) has been disrupted by the supernova (SN) reverse shock. Using Chandra X-ray observations acquired over a baseline of 15 years, we measure a pulsar velocity of 720 (+290/-215) km/s and a direction of motion of 14 +/- 22 degrees west of south. We use this measurement to constrain a hydrodynamical model for the evolution of this system and find that its morphology is well-described by an SNR expanding in an ambient density gradient that increases from east to west. The effect of the density gradient and the pulsar's motion is an asymmetric interaction between the SN reverse shock and the PWN that displaces the bulk of the PWN material away from the pulsar, towards the northeast. The simulation is consistent with an SNR age of 11,000 years, an SN ejecta mass of 10 solar masses, and an average surrounding density of 0.4 cm^-3. However, a combination of a higher SN ejecta mass and ambient density can produce a similar SNR morphology at a later age. △ Less

Submitted 24 November, 2017; originally announced November 2017.

Comments: Accepted to ApJ

Competitive Routing in Hybrid Communication Networks

Authors: Daniel Jung, Christina Kolb, Christian Scheideler, Jannik Sundermeier

Abstract: Routing is a challenging problem for wireless ad hoc networks, especially when the nodes are mobile and spread so widely that in most cases multiple hops are needed to route a message from one node to another. In fact, it is known that any online routing protocol has a poor performance in the worst case, in a sense that there is a distribution of nodes resulting in bad routing paths for that proto… ▽ More Routing is a challenging problem for wireless ad hoc networks, especially when the nodes are mobile and spread so widely that in most cases multiple hops are needed to route a message from one node to another. In fact, it is known that any online routing protocol has a poor performance in the worst case, in a sense that there is a distribution of nodes resulting in bad routing paths for that protocol, even if the nodes know their geographic positions and the geographic position of the destination of a message is known. The reason for that is that radio holes in the ad hoc network may require messages to take long detours in order to get to a destination, which are hard to find in an online fashion. In this paper, we assume that the wireless ad hoc network can make limited use of long-range links provided by a global communication infrastructure like a cellular infrastructure or a satellite in order to compute an abstraction of the wireless ad hoc network that allows the messages to be sent along near-shortest paths in the ad hoc network. We present distributed algorithms that compute an abstraction of the ad hoc network in $\mathcal{O}\left(\log ^2 n\right)$ time using long-range links, which results in $c$-competitive routing paths between any two nodes of the ad hoc network for some constant $c$ if the convex hulls of the radio holes do not intersect. We also show that the storage needed for the abstraction just depends on the number and size of the radio holes in the wireless ad hoc network and is independent on the total number of nodes, and this information just has to be known to a few nodes for the routing to work. △ Less

Submitted 28 February, 2018; v1 submitted 25 October, 2017; originally announced October 2017.

ACM Class: C.2.4

Self-Stabilizing Supervised Publish-Subscribe Systems

Authors: Michael Feldmann, Christina Kolb, Christian Scheideler, Thim Strothmann

Abstract: In this paper we present two major results: First, we introduce the first self-stabilizing version of a supervised overlay network by presenting a self-stabilizing supervised skip ring. Secondly, we show how to use the self-stabilizing supervised skip ring to construct an efficient self-stabilizing publish-subscribe system. That is, in addition to stabilizing the overlay network, every subscriber… ▽ More In this paper we present two major results: First, we introduce the first self-stabilizing version of a supervised overlay network by presenting a self-stabilizing supervised skip ring. Secondly, we show how to use the self-stabilizing supervised skip ring to construct an efficient self-stabilizing publish-subscribe system. That is, in addition to stabilizing the overlay network, every subscriber of a topic will eventually know all of the publications that have been issued so far for that topic. The communication work needed to processes a subscribe or unsubscribe operation is just a constant in a legitimate state, and the communication work of checking whether the system is still in a legitimate state is just a constant on expectation for the supervisor as well as any process in the system. △ Less

Submitted 22 December, 2017; v1 submitted 23 October, 2017; originally announced October 2017.

Evolution of a Pulsar Wind Nebula within a Composite Supernova Remnant

Authors: Christopher Kolb, John Blondin, Patrick Slane, Tea Temim

Abstract: The interaction between a pulsar wind nebula (PWN) and its host supernova remnant (SNR) can produce a vast array of observable structures. Asymmetry present within these structures derives from the complexity of the composite system, where many factors take turns playing a dominating hand throughout the stages of composite SNR evolution. Of particular interest are systems characterized by blastwav… ▽ More The interaction between a pulsar wind nebula (PWN) and its host supernova remnant (SNR) can produce a vast array of observable structures. Asymmetry present within these structures derives from the complexity of the composite system, where many factors take turns playing a dominating hand throughout the stages of composite SNR evolution. Of particular interest are systems characterized by blastwave expansion within a nonuniform interstellar medium (ISM), which contain an active pulsar having a substantial "kick" velocity (upward of 300 km s$^{-1}$ ), because these systems tend to produce complex morphologies. We present a numerical model that employs these and several other factors in an effort to generate asymmetry similar to that seen in various X-ray and radio observations. We find that the main parameters driving structure are ISM uniformity and total pulsar spin- down energy, with secondary contributions from factors such as pulsar trajectory and initial spin-down luminosity. We also investigate the dynamics behind PWN "tails," which may form to link active pulsars to a crushed, relic nebula as the reverse shock passes. We find that the directions of such tails are not good indicators of pulsar motion, but direction does reveal the flow of ejecta created by the passage of a reverse shock. △ Less

Submitted 19 July, 2017; originally announced July 2017.

Comments: 13 pages, 14 figures, accepted to ApJ

The SuperWASP catalogue of 4963 RR Lyr stars: identification of 983 Blazhko effect candidates

Authors: Paul A. Greer, Stan G. Payne, Andrew J. Norton, Pierre F. L. Maxted, Barry Smalley, Richard G. West, Peter J. Wheatley, Ulrich C. Kolb

Abstract: We set out to compile a catalogue of RRab pulsating variables in the SuperWASP archive and identify candidate Blazhko effect objects within this catalogue. We analysed their light curves and power spectra for correlations in their common characteristics to further our understanding of the phenomenon. Pulsation periods were found for each SWASP RRab object using PDM techniques. Low frequency period… ▽ More We set out to compile a catalogue of RRab pulsating variables in the SuperWASP archive and identify candidate Blazhko effect objects within this catalogue. We analysed their light curves and power spectra for correlations in their common characteristics to further our understanding of the phenomenon. Pulsation periods were found for each SWASP RRab object using PDM techniques. Low frequency periodic signals detected in the CLEAN power spectra of RRab stars were matched with modulation sidebands and combined with pairs of sidebands to produce a list of candidate Blazhko periods. A novel technique was used in an attempt to identify Blazhko effect stars by comparing scatter at different parts of the folded light curve. Pulsation amplitudes were calculated based on phase folded light curves. The SuperWASP RRab catalogue consists of 4963 objects of which 3397 are previously unknown. We discovered 983 distinct candidates for Blazhko effect objects, 613 of these being previously unknown in the literature as RR Lyrae stars, and 894 are previously unknown to be Blazhko effect stars. Correlations were investigated between the scatter of points on the light curve, the periods and amplitudes of the objects' pulsations, and those of the Blazhko effect. A statistical analysis has been performed on a large population of Blazhko effect stars from the wide-field SuperWASP survey. No correlations were found between the Blazhko period and other parameters including the Blazhko amplitude, although we confirmed a lower rate of occurrence of the Blazhko effect in long pulsation period objects. △ Less

Submitted 19 July, 2017; v1 submitted 7 July, 2017; originally announced July 2017.

Comments: Includes Note added in Press re: CRTS cross-match

Journal ref: A&A 607, A11 (2017)

Late-time Evolution of Composite Supernova Remnants: Deep Chandra Observations and Hydrodynamical Modeling of a Crushed Pulsar Wind Nebula in SNR G327.1-1.1

Authors: Tea Temim, Patrick Slane, Christopher Kolb, John Blondin, John P. Hughes, Niccolo Bucciantini

Abstract: In an effort to better understand the evolution of composite supernova remnants (SNRs) and the eventual fate of relativistic particles injected by their pulsars, we present a multifaceted investigation of the interaction between a pulsar wind nebula (PWN) and its host SNR G327.1-1.1. Our 350 ks Chandra X-ray observations of SNR G327.1-1.1 reveal a highly complex morphology; a cometary structure re… ▽ More In an effort to better understand the evolution of composite supernova remnants (SNRs) and the eventual fate of relativistic particles injected by their pulsars, we present a multifaceted investigation of the interaction between a pulsar wind nebula (PWN) and its host SNR G327.1-1.1. Our 350 ks Chandra X-ray observations of SNR G327.1-1.1 reveal a highly complex morphology; a cometary structure resembling a bow shock, prong-like features extending into large arcs in the SNR interior, and thermal emission from the SNR shell. Spectral analysis of the non-thermal emission offers clues about the origin of the PWN structures, while enhanced abundances in the PWN region provide evidence for mixing of supernova ejecta with PWN material. The overall morphology and spectral properties of the SNR suggest that the PWN has undergone an asymmetric interaction with the SNR reverse shock (RS) that can occur as a result of a density gradient in the ambient medium and/or a moving pulsar that displaces the PWN from the center of the remnant. We present hydrodynamical simulations of G327.1-1.1 that show that its morphology and evolution can be described by a ~ 17,000 yr old composite SNR that expanded into a density gradient with an orientation perpendicular to the pulsar's motion. We also show that the RS/PWN interaction scenario can reproduce the broadband spectrum of the PWN from radio to gamma-ray wavelengths. The analysis and modeling presented in this work have important implications for our general understanding of the structure and evolution of composite SNRs. △ Less

Submitted 9 June, 2015; originally announced June 2015.

Comments: 13 pages, 3 tables, 8 figures

The doubly eclipsing quintuple low-mass star system 1SWASP J093010.78+533859.5

Authors: M. E. Lohr, A. J. Norton, E. Gillen, R. Busuttil, U. C. Kolb, S. Aigrain, A. McQuillan, S. T. Hodgkin, E. González

Abstract: Our discovery of 1SWASP J093010.78+533859.5 as a probable doubly eclipsing quadruple system containing a contact binary with P~0.23 d and a detached binary with P~1.31 d was announced in 2013. Subsequently Koo et al. confirmed the detached binary spectroscopically and identified a fifth set of static spectral lines at its location, corresponding to a further non-eclipsing component of the system.… ▽ More Our discovery of 1SWASP J093010.78+533859.5 as a probable doubly eclipsing quadruple system containing a contact binary with P~0.23 d and a detached binary with P~1.31 d was announced in 2013. Subsequently Koo et al. confirmed the detached binary spectroscopically and identified a fifth set of static spectral lines at its location, corresponding to a further non-eclipsing component of the system. Here we present new spectroscopic and photometric observations, allowing confirmation of the contact binary and improved modelling of all four eclipsing components. The detached binary is found to contain components of masses 0.837(8) and 0.674(7) M_sol, with radii of 0.832(18) and 0.669(18) R_sol and effective temperatures of 5185(-20,+25) and 4325(-15,+20) K respectively, the contact system has masses 0.86(2) and 0.341(11) M_sol, radii of 0.79(4) and 0.52(5) R_sol respectively, and a common T_eff of 4700(50) K. The fifth star is of similar temperature and spectral type to the primaries in the two binaries. Long-term photometric observations indicate the presence of a spot on one component of the detached binary, moving at an apparent rate of approximately one rotation every two years. Both binaries have consistent system velocities around -11 to -12 km/s, which match the average radial velocity of the fifth star, consistent distance estimates for both subsystems of d=78(3) and d=73(4) pc are also found, and (with some further assumptions) of d=83(9) pc for the fifth star. These findings strongly support the claim that both binaries (and very probably all five stars) are gravitationally bound in a single system. The consistent angles of inclination found for the two binaries (88.2(3) and 86(4) degrees) may also indicate that they originally formed by fragmentation (~9-10 Gyr ago) from a single protostellar disk and subsequently remained in the same orbital plane. △ Less

Submitted 27 April, 2015; originally announced April 2015.

Comments: 13 pages, 23 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 578, A103 (2015)

XMM-Newton and Chandra observations of the M31 globular cluster black hole candidate XB135: a heavyweight contender cut down to size

Authors: R. Barnard, F. Primini, M. R. Garcia, U. C. Kolb, S. S. Murray

Abstract: CXOM31 J004252.030+413107.87 is one of the brightest X-ray sources within the D_25 region of M31, and associated with a globular cluster (GC) known as B135; we therefore call this X-ray source XB135. XB135 is a low mass X-ray binary (LMXB) that apparently exhibited hard state characteristics at 0.3--10 keV luminosities 4--6 E+38 erg/s, and the hard state is only observed below ~10% Eddington. If t… ▽ More CXOM31 J004252.030+413107.87 is one of the brightest X-ray sources within the D_25 region of M31, and associated with a globular cluster (GC) known as B135; we therefore call this X-ray source XB135. XB135 is a low mass X-ray binary (LMXB) that apparently exhibited hard state characteristics at 0.3--10 keV luminosities 4--6 E+38 erg/s, and the hard state is only observed below ~10% Eddington. If true, the accretor would be a high mass black hole (BH) (> ~50 M_Sun); such a BH may be formed from direct collapse of a metal-poor, high mass star, and the very low metalicity of B135 (0.015 Z_Sun) makes such a scenario plausible. We have obtained new XMM-Newton and Chandra HRC observations to shed light on the nature of this object. We find from the HRC observation that XB135 is a single point source located close to the center of B135. The new XMM-Newton spectrum is consistent with a rapidly spinning ~10--20 M_Sun BH in the steep power law or thermal dominant state, but inconsistent with the hard state that we previously assumed. We cannot formally reject three component emission models that have been associated with high luminosity neutron star LMXBs (known as Z-sources); however, we prefer a BH accretor. We note that deeper observation of XB135 could discriminate against a neutron star accretor. △ Less

Submitted 26 January, 2015; originally announced January 2015.

Comments: Accepted for publication in ApJ; 8 pages, 6 figures

Energy-dependent evolution in IC10 X-1: hard evidence for an extended corona, and implications

Authors: R. Barnard, J. F. Steiner, A. F. Prestwich, I. R. Stevens, J. S. Clark, U. C. Kolb

Abstract: We have analyzed a ~130 ks XMM-Newton observation of the dynamically confirmed black hole + Wolf-Rayet (BH+WR) X-ray binary (XB) IC10 X-1, covering ~1 orbital cycle. This system experiences periodic intensity dips every ~35 hours. We find that energy-independent evolution is rejected at a >5 sigma level. The spectral and timing evolution of IC10 X-1 are best explained by a compact disk blackbody a… ▽ More We have analyzed a ~130 ks XMM-Newton observation of the dynamically confirmed black hole + Wolf-Rayet (BH+WR) X-ray binary (XB) IC10 X-1, covering ~1 orbital cycle. This system experiences periodic intensity dips every ~35 hours. We find that energy-independent evolution is rejected at a >5 sigma level. The spectral and timing evolution of IC10 X-1 are best explained by a compact disk blackbody and an extended Comptonized component, where the thermal component is completely absorbed and the Comptonized component is partially covered during the dip. We consider three possibilities for the absorber: cold material in the outer accretion disk, as is well documented for Galactic neutron star (NS) XBs at high inclination; a stream of stellar wind that is enhanced by traveling through the L1 point; and a spherical wind. We estimated the corona radius (r_ADC) for IC10 X-1 from the dip ingress to be ~1 E+6 km, assuming absorption from the outer disk, and found it to be consistent with the relation between r_m ADC and 1--30 keV luminosity observed in Galactic NS XBs that spans 2 orders of magnitude. For the other two scenarios, the corona would be larger. Prior BH mass (M_BH) estimates range over 23--38 M_Sun, depending on the inclination and WR mass. For disk absorption, the inclination, i, is likely to be ~60--80 degrees, with M_BH ~24--41 M_Sun. Alternatively, the L1-enhanced wind requires i ~80 degrees, suggesting ~24--33 M_Sun. For a spherical absorber, i ~ 40 degrees, and M_BH ~50--65 M_Sun. △ Less

Submitted 21 July, 2014; originally announced July 2014.

Comments: Accepted for publication in ApJ; 10 pages, 4 figures

Period and amplitude variations in post-common-envelope eclipsing binaries observed with SuperWASP

Authors: M. E. Lohr, A. J. Norton, D. R. Anderson, A. Collier Cameron, F. Faedi, C. A. Haswell, C. Hellier, S. T. Hodgkin, K. Horne, U. C. Kolb, P. F. L. Maxted, D. Pollacco, I. Skillen, B. Smalley, R. G. West, P. J. Wheatley

Abstract: Period or amplitude variations in eclipsing binaries may reveal the presence of additional massive bodies in the system, such as circumbinary planets. Here, we have studied twelve previously-known eclipsing post-common-envelope binaries for evidence of such light curve variations, on the basis of multi-year observations in the SuperWASP archive. The results for HW Vir provided strong evidence for… ▽ More Period or amplitude variations in eclipsing binaries may reveal the presence of additional massive bodies in the system, such as circumbinary planets. Here, we have studied twelve previously-known eclipsing post-common-envelope binaries for evidence of such light curve variations, on the basis of multi-year observations in the SuperWASP archive. The results for HW Vir provided strong evidence for period changes consistent with those measured by previous studies, and help support a two-planet model for the system. ASAS J102322-3737.0 exhibited plausible evidence for a period increase not previously suggested; while NY Vir, QS Vir and NSVS 14256825 afforded less significant support for period change, providing some confirmation to earlier claims. In other cases, period change was not convincingly observed; for AA Dor and NSVS 07826147, previous findings of constant period were confirmed. This study allows us to present hundreds of new primary eclipse timings for these systems, and further demonstrates the value of wide-field high-cadence surveys like SuperWASP for the investigation of variable stars. △ Less

Submitted 23 May, 2014; originally announced May 2014.

Comments: 11(+19) pages, 4 figures, to be published in Astronomy & Astrophysics. Fig. 1 shown in reduced resolution to minimize file size

Parameters of two low-mass contact eclipsing binaries near the short-period limit

Authors: M. E. Lohr, S. T. Hodgkin, A. J. Norton, U. C. Kolb

Abstract: The two objects 1SWASP J150822.80-054236.9 and 1SWASP J160156.04+202821.6 were initially detected from their SuperWASP archived light curves as candidate eclipsing binaries with periods close to the short-period cut-off of the orbital period distribution of main sequence binaries, at ~0.2 d. Here, using INT spectroscopic data, we confirm them as double-lined spectroscopic and eclipsing binaries, i… ▽ More The two objects 1SWASP J150822.80-054236.9 and 1SWASP J160156.04+202821.6 were initially detected from their SuperWASP archived light curves as candidate eclipsing binaries with periods close to the short-period cut-off of the orbital period distribution of main sequence binaries, at ~0.2 d. Here, using INT spectroscopic data, we confirm them as double-lined spectroscopic and eclipsing binaries, in contact configuration. Following modelling of their visual light curves and radial velocity curves, we determine their component and system parameters to precisions between ~2 and 11%. The former system contains 1.07 and 0.55 M_sun components, with radii of 0.90 and 0.68 R_sun respectively; its primary exhibits pulsations with period 1/6 the orbital period of the system. The latter contains 0.86 and 0.57 M_sun components, with radii of 0.75 and 0.63R_sun respectively. △ Less

Submitted 3 February, 2014; originally announced February 2014.

Comments: 10 pages, 29 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: Astronomy & Astrophysics 563 (2014) A34

One, two or three stars? An investigation of an unusual eclipsing binary candidate undergoing dramatic period changes

Authors: M. E. Lohr, A. J. Norton, U. C. Kolb, D. R. S. Boyd

Abstract: We report our investigation of 1SWASP J234401.81-212229.1, a variable with a 18461.6 s period. After identification in a 2011 search of the SuperWASP archive for main-sequence eclipsing binary candidates near the distribution's short-period limit of approx. 0.20 d, it was measured to be undergoing rapid period decrease in our earlier work, though later observations supported a cyclic variation in… ▽ More We report our investigation of 1SWASP J234401.81-212229.1, a variable with a 18461.6 s period. After identification in a 2011 search of the SuperWASP archive for main-sequence eclipsing binary candidates near the distribution's short-period limit of approx. 0.20 d, it was measured to be undergoing rapid period decrease in our earlier work, though later observations supported a cyclic variation in period length. Spectroscopic data obtained in 2012 with the Southern African Large Telescope did not, however, support the interpretation of the object as a normal eclipsing binary. Here, we consider three possible explanations consistent with the data: a single-star oblique rotator model in which variability results from stable cool spots on opposite magnetic poles; a two-star model in which the secondary is a brown dwarf; and a three-star model involving a low-mass eclipsing binary in a hierarchical triple system. We conclude that the latter is the most likely model. △ Less

Submitted 6 September, 2013; originally announced September 2013.

Comments: 8 pages, 20 figures. Accepted for publication at Astronomy & Astrophysics. After publication, table 2 will appear at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/

Near-UV Absorption, Chromospheric Activity, and Star-Planet Interactions in the WASP-12 system

Authors: C. A. Haswell, L. Fossati, T. Ayres, K. France, C. S. Froning, S. Holmes, U. C. Kolb, R. Busuttil, R. A. Street, L. Hebb, A. Collier Cameron, B. Enoch, V. Burwitz, J. Rodriguez, R. G. West, D. Pollacco, P. J. Wheatley, A. Carter

Abstract: We observed the extreme close-in hot Jupiter system WASP-12 with HST. Near-UV transits up to three times deeper than the optical transit of WASP-12b reveal extensive diffuse gas, extending well beyond the Roche lobe. The distribution of absorbing gas varies between visits. The deepest NUV transits are at wavelength ranges with strong photospheric absorption, implying the absorbing gas may have tem… ▽ More We observed the extreme close-in hot Jupiter system WASP-12 with HST. Near-UV transits up to three times deeper than the optical transit of WASP-12b reveal extensive diffuse gas, extending well beyond the Roche lobe. The distribution of absorbing gas varies between visits. The deepest NUV transits are at wavelength ranges with strong photospheric absorption, implying the absorbing gas may have temperature and composition similar to the stellar photosphere. Our spectra reveal significantly enhanced absorption (greater than 3 σbelow the median) at ~200 wavelengths on each of two HST visits; 65 of these wavelengths are consistent between the two visits, using a strict criterion for velocity matching which excludes matches with velocity shifts exceeding ~20 km/s. Excess transit depths are robustly detected throughout the inner wings of the MgII resonance lines independently on both HST visits. We detected absorption in FeII 2586A, the heaviest species yet detected in an exoplanet transit. The MgII line cores have zero flux, emission cores exhibited by every other observed star of similar age and spectral type are conspicuously absent. WASP-12 probably produces normal MgII profiles, but the inner portions of these strong resonance lines are likely affected by extrinsic absorption. The required Mg+ column is an order of magnitude greater than expected from the ISM, though we cannot completely dismiss that possibility. A more plausible source of absorption is gas lost by WASP-12b. We show that planetary mass loss can produce the required column. Our Visit 2 NUV light curves show evidence for a stellar flare. We show that some of the possible transit detections in resonance lines of rare elements may be due instead to non-resonant transitions in common species. We present optical observations and update the transit ephemeris. △ Less

Submitted 25 January, 2013; v1 submitted 9 January, 2013; originally announced January 2013.

Comments: 79 pages, 23 figures, published in ApJ "replaced" version is identical except for corrected punctuation in the arXiv abstract

Journal ref: 2012, The Astrophysical Journal, 760, 79

Period and period change measurements for 143 SuperWASP eclipsing binary candidates near the short-period limit and discovery of a doubly eclipsing quadruple system

Authors: M. E. Lohr, A. J. Norton, U. C. Kolb, P. F. L. Maxted, I. Todd, R. G. West

Abstract: Building on previous work, a new search of the SuperWASP archive was carried out to identify eclipsing binary systems near the short-period limit. 143 candidate objects were detected with orbital periods between 16000 and 20000 s, of which 97 are new discoveries. Period changes significant at 1 sigma or more were detected in 74 of these objects, and in 38 the changes were significant at 3 sigma or… ▽ More Building on previous work, a new search of the SuperWASP archive was carried out to identify eclipsing binary systems near the short-period limit. 143 candidate objects were detected with orbital periods between 16000 and 20000 s, of which 97 are new discoveries. Period changes significant at 1 sigma or more were detected in 74 of these objects, and in 38 the changes were significant at 3 sigma or more. The significant period changes observed followed an approximately normal distribution with a half-width at half-maximum of ~0.1 s/yr. There was no apparent relationship between period length and magnitude or direction of period change. Amongst several interesting individual objects studied, 1SWASP J093010.78+533859.5 is presented as a new doubly eclipsing quadruple system, consisting of a contact binary with a 19674.575 s period and an Algol-type binary with a 112799.109 s period, separated by 66.1 AU, being the sixth known system of this type. △ Less

Submitted 21 November, 2012; v1 submitted 25 October, 2012; originally announced October 2012.

Comments: Main paper: 11 pages, 11 figures, with data table and 36-page additional figure to appear only in electronic version. Accepted by Astronomy & Astrophysics. Includes minor revisions made following referee's report

Period decrease in three SuperWASP eclipsing binary candidates near the short-period limit

Authors: M. E. Lohr, A. J. Norton, U. C. Kolb, D. R. Anderson, F. Faedi, R. G. West

Abstract: SuperWASP light curves for 53 W UMa-type eclipsing binary (EB) candidates, identified in previous work as being close to the contact binary short-period limit, were studied for evidence of period change. The orbital periods of most of the stars were confirmed, and period decrease, significant at more than 5 sigma, was observed in three objects: 1SWASP J174310.98+432709.6 (-0.055 \pm0.003 s/yr), 1S… ▽ More SuperWASP light curves for 53 W UMa-type eclipsing binary (EB) candidates, identified in previous work as being close to the contact binary short-period limit, were studied for evidence of period change. The orbital periods of most of the stars were confirmed, and period decrease, significant at more than 5 sigma, was observed in three objects: 1SWASP J174310.98+432709.6 (-0.055 \pm0.003 s/yr), 1SWASP J133105.91+121538.0 (-0.075 \pm0.013 s/yr) and 1SWASP J234401.81-212229.1 (-0.313 \pm0.019 s/yr). The magnitudes of the observed period changes cannot be explained by magnetic braking or gravitational radiation effects, and are most likely primarily due to unstable mass transfer from primary to secondary components, possibly accompanied by unstable mass and angular momentum loss from the systems. If these period decreases persist, the systems could merge on a relatively short timescale. △ Less

Submitted 8 May, 2012; originally announced May 2012.

Comments: 6 pages, 7 figures. Accepted for publication in A&A (3 May 2012)

Journal ref: Astronomy & Astrophysics (2012), 542, id.A124

NGC300 X-1 and IC10 X-1: a new breed of black hole binary?

Authors: R. Barnard, J. S. Clark, U. C. Kolb

Abstract: [ABRIDGED] IC10 X-1 has recently been confirmed as a black hole (BH) + Wolf-Rayet (WR) X-ray binary, and NGC300 X-1 is thought to be. IC10 X-1 and NGC300 X-1 have similar X-ray properties, with luminosities ~10^38 erg/s, and orbital periods ~30 hr. We investigate similarities between these two, as well as differences between them and the known Galactic BH binary systems. We have examined XMM-New… ▽ More [ABRIDGED] IC10 X-1 has recently been confirmed as a black hole (BH) + Wolf-Rayet (WR) X-ray binary, and NGC300 X-1 is thought to be. IC10 X-1 and NGC300 X-1 have similar X-ray properties, with luminosities ~10^38 erg/s, and orbital periods ~30 hr. We investigate similarities between these two, as well as differences between them and the known Galactic BH binary systems. We have examined XMM-Newton observations of NGC300 X-1 and IC10 X-1. We extracted lightcurves and spectra; power density spectra (PDS) were constructed from the lightcurves, and the X-ray emission spectra were modeled. Each source exhibits PDS that are characteristic of disc-accreting X-ray binaries (XBs) in the high state. In this state, Galactic XBs with known BH primaries have soft, thermal emission; however the emission spectra of our targets are predominantly non-thermal. Furthermore, the Observation 1 spectrum of NGC300 X-1 is strikingly similar to that of IC10 X-1. The remarkable similarity between the behaviour of NGC300 X-1 in Observation 1 and that of IC10 X-1 lends strong evidence for NGC300 X-1 being a (BH+WR) binary. The unusual spectra of NGC300 X-1 and IC10 X-1 may be due to these systems existing in a persistently high state, whereas all known BH LMXBs are transient. BH XBs in a persistent high state could retain their corona, and hence exhibit a large non-thermal component. LMC X-1 is a BH XB that has only been observed in the high state, and its spectrum is remarkably similar to those of our targets. We therefore classify NGC300 X-1, IC10 X-1 and perhaps LMC X-1 as a new breed of BH XB, defined by their persistently high accretion rates and consequent stable disc configuration and corona. This scenario may also explain the lack of ultraluminous X-ray sources in the canonical soft state. △ Less

Submitted 3 July, 2008; originally announced July 2008.

Comments: Approved for publication in A&A. 8 pages, 5 figures

New XMM-Newton analysis of three bright X-ray sources in M31 globular clusters, including a new black hole candidate

Authors: R. Barnard, H. Stiele, D. Hatzidimitriou, A. K. H. Kong, B. F. Williams, W. Pietsch, U. C. Kolb, F. Haberl, G. Sala

Abstract: We present detailed analysis of three globular cluster X-ray sources in the XMM-Newton extended survey of M31. The X-ray counterpart to the M31 globular cluster Bo 45 (XBo 45) was observed with XMM-Newton on 2006 December 26. Its combined pn+MOS 0.3--10 keV lightcurve exhibited a r.m.s variability of ~10%, and its 0.3--7.0 keV emission spectrum was well described by an absorbed power law with ph… ▽ More We present detailed analysis of three globular cluster X-ray sources in the XMM-Newton extended survey of M31. The X-ray counterpart to the M31 globular cluster Bo 45 (XBo 45) was observed with XMM-Newton on 2006 December 26. Its combined pn+MOS 0.3--10 keV lightcurve exhibited a r.m.s variability of ~10%, and its 0.3--7.0 keV emission spectrum was well described by an absorbed power law with photon index 1.44$\pm$0.12. Its variability and emission is characteristic of low mass X-ray binaries (LMXBs) in the low-hard state, whether the accretor is a neutron star or black hole. Such behaviour is typically observed at luminosities $\la$10% Eddington. However, XBo 45 exhibited this behaviour at an unabsorbed, 0.3--10 keV luminosity of 2.5$\pm0.2\times 10^{38}$ erg s$^{-1}$, or{~140%} Eddington for a 1.4 $M_{\odot}$ neutron star accreting hydrogen. Hence, we identify XBo 45 as a new candidate black hole LMXB. XBo 45 appears to have been consistently bright for ~30 years, consistent with theoretical prediction for a globular cluster black hole binary formed via tidal capture. Bo 375 was observed in the 2007, January 2 XMM-Newton observation, and has a two-component spectrum that is typical for a bright neutron star LMXB. Bo 135 was observed in the same field as Bo 45, and could contain either a black hole or neutron star. △ Less

Submitted 21 August, 2008; v1 submitted 1 July, 2008; originally announced July 2008.

Comments: Accepted by ApJ, 16 pages, 5 figures. This version includes the final changes made at the request of the referee

On the artificial nature of aperiodic variability in XMM-Newton observations of M31 X-ray sources and the ultraluminous X-ray source NGC 4559 ULX-7

Authors: R. Barnard, S. Trudolyubov, U. C. Kolb, C. A. Haswell, J. P. Osborne, W. C. Priedhorsky

Abstract: Context: Power density spectra (PDS) that are characteristic of low mass X-ray binaries (LMXBs) have been previously reported for M31 X-ray sources, observed by XMM-Newton. However, we have recently discovered that these PDS result from the improper addition/subtraction of non-simultaneous lightcurves. Aims: To understand the properties and origins of the artefact. Methods: We re-analysed our XM… ▽ More Context: Power density spectra (PDS) that are characteristic of low mass X-ray binaries (LMXBs) have been previously reported for M31 X-ray sources, observed by XMM-Newton. However, we have recently discovered that these PDS result from the improper addition/subtraction of non-simultaneous lightcurves. Aims: To understand the properties and origins of the artefact. Methods: We re-analysed our XMM-Newton observations of M31 with non-simultaneous and simultaneous lightcurves, then combined simulated lightcurves at various intensities with various offsets and found that the artefact is more dependent on the offset than the intensity. Results: The lightcurves produced by the XMM-Newton Science Analysis Software (SAS) are non-synchronised by default. This affects not only the combination of lightcurves from the three EPIC detectors (MOS1, MOS2 and pn), but also background subtraction in the same CCD. It is therefore imperative that all SAS-generated lightcurves are synchronised by time filtering, even if the whole observation is to be used. We also find that the reported timing behaviour for NGC 4559 ULX-7 was also contaminated by the artefact; there is no significant variability in the correctly-combined lightcurves of NGC 4559 ULX-7. Hence, the classification of this source as an intermediate-mass black hole is no longer justified. Conclusions: While previous timing results from M31 have been proven wrong, and also the broken power law PDS in NGC 4559 ULX-7, XMM-Newton was able to detect aperiodic variability in just 3 ks of observations of NGC 5408 ULX1. Hence XMM-Newton remains a viable tool for analysing variability in extra-galactic X-ray sources. △ Less

Submitted 6 March, 2007; originally announced March 2007.

Comments: Accepted for publication in A&A; 6 pages, 4 figures

Artificial variability in XMM-Newton observations of X-ray sources: M31 as a case study

Authors: R. Barnard, S. Trudolyubov, C. A. Haswell, U. C. Kolb, J. P. Osborne, W. H. Priedhorsky

Abstract: Power density spectra (PDS) that are characteristic of low mass X-ray binaries (LMXBs) have been previously reported for M31 X-ray sources observed by XMM-Newton. However, we have recently discovered that these PDS are false positives resulting from the improper manipulation of non-simultaneous lightcurves. The lightcurves produced by the XMM-Newton Science Analysis Software (SAS) are non-synchr… ▽ More Power density spectra (PDS) that are characteristic of low mass X-ray binaries (LMXBs) have been previously reported for M31 X-ray sources observed by XMM-Newton. However, we have recently discovered that these PDS are false positives resulting from the improper manipulation of non-simultaneous lightcurves. The lightcurves produced by the XMM-Newton Science Analysis Software (SAS) are non-synchronised by default. This affects not only the combination of lightcurves from the three EPIC detectors (MOS1, MOS2 and pn), but also background subtraction in the same CCD. It is therefore imperative that all SAS-generated lightcurves are synchronised by time filtering, even if the whole observation is to be used. We combined simulated lightcurves at various intensities with various offsets and found that the artefact is more dependent on the offset than the intensity. While previous timing results from M31 have been proven wrong, and also the broken power law PDS in NGC 4559 ULX-7, XMM-Newton was able to detect aperiodic variability in just 3 ks of observations of NGC 5408 ULX1. Hence XMM-Newton remains a viable tool for analysing variability in extra-galactic X-ray sources. △ Less

Submitted 2 October, 2006; originally announced October 2006.

Comments: 6 pages, 5 figures; to appear in proceedings for "The multicoloured landscape of compact objects and their explosive progenitors", Cefalu', Sicily, 2006 June 11-24, AIP, in press

Black hole hunting in the Andromeda Galaxy

Authors: R. Barnard, J. P. Osborne, U. C. Kolb, C. A. Haswell

Abstract: We present a new technique for identifying stellar mass black holes in low mass X-ray binaries (LMXBs), and apply it to XMM-Newton observations of M31. We examine X-ray time series variability seeking power density spectra (PDS) typical of LMXBs accreting at a low accretion rate (which we refer to as Type A PDS); these are very similar for black hole and neutron star LMXBs. Galactic neutron star… ▽ More We present a new technique for identifying stellar mass black holes in low mass X-ray binaries (LMXBs), and apply it to XMM-Newton observations of M31. We examine X-ray time series variability seeking power density spectra (PDS) typical of LMXBs accreting at a low accretion rate (which we refer to as Type A PDS); these are very similar for black hole and neutron star LMXBs. Galactic neutron star LMXBs exhibit Type A PDS at low luminosities (~10^36--10^37 erg/s) while black hole LMXBs can exhibit them at luminosities >10^38 erg/s. We propose that Type A PDS are confined to luminosities below a critical fraction of the Eddington limit, $l_c$ that is constant for all LMXBs; we have examined asample of black hole and neutron star LMXBs and find they are all consistent with $l_c$ = 0.10+/-0.04 in the 0.3--10 keV band. We present luminosity and PDS data from 167 observations of X-ray binaries in M31 that provide strong support for our hypothesis. Since the theoretical maximum mass for a neutron star is \~3.1 M_Sun, we therefore assert that any LMXB that exhibits a Type A PDS at a 0.3--10 keV luminosity greater than 4 x 10^37 erg/s is likely to contain a black hole primary. We have found eleven new black hole candidates in M31 using this method. We focus on XMM-Newton observations of RX J0042.4+4112, an X-ray source in M31 and find the mass of the primary to be 7+/-2 M_Sun, if our assumptions are correct. Furthermore, RX J0042.4+4112 is consistently bright in \~40 observations made over 23 years, and is likely to be a persistently bright LMXB; by contrast all known Galactic black hole LMXBs are transient. Hence our method may be used to find black holes in known, persistently bright Galactic LMXBs and also in LMXBs in other galaxies. △ Less

Submitted 6 September, 2004; originally announced September 2004.

Comments: 6 Pages, 6 figures. To appear in the conference proceedings of "Interacting Binaries: Accretion, Evolution and Outcomes" (Cefalu, July 4-10 2004)

The minimum period problem in cataclysmic variables

Authors: J. Barker, U. C. Kolb

Abstract: We investigate if consequential angular momentum losses (CAML) or an intrinsic deformation of the donor star in CVs could increase the CV bounce period from the canonical theoretical value ~65 min to the observed value $P_{min} \approx77$ min, and if a variation of these effects in a CV population could wash out the theoretically predicted accumulation of systems near the minimum period (the per… ▽ More We investigate if consequential angular momentum losses (CAML) or an intrinsic deformation of the donor star in CVs could increase the CV bounce period from the canonical theoretical value ~65 min to the observed value $P_{min} \approx77$ min, and if a variation of these effects in a CV population could wash out the theoretically predicted accumulation of systems near the minimum period (the period spike). We are able to construct suitably mixed CV model populations that a statisticial test cannot rule out as the parent population of the observed CV sample. However, the goodness of fit is never convincing, and always slightly worse than for a simple, flat period distribution. Generally, the goodness of fit is much improved if all CVs are assumed to form at long orbital periods. The weighting suggested by King, Schenker & Hameury (2002) does not constitute an improvment if a realistically shaped input period distribution is used. Put your abstract here. △ Less

Submitted 5 December, 2002; originally announced December 2002.

Comments: 10 pages, Latex, 13 postscript figures, Accepted for publication in MNRAS

XMM-Newton observations of variability in the X-ray binaries of M31: flares dips and a burst

Authors: R. Barnard, U. C. Kolb, J. P. Osborne

Abstract: The XMM-Newton satellite has made four observations of the core of M31 as part of a survey that will span the entire galaxy. The majority of the X-ray emission is contributed by point sources, chiefly X-ray Binaries (XB). Exciting early results from a survey of time variability in three XMM exposures of the M31 core reveal XB exhibiting Sco X-1 like flaring, dips and a burst. The flaring source… ▽ More The XMM-Newton satellite has made four observations of the core of M31 as part of a survey that will span the entire galaxy. The majority of the X-ray emission is contributed by point sources, chiefly X-ray Binaries (XB). Exciting early results from a survey of time variability in three XMM exposures of the M31 core reveal XB exhibiting Sco X-1 like flaring, dips and a burst. The flaring source in M31 also exhibits variability in intensity and spectral hardness that is characteristic of a Z-source. △ Less

Submitted 27 March, 2002; originally announced March 2002.

Comments: 4 pages, 7 figures, to appear in Proc. Symposium ' New Visions of the X-ray Universe in the XMM-Newton and Chandra Era' ESA SP-488 August 2002