Dark matter particles from the Galactic halo can be gravitationally trapped in the solar core or in external orbits. The enhanced density of dark matter particles either in the solar core or in external orbits can result in the annihilation of these particles producing gamma rays via long-lived intermediate states or directly outside the Sun, respectively. These processes would yield characteristic features in the energy spectrum of the subsequent gamma rays, i.e., a box-like or line-like shaped feature, respectively. We have performed a dedicated analysis using a 10-years sample of gamma-ray events from the Sun collected by the Fermi Large Area Telescope searching for spectral features in the energy spectrum as a signature of dark matter annihilation. In the scenario of gamma-ray production via long-lived mediators we have also evaluated the dark matter-nucleon spin-dependent and spin-independent scattering cross section constraints from the flux limits in a dark matter mass range from 3 GeV/c$^2$ up to about 1.8 TeV/c$^2$. In the mass range up to about 150 GeV/c$^2$ the limits are in the range $10^{-46} - 10^{-45}$ cm$^{2}$ for the spin-dependent scattering and in the range $10^{-48} - 10^{-47}$ cm$^{2}$ for the spin-independent case. The range of variation depends on the decay length of the mediator.

Calculating the abundance of thermally produced dark matter particles has become a standard procedure, with sophisticated methods guaranteeing a precision that matches the percent-level accuracy in the observational determination of the dark matter density. Here, we point out that one of the main assumptions in the commonly adopted formalism, namely local thermal equilibrium during the freeze-out of annihilating dark matter particles, does not have to be satisfied in general. We present two methods for how to deal with such situations, in which the kinetic decoupling of dark matter happens so early that it interferes with the chemical decoupling process: i) an approximate treatment in terms of a coupled system of differential equations for the leading momentum moments of the dark matter distribution, and ii) a full numerical solution of the Boltzmann equation in phase-space. For illustration, we apply these methods to the case of Scalar Singlet dark matter. We explicitly show that even in this simple model the prediction for the dark matter abundance can be affected by up to one order of magnitude compared to the traditional treatment.

We use 7 years of electron and positron Fermi-LAT data to search for a possible excess in the direction of the Sun in the energy range from 42 GeV to 2 TeV. In the absence of a positive signal we derive flux upper limits which we use to constrain two different dark matter (DM) models producing $e^+ e^-$ fluxes from the Sun. In the first case we consider DM model being captured by the Sun due to elastic scattering and annihilation into $e^+ e^-$ pairs via a long-lived light mediator that can escape the Sun. In the second case we consider instead a model where DM density is enhanced around the Sun through inelastic scattering and the DM annihilates directly into $e^+ e^-$ pairs. In both cases we perform an optimal analysis, searching specifically for the energy spectrum expected in each case, i.e., a box-like shaped and line-like shaped spectrum respectively. No significant signal is found and we can place limits on the spin-independent cross-section in the range from $10^{-46}~cm^2$ to $10^{-44}~cm^2$ and on the spin-dependent cross-section in the range from $10^{-43}~cm^2$ to $10^{-41}~cm^2$. In the case of inelastic scattering the limits on the cross-section are in the range from $10^{-43}~cm^2$ to $10^{-41}~cm^2$. The limits depend on the life time of the mediator (elastic case) and on the mass splitting value (inelastic case), as well as on the assumptions made for the size of the deflections of electrons and positrons in the interplanetary magnetic field.

Kinetic decoupling of dark matter typically happens much later than chemical freeze-out. In fact, local thermal equilibrium is an important assumption for the usual relic density calculations based on solving the Boltzmann equation (for its 0-th moment) describing the dark matter number density. But is this assumption always justified? Here we address this question and discuss the consequences of more accurate treatments. The first treatment is relying on the inclusion of higher moments of the Boltzmann equation and the second on solving the evolution of the phase-space distribution function fully numerically. For illustration, these methods are applied to the Scalar Singlet model, often referred to as the simplest WIMP DM possibility from a model-building perspective. It is explicitly shown that even in this simple model the prediction for the dark matter abundance can be affected by as much as one order of magnitude.

We explore the phenomenology of having a second epoch of dark matter annihilation into dark radiation long after the standard thermal freeze-out. Such a hidden reannihilation process could affect visible sectors only gravitationally. As a concrete realization we consider self-interacting dark matter (SIDM) with a light force mediator coupled to dark radiation. We demonstrate how resonantly Sommerfeld enhanced cross sections emerge to induce the reannihilation epoch. The effect is a temporally local modification of the Hubble expansion rate and we show that the Cosmic Microwave Background (CMB) measurements -- as well as other observations -- have a high sensitivity to observe this phenomenon. Special attention is given to the model region where late kinetic decoupling and strong self-interactions can alleviate several small-scale problems in the cold dark matter paradigm at the same time. Interestingly, we find that reannihilation might here also simultaneously lower the tension between CMB and low-redshift astronomical observations of $H_{0}$ and $\sigma_{8}$. Moreover, we identify reannihilation as a clear signature to discriminate between the phenomenologically otherwise almost identical vector and scalar mediator realizations of SIDM.

