We develop an analytic model for the power spectra of polarized filamentary structures as a way to study the Galactic polarization foreground to the Cosmic Microwave Background. Our approach is akin to the cosmological halo-model framework, and reproduces the main features of the Planck 353 GHz power spectra. We model the foreground as randomly-oriented, three-dimensional, spheroidal filaments, accounting for their projection onto the sky. The main tunable parameters are the distribution of filament sizes, the filament physical aspect ratio, and the dispersion of the filament axis around the local magnetic field direction. The abundance and properties of filaments as a function of size determine the slopes of the foreground power spectra, as we show via scaling arguments. The filament aspect ratio determines the ratio of $B$-mode power to $E$-mode power, and specifically reproduces the Planck-observed dust ratio of one-half when the short axis is roughly one-fourth the length of the long axis. Filament misalignment to the local magnetic field determines the $TE$ cross-correlation, and to reproduce Planck measurements, we need a (three-dimensional) misalignment angle with a root mean squared dispersion of about 50 degrees. These parameters are not sensitive to the particular filament density profile. By artificially skewing the distribution of the misalignment angle, this model can reproduce the Planck-observed (and parity-violating) $TB$ correlation. The skewing of the misalignment angle necessary to explain $TB$ will cause a yet-unobserved, positive $EB$ dust correlation, a possible target for future experiments.

Traditional cosmological hydrodynamics simulations fail to spatially resolve the circumgalatic medium (CGM), the reservoir of tenuous gas surrounding a galaxy and extending to its virial radius. We introduce the technique of Enhanced Halo Resolution (EHR), enabling more realistic physical modeling of the simulated CGM by consistently forcing gas refinement to smaller scales throughout the virial halo of a simulated galaxy. We investigate the effects of EHR in the Tempest simulations, a suite of Enzo-based cosmological zoom simulations following the evolution of an L* galaxy, resolving spatial scales of 500 comoving pc out to 100 comoving kpc in galactocentric radius. Among its many effects, EHR (1) changes the thermal balance of the CGM, increasing its cool gas content and decreasing its warm/hot gas content; (2) preserves cool gas structures for longer periods; and (3) enables these cool clouds to exist at progressively smaller size scales. Observationally, this results in a boost in "low ions" like H I and a drop in "high ions" like O VI throughout the CGM. These effects of EHR do not converge in the Tempest simulations, but extrapolating these trends suggests that the CGM in reality is a mist consisting of ubiquitous, small, long-lived, cool clouds suspended in a hot medium at the virial temperature of the halo. Additionally, we explore the physical mechanisms to explain why EHR produces the above effects, proposing that it works both by (1) better sampling the distribution of CGM phases enabling runaway cooling in the denser, cooler tail of the phase distribution; and (2) preventing cool gas clouds from artificially mixing with the ambient hot halo and evaporating. Evidence is found for both EHR mechanisms occurring in the Tempest simulations.

We show how new and upcoming advances in the age of time-domain and multi-wavelength astronomy will open up a new venue to probe the diversity of SN~Ia. We discuss this in the context of the ELT (ESO), as well as space based instrument such as James Webb Space Telescope (JWST). As examples we demonstrate how the power of very early observations, within hours to days after the explosion, and very late-time observations, such as light curves and mid-infrared spectra beyond 3 years, can be used to probe the link to progenitors and explosion scenarios. We identify the electron-capture cross sections of Cr, Mn, and Ni/Co as one of the limiting factors we will face in the future.

As the closest Type Ia supernova in decades, SN 2014J provides a unique opportunity for detailed investigation into observational signatures of the progenitor system and explosion mechanism in addition to burning product distribution. We present a late-time near-infrared spectral series from Gemini-N at $307-466$ days after the explosion. Following the $H$-band evolution probes the distribution of radioactive iron group elements, the extent of mixing, and presence of magnetic fields in the expanding ejecta. Comparing the isolated $1.6440$ $\mu$m [Fe II] emission line with synthetic models shows consistency with a Chandrasekhar-mass white dwarf of $\rho_c=0.7\times10^9$ g cm${}^{-3}$ undergoing a delayed detonation. The ratio of the flux in the neighboring $1.54$ $\mu$m emission feature to the flux in the $1.6440$ $\mu$m feature shows evidence of some limited mixing of stable and radioactive iron group elements in the central regions. Additionally, the evolution of the $1.6440$ $\mu$m line shows an intriguing asymmetry. When measuring line-width of this feature, the data show an increase in line width not seen in the evolution of the synthetic spectra, corresponding to $\approx1{,}000$ km s${}^{-1}$, which could be caused by a localized transition to detonation producing asymmetric ionization in the ejecta. Using the difference in width between the different epochs, an asymmetric component in the central regions, corresponding to approximately the inner $2\times10^{-4}$ of white dwarf mass suggests an off-center ignition of the initial explosion and hence of the kinematic center from the chemical center. Several additional models investigated, including a He detonation and a merger, have difficulty reproducing the features seen these spectra.

