The $\rho$ Oph molecular cloud is one of the best examples of spinning dust emission, first detected by the Cosmic Background Imager (CBI). Here we present 4.5 arcmin observations with CBI 2 that confirm 31 GHz emission from $\rho$ Oph W, the PDR exposed to B-type star HD 147889, and highlight the absence of signal from S1, the brightest IR nebula in the complex. In order to quantify an association with dust-related emission mechanisms, we calculated correlations at different angular resolutions between the 31 GHz map and proxies for the column density of IR emitters, dust radiance and optical depth templates. We found that the 31 GHz emission correlates best with the PAH column density tracers, while the correlation with the dust radiance improves when considering emission that is more extended (from the shorter baselines), suggesting that the angular resolution of the observations affects the correlation results. A proxy for the spinning dust emissivity reveals large variations within the complex, with a dynamic range of 25 at 3$\sigma$ and a variation by a factor of at least 23, at 3$\sigma$, between the peak in $\rho$ Oph W and the location of S1, which means that environmental factors are responsible for boosting spinning dust emissivities locally.

The sightline toward the luminous blue hypergiant Cyg OB2-12 is widely used in studying interstellar dust on account of its large extinction ($A_V \simeq 10$ mag) and the fact that this extinction appears to be dominated by dust typical of the diffuse interstellar medium. We present a new analysis of archival ISO-SWS and Spitzer IRS observations of Cyg OB2-12 using a model of the emission from the star and its stellar wind to determine the total extinction $A_\lambda$ from 2.4--37 $\mu$m. In addition to the prominent 9.7 and 18 $\mu$m silicate features, we robustly detect absorption features associated with polycyclic aromatic hydrocarbons (PAHs), including the first identification of the 7.7 $\mu$m feature in absorption. The 3.3 $\mu$m aromatic feature is found to be much broader in absorption than is typically seen in emission. The 3.4 and 6.85 $\mu$m aliphatic hydrocarbon features are observed with relative strengths consistent with observation of these features on sightlines toward the Galactic Center. We identify and characterize more than sixty spectral lines in this wavelength range, which may be useful in constraining models of the star and its stellar wind. Based on this analysis, we present an extinction curve $A_\lambda/A_{2.2 \mu m}$ that extrapolates smoothly to determinations of the mean Galactic extinction curve at shorter wavelengths and to dust opacities inferred from emission at longer wavelengths, providing a new constraint on models of interstellar dust in the mid-infrared.

What is the expected infrared output of elliptical galaxies? Here we report the latest findings obtained in this high time resolution (~10 years) and high spatial resolution (2.5 parsec at center) study. We add a set of grain physics to the MACER code, including (a) dust grains made in passive stellar evolution; (b) dust grain growth due to collision and sticking; (c) grain destruction due to thermal sputtering; (d) dust cooling of hot gas via inelastic collisions; and (e) radiation pressure on dust grains. The code improvements enable us to analyze the effects of dust on metal depletion and AGN obscuration, and also to assess the infrared output of the modeled galaxies. We simulate a representative massive elliptical galaxy of a central stellar velocity dispersion ~ 260 km/s and modest rotation. We find that: (1) the circumnuclear disk (of a size <~ 1 kpc) is dusty in its outer region where most of the metals are in dust grains, while in the inner disk most of the dust grains are destroyed by the AGN irradiation; (2) the dusty disk is optically thick to both the starlight within the disk and the radiation from the central AGN. Thus the AGN will be obscured behind the disk, and the latter is of a covering factor ~ 0.2; (3) the dust infrared emission is mainly due to the AGN irradiation. The median infrared luminosity is ~ 2e44 erg/s, and it can reach >~ 1e46 erg/s during outbursts; (4) the duty cycles of the AGN activities, star formation, and the dust infrared luminosity roughly match observations, e.g., in most of its lifetime, the simulated galaxy is a stereotypical "quiescent" elliptical galaxy with L_{IR} ~ 1e11*L_{solar}, while the star formation rate can exceed 250 M_{solar}/yr during central outbursts.