Despite being one of the weakest dimers in nature, low-spectral-resolution Voyager/IRIS observations revealed the presence of (H$_2$)$_2$ dimers on Jupiter and Saturn in the 1980s. However, the collision-induced H$_2$-H$_2$ opacity databases widely used in planetary science (Borysow et al., 1985; Orton et al., 2007; Richard et al., 2012) have thus far only included free-to-free transitions and have neglected the contributions of dimers. Dimer spectra have both fine-scale structure near the S$(0)$ and S$(1)$ quadrupole lines (354 and 587 cm$^{-1}$, respectively), and broad continuum absorption contributions up to $\pm50$ cm$^{-1}$ from the line centres. We develop a new ab initio model for the free-to-bound, bound-to-free and bound-to-bound transitions of the hydrogen dimer for a range of temperatures (40-400 K) and para-hydrogen fractions (0.25-1.0). The model is validated against low-temperature laboratory experiments, and used to simulate the spectra of the giant planets. The new collision-induced opacity database permits high-resolution (0.5-1.0 cm$^{-1}$) spectral modelling of dimer spectra near S$(0)$ and S$(1)$ in both Cassini Composite Infrared Spectrometer (CIRS) observations of Jupiter and Saturn, and in Spitzer Infrared Spectrometer (IRS) observations of Uranus and Neptune for the first time. Furthermore, the model reproduces the dimer signatures observed in Voyager/IRIS data near S$(0)$ (McKellar et al., 1984) on Jupiter and Saturn, and generally lowers the amount of para-H$_2$ (and the extent of disequilibrium) required to reproduce IRIS observations.

The nature of dark matter is a longstanding enigma of physics; it may consist of particles beyond the Standard Model that are still elusive to experiments. Among indirect search techniques, which look for stable products from the annihilation or decay of dark matter particles, or from axions coupling to high-energy photons, observations of the $\gamma$-ray sky have come to prominence over the last few years, because of the excellent sensitivity of the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope mission. The LAT energy range from 20 MeV to above 300 GeV is particularly well suited for searching for products of the interactions of dark matter particles. In this report we describe methods used to search for evidence of dark matter with the LAT, and review the status of searches performed with up to six years of LAT data. We also discuss the factors that determine the sensitivities of these searches, including the magnitudes of the signals and the relevant backgrounds, considering both statistical and systematic uncertainties. We project the expected sensitivities of each search method for 10 and 15 years of LAT data taking. In particular, we find that the sensitivity of searches targeting dwarf galaxies, which provide the best limits currently, will improve faster than the square root of observing time. Current LAT limits for dwarf galaxies using six years of data reach the thermal relic level for masses up to 120 GeV for the $b\bar{b}$ annihilation channel for reasonable dark matter density profiles. With projected discoveries of additional dwarfs, these limits could extend to about 250 GeV. With as much as 15 years of LAT data these searches would be sensitive to dark matter annihilations at the thermal relic cross section for masses to greater than 400 GeV (200 GeV) in the $b\bar{b}$ ($\tau^+ \tau^-$) annihilation channels.

We present a detailed study of the annihilation signals of the inert dark matter doublet model in its high mass regime. Concretely, we study the prospects to observe gamma-ray signals of the model in current and projected Cherenkov telescopes taking into account the Sommerfeld effect and including the contribution to the spectrum from gamma-ray lines as well as from internal bremsstrahlung. We show that present observations of the galactic center by the H.E.S.S. instrument are able to exclude regions of the parameter space that give the correct dark matter relic abundance. In particular, models with the charged and the neutral components of the inert doublet nearly degenerate in mass have strong gamma-ray signals. Furthermore, for dark matter particle masses above 1 TeV, we find that the non-observation of the continuum of photons generated by the hadronization of the annihilation products typically give stronger constraints on the model parameters than the sharp spectral features associated to annihilation into monochromatic photons and the internal bremsstrahlung process. Lastly, we also analyze the interplay between indirect and direct detection searches for this model, concluding that the prospects for the former are more promising. In particular, we find that the upcoming Cherenkov Telescope Array will be able to probe a significant part of the high mass regime of the model.

We report on a search for monochromatic $\gamma$-ray features in the spectra of galaxy clusters observed by the \emph{Fermi} Large Area Telescope. Galaxy clusters are the largest structures in the Universe that are bound by dark matter (DM), making them an important testing ground for possible self-interactions or decays of the DM particles. Monochromatic $\gamma$-ray lines provide a unique signature due to the absence of astrophysical backgrounds and are as such considered a smoking-gun signature for new physics. An unbinned joint likelihood analysis of the sixteen most promising clusters using five years of data at energies between 10 and 400 GeV revealed no significant features. For the case of self-annihilation, we set upper limits on the monochromatic velocity-averaged interaction cross section. These limits are compatible with those obtained from observations of the Galactic Center, albeit weaker due to the larger distance to the studied clusters.

To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays.

In the energy range from few TeV to 25 TeV, upper bounds on the dark matter decay rate into high energy monochromatic neutrinos have recently become comparable to those on monochromatic gamma-ray lines. This implies clear possibilities of a future double "smoking-gun" evidence for the dark matter particle, from the observation of both a gamma and a neutrino line at the same energy. In particular, we show that a scenario where both lines are induced from the same dark matter particle decay leads to correlations that can already be tested. We study this "double monochromatic" scenario by considering the complete list of lowest dimensional effective operators that could induce such a decay. Furthermore, we argue that, on top of lines from decays into two-body final states, three-body final states can also be highly relevant. In addition to producing a distinct hard photon spectrum, three-body final states also produce a line-like feature in the neutrino spectrum that can be searched for by neutrino telescopes.