This article presents the study of the effects of magnetic fields on non-distributed nuclear burning fronts as a possible solution to a fundamental problem for the thermonuclear explosion of a Chandrasekhar mass ($M_{Ch}$) white dwarf (WD), the currently favored scenario for the majority of Type Ia SNe (SNe~Ia). All existing 3D hydrodynamical simulations predict strong global mixing of the burning products due to Rayleigh-Taylor (RT) instabilities, which is in contradiction with observations. As a first step and to study the flame physics we present a set of computational magneto-hydrodynamic (MHD) models in rectangular flux tubes, resembling a small inner region of a WD. We consider initial magnetic fields up to $10^{12}\,\,\mathrm{G}$ of various orientations. We find an increasing suppression of RT instabilities starting at about $10^9\,\,\mathrm{G}$. The front speed tends to decrease with increasing magnitude up to about $10^{11}\,\,\mathrm{G}$. For even higher fields new small scale finger-like structures develop, which increase the burning speed by a factor of 3 to 4 above the field-free RT-dominated regime. We suggest that the new instability may provide sufficiently accelerated energy production during the distributed burning regime to go over the Chapman-Jougey limit and trigger a detonation. Finally we discuss the possible origins of high magnetic fields during the final stage of the progenitor evolution or the explosion.

The column density probability distribution function (N-PDF) of GMCs has been used as a diagnostic of star formation. Simulations and analytic predictions have suggested the N-PDF is composed of a low density lognormal component and a high density power-law component, tracing turbulence and gravitational collapse, respectively. In this paper, we study how various properties of the true 2D column density distribution create the shape, or "anatomy" of the PDF. We test our ideas and analytic approaches using both a real, observed, PDF based on Herschel observations of dust emission as well as a simulation that uses the ENZO code. Using a dendrogram analysis, we examine the three main components of the N-PDF: the lognormal component, the power-law component, and the transition point between these two components. We find that the power-law component of an N-PDF is the summation of N-PDFs of power-law substructures identified by the dendrogram algorithm. We also find that the analytic solution to the transition point between lognormal and power-law components proposed by Burkhart, Stalpes & Collins (2017) is applicable when tested on observations and simulations, within the uncertainties. We reconfirm and extend the results of Lombardi, Alves & Lada (2015), which stated that the lognormal component of the N-PDF is difficult to constrain due to the artificial choice of the map area. Based on the resulting anatomy of the N-PDF, we suggest avoiding analyzing the column density structures of a star forming region based solely on fits to the lognormal component of an N-PDF. We also suggest applying the N-PDF analysis in combination with the dendrogram algorithm, to obtain a more complete picture of the global and local environments and their effects on the density structures.

We study giant molecular cloud (GMC) collisions and their ability to trigger star cluster formation. We further develop our three dimensional magnetized, turbulent, colliding GMC simulations by implementing star formation sub-grid models. Two such models are explored: (1) "Density-Regulated," i.e., fixed efficiency per free-fall time above a set density threshold; (2) "Magnetically-Regulated," i.e., fixed efficiency per free-fall time in regions that are magnetically supercritical. Variations of parameters associated with these models are also explored. In the non-colliding simulations, the overall level of star formation is sensitive to model parameter choices that relate to effective density thresholds. In the GMC collision simulations, the final star formation rates and efficiencies are relatively independent of these parameters. Between non-colliding and colliding cases, we compare the morphologies of the resulting star clusters, properties of star-forming gas, time evolution of the star formation rate (SFR), spatial clustering of the stars, and resulting kinematics of the stars in comparison to the natal gas. We find that typical collisions, by creating larger amounts of dense gas, trigger earlier and enhanced star formation, resulting in 10 times higher SFRs and efficiencies. The star clusters formed from GMC collisions show greater spatial sub-structure and more disturbed kinematics.

