We outline the experimental concept and key scientific capabilities of AION (Atom Interferometer Observatory and Network), a proposed UK-based experimental programme using cold strontium atoms to search for ultra-light dark matter, to explore gravitational waves in the mid-frequency range between the peak sensitivities of the LISA and LIGO/Virgo/ KAGRA/INDIGO/Einstein Telescope/Cosmic Explorer experiments, and to probe other frontiers in fundamental physics. AION would complement other planned searches for dark matter, as well as probe mergers involving intermediate mass black holes and explore early universe cosmology. AION would share many technical features with the MAGIS experimental programme in the US, and synergies would flow from operating AION in a network with this experiment, as well as with other atom interferometer experiments such as MIGA, ZAIGA and ELGAR. Operating AION in a network with other gravitational wave detectors such as LIGO, Virgo and LISA would also offer many synergies.

Data from the Gaia satellite show that the solar neighbourhood of the Milky Way's stellar halo is imprinted with substructure from several accretion events. Evidence of these events is found in "the Shards", stars clustering with high significance in both action space and metallicity. Stars in the Shards share a common origin, likely as ancient satellite galaxies of the Milky Way, so will be embedded in dark matter (DM) counterparts. These "Dark Shards" contain two substantial streams (S1 and S2), as well as several retrograde, prograde and lower energy substructures. The retrograde stream S1 has a very high Earth-frame speed of $\sim 550$ km s$^{-1}$ while S2 moves on a prograde, but highly polar orbit and enhances the peak of the speed distribution at around $300$ km s$^{-1}$. The presence of the Dark Shards locally leads to modifications of many to the fundamental properties of experimental DM signals. The~S2 stream in particular gives rise to an array of effects in searches for axions and in the time dependence of nuclear recoils: shifting the peak day, inducing nonsinusoidal distortions, and increasing the importance of the gravitational focusing of DM by the Sun. Dark Shards additionally bring new features for directional signals, while also enhancing the DM flux towards Cygnus.

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

