Frequency domain multiplexing (FDM) is the baseline readout system for the X-ray Integral Field Unit (X-IFU) on board the Athena mission. Under the FDM scheme, TESs are coupled to a passive LC filter and biased with alternating current (AC bias) at MHz frequencies. Using high-quality factor LC filters and room temperature electronics developed at SRON and low-noise two-stage SQUID amplifiers provided by VTT, we have recently demonstrated good performance with the FDM readout of Mo/Au TES calorimeters with Au/Bi absorbers. We have achieved a performance requested for the demonstration model (DM) with the single pixel AC bias ($\Delta E=$1.8 eV) and 9 pixel multiplexing ($\Delta E=$2.6 eV) modes. We have also demonstrated 14-pixel multiplexing with an average energy resolution of 3.3 eV, which is limited by non-fundamental issues related to FDM readout in our lab setup.

The Micro-X sounding rocket flew for the first time on July 22, 2018, becoming the first program to fly Transition-Edge Sensors and multiplexing SQUID readout electronics in space. While a rocket pointing failure led to no time on-target, the success of the flight systems was demonstrated. The successful flight operation of the instrument puts the program in a position to modify the payload for indirect galactic dark matter searches. The payload modifications are motivated by the science requirements of this observation. Micro-X can achieve world-leading sensitivity in the keV regime with a single flight. Dark matter sensitivity projections have been updated to include recent observations and the expected sensitivity of Micro-X to these observed fluxes. If a signal is seen (as seen in the X-ray satellites), Micro-X can differentiate an atomic line from a dark matter signature.

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

We describe an on-sky demonstration of a microwave-multiplexing readout system in one of the receivers of the Keck Array, a polarimetry experiment observing the cosmic microwave background at the South Pole. During the austral summer of 2018-2019, we replaced the time-division multiplexing readout system with microwave-multiplexing components including superconducting microwave resonators coupled to radio-frequency superconducting quantum interference devices at the sub-Kelvin focal plane, coaxial-cable plumbing and amplification between room temperature and the cold stages, and a SLAC Microresonator Radio Frequency system for the warm electronics. In the range 5-6 GHz, a single coaxial cable reads out 528 channels. The readout system is coupled to transition-edge sensors, which are in turn coupled to 150-GHz slot-dipole phased-array antennas. Observations began in April 2019, and we report here on an initial characterization of the system performance.

Micro-X is a sounding rocket-borne instrument that uses a microcalorimeter array to perform high-resolution X-ray spectroscopy. This instrument flew for the first time on July 22nd, 2018 from the White Sands Missile Range, USA. This flight marks the first successful operation of a Transition-Edge Sensor array and its time division multiplexing read-out system in space. This launch was dedicated to the observation of the supernova remnant Cassiopeia A. A failure in the attitude control system prevented the rocket from pointing and led to no time on target. The on-board calibration source provided X-rays in flight, and it is used to compare detector performance during pre-flight integration, flight, and after the successful post-flight recovery. This calibration data demonstrates the capabilities of the detector in a space environment as well as its potential for future flights.

We report on the development of multi-absorber transition edge sensors (TESs), referred to as hydras. A hydra consists of multiple x-ray absorbers each with a different thermal conductance to a TES. Position information is encoded in the pulse shape. With some trade-off in performance, hydras enable very large format arrays without the prohibitive increase in bias and read-out components associated with arrays of individual TESs. Hydras are under development for the next generation of space telescope such as Lynx. Lynx is a NASA concept under study that will combine a < 1 arcsecond angular resolution optic with 100,000-pixel microcalorimeter array with energy resolution of deltaE_FWHM ~ 3 eV in the soft x-ray energy range. We present first results from hydras with 25-pixels for Lynx. Designs with absorbers on a 25 micron and 50 micron pitch are studied. Arrays incorporate, for the first time, microstrip buried wiring layers of suitable pitch and density required to readout a full-scale Lynx array. The resolution from the coadded energy histogram including all 25-pixels was deltaE_FWHM = 1.66+/-0.02 eV and 3.34+/-0.06 eV at an energy of 1.5 keV for the 25 micron and 50 micron absorber designs respectively. Position discrimination is demonstrated from parameterization of the rise-time.

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