We present an analysis of photometric observations of the eclipsing novalike variable DW UMa made by the CBA consortium between 1999 and 2015. Analysis of 372 new and 260 previously published eclipse timings reveals a 13.6 year period or quasi-period in the times of minimum light. The seasonal light curves show a complex spectrum of periodic signals: both positive and negative superhumps, likely arising from a prograde apsidal precession and a retrograde nodal precession of the accretion disc. These signals appear most prominently and famously as sidebands of the orbital frequency but the precession frequencies themselves, at 0.40 and 0.22 cycles per day, are also seen directly in the power spectrum. The superhumps are sometimes seen together and sometimes separately. The depth, width and skew of eclipses are all modulated in phase with both nodal and apsidal precession of the tilted and eccentric accretion disc. The superhumps, or more correctly the precessional motions which produce them, may be essential to understanding the mysterious SW Sextantis syndrome. Disc wobble and eccentricity can both produce Doppler signatures inconsistent with the true dynamical motions in the binary, and disc wobble might boost the mass-transfer rate by enabling the hot white dwarf to directly irradiate the secondary star.

We investigate giant molecular cloud (GMCs) collisions and their ability to induce gravitational instability and thus star formation. This mechanism may be a major driver of star formation activity in galactic disks. We carry out a series of three dimensional, magnetohydrodynamics (MHD), adaptive mesh refinement (AMR) simulations to study how cloud collisions trigger formation of dense filaments and clumps. Heating and cooling functions are implemented based on photo-dissociation region (PDR) models that span the atomic to molecular transition and can return detailed diagnostic information. The clouds are initialized with supersonic turbulence and a range of magnetic field strengths and orientations. Collisions at various velocities and impact parameters are investigated. Comparing and contrasting colliding and non-colliding cases, we characterize morphologies of dense gas, magnetic field structure, cloud kinematic signatures, and cloud dynamics. We present key observational diagnostics of cloud collisions, especially: relative orientations between magnetic fields and density structures, like filaments; $^{13}$CO($J$=2-1), $^{13}$CO($J$=3-2), and $^{12}$CO($J$=8-7) integrated intensity maps and spectra; and cloud virial parameters. We compare these results to observed Galactic clouds.

We derive an analytic expression for the transitional column density value ($s_t$) between the lognormal and power-law form of the probability distribution function (PDF) in star-forming molecular clouds. Our expression for $s_t$ depends on the mean column density, the variance of the lognormal portion of the PDF, and the slope of the power-law portion of the PDF. We show that $s_t$ can be related to physical quantities such as the sonic Mach number of the flow and the power-law index for a self-gravitating isothermal sphere. This implies that the transition point between the lognormal and power-law density/column density PDF represents the critical density where turbulent and thermal pressure balance, the so-called "post-shock density." We test our analytic prediction for the transition column density using dust PDF observations reported in the literature as well as numerical MHD simulations of self-gravitating supersonic turbulence with the Enzo code. We find excellent agreement between the analytic $s_t$ and the measured values from the numerical simulations and observations (to within 1.5 A$_V$). We discuss the utility of our expression for determining the properties of the PDF from unresolved low density material in dust observations, for estimating the post-shock density, and for determining the HI-H$_2$ transition in clouds.

Additional physics beyond standard hydrodynamics is needed to fully model the intracluster medium (ICM); however, as we move to more sophisticated models, it is important to consider the role of magnetic fields and the way the fluid approximation breaks down. This paper represents a first step towards developing a self-consistent model of the ICM by characterizing the statistical properties of magnetic fields in cosmological simulations of galaxy clusters. We find that plasma conditions are largely homogeneous across a range of cluster masses and relaxation states. We also find that the magnetic field length scales are resolution dependent and not based on any particular physical process. Energy transfer mechanisms and scales are also identified, and imply the existence of small scale dynamo action. The scales of the small scale dynamo are resolution limited and driven by numerical resistivity and viscosity.