Predicting signals in experiments to directly detect dark matter (DM) requires a form for the local DM velocity distribution. Hitherto, the standard halo model (SHM), in which velocities are isotropic and follow a truncated Gaussian law, has performed this job. New data, however, suggest that a substantial fraction of our stellar halo lies in a strongly radially anisotropic population, the 'Gaia Sausage'. Inspired by this recent discovery, we introduce an updated DM halo model, the SHM$^{++}$, which includes a `Sausage' component, thus better describing the known features of our galaxy. The SHM$^{++}$ is a simple analytic model with five parameters: the circular speed, local escape speed and local DM density, which we update to be consistent with the latest data, and two new parameters: the anisotropy and the density of DM in the Sausage. The impact of the SHM$^{++}$ on signal models for WIMPs and axions is rather modest since the multiple changes and updates have competing effects. The largest change occurs for directional detectors which have sensitivity to the full three-dimensional velocity distribution.

The recently discovered S1 stream passes through the Solar neighbourhood on a low inclination, counter-rotating orbit. The progenitor of S1 is a dwarf galaxy with a total mass comparable to the present-day Fornax dwarf spheroidal, so the stream is expected to have a significant DM component. We compute the effects of the S1 stream on WIMP and axion detectors as a function of the density of its unmeasured dark component. In WIMP detectors the S1 stream supplies more high energy nuclear recoils so will marginally improve DM detection prospects. We find that even if S1 comprises less than 10% of the local density, multi-ton xenon WIMP detectors can distinguish the S1 stream from the bulk halo in the relatively narrow mass range between 5 and 25 GeV. In directional WIMP detectors such as CYGNUS, S1 increases DM detection prospects more substantially since it enhances the anisotropy of the WIMP signal. Finally, we show that axion haloscopes possess by far the greatest potential sensitivity to the S1 stream. Once the axion mass has been discovered, the distinctive velocity distribution of S1 can easily be extracted from the axion power spectrum.

Dark matter (DM) particles with mass in the sub-GeV range are an attractive alternative to heavier weakly-interacting massive particles, but direct detection of such light particles is challenging. If however DM-nucleus scattering leads to ionisation of the recoiling atom, the resulting electron may be detected even if the nuclear recoil is unobservable. We demonstrate that including this effect significantly enhances direct detection sensitivity to sub-GeV DM. Existing experiments set world-leading limits, and future experiments may probe the cross sections relevant for thermal freeze-out.

One of the most promising strategies to identify the nature of dark matter consists in the search for new particles at accelerators and with so-called direct detection experiments. Working within the framework of simplified models, and making use of machine learning tools to speed up statistical inference, we address the question of what we can learn about dark matter from a detection at the LHC and a forthcoming direct detection experiment. We show that with a combination of accelerator and direct detection data, it is possible to identify newly discovered particles as dark matter, by reconstructing their relic density assuming they are weakly interacting massive particles (WIMPs) thermally produced in the early Universe, and demonstrating that it is consistent with the measured dark matter abundance. An inconsistency between these two quantities would instead point either towards additional physics in the dark sector, or towards a non-standard cosmology, with a thermal history substantially different from that of the standard cosmological model.

Existing xenon dark matter (DM) direct detection experiments can probe the DM-nucleon interaction of DM with a sub-GeV mass through a search for photon emission from the recoiling xenon atom. We show that LUX's constraints on sub-GeV DM, which utilise the scintillation (S1) and ionisation (S2) signals, are approximately three orders of magnitude more stringent than previous xenon constraints in this mass range, derived from the XENON10 and XENON100 S2-only searches. The new LUX constraints provide the most stringent direct detection constraints for DM particles with a mass below 0.5 GeV. In addition, the photon emission signal in LUX and its successor LZ maintain the discrimination between background and signal events so that an unambiguous discovery of sub-GeV DM is possible. We show that LZ has the potential to reconstruct the DM mass with 20% accuracy for particles lighter than 0.5 GeV.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

Dark matter detectors that utilize liquid xenon have now achieved tonne-scale targets, giving them sensitivity to all flavours of supernova neutrinos via coherent elastic neutrino-nucleus scattering. Considering for the first time a realistic detector model, we simulate the expected supernova neutrino signal for different progenitor masses and nuclear equations of state in existing and upcoming dual-phase liquid xenon experiments. We show that the proportional scintillation signal (S2) of a dual-phase detector allows for a clear observation of the neutrino signal and guarantees a particularly low energy threshold, while the backgrounds are rendered negligible during the supernova burst. XENON1T (XENONnT and LZ; DARWIN) experiments will be sensitive to a supernova burst up to 25 (35; 65) kpc from Earth at a significance of more than 5 sigma, observing approximately 35 (123; 704) events from a 27 Msun supernova progenitor at 10 kpc. Moreover, it will be possible to measure the average neutrino energy of all flavours, to constrain the total explosion energy, and to reconstruct the supernova neutrino light curve. Our results suggest that a large xenon detector such as DARWIN will be competitive with dedicated neutrino telescopes, while providing complementary information that is not otherwise accessible.

Over the past few decades, an anomalous 511 keV gamma-ray line has been observed from the centre of the Milky Way. Dark matter (DM) in the form of light weakly interacting massive particles (WIMPs) annihilating into electron-positron pairs has been one of the leading hypotheses of the observed emission. Given the small required cross section, a further coupling to lighter particles is required to produce the correct relic density. Here, we derive constraints from the Planck satellite on light WIMPs that were in equilibrium with either the neutrino or electron sector in the early universe. For the neutrino sector, we obtain a lower bound on the WIMP mass of 4 MeV for a real scalar and 10 MeV for a Dirac fermion DM particle, at 95% CL. For the electron sector, we find even stronger bounds of 7 and 11 MeV, respectively. Using these results, we show that, in the absence of additional ingredients such as dark radiation, the light thermally produced WIMP explanation of the 511 keV excess is strongly disfavoured by the latest cosmological data. This suggests an unknown astrophysical or more exotic DM source of the signal.

Dark matter can scatter and excite a nucleus to a low-lying excitation in a direct detection experiment. This signature is distinct from the canonical elastic scattering signal because the inelastic signal also contains the energy deposited from the subsequent prompt de-excitation of the nucleus. A measurement of the elastic and inelastic signal will allow a single experiment to distinguish between a spin-independent and spin-dependent interaction. For the first time, we characterise the inelastic signal for two-phase xenon detectors in which dark matter inelastically scatters off the Xe-129 or Xe-131 isotope. We do this by implementing a realistic simulation of a typical tonne-scale two-phase xenon detector and by carefully estimating the relevant background signals. With our detector simulation, we explore whether the inelastic signal from the axial-vector interaction is detectable with upcoming tonne-scale detectors. We find that two-phase detectors allow for some discrimination between signal and background so that it is possible to detect dark matter that inelastically scatters off either the Xe-129 or Xe-131 isotope for dark matter particles that are heavier than approximately 100 GeV. If, after two years of data, the XENON1T search for elastic scattering nuclei finds no evidence for dark matter, the possibility of ever detecting an inelastic signal from the axial-vector interaction will be almost entirely excluded.

Several groups have identified an extended excess of gamma rays over the modeled foreground and background emissions towards the Galactic center (GC) based on observations with the Fermi Large Area Telescope. This excess emission is compatible in morphology and spectrum with a telltale sign from dark matter (DM) annihilation. Here, we present a critical reassessment of DM interpretations of the GC signal in light of the foreground and background uncertainties that some of us recently outlaid in Calore et al. 2014. We find that a much larger number of DM models fits the gamma-ray data than previously noted. In particular: (1) In the case of DM annihilation into $\bar{b}b$, we find that even large DM masses up to $m_\chi \simeq$ 74 GeV are allowed with a $p$-value $> 0.05$. (2) Surprisingly, annihilation into non-relativistic hh gives a good fit to the data. (3) The inverse Compton emission from $\mu^+\mu^-$ with $m_\chi\sim$ 60-70 GeV can also account for the excess at higher latitudes, $|b|>2^\circ$, both in its spectrum and morphology. We also present novel constraints on a large number of mixed annihilation channels, including cascade annihilation involving hidden sector mediators. Finally, we show that the current limits from dwarf spheroidal observations are not in tension with a DM interpretation when uncertainties on the DM halo profile are accounted for.

Previous attempts at explaining the gamma-ray excess near the Galactic Centre have focussed on dark matter annihilation directly into Standard Model particles. This results in a preferred dark matter mass of 30-40 GeV (if the annihilation is into b quarks) or 10 GeV (if it is into leptons). Here we show that the gamma-ray excess is also consistent with heavier dark matter particles; in models of secluded dark matter, dark matter with mass up to 76 GeV provides a good fit to the data. This occurs if the dark matter first annihilates to an on-shell particle that subsequently decays to Standard Model particles through a portal interaction. This is a generic process that works in models with annihilation, semi-annihilation or both. We explicitly demonstrate this in a model of hidden vector dark matter with an SU(2) gauge group in the hidden sector.

We show that it is possible for WIMP dark matter to produce a large signal in indirect dark matter searches without producing signals elsewhere. We illustrate our point by fitting the Fermi-LAT extended galactic gamma-ray excess with a simple model of Dirac dark matter that annihilates primarily into b quarks via a pseudoscalar. Current collider constraints are weak while the 14 TeV LHC run will constrain a limited portion of the parameter space. No signal is expected in additional indirect searches or at future direct detection experiments. Our results emphasise the importance of fully understanding potential indirect signals of dark matter as they may provide the only information about the dark matter particle.

The Earth's velocity relative to the Sun in galactic coordinates is required in the rate calculation for direct detection experiments. We provide a rigorous derivation of this quantity to first order in the eccentricity of the Earth's orbit. We also discuss the effect of the precession of the equinoxes, which has hitherto received little explicit discussion. Comparing with other expressions in the literature, we confirm that the expression of Lee, Lisanti and Safdi is correct, while the expression of Lewin and Smith, the de facto standard expression, contains an error. For calculations of the absolute event rate, the leading order expression is sufficient while for modulation searches, an expression with the eccentricity is required for accurate predictions of the modulation phase.

The CDMS-II collaboration has reported 3 events in a Si detector, which are consistent with being nuclear recoils due to scattering of Galactic dark matter particles with a mass of about 8.6 GeV and a cross-section on neutrons of about 2 x 10^-41 cm^2. While a previous result from the XENON10 experiment has supposedly ruled out such particles as dark matter, we find by reanalysing the XENON10 data that this is not the case. Some tension remains however with the upper limit placed by the XENON100 experiment, independently of astrophysical uncertainties concerning the Galactic dark matter distribution. We explore possible ways of ameliorating this tension by altering the properties of dark matter interactions. Nevertheless, even with standard couplings, light dark matter is consistent with both CDMS and XENON10/100.

Motivated by recent claims of lines in the Fermi gamma-ray spectrum, we critically examine means of enhancing neutralino annihilation into neutral gauge bosons. The signal can be boosted while remaining consistent with continuum photon constraints if a new singlet-like pseudoscalar is present. We consider singlet extensions of the MSSM, focusing on the NMSSM, where a `well-tempered' neutralino can explain the lines while remaining consistent with current constraints. We adopt a complementary numerical and analytic approach throughout in order to gain intuition for the underlying physics. The scenario requires a rich spectrum of light neutralinos and charginos leading to characteristic phenomenological signatures at the LHC whose properties we explore. Future direct detection prospects are excellent, with sizeable spin-dependent and spin-independent cross-sections.

