• ### The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package(1801.02634)

Jan. 16, 2018 astro-ph.IM
The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed 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 inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.
• ### The DECam Plane Survey: Optical photometry of two billion objects in the southern Galactic plane(1710.01309)

Oct. 3, 2017 astro-ph.GA
The DECam Plane Survey is a five-band optical and near-infrared survey of the southern Galactic plane with the Dark Energy Camera at Cerro Tololo. The survey is designed to reach past the main-sequence turn-off at the distance of the Galactic center through a reddening E(B-V) of 1.5 mag. Typical single-exposure 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.
• ### Design of the LBNF Beamline Target Station(1612.07293)

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 sub-systems and shielding infrastructure. Protons hit the target producing mesons which are then focused by the horns into a helium-filled 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 physics-driven component optimization of the Target Station.
• ### Analogy between the "Hidden Order" and the Orbital Antiferromagnetism in URu$_{2-x}$Fe$_x$Si$_2$(1608.01748)

Nov. 28, 2016 cond-mat.str-el
We study URu$_{2-x}$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 pseudovector-like $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 degrees-of-freedom of the uranium-5$f$ electrons.
• ### Deciding Graph non-Hamiltonicity via a Closure Algorithm(1611.01710)

Nov. 5, 2016 cs.DS, cs.CC
We present a matching and LP based heuristic algorithm that decides graph non-Hamiltonicity. 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,n-1$, 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 stochastic-like $\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 non-Hamiltonian $G$, integer solutions cannot exist. We recognize non-Hamiltonicity 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 non-Hamiltonian. Over 100 non-Hamiltonian graphs (10 through 104 vertices) and 2000 randomized 31 vertex non-Hamiltonian graphs are tested and correctly decided non-Hamiltonian. 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.
• ### Observation of a Discrete Time Crystal(1609.08684)

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, non-equilibrium Floquet systems subject to a periodic drive can exhibit persistent time-correlations at an emergent sub-harmonic 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 many-body localization (MBL) conditions, and observe a sub-harmonic temporal response that is robust to external perturbations. Such a time crystal opens the door for studying systems with long-range spatial-temporal correlations and novel phases of matter that emerge under intrinsically non-equilibrium conditions.
• ### Power Spectral Density of Magnetization Dynamics Driven by a Jump-Noise Process(1608.04733)

Aug. 16, 2016 cond-mat.mes-hall
Random magnetization dynamics driven by a jump-noise 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.
• ### Collective excitations of dynamic Fermi surface deformations in BaFe$_2$(As$_{0.5}$P$_{0.5}$)$_2$(1607.06575)

We use electronic Raman scattering to study the low-energy 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 non-Fermi liquid to Femi liquid crossover determined by transport, suggesting that the Pomeranchuk instability plays an important role in the low-energy physics of the Fe-based 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}$ protons-on-target 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 kton-equivalent exposure of $2.74 \times 10^{20}$ protons-on-target 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.38-0.65, both at the 68% confidence level, with two statistically-degenerate best fit points at $\sin^2\theta_{23} =$ 0.43 and 0.60. Results for the inverted mass hierarchy are also presented.
• ### Development of readout electronics for POLARBEAR-2 Cosmic Microwave Background experiment(1512.07663)

Dec. 23, 2015 physics.ins-det, astro-ph.IM
The readout of transition-edge sensor (TES) bolometers with a large multiplexing factor is key for the next generation Cosmic Microwave Background (CMB) experiment, Polarbear-2, 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 non-linear 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).
• ### The POLARBEAR-2 and the Simons Array Experiment(1512.07299)

Dec. 22, 2015 astro-ph.CO, astro-ph.IM
We present an overview of the design and status of the \Pb-2 and the Simons Array experiments. \Pb-2 is a Cosmic Microwave Background polarimetry experiment which aims to characterize the arc-minute angular scale B-mode signal from weak gravitational lensing and search for the degree angular scale B-mode signal from inflationary gravitational waves. The receiver has a 365~mm diameter focal plane cooled to 270~milli-Kelvin. The focal plane is filled with 7,588 dichroic lenslet-antenna coupled polarization sensitive Transition Edge Sensor (TES) bolometric pixels that are sensitive to 95~GHz and 150~GHz bands simultaneously. The TES bolometers are read-out 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. \Pb-2 will deploy in 2016 in the Atacama desert in Chile. The Simons Array is a project to further increase sensitivity by deploying three \Pb-2 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 tensor-to-scalar 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 MOLLER Experiment: An Ultra-Precise Measurement of the Weak Mixing Angle Using M{\o}ller Scattering(1411.4088)