Line intensity mapping (LIM) is a promising approach to study star formation and the interstellar medium (ISM) in galaxies by measuring the aggregate line emission from the entire galaxy population. In this work, we develop a simple yet physically-motivated framework for modeling the line emission as would be observed in LIM experiments, which is done by building on analytic models of the Cosmic Infrared Background that connect total infrared luminosity of galaxies to their host dark matter halos. We present models of the HI 21cm, CO(1-0), [C II] 158$\mu$m, and [N II] 122 and 205$\mu$m lines consistent with current observational constraints. With four case studies of various combinations of these lines that probe different ISM phases, we demonstrate the potential for reliably extracting physical properties of the ISM, and the evolution of these properties with cosmic time, from auto and cross-correlation analysis of these lines as measured by future LIM experiments.

Recent analyses of 21-cm neutral hydrogen (HI) emission have demonstrated that HI gas is organized into linear filamentary structures that are preferentially aligned with the local magnetic field, and that the coherence of these structures in velocity space traces line-of-sight magnetic field tangling. On this basis, we introduce a paradigm for modeling the properties of magnetized, dusty regions of the interstellar medium, using the orientation of HI structure at different velocities to map "magnetically coherent" regions of space. We construct three-dimensional (position-position-velocity) Stokes parameter maps using HI4PI full-sky spectroscopic HI data. We compare these maps, integrated over the velocity dimension, to Planck maps of the polarized dust emission at 353 GHz. Without any free parameters governing the relation between HI intensity and dust emission, we find that our $Q$ and $U$ maps are highly correlated ($r > 0.75$) with the 353 GHz $Q$ and $U$ maps of polarized dust emission observed by Planck and reproduce many of its large-scale features. The $E/B$ ratio of the dust emission maps agrees well with the HI-derived maps at large angular scales ($\ell \lesssim 120$), supporting the interpretation that this asymmetry arises from the coupling of linear density structures to the Galactic magnetic field. We demonstrate that our 3D Stokes parameter maps constrain the 3D structure of the Galactic interstellar medium and the orientation of the interstellar magnetic field.

It is well known that aligned, aspherical dust grains emit polarized radiation and that the degree of polarization depends on the angle $\psi$ between the interstellar magnetic field and the line of sight. However, anisotropy of the dust absorption cross sections also modulates the $total\ intensity$ of the radiation as the viewing geometry changes. We report a detection of this effect in the high Galactic latitude $Planck$ data, finding that the 353 GHz dust intensity per $N_{\rm HI}$ is smaller when the Galactic magnetic field is mostly in the plane of the sky and larger when the field is mostly along the line of sight. These variations are of opposite sign and roughly equal magnitude as the changes in polarized intensity per $N_{\rm HI}$ with $\psi$, as predicted. In principle, the variation in intensity can be used in conjunction with the dust polarization angle to constrain the full 3D orientation of the Galactic magnetic field.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

This white paper addresses key challenges for the design of next-decade Cosmic Microwave Background (CMB) experiments, and for assessing their capability to extract cosmological information from CMB polarization. We focus here on the challenges posed by foreground emission, CMB lensing, and instrumental systematics to detect the signal that arises from gravitational waves sourced by inflation and parameterized by $r$, at the level of $r \sim 10^{-3}$ or lower, as proposed for future observational efforts. We argue that more accurate and robust analysis and simulation tools are required for these experiments to realize their promise. We are optimistic that the capability to simulate the joint impact of foregrounds, CMB lensing, and systematics can be developed to the level necessary to support the design of a space mission at $r \sim 10^{-4}$ in a few years. We make the case here for supporting such work. Although ground-based efforts present additional challenges (e.g., atmosphere, ground pickup), which are not addressed here, they would also benefit from these improved simulation capabilities.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

