
The Astropy project supports and fosters the development of opensource and
openlydeveloped Python packages that provide commonlyneeded functionality to
the astronomical community. A key element of the Astropy project is the core
package Astropy, which serves as the foundation for more specialized projects
and packages. In this article, we provide an overview of the organization of
the Astropy project and summarize key features in the core package as of the
recent major release, version 2.0. We then describe the project infrastructure
designed to facilitate and support development for a broader ecosystem of
interoperable packages. We conclude with a future outlook of planned new
features and directions for the broader Astropy project.

The DECam Plane Survey is a fiveband optical and nearinfrared survey of the
southern Galactic plane with the Dark Energy Camera at Cerro Tololo. The survey
is designed to reach past the mainsequence turnoff at the distance of the
Galactic center through a reddening E(BV) of 1.5 mag. Typical singleexposure
depths are 23.7, 22.8, 22.3, 21.9, and 21.0 mag in the grizY bands, with seeing
around 1 arcsecond. The footprint covers the Galactic plane with b < 4
degrees, 5 degrees > l > 120 degrees. The survey pipeline simultaneously
solves for the positions and fluxes of tens of thousands of sources in each
image, delivering positions and fluxes of roughly two billion stars with better
than 10 mmag precision. Most of these objects are highly reddened and deep in
the Galactic disk, probing the structure and properties of the Milky Way and
its interstellar medium. The full survey is publicly available.

The Long Baseline Neutrino Facility (LBNF) project will build a beamline
located at Fermilab to create and aim an intense neutrino beam of appropriate
energy range toward the DUNE detectors at the SURF facility in Lead, South
Dakota. Neutrino production starts in the Target Station, which consists of a
solid target, magnetic focusing horns, and the associated subsystems and
shielding infrastructure. Protons hit the target producing mesons which are
then focused by the horns into a heliumfilled decay pipe where they decay into
muons and neutrinos. The target and horns are encased in actively cooled steel
and concrete shielding in a chamber called the target chase. The reference
design chase is filled with air, but nitrogen and helium are being evaluated as
alternatives. A replaceable beam window separates the decay pipe from the
target chase. The facility is designed for initial operation at 1.2 MW, with
the ability to upgrade to 2.4 MW, and is taking advantage of the experience
gained by operating Fermilab's NuMI facility. We discuss here the design
status, associated challenges, and ongoing R&D and physicsdriven component
optimization of the Target Station.

We study URu$_{2x}$Fe$_x$Si$_2$, in which two types of staggered phases
compete at low temperature as the iron concentration $x$ is varied: the
nonmagnetic "hidden order" (HO) phase below the critical concentration $x_c$,
and unconventional antiferromagnetic (AF) phase above $x_c$. By using
polarization resolved Raman spectroscopy, we detect a collective mode of
pseudovectorlike $A_{2g}$ symmetry whose energy continuously evolves with
increasing $x$; it monotonically decreases in the HO phase until it vanishes at
$x=x_c$, and then reappears with increasing energy in the AF phase. The mode's
evolution provides direct evidence for unified order parameter for both
nonmagnetic and magnetic phases arising from the orbital degreesoffreedom of
the uranium5$f$ electrons.

We present a matching and LP based heuristic algorithm that decides graph
nonHamiltonicity. Each of the $n!$ Hamilton cycles in a complete directed
graph on $n+1$ vertices corresponds with each of the $n!$ $n$permutation
matrices $P$, such that $p_{u,i}=1$ if and only if the $i^{th}$ arc in a cycle
enters vertex $u$, starting and ending at vertex $n+1$. A graph instance ($G$)
is initially coded as exclusion set $E$, whose members are pairs of components
of $P$, $\{p_{u,i} ,p_{v,i+1}\}, i=1,n1$, for each arc $(u,v)$ not in $G$. For
each $\{p_{u,i} ,p_{v,i+1}\}\in E$, the set of $P$ satisfying
$p_{u,i}=p_{v,i+1}=1$ correspond with a set of cycles not in $G$. Accounting
for all arcs not in $G$, $E$ codes precisely the set of cycles not in $G$. A
doubly stochasticlike $\mathcal{O}$($n^4$) formulation of the Hamilton cycle
decision problem is then constructed. Each $\{p_{u,i} ,p_{v,j}\}$ is coded as
variable $q_{u,i,v,j}$ such that the set of integer extrema is the set of all
permutations. We model $G$ by setting each $q_{u,i,v,j}=0$ in correspondence
with each $\{p_{u,i} ,p_{v,j}\}\in E$ such that for nonHamiltonian $G$,
integer solutions cannot exist. We recognize nonHamiltonicity by iteratively
deducing additional $q_{u,i,v,j}$ that can be set zero and expanding $E$ until
the formulation becomes infeasible, in which case we recognize that no integer
solutions exists i.e. $G$ is decided nonHamiltonian. Over 100 nonHamiltonian
graphs (10 through 104 vertices) and 2000 randomized 31 vertex nonHamiltonian
graphs are tested and correctly decided nonHamiltonian. For Hamiltonian $G$,
the complement of $E$ provides information about covers of matchings, perhaps
useful in searching for cycles. We also present an example where the algorithm
fails to deduce any integral value for any $q_{u,i,v,j}$ i.e. $G$ is undecided.