For the first time in human history, technologies have matured sufficiently to enable a mission capable of discovering and characterizing habitable planets like Earth orbiting sunlike stars other than the Sun. At the same time, such a platform would enable unique science not possible from ground-based facilities. This science is broad and exciting, ranging from new investigations of our own solar system to a full range of astrophysics disciplines. The Habitable Exoplanet Observatory, or HabEx, is one of four studies currently being undertaken by NASA in preparation for the 2020 Astrophysics Decadal Survey. HabEx has been designed to be the Great Observatory of the 2030s, with community involvement through a competed and funded Guest Observer (GO) program. This interim report describes the HabEx baseline concept, which is a space-based 4-meter diameter telescope mission concept with ultraviolet (UV), optical, and near-infrared (near-IR) imaging and spectroscopy capabilities. More information on HabEx can be found at https://www.jpl.nasa.gov/habex

The next generation of cryogenic CMB and submillimeter cameras under development require densely instrumented sensor arrays to meet their science goals. The readout of large numbers ($\sim$10,000--100,000 per camera) of sub-Kelvin sensors, for instance as proposed for the CMB-S4 experiment, will require substantial improvements in cold and warm readout techniques. To reduce the readout cost per sensor and integration complexity, efforts are presently focused on achieving higher multiplexing density while maintaining readout noise subdominant to intrinsic detector noise. Highly-multiplexed cold readout technologies in active development include Microwave Kinetic Inductance Sensors (MKIDs) and microwave rf-SQUIDs. Both exploit the high quality factors of superconducting microwave resonators to densely channelize sub-Kelvin sensors into the bandwidth of a microwave transmission line. We present advancements in the development of a new warm readout system for microwave SQUID multiplexing, the SLAC Superconducting Microresonator RF electronics, or SMuRF. The SMuRF system is unique in its ability to track each tone, minimizing the total RF power required to read out each resonator, thereby significantly reducing the linearity requirements on the cold and warm readout. Here, we present measurements of the readout noise and linearity of the first full SMuRF system, including a demonstration of closed-loop tone tracking on a 528 channel cryogenic microwave SQUID multiplexer. SMuRF is being explored as a potential readout solution for several future CMB projects including Simons Observatory, BICEP Array, CCAT-prime, Ali-CPT, and CMB-S4. Parallel development of the platform is underway to adapt SMuRF to read out both MKID and fast X-ray TES calorimeter arrays.

The X-ray Integral Field Unit (X-IFU) is the high resolution X-ray spectrometer of the ESA Athena X-ray observatory. Over a field of view of 5' equivalent diameter, it will deliver X-ray spectra from 0.2 to 12 keV with a spectral resolution of 2.5 eV up to 7 keV on ~5 arcsecond pixels. The X-IFU is based on a large format array of super-conducting molybdenum-gold Transition Edge Sensors cooled at about 90 mK, each coupled with an absorber made of gold and bismuth with a pitch of 249 microns. A cryogenic anti-coincidence detector located underneath the prime TES array enables the non X-ray background to be reduced. A bath temperature of about 50 mK is obtained by a series of mechanical coolers combining 15K Pulse Tubes, 4K and 2K Joule-Thomson coolers which pre-cool a sub Kelvin cooler made of a 3He sorption cooler coupled with an Adiabatic Demagnetization Refrigerator. Frequency domain multiplexing enables to read out 40 pixels in one single channel. A photon interacting with an absorber leads to a current pulse, amplified by the readout electronics and whose shape is reconstructed on board to recover its energy with high accuracy. The defocusing capability offered by the Athena movable mirror assembly enables the X-IFU to observe the brightest X-ray sources of the sky (up to Crab-like intensities) by spreading the telescope point spread function over hundreds of pixels. Thus the X-IFU delivers low pile-up, high throughput (>50%), and typically 10 eV spectral resolution at 1 Crab intensities, i.e. a factor of 10 or more better than Silicon based X-ray detectors. In this paper, the current X-IFU baseline is presented, together with an assessment of its anticipated performance in terms of spectral resolution, background, and count rate capability. The X-IFU baseline configuration will be subject to a preliminary requirement review that is scheduled at the end of 2018.

