Showing 1–50 of 136 results for author: Barnes, G

OpenMedLM: Prompt engineering can out-perform fine-tuning in medical question-answering with open-source large language models

Authors: Jenish Maharjan, Anurag Garikipati, Navan Preet Singh, Leo Cyrus, Mayank Sharma, Madalina Ciobanu, Gina Barnes, Rahul Thapa, Qingqing Mao, Ritankar Das

Abstract: LLMs have become increasingly capable at accomplishing a range of specialized-tasks and can be utilized to expand equitable access to medical knowledge. Most medical LLMs have involved extensive fine-tuning, leveraging specialized medical data and significant, thus costly, amounts of computational power. Many of the top performing LLMs are proprietary and their access is limited to very few resear… ▽ More LLMs have become increasingly capable at accomplishing a range of specialized-tasks and can be utilized to expand equitable access to medical knowledge. Most medical LLMs have involved extensive fine-tuning, leveraging specialized medical data and significant, thus costly, amounts of computational power. Many of the top performing LLMs are proprietary and their access is limited to very few research groups. However, open-source (OS) models represent a key area of growth for medical LLMs due to significant improvements in performance and an inherent ability to provide the transparency and compliance required in healthcare. We present OpenMedLM, a prompting platform which delivers state-of-the-art (SOTA) performance for OS LLMs on medical benchmarks. We evaluated a range of OS foundation LLMs (7B-70B) on four medical benchmarks (MedQA, MedMCQA, PubMedQA, MMLU medical-subset). We employed a series of prompting strategies, including zero-shot, few-shot, chain-of-thought (random selection and kNN selection), and ensemble/self-consistency voting. We found that OpenMedLM delivers OS SOTA results on three common medical LLM benchmarks, surpassing the previous best performing OS models that leveraged computationally costly extensive fine-tuning. The model delivers a 72.6% accuracy on the MedQA benchmark, outperforming the previous SOTA by 2.4%, and achieves 81.7% accuracy on the MMLU medical-subset, establishing itself as the first OS LLM to surpass 80% accuracy on this benchmark. Our results highlight medical-specific emergent properties in OS LLMs which have not yet been documented to date elsewhere, and showcase the benefits of further leveraging prompt engineering to improve the performance of accessible LLMs for medical applications. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Permutation invariant tensor models and partition algebras

Authors: George Barnes, Adrian Padellaro, Sanjaye Ramgoolam

Abstract: Matrix models with continuous symmetry are powerful tools for studying quantum gravity and holography. Tensor models have also found applications in holographic quantum gravity. Matrix models with discrete permutation symmetry have been shown to satisfy large $N$ factorisation properties relevant to holography, while also having applications to the statistical analysis of ensembles of real-world m… ▽ More Matrix models with continuous symmetry are powerful tools for studying quantum gravity and holography. Tensor models have also found applications in holographic quantum gravity. Matrix models with discrete permutation symmetry have been shown to satisfy large $N$ factorisation properties relevant to holography, while also having applications to the statistical analysis of ensembles of real-world matrices. Here we develop 3-index tensor models in dimension $D$ with a discrete symmetry of permutations in the symmetric group $S_D$. We construct the most general permutation invariant Gaussian tensor model using the representation theory of symmetric groups and associated partition algebras. We define a representation basis for the 3-index tensors, where the two-point function is diagonalised. Inverting the change of basis gives an explicit formula for the two-point function in the tensor basis for general $D$. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: 28 pages + 11 pages appendices, 5 figures

Report number: QMUL-PH-23-35

Permutation invariant Gaussian matrix models for financial correlation matrices

Authors: George Barnes, Sanjaye Ramgoolam, Michael Stephanou

Abstract: We construct an ensemble of correlation matrices from high-frequency foreign exchange market data, with one matrix for every day for 446 days. The matrices are symmetric and have vanishing diagonal elements after subtracting the identity matrix. For this case, we construct the general permutation invariant Gaussian matrix model, which has 4 parameters characterised using the representation theory… ▽ More We construct an ensemble of correlation matrices from high-frequency foreign exchange market data, with one matrix for every day for 446 days. The matrices are symmetric and have vanishing diagonal elements after subtracting the identity matrix. For this case, we construct the general permutation invariant Gaussian matrix model, which has 4 parameters characterised using the representation theory of symmetric groups. The permutation invariant polynomial functions of the symmetric, diagonally vanishing matrices have a basis labelled by undirected loop-less graphs. Using the expectation values of the general linear and quadratic permutation invariant functions of the matrices in the dataset, the 4 parameters of the matrix model are determined. The model then predicts the expectation values of the cubic and quartic polynomials. These predictions are compared to the data to give strong evidence for a good overall fit of the permutation invariant Gaussian matrix model. The linear, quadratic, cubic and quartic polynomial functions are then used to define low-dimensional feature vectors for the days associated to the matrices. These vectors, with choices informed by the refined structure of small non-Gaussianities, are found to be effective as a tool for anomaly detection in market states: statistically significant correlations are established between atypical days as defined using these feature vectors, and days with significant economic events as recognized in standard foreign exchange economic calendars. They are also shown to be useful as a tool for ranking pairs of days in terms of their similarity, yielding a strongly statistically significant correlation with a ranking based on a higher dimensional proxy for visual similarity. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 49 pages + appendices, 8 figures

Report number: QMUL-PH-23-08

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

Implications of Different Solar Photospheric Flux-Transport Models for Global Coronal and Heliospheric Modeling

Authors: Graham Barnes, Marc L. DeRosa, Shaela I. Jones, Charles N. Arge, Carl J. Henney, Mark C. M. Cheung

Abstract: The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open… ▽ More The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open magnetic field footpoints, and the presence of coronal magnetic null points throughout various phases of a solar activity cycle, focusing on the months of April in even-numbered years between 2012 and 2020, inclusive. We find that there is a solar cycle dependence to each of the metrics considered, but there is not a single phase of the cycle in which all the metrics indicate good agreement between the models. The metrics also reveal large, transient differences between the models when a new active region is rotating into the assimilation window. The evolution of the surface flux is governed by a combination of large scale flows and comparatively small scale motions associated with convection. Because the latter flows evolve rapidly, there are intervals during which their impact on the surface flux can only be characterized in a statistical sense, thus their impact is modeled by introducing a random evolution that reproduces the typical surface flux evolution. We find that the differences between the predicted properties are dominated by differences in the model assumptions and implementation, rather than selection of a particular realization of the random evolution. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted for publication in The Astrophysical Journal

Properties of Flare-Imminent versus Flare-Quiet Active Regions from the Chromosphere through the Corona II: NonParametric Discriminant Analysis Results from the NWRA Classification Infrastructure (NCI)

Authors: KD Leka, Karin Dissauer, Graham Barnes, Eric L. Wagner

Abstract: A large sample of active-region-targeted time-series images from the Solar Dynamics Observatory / Atmospheric Imaging Assembly, the AIA Active Region Patch database ("AARPs", Paper I: Dissauer et al 2022) is used to investigate whether parameters describing the coronal, transition region, and chromospheric emission can differentiate a region that will imminently produce a solar flare from one that… ▽ More A large sample of active-region-targeted time-series images from the Solar Dynamics Observatory / Atmospheric Imaging Assembly, the AIA Active Region Patch database ("AARPs", Paper I: Dissauer et al 2022) is used to investigate whether parameters describing the coronal, transition region, and chromospheric emission can differentiate a region that will imminently produce a solar flare from one that will not. Parametrizations based on moment analysis of direct and running-difference images provide for physically-interpretable results from nonparametric discriminant analysis. Across four event definitions including both 24hr and 6hr validity periods, 160 image-based parameters capture the general state of the atmosphere, rapid brightness changes, and longer-term intensity evolution. We find top Brier Skill Scores in the 0.07--0.33 range, True Skill Statistics in the 0.68--0.82 range (both depending on event definition), and Receiver Operating Characteristic Skill Scores above 0.8. Total emission can perform notably as can steeply increasing or decreasing brightness, although mean brightness measures do not, demonstrating the well-known active-region-size/flare-productivity relation. Once a region is flare productive, the active-region coronal plasma appears to stay hot. The 94AA filter data provides the most parameters with discriminating power, with indications that it benefits from sampling multiple physical regimes. In particular, classification success using higher-order moments of running difference images indicate a propensity for flare-imminent regions to display short-lived small-scale brightening events. Parameters describing the evolution of the corona can provide flare-imminent indicators, but at no preference over "static" parameters. Finally, all parameters and NPDA-derived probabilities are available to the community for additional research. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Astrophysical Journal, in press

Large-Scale Spatial Cross-Calibration of Hinode/SOT-SP and SDO/HMI

Authors: David F. Fouhey, Richard E. L. Higgins, Spiro K. Antiochos, Graham Barnes, Marc L. DeRosa, J. Todd Hoeksema, K. D. Leka, Yang Liu, Peter W. Schuck, Tamas I. Gombosi

Abstract: We investigate the cross-calibration of the Hinode/SOT-SP and SDO/HMI instrument meta-data, specifically the correspondence of the scaling and pointing information. Accurate calibration of these datasets gives the correspondence needed by inter-instrument studies and learning-based magnetogram systems, and is required for physically-meaningful photospheric magnetic field vectors. We approach the p… ▽ More We investigate the cross-calibration of the Hinode/SOT-SP and SDO/HMI instrument meta-data, specifically the correspondence of the scaling and pointing information. Accurate calibration of these datasets gives the correspondence needed by inter-instrument studies and learning-based magnetogram systems, and is required for physically-meaningful photospheric magnetic field vectors. We approach the problem by robustly fitting geometric models on correspondences between images from each instrument's pipeline. This technique is common in computer vision, but several critical details are required when using scanning slit spectrograph data like Hinode/SOT-SP. We apply this technique to data spanning a decade of the Hinode mission. Our results suggest corrections to the published Level 2 Hinode/SOT-SP data. First, an analysis on approximately 2,700 scans suggests that the reported pixel size in Hinode/SOT-SP Level 2 data is incorrect by around 1%. Second, analysis of over 12,000 scans show that the pointing information is often incorrect by dozens of arcseconds with a strong bias. Regression of these corrections indicates that thermal effects have caused secular and cyclic drift in Hinode/SOT-SP pointing data over its mission. We offer two solutions. First, direct co-alignment with SDO/HMI data via our procedure can improve alignments for many Hinode/SOT-SP scans. Second, since the pointing errors are predictable, simple post-hoc corrections can substantially improve the pointing. We conclude by illustrating the impact of this updated calibration on derived physical data products needed for research and interpretation. Among other things, our results suggest that the pointing errors induce a hemispheric bias in estimates of radial current density. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: Under revisions at ApJS

Permutation symmetry in large N Matrix Quantum Mechanics and Partition Algebras

Authors: George Barnes, Adrian Padellaro, Sanjaye Ramgoolam