Dec. 3, 2014 hep-ph, nucl-ex
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 parity-violating asymmetry $A_{PV}$ in polarized electron-electron (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 CP-conserving process over the next decade, and yields a unique window to new physics at MeV and multi-TeV 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 parity-violating electron scattering experiments
• ### Experimental Realization of a Quantum Integer-Spin Chain with Controllable Interactions(1410.0937)

The physics of interacting integer-spin chains has been a topic of intense theoretical interest, particularly in the context of symmetry-protected topological phases. However, there has not been a controllable model system to study this physics experimentally. We demonstrate how spin-dependent forces on trapped ions can be used to engineer an effective system of interacting spin-1 particles. Our system evolves coherently under an applied spin-1 XY Hamiltonian with tunable, long-range 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 three-level 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 symmetry-protected order in spin-1 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 single-differential cross-sections d(sigma)/dQ2, d(sigma)/dx and d(sigma)/dy and the reduced cross-section were measured in the kinematic region Q2 > 185 GeV2 and y < 0.9, where Q2 is the four-momentum 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 pb-1 collected with the ZEUS detector in 2006 and 2007 at a centre-of-mass 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 e-p 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.
• ### Coherent Imaging Spectroscopy of a Quantum Many-Body Spin System(1401.5751)

Jan. 22, 2014 quant-ph
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 fully-connected quantum many-body 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 proof-of-principle demonstrations and exceed the resources of conventional computers.
• ### Beat note stabilization of mode-locked lasers for quantum information processing(1312.6846)

We stabilize a chosen radiofrequency beat note between two optical fields derived from the same mode-locked 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 two-photon stimulated Raman transitions between the hyperfine ground states of trapped 171-Yb+ quantum bits.
• ### Mission design of LiteBIRD(1311.2847)

Nov. 12, 2013 astro-ph.CO, astro-ph.IM
LiteBIRD is a next-generation satellite mission to measure the polarization of the cosmic microwave background (CMB) radiation. On large angular scales the B-mode 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 kilo-pixel superconducting detector array on a cryogenically cooled sub-Kelvin 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 TES-based polarimeter designed to achieve the target sensitivity of 2~$\mu$K$\cdot$arcmin over the frequency range 50 to 320~GHz.
• ### Dark Energy and CMB(1309.5386)

Sept. 24, 2013 astro-ph.CO
The American Physical Society's Division of Particles and Fields initiated a long-term planning exercise over 2012-13, with the goal of developing the community's long term aspirations. The sub-group "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.
• ### Time-dependent density functional theory for X-ray near-edge spectroscopy(1306.6571)

We derive a time-dependent density functional theory appropriate for calculating the near-edge X-ray absorption spectrum in molecules and condensed matter. The basic assumption is to increase the space of many-body 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 core-electron field. The equations can be solved numerically nearly as easily as the ordinary real-time time-dependent Kohn-Sham equations. We carry out the solution under conditions that permit comparison with the expected power-law behavior. Our extracted power-law 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 electron-electron 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 low-mass 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 pb-1 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.
• ### Emergence and Frustration of Magnetic Order with Variable-Range Interactions in a Trapped Ion Quantum Simulator(1210.0142)

Sept. 29, 2012 quant-ph, cond-mat.mtrl-sci
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 spin-spin interactions given by the geometry of the spin lattice or by the presence of long-range 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 small-scale experiments have moved in this direction. In this article, we perform a quantum simulation of a long-range 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 spatially-resolved 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.
• Inclusive-jet cross sections have been measured in the reaction ep->e+jet+X for photon virtuality Q2 < 1 GeV2 and gamma-p centre-of-mass energies in the region 142 < W(gamma-p) < 293 GeV with the ZEUS detector at HERA using an integrated luminosity of 300 pb-1. Jets were identified using the kT, anti-kT or SIScone jet algorithms in the laboratory frame. Single-differential 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 double-differential inclusive-jet cross sections are presented as functions of ETjet in different regions of etajet. Next-to-leading-order QCD calculations give a good description of the measurements, except for jets with low ETjet and high etajet. The influence of non-perturbative 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 energy-scale 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 anti-kT 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 pb-1. 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. Leading-logarithm parton-shower Monte Carlo calculations interfaced to the Lund string fragmentation model describe the data reasonably well in the whole range measured. Next-to-leading-order 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, anti-quarks and gluons yielding K0s and Lambda/bar Lambda strange hadrons.
• A search for single-top 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 single-top production via flavour changing neutral current transitions. The result was combined with a previous ZEUS result yielding a total luminosity of 0.50fb-1. A 95% credibility level upper limit of 0.13 pb was obtained for the cross section at the centre-of-mass energy of $\sqrt{s}=315\gev$.