Galactic dust emission is polarized at unexpectedly high levels, as revealed by Planck. The origin of the observed $\simeq 20\%$ polarization fractions can be identified by characterizing the properties of optical starlight polarization in a region with maximally polarized dust emission. We measure the R-band linear polarization of 22 stars in a region with a submillimeter polarization fraction of $\simeq 20%$. A subset of 6 stars is also measured in the B, V and I bands to investigate the wavelength dependence of polarization. We find that starlight is polarized at correspondingly high levels. Through multiband polarimetry we find that the high polarization fractions are unlikely to arise from unusual dust properties, such as enhanced grain alignment. Instead, a favorable magnetic field geometry is the most likely explanation, and is supported by observational probes of the magnetic field morphology. The observed starlight polarization exceeds the classical upper limit of $\left[p_V/E\left(B-V\right)\right]_{\rm max} = 9$%mag$^{-1}$ and is at least as high as 13%mag$^{-1}$ that was inferred from a joint analysis of Planck data, starlight polarization and reddening measurements. Thus, we confirm that the intrinsic polarizing ability of dust grains at optical wavelengths has long been underestimated.

Understanding the physics of how stars form is a highly-prioritized goal of modern Astrophysics, in part because star formation is linked to both galactic dynamics on large scales and to the formation of planets on small scales. It is well-known that stars form from the gravitational collapse of molecular clouds, which are in turn formed out of the turbulent interstellar medium. Star formation is highly inefficient, with one of the likely culprits being the regulation against gravitational collapse provided by magnetic fields. Measurement of the polarized emission from interstellar dust grains, which are partially aligned with the magnetic field, provides a key tool for understanding the role these fields play in the star formation process. Over the past decade, much progress has been made by the most recent generation of polarimeters operating over a range of wavelengths (from the far-infrared through the millimeter part of the spectrum) and over a range of angular resolutions (from less than an arcsecond through fractions of a degree). Future developments in instrument sensitivity for ground-based, airborne, and space-borne polarimeters operating over range of spatial scales are critical for enabling revolutionary steps forward in our understanding of the magnetized turbulence from which stars are formed.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

Polarized Galactic foregrounds are one of the primary sources of systematic error in measurements of the B-mode polarization of the Cosmic Microwave Background (CMB). Experiments are becoming increasingly sensitive to complexities in the foreground frequency spectra that are not captured by standard parametric models, potentially affecting our ability to efficiently separate out these components. Employing a suite of dust models encompassing a variety of physical effects, we simulate observations of a future seven-band CMB experiment to assess the impact of these complexities on parametric component separation. We identify configurations of frequency bands that minimize the `model errors' caused by fitting simple parametric models to more complex `true' foreground spectra, which bias the inferred CMB signal. We find that: (a) fits employing a simple two parameter modified blackbody (MBB) dust model tend to produce significant bias in the recovered polarized CMB signal in the presence of physically realistic dust foregrounds; (b) generalized MBB models with three additional parameters reduce this bias in most cases, but non-negligible biases can remain, and can be hard to detect; and (c) line of sight effects, which give rise to frequency decorrelation, and the presence of iron grains are the most problematic complexities in the dust emission for recovering the true CMB signal. More sophisticated simulations will be needed to demonstrate that future CMB experiments can successfully mitigate these more physically realistic dust foregrounds.