Abstract: We describe the implications of permutation symmetry for the state space and dynamics of quantum mechanical systems of matrices of general size $N$. We solve the general 11- parameter permutation invariant quantum matrix harmonic oscillator Hamiltonian and calculate the canonical partition function. The permutation invariant sector of the Hilbert space, for general Hamiltonians, can be described u… ▽ More We describe the implications of permutation symmetry for the state space and dynamics of quantum mechanical systems of matrices of general size $N$. We solve the general 11- parameter permutation invariant quantum matrix harmonic oscillator Hamiltonian and calculate the canonical partition function. The permutation invariant sector of the Hilbert space, for general Hamiltonians, can be described using partition algebra diagrams forming the bases of a tower of partition algebras $P_k(N)$. The integer $k$ is interpreted as the degree of matrix oscillator polynomials in the quantum mechanics. Families of interacting Hamiltonians are described which are diagonalised by a representation theoretic basis for the permutation invariant subspace which we construct for $ N \ge 2k $. These include Hamiltonians for which the low-energy states are permutation invariant and can give rise to large ground state degeneracies related to the dimensions of partition algebras. A symmetry-based mechanism for quantum many body scars discussed in the literature can be realised in these matrix systems with permutation symmetry. A mapping of the matrix index values to lattice sites allows a realisation of the mechanism in the context of modified Bose-Hubbard models. Extremal correlators analogous to those studied in AdS/CFT are shown to obey selection rules based on Clebsch-Gordan multiplicities (Kronecker coefficients) of symmetric groups. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 64 pages + Appendices, 2 figures

Report number: QMUL-PH-22-19

The eta Aquilae System: Radial Velocities and Astrometry in Search of eta Aql B

Authors: G. Fritz Benedict, Thomas G. Barnes III, Nancy R. Evans, William D. Cochran, Richard I. Anderson, Barbara E. McArthur, Thomas E. Harrison

Abstract: The classical Cepheid eta Aql was not included in past Leavitt Law work (Benedict et al. 2007) because of a presumed complicating orbit due to a known B9.8V companion. To determine the orbit of eta Aql B, we analyze a significant number of radial velocity measures (RV) from eight sources. With these we establish the RV variation due to Cepheid pulsation, using a twelve Fourier coefficient model, w… ▽ More The classical Cepheid eta Aql was not included in past Leavitt Law work (Benedict et al. 2007) because of a presumed complicating orbit due to a known B9.8V companion. To determine the orbit of eta Aql B, we analyze a significant number of radial velocity measures (RV) from eight sources. With these we establish the RV variation due to Cepheid pulsation, using a twelve Fourier coefficient model, while solving for velocity offsets required to bring the RV data sets into coincidence. RV residuals provide no evidence of orbital motion, suggesting either nearly face-on orientation or very long period. Reanalysis of Hubble Space Telescope Fine Guidance Sensor astrometry now includes reference star parallax and proper motion priors from Gaia EDR3. As modeling confirmation, we reanalyze zeta Gem in parallel, deriving zeta Gem parallax and proper motion values consistent with Gaia EDR3, and consistent with the Benedict 2007 Leavitt Law. In an effort to further characterize eta Aql B, we hypothesize that eta Aql residuals larger than those of the associated reference stars or a parallax inconsistent with EDR3 and the Benedict 2007 Leavitt Law indicate unmodeled orbital motion. Using the astrometric noise or parallax mismatch with EDR3 we estimate possible periods and mass for eta Aql B. Ascribing photocenter motion to the photometric variation of the Cepheid, eta Aql A, yields a plausible separation, consistent with a long period, explaining the lack of RV variation. None of these approaches yields an unassailable characterization of the eta Aql A-B system △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: to appear in The Astronomical Journal

Hidden symmetries and Large N factorisation for permutation invariant matrix observables

Authors: George Barnes, Adrian Padellaro, Sanjaye Ramgoolam

Abstract: Permutation invariant polynomial functions of matrices have previously been studied as the observables in matrix models invariant under $S_N$, the symmetric group of all permutations of $N$ objects. In this paper, the permutation invariant matrix observables (PIMOs) of degree $k$ are shown to be in one-to-one correspondence with equivalence classes of elements in the diagrammatic partition algebra… ▽ More Permutation invariant polynomial functions of matrices have previously been studied as the observables in matrix models invariant under $S_N$, the symmetric group of all permutations of $N$ objects. In this paper, the permutation invariant matrix observables (PIMOs) of degree $k$ are shown to be in one-to-one correspondence with equivalence classes of elements in the diagrammatic partition algebra $P_k(N)$. On a 4-dimensional subspace of the 13-parameter space of $S_N$ invariant Gaussian models, there is an enhanced $O(N)$ symmetry. At a special point in this subspace, is the simplest $O(N)$ invariant action. This is used to define an inner product on the PIMOs which is expressible as a trace of a product of elements in the partition algebra. The diagram algebra $P_k(N)$ is used to prove the large $N$ factorisation property of this inner product, which generalizes a familiar large $N$ factorisation for inner products of matrix traces invariant under continuous symmetries. △ Less

Submitted 9 April, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 35 pages, 1 figure, V2: typos corrected and refs added, V3: sage code added

Report number: QMUL-PH-21-46

Interference suppression techniques for OPM-based MEG: Opportunities and challenges

Authors: Robert A Seymour, Nicholas Alexander, Stephanie Mellor, George C O'Neill, Tim M Tierney, Gareth R Barnes, Eleanor A Maguire

Abstract: One of the primary technical challenges facing magnetoencephalography (MEG) is that the magnitude of neuromagnetic fields is several orders of magnitude lower than interfering signals. Recently, a new type of sensor has been developed - the optically pumped magnetometer (OPM). These sensors can be placed directly on the scalp and move with the head during participant movement, making them wearable… ▽ More One of the primary technical challenges facing magnetoencephalography (MEG) is that the magnitude of neuromagnetic fields is several orders of magnitude lower than interfering signals. Recently, a new type of sensor has been developed - the optically pumped magnetometer (OPM). These sensors can be placed directly on the scalp and move with the head during participant movement, making them wearable. This opens up a range of exciting experimental and clinical opportunities for OPM-based MEG experiments, including paediatric studies, and the incorporation of naturalistic movements into neuroimaging paradigms. However, OPMs face some unique challenges in terms of interference suppression, especially in situations involving mobile participants, and when OPMs are integrated with electrical equipment required for naturalistic paradigms, such as motion capture systems. Here we briefly review various hardware solutions for OPM interference suppression. We then outline several signal processing strategies aimed at increasing the signal from neuromagnetic sources. These include regression-based strategies, temporal filtering and spatial filtering approaches. The focus is on the practical application of these signal processing algorithms to OPM data. In a similar vein, we include two worked-through experiments using OPM data collected from a whole-head sensor array. These tutorial-style examples illustrate how the steps for suppressing external interference can be implemented, including the associated data and code so that researchers can try the pipelines for themselves. With the popularity of OPM-based MEG rising, there will be an increasing need to deal with interference suppression. We hope this practical paper provides a resource for OPM-based MEG researchers to build upon. △ Less

Submitted 29 November, 2021; v1 submitted 6 October, 2021; originally announced October 2021.

Comments: 56 pages, 19 figures, supplementary materials available on request

SynthIA: A Synthetic Inversion Approximation for the Stokes Vector Fusing SDO and Hinode into a Virtual Observatory

Authors: Richard E. L. Higgins, David F. Fouhey, Spiro K. Antiochos, Graham Barnes, Mark C. M. Cheung, J. Todd Hoeksema, KD Leka, Yang Liu, Peter W. Schuck, Tamas I. Gombosi

Abstract: Both NASA's Solar Dynamics Observatory (SDO) and the JAXA/NASA Hinode mission include spectropolarimetric instruments designed to measure the photospheric magnetic field. SDO's Helioseismic and Magnetic Imager (HMI) emphasizes full-disk high-cadence and good spatial resolution data acquisition while Hinode's Solar Optical Telescope Spectro-Polarimeter (SOT-SP) focuses on high spatial resolution an… ▽ More Both NASA's Solar Dynamics Observatory (SDO) and the JAXA/NASA Hinode mission include spectropolarimetric instruments designed to measure the photospheric magnetic field. SDO's Helioseismic and Magnetic Imager (HMI) emphasizes full-disk high-cadence and good spatial resolution data acquisition while Hinode's Solar Optical Telescope Spectro-Polarimeter (SOT-SP) focuses on high spatial resolution and spectral sampling at the cost of a limited field of view and slower temporal cadence. This work introduces a deep-learning system named SynthIA (Synthetic Inversion Approximation), that can enhance both missions by capturing the best of each instrument's characteristics. We use SynthIA to produce a new magnetogram data product, SynodeP (Synthetic Hinode Pipeline), that mimics magnetograms from the higher spectral resolution Hinode/SOT-SP pipeline, but is derived from full-disk, high-cadence, and lower spectral-resolution SDO/HMI Stokes observations. Results on held-out data show that SynodeP has good agreement with the Hinode/SOT-SP pipeline inversions, including magnetic fill fraction, which is not provided by the current SDO/HMI pipeline. SynodeP further shows a reduction in the magnitude of the 24-hour oscillations present in the SDO/HMI data. To demonstrate SynthIA's generality, we show the use of SDO/AIA data and subsets of the HMI data as inputs, which enables trade-offs between fidelity to the Hinode/SOT-SP inversions, number of observations used, and temporal artifacts. We discuss possible generalizations of SynthIA and its implications for space weather modeling. This work is part of the NASA Heliophysics DRIVE Science Center (SOLSTICE) at the University of Michigan under grant NASA 80NSSC20K0600E, and will be open-sourced. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Pulsational instability of pre-main sequence models from accreting protostars. I. Constraining the input physics for accretion with spectroscopic parameters and stellar pulsations

Authors: T. Steindl, K. Zwintz, T. G. Barnes, M. Muellner, E. I. Vorobyov