With its array of 3840 Transition Edge Sensors (TESs) operated at 90 mK, the X-Ray Integral Field Unit (X-IFU) on board the ESA L2 mission Athena will provide spatially resolved high-resolution spectroscopy (2.5 eV FWHM up to 7 keV) over the 0.2 to 12 keV bandpass. The in-flight performance of the X-IFU will be strongly affected by the calibration of the instrument. Uncertainties in the knowledge of the overall system, from the filter transmission to the energy scale, may introduce systematic errors in the data, which could potentially compromise science objectives - notably those involving line characterisation e.g. turbulence velocity measurements - if not properly accounted for. Defining and validating calibration requirements is therefore of paramount importance. In this paper, we put forward a simulation tool based on the most up-to-date configurations of the various subsystems (e.g. filters, detector absorbers) which allows us to estimate systematic errors related to uncertainties in the instrumental response. Notably, the effect of uncertainties in the energy resolution and of the instrumental quantum efficiency on X-IFU observations is assessed, by taking as a test case the measurements of the iron K complex in the hot gas surrounding clusters of galaxies. In-flight and ground calibration of the energy resolution and the quantum efficiency is also addressed. We demonstrate that provided an accurate calibration of the instrument, such effects should be low in both cases with respect to statistics during observations.

Large arrays of cryogenic sensors for various imaging applications ranging across x-ray, gamma-ray, Cosmic Microwave Background (CMB), mm/sub-mm, as well as particle detection increasingly rely on superconducting microresonators for high multiplexing factors. These microresonators take the form of microwave SQUIDs that couple to Transition-Edge Sensors (TES) or Microwave Kinetic Inductance Detectors (MKIDs). In principle, such arrays can be read out with vastly scalable software-defined radio using suitable FPGAs, ADCs and DACs. In this work, we share plans and show initial results for SLAC Microresonator Radio Frequency (SMuRF) electronics, a next-generation control and readout system for superconducting microresonators. SMuRF electronics are unique in their implementation of specialized algorithms for closed-loop tone tracking, which consists of fast feedback and feedforward to each resonator's excitation parameters based on transmission measurements. Closed-loop tone tracking enables improved system linearity, a significant increase in sensor count per readout line, and the possibility of overcoupled resonator designs for enhanced dynamic range. Low-bandwidth prototype electronics were used to demonstrate closed-loop tone tracking on twelve 300-kHz-wide microwave SQUID resonators, spaced at $\sim$6 MHz with center frequencies $\sim$5-6 GHz. We achieve multi-kHz tracking bandwidth and demonstrate that the noise floor of the electronics is subdominant to the noise intrinsic in the multiplexer.

New radio (MeerKAT and Parkes) and X-ray (XMM-Newton, Swift, Chandra, and NuSTAR) observations of PSR J1622-4950 indicate that the magnetar, in a quiescent state since at least early 2015, reactivated between 2017 March 19 and April 5. The radio flux density, while variable, is approximately 100x larger than during its dormant state. The X-ray flux one month after reactivation was at least 800x larger than during quiescence, and has been decaying exponentially on a 111+/-19 day timescale. This high-flux state, together with a radio-derived rotational ephemeris, enabled for the first time the detection of X-ray pulsations for this magnetar. At 5%, the 0.3-6 keV pulsed fraction is comparable to the smallest observed for magnetars. The overall pulsar geometry inferred from polarized radio emission appears to be broadly consistent with that determined 6-8 years earlier. However, rotating vector model fits suggest that we are now seeing radio emission from a different location in the magnetosphere than previously. This indicates a novel way in which radio emission from magnetars can differ from that of ordinary pulsars. The torque on the neutron star is varying rapidly and unsteadily, as is common for magnetars following outburst, having changed by a factor of 7 within six months of reactivation.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