We show that there exists only a quite limited number of higher dimensional operators which can naturally lead to a slow decay of dark matter particles into monochromatic photons. As each of these operators inevitably induces decays into particles other than photons, we show that the gamma-lines it induces are always accompanied by a continuum flux of cosmic rays. Hence constraints on cosmic-ray fluxes imply constraints on the intensity of gamma-lines and vice-versa. A comparison with up to date observational bounds shows the possibilities to observe or exclude cosmic rays associated to gamma-line emission, so that one could better determine the properties of the DM particle, possibly discriminating between some of the operators.

Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $\gamma$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for searching for very extended emission in the vicinity of nearby Virgo galaxy cluster. Our analysis reveals statistically significant extended emission which can be well characterized by a uniformly emitting disk profile with a radius of 3\deg that moreover is offset from the cluster center. We demonstrate that the significance of this extended emission strongly depends on the adopted interstellar emission model (IEM) and is most likely an artifact of our incomplete description of the IEM in this region. We also search for and find new point source candidates in the region. We then derive conservative upper limits on the velocity-averaged DM pair annihilation cross section from Virgo. We take into account the potential $\gamma$-ray flux enhancement due to DM sub-halos and its complex morphology as a merging cluster. For DM annihilating into $b\overline{b}$, assuming a conservative sub-halo model setup, we find limits that are between 1 and 1.5 orders of magnitude above the expectation from the thermal cross section for $m_{\mathrm{DM}}\lesssim100\,\mathrm{GeV}$. In a more optimistic scenario, we exclude $\langle \sigma v \rangle\sim3\times10^{-26}\,\mathrm{cm^{3}\,s^{-1}}$ for $m_{\mathrm{DM}}\lesssim40\,\mathrm{GeV}$ for the same channel. Finally, we derive upper limits on the $\gamma$-ray-flux produced by hadronic cosmic-ray interactions in the inter cluster medium. We find that the volume-averaged cosmic-ray-to-thermal pressure ratio is less than $\sim6\%$.

We search for evidence of dark matter (DM) annihilation in the isotropic gamma-ray background (IGRB) measured with 50 months of Fermi Large Area Telescope (LAT) observations. An improved theoretical description of the cosmological DM annihilation signal, based on two complementary techniques and assuming generic weakly interacting massive particle (WIMP) properties, renders more precise predictions compared to previous work. More specifically, we estimate the cosmologically-induced gamma-ray intensity to have an uncertainty of a factor ~20 in canonical setups. We consistently include both the Galactic and extragalactic signals under the same theoretical framework, and study the impact of the former on the IGRB spectrum derivation. We find no evidence for a DM signal and we set limits on the DM-induced isotropic gamma-ray signal. Our limits are competitive for DM particle masses up to tens of TeV and, indeed, are the strongest limits derived from Fermi LAT data at TeV energies. This is possible thanks to the new Fermi LAT IGRB measurement, which now extends up to an energy of 820 GeV. We quantify uncertainties in detail and show the potential this type of search offers for testing the WIMP paradigm with a complementary and truly cosmological probe of DM particle signals.

Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.

From data recently reported from the IceCube telescope, we derive new bounds on the monochromatic neutrino signal produced from dark matter particle decays. In the few TeV to tens of TeV energy range, these bounds turn out to be better than previous limits by more than an order of magnitude. As a result, intensity constraints on neutrino lines at energies above a few TeV are now comparable to those on gamma-ray lines. From the same data sample, we also perform a detailed search for a neutrino line, showing that there is no significant hint for such a signal.

The origin of the extragalactic $\gamma$-ray background (EGB) has been debated for some time. { The EGB comprises the $\gamma$-ray emission from resolved and unresolved extragalactic sources, such as blazars, star-forming galaxies and radio galaxies, as well as radiation from truly diffuse processes.} This letter focuses on the blazar source class, the most numerous detected population, and presents an updated luminosity function and spectral energy distribution model consistent with the blazar observations performed by the {\it Fermi} Large Area Telescope (LAT). We show that blazars account for 50$^{+12}_{-11}$\,\% of the EGB photons ($>$0.1\,GeV), and that {\it Fermi}-LAT has already resolved $\sim$70\,\% of this contribution. Blazars, and in particular low-luminosity hard-spectrum nearby sources like BL Lacs, are responsible for most of the EGB emission above 100\,GeV. We find that the extragalactic background light, which attenuates blazars' high-energy emission, is responsible for the high-energy cut-off observed in the EGB spectrum. Finally, we show that blazars, star-forming galaxies and radio galaxies can naturally account for the amplitude and spectral shape of the background in the 0.1--820\,GeV range, leaving only modest room for other contributions. This allows us to set competitive constraints on the dark-matter annihilation cross section.