Abstract: Context: The pre-main sequence evolution is often simplified by choosing classical initial models. These have large initial radii and sufficient uniform contraction to make them fully convective. Contrary to that, real stars are born as small protostellar seeds in collapsing molecular clouds and obtain their final mass by means of accretion. Aims: We aim to constrain the input physics of accretion… ▽ More Context: The pre-main sequence evolution is often simplified by choosing classical initial models. These have large initial radii and sufficient uniform contraction to make them fully convective. Contrary to that, real stars are born as small protostellar seeds in collapsing molecular clouds and obtain their final mass by means of accretion. Aims: We aim to constrain the input physics of accretion on protostellar seeds with observed spectroscopic parameters and stellar pulsations of young stellar objects and pre-main sequence stars. Methods: We conducted a literature search for spectroscopic samples of young stellar objects and pre-main sequence stars including all previously known pulsators. The sample size of pulsating pre-main sequence stars is increased by analysing TESS observations and presenting discoveries in CoRoT data. We employ MESA and GYRE to calculate evolutionary tracks of accreting protostellar seeds in a constant accretion scenario, the subsequent pre-main sequence evolution, and their pulsation properties. The results are then compared with observations to constrain the input physics. Results: We discuss 16 formerly unknown pulsating pre-main sequence stars and candidates that are either of SPB, $δ$ Scuti,$γ$ Doradus or $δ$ Scuti - $γ$ Doradus hybrid type. We find that evolutionary tracks with a mass accretion rate of $5\times10^{-6} M_\odot/{\rm yr}$ and fraction of injected accretion energy of $β=0.1$ provide the best results in enveloping the spectroscopic parameters of pre-main sequence stars in a constant accretion scenario. The calculated instability regions constrain the atmospheric boundary conditions to Eddington Gray atmospheres; we discuss the future potential for additional constraints by instability regions that are dependent on radial order. We present a possible candidate for pulsations in M-type young stellar objects. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 21 pages plus 12 for the appendix, 17 figures plus 27 for the appendix, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 654, A36 (2021)

A Radial Velocity Search for Binary RR Lyrae Variables

Authors: Thomas G. Barnes III, Elisabeth Guggenberger, Katrien Kolenberg

Abstract: We report 272 radial velocities for 19 RR Lyrae variables. For most of the stars we have radial velocities for the complete pulsation cycle. These data are used to determine robust center--of--mass radial velocities that have been compared to values from the literature in a search for evidence of binary systems. Center--of--mass velocities were determined for each star using Fourier Series and Tem… ▽ More We report 272 radial velocities for 19 RR Lyrae variables. For most of the stars we have radial velocities for the complete pulsation cycle. These data are used to determine robust center--of--mass radial velocities that have been compared to values from the literature in a search for evidence of binary systems. Center--of--mass velocities were determined for each star using Fourier Series and Template fits to the radial velocities. Our center--of--mass velocities have uncertainties from $\pm0.16$ km s$^{-1}$ to $\pm$2.5 km s$^{-1}$, with a mean uncertainty of $\pm$0.92 km s$^{-1}$. We combined our center--of--mass velocities with values from the literature to look for deviations from the mean center--of--mass velocity of each star. Fifteen RR Lyrae show no evidence of binary motion (BK And, CI And, Z CVn, DM Cyg, BK Dra, RR Gem, XX Hya, SZ Leo, BX Leo, TT Lyn, CN Lyr, TU Per, U Tri, RV UMa, and AV Vir). In most cases this conclusion is reached due to the sporadic sampling of the center--of--mass velocities over time. Three RR Lyrae show suspicious variation in the center--of--mass velocities that may indicate binary motion but do not prove it (SS Leo, ST Leo, and AO Peg). TU UMa was observed by us near a predicted periastron passage (at 0.14 in orbital phase) but the absence of additional center--of--mass velocities near periastron make the binary detection, based on radial velocities alone, uncertain. Two stars in our sample show $Hγ$ emission in phases 0.9--1.0: SS Leo and TU UMa. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 5 tables, 21 figures

AutoMeKin2021: An open-source program for automated reaction discovery

Authors: Emilio Martínez-Núñez, George L. Barnes, David R. Glowacki, Sabine Kopec, Daniel Peláez, Aurelio Rodríguez, Roberto Rodríguez-Fernández, Robin J. Shannon, James J. P. Stewart, Pablo G. Tahoces, Saulo A. Vazquez

Abstract: AutoMeKin2021 is an updated version of tsscds2018, a program for the automated discovery of reaction mechanisms (J. Comput. Chem. 2018, 39, 1922-1930). This release features a number of new capabilities: rare-event molecular dynamics simulations to enhance reaction discovery, extension of the original search algorithm to study van der Waals complexes, use of chemical knowledge, a new search algori… ▽ More AutoMeKin2021 is an updated version of tsscds2018, a program for the automated discovery of reaction mechanisms (J. Comput. Chem. 2018, 39, 1922-1930). This release features a number of new capabilities: rare-event molecular dynamics simulations to enhance reaction discovery, extension of the original search algorithm to study van der Waals complexes, use of chemical knowledge, a new search algorithm based on bond-order time series analysis, statistics of the chemical reaction networks, a web application to submit jobs, and other features. The source code, manual, installation instructions and the website link are available at: https://rxnkin.usc.es/index.php/AutoMeKin △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 25 pages, 9 figures and 2 tables

Permutation invariant Gaussian 2-matrix models

Authors: George Barnes, Adrian Padellaro, Sanjaye Ramgoolam

Abstract: We construct the general permutation invariant Gaussian 2-matrix model for matrices of arbitrary size $D$. The parameters of the model are given in terms of variables defined using the representation theory of the symmetric group $S_D$. A correspondence is established between the permutation invariant polynomial functions of the matrix variables (the observables of the model) and directed colored… ▽ More We construct the general permutation invariant Gaussian 2-matrix model for matrices of arbitrary size $D$. The parameters of the model are given in terms of variables defined using the representation theory of the symmetric group $S_D$. A correspondence is established between the permutation invariant polynomial functions of the matrix variables (the observables of the model) and directed colored graphs, which sheds light on stability properties in the large $D$ counting of these invariants. The refined counting of the graphs is given in terms of double cosets involving permutation groups defined by the local structure of the graphs. Linear and quadratic observables are transformed to an $S_D$ representation theoretic basis and are used to define the convergent Gaussian measure. The perturbative rules for the computation of expectation values of graph-basis observables of any degree are given in terms of the representation theoretic parameters. Explicit results for a number of observables of degree up to four are given along with a Sage programme that computes general expectation values. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 80 pages, 11 figures

Report number: QMUL-PH-21-19

Fast and Accurate Emulation of the SDO/HMI Stokes Inversion with Uncertainty Quantification

Authors: Richard E. L. Higgins, David F. Fouhey, Dichang Zhang, Spiro K. Antiochos, Graham Barnes, J. Todd Hoeksema, K. D. Leka, Yang Liu, Peter W. Schuck, Tamas I. Gombosi

Abstract: The Helioseismic and Magnetic Imager (HMI) onboard NASA's Solar Dynamics Observatory (SDO) produces estimates of the photospheric magnetic field which are a critical input to many space weather modelling and forecasting systems. The magnetogram products produced by HMI and its analysis pipeline are the result of a per-pixel optimization that estimates solar atmospheric parameters and minimizes dis… ▽ More The Helioseismic and Magnetic Imager (HMI) onboard NASA's Solar Dynamics Observatory (SDO) produces estimates of the photospheric magnetic field which are a critical input to many space weather modelling and forecasting systems. The magnetogram products produced by HMI and its analysis pipeline are the result of a per-pixel optimization that estimates solar atmospheric parameters and minimizes disagreement between a synthesized and observed Stokes vector. In this paper, we introduce a deep learning-based approach that can emulate the existing HMI pipeline results two orders of magnitude faster than the current pipeline algorithms. Our system is a U-Net trained on input Stokes vectors and their accompanying optimization-based VFISV inversions. We demonstrate that our system, once trained, can produce high-fidelity estimates of the magnetic field and kinematic and thermodynamic parameters while also producing meaningful confidence intervals. We additionally show that despite penalizing only per-pixel loss terms, our system is able to faithfully reproduce known systematic oscillations in full-disk statistics produced by the pipeline. This emulation system could serve as an initialization for the full Stokes inversion or as an ultra-fast proxy inversion. This work is part of the NASA Heliophysics DRIVE Science Center (SOLSTICE) at the University of Michigan, under grant NASA 80NSSC20K0600E, and has been open sourced. △ Less

Submitted 27 August, 2021; v1 submitted 31 March, 2021; originally announced March 2021.

On measuring divergence for magnetic field modeling

Authors: S. A. Gilchrist, K. D. Leka, G. Barnes, M. S. Wheatland, M. L. DeRosa

Abstract: A physical magnetic field has a divergence of zero. Numerical error in constructing a model field and computing the divergence, however, introduces a finite divergence into these calculations. A popular metric for measuring divergence is the average fractional flux $\langle |f_{i}| \rangle$. We show that $\langle |f_{i}| \rangle$ scales with the size of the computational mesh, and may be a poor me… ▽ More A physical magnetic field has a divergence of zero. Numerical error in constructing a model field and computing the divergence, however, introduces a finite divergence into these calculations. A popular metric for measuring divergence is the average fractional flux $\langle |f_{i}| \rangle$. We show that $\langle |f_{i}| \rangle$ scales with the size of the computational mesh, and may be a poor measure of divergence because it becomes arbitrarily small for increasing mesh resolution, without the divergence actually decreasing. We define a modified version of this metric that does not scale with mesh size. We apply the new metric to the results of DeRosa et al. (2015), who measured $\langle |f_{i}| \rangle$ for a series of Nonlinear Force-Free Field (NLFFF) models of the coronal magnetic field based on solar boundary data binned at different spatial resolutions. We compute a number of divergence metrics for the DeRosa et al. (2015) data and analyze the effect of spatial resolution on these metrics using a non-parametric method. We find that some of the trends reported by DeRosa et al. (2015) are due to the intrinsic scaling of $\langle |f_{i}| \rangle$. We also find that different metrics give different results for the same data set and therefore there is value in measuring divergence via several metrics. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: Accepted for publication in ApJ

Journal ref: ApJ 900 136 (2020)

A Comparison of Flare Forecasting Methods. IV. Evaluating Consecutive-Day Forecasting Patterns

Authors: Sung-Hong Park, K. D. Leka, Kanya Kusano, Jesse Andries, Graham Barnes, Suzy Bingham, D. Shaun Bloomfield, Aoife E. McCloskey, Veronique Delouille, David Falconer, Peter T. Gallagher, Manolis K. Georgoulis, Yuki Kubo, Kangjin Lee, Sangwoo Lee, Vasily Lobzin, JunChul Mun, Sophie A. Murray, Tarek A. M. Hamad Nageem, Rami Qahwaji, Michael Sharpe, Rob A. Steenburgh, Graham Steward, Michael Terkildsen