In order to constrain the evolutionary state of the red supergiant Betelgeuse, we have produced a suite of models with ZAMS masses from 15 to 25 Msun in intervals of 1 Msun including the effects of rotation. The models were computed with the stellar evolutionary code MESA. For non-rotating models we find results that are similar to other work. It is somewhat difficult to find models that agree within 1 sigma of the observed values of R, Teff and L, but modestly easy within 3 sigma uncertainty. Incorporating the nominal observed rotational velocity, ~15 km/s, yields significantly different, and challenging, constraints. This velocity constraint is only matched when the models first approach the base of the red supergiant branch (RSB), having crossed the Hertzsprung gap, but not yet having ascended the RSB and most violate even generous error bars on R, Teff and L. Models at the tip of the RSB typically rotate at only ~0.1 km/s, independent of any reasonable choice of initial rotation. We discuss the possible uncertainties in our modeling and the observations, including the distance to Betelgeuse, the rotation velocity, and model parameters. We summarize various options to account for the rotational velocity and suggest that one possibility is that Betelgeuse merged with a companion star of about 1 Msun as it ascended the RSB, in the process producing the ring structure observed at about 7' away. A past coalescence would complicate attempts to understand the evolutionary history and future of Betelgeuse.

We are developing the frequency domain multiplexing (FDM) read-out of transition-edge sensor (TES) microcalorimeters for the X-ray Integral Field Unit (X-IFU) instrument on board of the future European X-Ray observatory Athena. The X-IFU instrument consists of an array of $\sim$3840 TESs with a high quantum efficiency ($>$90 \%) and spectral resolution $\Delta E$=2.5 eV $@$ 7 keV ($E/\Delta E\sim$2800). FDM is currently the baseline readout system for the X-IFU instrument. Using high quality factor LC filters and room temperature electronics developed at SRON and low-noise two stage SQUID amplifiers provided by VTT, we have recently demonstrated good performance with the FDM readout of Mo/Au TES calorimeters with Au/Bi absorbers. An integrated noise equivalent power resolution of about 2.0 eV at 1.7 MHz has been demonstrated with a pixel from a new TES array from NASA/Goddard (GSFC-A2). We have achieved X-ray energy resolutions $\sim$2.5 eV at AC bias frequency at 1.7 MHz in the single pixel read-out. We have also demonstrated for the first time an X-ray energy resolution around 3.0 eV in a 6 pixel FDM read-out with TES array (GSFC-A1). In this paper we report on the single pixel performance of these microcalorimeters under MHz AC bias, and further results of the performance of these pixels under FDM.

A comprehensive study on the atmospheric neutrino flux in the energy region from sub-GeV up to several TeV using the Super-Kamiokande water Cherenkov detector is presented in this paper. The energy and azimuthal spectra of the atmospheric ${\nu}_e+{\bar{\nu}}_e$ and ${\nu}_{\mu}+{\bar{\nu}}_{\mu}$ fluxes are measured. The energy spectra are obtained using an iterative unfolding method by combining various event topologies with differing energy responses. The azimuthal spectra depending on energy and zenith angle, and their modulation by geomagnetic effects, are also studied. A predicted east-west asymmetry is observed in both the ${\nu}_e$ and ${\nu}_{\mu}$ samples at 8.0 {\sigma} and 6.0 {\sigma} significance, respectively, and an indication that the asymmetry dipole angle changes depending on the zenith angle was seen at the 2.2 {\sigma} level. The measured energy and azimuthal spectra are consistent with the current flux models within the estimated systematic uncertainties. A study of the long-term correlation between the atmospheric neutrino flux and the solar magnetic activity cycle is also performed, and a weak indication of a correlation was seen at the 1.1 {\sigma} level, using SK I-IV data spanning a 20 year period. For particularly strong solar activity periods known as Forbush decreases, no theoretical prediction is available, but a deviation below the typical neutrino event rate is seen at the 2.4 {\sigma} level.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

