-
HyperColorization: Propagating spatially sparse noisy spectral clues for reconstructing hyperspectral images
Authors:
M. Kerem Aydin,
Qi Guo,
Emma Alexander
Abstract:
Hyperspectral cameras face challenging spatial-spectral resolution trade-offs and are more affected by shot noise than RGB photos taken over the same total exposure time. Here, we present a colorization algorithm to reconstruct hyperspectral images from a grayscale guide image and spatially sparse spectral clues. We demonstrate that our algorithm generalizes to varying spectral dimensions for hype…
▽ More
Hyperspectral cameras face challenging spatial-spectral resolution trade-offs and are more affected by shot noise than RGB photos taken over the same total exposure time. Here, we present a colorization algorithm to reconstruct hyperspectral images from a grayscale guide image and spatially sparse spectral clues. We demonstrate that our algorithm generalizes to varying spectral dimensions for hyperspectral images, and show that colorizing in a low-rank space reduces compute time and the impact of shot noise. To enhance robustness, we incorporate guided sampling, edge-aware filtering, and dimensionality estimation techniques. Our method surpasses previous algorithms in various performance metrics, including SSIM, PSNR, GFC, and EMD, which we analyze as metrics for characterizing hyperspectral image quality. Collectively, these findings provide a promising avenue for overcoming the time-space-wavelength resolution trade-off by reconstructing a dense hyperspectral image from samples obtained by whisk or push broom scanners, as well as hybrid spatial-spectral computational imaging systems.
△ Less
Submitted 18 March, 2024;
originally announced March 2024.
-
Understanding Trap States in InP and GaP Quantum Dots Through Density Functional Theory
Authors:
Ezra Alexander,
Matthias Kick,
Alexandra McIsaac,
Troy Van Voorhis
Abstract:
The widespread application of III-V colloidal quantum dots (QDs) as non-toxic, highly tunable emitters is stymied by their high density of trap states. Here, we utilize density functional theory (DFT) to investigate trap state formation in a diverse set of realistically passivated core-only InP and GaP QDs. Through orbital localization techniques, we deconvolute the dense manifold of trap states t…
▽ More
The widespread application of III-V colloidal quantum dots (QDs) as non-toxic, highly tunable emitters is stymied by their high density of trap states. Here, we utilize density functional theory (DFT) to investigate trap state formation in a diverse set of realistically passivated core-only InP and GaP QDs. Through orbital localization techniques, we deconvolute the dense manifold of trap states to allow for detailed assignment of surface defects. We find that the three-coordinate species dominate trapping in III-V QDs and identify features in the geometry and charge environment of trap centers capable of deepening, or sometimes passivating, traps. Furthermore, we observe stark differences in surface reconstruction between InP and GaP, where the more labile InP reconstructs to passivate three-coordinate indium at the cost of distortion elsewhere. These results offer explanations for experimentally observed trapping behavior and suggest new avenues for controlling trap states in III-V QDs.
△ Less
Submitted 26 June, 2024; v1 submitted 1 March, 2024;
originally announced March 2024.
-
PySSED: an automated method of collating and fitting stellar spectral energy distributions
Authors:
Iain McDonald,
Albert A. Zijlstra,
Nick L. J. Cox,
Emma L. Alexander,
Alexander Csukai,
Ria Ramkumar,
Alexander Hollings
Abstract:
Stellar atmosphere modelling predicts the luminosity and temperature of a star, together with parameters such as the effective gravity and the metallicity, by reproducing the observed spectral energy distribution. Most observational data comes from photometric surveys, using a variety of passbands. We herein present the Python Stellar Spectral Energy Distribution (PySSED) routine, designed to comb…
▽ More
Stellar atmosphere modelling predicts the luminosity and temperature of a star, together with parameters such as the effective gravity and the metallicity, by reproducing the observed spectral energy distribution. Most observational data comes from photometric surveys, using a variety of passbands. We herein present the Python Stellar Spectral Energy Distribution (PySSED) routine, designed to combine photometry from disparate catalogues, fit the luminosity and temperature of stars, and determine departures from stellar atmosphere models such as infrared or ultraviolet excess. We detail the routine's operation, and present use cases on both individual stars, stellar populations, and wider regions of the sky. PySSED benefits from fully automated processing, allowing fitting of arbitrarily large datasets at the rate of a few seconds per star.
△ Less
Submitted 19 February, 2024;
originally announced February 2024.
-
Broad Yet Narrow: Super-resolution techniques to simulate electronic spectra of large molecular systems
Authors:
Matthias Kick,
Ezra Alexander,
Anton Beiersdorfer,
Troy Van Voorhis
Abstract:
An accurate treatment of electronic spectra in large systems with a technique such as time dependent density functional theory (TDDFT) is computationally challenging. Due to the Nyquist sampling theorem, direct real time simulations must be prohibitively long in order to recover a suitably sharp resolution in frequency space. Super-resolution techniques such as compressed sensing and MUSIC assume…
▽ More
An accurate treatment of electronic spectra in large systems with a technique such as time dependent density functional theory (TDDFT) is computationally challenging. Due to the Nyquist sampling theorem, direct real time simulations must be prohibitively long in order to recover a suitably sharp resolution in frequency space. Super-resolution techniques such as compressed sensing and MUSIC assume only a small number of excitations contribute to the spectrum, which fails in large molecular systems where the number of excitations is typically very large. We present a new approach that combines exact short time dynamics with approximate frequency space methods to capture large narrow features embedded in a dense manifold of smaller nearby peaks. We show that our approach can accurately capture narrow features and broad quasi-continuum of states at the same time - even when the features overlap in frequency. Our approach is able reduce the required simulation time by a factor of 20-40 with respect to standard Fourier analysis and shows promise for the accurate whole-spectrum prediction of large molecules and materials.
△ Less
Submitted 12 January, 2024;
originally announced January 2024.
-
Hydra II: Characterisation of Aegean, Caesar, ProFound, PyBDSF, and Selavy source finders
Authors:
M. M. Boyce,
A. M. Hopkins,
S. Riggi,
L. Rudnick,
M. Ramsay,
C. L. Hale,
J. Marvil,
M. Whiting,
P. Venkataraman,
C. P. O'Dea,
S. A. Baum,
Y. A. Gordon,
A. N. Vantyghem,
M. Dionyssiou,
H. Andernach,
J. D. Collier,
J. English,
B. S. Koribalski,
D. Leahy,
M. J. Michałowski,
S. Safi-Harb,
M. Vaccari,
E. Alexander,
M. Cowley,
A. D. Kapinska
, et al. (2 additional authors not shown)
Abstract:
We present a comparison between the performance of a selection of source finders using a new software tool called Hydra. The companion paper, Paper~I, introduced the Hydra tool and demonstrated its performance using simulated data. Here we apply Hydra to assess the performance of different source finders by analysing real observational data taken from the Evolutionary Map of the Universe (EMU) Pil…
▽ More
We present a comparison between the performance of a selection of source finders using a new software tool called Hydra. The companion paper, Paper~I, introduced the Hydra tool and demonstrated its performance using simulated data. Here we apply Hydra to assess the performance of different source finders by analysing real observational data taken from the Evolutionary Map of the Universe (EMU) Pilot Survey. EMU is a wide-field radio continuum survey whose primary goal is to make a deep ($20μ$Jy/beam RMS noise), intermediate angular resolution ($15^{\prime\prime}$), 1\,GHz survey of the entire sky south of $+30^{\circ}$ declination, and expecting to detect and catalogue up to 40 million sources. With the main EMU survey expected to begin in 2022 it is highly desirable to understand the performance of radio image source finder software and to identify an approach that optimises source detection capabilities. Hydra has been developed to refine this process, as well as to deliver a range of metrics and source finding data products from multiple source finders. We present the performance of the five source finders tested here in terms of their completeness and reliability statistics, their flux density and source size measurements, and an exploration of case studies to highlight finder-specific limitations.
△ Less
Submitted 27 April, 2023;
originally announced April 2023.
-
Hydra I: An extensible multi-source-finder comparison and cataloguing tool
Authors:
M. M. Boyce,
A. M. Hopkins,
S. Riggi,
L. Rudnick,
M. Ramsay,
C. L. Hale,
J. Marvil,
M. Whiting,
P. Venkataraman,
C. P. O'Dea,
S. A. Baum,
Y. A. Gordon,
A. N. Vantyghem,
M. Dionyssiou,
H. Andernach,
J. D. Collier,
J. English,
B. S. Koribalski,
D. Leahy,
M. J. Michałowski,
S. Safi-Harb,
M. Vaccari,
E. Alexander,
M. Cowley,
A. D. Kapinska
, et al. (2 additional authors not shown)
Abstract:
The latest generation of radio surveys are now producing sky survey images containing many millions of radio sources. In this context it is highly desirable to understand the performance of radio image source finder (SF) software and to identify an approach that optimises source detection capabilities. We have created Hydra to be an extensible multi-SF and cataloguing tool that can be used to comp…
▽ More
The latest generation of radio surveys are now producing sky survey images containing many millions of radio sources. In this context it is highly desirable to understand the performance of radio image source finder (SF) software and to identify an approach that optimises source detection capabilities. We have created Hydra to be an extensible multi-SF and cataloguing tool that can be used to compare and evaluate different SFs. Hydra, which currently includes the SFs Aegean, Caesar, ProFound, PyBDSF, and Selavy, provides for the addition of new SFs through containerisation and configuration files. The SF input RMS noise and island parameters are optimised to a 90\% ''percentage real detections'' threshold (calculated from the difference between detections in the real and inverted images), to enable comparison between SFs. Hydra provides completeness and reliability diagnostics through observed-deep ($\mathcal{D}$) and generated-shallow ($\mathcal{S}$) images, as well as other statistics. In addition, it has a visual inspection tool for comparing residual images through various selection filters, such as S/N bins in completeness or reliability. The tool allows the user to easily compare and evaluate different SFs in order to choose their desired SF, or a combination thereof. This paper is part one of a two part series. In this paper we introduce the Hydra software suite and validate its $\mathcal{D/S}$ metrics using simulated data. The companion paper demonstrates the utility of Hydra by comparing the performance of SFs using both simulated and real images.
△ Less
Submitted 27 April, 2023;
originally announced April 2023.
-
Radio Galaxy Zoo EMU: Towards a Semantic Radio Galaxy Morphology Taxonomy
Authors:
Micah Bowles,
Hongming Tang,
Eleni Vardoulaki,
Emma L. Alexander,
Yan Luo,
Lawrence Rudnick,
Mike Walmsley,
Fiona Porter,
Anna M. M. Scaife,
Inigo Val Slijepcevic,
Elizabeth A. K. Adams,
Alexander Drabent,
Thomas Dugdale,
Gülay Gürkan,
Andrew M. Hopkins,
Eric F. Jimenez-Andrade,
Denis A. Leahy,
Ray P. Norris,
Syed Faisal ur Rahman,
Xichang Ouyang,
Gary Segal,
Stanislav S. Shabala,
O. Ivy Wong
Abstract:
We present a novel natural language processing (NLP) approach to deriving plain English descriptors for science cases otherwise restricted by obfuscating technical terminology. We address the limitations of common radio galaxy morphology classifications by applying this approach. We experimentally derive a set of semantic tags for the Radio Galaxy Zoo EMU (Evolutionary Map of the Universe) project…
▽ More
We present a novel natural language processing (NLP) approach to deriving plain English descriptors for science cases otherwise restricted by obfuscating technical terminology. We address the limitations of common radio galaxy morphology classifications by applying this approach. We experimentally derive a set of semantic tags for the Radio Galaxy Zoo EMU (Evolutionary Map of the Universe) project and the wider astronomical community. We collect 8,486 plain English annotations of radio galaxy morphology, from which we derive a taxonomy of tags. The tags are plain English. The result is an extensible framework which is more flexible, more easily communicated, and more sensitive to rare feature combinations which are indescribable using the current framework of radio astronomy classifications.
△ Less
Submitted 14 April, 2023;
originally announced April 2023.
-
Thermal Spread Functions (TSF): Physics-guided Material Classification
Authors:
Aniket Dashpute,
Vishwanath Saragadam,
Emma Alexander,
Florian Willomitzer,
Aggelos Katsaggelos,
Ashok Veeraraghavan,
Oliver Cossairt
Abstract:
Robust and non-destructive material classification is a challenging but crucial first-step in numerous vision applications. We propose a physics-guided material classification framework that relies on thermal properties of the object. Our key observation is that the rate of heating and cooling of an object depends on the unique intrinsic properties of the material, namely the emissivity and diffus…
▽ More
Robust and non-destructive material classification is a challenging but crucial first-step in numerous vision applications. We propose a physics-guided material classification framework that relies on thermal properties of the object. Our key observation is that the rate of heating and cooling of an object depends on the unique intrinsic properties of the material, namely the emissivity and diffusivity. We leverage this observation by gently heating the objects in the scene with a low-power laser for a fixed duration and then turning it off, while a thermal camera captures measurements during the heating and cooling process. We then take this spatial and temporal "thermal spread function" (TSF) to solve an inverse heat equation using the finite-differences approach, resulting in a spatially varying estimate of diffusivity and emissivity. These tuples are then used to train a classifier that produces a fine-grained material label at each spatial pixel. Our approach is extremely simple requiring only a small light source (low power laser) and a thermal camera, and produces robust classification results with 86% accuracy over 16 classes.
△ Less
Submitted 2 April, 2023;
originally announced April 2023.
-
Galaxy Image Deconvolution for Weak Gravitational Lensing with Unrolled Plug-and-Play ADMM
Authors:
Tianao Li,
Emma Alexander
Abstract:
Removing optical and atmospheric blur from galaxy images significantly improves galaxy shape measurements for weak gravitational lensing and galaxy evolution studies. This ill-posed linear inverse problem is usually solved with deconvolution algorithms enhanced by regularisation priors or deep learning. We introduce a so-called "physics-informed deep learning" approach to the Point Spread Function…
▽ More
Removing optical and atmospheric blur from galaxy images significantly improves galaxy shape measurements for weak gravitational lensing and galaxy evolution studies. This ill-posed linear inverse problem is usually solved with deconvolution algorithms enhanced by regularisation priors or deep learning. We introduce a so-called "physics-informed deep learning" approach to the Point Spread Function (PSF) deconvolution problem in galaxy surveys. We apply algorithm unrolling and the Plug-and-Play technique to the Alternating Direction Method of Multipliers (ADMM), in which a neural network learns appropriate hyperparameters and denoising priors from simulated galaxy images. We characterise the time-performance trade-off of several methods for galaxies of differing brightness levels as well as our method's robustness to systematic PSF errors and network ablations. We show an improvement in reduced shear ellipticity error of 38.6% (SNR=20)/45.0% (SNR=200) compared to classic methods and 7.4% (SNR=20)/33.2% (SNR=200) compared to modern methods.
△ Less
Submitted 13 March, 2023; v1 submitted 2 November, 2022;
originally announced November 2022.
-
A New Task: Deriving Semantic Class Targets for the Physical Sciences
Authors:
Micah Bowles,
Hongming Tang,
Eleni Vardoulaki,
Emma L. Alexander,
Yan Luo,
Lawrence Rudnick,
Mike Walmsley,
Fiona Porter,
Anna M. M. Scaife,
Inigo Val Slijepcevic,
Gary Segal
Abstract:
We define deriving semantic class targets as a novel multi-modal task. By doing so, we aim to improve classification schemes in the physical sciences which can be severely abstracted and obfuscating. We address this task for upcoming radio astronomy surveys and present the derived semantic radio galaxy morphology class targets.
We define deriving semantic class targets as a novel multi-modal task. By doing so, we aim to improve classification schemes in the physical sciences which can be severely abstracted and obfuscating. We address this task for upcoming radio astronomy surveys and present the derived semantic radio galaxy morphology class targets.
△ Less
Submitted 27 October, 2022; v1 submitted 26 October, 2022;
originally announced October 2022.
-
Word Clouds in the Wild
Authors:
Rebecca M. M. Hicke,
Maanya Goenka,
Eric Alexander
Abstract:
Word clouds are frequently used to analyze and communicate text data in many domains. In order to help guide research on improving the legibility of word clouds, we have conducted a survey of their usage in Digital Humanities academia and journalism. Using a modified grounded theory approach, we sought to identify the most common purposes for which word clouds were employed and the most common vis…
▽ More
Word clouds are frequently used to analyze and communicate text data in many domains. In order to help guide research on improving the legibility of word clouds, we have conducted a survey of their usage in Digital Humanities academia and journalism. Using a modified grounded theory approach, we sought to identify the most common purposes for which word clouds were employed and the most common visual encodings they contained. Our findings indicate that font size, color, and word placement dominate as the primary data-encoding channels, as we hypothesized. Perhaps more surprisingly, we found that asking viewers to perform analytical tasks with word clouds was relatively common, especially in DH sources. This suggests that research into the interactions of these visual encoding channels (particularly in regards to legibility) is warranted.
△ Less
Submitted 14 October, 2022;
originally announced October 2022.
-
Identifying anomalous radio sources in the EMU Pilot Survey using a complexity-based approach
Authors:
Gary Segal,
David Parkinson,
Ray Norris,
Andrew M. Hopkins,
Heinz Andernach,
Emma L. Alexander,
Ettore Carretti,
Bärbel S. Koribalski,
Letjatji S. Legodi,
Sarah Leslie,
Yan Luo,
Jonathon C. S. Pierce,
Hongming Tang,
Eleni Vardoulaki,
Tessa Vernstrom
Abstract:
The Evolutionary Map of the Universe (EMU) large-area radio continuum survey will detect tens of millions of radio galaxies, giving an opportunity for the detection of previously unknown classes of objects. To maximise the scientific value and make new discoveries, the analysis of this data will need to go beyond simple visual inspection. We propose the coarse-grained complexity, a simple scalar q…
▽ More
The Evolutionary Map of the Universe (EMU) large-area radio continuum survey will detect tens of millions of radio galaxies, giving an opportunity for the detection of previously unknown classes of objects. To maximise the scientific value and make new discoveries, the analysis of this data will need to go beyond simple visual inspection. We propose the coarse-grained complexity, a simple scalar quantity relating to the minimum description length of an image, that can be used to identify unusual structures. The complexity can be computed without reference to the broader sample or existing catalogue data, making the computation efficient on new surveys at very large scales (such as the full EMU survey). We apply our coarse-grained complexity measure to data from the EMU Pilot Survey to detect and confirm anomalous objects in this data set and produce an anomaly catalogue. Rather than work with existing catalogue data using a specific source detection algorithm, we perform a blind scan of the area, computing the complexity using a sliding square aperture. The effectiveness of the complexity measure for identifying anomalous objects is evaluated using crowd-sourced labels generated via the Zooniverse.org platform. We find that the complexity scan identifies unusual sources, such as odd radio circles, by partitioning on complexity. We achieve partitions where 5\% of the data is estimated to be 86\% complete, and 0.5\% is estimated to be 94\% pure, with respect to anomalies and use this to produce an anomaly catalogue.
△ Less
Submitted 3 March, 2023; v1 submitted 29 June, 2022;
originally announced June 2022.
-
CS-ROMER: A novel compressed sensing framework for Faraday depth reconstruction
Authors:
Miguel Cárcamo,
Anna M. M. Scaife,
Emma L. Alexander,
J. Patrick Leahy
Abstract:
The reconstruction of Faraday depth structure from incomplete spectral polarization radio measurements using the RM Synthesis technique is an under-constrained problem requiring additional regularisation. In this paper we present cs-romer: a novel object-oriented compressed sensing framework to reconstruct Faraday depth signals from spectro-polarization radio data. Unlike previous compressed sensi…
▽ More
The reconstruction of Faraday depth structure from incomplete spectral polarization radio measurements using the RM Synthesis technique is an under-constrained problem requiring additional regularisation. In this paper we present cs-romer: a novel object-oriented compressed sensing framework to reconstruct Faraday depth signals from spectro-polarization radio data. Unlike previous compressed sensing applications, this framework is designed to work directly with data that are irregularly sampled in wavelength-squared space and to incorporate multiple forms of compressed sensing regularisation. We demonstrate the framework using simulated data for the VLA telescope under a variety of observing conditions, and we introduce a methodology for identifying the optimal basis function for reconstruction of these data, using an approach that can also be applied to datasets from other telescopes and over different frequency ranges. In this work we show that the delta basis function provides optimal reconstruction for VLA L-band data and we use this basis with observations of the low-mass galaxy cluster Abell 1314 in order to reconstruct the Faraday depth of its constituent cluster galaxies. We use the cs-romer framework to de-rotate the Galactic Faraday depth contribution directly from the wavelength-squared data and to handle the spectral behaviour of different radio sources in a direction-dependent manner. The results of this analysis show that individual galaxies within Abell 1314 deviate from the behaviour expected for a Faraday-thin screen such as the intra-cluster medium and instead suggest that the Faraday rotation exhibited by these galaxies is dominated by their local environments.
△ Less
Submitted 26 October, 2022; v1 submitted 3 May, 2022;
originally announced May 2022.
-
Isoperimetric 3- and 4-bubble results on $\mathbb{R}$ with density $|x|$
Authors:
Evan Alexander,
Emily Burns,
John Ross,
Jesse Stovall,
Zariah Whyte
Abstract:
We study the isoperimetric problem on $\mathbb{R}^1$ with a prescribed density function $f(x) = |x|$. Under these conditions, we find that isoperimetric $3$-bubble and $4$-bubble results satisfy a regular structure. As our regions increase in size, the intervals that form them alternate back-and-forth across the origin, with the smaller regions closer to the origin. This expands on previously know…
▽ More
We study the isoperimetric problem on $\mathbb{R}^1$ with a prescribed density function $f(x) = |x|$. Under these conditions, we find that isoperimetric $3$-bubble and $4$-bubble results satisfy a regular structure. As our regions increase in size, the intervals that form them alternate back-and-forth across the origin, with the smaller regions closer to the origin. This expands on previously known observations about the single- and double-bubble results on $\mathbb{R}$ with density $|x|^p$.
△ Less
Submitted 5 January, 2022;
originally announced January 2022.
-
Automated Multi-Process CTC Detection using Deep Learning
Authors:
Elena Alexander,
Kam W. Leong,
Andrew F. Laine
Abstract:
Circulating Tumor Cells (CTCs) bear great promise as biomarkers in tumor prognosis. However, the process of identification and later enumeration of CTCs require manual labor, which is error-prone and time-consuming. The recent developments in object detection via Deep Learning using Mask-RCNNs and wider availability of pre-trained models have enabled sensitive tasks with limited data of such to be…
▽ More
Circulating Tumor Cells (CTCs) bear great promise as biomarkers in tumor prognosis. However, the process of identification and later enumeration of CTCs require manual labor, which is error-prone and time-consuming. The recent developments in object detection via Deep Learning using Mask-RCNNs and wider availability of pre-trained models have enabled sensitive tasks with limited data of such to be tackled with unprecedented accuracy. In this report, we present a novel 3-stage detection model for automated identification of Circulating Tumor Cells in multi-channel darkfield microscopic images comprised of: RetinaNet based identification of Cytokeratin (CK) stains, Mask-RCNN based cell detection of DAPI cell nuclei and Otsu thresholding to detect CD-45s. The training dataset is composed of 46 high variance data points, with 10 Negative and 36 Positive data points. The test set is composed of 420 negative data points. The final accuracy of the pipeline is 98.81%.
△ Less
Submitted 26 September, 2021;
originally announced September 2021.
-
Prediction of the Existence of LiCH, a Carbene-like Organometallic Molecule
Authors:
Jason M. Montgomery,
Ezra Alexander,
David A. Mazziotti
Abstract:
Carbenes comprise a well-known class of organometallic compounds consisting of a neutral, divalent carbon and two unshared electrons. Carbenes can have singlet or triplet ground states, each giving rise to a distinct reactivity. Methylene, CH$_2$, the parent hydride, is well-known to be bent in its triplet ground state. Here we predict the existence of LiCH, a carbene-like organometallic molecule.…
▽ More
Carbenes comprise a well-known class of organometallic compounds consisting of a neutral, divalent carbon and two unshared electrons. Carbenes can have singlet or triplet ground states, each giving rise to a distinct reactivity. Methylene, CH$_2$, the parent hydride, is well-known to be bent in its triplet ground state. Here we predict the existence of LiCH, a carbene-like organometallic molecule. Computationally, we treat the electronic structure with parametric and variational two-electron reduced density matrix (2-RDM) methods, which are capable of capturing multireference correlation typically associated with the singlet state of a diradical. Similar to methylene, LiCH is a triplet ground state with a predicted 15.8 kcal/mol singlet-triplet gap. However, unlike methylene, LiCH is linear in both the triplet state and the lowest excited singlet state. Furthermore, the singlet state is found to exhibit strong electron correlation as a diradical. In comparison to dissociation channels Li + CH and Li$^+$ + CH$^-$, the LiCH was found to be stable by approximately 77 kcal/mol.
△ Less
Submitted 4 August, 2020;
originally announced August 2020.
-
The ASKAP-EMU Early Science Project:Radio Continuum Survey of the Small Magellanic Cloud
Authors:
T. D. Joseph,
M. D. Filipović,
E. J. Crawford,
I. Bojičić,
E. L. Alexander,
G. F. Wong,
H. Andernach,
H. Leverenz,
R. P. Norris,
R. Z. E. Alsaberi,
C. Anderson,
L. A. Barnes,
L. M. Bozzetto,
F. Bufano,
J. D. Bunton,
F. Cavallaro,
J. D. Collier,
H. Dénes,
Y. Fukui,
T. Galvin,
F. Haberl,
A. Ingallinera,
A. D. Kapinska,
B. S. Koribalski,
R. Kothes
, et al. (31 additional authors not shown)
Abstract:
We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observ…
▽ More
We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0"x30.0" and 16.3"x15.1", respectively. The median Root Mean Squared (RMS) noise values are 186$μ$Jy beam$^{-1}$ (960 MHz) and 165$μ$Jy beam$^{-1}$ (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg$^2$. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies.
△ Less
Submitted 23 September, 2019;
originally announced September 2019.
-
A new functional RG flow: regulator-sourced 2PI versus average 1PI
Authors:
Elizabeth Alexander,
Peter Millington,
Jordan Nursey,
Paul M. Saffin
Abstract:
We derive the renormalization group evolution of the quartic scalar theory with spontaneous symmetry breaking from an alternative flow equation, obtained within the externally sourced two-particle irreducible framework due to Garbrecht and Millington. In order to make a straightforward comparison with the evolution from the standard Wetterich-Morris-Ellwanger equation, we employ the Litim regulato…
▽ More
We derive the renormalization group evolution of the quartic scalar theory with spontaneous symmetry breaking from an alternative flow equation, obtained within the externally sourced two-particle irreducible framework due to Garbrecht and Millington. In order to make a straightforward comparison with the evolution from the standard Wetterich-Morris-Ellwanger equation, we employ the Litim regulator, work to lowest order in the derivative expansion and neglect anomalous scaling. By this means, we illustrate the leading differences between analytic expressions for the resulting threshold and (non-perturbative) beta functions. In four dimensions, we find that the positions of the potential minima and the cosmological constant evolve more rapidly with scale compared to the standard approach, whereas the quartic coupling evolves more slowly, albeit by a small amount. These differences may have implications for the asymptotic safety programme, as well as our understanding of the non-perturbative scale evolution of the Standard Model Higgs sector.
△ Less
Submitted 6 August, 2019;
originally announced August 2019.
-
Alternative flow equation for the functional renormalization group
Authors:
Elizabeth Alexander,
Peter Millington,
Jordan Nursey,
Paul M. Saffin
Abstract:
We derive an alternative to the Wetterich-Morris-Ellwanger equation by means of the two-particle irreducible (2PI) effective action, exploiting the method of external sources due to Garbrecht and Millington. The latter allows the two-point source of the 2PI effective action to be associated consistently with the regulator of the renormalization group flow. We show that this procedure leads to a fl…
▽ More
We derive an alternative to the Wetterich-Morris-Ellwanger equation by means of the two-particle irreducible (2PI) effective action, exploiting the method of external sources due to Garbrecht and Millington. The latter allows the two-point source of the 2PI effective action to be associated consistently with the regulator of the renormalization group flow. We show that this procedure leads to a flow equation that differs from that obtained in the standard approach based on the average one-particle irreducible effective action.
△ Less
Submitted 28 July, 2021; v1 submitted 15 July, 2019;
originally announced July 2019.
-
Spectroscopy with the Engineering Development Array: cold H$^{+}$ at 63 MHz towards the Galactic Centre
Authors:
J. B. R. Oonk,
E. L. Alexander,
J. W. Broderick,
M. Sokolowski,
R. Wayth
Abstract:
The Engineering Development Array (EDA) is a single test station for Square Kilometre Array (SKA) precursor technology. We have used the EDA to detect low-frequency radio recombination lines (RRLs) from the Galactic Centre region. Low-frequency RRLs are an area of interest for future low-frequency SKA work as these lines provide important information on the physical properties of the cold neutral…
▽ More
The Engineering Development Array (EDA) is a single test station for Square Kilometre Array (SKA) precursor technology. We have used the EDA to detect low-frequency radio recombination lines (RRLs) from the Galactic Centre region. Low-frequency RRLs are an area of interest for future low-frequency SKA work as these lines provide important information on the physical properties of the cold neutral medium. In this project we investigate the EDA, its bandpass and the radio frequency interference environment for low-frequency spectroscopy. We present line spectra from 30 to 325 MHz for the Galactic Centre region. The decrease in sensitivity for the EDA at the low end of the receiver prevents carbon and hydrogen RRLs to be detected below 40 and 60 MHz respectively. RFI strongly affects frequencies in the range 276-292, 234-270, 131-138, 95-102 and below 33 MHz. Cn$α$ RRLs were detected in absorption for quantum levels n = 378 to 550 (39-121 MHz) and in emission for n = 272 to 306 (228-325 MHz). Cn$β$ lines were detected in absorption for n = 387 to 696 (39-225 MHz). Hn$α$ RRLs were detected in emission for n = 272 to 480 (59-325 MHz). Hn$β$ lines were detected for n = 387 to 453 (141-225 MHz). The stacked Hn$α$ detection at 63 MHz is the lowest frequency detection made for hydrogen RRLs and shows that a cold (partially) ionized medium exists along the line of sight to the Galactic Centre region. The size and velocity of this cold H$^{+}$ gas indicates that it is likely associated with the nearby Riegel-Crutcher cloud.
△ Less
Submitted 6 July, 2019;
originally announced July 2019.
-
Prediction of a red nova outburst in KIC 9832227
Authors:
Lawrence A. Molnar,
Daniel M. Van Noord,
Karen Kinemuchi,
Jason P. Smolinski,
Cara E. Alexander,
Evan M. Cook,
Byoungchan Jang,
Henry A. Kobulnicky,
Christopher J. Spedden,
Steven D. Steenwyk
Abstract:
We present the first identification of a candidate precursor for an imminent red nova. Our prediction is based on the example of the precursor to the red nova V1309 Sco, which was retrospectively found to be a contact binary with an exponentially decreasing period. We explore the use of this distinctive timing signature to identify precursors, developing the observational and analysis steps needed…
▽ More
We present the first identification of a candidate precursor for an imminent red nova. Our prediction is based on the example of the precursor to the red nova V1309 Sco, which was retrospectively found to be a contact binary with an exponentially decreasing period. We explore the use of this distinctive timing signature to identify precursors, developing the observational and analysis steps needed. We estimate that our Galaxy has roughly 1-10 observable precursors. Specifically, we lay out the observational case for KIC 9832227, which we identified as a tentative candidate two years ago (Molnar et al. 2015, AAS Meeting Abstracts 415.05). Orbital timing over the past two years has followed the tentative exponential fit. As of late 2015, the period time derivative went beyond the range found in other systems (dP/dt < 1x10^{-8}), a necessary criterion for a serious candidate. We estimate time of merger is the year 2022.2(7). Double absorption line spectra confirm directly the 0.458 d light curve period is a contact binary system and yield a mass ratio m_B/m_A = 0.228(3). Closer analysis of the Kepler timing data shows evidence of a component C with orbital period P_C = 590(8) days and m_C x sin i_C = 0.11 solar masses. An alternative interpretation of the long term timing trend, light travel time delay due to orbit around a distant component D, is ruled out by the spectroscopic data for any nondegenerate star. Additional measurements are needed to test further the merging hypothesis and to utilize fully this fortuitous opportunity.
△ Less
Submitted 18 April, 2017;
originally announced April 2017.
-
Can Large Time Delays Observed in Light Curves of Coronal Loops be Explained by Impulsive Heating?
Authors:
Roberto Lionello,
Caroline E. Alexander,
Amy R. Winebarger,
Jon A. Linker,
Zoran Mikić
Abstract:
The light curves of solar coronal loops often peak first in channels associated with higher temperatures and then in those associated with lower. The time delays between the different narrowband EUV channels have been measured for many individual loops and recently for every pixel of an active region observation. Time delays between channels for an active region exhibit a wide range of values, wit…
▽ More
The light curves of solar coronal loops often peak first in channels associated with higher temperatures and then in those associated with lower. The time delays between the different narrowband EUV channels have been measured for many individual loops and recently for every pixel of an active region observation. Time delays between channels for an active region exhibit a wide range of values, with maxima $>$ 5,000\,s. These large time delays make up 3-26\% (depending on the channel pair) of the pixels where a significant, positive time delay is measured. It has been suggested that time delays can be explained by impulsive heating. In this paper, we investigate whether the largest observed time delays can be explained by this hypothesis by simulating a series of coronal loops with different heating rates, loop lengths, abundances, and geometries to determine the range of expected time delays between a set of four EUV channels. We find that impulsive heating cannot address the largest time delays observed in two of the channel pairs and that the majority of the large time delays can only be explained by long, expanding loops with photospheric abundances. Additional observations may rule out these simulations as an explanation for the long time delays. We suggest that either the time delays found in this manner may not be representative of real loop evolution, or that the impulsive heating and cooling scenario may be too simple to explain the observations and other heating scenarios must be explored.
△ Less
Submitted 18 December, 2015;
originally announced December 2015.
-
The initial conditions of observed star clusters - I. Method description and validation
Authors:
J. T. Pijloo,
S. F. Portegies Zwart,
P. E. R. Alexander,
M. Gieles,
S. S. Larsen,
P. J. Groot,
B. Devecchi
Abstract:
We have coupled a fast, parametrized star cluster evolution code to a Markov Chain Monte Carlo code to determine the distribution of probable initial conditions of observed star clusters, which may serve as a starting point for future $N$-body calculations. In this paper we validate our method by applying it to a set of star clusters which have been studied in detail numerically with $N$-body simu…
▽ More
We have coupled a fast, parametrized star cluster evolution code to a Markov Chain Monte Carlo code to determine the distribution of probable initial conditions of observed star clusters, which may serve as a starting point for future $N$-body calculations. In this paper we validate our method by applying it to a set of star clusters which have been studied in detail numerically with $N$-body simulations and Monte Carlo methods: the Galactic globular clusters M4, 47 Tucanae, NGC 6397, M22, $ω$ Centauri, Palomar 14 and Palomar 4, the Galactic open cluster M67, and the M31 globular cluster G1. For each cluster we derive a distribution of initial conditions that, after evolution up to the cluster's current age, evolves to the currently observed conditions. We find that there is a connection between the morphology of the distribution of initial conditions and the dynamical age of a cluster and that a degeneracy in the initial half-mass radius towards small radii is present for clusters which have undergone a core collapse during their evolution. We find that the results of our method are in agreement with $N$-body and Monte Carlo studies for the majority of clusters. We conclude that our method is able to find reliable posteriors for the determined initial mass and half-mass radius for observed star clusters, and thus forms an suitable starting point for modeling an observed cluster\rq{}s evolution.
△ Less
Submitted 15 July, 2015;
originally announced July 2015.
-
Trigger Mechanism of Solar Subflares in a Braided Coronal Magnetic Structure
Authors:
Sanjiv K. Tiwari,
Caroline E. Alexander,
Amy R. Winebarger,
Ronald L. Moore
Abstract:
Fine-scale braiding of coronal magnetic loops by continuous footpoint motions may power coronal heating via nanoflares, which are spontaneous fine-scale bursts of internal reconnection. An initial nanoflare may trigger an avalanche of reconnection of the braids, making a microflare or larger subflare. In contrast to this internal triggering of subflares, we observe external triggering of subflares…
▽ More
Fine-scale braiding of coronal magnetic loops by continuous footpoint motions may power coronal heating via nanoflares, which are spontaneous fine-scale bursts of internal reconnection. An initial nanoflare may trigger an avalanche of reconnection of the braids, making a microflare or larger subflare. In contrast to this internal triggering of subflares, we observe external triggering of subflares in a braided coronal magnetic field observed by the {\it High-resolution Coronal Imager (Hi-C)}. We track the development of these subflares using 12 s cadence images acquired by {\it SDO}/AIA in 1600, 193, 94 Å, and registered magnetograms of {\it SDO}/HMI, over four hours centered on the {\it Hi-C} observing time. These data show numerous recurring small-scale brightenings in transition-region emission happening on polarity inversion lines where flux cancellation is occurring. We present in detail an example of an apparent burst of reconnection of two loops in the transition region under the braided coronal field, appropriate for releasing a short reconnected loop downward and a longer reconnected loop upward. The short loop presumably submerges into the photosphere, participating in observed flux cancellation. A subflare in the overlying braided magnetic field is apparently triggered by the disturbance of the braided field by the reconnection-released upward loop. At least 10 subflares observed in this braided structure appear to be triggered this way. How common this external trigger mechanism for coronal subflares is in other active regions, and how important it is for coronal heating in general, remain to be seen.
△ Less
Submitted 15 October, 2014;
originally announced October 2014.
-
The Effects of Regional Vaccination Heterogeneity on Measles Outbreaks with France as a Case Study
Authors:
Elaine T. Alexander,
Savanah D. McMahon,
Nicholas Roberts,
Emilio Sutti,
Daniel Burkow,
Miles Manning,
Kamuela E. Yong,
Sergei Suslov
Abstract:
The rubeola virus, commonly known as measles, is one of the major causes of vaccine-preventable deaths among children worldwide. This is the case despite the fact that an effective vaccine is widely available. Even in developed countries elimination efforts have fallen short as seen by recent outbreaks in Europe, which had over 30,000 cases reported in 2010. The string of measles outbreaks in Fran…
▽ More
The rubeola virus, commonly known as measles, is one of the major causes of vaccine-preventable deaths among children worldwide. This is the case despite the fact that an effective vaccine is widely available. Even in developed countries elimination efforts have fallen short as seen by recent outbreaks in Europe, which had over 30,000 cases reported in 2010. The string of measles outbreaks in France from 2008-2011 is of particular interest due to the documented disparity in regional vaccination coverage. The impact of heterogeneous vaccine coverage on disease transmission is a broad interest and the focus of this study. A Susceptible-Exposed-Infectious-Recovered (SEIR) multi-patch epidemiological model capturing the regional differences in vaccination rates and mixing is introduced. The mathematical analysis of a two-patch system is carried out to help our understanding of the behavior of multi-patch systems. Numerical simulations are generated to aid the study of the system's qualitative dynamics. Data from the recent French outbreaks were used to generate parameter values and to help connect theory with application. Our findings show that heterogeneous vaccination coverage increases controlled reproduction number compared to comparable homogeneous coverage.
△ Less
Submitted 31 July, 2014;
originally announced August 2014.
-
Sparkling EUV bright dots observed with Hi-C
Authors:
S. Regnier,
C. E. Alexander,
R. W. Walsh,
A. R. Winebarger,
J. Cirtain,
L. Golub,
K. E. Korreck,
N. Mitchell,
S. Platt,
M. Weber,
B. De Pontieu,
A. Title,
K. Kobayashi,
S. Kuzin,
C. E. DeForest
Abstract:
Observing the Sun at high time and spatial scales is a step towards understanding the finest and fundamental scales of heating events in the solar corona. The Hi-C instrument has provided the highest spatial and temporal resolution images of the solar corona in the EUV wavelength range to date. Hi-C observed an active region on 11 July 2012, which exhibits several interesting features in the EUV l…
▽ More
Observing the Sun at high time and spatial scales is a step towards understanding the finest and fundamental scales of heating events in the solar corona. The Hi-C instrument has provided the highest spatial and temporal resolution images of the solar corona in the EUV wavelength range to date. Hi-C observed an active region on 11 July 2012, which exhibits several interesting features in the EUV line at 193Å: one of them is the existence of short, small brightenings ``sparkling" at the edge of the active region; we call these EUV Bright Dots (EBDs). Individual EBDs have a characteristic duration of 25s with a characteristic length of 680 km. These brightenings are not fully resolved by the SDO/AIA instrument at the same wavelength, however, they can be identified with respect to the Hi-C location of the EBDs. In addition, EBDs are seen in other chromospheric/coronal channels of SDO/AIA suggesting a temperature between 0.5 and 1.5 MK. Based on their frequency in the Hi-C time series, we define four different categories of EBDs: single peak, double peak, long duration, and bursty EBDs. Based on a potential field extrapolation from an SDO/HMI magnetogram, the EBDs appear at the footpoints of large-scale trans-equatorial coronal loops. The Hi-C observations provide the first evidence of small-scale EUV heating events at the base of these coronal loops, which have a free magnetic energy of the order of 10$^{26}$ erg.
△ Less
Submitted 11 February, 2014;
originally announced February 2014.
-
Anti-parallel EUV flows observed along active region filament threads with Hi-C
Authors:
Caroline E. Alexander,
Robert W. Walsh,
Stephane Regnier,
Jonathan Cirtain,
Amy R. Winebarger,
Leon Golub,
Ken Kobayashi,
Simon Platt,
Nick Mitchell,
Kelly Korreck,
Bart DePontieu,
Craig DeForest,
Mark Weber,
Alan Title,
Sergey Kuzin
Abstract:
Plasma flows within prominences/filaments have been observed for many years and hold valuable clues concerning the mass and energy balance within these structures. Previous observations of these flows primarily come from H-alpha and cool EUV lines (e.g., 304A) where estimates of the size of the prominence threads has been limited by the resolution of the available instrumentation. Evidence of `cou…
▽ More
Plasma flows within prominences/filaments have been observed for many years and hold valuable clues concerning the mass and energy balance within these structures. Previous observations of these flows primarily come from H-alpha and cool EUV lines (e.g., 304A) where estimates of the size of the prominence threads has been limited by the resolution of the available instrumentation. Evidence of `counter-steaming' flows has previously been inferred from these cool plasma observations but now, for the first time, these flows have been directly imaged along fundamental filament threads within the million degree corona (at 193A). In this work we present observations of an active region filament observed with Hi-C that exhibits anti-parallel flows along adjacent filament threads. Complementary data from SDO/AIA and HMI are presented. The ultra-high spatial and temporal resolution of Hi-C allow the anti-parallel flow velocities to be measured (70-80 km/s) and gives an indication of the resolvable thickness of the individual strands (0.8'' +/- 0.1''). The temperature distribution of the plasma flows was estimated to be log T(K) = 5.45 +/- 0.10 using EM loci analysis. We find that SDO/AIA cannot clearly observe these anti-parallel flows nor measure their velocity or thread width due to its larger pixel size. We suggest that
anti-parallel/counter-streaming flows are likely commonplace within all filaments and are currently not observed in EUV due to current instrument spatial resolution.
△ Less
Submitted 21 June, 2013;
originally announced June 2013.
-
Constraining the initial conditions of globular clusters using their radius distribution
Authors:
Poul E. R. Alexander,
Mark Gieles
Abstract:
Studies of extra-galactic globular clusters have shown that the peak size of the globular cluster (GC) radius distribution (RD) depends only weakly on galactic environment, and can be used as a standard ruler. We model RDs of GC populations using a simple prescription for a Hubble time of relaxation driven evolution of cluster mass and radius, and explore the conditions under which the RD can be u…
▽ More
Studies of extra-galactic globular clusters have shown that the peak size of the globular cluster (GC) radius distribution (RD) depends only weakly on galactic environment, and can be used as a standard ruler. We model RDs of GC populations using a simple prescription for a Hubble time of relaxation driven evolution of cluster mass and radius, and explore the conditions under which the RD can be used as a standard ruler. We consider a power-law cluster initial mass function (CIMF) with and without an exponential truncation, and focus in particular on a flat and a steep CIMF (power-law indices of 0 and -2, respectively). For the initial half-mass radii at birth we adopt either Roche-lobe filling conditions ('filling',meaning that the ratio of half-mass to Jacobi radius is approximately rh/rJ ~ 0.15) or strongly Roche-lobe under-filling conditions ('under-filling', implying that initially rh/rJ << 0.15). Assuming a constant orbital velocity about the galaxy centre we find for a steep CIMF that the typical half-light radius scales with galactocentric radius RG as RG^1/3. This weak scaling is consistent with observations, but this scenario has the (well known) problem that too many low-mass clusters survive. A flat CIMF with 'filling' initial conditions results in the correct mass function at old ages, but with too many large (massive) clusters at large RG. An 'underfilling' GC population with a flat CIMF also results in the correct mass function, and can also successfully reproduce the shape of the RD, with a peak size that is (almost) independent of RG. In this case, the peak size depends (almost) only on the peak mass of the GC mass function. The (near) universality of the GC RD is therefore because of the (near) universality of the CIMF. There are some extended GCs in the outer halo of the Milky Way that cannot be explained by this model.
△ Less
Submitted 13 February, 2013;
originally announced February 2013.
-
A prescription and fast code for the long-term evolution of star clusters
Authors:
Poul E. R. Alexander,
Mark Gieles
Abstract:
We introduce the star cluster evolution code Evolve Me A Cluster of StarS (EMACSS), a simple yet physically motivated computational model that describes the evolution of some fundamental properties of star clusters in static tidal fields. We base our prescription upon the flow of energy within the cluster, which is a constant fraction of the total energy per half-mass relaxation time. According to…
▽ More
We introduce the star cluster evolution code Evolve Me A Cluster of StarS (EMACSS), a simple yet physically motivated computational model that describes the evolution of some fundamental properties of star clusters in static tidal fields. We base our prescription upon the flow of energy within the cluster, which is a constant fraction of the total energy per half-mass relaxation time. According to Henon's predictions, this flow is independent of the precise mechanisms for energy production within the core, and therefore does not require a complete description of the many-body interactions therein. For a cluster of equal-mass stars, we thence use dynamical theory and analytic descriptions of escape mechanisms to construct a series of coupled differential equations expressing the time evolution of cluster mass and radius. These equations are numerically solved using a fourth-order Runge-Kutta integration kernel, and the results benchmarked against a data base of direct N-body simulations. We use simulations containing a modest initial number of stars (1024 <= N <= 65 536) and point-mass tidal fields of various strengths. Our prescription is publicly available and reproduces the N-body results to within ~10 per cent accuracy for the entire post-collapse evolution of star clusters.
△ Less
Submitted 21 March, 2012;
originally announced March 2012.
-
A new look at a polar crown cavity as observed by SDO/AIA
Authors:
S. Regnier,
R. W. Walsh,
C. E. Alexander
Abstract:
The Solar Dynamics Observatory (SDO) was launched in February 2010 and is now providing an unprecedented view of the solar activity at high spatial resolution and high cadence covering a broad range of temperature layers of the atmosphere. We aim at defining the structure of a polar crown cavity and describing its evolution during the erupting process. We use high cadence time series of SDO/AIA ob…
▽ More
The Solar Dynamics Observatory (SDO) was launched in February 2010 and is now providing an unprecedented view of the solar activity at high spatial resolution and high cadence covering a broad range of temperature layers of the atmosphere. We aim at defining the structure of a polar crown cavity and describing its evolution during the erupting process. We use high cadence time series of SDO/AIA observations at 304Å(50000 K) and 171Å(0.6 MK) to determine the structure of the polar crown cavity and its associated plasma as well as the evolution of the cavity during the different phases of the eruption. We report on the observations recorded on 13 June 2010 located on the North-West limb. We observe coronal plasma shaped by magnetic field lines with a negative curvature (U-shape) sitting at the bottom of a cavity. The cavity is located just above the polar crown filament material. We thus observe the inner part of the cavity above the filament as depicted in the classical three part Coronal Mass Ejection (CME) model composed of a filament, a cavity and a CME front. The filament (in this case a polar crown filament) is part of the cavity and makes a continuous structuring from the filament to the CME front depicted by concentric ellipses (in a 2D cartoon). We propose to define a polar crown cavity as a density depletion sitting above denser polar crown filament plasma drained down the cavity due to gravity. As part of the polar crown filament, plasma at different temperatures (ranging from 50000K to 0.6 MK) is observed at the same location on the cavity dips and sustained by a competition between the gravity and the curvature of magnetic field lines. The eruption of the polar crown cavity as a solid body can be decomposed into two phases: a slow rise at a speed of 0.6 km/s, and an acceleration phase at a mean speed of 25 km/s.
△ Less
Submitted 18 July, 2011;
originally announced July 2011.
-
Generating and solving the mean field and pair approximation equations in epidemiological models
Authors:
Murray E. Alexander,
Randy Kobes
Abstract:
The pair approximation is a simple, low-order method to incorporate effects of local spatial structure in epidemiological models. However, since for K state variables in a model there are K(K+1)/2 equations in the pair approximation, generating these equations, although straightforward, can become tedious. In this paper we describe two programs written in Perl to simplify this process - one to con…
▽ More
The pair approximation is a simple, low-order method to incorporate effects of local spatial structure in epidemiological models. However, since for K state variables in a model there are K(K+1)/2 equations in the pair approximation, generating these equations, although straightforward, can become tedious. In this paper we describe two programs written in Perl to simplify this process - one to construct the equations, and the other to generate Matlab/Octave functions to numerically integrate the equations using a positivity-preserving method. A third utility program is also included which generates a Maple file that can be used within the Maple symbolic manipulation program to simplify algebraically some of the terms in the generated file describing the equations, as well as to check that the usual combination of equations sum up to zero, which is expected in cases where the total population sizes are conserved.
△ Less
Submitted 16 July, 2010;
originally announced July 2010.
-
Expansion of the conditional probability function in a network with nearest-neighbour degree correlations
Authors:
Murray E. Alexander,
Randy Kobes
Abstract:
A useful property of a network that can be used to characterize many systems is the degree distribution. However, many complex networks exhibit higher--order degree correlations that must be studied through other means, such as clustering coefficients, the Newman r factor, and the average nearest neighbour degree (ANND). In this paper we develop an expansion of the conditional probability that can…
▽ More
A useful property of a network that can be used to characterize many systems is the degree distribution. However, many complex networks exhibit higher--order degree correlations that must be studied through other means, such as clustering coefficients, the Newman r factor, and the average nearest neighbour degree (ANND). In this paper we develop an expansion of the conditional probability that can be used to parameterize such degree correlations. The measures of degree correlations associated with this expansion can be used to signal the presence of non--linear correlations.
△ Less
Submitted 5 July, 2010;
originally announced July 2010.