Abstract: A crucial challenge to successful flare prediction is forecasting periods that transition between "flare-quiet" and "flare-active". Building on earlier studies in this series (Barnes et al. 2016; Leka et al. 2019a,b) in which we describe methodology, details, and results of flare forecasting comparison efforts, we focus here on patterns of forecast outcomes (success and failure) over multi-day per… ▽ More A crucial challenge to successful flare prediction is forecasting periods that transition between "flare-quiet" and "flare-active". Building on earlier studies in this series (Barnes et al. 2016; Leka et al. 2019a,b) in which we describe methodology, details, and results of flare forecasting comparison efforts, we focus here on patterns of forecast outcomes (success and failure) over multi-day periods. A novel analysis is developed to evaluate forecasting success in the context of catching the first event of flare-active periods, and conversely, of correctly predicting declining flare activity. We demonstrate these evaluation methods graphically and quantitatively as they provide both quick comparative evaluations and options for detailed analysis. For the testing interval 2016-2017, we determine the relative frequency distribution of two-day dichotomous forecast outcomes for three different event histories (i.e., event/event, no-event/event and event/no-event), and use it to highlight performance differences between forecasting methods. A trend is identified across all forecasting methods that a high/low forecast probability on day-1 remains high/low on day-2 even though flaring activity is transitioning. For M-class and larger flares, we find that explicitly including persistence or prior flare history in computing forecasts helps to improve overall forecast performance. It is also found that using magnetic/modern data leads to improvement in catching the first-event/first-no-event transitions. Finally, 15% of major (i.e., M-class or above) flare days over the testing interval were effectively missed due to a lack of observations from instruments away from the Earth-Sun line. △ Less

Submitted 21 January, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 33 pages, 13 figures, accepted for publication in ApJ

A Comparison of Flare Forecasting Methods. III. Systematic Behaviors of Operational Solar Flare Forecasting Systems

Authors: K. D. Leka, Sung-Hong Park, Kanya Kusano, Jesse Andries, Graham Barnes, Suzy Bingham, D. Shaun Bloomfield, Aoife E. McCloskey, Veronique Delouille, David Falconer, Peter T. Gallagher, Manolis K. Georgoulis, Yuki Kubo, Kangjin Lee, Sangwoo Lee, Vasily Lobzin, JunChul Mun, Sophie A. Murray, Tarek A. M. Hamad Nageem, Rami Qahwaji, Michael Sharpe, Rob Steenburgh, Graham Steward, Michael Terkildsen

Abstract: A workshop was recently held at Nagoya University (31 October - 02 November 2017), sponsored by the Center for International Collaborative Research, at the Institute for Space-Earth Environmental Research, Nagoya University, Japan, to quantitatively compare the performance of today's operational solar flare forecasting facilities. Building upon Paper I of this series (Barnes et al. 2016), in Paper… ▽ More A workshop was recently held at Nagoya University (31 October - 02 November 2017), sponsored by the Center for International Collaborative Research, at the Institute for Space-Earth Environmental Research, Nagoya University, Japan, to quantitatively compare the performance of today's operational solar flare forecasting facilities. Building upon Paper I of this series (Barnes et al. 2016), in Paper II (Leka et al. 2019) we described the participating methods for this latest comparison effort, the evaluation methodology, and presented quantitative comparisons. In this paper we focus on the behavior and performance of the methods when evaluated in the context of broad implementation differences. Acknowledging the short testing interval available and the small number of methods available, we do find that forecast performance: 1) appears to improve by including persistence or prior flare activity, region evolution, and a human "forecaster in the loop"; 2) is hurt by restricting data to disk-center observations; 3) may benefit from long-term statistics, but mostly when then combined with modern data sources and statistical approaches. These trends are arguably weak and must be viewed with numerous caveats, as discussed both here and in Paper II. Following this present work, we present in Paper IV a novel analysis method to evaluate temporal patterns of forecasting errors of both types (i.e., misses and false alarms; Park et al. 2019). Hence, most importantly, with this series of papers we demonstrate the techniques for facilitating comparisons in the interest of establishing performance-positive methodologies. △ Less

Submitted 5 July, 2019; originally announced July 2019.

Comments: 23 pages, 6 figures, accepted for publication in The Astrophysical Journal

A Comparison of Flare Forecasting Methods. II. Benchmarks, Metrics and Performance Results for Operational Solar Flare Forecasting Systems

Authors: K. D. Leka, Sung-Hong Park, Kanya Kusano, Jesse Andries, Graham Barnes, Suzy Bingham, D. Shaun Bloomfield, Aoife E. McCloskey, Veronique Delouille, David Falconer, Peter T. Gallagher, Manolis K. Georgoulis, Yuki Kubo, Kangjin Lee, Sangwoo Lee, Vasily Lobzin, JunChul Mun, Sophie A. Murray, Tarek A. M. Hamad Nageem, Rami Qahwaji, Michael Sharpe, Rob Steenburgh, Graham Steward, Michael Terkildsen

Abstract: Solar flares are extremely energetic phenomena in our Solar System. Their impulsive, often drastic radiative increases, in particular at short wavelengths, bring immediate impacts that motivate solar physics and space weather research to understand solar flares to the point of being able to forecast them. As data and algorithms improve dramatically, questions must be asked concerning how well the… ▽ More Solar flares are extremely energetic phenomena in our Solar System. Their impulsive, often drastic radiative increases, in particular at short wavelengths, bring immediate impacts that motivate solar physics and space weather research to understand solar flares to the point of being able to forecast them. As data and algorithms improve dramatically, questions must be asked concerning how well the forecasting performs; crucially, we must ask how to rigorously measure performance in order to critically gauge any improvements. Building upon earlier-developed methodology (Barnes et al, 2016, Paper I), international representatives of regional warning centers and research facilities assembled in 2017 at the Institute for Space-Earth Environmental Research, Nagoya University, Japan to - for the first time - directly compare the performance of operational solar flare forecasting methods. Multiple quantitative evaluation metrics are employed, with focus and discussion on evaluation methodologies given the restrictions of operational forecasting. Numerous methods performed consistently above the "no skill" level, although which method scored top marks is decisively a function of flare event definition and the metric used; there was no single winner. Following in this paper series we ask why the performances differ by examining implementation details (Leka et al. 2019, Paper III), and then we present a novel analysis method to evaluate temporal patterns of forecasting errors in (Park et al. 2019, Paper IV). With these works, this team presents a well-defined and robust methodology for evaluating solar flare forecasting methods in both research and operational frameworks, and today's performance benchmarks against which improvements and new methods may be compared. △ Less

Submitted 5 July, 2019; originally announced July 2019.

Comments: 26 pages, 5 figures, accepted for publication in the Astrophysical Journal Supplement Series

The effect of metallicity on Cepheid Period-Luminosity relations from a Baade-Wesselink analysis of Cepheids in the Milky Way and Magellanic Clouds

Authors: W. Gieren, J. Storm, P. Konorski, M. Górski, B. Pilecki, I. Thompson, G. Pietrzyński, D. Graczyk, T. G. Barnes, P. Fouqué, N. Nardetto, A. Gallenne, P. Karczmarek, K. Suchomska, P. Wielgórski, M. Taormina, B. Zgirski