PASIPHAE (the Polar-Areas Stellar Imaging in Polarization High-Accuracy Experiment) is an optopolarimetric survey aiming to measure the linear polarization from millions of stars, and use these to create a three-dimensional tomographic map of the magnetic field threading dust clouds within the Milky Way. This map will provide invaluable information for future CMB B-mode experiments searching for inflationary gravitational waves, providing unique information regarding line-of-sight integration effects. Optical polarization observations of a large number of stars at known distances, tracing the same dust that emits polarized microwaves, can map the magnetic field between them. The Gaia mission is measuring distances to a billion stars, providing an opportunity to produce a tomographic map of Galactic magnetic field directions, using optical polarization of starlight. Such a map will not only boost CMB polarization foreground removal, but it will also have a profound impact in a wide range of astrophysical research, including interstellar medium physics, high-energy astrophysics, and evolution of the Galaxy. Taking advantage of the novel technology implemented in our high-accuracy Wide-Area Linear Optical Polarimeters (WALOPs) currently under construction at IUCAA, India, we will engage in a large-scale optopolarimetric program that can meet this challenge: a survey of both northern and southern Galactic polar regions targeted by CMB experiments, covering over 10,000 square degrees, which will measure linear optical polarization of over 360 stars per square degree (over 3.5 million stars, a 1000-fold increase over the state of the art). The survey will be conducted concurrently from the South African Astronomical Observatory in Sutherland, South Africa in the southern hemisphere, and the Skinakas Observatory in Crete, Greece, in the north.

We discuss the nature of a discrete, compact radio source (NGC 4725 B) located $\approx$1.9 kpc from the nucleus in the nearby star-forming galaxy NGC 4725, which we believe to be a new detection of extragalactic Anomalous Microwave Emission (AME). Based on detections at 3, 15, 22, 33, and 44 GHz, NGC 4725 B is a $\mu$Jy radio source peaking at $\approx$33 GHz. While the source is not identified in $BVRI$ photometry, we detect counterparts in the mid-infrared $Spitzer$/IRAC bands (3.6, 4.5, 5.8, 8.0 $\mu$m) that appear to be associated with dust emission in the central region of NGC 4725. Consequently, we conclude that NGC 4725 B is a new detection of AME, and very similar to a recent detection of AME in an outer-disk star-forming region in NGC 6946. We find that models of electric dipole emission from rapidly rotating ultra-small grains are able to reproduce the radio spectrum for reasonable interstellar medium conditions. Given the lack of an optical counterpart and the shape of the radio spectrum, NGC 4725 B appears consistent with a nascent star-forming region in which young ($\lesssim 3$ Myr) massive stars are still highly enshrouded by their natal cocoons of gas and dust with insufficient supernovae occurring to produce a measurable amount of synchrotron emission.

Anomalous Microwave Emission (AME) is a component of diffuse Galactic radiation observed at frequencies in the range $\approx 10$-60 GHz. AME was first detected in 1996 and recognised as an additional component of emission in 1997. Since then, AME has been observed by a range of experiments and in a variety of environments. AME is spatially correlated with far-IR thermal dust emission but cannot be explained by synchrotron or free-free emission mechanisms, and is far in excess of the emission contributed by thermal dust emission with the power-law opacity consistent with the observed emission at sub-mm wavelengths. Polarization observations have shown that AME is very weakly polarized ($\lesssim 1$%). The most natural explanation for AME is rotational emission from ultra-small dust grains ("spinning dust"), first postulated in 1957. Magnetic dipole radiation from thermal fluctuations in the magnetization of magnetic grain materials may also be contributing to the AME, particularly at higher frequencies ($\gtrsim 50$ GHz). AME is also an important foreground for Cosmic Microwave Background analyses. This paper presents a review and the current state-of-play in AME research, which was discussed in an AME workshop held at ESTEC, The Netherlands, June 2016.

Although interstellar grains are known to be aspherical, their actual shapes remain poorly constrained. We assess whether three distributions of ellipsoidal shapes from the literature, employed principally for their convenient mathematical properties, are suitable for describing the shapes of interstellar grains. Randomly-selected shapes from each distribution are shown as illustrations. The often-used BHCDE distribution includes a very large fraction of extreme shapes: fully 10% of random draws have axial ratio a_3/a_1 > 19.7, and 5% have a_3/a_1 > 33. Of the three distributions considered, the CDE2 appears to be most realistic. For each of the three CDEs considered, we derive shape-averaged cross sections for extinction and polarization. Finally, we describe a method for "synthesizing" a dielectric function for an assumed shape or shape distribution if the actual absorption cross sections per grain volume are known from observations. This synthetic dielectric function predicts the wavelength dependence of polarization, which can then be compared to observations to constrain the grain shape.