Spontaneous symmetry breaking is a fundamental concept in many areas of
physics, ranging from cosmology and particle physics to condensed matter. A
prime example is the breaking of spatial translation symmetry, which underlies
the formation of crystals and the phase transition from liquid to solid.
Analogous to crystals in space, the breaking of translation symmetry in time
and the emergence of a "time crystal" was recently proposed, but later shown to
be forbidden in thermal equilibrium. However, nonequilibrium Floquet systems
subject to a periodic drive can exhibit persistent timecorrelations at an
emergent subharmonic frequency. This new phase of matter has been dubbed a
"discrete time crystal" (DTC). Here, we present the first experimental
observation of a discrete time crystal, in an interacting spin chain of trapped
atomic ions. We apply a periodic Hamiltonian to the system under manybody
localization (MBL) conditions, and observe a subharmonic temporal response
that is robust to external perturbations. Such a time crystal opens the door
for studying systems with longrange spatialtemporal correlations and novel
phases of matter that emerge under intrinsically nonequilibrium conditions.

Random magnetization dynamics driven by a jumpnoise process is reduced to
stochastic magnetic energy dynamics on specific graphs using an averaging
technique. An approach to analyzing stochastic energy dynamics on graphs is
presented and applied to the calculation of power spectral density of random
magnetization dynamics. An eigenvalue technique for computing the power
spectral density under specific cases is also presented and illustrated by
numerical results.

We use electronic Raman scattering to study the lowenergy excitations in
BaFe$_2$(As$_{0.5}$P$_{0.5}$)$_2$ ($T_c \approx 16$ K) samples. In addition to
a superconducting pair breaking peak (2$\Delta=6.7$ meV) in the A$_{1g}$
channel with a linear tail towards zero energy, suggesting a nodal gap
structure, we detect spectral features associated to Pomeranchuk oscillations
in the A$_{1g}$, B$_{1g}$ and B$_{2g}$ channels. We argue that the small Fermi
energy of the system is an essential condition for these Pomeranchuk
oscillations to be underdamped. The Pomeranchuk oscillations have the same
frequencies in the B$_{1g}$ and B$_{2g}$ channels, which we explain by the
mixing of these symmetries resulting from the removal of the $\sigma_v$ and
$\sigma_v$ symmetry planes due to a large As/P disorder. Interestingly, we show
that the temperature at which the peaks corresponding to the Pomeranchuk
oscillations get underdamped is consistent with the nonFermi liquid to Femi
liquid crossover determined by transport, suggesting that the Pomeranchuk
instability plays an important role in the lowenergy physics of the Febased
superconductors.

We report results from the first search for $\nu_\mu\to\nu_e$ transitions by
the NOvA experiment. In an exposure equivalent to $2.74\times10^{20}$
protonsontarget in the upgraded NuMI beam at Fermilab, we observe 6 events in
the Far Detector, compared to a background expectation of $0.99\pm0.11$ (syst.)
events based on the Near Detector measurement. A secondary analysis observes 11
events with a background of $1.07\pm0.14$ (syst.). The $3.3\sigma$ excess of
events observed in the primary analysis disfavors $0.1\pi < \delta_{CP} <
0.5\pi$ in the inverted mass hierarchy at the 90% C.L.

This paper reports the first measurement using the NOvA detectors of
$\nu_\mu$ disappearance in a $\nu_\mu$ beam. The analysis uses a 14
ktonequivalent exposure of $2.74 \times 10^{20}$ protonsontarget from the
Fermilab NuMI beam. Assuming the normal neutrino mass hierarchy, we measure
$\Delta m^{2}_{32}=(2.52^{+0.20}_{0.18})\times 10^{3}$ eV$^{2}$ and
$\sin^2\theta_{23}$ in the range 0.380.65, both at the 68% confidence level,
with two statisticallydegenerate best fit points at $\sin^2\theta_{23} = $
0.43 and 0.60. Results for the inverted mass hierarchy are also presented.