Recent work on increasing the effective number of neutrino species (Neff) in the early universe has focussed on introducing extra relativistic species (`dark radiation'). We draw attention to another possibility: a new particle of mass less than 10 MeV that remains in thermal equilibrium with neutrinos until it becomes non-relativistic increases the neutrino temperature relative to the photons. We demonstrate that this leads to a value of Neff that is greater than three and that Neff at CMB formation is larger than at BBN. We investigate the constraints on such particles from the primordial abundance of helium and deuterium created during BBN and from the CMB power spectrum measured by ACT and SPT and find that they are presently relatively unconstrained. We forecast the sensitivity of the Planck satellite to this scenario: in addition to dramatically improving constraints on the particle mass, in some regions of parameter space it can discriminate between the new particle being a real or complex scalar.

We study the impact of the assumed velocity distribution of galactic dark matter particles on the interpretation of results from nuclear recoil detectors. By converting experimental data to variables that make the astrophysical unknowns explicit, different experiments can be compared without implicit assumptions concerning the dark matter halo. We extend this framework to include the annual modulation signal, as well as multiple target elements. Recent results from DAMA, CoGeNT and CRESST-II can be brought into agreement if the velocity distribution is very anisotropic and thus allows a large modulation fraction. However constraints from CDMS and XENON cannot be evaded by appealing to such astrophysical uncertainties alone.

Motivated by the recent interest in light WIMPs of mass ~O(10 GeV), an extension of the elastic, spin-independent WIMP-nucleon cross-section limits resulting from the CRESST-II commissioning run (2007) are presented. Previously, these data were used to set cross-section limits from 1000 GeV down to ~17 GeV, using tungsten recoils, in 47.9 kg-days of exposure of calcium tungstate. Here, the overlap of the oxygen and calcium bands with the acceptance region of the commissioning run data set is reconstructed using previously published quenching factors. The resulting elastic WIMP cross section limits, accounting for the additional exposure of oxygen and calcium, are presented down to 5 GeV.

The CoGeNT collaboration has reported evidence of an annual modulation in its first fifteen months of data. Here we compare the amplitude and phase of this signal to the modulation observed by the DAMA collaboration, assuming that both arise due to elastically scattering dark matter (DM). We directly map the CoGeNT signal to the DAMA detector without specifying any astrophysical parameters and compare this with the signal measured by DAMA. We also compare with constraints from CDMS II and XENON10. We find that DM of mass 5-14 GeV that couples equally to protons and neutrons is strongly disfavoured. Isospin-violating DM fares better but requires a boosted modulation fraction.

DAMA observes an annual modulation in their event rate, as might be expected from dark matter scatterings, while CoGeNT has reported evidence for a similar modulation. The simplest interpretation of these findings in terms of dark matter-nucleus scatterings is excluded by other direct detection experiments. We consider the robustness of these exclusions with respect to assumptions regarding the scattering and find that isospin-violating inelastic dark matter helps alleviate this tension and allows marginal compatibility between experiments. Isospin-violation can significantly weaken the XENON constraints, while inelasticity enhances the annual modulation fraction of the signal, bringing the CoGeNT and CDMS results into better agreement.

The effects of astrophysical uncertainties on the exclusion limits at dark matter direct detection experiments are investigated for three scenarios: elastic, momentum dependent and inelastically scattering dark matter. We find that varying the dark matter galactic escape velocity and the Sun's circular velocity can lead to significant variations in the exclusion limits for light ($\lesssim10$ GeV) elastic and inelastic scattering dark matter. We also calculate the limits using one hundred velocity distributions extracted from the Via Lactea II and GHALO N-body simulations and find that a Maxwell-Boltzmann distribution with the same astrophysical parameters generally sets less constraining limits. The elastic and momentum dependent limits remain robust for masses $\gtrsim50$ GeV under variations of the astrophysical parameters and the form of the velocity distribution.

Aims. The brown dwarf (BD) formation process has not yet been completely understood. To shed more light on the differences and similarities between star and BD formation processes, we study and compare the disk fraction among both kinds of objects over a large angular region in the Taurus cloud. In addition, we examine the spatial distribution of stars and BD relative to the underlying molecular gas Methods. In this paper, we present new and updated photometry data from the Infrared Array Camera (IRAC) aboard the Spitzer Space Telescope on 43 BDs in the Taurus cloud, and recalculate of the BD disk fraction in this region. We also useed recently available CO mm data to study the spatial distribution of stars and BDs relative to the cloud's molecular gas. Results. We find that the disk fraction among BDs in the Taurus cloud is 41 \pm 12%, a value statistically consistent with the one among TTS (58 \pm 9%). We find that BDs in transition from a state where they have a disk to a diskless state are rare, and we study one isolated example of a transitional disk with an inner radius of \approx 0.1 AU (CFHT BD Tau 12, found via its relatively small mid-IR excess compared to most members of Taurus that have disks. We find that BDs are statistically found in regions of similar molecular gas surface density to those associated with stars. Furthermore, we find that the gas column density distribution is almost identical for stellar and substellar objects with and without disks.

Using data from the SCUBA Legacy Catalogue (850 um) and Spitzer Space Telescope (3.6 - 70 um), we explore dense cores in the Ophiuchus, Taurus, Perseus, Serpens, and Orion molecular clouds. We develop a new method to discriminate submillimeter cores found by SCUBA as starless or protostellar, using point source photometry from Spitzer wide field surveys. First, we identify infrared sources with red colors associated with embedded young stellar objects (YSOs). Second, we compare the positions of these YSO-candidates to our submillimeter cores. With these identifications, we construct new, self-consistent starless and protostellar core mass functions (CMFs) for the five clouds. We find best fit slopes to the high-mass end of the CMFs of -1.26 +/- 0.20, -1.22 +/- 0.06, -0.95 +/- 0.20, and -1.67 +/- 0.72 for Ophiuchus, Taurus, Perseus, and Orion, respectively. Broadly, these slopes are each consistent with the -1.35 power-law slope of the Salpeter IMF at higher masses, but suggest some differences. We examine a variety of trends between these CMF shapes and their parent cloud properties, potentially finding a correlation between the high-mass slope and core temperature. We also find a trend between core mass and effective size, but we are very limited by sensitivity. We make similar comparisons between core mass and size with visual extinction (for A_V >= 3) and find no obvious trends. We also predict the numbers and mass distributions of cores that future surveys with SCUBA-2 may detect in each of these clouds.

A simple theory of supersymmetric dark matter (DM) naturally linked to neutrino flavour physics is studied. The DM sector comprises a spectrum of mixed lhd-rhd sneutrino states where both the sneutrino flavour structure and mass splittings are determined by the associated neutrino masses and mixings. Prospects for indirect detection from solar capture are good due to a large sneutrino-nucleon cross-section afforded by the inelastic splitting (solar capture limits exclude an explanation of DAMA/LIBRA). We find parameter regions where all heavier states will have decayed, leaving only one flavour mixture of sneutrino as the candidate DM. Such regions have a unique `smoking gun' signature--sneutrino annihilation in the Sun produces a pair of neutrino mass eigenstates free from vacuum oscillations, with the potential for detection at neutrino telescopes through the observation of a hard spectrum of nu_mu and nu_tau (for a normal neutrino hierarchy). Next generation direct detection experiments can explore much of the parameter space through both elastic and inelastic scattering. We show in detail that the observed neutrino masses and mixings can arise as a consequence of supersymmetry breaking effects in the sneutrino DM sector, consistent with all experimental constraints.

We present new high spatial resolution (<~ 0.1") 1-5 micron adaptive optics images, interferometric 1.3 mm continuum and 12CO 2-1 maps, and 350 micron, 2.8 and 3.3 mm fluxes measurements of the HV Tau system. Our adaptive optics images reveal an unusually slow orbital motion within the tight HV Tau AB pair that suggests a highly eccentric orbit and/or a large deprojected physical separation. Scattered light images of the HV Tau C edge-on protoplanetary disk suggest that the anisotropy of the dust scattering phase function is almost independent of wavelength from 0.8 to 5 micron, whereas the dust opacity decreases significantly over the same range. The images further reveal a marked lateral asymmetry in the disk that does not vary over a timescale of 2 years. We further detect a radial velocity gradient in the disk in our 12CO map that lies along the same position angle as the elongation of the continuum emission, which is consistent with Keplerian rotation around an 0.5-1 Msun central star, suggesting that it could be the most massive component in the triple system. We use a powerful radiative transfer model to compute synthetic disk observations and use a Bayesian inference method to extract constraints on the disk properties. Each individual image, as well as the spectral energy distribution, of HV Tau C can be well reproduced by our models with fully mixed dust provided grain growth has already produced larger-than-interstellar dust grains. However, no single model can satisfactorily simultaneously account for all observations. We suggest that future attempts to model this source include more complex dust properties and possibly vertical stratification. (Abridged)

We report on the properties of pre-main-sequence objects in the Taurus molecular clouds as observed in 7 mid- and far-infrared bands with the Spitzer Space Telescope. There are 215 previously-identified members of the Taurus star-forming region in our ~44 square degree map; these members exhibit a range of Spitzer colors that we take to define young stars still surrounded by circumstellar dust (noting that ~20% of the bonafide Taurus members exhibit no detectable dust excesses). We looked for new objects in the survey field with similar Spitzer properties, aided by extensive optical, X-ray, and ultraviolet imaging, and found 148 candidate new members of Taurus. We have obtained follow-up spectroscopy for about half the candidate sample, thus far confirming 34 new members, 3 probable new members, and 10 possible new members, an increase of 15-20% in Taurus members. Of the objects for which we have spectroscopy, 7 are now confirmed extragalactic objects, and one is a background Be star. The remaining 93 candidate objects await additional analysis and/or data to be confirmed or rejected as Taurus members. Most of the new members are Class II M stars and are located along the same cloud filaments as the previously-identified Taurus members. Among non-members with Spitzer colors similar to young, dusty stars are evolved Be stars, planetary nebulae, carbon stars, galaxies, and AGN.

Circumstellar discs are expected to be the nursery of planets. Grain growth within such discs is the first step in the planet formation process. The Bok globule CB 26 harbours such a young disc. We present a detailed model of the edge-on circumstellar disc and its envelope in the Bok globule CB 26. The model is based on HST near-infrared maps in the I, J, H, and K bands, OVRO and SMA radio maps at 1.1mm, 1.3mm and 2.7mm, and the spectral energy distribution (SED) from 0.9 microns to 3mm. New photometric and spectroscopic data from the Spitzer Space Telescope and the Caltech Submilimeter Observatory have been obtained and are part of our analysis. Using the self-consistent radiative transfer code MC3D, the model we construct is able to discriminate parameter sets and dust properties of both its parts, namely envelope and disc. We find that the disc has an inner hole with a radius of 45 +/- 5 AU. Based on a dust model including silicate and graphite the maximum grain size needed to reproduce the spectral millimetre index is 2.5 microns. Features seen in the near-infrared images, dominated by scattered light, can be described as a result of a rotating envelope. Successful employment of ISM dust in both the disc and envelope hint that grain growth may not yet play a significant role for the appearance of this system. A larger inner hole gives rise to the assumption that CB 26 is a circumbinary disc.

We present Spitzer images of the Taurus Complex (TC) and take advantage of the sensitivity and spatial resolution of the observations to characterize the diffuse IR emission across the cloud. This work highlights evidence of dust evolution within the translucent sections of the archetype reference for studies of quiescent molecular clouds. We combine Spitzer 160 um and IRAS 100 um observations to produce a dust temperature map and a far-IR dust opacity map at 5' resolution. The average dust temperature is about 14.5K with a dispersion of +/-1K across the cloud. The far-IR dust opacity is a factor 2 larger than the average value for the diffuse ISM. This opacity increase and the attenuation of the radiation field (RF) both contribute to account for the lower emission temperature of the large grains. The structure of the TC significantly changes in the mid-IR images that trace emission from PAHs and VSGs. We focus our analysis of the mid-IR emission to a range of ecliptic latitudes where the zodiacal light residuals are small. Within this cloud area, there are no 8 and 24 um counterparts to the brightest 160 um emission features. Conversely, the 8 and 24 um images reveal filamentary structure that is strikingly inconspicuous in the 160 um and extinction maps. The IR colors vary over sub-parsec distances across this filamentary structure. We compare the observed colors with model calculations quantifying the impact of the RF intensity and the abundance of stochastically heated particles on the dust SED. To match the range of observed colors, we have to invoke variations by a factor of a few of both the interstellar RF and the abundance of PAHs and VSGs. We conclude that within this filamentary structure a significant fraction of the dust mass cycles in and out the small size end of the dust size distribution.

We combine the results from several multiplicity surveys of pre-main-sequence stars located in four nearby star-forming regions with Spitzer data from three different Legacy Projects. This allows us to construct a sample of 349 targets, including 125 binaries, which we use to to investigate the effect of companions on the evolution of circumstellar disks. We find that the distribution of projected separations of systems with Spitzer excesses is significantly different (P ~2.4e-5, according to the KS test for binaries with separations < 400 AU) from that of systems lacking evidence for a disk. As expected, systems with projected separations < 40 AU are half as likely to retain at least one disk than are systems with projected separations in the 40-400 AU range. These results represent the first statistically significant evidence for a correlation between binary separation and the presence of an inner disk (r ~ 1 AU). Several factors (e.g., the incompleteness of the census of close binaries, the use of unresolved disk indicators, and projection effects) have previously masked this correlation in smaller samples. We discuss the implications of our findings for circumstellar disk lifetimes and the formation of planets in multiple systems.

The DAMA collaboration have claimed to detect particle dark matter (DM) via an annual modulation in their observed recoil event rate. This appears to be in strong disagreement with the null results of other experiments if interpreted in terms of elastic DM scattering, while agreement for a small region of parameter space is possible for inelastic DM (iDM) due to the altered kinematics of the collision. To date most analyses assume a simple galactic halo DM velocity distribution, the Standard Halo Model, but direct experimental support for the SHM is severely lacking and theoretical studies indicate possible significant differences. We investigate the dependence of DAMA and the other direct detection experiments on the local DM velocity distribution, utilizing the results of the Via Lactea and Dark Disc numerical simulations. We also investigate effects of varying the solar circular velocity, the DM escape velocity, and the DAMA quenching factor within experimental limits. Our data set includes the latest ZEPLIN-III results, as well as full publicly available data sets. Due to the more sensitive dependence of the inelastic cross section on the velocity distribution, we find that with Via Lactea the DAMA results can be consistent with all other experiments over an enlarged region of iDM parameter space, with higher mass particles being preferred, while Dark Disc does not lead to an improvement. A definitive test of DAMA for iDM requires heavy element detectors.

We investigate the environment of the very low-luminosity object L1521F IRS using data from the Taurus Spitzer Legacy Survey. The MIPS 160 micron image shows both extended emission from the Taurus cloud as well as emission from multiple cold cores over a 1 X 2 deg region. Analysis shows that the cloud dust temperature is 14.2 +- 0.4 K and the extinction ratio is A_160/A_K = 0.010 +- 0.001 up to A_V ~ 4 mag. We find kappa_160 = 0.23 +- 0.046 cm^2/g for the specific opacity of the gas-dust mixture. Therefore, for dust in the Taurus cloud we find the 160 um opacity is significantly higher than that measured for the diffuse ISM, but not too different from dense cores, even at modest extinction values. Furthermore, the 160 um image shows features that do not appear in the IRAS 100 um image. We identify six regions as cold cores, i.e. colder than 14.2 K, all of which have counterparts in extinction maps or C18O maps. We compare the effects of L1521F IRS on its natal core and find there is no evidence for dust heating at 160 or 100 um by the embedded source. From the infrared luminosity L_TIR = 0.024 Lo we find L_bol = 0.034 - 0.046L_o, thus confirming the source's low-luminosity. Comparison of L1521F IRS with theoretical simulations for the very early phases of star formation appears to rule out the first core collapse phase. The evolutionary state appears similar to or younger than the class 0 phase, and the estimated mass is likely to be substellar.

We present the Chandra High Energy Transmission Grating Spectrometer (HETGS) and Keck observations of HDE 245059, a young weak-lined T Tauri star (WTTS), member of the pre-main sequence group in the Lambda Orionis Cluster. Our high spatial resolution, near-infrared observations with Keck reveal that HDE 245059 a binary separated by 0.87". Based on this new information we have obtained an estimate of the masses of the binary components; 3M_{sun} and 2.5M_{sun} for the north and south components, respectively. We have estimated the age of the system to be ~2-3 Myr. We detect both components of the binary in the zeroth order Chandra image and in the grating spectra. Our fits to the spectrum of the binary have shown that the emission is dominated by a plasma between 8 and 15 MK, a soft component at 4 MK and a hard component at 50 MK are also detected. The value of the hydrogen column density was low, 8 x 10^{19} cm^{-2}, likely due to the clearing of the inner region of the Lambda Orionis cloud. The abundance pattern shows an inverse First Ionization Potential (FIP) effect for all elements from O to Fe, the only exception being Ca. A 3-T model was fitted to the individual zeroth order spectra using the abundances derived for the binary. We have also obtained several lines fluxes from the grating spectra. The fits to the triplets show no evidence of high densities. We conclude that the X-ray properties of the weak-lined T Tau binary HDE 245059 are similar to those generally observed in other weak-lined T Tau stars. Although its accretion history may have been affected by the clearing of the interstellar material around Lambda Ori, its coronal properties appears not to have been strongly modified.

We present maps of 14.4 deg^2 of the Ophiuchus dark clouds observed by the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). These high quality maps depict both numerous point sources as well as extended dust emission within the star-forming and non-star-forming portions of these clouds. Using PSF-fitting photometry, we detect 5779 sources at 24 um and 81 sources at 70 um at the 10 sigma level of significance. Three hundred twenty-three candidate young stellar objects (YSOs) were identified according to their positions on the MIPS/2MASS K versus K$-$[24] color-magnitude diagrams as compared to 24 um detections in the SWIRE extragalactic survey. We find that more than half of the YSO candidates, and almost all the ones with protostellar Class I spectral energy distributions, are confined to the known cluster and aggregates.

We present the results of a five year monitoring campaign of the close binary TWA 5Aab in the TW Hydrae association, using speckle and adaptive optics on the W.M. Keck 10 m telescopes. These measurements were taken as part of our ongoing monitoring of pre-main sequence (PMS) binaries in an effort to increase the number of dynamically determined PMS masses and thereby calibrate the theoretical PMS evolutionary tracks. Our observations have allowed us to obtain the first determination of this system's astrometric orbit. We find an orbital period of 5.94 +- 0.09 years and a semi-major axis of 0.066" +- 0.005". Combining these results with a kinematic distance, we calculate a total mass of 0.71 +- 0.14 M_sun (D/44 pc)^3. for this system. This mass measurement, as well as the estimated age of this system, are consistent to within 2$\sigma$ of all theoretical models considered. In this analysis, we properly account for correlated uncertainties, and show that while these correlations are generally ignored, they increase the formal uncertainties by up to a factor of five and therefore are important to incorporate. With only a few more years of observation, this type of measurement will allow the theoretical models to be distinguished.

The formation of a binary system surrounded by disks is the most common outcome of stellar formation. Hence studying and understanding the formation and the evolution of binary systems and associated disks is a cornerstone of star formation science. Moreover, since the components within binary systems are coeval and the sizes of their disks are fixed by the tidal truncation of their companion, binary systems provide an ideal "laboratory" in which to study disk evolution under well defined boundary conditions. In this paper, we review observations of several inner disk diagnostics in multiple systems, including hydrogen emission lines (indicative of ongoing accretion), $K-L$ and $K-N$ color excesses (evidence of warm inner disks), and polarization (indicative of the relative orientations of the disks around each component). We examine to what degree these properties are correlated within binary systems and how this degree of correlation depends on parameters such as separation and binary mass ratio. These findings will be interpreted both in terms of models that treat each disk as an isolated reservoir and those in which the disks are subject to re-supply from some form of circumbinary reservoir, the observational evidence for which we will also critically review. The planet forming potential of multiple star systems is discussed in terms of the relative lifetimes of disks around single stars, binary primaries and binary secondaries. Finally, we summarize several potentially revealing observational problems and future projects that could provide further insight into disk evolution in the coming decade

We present a high spatial resolution, 10-20 micron survey of 65 T Tauri binary stars in Taurus, Ophiuchus, and Corona Australis using the Keck 10 m telescopes. Designed to probe the inner ~1 AU region of the circumstellar disks around the individual stellar components in these binary systems, this study increases the number of binaries with spatially resolved measurements at 10 micron by a factor of ~5. Combined with resolved near-infrared photometry and spectroscopic accretion diagnostics, we find that ~10% of stars with a mid-infrared excess do not appear to be accreting. In contrast to an actively accreting disk system, these passive disks have significantly lower near-infrared colors that are, in most cases, consistent with photospheric emission, suggesting the presence of an inner disk hole. In addition, there appears to be a spectral type/mass dependence associated with the presence of a passive disk, with all passive disks occurring around M type stars. The possibility that the passive disks are caused by the presence of an as yet undetected companion at a small separation (0.3-3 AU) is possible for any individual system, however, it cannot account for the spectral type dependence of the passive disk sample as a whole. We propose that these passive disks represent a subset of T Tauri stars that are undergoing significant disk evolution. The fraction of observed passive disks and the observed spectral type dependence can both be explained by models of disk evolution that include disk photoevaporation from the central star. (abridged).

We have obtained the first high spatial (0.05'') and spectral (R~35000) resolution 2 micron spectrum of the T Tau S tight binary system using adaptive optics on the Keck II telescope. We have also obtained the first 3.8 and 4.7 micron images that resolve the three components of the T Tau multiple system, as well as new 1.6 and 2.2 micron images. Together with its very red near-infrared colors, the spectrum of T Tau Sb shows that this T Tauri star is extincted by a roughly constant extinction of Av~15 mag, which is probably the 0.7''x0.5'' circumbinary structure recently observed in absorption in the ultraviolet. T Tau Sa, which is also observed through this screen and is actively accreting, further possesses a small edge-on disk that is evidenced by warm (390 K), narrow overtone CO rovibrational absorption features in our spectrum. We find that T Tau Sa is most likely an intermediate-mass star surrounded by a semi-transparent 2-3 AU-radius disk whose asymmetries and short Keplerian rotation explain the large photometric variability of the source on relatively short timescales. We also show that molecular hydrogen emission exclusively arises from the gas that surrounds T Tau S and that its spatial and kinematic structure, while providing suggestive evidence for a jet-like structure, is highly complex.

The explosion that results in a cosmic gamma-ray burst (GRB) is thought to produce emission from two physical processes -- the activity of the central engine gives rise to the high-energy emission of the burst through internal shocking and the subsequent interaction of the flow with the external environment produces long-wavelength afterglow. While afterglow observations continue to refine our understanding of GRB progenitors and relativistic shocks, gamma-ray observations alone have not yielded a clear picture of the origin of the prompt emission nor details of the central engine. Only one concurrent visible-light transient has been found and was associated with emission from an external shock. Here we report the discovery of infrared (IR) emission contemporaneous with a GRB, beginning 7.2 minutes after the onset of GRB 041219a. Our robotic telescope acquired 21 images during the active phase of the burst, yielding the earliest multi-colour observations of any long-wavelength emission associated with a GRB. Analysis of an initial IR pulse suggests an origin consistent with internal shocks. This opens a new possibility to study the central engine of GRBs with ground-based observations at long wavelengths.

We present here the results of astrometric, photometric and spectroscopic observations leading to the determination of the orbit and dynamical masses of the binary L dwarf 2MASSW J0746425+2000321. High angular resolution observations spread over almost 4 years and obtained with the Hubble Space Telescope (HST), the ESO Very Large Telescope (VLT), and a the W. M. Keck Observatory (Keck) allow us to cover 36% of the period, corresponding to 60% of the orbit, and, for the first time, to derive a precise estimate of the total and individual masses of such a late-type object. We find an orbital period of 3850.9$^{+904}_{-767}$ days. The corresponding total mass is 0.146$^{+0.016}_{-0.006}$ M$_{\sun}$, with uncertainties depending on the distance. Spatially resolved low resolution optical (550--1025 nm) spectra have been obtained with HST/STIS, allowing us to measure the spectral types of the two components (L0$\pm$0.5 for the primary and L1.5$\pm$0.5 for the secondary). We also present precise photometry of the individual components measured on the high angular resolution images obtained with HST/ACS and WFPC2 (visible), VLT/NACO (J, H and Ks bands) and Keck I (Ks) band). These spectral and photometric measurements enable us to estimate their effective temperatures and mass ratio, and to place the object accurately in a H--R diagram. The binary system is most likely formed by a primary with a mass of 0.085$\pm$0.010 M$_{\sun}$ and a secondary with a mass of 0.066$\pm$0.006 M$_{\sun}$, thus clearly substellar, for an age of approximately 300$\pm$150 Myr. H$\alpha$ variability indicates chromospheric and/or magnetic activity.

