We study the projected spatial offset between the ultraviolet continuum and Ly$\alpha$ emission for a sample of 65 lensed and unlensed galaxies in the Epoch of Reionization (5<z<7), the first such study at these redshifts, in order to understand the potential for these offsets to confuse estimates of the Ly$\alpha$ properties of galaxies observed in slit spectroscopy. While we find that ~40% of galaxies in our sample show significant projected spatial offsets ($|\Delta_{Ly\alpha-UV}|$), we find a modest average offset of 0.61$\pm$0.08 kpc. A small fraction of our sample, ~10%, exhibits offsets of 2-4 kpc, sizes that are larger than the effective radii of typical galaxies at these redshifts. An internal comparison and a comparison to studies at lower redshift yielded no significant evidence of evolution of $|\Delta_{Ly\alpha-UV}|$ with redshift. In our own sample, UV-bright galaxies showed offsets a factor of three greater than their fainter counterparts, 0.89$\pm$0.18 vs. 0.27$\pm$0.05 kpc, respectively. We argue that offsets are likely not the result of merging processes, but are rather due to anisotropic internal processes internal to the galaxies as a result of stellar feedback, which facilitates Ly$\alpha$ fluorescence and/or backscattering from nearby or outflowing gas. The reduction in the Ly$\alpha$ flux due to offset effects for various observational setups was quantified through mock observations of simple simulations. It was found that the loss of Ly$\alpha$ photons for galaxies with average offsets is not, if corrected for, a limiting factor for all but the narrowest slit widths (<0.4''). However, for the largest offsets, if such offsets are mostly perpendicular to the slit major axis, slit losses were found to be extremely severe in cases where slit widths of $\le$1'' were employed, such as those planned for James Webb Space Telescope/NIRSpec observations. (abridged)

We present results from the KMOS Lens-Amplified Spectroscopic Survey (KLASS), an ESO Very Large Telescope (VLT) large program using gravitational lensing to study the spatially resolved kinematics of 44 star-forming galaxies at 0.6<z<2.3 with a stellar mass of 8.1<log(M$_\star$/M$_{\odot}$)<11.0. These galaxies are located behind six galaxy clusters selected from the HST Grism Lens-Amplified Survey from Space (GLASS). We find that the majority of the galaxies show a rotating disk, but most of the rotation-dominated galaxies only have a low $\upsilon_{rot}/\sigma_0$ ratio (median of $\upsilon_{rot}/\sigma_0\sim2.5$). We explore the Tully-Fisher relation by adopting the circular velocity, $V_{circ}=(\upsilon_{rot}^2+3.4\sigma_0^2)^{1/2}$, to account for pressure support. We find that our sample follows a Tully-Fisher relation with a positive zero-point offset of +0.18 dex compared to the local relation, consistent with more gas-rich galaxies that still have to convert most of their gas into stars. We find a strong correlation between the velocity dispersion and stellar mass in the KLASS sample. When combining our data to other surveys from the literature, we also see an increase of the velocity dispersion with stellar mass at all redshift. We obtain an increase of $\upsilon_{rot}/\sigma_0$ with stellar mass at 0.5<z<1.0. This could indicate that massive galaxies settle into regular rotating disks before the low-mass galaxies. For higher redshift (z>1), we find a weak increase or flat trend. We investigate the relation between the rest-frame UV clumpiness of galaxies and their global kinematic properties. We find no clear trend between the clumpiness and the velocity dispersion and $\upsilon_{rot}/\sigma_0$. This could suggest that the kinematic properties of galaxies evolve after the clumps formed in the galaxy disk or that the clumps can form in different physical conditions.

Recent zoom-in cosmological simulations have shown that stellar feedback can flatten the inner density profile of the dark matter halo in low-mass galaxies. A correlation between the stellar/gas velocity dispersion ($\sigma_{star}$, $\sigma_{gas}$) and the specific star formation rate (sSFR) is predicted as an observational test of the role of stellar feedback in re-shaping the dark matter density profile. In this work we test the validity of this prediction by studying a sample of star-forming galaxies at $0.6<z<1.0$ from the LEGA-C survey, which provides high signal-to-noise measurements of stellar and gas kinematics. We find that a weak but significant correlation between $\sigma_{star}$ (and $\sigma_{gas}$) and sSFR indeed exists for galaxies in the lowest mass bin (M$_\ast\sim10^{10}\,$M$_\odot$). This correlation, albeit with a $\sim$35% scatter, holds for different tracers of star formation, and becomes stronger with redshift. This result generally agrees with the picture that at higher redshifts star formation rate was generally higher, and galaxies at M$_\ast\lesssim10^{10}\,$M$_\odot$ have not yet settled into a disk. As a consequence, they have shallower gravitational potentials more easily perturbed by stellar feedback. The observed correlation between $\sigma_{star}$ (and $\sigma_{gas}$) and sSFR supports the scenario predicted by cosmological simulations, in which feedback-driven outflows cause fluctuations in the gravitation potential that flatten the density profiles of low-mass galaxies.

Physical properties of galaxies at z>7 are of interest for understanding both the early phases of star formation and the process of cosmic reionization. Chemical abundance measurements offer valuable information on the integrated star formation history, and hence ionizing photon production, as well as the rapid gas accretion expected at such high redshifts. We use reported measurements of [O III] 88$\mu$m emission and star formation rate to estimate gas-phase oxygen abundances in five galaxies at z=7.1-9.1 using the direct T_e method. We find typical abundances 12+log(O/H) = 7.9 ($\sim$0.2 times the solar value) and an evolution of 0.9$\pm$0.5 dex in oxygen abundance at fixed stellar mass from z$\simeq$8 to 0. These results are compatible with theoretical predictions, albeit with large (conservative) uncertainties in both mass and metallicity. We assess both statistical and systematic uncertainties to identify promising means of improvement with the Atacama Large Millimeter Array (ALMA) and the James Webb Space Telescope (JWST). In particular we highlight [O III] 52$\mu$m as a valuable feature for robust metallicity measurements. Precision of 0.1-0.2 dex in T_e-based O/H abundance can be reasonably achieved for galaxies at z$\approx$5-8 by combining [O III] 52$\mu$m with rest-frame optical strong lines. It will also be possible to probe gas mixing and mergers via resolved T_e-based abundances on kpc scales. With ALMA and JWST, direct metallicity measurements will thus be remarkably accessible in the reionization epoch.

Despite significant progress both observationally and theoretically, the origin of high-ionization nebular HeII emission in galaxies dominated by stellar photoionization remains unclear. Accretion-powered radiation from high-mass X-ray binaries (HMXBs) is still one of the leading proposed explanations for the missing $\mathrm{He^+}$-ionizing photons, but this scenario has yet to be conclusively tested. In this paper, we present nebular line predictions from a grid of photoionization models with input SEDs containing the joint contribution of both stellar atmospheres and a multi-color disk model for HMXBs. This grid demonstrates that HMXBs are inefficient producers of the photons necessary to power HeII, and can only boost this line substantially in galaxies with HMXB populations large enough to power X-ray luminosities of $10^{42}$ erg/s per unit star formation rate (SFR). To test this, we assemble a sample of eleven low-redshift star-forming galaxies with high-quality constraints on both X-ray emission from Chandra and HeII emission from deep optical spectra, including new observations with the MMT. These data reveal that the HMXB populations of these nearby systems are insufficient to account for the observed HeII strengths, with typical X-ray luminosities or upper limits thereon of only $10^{40}$-$10^{41}$ erg/s per SFR. This indicates that HMXBs are not the dominant source of $\mathrm{He^+}$ ionization in these metal-poor star-forming galaxies. We suggest that the solution may instead reside in revisions to stellar wind predictions, softer X-ray sources, or very hot products of binary evolution at low metallicity.

Galaxy evolution is thought to be driven in large part by the flow of gas between galaxies and the circumgalactic medium (CGM), a halo of metal-enriched gas extending out to $\gtrsim100$ kpc from each galaxy. Studying the spatial structure of the CGM holds promise for understanding these gas flow mechanisms; however, the common method using background quasar sightlines provides minimal spatial information. Recent works have shown the utility of extended background sources such as giant gravitationally lensed arcs. Using background lensed arcs from the CSWA 38 lens system, we continuously probed, at a resolution element of about 15 kpc$^2$, the spatial and kinematic distribution of MgII absorption in a star-forming galaxy at $z=0.77$ (stellar mass $\approx 10^{9.7}$ M$_\odot$, star formation rate $\approx 10$ M$_\odot$ yr$^{-1}$) at impact parameters $D=5-40$ kpc. Our results present an anisotropic, optically thick medium whose absorption strength decreases with increasing impact parameter, in agreement with the statistics towards quasars and other gravitational arcs. Furthermore, we detect low line-of-sight velocities ($v\approx-25-80$ km s$^{-1}$) and relatively high velocity dispersion ($\sigma\approx50\pm30$ km s$^{-1}$) in the MgII gas. These measures provide evidence of a mainly pressure-supported, metal-enriched diffuse gas recycling through the CGM rather than an active galactic outflow.

We present the results of deep Chandra and XMM-Newton observations of a complex merging galaxy cluster Abell 2256 (A2256) that hosts a spectacular radio relic. The temperature and metallicity maps show clear evidence of a merger between the western subcluster and the primary cluster. We detect five X-ray surface brightness edges. Three of them near the cluster center are cold fronts (CFs): CF1 is associated with the infalling subcluster; CF2 is located in the east of the primary cluster and CF3 is to the west of the primary cluster core. The other two edges at cluster outskirts are shock fronts (SFs): SF1 near the radio relic in the NW has Mach numbers derived from the temperature and the density jumps respectively of $M_T=1.62\pm0.12$ and $M_\rho=1.23\pm0.06$; SF2 in the SE has $M_T=1.54\pm0.05$ and $M_\rho=1.16\pm0.13$. In the region of the radio relic there is no evidence for the correlation between X-ray and radio substructures, from which we estimate an upper limit for the inverse Compton emission, and therefore set a lower limit on the magnetic field ($\sim$ 450 kpc from primary cluster center) of $B>1.0\ \mu$G for a single power-law electron spectrum or $B>0.4\ \mu$G for a broken power-law electron spectrum. We propose a merger scenario including a primary cluster, a subcluster, and a group. Our merger scenario accounts for the X-ray edges, diffuse radio features, and galaxy kinematics, as well as projection effects.

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

We present the status of on-going detector development efforts for our joint NASA/CNES balloon-borne UV multi-object spectrograph, the Faint Intergalactic Redshifted Emission Balloon (FIREBall-2; FB-2). FB-2 demonstrates a new UV detector technology, the delta-doped Electron Multiplying CCD (EMCCD), in a low risk suborbital environment, to prove the performance of EMCCDs for future space missions and Technology Readiness Level (TRL) advancement. EMCCDs can be used in photon counting (PC) mode to achieve extremely low readout noise ($<$1 electron). Our testing has focused on reducing clock-induced-charge (CIC) through wave shaping and well depth optimization with a \nuvu V2 CCCP Controller, measuring CIC at 0.001 e$^{-}$/pixel/frame. This optimization also includes methods for reducing dark current, via cooling, and substrate voltage levels. We discuss the challenges of removing cosmic rays, which are also amplified by these detectors, as well as a data reduction pipeline designed for our noise measurement objectives. FB-2 flew in 2018, providing the first time an EMCCD was used for UV observations in the stratosphere. FB-2 is currently being built up to fly again in 2020, and improvements are being made to the EMCCD to continue optimizing its performance for better noise control.

The degree of turbulent pressure support by residual gas motions in galaxy clusters is not well known. Mass modelling of combined X-ray and Sunyaev Zel'dovich observations provides an estimate of turbulent pressure support in the outer regions of several galaxy clusters. Here, we test two different filtering techniques to disentangle bulk from turbulent motions in non-radiative high-resolution cosmological simulations of galaxy clusters using the cosmological hydro code ENZO. We find that the radial behavior of the ratio of non-thermal pressure to total gas pressure as a function of cluster-centric distance can be described by a simple polynomial function. The typical non-thermal pressure support in the centre of clusters is $\sim$5%, increasing to $\sim$15% in the outskirts, in line with the pressure excess found in recent X-ray observations. While the complex dynamics of the ICM makes it impossible to reconstruct a simple correlation between turbulent motions and hydrostatic bias, we find that a relation between them can be established using the median properties of a sample of objects. Moreover, we estimate the contribution of radial accelerations to the non-thermal pressure support and conclude that it decreases moving outwards from 40% (in the core) to 15% (in the cluster's outskirts). Adding this contribution to one provided by turbulence, we show that it might account for the entire observed hydrostatic bias in the innermost regions of the clusters, and for less than 80% of it at $r > 0.8 r_{200, m}$.

We present high spatial resolution LBTI/NOMIC $9-12$ $\mu m$ images of VY CMa and its massive outflow feature, the Southwest (SW) Clump. Combined with high-resolution imaging from HST ($0.4-1$ $\mu m$) and LBT/LMIRCam ($1-5$ $\mu m$), we isolate the spectral energy distribution (SED) of the clump from the star itself. Using radiative-transfer code DUSTY, we model both the scattered light from VY CMa and the thermal emission from the dust in the clump to estimate the optical depth, mass, and temperature of the SW Clump. The SW Clump is optically thick at 8.9 $\mu m$ with a brightness temperature of $\sim$200 K. With a dust chemistry of equal parts silicates and metallic iron, as well as assumptions on grain size distribution, we estimate a dust mass of $5.4\times10^{-5}\,M_\odot$. For a gas--to--dust ratio of 100, this implies a total mass of $5.4\times10^{-3}\,M_\odot$. Compared to the typical mass-loss rate of VY CMa, the SW Clump represents an extreme, localized mass-loss event from $\lesssim300$ years ago.

Magnetic fields play vital roles in intracluster media (ICMs), but estimating their strengths and distributions from observations is a major challenge. Faraday rotation measures (RMs) are widely applied to this task, so it is critical to understand inherent uncertainties in RM analysis. In this paper, we seek to characterize those uncertainties given the types of information available today, independent of the specific technique used. We conduct synthetic RM observations through the ICM of a galaxy cluster drawn from an MHD cosmological simulation in which the magnetic field is known. We analyze the synthetic RM observations using an analytical formalism based on commonly used model assumptions allowing us to relate model physical variables to outcome uncertainties. Despite the simplicity of some assumptions, and unknown physical parameters, we are able to extract an approximate magnitude of the central magnetic field within an apparently irreducible uncertain factor approximately 3. Principal, largely irreducible, uncertainties come from the unknown depth along the line of sight of embedded polarized sources, the lack of robust coherence lengths from area-constrained polarization sampling, and the unknown scaling between ICM electron density and magnetic field strength. The RM-estimated central magnetic field strengths span more than an order of magnitude including the full range of synthetic experiments.

We present measurements of [SIII]$\lambda\lambda$9069,9531 for a sample of $z\sim1.5$ star-forming galaxies, the first sample with measurements of these lines at z>0.1. We employ the line ratio S$_{32}$$\equiv$[SIII]$\lambda\lambda$9069,9531/[SII]$\lambda\lambda$6716,6731 as a novel probe of evolving ISM conditions. Since this ratio includes the low-ionization line [SII], it is crucial that the effects of diffuse ionized gas (DIG) on emission-line ratios be accounted for in $z\sim0$ integrated galaxy spectra, or else that comparisons be made to samples of local HII regions in which DIG emission is not present. We find that S$_{32}$ decreases with increasing stellar mass at both $z\sim1.5$ and $z\sim0$, but that the dependence is weak suggesting S$_{32}$ has a very shallow anticorrelation with metallicity, in contrast with O$_{32}$ that displays a strong metallicity dependence. As a result, S$_{32}$ only mildly evolves with redshift at fixed stellar mass. The $z\sim1.5$ sample is systematicallty offset towards lower S$_{32}$ and higher [SII]/H$\alpha$ at fixed [OIII]/H$\beta$ relative to $z=0$ HII regions. By comparing to photoionization model grids, we find that such trends can be explained by a scenario in which the ionizing spectrum is harder at fixed O/H with increasing redshift, but are inconsistent with an increase in ionization parameter at fixed O/H. This analysis demonstrates the advantages of expanding beyond the strongest rest-optical lines for evolutionary studies, and the particular utility of [SIII] for characterizing evolving ISM conditions and stellar compositions. These measurements provide a basis for estimating [SIII] line strengths for high-redshift galaxies, a line that the James Webb Space Telescope will measure out to z~5.5.

We present hitherto the largest sample of gas-phase metallicity radial gradients measured at sub-kiloparsec resolution in star-forming galaxies in the redshift range of $z\in[1.2, 2.3]$. These measurements are enabled by the synergy of slitless spectroscopy from the Hubble Space Telescope near-infrared channels and the lensing magnification from foreground galaxy clusters. Our sample consists of 79 galaxies with stellar mass ranging from 10$^7$ to 10$^{10}$ $M_\odot$, instantaneous star-formation rate in the range of [1, 100] $M_\odot$/yr, and global metallicity [$\frac{1}{12}$, 2] solar. At 2-$\sigma$ confidence level, 14/79 galaxies in our sample show negative radial gradients, whereas 7/79 show inverted gradients. Combining ours and all other metallicity gradients obtained at similar resolution currently available in the literature, we measure a negative mass dependence of $\Delta\log({\rm O/H})/\Delta r~ [\mathrm{dex~kpc^{-1}}] = (-0.020\pm0.007) + (-0.014\pm0.008) \log(M_\ast/10^{9.4}M_\odot)$ with the intrinsic scatter being $\sigma=0.060\pm0.006$ over four orders of magnitude in stellar mass. Our result is consistent with strong feedback, not secular processes, being the primary governor of the chemo-structural evolution of star-forming galaxies during the disk mass assembly at cosmic noon. We also find that the intrinsic scatter of metallicity gradients increases with decreasing stellar mass and increasing specific star-formation rate. This increase in the intrinsic scatter is likely caused by the combined effect of cold-mode gas accretion and merger-induced starbursts, with the latter more predominant in the dwarf mass regime of $M_\ast\lesssim10^9 M_\odot$.

We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $\mu$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $\mu$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um.

We report near simultaneous imaging using LMIRCam on the LBTI of the quadruply imaged lensed quasar HS 0810+2554 at wavelengths of 2.16, 3.7 and $4.78~\mu$m with a Full Width Half Max (FWHM) spatial resolution of $0^{\prime\prime}\!\!.13$, $0^{\prime\prime}\!\!.12$ and $0^{\prime\prime}\!\!.15$ respectively, comparable to HST optical imaging. In the $\rm{z} = 1.5$ rest frame of the quasar, the observed wavelengths correspond to 0.86, 1.48, and $1.91~\mu$m respectively. The two brightest images in the quad, A and B, are clearly resolved from each other with a separation of $0.187^{\prime\prime}$. The flux ratio of these two images (A/B) trends from 1.79 to 1.23 from 2.16 to $4.78~\mu$m. The trend in flux ratio is consistent with the $2.16~\mu$m flux originating from a small sized accretion disk in the quasar that experiences only microlensing. The excess flux above the contribution from the accretion disk at the two longer wavelengths originates from a larger sized region that experiences no microlensing. A simple model employing multiplicative factors for image B due to stellar microlensing $(m)$ and sub-structure millilensing $(M)$ is presented. The result is tightly constrained to the product $m\times M=1.79$. Given the observational errors, the 60\% probability contour for this product stretches from $m= 2.6$, $M = 0.69$ to $m= 1.79$, $M = 1.0$, where the later is consistent with microlensing only.

We present a numerical study of the interactions between the elongated AGN outflows representing an evolved, narrow-angle tail (NAT) radio galaxy and planar, transverse ICM shock fronts characteristic of those induced by galaxy cluster mergers (incident Mach numbers 2 - 4). The simulated NAT formation was reported previously in \cite{on19a}. Our simulations utilize a three-dimensional, Eulerian magnetohydrodynamic code along with energy-dependent Eulerian transport of passive cosmic ray electrons. Our analysis of the shock/NAT interaction applies a Riemann problem-based theoretical model to interpret complex shock front behavior during passage through the highly heterogeneous structures of the simulated NAT tails. In addition to shock compression, shock-induced vortical motions are observed within the tails that contribute to coherent turbulent dynamo processes that continue to amplify the magnetic fields in the tails well after initial shock compression. We analyze synthetic radio observations spanning the NAT-shock interaction period, and examine the brightness, spectral and polarization properties of our shock-rejuvenated radio tails, as well as the extent to which the pre-shock states of the plasma and particle populations in our tails influence post-shock observations. Finally, we evaluate our findings in the possible context of a physical analogy to our simulated NAT providing the precursor to a cluster ``radio relic'' associated with an impacting ICM shock.

We report a 3D MHD simulation study of the interactions between radio galaxies and galaxy-cluster-media shocks in which the incident shock normals are orthogonal to the bipolar AGN jets. Before shock impact, light, supersonic jets inflate lobes (cavities) in a static, uniform ICM. We examine three AGN activity scenarios: 1) continued, steady jet activity; 2) jet source cycled off coincident with shock/radio lobe impact; 3) jet activity ceased well before shock arrival (a "radio phoenix" scenario). The simulations follow relativistic electrons (CRe) introduced by the jets, enabling synthetic radio synchrotron images and spectra. Such encounters can be decomposed into an abrupt shock transition and a subsequent long term post shock wind. Shock impact disrupts the pre-formed, low density RG cavities into two ring vortices embedded in the post shock wind. Dynamical processes cause the vortex pair to merge as they propagate downwind somewhat faster than the wind itself. When the AGN jets remain active ram pressure bends the jets downwind, generating a narrow angle tail morphology aligned with the axis of the vortex ring. The deflected jets do not significantly alter dynamical evolution of the vortex ring. However, active jets and their associated tails do dominate the synchrotron emission, compromising the observability of the vortex structures. Downwind-directed momentum concentrated by the jets impacts and alters the post-encounter shock. In the "radio phoenix" scenario, no DSA of the fossil electron population is required to account for the observed brightening and flattening of the spectra, adiabatic compression effects are sufficient.

We report J-band MOSFIRE spectroscopy of a low-mass (log$(M_*/M_\odot)=8.62^{+0.10}_{-0.06}$) star-forming galaxy at $z=0.997$ showing the detection of [NII] and [SII] alongside a strong H$\alpha$ line. We derive a gas-phase metallicity of log$(\text{O}/\text{H})=7.99^{+0.13}_{-0.23}$, placing this object in a region of $M_* - Z$ space that is sparsely populated at this redshift. Furthermore, many existing metallicity measurements in this $M_* - z$ regime are derived from only [NII]/H$\alpha$ (N2), a diagnostic widely used in high-redshift metallicity studies despite the known strong degeneracy with the ionization parameter and resulting large systematic uncertainty. We demonstrate that even in a regime where [NII] and [SII] are at the detection limit and the measurement uncertainty associated with the [NII]/[SII] ratio is high (S/N~3), the more sophisticated Dopita et al. diagnostic provides an improved constraint compared to N2 by reducing the systematic uncertainty due to the ionization parameter. This approach does not, however, dispel uncertainty associated with stochastic or systematic variations in the nitrogen-to-oxygen abundance ratio. While this approach improves upon N2, future progress in extending metallicity studies into this low-mass regime will require larger samples to allow for stochastic variations, as well as careful consideration of the global trends among dwarf galaxies in all physical parameters, not just metallicity.

We present a 3D magnetohydrodynamic (MHD) study of narrow-angle tail (NAT) radio galaxy (RG) dynamics, including passive cosmic ray electrons. We follow evolution of a bipolar-jet RG in a persistent crosswind through hundreds of Myr. We confirm traditional jet-bending models, while noting that our NAT exhibits a transitional morphology reminiscent of wide-angle radio tails. Once deflected, jets remain internally stable, but are intermittently disrupted by external disturbances induced by the NAT dynamics itself. The disruptions enhance jet and tail magnetic fields. Disrupted jet plasma is heterogeneously mixed with denser wind plasma, yielding patchy, filamentary tails that grow longer at a rate exceeding the wind speed. Such fast tail extension could, for example, allow NAT tails to overtake extraneous ICM features, such as shocks and shear layers downwind of where the tails first form. Those events, in turn, could generate enhanced radio emissions within the ICM features themselves that do not follow the geometrical extension of the tails past the encounter. Analysis of synthetic radio observations reveals an extended time period once the NAT has developed in which it displays a nearly steady-state morphology with integrated fluxes that are roughly constant, along with a self-similar, curved integrated spectrum. In an appendix, we outline a simple analytic jet trajectory formalism with one adjustable parameter, using it to illustrate explicit trajectories that extend the classic bending model to arbitrary jet-wind orientations.

Ultraviolet (UV) observations of local star-forming galaxies have begun to establish an empirical baseline for interpreting the rest-UV spectra of reionization-era galaxies. However, existing high-ionization emission line measurements at $z>6$ ($\mathrm{W_{C IV, 0}} \gtrsim 20$ {\AA}) are uniformly stronger than observed locally ($\mathrm{W_{C IV, 0}} \lesssim 2$ {\AA}), likely due to the relatively high metallicities ($Z/Z_\odot > 0.1$) typically probed by UV surveys of nearby galaxies. We present new HST/COS spectra of six nearby ($z<0.01$) extremely metal-poor galaxies (XMPs, $Z/Z_\odot \lesssim 0.1$) targeted to address this limitation and provide constraints on the highly-uncertain ionizing spectra powered by low-metallicity massive stars. Our data reveal a range of spectral features, including one of the most prominent nebular C IV doublets yet observed in local star-forming systems and strong He II emission. Using all published UV observations of local XMPs to-date, we find that nebular C IV emission is ubiquitous in very high specific star formation rate systems at low metallicity, but still find equivalent widths smaller than those measured in individual lensed systems at $z>6$. Our moderate-resolution HST/COS data allow us to conduct an analysis of the stellar winds in a local nebular C IV emitter, which suggests that some of the tension with $z>6$ data may be due to existing local samples not yet probing sufficiently high $\mathrm{\alpha/Fe}$ abundance ratios. Our results indicate that C IV emission can play a crucial role in the JWST and ELT era by acting as an accessible signpost of very low metallicity ($Z/Z_\odot < 0.1$) massive stars in assembling reionization-era systems.

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.

We report the first sub-kiloparsec spatial resolution measurements of strongly inverted gas-phase metallicity gradients in two dwarf galaxies at $z$$\sim$2. The galaxies have stellar masses $\sim$$10^9M_\odot$, specific star-formation rate $\sim$20 Gyr$^{-1}$, and global metallicity $12+\log({\rm O/H})\sim8.1$ (1/4 solar), assuming the Maiolino et al. (2008) strong line calibrations of OIII/Hb and OII/Hb. Their metallicity radial gradients are measured to be highly inverted, i.e., 0.122$\pm$0.008 and 0.111$\pm$0.017 dex/kpc, which is hitherto unseen at such small masses in similar redshift ranges. From the Hubble Space Telescope observations of the source nebular emission and stellar continuum, we present the 2-dimensional spatial maps of star-formation rate surface density, stellar population age, and gas fraction, which show that our galaxies are currently undergoing rapid mass assembly via disk inside-out growth. More importantly, using a simple chemical evolution model, we find that the gas fractions for different metallicity regions cannot be explained by pure gas accretion. Our spatially resolved analysis based on a more advanced gas regulator model results in a spatial map of net gaseous outflows, triggered by active central starbursts, that potentially play a significant role in shaping the spatial distribution of metallicity by effectively transporting stellar nucleosynthesis yields outwards. The relation between wind mass loading factors and stellar surface densities measured in different regions of our galaxies shows that a single type of wind mechanism, driven by either energy or momentum conservation, cannot explain the entire galaxy. These sources present a unique constraint on the effects of gas flows on the early phase of disk growth from the perspective of spatially resolved chemical evolution within individual systems.

We report on polarimetric maps made with HAWC+/SOFIA toward Rho Oph A, the densest portion of the Rho Ophiuchi molecular complex. We employed HAWC+ bands C (89 $\mu$m) and D (154 $\mu$m). The slope of the polarization spectrum was investigated by defining the quantity R_DC = p_D/p_C, where p_C and p_D represent polarization degrees in bands C and D, respectively. We find a clear correlation between R_DC and the molecular hydrogen column density across the cloud. A positive slope (R_DC > 1) dominates the lower density and well illuminated portions of the cloud, that are heated by the high mass star Oph S1, whereas a transition to a negative slope (R_DC < 1) is observed toward the denser and less evenly illuminated cloud core. We interpret the trends as due to a combination of: (1) Warm grains at the cloud outskirts, which are efficiently aligned by the abundant exposure to radiation from Oph S1, as proposed in the radiative torques theory; and (2) Cold grains deep in the cloud core, which are poorly aligned due to shielding from external radiation. To assess this interpretation, we developed a very simple toy model using a spherically symmetric cloud core based on Herschel data, and verified that the predicted variation of R_DC is consistent with the observations. This result introduces a new method that can be used to probe the grain alignment efficiency in molecular clouds, based on the analysis of trends in the far-infrared polarization spectrum.

Observations have revealed massive (logM*/Msun>11) galaxies that were already dead when the universe was only ~2 Gyr. Given the short time before these galaxies were quenched, their past histories and quenching mechanism(s) are of particular interest. In this paper, we study star formation histories (SFHs) of 24 massive galaxies at 1.6<z<2.5. A deep slitless spectroscopy + imaging data set collected from multiple Hubble Space Telescope surveys allows robust determination of their spectral energy distributions and SFHs with no functional assumption on their forms. We find that most of our massive galaxies had formed > 50% of their extant masses by ~1.5 Gyr before the time of observed redshifts, with a trend where more massive galaxies form earlier. Their stellar-phase metallicities are already compatible with those of local early-type galaxies, with a median value of logZ*/Zsun=0.25 and scatter of ~0.15dex. In combination with the reconstructed SFHs, we reveal their rapid metallicity evolution from z~5.5 to ~2.2 at a rate of ~0.2dex/Gyr in log Z*/Zsun. Interestingly, the inferred stellar-phase metallicities are, when compared at half-mass time, ~0.25dex higher than observed gas-phase metallicities of star forming galaxies. While systematic uncertainties remain, this may imply that these quenched galaxies have continued low-level star formation, rather than abruptly terminating their star formation activity, and kept enhancing their metallicity until recently.

We constrain the distribution of spatially offset Lyman-alpha emission (Ly$\alpha$) relative to rest-frame ultraviolet emission in $\sim300$ high redshift ($3<z<5.5$) Lyman-break galaxies (LBGs) exhibiting Ly$\alpha$ emission from VANDELS, a VLT/VIMOS slit-spectroscopic survey of the CANDELS Ultra Deep Survey and Chandra Deep Field South fields (${\simeq0.2}~\mathrm{deg}^2$ total). Because slit spectroscopy compresses two-dimensional spatial information into one spatial dimension, we use Bayesian inference to recover the underlying Ly$\alpha$ spatial offset distribution. We model the distribution using a 2D circular Gaussian, defined by a single parameter $\sigma_{r,\mathrm{Ly}\alpha}$, the standard deviation expressed in polar coordinates. Over the entire redshift range of our sample ($3<z<5.5$), we find $\sigma_{r,\mathrm{Ly}\alpha}=1.70^{+0.09}_{-0.08}$ kpc ($68\%$ conf.), corresponding to $\sim0.25$ arcsec at $\langle z\rangle=4.5$. We also find that $\sigma_{r,\mathrm{Ly}\alpha}$ decreases significantly with redshift. Because Ly$\alpha$ spatial offsets can cause slit-losses, the decrease in $\sigma_{r,\mathrm{Ly}\alpha}$ with redshift can partially explain the increase in the fraction of Ly$\alpha$ emitters observed in the literature over this same interval, although uncertainties are still too large to reach a strong conclusion. If $\sigma_{r,\mathrm{Ly}\alpha}$ continues to decrease into the reionization epoch, then the decrease in Ly$\alpha$ transmission from galaxies observed during this epoch might require an even higher neutral hydrogen fraction than what is currently inferred. Conversely, if spatial offsets increase with the increasing opacity of the IGM, slit losses may explain some of the drop in Ly$\alpha$ transmission observed at $z>6$. Spatially resolved observations of Ly$\alpha$ and UV continuum at $6<z<8$ are needed to settle the issue.

