Identifying invertebrate invasions using morphological and molecular analyses: North American Daphnia ‘pulex’ in New Zealand fresh waters

We used a DNA barcoding approach to identify specimens of the Daphnia pulex complex occurring in New Zealand lakes, documenting the establishment of non-indigenous North American Daphnia 'pulex'. Morphological delineation of species in this complex is problematic due to a lack of good morphological traits to distinguish the species, as there is a relatively high degree of morphological stasis within the group through evolutionary time. Accordingly, genetic analyses were used to determine the specific identity and likely geographic origin of this species. Morphologically, individuals most closely resembled Daphnia pulicaria or Daphnia pulex sensu lato, which cannot be separated morphologically. Furthermore, each of these taxa comprises separate species in North America and Europe, despite carrying the same names. We identified individuals using a 658 bp nucleotide portion of the mitochondrial cytochrome c oxidase subunit 1 gene (COI) as North American Daphnia 'pulex', being distinct from European Daphnia pulex sensu stricto and D. pulicaria from Europe or North America. Cellulose allozyme electrophoresis was used to confirm that individuals were not hybrids with D. pulicaria. North American Daphnia 'pulex' in New Zealand were first recorded in New Zealand from South Island lakes that are popular for overseas recreational fishers, indicating a possible source of introduction for this species (e.g. on/in fishing gear). Our study provides an additional example of how genetic techniques can be used for the accurate identification of non-indigenous taxa, particularly when morphological species determination is not possible. The growth of global databases such as GenBank and Barcode of Life Datasystems (BOLD) will further enhance this identification capacity

Milky Way Mass and Potential Recovery Using Tidal Streams in a Realistic Halo

We present a new method for determining the Galactic gravitational potential based on forward modeling of tidal stellar streams. We use this method to test the performance of smooth and static analytic potentials in representing realistic dark matter halos, which have substructure and are continually evolving by accretion. Our FAST-FORWARD method uses a Markov Chain Monte Carlo algorithm to compare, in 6D phase space, an "observed" stream to models created in trial analytic potentials. We analyze a large sample of streams evolved in the Via Lactea II (VL2) simulation, which represents a realistic Galactic halo potential. The recovered potential parameters are in agreement with the best fit to the global, present-day VL2 potential. However, merely assuming an analytic potential limits the dark matter halo mass measurement to an accuracy of 5 to 20%, depending on the choice of analytic parametrization. Collectively, mass estimates using streams from our sample reach this fundamental limit, but individually they can be highly biased. Individual streams can both under- and overestimate the mass, and the bias is progressively worse for those with smaller perigalacticons, motivating the search for tidal streams at galactocentric distances larger than 70 kpc. We estimate that the assumption of a static and smooth dark matter potential in modeling of the GD-1 and Pal5-like streams introduces an error of up to 50% in the Milky Way mass estimates.Comment: 12 pages, 6 figures, submitted to ApJ; more information on our stream sample and a movie of the potential recovery method used can be found at http://www.astro.yale.edu/abonaca/research/potential_recovery.htm

Steady-state simulations using weighted ensemble path sampling

We extend the weighted ensemble (WE) path sampling method to perform rigorous statistical sampling for systems at steady state. The straightforward steady-state implementation of WE is directly practical for simple landscapes, but not when significant metastable intermediates states are present. We therefore develop an enhanced WE scheme, building on existing ideas, which accelerates attainment of steady state in complex systems. We apply both WE approaches to several model systems confirming their correctness and efficiency by comparison with brute-force results. The enhanced version is significantly faster than the brute force and straightforward WE for systems with WE bins that accurately reflect the reaction coordinate(s). The new WE methods can also be applied to equilibrium sampling, since equilibrium is a steady state

A deductive statistical mechanics approach for granular matter

We introduce a deductive statistical mechanics approach for granular materials which is formally built from few realistic physical assumptions. The main finding is an universal behavior for the distribution of the density fluctuations. Such a distribution is the equivalent of the Maxwell-Boltzmann's distribution in the kinetic theory of gasses. The comparison with a very extensive set of experimental and simulation data for packings of monosized spherical grains, reveals a remarkably good quantitative agreement with the theoretical predictions for the density fluctuations both at the grain level and at the global system level. Such agreement is robust over a broad range of packing fractions and it is observed in several distinct systems prepared by using different methods. The equilibrium distributions are characterized by only one parameter ($k$) which is a quantity very sensitive to changes in the structural organization. The thermodynamical equivalent of $k$ and its relation with the `granular temperature' are also discussed.Comment: 15 pages, 6 figure