We present a new map of interstellar reddening, covering the 39\% of the sky with low {\rm HI} column densities ($N_{\rm HI} < 4\times10^{20}\,\rm cm^{-2}$ or $E(B-V)\approx 45\rm\, mmag$) at $16\overset{'}{.}1$ resolution, based on all-sky observations of Galactic HI emission by the HI4PI Survey. In this low column density regime, we derive a characteristic value of $N_{\rm HI}/E(B-V) = 8.8\times10^{21}\, \rm\, cm^{2}\, mag^{-1}$ for gas with $|v_{\rm LSR}| < 90\,\rm km\, s^{-1}$ and find no significant reddening associated with gas at higher velocities. We compare our HI-based reddening map with the Schlegel, Finkbeiner, and Davis (1998, SFD) reddening map and find them consistent to within a scatter of $\simeq 5\,\rm mmag$. Further, the differences between our map and the SFD map are in excellent agreement with the low resolution ($4\overset{\circ}{.}5$) corrections to the SFD map derived by Peek and Graves (2010) based on observed reddening toward passive galaxies. We therefore argue that our HI-based map provides the most accurate interstellar reddening estimates in the low column density regime to date. Our reddening map is made publicly available (http://dx.doi.org/10.7910/DVN/AFJNWJ).

Interstellar iron in the form of metallic iron nanoparticles may constitute a component of the interstellar dust. We compute the stability of iron nanoparticles to sublimation in the interstellar radiation field, finding that iron clusters can persist down to a radius of $\simeq 4.5\,$\AA, and perhaps smaller. We employ laboratory data on small iron clusters to compute the photoelectric yields as a function of grain size and the resulting grain charge distribution in various interstellar environments, finding that iron nanoparticles can acquire negative charges particularly in regions with high gas temperatures and ionization fractions. If $\gtrsim 10\%$ of the interstellar iron is in the form of ultrasmall iron clusters, the photoelectric heating rate from dust may be increased by up to tens of percent relative to dust models with only carbonaceous and silicate grains.

In light of recent observational results indicating an apparent lack of correlation between the Anomalous Microwave Emission (AME) and mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs), we assess whether rotational emission from spinning silicate and/or iron nanoparticles could account for the observed AME without violating observational constraints on interstellar abundances, ultraviolet extinction, and infrared emission. By modifying the SpDust code to compute the rotational emission from these grains, we find that nanosilicate grains could account for the entirety of the observed AME, whereas iron grains could be responsible for only a fraction, even for extreme assumptions on the amount of interstellar iron concentrated in ultrasmall iron nanoparticles. Given the added complexity of contributions from multiple grain populations to the total spinning dust emission, as well as existing uncertainties due to the poorly-constrained grain size, charge, and dipole moment distributions, we discuss generic, carrier-independent predictions of spinning dust theory and observational tests that could help identify the AME carrier(s).

The quantization of energy levels in very nanoparticles suppresses dissipative processes that convert grain rotational kinetic energy into heat. For grains small enough to have GHz rotation rates, the suppression of dissipation can be extreme. As a result, alignment of such grains is suppressed. This applies both to alignment of the grain body with its angular momentum J, and to alignment of J with the local magnetic field B_0. If the anomalous microwave emission is rotational emission from spinning grains, it will be negligibly polarized at GHz frequencies, with P < 10^{-6} at frequencies above 10 GHz.

We employ an all-sky map of the anomalous microwave emission (AME) produced by component separation of the microwave sky to study correlations between the AME and Galactic dust properties. We find that while the AME is highly correlated with all tracers of dust emission, the best predictor of the AME strength is the dust radiance. Fluctuations in the AME intensity per dust radiance are uncorrelated with fluctuations in the emission from polycyclic aromatic hydrocarbons (PAHs), casting doubt on the association between AME and PAHs. The PAH abundance is strongly correlated with the dust optical depth and dust radiance, consistent with PAH destruction in low density regions. We find that the AME intensity increases with increasing radiation field strength, at variance with predictions from the spinning dust hypothesis. Finally, the temperature-dependence of the AME per dust radiance disfavors the interpretation of the AME as thermal emission. A reconsideration of other AME carriers, such as ultrasmall silicates, and other emission mechanisms, such as magnetic dipole emission, is warranted.