We introduce the OSIRIS Lens-Amplified Survey (OLAS), a kinematic survey of gravitationally lensed galaxies at cosmic noon taken with Keck adaptive optics. In this paper we present spatially resolved spectroscopy and nebular emission kinematic maps for 17 star forming galaxies with stellar masses 8 < log($M_*$/$M_{\odot}$) < 9.8 and redshifts 1.2 < z < 2.3. OLAS is designed to probe the stellar mass ($M_*$) and specific star formation rate (sSFR) range where simulations suggest that stellar feedback is most effective at driving gaseous outflows that create galaxy-wide potential fluctuations which can generate dark matter cores. We compare our kinematic data with the trend between sSFR, $M_*$ and H$\alpha$ velocity dispersion, $\sigma$, from the Feedback In Realistic Environments (FIRE) simulations. Our observations reveal a correlation between sSFR and sigma at fixed $M_*$ that is similar to the trend predicted by simulations: feedback from star formation drives star-forming gas and newly formed stars into more dispersion dominated orbits. The observed magnitude of this effect is in good agreement with the FIRE simulations, in which feedback alters the central density profiles of low mass galaxies, converting dark matter cusps into cores over time. Our data support the scenario that stellar feedback drives gaseous outflows and potential fluctuations, which in turn drive dark matter core formation in dwarf galaxies.

We report from a study utilizing 3D MHD simulations, including cosmic-ray electrons, of the interactions between radio galaxies (RGs) and dynamically active ICMs. Here we consider interactions involving plane ICM shocks having Mach numbers 2--4 and their normals aligned with steady, active bipolar RG jets penetrating uniform, stationary ICMs. Shock impact disrupts the pre-formed RG jet cocoons into ring vortex structures. Sufficiently strong post-shock winds can stop and even reverse the upwind jet, and strip jets to virtually "naked" states, leaving them without a surrounding cocoon. Strong shock-induced vorticity can also disrupt the downwind jet, so that the ring vortex remnant of the cocoons appears ahead of that jet's visible terminus. Magnetic field amplification in the ring vortex can significantly enhance its synchrotron emissions well after the vortex becomes isolated from the RG and its fresh CRe supply. We examine these dynamics and their observable consequences in detail.

HST/STIS spectra of the small clumps and filaments closest to the central star in VY CMa reveal that the very strong K I emission and TiO and VO molecular emission, long thought to form in a dusty circumstellar shell, actually originate in a few small clumps 100's of AU from the star. The K I lines are 10 to 20 times stronger in these nearest ejecta than on the star. The observations also confirm VO as a circumstellar molecule. In this letter we discuss the spectra of the features, their motions and ages, and the identification of the molecular emission. The strength of the atomic and molecular features in the small clumps present an astrophysical problem for the excitation process. We show that the clumps must have a nearly clear line of sight to the star's radiation.

Galaxy clusters are the largest and most massive bound objects resulting from cosmic hierarchical structure formation. Baryons account for somewhat more than 10% of that mass, with roughly 90% of the baryonic matter distributed throughout the clusters as hot ($T>1$ keV), high-$\beta$, very weakly collisional plasma; the so-called "intracluster medium" (ICM). Cluster mergers, close gravitational encounters and accretion, along with violent feedback from galaxies and relativistic jets from active galactic nuclei, drive winds, gravity waves, turbulence and shocks within the ICM. Those dynamics, in turn, generate cluster-scale magnetic fields and accelerate and mediate the transport of high-energy charged particles. Kinetic-scale, collective plasma processes define the basic character and fundamental signatures of these ICM phenomena, which are observed primarily by X-ray and radio astronomers.

We study the impact of cluster environment on the morphology of a sample of 4304 extended radio galaxies from Radio Galaxy Zoo. A total of 87% of the sample lies within a projected 15 Mpc of an optically identified cluster. Brightest cluster galaxies (BCGs) are more likely than other cluster members to be radio sources, and are also moderately bent. The surface density as a function of separation from cluster center of non-BCG radio galaxies follows a power law with index $-1.10\pm 0.03$ out to $10~r_{500}$ ($\sim 7~$Mpc), which is steeper than the corresponding distribution for optically selected galaxies. Non-BCG radio galaxies are statistically more bent the closer they are to the cluster center. Within the inner $1.5~r_{500}$ ($\sim 1~$Mpc) of a cluster, non-BCG radio galaxies are statistically more bent in high-mass clusters than in low-mass clusters. Together, we find that non-BCG sources are statistically more bent in environments that exert greater ram pressure. We use the orientation of bent radio galaxies as an indicator of galaxy orbits and find that they are preferentially in radial orbits. Away from clusters, there is a large population of bent radio galaxies, limiting their use as cluster locators; however, they are still located within statistically overdense regions. We investigate the asymmetry in the tail length of sources that have their tails aligned along the radius vector from the cluster center, and find that the length of the inward-pointing tail is weakly suppressed for sources close to the center of the cluster.

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.

Detections and non-detections of Lyman alpha (Ly$\alpha$) emission from $z>6$ galaxies ($<1$ Gyr after the Big Bang) can be used to measure the timeline of cosmic reionization. Of key interest to measuring reionization's mid-stages, but also increasing observational challenge, are observations at z > 7, where Ly$\alpha$ redshifts to near infra-red wavelengths. Here we present a search for z > 7.2 Ly$\alpha$ emission in 53 intrinsically faint Lyman Break Galaxy candidates, gravitationally lensed by massive galaxy clusters, in the KMOS Lens-Amplified Spectroscopic Survey (KLASS). With integration times of ~7-10 hours, we detect no Ly$\alpha$ emission with S/N>5 in our sample. We determine our observations to be 80% complete for 5$\sigma$ spatially and spectrally unresolved emission lines with integrated line flux $>5.7\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$. We define a photometrically selected sub-sample of 29 targets at $z=7.9\pm0.6$, with a median 5$\sigma$ Ly$\alpha$ EW limit of 58A. We perform a Bayesian inference of the average intergalactic medium (IGM) neutral hydrogen fraction using their spectra. Our inference accounts for the wavelength sensitivity and incomplete redshift coverage of our observations, and the photometric redshift probability distribution of each target. These observations, combined with samples from the literature, enable us to place a lower limit on the average IGM neutral hydrogen fraction of $> 0.76 \; (68\%), \; > 0.46 \; (95\%)$ at z ~ 8, providing further evidence of rapid reionization at z~6-8. We show that this is consistent with reionization history models extending the galaxy luminosity function to $M_\textrm{UV} \lesssim -12$, with low ionizing photon escape fractions, $f_\textrm{esc} \lesssim 15\%$.

We implement a steady, one-dimensional flow model for the X-ray jet of Centaurus A in which entrainment of stellar mass loss is the primary cause of dissipation. Using over 260 ks of new and archival Chandra/ACIS data, we have constrained the temperature, density and pressure distributions of gas in the central regions of the host galaxy of Centaurus A, and so the pressure throughout the length of its jet. The model is constrained by the observed profiles of pressure and jet width, and conserves matter and energy, enabling us to estimate jet velocities, and hence all the other flow properties. Invoking realistic stellar populations within the jet, we find that the increase in its momentum flux exceeds the net pressure force on the jet unless only about one half of the total stellar mass loss is entrained. For self-consistent models, the bulk speed only falls modestly, from ~0.67c to ~0.52c over the range of 0.25-5.94 kpc from the nucleus. The sonic Mach number varies between ~5.3 and 3.6 over this range.

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.

We quantify the star formation (SF) in the inner cores ($\mathcal{R}$/$R_{200}$$\leq$0.3) of 24 massive galaxy clusters at 0.2$\lesssim$$z$$\lesssim$0.9 observed by the $Herschel$ Lensing Survey and the Cluster Lensing and Supernova survey with $Hubble$. These programmes, covering the rest-frame ultraviolet to far-infrared regimes, allow us to accurately characterize stellar mass-limited ($\mathcal{M}_{*}$$>$$10^{10}$ $M_{\odot}$) samples of star-forming cluster members (not)-detected in the mid- and/or far-infrared. We release the catalogues with the photometry, photometric redshifts, and physical properties of these samples. We also quantify the SF displayed by comparable field samples from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. We find that in intermediate-$z$ cluster cores, the SF activity is suppressed with respect the field in terms of both the fraction ($\mathcal{F}$) of star-forming galaxies (SFG) and the rate at which they form stars ($\mathcal{SFR}$ and $s\mathcal{SFR} = \mathcal{SFR}/\mathcal{M}_{*}$). On average, the $\mathcal{F}$ of SFGs is a factor $\sim$$2$ smaller in cluster cores than in the field. Furthermore, SFGs present average $\mathcal{SFR}$ and $s\mathcal{SFR}$ typically $\sim$0.3 dex smaller in the clusters than in the field along the whole redshift range probed. Our results favour long time-scale quenching physical processes as the main driver of SF suppression in the inner cores of clusters since $z$$\sim$0.9, with shorter time-scale processes being very likely responsible for a fraction of the missing SFG population.

We present Far-Infrared polarimetry observations of M82 at 53 and $154~\mu \rm{m}$ and NGC 253 at $89~\mu \rm{m}$, which were taken with HAWC+ in polarimetry mode on the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization of M82 at $53~\mu \rm{m}$ clearly shows a magnetic field geometry perpendicular to the disk in the hot dust emission. For M82 the polarization at $154~\mu \rm{m}$ shows a combination of field geometry perpendicular to the disk in the nuclear region, but closer to parallel to the disk away from the nucleus. The fractional polarization at $53~\mu \rm{m}$ $(154~\mu \rm{m})$ ranges from 7% (3%) off nucleus to 0.5% (0.3%) near the nucleus. A simple interpretation of the observations of M82 invokes a massive polar outflow, dragging the field along, from a region $\sim 700$~pc in diameter that has entrained some of the gas and dust, creating a vertical field geometry seen mostly in the hotter $(53~\mu \rm{m})$ dust emission. This outflow sits within a larger disk with a more typical planar geometry that more strongly contributes to the cooler $(154~\mu \rm{m})$ dust emission. For NGC 253, the polarization at $89~\mu \rm{m}$ is dominated by a planar geometry in the tilted disk, with weak indication of a vertical geometry above and below the plane from the nucleus. The polarization observations of NGC 253 at $53~\mu \rm{m}$ were of insufficient S/N for detailed analysis.

Due to increase in computing power, high-order Eulerian schemes will likely become instrumental for the simulations of turbulence and magnetic field amplification in astrophysical fluids in the next years. We present the implementation of a fifth order weighted essentially non-oscillatory scheme for constrained-transport magnetohydrodynamics into the code WOMBAT. We establish the correctness of our implementation with an extensive number tests. We find that the fifth order scheme performs as accurately as a common second order scheme at half the resolution. We argue that for a given solution quality the new scheme is more computationally efficient than lower order schemes in three dimensions. We also establish the performance characteristics of the solver in the WOMBAT framework. Our implementation fully vectorizes using flattened arrays in thread-local memory. It performs at about 0.6 Million zones per second per node on Intel Broadwell. We present scaling tests of the code up to 98 thousand cores on the Cray XC40 machine "Hazel Hen", with a sustained performance of about 5 percent of peak at scale.

In galaxy clusters, modern radio interferometers observe non-thermal radio sources with unprecedented spatial and spectral resolution. For the first time, the new data allows to infer the structure of the intra-cluster magnetic fields on small scales via Faraday tomography. This leap forward demands new numerical models for the amplification of magnetic fields in cosmic structure formation - the cosmological magnetic dynamo. Here we present a novel numerical approach to astrophyiscal MHD simulations aimed to resolve this small-scale dynamo in future cosmological simulations. As a first step, we implement a fifth order WENO scheme in the new code WOMBAT. We show that this scheme doubles the effective resolution of the simulation and is thus less expensive than common second order schemes. WOMBAT uses a novel approach to parallelization and load balancing developed in collaboration with performance engineers at Cray Inc. This will allow us scale simulation to the exaflop regime and achieve kpc resolution in future cosmological simulations of galaxy clusters. Here we demonstrate the excellent scaling properties of the code and argue that resolved simulations of the cosmological small scale dynamo within the whole virial radius are possible in the next years.

Galactic-scale outflows regulate the stellar mass growth and chemical enrichment of galaxies, yet key outflow properties such as the chemical composition and mass loss rate remain largely unknown. We address these properties with Keck/ESI echellete spectra of nine gravitationally lensed z=2-3 star forming galaxies, probing a range of absorption transitions. Interstellar absorption in our sample is dominated by outflowing material, with typical velocities -150 km/s. Approximately 80% of the total column density is associated with a net outflow. Mass loss rates in the low ionization phase are comparable to or in excess of the star formation rate, with total outflow rates likely higher when accounting for ionized gas. Of order half of the heavy element yield from star formation is ejected in the low ionization phase, confirming that outflows play a critical role in regulating galaxy chemical evolution. Covering fractions vary and are in general non-uniform, with most galaxies having incomplete covering by the low ions across all velocities. Low ion abundance patterns show remarkably little scatter, revealing a distinct "chemical fingerprint" of outflows. Gas phase Si/Fe abundances are significantly super-solar ([Si/Fe]$\gtrsim$0.4) indicating a combination of $\alpha$-enhancement and dust depletion. Derived properties are comparable to the most kinematically broad, metal-rich, and depleted intergalactic absorption systems at similar redshifts, suggesting that these extreme systems are associated with galactic outflows at impact parameters conservatively within a few tens of kpc. We discuss implications of the abundance patterns in z=2-3 galaxies and the role of outflows at this epoch.

Due to their late formation in cosmic history, clusters of galaxies are not fully in hydrostatic equilibrium and the gravitational pull of their mass at a given radius is expected not to be entirely balanced by the thermal gas pressure. Turbulence may supply additional pressure, and recent (X-ray and SZ) hydrostatic mass reconstructions claim a pressure support of $\sim 5-15\%$ of the total pressure at $R_{\rm 200}$. In this work we show that, after carefully disentangling bulk from small-scale turbulent motions in high-resolution simulations of galaxy clusters, we can constrain which fraction of the gas kinetic energy effectively provides pressure support in the cluster's gravitational potential. While the ubiquitous presence of radial inflows in the cluster can lead to significant bias in the estimate of the non-thermal pressure support, we report that only a part of this energy effectively acts as a source of pressure, providing a support of the order of $\sim 10\%$ of the total pressure at $R_{\rm 200}$.

We present spatially resolved analysis of a lensed galaxy, SDSS1958+5950 at $z = 2.225$, from the Cambridge Sloan Survey of Wide Arcs in the Sky (CASSOWARY). We use our new high resolution imaging data to construct a robust lens model for the galaxy group at $z = 0.214$. We employ the updated lens model to combine the Integral Field Spectrographic observations on two highly distorted images of the lensed target. We adopt a forward-modeling approach to deconvolve the effects of point spread function from the combined source-plane reconstruction. The approach is adapted to the lens model magnification and enables a resolution of $\sim$170 pc in the galaxy-source plane. We propose an ongoing merger as the origin of the lensed system on the basis of its source-plane morphology, kinematics and rest-frame emission line ratios. Using our novel technique of adaptive coadded source plane reconstruction, we are able to detect different components in the velocity gradient that were not seen in previous studies of this object, plausibly belonging to different components in the merging system.

EBEX was a long-duration balloon-borne experiment to measure the polarization of the cosmic microwave background. The experiment had three frequency bands centered at 150, 250, and 410 GHz and was the first to use a kilo-pixel array of transition edge sensor (TES) bolometers aboard a balloon platform; shortly after reaching float we operated 504, 342, and 109 TESs at each of the bands, respectively. We describe the design and characterization of the array and the readout system. We give the distributions of measured thermal conductances, normal resistances, and transition temperatures. With the exception of the thermal conductance at 150 GHz. We measured median low-loop-gain time constants $\tau_{0}=$ 88, 46, and 57 ms and compare them to predictions. Two measurements of bolometer absorption efficiency show high ($\sim$0.9) efficiency at 150 GHz and medium ($\sim$0.35, and $\sim$0.25) at the two higher bands, respectively. We measure a median total optical load of 3.6, 5.3 and 5.0 pW absorbed at the three bands, respectively. EBEX pioneered the use of the digital version of the frequency domain multiplexing (FDM) system which multiplexed the bias and readout of 16 bolometers onto two wires. We present accounting of the measured noise equivalent power. The median per-detector noise equivalent temperatures referred to a black body with a temperature of 2.725 K are 400, 920, and 14500 $\mu$K$\sqrt{s}$ for the three bands, respectively. We compare these values to our pre-flight predictions and to a previous balloon payload, discuss the sources of excess noise, and the path for a future payload to make full use of the balloon environment.

Measurements of the galaxy UV luminosity function at z>6 suggest that young stars hosted in low-mass star-forming galaxies produced the bulk of hydrogen-ionizing photons necessary to reionize the intergalactic medium (IGM) by redshift z~6. Whether star-forming galaxies dominated cosmic reionization, however, also depends on their stellar populations and interstellar medium properties, which set, among other things, the production rate of H-ionizing photons, $\xi_\text{ion}^\star$, and the fraction of these escaping into the IGM. Given the difficulty of constraining with existing observatories the physical properties of z>6 galaxies, in this work we focus on a sample of ten nearby objects showing UV spectral features comparable to those observed at z>6. We use the new-generation Beagle tool to model the UV-to-optical photometry and UV/optical emission lines of these Local 'analogues' of high-redshift galaxies, finding that our relatively simple, yet fully self-consistent, physical model can successfully reproduce the different observables considered. Our galaxies span a broad range of metallicities and are characterised by high ionization parameters, low dust attenuation, and very young stellar populations. Through our analysis, we derive a novel diagnostic of the production rate of H-ionizing photons per unit UV luminosity, $\xi_\text{ion}^\star$, based on the equivalent width of the bright $[\text{OIII}] \lambda 4959,5007$ doublet, which does not require measurements of H-recombination lines. This new diagnostic can be used to estimate $\xi_\text{ion}^\star$ from future direct measurements of the $[\text{OIII}] \lambda 4959,5007$ line using JWST/NIRSpec (out to z~9.5), and by exploiting the contamination by $\text{H}\beta + [\text{OIII}] \lambda 4959,5007$ of photometric observations of distant galaxies, for instance from existing Spitzer/IRAC data and from future ones with JWST/NIRCam.

We present the first comprehensive evolutionary analysis of the rest-frame UV spectroscopic properties of star-forming galaxies at z ~ 2-4. We match samples at different redshifts in UV luminosity and stellar mass, and perform systematic measurements of spectral features and stellar population modeling. By creating composite spectra grouped according to Ly$\alpha$ equivalent width (EW), and various galaxy properties, we study the evolutionary trends among Ly$\alpha$, low- and high-ionization interstellar (LIS and HIS) absorption features, and integrated galaxy properties. We also examine the redshift evolution of Ly$\alpha$ and LIS absorption kinematics, and fine-structure emission EWs. The connections among the strengths of Ly$\alpha$, LIS lines, and dust extinction are redshift-independent, as is the decoupling of Ly$\alpha$ and HIS line strengths, and the bulk outflow kinematics as traced by LIS lines. Stronger Ly$\alpha$ emission is observed at higher redshift at fixed UV luminosity, stellar mass, SFR, and age. Much of this variation in average Ly$\alpha$ strength with redshift, and the variation in Ly$\alpha$ strength at fixed redshift, can be explained in terms of variations in neutral gas covering fraction and/or dust content in the ISM and CGM. However, based on the connection between Ly$\alpha$ and CIII] emission strengths, we additionally find evidence for variations in the intrinsic production rate of Ly$\alpha$ photons at the highest Ly$\alpha$ EWs. The challenge now is to understand the observed evolution in neutral gas covering fraction and dust extinction within a coherent model for galaxy formation, and make robust predictions for the escape of ionizing radiation at z > 6.

New MMT/MIRAC (9-11 {\mu}m), SOFIA/FORCAST (11-37 {\mu}m), and Herschel/PACS (70 and 160 {\mu}m) infrared (IR) imaging and photometry is presented for three famous OH/IR red supergiants (NML Cyg, VX Sgr, and S Per) and two normal red supergiants (RS Per and T Per). We model the observed spectral energy distributions (SEDs) using radiative transfer code DUSTY. Azimuthal average profiles from the SOFIA/FORCAST imaging, in addition to dust mass distribution profiles from DUSTY, constrain the mass-loss histories of these supergiants. For all of our observed supergiants, the DUSTY models suggest that constant mass-loss rates do not produce enough dust to explain the observed infrared emission in the stars' SEDs. Combining our results with Shenoy et al. (2016) (Paper I) we find mixed results with some red supergiants showing evidence for variable and high mass-loss events while others have constant mass loss over the past few thousand years.

Local massive early-type galaxies are believed to have completed most of their star formation $\sim10$Gyr ago and evolved without having substantial star formation since. If so, their progenitors should have roughly solar stellar metallicities ($Z_*$), comparable to their values today. We report the discovery of two lensed massive ($\log M_*/M_\odot\sim11$), $z\sim2.2$ dead galaxies, that appear markedly metal deficient given this scenario. Using 17-band $HST$+$K_{s}$+$Spitzer$ photometry and deep $HST$ grism spectra from the GLASS and SN Refsdal follow-up campaigns covering features near $\lambda_{\rm rest}\sim4000$\AA, we find these systems to be dominated by A-type stars with $\log Z_*/Z_\odot=-0.40\pm0.02$ and $-0.49\pm0.03$ ($30$-$40\%$ solar) under standard assumptions. The second system's lower metallicity is robust to isochrone changes, though this choice can drive the first system's from $\log Z_*/Z_\odot=-0.6$ to 0.1. If these two galaxies are representative of larger samples, this finding suggests that evolutionary paths other than dry minor-merging are required for these massive galaxies. Future analyses with direct metallicity measurements-e.g., by the $James\ Webb\ Space\ Telescope$-will provide critical insight into the nature of such phenomena.

We discuss a systematic effect associated with measuring polarization with a continuously rotating half-wave plate. The effect was identified with the data from the E and B Experiment (EBEX), which was a balloon-borne instrument designed to measure the polarization of the CMB as well as that from Galactic dust. The data show polarization fraction larger than 10\% while less than 3\% were expected from instrumental polarization. We give evidence that the excess polarization is due to detector non-linearity in the presence of a continuously rotating HWP. The non-linearity couples intensity signals into polarization. We develop a map-based method to remove the excess polarization. Applying this method for the 150 (250) GHz bands data we find that 81\% (92\%) of the excess polarization was removed. Characterization and mitigation of this effect is important for future experiments aiming to measure the CMB B-modes with a continuously rotating HWP.

We present deep spectroscopic observations of a Lyman-break galaxy candidate (hereafter MACS1149-JD) at $z\sim9.5$ with the $\textit{Hubble}$ Space Telescope ($\textit{HST}$) WFC3/IR grisms. The grism observations were taken at 4 distinct position angles, totaling 34 orbits with the G141 grism, although only 19 of the orbits are relatively uncontaminated along the trace of MACS1149-JD. We fit a 3-parameter ($z$, F160W mag, and Ly$\alpha$ equivalent width) Lyman-break galaxy template to the three least contaminated grism position angles using an MCMC approach. The grism data alone are best fit with a redshift of $z_{\mathrm{grism}}=9.53^{+0.39}_{-0.60}$ ($68\%$ confidence), in good agreement with our photometric estimate of $z_{\mathrm{phot}}=9.51^{+0.06}_{-0.12}$ ($68\%$ confidence). Our analysis rules out Lyman-alpha emission from MACS1149-JD above a $3\sigma$ equivalent width of 21 \AA{}, consistent with a highly neutral IGM. We explore a scenario where the red $\textit{Spitzer}$/IRAC $[3.6] - [4.5]$ color of the galaxy previously pointed out in the literature is due to strong rest-frame optical emission lines from a very young stellar population rather than a 4000 \AA{} break. We find that while this can provide an explanation for the observed IRAC color, it requires a lower redshift ($z\lesssim9.1$), which is less preferred by the $\textit{HST}$ imaging data. The grism data are consistent with both scenarios, indicating that the red IRAC color can still be explained by a 4000 \AA{} break, characteristic of a relatively evolved stellar population. In this interpretation, the photometry indicate that a $340^{+29}_{-35}$ Myr stellar population is already present in this galaxy only $\sim500~\mathrm{Myr}$ after the Big Bang.

Nearby dwarf galaxies provide a unique laboratory in which to test stellar population models below $Z_\odot/2$. Such tests are particularly important for interpreting the surprising high-ionization UV line emission detected at $z>6$ in recent years. We present HST/COS ultraviolet spectra of ten nearby metal-poor star-forming galaxies selected to show He II emission in SDSS optical spectra. The targets span nearly a dex in gas-phase oxygen abundance ($7.8<12+\log\mathrm{O/H}<8.5$) and present uniformly large specific star formation rates (sSFR $\sim 10^2$ $\mathrm{Gyr}^{-1}$). The UV spectra confirm that metal-poor stellar populations can power extreme nebular emission in high-ionization UV lines, reaching C III] equivalent widths comparable to those seen in systems at $z\sim 6-7$. Our data reveal a marked transition in UV spectral properties with decreasing metallicity, with systems below $12+\log\mathrm{O/H}\lesssim 8.0$ ($Z/Z_\odot \lesssim 1/5$) presenting minimal stellar wind features and prominent nebular emission in He II and C IV. This is consistent with nearly an order of magnitude increase in ionizing photon production beyond the $\mathrm{He^+}$-ionizing edge relative to H-ionizing flux as metallicity decreases below a fifth solar, well in excess of standard stellar population synthesis predictions. Our results suggest that often neglected sources of energetic radiation such as stripped binary products and very massive O-stars produce a sharper change in the ionizing spectrum with decreasing metallicity than expected. Consequently, nebular emission in C IV and He II powered by these stars may provide useful metallicity constraints in the reionization era.

Turbulence generated by large-scale motions during structure formation affects the evolution of the thermal and non-thermal components of the intracluster medium. As enstrophy is a measure of the magnitude of vorticity, we study the generation and evolution of turbulence by analysing the Lagrangian history of enstrophy. For this purpose we combine cosmological simulations carried out with the ENZO-code with our Lagrangian post-processing tool CRaTer. This way we are able to quantify the individual source terms of enstrophy in the course of the accretion of groups onto galaxy clusters. Here we focus on the redshift range from $z=1$ to $z=0$. Finally, we measure the rate of dissipation of turbulence and estimate the resulting amplification of intracluster magnetic fields. We find that compressive and baroclinic motions are the main sources of enstrophy, while stretching motions and dissipation affect most of the ensuing enstrophy evolution. The rate of turbulent dissipation is able to sustain the amplification of intracluster magnetic fields to observed levels.

According to structure formation simulations, weak shocks with typical Mach number, $M_{\rm s}\lesssim 3$, are expected to form in merging galaxy clusters. The presence of such shocks has been indicated by X-ray and radio observations of many merging clusters. In particular, diffuse radio sources known as radio relics could be explained by synchrotron-emitting electrons accelerated via diffusive shock acceleration (Fermi I) at quasi-perpendicular shocks. Here we also consider possible roles of stochastic acceleration (Fermi II) by compressive MHD turbulence downstream of the shock. Then we explore a puzzling discrepancy that for some radio relics, the shock Mach number inferred from the radio spectral index is substantially larger than that estimated from X-ray observations. This problem could be understood, if shock surfaces associated with radio relics consist of multiple shocks with different strengths.In that case, X-ray observations tend to pick up the part of shocks with lower Mach numbers and higher kinetic energy flux, while radio emissions come preferentially from the part of shocks with higher Mach numbers and higher cosmic ray (CR) production. We also show that the Fermi I reacceleration model with preexisting fossil electrons supplemented by Fermi II acceleration due to postshock turbulence could reproduce observed profiles of radio flux densities and integrated radio spectra of two giant radio relics. This study demonstrates the CR electrons can be accelerated at collisionless shocks in galaxy clusters just like supernova remnant shock in the interstellar medium and interplanetary shocks in the solar wind.

Although many of the observed properties of giant radio relics detected in the outskirts of galaxy clusters can be explained by relativistic electrons accelerated at merger-driven shocks, significant puzzles remain. In the case of the so-called Toothbrush relic, the shock Mach number estimated from X-ray observations ($M_{\rm X}\approx1.2-1.5$) is substantially weaker than that inferred from the radio spectral index ($M_{\rm rad}\approx2.8$).Toward understanding such a discrepancy, we here consider the following diffusive shock acceleration (DSA) models:(1) weak-shock models with $M_{\rm s}\lesssim 2$ and a preexisting population of cosmic-ray electrons (CRe) with a flat energy spectrum,and (2) strong-shock models with $M_{\rm s}\approx3$ and either shock-generated suprathermal electrons or preexisting fossil CRe. We calculate the synchrotron emission from the accelerated CRe, following the time evolution of the electron DSA, and subsequent radiative cooling and postshock turbulent acceleration (TA). We find that both models could reproduce reasonably well the observed integrated radio spectrum of the Toothbrush relic, but the observed broad transverse profile requires the stochastic acceleration by downstream turbulence, which we label "turbulent acceleration" or TA to distinguish it from DSA. Moreover, to account for the almost uniform radio spectral index profile along the length of the relic, the weak-shock models require a preshock region over 400~kpc with a uniform population of preexisting CRe with a high cutoff energy ($\gtrsim 40$ GeV). Due to the short cooling time, it is challenging to explain the origin of such energetic electrons. Therefore, we suggest the strong-shock models with low-energy seed CRe ($\lesssim 150$~MeV) are preferred for the radio observations of this relic.

The E and B Experiment (EBEX) was a long-duration balloon-borne cosmic microwave background polarimeter that flew over Antarctica in 2013. We describe the experiment's optical system, receiver, and polarimetric approach, and report on their in-flight performance. EBEX had three frequency bands centered on 150, 250, and 410 GHz. To make efficient use of limited mass and space we designed a 115 cm$^{2}$sr high throughput optical system that had two ambient temperature mirrors and four anti-reflection coated polyethylene lenses per focal plane. All frequency bands shared the same optical train. Polarimetry was achieved with a continuously rotating achromatic half-wave plate (AHWP) that was levitated with a superconducting magnetic bearing (SMB). Rotation stability was 0.45 % over a period of 10 hours, and angular position accuracy was 0.01 degrees. This is the first use of a SMB in astrophysics. The measured modulation efficiency was above 90 % for all bands. To our knowledge the 109 % fractional bandwidth of the AHWP was the broadest implemented to date. The receiver that contained one lens and the AHWP at a temperature of 4 K, the polarizing grid and other lenses at 1 K, and the two focal planes at 0.25 K performed according to specifications giving focal plane temperature stability with fluctuation power spectrum that had $1/f$ knee at 2 mHz. EBEX was the first balloon-borne instrument to implement technologies characteristic of modern CMB polarimeters including high throughput optical systems, and large arrays of transition edge sensor bolometric detectors with mutiplexed readouts.