Some Systematics of Galactic Globular Clusters

The global properties of all known Galactic globular clusters are examined. The relationship between the luminosities and the metallicities of Galactic globular clusters is found to be complex. Among luminous clusters there is a correlation in the sense that the oldest clusters are slightly more metal deficient than are younger clusters. However, no such clear-cut relationship is found among the faintest globular clusters. The central concentration index C of globular clusters is seen to be independent of metallicity. The dependence of the half-light radii of globular clusters on their Galactocentric distances can be approximated by the relation $R_h \alpha R^{2/3}_{gc}$. Clusters with collapsed cores are mostly situated close to the Galactic nucleus. For $R_{gc} < 10$ kpc the luminosities and the radii of clusters appear to be uncorrelated. The Galaxy differs from the LMC and the SMC in that it appears to lack highly flattened luminous clusters. Galactic globular clusters with ages $\geq$ 13.0 Gyr are all of Oosterhoff type II, whereas almost all of those with ages $<$ 13.0 Gyr have been assigned to Oosterhoff type I. Globular clusters with ages $<$11.5 Gyr are all located in the outer Galactic halo, have below-average luminosities and above-average radii. On the other hand the very old globular cluster NGC 6522 is situated close to the Galactic nucleus.Comment: PASP, in pres

Single-shot single-gate RF spin readout in silicon

For solid-state spin qubits, single-gate RF readout can help minimise the number of gates required for scale-up to many qubits since the readout sensor can integrate into the existing gates required to manipulate the qubits (Veldhorst 2017, Pakkiam 2018). However, a key requirement for a scalable quantum computer is that we must be capable of resolving the qubit state within single-shot, that is, a single measurement (DiVincenzo 2000). Here we demonstrate single-gate, single-shot readout of a singlet-triplet spin state in silicon, with an average readout fidelity of $82.9\%$ at a $3.3~\text{kHz}$ measurement bandwidth. We use this technique to measure a triplet $T_-$ to singlet $S_0$ relaxation time of $0.62~\text{ms}$ in precision donor quantum dots in silicon. We also show that the use of RF readout does not impact the maximum readout time at zero detuning limited by the $S_0$ to $T_-$ decay, which remained at approximately $2~\text{ms}$. This establishes single-gate sensing as a viable readout method for spin qubits

Globular Cluster Streams as Galactic High-Precision Scales - The Poster Child Palomar 5

Using the example of the tidal stream of the Milky Way globular cluster Palomar 5 (Pal 5), we demonstrate how observational data on streams can be efficiently reduced in dimensionality and modeled in a Bayesian framework. Our approach combines detection of stream overdensities by a Difference-of-Gaussians process with fast streakline models, a continuous likelihood function built from these models, and inference with MCMC. By generating $\approx10^7$ model streams, we show that the geometry of the Pal 5 debris yields powerful constraints on the solar position and motion, the Milky Way and Pal 5 itself. All 10 model parameters were allowed to vary over large ranges without additional prior information. Using only SDSS data and a few radial velocities from the literature, we find that the distance of the Sun from the Galactic Center is $8.30\pm0.25$ kpc, and the transverse velocity is $253\pm16$ km/s. Both estimates are in excellent agreement with independent measurements of these quantities. Assuming a standard disk and bulge model, we determine the Galactic mass within Pal 5's apogalactic radius of 19 kpc to be $(2.1\pm0.4)\times10^{11}$ M$_\odot$. Moreover, we find the potential of the dark halo with a flattening of $q_z = 0.95^{+0.16}_{-0.12}$ to be essentially spherical within the radial range that is effectively probed by Pal 5. We also determine Pal 5's mass, distance and proper motion independently from other methods, which enables us to perform vital cross-checks. We conclude that with more observational data and by using additional prior information, the precision of this method can be significantly increased.Comment: 28 pages, 14 figures, submitted to ApJ (revised version), comments welcom