Abstract: The extragalactic distance scale builds on the Cepheid period-luminosity (PL) relation. In this paper, we want to carry out a strictly differential comparison of the absolute PL relations obeyed by classical Cepheids in the Milky Way (MW), LMC and SMC galaxies. Taking advantage of the substantial metallicity difference among the Cepheid populations in these three galaxies, we want to establish a p… ▽ More The extragalactic distance scale builds on the Cepheid period-luminosity (PL) relation. In this paper, we want to carry out a strictly differential comparison of the absolute PL relations obeyed by classical Cepheids in the Milky Way (MW), LMC and SMC galaxies. Taking advantage of the substantial metallicity difference among the Cepheid populations in these three galaxies, we want to establish a possible systematic trend of the PL relation absolute zero point as a function of metallicity, and determine the size of such an effect in optical and near-infrared photometric bands. We are using the IRSB Baade-Wesselink type method as calibrated by Storm et al. to determine individual distances to the Cepheids in our samples in MW, LMC and SMC. For our analysis, we use a greatly enhanced sample of Cepheids in the SMC (31 stars) as compared to the small sample (5 stars) available in our previous work. We use the distances to determine absolute Cepheid PL relations in optical and near-infrared bands in each of the three galaxies.} {Our distance analysis of 31 SMC Cepheids with periods from 4-69 days yields tight PL relations in all studied bands, with slopes consistent with the corresponding LMC and MW relations. Adopting the very accurately determined LMC slopes for the optical and near-infrared bands, we determine the zero point offsets between the corresponding absolute PL relations in the 3 galaxies. We find that in all bands the metal-poor SMC Cepheids are intrinsically fainter than their more metal-rich counterparts in the LMC and MW. In the $K$ band the metallicity effect is $-0.23\pm0.06$~mag/dex while in the $V,(V-I)$ Wesenheit index it is slightly stronger, $-0.34\pm0.06$~mag/dex. We find some evidence that the PL relation zero point-metallicity relation might be nonlinear, becoming steeper for lower metallicities. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 20 pages, accepted for publication by A&A

Journal ref: A&A 620, A99 (2018)

Immersive Virtual Reality Experiences for All-Sky Data

Authors: C. J. Fluke, D. G. Barnes

Abstract: Spherical coordinate systems, which are ubiquitous in astronomy, cannot be shown without distortion on flat, two-dimensional surfaces. This poses challenges for the two complementary phases of visual exploration -- making discoveries in data by looking for relationships, patterns or anomalies -- and publication -- where the results of an exploration are made available for scientific scrutiny or co… ▽ More Spherical coordinate systems, which are ubiquitous in astronomy, cannot be shown without distortion on flat, two-dimensional surfaces. This poses challenges for the two complementary phases of visual exploration -- making discoveries in data by looking for relationships, patterns or anomalies -- and publication -- where the results of an exploration are made available for scientific scrutiny or communication. This is a long-standing problem, and many practical solutions have been developed. Our allskyVR approach provides a workflow for experimentation with commodity virtual reality head-mounted displays. Using the free, open source S2PLOT programming library, and the A-Frame WebVR browser-based framework, we provide a straightforward way to visualise all-sky catalogues on a user-centred, virtual celestial sphere. The allskyVR distribution contains both a quickstart option, complete with a gaze-based menu system, and a fully customisable mode for those who need more control of the immersive experience. The software is available for download from: https://github.com/cfluke/allskyVR △ Less

Submitted 8 May, 2018; originally announced May 2018.

Comments: 11 pages, 3 figures, accepted for publication in Publications of the Astronomical Society of Australia

The NWRA Classification Infrastructure: Description and Extension to the Discriminant Analysis Flare Forecasting System (DAFFS)

Authors: K. D. Leka, Graham Barnes, Eric L. Wagner

Abstract: A classification infrastructure built upon Discriminant Analysis has been developed at NorthWest Research Associates for examining the statistical differences between samples of two known populations. Originating to examine the physical differences between flare-quiet and flare-imminent solar active regions, we describe herein some details of the infrastructure including: parametrization of large… ▽ More A classification infrastructure built upon Discriminant Analysis has been developed at NorthWest Research Associates for examining the statistical differences between samples of two known populations. Originating to examine the physical differences between flare-quiet and flare-imminent solar active regions, we describe herein some details of the infrastructure including: parametrization of large datasets, schemes for handling "null" and "bad" data in multi-parameter analysis, application of non-parametric multi-dimensional Discriminant Analysis, an extension through Bayes' theorem to probabilistic classification, and methods invoked for evaluating classifier success. The classifier infrastructure is applicable to a wide range of scientific questions in solar physics. We demonstrate its application to the question of distinguishing flare-imminent from flare-quiet solar active regions, updating results from the original publications that were based on different data and much smaller sample sizes. Finally, as a demonstration of "Research to Operations" efforts in the space-weather forecasting context, we present the Discriminant Analysis Flare Forecasting System (DAFFS), a near-real-time operationally-running solar flare forecasting tool that was developed from the research-directed infrastructure. △ Less

Submitted 19 February, 2018; originally announced February 2018.

Comments: J. Space Weather Space Climate: Accepted / in press; access supplementary materials through journal; some figures are less than full resolution for arXiv

Does Nearby Open Flux Affect the Eruptivity of Solar Active Regions?

Authors: Marc L. DeRosa, Graham Barnes

Abstract: The most energetic solar flares are typically associated with the ejection of a cloud of coronal material into the heliosphere in the form of a coronal mass ejection (CME). However, there exist large flares which are not accompanied by a CME. The existence of these non-eruptive flares raises the question of whether such flares suffer from a lack of access to nearby open fields in the vicinity abov… ▽ More The most energetic solar flares are typically associated with the ejection of a cloud of coronal material into the heliosphere in the form of a coronal mass ejection (CME). However, there exist large flares which are not accompanied by a CME. The existence of these non-eruptive flares raises the question of whether such flares suffer from a lack of access to nearby open fields in the vicinity above the flare (reconnection) site. In this study, we use a sample of 56 flares from Sunspot Cycles 23 and 24 to test whether active regions that produce eruptive X-class flares are preferentially located near coronal magnetic field domains that are open to the heliosphere, as inferred from a potential field source surface model. The study shows that X-class flares having access to open fields are eruptive at a higher rate than those for which access is lacking. The significance of this result should be moderated due to the small number of non-eruptive X-class flares in the sample, based on the associated Bayes factor. △ Less

Submitted 30 May, 2018; v1 submitted 4 February, 2018; originally announced February 2018.

Comments: ApJ accepted

A cautionary tale of interpreting O-C diagrams: period instability in a classical RR Lyr Star Z CVn mimicking as a distant companion

Authors: M. Skarka, J. Liška, R. Dřevěný, E. Guggenberger, Á. Sódor, T. G. Barnes, K. Kolenberg

Abstract: We present a comprehensive study of Z CVn, an RR Lyrae star that shows long-term cyclic variations of its pulsation period. A possible explanation suggested from the shape of the O-C diagram is the light travel-time effect, which we thoroughly examine. We used original photometric and spectroscopic measurements and investigated the period evolution using available maximum times spanning more than… ▽ More We present a comprehensive study of Z CVn, an RR Lyrae star that shows long-term cyclic variations of its pulsation period. A possible explanation suggested from the shape of the O-C diagram is the light travel-time effect, which we thoroughly examine. We used original photometric and spectroscopic measurements and investigated the period evolution using available maximum times spanning more than one century. If the binary hypothesis is valid, Z CVn orbits around a black hole with minimal mass of $56.5$ $\mathfrak{M}_{\odot}$ on a very wide ($P_{\rm orbit}=78.3$ years) and eccentric orbit ($e=0.63$). We discuss the probability of a formation of a black hole-RR Lyrae pair and, although we found it possible, there is no observational evidence of the black hole in the direction to Z CVn. However, the main objection against the binary hypothesis is the comparison of the systemic radial velocity curve model and spectroscopic observations that clearly show that Z CVn cannot be bound in such a binary. Therefore, the variations of pulsation period are likely intrinsic to the star. This finding represents a discovery/confirmation of a new type of cyclic period changes in RR Lyrae stars. By the analysis of our photometric data, we found that the Blazhko modulation with period of 22.931 d is strongly dominant in amplitude. The strength of the phase modulation varies and is currently almost undetectable. We also estimated photometric physical parameters of Z CVn and investigated their variations during the Blazhko cycle using the Inverse Baade-Wesselink method. △ Less

Submitted 18 October, 2017; originally announced October 2017.

Comments: 15 pages, 8 tables, 9 figures, accepted for publication in MNRAS

Real-time colouring and filtering with graphics shaders

Authors: Dany Vohl, Christopher J. Fluke, David G. Barnes, Amr H. Hassan

Abstract: Despite the popularity of the Graphics Processing Unit (GPU) for general purpose computing, one should not forget about the practicality of the GPU for fast scientific visualisation. As astronomers have increasing access to three dimensional (3D) data from instruments and facilities like integral field units and radio interferometers, visualisation techniques such as volume rendering offer means t… ▽ More Despite the popularity of the Graphics Processing Unit (GPU) for general purpose computing, one should not forget about the practicality of the GPU for fast scientific visualisation. As astronomers have increasing access to three dimensional (3D) data from instruments and facilities like integral field units and radio interferometers, visualisation techniques such as volume rendering offer means to quickly explore spectral cubes as a whole. As most 3D visualisation techniques have been developed in fields of research like medical imaging and fluid dynamics, many transfer functions are not optimal for astronomical data. We demonstrate how transfer functions and graphics shaders can be exploited to provide new astronomy-specific explorative colouring methods. We present 12 shaders, including four novel transfer functions specifically designed to produce intuitive and informative 3D visualisations of spectral cube data. We compare their utility to classic colour mapping. The remaining shaders highlight how common computation like filtering, smoothing and line ratio algorithms can be integrated as part of the graphics pipeline. We discuss how this can be achieved by utilising the parallelism of modern GPUs along with a shading language, letting astronomers apply these new techniques at interactive frame rates. All shaders investigated in this work are included in the open source software shwirl (Vohl, 2017). △ Less

Submitted 3 July, 2017; originally announced July 2017.

Comments: 24 pages, 19 figures, 14 algorithms, 1 table. Accepted for publication by Monthly Notices of the Royal Astronomical Society

A comprehensive study of young B stars in NGC 2264: I. Space photometry and asteroseismology

Authors: K. Zwintz, E. Moravveji, P. I. Papics, A. Tkachenko, N. Przybilla, M. -F. Nieva, R. Kuschnig, V. Antoci, D. Lorenz, N. Themessl, L. Fossati, T. G. Barnes

Abstract: Space photometric time series of the most massive members of the young open cluster NGC 2264 allow us to study their different sources of variability down to the millimagnitude level and permits a search for Slowly Pulsating B (SPB) type pulsation among objects that are only a few million years old. Our goal is to conduct a homogeneous study of young B type stars in the cluster NGC 2264 using phot… ▽ More Space photometric time series of the most massive members of the young open cluster NGC 2264 allow us to study their different sources of variability down to the millimagnitude level and permits a search for Slowly Pulsating B (SPB) type pulsation among objects that are only a few million years old. Our goal is to conduct a homogeneous study of young B type stars in the cluster NGC 2264 using photometric time series from space in combination with high-resolution spectroscopy and spectropolarimetry obtained from the ground. The latter will be presented in a separate follow-up article. We performed frequency analyses for eleven B stars in the field of the young cluster NGC 2264 using photometric time series from the MOST, CoRoT and Spitzer space telescopes and the routines Period04 and SigSpec. We employ the MESA stellar evolution code in combination with the oscillation code GYRE to identify the pulsation modes for two SPB stars which exhibit short period spacing series. From our analysis we identify four objects that show SPB pulsations, five stars that show rotational modulation of their light curves caused by spots, one star that is identified to be a binary, and one object in the field of the cluster that is found to be a non-member Be star. In two SPB stars we detect a number of regularly spaced pulsation modes that are compatible with being members of a g mode period series. Despite NGC 2264's young age, our analysis illustrates that its B type members have already arrived on the zero-age main sequence (ZAMS). Our asteroseismic analysis yields masses between 4 and 6 Msun and ages between 1 and 6 million years, which agree well to the overall cluster age. △ Less

Submitted 19 March, 2017; originally announced March 2017.

Comments: 19 pages, 20 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 601, A101 (2017)

Evaluating (and Improving) Estimates of the Solar Radial Magnetic Field Component from Line-of-Sight Magnetograms

Authors: K. D. Leka, G. Barnes, E. L. Wagner

Abstract: Although for many solar physics problems the desirable or meaningful boundary is the radial component of the magnetic field $B_{\rm r}$, the most readily available measurement is the component of the magnetic field along the line-of-sight to the observer, $B_{\rm los}$. As this component is only equal to the radial component where the viewing angle is exactly zero, some approximation is required t… ▽ More Although for many solar physics problems the desirable or meaningful boundary is the radial component of the magnetic field $B_{\rm r}$, the most readily available measurement is the component of the magnetic field along the line-of-sight to the observer, $B_{\rm los}$. As this component is only equal to the radial component where the viewing angle is exactly zero, some approximation is required to estimate $B_{\rm r}$ at all other observed locations. In this study, a common approximation known as the "$μ$-correction", which assumes all photospheric field to be radial, is compared to a method which invokes computing a potential field that matches the observed $B_{\rm los}$, from which the potential field radial component, $B_{\rm r}^{\rm pot}$ is recovered. We demonstrate that in regions that are truly dominated by radially-oriented field at the resolution of the data employed, the $μ$-correction performs acceptably if not better than the potential-field approach. However, it is also shown that for any solar structure which includes horizontal fields, i.e. active regions, the potential-field method better recovers both the strength of the radial field and the location of magnetic neutral line. △ Less

Submitted 17 January, 2017; originally announced January 2017.

Comments: Accepted for publication in Solar Physics

Collaborative visual analytics of radio surveys in the Big Data era

Authors: Dany Vohl, Christopher J. Fluke, Amr H. Hassan, David G. Barnes, Virginia A. Kilborn

Abstract: Radio survey datasets comprise an increasing number of individual observations stored as sets of multidimensional data. In large survey projects, astronomers commonly face limitations regarding: 1) interactive visual analytics of sufficiently large subsets of data; 2) synchronous and asynchronous collaboration; and 3) documentation of the discovery workflow. To support collaborative data inquiry,… ▽ More Radio survey datasets comprise an increasing number of individual observations stored as sets of multidimensional data. In large survey projects, astronomers commonly face limitations regarding: 1) interactive visual analytics of sufficiently large subsets of data; 2) synchronous and asynchronous collaboration; and 3) documentation of the discovery workflow. To support collaborative data inquiry, we present encube, a large-scale comparative visual analytics framework. Encube can utilise advanced visualization environments such as the CAVE2 (a hybrid 2D and 3D virtual reality environment powered with a 100 Tflop/s GPU-based supercomputer and 84 million pixels) for collaborative analysis of large subsets of data from radio surveys. It can also run on standard desktops, providing a capable visual analytics experience across the display ecology. Encube is composed of four primary units enabling compute-intensive processing, advanced visualisation, dynamic interaction, parallel data query, along with data management. Its modularity will make it simple to incorporate astronomical analysis packages and Virtual Observatory capabilities developed within our community. We discuss how encube builds a bridge between high-end display systems (such as CAVE2) and the classical desktop, preserving all traces of the work completed on either platform -- allowing the research process to continue wherever you are. △ Less