We present the first 3.8 micron image of the dusty ring surrounding the young binary system GG Tau, obtained with the W. M. Keck II 10m telescope's adaptive optics system. This is the longest wavelength at which the ring has been detected in scattered light so far, allowing a multi-wavelength analysis of the scattering properties of the dust grains present in this protoplanetary disk in combination with previous, shorter wavelengths, HST images. We find that the scattering phase function of the dust grains in the disk is only weakly dependent on the wavelength. This is inconsistent with dust models inferred from observations of the interstellar medium or dense molecular clouds. In particular, the strongly forward-throwing scattering phase function observed at 3.8 micron implies a significant increase in the population of large (>~1 micron) grains, which provides direct evidence for grain growth in the ring. However, the grain size distribution required to match the 3.8 micron image of the ring is incompatible with its published 1 micron polarization map, implying that the dust population is not uniform throughout the ring. We also show that our 3.8 micron scattered light image probes a deeper layer of the ring than previous shorter wavelength images, as demonstrated by a shift in the location of the inner edge of the disk's scattered light distribution between 1 and 3.8 micron. We therefore propose a stratified structure for the ring in which the surface layers, located ~50 AU above the ring midplane, contain dust grains that are very similar to those found in dense molecular clouds, while the region of the ring located ~25 AU from the midplane contains significantly larger grains. This stratified structure is likely the result of vertical dust settling and/or preferred grain growth in the densest parts of the ring.