The {\gamma}-ray sky can be decomposed into individually detected sources, diffuse emission attributed to the interactions of Galactic cosmic rays with gas and radiation fields, and a residual all-sky emission component commonly called the isotropic diffuse {\gamma}-ray background (IGRB). The IGRB comprises all extragalactic emissions too faint or too diffuse to be resolved in a given survey, as well as any residual Galactic foregrounds that are approximately isotropic. The first IGRB measurement with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) used 10 months of sky-survey data and considered an energy range between 200 MeV and 100 GeV. Improvements in event selection and characterization of cosmic-ray backgrounds, better understanding of the diffuse Galactic emission, and a longer data accumulation of 50 months, allow for a refinement and extension of the IGRB measurement with the LAT, now covering the energy range from 100 MeV to 820 GeV. The IGRB spectrum shows a significant high-energy cutoff feature, and can be well described over nearly four decades in energy by a power law with exponential cutoff having a spectral index of $2.32\pm0.02$ and a break energy of $(279\pm52)$ GeV using our baseline diffuse Galactic emission model. The total intensity attributed to the IGRB is $(7.2\pm0.6) \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above 100 MeV, with an additional $+15$%/$-30$% systematic uncertainty due to the Galactic diffuse foregrounds.

A minimal extension of the standard model to naturally generate small neutrino masses and provide a dark matter candidate is proposed. The dark matter particle is part of a new scalar doublet field that plays a crucial role in radiatively generating neutrino masses. The symmetry that stabilizes the dark matter also suppresses neutrino masses to appear first at three-loop level. Without the need of right-handed neutrinos or other very heavy new fields, this offers an attractive explanation of the hierarchy between the electroweak and neutrino mass scales. The model has distinct verifiable predictions for the neutrino masses, flavor mixing angles, colliders and dark matter signals.

We re-evaluate the extragalactic gamma-ray flux prediction from dark matter annihilation in the approach of integrating over the nonlinear matter power spectrum, extrapolated to the free-streaming scale. We provide an estimate of the uncertainty based entirely on available N-body simulation results and minimal theoretical assumptions. We illustrate how an improvement in the simulation resolution, exemplified by the comparison between the Millennium and Millennium II simulations, affects our estimate of the flux uncertainty and we provide a "best guess" value for the flux multiplier, based on the assumption of stable clustering for the dark matter perturbations described as a collision-less fluid. We achieve results comparable to traditional Halo Model calculations, but with a much simpler procedure and a more general approach, as it relies only on one, directly measurable quantity. In addition we discuss the extension of our calculation to include baryonic effects as modeled in hydrodynamical cosmological simulations and other possible sources of uncertainty that would in turn affect indirect dark matter signals. Upper limit on the integrated power spectrum from supernovae lensing magnification are also derived and compared with theoretical expectations.

Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $\gamma$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $\gamma$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7\sigma$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $\gamma$-ray flux from individual clusters in our sample.

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $\gamma$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

The dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter-dominated objects known. Due to their proximity, high dark matter content, and lack of astrophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most promising targets for the indirect detection of dark matter via gamma rays. Here we report on gamma-ray observations of 25 Milky Way dwarf spheroidal satellite galaxies based on 4 years of Fermi Large Area Telescope (LAT) data. None of the dwarf galaxies are significantly detected in gamma rays, and we present gamma-ray flux upper limits between 500 MeV and 500 GeV. We determine the dark matter content of 18 dwarf spheroidal galaxies from stellar kinematic data and combine LAT observations of 15 dwarf galaxies to constrain the dark matter annihilation cross section. We set some of the tightest constraints to date on the the annihilation of dark matter particles with masses between 2 GeV and 10 TeV into prototypical Standard Model channels. We find these results to be robust against systematic uncertainties in the LAT instrument performance, diffuse gamma-ray background modeling, and assumed dark matter density profile.

A distinct signature for dark matter in the form of weakly interacting massive particles (WIMPs) would be the detection of a monochromatic spectral line in the gamma-ray sky. The Fermi-LAT collaboration has searched for such a line in the energy range from 5 to 300 GeV in five sky regions around the Galactic centre. No globally significant line is detected, and 95% CL upper limits on monochromatic-line strengths are presented. The smallest search region reveals a line-like structure at 133 GeV with a local significance of 2.9 sigma after 4.4 years of data, which translates to less than 1 sigma global significance from a trial factor of around 200.

We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope in the Milky Way Halo region searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e+e- produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