The E and B Experiment (EBEX) was a long-duration balloon-borne instrument designed to measure the polarization of the cosmic microwave background (CMB) radiation. EBEX was the first balloon-borne instrument to implement a kilo-pixel array of transition edge sensor (TES) bolometric detectors and the first CMB experiment to use the digital version of the frequency domain multiplexing system for readout of the TES array. The scan strategy relied on 40 s peak-to-peak constant velocity azimuthal scans. We discuss the unique demands on the design and operation of the payload that resulted from these new technologies and the scan strategy. We describe the solutions implemented including the development of a power system designed to provide a total of at least 2.3 kW, a cooling system to dissipate 590 W consumed by the detectors' readout system, software to manage and handle the data of the kilo-pixel array, and specialized attitude reconstruction software. We present flight performance data showing faultless management of the TES array, adequate powering and cooling of the readout electronics, and constraint of attitude reconstruction errors such that the spurious B-modes they induced were less than 10% of CMB B-mode power spectrum with $r=0.05$.

We use Atacama Large Millimeter/submillimeter Array Band 5 science verification observations of the red supergiant VY CMa to study the polarization of SiO thermal/masers lines and dust continuum at ~1.7 mm wavelength. We analyse both linear and circular polarization and derive the magnetic field strength and structure, assuming the polarization of the lines originates from the Zeeman effect, and that of the dust originates from aligned dust grains. We also discuss other effects that could give rise to the observed polarization. We detect, for the first time, significant polarization (~3%) of the circumstellar dust emission at millimeter wavelengths. The polarization is uniform with an electric vector position angle of $\sim8^\circ$. Varying levels of linear polarization are detected for the J=4-3 28SiO v=0, 1, 2, and 29SiO v=0, 1 lines, with the strongest polarization fraction of ~30% found for the 29SiO v=1 maser. The linear polarization vectors rotate with velocity, consistent with earlier observations. We also find significant (up to ~1%) circular polarization in several lines, consistent with previous measurements. We conclude that the detection is robust against calibration and regular instrumental errors, although we cannot yet fully rule out non-standard instrumental effects. Emission from magnetically aligned grains is the most likely origin of the observed continuum polarization. This implies that the dust is embedded in a magnetic field >13 mG. The maser line polarization traces the magnetic field structure. The magnetic field in the gas and dust is consistent with an approximately toroidal field configuration, but only higher angular resolution observations will be able to reveal more detailed field structure. If the circular polarization is due to Zeeman splitting, it indicates a magnetic field strength of ~1-3 Gauss, consistent with previous maser observations.

The CIII] and CIV rest-frame UV emission lines are powerful probes of the ionizations states of galaxies. They have furthermore been suggested as alternatives for spectroscopic redshift confirmation of objects at the epoch of reionization ($z>6$), where the most frequently used redshift indicator, Ly$\alpha$, is attenuated by the high fraction of neutral hydrogen in the inter-galactic medium. However, currently only very few confirmations of carbon UV lines at these high redshifts exist, making it challenging to quantify these claims. Here, we present the detection of CIV$\lambda\lambda$1548,1551\AA\ in \HST\ slitless grism spectroscopy obtained by GLASS of a Ly$\alpha$ emitter at $z=6.11$ multiply imaged by the massive foreground galaxy cluster RXJ2248. The CIV emission is detected at the 3--5$\sigma$ level in two images of the source, with marginal detection in two other images. We do not detect significant CIII]$\lambda\lambda$1907,1909\AA\ emission implying an equivalent width EW$_\textrm{CIII]}<20$\AA\ (1$\sigma$) and $\textrm{CIV/CIII}>0.7$ (2$\sigma$). Combined with limits on the rest-frame UV flux from the HeII$\lambda$1640\AA\ emission line and the OIII]$\lambda\lambda$1661,1666\AA\ doublet, we put constraints on the metallicity and the ionization state of the galaxy. The estimated line ratios and equivalent widths do not support a scenario where an AGN is responsible for ionizing the carbon atoms. SED fits including nebular emission lines imply a source with a mass of log(M/M$_\odot)\sim9$, SFR of around 10M$_\odot$/yr, and a young stellar population $<50$Myr old. The source shows a stronger ionizing radiation field than objects with detected CIV emission at $z<2$ and adds to the growing sample of low-mass (log(M/M$_\odot)\lesssim9$) galaxies at the epoch of reionization with strong radiation fields from star formation.

We present the first results of the KMOS Lens-Amplified Spectroscopic Survey (KLASS), a new ESO Very Large Telescope (VLT) large program, doing multi-object integral field spectroscopy of galaxies gravitationally lensed behind seven galaxy clusters selected from the HST Grism Lens-Amplified Survey from Space (GLASS). Using the power of the cluster magnification we are able to reveal the kinematic structure of 25 galaxies at $0.7 \lesssim z \lesssim 2.3$, in four cluster fields, with stellar masses $8 \lesssim \log{(M_\star/M_\odot)} \lesssim 11$. This sample includes 5 sources at $z>1$ with lower stellar masses than in any previous kinematic IFU surveys. Our sample displays a diversity in kinematic structure over this mass and redshift range. The majority of our kinematically resolved sample is rotationally supported, but with a lower ratio of rotational velocity to velocity dispersion than in the local universe, indicating the fraction of dynamically hot disks changes with cosmic time. We find no galaxies with stellar mass $<3 \times 10^9 M_\odot$ in our sample display regular ordered rotation. Using the enhanced spatial resolution from lensing, we resolve a lower number of dispersion dominated systems compared to field surveys, competitive with findings from surveys using adaptive optics. We find that the KMOS IFUs recover emission line flux from HST grism-selected objects more faithfully than slit spectrographs. With artificial slits we estimate slit spectrographs miss on average 60% of the total flux of emission lines, which decreases rapidly if the emission line is spatially offset from the continuum.

We report the discovery of the quadruply lensed quasar J1433+6007, mined in the SDSS DR12 photometric catalogues using a novel outlier-selection technique, without prior spectroscopic or UV excess information. Discovery data obtained at the Nordic Optical telescope (NOT, La Palma) show nearly identical quasar spectra at $z_s=2.74$ and four quasar images in a fold configuration, one of which sits on a blue arc. The deflector redshift is $z_{l}=0.407,$ from Keck-ESI spectra. We describe the selection procedure, discovery and follow-up, image positions and $BVRi$ magnitudes, and first results and forecasts from simple lens models.

We present a new code for astrophysical magneto-hydrodynamics specifically designed and optimized for high performance and scaling on modern and future supercomputers. We describe a novel hybrid OpenMP/MPI programming model that emerged from a collaboration between Cray, Inc. and the University of Minnesota. This design utilizes MPI-RMA optimized for thread scaling, which allows the code to run extremely efficiently at very high thread counts ideal for the latest generation of the multi-core and many-core architectures. Such performance characteristics are needed in the era of "exascale" computing. We describe and demonstrate our high-performance design in detail with the intent that it may be used as a model for other, future astrophysical codes intended for applications demanding exceptional performance.

We present the results of ALMA spectroscopic follow-up of a $z=6.765$ Lyman-$\alpha$ emitting galaxy behind the cluster RXJ1347-1145. We report the detection of [CII]158$\mu$m line fully consistent with the Lyman-$\alpha$ redshift and with the peak of the optical emission. Given the magnification of $\mu=5.0 \pm 0.3$ the intrinsic (corrected for lensing) luminosity of the [CII] line is $L_{[CII]} =1.4^{+0.2}_{-0.3} \times 10^7L_{\odot}$, which is ${\sim}5$ times fainter than other detections of $z\sim 7$ galaxies. The result indicates that low $L_{[CII]}$ in $z\sim 7$ galaxies compared to the local counterparts might be caused by their low metallicities and/or feedback. The small velocity off-set ($\Delta v = 20_{-40}^{+140} \rm km/s$) between the Lyman-$\alpha$ and [CII] line is unusual, and may be indicative of ionizing photons escaping.

A number of exciting concepts are under development for Flagship, Probe class, Explorer class, and Suborbital class NASA missions in the ultraviolet/optical spectral ranges. These missions will depend on high performance silicon detector arrays being delivered affordably and in high numbers. In a focused effort we have advanced delta-doping technology to high throughput and high yield wafer-scale processing, encompassing a multitude of state-of-the-art silicon-based detector formats and designs. As part of this technology advancement and in preparation for upcoming missions, we have embarked on a number of field observations, instrument integrations, and independent evaluations of delta-doped arrays. In this paper, we present recent data and innovations from the Advanced Detectors and Systems program at JPL, including two-dimensional doping technology; our end-to-end post-fabrication processing of high performance UV/Optical/NIR arrays; and advanced coatings for detectors and optical elements. Additionally, we present examples of past, in-progress, and planned observations and deployments of delta-doped arrays.

We report on multilayer high efficiency antireflection coating (ARC) design and development for use at UV wavelengths on CCDs and other Si-based detectors. We have previously demonstrated a set of single-layer coatings, which achieve >50% quantum efficiency (QE) in four bands from 130 to 300 nm. We now present multilayer coating designs that significantly outperform our previous work between 195 and 215 nm. Using up to 11 layers, we present several model designs to reach QE above 80%. We also demonstrate the successful performance of 5 and 11 layer ARCs on silicon and fused silica substrates. Finally, we present a five-layer coat- ing deposited onto a thinned, delta-doped CCD and demonstrate external QE greater than 60% between 202 and 208 nm, with a peak of 67.6% at 206 nm.

In this paper we lay out a simple set of relationships connecting the dynamics of fast plasma jets to the dynamical state of their ambient media. The objective is to provide a tool kit that can be used to connect the morphologies of radio AGNs in galaxy clusters to the dynamical state of the local ICM. The formalism is intended to apply to jets whether they are relativistic or non-relativistic. Special attention is paid to interactions involving ICM shocks, although the results can be applied more broadly. Our formalism emphasizes the importance of the relative Mach number of the impacting ICM flow and the internal Mach number of the AGN jet in determining how the AGN outflows evolve.

We present images and spectral energy distributions (SEDs) of massive young stellar objects (YSOs) in three star-forming H II regions of the Large Magellanic Cloud: N159A, N159 Papillon, and N160. We use photometry from SOFIA/FORCAST at 25.3--37.1 um to constrain model fits to the SEDs and determine luminosities, ages, and dust content of the embedded YSOs and their local environments. By placing these sources on mid-infrared color-magnitude and color-color diagrams, we analyze their dust properties and consider their evolutionary status. Since each object in the FORCAST images has an obvious bright near-infrared counterpart in Spitzer Space Telescope images, we do not find any evidence for new, very cool, previously-undiscovered Class 0 YSOs. Additionally, based on its mid-infrared colors and model parameters, N159A is younger than N160 and the Papillon. The nature of the first extragalactic protostars in N159, P1 and P2, is also discussed.

(Abridged) We combine deep HST grism spectroscopy with a new Bayesian method to derive maps of gas-phase metallicity, nebular dust extinction, and star-formation rate for 10 star-forming galaxies at high redshift ($1.2<z<2.3$). Exploiting lensing magnification by the foreground cluster MACS1149.6+2223, we reach sub-kpc spatial resolution and push the stellar mass limit associated with such high-z spatially resolved measurements below $10^8M_\odot$ for the first time. Our maps exhibit diverse morphologies, indicative of various effects such as efficient radial mixing from tidal torques, rapid accretion of low-metallicity gas, etc., which can affect the gas and metallicity distributions in individual galaxies. Based upon an exhaustive sample of all existing sub-kpc metallicity gradients at high-z, we find that predictions given by analytical chemical evolution models assuming a relatively extended star-formation profile in the early disk formation phase can explain the majority of observed gradients, without involving galactic feedback or radial outflows. We observe a tentative correlation between stellar mass and metallicity gradient, consistent with the downsizing galaxy formation picture that more massive galaxies are more evolved into a later phase of disk growth, where they experience more coherent mass assembly at all radii and thus show shallower metallicity gradients. In addition, we compile a sample of homogeneously cross-calibrated integrated metallicity measurements spanning three orders of magnitude in stellar mass at $z\sim1.8$. We use this sample to study the mass-metallicity relation (MZR) and test the fundamental metallicity relation (FMR). The slope of the observed MZR can rule out the momentum-driven wind model at 3-$\sigma$ confidence level. We find no significant offset with respect to the FMR, taking into account the intrinsic scatter and measurement uncertainties.