We summarize here the observed properties of ``infrared companions'' to T Tauri stars and argue that their observational properties are identical to those of Class I sources. They may therefore be embedded protostars in a much earlier evolutionary phase than T Tauri stars, in which case these multiple systems are significantly non-coeval as opposed to the majority of young binary systems. They would have formed through a different mechanism than core fragmentation. The only distinction between IRCs and Class I sources is that they lie within a few tens of AU of a T Tauri star, and so they cannot be at the center of a vast optically thick envelope as is believed to be the case for protostars. We discuss whether systems with an IRC are really candidates for non-coeval multiple systems.

We report spatially resolved 11.8 micron images, obtained at the W. M. Keck 10 m telescope, of the protoplanetary disk around the pre--main-sequence star HK Tau B. The mid-infrared morphology and astrometry of HK Tau B with respect to HK Tau A indicate that the flux observed in the mid-infrared from HK Tau B has been scattered off the upper surface of its nearly edge-on disk. This is the first example of a protoplanetary disk observed in scattered light at mid-infrared wavelengths. Monte Carlo simulations of this disk show that the extent (FWHM =0."5, or 70 AU) of the scattered light nebula in the mid-infrared is very sensitive to the dust size distribution. The 11.8 micron measurement can be best modelled by a dust grain population that contains grains on the order of 1.5-3 micron in size; grain populations with exclusively sub-micron grain sizes or power law size distributions that extend beyond 5 micron cannot reproduce the observed morphology. These grains are significantly larger than those expected in the ISM implying that grain growth has occurred; whether this growth is a result of dust evolution within the disk itself or had originally occurred within the dark cloud remains an open question.

We present new multi-epoch near-infrared and optical high-angular images of the V773 Tau pre-main sequence triple system, a weak-line T Tauri (WTTS) system in which the presence of an evolved, ``fossil'' protoplanetary disk has been inferred on the basis of a significant infrared excess. Our images reveal a fourth object bound to the system, V773 Tau D. While it is much fainter than all other components at 2 micron, it is the brightest source in the system at 4.7 micron. We also present medium-resolution K band adaptive optics spectroscopy of this object, which is featureless with the exception of a weak Br gamma emission line. Based on this spectrum and on the spectral energy distribution of the system, we show that V773 Tau D is another member of the small class of ``infrared companions'' (IRCs) to T Tauri stars. It is the least luminous, and probably the least massive, component of the system, as opposed to most other IRCs, which suggests that numerous low-luminosity IRCs such as V773 Tau D may still remain to be discovered. Furthermore, it is the source of the strong IR excess in the system. We therefore reject the interpretation of this excess as the signature of a fossil (or ``passive'') disk and further suggest that these systems may be much less frequent than previously thought. We further show that V773 Tau C is a variable classical T Tauri star (CTTS) and that its motion provides a well constrained orbital model. We show that V773 Tau D can be dynamically stable within this quadruple system if its orbit is highly inclined. Finally, V773 Tau is the first multiple system to display such a variety of evolutionary states (WTTS, CTTS, IRC), which may be the consequence of the strong star-star interactions in this compact quadruple system.

A combination of high-resolution and wide-field imaging reveals two binary stars and one triple star system among the sample of the first 11 stars with planets detected by radial velocity variations. High resolution speckle or adaptive optics (AO) data probe subarcsecond scales down to the diffraction limit of the Keck 10-m or Lick 3-m, and direct images or AO images are sensitive to a wider field, extending to 10" or 38", depending upon the camera. One of the binary systems -- HD 114762 -- was not previously known to be a spatially resolved multiple system; additional data taken with the combination of Keck adaptive optics and NIRSPEC are used to characterize the new companion. The second binary system -- Tau Boo -- was a known multiple with two conflicting orbital solutions; the current data will help constrain the discrepant estimates of periastron time and separation. Another target -- 16 Cyg B -- was a known common proper motion binary, but the current data resolve a new third component, close to the wide companion 16 Cyg A. Both the HD 114762 and 16 Cyg B systems harbor planets in eccentric orbits, while the Tau Boo binary contains an extremely close planet in a tidally-circularized orbit. Although the sample is currently small, the proportion of binary systems is comparable to that measured in the field over a similar separation range. Incorporating the null result from another companion search project lowers the overall fraction of planets in binary systems, but the detections in our survey reveal that planets can form in binaries separated by less than 50 AU.

We present deep, near-infrared images of the circumbinary disk surrounding the pre-main-sequence binary star, GG Tau A, obtained with NICMOS aboard the Hubble Space Telescope. The spatially resolved proto-planetary disk scatters roughly 1.5% of the stellar flux, with a near-to-far side flux ratio of ~1.4, independent of wavelength, and colors that are comparable to the central source; all of these properties are significantly different from the earlier ground-based observations. New Monte Carlo scattering simulations of the disk emphasize that the general properties of the disk, such as disk flux, near side to far side flux ratio and integrated colors, can be approximately reproduced using ISM-like dust grains, without the presence of either circumstellar disks or large dust grains, as had previously been suggested. A single parameter phase function is fitted to the observed azimuthal variation in disk flux, providing a lower limit on the median grain size of 0.23 micron. Our analysis, in comparison to previous simulations, shows that the major limitation to the study of grain growth in T Tauri disk systems through scattered light lies in the uncertain ISM dust grain properties. Finally, we use the 9 year baseline of astrometric measurements of the binary to solve the complete orbit, assuming that the binary is coplanar with the circumbinary ring. We find that the estimated 1 sigma range on disk inner edge to semi-major axis ratio, 3.2 < Rin/a < 6.7, is larger than that estimated by previous SPH simulations of binary-disk interactions.

We obtained new near-infrared images of the prototypical pre-main sequence triple system T Tau, as well as the first resolved medium-resolution spectra of the close pair T Tau S. At the time of our observations, the tight binary had a 13 AU projected separation and showed significant motion since its discovery, three years before. The orbit cannot be strongly constrained yet, but the observed motion of T Tau Sb with respect to T Tau Sa suggests that the system is at least as massive as T Tau N itself. This may indicate that T Tau N is not the most massive star in the system. The spectrum of T Tau Sa, which is totally featureless except for a strong Br gamma emission line, identifies this component with the ``infrared companion'', whose exact nature remains obscure but may be the consequence of it being the most massive component of the system. Contrasting sharply with T Tau Sa, the spectrum of T Tau Sb shows numerous photospheric features consistent with an early-M spectral type. The presence of a strong Br gamma emission line and of a significant veiling continuum classifies this object as a deeply embedded T Tauri star. From these observations, we conclude that both components of T Tau S are embedded in their own dense circumstellar cocoon of material, which are probably fed by a much more extended structure.