Using 1 cm and 3 mm CARMA and 2 mm GISMO observations, we follow up the first extragalactic detection of anomalous microwave emission (AME) reported by Murphy et al. 2010 in an extranuclear region (Enuc. 4) of the nearby face-on spiral galaxy NGC 6946. We find the spectral shape and peak frequency of AME in this region to be consistent with models of spinning dust emission. However, the strength of the emission far exceeds the Galactic AME emissivity given the abundance of polycyclic aromatic hydrocarbons (PAHs) in that region. Using our galaxy-wide 1 cm map (21" resolution), we identify a total of eight 21"x21" regions in NGC 6946 that harbour AME at >95% significance at levels comparable to that observed in Enuc. 4. The remainder of the galaxy has 1 cm emission consistent with or below the observed Galactic AME emissivity per PAH surface density. We probe relationships between the detected AME and dust surface density, PAH emission, and radiation field, though no environmental property emerges to delineate regions with strong versus weak or non-existent AME. On the basis of these data and other AME observations in the literature, we determine that the AME emissivity per unit dust mass is highly variable. We argue that the spinning dust hypothesis, which predicts the AME power to be approximately proportional to the PAH mass, is therefore incomplete.

Photometric observations of galaxies at submillimetre to millimetre wavelengths (50 - 1000 GHz) are susceptible to spatial variations in both the background CMB temperature and CIB emission that can be comparable to the flux from the target galaxy. We quantify the residual uncertainty when background emission inside a circular aperture is estimated by the mean flux in a surrounding annular region, assumed to have no contribution from the source of interest. We present simple formulae to calculate this uncertainty as a function of wavelength and aperture size. Drawing on examples from the literature, we illustrate the use of our formalism in practice and highlight cases in which uncertainty in the background subtraction needs to be considered in the error analysis. We make the code used to calculate the uncertainties publicly available on the web.

We study the effects of a detailed dust treatment on the properties and evolution of early-type galaxies containing central black holes, as determined by AGN feedback. We find that during cooling flow episodes, radiation pressure on the dust in and interior to infalling shells of cold gas can greatly impact the amount of gas able to be accreted and therefore the frequency of AGN bursts. However, the overall hydrodynamic evolution of all models, including mass budget, is relatively robust to the assumptions on dust. We find that IR re-emission from hot dust can dominate the bolometric luminosity of the galaxy during the early stages of an AGN burst, reaching values in excess of $10^{46}$ erg/s. The AGN-emitted UV is largely absorbed, but the optical depth in the IR does not exceed unity, so the radiation momentum input never exceeds $L_{\rm BH}/c$. We constrain the viability of our models by comparing the AGN duty cycle, broadband luminosities, dust mass, black hole mass, and other model predictions to current observations. These constraints force us to models wherein the dust to metals ratios are $\simeq 1%$ of the Galactic value, and only models with a dynamic dust to gas ratio are able to produce both quiescent galaxies consistent with observations and high obscured fractions during AGN "on" phases. During AGN outbursts, we predict that a large fraction of the FIR luminosity can be attributed to warm dust emission ($\simeq100$ K) from dense dusty gas within $\leq 1$ kpc reradiating the AGN UV emission.

The total anisotropy of a diffuse background composed of two or more sources, such as the Fermi-LAT--measured gamma-ray background, is set by the anisotropy of each source population and the contribution of each population to the total intensity. The total anisotropy as a function of energy (the anisotropy energy spectrum) will modulate as the relative contributions of the sources change, implying that the anisotropy energy spectrum also encodes the intensity spectrum of each source class. We develop techniques, applicable to any such diffuse background, for unraveling the intensity spectrum of each component source population given a measurement of the total intensity spectrum and the total anisotropy energy spectrum, without introducing \emph{a priori} assumptions about the spectra of the source classes. We demonstrate the potential of these methods by applying them to example scenarios for the composition of the Fermi-LAT gamma-ray background consistent with current data and feasible within 10 years of observation.

The presence of ferromagnetic or ferrimagnetic nanoparticles in the interstellar medium would give rise to magnetic dipole radiation at microwave and submm frequencies. Such grains may account for the strong mm-wavelength emission observed from a number of low-metallicity galaxies, including the Small Magellanic Cloud. We show how to calculate the absorption and scattering cross sections for such grains, with particular attention to metallic Fe, magnetite Fe3O4, and maghemite gamma-Fe2O3, all potentially present in the interstellar medium. The rate of Davis-Greenstein alignment by magnetic dissipation is also estimated. We determine the temperature of free-flying magnetic grains heated by starlight and we calculate the polarization of the magnetic dipole emission from both free-fliers and inclusions. For inclusions, the magnetic dipole emission is expected to be polarized orthogonally relative to the normal electric dipole radiation. Finally, we present self-consistent dielectric functions for metallic Fe, magnetite Fe3O4, and maghemite gamma-Fe2O3, enabling calculation of absorption and scattering cross sections from microwave to X-ray wavelengths.

The Small Magellanic Cloud (SMC) has surprisingly strong submm and mm-wavelength emission that is inconsistent with standard dust models, including those with emission from spinning dust. Here we show that the emission from the SMC may be understood if the interstellar dust mixture includes magnetic nanoparticles, emitting magnetic dipole radiation resulting from thermal fluctuations in the magnetization. The magnetic grains can be metallic iron, magnetite Fe3O4, or maghemite gamma-Fe2O3. The required mass of iron is consistent with elemental abundance constraints. The magnetic dipole emission is predicted to be polarized orthogonally to the normal electric dipole radiation if the nanoparticles are inclusions in larger grains. We speculate that other low-metallicity galaxies may also have a large fraction of the interstellar Fe in magnetic materials.

The energy-dependence of the anisotropy (the anisotropy energy spectrum) of the large-scale diffuse gamma-ray background can reveal the presence of multiple source populations. Annihilating dark matter in the substructure of the Milky Way halo could give rise to a modulation in the anisotropy energy spectrum of the diffuse gamma-ray emission measured by Fermi, enabling the detection of a dark matter signal. We determine the detectability of a dark-matter-induced modulation for scenarios in which unresolved blazars are the primary contributor to the measured emission above ~1 GeV and find that in some scenarios pair-annihilation cross sections of order the value expected for thermal relic dark matter can produce a detectable feature. We anticipate that the sensitivity of this technique to specific dark matter models could be improved by tailored likelihood analysis methods.

Exploration of time domain is now a vibrant area of research in astronomy, driven by the advent of digital synoptic sky surveys. While panoramic surveys can detect variable or transient events, typically some follow-up observations are needed; for short-lived phenomena, a rapid response is essential. Ability to automatically classify and prioritize transient events for follow-up studies becomes critical as the data rates increase. We have been developing such methods using the data streams from the Palomar-Quest survey, the Catalina Sky Survey and others, using the VOEventNet framework. The goal is to automatically classify transient events, using the new measurements, combined with archival data (previous and multi-wavelength measurements), and contextual information (e.g., Galactic or ecliptic latitude, presence of a possible host galaxy nearby, etc.); and to iterate them dynamically as the follow-up data come in (e.g., light curves or colors). We have been investigating Bayesian methodologies for classification, as well as discriminated follow-up to optimize the use of available resources, including Naive Bayesian approach, and the non-parametric Gaussian process regression. We will also be deploying variants of the traditional machine learning techniques such as Neural Nets and Support Vector Machines on datasets of reliably classified transients as they build up.