Submitted 2 December, 2016; originally announced December 2016.

Comments: Accepted for Proceedings of IAU Symposium 325

An interactive, comparative and quantitative 3D visualization system for large-scale spectral-cube surveys using CAVE2

Authors: Dany Vohl, Christopher J. Fluke, Amr H. Hassan, David G. Barnes

Abstract: As the quantity and resolution of spectral-cubes from optical/infrared and radio surveys increase, desktop-based visualization and analysis solutions must adapt and evolve. Novel immersive 3D environments such as the CAVE2 at Monash University can overcome personal computer's visualization limitations. CAVE2 is part advanced 2D/3D visualization space (80 stereo-capable screens providing a total of… ▽ More As the quantity and resolution of spectral-cubes from optical/infrared and radio surveys increase, desktop-based visualization and analysis solutions must adapt and evolve. Novel immersive 3D environments such as the CAVE2 at Monash University can overcome personal computer's visualization limitations. CAVE2 is part advanced 2D/3D visualization space (80 stereo-capable screens providing a total of 84 million pixels) and part supercomputer ($\sim100$ TFLOPS of integrated GPU-based processing power). We present a novel visualization system enabling simultaneous 3D comparative visualization of $\sim100$ spectral-cubes. With CAVE2 augmented by our newly developed web-based controller interface, astronomers can easily organise spectral-cubes on the different display panels, apply real-time transforms to one or many spectral cubes, and request quantitative information about the displayed data. We also discuss how such a solution can help accelerate the discovery rate in varied research scenarios. △ Less

Submitted 3 October, 2016; originally announced October 2016.

Comments: 4 pages, 2 figures, to appear in Proceedings of ADASS XXV (Sydney, Australia; 2015), edited by N. P. F. Lorente & K. Shortridge (San Francisco: ASP), ASP Conf. Series

Large-scale comparative visualisation of sets of multidimensional data

Authors: Dany Vohl, David G. Barnes, Christopher J. Fluke, Govinda Poudel, Nellie Georgiou-Karistianis, Amr H. Hassan, Yuri Benovitski, Tsz Ho Wong, Owen Kaluza, Toan D. Nguyen, C. Paul Bonnington

Abstract: We present encube $-$ a qualitative, quantitative and comparative visualisation and analysis system, with application to high-resolution, immersive three-dimensional environments and desktop displays. encube extends previous comparative visualisation systems by considering: 1) the integration of comparative visualisation and analysis into a unified system; 2) the documentation of the discovery pro… ▽ More We present encube $-$ a qualitative, quantitative and comparative visualisation and analysis system, with application to high-resolution, immersive three-dimensional environments and desktop displays. encube extends previous comparative visualisation systems by considering: 1) the integration of comparative visualisation and analysis into a unified system; 2) the documentation of the discovery process; and 3) an approach that enables scientists to continue the research process once back at their desktop. Our solution enables tablets, smartphones or laptops to be used as interaction units for manipulating, organising, and querying data. We highlight the modularity of encube, allowing additional functionalities to be included as required. Additionally, our approach supports a high level of collaboration within the physical environment. We show how our implementation of encube operates in a large-scale, hybrid visualisation and supercomputing environment using the CAVE2 at Monash University, and on a local desktop, making it a versatile solution. We discuss how our approach can help accelerate the discovery rate in a variety of research scenarios. △ Less

Submitted 3 October, 2016; originally announced October 2016.

Comments: 26 pages, 11 figures, 5 tables. Accepted for publication in PeerJ Computer Science

A Fixed-point Scheme for the Numerical Construction of Magnetohydrostatic Atmospheres in Three Dimensions

Authors: S. A. Gilchrist, D. C. Braun, G. Barnes

Abstract: Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three d… ▽ More Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three dimensions. Our method is a fixed-point iteration scheme that extends the method of Grad and Rubin (Proc. 2nd Int. Conf. on Peaceful Uses of Atomic Energy 31, 190, 1958) to include a finite gravity force. We apply the method to a test case to demonstrate the method in general and our implementation in code in particular. △ Less

Submitted 27 October, 2016; v1 submitted 2 September, 2016; originally announced September 2016.

Comments: Accepted for publication in Solar Physics

A Comparison of Flare Forecasting Methods, I: Results from the "All-Clear" Workshop

Authors: G. Barnes, K. D. Leka, C. J. Schrijver, T. Colak, R. Qahwaji, O. W. Ashamari, Y. Yuan, J. Zhang, R. T. J. McAteer, D. S. Bloomfield, P. A. Higgins, P. T. Gallagher, D. A. Falconer, M. K. Georgoulis, M. S. Wheatland, C. Balch, T. Dunn, E. L. Wagner

Abstract: Solar flares produce radiation which can have an almost immediate effect on the near-Earth environment, making it crucial to forecast flares in order to mitigate their negative effects. The number of published approaches to flare forecasting using photospheric magnetic field observations has proliferated, with varying claims about how well each works. Because of the different analysis techniques a… ▽ More Solar flares produce radiation which can have an almost immediate effect on the near-Earth environment, making it crucial to forecast flares in order to mitigate their negative effects. The number of published approaches to flare forecasting using photospheric magnetic field observations has proliferated, with varying claims about how well each works. Because of the different analysis techniques and data sets used, it is essentially impossible to compare the results from the literature. This problem is exacerbated by the low event rates of large solar flares. The challenges of forecasting rare events have long been recognized in the meteorology community, but have yet to be fully acknowledged by the space weather community. During the interagency workshop on "all clear" forecasts held in Boulder, CO in 2009, the performance of a number of existing algorithms was compared on common data sets, specifically line-of-sight magnetic field and continuum intensity images from MDI, with consistent definitions of what constitutes an event. We demonstrate the importance of making such systematic comparisons, and of using standard verification statistics to determine what constitutes a good prediction scheme. When a comparison was made in this fashion, no one method clearly outperformed all others, which may in part be due to the strong correlations among the parameters used by different methods to characterize an active region. For M-class flares and above, the set of methods tends towards a weakly positive skill score (as measured with several distinct metrics), with no participating method proving substantially better than climatological forecasts. △ Less

Submitted 22 August, 2016; originally announced August 2016.

Comments: Accepted for publication in ApJ

Mapping spaces and automorphism groups of toric noncommutative spaces

Authors: Gwendolyn E. Barnes, Alexander Schenkel, Richard J. Szabo

Abstract: We develop a sheaf theory approach to toric noncommutative geometry which allows us to formalize the concept of mapping spaces between two toric noncommutative spaces. As an application we study the `internalized' automorphism group of a toric noncommutative space and show that its Lie algebra has an elementary description in terms of braided derivations. We develop a sheaf theory approach to toric noncommutative geometry which allows us to formalize the concept of mapping spaces between two toric noncommutative spaces. As an application we study the `internalized' automorphism group of a toric noncommutative space and show that its Lie algebra has an elementary description in terms of braided derivations. △ Less

Submitted 10 April, 2017; v1 submitted 15 June, 2016; originally announced June 2016.

Comments: 30 pages. v2: Final version published in Letters in Mathematical Physics

Report number: EMPG-16-14 MSC Class: 16T05; 18F20; 53D55; 81R60

Journal ref: Letters in Mathematical Physics, 107(9), 1591-1628 (2017)

Working with Nonassociative Geometry and Field Theory

Authors: Gwendolyn E. Barnes, Alexander Schenkel, Richard J. Szabo

Abstract: We review aspects of our formalism for differential geometry on noncommutative and nonassociative spaces which arise from cochain twist deformation quantization of manifolds. We work in the simplest setting of trivial vector bundles and flush out the details of our approach providing explicit expressions for all bimodule operations, and for connections and curvature. As applications, we describe t… ▽ More We review aspects of our formalism for differential geometry on noncommutative and nonassociative spaces which arise from cochain twist deformation quantization of manifolds. We work in the simplest setting of trivial vector bundles and flush out the details of our approach providing explicit expressions for all bimodule operations, and for connections and curvature. As applications, we describe the constructions of physically viable action functionals for Yang-Mills theory and Einstein-Cartan gravity on noncommutative and nonassociative spaces, as first steps towards more elaborate models relevant to non-geometric flux deformations of geometry in closed string theory. △ Less

Submitted 13 February, 2016; v1 submitted 27 January, 2016; originally announced January 2016.

Comments: 20 pages; v2: Reference added; Contribution to the proceedings of the Corfu Summer Institute on Elementary Particle Physics and Gravity, September 1-26, 2015, Corfu, Greece; Final version published in Proceedings of Science

Report number: EMPG-16-02

Journal ref: PoS CORFU2015 (2015) 081

The Ultimate Display

Authors: C. J. Fluke, D. G. Barnes

Abstract: Astronomical images and datasets are increasingly high-resolution and multi-dimensional. The vast majority of astronomers perform all of their visualisation and analysis tasks on low-resolution, two-dimensional desktop monitors. If there were no technological barriers to designing the ultimate stereoscopic display for astronomy, what would it look like? What capabilities would we require of our co… ▽ More Astronomical images and datasets are increasingly high-resolution and multi-dimensional. The vast majority of astronomers perform all of their visualisation and analysis tasks on low-resolution, two-dimensional desktop monitors. If there were no technological barriers to designing the ultimate stereoscopic display for astronomy, what would it look like? What capabilities would we require of our compute hardware to drive it? And are existing technologies even close to providing a true 3D experience that is compatible with the depth resolution of human stereoscopic vision? We consider the CAVE2 (an 80 Megapixel, hybrid 2D and 3D virtual reality environment directly integrated with a 100 Tflop/s GPU-powered supercomputer) and the Oculus Rift (a low- cost, head-mounted display) as examples at opposite financial ends of the immersive display spectrum. △ Less

Submitted 13 January, 2016; originally announced January 2016.

Comments: 4 pages, 2 figures, to appear in Proceedings of ADASS XXV (Sydney, Australia; 2015), edited by N. P. F. Lorente & K. Shortridge (San Francisco: ASP), ASP Conf. Series

New systemic radial velocities of suspected RR Lyrae binary stars

Authors: Elisabeth Guggenberger, Thomas G. Barnes, Katrien Kolenberg