We present high-angular (0.17$-$0.35 arcsec) resolution imaging polarimetric observations of Mrk 231 in the 3.1 $\mu$m filter using MMT-Pol on the 6.5-m MMT, and in the 8.7 $\mu$m, 10.3 $\mu$m, and 11.6 $\mu$m filters using CanariCam on the 10.4-m Gran Telescopio CANARIAS. In combination with already published observations, we compile the 1$-$12 $\mu$m total and polarized nuclear spectral energy distribution (SED). The total flux SED in the central 400 pc is explained as the combination of 1) a hot (731 $\pm$ 4 K) dusty structure, directly irradiated by the central engine, which is at 1.6 $\pm$ 0.1 pc away and attributed to be in the pc-scale polar region, 2) an optically-thick, smooth and disk-like dusty structure (`torus') with an inclination of 48 $\pm$ 23$^{\circ}$ surrounding the central engine, and 3) an extinguished (A$_{\mbox{V}} =$ 36 $\pm$ 5 mag) starburst component. The polarized SED decreases from 0.77 $\pm$ 0.14 per cent at 1.2 $\mu$m to 0.31 $\pm$ 0.15 per cent at 11.6 $\mu$m and follows a power-law function, $\lambda^{\sim0.57}$. The polarization angle remains constant ($\sim$108$^{\circ}$) in the 1$-$12 $\mu$m wavelength range. The dominant polarization mechanism is explained as scattering off hot dust grains in the pc-scale polar regions.

Turbulence is a key ingredient for the evolution of the intracluster medium, whose properties can be predicted with high resolution numerical simulations. We present initial results on the generation of solenoidal and compressive turbulence in the intracluster medium during the formation of a small-size cluster using highly resolved, non-radiative cosmological simulations, with a refined monitoring in time. In this first of a series of papers, we closely look at one simulated cluster whose formation was distinguished by a merger around $z \sim 0.3$. We separate laminar gas motions, turbulence and shocks with dedicated filtering strategies and distinguish the solenoidal and compressive components of the gas flows using Hodge-Helmholtz decomposition. Solenoidal turbulence dominates the dissipation of turbulent motions ($\sim 95\%$) in the central cluster volume at all epochs. The dissipation via compressive modes is found to be more important ($\sim 30 \%$ of the total) only at large radii ($\geq 0.5 ~r_{\rm vir}$) and close to merger events. We show that enstrophy (vorticity squared) is good proxy of solenoidal turbulence. All terms ruling the evolution of enstrophy (i.e. baroclinic, compressive, stretching and advective terms) are found to be significant, but in amounts that vary with time and location. Two important trends for the growth of enstrophy in our simulation are identified: first, enstrophy is continuously accreted into the cluster from the outside, and most of that accreted enstrophy is generated near the outer accretion shocks by baroclinic and compressive processes. Second, in the cluster interior vortex stretching is dominant, although the other terms also contribute substantially.

The fraction of ionizing photons escaping from high-redshift star-forming galaxies remains a key obstacle in evaluating whether galaxies were the primary agents of cosmic reionization. We previously proposed using the covering fraction of low-ionization gas, measured via deep absorption line spectroscopy, as a proxy. We now present a significant update, sampling seven gravitationally-lensed sources at $4<z<5$. We show that the absorbing gas in our sources is spatially inhomogeneous with a median covering fraction of 66\%. Correcting for reddening according to a dust-in-cloud model, this implies an estimated absolute escape fraction of $\simeq19\pm6$\%. With possible biases and uncertainties, collectively we find the average escape fraction could be reduced to no less than 11\%, excluding the effect of spatial variations. For one of our lensed sources, we have sufficient signal/noise to demonstrate the presence of such spatial variations and scatter in its dependence on the Ly$\alpha$ equivalent width, consistent with recent simulations. If this source is typical, our lower limit to the escape fraction could be reduced by a further factor $\simeq$2. Across our sample, we find a modest anti-correlation between the inferred escape fraction and the local star formation rate, consistent with a time delay between a burst and leaking Lyman continuum photons. Our analysis demonstrates considerable variations in the escape fraction consistent with being governed by the small-scale behavior of star-forming regions, whose activities fluctuate over short timescales. This supports the suggestion that the escape fraction may increase toward the reionization era when star formation becomes more energetic and burst-like.

We have discovered a previously unreported poor cluster of galaxies (RGZ-CL J0823.2+0333) through an unusual giant wide-angle tail radio galaxy found in the Radio Galaxy Zoo project. We obtained a spectroscopic redshift of $z=0.0897$ for the E0-type host galaxy, 2MASX J08231289+0333016, leading to M$_r = -22.6$ and a $1.4\,$GHz radio luminosity density of $L_{\rm 1.4} = 5.5\times10^{24}$ W Hz$^{-1}$. These radio and optical luminosities are typical for wide-angle tailed radio galaxies near the borderline between Fanaroff-Riley (FR) classes I and II. The projected largest angular size of $\approx8\,$arcmin corresponds to $800\,$kpc and the full length of the source along the curved jets/trails is $1.1\,$Mpc in projection. X-ray data from the XMM-Newton archive yield an upper limit on the X-ray luminosity of the thermal emission surrounding RGZ J082312.9+033301,at $1.2-2.6\times10^{43}$ erg s$^{-1}$ for assumed intra-cluster medium temperatures of $1.0-5.0\,$keV. Our analysis of the environment surrounding RGZ J082312.9+033301 indicates that RGZ J082312.9+033301 lies within a poor cluster. The observed radio morphology suggests that (a) the host galaxy is moving at a significant velocity with respect to an ambient medium like that of at least a poor cluster, and that (b) the source may have had two ignition events of the active galactic nucleus with $10^7\,$yrs in between. This reinforces the idea that an association between RGZ J082312.9+033301, and the newly discovered poor cluster exists.

The supernova impostor PSN J09132750+7627410 in NGC 2748 reached a maximum luminosity of approximately -14 mag. It was quickly realized that its was not a true supernova, but another example of a non-terminal giant eruption. PSN J09132750+7627410 is distinguished by multiple P Cygni absorption minima in the Balmer emission lines that correspond to outflow velocities of -400, -1100, and -1600 km/s. Multiple outflows have been observed in only a few other objects. In this paper we describe the evolution of the spectrum and the P Cygni profiles for three months past maximum, the post-maximum formation of a cool, dense wind, and the identification of a possible progenitor. One of the possible progenitors is an infrared source. Its pre-eruption spectral energy distribution suggests a bolometric luminosity of -8.3 mag and a dust temperature of 780 degrees K. If it is the progenitor it is above the AGB limit unlike the intermediate luminosity red transients. The three P Cygni profiles could be due to ejecta from the current eruption, the wind of the progenitor, or previous mass loss events. We suggest that they were all formed as part of the same high mass loss event and are due to material ejected at different velocities or energies. We also suggest that multiple outflows during giant eruptions may be more common than reported.

Recent supernova and transient surveys have revealed an increasing number of non-terminal stellar eruptions. Though the progenitor class of these eruptions includes the most luminous stars, little is known of the pre-supernova mechanics of massive stars in their most evolved state, thus motivating a census of possible progenitors. From surveys of evolved and unstable luminous star populations in nearby galaxies, we select a sample of yellow and red supergiant candidates in M31 and M33 for review of their spectral characteristics and spectral energy distributions. Since the position of intermediate and late-type supergiants on the color-magnitude diagram can be heavily contaminated by foreground dwarfs, we employ spectral classification and multi-band photometry from optical and near-infrared surveys to confirm membership. Based on spectroscopic evidence for mass loss and the presence of circumstellar dust in their SEDs, we find that $30-40\%$ of the yellow supergiants are likely in a post-red supergiant state. Comparison with evolutionary tracks shows that these mass-losing, post-RSGs have initial masses between $20-40\,M_{\odot}$. More than half of the observed red supergiants in M31 and M33 are producing dusty circumstellar ejecta. We also identify two new warm hypergiants in M31, J004621.05+421308.06 and J004051.59+403303.00, both of which are likely in a post-RSG state.

We discuss spatially resolved emission line spectroscopy secured for a total sample of 15 gravitationally lensed star-forming galaxies at a mean redshift of $z\simeq2$ based on Keck laser-assisted adaptive optics observations undertaken with the recently-improved OSIRIS integral field unit (IFU) spectrograph. By exploiting gravitationally lensed sources drawn primarily from the CASSOWARY survey, we sample these sub-L$^{\ast}$ galaxies with source-plane resolutions of a few hundred parsecs ensuring well-sampled 2-D velocity data and resolved variations in the gas-phase metallicity. Such high spatial resolution data offers a critical check on the structural properties of larger samples derived with coarser sampling using multiple-IFU instruments. We demonstrate how kinematic complexities essential to understanding the maturity of an early star-forming galaxy can often only be revealed with better sampled data. Although we include four sources from our earlier work, the present study provides a more representative sample unbiased with respect to emission line strength. Contrary to earlier suggestions, our data indicates a more diverse range of kinematic and metal gradient behavior inconsistent with a simple picture of well-ordered rotation developing concurrently with established steep metal gradients in all but merging systems. Comparing our observations with the predictions of hydrodynamical simulations suggests that gas and metals have been mixed by outflows or other strong feedback processes, flattening the metal gradients in early star-forming galaxies.

The E and B Experiment (EBEX) is a balloon-borne polarimeter designed to measure the polarization of the cosmic microwave background radiation and to characterize the polarization of galactic dust. EBEX was launched December 29, 2012 and circumnavigated Antarctica observing $\sim$6,000 square degrees of sky during 11 days at three frequency bands centered around 150, 250 and 410 GHz. EBEX was the first experiment to operate a kilo-pixel array of transition-edge sensor bolometers and a continuously rotating achromatic half-wave plate aboard a balloon platform. It also pioneered the use of detector readout based on digital frequency domain multiplexing. We describe the temperature calibration of the experiment. The gain response of the experiment is calibrated using a two-step iterative process. We use signals measured on passes across the Galactic plane to convert from readout-system counts to power. The effective smoothing scale of the EBEX optics and the star camera-to-detector offset angles are determined through \c{hi}2 minimization using the compact HII region RCW 38. This two-step process is initially performed with parameters measured before the EBEX 2013 flight and then repeated until the calibration factor and parameters converge.

We present a follow-up study to the imaging polarimetry performed by Hayes et. al. 2011 on LAB1 in the SSA22 protocluster region. Arguably the most well-known Lyman-$\alpha$ "blob", this radio-quiet emission-line nebula likely hosts a galaxy which is either undergoing significant star formation or hosts an AGN, or both. We obtain deep, spatially resolved spectro-polarimetry of the Ly$\alpha$ emission and detect integrated linear polarization of $9$-$13\%\pm2$-$3\%$ at a distance of approximately 15 kpc north and south of the peak of the Lyman-$\alpha$ surface brightness with polarization vectors lying tangential to the galactic central source. In these same regions, we also detect a wavelength dependence in the polarization which is low at the center of the Ly$\alpha$ line profile and rises substantially in the wings of the profile. These polarization signatures are easily explained by a weak out-flowing shell model. The spectral dependence of the polarization presented here provide a framework for future observations and interpretations of the southern portion of LAB1 in that any model for this system must be able to reproduce this particular spectral dependence. However, questions still remain for the northern-most spur of LAB1. In this region we detect total linear polarization of between $3$ and $20\%$ at the $5\%$ significance level. Simulations predict that polarization should increase with radius for a symmetric geometry. That the northern spur does not suggests either that this region is not symmetric (which is likely) and exhibits variations in columns density, or that it is kinematically distinct from the rest of LAB1 and powered by another mechanism altogether.

We present mid- and far- IR imaging of four famous hypergiant stars: the red supergiants $\mu$ Cep and VY CMa, and the warm hypergiants IRC +10420 and $\rho$ Cas. Our 11 to 37 $\mu$m SOFIA/FORCAST imaging probes cool dust not detected in visual and near-IR imaging studies. Adaptive optics (AO) 8 - 12 $\mu$m imaging of $\mu$ Cep and IRC +10420 with MMT/MIRAC reveals extended envelopes that are the likely sources of these stars' strong silicate emission features. We find $\mu$ Cep's mass-loss rate to have declined by about a factor of 5 over a 13,000 history, ranging from 5 $\times$ 10$^{-6}$ down to $\sim$1 $\times$ 10$^{-6}$ $M_{\odot}$ yr$^{-1}$. The morphology of VY CMa indicates a cooler dust component coincident with the highly asymmetric reflection nebulae seen in the visual and near-IR. The lack of cold dust at greater distances around VY CMa indicates its mass-loss history is limited to the last $\sim$1200 years, with an average rate of 6 $\times$ 10$^{-4}$ $M_{\odot}$ yr$^{-1}$. We find two distinct periods in the mass-loss history of IRC +10420 with a high rate of 2 $\times$ 10$^{-3}$ $M_{\odot}$ yr$^{-1}$ until approximately 2000 yr ago, followed by an order of magnitude decrease in the recent past. We interpret this change as evidence of its evolution beyond the RSG stage. Our new infrared photometry of $\rho$ Cas is consistent with emission from the expanding dust shell ejected in its 1946 eruption, with no evidence of newer dust formation from its more recent events.

[abbreviated] We present a census of Ly\alpha\ emission at $z\gtrsim7$ utilizing deep near infrared HST grism spectroscopy from the first six completed clusters of the Grism Lens-Amplified Survey from Space (GLASS). In 24/159 photometrically selected galaxies we detect emission lines consistent with Ly\alpha\ in the GLASS spectra. Based on the distribution of signal-to-noise ratios and on simulations we expect the completeness and the purity of the sample to be 40-100% and 60-90%, respectively. For the objects without detected emission lines we show that the observed (not corrected for lensing magnification) 1$\sigma$ flux limits reaches $5\times10^{-18}$erg/s/cm$^{2}$ per position angle over the full wavelength range of GLASS (0.8-1.7$\mu$m). Based on the conditional probability of Ly\alpha\ emission measured from the ground at $z\sim7$ we would have expected 12-18 Ly\alpha\ emitters. This is consistent with the number of detections, within the uncertainties, confirming the drop in Ly\alpha\ emission with respect to $z\sim6$. These candidates include a promising source at $z=8.1$. The spatial extent of Ly\alpha\ in a deep stack of the most convincing Ly\alpha\ emitters with $\langle z\rangle=7.2$ is consistent with that of the rest-frame UV continuum. Extended Ly$\alpha$ emission, if present, has a surface brightness below our detection limit, consistent with the properties of lower redshift comparison samples. From the stack we estimate upper limits on rest-frame UV emission line ratios and find $f_\textrm{CIV} / f_\textrm{Ly${\alpha}$} \lesssim 0.32$ and $f_\textrm{CIII]} / f_\textrm{Ly$\alpha$} \lesssim 0.23$ in good agreement with other values published in the literature.

Radio relics are patches of diffuse synchrotron radio emission that trace shock waves. Relics are thought to form when intra-cluster medium electrons are accelerated by cluster merger induced shock waves through the diffusive shock acceleration mechanism. In this paper, we present observations spanning 150 MHz to 30 GHz of the `Sausage' and `Toothbrush' relics from the Giant Metrewave and Westerbork telescopes, the Karl G. Jansky Very Large Array, the Effelsberg telescope, the Arcminute Microkelvin Imager and Combined Array for Research in Millimeter-wave Astronomy. We detect both relics at 30 GHz, where the previous highest frequency detection was at 16 GHz. The integrated radio spectra of both sources clearly steepen above 2 GHz, at the >6$\sigma$ significance level, supports the spectral steepening previously found in the `Sausage' and the Abell 2256 relic. Our results challenge the widely adopted simple formation mechanism of radio relics and suggest more complicated models have to be developed that, for example, involve re-acceleration of aged seed electrons.

We give an overview of the Grism Lens Amplified Survey from Space (GLASS), a large Hubble Space Telescope program aimed at obtaining grism spectroscopy of the fields of ten massive clusters of galaxies at redshift z=0.308-0.686, including the Hubble Frontier Fields (HFF). The Wide Field Camera 3 yields near infrared spectra of the cluster cores, covering the wavelength range 0.81-1.69mum through grisms G102 and G141, while the Advanced Camera for Surveys in parallel mode provides G800L spectra of the infall regions of the clusters. The WFC3 spectra are taken at two almost orthogonal position angles in order to minimize the effects of confusion. After summarizing the scientific drivers of GLASS, we describe the sample selection as well as the observing strategy and data processing pipeline. We then utilize MACSJ0717.5+3745, a HFF cluster and the first one observed by GLASS, to illustrate the data quality and the high-level data products. Each spectrum brighter than H_AB=23 is visually inspected by at least two co-authors and a redshift is measured when sufficient information is present in the spectra. Furthermore, we conducted a thorough search for emission lines through all the GLASS WFC3 spectra with the aim of measuring redshifts for sources with continuum fainter than H_AB=23. We provide a catalog of 139 emission-line based spectroscopic redshifts for extragalactic sources, including three new redshifts of multiple image systems (one probable, two tentative). In addition to the data itself we also release software tools that are helpful to navigate the data.

We present the first direct calibration of strong-line metallicity diagnostics at significant cosmological distances using a sample at z=0.8 drawn from the DEEP2 Galaxy Redshift Survey. Oxygen and neon abundances are derived from measurements of electron temperature and density. We directly compare various commonly used relations between gas-phase metallicity and strong line ratios of O, Ne, and H at z=0.8 and z=0. There is no evolution with redshift at high precision ($\Delta \log{\mathrm{O/H}} = -0.01\pm0.03$, $\Delta \log{\mathrm{Ne/O}} = 0.01 \pm 0.01$). O, Ne, and H line ratios follow the same locus at z=0.8 as at z=0 with $\lesssim$0.02 dex evolution and low scatter ($\lesssim$0.04 dex). This suggests little or no evolution in physical conditions of HII regions at fixed oxygen abundance, in contrast to models which invoke more extreme properties at high redshifts. We speculate that offsets observed in the [N II]/H$\alpha$ versus [O III]/H$\beta$ diagram at high redshift are therefore due to [NII] emission, likely as a result of relatively high N/O abundance. If this is indeed the case, then nitrogen-based metallicity diagnostics suffer from systematic errors at high redshift. Our findings indicate that locally calibrated abundance diagnostics based on alpha-capture elements can be reliably applied at z$\simeq$1 and possibly at much higher redshifts. This constitutes the first firm basis for the widespread use of empirical calibrations in high redshift metallicity studies.

We present a strong and weak lensing reconstruction of the massive cluster Abell 2744, the first cluster for which deep Hubble Frontier Field (HFF) images and spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS) are available. By performing a targeted search for emission lines in multiply imaged sources using the GLASS spectra, we obtain 5 high-confidence spectroscopic redshifts and 2 tentative ones. We confirm 1 strongly lensed system by detecting the same emission lines in all 3 multiple images. We also search for additional line emitters blindly and use the full GLASS spectroscopic catalog to test reliability of photometric redshifts for faint line emitters. We see a reasonable agreement between our photometric and spectroscopic redshift measurements, when including nebular emission in photometric redshift estimations. We introduce a stringent procedure to identify only secure multiple image sets based on colors, morphology, and spectroscopy. By combining 7 multiple image systems with secure spectroscopic redshifts (at 5 distinct redshift planes) with 18 multiple image systems with secure photometric redshifts, we reconstruct the gravitational potential of the cluster pixellated on an adaptive grid, using a total of 72 images. The resulting mass map is compared with a stellar mass map obtained from the deep Spitzer Frontier Fields data to study the relative distribution of stars and dark matter in the cluster. We find that the stellar to total mass ratio varies substantially across the cluster field, suggesting that stars do not trace exactly the total mass in this interacting system. The maps of convergence, shear, and magnification are made available in the standard HFF format.

We analyze high resolution simulations of compressible, MHD turbulence with properties resembling conditions in galaxy clusters. The flow is driven to turbulence Mach number $\mathcal{M}_t \sim 1/2$ in an isothermal medium with an initially very weak, uniform seed magnetic field ($\beta = P_g/P_B = 10^6$). Since cluster turbulence is likely to result from a mix of sheared (solenoidal) and compressive forcing processes, we examine the distinct turbulence properties for both cases. In one set of simulations velocity forcing is entirely solenoidal ($\nabla\cdot \delta {\vec u} = 0$), while in the other it is entirely compressive ($\nabla\times \delta {\vec u} = 0$). Both cases develop a mixture of solenoidal and compressive turbulent motions, since each generates the other. The development of compressive turbulent motions leads to shocks, even when the turbulence is solenoidally forced and subsonic. Shocks, in turn, produce and amplify vorticity, which is especially important in compressively forced turbulence. To clarify those processes we include a pair of appendices that look in detail at vorticity evolution in association with shocks. From our simulation analyses we find that magnetic fields amplified to near saturation levels in predominantly solenoidal turbulence can actually enhance vorticity on small scales by concentrating and stabilizing shear. The properties, evolution rates and relative contributions of the kinetic and magnetic turbulent elements depend strongly on the character of the forcing. Specifically, shocks are stronger, but vorticity evolution and magnetic field amplification are slower and weaker when the turbulence is compressively forced. We identify a simple relation to estimate characteristic shock strengths in terms of the turbulence Mach number and the character of the forcing. Our results will be helpful in understanding flow motions in galaxy clusters.

We present results from the first twelve months of operation of Radio Galaxy Zoo, which upon completion will enable visual inspection of over 170,000 radio sources to determine the host galaxy of the radio emission and the radio morphology. Radio Galaxy Zoo uses $1.4\,$GHz radio images from both the Faint Images of the Radio Sky at Twenty Centimeters (FIRST) and the Australia Telescope Large Area Survey (ATLAS) in combination with mid-infrared images at $3.4\,\mu$m from the {\it Wide-field Infrared Survey Explorer} (WISE) and at $3.6\,\mu$m from the {\it Spitzer Space Telescope}. We present the early analysis of the WISE mid-infrared colours of the host galaxies. For images in which there is $>\,75\%$ consensus among the Radio Galaxy Zoo cross-identifications, the project participants are as effective as the science experts at identifying the host galaxies. The majority of the identified host galaxies reside in the mid-infrared colour space dominated by elliptical galaxies, quasi-stellar objects (QSOs), and luminous infrared radio galaxies (LIRGs). We also find a distinct population of Radio Galaxy Zoo host galaxies residing in a redder mid-infrared colour space consisting of star-forming galaxies and/or dust-enhanced non star-forming galaxies consistent with a scenario of merger-driven active galactic nuclei (AGN) formation. The completion of the full Radio Galaxy Zoo project will measure the relative populations of these hosts as a function of radio morphology and power while providing an avenue for the identification of rare and extreme radio structures. Currently, we are investigating candidates for radio galaxies with extreme morphologies, such as giant radio galaxies, late-type host galaxies with extended radio emission, and hybrid morphology radio sources.

We present J' and K' imaging linear polarimetric adaptive optics observations of NGC 1068 using MMT-Pol on the 6.5-m MMT. These observations allow us to study the torus from a magnetohydrodynamical (MHD) framework. In a 0.5" (30 pc) aperture at K', we find that polarisation arising from the passage of radiation from the inner edge of the torus through magnetically aligned dust grains in the clumps is the dominant polarisation mechanism, with an intrinsic polarisation of 7.0%$\pm$2.2%. This result yields a torus magnetic field strength in the range of 4$-$82 mG through paramagnetic alignment, and 139$^{+11}_{-20}$ mG through the Chandrasekhar-Fermi method. The measured position angle (P.A.) of polarisation at K$'$ is found to be similar to the P.A. of the obscuring dusty component at few parsec scales using infrared interferometric techniques. We show that the constant component of the magnetic field is responsible for the alignment of the dust grains, and aligned with the torus axis onto the plane of the sky. Adopting this magnetic field configuration and the physical conditions of the clumps in the MHD outflow wind model, we estimate a mass outflow rate $\le$0.17 M$_{\odot}$ yr$^{-1}$ at 0.4 pc from the central engine for those clumps showing near-infrared dichroism. The models used were able to create the torus in a timescale of $\geq$10$^{5}$ yr with a rotational velocity of $\leq$1228 km s$^{-1}$ at 0.4 pc. We conclude that the evolution, morphology and kinematics of the torus in NGC 1068 can be explained within a MHD framework.

The Keck science community is entering an era of unprecedented change. Powerful new instrument like ZTF, JWST, LSST, and the ELTs will catalyze this change, and we must be ready to take full advantage to maintain our position of scientific leadership. The best way to do this is to continue the UC and Caltech tradition of technical excellence in instrumentation. In this whitepaper we describe a new instrument called KRAKENS to help meet these challenges. KRAKENS uses a unique detector technology (MKIDs) to enable groundbreaking science across a wide range of astrophysical research topics. This document will lay out the detailed expected science return of KRAKENS.

We present 2 - 5 micron adaptive optics (AO) imaging and polarimetry of the famous hypergiant stars IRC +10420 and VY Canis Majoris. The imaging polarimetry of IRC +10420 with MMT-Pol at 2.2 micron resolves nebular emission with intrinsic polarization of 30%, with a high surface brightness indicating optically thick scattering. The relatively uniform distribution of this polarized emission both radially and azimuthally around the star confirms previous studies that place the scattering dust largely in the plane of the sky. Using constraints on scattered light consistent with the polarimetry at 2.2 micron, extrapolation to wavelengths in the 3 - 5 micron band predicts a scattered light component significantly below the nebular flux that is observed in our LBT/LMIRCam 3 - 5 micron AO imaging. Under the assumption this excess emission is thermal, we find a color temperature of ~ 500 K is required, well in excess of the emissivity-modified equilibrium temperature for typical astrophysical dust. The nebular features of VY CMa are found to be highly polarized (up to 60%) at 1.3 micron, again with optically thick scattering required to reproduce the observed surface brightness. This star's peculiar nebular feature dubbed the "Southwest Clump" is clearly detected in the 3.1 micron polarimetry as well, which, unlike IRC+10420, is consistent with scattered light alone. The high intrinsic polarizations of both hypergiants' nebulae are compatible with optically thick scattering for typical dust around evolved dusty stars, where the depolarizing effect of multiple scatters is mitigated by the grains' low albedos.

We present adaptive optics-assisted integral field spectroscopy around the Ha or Hb lines of 12 gravitationally lensed galaxies obtained with VLT/SINFONI, Keck/OSIRIS and Gemini/NIFS. We combine these data with previous observations and investigate the dynamics and star formation properties of 17 lensed galaxies at z = 1-4. Thanks to gravitational magnification of 1.4 - 90x by foreground clusters, effective spatial resolutions of 40 - 700 pc are achieved. The magnification also allows us to probe lower star formation rates and stellar masses than unlensed samples; our target galaxies feature dust-corrected SFRs derived from Ha or Hb emission of 0.8 - 40Msol/yr, and stellar masses M* ~ 4e8 - 6e10 Msol. All of the galaxies have velocity gradients, with 59% consistent with being rotating discs and a likely merger fraction of 29%, with the remaining 12% classed as 'undetermined.' We extract 50 star-forming clumps with sizes in the range 60pc - 1kpc from the Ha (or Hb) maps, and find that their surface brightnesses and their characteristic luminosities evolve to higher luminosities with redshift. We show that this evolution can be described by fragmentation on larger scales in gas-rich discs, and is likely to be driven by evolving gas fractions.

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.

Our objectives are to determine the properties of the interstellar medium (ISM) and of star-formation in typical star-forming galaxies at high redshift. Following up on our previous multi-wavelength observations with HST, Spitzer, Herschel, and the Plateau de Bure Interferometer (PdBI), we have studied a strongly lensed z=2.013 galaxy, the arc behind the galaxy cluster MACS J0451+0006, with ALMA to measure the [CII] 158 micron emission line, one of the main coolants of the ISM. [CII] emission from the southern part of this galaxy is detected at 10 $\sigma$. Taking into account strong gravitational lensing, which provides a magnification of $\mu=49$, the intrinsic lensing-corrected [CII]158 micron luminosity is $L(CII)=1.2 \times 10^8 L_\odot$. The observed ratio of [CII]-to-IR emission, $L(CII)/L(FIR) \approx (1.2-2.4) \times 10^{-3}$, is found to be similar to that in nearby galaxies. The same also holds for the observed ratio $L(CII)/L(CO)=2.3 \times 10^3$, which is comparable to that of star-forming galaxies and active galaxy nuclei (AGN) at low redshift. We utilize strong gravitational lensing to extend diagnostic studies of the cold ISM to an order of magnitude lower luminosity ($L(IR) \sim (1.1-1.3) \times 10^{11} L_\odot$) and SFR than previous work at high redshift. While larger samples are needed, our results provide evidence that the cold ISM of typical high redshift galaxies has physical characteristics similar to normal star forming galaxies in the local Universe.

We present spatially resolved gas-phase metallicity for a system of three galaxies at z=1.85 detected in the Grism Lens-Amplified Survey from Space (GLASS). The combination of HST's diffraction limit and strong gravitational lensing by the cluster MACS J0717+3745 results in a spatial resolution of ~200-300 pc, enabling good spatial sampling despite the intrinsically small galaxy sizes. The galaxies in this system are separated by 50-200 kpc in projection and are likely in an early stage of interaction, evidenced by relatively high specific star formation rates. Their gas-phase metallicities are consistent with larger samples at similar redshift, star formation rate, and stellar mass. We obtain a precise measurement of the metallicity gradient for one galaxy and find a shallow slope compared to isolated galaxies at high redshift, consistent with a flattening of the gradient due to gravitational interaction. An alternative explanation for the shallow metallicity gradient and elevated star formation rate is rapid recycling of metal-enriched gas, but we find no evidence for enhanced gas-phase metallicities which should result from this effect. Notably, the measured stellar masses log(M/Msun) = 7.2-9.1 probe to an order of magnitude below previous mass-metallicity studies at this redshift. The lowest mass galaxy has properties similar to those expected for Fornax at this redshift, indicating that GLASS is able to directly study the progenitors of local group dwarf galaxies on spatially resolved scales. Larger samples from the full GLASS survey will be ideal for studying the effects of feedback, and the time evolution of metallicity gradients. These initial results demonstrate the utility of HST spectroscopy combined with gravitational lensing for characterizing resolved physical properties of galaxies at high redshift.

For the Class 0 protostar, L1527, we compare 131 polarization vectors from SCUPOL/JCMT, SHARP/CSO and TADPOL/CARMA observations with the corresponding model polarization vectors of four ideal-MHD, non-turbulent, cloud core collapse models. These four models differ by their initial magnetic fields before collapse; two initially have aligned fields (strong and weak) and two initially have orthogonal fields (strong and weak) with respect to the rotation axis of the L1527 core. Only the initial weak orthogonal field model produces the observed circumstellar disk within L1527. This is a characteristic of nearly all ideal-MHD, non-turbulent, core collapse models. In this paper we test whether this weak orthogonal model also has the best agreement between its magnetic field structure and that inferred from the polarimetry observations of L1527. We found that this is not the case; based on the polarimetry observations the most favored model of the four is the weak aligned model. However, this model does not produce a circumstellar disk, so our result implies that a non-turbulent, ideal-MHD global collapse model probably does not represent the core collapse that has occurred in L1527. Our study also illustrates the importance of using polarization vectors covering a large area of a cloud core to determine the initial magnetic field orientation before collapse; the inner core magnetic field structure can be highly altered by a collapse and so measurements from this region alone can give unreliable estimates of the initial field configuration before collapse.

We present near infrared polarimetry data of background stars shining through a selection of starless cores taken in the $K$ band, probing visual extinctions up to $A_{V} \sim 48$. We find that $P_K/{\tau _K}$ continues to decline with increasing $A_{V}$ with a power law slope of roughly -0.5. Examination of published submillimeter (submm) polarimetry of starless cores suggests that by $A_{V} \gtrsim 20$ the slope for $P$ vs. $\tau$ becomes $\sim -1$, indicating no grain alignment at greater optical depths. Combining these two data sets, we find good evidence that, in the absence of a central illuminating source, the dust grains in dense molecular cloud cores with no internal radiation source cease to become aligned with the local magnetic field at optical depths greater than $A_V \sim 20$. A simple model relating the alignment efficiency to the optical depth into the cloud reproduces the observations well.

EBEX is a balloon-borne telescope designed to measure the polarization of the cosmic microwave background radiation. During its eleven day science flight in the Austral Summer of 2012, it operated 955 spider-web transition edge sensor (TES) bolometers separated into bands at 150, 250 and 410 GHz. This is the first time that an array of TES bolometers has been used on a balloon platform to conduct science observations. Polarization sensitivity was provided by a wire grid and continuously rotating half-wave plate. The balloon implementation of the bolometer array and readout electronics presented unique development requirements. Here we present an outline of the readout system, the remote tuning of the bolometers and Superconducting QUantum Interference Device (SQUID) amplifiers, and preliminary current noise of the bolometer array and readout system.

We present results from imaging of the radio filaments in the southern giant lobe of Centaurus A using data from GMRT observations at 325 and 235 MHz, and outcomes from filament modelling. The observations reveal a rich filamentary structure, largely matching the morphology at 1.4 GHz. We find no clear connection of the filaments to the jet. We seek to constrain the nature and origin of the vertex and vortex filaments associated with the lobe and their role in high-energy particle acceleration. We deduce that these filaments are at most mildly overpressured with respect to the global lobe plasma showing no evidence of large-scale efficient Fermi I-type particle acceleration, and persist for ~ 2-3 Myr. We demonstrate that the dwarf galaxy KK 196 (AM 1318-444) cannot account for the features, and that surface plasma instabilities, the internal sausage mode and radiative instabilities are highly unlikely. An internal tearing instability and the kink mode are allowed within the observational and growth time constraints and could develop in parallel on different physical scales. We interpret the origin of the vertex and vortex filaments in terms of weak shocks from transonic MHD turbulence or from a moderately recent jet activity of the parent AGN, or an interplay of both.

We consider a phenomenological model for the thermal leakage injection in the diffusive shock acceleration (DSA) process, in which suprathermal protons and electrons near the shock transition zone are assumed to have the so-called $\kappa$-distributions produced by interactions of background thermal particles with pre-existing and/or self-excited plasma/MHD waves or turbulence. The $\kappa$-distribution has a power-law tail, instead of an exponential cutoff, well above the thermal peak momentum. So there are a larger number of potential seed particles with momentum, above that required for participation in the DSA process. As a result, the injection fraction for the $\kappa$-distribution depends on the shock Mach number much less severely compared to that for the Maxwellian distribution. Thus, the existence of $\kappa$-like suprathermal tails at shocks would ease the problem of extremely low injection fractions, especially for electrons and especially at weak shocks such as those found in the intracluster medium. We suggest that the injection fraction for protons ranges $10^{-4}-10^{-3}$ for a $\kappa$-distribution with $10 < \kappa_p < 30$ at quasi-parallel shocks, while the injection fraction for electrons becomes $10^{-6}-10^{-5}$ for a $\kappa$-distribution with $\kappa_e < 2$ at quasi-perpendicular shocks. For such $\kappa$ values the ratio of cosmic ray electrons to protons naturally becomes $K_{e/p}\sim 10^{-3}-10^{-2}$, which is required to explain the observed ratio for Galactic cosmic rays.

Radio observations prove the existence of relativistic particles and magnetic field associated with the intra-cluster-medium (ICM) through the presence of extended synchrotron emission in the form of radio halos and peripheral relics. This observational evidence has fundamental implications on the physics of the ICM. Non-thermal components in galaxy clusters are indeed unique probes of very energetic processes operating within clusters that drain gravitational and electromagnetic energy into cosmic rays and magnetic fields. These components strongly affect the (micro-)physical properties of the ICM, including viscosity and electrical conductivities, and have also potential consequences on the evolution of clusters themselves. The nature and properties of cosmic rays in galaxy clusters, including the origin of the observed radio emission on cluster-scales, have triggered an active theoretical debate in the last decade. Only recently we can start addressing some of the most important questions in this field, thanks to recent observational advances, both in the radio and at high energies. The properties of cosmic rays and of cluster non-thermal emissions depend on the dynamical state of the ICM, the efficiency of particle acceleration mechanisms in the ICM and on the dynamics of these cosmic rays. In this review we discuss in some detail the acceleration and transport of cosmic rays in galaxy clusters and the most relevant observational milestones that have provided important steps on our understanding of this physics. Finally, looking forward to the possibilities from new generations of observational tools, we focus on what appear to be the most important prospects for the near future from radio and high-energy observations.

The Grism Lens-Amplified Survey from Space (GLASS) is a Hubble Space Telescope (HST) Large Program, which will obtain 140 orbits of grism spectroscopy of the core and infall regions of 10 galaxy clusters, selected to be among the very best cosmic telescopes. Extensive HST imaging is available from many sources including the CLASH and Frontier Fields programs. We introduce the survey by analyzing spectra of faint multiply-imaged galaxies and $z\gtrsim6$ galaxy candidates obtained from the first seven orbits out of fourteen targeting the core of the Frontier Fields cluster MACS0717.5+3745. Using the G102 and G141 grisms to cover the wavelength range 0.8-1.7$\mu$m, we confirm 4 strongly lensed systems by detecting emission lines in each of the images. For the 9 $z\gtrsim6$ galaxy candidates clear from contamination, we do not detect any emission lines down to a seven-orbit 1$\sigma$ noise level of $\sim$5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. Taking lensing magnification into account, our flux sensitivity reaches $\sim$0.2-5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. These limits over an uninterrupted wavelength range rule out the possibility that the high-$z$ galaxy candidates are instead strong line emitters at lower redshift. These results show that by means of careful modeling of the background - and with the assistance of lensing magnification - interesting flux limits can be reached for large numbers of objects, avoiding pre-selection and the wavelength restrictions inherent to ground-based multi-slit spectroscopy. These observations confirm the power of slitless HST spectroscopy even in fields as crowded as a cluster core.

We present Submillimeter Array (SMA) [CII] 158um and Jansky Very Large Array (JVLA) $^{12}$CO(1-0) line emission maps for the bright, lensed, submillimeter source at $z=5.2430$ behind Abell 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple $^{12}$CO rotational transitions, [CI], water and [NII], providing some of the best constraints on the properties of the interstellar medium (ISM) in a galaxy at $z>5$. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8-1000um) = (1.6$\pm$0.1)x10^14 L_sun/mu, where the total magnification mu_total = 8.9$\pm$1.9, via a new lens model from the [CII] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components which correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source plane separation <1 kpc. The reddest dominates the continuum map (de-magnified L_FIR = (1.1$\pm$0.2)x10^13 L_sun), and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (dV ~ 500 km/s) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region, and photo-dissociation region (PDR) analysis suggests that it is lower density, cooler and forming stars less vigorously than the other components. Strikingly, it has very strong [NII] emission which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity.

We examine how the total X-ray luminosity correlates with stellar mass, stellar population, and redshift for a K-band limited sample of ~3500 galaxies at 0.5<z<2.0 from the NEWFIRM Medium Band Survey in the COSMOS field. The galaxy sample is divided into 32 different galaxy types, based on similarities between the spectral energy distributions. For each galaxy type, we further divide the sample into bins of redshift and stellar mass, and perform an X-ray stacking analysis using the Chandra COSMOS (C-COSMOS) data. We find that full band X-ray luminosity is primarily increasing with stellar mass, and at similar mass and spectral type is higher at larger redshifts. When comparing at the same stellar mass, we find that the X-ray luminosity is slightly higher for younger galaxies (i.e., weaker 4000\AA breaks), but the scatter in this relation is large. We compare the observed X-ray luminosities to those expected from low and high mass X-ray binaries (XRBs). For blue galaxies, XRBs can almost fully account for the observed emission, while for older galaxies with larger 4000\AA breaks, active galactic nuclei (AGN) or hot gas dominate the measured X-ray flux. After correcting for XRBs, the X-ray luminosity is still slightly higher in younger galaxies, although this correlation is not significant. AGN appear to be a larger component of galaxy X-ray luminosity at earlier times, as the hardness ratio increases with redshift. Together with the slight increase in X-ray luminosity this may indicate more obscured AGNs or higher accretion rates at earlier times.

We explore nonlinear effects of wave-particle interactions on the diffusive shock acceleration (DSA) process in Type Ia-like, SNR blast waves, by implementing phenomenological models for magnetic field amplification, Alfv'enic drift, and particle escape in time-dependent numerical simulations of nonlinear DSA. For typical SNR parameters the CR protons can be accelerated to PeV energies only if the region of amplified field ahead of the shock is extensive enough to contain the diffusion lengths of the particles of interest. Even with the help of Alfv'enic drift, it remains somewhat challenging to construct a nonlinear DSA model for SNRs in which order of 10 % of the supernova explosion energy is converted to the CR energy and the magnetic field is amplified by a factor of 10 or so in the shock precursor, while, at the same time, the energy spectrum of PeV protons is steeper than E^{-2}. To explore the influence of these physical effects on observed SNR emissions, we also compute resulting radio-to-gamma-ray spectra. Nonthermal emission spectra, especially in X-ray and gamma-ray bands,depend on the time dependent evolution of CR injection process, magnetic field amplification, and particle escape, as well as the shock dynamic evolution. This result comes from the fact that the high energy end of the CR spectrum is composed of the particles that are injected in the very early stages of blast wave evolution. Thus it is crucial to understand better the plasma wave-particle interactions associated with collisionless shocks in detail modeling of nonthermal radiation from SNRs.

We present adaptive optics images of the extreme red supergiant VY Canis Majoris in the Ks, L' and M bands (2.15 to 4.8 micron) made with LMIRCam on the Large Binocular Telescope (LBT). The peculiar "Southwest Clump" previously imaged from 1 to 2.2 micron appears prominently in all three filters. We find its brightness is due almost entirely to scattering, with the contribution of thermal emission limited to at most 25%. We model its brightness as optically thick scattering from silicate dust grains using typical size distributions. We find a lower limit mass for this single feature of 5E-03 Msun to 2.5E-02 Msun depending on the assumed gas-to-dust ratio. The presence of the Clump as a distinct feature with no apparent counterpart on the other side of the star is suggestive of an ejection event from a localized region of the star and is consistent with VY CMa's history of asymmetric high mass loss events.

We present a census of z(abs) < 2, intrinsic (those showing partial coverage) and associated [z(abs) ~ z(em)] quasar absorption-line systems detected in the Hubble Space Telescope archive of Space Telescope Imaging Spectrograph echelle spectra. This work complements the Misawa et al. (2007) survey of 2 < z(em) < 4 quasars that selects systems using similar techniques. We confirm the existence of so-called "strong N V" intrinsic systems (where the equivalent width of H I Ly alpha is small compared to N V 1238) presented in that work, but find no convincing cases of "strong C IV" intrinsic systems at low redshift/luminosity. Moreover, we also report on the existence of "strong O VI" systems. From a comparison of partial coverage results as a function of ion, we conclude that systems selected by the N V ion have the highest probability of being intrinsic. By contrast, the C IV and O VI ions are poor selectors. Of the 30 O VI systems tested, only two of the systems in the spectrum on 3C 351 show convincing evidence for partial coverage. However, there is a 3-sigma excess in the number of absorbers near the quasar redshift (|Delta v| <= 5000 km/s) over absorbers at large redshift differences. In at least two cases, the associated O VI systems are known not to arise close to the accretion disk of the quasar.

We present rest-frame optical spectra for a sample of 9 low-mass star-forming galaxies in the redshift range 1.5 < z < 3 which are gravitationally lensed by foreground clusters. We used Triplespec, an echelle spectrograph at the Palomar 200-inch telescope that is very effective for this purpose as it samples the entire near-infrared spectrum simultaneously. By measuring the flux of nebular emission lines we derive gas phase metallicities and star formation rates, and by fitting the optical to infrared spectral energy distributions we obtain stellar masses. Taking advantage of the high magnification due to strong lensing we are able to probe the physical properties of galaxies with stellar masses in the range 7.8 < log M/Msun < 9.4 whose star formation rates are similar to those of typical star-forming galaxies in the local universe. We compare our results with the locally determined relation between stellar mass, gas metallicity and star formation rate. Our data are in excellent agreement with this relation, with an average offset <Delta log O/H> = 0.01 +/- 0.08, suggesting a universal relationship. Remarkably, the scatter around the fundamental metallicity relation is only 0.24 dex, smaller than that observed locally at the same stellar masses, which may provide an important additional constraint for galaxy evolution models.

The fraction of ionizing photons that escape from young star-forming galaxies is one of the largest uncertainties in determining the role of galaxies in cosmic reionization. Yet traditional techniques for measuring this fraction are inapplicable at the redshifts of interest due to foreground screening by the Lyman alpha forest. In an earlier study, we demonstrated a reduction in the equivalent width of low-ionization absorption lines in composite spectra of Lyman break galaxies at z=4 compared to similar measures at z=3. This might imply a lower covering fraction of neutral gas and hence an increase with redshift in the escape fraction of ionizing photons. However, our spectral resolution was inadequate to differentiate between several alternative explanations, including changes with redshift in the outflow kinematics. Here we present higher quality spectra of 3 gravitationally lensed Lyman break galaxies at z=4 with a spectral resolution sufficient to break this degeneracy of interpretation. We present a method for deriving the covering fraction of low-ionization gas as a function of outflow velocity and compare the results with similar quality data taken for galaxies at lower redshift. We find a significant trend of lower covering fractions of low-ionization gas for galaxies with strong \Lya emission. In combination with the demographic trends of \Lya emission with redshift from our earlier work, our results provide new evidence for a reduction in the average H I covering fraction, and hence an increase in the escape fraction of ionizing radiation from Lyman break galaxies, with redshift.

An optically and geometrically thick torus obscures the central engine of Active Galactic Nuclei (AGN) from some lines of sight. From a magnetohydrodynamical framework, the torus can be considered to be a particular region of clouds surrounding the central engine where the clouds are dusty and optically thick. In this framework, the magnetic field plays an important role in the creation, morphology and evolution of the torus. If the dust grains within the clouds are assumed to be aligned by paramagnetic alignment, then the ratio of the intrinsic polarisation and visual extinction, P(%)/Av, is a function of the magnetic field strength. To estimate the visual extinction through the torus and constrain the polarisation mechanisms in the nucleus of AGN, we developed a polarisation model to fit both the total and polarised flux in a 1.2" (~263pc) aperture of the type 2 AGN, IC5063. The polarisation model is consistent with the nuclear polarisation observed at K being produced by dichroic absorption from aligned dust grains with a visual extinction through the torus of 48$\pm$2 mag. We estimated the intrinsic polarisation arising from dichroic absorption to be P$_{K}^{dic}$=12.5$\pm$2.7%. We consider the physical conditions and environment of the gas and dust for the torus of IC5063. Then, through paramagnetic alignment, we estimate a magnetic field strength in the range of 12-128mG in the NIR emitting regions of the torus of IC5063. Alternatively, we estimate the magnetic field strength in the plane of the sky using the Chandrasekhar-Fermi method. The minimum magnetic field strength in the plane of the sky is estimated to be 13 and 41 mG depending of the conditions within the torus of IC5063. These techniques afford the chance to make a survey of AGN, to investigate the effects of magnetic field strength on the torus, accretion, and interaction to the host galaxy.

Bright gravitationally lensed galaxies provide our most detailed view of galaxies at high redshift. Yet as a result of the small number of ultra-bright z~2 lensed systems with confirmed redshifts, most detailed spectroscopic studies have been limited in their scope. With the goal of increasing the number of bright lensed galaxies available for detailed follow-up, we have undertaken a spectroscopic campaign targeting wide separation (>3 arcsec) galaxy-galaxy lens candidates within the Sloan Digital Sky Survey (SDSS). Building on the earlier efforts of our CASSOWARY survey, we target a large sample of candidate galaxy-galaxy lens systems in SDSS using a well-established search algorithm which identifies blue arc-like structures situated around luminous red galaxies. In this paper, we present a new redshift catalog containing 25 lensed sources in SDSS confirmed through spectroscopic follow-up of candidate galaxy-galaxy lens systems. Included in this new sample are two of the brightest galaxies (r=19.6 and 19.7) galaxies known at z~2, a low metallicity (12 + log (O/H)~8.0) extreme nebular line emitting galaxy at z=1.43, and numerous systems for which detailed follow-up will be possible. The source redshifts span 0.9<z<2.5 (median redshift of 1.9), and their optical magnitudes are in the range 19.6<r<22.3. We present a brief source-by-source discussion of the spectroscopic properties extracted from our confirmatory spectra and discuss some initial science results. With more than 50 gravitationally lensed z>1 galaxies now confirmed within SDSS, it will soon be possible for the first time to develop generalized conclusions from detailed spectroscopic studies of the brightest lensed systems at high redshift.

We present a 3-5um LBT/MMT adaptive optics imaging study of three Upper Scorpius stars with brown dwarf (BD) companions with very low-masses/mass ratios (M_BD < 25M_Jup; M_BD / M_star ~ 1-2%), and wide separations (300-700 AU): GSC 06214, 1RXS 1609, and HIP 78530. We combine these new thermal IR data with existing 1-4um and 24um photometry to constrain the properties of the BDs and identify evidence for circumprimary/secondary disks in these unusual systems. We confirm that GSC 06214B is surrounded by a disk, further showing this disk produces a broadband IR excess due to small dust near the dust sublimation radius. An unresolved 24um excess in the system may be explained by the contribution from this disk. 1RXS 1609B exhibits no 3-4um excess, nor does its primary; however, the system as a whole has a modest 24um excess, which may come from warm dust around the primary and/or BD. Neither object in the HIP 78530 system exhibits near- to mid-IR excesses. We additionally find that the 1-4um colors of HIP 78530B match a spectral type of M3+-2, inconsistent with the M8 spectral type assigned based on its near-IR spectrum, indicating it may be a low-mass star rather than a BD. We present new upper limits on additional low-mass companions in the system (<5M_Jup beyond 175AU). Finally, we examine the utility of circumsecondary disks as probes of the formation histories of wide BD companions, finding that the presence of a disk may disfavor BD formation near the primary with subsequent outward scattering.

SN2011ht has been described both as a true supernova and as an impostor. In this paper, we conclude that it does not match some basic expectations for a core-collapse event. We discuss SN2011ht's spectral evolution from a hot dense wind to a cool dense wind, followed by the post-plateau appearance of a faster low density wind during a rapid decline in luminosity. We identify a slow dense wind expanding at only 500--600 km/s, present throughout the eruption. A faster wind speed V ~ 900 km/s may be identified with a second phase of the outburst. There is no direct or significant evidence for any flow speed above 1000 km/s; the broad asymmetric wings of Balmer emission lines in the hot wind phase were due to Thomson scattering, not bulk motion. We estimate a mass loss rate of order 0.04 Msun/yr during the hot dense wind phase of the event. There is no evidence that the kinetic energy substantially exceeded the luminous energy, roughly 2 X 10^49 ergs; so the total energy was far less than a true SN. We suggest that SN2011ht was a giant eruption driven by super-Eddington radiation pressure, perhaps beginning about 6 months before the discovery. A strongly non-spherical SN might also account for the data, at the cost of more free parameters.

We present HST/WFC3 narrowband imaging of the H-alpha emission in a sample of eight gravitationally-lensed galaxies at z = 1 - 1.5. The magnification caused by the foreground clusters enables us to obtain a median source plane spatial resolution of 360pc, as well as providing magnifications in flux ranging from ~10x to ~50x. This enables us to identify resolved star-forming HII regions at this epoch and therefore study their H-alpha luminosity distributions for comparisons with equivalent samples at z ~ 2 and in the local Universe. We find evolution in the both luminosity and surface brightness of HII regions with redshift. The distribution of clump properties can be quantified with an HII region luminosity function, which can be fit by a power law with an exponential break at some cut-off, and we find that the cut-off evolves with redshift. We therefore conclude that `clumpy' galaxies are seen at high redshift because of the evolution of the cut-off mass; the galaxies themselves follow similar scaling relations to those at z = 0, but their HII regions are larger and brighter and thus appear as clumps which dominate the morphology of the galaxy. A simple theoretical argument based on gas collapsing on scales of the Jeans mass in a marginally unstable disk shows that the clumpy morphologies of high-z galaxies are driven by the competing effects of higher gas fractions causing perturbations on larger scales, partially compensated by higher epicyclic frequencies which stabilise the disk.

We present a study of the evolution of several classes of MgII absorbers, and their corresponding FeII absorption, over a large fraction of cosmic history: 2.3 to 8.7 Gyrs from the Big Bang. Our sample consists of 87 strong (Wr(MgII)>0.3 A) MgII absorbers, with redshifts 0.2<z<2.5, measured in 81 quasar spectra obtained from the Very Large Telescope(VLT)/Ultraviolet and Visual Echelle Spectrograph(UVES) archives of high-resolution spectra (R \sim 45,000). No evolutionary trend in Wr(FeII)/Wr(MgII) is found for moderately strong MgII absorbers (0.3<Wr(MgII)<1.0 A). However, at lower z we find an absence of very strong MgII absorbers (those with Wr(MgII)>1 A) with small ratios of equivalent widths of FeII to MgII. At high z, very strong MgII absorbers with both small and large Wr(FeII)/Wr(MgII) values are present. We compare our findings to a sample of 100 weak MgII absorbers (Wr(MgII)<0.3 A) found in the same quasar spectra by Narayanan et al. (2007). The main effect driving the evolution of very strong MgII systems is the difference between the kinematic profiles at low and high redshifts. At high z, we observe that, among the very strong MgII absorbers, all of the systems with small ratios of Wr(FeII)/Wr(MgII) have relatively large velocity spreads, resulting in less saturated profiles. At low z, such kinematically spread systems are absent, and both FeII and MgII are saturated, leading to Wr(FeII)/Wr(MgII) values that are all close to 1. The high redshift, small Wr(FeII)/Wr(MgII) systems could correspond to sub-DLA systems, many of which have large velocity spreads and are possibly linked to superwinds in star forming galaxies. In addition to the change in saturation due to kinematic evolution, the smaller Wr(FeII)/Wr(MgII) values could be due to a lower abundance of Fe at high z, which would indicate relatively early stages of star formation in those environments.

We present and discuss measurements of the gas-phase metallicity gradient in gravitationally lensed galaxies at z=2.0-2.4 based on adaptive optics-assisted imaging spectroscopy with the Keck II telescope. Through deep exposures we have secured high signal to noise data for four galaxies with well-understood kinematic properties. Three galaxies with well-ordered rotation reveal metallicity gradients in the sense of having lower gas-phase metallicities at larger galactocentric radii. Two of these display gradients much steeper than found locally, while a third has one similar to that seen in local disk galaxies. The fourth galaxy exhibits complex kinematics indicative of an ongoing merger and reveals an "inverted" gradient with lower metallicity in the central regions. By comparing our sample to similar data in the literature for lower redshift galaxies, we determine that, on average, metallicity gradients must flatten by a factor of 2.6 +/- 0.9 between z=2.2 and the present epoch. This factor is in rough agreement with the size growth of massive galaxies suggesting that inside-out growth can account for the evolution of metallicity gradients. Since the addition of our new data provides the first indication of a coherent picture of this evolution, we develop a simple model of chemical evolution to explain the collective data. We find that metallicity gradients and their evolution can be explained by the inward radial migration of gas together with a radial variation in the mass loading factor governing the ratio of outflowing gas to the local star formation rate. Average mass loading factors of \lsim 2 are inferred from our model in good agreement with direct measurements of outflowing gas in z \simeq 2 galaxies.

We present diffraction-limited \ks band and \lprime adaptive optics images of the edge-on debris disk around the nearby F2 star HD 15115, obtained with a single 8.4 m primary mirror at the Large Binocular Telescope. At \ks band the disk is detected at signal-to-noise per resolution element (SNRE) \about 3-8 from \about 1-2\fasec 5 (45-113 AU) on the western side, and from \about 1.2-2\fasec 1 (63-90 AU) on the east. At \lprime the disk is detected at SNRE \about 2.5 from \about 1-1\fasec 45 (45-90 AU) on both sides, implying more symmetric disk structure at 3.8 \microns . At both wavelengths the disk has a bow-like shape and is offset from the star to the north by a few AU. A surface brightness asymmetry exists between the two sides of the disk at \ks band, but not at \lprime . The surface brightness at \ks band declines inside 1\asec (\about 45 AU), which may be indicative of a gap in the disk near 1\asec. The \ks - \lprime disk color, after removal of the stellar color, is mostly grey for both sides of the disk. This suggests that scattered light is coming from large dust grains, with 3-10 \microns -sized grains on the east side and 1-10 \microns dust grains on the west. This may suggest that the west side is composed of smaller dust grains than the east side, which would support the interpretation that the disk is being dynamically affected by interactions with the local interstellar medium.

Recent radio observations have identified a class of structures, so-called radio relics, in clusters of galaxies. The radio emission from these sources is interpreted as synchrotron radiation from GeV electrons gyrating in microG-level magnetic fields. Radio relics, located mostly in the outskirts of clusters, seem to associate with shock waves, especially those developed during mergers. In fact, they seem to be good structures to identify and probe such shocks in intracluster media (ICMs), provided we understand the electron acceleration and re-acceleration at those shocks. In this paper, we describe time-dependent simulations for diffusive shock acceleration at weak shocks that are expected to be found in ICMs. Freshly injected as well as pre-existing populations of cosmic-ray (CR) electrons are considered, and energy losses via synchrotron and inverse Compton are included. We then compare the synchrotron flux and spectral distributions estimated from the simulations with those in two well-observed radio relics in CIZA J2242.8+5301 and ZwCl0008.8+5215. Considering that the CR electron injection is rather inefficient at weak shocks with Mach number M <~ a few, the existence of radio relics could indicate the pre-existing population of low-energy CR electrons in ICMs. The implication of our results on the merger shock scenario of radio relics is discussed.

As the only directly imaged multiple planet system, HR 8799 provides a unique opportunity to study the physical properties of several planets in parallel. In this paper, we image all four of the HR 8799 planets at H-band and 3.3 microns with the new LBT adaptive optics system, PISCES, and LBTI/LMIRCam. Our images offer an unprecedented view of the system, allowing us to obtain H and 3.3$ micron photometry of the innermost planet (for the first time) and put strong upper-limits on the presence of a hypothetical fifth companion. We find that all four planets are unexpectedly bright at 3.3 microns compared to the equilibrium chemistry models used for field brown dwarfs, which predict that planets should be faint at 3.3 microns due to CH4 opacity. We attempt to model the planets with thick-cloudy, non-equilibrium chemistry atmospheres, but find that removing CH4 to fit the 3.3 micron photometry increases the predicted L' (3.8 microns) flux enough that it is inconsistent with observations. In an effort to fit the SED of the HR 8799 planets, we construct mixtures of cloudy atmospheres, which are intended to represent planets covered by clouds of varying opacity. In this scenario, regions with low opacity look hot and bright, while regions with high opacity look faint, similar to the patchy cloud structures on Jupiter and L/T transition brown-dwarfs. Our mixed cloud models reproduce all of the available data, but self-consistent models are still necessary to demonstrate their viability.

We present a pair of 3-d magnetohydrodynamical simulations of intermittent jets from a central active galactic nucleus (AGN) in a galaxy cluster extracted from a high resolution cosmological simulation. The selected cluster was chosen as an apparently relatively relaxed system, not having undergone a major merger in almost 7 Gyr. Despite this characterization and history, the intra-cluster medium (ICM) contains quite active "weather". We explore the effects of this ICM weather on the morphological evolution of the AGN jets and lobes. The orientation of the jets is different in the two simulations so that they probe different aspects of the ICM structure and dynamics. We find that even for this cluster that can be characterized as relaxed by an observational standard, the large-scale, bulk ICM motions can significantly distort the jets and lobes. Synthetic X-ray observations of the simulations show that the jets produce complex cavity systems, while synthetic radio observations reveal bending of the jets and lobes similar to wide-angle tail (WAT) radio sources. The jets are cycled on and off with a 26 Myr period using a 50% duty cycle. This leads to morphological features similar to those in "double-double" radio galaxies. While the jet and ICM magnetic fields are generally too weak in the simulations to play a major role in the dynamics, Maxwell stresses can still become locally significant.

We present and discuss the mean rest-frame ultraviolet spectrum for a sample of 81 Lyman Break Galaxies (LBGs) selected to be B-band dropouts with a mean redshift of z=3.9 and apparent magnitudes z_AB<26. Most of the individual spectra are drawn from our ongoing survey in the GOODS fields with the Keck DEIMOS spectrograph, and we have augmented our sample with published data taken with FORS2 on the VLT. In general we find similar trends in the spectral diagnostics to those found in the earlier, more extensive survey of LBGs at z=3 undertaken by Shapley et al (2003). Specifically, we find low-ionization absorption lines which trace the presence of neutral outflowing gas are weaker in galaxies with stronger Lyman-alpha emission, bluer UV spectral slopes, lower stellar masses, lower UV luminosities and smaller half-light radii. This is consistent with a physical picture whereby star formation drives outflows of neutral gas which scatters Lyman-alpha and gives rise to strong low-ionization absorption lines, while increasing the stellar mass, size, metallicity, and dust content of galaxies. Typical galaxies are thus expected to have stronger Lyman-alpha emission and weaker low-ionization absorption at earlier times (higher redshifts). Indeed, our mean spectrum at z=4 shows somewhat weaker low-ionization absorption lines than at z=3 and available data at higher redshift indicates a rapid decrease in low-ionization absorption strength with redshift. We argue that the reduced low-ionization absorption is likely caused by a decrease in the covering fraction and/or velocity range of outflowing neutral gas at earlier epochs. Our continuing survey will enable us to extend these diagnostics more reliably to higher redshift and determine the implications for the escape fraction of ionizing photons which governs the role of early galaxies in cosmic reionization. [Abridged]

New visible and K-band polarization measurements on stars surrounding molecular clouds in Orion A and stars in the BN vicinity are presented. Our results confirm that magnetic fields located inside the Orion A molecular clouds and in their close neighborhood are spatially connected. On and around the BN object, we measured the angular offsets between the K-band polarization data and available submm data. We find high values of the polarization degree, P_{K}, and of the optical depth, \tau_{K}, close to an angular offset position of 90^{\circ} whereas lower values of P_{K} and \tau_{K} are observed for smaller angular offsets. We interpret these results as evidence for the presence of various magnetic field components toward lines of sight in the vicinity of BN. On a larger scale, we measured the distribution of angular offsets between available H-band polarization data and the same submm data set. Here we find an increase of <P_{H}> with angular offset which we interpret as a rotation of the magnetic field by \lesssim 60^{\circ}. This trend generalizes previous results on small scale toward and around lines of sight to BN and is consistent with a twist of the magnetic field on a larger scale towards OMC-1. A comparison of our results with several other studies suggests that a two-component magnetic field, maybe helical, could be wrapping the OMC-1 filament.

The gaseous media of galaxy clusters and cosmic filaments, which constitute most of the baryonic matter in the universe, is highly dynamic. It is also probably turbulent, although the turbulence properties are poorly known. The gas is highly rarefied, essentially fully ionized plasma. Observational evidence suggests intracluster media (ICMs) are magnetized at some level. There are several possible origins for ICM seed fields; the observed fields are likely the result of turbulence in the ICM. We are engaged in a simulation study designed to understand in this context how very weak initial magnetic fields evolve in driven turbulence. We find that the magnetic fields eventually evolve towards equipartition levels with the vortical, solenoidal kinetic energy in the turbulence. As they do so the topology of the field structures transition from filamentary forms into ribbon-like structures in which the field orientations are laminated with vorticity structures.

We present calculations of expected continuum emissions from Sedov-Taylor phase Type Ia supernova remnants (SNRs), using the energy spectra of cosmic ray (CR) electrons and protons from nonlinear diffusive shock acceleration (DSA) simulations. A new, general-purpose radiative process code, Cosmicp, was employed to calculate the radiation expected from CR electrons and protons and their secondary products. These radio, X-ray and gamma-ray emissions are generally consistent with current observations of Type Ia SNRs. The emissions from electrons in these models dominate the radio through X-ray bands. Decays of \pi^0 s from p-p collisions mostly dominate the gamma-ray range, although for a hot, low density ISM case (n_{ISM}=0.003 cm^{-3}), the pion decay contribution is reduced sufficiently to reveal the inverse Compton contribution to TeV gamma-rays. In addition, we present simple scalings for the contributing emission processes to allow a crude exploration of model parameter space, enabling these results to be used more broadly. We also discuss the radial surface brightness profiles expected for these model SNRs in the X-ray and gamma-ray bands.

We report optical imaging, optical and near-infrared polarimetry, and Spitzer mid-infrared spectroscopy of comet C/2007 N3 (Lulin). Polarimetric observations were obtained in R (0.676 micron) at phase angles from 0.44 degrees to 21 degrees with simultaneous observations in H (1.65 micron) at 4.0 degrees, exploring the negative branch in polarization. Comet C/2007 N3 (Lulin) shows typical negative polarization in the optical as well as a similar negative branch near-infrared wavelengths. The 10 micron silicate feature is only weakly in emission and according to our thermal models, is consistent with emission from a mixture of silicate and carbon material. We argue that large, low-porosity (akin to Ballistic Particle Cluster Aggregates) rather absorbing aggregate dust particles best explain both the polarimetric and the mid-infrared spectral energy distribution.

In this paper we present our system design and methodology for making absolute quantum efficiency (QE) measurements through the vacuum ultraviolet (VUV) and verify the system with delta-doped silicon CCDs. Delta-doped detectors provide an excellent platform to validate measurements through the VUV due to their enhanced UV response. The requirements for measuring QE through the VUV are more strenuous than measurements in the near UV and necessitate, among other things, the use of a vacuum monochromator, good dewar chamber vacuum to prevent on-chip condensation, and more stringent handling requirements.

We have used Molecular Beam Epitaxy (MBE)-based delta doping technology to demonstrate near 100% internal quantum efficiency (QE) on silicon electron-multiplied Charge Coupled Devices (EMCCDs) for single photon counting detection applications. Furthermore, we have used precision techniques for depositing antireflection (AR) coatings by employing Atomic Layer Deposition (ALD) and demonstrated over 50% external QE in the far and near-ultraviolet in megapixel arrays. We have demonstrated that other device parameters such as dark current are unchanged after these processes. In this paper, we report on these results and briefly discuss the techniques and processes employed.

We have studied the properties of giant star forming clumps in five z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. The clumps reside in disk regions where the Toomre Q-parameter is below unity, consistent with their being bound and having formed from gravitational instability. Broad H{\alpha}/[NII] line wings demonstrate that the clumps are launching sites of powerful outflows. The inferred outflow rates are comparable to or exceed the star formation rates, in one case by a factor of eight. Typical clumps may lose a fraction of their original gas by feedback in a few hundred million years, allowing them to migrate into the center. The most active clumps may lose much of their mass and disrupt in the disk. The clumps leave a modest imprint on the gas kinematics. Velocity gradients across the clumps are 10-40 km/s/kpc, similar to the galactic rotation gradients. Given beam smearing and clump sizes, these gradients may be consistent with significant rotational support in typical clumps. Extreme clumps may not be rotationally supported; either they are not virialized, or they are predominantly pressure supported. The velocity dispersion is spatially rather constant and increases only weakly with star formation surface density. The large velocity dispersions may be driven by the release of gravitational energy, either at the outer disk/accreting streams interface, and/or by the clump migration within the disk. Spatial variations in the inferred gas phase oxygen abundance are broadly consistent with inside-out growing disks, and/or with inward migration of the clumps.

We present a consistent 3D model for the head-tail radio galaxy NGC 1265 that explains the complex radio morphology and spectrum by a past passage of the galaxy and radio bubble through a shock wave. Using analytical solutions to the full Riemann problem and hydrodynamical simulations, we study how this passage transformed the plasma bubble into a toroidal vortex ring. Adiabatic compression of the aged electron population causes it to be energized and to emit low-surface brightness and steep-spectrum radio emission. The large infall velocity of NGC 1265 and the low Faraday rotation measure values and variance of the jet strongly argue that this transformation was due to the accretion shock onto Perseus situated roughly at R_200. Estimating the volume change of the radio bubble enables inferring a shock Mach number of M = 4.2_{-1.2}^{+0.8}, a density jump of 3.4_{-0.4}^{+0.2}, a temperature jump of 6.3_{-2.7}^{+2.5}, and a pressure jump of 21.5 +/- 10.5 while allowing for uncertainties in the equation of state of the radio plasma and volume of the torus. Extrapolating X-ray profiles, we obtain upper limits on the gas temperature and density in the infalling warm-hot intergalactic medium of kT < 0.4 keV and n < 5e-5 / cm^3. The orientation of the ellipsoidally shaped radio torus in combination with the direction of the galaxy's head and tail in the plane of the sky is impossible to reconcile with projection effects. Instead, this argues for post-shock shear flows that have been caused by curvature in the shock surface with a characteristic radius of 850 kpc. The energy density of the shear flow corresponds to a turbulent-to-thermal energy density of 14%. The shock-injected vorticity might be important in generating and amplifying magnetic fields in galaxy clusters. Future LOFAR observations of head-tail galaxies can be complementary probes of accretion shocks onto galaxy clusters.

Observations of X-ray cavities formed by powerful jets from AGN in galaxy cluster cores are widely used to estimate the energy output of the AGN. Using methods commonly applied to observations of clusters, we conduct synthetic X-ray observations of 3D MHD simulated jet-ICM interactions to test the reliability of measuring X-ray cavity power. These measurements are derived from empirical estimates of the enthalpy content of the cavities and their implicit ages. We explore how such physical factors as jet intermittency and observational conditions such as orientation of the jets with respect to the line of sight impact the reliability of observational measurements of cavity enthalpy and age. An estimate of the errors in these quantities can be made by directly comparing "observationally" derived values with "actual" values from the simulations. In our tests, cavity enthalpy derived from observations was typically within a factor of two of the simulation values. Cavity age and, therefore, cavity power are sensitive to the accuracy of the estimated inclination angle of the jets. Cavity age and power estimates within a factor of two of the actual values are possible given an accurate inclination angle.

Cluster media are dynamical, not static; observational evidence suggests they are turbulent. High-resolution simulations of the intracluster media (ICMs) and of idealized, similar media help us understand the complex physics and astrophysics involved. We present a brief overview of the physics behind ICM turbulence and outline the processes that control its development. High-resolution, compressible, isothermal MHD simulations are used to illustrate important dynamical properties of turbulence that develops in media with initially very weak magnetic fields. The simulations follow the growth of magnetic fields and reproduce the characteristics of turbulence. These results are also compared with full cluster simulations that have examined the properties of ICM turbulence.

We present and analyse near-infrared spectroscopy for a sample of 28 gravitationally- lensed star-forming galaxies in the redshift range 1.5 < z < 5, observed mostly with the Keck II telescope. With typical magnifications of ~1.5-4 magnitudes, our survey provides a valuable census of star formation rates, gas-phase metallicities and dynamical masses for a representative sample of low luminosity galaxies seen at a formative period in cosmic history. We find less evolution in the mass-metallicity relation compared to earlier work that focused on more luminous systems with z - 2-3, especially in the low mass (- 10^9 Msol) where our sample is - 0.25 dex more metal-rich. We interpret this offset as a result of the lower star formation rates (typically a factor of -10 lower) for a given stellar mass in our sub-luminous systems. Taking this effect into account, we conclude our objects are consistent with a fundamental metallicity relation recently proposed from unlensed observations.

We present near-infrared imaging spectroscopy of the strongly-lensed z=2.00 galaxy SDSS J120601.69+514227.8 (`the Clone arc'). Using OSIRIS on the Keck 2 telescope with laser guide star adaptive optics, we achieve resolved spectroscopy with 0.20 arcsecond FWHM resolution in the diagnostic emission lines [O III], Halpha, and [N II]. The lensing magnification allows us to map the velocity and star formation from Halpha emission at a physical resolution of ~300 pc in the galaxy source plane. With an integrated star formation rate of ~50 Msun/yr, the galaxy is typical of sources similarly studied at this epoch. It is dispersion-dominated with a velocity gradient of +/- 80 km/s and average dispersion sigma = 85 km/s; the dynamical mass is 2.4 \times 10^{10} Msun within a half-light radius of 2.9 kpc. Robust detection of [N II] emission across the entire OSIRIS field of view enables us to trace the gas-phase metallicity distribution with 500 pc resolution. We find a strong radial gradient in both the [N II]/Halpha and [O III]/Halpha ratios indicating a metallicity gradient of -0.27 +/- 0.05 dex/kpc with central metallicity close to solar. We demonstrate that the gradient is seen independently in two multiple images. While the physical gradient is considerably steeper than that observed in local galaxies, in terms of the effective radius at that epoch, the gradient is similar. This suggests that subsequent growth occurs in an inside-out manner with the inner metallicity gradient diminished over time due to radial mixing and enrichment from star formation.

We provide a both qualitative and quantitative comparison among different approaches aimed to solve the problem of non-linear diffusive acceleration of particles at shocks. In particular, we show that state-of-the-art models (numerical, Monte Carlo and semi-analytical), even if based on different physical assumptions and implementations, for typical environmental parameters lead to very consistent results in terms of shock hydrodynamics, cosmic ray spectrum and also escaping flux spectrum and anisotropy. Strong points and limits of each approach are also discussed, as a function of the problem one wants to study.

EBEX is a NASA-funded balloon-borne experiment designed to measure the polarization of the cosmic microwave background (CMB). Observations will be made using 1432 transition edge sensor (TES) bolometric detectors read out with frequency multiplexed SQuIDs. EBEX will observe in three frequency bands centered at 150, 250, and 410 GHz, with 768, 384, and 280 detectors in each band, respectively. This broad frequency coverage is designed to provide valuable information about polarized foreground signals from dust. The polarized sky signals will be modulated with an achromatic half wave plate (AHWP) rotating on a superconducting magnetic bearing (SMB) and analyzed with a fixed wire grid polarizer. EBEX will observe a patch covering ~1% of the sky with 8' resolution, allowing for observation of the angular power spectrum from \ell = 20 to 1000. This will allow EBEX to search for both the primordial B-mode signal predicted by inflation and the anticipated lensing B-mode signal. Calculations to predict EBEX constraints on r using expected noise levels show that, for a likelihood centered around zero and with negligible foregrounds, 99% of the area falls below r = 0.035. This value increases by a factor of 1.6 after a process of foreground subtraction. This estimate does not include systematic uncertainties. An engineering flight was launched in June, 2009, from Ft. Sumner, NM, and the long duration science flight in Antarctica is planned for 2011. These proceedings describe the EBEX instrument and the North American engineering flight.

Images of the circumstellar ejecta associated with the post-red supergiant IRC +10420 show a complex ejecta with visual evidence for episodic mass loss. In this paper we describe the transverse motions of numerous knots, arcs and condensations in the inner ejecta measured from second epoch {\it HST/WFPC2} images. When combined with the radial motions for several of the features, the total space motion and direction of the outflows show that they were ejected at different times, in different directions, and presumably from separate regions on the surface of the star. These discrete structures in the ejecta are kinematically distinct from the general expansion of the nebula and their motions are dominated by their transverse velocities. They are apparently all moving within a few degrees of the plane of the sky. We are thus viewing IRC +10420 nearly pole-on and looking nearly directly down onto its equatorial plane. We also discuss the role of surface activity and magnetic fields on IRC +10420's recent mass loss history.

We describe the cryostat and supporting electronics for the EBEX experiment. EBEX is a balloon-borne polarimeter designed to measure the B-mode polarization of the cosmic microwave background radiation. The instrument includes a 1.5 meter Gregorian-type telescope and 1432 bolometric transition edge sensor detectors operating at 0.3 K. Electronics for monitoring temperatures and controlling cryostat refrigerators is read out over CANbus. A timing system ensures the data from all subsystems is accurately synchronized. EBEX completed an engineering test flight in June 2009 during which the cryogenics and supporting electronics performed according to predictions. The temperatures of the cryostat were stable, and an analysis of a subset of the data finds no scan synchronous signal in the cryostat temperatures. Preparations are underway for an Antarctic flight.

Several neutral species (MgI, SiI, CaI, FeI) have been detected in a weak MgII absorption line system (W_r(2796)~0.15 Angstroms) at z~0.45 along the sightline toward HE0001-2340. These observations require extreme physical conditions, as noted in D'Odorico (2007). We place further constraints on the properties of this system by running a wide grid of photoionization models, determining that the absorbing cloud that produces the neutral absorption is extremely dense (~100-1000/cm^3), cold (<100 K), and has significant molecular content (~72-94%). Structures of this size and temperature have been detected in Milky Way CO surveys, and have been predicted in hydrodynamic simulations of turbulent gas. In order to explain the observed line profiles in all neutral and singly ionized chemical transitions, the lines must suffer from unresolved saturation and/or the absorber must partially cover the broad emission line region of the background quasar. In addition to this highly unusual cloud, three other ordinary weak MgII clouds (within densities of ~0.005/cm^3 and temperatures of ~10000K) lie within 500 km/s along the same sightline. We suggest that the "bare molecular cloud", which appears to reside outside of a galaxy disk, may have had in situ star formation and may evolve into an ordinary weak MgII absorbing cloud.

We report on a series of three-dimensional magnetohydrodynamic simulations of active galactic nucleus (AGN) jet propagation in realistic models of magnetized galaxy clusters. We are primarily interested in the details of energy transfer between jets and the intracluster medium (ICM) to help clarify what role such flows could have in the reheating of cluster cores. Our simulated jets feature a range of intermittency behaviors, including intermittent jets that periodically switch on and off and one model jet that shuts down completely, naturally creating a relic plume. The ICM into which these jets propagate incorporates tangled magnetic field geometries and density substructure designed to mimic some likely features of real galaxy clusters. We find that our jets are characteristically at least 60% efficient at transferring thermal energy to the ICM. Irreversible heat energy is not uniformly distributed, however, instead residing preferentially in regions very near the jet/cocoon boundaries. While intermittency affects the details of how, when, and where this energy is deposited, all of our models generically fail to heat the cluster cores uniformly. Both the detailed density structure and nominally weak magnetic fields in the ICM play interesting roles in perturbing the flows, particularly when the jets are non-steady. Still, this perturbation is never sufficient to isotropize the jet energy deposition, suggesting that some other ingredient is required for AGN jets to successfully reheat cluster cores.

We propose a model for Diffusive Shock Acceleration (DSA) in which stochastic magnetic fields in the shock precursor are generated through purely fluid mechanisms of a so-called small-scale dynamo. This contrasts with previous DSA models that considered magnetic fields amplified through cosmic ray streaming instabilities; i.e., either by way of individual particles resonant scattering in the magnetic fields, or by macroscopic electric currents associated with large-scale cosmic ray streaming. Instead, in our picture, the solenoidal velocity perturbations that are required for the dynamo to work are produced through the interactions of the pressure gradient of the cosmic ray precursor and density perturbations in the inflowing fluid. Our estimates show that this mechanism provides fast growth of magnetic field and is very generic. We argue that for supernovae shocks the mechanism is capable of generating upstream magnetic fields that are sufficiently strong for accelerating cosmic rays up to around 10^16 eV. No action of any other mechanism is necessary.

We present spatially-resolved dynamics for six strongly lensed star-forming galaxies at z=1.7-3.1, each enlarged by a linear magnification factor ~8. Using the Keck laser guide star AO system and the OSIRIS integral field unit spectrograph we resolve kinematic and morphological detail in our sample with an unprecedented fidelity, in some cases achieving spatial resolutions of ~100 pc. With one exception our sources have diameters ranging from 1-7 kpc, star formation rates of 2-40 Msun/yr (uncorrected for extinction) and dynamical masses of 10^(9.7-10.3) Msun. With this exquisite resolution we find that four of the six galaxies display coherent velocity fields consistent with a simple rotating disk model, which can only be recovered with the considerably improved spatial resolution and sampling from the combination of adaptive optics and strong gravitational lensing. Our model fits imply ratios for the systemic to random motion, V sin(i)/sigma, ranging from 0.5-1.3 and Toomre disk parameters Q<1. The large fraction of well-ordered velocity fields in our sample is consistent with data analyzed for larger, more luminous sources at this redshift. Our high resolution data further reveal that all six galaxies contain multiple giant star-forming HII regions whose resolved diameters are in the range 300 pc - 1.0 kpc, consistent with the Jeans length expected in the case of dispersion support. The density of star formation in these regions is ~100 times higher than observed in local spirals; such high values are only seen in the most luminous local starbursts. The global dynamics and demographics of star formation in these HII regions suggest that vigorous star formation is primarily governed by gravitational instability in primitive rotating disks.

Observations of X-ray cavities formed by powerful jets from AGN in galaxy cluster cores are commonly used to estimate the mechanical luminosity of these sources. We test the reliability of observationally measuring this power with synthetic X-ray observations of 3-D MHD simulations of jets in a galaxy cluster environment. We address the role that factors such as jet intermittency and orientation of the jets on the sky have on the reliability of observational measurements of cavity enthalpy and age. An estimate of the errors in these quantities can be made by directly comparing ``observationally'' derived values with values from the simulations. In our tests, cavity enthalpy, age and mechanical luminosity derived from observations are within a factor of two of the simulation values.

We present a detailed study of the spatially resolved kinematics, star-formation and stellar mass in a highly amplified galaxy at z=4.92 behind the lensing cluster MS1358+62. We use the observed optical, near- and mid-infrared imaging from HST ACS & NICMOS and Spitzer IRAC to derive the stellar mass and the Gemini/NIFS IFU to investigate the velocity structure of the galaxy from the nebular [OII] emission. Using a detailed gravitational lens model, we account for lensing amplification factor 12.+/-2.0 and find that this intrinsically L* galaxy has a stellar mass of M*=7+/-2x10^8Mo, a dynamical mass of Mdyn=3+/-1x10^9csc^2(i)Mo (within of 2kpc) and a star-formation rate of 42+/-8Mo/yr. The source-plane UV/optical morphology of this galaxy is dominated by five discrete star-forming regions. Exploiting the dynamical information we derive masses for individual star-forming regions of Mcl~10^(8-9)Mo with sizes of ~200pc. We find that, at a fixed size, the star-formation rate density within these HII regions is approximately two orders of magnitude greater than those observed in local spiral/starburst galaxies, but consistent with the most massive HII regions in the local Universe such as 30Dor. Finally, we compare the spatially resolved nebular emission line velocity with the Ly-alpha and UV-ISM lines and find that this galaxy is surrounded by a galactic scale outflow in which the Ly-alpha appears redshifted by ~150km/s and the UV-ISM lines blue-shifted by -200km/s from the (systemic) nebular emission. The velocity structure of the outflow mirrors that of the nebular emission suggesting the outflow is young (~15Myr), and has yet to burst out of the system. Taken together, these results suggest that this young galaxy is undergoing its first major epoch of mass assembly.

We outline the challenges associated with the development and construction of large focal plane arrays for use both on the ground and in space. Using lessons learned from existing JPL-led and ASU/JPL partnership efforts to develop technology for, and design such arrays and imagers for large focal planes, we enumerate here the remaining problems that need to be solved to make such a venture viable. Technologies we consider vital for further development include: (1) architectures, processes, circuits, and readout solutions for production and integration of four-side buttable, low-cost, high-fidelity, high-performance, and high-reliability CCD and CMOS imagers; (2) modular, four-side buttable packaging of CCD/CMOS imagers; (3) techniques and hardware to test and characterize the large number of chips required to produce the hundreds of flight-grade detectors needed for large focal-plane missions being conceived at this time; (4) ground based testbed needs, such as a large format camera mounted on a ground-based telescope, to field test the detectors and the focal plane technology solutions; and (5) validation of critical sub-components of the design on a balloon mission to ensure their flight-readiness. This paper outlines the steps required to provide a mature solution to the astronomical community with a minimal investment, building on years of planning and investments already completed at JPL.

A description is given of the algorithms implemented in the AstroBEAR adaptive mesh refinement code for ideal magnetohydrodynamics. The code provides several high resolution, shock capturing schemes which are constructed to maintain conserved quantities of the flow in a finite volume sense. Divergence free magnetic field topologies are maintained to machine precision by collating the components of the magnetic field on a cell-interface staggered grid and utilizing the constrained transport approach for integrating the induction equations. The maintenance of magnetic field topologies on adaptive grids is achieved using prolongation and restriction operators which preserve the divergence and curl of the magnetic field across co-located grids of different resolution. The robustness and correctness of the code is demonstrated by comparing the numerical solution of various tests with analytical solutions or previously published numerical solutions obtained by other codes.

The startling discovery of Prochter et al. (2006) that the frequency of very strong (W_r(2796)>1 A) MgII absorbers along gamma-ray burst (GRB) lines of sight ([dN/dz]_{GRB} = 0.90) is more than three times the frequency along quasar lines of sight ([dN/dz]_{QSO} = 0.24), over similar redshift ranges, has yet to be understood. We reconsider the possibility that the excess of very strong MgII absorbers toward GRBs is intrinsic either to the GRBs themselves or to their immediate environment, and associated with bulk outflows with velocities as large as v_{max} ~ 0.3c. In order to examine this hypothesis, we accumulate a sample of 27 W_r(2796) > 1 A absorption systems found toward 81 quasars, and compare their properties to those of 9 W_r(2796)>1 A absorption systems found toward 6 GRBs; all systems have been observed at high spectral resolution (R = 45,000) using the Ultraviolet and Visual Echelle Spectrograph on the Very Large Telescope. We make multiple comparisons of the absorber properties across the two populations, testing for differences in metallicity, ionization state, abundance patterns, dust abundance, kinematics, and phase structure. We find no significant differences between the two absorber populations using any of these metrics, implying that, if the excess absorbers toward GRB lines of sight are indeed intrinsic, they must be produced in a process which has strong similarities to the processes yielding strong MgII systems in association with intervening galaxies. Although this may seem a priori unlikely, given the high outflow velocities required for any intrinsic model, we note that the same conclusion was reached, recently, with respect to the narrow absorption line systems seen in some quasars.

For over four decades, synchrotron-radiating sources have played a series of pathfinding roles in the study of galaxy clusters and large scale structure. Such sources are uniquely sensitive to the turbulence and shock structures of large-scale environments, and their cosmic rays and magnetic fields often play important dynamic and thermodynamic roles. They provide essential complements to studies at other wavebands. Over the next decade, they will fill essential gaps in both cluster astrophysics and the cosmological growth of structure in the universe, especially where the signatures of shocks and turbulence, or even the underlying thermal plasma itself, are otherwise undetectable. Simultaneously, synchrotron studies offer a unique tool for exploring the fundamental question of the origins of cosmic magnetic fields. This work will be based on the new generation of m/cm-wave radio telescopes now in construction, as well as major advances in the sophistication of 3-D MHD simulations.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the fields of stars and stellar evolution. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of Galactic science. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the fields of planetary systems and star formation. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of extragalactic astrophysics. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of cosmology. Community-based White Paper to Astro2010 in response to the call for such papers.

We present results from and the analysis of data from MAXIPOL, a balloon-borne experiment designed to measure the polarization in the Cosmic Microwave Background (CMB). MAXIPOL is the first CMB experiment to obtain results using a rotating half-wave plate as a rapid polarization modulator. We report results from observations of a sky area of 8 deg^2 with 10-arcmin resolution, providing information up to l~700. We use a maximum-likelihood method to estimate maps of the Q and U Stokes parameters from the demodulated time streams, and then both Bayesian and frequentist approaches to compute the EE, EB, and BB power spectra. Detailed formalisms of the analyses are given. A variety of tests show no evidence for systematic errors. The Bayesian analysis gives weak evidence for an EE signal. The EE power is 55^{+51}_{-45} \mu K^2 at the 68% confidence level for l=151-693. Its likelihood function is asymmetric and skewed positive such that with a uniform prior the probability of a positive EE power is 96%. The powers of EB and BB signals at the 68% confidence level are 18^{+27}_{-34} \mu K^2 and -31^{+31}_{-19} \mu K^2 respectively and thus consistent with zero. The upper limit of the BB-mode at the 95% confidence level is 9.5 \mu K. Results from the frequentist approach are in agreement within statistical errors. These results are consistent with the current concordance LCDM model.

We use kinetic simulations of diffusive shock acceleration (DSA) to study the time-dependent evolution of plane, quasi-parallel, cosmic-ray (CR) modified shocks. Thermal leakage injection of low energy CRs and finite Alfv\'en wave propagation and dissipation are included. Bohm diffusion as well as the diffusion with the power-law momentum dependence are modeled. As long as the acceleration time scale to relativistic energies is much shorter than the dynamical evolution time scale of the shocks, the precursor and subshock transition approach the time-asymptotic state, which depends on the shock sonic and Alfv\'enic Mach numbers and the CR injection efficiency. For the diffusion models we employ, the shock precursor structure evolves in an approximately self-similar fashion, depending only on the similarity variable, x/(u_s t). During this self-similar stage, the CR distribution at the subshock maintains a characteristic form as it evolves: the sum of two power-laws with the slopes determined by the subshock and total compression ratios with an exponential cutoff at the highest accelerated momentum, p_{max}(t). Based on the results of the DSA simulations spanning a range of Mach numbers, we suggest functional forms for the shock structure parameters, from which the aforementioned form of CR spectrum can be constructed. These analytic forms may represent approximate solutions to the DSA problem for astrophysical shocks during the self-similar evolutionary stage as well as during the steady-state stage if p_{max} is fixed.

We report results of 3D MHD simulations of the dynamics of buoyant bubbles in magnetized galaxy cluster media. The simulations are three dimensional extensions of two dimensional calculations reported by Jones & De Young (2005). Initially spherical bubbles and briefly inflated spherical bubbles all with radii a few times smaller than the intracluster medium (ICM) scale height were followed as they rose through several ICM scale heights. Such bubbles quickly evolve into a toroidal form that, in the absence of magnetic influences, is stable against fragmentation in our simulations. This ring formation results from (commonly used) initial conditions that cause ICM material below the bubbles to drive upwards through the bubble, creating a vortex ring; that is, hydrostatic bubbles develop into "smoke rings", if they are initially not very much smaller or very much larger than the ICM scale height. Even modest ICM magnetic fields with beta = P_gas/P_mag ~ 10^3 can influence the dynamics of the bubbles, provided the fields are not tangled on scales comparable to or smaller than the size of the bubbles. Quasi-uniform, horizontal fields with initial beta ~ 10^2 bifurcated our bubbles before they rose more than about a scale height of the ICM, and substantially weaker fields produced clear distortions. On the other hand, tangled magnetic fields with similar, modest strengths are generally less easily amplified by the bubble motions and are thus less influential in bubble evolution. Inclusion of a comparably strong, tangled magnetic field inside the initial bubbles had little effect on our bubble evolution, since those fields were quickly diminished through expansion of the bubble and reconnection of the initial field.

This is a report on the findings of the extragalactic science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses extragalactic science topics including active galactic nuclei, cosmic ray acceleration in galaxies, galaxy clusters and large scale structure formation shocks, and the study of the extragalactic infrared and optical background radiation. The scientific potential of ground based gamma-ray observations of Gamma-Ray Bursts and dark matter annihilation radiation is covered in other sections of the white paper.

This is a report on the findings of the SNR/cosmic-ray working group for the white paper on the status and future of ground-based gamma-ray astronomy. The white paper is an APS commissioned document, and the overall version has also been released and can be found on astro-ph. This detailed section of the white paper discusses the status of past and current attempts to observe shell-type supernova remnants and diffuse emission from cosmic rays at GeV-TeV energies. We concentrate on the potential of future ground-based gamma-ray experiments to study the acceleration of relativistic charged particles which is one of the main unsolved, yet fundamental, problems in modern astrophysics. The acceleration of particles relies on interactions between energetic particles and magnetic turbulence. In the case of SNRs we can perform spatially resolved studies in systems with known geometry, and the plasma physics deduced from these observations will help us to understand other systems where rapid particle acceleration is believed to occur and where observations as detailed as those of SNRs are not possible.

We study the performance of achromatic half-wave plates (AHWP) as a function of their construction parameters, the detection bandwidth of a power detector operating in the millimeter wave band, and the spectral shape of the incident radiation. We focus particular attention on the extraction of the degree of incident polarization and its orientation angle from the intensity measured as a function of AHWP rotation angle, which we call the IVA (intensity versus angle). We quantify the phase offset of the IVA and point to potential systematic errors in the extraction of this offset in cases where the incident spectrum is not sufficiently well known. We show how the phase offset and modulation efficiency of the AHWP depend on the relative angles between the plates in the stack and find that high modulation efficiency can be achieved with alignment accuracy of few degrees.

We report imaging polarimetry of segments B and C of the Jupiter-family Comet 73P/Schwassmann-Wachmann 3 in the I and H bandpasses at solar phase angles of approximately 35 and 85deg. The level of polarization was typical for active comets, but larger than expected for a Jupiter-family comet. The polarimetric color was slightly red (dP/dL = +1.2 +/- 0.4) at a phase angle of ~ 35deg and either neutral or slightly blue at a phase angle of ~ 85deg. Observations during the closest approach from 2006 May 11-13 achieved a resolution of 35 km at the nucleus. Both segments clearly depart from a 1/rho surface brightness for the first 50 - 200 km from the nucleus. Simulations of radiation driven dust dynamics can reproduce some of the observed coma morphology, but only with a wide distribution of initial dust velocities (at least a factor of 10) for a given grain radius. Grain aggregate breakup and fragmentation are able to reproduce the observed profile perpendicular to the Sun-Comet axis, but fit the observations less well along this axis (into the tail). The required fragmentation is significant, with a reduction in the mean grain aggregate size by about a factor of 10. A combination of the two processes could possibly explain the surface brightness profile of the comet.

Multi-wavelength imaging polarimetry at far-infrared wavelengths has proven to be an excellent tool for studying the physical properties of dust, molecular clouds, and magnetic fields in the interstellar medium. Although these wavelengths are only observable from airborne or space-based platforms, no first-generation instrument for the Stratospheric Observatory for Infrared Astronomy (SOFIA) is presently designed with polarimetric capabilities. We study several options for upgrading the High-resolution Airborne Wideband Camera (HAWC) to a sensitive FIR polarimeter. HAWC is a 12 x 32 pixel bolometer camera designed to cover the 53 - 215 micron spectral range in 4 colors, all at diffraction-limited resolution (5 - 21 arcsec). Upgrade options include: (1) an external set of optics which modulates the polarization state of the incoming radiation before entering the cryostat window; (2) internal polarizing optics; and (3) a replacement of the current detector array with two state-of-the-art superconducting bolometer arrays, an upgrade of the HAWC camera as well as polarimeter. We discuss a range of science studies which will be possible with these upgrades including magnetic fields in star-forming regions and galaxies and the wavelength-dependence of polarization.

Chang et al. (2006) reported millisecond duration dips in the X-ray intensity of Sco X-1 and attributed them to occultations of the source by small trans-Neptunian objects (TNOs). We have found multiple lines of evidence that these dips are not astronomical in origin, but rather the result of high-energy charged particle events in the RXTE PCA detectors. Our analysis of the RXTE data indicates that at most 10% of the observed dips in Sco X-1 could be due to occultations by TNOs, and, furthermore, we find no positive or supporting evidence for any of them being due to TNOs. We therefore believe that it is a mistake to conclude that any TNOs have been detected via occultation of Sco X-1.

Magnetic fields probably play a central role in the dynamics and thermodynamics of ICMs and their interactions with AGNs, despite the fact that the fields usually contribute relatively little pressure; i.e., the ICM is a ``high-$\beta$'' plasma. More typically, the roles of magnetic fields come through ``microscopic'' influences on charged particle behaviors, and through magnetic tension, which can still be significant in subsonic, high-$\beta$ flows. I briefly review these issues, while exploring the underlying question of using the commonly-applied magnetohydrodynamics model in the ICM when Coulomb scattering mean free paths can sometimes exceed tens of kiloparsecs.

We present results of time evolution of oblique MHD plane shocks including diffusive cosmic ray acceleration with backreaction on the plasma flows. The simulations include self-consistent effects of finite Alfven wave propagation and dissipation. From the computed cosmic ray particle phase space distributions we calculate expected leptonic and hadronic emissions resulting from interactions between the cosmic rays, magnetic fields, the thermal particle population and relevant astrophysical photon fields.

Using an improved version of the previously introduced CRASH (Cosmic Ray Acceleration SHock) code, we have calculated the time evolution of cosmic-ray (CR) modified quasi-parallel plane shocks for Bohm-like diffusion, including self-consistent models of Alfven wave drift and dissipation, along with thermal leakage injection of CRs. The new simulations follow evolution of the CR distribution to much higher energies than our previous study, providing a better examination of evolutionary and asymptotic behaviors. The postshock CR pressure becomes constant after quick initial adjustment, since the evolution of the CR partial pressure expressed in terms of a momentum similarity variable is self-similar. The shock precursor, which scales as the diffusion length of the highest energy CRs, subsequently broadens approximately linearly with time, independent of diffusion model, so long as CRs continue to be accelerated to ever-higher energies. This means the nonlinear shock structure can be described approximately in terms of the similarity variable, x/(u_s t), where u_s is the shock speed once the postshock pressure reaches an approximate time asymptotic state. As before, the shock Mach number is the key parameter determining the evolution and the CR acceleration efficiency, although finite Alfven wave drift and wave energy dissipation in the shock precursor reduce the effective velocity change experienced by CRs, so reduce acceleration efficiency noticeably, thus, providing a second important parameter at low and moderate Mach numbers.

We use imaging polarimetry taken with the HST/ACS/HRC to explore the three dimensional structure of the circumstellar dust distribution around the red supergiant VY Canis Majoris. The polarization vectors of the nebulosity surrounding VY CMa show a strong centro-symmetric pattern in all directions except directly East and range from 10% - 80% in fractional polarization. In regions that are optically thin, and therefore likely have only single scattering, we use the fractional polarization and photometric color to locate the physical position of the dust along the line-of-sight. Most of the individual arc-like features and clumps seen in the intensity image are also features in the fractioanl polarization map. These features must be distinct geometric objects. If they were just local density enhancements, the fractional polarization would not change so abruptly at the edge of the feature. The location of these features in the ejecta of VY CMa using polarimetry provides a determination of their 3D geometry independent of, but in close agreement with, the results from our study of their kinematics (Paper I).

Images of the complex circumstellar nebula associated with the famous red supergiant VY CMa show evidence for multiple and asymmetric mass loss events over the past 1000 yrs. Doppler velocities of the arcs and knots in the ejecta showed that they are not only spatially distinct but also kinematically separate from the surrounding diffuse material. In this paper we describe second epoch HST/WFPC2 images to measure the transverse motions which when combined with the radial motions provide a complete picture of the kinematics of the ejecta including the total space motions and directions of the outflows. Our results show that the arcs and clumps of knots are moving at different velocities, in different directions, and at different angles relative to the plane of the sky and to the star, confirming their origin from eruptions at different times and from physically separate regions on the star. We conclude that the morphology and kinematics of the arcs and knots are consistent with a history of mass ejections not aligned with any presumed axis of symmetry. The arcs and clumps represent relatively massive outflows and ejections of gas very likely associated with large -- scale convective activity and magnetic fields.

Chang et al. (2006) reported millisecond duration dips in the X-ray intensity of Sco X-1 and attributed them to occultations of the source by small trans-Neptunian objects (TNOs). We have found evidence that these dips are in fact not astronomical in origin, but rather the result of high-energy charged particle events in the RXTE PCA detectors.

We discuss MAXIPOL, a bolometric balloon-borne experiment designed to measure the E-mode polarization of the cosmic microwave background radiation (CMB). MAXIPOL is the first bolometric CMB experiment to observe the sky using rapid polarization modulation. To build MAXIPOL, the CMB temperature anisotropy experiment MAXIMA was retrofitted with a rotating half-wave plate and a stationary analyzer. We describe the instrument, the observations, the calibration and the reduction of data collected with twelve polarimeters operating at 140 GHz and with a FWHM beam size of 10 arcmin. We present maps of the Q and U Stokes parameters of an 8 deg^2 region of the sky near the star Beta Ursae Minoris. The power spectra computed from these maps give weak evidence for an EE signal. The maximum-likelihood amplitude of l(l+1)C^{EE}_{l}/(2 pi) is 55_{-45}^{+51} uK^2 (68%), and the likelihood function is asymmetric and skewed positive such that with a uniform prior the probability that the amplitude is positive is 96%. This result is consistent with the expected concordance LCDM amplitude of 14 uK^2. The maximum likelihood amplitudes for l(l+1)C^{BB}_{l}/(2 pi) and $\ell(\ell+1)C^{EB}_{\ell}/2\pi$ are -31_{-19}^{+31} and 18_{-34}^{+27} uK^2 (68%), respectively, which are consistent with zero. All of the results are for one bin in the range 151 < l < 693. Tests revealed no residual systematic errors in the time or map domain. A comprehensive discussion of the analysis of the data is presented in a companion paper.

We have developed a cosmic ray (CR) shock code in one dimensional spherical geometry with which the particle distribution, the gas flow and their nonlinear interaction can be followed numerically in a frame comoving with an expanding shock. In order to accommodate a very wide dynamic range of diffusion length scales in the CR shock problem, we have incorporated subzone shock tracking and adaptive mesh refinement techniques. We find the spatial grid resolution required for numerical convergence is less stringent in this code compared to typical, fixed-grid Eulerian codes. The improved convergence behavior derives from maintaining the shock discontinuity inside the same grid zone in the comoving code. That feature improves numerical estimates of the compression rate experienced by CRs crossing the subshock compared to codes that allow the subshock to drift on the grid. Using this code with a Bohm-like diffusion model we have calculated the CR acceleration and the nonlinear feedback at supernova remnant shocks during the Sedov-Taylor stage. Similarly to plane-parallel shocks, with an adopted thermal leakage injection model, about 10^{-3} of the particles that pass through the shock and up to 60 % of the explosion energy are transferred to the CR component. These results are in good agreement with previous nonlinear spherical CR shock calculations of Berezhko and collaborators.

Variable A in M33 is a member of a rare class of highly luminous, evolved stars near the upper luminosity boundary that show sudden and dramatic shifts in apparent temperature due to the formation of optically thick winds in high mass loss episodes. Recent optical and infrared spectroscopy and imaging reveal that its ``eruption'' begun in $\sim$1950 has ended, {\it lasting $\approx$ 45 yrs}. Our current observations show major changes in its wind from a cool, dense envelope to a much warmer state surrounded by low density gas with rare emission lines of Ca II, [Ca II] and K I. Its spectral energy distribution has unexpectedly changed, especially at the long wavelengths, with a significant decrease in its apparent flux, while the star remains optically obscured. We conclude that much of its radiation is now escaping out of our line of sight. We attribute this to the changing structure and distribution of its circumstellar ejecta corresponding to the altered state of its wind as the star recovers from a high mass loss event.

We present a set of high-resolution 3D MHD simulations of steady light, supersonic jets, exploring the influence of jet Mach number and the ambient medium on jet propagation and energy deposition over long distances. The results are compared to simple self-similar scaling relations for the morphological evolution of jet-driven structures and to previously published 2D simulations. For this study we simulated the propagation of light jets with internal Mach numbers 3 and 12 to lengths exceeding 100 initial jet radii in both uniform and stratified atmospheres. The propagating jets asymptotically deposit approximately half of their energy flux as thermal energy in the ambient atmosphere, almost independent of jet Mach number or the external density gradient. Nearly one-quarter of the jet total energy flux goes directly into dissipative heating of the ICM, supporting arguments for effective feedback from AGNs to cluster media. The remaining energy resides primarily in the jet and cocoon structures. Despite having different shock distributions and magnetic field features, global trends in energy flow are similar among the different models. As expected the jets advance more rapidly through stratified atmospheres than uniform environments. The asymptotic head velocity in King-type atmospheres shows little or no deceleration. This contrasts with jets in uniform media with heads that are slowed as they propagate. This suggests that the energy deposited by jets of a given length and power depends strongly on the structure of the ambient medium. While our low-Mach jets are more easily disrupted, their cocoons obey evolutionary scaling relations similar to the high-Mach jets.

We have developed a new, very efficient numerical scheme to solve the CR diffusion convection equation that can be applied to the study of the nonlinear time evolution of CR modified shocks for arbitrary spatial diffusion properties. The efficiency of the scheme derives from its use of coarse-grained finite momentum volumes. This approach has enabled us, using $\sim 10 - 20$ momentum bins spanning nine orders of magnitude in momentum, to carry out simulations that agree well with results from simulations of modified shocks carried out with our conventional finite difference scheme requiring more than an order of magnitude more momentum points. The coarse-grained, CGMV scheme reduces execution times by a factor approximately half the ratio of momentum bins used in the two methods. Depending on the momentum dependence of the diffusion, additional economies in required spatial and time resolution can be utilized in the CGMV scheme, as well. These allow a computational speed-up of at least an order of magnitude in some cases.

In order to better understand the origin and evolution of relic radio bubbles in clusters of galaxies, we report on an extensive set of 2D MHD simulations of hot buoyant bubbles evolving in a realistic intracluster medium. Our bubbles are inflated near the base of the ICM over a finite time interval from a region whose magnetic field is isolated from the ICM. We confirm both the early conjecture from linear analysis and the later results based on preformed MHD bubbles; namely, that very modest ICM magnetic fields can stabilize the rising bubbles against disruption by Rayleigh-Taylor and Kelvin-Helmholtz instabilities. We find in addition that amplification of the ambient fields as they stretch around the bubbles can be sufficient to protect the bubbles or their initial fragments even if the fields are initially much too weak to play a significant role early in the evolution of the bubbles. Indeed, even with initial fields less than a micro-Gauss and values of $\beta = P_g/P_b$ approaching $10^5$, magnetic stresses in our simulations eventually became large enough to influence the bubble evolution. Magnetic field influence also depends significantly on the geometry of the ICM field and on the topology of the field at the bubble/ICM interface. For example, reconnection of anti-parallel fields across the bubble top greatly reduced the ability of the magnetic field to inhibit disruptive instabilities. Our results confirm earlier estimates of $10^8$ yr for relic radio bubble lifetimes and show that magnetic fields can account for the long term stability of these objects against disruption by surface instabilities. In addition these calculations show that lifting and mixing of the ambient ICM may be a critical function of field geometries in both the ICM and in the bubble interior.

EBEX is a balloon-borne polarimeter designed to measure the intensity and polarization of the cosmic microwave background radiation. The measurements would probe the inflationary epoch that took place shortly after the big bang and would significantly improve constraints on the values of several cosmological parameters. EBEX is unique in its broad frequency coverage and in its ability to provide critical information about the level of polarized Galactic foregrounds which will be necessary for all future CMB polarization experiments. EBEX consists of a 1.5 m Dragone-type telescope that provides a resolution of less than 8 arcminutes over four focal planes each of 4 degree diffraction limited field of view at frequencies up to 450 GHz. The experiment is designed to accommodate 330 transition edge bolometric detectors per focal plane, for a total of up to 1320 detectors. EBEX will operate with frequency bands centered at 150, 250, 350, and 450 GHz. Polarimetry is achieved with a rotating achromatic half-wave plate. EBEX is currently in the design and construction phase, and first light is scheduled for 2008.

We estimate the Mass Function of void galaxies in the second public data release of the Sloan Digital Sky Survey from a sample of 1000 galaxies with local density contrasts of delta_v < -0.6. The galaxy sample is split into ellipticals and spirals using a color-Sersic index criteria. We estimate the virial masses of ellipticals using the measured spectral line-widths along with the observed size. Projection effects and uncertainties in halo properties make mass estimates of spirals more difficult. We use an inversion of the Tully-Fisher relation to estimate the isothermal rotational velocity, and introduce a scaling factor to estimate the halo extent. We then fit the measured mass function against a theoretical Press-Schechter model, and find that the distribution of galaxies in voids appears to be nearly unbiased compared to the mass.

We have calculated the evolution of cosmic ray (CR) modified astrophysical shocks for a wide range of shock Mach numbers and shock speeds through numerical simulations of diffusive shock acceleration (DSA) in 1D quasi- parallel plane shocks. The simulations include thermal leakage injection of seed CRs, as well as pre-existing, upstream CR populations. Bohm-like diffusion is assumed. We model shocks similar to those expected around cosmic structure pancakes as well as other accretion shocks driven by flows with upstream gas temperatures in the range $T_0=10^4-10^{7.6}$K and shock Mach numbers spanning $M_s=2.4-133$. We show that CR modified shocks evolve to time-asymptotic states by the time injected particles are accelerated to moderately relativistic energies ($p/mc \gsim 1$), and that two shocks with the same Mach number, but with different shock speeds, evolve qualitatively similarly when the results are presented in terms of a characteristic diffusion length and diffusion time. For these models the time asymptotic value for the CR acceleration efficiency is controlled mainly by shock Mach number. The modeled high Mach number shocks all evolve towards efficiencies $\sim 50$%, regardless of the upstream CR pressure. On the other hand, the upstream CR pressure increases the overall CR energy in moderate strength shocks ($M_s \sim {\rm a few}$). (abridged)

We present observations of the giant HII region complex N159 in the LMC using IRAC on the {\it Spitzer Space Telescope}. One of the two objects previously identified as protostars in N159 has an SED consistent with classification as a Class I young stellar object (YSO) and the other is probably a Class I YSO as well, making these two stars the youngest stars known outside the Milky Way. We identify two other sources that may also be Class I YSOs. One component, N159AN, is completely hidden at optical wavelengths, but is very prominent in the infrared. The integrated luminosity of the entire complex is L $\approx 9\times10^6$L$_{\odot}$, consistent with the observed radio emission assuming a normal Galactic initial mass function (IMF). There is no evidence for a red supergiant population indicative of an older burst of star formation. The N159 complex is 50 pc in diameter, larger in physical size than typical HII regions in the Milky Way with comparable luminosity. We argue that all of the individual components are related in their star formation history. The morphology of the region is consistent with a wind blown bubble $\approx 1-2Myr-old that has initiated star formation now taking place at the rim. Other than its large physical size, star formation in N159 appears to be indistinguishable from star formation in the Milky Way.

e present Near Infrared imaging polarimetry of three regions of massive star formation, G$192.16 - 3.82$, Cepheus A, and W42. In W42 we have discovered a new bipolar nebula located at the far side of the HII region behind the visible cluster of exciting stars. The axis of this new nebula is aligned with the magnetic field threading the entire cluster region. Polarization in the bipolar outflow nebulosity associated with G192.16 is consistent with a single illuminating source, too faint to be detected at $2 \micron$. Polarization in the reflection nebulosity associated with Ceph A requires more than one illuminating source, although HW2 is clearly dominant. In all three objects, the magnetic field in the outflow at distances greater than $\sim 0.2$ pc is radial. In G192.16 the magnetic field geometry closer than $\sim 0.2$ pc to the embedded star appears chaotic. For G192.16 the outflow is not aligned with the surrounding magnetic field, which lies in the galactic plane. In Ceph A, the outflow axis could be interpreted as being aligned with the galactic plane, but the magnetic field threading the region is not. Only in the case of W42 is the magnetic field threading the HII region aligned with the mean field in the surrounding galactic plane.

We present high-resolution X-ray proper motion measurements of Cassiopeia A using Chandra observations from 2000 and 2002. We separate the emission into four spectrally distinct classes: Si-dominated, Fe-dominated, low-energy-enhanced, and continuum-dominated. These classes also represent distinct spatial and kinematic components. The Si- and Fe-dominated classes are ejecta and have a mean expansion rate of 0.2%/yr. This is the same as for the forward shock filaments but less than the 0.3%/yr characteristic of optical ejecta. The low-energy-enhanced spectral class possibly illuminates a clumpy circumstellar component and has a mean expansion rate of 0.05%/yr. The continuum-dominated emission likely represents the forward shock and consists of diffuse circumstellar material which is seen as a circular ring around the periphery of the remnant as well as projected across the center.

Observations with the Extreme Ultraviolet Explorer (EUVE) have shown the Coma Cluster to be a source of EUV emission in excess of that produced by X-ray gas in the cluster. We have re-examined the EUVE data on this cluster in an attempt to obtain clues as to the origin of this emission. We find two important new results. First, the ratio between the azimuthally averaged EUV excess emission and the ROSAT hard X-ray flux is constant as a function of distance from the cluster center outward. Second, a correlation analysis between the EUV excess emission and the X-ray emission shows that on a detailed level the EUV excess is spatially closely related to the X-ray emission. These findings contradict previous suggestions as to the underlying source of the diffuse EUV emission in Coma and provide important information in regards to the true source of this emission. We propose a new explanation for the source of this emission: inverse Compton scattering of microwave background photons by secondary electrons and positrons. We explore this possibility in some detail and show that it is consistent with all of the available observational evidence. The parent cosmic ray protons may have been produced by any of a number of sources, including supernovae, active galaxies, galactic winds, and cluster formation shocks, but we believe that the most likely source is cluster formation shocks. If the EUV emission in the Coma Cluster is, in fact, the result of secondary electrons, this may be the only direct evidence for secondary electrons in the intracluster medium of a cluster of galaxies, since recent work suggests that secondary electrons may not be the cause of radio halos.

We present an extensive synthetic observational analysis of numerically- simulated radio galaxies designed to explore the effectiveness of conventional observational analyses at recovering physical source properties. These are the first numerical simulations with sufficient physical detail to allow such a study. The present paper focuses on extraction of magnetic field properties from nonthermal intensity information. Synchrotron and inverse-Compton intensities provided meaningful information about distributions and strengths of magnetic fields, although considerable care was called for. Correlations between radio and X-ray surface brightness correctly revealed useful dynamical relationships between particles and fields. Magnetic field strength estimates derived from the ratio of X-ray to radio intensity were mostly within about a factor of two of the RMS field strength along a given line of sight. When emissions along a given line of sight were dominated by regions close to the minimum energy/equipartition condition, the field strengths derived from the standard power-law-spectrum minimum energy calculation were also reasonably close to actual field strengths, except when spectral aging was evident. Otherwise, biases in the minimum- energy magnetic field estimation mirrored actual differences from equipartition. The ratio of the inverse-Compton magnetic field to the minimum-energy magnetic field provided a rough measure of the actual total energy in particles and fields in most instances, within an order of magnitude. This may provide a practical limit to the accuracy with which one may be able to establish the internal energy density or pressure of optically thin synchrotron sources.

We discuss MAXIPOL, a bolometric balloon-borne experiment designed to measure the E-mode polarization anisotropy of the cosmic microwave background radiation (CMB) on angular scales of 10 arcmin to 2 degrees. MAXIPOL is the first CMB experiment to collect data with a polarimeter that utilizes a rotating half-wave plate and fixed wire-grid polarizer. We present the instrument design, elaborate on the polarimeter strategy and show the instrument performance during flight with some time domain data. Our primary data set was collected during a 26 hour turnaround flight that was launched from the National Scientific Ballooning Facility in Ft. Sumner, New Mexico in May 2003. During this flight five regions of the sky were mapped. Data analysis is in progress.

Interstellar polarization from aligned dust grains can be measured both in transmission at visible and near-infrared wavelengths and in emission at far-infrared and sub-mm wavelengths. These observations can help predict the behavior of foreground contamination of CMB polarimetry by dust in the Milky Way. Fractional polarization in emission from aligned dust grains will be at the higher range of currently observed values of 4-10%. Away from the galactic plane, fluctuations in Q and U will be dominated by fluctuations in intensity, and less influenced by fluctuations in fractional polarization and position angle.

We have studied forced turbulence of compressible magnetohydrodynamic (MHD) flows through two-dimensional simulations with different numerical resolutions. First, hydrodynamic turbulence with Mach number $<M_s >_{\rm init} \equiv < v >_{\rm rms}/ c_s = 1$ and density compression ${< \delta\rho / \rho >}_{\rm rms} \simeq 0.45$ was generated by enforcing a random force. Then, initial, uniform magnetic fields of various strengths were added with Alfv\'enic Mach number $<M_A >_{\rm init} \equiv < v >_{\rm rms} / c_{A, {\rm init}} \gg 1$. An isothermal equation of state was employed, and no explicit dissipation was included. After the MHD turbulence is saturated, the resulting flows are categorized as very weak field (VWF), weak field (WF), and strong field (SF) classes, which have $<M_A > \equiv < v >_{\rm rms} / < c_A >_{\rm rms} \gg 1$, $<M_A > > 1$, and $<M_A > \sim 1$, respectively. Not only in the SF regime but also in the WF regime, turbulent transport is suppressed by the magnetic field. In the SF cases, the energy power spectra in the inertial range, although no longer power-law, exhibit a range with slopes close to $\sim1.5$, hinting the Iroshnikov-Kraichnan spectrum. Our simulations were able to produce the SF class behaviors only with high resolution of at least $1024^2$ grid cells. The specific requirements for the simulation of the SF class should depend on the code (and the numerical scheme) as well as the initial setup, but our results do indicate that very high resolution would be required for converged results in simulation studies of MHD turbulence.

We study the properties of cosmological shock waves identified in high-resolution, N-body/hydrodynamic simulations of a $\Lambda$CDM universe and their role on thermalization of gas and acceleration of nonthermal, cosmic ray (CR) particles. External shocks form around sheets, filaments and knots of mass distribution when the gas in void regions accretes onto them. Within those nonlinear structures, internal shocks are produced by infall of previously shocked gas to filaments and knots, and during subclump mergers, as well as by chaotic flow motions. Due to the low temperature of the accreting gas, the Mach number of external shocks is high, extending up to $M\sim 100$ or higher. In contrast, internal shocks have mostly low Mach numbers. For all shocks of $M\ge1.5$ the mean distance between shock surfaces over the entire computed volume is $\sim4 h^{-1}$ Mpc at present, or $\sim 1 h^{-1}$ Mpc for internal shocks within nonlinear structures. Identified external shocks are more extensive, with their surface area $\sim2$ times larger than that of identified internal shocks at present. However, especially because of higher preshock densities, but also due to higher shock speeds, internal shocks dissipate more energy. Hence, the internal shocks are mainly responsible for gas thermalization as well as CR acceleration. In fact, internal shocks with $2 \la M \la 4$ contribute $\sim 1/2$ of the total dissipation. Using a nonlinear diffusive shock acceleration model for CR protons, we estimate the ratio of CR energy to gas thermal energy dissipated at cosmological shock waves to be $\sim1/2$ through the history of the universe. Our result supports scenarios in which the intracluster medium contains energetically significant populations of CRs.

Measurements of the polarization of the cosmic microwave background (CMB) radiation are expected to significantly increase our understanding of the early universe. We present a design for a CMB polarimeter in which a cryogenically cooled half wave plate rotates by means of a high-temperature superconducting (HTS) bearing. The design is optimized for implementation in MAXIPOL, a balloon-borne CMB polarimeter. A prototype bearing, consisting of commercially available ring-shaped permanent magnet and an array of YBCO bulk HTS material, has been constructed. We measured the coefficient of friction as a function of several parameters including temperature between 15 and 80 K, rotation frequency between 0.3 and 3.5 Hz, levitation distance between 6 and 10 mm, and ambient pressure between 10^{-7} and 1 torr. The low rotational drag of the HTS bearing allows rotations for long periods of time with minimal input power and negligible wear and tear thus making this technology suitable for a future satellite mission.

The Parker instability, which has been considered as a process governing the structure of the interstellar medium, is induced by the buoyancy of magnetic field and cosmic rays. In previous studies, while the magnetic field has been fully incorporated in the context of isothermal magnetohydrodynamics, cosmic rays have been normally treated with the simplifying assumption of infinite diffusion along magnetic field lines but no diffusion across them. The cosmic ray diffusion is, however, finite. In this work, we take into account fully the diffusion process of cosmic rays in a linear stability analysis of the Parker instability. Cosmic rays are described with the diffusion-convection equation. With realistic values of cosmic ray diffusion coefficients expected in the interstellar medium, we show that the result of previous studies with the simplifying assumption on cosmic ray diffusion applies well. Finiteness of parallel diffusion decreases the growth rate of the Parker instability, while the relatively smaller perpendicular diffusion has no significant effect. We discuss the implication of our result on the role of the Parker instability in the interstellar medium.

We present broadband imaging polarimetry of DR21 at 2.2$\mu$m. Background stars shining through the lobes of the bipolar outflow show polarization aligned with the long axis of the outflow, indicating a magnetic field geometry oriented along the flow axis. There is no indication of a spiral or turbulent magnetic field geometry in the lobes. The polarization of stars in the central cluster has a different position angle than the lobes and is in good agreement with millimeter polarimetry. The nebulosity in the Eastern lobe has moderate to high polarization consistent with scattering of continuum light from the central cluster. We were unable to detect polarization of the nebulosity in the Western lobe at the 4.2% (3sigma) level.

We report the discovery of polarized flux at 2.2 micron from the bright shell of the approximately 320 year old supernova remnant Cas A. The fractional polarizations are comparable at 6 cm and 2.2 micron, and the polarization angles are similar, demonstrating that synchrotron radiation from the same relativistic plasma is being observed at these widely separated wavebands. The relativistic electrons radiating at 2.2 micron have an energy of ~ 150 GeV, (gamma ~ 3e5), assuming an ~500 microGauss magnetic field. The total intensity at 2.2 micron lies close to the power law extrapolation from radio frequencies, showing that relativistic particle acceleration is likely an ongoing process; the infrared emitting electrons were accelerated no longer than ~80 years ago. There is a small but significant concave curvature to the spectrum, as expected if the accelerating shocks have been modified by the back pressure of the cosmic rays; given calibration uncertainties, this conclusion must be considered tentative at present. The 2.2 micron polarization angles and the emission-line filaments observed by HST are both offset from the local radial direction by 10 - 20 degrees, providing evidence that the magnetic fields in Cas A are generated by Rayleigh-Taylor instabilities in the decelerating ejecta.

Cosmological hydrodynamic simulations of large scale structure in the universe have shown that accretion shocks and merger shocks form due to flow motions associated with the gravitational collapse of nonlinear structures. Estimated speed and curvature radius of these shocks could be as large as a few 1000 km/s and several Mpc, respectively. According to the diffusive shock acceleration theory, populations of cosmic-ray particles can be injected and accelerated to very high energy by astrophysical shocks in tenuous plasmas. In order to explore the cosmic ray acceleration at the cosmic shocks, we have performed nonlinear numerical simulations of cosmic ray (CR) modified shocks with the newly developed CRASH (Cosmic Ray Amr SHock) numerical code. We adopted the Bohm diffusion model for CRs, based on the hypothesis that strong Alfv\'en waves are self-generated by streaming CRs. The shock formation simulation includes a plasma-physics-based ``injection'' model that transfers a small proportion of the thermal proton flux through the shock into low energy CRs for acceleration there. We found that, for strong accretion shocks, CRs can absorb most of shock kinetic energy and the accretion shock speed is reduced up to 20 %, compared to pure gas dynamic shocks. For merger shocks with small Mach numbers, however, the energy transfer to CRs is only about 10-20 % with an associated CR particle fraction of $10^{-3}$. Nonlinear feedback due to the CR pressure is insignificant in the latter shocks. Although detailed results depend on models for the particle diffusion and injection, these calculations show that cosmic shocks in large scale structure could provide acceleration sites of extragalactic cosmic rays of the highest energy.

In this lecture I outline some of the underlying physics issues associated with accelerators plausibly capable of explaining the UHECRs up to ZeV energies. I concentrate on the concentrate on mechanisms and their constraints, but provide a brief background on on observations and the constraints they supply, as well.

A numerical scheme that incorporates a thermal leakage injection model into a combined gas dynamics and cosmic ray (CR, hereafter) diffusion-convection code has been developed. The particle injection is followed numerically by filtering the diffusive flux of suprathermal particles across the shock to the upstream region according to a velocity-dependent transparency function that controls the fraction of leaking particles. We have studied CR injection and acceleration efficiencies during the evolution of CR modified planar shocks for a wide range of initial shock Mach numbers, $M_0$, assuming a Bohm-like diffusion coefficient. The injection process is very efficient when the subshock is strong, leading to fast and significant modification of the shock structure. As the CR pressure increases, the subshock weakens and the injection rate decreases accordingly, so that the subshock does not disappear. Although some fraction of the particles injected early in the evolution continue to be accelerated to ever higher energies, the postshock CR pressure reaches an approximate time-asymptotic value due to a balance between fresh injection/acceleration and advection/diffusion of the CR particles away from the shock. We conclude that the injection rates in strong parallel shocks are sufficient to lead to rapid nonlinear modifications to the shock structures and that self-consistent injection and time-dependent simulations are crucial to understanding the non-linear evolution of CR modified shocks.

Shocks are a ubiquitous consequence of cosmic structure formation, and they play an essential role in heating galaxy cluster media. Virtually all of the gas in clusters has been processed by one or more shocks of at least moderate strength. These are collisionless shocks, so likely sites for diffusive shock acceleration of high energy particles. We have carried out numerical simulations of cosmic structure formation that directly include acceleration and transport of nonthermal protons, as well as primary and secondary electrons. Nonthermal emissions have also been computed from the resulting particle spatial and energy distributions. Here we outline some of our current findings, showing that nonthermal protons may contribute a significant pressure in cluster media, and that expected radio, X-ray and $\gamma$-ray emissions from these populations should be important cluster diagnostics.

The theory of shock acceleration predicts the maximum particle energy to be limited only by the acceleration time and the size (geometry) of the shock. This led to optimistic estimates for the galactic cosmic ray energy achievable in the SNR shocks. The estimates imply that the accelerated particles, while making no strong impact on the shock structure (test particle approach) are nevertheless scattered by strong self-generated Alfven waves (turbulent boost) needed to accelerate them quickly. We demonstrate that these two assumptions are in conflict when applied to SNRs of the age required for cosmic ray acceleration to the ``knee'' energy. We study the combined effect of acceleration nonlinearity (shock modification by accelerated particles) and wave generation on the acceleration process. We show that the refraction of self-generated waves resulting from the deceleration of the plasma flow by the pressure of energetic particles causes enhanced losses of these particles. This effect slows down the acceleration and changes the shape of particle spectrum near the cut-off. The implications for observations of TeV emission from SNR remnants are also discussed.

We study the generation and distribution of high energy electrons in cosmic environment and their observational consequences by carrying out the first cosmological simulation that includes directly cosmic ray (CR) particles. Starting from cosmological initial conditions we follow the evolution of primary and secondary electrons (CRE), CR ions (CRI) and a passive magnetic field. CRIs and primary CREs are injected and accelerated at large scale structure shocks. Secondary CREs are continuously generated through inelastic p-p collisions. We include spatial transport, adiabatic expansion/compression, Coulomb collisions, bremsstrahlung, synchrotron (SE)and inverse Compton (IC) emission. We find that, from the perspective of cosmic shock energy and acceleration efficiency, the few detections of hard X-ray radiation excess could be explained in the framework of IC emission of primary CREs in clusters undergoing high accretion/merger phase. Instead, IC emission from both primary and secondary CREs accounts at most for a small fraction of the radiation excesses detected in the extreme-UV (except for the Coma cluster as reported by Bowyer et al.1999). Next, we calculate the SE after normalizing the magnetic field so that for a Coma-like cluster <B^2>^1/2~3 \muG. Our results indicate that the SE from secondary CREs reproduces several general properties of radio halos, including the recently found P_1.4GHz vs T relation, the morphology and polarization of the emitting region and, to some extent, the spectral index. Moreover, SE from primary CREs turns out sufficient to power extended regions resembling radio relics observed at the outskirts of clusters. Again we find striking resemblance between morphology, polarization and spectral index of our synthetic maps and those reported in the literature.

The linear theory of shock acceleration predicts the maximum particle energy to be limited only by the acceleration time and the size of the shock. We study the combined effect of acceleration nonlinearity (shock modification by accelerated particles, that must be present in strong astrophysical shocks) and propagation of Alfven waves that are responsible for particle confinement to the shock front. We show that wave refraction to larger wave numbers in the nonlinearly modified flow causes enhanced losses of particles in the momentum range p_max/R < p < p_max, where R>1 is the nonlinear pre-compression of the flow and p_max is a conventional maximum momentum, that could be reached if there was no refraction.

We investigate the production of cosmic ray (CR) protons at cosmological shocks by performing, for the first time, numerical simulations of large scale structure formation that include directly the acceleration, transport and energy losses of the high energy particles. CRs are injected at shocks according to the thermal leakage model and, thereafter, accelerated to a power-law distribution as indicated by the test particle limit of the diffusive shock acceleration theory. The evolution of the CR protons accounts for losses due to adiabatic expansion/compression, Coulomb collisions and inelastic p-p scattering. Our results suggest that CR protons produced at shocks formed in association with the process of large scale structure formation could amount to a substantial fraction of the total pressure in the intra-cluster medium. Their presence should be easily revealed by GLAST through detection of gamma-ray flux from the decay of neutral pions produced in inelastic p-p collisions of such CR protons with nuclei of the intra-cluster gas. This measurement will allow a direct determination of the CR pressure contribution in the intra-cluster medium. We also find that the spatial distribution of CR is typically more irregular than that of the thermal gas because it is more influenced by the underlying distribution of shocks. This feature is reflected in the appearance of our gamma-ray synthetic images. Finally, the average CR pressure distribution appears statistically slightly more extended than the thermal pressure.

We present the first three-dimensional MHD radio galaxy simulations that explicitly model transport of relativistic electrons, including diffusive acceleration at shocks as well as radiative and adiabatic cooling in smooth flows. We discuss three simulations of light Mach 8 jets, designed to explore the effects of shock acceleration and radiative aging on the nonthermal particle populations that give rise to synchrotron and inverse-Compton radiations. We also conduct detailed synthetic radio observations of our simulated objects. We have gained several key insights from this approach: 1. The jet head in these multidimensional simulations is extremely complex. The classical jet termination shock is often absent, but motions of the jet terminus spin a ``shock-web complex'' within the backflowing jet material of the head. 2. Understanding the spectral distribution of energetic electrons in these simulations relies partly upon understanding the shock-web complex, for it can give rise to distributions that confound interpretation in terms of the standard model for radiative aging of radio galaxies. 3. The magnetic field outside of the jet itself becomes very intermittent and filamentary in these simulations, yet adiabatic expansion causes most of the cocoon volume to be occupied by field strengths considerably diminished below the nominal jet value. Thus population aging rates vary considerably from point to point.

We investigate the effects of rotation on the evolution of the Parker instability in an exponentially-stratified and uniformly-rotating magnetized gas disk by carrying out three-dimensional numerical simulations with an isothermal magnetohydrodynamic code. The instability has been initialized by random velocity perturbations. In the linear stage, the evolution is not much different from that without rotation and the mixed (undular + interchange) mode regulates the system. The interchange mode induces alternating dense and rarefied regions with small radial wavelengths, while the undular mode bends the magnetic field lines in the plane of azimuthal and vertical directions. In the nonlinear stage, flow motion overall becomes chaotic as in the case without rotation. However, as the gas in higher positions slides down along field lines forming supersonic flows, the Coriolis force becomes important. As oppositely directed flows fall into valleys along both sides of magnetic field lines, they experience the Coriolis force toward opposite directions, which twists magnetic field lines there. Hence, we suggest that the Coriolis force plays a role in randomizing magnetic field. The three-dimensional density structure formed by the instability is still sheet-like with the short dimension along the radial direction, as in the case without rotation. However, the long dimension is now slightly tilted with respect to the mean field direction. The shape of high density regions is a bit rounder. The maximum enhancement factor of the vertical column density relative to its initial value is about 1.5, which is smaller than that in the case without rotation. We conclude that uniform rotation doesn't change our point of view that the Parker instability alone is not a viable mechanism for the formation of giant molecular clouds.

We report the localization of the forward and reversed shock fronts in the young supernova remnant Cas-A using X-ray data obtained with the Chandra Observatory. High resolution X-ray maps resolve a previously unseen X-ray feature encompassing the extremity of the remnant. This feature consists of thin, tangential wisps of emission bordering the outer edge of the thermal X-ray and radio remnant, forming a circular rim, approx. 2.7 in radius. Radio images show a sharp rise in brightness at this X-ray rim, along with a large jump in the synchrotron polarization angle. These characteristics suggest that these wisps are the previously unresolved signature of the forward, or outer, shock. Similarly, we identify the sharp rise in emissivity of the bright shell for both the radio and X-ray line emission associated with the reverse shock. The derived ratio of the averaged forward and reverse shock radii of approx. 3:2 constrains the remnant to have swept up roughly the same amount of mass as was ejected; this suggests that Cas-A is just entering the Sedov phase. Comparison of the X-ray spectra from the two shock regions shows that the equivalent widths of prominent emission lines are significantly lower exterior to the bright shell, as expected if they are respectively identified with the shocked circumstellar material and shocked ejecta. Furthermore, the spectrum of the outer rim itself is dominated by power-law emission, likely the counterpart of the non-thermal component previously seen at energies above 10 keV.

Cosmic-rays are ubiquitous, but their origins are surprisingly difficult to understand. A review is presented of some of the basic issues common to cosmic particle accelerators and arguments leading to the likely importance of diffusive shock acceleration as a general explanation. The basic theory of diffusive shock acceleration is outlined, followed by a discussion of some of the key issues that still prevent us from a full understanding of its outcomes. Some recent insights are mentioned at the end that may help direct ultimate resolution of our uncertainties.

We report on an effort to study the connections between dynamics in simulated radio galaxy plasma flows and the properties of nonthermal electron populations carried in those flows. To do this we have introduced a new numerical scheme for electron transport that allows a much more detailed look at this problem than has been possible before. Especially when the dynamics is fully three dimensional the flows are generally chaotic in the cocoon, and the jet itself can flail about violently. The bending jet can pinch itself off and redirect itself to enhance its penetration of the ambient medium. These behaviors often eliminate the presence of a strong jet termination shock, which is assumed present in all modern cartoon models of the RG phenomenon. Instead a much more complex ``shock web'' forms near the end of the jet that leads to a far less predictable pattern of particle acceleration. Similarly, the magnetic fields in these flows are highly filamented, as well as spatially and temporally intermittent. This leads to a very localized and complex pattern of synchrotron aging for relativistic electron populations, which makes it difficult to use properties of the electron spectrum to infer the local rate of aging.

We perform a series of so-called ``synthetic observations'' on a set of 3D MHD jet simulations which explicitly include energy-dependent transport of relativistic electrons, as described in the companion paper by Jones, Tregillis, & Ryu. Analyzing them in light of the complex source dynamics and energetic particle distributions described in that paper, we find that the standard model for radiative aging in radio galaxies does not always adequately reflect the detailed source structure.

We have developed a novel computer code designed to follow the evolution of cosmic-ray modified shocks, including the full momentum dependence of the particles for a realistic diffusion coefficient model. In this form the problem is technically very difficult, because one needs to cover a wide range of diffusive scales, beginning with those slightly larger than the physical shock thickness. With most finite difference schemes for Euler's equations the numerical shock thickness is at least one zone across, so this provides a lower bound on the physical scale for diffusive transport computation. Our code uses sub-zone shock tracking and multi-level adaptive mesh refinement to provide enhanced spatial resolution around shocks at modest cost compared to the coarse grid and vastly improved cost effectiveness compared to a uniform, highly refined grid. We present and discuss the implications from our initial results.

One of the key questions to understanding the efficiency of diffusive shock acceleration of the cosmic rays (CRs) is the injection process from thermal particles. A self-consistent injection model based on the interactions of the suprathermal particles with self-generated magneto-hydrodynamic waves has been developed recently by Malkov (1998, Phys. Rev. E 58, 4911). By adopting this analytic solution, a numerical treatment of the plasma-physical injection model at a strong quasi-parallel shock has been devised and incorporated into the combined gas dynamics and the CR diffusion-convection code. In order to investigate self-consistently the injection and acceleration efficiencies, we have applied this code to the CR modified shocks of both high and low Mach numbers (M=30 and M=2.24) with a Bohm type diffusion model. We find the injection process is self-regulated in such a way that the injection rate reaches and stays at a nearly stable value after quick initial adjustment. For both shocks about 0.1% of the incoming thermal particles are injected into the CRs. For the weak shock, the shock has reached a steady state within our integration time and about 10% of the total available shock energy is transfered into the CR energy density. The strong shock has achieved a higher acceleration efficiency of about 20% by the end of our simulation, but has not yet reached a steady-state.

We review the possible roles of large scale shocks as particle accelerators in clusters of galaxies. Recent observational and theoretical work has suggested that high energy charged particles may constitute a substantial pressure component in clusters. If true that would alter the expected dynamical evolution of clusters and increase the dynamical masses consistent with hydrostatic equilibrium. Moderately strong shocks are probably common in clusters, through the actions of several agents. The most obvious of these agents include winds from galaxies undergoing intense episodes of starbursts, active galaxies and cosmic inflows, such as accretion and cluster mergers. We describe our own work derived from simulations of large scale structure formation, in which we have, for the first time, explicitly included passive components of high energy particles. We find, indeed that shocks associated with these large scale flows can lead to nonthermal particle pressures big enough to influence cluster dynamics. These same simulations allow us also to compute nonthermal emissions from the clusters. Here we present resulting predictions of gamma-ray fluxes.

Imaging polarimetry at 1.65 and 2.2 $\mu $m is presented for the classic starburst galaxy M82 and the advanced merger system Cen A. Polarimetry at near IR wavelengths allows the magnetic field geometry in galaxies to be probed much deeper into dusty regions than optical polarimetry. In M82, the magnetic field throughout the nucleus has a polar geometry, presumably due to the massive vertical flow that is a result of the intense star formation there. Fully two thirds of the line of sight dust through to the center of M82 contains a vertical magnetic field. In Cen A, the prominent dust lane shows a normal planar field geometry. There is no indication of significant disturbance in the field geometry in the dust lane and the polarization strength is near normal for the amount of extinction. Either the magnetic field geometry was well maintained during the merger, or it reestablished itself very easily.

We investigate through high resolution 3D simulations the nonlinear evolution of compressible magnetohydrodynamic flows subject to the Kelvin-Helmholtz instability. We confirm in 3D flows the conclusion from our 2D work that even apparently weak magnetic fields embedded in Kelvin-Helmholtz unstable plasma flows can be fundamentally important to nonlinear evolution of the instability. In fact, that statement is strengthened in 3D by this work, because it shows how field line bundles can be stretched and twisted in 3D as the quasi-2D Cat's Eye vortex forms out of the hydrodynamical motions. In our simulations twisting of the field may increase the maximum field strength by more than a factor of two over the 2D effect. If, by these developments, the Alfv\'en Mach number of flows around the Cat's Eye drops to unity or less, our simulations suggest magnetic stresses will eventually destroy the Cat's Eye and cause the plasma flow to self-organize into a relatively smooth and apparently stable flow that retains memory of the original shear. For our flow configurations the regime in 3D for such reorganization is $4\lesssim M_{Ax} \lesssim 50$, expressed in terms of the Alfv\'en Mach number of the original velocity transition and the initial Alfv\'en speed projected to the flow plan. For weaker fields the instability remains essentially hydrodynamic in early stages, and the Cat's Eye is destroyed by the hydrodynamic secondary instabilities of a 3D nature. Then, the flows evolve into chaotic structures that approach decaying isotropic turbulence. In this stage, there is considerable enhancement to the magnetic energy due to stretching, twisting, and turbulent amplification, which is retained long afterwards. The magnetic energy eventually catches up to the kinetic energy, and the nature of flows become magnetohydrodynamic.

We present results from three-dimensional (3-D) numerical simulations investigating the magnetohydrodynamics of cloud-wind interactions. The initial cloud is spherical while the magnetic field is uniform and transverse to the cloud motion. A simplified analytical model that describes the magnetic energy evolution in front of the cloud is developed and compared with simulation results. In addition, it is found the interaction of the cloud with a magnetized interstellar medium (ISM) results in the formation of a highly structured magnetotail. The magnetic flux in the wake of the cloud organizes into flux ropes and a reconnection, current sheet is developed, as field lines of opposite polarity are brought close together near the symmetry axis. At the same time, magnetic pressure is strongly enhanced at the leading edge of the cloud from the stretching of the field lines that occurs there. This has an important dynamical effect on the subsequent evolution of the cloud, since some unstable modes tend to be strongly enhanced.

We investigate the dynamical importance of a newly recognized possible source of significant feedback generated during structure formation; namely cosmic ray (CR) pressure. We present evidence for the existence of numerous shocks in the hot gas of galaxy clusters (GCs). We employ for the first time an explicit numerical treatment of CR acceleration and transport in hydro simulations of structure formation. According to our results, CRs provide an important fraction of the total pressure inside GCs, up to several tenths. This was true even at high redshift (z=2), meaning that such non-thermal component could affect the evolution of structure formation.

We have examined the properties of shock waves in simulations of large scale structure formation for two cosmological scenarios (a SCDM and a LCDM with Omega =1). Large-scale shocks result from accretion onto sheets, filaments and Galaxy Clusters (GCs) on a scale of circa 5 Mpc/h in both cases. Energetic motions, both residual of past accretion history and due to current asymmetric inflow along filaments, generate additional, common shocks on a scale of about 1 Mpc/h, which penetrate deep inside GCs. Also collisions between substructures inside GCs form merger shocks. Consequently, the topology of the shocks is very complex and highly connected. During cosmic evolution the comoving shock surface density decreases, reflecting the ongoing structure merger process in both scenarios. Accretion shocks have very high Mach numbers (10-10^3), when photo-heating of the pre-shock gas is not included. The typical shock speed is of order v_{sh}(z) =H(z)lambda_{NL}(z), with lambda_{NL}(z) the wavelength scale of the nonlinear perturbation at the given epoch. However, the Mach number for shocks occuring within clusters is usually smaller (3-10), due to the fact that the intracluster gas is already hot. Statistical fits of shock speed around GCs as a function of GCs temperature give power-law's in accord with 1-D predictions. However, a very different result is obtained for fits of the shock radius, reflecting the very complex shock structures forming in 3-D simulations. The in-flowing kinetic energy across such shocks, giving the power available for cosmic-ray acceleration, is comparable to the cluster X-ray luminosity emitted from a central region of radius 0.5 Mpc/h. Considering their large size and long lifetimes, those shocks are potentially interesting sites for cosmic-ray acceleration, if modest magnetic fields exist within them.

Numerical results for particle acceleration at multiple oblique shocks are presented. We calculate the steady state spectral slope of test particles accelerated by the first order Fermi process. The results are compared to analytical treatments, for parameters, where the diffusion approximation does apply. Effects of injection and finite shock extend are included phenomenologically. We find the spectrum of accelerated particles to harden substantially at multiple oblique shocks and discuss the influence of the number of shocks compared to the obliquity itself.

We present simulations of the propagation of magnetized jets. This work differs from previous studies in that the cross-sectional distributions of the jets's state variables are derived from analytical models for magneto-centrifugal launching. The source is a magnetized rotator whose properties are specfied as boundary conditions. The jets in these simulations are considerably more complex than the ``top-hat''constant density etc. profiles used in previous work. We find that density and magnetic field stratification (with radius) in the jet leads to new behavior including the separation of an inner jet core from a low density collar. We find this {\it jet within a jet} structure, along with the magnetic stresses, leads to propagation behaviors not observed in previous simulation studies. Our methodology allows us to compare MHD jets from different types of sources whose properties could ultimately be derived from the behavior of the propagating jets.

A cosmological hydrodynamic code is described, which includes a routine to compute cosmic ray acceleration and transport in a simplified way. The routine was designed to follow explicitly diffusive acceleration at shocks, and second-order Fermi acceleration and adiabatic loss in smooth flows. Synchrotron cooling of the electron population can also be followed. The updated code is intended to be used to study the properties of nonthermal synchrotron emission and inverse Compton scattering from electron cosmic rays in clusters of galaxies, in addition to the properties of thermal bremsstrahlung emission from hot gas. The results of a test simulation using a grid of $128^3$ cells are presented, where cosmic rays and magnetic field have been treated passively and synchrotron cooling of cosmic ray electrons has not been included.

We investigate the nonlinear effects of uniform rotation on the Parker instability in an exponentially-stratified disk through high-resolution simulations. During the linear stage, the speed of gas motion is subsonic and the evolution with the rotation is not much different from that without the rotation. This is because the Coriolis force is small. During the nonlinear stage, oppositely-directed supersonic flows near a magnetic valley are under the influence of the Coriolis force with different directions, resulting in twisted magnetic field lines near the valley. Sheet-like structures, which are tilted with respect to the initial field direction, are formed with an 1.5 enhancement of column density with respect to its initial value. Even though uniform rotation doesn't give much impact on density enhancement, it generates helically twisted field lines, which may become an additional support mechanism of clouds.

We present results from the first three-dimensional numerical simulations of moderately supersonic cloud motion through a tenuous, magnetized medium. We show that the interaction of the cloud with a magnetic field perpendicular to its motion has a great dynamical impact on the development of instabilities at the cloud surface. Even for initially spherical clouds, magnetic field lines become trapped in surface deformations and undergo stretching. The consequent field amplification that occurs there and particularly its variation across the cloud face then dramatically enhance the growth rate of Rayleigh-Taylor unstable modes, hastening the cloud disruption.

We have carried out simulations of the nonlinear evolution of the magnetohydrodynamic (MHD) Kelvin-Helmholtz (KH) instability for compressible fluids in $2\frac{1}{2}$-dimensions, extending our previous work by Frank et al (1996) and Jones \etal (1997). In the present work we have simulated flows in the x-y plane in which a ``sheared'' magnetic field of uniform strength ``smoothly'' rotates across a thin velocity shear layer from the z direction to the x direction, aligned with the flow field. We focus on dynamical evolution of fluid features, kinetic energy dissipation, and mixing of the fluid between the two layers, considering their dependence on magnetic field strength for this geometry. The introduction of magnetic shear can allow a Cat's Eye-like vortex to form, even when the field is stronger than the nominal linear instability limit given above. For strong fields that vortex is asymmetric with respect to the preliminary shear layer, however, so the subsequent dissipation is enhanced over the uniform field cases of comparable field strength. In fact, so long as the magnetic field achieves some level of dynamical importance during an eddy turnover time, the asymmetries introduced through the magnetic shear will increase flow complexity, and, with that, dissipation and mixing. The degree of the fluid mixing between the two layers is strongly influenced by the magnetic field strength. Mixing of the fluid is most effective when the vortex is disrupted by magnetic tension during transient reconnection, through local chaotic behavior that follows.

We have applied an effective numerical scheme for cosmic-ray transport to 3D MHD simulations of jet flow in radio galaxies (see the companion paper by Jones et al. 1999). The marriage of relativistic particle and 3D magnetic field information allows us to construct a rich set of ``synthetic observations'' of our simulated objects. The information is sufficient to calculate the ``true'' synchrotron emissivity at a given frequency using explicit information about the relativistic electrons. This enables us to produce synchrotron surface-brightness maps, including polarization. Inverse-Compton X-ray surface-brightness maps may also be produced. First results intended to explore the connection between jet dynamics and electron transport in radio lobes are discussed. We infer lobe magnetic field values by comparison of synthetically observed X-ray and synchrotron fluxes, and find these ``inverse-Compton'' fields to be quite consistent with the actual RMS field averaged over the lobe. The simplest minimum energy calculation from the synthetic observations also seems to agree with the actual simulated source properties.

We have developed an economical, effective numerical scheme for cosmic-ray transport suitable for treatment of electrons up to a few hundreds of GeV in multidimensional simulations of radio galaxies. The method follows the electron population in sufficient detail to allow computation of synthetic radio and X-ray observations of the simulated sources, including spectral properties (see the companion paper by Tregillis et al. 1999). The cosmic-ray particle simulations can follow the effects of shock acceleration, second-order Fermi acceleration as well as radiative and adiabatic energy losses. We have applied this scheme to 2-D and 3-D MHD simulations of jet-driven flows and have begun to explore links between dynamics and the properties of high energy electron populations in radio lobes. The key initial discovery is the great importance to the high energy particle population of the very unsteady and inhomogeneous flows, especially near the end of the jet. Because of this, in particular, our simulations show that a large fraction of the particle population flowing from the jet into the cocoon never passes through strong shocks. The shock strengths encountered are not simply predicted by 1-D models, and are quite varied. Consequently, the emergent electron spectra are highly heterogeneous. Rates of synchrotron aging in "hot-spots" seem similarly to be very uneven, enhancing complexity in the spectral properties of electrons as they emerge into the lobes and making more difficult the task of comparing dynamical and radiative ages.

The properties of cosmic shock waves are studied through numerical simulations in two cosmological scenarios (SCDM and LCDM). The scaling relations for the average radius and velocity associated with the accretion shocks are somewhat different, yet qualitatively similar to the self similar solutions for a flat Omega_M=1 universe. The energy supplied by infalling gas at accretion shock waves is large enough to sustain production of abundant cosmic ray populations if a viable acceleration mechanism can take place there. Finally, in addition to shocks created by the encounter of the merging ICMs of two clusters of galaxies, accretion shocks associated with the merging clusters generate strong ``relic'' shocks which propagate through the ICM producing additional heating of the ICM, and associated CR acceleration.

The theory of diffusive particle acceleration explains the spectral properties of the cosmic rays below energies of approx. 10^6 GeV as produced at strong shocks in supernova remnants (SNR's). To supply the observed flux of cosmic rays, a significant fraction of the energy released by a supernova has to be transfered to cosmic rays. The key to the question of the efficiency of SNR's in producing cosmic rays is the injection process from thermal energies. A self-consistent model has to take into account the interaction of the accelerated particles with magneto-hydrodynamic waves, which generate the particle diffusion, a requisite for the acceleration process. Such a nonlinear model of the turbulent background plasma has been developed recently (Malkov, 1998, Phys. Rev. E 58, 4911). We use this model for the first numerical treatment of the gas dynamics and the diffusion-convection equation at a quasi-parallel strong shock, which incorporates a plasma-physical injection model to investigate the cosmic ray production.

We investigate, through 2D MHD numerical simulations, the interaction of a uniform magnetic field oblique to a moving interstellar cloud. In particular we explore the transformation of cloud kinetic energy into magnetic energy as a result of field line stretching. Some previous simulations have emphasized the possible dynamical importance of a ``magnetic shield'' formed around clouds when the magnetic field is perpendicular to the cloud motion (Jones et al. 1996, Miniati et al. 1998). It was not clear, however, how dependent those findings were to the assumed field configuration and cloud properties. To expand our understanding of this effect, we examine several new cases by varing the magnetic field orientation angle with respect to the cloud motion (\theta), the cloud-background density contrast, and the cloud Mach number. We show that in 2D and with \theta large enough, the magnetic field tension can become dominant in the dynamics of the motion of high density contrast, low Mach number clouds. In such cases a significant fraction of cloud kinetic energy can be transformed into magnetic energy with the magnetic pressure at the cloud nose exceeding the ram pressure of the impinging flow. We derive a characteristic timescale for this process of energy ``conversion''. We find also that unless the cloud motion is highly aligned to the magnetic field, reconnection through tearing mode instabilities in the cloud wake limit the formation of a strong flux rope feature following the cloud. Finally we attempt to interpret some observational properties of the magnetic field in view of our results.

We present a multi-dimensional numerical code to solve isothermal magnetohydrodynamic (IMHD) equations for use in modeling astrophysical flows. First, we have built a one-dimensional code which is based on an explicit finite-difference method on an Eulerian grid, called the total variation diminishing (TVD) scheme. Recipes for building the one-dimensional IMHD code, including the normalized right and left eigenvectors of the IMHD Jacobian matrix, are presented. Then, we have extended the one-dimensional code to a multi-dimensional IMHD code through a Strang-type dimensional splitting. In the multi-dimensional code, an explicit cleaning step has been included to eliminate non-zero $\nabla\cdot B$ at every time step. To estimate the proformance of the code, one- and two-dimensional IMHD shock tube tests, and the decay test of a two-dimensional Alfv\'{e}n wave have been done. As an example of astrophysical applications, we have simulated the nonlinear evolution of the two-dimensional Parker instability under a uniform gravity.

We introduce a simple and economical but effective method for including relativistic electron transport in multi-dimensional simulations of radio galaxies. The method is designed to follow explicitly diffusive acceleration at shocks, and, in smooth flows 2nd order Fermi acceleration plus adiabatic and synchrotron cooling. We are able to follow both the spatial and energy distributions of the electrons, so that direct synchrotron emission properties can be modeled in time-dependent flows for the first time. Here we present first results in the form of some axis-symmetric MHD simulations of Mach 20 light jet flows. These show clearly the importance of nonsteady terminal shocks that develop in such flows even when the jet inflow is steady. As a result of this and other consequences of the fundamentally driven character of jets, we find complex patterns of emissivities and synchrotron spectra, including steep spectral gradients in hot spots, islands of distinct spectra electrons within the lobes and spectral gradients coming from the dynamical histories of a given flow element rather than from synchrotron aging of the embedded electrons. In addition, spectral aging in the lobes tends to proceed more slowly than one would estimate from regions of high emissivity.

We present two-dimensional MHD simulations of the evolution of a young Type Ia supernova remnant during its interaction with an interstellar cloud of comparable size at impact. We include for the first time in such simulations explicit relativistic electron transport, including spectral information using a simple but effective scheme that follows their acceleration at shocks and subsequent transport. From this information we also model radio synchrotron emission, including spectra. The principal conclusions from these experiments are: 1) Independent of the cloud interaction, the SNR reverse shock can be an efficient site for particle acceleration in a young SNR. 2) At these early times the synchrotron spectral index due to electrons accelerated at the primary shocks should be close to 0.5 unless those shocks are modified by cosmic-ray pressures. However, interaction with the cloud generates regions of distinctly steeper spectra, which may complicate interpretation in terms of global dynamical models for SNR evolution. 3) The internal motions within the SNR become highly turbulent following the cloud interaction. 4) An initially uniform interstellar magnetic field is preferentially amplified along the magnetic equator of the SNR, primarily due to biased amplification by instabilities. Independent of the external field configuration, there is a net radial direction to this field inside the SNR. 5) Filamentary radio structures correlate well with magnetic filaments, while diffuse emission follows the electron distribution. 6) Interaction with the cloud enhances both the electron population and the radio emission.

Using an isothermal MHD code, we have performed three-dimensional, high-resolution simulations of the Parker instability. The initial equilibrium system is composed of exponentially-decreasing isothermal gas and magnetic field (along the azimuthal direction) under a uniform gravity. The evolution of the instability can be divided into three phases: linear, nonlinear, and relaxed. During the linear phase, the perturbations grow exponentially with a preferred scale along the azimuthal direction but with smallest possible scale along the radial direction, as predicted from linear analyses. During the nonlinear phase, the growth of the instability is saturated and flow motion becomes chaotic. Magnetic reconnection occurs, which allows gas to cross field lines. This, in turn, results in the redistribution of gas and magnetic field. The system approaches a new equilibrium in the relaxed phase, which is different from the one seen in two-dimensional works. The structures formed during the evolution are sheet-like or filamentary, whose shortest dimension is radial. Their maximum density enhancement factor relative to the initial value is less than 2. Since the radial dimension is too small and the density enhancement is too low, it is difficult to regard the Parker instability alone as a viable mechanism for the formation of giant molecular clouds.

We extend previous studies of the physics of interstellar cloud collisions by beginning investigation of the role of magnetic fields through 2D magnetohydrodynamic (MHD) numerical simulations. We study head-on collisions between equal mass, mildly supersonic diffuse clouds. We include a moderate magnetic field and two limiting field geometries, with the field lines parallel (aligned) and perpendicular (transverse) to the colliding cloud motion. We explore both adiabatic and radiative cases, as well as symmetric and asymmetric ones. We also compute collisions between clouds evolved through prior motion in the intercloud medium and compare with unevolved cases. We find that: In the (i) aligned case, adiabatic collisions, like their HD counterparts, are very disruptive, independent of the cloud symmetry. However, when radiative processes are taken into account, partial coalescence takes place even in the asymmetric case, unlike the HD calculations. In the (ii) transverse case, collisions between initially adjacent unevolved clouds are almost unaffected by magnetic fields. However, the interaction with the magnetized intercloud gas during the pre-collision evolution produces a region of very high magnetic energy in front of the cloud. In collisions between evolved clouds with transverse field geometry, this region acts like a ``bumper'', preventing direct contact between the clouds, and eventually reverses their motion. The ``elasticity'', defined as the ratio of the final to the initial kinetic energy of each cloud, is about 0.5-0.6 in the cases we considered. This behavior is found both in adiabatic and radiative cases.

A description is given for preserving ${\bmsy\nabla}\cdot{\vec B}=0$ in a magnetohydrodynamic (MHD) code that employs the upwind, Total Variation Diminishing (TVD) scheme and the Strang-type operator splitting for multi-dimensionality. The method is based on the staggered mesh technique to constrain the transport of magnetic field: the magnetic field components are defined at grid interfaces with their advective fluxes on grid edges, while other quantities are defined at grid centers. The magnetic field at grid centers for the upwind step is calculated by interpolating the values from grid interfaces. The advective fluxes on grid edges for the magnetic field evolution are calculated from the upwind fluxes at grid interfaces. Then, the magnetic field can be maintained with ${\bmsy\nabla}\cdot{\vec B}=0$ exactly, if this is so initially, while the upwind scheme is used for the update of fluid quantities. The correctness of the code is demonstrated through tests comparing numerical solutions either with analytic solutions or with numerical solutions from the code using an explicit divergence-cleaning method. Also the robustness is shown through tests involving realistic astrophysical problems.

We have studied head-on collisions between equal-mass, mildly supersonic (Mach number 1.5) HI clouds, in a standard Two-phase ISM (T_cl = 74 K, n_cl = 22 cm^-3, \chi = 100). We explore the role of various factors, including the radiative cooling parameter \eta = t_rad/t_coll (t_coll=R_c/v_c), evolutionary modifications on the cloud structure (by colliding clouds ``evolved'' through independent motion within the intercloud medium (WIM)), and the symmetry of the problem (by colliding initially identical clouds, evolved to different ages before impact). The presence of bow shocks and ram pressure from material in the cloud wake, developed during such evolution through the WIM, significantly alters these interactions with respect to the standard case of non-evolved clouds. In general, in our adiabatic collisions the clouds are disrupted and convert their gas into a few low density contrast clumps. By contrast, for symmetric radiative cases we find that the two clouds coalesce, with almost all the initial kinetic energy radiated away. On the other hand, for both adiabatic and radiative collisions, asymmetric collisions have a much greater tendency to disrupt the two clouds. Fragmentation of the clouds may occur, and instabilities are in general enhanced. In addition, radiative cooling is less efficient in our asymmetric interactions, so that those parts of the clouds that initially seem to merge are more likely to re-expand and fade into the WIM. Since the majority of real cloud collisions should be asymmetric for one reason or another, we conclude that most gasdynamical diffuse cloud collisions will be disruptive, at least in the absence of significant self-gravity or of a significant magnetic field.

We present multi-dimensional simulations of magnetized radiative jets appropriate to Young Stellar Objects. Magnetized jets subject to collisionally excited radiative losses have not, as yet, received extensive scrutiny. The purpose of this letter is to articulate the propagation dynamics of radiative MHD jets in the context of the extensive jet literature. Most importantly, we look for morphological and kinematic diagnostics that may distinguish hydrodynamic protostellar jets from their magnetically dominated cousins. Our simulations are axisymmetric (2.5-D). A toroidal B_phi field geometry is used. Our models have high sonic Mach numbers (Ms approx 10), but lower fast mode Mach number (Mf approx 5). This is approximately the case for jets formed via disk-wind or X-wind models - currently the consensus choice for launching and collimating YSO jets. Time-dependent radiative losses are included via a coronal cooling curve. Our results demonstrate that the morphology and propagation characteristics of strongly magnetized radiative jets can differ significantly from jets with weak fields. In particular the formation of nose-cones via post-shock hoop stresses leads to narrow bow shocks and enhanced bow shock speeds. In addition, the hoop stresses produce strong shocks in the jet beam which constrasts with the relatively unperturbed beam in radiative hydrodynamic jets. Our simulations show that pinch modes produced by magnetic tension can strongly effect magnetized protostellar jets.

This paper reports on a workshop hosted by the University of Minnesota, March 23-26, 1997. It addressed fundamental dynamical issues associated with the evolution of shell supernova remnants and the relationships between supernova remnants and their environments. The workshop considered, in addition to classical shell SNRs, dynamical issues involving X-ray filled composite remnants and pulsar driven shells, such as that in the Crab Nebula. Approximately 75 participants with wide ranging interests attended the workshop. An even larger community helped through extensive on-line debates prior to the meeting. Each of the several sessions, organized mostly around chronological labels, also addressed some underlying, general physical themes: How are SNR dynamics and structures modified by the character of the CSM and the ISM and vice versa? How are magnetic fields generated in SNRs and how do magnetic fields influence SNRs? Where and how are cosmic-rays (electrons and ions) produced in SNRs and how does their presence influence or reveal SNR dynamics? How does SNR blast energy partition into various components over time and what controls conversion between components? In lieu of a proceedings volume, we present here a synopsis of the workshop in the form of brief summaries of the workshop sessions. The sharpest impressions from the workshop were the crucial and under-appreciated roles that environments have on SNR appearance and dynamics and the critical need for broad-based studies to understand these beautiful, but enigmatic objects. \\

We have carried out high resolution MHD simulations of the nonlinear evolution of Kelvin-Helmholtz unstable flows in 2 1/2 dimensions. The modeled flows and fields were initially uniform except for a thin shear layer with a hyperbolic tangent velocity profile and a small, normal mode perturbation. The calculations consider periodic sections of flows containing magnetic fields parallel to the shear layer, but projecting over a full range of angles with respect to the flow vectors. They are intended as preparation for fully 3D calculations and to address two specific questions raised in earlier work: 1) What role, if any, does the orientation of the field play in nonlinear evolution of the MHD Kelvin-Helmholtz instability in 2 1/2 D. 2) Given that the field is too weak to stabilize against a linear perturbation of the flow, how does the nonlinear evolution of the instability depend on strength of the field. The magnetic field component in the third direction contributes only through minor pressure contributions, so the flows are essentially 2D. Even a very weak field can significantly enhance the rate of energy dissipation. In all of the cases we studied magnetic field amplification by stretching in the vortex is limited by tearing mode, ``fast'' reconnection events that isolate and then destroy magnetic flux islands within the vortex and relax the fields outside the vortex. If the magnetic tension developed prior to reconnection is comparable to Reynolds stresses in the flow, that flow is reorganized during reconnection. Otherwise, the primary influence on the plasma is generation of entropy. The effective expulsion of flux from the vortex is very similar to that shown by Weiss for passive fields in idealized vortices with large magnetic Reynolds numbers. We demonstrated that this

We examine the formation and evolution of the density enhancement (density spike) that appears downstream of strong, cosmic-ray-modified shocks. This feature results from temporary overcompression of the flow by the combined cosmic-ray shock precursor/gas subshock. Formation of the density spike is expected whenever shock modification by cosmic-ray pressure increases strongly. That occurence may be anticipated for newly generated strong shocks or for cosmic-ray-modified shocks encountering a region of higher external density, for example. The predicted mass density within the spike increases with the shock Mach number and with shocks more dominated by cosmic-ray pressure. We find this spike to be linearly unstable under a modified Rayleigh-Taylor instability criterion at the early stage of its formation. We confirm this instability numerically using two independent codes based on the two-fluid model for cosmic-ray transport. These two-dimensional simulations show that the instability grows impulsively at early stages and then slows down as the gradients of total pressure and gas density decrease. Observational discovery of this unstable density spike behind shocks, possibly through radio emission enhanced by the amplified magnetic fields would provide evidence for the existence of strongly cosmic-ray modified shock structures.

We report simulations of diffusive particle acceleration in oblique magnetohydrodynamical (MHD) shocks. These calculations are based on extension to oblique shocks of a numerical model for ``thermal leakage'' injection of particles at low energy into the cosmic-ray population. That technique, incorporated into a fully dynamical diffusion-convection formalism, was recently introduced for parallel shocks by Kang \& Jones (1995). Here, we have compared results of time dependent numerical simulations using our technique with Monte Carlo simulations by Ellison, Baring \& Jones 1995 and with {\it in situ} observations from the Ulysses spacecraft of oblique interplanetary shocks discussed by Baring \etal (1995). Through the success of these comparisons we have demonstrated that our {diffusion-convection} method and injection techniques provide a practical tool to capture essential physics of the injection process and particle acceleration at oblique MHD shocks. In addition to the diffusion-convection simulations, we have included time dependent two-fluid simulations for a couple of the shocks to demonstrate the basic validity of that formalism in the oblique shock context. Using simple models for the two-fluid closure parameters based on test-particle considerations, we find good agreement with the dynamical properties of the more detailed diffusion-convection results. We emphasize, however, that such two-fluid results can be sensitive to the properties of these closure parameters when the flows are not truly steady. Furthermore, we emphasize through example how the validity of the two-fluid formalism does not necessarily mean that {\it steady-state} two-fluid models provide a reliable tool for predicting the efficiency of particle acceleration in real shocks.

We present 3-D numerical hydrodynamical simulations of precessing supersonic heavy jets to explore how well they serve as a model for generating molecular outflows from Young Stellar Objects. The dynamics are studied with a number of high resolution simulations on a Cartesian grid (128x128x128 zones) using a high order finite difference method. A range of cone angles and precession rates were included in the study. Two higher resolution runs (256x256x256 zones) were made for comparison in order to confirm numerical convergence of global flow characteristics. Morphological, kinematical and dynamical characteristics of precessing jets are described and compared to important properties of straight jets and also to observations of YSOs. In order to examine the robustness of precessing jets as a mean to produce molecular outflows around Young Stellar Objects, ``synthetic observations'' of the momentum distributions of the simulated precessing jets are compared to observations of molecular outflows. It is found that precessing jets match better the morphology, highly forward driven momentum and momentum distributions along the long axis of molecular outflows than do wind-driven or straight jet-driven flow models.

We report an extensive set of two-dimensional MHD simulations exploring the role and evolution of magnetic fields in the dynamics of supersonic plasma clumps. We examine the influence of both ambient field strength and orientation on the problem. Of those two characteristics, field orientation is far more important in the cases we have considered with $\beta_0 = p_g/p_b \ge 1$. That is due to the geometry-sensitivity of field stretching/amplification from large-scale shearing motions around the bullet. When the ambient magnetic field is transverse to the bullet motion, even a very modest field, well below equipartition strength, can be amplified by field line stretching around the bullet within a couple of bullet crushing times so that Maxwell stresses become comparable to the ram pressure associated with the bullet motion. The possibility is discussed that those situations might lead to large, induced electric potentials capable of accelerating charged particles. When the ambient field is aligned to the bullet motion, on the other hand, reconnection- prone topologies develop that shorten the stretched field and release much of the excess energy it contains. In this geometry, the Maxwell stresses on the bullet never approach the ram pressure level. In both cases, however, the presence of a field with even moderate initial strength acts to help the flow realign itself around the bullet into a smoother, more laminar form. That reduces bullet fragmentation tendencies caused by destructive instabilities. Eddies seem less effective at field amplification than flows around the bullet, because fields within eddies tend to be expelled to the eddy perimeters. Similar effects cause the magnetic field within the bullet itself to be reduced below its initial value over time.

We discover a new dynamical mechanism that significantly enhances the growth of Rayleigh-Taylor fingers developed near the contact interface between the supernova ejecta and swept-up ambient gas in young supernova remnants if the supernova remnant expands into a clumpy (cloudy) circumstellar medium. Our numerical simulation demonstrates that large Rayleigh-Taylor fingers can obtain a sufficient terminal velocity to protrude through the forward shock front by taking extra kinetic energy from vorticies generated by shock-cloud interactions. We suggest this mechanism as a means to generate the aspherical expansion of the supernova ejecta. Ambient magnetic fields are stretched and amplified as the Rayleigh-Taylor fingers protrude, possibly leading to strongly enhanced radio emission. The material in the protrusions originates from the ejected stellar material with greatly enhanced heavy elements. Therefore, it can be a strong X-ray emitter. The timescale for the Rayleigh-Taylor fingers to reach the forward shock depends on the size, mass density and distribution of clouds being engulfed by the supernova shock, although the details will require further numerical investigation.

Using a new numerical code we have carried out two-dimensional simulations of the nonlinear evolution of unstable sheared magnetohydrodynamic flows. We considered two cases: a strong magnetic field (Alfven Mach number, M_a = 2.5) and a weak field (M_a =5). Each flow rapidly evolves until it reaches a nearly steady condition, which is fundamentally different from the analogous gasdynamic state. Both MHD flows relax to a stable, laminar flow on timescales less than or of the order of 15 linear growth times, measured from saturation of the instability. That timescale is several orders of magnitude less than the nominal dissipation time for these simulated flows, so this condition represents an quasi-steady relaxed state. The strong magnetic field case reaches saturation as magnetic tension in the displaced flow boundary becomes sufficient to stabilize it. That flow then relaxes in a straightforward way to the steady, laminar flow condition. The weak magnetic field case, on the other hand, begins development of the vortex expected for gasdynamics, but that vortex is destroyed by magnetic stresses that locally become strong. Magnetic topologies lead to reconnection and dynamical alignment between magnetic and velocity fields. Together these processes produce a sequence of intermittent vortices and subsequent relaxation to a nearly laminar flow condition in which the magnetic cross helicity is nearly maximized. Remaining irregularities consist of a pair of flux tubes straddling the shear layer. Fluctuations within those features are closely aligned, representing Alfv\'en waves propagating locally downstream.

Three-dimensional hydrodynamic simulations of large scale structure in the Universe have shown that accretion shocks form during the gravitational collapse of one-dimensional caustics, and that clusters of galaxies formed at intersections of the caustics are surrounded by these accretion shocks. Estimated speed and curvature radius of the shocks are 1000-3000 \kms and about 5 Mpc, respectively, in the $\Omega=1$ CDM universe. Assuming that energetic protons are accelerated by these accretion shocks via the first-order Fermi process and modeling particle transport around the shocks through Bohm diffusion, we suggest that protons can be accelerated up to the {\it Greisen cutoff energy} near $6\times 10^{19}$ eV, provided the mean magnetic field strength in the region around the shocks is at least of order a microgauss. We have also estimated the proton flux at earth from the Virgo cluster. Assuming a few (1-10) \% of the ram pressure of the infalling matter would be transferred to the cosmic-rays, the estimated flux for $ E \sim 10^{19}$eV is consistent with observations, so that such clusters could be plausible sources of the UHE CRs.

We present for astrophysical use a multi-dimensional numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference method on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. The constraint of a divergence-free field is enforced exactly by calculating a correction via a gauge transformation in each time step. Results from two-dimensional shock tube tests show that the code captures correctly discontinuities in all three MHD waves families as well as contact discontinuities. The numerical viscosities and resistivity in the code, which are useful in order to understand simulations involving turbulent flows, are estimated through the decay of two-dimensional linear waves. Finally, the robustness of the code in two-dimensions is demonstrated through calculations of the Kelvin-Helmholtz instability and the Orszag-Tang vortex.

Direct comparisons of diffusive particle acceleration numerical simulations have been made against Monte Carlo and hybrid plasma simulations by Ellison {\it et. al.} (1993) and against observations at the earth's bow shock presented by Ellison {\it et. al.} (1990). Toward this end we have introduced a new numerical scheme for injection of cosmic-ray particles out of the thermal plasma, modeled by way of the diffusive scattering process itself; that is, the diffusion and acceleration across the shock front of particles out of the suprathermal tail of the Maxwellian distribution. We find that all of these computational methods (diffusion-advection, two-fluid, Monte Carlo and hybrid) are in substantial agreement on the issues they can simultaneously address, so that the essential physics of diffusive particle acceleration is adequately contained within each. This is despite the fact that each makes what appear to be very different assumptions or approximations.

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic ray (CR) modified plane shock evolution using the two-fluid model for diffusive particle acceleration. The numerical shocks evolve to published analytical steady state properties. In order to probe the dynamical role of magnetic fields we have explored for these time asymptotic states the parameter space of upstream fast mode Mach number, $M_f$, and plasma $\beta$, compiling the results into maps of dynamical steady state CR acceleration efficiency, $\epsilon_c$. These maps, along with additional numerical experiments, show that $\epsilon_c$ is reduced through the action of compressive work on tangential magnetic fields in CR-MHD shocks. Thus $\epsilon_c$ in low $\beta$, moderate $M_f$ shocks tends to be smaller in quasi perpendicular shocks than it would be high $\beta$ shocks of the same $M_f$. This result supports earlier conclusions that strong, oblique magnetic fields inhibit diffusive shock acceleration. For quasi parallel shocks with $\beta < 1$, on the other hand, $\epsilon_c$ seems to be increased at a given $M_f$ when compared to high $\beta$ shocks. The apparent contradiction to the first conclusion results, however, from the fact that for small $\beta$ quasi parallel shocks, the fast mode Mach number is not a good measure of compression through the shock. That is better reflected in the sonic Mach number, which is greater. Acceleration efficiencies for high and low $\beta$ having comparable sonic Mach numbers are more similar. Time evolution of CR-MHD shocks is qualitatively similar to CR-gasdynamical shocks. However, several potentially interesting differences are apparent.

We describe a numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference scheme on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. We also describe a nonlinear Riemann solver for ideal MHD, which includes rarefactions as well as shocks and produces exact solutions for two-dimensional magnetic field structures as well as for the three-dimensional ones. The numerical code and the Riemann solver have been used to test each other. Extensive tests encompassing all the possible ideal MHD structures with planar symmetries (\ie ~one-dimensional flows) are presented. These include those for which the field structure is two-dimensional ({\it i.e.}, those flows often called ``$1 + 1/2$ dimensional'') as well as those for which the magnetic field plane rotates ({\it i.e.,}, those flows often called ``$1 + 1/2 + 1/2$ dimensional''). Results indicate that the code can resolve strong fast, slow, and magnetosonic shocks within 2-4 cells while more cells are required if shocks become weak. With proper stiffening, rotational discontinuities are resolved within 3-5 cells. Contact discontinuities are also resolved within 3-5 cells with stiffening and 6-8 cells without stiffening, while the stiffening on contact discontinuities in some cases generates numerical oscillations. Tangential discontinuities spread over more than 10 cells. Our tests confirm that slow compound structures with two-dimensional magnetic field are composed of intermediate shocks (so called ``2-4''