The nearby cloud L1642 is one of only two known very high latitude (|b| > 30 deg) clouds actively forming stars. It is a rare example of star formation in isolated conditions, and can reveal important details of star formation in general, e.g., of the effect of magnetic fields. We compare Herschel dust emission structures and magnetic field orientation revealed by Planck polarization maps in L1642. The high-resolution ($\sim20"$) Herschel data reveal a complex structure including a dense, compressed central clump, and low density striations. The Planck polarization data (at 10$'$ resolution) reveal an ordered magnetic field pervading the cloud and aligned with the surrounding striations. There is a complex interplay between the cloud structure and large scale magnetic field. This suggests that the magnetic field is closely linked to the formation and evolution of the cloud. CO rotational emission confirms that the striations are connected with the main clumps and likely to contain material either falling into or flowing out of the clumps. There is a clear transition from aligned to perpendicular structures approximately at a column density of $N_{\rm{H}} = 1.6 \times 10^{21}\, {{\rm cm}}^{-2}$. Comparing the Herschel maps with the Planck polarization maps shows the close connection between the magnetic field and cloud structure even in the finest details of the cloud.

We report the results of a long campaign of time-series photometry on the nova-like variable UX Ursae Majoris during 2015. It spanned 150 nights, with ~1800 hours of coverage on 121 separate nights. The star was in its normal `high state' near magnitude V=13, with slow waves in the light curve and eclipses every 4.72 hours. Remarkably, the star also showed a nearly sinusoidal signal with a full amplitude of 0.44 mag and a period of 3.680 +/- 0.007 d. We interpret this as the signature of a retrograde precession (wobble) of the accretion disc. The same period is manifest as a +/-33 s wobble in the timings of mid-eclipse, indicating that the disc's centre of light moves with this period. The star also showed strong `negative superhumps' at frequencies w_orb+N and 2w_orb+N, where w_orb and N are respectively the orbital and precession frequencies. It is possible that these powerful signals have been present, unsuspected, throughout the more than 60 years of previous photometric studies.

Classical Sweet-Parker models of reconnection predict that reconnection rates depend inversely on the resistivity, usually parameterized using the dimensionless Lundquist number ($\Lund$). We describe magnetohydrodynamic (MHD) simulations using a static, nested grid that show the development of a three-dimensional instability in the plane of a current sheet between reversing field lines without a guide field. The instability leads to rapid reconnection of magnetic field lines at a rate independent of $\Lund$ over at least the range $3.2\times 10^3 \lesssim \Lund \lesssim 3.2 \times 10^5$ resolved by the simulations. We find that this instability occurs even for cases with $\Lund \lesssim 10^4$ that in our models appear stable to the recently described, two-dimensional, plasmoid instability. Our results suggest that three-dimensional, MHD processes alone produce fast (resistivity independent) reconnection without recourse to kinetic effects or external turbulence. The unstable reconnection layers provide a self-consistent environment in which the extensively studied turbulent reconnection process can occur.

We study the observable signatures of self-gravitating MHD turbulence by applying the probability density functions (PDFs) and the spatial density power spectrum to synthetic column density maps. We find that there exists three characterizable stages of the evolution of the collapsing cloud which we term "early," "intermediate," and "advanced." At early times, i.e. $t<0.15t_{ff}$, the column density has a power spectral slope similar to nongravitating supersonic turbulence and a lognormal distribution. At an intermediate stage, i.e. $0.15t_{ff}< t \leq 0.35t_{ff}$, there exists signatures of the prestellar cores in the shallower PDF and power spectrum power law slopes. The column density PDF power law tails at these times have line of sight averaged slopes ranging from -2.5 to -1.5 with shallower values belonging to simulations with lower magnetic field strength. The density power spectrum slope becomes shallow and can be characterized by $P(k)=A_1k^{\beta_2}e^{-k/k_c}$, where $A_1$ describes the amplitude, $k^{\beta_2}$ describes the classical power law behavior and the scale $k_c$ characterizes the turn over from turbulence dominated to self-gravity dominated. At advanced stages of collapse, i.e. $\approx t>0.35t_{ff}$, the power spectral slope is positive valued, and a dramatic increase is observed in the PDF moments and the Tsallis incremental PDF parameters, which gives rise to deviations between PDF-sonic Mach number relations. Finally, we show that the imprint of gravity on the density power spectrum can be replicated in non-gravitating turbulence by introducing a delta-function with amplitude equivalent to the maximum valued point in a given self-gravitating map. We find that the turbulence power spectrum restored through spatial filtering of the high density material.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

Comparisons of the integrated thermal pressure support of gas against its gravitational potential energy lead to critical mass scales for gravitational instability such as the Jeans and the Bonnor-Ebert masses, which play an important role in analysis of many physical systems, including the heuristics of numerical simulations. In a strict theoretical sense, however, neither the Jeans nor the Bonnor-Ebert mass are meaningful when applied locally to substructure in a self-gravitating turbulent medium. For this reason, we investigate the local support by thermal pressure, turbulence, and magnetic fields against gravitational compression through an approach that is independent of these concepts. At the centre of our approach is the dynamical equation for the divergence of the velocity field. We carry out a statistical analysis of the source terms of the local compression rate (the negative time derivative of the divergence) for simulations of forced self-gravitating turbulence in periodic boxes with zero, weak, and moderately strong mean magnetic fields (measured by the averages of the magnetic and thermal pressures). We also consider the amplification of the magnetic field energy by shear and by compression. Thereby, we are able to demonstrate that the support against gravity is dominated by thermal pressure fluctuations, although magnetic pressure also yields a significant contribution. The net effect of turbulence in the highly supersonic regime, however, is to enhance compression rather than supporting overdense gas even if the vorticity is very high. This is incommensurate with the support of the highly dynamical substructures in magneto-turbulent fluids being determined by local virial equilibria of volume energies without surface stresses.

We present the discovery and timing solutions of five new pulsars by students involved in the Pulsar Search Collaboratory (PSC), a NSF-funded joint program between the National Radio Astronomy Observatory and West Virginia University designed to excite and engage high-school students in Science, Technology, Engineering, and Mathematics (STEM) and related fields. We encourage students to pursue STEM fields by apprenticing them within a professional scientific community doing cutting edge research, specifically by teaching them to search for pulsars. The students are analyzing 300 hours of drift-scan survey data taken with the Green Bank Telescope at 350 MHz. These data cover 2876 square degrees of the sky. Over the course of five years, more than 700 students have inspected diagnostic plots through a web-based graphical interface designed for this project. The five pulsars discovered in the data have spin periods ranging from 3.1 ms to 4.8 s. Among the new discoveries are - PSR J1926-1314, a long period, nulling pulsar; PSR J1821+0155, an isolated, partially recycled 33-ms pulsar; and PSR J1400-1438, a millisecond pulsar in a 9.5-day orbit whose companion is likely a white dwarf star.

Non-thermal radio emission from cosmic ray electrons in the vicinity of merging galaxy clusters is an important tracer of cluster merger activity, and is the result of complex physical processes that involve magnetic fields, particle acceleration, gas dynamics, and radiation. In particular, objects known as radio relics are thought to be the result of shock-accelerated electrons that, when embedded in a magnetic field, emit synchrotron radiation in the radio wavelengths. In order to properly model this emission, we utilize the adaptive mesh refinement simulation of the magnetohydrodynamic evolution of a galaxy cluster from cosmological initial conditions. We locate shock fronts and apply models of cosmic ray electron acceleration that are then input into radio emission models. We have determined the thermodynamic properties of this radio-emitting plasma and constructed synthetic radio observations to compare to observed galaxy clusters. We find a significant dependence of the observed morphology and radio relic properties on the viewing angle of the cluster, raising concerns regarding the interpretation of observed radio features in clusters. We also find that a given shock should not be characterized by a single Mach number. We find that the bulk of the radio emission comes from gas with T>5x10^7, \rho~10^(-28)-10^(-27) g/cm^3, with magnetic field strengths of 0.1-1.0 \mu G and shock Mach numbers of M~3-6. We present an analysis of the radio spectral index which suggests that the spatial variation of the spectral index can mimic synchrotron aging. Finally, we examine the polarization fraction and position angle of the simulated radio features, and compare to observations.

Radio observations of galaxy clusters show that there are $\mu$G magnetic fields permeating the intra-cluster medium (ICM), but it is hard to accurately constrain the strength and structure of the magnetic fields without the help of advanced computer simulations. We present qualitative comparisons of synthetic VLA observations of simulated galaxy clusters to radio observations of Faraday Rotation Measure (RM) and radio halos. The cluster formation is modeled using adaptive mesh refinement (AMR) magneto-hydrodynamic (MHD) simulations with the assumption that the initial magnetic fields are injected into the ICM by active galactic nuclei (AGNs) at high redshift. In addition to simulated clusters in Xu et al. (2010, 2011), we present a new simulation with magnetic field injections from multiple AGNs. We find that the cluster with multiple injection sources is magnetized to a similar level as in previous simulations with a single AGN. The RM profiles from simulated clusters, both $|RM|$ and the dispersion of RM ($\sigma_{RM}$), are consistent at a first-order with the radial distribution from observations. The correlations between the $\sigma_{RM}$ and X-ray surface brightness from simulations are in a broad agreement with the observations, although there is an indication that the simulated clusters could be slightly over-dense and less magnetized with respect to those in the observed sample. In addition, the simulated radio halos agree with the observed correlations between the radio power versus the cluster X-ray luminosity and between the radio power versus the radio halo size. These studies show that the cluster wide magnetic fields that originate from AGNs and are then amplified by the ICM turbulence (Xu et al. 2010) match observations of magnetic fields in galaxy clusters.

We examine the effects of self-gravity and magnetic fields on supersonic turbulence in isothermal molecular clouds with high resolution simulations and adaptive mesh refinement. These simulations use large root grids (512^3) to capture turbulence and four levels of refinement to capture high density, for an effective resolution of 8,196^3. Three Mach 9 simulations are performed, two super-Alfv\'enic and one trans-Alfv\'enic. We find that gravity splits the clouds into two populations, one low density turbulent state and one high density collapsing state. The low density state exhibits properties similar to non-self-gravitating in this regime, and we examine the effects of varied magnetic field strength on statistical properties: the density probability distribution function is approximately lognormal; velocity power spectral slopes decrease with field strength; alignment between velocity and magnetic field increases with field; the magnetic field probability distribution can be fit to a stretched exponential. The high density state is characterized by self-similar spheres; the density PDF is a power-law; collapse rate decreases with increasing mean field; density power spectra have positive slopes, P({\rho},k) \propto k; thermal-to-magnetic pressure ratios are unity for all simulations; dynamic-to-magnetic pressure ratios are larger than unity for all simulations; magnetic field distribution is a power-law. The high Alfv\'en Mach numbers in collapsing regions explain recent observations of magnetic influence decreasing with density. We also find that the high density state is found in filaments formed by converging flows, consistent with recent Herschel observations. Possible modifications to existing star formation theories are explored.

Theory and simulations suggest that magnetic fields from radio jets and lobes powered by their central super massive black holes can be an important source of magnetic fields in the galaxy clusters. This is paper II in a series of studies where we present self-consistent high-resolution adaptive mesh refinement cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an active galactic nucleus (AGN). We studied 12 different galaxy clusters with virial masses ranging from 1 $\times$ 10$^{14}$ to 2 $\times$ 10$^{15}$ M$_{\odot}$. In this work we examine the effects of the mass and merger history on the final magnetic properties. We find that the evolution of magnetic fields is qualitatively similar to those of previous studies. In most clusters, the injected magnetic fields can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process with hierarchical mergers, while the amplification history and the magnetic field distribution depend on the cluster formation and magnetism history. This can be very different for different clusters. The total magnetic energies in these clusters are between 4 $\times$ 10$^{57}$ and $10^{61}$ erg, which is mainly decided by the cluster mass, scaling approximately with the square of the total mass. Dynamically older relaxed clusters usually have more magnetic fields in their ICM. The dynamically very young clusters may be magnetized weakly since there is not enough time for magnetic fields to be amplified.

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

We study the reliability of mass estimates obtained for molecular cloud cores using sub-millimetre and infrared dust emission. We use magnetohydrodynamic simulations and radiative transfer to produce synthetic observations with spatial resolution and noise levels typical of Herschel surveys. We estimate dust colour temperatures using different pairs of intensities, calculate column densities and compare the estimated masses with the true values. We compare these results to the case when all five Herschel wavelengths are available. We investigate the effects of spatial variations of dust properties and the influence of embedded heating sources. Wrong assumptions of dust opacity and its spectral index beta can cause significant systematic errors in mass estimates. These are mainly multiplicative and leave the slope of the mass spectrum intact, unless cores with very high optical depth are included. Temperature variations bias colour temperature estimates and, in quiescent cores with optical depths higher than for normal stable cores, masses can be underestimated by up to one order of magnitude. When heated by internal radiation sources the observations recover the true mass spectra. The shape, although not the position, of the mass spectrum is reliable against observational errors and biases introduced in the analysis. This changes only if the cores have optical depths much higher than expected for basic hydrostatic equilibrium conditions. Observations underestimate the value of beta whenever there are temperature variations along the line of sight. A bias can also be observed when the true beta varies with wavelength. Internal heating sources produce an inverse correlation between colour temperature and beta that may be difficult to separate from any intrinsic beta(T) relation of the dust grains. This suggests caution when interpreting the observed mass spectra and the spectral indices.

We report unfiltered CCD photometry of the first confirmed superoutburst of the recently discovered dwarf nova, SDSS J083931.35+282824.0 in April 2010. From a quiescence magnitude of ~19.8 it rose to 14.0, an outburst amplitude of at least 5.8 magnitudes. Only the plateau phase of the outburst was observed during which superhumps with peak-to-peak amplitude of up to 0.28 magnitudes were present, confirming this to be an SU UMa type dwarf nova. The mean superhump period was Psh = 0.07836(2) during the first 3 days and this subsequently decreased to 0.07800(3) d. Analysis of the data revealed tentative evidence for an orbital period Porb = 0.07531(25) d. The fractional superhump period excess was epsilon = 0.039(6), which is consistent with other dwarf novae of similar orbital period.

We present a series of cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an Active Galactic Nucleus (AGN). Specifically, we investigate the influence of both the epoch of AGN (z $\sim$ 3-0.5) and the AGN energy ($\sim$ 3 $\times$ 10$^{57}$ - 2 $\times$ 10$^{60}$ ergs)on the final magnetic field distribution in a relatively massive cluster (M$_{vir}$ $\sim$10$^{15}$ M$_\odot$). We find that as long as the AGN magnetic fields are ejected before the major mergers in the cluster formation history, magnetic fields can be transported throughout the cluster and can be further amplified by the intra-cluster medium (ICM) turbulence cause by hierarchical mergers during the cluster formation process. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs, with micro Gauss fields distributed over $\sim$ Mpc scale. The amplification of the total magnetic energy by the ICM turbulence can be significant, up to $\sim$1000 times in some cases. Therefore even weak magnetic fields from AGNs can be used to magnetize the cluster to the observed level. The final magnetic energy in the ICM is determined by the ICM turbulent energy, with a weak dependence on the AGN injection energy. We discuss the properties of magnetic fields throughout the cluster and the synthetic Faraday rotation measure maps they produce. We also show that high spatial resolution over most of the magnetic regions of the cluster is very important to capture the small scale dynamo process and maintain the magnetic field structure in our simulations.

In this work, we present the mass and magnetic distributions found in a recent Adaptive Mesh Refinement (AMR) MHD simulation of supersonic, \sa, self gravitating turbulence. Powerlaw tails are found in both volume density and magnetic field probability density functions, with $P(\rho) \propto \rho^{-1.67}$ and $P(B)\propto B^{-2.74}$. A power law is also found between magnetic field strength and density, with $B\propto \rho^{0.48}$, throughout the collapsing gas. The mass distribution of gravitationally bound cores is shown to be in excellent agreement with recent observation of prestellar cores. The mass to flux distribution of cores is also found to be in excellent agreement with recent Zeeman splitting measurements.

Power law mass-size and linewidth-size correlations, two of "Larson's laws," are often studied to assess the dynamical state of clumps within molecular clouds. Using the result of a hydrodynamic simulation of a molecular cloud, we investigate how geometric projection may affect the derived Larson relationships. We find that large scale structures in the column density map have similar masses and sizes to those in the 3D simulation (PPP). Smaller scale clumps in the column density map are measured to be more massive than the PPP clumps, due to the projection of all emitting gas along lines of sight. Further, due to projection effects, structures in a synthetic spectral observation (PPV) may not necessarily correlate with physical structures in the simulation. In considering the turbulent velocities only, the linewidth-size relationship in the PPV cube is appreciably different from that measured from the simulation. Including thermal pressure in the simulated linewidths imposes a minimum linewidth, which results in a better agreement in the slopes of the linewidth-size relationships, though there are still discrepancies in the offsets, as well as considerable scatter. Employing commonly used assumptions in a virial analysis, we find similarities in the computed virial parameters of the structures in the PPV and PPP cubes. However, due to the discrepancies in the linewidth- and mass- size relationships in the PPP and PPV cubes, we caution that applying a virial analysis to observed clouds may be misleading due to geometric projection effects. We speculate that consideration of physical processes beyond kinetic and gravitational pressure would be required for accurately assessing whether complex clouds, such as those with highly filamentary structure, are bound.

We present RXTE observations covering six normal outbursts of the dwarf nova SU UMa, the prototype of its class. The outbursts showed consistent X-ray behaviour with the X-ray count rate dropping suddenly by a factor of our, and with evidence for a half day delay between the optical rise and the X-ray suppression. In contrast to SS Cyg, an X-ray flux increase at the beginning of outburst was not observed, although it is expected from boundary layer models. The X-ray flux was high and decreasing during quiescence, in conflict with the disc instability model. The X-ray spectrum of SU UMa was softer in outburst than during quiescence, and it was consistent with constant reflection.

BZ UMa is a cataclysmic variable star whose specific classification has eluded researchers since its discovery in 1968. It has outburst and spectral properties consistent with both U Gem class dwarf novae and intermediate polars. We present new photometric and polarimetric measurements of recent outbursts, including the first detected superoutburst of the system. Statistical analysis of these and archival data from outbursts over the past 40 years present a case for BZ UMa as a non-magnetic, U Gem class, SU-UMa subclass dwarf novae.

We present self-consistent cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and the magnetic field ejection by an active galactic nucleus (AGN). We find that the magnetic fields ejected by the AGNs, though initially distributed in relatively small volumes, can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process. The ICM turbulence is shown to be generated and sustained by the frequent mergers of smaller halos. Furthermore, a cluster-wide dynamo process is shown to exist in the ICM and amplify the magnetic field energy and flux. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs while micro Gauss ($\mu$G) fields can distribute over $\sim$ Mpc scales throughout the whole cluster. This finding shows that magnetic fields from AGNs, being further amplified by the ICM turbulence through small-scale dynamo processes, can be the origin of cluster-wide magnetic fields.

We present photometry of the first reported superoutburst of the dwarf nova VSX J074727.6+065050 during 2008 January and February. At its brightest the star reached magnitude 11.4 and this was followed by a slow decline at 0.09 mag/d for 19 days, corresponding to the plateau phase. There was then a rapid decline at 1.66 mag/d to a temporary minimum at magnitude 16.6 where it stayed for 2 to 3 days after which there were six remarkable echo outbursts before the star gradually faded back towards quiescence at ~magnitude 19.5. The overall outburst amplitude was at least 8 magnitudes and it lasted more than 80 days. During the plateau phase we observed common superhumps with Psh = 0.06070(6) d, but the period increased to Psh = 0.06151(5) d coinciding with the end of the plateau phase and the onset of the rapid decline. This corresponds to a continuous period change with P^dot = +4.4(9) x 10-5. During the echo outbursts there was a superhump regime with Psh = 0.06088(49) d. Evidence is presented which is consistent with the star being a member of the WZ Sge family of dwarf novae.

In this work, we present MHDEnzo, the extension of the cosmological code Enzo to include the effects magnetic fields through the ideal MHD approximation. We use a higher order Godunov Riemann solver for the computation of interface fluxes. We use two constrained transport methods to compute the electric field from those interface fluxes, which simultaneously advances the induction equation and maintains the divergence of the magnetic field. A third order divergence free reconstruction technique is used to interpolate the magnetic fields in the block structured AMR framework already extant in Enzo. This reconstruction also preserves the divergence of the magnetic field to machine precision. We use operator splitting to include gravity and cosmological expansion. We then present a series of cosmological and non cosmological tests problems to demonstrate the quality of solution resulting from this combination of solvers.

We present cosmological magnetohydrodynamic simulations of the formation of a galaxy cluster with magnetic energy feedback from an active galactic nuclei (AGN). We demonstrate that X-ray cavities can be produced by the magnetically dominated jet-lobe system that is supported by a central axial current. The cavities are magnetically dominated and their morphology is determined jointedly by the magnetic fields and the background cluster pressure profile. The expansion and motion of the cavities are driven initially by the Lorentz force of the magnetic fields, and the cavities only become buoyant at late stages ($> 500$ Myr). We find that up to $80%-90%$ of the injected magnetic energy goes into doing work against the hot cluster medium, heating it, and lifting it in the cluster potential.

We report the results of the first self-consistent three-dimensional adaptive mesh refinement magnetohydrodynamical simulations of Population III star formation including the Biermann Battery effect. We find that the Population III stars formed including this effect are both qualitatively and quantitatively similar to those from hydrodynamics-only (non-MHD) cosmological simulations. We observe peak magnetic fields of ~10^-9 G in the center of our star-forming halo at z ~ 17.55. The magnetic fields created by the Biermann Battery effect are predominantly formed early in the evolution of the primordial halo at low density and large spatial scales, and then grow through compression and by shear flows. The fields seen in this calculation are never large enough to be dynamically important (with beta >= 10^{15} at all times), and should be considered the minimum possible fields in existence during Population III star formation, and may be seed fields for the stellar dynamo or the magnetorotational instability at higher densities and smaller spatial scales.

Magnetic fields play an important role in almost all astrophysical phenomena including star formation. But due to the difficulty in analytic modeling and observation, magnetic fields are still poorly studied and numerical simulation has become a major tool. We have implemented a cosmological magnetohydrodynamics package for Enzo which is an AMR hydrodynamics code designed to simulate structure formation. We use the TVD solver developed by S. Li as the base solver. In addition, we employ the constrained transport(CT) algorithm as described by D. Balsara. For interpolation magnetic fields to fine grids we used a divergence free quadratic reconstruction, also described by Balsara. We present results from several test problems including MHD caustics, MHD pancake and galaxy cluster formation with magnetic fields. We also discuss possible applications of our AMR MHD code to first star research.