The readout of transitionedge sensor (TES) bolometers with a large
multiplexing factor is key for the next generation Cosmic Microwave Background
(CMB) experiment, Polarbear2, having 7,588 TES bolometers. To enable the large
arrays, we have been developing a readout system with a multiplexing factor of
40 in the frequency domain. Extending that architecture to 40 bolometers
requires an increase in the bandwidth of the SQUID electronics above 4 MHz.
This paper focuses on cryogenic readout and shows how it affects cross talk and
the responsivity of the TES bolometers. A series resistance, such as equivalent
series resistance (ESR) of capacitors for LC filters, leads to nonlinear
response of the bolometers. A wiring inductance modulates a voltage across the
bolometers and causes cross talk. They should be controlled well to reduce
systematic errors in CMB observations. We have been developing a cryogenic
readout with a low series impedance and have tuned bolometers in the middle of
their transition at a high frequency (> 3 MHz).

We present an overview of the design and status of the \Pb2 and the Simons
Array experiments. \Pb2 is a Cosmic Microwave Background polarimetry
experiment which aims to characterize the arcminute angular scale Bmode
signal from weak gravitational lensing and search for the degree angular scale
Bmode signal from inflationary gravitational waves. The receiver has a 365~mm
diameter focal plane cooled to 270~milliKelvin. The focal plane is filled with
7,588 dichroic lensletantenna coupled polarization sensitive Transition Edge
Sensor (TES) bolometric pixels that are sensitive to 95~GHz and 150~GHz bands
simultaneously. The TES bolometers are readout by SQUIDs with 40 channel
frequency domain multiplexing. Refractive optical elements are made with high
purity alumina to achieve high optical throughput. The receiver is designed to
achieve noise equivalent temperature of 5.8~$\mu$K$_{CMB}\sqrt{s}$ in each
frequency band. \Pb2 will deploy in 2016 in the Atacama desert in Chile. The
Simons Array is a project to further increase sensitivity by deploying three
\Pb2 type receivers. The Simons Array will cover 95~GHz, 150~GHz and 220~GHz
frequency bands for foreground control. The Simons Array will be able to
constrain tensortoscalar ratio and sum of neutrino masses to $\sigma(r) =
6\times 10^{3}$ at $r = 0.1$ and $\sum m_\nu (\sigma =1)$ to 40 meV.

The physics case and an experimental overview of the MOLLER (Measurement Of a
Lepton Lepton Electroweak Reaction) experiment at the 12 GeV upgraded Jefferson
Lab are presented. A highlight of the Fundamental Symmetries subfield of the
2007 NSAC Long Range Plan was the SLAC E158 measurement of the parityviolating
asymmetry $A_{PV}$ in polarized electronelectron (M{\o}ller) scattering. The
proposed MOLLER experiment will improve on this result by a factor of five,
yielding the most precise measurement of the weak mixing angle at low or high
energy anticipated over the next decade. This new result would be sensitive to
the interference of the electromagnetic amplitude with new neutral current
amplitudes as weak as $\sim 10^{3}\cdot G_F$ from as yet undiscovered dynamics
beyond the Standard Model. The resulting discovery reach is unmatched by any
proposed experiment measuring a flavor and CPconserving process over the next
decade, and yields a unique window to new physics at MeV and multiTeV scales,
complementary to direct searches at high energy colliders such as the Large
Hadron Collider (LHC). The experiment takes advantage of the unique opportunity
provided by the upgraded electron beam energy, luminosity, and stability at
Jefferson Laboratory and the extensive experience accumulated in the community
after a round of recent successfully completed parityviolating electron
scattering experiments