The X-ray Integral Field Unit (X-IFU) on board the Advanced Telescope for High-ENergy Astrophysics (Athena) will provide spatially resolved high-resolution X-ray spectroscopy from 0.2 to 12 keV, with 5 arc second pixels over a field of view of 5 arc minute equivalent diameter and a spectral resolution of 2.5 eV up to 7 keV. In this paper, we first review the core scientific objectives of Athena, driving the main performance parameters of the X-IFU, namely the spectral resolution, the field of view, the effective area, the count rate capabilities, the instrumental background. We also illustrate the breakthrough potential of the X-IFU for some observatory science goals. Then we briefly describe the X-IFU design as defined at the time of the mission consolidation review concluded in May 2016, and report on its predicted performance. Finally, we discuss some options to improve the instrument performance while not increasing its complexity and resource demands (e.g. count rate capability, spectral resolution). The X-IFU will be provided by an international consortium led by France, The Netherlands and Italy, with further ESA member state contributions from Belgium, Finland, Germany, Poland, Spain, Switzerland and two international partners from the United States and Japan.

We present a real-time supernova neutrino burst monitor at Super-Kamiokande (SK). Detecting supernova explosions by neutrinos in real time is crucial for giving a clear picture of the explosion mechanism. Since the neutrinos are expected to come earlier than light, a fast broadcasting of the detection may give astronomers a chance to make electromagnetic radiation observations of the explosions right at the onset. The role of the monitor includes a fast announcement of the neutrino burst detection to the world and a determination of the supernova direction. We present the online neutrino burst detection system and studies of the direction determination accuracy based on simulations at SK.

We determine the neutrino spectra arising from low-mass (4-10 GeV) dark matter annihilating in the sun. We also determine the low-mass dark matter capture rates (element by element in the sun), assuming dark matter interacts either through elastic contact interactions, elastic long-range interactions, or inelastic contact interactions. These are the non-detector-specific data needed for determining the sensitivity of a neutrino detector to dark matter annihilating in the sun. As an application, we estimate the sensitivity of a one kiloton liquid scintillation neutrino detector (such as KamLAND) and LBNE (LAr-based) to low-mass dark matter with long-range interactions and compare this to the expected CDMS sensitivity. It is found that KamLAND's sensitivity can exceed that obtainable from the current CDMS data set by up to two orders of magnitude.

The "Airships: A New Horizon for Science" study at the Keck Institute for Space Studies investigated the potential of a variety of airships currently operable or under development to serve as observatories and science instrumentation platforms for a range of space, atmospheric, and Earth science. The participants represent a diverse cross-section of the aerospace sector, NASA, and academia. Over the last two decades, there has been wide interest in developing a high altitude, stratospheric lighter-than-air (LTA) airship that could maneuver and remain in a desired geographic position (i.e., "station-keeping") for weeks, months or even years. Our study found considerable scientific value in both low altitude (< 40 kft) and high altitude (> 60 kft) airships across a wide spectrum of space, atmospheric, and Earth science programs. Over the course of the study period, we identified stratospheric tethered aerostats as a viable alternative to airships where station-keeping was valued over maneuverability. By opening up the sky and Earth's stratospheric horizon in affordable ways with long-term flexibility, airships allow us to push technology and science forward in a project-rich environment that complements existing space observatories as well as aircraft and high-altitude balloon missions.

Procedures and results on hardware level detector calibration in Super-Kamiokande (SK) are presented in this paper. In particular, we report improvements made in our calibration methods for the experimental phase IV in which new readout electronics have been operating since 2008. The topics are separated into two parts. The first part describes the determination of constants needed to interpret the digitized output of our electronics so that we can obtain physical numbers such as photon counts and their arrival times for each photomultiplier tube (PMT). In this context, we developed an in-situ procedure to determine high-voltage settings for PMTs in large detectors like SK, as well as a new method for measuring PMT quantum efficiency and gain in such a detector. The second part describes the modeling of the detector in our Monte Carlo simulation, including in particular the optical properties of its water target and their variability over time. Detailed studies on the water quality are also presented. As a result of this work, we achieved a precision sufficient for physics analysis over a wide energy range (from a few MeV to above a TeV). For example, the charge determination was understood at the 1% level, and the timing resolution was 2.1 nsec at the one-photoelectron charge level and 0.5 nsec at the 100-photoelectron charge level.