Abstract: Among the tens of thousands of known RR Lyrae stars there are only a handful that show indications of possible binarity. The question why this is the case is still unsolved, and has recently sparked several studies dedicated to the search for additional RR Lyraes in binary systems. Such systems are particularly valuable because they might allow to constrain the stellar mass. Most of the recent stu… ▽ More Among the tens of thousands of known RR Lyrae stars there are only a handful that show indications of possible binarity. The question why this is the case is still unsolved, and has recently sparked several studies dedicated to the search for additional RR Lyraes in binary systems. Such systems are particularly valuable because they might allow to constrain the stellar mass. Most of the recent studies, however, are based on photometry by finding a light time effect in the timings of maximum light. This approach is a very promising and successful one, but it has a major drawback: by itself, it cannot serve as a definite proof of binarity, because other phenomena such as the Blazhko effect or intrinsic period changes could lead to similar results. Spectroscopic radial velocity measurements, on the other hand, can serve as definite proof of binarity. We have therefore started a project to study spectroscopically RR Lyrae stars that are suspected to be binaries. We have obtained radial velocity (RV) curves with the 2.1m telescope at McDonald observatory. From these we derive systemic RVs which we will compare to previous measurements in order to find changes induced by orbital motions. We also construct templates of the RV curves that can facilitate future studies. We also observed the most promising RR Lyrae binary candidate, TU UMa, as no recent spectroscopic measurements were available. We present a densely covered pulsational RV curve, which will be used to test the predictions of the orbit models that are based on the O-C variations. △ Less

Submitted 2 December, 2015; originally announced December 2015.

Comments: 4 pages, 1 figure, to appear in the No. 105 issue of the Communications from the Konkoly Observatory of the Hungarian Academy of Sciences

Achieving Consistent Doppler Measurements from SDO/HMI Vector Field Inversions

Authors: P. W. Schuck, Spiro Antiochos, K. D. Leka, Graham Barnes

Abstract: NASA's Solar Dynamics Observatory is delivering vector field observations of the full solar disk with unprecedented temporal and spatial resolution; however, the satellite is in a highly inclined geostationary orbit. The relative spacecraft-Sun velocity varies by $\pm3$~km/s over a day which introduces major orbital artifacts in the Helioseismic Magnetic Imager data. We demonstrate that the orbita… ▽ More NASA's Solar Dynamics Observatory is delivering vector field observations of the full solar disk with unprecedented temporal and spatial resolution; however, the satellite is in a highly inclined geostationary orbit. The relative spacecraft-Sun velocity varies by $\pm3$~km/s over a day which introduces major orbital artifacts in the Helioseismic Magnetic Imager data. We demonstrate that the orbital artifacts contaminate all spatial and temporal scales in the data. We describe a newly-developed three stage procedure for mitigating these artifacts in the Doppler data derived from the Milne-Eddington inversions in the HMI Pipeline. This procedure was applied to full disk images of AR11084 to produce consistent Dopplergrams. The data adjustments reduce the power in the orbital artifacts by 31dB. Furthermore, we analyze in detail the corrected images and show that our procedure greatly improve the temporal and spectral properties of the data without adding any new artifacts. We conclude that this new and easily implemented procedure makes a dramatic improvement in the consistency of the HMI data and in its usefulness for precision scientific studies. △ Less

Submitted 20 November, 2015; originally announced November 2015.

Comments: 58 pages, 19 figures, submitted to ApJ

Quantum Thermodynamics, Entropy of the Universe, Free Energy, and the Second Law

Authors: George L. Barnes, Phillip C. Lotshaw, Michael E. Kellman

Abstract: We take the view that the standard von Neumann definition, in which the entropy $S^{vN}$ of a pure state is zero, is in evident conflict with the statement of the second law that the entropy of the universe $S_{univ}$ increases in spontaneous processes, $ ΔS_{univ} > 0$. Here we seek an alternative entropy of the universe $S_{univ}$ that is in accord with the second law, in a spirit not dissimilar… ▽ More We take the view that the standard von Neumann definition, in which the entropy $S^{vN}$ of a pure state is zero, is in evident conflict with the statement of the second law that the entropy of the universe $S_{univ}$ increases in spontaneous processes, $ ΔS_{univ} > 0$. Here we seek an alternative entropy of the universe $S_{univ}$ that is in accord with the second law, in a spirit not dissimilar to von Neumann himself in lesser-known work. We perform simulations of time dependent dynamics for a previously developed \cite{polyadbath} model quantum system becoming entangled with a quantum environment. We test the new definition of the entropy of the system-environment "universe" against the standard thermodynamic relation $ ΔF_{sys} = - T ΔS_{univ} $, calculating system properties using the reduced density matrix and standard von Neumann entropy. Good agreement is obtained, showing the compatibility of an entropy for a pure state of a universe with the statement of the second law and the concept of free energy. Interesting deviation from microcanonical behavior within the zero order energy shell is observed in a context of effectively microcanonical behavior within the much larger total basis of the time dependent universe. △ Less

Submitted 23 October, 2018; v1 submitted 19 November, 2015; originally announced November 2015.

Comments: 24 pages, 4 figures, and 1 table. In v3 we have added Phillip Lotshaw as a coauthor. We have extensively revised the discussion of the meaning of the entropy introduced here, and related this work considerably more to other contemporary work on foundations of quantum thermodynamics and statistical mechanics

The Influence of Spatial Resolution on Nonlinear Force-Free Modeling

Authors: M. L. DeRosa, M. S. Wheatland, K. D. Leka, G. Barnes, T. Amari, A. Canou, S. A. Gilchrist, J. K. Thalmann, G. Valori, T. Wiegelmann, C. J. Schrijver, A. Malanushenko, X. Sun, S. Régnier

Abstract: The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLF… ▽ More The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting. △ Less

Submitted 21 August, 2015; originally announced August 2015.

Comments: Accepted to ApJ; comments/corrections to this article are welcome via e-mail, even after publication

Nonassociative geometry in quasi-Hopf representation categories II: Connections and curvature

Authors: Gwendolyn E. Barnes, Alexander Schenkel, Richard J. Szabo

Abstract: We continue our systematic development of noncommutative and nonassociative differential geometry internal to the representation category of a quasitriangular quasi-Hopf algebra. We describe derivations, differential operators, differential calculi and connections using universal categorical constructions to capture algebraic properties such as Leibniz rules. Our main result is the construction of… ▽ More We continue our systematic development of noncommutative and nonassociative differential geometry internal to the representation category of a quasitriangular quasi-Hopf algebra. We describe derivations, differential operators, differential calculi and connections using universal categorical constructions to capture algebraic properties such as Leibniz rules. Our main result is the construction of morphisms which provide prescriptions for lifting connections to tensor products and to internal homomorphisms. We describe the curvatures of connections within our formalism, and also the formulation of Einstein-Cartan geometry as a putative framework for a nonassociative theory of gravity. △ Less

Submitted 2 May, 2016; v1 submitted 10 July, 2015; originally announced July 2015.

Comments: 29 pages. v2: Final version published in Journal of Geometry and Physics

Report number: EMPG-15-10 MSC Class: 16T05; 17B37; 46L87; 53D55

Journal ref: Journal of Geometry and Physics, Volume 106, August 2016, Pages 234-255

A Digital-Receiver for the Murchison Widefield Array

Authors: Thiagaraj Prabu, K. S. Srivani, D. Anish Roshi, P. A. Kamini, S. Madhavi, David Emrich, Brian Crosse, Andrew J. Williams, Mark Waterson, Avinash A. Deshpande, N. Udaya Shankar, Ravi Subrahmanyan, Frank H. Briggs, Robert F. Goeke, Steven J. Tingay, Melanie Johnston-Hollitt, Gopalakrishna M R, Edward H. Morgan, Joseph Pathikulangara, John D. Bunton, Grant Hampson, Christopher Williams, Stephen M. Ord, Randall B. Wayth, Deepak Kumar , et al. (33 additional authors not shown)

Abstract: An FPGA-based digital-receiver has been developed for a low-frequency imaging radio interferometer, the Murchison Widefield Array (MWA). The MWA, located at the Murchison Radio-astronomy Observatory (MRO) in Western Australia, consists of 128 dual-polarized aperture-array elements (tiles) operating between 80 and 300\,MHz, with a total processed bandwidth of 30.72 MHz for each polarization. Radio-… ▽ More An FPGA-based digital-receiver has been developed for a low-frequency imaging radio interferometer, the Murchison Widefield Array (MWA). The MWA, located at the Murchison Radio-astronomy Observatory (MRO) in Western Australia, consists of 128 dual-polarized aperture-array elements (tiles) operating between 80 and 300\,MHz, with a total processed bandwidth of 30.72 MHz for each polarization. Radio-frequency signals from the tiles are amplified and band limited using analog signal conditioning units; sampled and channelized by digital-receivers. The signals from eight tiles are processed by a single digital-receiver, thus requiring 16 digital-receivers for the MWA. The main function of the digital-receivers is to digitize the broad-band signals from each tile, channelize them to form the sky-band, and transport it through optical fibers to a centrally located correlator for further processing. The digital-receiver firmware also implements functions to measure the signal power, perform power equalization across the band, detect interference-like events, and invoke diagnostic modes. The digital-receiver is controlled by high-level programs running on a single-board-computer. This paper presents the digital-receiver design, implementation, current status, and plans for future enhancements. △ Less

Submitted 13 February, 2015; originally announced February 2015.

Comments: 14 pages, 7 figures

An In-Depth Spectroscopic Analysis of RR Lyr Variations over the Pulsation Cycle

Authors: L. Fossati, K. Kolenberg, D. V. Shulyak, A. Elmasli, V. Tsymbal, T. G. Barnes, E. Guggenberger, O. Kochukhov

Abstract: The stellar parameters of RR Lyrae stars vary considerably over a pulsation cycle, and their determination is crucial for stellar modelling. We present a detailed spectroscopic analysis of the pulsating star RR Lyr, the prototype of its class, over a complete pulsation cycle, based on high-resolution spectra collected at the 2.7-m telescope of McDonald Observatory. We used simultaneous photometry… ▽ More The stellar parameters of RR Lyrae stars vary considerably over a pulsation cycle, and their determination is crucial for stellar modelling. We present a detailed spectroscopic analysis of the pulsating star RR Lyr, the prototype of its class, over a complete pulsation cycle, based on high-resolution spectra collected at the 2.7-m telescope of McDonald Observatory. We used simultaneous photometry to determine the accurate pulsation phase of each spectrum and determined the effective temperature, the shape of the depth-dependent microturbulent velocity, and the abundance of several elements, for each phase. The surface gravity was fixed to 2.4. Element abundances resulting from our analysis are stable over the pulsation cycle. However, a variation in ionisation equilibrium is observed around minimum radius. We attribute this mostly to a dynamical acceleration contributing to the surface gravity. Variable turbulent convection on time scales longer than the pulsation cycle has been proposed as a cause for the Blazhko effect. We test this hypothesis to some extent by using the derived variable depth-dependent microturbulent velocity profiles to estimate their effect on the stellar magnitude. These effects turn out to be wavelength-dependent and much smaller than the observed light variations over the Blazhko cycle: if variations in the turbulent motions are entirely responsible for the Blazhko effect, they must surpass the scales covered by the microturbulent velocity. This work demonstrates the possibility of a self-consistent spectroscopic analysis over an entire pulsation cycle using static atmosphere models, provided one takes into account certain features of a rapidly pulsating atmosphere. △ Less