The physics of interacting integerspin chains has been a topic of intense
theoretical interest, particularly in the context of symmetryprotected
topological phases. However, there has not been a controllable model system to
study this physics experimentally. We demonstrate how spindependent forces on
trapped ions can be used to engineer an effective system of interacting spin1
particles. Our system evolves coherently under an applied spin1 XY Hamiltonian
with tunable, longrange couplings, and all three quantum levels at each site
participate in the dynamics. We observe the time evolution of the system and
verify its coherence by entangling a pair of effective threelevel particles
(`qutrits') with 86% fidelity. By adiabatically ramping a global field, we
produce ground states of the XY model, and we demonstrate an instance where the
ground state cannot be created without breaking the same symmetries that
protect the topological Haldane phase. This experimental platform enables
future studies of symmetryprotected order in spin1 systems and their use in
quantum applications.

Measurements of neutral current cross sections for deep inelastic scattering
in e+p collisions at HERA with a longitudinally polarised positron beam are
presented. The singledifferential crosssections d(sigma)/dQ2, d(sigma)/dx and
d(sigma)/dy and the reduced crosssection were measured in the kinematic region
Q2 > 185 GeV2 and y < 0.9, where Q2 is the fourmomentum transfer squared, x
the Bjorken scaling variable, and y the inelasticity of the interaction. The
measurements were performed separately for positively and negatively polarised
positron beams. The measurements are based on an integrated luminosity of 135.5
pb1 collected with the ZEUS detector in 2006 and 2007 at a centreofmass
energy of 318 GeV. The structure functions F3 and F3(gamma)Z were determined by
combining the e+p results presented in this paper with previously published ep
neutral current results. The asymmetry parameter A+ is used to demonstrate the
parity violation predicted in electroweak interactions. The measurements are
well described by the predictions of the Standard Model.

Quantum simulators, in which well controlled quantum systems are used to
reproduce the dynamics of less understood ones, have the potential to explore
physics that is inaccessible to modeling with classical computers. However,
checking the results of such simulations will also become classically
intractable as system sizes increase. In this work, we introduce and implement
a coherent imaging spectroscopic technique to validate a quantum simulation,
much as magnetic resonance imaging exposes structure in condensed matter. We
use this method to determine the energy levels and interaction strengths of a
fullyconnected quantum manybody system. Additionally, we directly measure the
size of the critical energy gap near a quantum phase transition. We expect this
general technique to become an important verification tool for quantum
simulators once experiments advance beyond proofofprinciple demonstrations
and exceed the resources of conventional computers.

We stabilize a chosen radiofrequency beat note between two optical fields
derived from the same modelocked laser pulse train, in order to coherently
manipulate quantum information. This scheme does not require access or active
stabilization of the laser repetition rate. We implement and characterize this
external lock, in the context of twophoton stimulated Raman transitions
between the hyperfine ground states of trapped 171Yb+ quantum bits.

LiteBIRD is a nextgeneration satellite mission to measure the polarization
of the cosmic microwave background (CMB) radiation. On large angular scales the
Bmode polarization of the CMB carries the imprint of primordial gravitational
waves, and its precise measurement would provide a powerful probe of the epoch
of inflation. The goal of LiteBIRD is to achieve a measurement of the
characterizing tensor to scalar ratio $r$ to an uncertainty of $\delta
r=0.001$. In order to achieve this goal we will employ a kilopixel
superconducting detector array on a cryogenically cooled subKelvin focal plane
with an optical system at a temperature of 4~K. We are currently considering
two detector array options; transition edge sensor (TES) bolometers and
microwave kinetic inductance detectors (MKID). In this paper we give an
overview of LiteBIRD and describe a TESbased polarimeter designed to achieve
the target sensitivity of 2~$\mu$K$\cdot$arcmin over the frequency range 50 to
320~GHz.

The American Physical Society's Division of Particles and Fields initiated a
longterm planning exercise over 201213, with the goal of developing the
community's long term aspirations. The subgroup "Dark Energy and CMB" prepared
a series of papers explaining and highlighting the physics that will be studied
with large galaxy surveys and cosmic microwave background experiments. This
paper summarizes the findings of the other papers, all of which have been
submitted jointly to the arXiv.

We derive a timedependent density functional theory appropriate for
calculating the nearedge Xray absorption spectrum in molecules and condensed
matter. The basic assumption is to increase the space of manybody wave
functions from one Slater determinant to two. The equations of motion derived
from Dirac's variational principle provide an exact solution for the linear
response when the interaction Hamiltonian has only a coreelectron field. The
equations can be solved numerically nearly as easily as the ordinary realtime
timedependent KohnSham equations. We carry out the solution under conditions
that permit comparison with the expected powerlaw behavior. Our extracted
powerlaw exponents are similar to those derived by Nozieres and DeDominicis,
but are not in quantitative agreement. We argue that our calculational method
can be readily generalized to density functionals that take into account the
more general electronelectron interactions that are needed for treating
dynamic effects such as plasmon excitations.

The production of Z0 bosons in the reaction ep > eZ0p*, where p* stands for
a proton or a lowmass nucleon resonance, has been studied in ep collisions at
HERA using the ZEUS detector. The analysis is based on a data sample collected
between 1996 and 2007, amounting to 496 pb1 of integrated luminosity. The Z0
was measured in the hadronic decay mode. The elasticity of the events was
ensured by a cut on eta_max < 3.0, where eta_max is the maximum pseudorapidity
of energy deposits in the calorimeter defined with respect to the proton beam
direction. A signal was observed at the Z0 mass. The cross section of the
reaction ep > eZ0p* was measured to be sigma(ep > eZ0p*) = 0.13 +/ 0.06
(stat.) +/ 0.01 (syst.) pb, in agreement with the Standard Model prediction of
0.16 pb. This is the first measurement of Z0 production in ep collisions.

Frustration, or the competition between interacting components of a network,
is often responsible for the complexity of many body systems, from social and
neural networks to protein folding and magnetism. In quantum magnetic systems,
frustration arises naturally from competing spinspin interactions given by the
geometry of the spin lattice or by the presence of longrange antiferromagnetic
couplings. Frustrated magnetism is a hallmark of poorly understood systems such
as quantum spin liquids, spin glasses and spin ices, whose ground states are
massively degenerate and can carry high degrees of quantum entanglement. The
controlled study of frustrated magnetism in materials is hampered by short
dynamical time scales and the presence of impurities, while numerical modeling
is generally intractable when dealing with dynamics beyond N~30 particles.
Alternatively, a quantum simulator can be exploited to directly engineer
prescribed frustrated interactions between controlled quantum systems, and
several smallscale experiments have moved in this direction. In this article,
we perform a quantum simulation of a longrange antiferromagnetic quantum Ising
model with a transverse field, on a crystal of up to N = 16 trapped Yb+ atoms.
We directly control the amount of frustration by continuously tuning the range
of interaction and directly measure spin correlation functions and their
dynamics through spatiallyresolved spin detection. We find a pronounced
dependence of the magnetic order on the amount of frustration, and extract
signatures of quantum coherence in the resulting phases.

Inclusivejet cross sections have been measured in the reaction ep>e+jet+X
for photon virtuality Q2 < 1 GeV2 and gammap centreofmass energies in the
region 142 < W(gammap) < 293 GeV with the ZEUS detector at HERA using an
integrated luminosity of 300 pb1. Jets were identified using the kT, antikT
or SIScone jet algorithms in the laboratory frame. Singledifferential cross
sections are presented as functions of the jet transverse energy, ETjet, and
pseudorapidity, etajet, for jets with ETjet > 17 GeV and 1 < etajet < 2.5. In
addition, measurements of doubledifferential inclusivejet cross sections are
presented as functions of ETjet in different regions of etajet.
Nexttoleadingorder QCD calculations give a good description of the
measurements, except for jets with low ETjet and high etajet. The influence of
nonperturbative effects not related to hadronisation was studied. Measurements
of the ratios of cross sections using different jet algorithms are also
presented; the measured ratios are well described by calculations including up
to O(alphas2) terms. Values of alphas(Mz) were extracted from the measurements
and the energyscale dependence of the coupling was determined. The value of
alphas(Mz) extracted from the measurements based on the kT jet algorithm is
alphas(Mz) = 0.1206 +0.0023 0.0022 (exp.) +0.0042 0.0035 (th.); the results
from the antikT and SIScone algorithms are compatible with this value and have
a similar precision.

Scaled momentum distributions for the strange hadrons K0s and Lambda/bar
Lambda were measured in deep inelastic ep scattering with the ZEUS detector at
HERA using an integrated luminosity of 330 pb1. The evolution of these
distributions with the photon virtuality, Q2, was studied in the kinematic
region 10<Q2<40000 GeV2 and 0.001<x<0.75, where x is the Bjorken scaling
variable. Clear scaling violations are observed. Predictions based on different
approaches to fragmentation were compared to the measurements.
Leadinglogarithm partonshower Monte Carlo calculations interfaced to the Lund
string fragmentation model describe the data reasonably well in the whole range
measured. Nexttoleadingorder QCD calculations based on fragmentation
functions, FFs, extracted from e+e data alone, fail to describe the
measurements. The calculations based on FFs extracted from a global analysis
including e+e, ep and pp data give an improved description. The measurements
presented in this paper have the potential to further constrain the FFs of
quarks, antiquarks and gluons yielding K0s and Lambda/bar Lambda strange
hadrons.

A search for singletop production, $ep \rightarrow etX$, has been performed
with the ZEUS detector at HERA using data corresponding to an integrated
luminosity of $0.37\fbi$. No evidence for top production was found, consistent
with the expectation from the Standard Model. Limits were computed for
singletop production via flavour changing neutral current transitions. The
result was combined with a previous ZEUS result yielding a total luminosity of
0.50fb1. A 95% credibility level upper limit of 0.13 pb was obtained for the
cross section at the centreofmass energy of $\sqrt{s}=315\gev$.