Recent detections of the starburst galaxies M82 and NGC 253 by gamma-ray telescopes suggest that galaxies rapidly forming massive stars are more luminous at gamma-ray energies compared to their quiescent relatives. Building upon those results, we examine a sample of 69 dwarf, spiral, and luminous and ultraluminous infrared galaxies at photon energies 0.1-100 GeV using 3 years of data collected by the Large Area Telescope (LAT) on the \textit{Fermi Gamma-ray Space Telescope} (\textit{Fermi}). Measured fluxes from significantly detected sources and flux upper limits for the remaining galaxies are used to explore the physics of cosmic rays in galaxies. We find further evidence for quasi-linear scaling relations between gamma-ray luminosity and both radio continuum luminosity and total infrared luminosity which apply both to quiescent galaxies of the Local Group and low-redshift starburst galaxies (conservative $P$-values $\lesssim0.05$ accounting for statistical and systematic uncertainties). The normalizations of these scaling relations correspond to luminosity ratios of $\log(L_{0.1-100 \rm{GeV}}/L_{1.4 \rm{GHz}}) = 1.7 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ and $\log(L_{0.1-100 \rm{GeV}}/L_{8-1000 \mu\rm{m}}) = -4.3 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ for a galaxy with a star formation rate of 1 $M_{\odot}$ yr$^{-1}$, assuming a Chabrier initial mass function. Using the relationship between infrared luminosity and gamma-ray luminosity, the collective intensity of unresolved star-forming galaxies at redshifts $0<z<2.5$ above 0.1 GeV is estimated to be 0.4-2.4 $\times 10^{-6}$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (4-23% of the intensity of the isotropic diffuse component measured with the LAT). We anticipate that $\sim10$ galaxies could be detected by their cosmic-ray induced gamma-ray emission during a 10-year \textit{Fermi} mission.

Dark matter particle annihilation or decay can produce monochromatic gamma-ray lines and contribute to the diffuse gamma-ray background. Flux upper limits are presented for gamma-ray spectral lines from 7 to 200 GeV and for the diffuse gamma-ray background from 4.8 GeV to 264 GeV obtained from two years of Fermi Large Area Telescope data integrated over most of the sky. We give cross section upper limits and decay lifetime lower limits for dark matter models that produce gamma-ray lines or contribute to the diffuse spectrum, including models proposed as explanations of the PAMELA and Fermi cosmic-ray data.

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy.

The fast outflow emerging from a region associated with massive star formation in the Orion Molecular Cloud 1 (OMC-1), located behind the Orion Nebula, appears to have been set in motion by an explosive event. Here we study the structure and dynamics of outflows in OMC-1. We combine radial velocity and proper motion data for near-IR emission of molecular hydrogen to obtain the first 3-dimensional (3D) structure of the OMC-1 outflow. Our work illustrates a new diagnostic tool for studies of star formation that will be exploited in the near future with the advent of high spatial resolution spectro-imaging in particular with data from the Atacama Large Millimeter Array (ALMA). We use published radial and proper motion velocities obtained from the shock-excited vibrational emission in the H2 v=1-0 S(1) line at 2.122 $\mu$m obtained with the GriF instrument on the Canada-France-Hawaii Telescope, the Apache Point Observatory, the Anglo-Australian Observatory and the Subaru Telescope. These data give the 3D velocity of ejecta yielding a 3D reconstruction of the outflows. This allows one to view the material from different vantage points in space giving considerable insight into the geometry. Our analysis indicates that the ejection occurred <720 years ago from a distorted ring-like structure of ~15" (6000 AU) in diameter centered on the proposed point of close encounter of the stars BN, source I and maybe also source n. We propose a simple model involving curvature of shock trajectories in magnetic fields through which the origin of the explosion and the centre defined by extrapolated proper motions of BN, I and n may be brought into spatial coincidence.

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/<I>^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

We revisit the computation of the extragalactic gamma-ray signal from cosmological dark matter annihilations. The prediction of this signal is notoriously model dependent, due to different descriptions of the clumpiness of the dark matter distribution at small scales, responsible for an enhancement with respect to the smoothly distributed case. We show how a direct computation of this "flux multiplier" in terms of the nonlinear power spectrum offers a conceptually simpler approach and may ease some problems, such as the extrapolation issue. In fact very simple analytical recipes to construct the power spectrum yield results similar to the popular Halo Model expectations, with a straightforward alternative estimate of errors. For this specific application, one also obviates to the need of identifying (often literature-dependent) concepts entering the Halo Model, to compare different simulations.

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

The single Higgs doublet in the standard model (SM) may be the simplest way of introducing electroweak symmetry breaking, but SM extensions with more scalar doublets are not excluded. A special case of the two Higgs doublet models is the inert doublet model - a minimalistic version with interesting phenomenology. These proceedings reviews the inert doublet model's theoretical setup, constraints, collider prospects and its dark matter phenomenology.

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

Dark matter (DM) particle annihilation or decay can produce monochromatic $\gamma$-rays readily distinguishable from astrophysical sources. $\gamma$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $\gamma$-ray line analysis, and integrated over most of the sky. We obtain $\gamma$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)

The inert doublet model is a minimal extension of the standard model introducing an additional SU(2) doublet with new scalar particles that could be produced at accelerators. While there exists no LEP II analysis dedicated for these inert scalars, the absence of a signal within searches for supersymmetric neutralinos can be used to constrain the inert doublet model. This translation however requires some care because of the different properties of the inert scalars and the neutralinos. We investigate what restrictions an existing DELPHI collaboration study of neutralino pair production can put on the inert scalars and discuss the result in connection with dark matter. We find that although an important part of the inert doublet model parameter space can be excluded by the LEP II data, the lightest inert particle still constitutes a valid dark matter candidate.