We present the result of an indirect search for high energy neutrinos from WIMP annihilation in the Sun using upward-going muon (upmu) events at Super-Kamiokande. Datasets from SKI-SKIII (3109.6 days) were used for the analysis. We looked for an excess of neutrino signal from the Sun as compared with the expected atmospheric neutrino background in three upmu categories: stopping, non-showering, and showering. No significant excess was observed. The 90% C.L. upper limits of upward-going muon flux induced by WIMPs of 100 GeV/c$^2$ were 6.4$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ and 4.0$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ for the soft and hard annihilation channels, respectively. These limits correspond to upper limits of 4.5$\times10^{-39}$ cm$^{-2}$ and 2.7$\times10^{-40}$ cm$^{-2}$ for spin-dependent WIMP-nucleon scattering cross sections in the soft and hard annihilation channels, respectively.

We present a modern ephemeris and modern light curve of the first-discovered Cepheid variable in M31, Edwin Hubble's M31-V1. Observers of the American Association of Variable Star Observers undertook these observations during the latter half of 2010. The observations were in support of an outreach program by the Space Telescope Science Institute's Hubble Heritage project, but the resulting data are the first concentrated observations of M31-V1 made in modern times. AAVSO observers obtained 214 V-band, Rc-band, and unfiltered observations from which a current ephemeris was derived. The ephemeris derived from these observations is JD(Max) = 2455430.5(+/-0.5) + 31.4 (+/-0.1) E. The period derived from the 2010 data are in agreement with the historic values of the period, but the single season of data precludes a more precise determination of the period or measurement of the period change using these data alone. However, using an ephemeris based upon the period derived by Baade and Swope we are able to fit all of the observed data acceptably well. Continued observations in the modern era will be very valuable in linking these modern data with data from the 1920s-30s and 1950s, and will enable us to measure period change in this historic Cepheid. In particular, we strongly encourage intensive observations of this star around predicted times of maximum to constrain the date of maximum to better than 0.5 days.

We have recently shown that normal-metal/superconductor (N/S) bilayer TESs (superconducting Transition-Edge Sensors) exhibit weak-link behavior.1 Here we extend our understanding to include TESs with added noise-mitigating normal-metal structures (N structures). We find TESs with added Au structures also exhibit weak-link behavior as evidenced by exponential temperature dependence of the critical current and Josephson-like oscillations of the critical current with applied magnetic field. We explain our results in terms of an effect converse to the longitudinal proximity effect (LoPE)1, the lateral inverse proximity effect (LaiPE), for which the order parameter in the N/S bilayer is reduced due to the neighboring N structures. Resistance and critical current measurements are presented as a function of temperature and magnetic field taken on square Mo/Au bilayer TESs with lengths ranging from 8 to 130 {\mu}m with and without added N structures. We observe the inverse proximity effect on the bilayer over in-plane distances many tens of microns and find the transition shifts to lower temperatures scale approximately as the inverse square of the in- plane N-structure separation distance, without appreciable broadening of the transition width. We also present evidence for nonequilbrium superconductivity and estimate a quasiparticle lifetime of 1.8 \times 10-10 s for the bilayer. The LoPE model is also used to explain the increased conductivity at temperatures above the bilayer's steep resistive transition.

We consider the prospects for liquid scintillation experiments (with a focus on KamLAND) to detect the flux of electron neutrinos arising from dark matter annihilation in the core of the sun. We show that, with data already taken, KamLAND can provide the greatest sensitivity to the dark matter-proton spin-dependent scattering cross-section for dark matter lighter than 20 GeV. It is also possible to probe the dark matter-nucleon spin-independent scattering cross-section for isospin-violating dark matter lighter than 10 GeV. KamLAND can thus potentially confirm the dark matter interpretation of the DAMA and CoGeNT signals, utilizing data already taken.