Submitted 16 October, 2014; originally announced October 2014.

Comments: Accepted for publication on MNRAS

Nonassociative geometry in quasi-Hopf representation categories I: Bimodules and their internal homomorphisms

Authors: Gwendolyn E. Barnes, Alexander Schenkel, Richard J. Szabo

Abstract: We systematically study noncommutative and nonassociative algebras A and their bimodules as algebras and bimodules internal to the representation category of a quasitriangular quasi-Hopf algebra. We enlarge the morphisms of the monoidal category of A-bimodules by internal homomorphisms, and describe explicitly their evaluation and composition morphisms. For braided commutative algebras A the full… ▽ More We systematically study noncommutative and nonassociative algebras A and their bimodules as algebras and bimodules internal to the representation category of a quasitriangular quasi-Hopf algebra. We enlarge the morphisms of the monoidal category of A-bimodules by internal homomorphisms, and describe explicitly their evaluation and composition morphisms. For braided commutative algebras A the full subcategory of symmetric A-bimodule objects is a braided closed monoidal category, from which we obtain an internal tensor product operation on internal homomorphisms. We describe how these structures deform under cochain twisting of the quasi-Hopf algebra, and apply the formalism to the example of deformation quantization of equivariant vector bundles over a smooth manifold. Our constructions set up the basic ingredients for the systematic development of differential geometry internal to the quasi-Hopf representation category, which will be tackled in the sequels to this paper, together with applications to models of noncommutative and nonassociative gravity such as those anticipated from non-geometric string theory. △ Less

Submitted 6 February, 2015; v1 submitted 22 September, 2014; originally announced September 2014.

Comments: v1: 56 pages. v2: 57 pages. Introduction expanded. Published in Journal of Geometry and Physics

Report number: EMPG-14-17 MSC Class: 16T05; 17B37; 46L87; 53D55

Journal ref: Journal of Geometry and Physics, Volume 89, March 2015, Pages 111-152

Using Regression Kernels to Forecast A Failure to Appear in Court

Authors: Richard Berk, Justin Bleich, Adam Kapelner, Jaime Henderson, Geoffrey Barnes, Ellen Kurtz

Abstract: Forecasts of prospective criminal behavior have long been an important feature of many criminal justice decisions. There is now substantial evidence that machine learning procedures will classify and forecast at least as well, and typically better, than logistic regression, which has to date dominated conventional practice. However, machine learning procedures are adaptive. They "learn" inductivel… ▽ More Forecasts of prospective criminal behavior have long been an important feature of many criminal justice decisions. There is now substantial evidence that machine learning procedures will classify and forecast at least as well, and typically better, than logistic regression, which has to date dominated conventional practice. However, machine learning procedures are adaptive. They "learn" inductively from training data. As a result, they typically perform best with very large datasets. There is a need, therefore, for forecasting procedures with the promise of machine learning that will perform well with small to moderately-sized datasets. Kernel methods provide precisely that promise. In this paper, we offer an overview of kernel methods in regression settings and compare such a method, regularized with principle components, to stepwise logistic regression. We apply both to a timely and important criminal justice concern: a failure to appear (FTA) at court proceedings following an arraignment. A forecast of an FTA can be an important factor is a judge's decision to release a defendant while awaiting trial and can influence the conditions imposed on that release. Forecasting accuracy matters, and our kernel approach forecasts far more accurately than stepwise logistic regression. The methods developed here are implemented in the R package kernReg currently available on CRAN. △ Less

Submitted 5 September, 2014; originally announced September 2014.

Comments: 43 pages, 5 figures, 2 tables, 2 appendices

Echography of young stars reveals their evolution

Authors: K. Zwintz, L. Fossati, T. Ryabchikova, D. Guenther, C. Aerts, T. G. Barnes, N. Themessl, D. Lorenz, C. Cameron, R. Kuschnig, S. Pollack-Drs, E. Moravveji, A. Baglin, J. M. Matthews, A. F. J. Moffat, E. Poretti, M. Rainer, S. M. Rucinski, D. Sasselov, W. W. Weiss

Abstract: We demonstrate that a seismic analysis of stars in their earliest evolutionary phases is a powerful method to identify young stars and distinguish their evolutionary states. The early star that is born from the gravitational collapse of a molecular cloud reaches at some point sufficient temperature, mass and luminosity to be detected. Accretion stops and the pre-main sequence star that emerges is… ▽ More We demonstrate that a seismic analysis of stars in their earliest evolutionary phases is a powerful method to identify young stars and distinguish their evolutionary states. The early star that is born from the gravitational collapse of a molecular cloud reaches at some point sufficient temperature, mass and luminosity to be detected. Accretion stops and the pre-main sequence star that emerges is nearly fully convective and chemically homogeneous. It will continue to contract gravitationally until the density and temperature in the core are high enough to start nuclear burning of hydrogen. We show that there is a relationship between detected pulsation properties for a sample of young stars and their evolutionary status illustrating the potential of asteroseismology for the early evolutionary phases. △ Less

Submitted 18 July, 2014; originally announced July 2014.

Comments: 23 pages, 9 figures, includes Supplementary Material. Science, published in Science Express on July 3, 2014

The Helioseismic and Magnetic Imager (HMI) Vector Magnetic Field Pipeline: Overview and Performance

Authors: J. Todd Hoeksema, Yang Liu, Keiji Hayashi, Xudong Sun, Jesper Schou, Sebastien Couvidat, Aimee Norton, Monica Bobra, Rebecca Centeno, K. D. Leka, Graham Barnes, Michael J. Turmon

Abstract: The Helioseismic and Magnetic Imager (HMI) began near-continuous full-disk solar measurements on 1 May 2010 from the Solar Dynamics Observatory (SDO). An automated processing pipeline keeps pace with observations to produce observable quantities, including the photospheric vector magnetic field, from sequences of filtergrams. The primary 720s observables were released in mid 2010, including Stokes… ▽ More The Helioseismic and Magnetic Imager (HMI) began near-continuous full-disk solar measurements on 1 May 2010 from the Solar Dynamics Observatory (SDO). An automated processing pipeline keeps pace with observations to produce observable quantities, including the photospheric vector magnetic field, from sequences of filtergrams. The primary 720s observables were released in mid 2010, including Stokes polarization parameters measured at six wavelengths as well as intensity, Doppler velocity, and the line-of-sight magnetic field. More advanced products, including the full vector magnetic field, are now available. Automatically identified HMI Active Region Patches (HARPs) track the location and shape of magnetic regions throughout their lifetime. The vector field is computed using the Very Fast Inversion of the Stokes Vector (VFISV) code optimized for the HMI pipeline; the remaining 180 degree azimuth ambiguity is resolved with the Minimum Energy (ME0) code. The Milne-Eddington inversion is performed on all full-disk HMI observations. The disambiguation, until recently run only on HARP regions, is now implemented for the full disk. Vector and scalar quantities in the patches are used to derive active region indices potentially useful for forecasting; the data maps and indices are collected in the SHARP data series, hmi.sharp_720s. Patches are provided in both CCD and heliographic coordinates. HMI provides continuous coverage of the vector field, but has modest spatial, spectral, and temporal resolution. Coupled with limitations of the analysis and interpretation techniques, effects of the orbital velocity, and instrument performance, the resulting measurements have a certain dynamic range and sensitivity and are subject to systematic errors and uncertainties that are characterized in this report. △ Less

Submitted 7 April, 2014; originally announced April 2014.

Comments: 42 pages, 19 figures, accepted to Solar Physics

The Helioseismic and Magnetic Imager (HMI) Vector Magnetic Field Pipeline: SHARPs -- Space-weather HMI Active Region Patches

Authors: Monica G. Bobra, Xudong Sun, J. Todd Hoeksema, Michael J. Turmon, Yang Liu, Keiji Hayashi, Graham Barnes, K. D. Leka

Abstract: A new data product from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) called Space-weather HMI Active Region Patches (SHARPs) is now available. SDO/HMI is the first space-based instrument to map the full-disk photospheric vector magnetic field with high cadence and continuity. The SHARP data series provide maps in patches that encompass automatically track… ▽ More A new data product from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) called Space-weather HMI Active Region Patches (SHARPs) is now available. SDO/HMI is the first space-based instrument to map the full-disk photospheric vector magnetic field with high cadence and continuity. The SHARP data series provide maps in patches that encompass automatically tracked magnetic concentrations for their entire lifetime; map quantities include the photospheric vector magnetic field and its uncertainty, along with Doppler velocity, continuum intensity, and line-of-sight magnetic field. Furthermore, keywords in the SHARP data series provide several parameters that concisely characterize the magnetic-field distribution and its deviation from a potential-field configuration. These indices may be useful for active-region event forecasting and for identifying regions of interest. The indices are calculated per patch and are available on a twelve-minute cadence. Quick-look data are available within approximately three hours of observation; definitive science products are produced approximately five weeks later. SHARP data are available at http://jsoc.stanford.edu and maps are available in either of two different coordinate systems. This article describes the SHARP data products and presents examples of SHARP data and parameters. △ Less

Submitted 7 April, 2014; originally announced April 2014.

Comments: 27 pages, 7 figures. Accepted to Solar Physics

The Helioseismic and Magnetic Imager (HMI) Vector Magnetic Field Pipeline: Optimization of the Spectral Line Inversion Code

Authors: R. Centeno, J. Schou, K. Hayashi, A. Norton, J. T. Hoeksema, Y. Liu, K. D. Leka, G. Barnes

Abstract: The Very Fast Inversion of the Stokes Vector (VFISV) is a Milne-Eddington spectral line inversion code used to determine the magnetic and thermodynamic parameters of the solar photosphere from observations of the Stokes vector in the 6173 A Fe I line by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We report on the modifications made to the original VFISV… ▽ More The Very Fast Inversion of the Stokes Vector (VFISV) is a Milne-Eddington spectral line inversion code used to determine the magnetic and thermodynamic parameters of the solar photosphere from observations of the Stokes vector in the 6173 A Fe I line by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We report on the modifications made to the original VFISV inversion code in order to optimize its operation within the HMI data pipeline and provide the smoothest solution in active regions. The changes either sped up the computation or reduced the frequency with which the algorithm failed to converge to a satisfactory solution. Additionally, coding bugs which were detected and fixed in the original VFISV release, are reported here. △ Less

Submitted 14 March, 2014; originally announced March 2014.

Comments: Accepted for publication in Solar Physics