We report the detection of quiescent H2 emission in a spatially resolved ring-like structure within 100 AU of T Tau N. We present evidence to show that the emission most likely arises from shocks in the atmosphere of a nearly face-on disk around T Tau N. Using high spatial resolution 3D spectroscopic K-band data, we trace the spatial distribution of several H2 NIR rovibrational lines in the vicinity of T Tau N. We detect weak H2 emission from the v=1-0 S(0), S(1), Q(1) lines and the v=2-1 S(1) line in a ring-like structure around T Tau N between 0.1'' (~15 AU) and 0.7'' (~100AU) from the star. The v=1-0 S(0) and v=2-1 S(1) lines are detected only in the outer parts of the ring structure. Closer to the star, the strong continuum limits our sensitivity to these lines. The total flux of the v=1-0 S(1) line is 1.8 *10^{-14} ergs s^{-1}cm^{-2}, similar to previous measurements of H2 in circumstellar disks. The velocity of the H2 emitting gas around T Tau N is consistent with the rest velocity of the star, and the H2 does not seem to be part of a collimated outflow. Both shocks impinging on the surface of a disk and irradiation of a disk by UV-photons and X-rays from the central star are plausible candidates for the H2 excitation mechanism. However, irradiation should not create a large degree of excitation at radii larger than 20 AU. Most likely the H2 emission arises in the atmosphere of a flared disk with radius 85-100 AU and mass 0.005-0.5Msun, where the gas is excited by shocks created when a wide-angle wind impinges on the disk. The H2 emission could also originate from shock excitation in the cavity walls of an envelope, but this requires an unusually high velocity of the wide-angle wind from T Tau N.

One way to unambiguously confirm the existence of particle dark matter and determine its mass would be to detect its annihilation into monochromatic gamma-rays in upcoming telescopes. One of the most minimal models for dark matter is the inert doublet model, obtained by adding another Higgs doublet with no direct coupling to fermions. For a mass between 40 and 80 GeV, the lightest of the new inert Higgs particles can give the correct cosmic abundance of cold dark matter in agreement with current observations. We show that for this scalar dark matter candidate, the annihilation signal of monochromatic \gamma\gamma and Z\gamma final states would be exceptionally strong. The energy range and rates for these gamma-ray line signals make them ideal to search for with the soon upcoming GLAST satellite.

An extensive set of optical and NIR photometry and low-resolution spectra the Type Ia supernova (SN Ia) 2003du was obtained using a number of facilities. The observations started 13 days before B-band maximum light and continued for 480 days with exceptionally good time sampling. The optical photometry was calibrated through the S-correction technique. The UBVRIJHK light curves and the color indices of SN 2003du closely resemble those of normal SNe Ia. SN 2003du reached a B-band maximum of 13.49 (+/-0.02) mag on JD2452766.38 (+/-0.5). We derive a B-band stretch parameter of 0.988 (+/-0.003), which corresponds to dM15=1.02 (+/-0.05), indicative of a SN Ia of standard luminosity. The reddening in the host galaxy was estimated by three methods, and was consistently found to be negligible. We estimate a distance modulus of 32.79 (+/-0.15) mag to the host galaxy, UGC 9391. The peak UVOIR bolometric luminosity of 1.35(+/-0.20) 10^43 erg/s and Arnett's rule implies that M(Ni56)=0.68 (+/-0.14)M_sun of Ni56 was synthesized during the explosion. Modeling of the UVOIR bolometric light curve also indicates M(Ni56) in the range 0.6-0.8 M_sun. The spectral evolution of SN 2003du at both optical and NIR wavelengths also closely resembles normal SNe Ia. In particular, the Si II ratio at maximum R(Si II)=0.22 (+/-0.02) and the time evolution of the blueshift velocities of the absorption line minima are typical. The pre-maximum spectra of SN 2003du showed conspicuous high-velocity features in the Ca II H&K doublet and infrared triplet, and possibly in Si II 6355, lines. We compare the time evolution of the profiles of these lines with other well-observed SNe Ia and we suggest that the peculiar pre-maximum evolution of Si II 6355 line in many SNe Ia is due to the presence of two blended absorption components.

High resolution cosmological N-body simulations of four galaxy-scale dark matter halos are compared to corresponding N-body/hydrodynamical simulations containing dark matter, stars and gas. The simulations without baryons share features with others described in the literature in that the dark matter density slope continuously decreases towards the center, with a density \rho~r^{-1.3+/-0.2}, at about 1% of the virial radius for our Milky Way sized galaxies. The central cusps in the simulations which also contain baryons steepen significantly, to \rho~r^{-1.9+/-0.2}, with an indication of the inner logarithmic slope converging. Models of adiabatic contraction of dark matter halos due to the central build-up of stellar/gaseous galaxies are examined. The simplest and most commonly used model, by Blumenthal et al., is shown to overestimate the central dark matter density considerably. A modified model proposed by Gnedin et al. is shown to be a considerable improvement, but not perfect. Moreover it is found that the contraction parameters not only depend on the orbital structure of the dark-matter-only halos but also on the stellar feedback prescription which is most relevant for the baryonic distribution. Implications for dark matter annihilation at the galactic center are discussed and it is found that although our simulations show a considerable reduced halo contraction as compared to the Blumenthal et al. model, the fluxes from dark matter annihilation is still expected to be enhanced by at least a factor of a hundred as compared to dark-matter-only halos. Finally, it is shown that while dark-matter-only halos are typically prolate, the dark matter halos containing baryons are mildly oblate with minor-to-major axis ratios of c/a=0.73+/-0.11, with their flattening aligned with the central baryonic disks.

