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

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.Comment: J. Space Weather Space Climate: Accepted / in press; access supplementary materials through journal; some figures are less than full resolution for arXi

Achieving Consistent Doppler Measurements from SDO/HMI Vector Field Inversions

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.Comment: 58 pages, 19 figures, submitted to Ap

BB Intermeson Potentials in the Quark Model

In this paper we derive quark model results for scattering amplitudes and equivalent low energy potentials for heavy meson pairs, in which each meson contains a heavy quark. This "BB" system is an attractive theoretical laboratory for the study of the nuclear force between color singlets; the hadronic system is relatively simple, and there are lattice gauge theory (LGT) results for V_BB(r) which may be compared to phenomenological models. We find that the quark model potential (after lattice smearing) has qualitative similarities to the LGT potential in the two B*B* channels in which direct comparison is possible, although there is evidence of a difference in length scales. The quark model prediction of equal magnitude but opposite sign for I=0 and I=1 potentials also appears similar to LGT results at intermediate r. There may however be a discrepancy between the LGT and quark model I=1 BB potentials. A numerical study of the two-meson Schrodinger equations in the (bqbar)(bqbar) and (cqbar)(cqbar) sectors with the quark model potentials finds a single "molecule", in the I=0 BB* sector. Binding in other channels might occur if the quark model forces are augmented by pion exchange.Comment: 30 pages, 5 figures, revtex and epsfig. Submitted to Phys. Rev.

Helioseismology of Pre-Emerging Active Regions II: Average Emergence Properties

We report on average subsurface properties of pre-emerging active regions as compared to areas where no active region emergence was detected. Helioseismic holography is applied to samples of the two populations (pre-emergence and without emergence), each sample having over 100 members, which were selected to minimize systematic bias, as described in Leka et al. We find that there are statistically significant signatures (i.e., difference in the means of more than a few standard errors) in the average subsurface flows and the apparent wave speed that precede the formation of an active region. The measurements here rule out spatially extended flows of more than about 15 m/s in the top 20 Mm below the photosphere over the course of the day preceding the start of visible emergence. These measurements place strong constraints on models of active region formation.Comment: 15 pages, 10 figures, ApJ (published

Spectroscopic membership for the populous 300 Myr-old open cluster NGC 3532

NGC 3532 is an extremely rich open cluster embedded in the Galactic disc, hitherto lacking a comprehensive, documented membership list. We provide membership probabilities from new radial velocity observations of solar-type and low-mass stars in NGC 3532, in part as a prelude to a subsequent study of stellar rotation in the cluster. Using extant optical and infra-red photometry we constructed a preliminary photometric membership catalogue, consisting of 2230 dwarf and turn-off stars. We selected 1060 of these for observation with the AAOmega spectrograph at the Anglo-Australian Telescope and 391 stars for observations with the Hydra-South spectrograph at the Victor Blanco Telescope, obtaining spectroscopic observations over a decade for 145 stars. We measured radial velocities for our targets through cross-correlation with model spectra and standard stars, and supplemented them with radial velocities for 433 additional stars from the literature. We also measured log g, Teff, and [Fe/H] from the AAOmega spectra. Together with proper motions from Gaia DR2 we find 660 exclusive members. The members are distributed across the whole cluster sequence, from giant stars to M dwarfs, making NGC 3532 one of the richest Galactic open clusters known to date, on par with the Pleiades. From further spectroscopic analysis of 153 dwarf members we find the metallicity to be marginally sub-solar, with [Fe/H]=-0.07. Exploiting trigonometric parallax measurements from Gaia DR2 we find a distance of $484^{+35}_{-30}$ pc. Based on the membership we provide an empirical cluster sequence in multiple photometric passbands. A comparison of the photometry of the measured cluster members with several recent model isochrones enables us to confirm the 300 Myr cluster age. However, all of the models evince departures from the cluster sequence in particular regions, especially in the lower mass range. (abridged)Comment: Accepted for publication in A&A. 19 pages, 18 Figures, and 6 Table

Dissociation cross sections of ground-state and excited charmonia with light mesons in the quark model

We present numerical results for the dissociation cross sections of ground-state, orbitally- and radially-excited charmonia in collisions with light mesons. Our results are derived using the nonrelativistic quark model, so all parameters are determined by fits to the experimental meson spectrum. Examples of dissociation into both exclusive and inclusive final states are considered. The dissociation cross sections of several C=(+) charmonia may be of considerable importance for the study of heavy ion collisions, since these states are expected to be produced more copiously than the J/psi. The relative importance of the productions of ground-state and orbitally-excited charmed mesons in a pion-charmonium collision is demonstrated through the $\sqrt {s}$-dependent charmonium dissociation cross sections.Comment: 9 pages, 8 figure

An HI census of Loose Groups of Galaxies

We present results from our Parkes Multibeam HI survey of 3 loose groups of galaxies that are analogous to the Local Group. This is a survey of groups containing only spiral galaxies with mean separations of a few hundred kpc, and total areas of approximately 1 sq. Mpc; groups similar to our own Local Group. We present a census of the HI-rich objects in these groups down to an M(HI), 1-sigma sensitivity ~7x10^5 M(sun), as well as the detailed properties of these detections from follow-up Compact Array observations. We found 7 new HI-rich members in the 3 groups, all of which have stellar counterparts and are, therefore, typical dwarf galaxies. The ratio of low-mass to high-mass gas-rich galaxies in these groups is less than in the Local Group meaning that the ``missing satellite'' problem is not unique. No high-velocity cloud analogs were found in any of the groups. If HVCs in these groups are the same as in the Local Group, this implies that HVCs must be located within ~300-400 kpc of the Milky Way.Comment: 6 pages, to appear in the ASP proceedings of IAU Symposium 217, "Recycling intergalactic and interstellar matter", eds. Pierre-Alain Duc, Jonathan Braine, Elias Brink

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

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.Comment: Accepted for publication in Solar Physic

Bi-layer splitting in overdoped high TcT_{c} cuprates

Recent angle-resolved photoemission data for overdoped Bi2212 are explained. Of the peak-dip-hump structure, the peak corresponds the $\vec q =0$ component of a hole condensate which appears at $T_c$. The fluctuating part of this same condensate produces the hump. The bilayer splitting is large enough to produce a bonding hole and an electron antibonding quasiparticle Fermi surface. Smaller bilayer splittings observed in some experiments reflect the interaction of the peak structure with quasiparticle states near, but not at, the Fermi surface.Comment: 4 pages with 2 figures - published versio