We report on TeV gamma-ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV gamma-ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hrs of VERITAS and 96 hrs of Whipple 10m data were acquired between January 2006 and June 2008. We present the results of a study of the TeV gamma-ray energy spectra as a function of time, and for different flux levels. On May 2nd and 3rd, 2008, bright TeV gamma-ray flares were detected with fluxes reaching the level of 10 Crab. The TeV gamma-ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio waveband. The combination of the RXTE and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV gamma-rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone SSC model.

We present the science case for a broadband X-ray imager with high-resolution spectroscopy, including simulations of X-ray spectral diagnostics of both active regions and solar flares. This is part of a trilogy of white papers discussing science, instrument (Bandler et al. 2010), and missions (Bookbinder et al. 2010) to exploit major advances recently made in transition-edge sensor (TES) detector technology that enable resolution better than 2 eV in an array that can handle high count rates. Combined with a modest X-ray mirror, this instrument would combine arcsecondscale imaging with high-resolution spectra over a field of view sufficiently large for the study of active regions and flares, enabling a wide range of studies such as the detection of microheating in active regions, ion-resolved velocity flows, and the presence of non-thermal electrons in hot plasmas. It would also enable more direct comparisons between solar and stellar soft X-ray spectra, a waveband in which (unusually) we currently have much better stellar data than we do of the Sun.

We have found experimentally that the critical current of a square superconducting transition-edge sensor (TES) depends exponentially upon the side length L and the square root of the temperature T. As a consequence, the effective transition temperature Tc of the TES is current-dependent and at fixed current scales as 1/L^2. We also have found that the critical current can show clear Fraunhofer-like oscillations in an applied magnetic field, similar to those found in Josephson junctions. The observed behavior has a natural theoretical explanation in terms of longitudinal proximity effects if the TES is regarded as a weak link between superconducting leads. We have observed the proximity effect in these devices over extraordinarily long lengths exceeding 100 um.

The Advanced Technology Large-Aperture Space Telescope (ATLAST) is a set of mission concepts for the next generation of UVOIR space observatory with a primary aperture diameter in the 8-m to 16-m range that will allow us to perform some of the most challenging observations to answer some of our most compelling questions, including "Is there life elsewhere in the Galaxy?" We have identified two different telescope architectures, but with similar optical designs, that span the range in viable technologies. The architectures are a telescope with a monolithic primary mirror and two variations of a telescope with a large segmented primary mirror. This approach provides us with several pathways to realizing the mission, which will be narrowed to one as our technology development progresses. The concepts invoke heritage from HST and JWST design, but also take significant departures from these designs to minimize complexity, mass, or both. Our report provides details on the mission concepts, shows the extraordinary scientific progress they would enable, and describes the most important technology development items. These are the mirrors, the detectors, and the high-contrast imaging technologies, whether internal to the observatory, or using an external occulter. Experience with JWST has shown that determined competitors, motivated by the development contracts and flight opportunities of the new observatory, are capable of achieving huge advances in technical and operational performance while keeping construction costs on the same scale as prior great observatories.

Since September 2005, the Whipple 10m Gamma-ray Telescope has been operated primarily as a blazar monitor. The five Northern Hemisphere blazars that have already been detected at the Whipple Observatory, Markarian 421, H1426+428, Markarian 501, 1ES 1959+650 and 1ES 2344+514, are monitored routinely each night that they are visible. We report on the Markarian 421 observations taken from November 2005 to June 2006 in the gamma-ray, X-ray, optical and radio bands. During this time, Markarian 421 was found to be variable at all wavelengths probed. Both the variability and the correlations among different energy regimes are studied in detail here. A tentative correlation, with large spread, was measured between the X-ray and gamma-ray bands, while no clear correlation was evident among the other energy bands. In addition to this, the well-sampled spectral energy distribution of Markarian 421 (1101+384) is presented for three different activity levels. The observations of the other blazar targets will be reported separately.