The extra-dimensional origin of dark matter is a fascinating and nowadays often discussed possibility. Here, we present the gamma-ray signatures that are expected from the self-annihilation of Kaluza-Klein dark matter particles. For comparison, we contrast this with the case of supersymmetry, where the neutralino annihilation spectra take a very different form. In both cases we find pronounced spectral signatures that could in principle be used to distinguish between these two types of dark matter candidates already with today's detector resolutions.

We investigate the internal consistency of the halo dark matter model which has been proposed by de Boer et al. to explain the excess of diffuse galactic gamma rays observed by the EGRET experiment. Any model based on dark matter annihilation into quark jets, such as the supersymmetric model proposed by de Boer et al., inevitably also predicts a primary flux of antiprotons from the same jets. Since propagation of the antiprotons in the unconventional, disk-dominated type of halo model used by de Boer et al. is strongly constrained by the measured ratio of boron to carbon nuclei in cosmic rays, we investigate the viability of the model using the DarkSUSY package to compute the gamma-ray and antiproton fluxes. We are able to show that their model is excluded by a wide margin from the measured flux of antiprotons. We therefore find that a model of the type suggested by Moskalenko et al., where the intensities of protons and electrons in the cosmic rays vary with galactic position, is far more plausible to explain the gamma excess.

Context. In order to understand the evolution of molecular clouds it is important to identify the departures from self-similarity associated with the scales of self-gravity and the driving of turbulence. Aims. A method is described based on structure functions for determining whether a region of gas, such as a molecular cloud, is fractal or contains structure with characteristic scale sizes. Methods. Using artificial data containing structure it is shown that derivatives of higher order structure functions provide a powerful way to detect the presence of characteristic scales should any be present and to estimate the size of such structures. The method is applied to observations of hot H2 in the Kleinman-Low nebula, north of the Trapezium stars in the Orion Molecular Cloud, including both brightness and velocity data. The method is compared with other techniques such as Fourier transform and histogram techniques. Results. It is found that the density structure, represented by H2 emission brightness in the K-band (2-2.5micron), exhibits mean characteristic sizes of 110, 550, 1700 and 2700AU. The velocity data show the presence of structure at 140, 1500 and 3500AU. Compared with other techniques such as Fourier transform or histogram, the method appears both more sensitive to characteristic scales and easier to interpret.

Aim: To study turbulence in the Orion Molecular Cloud (OMC1) by comparing observed and simulated characteristics of the gas motions. Method: Using a dataset of vibrationally excited H2 emission in OMC1 containing radial velocity and brightness which covers scales from 70AU to 30000AU, we present the transversal structure functions and the scaling of the structure functions with their order. These are compared with the predictions of two-dimensional projections of simulations of supersonic hydrodynamic turbulence. Results: The structure functions of OMC1 are not well represented by power laws, but show clear deviations below 2000AU. However, using the technique of extended self-similarity, power laws are recovered at scales down to 160AU. The scaling of the higher order structure functions with order deviates from the standard scaling for supersonic turbulence. This is explained as a selection effect of preferentially observing the shocked part of the gas and the scaling can be reproduced using line-of-sight integrated velocity data from subsets of supersonic turbulence simulations. These subsets select regions of strong flow convergence and high density associated with shock structure. Deviations of the structure functions in OMC1 from power laws cannot however be reproduced in simulations and remains an outstanding issue.

Observations are reported of H2 IR emission in the S(1) v=1-0 line at 2.121 microns in the Orion Molecular Cloud, OMC1, using the GriF instrument on the Canada-France-Hawaii Telescope. GriF is a combination of adaptive optics and Fabry-Perot interferometry, yielding a spatial resolution of 0.15" to 0.18" and a velocity discrimination as high as 1 km/s. Thanks to the high spatial and velocity resolution of the GriF data, 193 bright H2 emission regions can be identified in OMC1. The general characteristics of these features are described in terms of radial velocities, brightness and spatial displacement of maxima of velocity and brightness, the latter to yield the orientation of flows in the plane of the sky. Strong spatial correlation between velocity and bright H2 emission is found and serves to identify many features as shocks. Important results are: (i) velocities of the excited gas illustrate the presence of a zone to the south of BN-IRc2 and Peak 1, and the west of Peak 2, where there is a powerful blue-shifted outflow with an average velocity of -18 km/s. This is shown to be the NIR counterpart of an outflow identified in the radio from source I, a very young O-star. (ii) There is a band of weak velocity features (<5 km/s) in Peak 1 which may share a common origin through an explosive event, in the BN-IRc2 region, with the fast-moving fingers (or bullets) to the NW of OMC1. (iii) A proportion of the flows are likely to represent sites of low mass star formation and several regions show multiple outflows, probably indicative of multiple star formation within OMC1. The high spatial and velocity resolution of the GriF data show these and other features in more detail than has previously been possible.

We consider the gamma-ray spectrum from neutralino dark matter annihilations and show that internal bremsstrahlung of W pair final states gives a previously neglected source of photons at energies near the mass of the neutralino. For masses larger than about 1 TeV, and for present day detector resolutions, this results in a characteristic signal that may dominate not only over the continuous spectrum from W fragmentation, but also over the \gamma-\gamma and \gamma-Z line signals which are known to give large rates for heavy neutralinos. Observational prospects thus seem promising.

A statistical analysis is presented of the turbulent velocity structure in the Orion Molecular Cloud at scales ranging from 70 AU to 30000 AU. Results are based on IR Fabry-Perot interferometric observations of shock and photon-excited H2 in the K-band S(1) v=1-0 line at 2.121micron and refer to the dynamical characteristics of warm perturbed gas. Observations establish that the Larson size-linewidth relation is obeyed to the smallest scales studied here extending the range of validity of this relationship by nearly 2 orders of magnitude. The velocity probability distribution function (PDF) is constructed showing extended exponential wings, providing evidence of intermittency, further supported by the skewness and kurtosis of the velocity distribution. Variance and kurtosis of the PDF of velocity differences are constructed as a function of lag. The variance shows an approximate power law dependence on lag, with exponent significantly lower than the Kolmogorov value, and with deviations below 2000AU which are attributed to outflows and possibly disk structures associated with low mass star formation within OMC1. The kurtosis shows strong deviation from a gaussian velocity field, providing evidence of velocity correlations at small lags. Results agree accurately with semi-empirical simulations in Eggers & Wang (1998). In addition, 170 individual H2 emitting clumps have been analysed with sizes between 500 and 2200 AU. These show considerable diversity with regard to PDFs and variance functions. Our analysis constitutes the first characterization of the turbulent velocity field at the scale of star formation and provide a dataset which models of star-forming regions should aim to reproduce.

A TeV gamma-ray signal from the direction of the Galactic center (GC) has been detected by the H.E.S.S. experiment. Here, we investigate whether Kaluza-Klein (KK) dark matter annihilations near the GC can be the explanation. Including the contributions from internal bremsstrahlung as well as subsequent decays of quarks and tau leptons, we find a very flat gamma-ray spectrum which drops abruptly at the dark matter particle mass. For a KK mass of about 1 TeV, this gives a good fit to the H.E.S.S. data below 1 TeV. A similar model, with gauge coupling roughly three times as large and a particle mass of about 10 TeV, would give both the correct relic density and a photon spectrum that fits the complete range of data.

We investigate the fermionic one-loop cross section for the two photon annihilation of Kaluza-Klein (KK) dark matter particles in a model of universal extra dimensions (UED). This process gives a nearly mono-energetic gamma-ray line with energy equal to the KK dark matter particle mass. We find that the cross section is large enough that if a continuum signature is detected, the energy distribution of gamma-rays should end at the particle mass with a peak that is visible for an energy resolution of the detector at the percent level. This would give an unmistakable signature of a dark matter origin of the gamma-rays, and a unique determination of the dark matter particle mass, which in the case studied should be around 800 GeV. Unlike the situation for supersymmetric models where the two-gamma peak may or may not be visible depending on parameters, this feature seems to be quite robust in UED models, and should be similar in other models where annihilation into fermions is not helicity suppressed. The observability of the signal still depends on largely unknown astrophysical parameters related to the structure of the dark matter halo. If the dark matter near the galactic center is adiabatically contracted by the central star cluster, or if the dark matter halo has substructure surviving tidal effects, prospects for detection look promising.

In order not to be in conflict with observations it is crucial that extra dimensions, if they exist, are stable. It is shown that in the context of homogeneous extra dimensions, this can easily be achieved during both (4D) radiation and vacuum energy dominated eras of the cosmological evolution. During matter domination, however, there is no such possibility even for a very general class of stabilization mechanisms. Even if extra dimensions could be stabilized during matter domination, it is argued that they are generically time-varying during any transition period, such as the one from radiation to matter domination.

Near-infrared observations of line emission from excited H2 and in the continuum are reported in the direction of the Orion molecular cloud OMC1, using the European Southern Observatory Very Large Telescope UT4, equipped with the NAOS adaptive optics system and the CONICA infrared array camera. Spatial resolution has been achieved at close to the diffraction limit of the telescope (0.08" - 0.12") and images show a wealth of morphological detail. Structure is not fractal but shows two preferred scale sizes of 2.4" (1100 AU) and 1.2" (540 AU), where the larger scale may be associated with star formation.

The lightest Kaluza-Klein particle appearing in models with universal extra dimensions has recently been proposed as a viable dark matter candidate when the extra dimensions are compactified on a scale of the order of 1 TeV. Underlying assumptions of this proposal are that the size of the extra dimensions stays constant and that the evolution of the universe is given by standard cosmology. Here we investigate, both analytically and numerically, whether this is possible without introducing an explicit stabilization mechanism. By analyzing Einstein's field equations for a (3+n+1)-dimensional homogeneous, but in general anisotropic universe, we find that approximately static extra dimensions arise naturally during radiation domination. For matter domination, however, there are no solutions to the field equations that allow static extra dimensions or the usual behavior of the scale factor for ordinary three-dimensional space. We conclude that an explicit mechanism is needed in order to stabilize the extra dimensions and reproduce standard cosmology as we know it.