The effect of weak gravitational lensing on the cosmic microwave background (CMB) temperature anisotropies and polarization will provide access to cosmological information that cannot be obtained from the primary anisotropies alone. We compare the information content of the lensed B-mode polarization power spectrum, properly accounting for the non-Gaussian correlations between the power on different scales, with that of the unlensed CMB fields and the lensing potential. The latter represent the products of an (idealised) optimal analysis that exploits the lens-induced non-Gaussianity to reconstruct the fields. Compressing the non-Gaussian lensed CMB into power spectra is wasteful and leaves a tight degeneracy between the equation of state of dark energy and neutrino mass that is much stronger than in the more optimal analysis. Despite this, a power spectrum analysis will be a useful first step in analysing future B-mode polarization data. For this reason, we also consider how to extract accurate parameter constraints from the lensed B-mode power spectrum. We show with simulations that for cosmic-variance-limited measurements of the lensed B-mode power, including the non-Gaussian correlations in existing likelihood approximations gives biased parameter results. We develop a more refined likelihood approximation that performs significantly better. This new approximation should also be of more general interest in the wider context of parameter estimation from Gaussian CMB data.

A simple method is presented for the rapid simulation of statistically-isotropic non-Gaussian maps of CMB temperature fluctuations with a given power spectrum and analytically-calculable bispectrum and higher-order polyspectra. The $n$th-order correlators of the pixel values may also be calculated analytically. The cumulants of the simulated map may be used to obtain an expression for the probability density function of the pixel temperatures. The statistical properties of the simulated map are determined by the univariate non-Gaussian distribution from which pixel values are drawn independently in the first stage of the simulation process. We illustrate the method using a non-Gaussian distribution derived from the wavefunctions of the harmonic oscillator. The basic simulation method is easily extended to produce non-Gaussian maps with a given power spectrum and diagonal bispectrum.

We estimate the bispectrum of the Very Small Array data from the compact and extended configuration observations released in December 2002, and compare our results to those obtained from Gaussian simulations. There is a slight excess of large bispectrum values for two individual fields, but this does not appear when the fields are combined. Given our expected level of residual point sources, we do not expect these to be the source of the discrepancy. Using the compact configuration data, we put an upper limit of 5400 on the value of f_NL, the non-linear coupling parameter, at 95 per cent confidence. We test our bispectrum estimator using non-Gaussian simulations with a known bispectrum, and recover the input values.

We have carried out 3-D numerical simulations of the dynamical bar instability in a rotating star and the resulting gravitational radiation using both an Eulerian code written in cylindrical coordinates and a smooth particle hydrodynamics (SPH) code. The star is modeled initially as a polytrope with index $n = 3/2$ and $T_{\rm rot}/|W| \approx 0.30$, where $T_{\rm rot}$ is the rotational kinetic energy and $|W|$ is the gravitational potential energy. In both codes the gravitational field is purely Newtonian, and the gravitational radiation is calculated in the quadrupole approximation. We have run 3 simulations with the Eulerian code, varying the number of angular zones and the treatment of the boundary between the star and the vacuum. Using the SPH code we did 7 runs, varying the number of particles, the artificial viscosity, and the type of initial model. We compare the growth rate and rotation speed of the bar, the mass and angular momentum distributions, and the gravitational radiation quantities. We highlight the successes and difficulties of both methods, and make suggestions for future improvements.

We present the first calculations of the gravitational radiation produced by nonaxisymmetric dynamical instability in a rapidly rotating compact star. The star deforms into a bar shape, shedding $\sim 4\%$ of its mass and $\sim 17\%$ of its angular momentum. The gravitational radiation is calculated in the quadrupole approximation. For a mass $M \sim 1.4$ M$_{\odot}$ and radius $R \sim 10$ km, the gravitational waves have frequency $\sim 4$ kHz and amplitude $h \sim 2 \times 10^{-22}$ at the distance of the Virgo Cluster. They carry off energy $\Delta E/M \sim 0.1\%$ and radiate angular momentum $\Delta J/J \sim 0.7\%$.

We have developed a 3-D Eulerian hydrodynamics code to model sources of gravitational radiation. The code is written in cylindrical coordinates $(\varpi,z,\varphi)$ and has moving grids in the $\varpi$ and $z$-directions. We use Newtonian gravity and calculate the gravitational radiation in the quadrupole approximation. This code has been tested on a variety of problems to verify its accuracy and stability, and the results of these tests are reported here.