The outer crust properties of cold non-accreting neutron stars are studied within the framework of the quark-meson coupling (QMC) model, which includes the effects of modifications of the quark structure inside individual nucleons when they are within a high-density nuclear medium. With a unique set of five well-constrained adjustable parameters, which have a clear physical basis, the QMC model gives predictions for the ground state observables of even-even nuclei which agree with experiment as well as traditional models. Furthermore, it gives improved theoretical values for nuclei thought to play a role in the outer crusts of neutron stars but for which experimental data is not available. Using the latest experimental data tables wherever possible but otherwise the predictions from the QMC model, we construct an equation of state for the outer crust which is then used within stellar model calculations to obtain an equilibrium sequence of crustal layers, each characterized by a particular neutron rich nuclei. Various properties of the layers are calculated for a range of neutron-star masses, and comparisons are made with alternative equations of state from the literature. This leads to the conclusion that the QMC model successfully predicts the outer crust properties and is fully comparable with the more traditional mass models, which all depend on a larger number of parameters.

Understanding intermittent dissipation in the interstellar medium is integral to understanding the chemical abundances and line luminosities observed therein. In the cold neutral medium, continuous heating processes due to cosmic rays and photoelectric emission from dust grains determine a unique equilibrium temperature varying with density. Excursions from this are due to intermittent processes including shocks, viscous heating, and ambipolar diffusion. The high-temperature excursions are thought to explain the enhanced abundance of CH$^+$ observed along diffuse molecular sight-lines. Intermittent high temperatures should also have an impact on H$_2$ line luminosities. We carry out simulations of MHD turbulence in molecular clouds including heating and cooling, and post-process them to study H$_2$ line emission and hot-gas chemistry, particularly the formation of CH$^+$. We explore multiple magnetic field strengths and equations of state. We use a new H$_2$ cooling function for $n_{\rm H} \leq 10^5\,{\rm cm}^{-3}$, $T\leq 5000\,{\rm K}$, and variable H$_2$ fraction. Our models produce H$_2$ emission lines in accord with many observations, although extra excitation mechanisms are required in some clouds. For realistic r.m.s. magnetic field strengths ($\approx 10$ $\mu$G) and velocity dispersions, we reproduce observed CH$^+$ abundances. Comparison of predicted dust polarization with observations by ${\it Planck}$ suggests that the mean field $\gtrsim 5\,\mu$G, so that the turbulence is sub-Alfv\'enic. We recommend future work treating ions and neutrals as separate fluids to more accurately capture the effects of ambipolar diffusion on CH$^+$ abundance.

We present 3D general relativistic magnetohydrodynamic (GRMHD) simulations of the accretion flow surrounding Sagittarius A* that are initialized using larger-scale MHD simulations of the $\sim$ 30 Wolf--Rayet (WR) stellar winds in the Galactic center. The properties of the resulting accretion flow on horizon scales are set not by ad hoc initial conditions but by the observationally constrained properties of the WR winds with limited free parameters. For this initial study we assume a non-spinning black hole. Our simulations naturally produce a $\sim 10^{-8} M_\odot$ yr$^{-1}$ accretion rate, consistent with previous phenomenological estimates. We find that a magnetically arrested flow is formed by the continuous accretion of coherent magnetic field being fed from large radii. Near the event horizon, the magnetic field is so strong that it tilts the gas with respect to the initial angular momentum and concentrates the originally quasi-spherical flow to a narrow disk-like structure. We also present 230 GHz images calculated from our simulations where the inclination angle and physical accretion rate are not free parameters but are determined by the properties of the WR stellar winds. The image morphology is highly time variable. Linear polarization on horizon scales is coherent with weak internal Faraday rotation.

The design and implementation of a new framework for adaptive mesh refinement (AMR) calculations is described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic execution model based on a simple design called a "task list" that improves parallel performance by overlapping communication and computation, simplifies the inclusion of a diverse range of physics, and even enables multiphysics models involving different physics in different regions of the calculation. We describe physics modules implemented in this framework for both non-relativistic and relativistic magnetohydrodynamics (MHD). These modules adopt mature and robust algorithms originally developed for the Athena MHD code and incorporate new extensions: support for curvilinear coordinates, higher-order time integrators, more realistic physics such as a general equation of state, and diffusion terms that can be integrated with super-time-stepping algorithms. The modules show excellent performance and scaling, with well over 80% parallel efficiency on over half a million threads. The source code has been made publicly available.

Scattered light imaging has revealed nearly a dozen circumstellar disks around young Herbig Ae/Be stars$-$enabling studies of structures in the upper disk layers as potential signs of on-going planet formation. We present the first images of the disk around the variable Herbig Ae star PDS 201 (V* V351 Ori), and an analysis of the images and spectral energy distribution through 3D Monte-Carlo radiative transfer simulations and forward modelling. The disk is detected in three datasets with LBTI/LMIRCam at the LBT, including direct observations in the $Ks$ and $L'$ filters, and an $L'$ observation with the 360$^\circ$ vector apodizing phase plate coronagraph. The scattered light disk extends to a very large radius of $\sim$250 au, which places it among the largest of such disks. Exterior to the disk, we establish detection limits on substellar companions down to $\sim$5 M$_{Jup}$ at $\gtrsim$1.5" ($\gtrsim$500 au), assuming the Baraffe et al. (2015) models. The images show a radial gap extending to $\sim$0.4" ($\sim$140 au at a distance of 340 pc) that is also evident in the spectral energy distribution. The large gap is a possible signpost of multiple high-mass giant planets at orbital distances ($\sim$60-100 au) that are unusually massive and widely-separated compared to those of planet populations previously inferred from protoplanetary disk substructures.

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.

The Large Binocular Telescope Interferometer (LBTI) enables nulling interferometric observations across the N band (8 to 13 um) to suppress a star's bright light and probe for faint circumstellar emission. We present and statistically analyze the results from the LBTI/HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey for exozodiacal dust. By comparing our measurements to model predictions based on the Solar zodiacal dust in the N band, we estimate a 1 sigma median sensitivity of 23 zodis for early type stars and 48 zodis for Sun-like stars, where 1 zodi is the surface density of habitable zone (HZ) dust in the Solar system. Of the 38 stars observed, 10 show significant excess. A clear correlation of our detections with the presence of cold dust in the systems was found, but none with the stellar spectral type or age. The majority of Sun-like stars have relatively low HZ dust levels (best-fit median: 3 zodis, 1 sigma upper limit: 9 zodis, 95% confidence: 27 zodis based on our N band measurements), while ~20% are significantly more dusty. The Solar system's HZ dust content is consistent with being typical. Our median HZ dust level would not be a major limitation to the direct imaging search for Earth-like exoplanets, but more precise constraints are still required, in particular to evaluate the impact of exozodiacal dust for the spectroscopic characterization of imaged exo-Earth candidates.

FU Ori is the prototype of FU Orionis systems which are outbursting protoplanetary disks. Magnetic fields in FU Ori's accretion disks have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori's inner accretion disk. We find that (1) when the disk is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z$\sim$R, similar to the "surface accretion" mechanism in previous locally-isothermal MHD simulations. (2) A moderate disk wind is launched in the vertical field simulations with a terminal speed of $\sim$300-500 km/s and a mass loss rate of 1-10\% the disk accretion rate, which is consistent with observations. Disk wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disk photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems may reveal very different magnetic field structures. (4) Due to energy loss by the disk wind, the disk photosphere temperature is lower than that predicted by the thin disk theory, and the previously inferred disk accretion rate may be lower than the real accretion rate by a factor of $\sim$2-3.

We study the flow structure in 3D magnetohydrodynamic (MHD) simulations of accretion onto Sagittarius A* via the magnetized winds of the orbiting Wolf-Rayet stars. These simulations cover over 3 orders of magnitude in radius to reach $\approx$ 300 gravitational radii, with only one poorly constrained parameter (the magnetic field in the stellar winds). Even for winds with relatively weak magnetic fields (e.g., plasma $\beta$ $\sim$ $10^6$), flux freezing/compression in the inflowing gas amplifies the field to $\beta$ $\sim$ few well before it reaches the event horizon. Overall, the dynamics, accretion rate, and spherically averaged flow profiles (e.g., density, velocity) in our MHD simulations are remarkably similar to analogous hydrodynamic simulations. We attribute this to the broad distribution of angular momentum provided by the stellar winds, which sources accretion even absent much angular momentum transport. We find that the magneto-rotational instability is not important because of i) strong magnetic fields that are amplified by flux freezing/compression, and ii) the rapid inflow/outflow times of the gas and inefficient radiative cooling preclude circularization. The primary effect of magnetic fields is that they drive a polar outflow that is absent in hydrodynamics. The dynamical state of the accretion flow found in our simulations is unlike the rotationally supported tori used as initial conditions in horizon scale simulations, which could have implications for models being used to interpret Event Horizon Telescope and GRAVITY observations of Sgr A*.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

Recent high resolution observations unveil ring structures in circumstellar disks. The origin of these rings has been widely investigated under various theoretical scenarios. In this work we perform global 3D non-ideal MHD simulations including effects from both Ohmic resistivity and ambipolar diffusion (AD) to model the HL Tau disk. The non-ideal MHD diffusion profiles are calculated based on the global dust evolution calculation including sintering effects. Disk ionization structure changes dramatically across the snow line due to the change of dust size distribution close to snow line of major volatiles. We find that accretion is mainly driven by disk wind. Gaps and rings can be quickly produced from different accretion rates across snow line. Furthermore, ambipolar diffusion (AD) leads to highly preferential accretion at midplane, followed by magnetic reconnection. This results a local zone of decretion that drains of mass in the field reconnection area, which leaves a gap and an adjacent ring just outside it. Overall, under the favorable condition, both snow lines and non-ideal MHD effects can lead to gaseous gaps and rings in protoplanetary disks.

We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being the $M_\bullet-T_{\rm x}$, followed by $M_\bullet-L_{\rm x}$. The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the $M_\bullet-\sigma_{\rm e}$) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, $M_\bullet$ correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondi-like models have inconsistent anti-correlation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as sub-dominant over most of the cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion (CCA), the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev-Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions).

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

The Large Binocular Telescope Interferometer (LBTI) has the longest baseline in the world, 22.7 m, for performing astronomical interferometry in Fizeau mode, which involves beam combination in a focal plane and preserves a wide field-of-view. LBTI can operate in this mode at wavelengths of 1.2 to 5 and 8 to 12 {\mu}m, making it a unique platform for carrying out high-resolution imaging of circumstellar disks, evolved stars, solar system objects, and possibly searches for planets, in the thermal infrared. Over the past five years, LBTI has carried out a considerable number of interferometric observations by combining the beams near a pupil plane to carry out nulling interferometry. This mode is useful for measuring small luminosity level offsets, such as those of exozodiacal dust disks. The Fizeau mode, by contrast, is more useful for generating an image of the target because it has more (u, v) (Fourier) plane coverage. However, the Fizeau mode is still in an ongoing process of commissioning. Sensitive Fizeau observations require active phase control, increased automation, and the removal of non-common-path aberrations (NCPA) between the science and phase beams. This increased level of control will increase the fringe contrast, enable longer integrations, and reduce time overheads. We are in the process of writing a correction loop to remove NCPA, and have carried out tests on old and synthetic data. We have also carried out on-sky Fizeau engineering tests in fall 2018 and spring 2019. In this article, we share lessons learned and strategies developed as a result of these tests.

Theoretical studies suggest that a giant planet around the young star MWC 758 could be responsible for driving the spiral features in its circumstellar disk. Here, we present a deep imaging campaign with the Large Binocular Telescope with the primary goal of imaging the predicted planet. We present images of the disk in two epochs in the $L^{\prime}$ filter (3.8 $\mu m$) and a third epoch in the $M^{\prime}$ filter (4.8 $\mu m$). The two prominent spiral arms are detected in each observation, which constitute the first images of the disk at $M^\prime$, and the deepest yet in $L^\prime$ ($\Delta L^\prime=$12.1 exterior to the disk at 5$\sigma$ significance). We report the detection of a S/N$\sim$3.9 source near the end of the Sourthern arm, and, from the source's detection at a consistent position and brightness during multiple epochs, we establish a $\sim$90% confidence-level that the source is of astrophysical origin. We discuss the possibilities that this feature may be a) an unresolved disk feature, and b) a giant planet responsible for the spiral arms, with several arguments pointing in favor of the latter scenario. We present additional detection limits on companions exterior to the spiral arms, which suggest that a $\lesssim$4 M$_{Jup}$ planet exterior to the spiral arms could have escaped detection. Finally, we do not detect the companion candidate interior to the spiral arms reported recently by Reggiani et al. (2018), although forward modelling suggests that such a source would have likely been detected.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

We use 2D (axisymmetric) and 3D hydrodynamic simulations to study Bondi-Hoyle-Lyttleton (BHL) accretion with and without transverse upstream gradients. We mainly focus on the regime of high (upstream) Mach number, weak upstream gradients and small accretor size, which is relevant to neutron star (NS) accretion in wind-fed Supergiant X-ray binaries (SgXBs). We present a systematic exploration of the flow in this regime. When there are no upstream gradients, the flow is always stable regardless of accretor size or Mach number. For finite upstream gradients, there are three main types of behavior: stable flow (small upstream gradient), turbulent unstable flow without a disk (intermediate upstream gradient), and turbulent flow with a disk-like structure (relatively large upstream gradient). When the accretion flow is turbulent, the accretion rate decreases non-convergently as the accretor size decreases. The flow is more prone to instability and the disk is less likely to form than previously expected; the parameters of most observed SgXBs place them in the regime of a turbulent, disk-less accretion flow. Among the SgXBs with relatively well-determined parameters, we find OAO 1657-415 to be the only one that is likely to host a persistent disk (or disk-like structure); this finding is consistent with observations.

We report the first results of the extension of the QMC model for asymmetric dense matter at finite temperature. The effects of temperature on particle composition (including the full baryon octet content) of the core of (proto-)neutron stars, as well as on the equation of state, are studied. We consider both dense matter in chemical equilibrium and matter in which neutrinos are trapped. In order to simulate stellar temperature profiles that increase with density and stellar radius, the entropy per baryon is fixed. Under these conditions, the model predicts that proto-neutron stars are already born with hyperons present at about the threshold density for their appearance in cold neutron stars, reaching about 20% of the baryon content in the center of the most massive star produced.

In this work, we study matter in the cores of proto-neutron stars, focusing on the impact of their composition on the stellar structure. We begin by examining the effects of finite temperature (through a fixed entropy per baryon) and lepton fraction on purely nucleonic matter by making use of the DSH model . We then turn our attention to a relativistic mean-field model containing exotic degrees of freedom, the Chiral Mean Field (CMF) model, again, under the conditions of finite temperature and trapped neutrinos. In the latter, since both hyperons and quarks are found in the cores of large-mass stars, their interplay and the possibility of mixtures of phases is taken into account and analyzed. Finally, we discuss how stellar rotation can affect our results.

We present 2.9-4.1 micron integral field spectroscopy of the L4+L4 brown dwarf binary HD 130948BC, obtained with the Arizona Lenslets for Exoplanet Spectroscopy (ALES) mode of the Large Binocular Telescope Interferometer (LBTI). The HD 130948 system is a hierarchical triple system, in which the G2V primary is joined by two co-orbiting brown dwarfs. By combining the age of the system with the dynamical masses and luminosities of the substellar companions, we can test evolutionary models of cool brown dwarfs and extra-solar giant planets. Previous near-infrared studies suggest a disagreement between HD 130948BC luminosities and those derived from evolutionary models. We obtained spatially-resolved, low-resolution (R~20) L-band spectra of HD 130948B and C to extend the wavelength coverage into the thermal infrared. Jointly using JHK photometry and ALES L-band spectra for HD 130948BC, we derive atmospheric parameters that are consistent with parameters derived from evolutionary models. We leverage the consistency of these atmospheric quantities to favor a younger age (0.50 \pm 0.07 Gyr) of the system compared to the older age (0.79 \pm 0.22 Gyr) determined with gyrochronology in order to address the luminosity discrepancy.

An isovector-scalar meson is incorporated self-consistently into the quark-meson coupling description of nuclear matter and its most prominent effects on the structure of neutron stars are investigated. The recent measurement of GW170817 is used to constrain the strength of the isovector-scalar channel. With the imminent measurements of the neutron star radii in the NICER mission it is particularly notable that the inclusion of the isovector-scalar force has a significant impact. Indeed, the effect of this interaction on the neutron star radii and masses is larger than the uncertainties introduced by variations in the parameters of symmetric nuclear matter at saturation, namely the density, binding energy per nucleon and the symmetry energy. In addition, since the analysis of GW170817 has provided constraints on the binary tidal deformability of merging neutron stars, the predictions for this parameter within the quark-meson coupling model are explored, as well as the moment of inertia and the quadrupole moment of slowly rotating neutron stars.

We use global three dimensional radiation magneto-hydrodynamic simulations to study the properties of inner regions of accretion disks around a 5\times 10^8 solar mass black hole with mass accretion rates reaching 7% and 20% of the Eddington value. This region of the disk is supported by magnetic pressure with surface density significantly smaller than the values predicted by the standard thin disk model but with a much larger disk scale height. The disks do not show any sign of thermal instability over many thermal time scales. More than half of the accretion is driven by radiation viscosity in the optically thin corona region for the lower accretion rate case, while accretion in the optically thick part of the disk is driven by the Maxwell and Reynolds stresses from MRI turbulence. Coronae with gas temperatures > 10^8 K are generated only in the inner \approx 10 gravitational radii in both simulations, being more compact in the higher accretion rate case. In contrast to the thin disk model, surface density increases with increasing mass accretion rate, which causes less dissipation in the optically thin region and a relatively weaker corona. The simulation results may explain the formation of X-ray coronae in Active Galactic Nuclei (AGNs), the compact size of such coronae, and the observed trend of optical to X-ray luminosity with Eddington ratio for many AGNs.

This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed.

We investigate numerical convergence in simulations of magnetically arrested disks around spinning black holes. Using the general-relativistic magnetohydrodynamics code Athena++, we study the same system at four resolutions (up to an effective 512 by 256 by 512 cells) and with two different spatial reconstruction algorithms. The accretion rate and general large-scale structure of the flow agree across the simulations. This includes the amount of magnetic flux accumulated in the saturated state and the ensuing suppression of the magnetorotational instability from the strong field. The energy of the jet and the efficiency with which spin energy is extracted via the Blandford-Znajek process also show convergence. However the spatial structure of the jet shows variation across the set of grids employed, as do the Lorentz factors. Small-scale features of the turbulence, as measured by correlation lengths, are not fully converged. Despite convergence of a number of aspects of the flow, modeling of synchrotron emission shows that variability is not converged and decreases with increasing resolution even at our highest resolutions.

Gaseous disks have been proposed as a mechanism for facilitating mergers of binary black holes. We explore circumbinary disk systems to determine the evolution of the central binary. To do so, we perform 3D, hydrodynamic, locally isothermal simulations of circumbinary disks on a Cartesian grid. We focus on binaries of equal mass ratios on fixed circular orbits. To investigate the orbital evolution of the binary, we examine the various torques exerted on the system. For the case where the disk plane and binary orbital plane are aligned, we find that the total torque is positive so that the semi-major axis of the binary increases. For the misaligned case, we run simulations with the binary orbital plane and disk midplane misaligned by 45\degree and find the same results - the binary grows. The timescale for the circumbinary disk to realign to the plane of the binary is consistent with the global viscous timescale of the disk.

The gyro-resonant cosmic-ray (CR) streaming instability is believed to play a crucial role in CR transport, leading to growth of Alfv\'en waves at small scales that scatter CRs, and impacts the interaction of CRs with the ISM on large scales. However, extreme scale separation ($\lambda \ll \rm pc$), low cosmic ray number density ($n_{\rm CR}/n_{\rm ISM} \sim 10^{-9}$), and weak CR anisotropy ($\sim v_A/c$) pose strong challenges for proper numerical studies of this instability on a microphysical level. Employing the recently developed magnetohydrodynamic-particle-in-cell (MHD-PIC) method, which has unique advantages to alleviate these issues, we conduct one-dimensional simulations that quantitatively demonstrate the growth and saturation of the instability in the parameter regime consistent with realistic CR streaming in the large-scale ISM. Our implementation of the $\delta f$ method dramatically reduces Poisson noise and enables us to accurately capture wave growth over a broad spectrum, equally shared between left and right handed Alfv\'en modes. We are also able to accurately follow the quasi-linear diffusion of CRs subsequent to wave growth, which is achieved by employing phase randomization across periodic boundaries. Full isotropization of the CRs in the wave frame requires pitch angles of most CRs to efficiently cross $90^\circ$, and can be captured in simulations with relatively high wave amplitude and/or high spatial resolution. We attribute this crossing to non-linear wave-particle interaction (rather than mirror reflection) by investigating individual CR trajectories. We anticipate our methodology will open up opportunities for future investigations that incorporate additional physics.

The prevalence of binary stars at close separations implies that many of these systems will interact or merge during the binary's lifetime. This paper presents hydrodynamic simulations of the scenario of binary coalescence through unstable mass transfer, which drives the pair to closer separations. When the donor star does not rotate synchronously with respect to the orbit, dynamical tidal waves are excited in its envelope. We show that resonance crossings with high azimuthal-order $(m\sim3$ to $6$) fundamental modes induce a visible "polygram" distortion to the star. As the binary orbit tightens, the system sweeps through resonance with modes of decreasing azimuthal order, which are selectively excited. We compare our hydrodynamic simulations to predictions from linear theory of resonant mode excitation. The linear theory provides an estimate of mode amplitudes to within a factor of two, even as the oscillations become quite non-linear as the stars coalesce. We estimate that a wave with 10% radial amplitude generates approximately 1% photometric variability; this may be detectible if such a binary coalescence is caught in action by future photometric all-sky surveys.

We present the results of the largest $L^{\prime}$ ($3.8~\mu$m) direct imaging survey for exoplanets to date, the Large Binocular Telescope Interferometer (LBTI) Exozodi Exoplanet Common Hunt (LEECH). We observed 98 stars with spectral types from B to M. Cool planets emit a larger share of their flux in $L^{\prime}$ compared to shorter wavelengths, affording LEECH an advantage in detecting low-mass, old, and cold-start giant planets. We emphasize proximity over youth in our target selection, probing physical separations smaller than other direct imaging surveys. For FGK stars, LEECH outperforms many previous studies, placing tighter constraints on the hot-start planet occurrence frequency interior to $\sim20$ au. For less luminous, cold-start planets, LEECH provides the best constraints on giant-planet frequency interior to $\sim20$ au around FGK stars. Direct imaging survey results depend sensitively on both the choice of evolutionary model (e.g., hot- or cold-start) and assumptions (explicit or implicit) about the shape of the underlying planet distribution, in particular its radial extent. Artificially low limits on the planet occurrence frequency can be derived when the shape of the planet distribution is assumed to extend to very large separations, well beyond typical protoplanetary dust-disk radii ($\lesssim50$ au), and when hot-start models are used exclusively. We place a conservative upper limit on the planet occurrence frequency using cold-start models and planetary population distributions that do not extend beyond typical protoplanetary dust-disk radii. We find that $\lesssim90\%$ of FGK systems can host a 7 to 10 $M_{\mathrm{Jup}}$ planet from 5 to 50 au. This limit leaves open the possibility that planets in this range are common.

On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

The observed rotation measures (RMs) towards the galactic centre magnetar and towards Sagittarius A* provide a strong constraint on MHD models of the galactic centre accretion flow, probing distances from the black hole separated by many orders of magnitude. We show, using 3D simulations of accretion via magnetized stellar winds of the Wolf-Rayet stars orbiting the black hole, that the large, time-variable RM observed for the pulsar PSR J1745-2900 can be explained by magnetized wind-wind shocks of nearby stars in the clockwise stellar disc. In the same simulation, both the total X-ray luminosity integrated over 2-10$''$, the time variability of the magnetar's dispersion measure, and the RM towards Sagittarius A* are consistent with observations. We argue that (in order for the large RM of the pulsar to not be a priori unlikely) the pulsar should be on an orbit that keeps it near the clockwise disc of stars. We present a 2D RM map of the central 1/2 parsec of the galactic centre that can be used to test our models. Our simulations predict that Sgr A* is typically accreting a significantly ordered magnetic field that ultimately could result in a strongly magnetized flow with flux threading the horizon at $\sim$ 10$\%$ of the magnetically arrested limit.

We present detections of methane in R of $\sim$1300, L band spectra of VHS 1256 b and PSO 318.5, two low gravity, red, late L dwarfs that share the same colors as the HR 8799 planets. These spectra reveal shallow methane features, which indicate VHS 1256 b and PSO 318.5 have photospheres that are out of chemical equilibrium. Directly imaged exoplanets usually have redder near infrared colors than the field-age population of brown dwarfs on a color magnitude diagram. These objects along the L to T transition show reduced methane absorption and evidence of photospheric clouds. Compared to the H and K bands, L band (3 micron - 4 micron) spectroscopy provides stronger constraints on the methane abundances of brown dwarfs and directly imaged exoplanets that have similar effective temperatures as L to T transition objects. When combined with near infrared spectra, the L band extends our conventional wavelength coverage, increasing our understanding of atmospheric cloud structure. Our model comparisons show relatively strong vertical mixing and photospheric clouds can explain the molecular absorption features and continua of VHS 1256 b and PSO 318.5. We also discuss the implications of this work for future exoplanet focused instruments and observations with the James Webb Space Telescope.

Recent observations have revealed that the remnants of stellar-coalescence transients are bipolar. This raises the questions of how these bipolar morphologies arise and what they teach us about the mechanisms of mass ejection during stellar mergers and common-envelope phases. In this paper, we analyze hydrodynamic simulations of the lead-in to binary coalescence, a phase of unstable Roche lobe overflow that takes the binary from the Roche limit separation to the engulfment of the more compact accretor within the envelope of the extended donor. As mass transfer runs away at increasing rates, gas trails away from the binary. Contrary to previous expectations, early mass loss from the system remains bound to the binary and forms a circumbinary torus. Later ejecta, generated as the accretor grazes the surface of the donor, have very different morphologies and are unbound. These two components of mass loss from the binary interact as later, higher-velocity ejecta collide with the circumbinary torus formed by earlier mass loss. Unbound ejecta are redirected toward the poles, and escaping material creates a bipolar outflow. Our findings show that the transition from bound to unbound ejecta from coalescing binaries can explain the bipolar nature of their remnants, with implications for our understanding of the origin of bipolar remnants of stellar-coalescence transients and, perhaps, some preplanetary nebulae.

The X-ray signal from hydrogen-rich supernovae (SNe II) in the first tens to hundreds of days after the shock breakout encodes important information about the circumstellar material (CSM) surrounding their progenitors before explosion. In this study, we describe a way to generate the SN II X-ray light curves from hydrodynamical simulations performed with the code Athena++, using the X-ray package XSPEC. In addition, we employ a radiation diffusion hydrodynamic code SNEC for generating the optical light curves in different bands. In this numerical setup, we model the X-ray and optical emission from a set of progenitor models, consisting of either two (red supergiant + low density steady wind), or three (red supergiant + dense CSM + low density steady wind) components. We vary the density in the wind and the slope in the CSM to see how these parameters influence the resulting X-ray and optical light curves. Among our models, we identify one that is able to roughly reproduce both optical and X-ray data of the well observed SN 2013ej. In order to achieve this, the slope of the dense CSM in this model should be steeper than the one of a steady wind ($\rho\propto r^{-2}$), and closer to $\rho\propto r^{-5}$. On the other hand, we show that too steep and extended CSM profiles may produce excessive X-ray emission in the first few tens of days, up to a few orders of magnitude larger than observed. We conclude that ability to reproduce the observed X-ray signal from SNe~II together with their optical light curves is crucial in establishing the validity of different CSM models.

(Abridged) Stars more massive than $20-30M_{\odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff} \approx 2 - 4 \times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff} \approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5} M_{\odot}/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars.

[abridged] The presence of large amounts of dust in the habitable zones of nearby stars is a significant obstacle for future exo-Earth imaging missions. We executed an N band nulling interferometric survey to determine the typical amount of such exozodiacal dust around a sample of nearby main sequence stars. The majority of our data have been analyzed and we present here an update of our ongoing work. We find seven new N band excesses in addition to the high confidence confirmation of three that were previously known. We find the first detections around Sun-like stars and around stars without previously known circumstellar dust. Our overall detection rate is 23%. The inferred occurrence rate is comparable for early type and Sun-like stars, but decreases from 71% [+11%/-20%] for stars with previously detected mid- to far-infrared excess to 11% [+9%/-4%] for stars without such excess, confirming earlier results at high confidence. For completed observations on individual stars, our sensitivity is five to ten times better than previous results. Assuming a lognormal luminosity function of the dust, we find upper limits on the median dust level around all stars without previously known mid to far infrared excess of 11.5 zodis at 95% confidence level. The corresponding upper limit for Sun-like stars is 16 zodis. An LBTI vetted target list of Sun-like stars for exo-Earth imaging would have a corresponding limit of 7.5 zodis. We provide important new insights into the occurrence rate and typical levels of habitable zone dust around main sequence stars. Exploiting the full range of capabilities of the LBTI provides a critical opportunity for the detailed characterization of a sample of exozodiacal dust disks to understand the origin, distribution, and properties of the dust.

We present a fourth-order accurate finite volume method for the solution of ideal magnetohydrodynamics (MHD). The numerical method combines high-order quadrature rules in the solution of semi-discrete formulations of hyperbolic conservation laws with the upwind constrained transport (UCT) framework to ensure that the divergence-free constraint of the magnetic field is satisfied. A novel implementation of UCT that uses the piecewise parabolic method (PPM) for the reconstruction of magnetic fields at cell corners in 2D is introduced. The resulting scheme can be expressed as the extension of the second-order accurate constrained transport (CT) Godunov-type scheme that is currently used in the Athena astrophysics code. After validating the base algorithm on a series of hydrodynamics test problems, we present the results of multidimensional MHD test problems which demonstrate formal fourth-order convergence for smooth problems, robustness for discontinuous problems, and improved accuracy relative to the second-order scheme.

The Arizona Lenslets for Exoplanet Spectroscopy (ALES) is the world's first AO-fed thermal infrared integral field spectrograph, mounted inside the Large Binocular Telescope Interferometer (LBTI) on the LBT. An initial mode of ALES allows 3-4 micron spectra at R~20 with 0.026'' spaxels over a 1''x1'' field-of-view. We are in the process of upgrading ALES with additional wavelength ranges, spectral resolutions, and plate scales allowing a broad suite of science that takes advantage of ALES's unique ability to work at wavelengths >2 microns, and at the diffraction limit of the LBT's full 23.8 meter aperture.

Magnetohydrodynamic instabilities play an important role in accretion disks systems. Besides the well-known effects of the magnetorotational instability (MRI), the Parker-Rayleigh-Taylor instability (PRTI) also arises as an important mechanism to help in the formation of the coronal region around an accretion disk and in the production of magnetic reconnection events similar to those occurring in the solar corona. In this work, we have performed three-dimensional magnetohydrodynamical (3D-MHD) shearing-box numerical simulations of accretion disks with an initial stratified density distribution and a strong azimuthal magnetic field with a ratio between the thermal and magnetic pressures of the order of unity. This study aimed at verifying the role of these instabilities in driving fast magnetic reconnection in turbulent accretion disk/corona systems. As we expected, the simulations showed an initial formation of large-scale magnetic loops due to the PRTI followed by the development of a nearly steady-state turbulence driven by both instabilities. In this turbulent environment, we have employed an algorithm to identify the presence of current sheets produced by the encounter of magnetic flux ropes of opposite polarity in the turbulent regions of both the corona and the disk. We computed the magnetic reconnection rates in these locations obtaining average reconnection velocities in Alfv\'en speed units of the order of $0.13 \pm 0.09$ in the accretion disk and $0.17 \pm 0.10$ in the coronal region (with mean peak values of order of $0.2$), which are consistent with the predictions of the theory of turbulence-induced fast reconnection.

We present the data reduction pipeline, MEAD, for Arizona Lenslets for Exoplanet Spectroscopy (ALES), the first thermal infrared integral field spectrograph designed for high-contrast imaging. ALES is an upgrade of LMIRCam, the $1-5\,\mu$m imaging camera for the Large Binocular Telescope, capable of observing astronomical objects in the thermal infrared ($3-5\,\mu$m) to produce simultaneous spatial and spectral data cubes. The pipeline is currently designed to perform $L$-band ($2.8-4.2\,\mu$m) data cube reconstruction, relying on methods used extensively by current near-infrared integral field spectrographs. ALES data cube reconstruction on each spectra uses an optimal extraction method. The calibration unit comprises a thermal infrared source, a monochromator and an optical diffuser designed to inject specific wavelengths of light into LBTI to evenly illuminate the pupil plane and ALES lenslet array with monochromatic light. Not only does the calibration unit facilitate wavelength calibration for ALES and LBTI, but it also provides images of monochromatic point spread functions (PSFs). A linear combination of these monochromatic PSFs can be optimized to fit each spectrum in the least-square sense via $\chi^2$ fitting.

The integral field spectrograph configuration of the LMIRCam science camera within the Large Binocular Telescope Interferometer (LBTI) facilitates 2 to 5 $\mu$m spectroscopy of directly imaged gas-giant exoplanets. The mode, dubbed ALES, comprises magnification optics, a lenslet array, and direct-vision prisms, all of which are included within filter wheels in LMIRCam. Our observing approach includes manual adjustments to filter wheel positions to optimize alignment, on/off nodding to track sky-background variations, and wavelength calibration using narrow band filters in series with ALES optics. For planets with separations outside our 1"x1" field of view, we use a three-point nod pattern to visit the primary, secondary and sky. To minimize overheads we select the longest exposure times and nod periods given observing conditions, especially sky brightness and variability. Using this strategy we collected several datasets of low-mass companions to nearby stars.

Luminous red nova transients, presumably from stellar coalescence, exhibit long-term precursor emission over hundreds of binary orbits, leading to impulsive outbursts with durations similar to a single orbital period. In an effort to understand these signatures, we present and analyze a hydrodynamic model of unstable mass transfer from a giant-star donor onto a more compact accretor in a binary system. Our simulation begins with mass transfer at the Roche limit separation and traces a phase of runaway decay leading up to the plunge of the accretor within the envelope of the donor. We characterize the fluxes of mass and angular momentum through the system and show that the orbital evolution can be reconstructed from measurements of these quantities. The morphology of outflow from the binary changes significantly as the binary orbit tightens. At wide separations, a thin stream of relatively high-entropy gas trails from the outer Lagrange points. As the orbit tightens, the orbital motion desynchronizes from the donor's rotation, and low-entropy ejecta trace a broad fan of largely ballistic trajectories. An order-of-magnitude increase in mass ejection rate accompanies the plunge of the accretor with the envelope of the donor. We argue that this transition marks the precursor-to-outburst transition observed in stellar coalescence transients.

Monitoring the long-term radial velocity (RV) and acceleration of nearby stars has proven an effective method for directly detecting binary and substellar companions. Some fraction of nearby RV trend systems are expected to be comprised of compact objects that likewise induce a systemic Doppler signal. In this paper, we report the discovery of a white dwarf companion found to orbit the nearby ($\pi = 28.297 \pm 0.066$ mas) G9 V star HD 169889. High-contrast imaging observations using NIRC2 at Keck and LMIRCam at the LBT uncover the ($\Delta H = 9.76 \pm 0.16$, $\Delta L' = 9.60 \pm 0.03$) companion at an angular separation of 0.8'' (28 au). Thirteen years of precise Doppler observations reveal a steep linear acceleration in RV time series and place a dynamical constraint on the companion mass of $M \geq 0.369 \pm 0.010 M_{\odot}$. This "Sirius-like" system adds to the census of white dwarf companions suspected to be missing in the solar neighborhood.

We present Athena++ grid-based, hydrodynamic simulations of accretion onto Sagittarius A* via the stellar winds of the $\sim 30$ Wolf-Rayet stars within the central parsec of the galactic center. These simulations span $\sim$ 4 orders of magnitude in radius, reaching all the way down to 300 gravitational radii of the black hole, $\sim 32$ times further in than in previous work. We reproduce reasonably well the diffuse thermal X-ray emission observed by Chandra in the central parsec. The resulting accretion flow at small radii is a superposition of two components: 1) a moderately unbound, sub-Keplerian, thick, pressure-supported disc that is at most (but not all) times aligned with the clockwise stellar disc, and 2) a bound, low-angular momentum inflow that proceeds primarily along the southern pole of the disc. We interpret this structure as a natural consequence of a few of the innermost stellar winds dominating accretion, which produces a flow with a broad distribution of angular momentum. Including the star S2 in the simulation has a negligible effect on the flow structure. Extrapolating our results from simulations with different inner radii, we find an accretion rate of $\sim$ a few $\times 10^{-8} M_\odot$/yr at the horizon scale, consistent with constraints based on modeling the observed emission of Sgr A*. The flow structure found here can be used as more realistic initial conditions for horizon scale simulations of Sgr A*.

The Large Binocular Telescope Interferometer (LBTI) can perform Fizeau interferometry in the focal plane, which accesses spatial information out to the LBT's full 22.7-m edge-to-edge baseline. This mode has previously been used to obtain science data, but has been limited to observations where the optical path difference (OPD) between the two beams is not controlled, resulting in unstable fringes on the science detectors. To maximize the science return, we are endeavoring to stabilize the OPD and tip-tilt variations and make the LBTI Fizeau mode optimized and routine. Here we outline the optical configuration of LBTI's Fizeau mode and our strategy for commissioning this observing mode.

We report results from very long baseline interferometric (VLBI) observations of the supermassive black hole in the Galactic center, Sgr A*, at 1.3 mm (230 GHz). The observations were performed in 2013 March using six VLBI stations in Hawaii, California, Arizona, and Chile. Compared to earlier observations, the addition of the APEX telescope in Chile almost doubles the longest baseline length in the array, provides additional {\it uv} coverage in the N-S direction, and leads to a spatial resolution of $\sim$30 $\mu$as ($\sim$3 Schwarzschild radii) for Sgr A*. The source is detected even at the longest baselines with visibility amplitudes of $\sim$4-13% of the total flux density. We argue that such flux densities cannot result from interstellar refractive scattering alone, but indicate the presence of compact intrinsic source structure on scales of $\sim$3 Schwarzschild radii. The measured nonzero closure phases rule out point-symmetric emission. We discuss our results in the context of simple geometric models that capture the basic characteristics and brightness distributions of disk- and jet-dominated models and show that both can reproduce the observed data. Common to these models are the brightness asymmetry, the orientation, and characteristic sizes, which are comparable to the expected size of the black hole shadow. Future 1.3 mm VLBI observations with an expanded array and better sensitivity will allow a more detailed imaging of the horizon-scale structure and bear the potential for a deep insight into the physical processes at the black hole boundary.

The HOSTS (Hunt for Observable Signatures of Terrestrial Systems) survey searches for dust near the habitable zones (HZs) around nearby, bright main sequence stars. We use nulling interferometry in N band to suppress the bright stellar light and to probe for low levels of HZ dust around the 30 stars observed so far. Our overall detection rate is 18%, including four new detections, among which are the first three around Sun-like stars and the first two around stars without any previously known circumstellar dust. The inferred occurrence rates are comparable for early type and Sun-like stars, but decrease from 60 (+16/-21)% for stars with previously detected cold dust to 8 (+10/-3)% for stars without such excess, confirming earlier results at higher sensitivity. For completed observations on individual stars, our sensitivity is five to ten times better than previous results. Assuming a lognormal excess luminosity function, we put upper limits on the median HZ dust level of 13 zodis (95% confidence) for a sample of stars without cold dust and of 26 zodis when focussing on Sun-like stars without cold dust. However, our data suggest that a more complex luminosity function may be more appropriate. For stars without detectable LBTI excess, our upper limits are almost reduced by a factor of two, demonstrating the strength of LBTI target vetting for future exo-Earth imaging missions. Our statistics are so far limited and extending the survey is critical to inform the design of future exo-Earth imaging surveys.

We report the results of a neutrino search in Super-Kamiokande for coincident signals with the first detected gravitational wave produced by a binary neutron star merger, GW170817, which was followed by a short gamma-ray burst, GRB170817A, and a kilonova/macronova. We searched for coincident neutrino events in the range from 3.5 MeV to $\sim$100 PeV, in a time window $\pm$500 seconds around the gravitational wave detection time, as well as during a 14-day period after the detection. No significant neutrino signal was observed for either time window. We calculated 90% confidence level upper limits on the neutrino fluence for GW170817. From the upward-going-muon events in the energy region above 1.6 GeV, the neutrino fluence limit is $16.0^{+0.7}_{-0.6}$ ($21.3^{+1.1}_{-0.8}$) cm$^{-2}$ for muon neutrinos (muon antineutrinos), with an error range of $\pm5^{\circ}$ around the zenith angle of NGC4993, and the energy spectrum is under the assumption of an index of $-2$. The fluence limit for neutrino energies less than 100 MeV, for which the emission mechanism would be different than for higher-energy neutrinos, is also calculated. It is $6.6 \times 10^7$ cm$^{-2}$ for anti-electron neutrinos under the assumption of a Fermi-Dirac spectrum with average energy of 20 MeV.

We perform relativistic hydrodynamic simulations of internal shocks formed in microquasar jets by continuous variation of the bulk Lorentz factor, in order to investigate the internal shock model. We consider one-, two-, and flicker noise 20-mode variability. We observe emergence of a forward-reverse shock structure for each peak of the Lorentz factor modulation. The high pressure in the shocked layer launches powerful outflows perpendicular to the jet beam into the ambient medium. These outflows dominate the details of the jet's kinetic energy thermalization. They are responsible for mixing between the jet and surrounding medium and generate powerful shocks in the latter. These results do not concur with the popular picture of well-defined internal shells depositing energy as they collide within the confines of the jet, in fact collisions between internal shells themselves are quite rare in our continuous formulation of the problem. For each of our simulations, we calculate the internal energy deposited in the system, the "efficiency" of this deposition (defined as the ratio of internal to total flow energy), and the maximum temperature reached in order to make connections to emission mechanisms. We probe the dependence of these diagnostics on the Lorentz factor variation amplitudes, modulation frequencies, as well as the initial density ratio between the jet and the ambient medium.

We report new global ideal MHD simulations for thin accretion disks (with thermal scale height H/R=0.1 and 0.05) threaded by net vertical magnetic fields. Our computations span three orders of magnitude in radius, extend all the way to the pole, and are evolved for more than one viscous time over the inner decade in radius. Static mesh refinement is used to properly resolve MRI. We find that:(1) inward accretion occurs mostly in the upper magnetically dominated regions of the disk, similar to the predictions from some previous analytical work and the "coronal accretion" in previous GRMHD simulations. Rapid inflow in the upper layers combined with slow outflow at the midplane creates strong $R\phi$ and $z\phi$ stresses in the mean field; the vertically integrated $\alpha\sim 0.5-1$ when the initial field has $\beta_{0}=10^3$ at the midplane. (2) A quasi-static global field geometry is established in which flux transport by inflows at the surface is balanced by turbulent diffusion. The field is strongly pinched inwards at the surface. A steady-state advection-diffusion model, with turbulent magnetic Prandtl number of order unity, reproduces this geometry well. (3) Weak unsteady disk winds are launched at $z/R\sim1$ with the Alfven radius $R_{A}/R_{0}\sim3$. Although the wind is episodic, the time averaged properties are well described by steady wind theory. Wind is not efficient at transporting angular momentum. Even with $\beta_{0}=10^3$, only 5% of the angular momentum transport is driven by torque from the wind, and the wind mass flux from the inner decade of radius is only $\sim$ 0.4% of the mass accretion rate. With weaker fields or thinner disks, the wind contributes even less. (4) Most of the disk accretion is driven by the $R\phi$ stress from the MRI and global magnetic fields. Our simulations have many applications to astrophysical accretion disk systems.

We propose a novel method to constrain turbulence and bulk motions in massive galaxies, groups and clusters, exploring both simulations and observations. As emerged in the recent picture of the top-down multiphase condensation, the hot gaseous halos are tightly linked to all other phases in terms of cospatiality and thermodynamics. While hot halos (10^7 K) are perturbed by subsonic turbulence, warm (10^4 K) ionized and neutral filaments condense out of the turbulent eddies. The peaks condense into cold molecular clouds (< 100 K) raining in the core via chaotic cold accretion (CCA). We show all phases are tightly linked via the ensemble (wide-aperture) velocity dispersion along the line of sight. The correlation arises in complementary long-term AGN feedback simulations and high-resolution CCA runs, and is corroborated by the combined Hitomi and new IFU measurements in Perseus cluster. The ensemble multiphase gas distributions are characterized by substantial spectral line broadening (100-200 km/s) with mild line shift. On the other hand, pencil-beam detections sample the small-scale clouds displaying smaller broadening and significant line shift up to several 100 km/s, with increased scatter due to the turbulence intermittency. We present new ensemble sigma_v of the warm Halpha+[NII] gas in 72 observed cluster/group cores: the constraints are consistent with the simulations and can be used as robust proxies for the turbulent velocities, in particular for the challenging hot plasma (otherwise requiring extremely long X-ray exposures). We show the physically motivated criterion C = t_cool/t_eddy ~ 1 best traces the condensation extent region and presence of multiphase gas in observed clusters/groups. The ensemble method can be applied to many available datasets and can substantially advance our understanding of multiphase halos in light of the next-generation multiwavelength missions.

The interaction of a satellite with a gaseous disk results in the excitation of spiral density waves which remove angular momentum from the orbit. In addition, if the orbit is not coplanar with the disk, three-dimensional effects will excite bending and eccentricity waves. We perform three-dimensional hydrodynamic simulations to study nonlinear disk-satellite interaction in inviscid protoplanetary disks for a variety of orbital inclinations from $0^\circ$ to $180^\circ$. It is well known that three-dimensional effects are important even for zero inclination. In this work we (1) show that for planets with small inclinations (as in the Solar system), effects such as the total torque and migration rate strongly depend on the inclination and are significantly different (about 2.5 times smaller) from the two-dimensional case, (2) give formulae for the migration rate, inclination damping, and precession rate of planets with different inclination angles in disk with different scale heights, and (3) present the observational signatures of a planet on an inclined orbit with respect to the protoplanetary disk. For misaligned planets we find good agreement with linear theory in the limit of small inclinations, and with dynamical friction estimates for intermediate inclinations. We find that in the latter case, the dynamical friction force is not parallel to the relative planetary velocity. Overall, the derived formulae will be important for studying exoplanets with obliquity.

We perform two-dimensional and three-dimensional radiation hydrodynamic simulations to study cold clouds accelerated by radiation pressure on dust in the environment of rapidly star-forming galaxies dominated by infrared flux. We utilize the reduced speed of light approximation to solve the frequency-averaged, time-dependent radiative transfer equation. We find that radiation pressure is capable of accelerating the clouds to hundreds of kilometers per second while remaining dense and cold, consistent with observations. We compare these results to simulations where acceleration is provided by entrainment in a hot wind, where the momentum injection of the hot flow is comparable to the momentum in the radiation field. We find that the survival time of the cloud accelerated by the radiation field is significantly longer than that of a cloud entrained in a hot outflow. We show that the dynamics of the irradiated cloud depends on the initial optical depth, temperature of the cloud, and the intensity of the flux. Additionally, gas pressure from the background may limit cloud acceleration if the density ratio between the cloud and background is $\lesssim 10^{2}$. In general, a 10 pc-scale optically thin cloud forms a pancake structure elongated perpendicular to the direction of motion, while optically thick clouds form a filamentary structure elongated parallel to the direction of motion. The details of accelerated cloud morphology and geometry can also be affected by other factors, such as the cloud lengthscale, the reduced speed of light approximation, spatial resolution, initial cloud structure, and the dimensionality of the run, but these have relatively little affect on the cloud velocity or survival time.

We present the results of a search in the Super-Kamiokande (SK) detector for excesses of neutrinos with energies above a few GeV that are in the direction of the track events reported in IceCube. Data from all SK phases (SK-I through SK-IV) were used, spanning a period from April 1996 to April 2016 and corresponding to an exposure of 225 kilotonne-years . We considered the 14 IceCube track events from a data set with 1347 livetime days taken from 2010 to 2014. We use Poisson counting to determine if there is an excess of neutrinos detected in SK in a 10 degree search cone (5 degrees for the highest energy data set) around the reconstructed direction of the IceCube event. No significant excess was found in any of the search directions we examined. We also looked for coincidences with a recently reported IceCube multiplet event. No events were detected within a $\pm$ 500 s time window around the first detected event, and no significant excess was seen from that direction over the lifetime of SK.

We use global three dimensional radiation magneto-hydrodynamical simulations to study accretion disks onto a $5\times 10^8M_{\odot}$ black hole with accretion rates varying from $\sim 250L_{Edd}/c^2$ to $1500 L_{Edd}/c^2$. We form the disks with torus centered at $50-80$ gravitational radii with self-consistent turbulence initially generated by the magneto-rotational instability. We study cases with and without net vertical magnetic flux. The inner regions of all disks have radiation pressure $\sim 10^4-10^6$ times the gas pressure. Non-axisymmetric density waves that steepen into spiral shocks form as gas flows towards the black hole. In simulations without net vertical magnetic flux, Reynolds stress generated by the spiral shocks are the dominant mechanism to transfer angular momentum. Maxwell stress from MRI turbulence can be larger than the Reynolds stress only when net vertical magnetic flux is sufficiently large. Outflows are formed with speed $\sim 0.1-0.4c$. When the accretion rate is smaller than $\sim 500 L_{Edd}/c^2$, outflows start around $10$ gravitational radii and the radiative efficiency is $\sim 5\%-7\%$ with both magnetic field configurations. With accretion rate reaching $1500 L_{Edd}/c^2$, most of the funnel region close to the rotation axis becomes optically thick and the outflow only develops beyond $50$ gravitational radii. The radiative efficiency is reduced to $1\%$. We always find the kinetic energy luminosity associated with the outflow is only $\sim 15\%-30\%$ of the radiative luminosity. The mass flux lost in the outflow is $\sim 15\%-50\%$ of the net mass accretion rates. We discuss implications of our simulation results on the observational properties of these disks.

The early optical afterglow emission of several gamma-ray bursts (GRBs) shows a high linear polarization degree (PD) of tens of percent, suggesting an ordered magnetic field in the emission region. The light curves are consistent with being of a reverse shock (RS) origin. However, the magnetization parameter, $\sigma$, of the outflow is unknown. If $\sigma$ is too small, an ordered field in the RS may be quickly randomized due to turbulence driven by various perturbations so that the PD may not be as high as observed. Here we use the "Athena++" relativistic MHD code to simulate a relativistic jet with an ordered magnetic field propagating into a clumpy ambient medium, with a focus on how density fluctuations may distort the ordered magnetic field and reduce PD in the RS emission for different $\sigma$ values. For a given density fluctuation, we discover a clear power-law relationship between the relative PD reduction and the $\sigma$ value of the outflow. Such a relation may be applied to estimate $\sigma$ of the GRB outflows using the polarization data of early afterglows.

We perform global three-dimensional MHD simulations of unstratified accretion disks in cataclysmic variables (CVs). By including mass inflow via an accretion stream, we are able to evolve the disk to a steady state. We investigate the relative importance of spiral shocks and the magnetorotational instability (MRI) in driving angular momentum transport and how each depend on the geometry and strength of seed magnetic field and the Mach number of the disk (where Mach number is ratio of the azimuthal velocity and the sound speed of gas). We use a locally isothermal equation of state and adopt temperature profiles that are consistent with CV disk observations. Our results indicate that the relative importance of spiral shocks and MRI in driving angular momentum transport is controlled by the gas Mach number and the seed magnetic field strength. MRI and spiral shocks provide comparable efficiency of angular momentum transport when the disk Mach number is around 10 and the seed magnetic field has plasma $\beta=400$ (where $\beta$ is ratio of gas pressure and magnetic pressure). The MRI dominates whenever the seed field strength, or the disk Mach number, is increased. Among all of our simulations, the effective viscosity parameter $\alpha_{eff} \sim 0.016-0.1$ after MRI saturates and the disk reaches steady state. Larger values of $\alpha_{eff}$ are favored when the seed magnetic field has vertical components or the flow has stronger magnetization ($1/\beta$). Our models all indicate that the role of MRI in driving angular momentum transport thus mass accretion in CV disks is indispensable, especially in cool disks with weak spiral shocks.

Gravitationally focused wind accretion in binary systems consisting of an evolved star with a gaseous envelope and a compact accreting companion is a possible mechanism to explain mass transfer in symbiotic binaries. We study the mass accretion around the secondary caused by the strong wind from the primary late-type component using global three-dimensional hydrodynamic numerical simulations during quiescence and outburst stages. In particular, the dependence of the mass accretion rate on the mass-loss rate, wind parameters and phases of wind outburst development is considered. For a typical wind from an asymptotic giant branch star with a mass-loss rate of 1e-6 Msun/year and wind speeds of 20-50 km/s, the mass transfer through a focused wind results in efficient infall on to the secondary. Accretion rates onto the secondary of 5-20 per cent of the mass-loss from the primary are obtained during quiescence and outburst periods where the wind velocity and mass-loss rates are varied, about 20-50 per cent larger than in the standard Bondi-Hoyle-Lyttleton approximation. This mechanism could be an important method for explaining observed accretion luminosities and periodic modulations in the accretion rates for a broad range of interacting binary systems.

The Alice instrument on NASA's New Horizons spacecraft observed an ultraviolet solar occultation by Pluto's atmosphere on 2015 July 14. The transmission vs. altitude was sensitive to the presence of N2, CH4, C2H2, C2H4, C2H6, and haze. We derived line-of-sight abundances and local number densities for the 5 molecular species, and line-of-sight optical depth and extinction coefficients for the haze. We found the following major conclusions: 1) We confirmed temperatures in Pluto's upper atmosphere that were colder than expected before the New Horizons flyby, with upper atmospheric temperatures near 65-68 K. The inferred enhanced Jeans escape rates were (3e22-7e22) N2/s and (4e25-8e25) CH4/s at the exobase (at a radius of ~2900 km, or an altitude of ~1710 km). 2) We measured CH4 abundances from 80 to 1200 km above the surface. A joint analysis of the Alice CH4 and Alice and REX N2 measurements implied a very stable lower atmosphere with a small eddy diffusion coefficient, most likely between 550 and 4000 cm2/s. Such a small eddy diffusion coefficient placed the homopause within 12 km of the surface, giving Pluto a small planetary boundary layer. The inferred CH4 surface mixing ratio was ~0.28-0.35%. 3) The abundance profiles of the C2Hx hydrocarbons (C2H2, C2H4, C2H6) were not simply exponential with altitude. We detected local maxima in line-of-sight abundance near 410 km altitude for C2H4, near 320 km for C2H2, and an inflection point or the suggestion of a local maximum at 260 km for C2H6. We also detected local minima near 200 km altitude for C2H4, near 170 km for C2H2, and an inflection point or minimum near 170-200 km for C2H6. These compared favorably with models for hydrocarbon production near 300-400 km and haze condensation near 200 km, especially for C2H2 and C2H4 (Wong et al. 2017). 4) We found haze that had an extinction coefficient approximately proportional to N2 density.

Circumplanetary disks (CPDs) control the growth of planets, supply material for satellites to form, and provide observational signatures of young forming planets. We have carried out two dimensional hydrodynamical simulations with radiative cooling to study CPDs, and suggested a new mechanism to drive the disk accretion. Two spiral shocks are present in CPDs, excited by the central star. We find that spiral shocks can at least contribute to, if not dominate the angular momentum transport and energy dissipation in CPDs. Meanwhile, dissipation and heating by spiral shocks have a positive feedback on shock-driven accretion itself. As the disk is heated up by spiral shocks, the shocks become more open, leading to more efficient angular momentum transport. This shock driven accretion is, on the other hand, unsteady on a timescale of months/years due to production and destruction of vortices in disks. After being averaged over time, a quasi-steady accretion is reached from the planet's Hill radius all the way to the planet surface, and the disk $\alpha$-coefficient characterizing angular momentum transport due to spiral shocks is $\sim$0.001-0.02. The disk surface density ranges from 10 to 1000 g cm$^{-2}$ in our simulations, which is at least 3 orders of magnitude smaller than the "minimum mass sub-nebula" model used to study satellite formation; instead it is more consistent with the "gas-starved" satellite formation model. Finally, we calculate the millimeter flux emitted by CPDs at ALMA and EVLA wavelength bands and predict the flux for several recently discovered CPD candidates, which suggests that ALMA is capable of discovering these accreting CPDs.

The global evolution of protoplanetary disks (PPDs) has recently been shown to be largely controlled by the amount of poloidal magnetic flux threading the disk, which is further controlled by the poorly understood process of magnetic flux transport. In weakly ionized gas as in PPDs, magnetic flux is largely frozen in the electron fluid, except when resistivity is large. When the disk is largely laminar, we show that the relative drift between the electrons and ions (the Hall-drift), and the ions and neutral fluids (ambipolar-drift) can play a dominant role on the transport of magnetic flux. Using two-dimensional simulations that incorporate the Hall effect and ambipolar diffusion (AD) with prescribed diffusivities, we show that when large-scale poloidal field is aligned with disk rotation, the Hall effect rapidly drags magnetic flux inward at the midplane region, while it slowly pushes flux outward above/below the midplane. This leads to a highly radially elongated field configuration as a global manifestation of the Hall-shear instability. This field configuration further promotes rapid outward flux transport by AD at the midplane, leading to instability saturation. In quasi-steady state, magnetic flux is transported outward at approximately the same rate at all heights, and the rate is comparable to the Hall-free case. For anti-aligned field polarity, the Hall effect consistently transports magnetic flux outward, leading to a largely vertical field configuration in the midplane region. The field lines in the upper layer first bend radially inward and then outward to launch a disk wind. Overall, the net rate of outward flux transport is about twice faster than the aligned case. In addition, the rate of flux transport increases with increasing disk magnetization. The absolute rate of transport is sensitive to disk microphysics which remains to be explored in future studies.

We present results from the first 3D kinetic numerical simulation of magnetorotational turbulence and dynamo, using the local shearing-box model of a collisionless accretion disc. The kinetic magnetorotational instability grows from a subthermal magnetic field having zero net flux over the computational domain to generate self-sustained turbulence and outward angular-momentum transport. Significant Maxwell and Reynolds stresses are accompanied by comparable viscous stresses produced by field-aligned ion pressure anisotropy, which is regulated primarily by the mirror and ion-cyclotron instabilities through particle trapping and pitch-angle scattering. The latter endow the plasma with an effective viscosity that is biased with respect to the magnetic-field direction and spatio-temporally variable. Energy spectra suggest an Alfv\'en-wave cascade at large scales and a kinetic-Alfv\'en-wave cascade at small scales, with strong small-scale density fluctuations and weak non-axisymmetric density waves. Ions undergo non-thermal particle acceleration, their distribution accurately described by a kappa distribution. These results have implications for the properties of low-collisionality accretion flows, such as that near the black hole at the Galactic center.

A comprehensive study on the atmospheric neutrino flux in the energy region from sub-GeV up to several TeV using the Super-Kamiokande water Cherenkov detector is presented in this paper. The energy and azimuthal spectra of the atmospheric ${\nu}_e+{\bar{\nu}}_e$ and ${\nu}_{\mu}+{\bar{\nu}}_{\mu}$ fluxes are measured. The energy spectra are obtained using an iterative unfolding method by combining various event topologies with differing energy responses. The azimuthal spectra depending on energy and zenith angle, and their modulation by geomagnetic effects, are also studied. A predicted east-west asymmetry is observed in both the ${\nu}_e$ and ${\nu}_{\mu}$ samples at 8.0 {\sigma} and 6.0 {\sigma} significance, respectively, and an indication that the asymmetry dipole angle changes depending on the zenith angle was seen at the 2.2 {\sigma} level. The measured energy and azimuthal spectra are consistent with the current flux models within the estimated systematic uncertainties. A study of the long-term correlation between the atmospheric neutrino flux and the solar magnetic activity cycle is also performed, and a weak indication of a correlation was seen at the 1.1 {\sigma} level, using SK I-IV data spanning a 20 year period. For particularly strong solar activity periods known as Forbush decreases, no theoretical prediction is available, but a deviation below the typical neutrino event rate is seen at the 2.4 {\sigma} level.

We report the results from a search in Super-Kamiokande for neutrino signals coincident with the first detected gravitational wave events, GW150914 and GW151226, using a neutrino energy range from 3.5 MeV to 100 PeV. We searched for coincident neutrino events within a time window of $\pm$500 seconds around the gravitational wave detection time. Four neutrino candidates are found for GW150914 and no candidates are found for GW151226. The remaining neutrino candidates are consistent with the expected background events. We calculated the 90\% confidence level upper limits on the combined neutrino fluence for both gravitational wave events, which depends on event energy and topologies. Considering the upward going muon data set (1.6 GeV - 100 PeV) the neutrino fluence limit for each gravitational wave event is 14 - 37 (19 - 50) cm$^{-2}$ for muon neutrinos (muon antineutrinos), depending on the zenith angle of the event. In the other data sets, the combined fluence limits for both gravitational wave events range from 2.4$\times 10^{4}$ to 7.0$\times 10^{9}$ cm$^{-2}$.

L-band spectroscopy is a powerful probe of cool low-gravity atmospheres: The P, Q, and R branch fundamental transitions of methane near 3.3 $\mu$m provide a sensitive probe of carbon chemistry; cloud thickness modifies the spectral slope across the band; and H$_{3}^{+}$ opacity can be used to detect aurorae. Many directly imaged gas-giant companions to nearby young stars exhibit L-band fluxes distinct from the field population of brown dwarfs at the same effective temperature. Here we describe commissioning the L-band spectroscopic mode of Clio2, the 1-5 $\mu$m instrument behind the Magellan adaptive-optics system. We use this system to measure L-band spectra of directly imaged companions. Our spectra are generally consistent with the parameters derived from previous near-infrared spectra for these late M to early L type objects. Therefore, deviations from the field sequence are constrained to occur below 1500 K. This range includes the L-T transition for field objects and suggests that observed discrepancies are due to differences in cloud structure and CO/CH$_{4}$ chemistry.

Isolated pulsars are rotating neutron stars with accurately measured angular velocities $\Omega$, and their time derivatives that show unambiguously that the pulsars are slowing down. The commonly accepted view is that it arises through emission of magnetic dipole radiation (MDR) from a rotating magnetized body. The calculated energy loss by a rotating pulsar with a constant moment of inertia is assumed proportional to a model dependent power of $\Omega$. This relation leads to the power law $\dot{\Omega}$ = -K $\Omega^{\rm n}$ where $n$ is called the braking index. The MDR model predicts $n$ exactly equal to 3. Selected observations of isolated pulsars provide rather precise values of $n$, individually accurate to a few percent or better, in the range 1$ <$ n $ < $ 2.8, which is consistently less than the predictions of the MDR model. In spite of an extensive investigation of various modifications of the MDR model, no satisfactory explanation of observation has been found yet. The aim of this work is to determine the deviation of the value of $n$ from the canonical $n = 3$ for a star with a frequency dependent moment of inertia in the region of frequencies from zero to the Kepler velocity (onset of mass shedding by a rotating deformed star), in the macroscopic MDR model. For the first time, we use microscopic realistic equations of state (EoS) of the star to determine its behavior and structure. In addition, we examine the effects of the baryonic mass M$_{\rm B}$ of the star, and possible core superfluidity, on the value of the braking index within the MDR model.

The magnetic dipole radiation (MDR) model is currently the best approach we have to explain pulsar radiation. However a most characteristic parameter of the observed radiation, the braking index n$_{\rm obs}$ shows deviations for all the eight best studied isolated pulsars, from the simple model prediction n$_{\rm dip}$ = 3. The index depends upon the rotational frequency and its first and second time derivatives, but also on the assumption of that the magnetic dipole moment and inclination angle, and the moment of inertia of the pulsar are constant in time. In a recent paper [Phys. Rev. D 91, 063007 (2015)] we showed conclusively that changes in the moment of inertia with frequency alone, cannot explain the observed braking indices. Possible observational evidence for the magnetic dipole moment migrating away from the rotational axis at a rate $\dot\alpha$ $\sim$ 0.6$^{\circ}$ per 100 years over the life time of the Crab pulsar has been recently suggested by Lyne et al. In this paper, we explore the MDR model with constant moment of inertia and magnetic dipole moment but variable inclination angle $\alpha$. We first discuss the effect of the variation of $\alpha$ on the observed braking indices and show they all can be understood. However, no explanation for the origin of the change in $\alpha$ is provided. After discussion of the possible source(s) of magnetism in pulsars we propose a simple mechanism for the change in $\alpha$ based on a toy model in which the magnetic structure in pulsars consists of two interacting dipoles. We show that such a system can explain the Crab observation and the measured braking indices.

The Large Binocular Telescope Interferometer uses a near-infrared camera to measure the optical path length variations between the two AO-corrected apertures and provide high-angular resolution observations for all its science channels (1.5-13 $\mu$m). There is however a wavelength dependent component to the atmospheric turbulence, which can introduce optical path length errors when observing at a wavelength different from that of the fringe sensing camera. Water vapor in particular is highly dispersive and its effect must be taken into account for high-precision infrared interferometric observations as described previously for VLTI/MIDI or the Keck Interferometer Nuller. In this paper, we describe the new sensing approach that has been developed at the LBT to measure and monitor the optical path length fluctuations due to dry air and water vapor separately. After reviewing the current performance of the system for dry air seeing compensation, we present simultaneous H-, K-, and N-band observations that illustrate the feasibility of our feedforward approach to stabilize the path length fluctuations seen by the LBTI nuller.

We present a new general relativistic magnetohydrodynamics (GRMHD) code integrated into the Athena++ framework. Improving upon the techniques used in most GRMHD codes, ours allows the use of advanced, less diffusive Riemann solvers, in particular HLLC and HLLD. We also employ a staggered-mesh constrained transport algorithm suited for curvilinear coordinate systems in order to maintain the divergence-free constraint of the magnetic field. Our code is designed to work with arbitrary stationary spacetimes in one, two, or three dimensions, and we demonstrate its reliability in a number of tests. We also report on its promising performance and scalability.

To understand the conditions under which dense, molecular gas is able to form within a galaxy, we post-process a series of three-dimensional galactic-disk-scale simulations with ray-tracing based radiative transfer and chemical network integration to compute the equilibrium chemical and thermal state of the gas. In performing these simulations we vary a number of parameters, such as the ISRF strength, vertical scale height of stellar sources, cosmic ray flux, to gauge the sensitivity of our results to these variations. Self-shielding permits significant molecular hydrogen (H2) abundances in dense filaments around the disk midplane, accounting for approximately ~10-15% of the total gas mass. Significant CO fractions only form in the densest, n>~10^3 cm^-3, gas where a combination of dust, H2, and self-shielding attenuate the FUV background. We additionally compare these ray-tracing based solutions to photochemistry with complementary models where photo-shielding is accounted for with locally computed prescriptions. With some exceptions, these local models for the radiative shielding length perform reasonably well at reproducing the distribution and amount of molecular gas as compared with a detailed, global ray tracing calculation. Specifically, an approach based on the Jeans Length with a T=40K temperature cap performs the best in regards to a number of different quantitative measures based on the H2 and CO abundances.

Variable emission from Sgr~A*, the luminous counterpart to the super-massive black hole at the center of our Galaxy, arises from the innermost portions of the accretion flow. Better characterization of the variability is important for constraining models of the low-luminosity accretion mode powering Sgr~A*, and could further our ability to use variable emission as a probe of the strong gravitational potential in the vicinity of the $4\times10^{6}\mathrm{M}_{\odot}$ black hole. We use the \textit{Herschel} Spectral and Photometric Imaging Receiver (SPIRE) to monitor Sgr~A* at wavelengths that are difficult or impossible to observe from the ground. We find highly significant variations at 0.25, 0.35, and 0.5 mm, with temporal structure that is highly correlated across these wavelengths. While the variations correspond to $<$1% changes in the total intensity in the \textit{Herschel} beam containing Sgr~A*, comparison to independent, simultaneous observations at 0.85 mm strongly supports the reality of the variations. The lowest point in the light curves, $\sim$0.5 Jy below the time-averaged flux density, places a lower bound on the emission of Sgr~A* at 0.25 mm, the first such constraint on the THz portion of the SED. The variability on few hour timescales in the SPIRE light curves is similar to that seen in historical 1.3 mm data, where the longest time series is available, but the distribution of variations in the sub-mm do not show a tail of large-amplitude variations seen at 1.3 mm. Simultaneous X-ray photometry from XMM-Newton shows no significant variation within our observing period, which may explain the lack of very large variations if X-ray and submillimeter flares are correlated.

We present a real-time supernova neutrino burst monitor at Super-Kamiokande (SK). Detecting supernova explosions by neutrinos in real time is crucial for giving a clear picture of the explosion mechanism. Since the neutrinos are expected to come earlier than light, a fast broadcasting of the detection may give astronomers a chance to make electromagnetic radiation observations of the explosions right at the onset. The role of the monitor includes a fast announcement of the neutrino burst detection to the world and a determination of the supernova direction. We present the online neutrino burst detection system and studies of the direction determination accuracy based on simulations at SK.

We present results from the first global 3D MHD simulations of accretion disks in Cataclysmic Variable (CV) systems in order to investigate the relative importance of angular momentum transport via turbulence driven by the magnetorotational instability (MRI) compared to that driven by spiral shock waves. Remarkably, we find that even with vigorous MRI turbulence, spiral shocks are an important component to the overall angular momentum budget, at least when temperatures in the disk are high (so that Mach numbers are low). In order to understand the excitation, propagation, and damping of spiral density waves in our simulations more carefully, we perform a series of 2D global hydrodynamical simulations with various equation of states and both with and without mass inflow via the Lagrangian point (L1). Compared with previous similar studies, we find the following new results. 1) Linear wave dispersion relation fits the pitch angles of spiral density waves very well. 2) We demonstrate explicitly that mass accretion is driven by the deposition of negative angular momentum carried by the waves when they dissipate in shocks. 3) Using Reynolds stress scaled by gas pressure to represent the effective angular momentum transport rate alpha_{eff} is not accurate when mass accretion is driven by non-axisymmetric shocks. 4) Using the mass accretion rate measured in our simulations to directly measure alpha defined in standard thin-disk theory, we find 0.02 < alpha_{eff} < 0.05 for CV disks, consistent with observed values in quiescent states of dwarf novae (DNe). In this regime the disk may be too cool and neutral for the MRI to operate and spiral shocks are a possible accretion mechanism. However, we caution that our simulations use unrealistically low Mach numbers in this regime, and therefore future models with more realistic thermodynamics and non-ideal MHD are warranted.

The dynamics of solid bodies in protoplanetary disks are subject to the properties of any underlying gas turbulence. Turbulence driven by disk self-gravity shows features distinct from those driven by the magnetorotational instability (MRI). We study the dynamics of solids in gravito-turbulent disks with two-dimensional (in the disk plane), hybrid (particle and gas) simulations. Gravito-turbulent disks can exhibit stronger gravitational stirring than MRI-active disks, resulting in greater radial diffusion and larger eccentricities and relative speeds for large particles (those with dimensionless stopping times $t_{stop} \Omega > 1$, where $\Omega$ is the orbital frequency). The agglomeration of large particles into planetesimals by pairwise collisions is therefore disfavored in gravito-turbulent disks. However, the relative speeds of intermediate-size particles $t_{stop} \Omega \sim 1$ are significantly reduced as such particles are collected by gas drag and gas gravity into coherent filament-like structures with densities high enough to trigger gravitational collapse. First-generation planetesimals may form via gravitational instability of dust in marginally gravitationally unstable gas disks.

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

Accretion disks around supermassive black holes have regions where the Rosseland mean opacity can be much larger than the electron scattering opacity primarily due to the large number of bound-bound transitions in iron. We study the effects of this iron opacity "bump" on the thermal stability and vertical structure of radiation pressure dominated accretion disks, utilizing three dimensional radiation magneto-hydrodynamic simulations in the local shearing box approximation. The simulations self-consistently calculate the heating due to MHD turbulence caused by magneto-rotational instability and radiative cooling by using the radiative transfer module based on a variable Eddington tensor in Athena. For a $5\times 10^8$ solar mass black hole with $\sim 3\%$ of the Eddington luminosity, a model including the iron opacity bump maintains its structure for more than $10$ thermal times without showing significant signs of thermal runaway. In contrast, if only electron scattering and free-free opacity are included as in the standard thin disk model, the disk collapses on the thermal time scale. The difference is caused by a combination of 1) an anti-correlation between the total optical depth and the temperature, and 2) enhanced vertical advective energy transport. These results suggest that the iron opacity bump may have a strong impact on the stability and structure of AGN accretion disks, and may contribute to a dependence of AGN properties on metallicity. Since this opacity is relevant primarily in UV emitting regions of the flow, it may help to explain discrepancies between observation and theory that are unique to AGNs.

Recently, Gauza et al. (2015) reported the discovery of a companion to the late M-dwarf, VHS J125601.92-125723.9 (VHS 1256-1257). The companion's absolute photometry suggests its mass and atmosphere are similar to the HR 8799 planets. However, as a wide companion to a late-type star, it is more accessible to spectroscopic characterization. We discovered that the primary of this system is an equal-magnitude binary. For an age $\sim300$ Myr the A and B components each have a mass of $64.6^{+0.8}_{-2.0}~M_{\mathrm{Jup}}$, and the b component has a mass of $11.2^{+9.7}_{-1.8}$, making VHS 1256-1257 only the third brown dwarf triple system. There exists some tension between the spectrophotometric distance of $17.2\pm2.6$ pc and the parallax distance of $12.7\pm1.0$ pc. At 12.7 pc VHS1256-1257 A and B would be the faintest known M7.5 objects, and are even faint outliers among M8 types. If the larger spectrophotmetric distance is more accurate than the parallax, then the mass of each component increases. In particular, the mass of the b component increases well above the deuterium burning limit to $\sim35~M_{\mathrm{Jup}}$ and the mass of each binary component increases to $73^{+20}_{-17}~M_{\mathrm{Jup}}$. At 17.1 pc, the UVW kinematics of the system are consistent with membership in the AB~Dor moving group. The architecture of the system resembles a hierarchical stellar multiple suggesting it formed via an extension of the star-formation process to low masses. Continued astrometric monitoring will resolve this distance uncertainty and will provide dynamical masses for a new benchmark system.

Previous studies of the nonlinear regime of the magnetorotational instability in one particular type of shearing box model -- unstratified with no net magnetic flux -- find that without explicit dissipation (viscosity and resistivity) the saturation amplitude decreases with increasing numerical resolution. We show that this result is strongly dependent on the vertical aspect ratio of the computational domain $L_z/L_x$. When $L_z/L_x\lesssim 1$, we recover previous results. However, when the vertical domain is extended $L_z/L_x \gtrsim 2.5$, we find the saturation level of the stress is greatly increased (giving a ratio of stress to pressure $\alpha \gtrsim 0.1$), and moreover the results are independent of numerical resolution. Consistent with previous results, we find that saturation of the MRI in this regime is controlled by a cyclic dynamo which generates patches of strong toroidal field that switches sign on scales of $L_x$ in the vertical direction. We speculate that when $L_z/L_x\lesssim 1$, the dynamo is inhibited by the small size of the vertical domain, leading to the puzzling dependence of saturation amplitude on resolution. We show that previous toy models developed to explain the MRI dynamo are consistent with our results, and that the cyclic pattern of toroidal fields observed in stratified shearing box simulations (leading to the so-called butterfly diagram) may also be related. In tall boxes the saturation amplitude is insensitive to whether or not explicit dissipation is included in the calculations, at least for large magnetic Reynolds and Prandtl number. Finally, we show MRI turbulence in tall domains has a smaller critical $\rm{Pm}_c$, and an extended lifetime compared to $L_z/L_x\lesssim 1$ boxes.

As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the first generation of directly-imaged exoplanets were all L-type. Recently, Kuzuhara et al. (2013) announced the discovery of GJ 504 b, the first T dwarf exoplanet. GJ 504 b provides a unique opportunity to study the atmosphere of a new type of exoplanet with a ~500 K temperature that bridges the gap between the first directly imaged planets (~1000 K) and our own Solar System's Jupiter (~130 K). We observed GJ 504 b in three narrow L-band filters (3.71, 3.88, and 4.00 microns), spanning the red end of the broad methane fundamental absorption feature (3.3 microns) as part of the LEECH exoplanet imaging survey. By comparing our new photometry and literature photometry to a grid of custom model atmospheres, we were able to fit GJ 504 b's unusual spectral energy distribution for the first time. We find that GJ 504 b is well-fit by models with the following parameters: T_eff=544+/-10 K, g<600 m/s^2, [M/H]=0.60+/-0.12, cloud opacity parameter of f_sed=2-5, R=0.96+/-0.07 R_Jup, and log(L)=-6.13+/-0.03 L_Sun, implying a hot start mass of 3-30 M_jup for a conservative age range of 0.1-6.5 Gyr. Of particular interest, our model fits suggest that GJ 504 b has a super-stellar metallicity. Since planet formation can create objects with non-stellar metallicities, while binary star formation cannot, this result suggests that GJ 504 b formed like a planet, not like a binary companion.

Disk accretion at high rate onto a white dwarf or a neutron star has been suggested to result in the formation of a spreading layer (SL) - a belt-like structure on the object's surface, in which the accreted matter steadily spreads in the poleward (meridional) direction while spinning down. To assess its basic characteristics we perform two-dimensional hydrodynamic simulations of supersonic SLs in the relevant morphology with a simple prescription for cooling. We demonstrate that supersonic shear naturally present at the base of the SL inevitably drives sonic instability that gives rise to large scale acoustic modes governing the evolution of the SL. These modes dominate the transport of momentum and energy, which is intrinsically global and cannot be characterized via some form of local effective viscosity (e.g. $\alpha$-viscosity). The global nature of the wave-driven transport should have important implications for triggering Type I X-ray bursts in low mass X-ray binaries. The nonlinear evolution of waves into a system of shocks drives effective re-arrangement (sensitively depending on thermodynamical properties of the flow) and deceleration of the SL, which ultimately becomes transonic and susceptible to regular Kelvin-Helmholtz instability. We interpret this evolution in terms of the global structure of the SL and suggest that mixing of the SL material with the underlying stellar fluid should become effective only at intermediate latitudes on the accreting object's surface, where the flow has decelerated appreciably. In the near-equatorial regions the transport is dominated by acoustic waves and mixing is less efficient. We speculate that this latitudinal non-uniformity of mixing in accreting white dwarfs may be linked to the observed bipolar morphology of classical novae ejecta.

We explore the equation of state for nuclear matter in the quark-meson coupling model, including full Fock terms. The comparison with phenomenological constraints can be used to restrict the few additional parameters appearing in the Fock terms which are not present at Hartree level. Because the model is based upon the in-medium modification of the quark structure of the bound hadrons, it can be applied without additional parameters to include hyperons and to calculate the equation of state of dense matter in beta-equilibrium. This leads naturally to a study of the properties of neutron stars, including their maximum mass, their radii and density profiles.

Direct imaging observations have revealed spiral structures in protoplanetary disks. Previous studies have suggested that planet-induced spiral arms cannot explain some of these spiral patterns, due to the large pitch angle and high contrast of the spiral arms in observations. We have carried out three dimensional (3-D) hydrodynamical simulations to study spiral wakes/shocks excited by young planets. We find that, in contrast with linear theory, the pitch angle of spiral arms does depend on the planet mass, which can be explained by the non-linear density wave theory. A secondary (or even a tertiary) spiral arm, especially for inner arms, is also excited by a massive planet. With a more massive planet in the disk, the excited spiral arms have larger pitch angle and the separation between the primary and secondary arms in the azimuthal direction is also larger. We also find that although the arms in the outer disk do not exhibit much vertical motion, the inner arms have significant vertical motion, which boosts the density perturbation at the disk atmosphere. Combining hydrodynamical models with Monte-Carlo radiative transfer calculations, we find that the inner spiral arms are considerably more prominent in synthetic near-IR images using full 3-D hydrodynamical models than images based on 2-D models assuming vertical hydrostatic equilibrium, indicating the need to model observations with full 3-D hydrodynamics. Overall, companion-induced spiral arms not only pinpoint the companion's position but also provide three independent ways (pitch angle, separation between two arms, and contrast of arms) to constrain the companion's mass.

The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack of a reference solution has led various authors to assert the accuracy of their simulations based on ad-hoc proxies, e.g., the existence of small-scale structures. This paper proposes well-posed Kelvin-Helmholtz problems with smooth initial conditions and explicit diffusion. We show that in many cases numerical errors/noise can seed spurious small-scale structure in Kelvin-Helmholtz problems. We demonstrate convergence to a reference solution using both Athena, a Godunov code, and Dedalus, a pseudo-spectral code. Problems with constant initial density throughout the domain are relatively straightforward for both codes. However, problems with an initial density jump (which are the norm in astrophysical systems) exhibit rich behavior and are more computationally challenging. In the latter case, Athena simulations are prone to an instability of the inner rolled-up vortex; this instability is seeded by grid-scale errors introduced by the algorithm, and disappears as resolution increases. Both Athena and Dedalus exhibit late-time chaos. Inviscid simulations are riddled with extremely vigorous secondary instabilities which induce more mixing than simulations with explicit diffusion. Our results highlight the importance of running well-posed test problems with demonstrated convergence to a reference solution. To facilitate future comparisons, we include the resolved, converged solutions to the Kelvin-Helmholtz problems in this paper in machine-readable form.

The Large Binocular Telescope Interferometer (LBTI) is a strategic instrument of the LBT designed for high-sensitivity, high-contrast, and high-resolution infrared (1.5-13 $\mu$m) imaging of nearby planetary systems. To carry out a wide range of high-spatial resolution observations, it can combine the two AO-corrected 8.4-m apertures of the LBT in various ways including direct (non-interferometric) imaging, coronagraphy (APP and AGPM), Fizeau imaging, non-redundant aperture masking, and nulling interferometry. It also has broadband, narrowband, and spectrally dispersed capabilities. In this paper, we review the performance of these modes in terms of exoplanet science capabilities and describe recent instrumental milestones such as first-light Fizeau images (with the angular resolution of an equivalent 22.8-m telescope) and deep interferometric nulling observations.

Integral field spectrographs are an important technology for exoplanet imaging, due to their ability to take spectra in a high-contrast environment, and improve planet detection sensitivity through spectral differential imaging. ALES is the first integral field spectrograph capable of imaging exoplanets from 3-5$\mu$m, and will extend our ability to characterize self-luminous exoplanets into a wavelength range where they peak in brightness. ALES is installed inside LBTI/LMIRcam on the Large Binocular Telescope, taking advantage of existing AO systems, camera optics, and a HAWAII-2RG detector. The new optics that comprise ALES are a Keplerian magnifier, a silicon lenslet array with diffraction suppressing pinholes, a direct vision prism, and calibration optics. All of these components are installed in filter wheels making ALES a completely modular design. ALES saw first light at the LBT in June 2015.

Using 3D global hydro simulations coupled with radiative transfer calculations, we study the appearance of density waves induced by giant planets in direct imaging observations at near infrared wavelengths. We find that a 6 MJ planet in a typical disk around a 1 M_sun star can produce prominent and detectable spiral arms both interior and exterior to its orbit. The inner arms have (1) two well separated arms in roughly m=2 symmetry, (2) exhibit ~10-15 degrees pitch angles, (3) ~180-270 degrees extension in the azimuthal direction, and (4) ~150% surface brightness enhancement, all broadly consistent with observed spiral arms in the SAO 206462 and MWC 758 systems. The outer arms cannot explain observations as they are too tightly wound given typical disk scale height. We confirm previous results that the outer density waves excited by a 1 MJ planet exhibit low contrast in the IR and are practically not detectable. We also find that 3D effects of the waves are important. Compared to isothermal models, density waves in adiabatic disks exhibit weaker contrast in surface density but stronger contrast in scattered light images, due to a more pronounced vertical structure in the former caused by shock heating. To drive observed pairs of arms with an external companion on a circular orbit, a massive planet, possibly a brown dwarf, is needed at around [r~0.7", PA~10 degrees] (position angle PA from north to east) in SAO 206462 and [r~0.6", PA~10 degrees] in MWC 758. Their existence may be confirmed by direct imaging planet searches.

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

We have studied the impact of dust feedback on the survival and structure of vortices in protoplanetary discs using 2-D shearing box simulations with Lagrangian dust particles. We consider dust with a variety of sizes (stopping time $t_s = 10^{-2}\Omega^{-1} - 10^{2}\Omega^{-1}$, from fully coupled with the gas to the decoupling limit. We find that a vortex is destroyed by dust feedback when the total dust-to-gas mass ratio within the vortex is larger than 30-50%, independent of the dust size. The dust distribution can still be asymmetric in some cases after the vortex has been destroyed. With smaller amounts of dust, a vortex can survive for at least 100 orbits, and the maximum dust surface density within the vortex can be more than 100 times larger than the gas surface density, potentially facilitating planetesimal formation. On the other hand, in these stable vortices, small ($t_s < \Omega^{-1}$) and large ($t_s > \Omega^{-1}$) dust grains concentrate differently and affect the gas dynamics in different ways. The distribution of large dust is more elongated than that of small dust. Large dust ($t_s > \Omega^{-1}$) concentrates in the centre of the vortex and feedback leads to turn-over in vorticity towards the centre, forming a quiescent region within an anticyclonic vortex. Such a turn-over is absent if the vortex is loaded with small grains. We demonstrate that, in protoplanetary discs where both large and small dust grains are present and under the right condition, the concentration of large dust towards the vortex centre can lead to a quiescent centre, repelling the small dust and forming a small dust ring around the vortex centre. Such anticorrelations between small and large dust within vortices may explain the discrepancy between ALMA and near-IR scattered light observations in the asymmetric region of transitional discs.

Isolated pulsars are rotating neutron stars with accurately measured angular velocities $\Omega$, and their time derivatives which show unambiguously that the pulsars are slowing down. Although the exact mechanism of the spin-down is a question of debate in detail, the commonly accepted view is that it arises through emission of magnetic dipole radiation (MDR) from a rotating magnetized body. Other processes, including the emission of gravitational radiation, and of relativistic particles (pulsar wind), are also being considered. The calculated energy loss by a rotating pulsar with a constant moment of inertia is assumed proportional to a model dependent power of $\Omega$. This relation leads to the power law $\dot{\Omega}$ = -K $\Omega^{\rm n}$ where $n$ is called the braking index. The MDR model predicts $n$ exactly equal to 3. Selected observations of isolated pulsars provide rather precise values of $n$, individually accurate to a few percent or better, in the range 1$ <$ n $ < $ 2.8, which is consistently less than the predictions of the MDR model. In spite of an extensive investigation of various modifications of the MDR model, no satisfactory explanation of observation has been found as yet. We employ four physically realistic equations of state, and two computational codes, to model the dynamical effects of rotation on the braking index in the MDR model. In addition to this we simulate an effect on moment of inertia where we assume a certain amount of superfluid matter has manifest between the crust and core region of the star thus essentially eliminating momentum transfer between the two regions. We find that the effects of rotation on braking index are significant at high frequencies, but have little effect at frequencies consistent with the most accurately measured pulsars to-date.

Cold fronts -- contact discontinuities in the intracluster medium (ICM) of galaxy clusters -- should be disrupted by Kelvin-Helmholtz (K-H) instabilities due to the associated shear velocity. However, many observed cold fronts appear stable. This opens the possibility to place constraints on microphysical mechanisms that stabilize them, such as the ICM viscosity and/or magnetic fields. We performed exploratory high-resolution simulations of cold fronts arising from subsonic gas sloshing in cluster cores using the grid-based Athena MHD code, comparing the effects of isotropic Spitzer and anisotropic Braginskii viscosity (expected in a magnetized plasma). Magnetized simulations with full Braginskii viscosity or isotropic Spitzer viscosity reduced by a factor f ~ 0.1 are both in qualitative agreement with observations in terms of suppressing K-H instabilities. The RMS velocity of turbulence within the sloshing region is only modestly reduced by Braginskii viscosity. We also performed unmagnetized simulations with and without viscosity and find that magnetic fields have a substantial effect on the appearance of the cold fronts, even if the initial field is weak and the viscosity is the same. This suggests that determining the dominant suppression mechanism of a given cold front from X-ray observations (e.g. viscosity or magnetic fields) by comparison with simulations is not straightforward. Finally, we performed simulations including anisotropic thermal conduction, and find that including Braginskii viscosity in these simulations does not significant affect the evolution of cold fronts; they are rapidly smeared out by thermal conduction, as in the inviscid case.

We study dust transport in turbulent protoplanetary disks using three-dimensional global unstratified magnetohydrodynamic (MHD) simulations including Lagrangian dust particles. The turbulence is driven by the magnetorotational instability (MRI) with either ideal or non-ideal MHD that includes ambipolar diffusion (AD). In ideal MHD simulations, the surface density evolution (except for dust that drifts fastest), turbulent diffusion, and vertical scale height of dust can all be reproduced by simple one-dimensoinal and/or analytical models. However, in AD dominated simulations which simulate protoplanetary disks beyond 10s of AU, the vertical scale height of dust is larger than previously predicted. To understand this anomaly in more detail, we carry out both unstratified and stratified local shearing box simulations with Lagrangian particles, and find that turbulence in AD dominated disks has very different properties (e.g., temporal autocorrelation functions and power spectra) than turbulence in ideal MHD disks, which leads to quite different particle diffusion efficiency. For example, MRI turbulence with AD has a longer correlation time for the vertical velocity, which causes significant vertical particle diffusion and large dust scale height. In ideal MHD the Schmidt numbers ($Sc$) for radial and vertical turbulent diffusion are $Sc_{r}\sim 1$ and $Sc_{z}\gtrsim 3$, but in the AD dominated regime both $Sc_{r}$ and $Sc_{z}$ are $\lesssim 1$. Particle concentration in pressure bumps induced by MRI turbulence has also been studied. Since non-ideal MHD effects dominate most regions in protoplanetary disks, our study suggests that modeling dust transport in turbulence driven by MRI with non-ideal MHD effects is important for understanding dust transport in realistic protoplanetary disks.

We present multiple epochs of near-IR spectroscopy for a sample of 25 young stars, including T Tauri, Herbig Ae/Be, and FU Ori objects. Using the FSPEC instrument on the Bok 90-inch telescope, we obtained K-band spectra of the BrGamma transition of hydrogen, with a resolution of ~3500. Epochs were taken over a span of >1 year, sampling time-spacings of roughly one day, one month, and one year. The majority of our targets show BrGamma emission, and in some cases these are the first published detections. Time-variability is seen in approximately half of the targets showing BrGamma emission. We compare the observed variability with expectations for rotationally-modulated accretion onto the central stars and time-variable continuum emission or extinction from matter in the inner disk. Our observations are not entirely consistent with models of rotationally-modulated magnetospheric accretion. Further monitoring, over a larger number of epochs, will facilitate more quantitative constraints on variability timescales and amplitudes, and a more conclusive comparison with theoretical models.

Relativistic mean-field (RMF) models have been widely used in the study of many hadronic frameworks because of several important aspects not always present in nonrelativistic models, such as intrinsic Lorentz covariance, automatic inclusion of spin, appropriate saturation mechanism for nuclear matter, causality and, therefore, no problems related to superluminal speed of sound. With the aim of identifying the models which best satisfy well known properties of nuclear matter, we have analyzed $263$ parameterizations of seven different types of RMF models under three different sets of constraints related to symmetric nuclear matter, pure neutron matter, symmetry energy, and its derivatives. One of these (SET1) is formed of the same constraints used in a recent work [M. Dutra et al., Phys. Rev. C 85, 035201 (2012)] in which we analyzed $240$ Skyrme parameterizations. The results pointed to $2$ models consistent with all constraints. By using another set of constraints, namely, SET2a, formed by the updated versions of the previous one, we found $4$ models approved simultaneously. Finally, in the third set, named SET2b, in which the values of the constraints are more restrictive, we found $3$ consistent models. Another interesting feature of our analysis is that the results change dramatically if we do not consider the constraint regarding the volume part of the isospin incompressibility ($K_{\tau,\rm v}$). In this case, we have $35$ approved models in SET2a and $30$ in SET2b.

We study super-Eddington accretion flows onto black holes using a global three dimensional radiation magneto-hydrodynamical simulation. We solve the time dependent radiative transfer equation for the specific intensities to accurately calculate the angular distribution of the emitted radiation. Turbulence generated by the magneto-rotational instability provides self-consistent angular momentum transfer. The simulation reaches inflow equilibrium with an accretion rate ~220L_edd/c^2 and forms a radiation driven outflow along the rotation axis. The mechanical energy flux carried by the outflow is ~20% of the radiative energy flux. The total mass flux lost in the outflow is about 29% of the net accretion rate. The radiative luminosity of this flow is ~10L_edd. This yields a radiative efficiency ~4.5%, which is comparable to the value in a standard thin disk model. In our simulation, vertical advection of radiation caused by magnetic buoyancy transports energy faster than photon diffusion, allowing a significant fraction of the photons to escape from the surface of the disk before being advected into the black hole. We contrast our results with the lower radiative efficiencies inferred in most models, such as the slim disk model, which neglect vertical advection. Our inferred radiative efficiencies also exceed published results from previous global numerical simulations, which did not attribute a significant role to vertical advection. We briefly discuss the implications for the growth of supermassive black holes in the early universe and describe how these results provided a basis for explaining the spectrum and population statistics of ultraluminous X-ray sources.

Accretion disks are likely threaded by external vertical magnetic flux, which enhances the level of turbulence via the magnetorotational instability (MRI). Using shearing-box simulations, we find that such external magnetic flux also strongly enhances the amplitude of banded radial density variations known as zonal flows. Moreover, we report that vertical magnetic flux is strongly concentrated toward low-density regions of the zonal flow. Mean vertical magnetic field can be more than doubled in low-density regions, and reduced to nearly zero in high density regions in some cases. In ideal MHD, the scale on which magnetic flux concentrates can reach a few disk scale heights. In the non-ideal MHD regime with strong ambipolar diffusion, magnetic flux is concentrated into thin axisymmetric shells at some enhanced level, whose size is typically less than half a scale height. We show that magnetic flux concentration is closely related to the fact that the magnetic diffusivity of the MRI turbulence is anisotropic. In addition to a conventional Ohmic-like turbulent resistivity, we find that there is a correlation between the vertical velocity and horizontal magnetic field fluctuations that produces a mean electric field that acts to anti-diffuse the vertical magnetic flux. The anisotropic turbulent diffusivity has analogies to the Hall effect, and may have important implications for magnetic flux transport in accretion disks. The physical origin of magnetic flux concentration may be related to the development of channel flows followed by magnetic reconnection, which acts to decrease the mass-to-flux ratio in localized regions. The association of enhanced zonal flows with magnetic flux concentration may lead to global pressure bumps in protoplanetary disks that helps trap dust particles and facilitates planet formation.

We study particle trapping at the edge of a gap opened by a planet in a protoplanetary disk. In particular, we explore the effects of turbulence driven by the magnetorotational instability on particle trapping, using global three-dimensional magnetohydrodynamic (MHD) simulations including Lagrangian dust particles. We study disks either in the ideal MHD limit or dominated by ambipolar diffusion (AD) that plays an essential role at the outer regions of a protoplanetary disk. With ideal MHD, strong turbulence (the equivalent viscosity parameter $\alpha\sim 10^{-2}$) in disks prevents vortex formation at the edge of the gap opened by a 9 $M_{J}$ planet, and most particles (except the particles that drift fastest) pile up at the outer gap edge almost axisymmetrically. When AD is considered, turbulence is significantly suppressed ($\alpha\lesssim 10^{-3}$), and a large vortex forms at the edge of the planet induced gap, which survives $\sim$ 1000 orbits. The vortex can efficiently trap dust particles that span 3 orders of magnitude in size within 100 planetary orbits. We have also carried out two-dimensional hydrodynamical simulations using viscosity as an approximation to MHD turbulence. These hydrodynamical simulations can reproduce vortex generation at the gap edge as seen in MHD simulations. Finally, we use our simulation results to generate synthetic images for ALMA dust continuum observations on Oph IRS 48 and HD 142527, which show good agreement with existing observations. Predictions for future ALMA cycle 2 observations have been made. We conclude that the asymmetry in ALMA observations can be explained by dust trapping vortices and the existence of vortices could be the evidence that the outer protoplanetary disks are dominated by ambipolar diffusion with $\alpha<10^{-3}$ at the disk midplane.

One process affecting gas-rich cluster galaxies is ram pressure stripping, i.e. the removal of galactic gas through direct interaction with the intracluster medium. Galactic magnetic fields may have an important impact on the stripping rate and tail structure. We run the first magnetohydrodynamic simulations of ram pressure stripping that include a galactic magnetic field, using 159 pc resolution throughout our entire domain in order to resolve mixing throughout the tail. We find very little difference in the total amount of gas removed from the unmagnetized and magnetized galaxies, although a magnetic field with a radial component will initially accelerate stripped gas more quickly. In general, we find that magnetic fields in the disk lead to slower velocities in the stripped gas near the disk and faster velocities farther from the disk. We also find that magnetic fields in the galactic gas lead to larger unmixed structures in the tail. Finally, we discuss whether ram pressure stripped tails can magnetize the ICM. We find that the total magnetic energy density grows as the tail lengthens, likely through turbulence. There are microGauss-strength fields in the tail in all of our MHD runs, which survive to at least 100 kpc from the disk (the edge of our simulated region), indicating that the area-filling factor of magnetized tails in a cluster could be large.

We present new K-band spectroscopy of the UY Aur binary star system. Our data are the first to show H$_{2}$ emission in the spectrum of UY Aur A and the first to spectrally resolve the Br{\gamma} line in the spectrum of UY Aur B. We see an increase in the strength of the Br{\gamma} line in UY Aur A and a decrease in Br{\gamma} and H$_{2}$ line luminosity for UY Aur B compared to previous studies. Converting Br{\gamma} line luminosity to accretion rate, we infer that the accretion rate onto UY Aur A has increased by $2 \times 10^{-9}$ M$_{\odot}$ yr$^{-1}$ per year since a rate of zero was observed in 1994. The Br{\gamma} line strength for UY Aur B has decreased by a factor of 0.54 since 1994, but the K-band flux has increased by 0.9 mags since 1998. The veiling of UY Aur B has also increased significantly. These data evince a much more luminous disk around UY Aur B. If the lower Br{\gamma} luminosity observed in the spectrum of UY Aur B indicates an intrinsically smaller accretion rate onto the star, then UY Aur A now accretes at a higher rate than UY Aur B. However, extinction at small radii or mass pile-up in the circumstellar disk could explain decreased Br{\gamma} emission around UY Aur B even when the disk luminosity implies an increased accretion rate. In addition to our scientific results for the UY Aur system, we discuss a dedicated pipeline we have developed for the reduction of echelle-mode data from the ARIES spectrograph.

We present Keck Interferometer observations of T Tauri and Herbig Ae/Be stars with a spatial resolution of a few milliarcseconds and a spectral resolution of ~2000. Our observations span the K-band, and include the Br gamma transition of Hydrogen and the v=2-0 and v=3-1 transitions of carbon monoxide. For several targets we also present data from Keck/NIRSPEC that provide higher spectral resolution, but a seeing-limited spatial resolution, of the same spectral features. We analyze the Br gamma emission in the context of both disk and infall/outflow models, and conclude that the Br gamma emission traces gas at very small stellocentric radii, consistent with the magnetospheric scale. However some Br gamma-emitting gas also seems to be located at radii of >0.1 AU, perhaps tracing the inner regions of magnetically launched outflows. CO emission is detected from several objects, and we generate disk models that reproduce both the KI and NIRSPEC data well. We infer the CO spatial distribution to be coincident with the distribution of continuum emission in most cases. Furthermore the Br gamma emission in these objects is roughly coincident with both the CO and continuum emission. We present potential explanations for the spatial coincidence of continuum, Br gamma, and CO overtone emission, and explore the implications for the low occurrence rate of CO overtone emission in young stars. Finally, we provide additional discussion of V1685 Cyg, which is unusual among our sample in showing large differences in emitting region size and spatial position as a function of wavelength.

Hybrid-kinetic numerical simulations of firehose and mirror instabilities in a collisionless plasma are performed in which pressure anisotropy is driven as the magnetic field is changed by a persistent linear shear $S$. For a decreasing field, it is found that mostly oblique firehose fluctuations grow at ion Larmor scales and saturate with energies $\sim$$S^{1/2}$; the pressure anisotropy is pinned at the stability threshold by particle scattering off microscale fluctuations. In contrast, nonlinear mirror fluctuations are large compared to the ion Larmor scale and grow secularly in time; marginality is maintained by an increasing population of resonant particles trapped in magnetic mirrors. After one shear time, saturated order-unity magnetic mirrors are formed and particles scatter off their sharp edges. Both instabilities drive sub-ion-Larmor--scale fluctuations, which appear to be kinetic-Alfv\'{e}n-wave turbulence. Our results impact theories of momentum and heat transport in astrophysical and space plasmas, in which the stretching of a magnetic field by shear is a generic process.

The incompressibility (compression modulus) $K_{\rm 0}$ of infinite symmetric nuclear matter at saturation density has become one of the major constraints on mean-field models of nuclear many-body systems as well as of models of high density matter in astrophysical objects and heavy-ion collisions. We present a comprehensive re-analysis of recent data on GMR energies in even-even $^{\rm 112-124}$Sn and $^{\rm 106,100-116}$Cd and earlier data on 58 $\le$ A $\le$ 208 nuclei. The incompressibility of finite nuclei $K_{\rm A}$ is expressed as a leptodermous expansion with volume, surface, isospin and Coulomb coefficients $K_{\rm vol}$, $K_{\rm surf}$, $K_\tau$ and $K_{\rm coul}$. \textit{Assuming} that the volume coefficient $K_{\rm vol}$ is identified with $K_{\rm 0}$, the $K_{\rm coul}$ = -(5.2 $\pm$ 0.7) MeV and the contribution from the curvature term K$_{\rm curv}$A$^{\rm -2/3}$ in the expansion is neglected, compelling evidence is found for $K_{\rm 0}$ to be in the range 250 $ < K_{\rm 0} < $ 315 MeV, the ratio of the surface and volume coefficients $c = K_{\rm surf}/K_{\rm vol}$ to be between -2.4 and -1.6 and $K_{\rm \tau}$ between -840 and -350 MeV. We show that the generally accepted value of $K_{\rm 0}$ = (240 $\pm$ 20) MeV can be obtained from the fits provided $c \sim$ -1, as predicted by the majority of mean-field models. However, the fits are significantly improved if $c$ is allowed to vary, leading to a range of $K_{\rm 0}$, extended to higher values. A self-consistent simple (toy) model has been developed, which shows that the density dependence of the surface diffuseness of a vibrating nucleus plays a major role in determination of the ratio K$_{\rm surf}/K_{\rm vol}$ and yields predictions consistent with our findings.

(Abridged) We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics when the light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the Athena code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers. The method is ideally suited for problems in which characteristic velocities are non-relativistic, but still within a few percent or more of the speed of light. The method is an intermediate step towards algorithms for fully relativistic flows.

We use our variable Eddington tensor (VET) radiation hydrodynamics code to perform two-dimensional simulations to study the impact of radiation forces on atmospheres composed of dust and gas. Our setup closely follows that of Krumholz & Thompson, assuming that dust and gas are well-coupled and that the radiation field is characterized by blackbodies with temperatures >~ 80 K, as might be found in ultraluminous infrared galaxies. In agreement with previous work, we find that Rayleigh-Taylor instabilities develop in radiation supported atmospheres, leading to inhomogeneities that limit momentum exchange between radiation and dusty gas, and eventually providing a near balance of the radiation and gravitational forces. However, the evolution of the velocity and spatial distributions of the gas differs significantly from previous work, which utilized a less accurate flux-limited diffusion (FLD) method. Our VET simulations show continuous net acceleration of the gas, with no steady-state reached by the end of the simulation. In contrast, FLD results show little net acceleration of the gas and settle in to a quasi-steady, turbulent state with low velocity dispersion. The discrepancies result primarily from the inability of FLD to properly model the variation of the radiation field around structures that are less than a few optical depths across. We conclude that radiation feedback remains a viable mechanism for driving high-Mach number turbulence. We discuss implications for observed systems and global numerical simulations of feedback, but more realistic setups are needed to make robust observational predictions and assess the prospect of launching outflows with radiation.

A new mechanism to form a magnetic pressure supported, high temperature corona above the photosphere of an accretion disk is explored using three dimensional radiation magneto-hydrodynamic (MHD) simulations. The thermal properties of the disk are calculated self-consistently by balancing radiative cooling through the surfaces of the disk with heating due to dissipation of turbulence driven by magneto-rotational instability (MRI). As has been noted in previous work, we find the dissipation rate per unit mass increases dramatically with height above the mid-plane, in stark contrast to the alpha-disk model which assumes this quantity is a constant. Thus, we find that in simulations with a low surface density (and therefore a shallow photosphere), the fraction of energy dissipated above the photosphere is significant (about 3.4% in our lowest surface density model), and this fraction increases as surface density decreases. When a significant fraction of the accretion energy is dissipated in the optically thin photosphere, the gas temperature increases substantially and a high temperature, magnetic pressure supported corona is formed. The volume-averaged temperature in the disk corona is more than 10 times larger than at the disk mid-plane. Moreover, gas temperature in the corona is strongly anti-correlated with gas density, which implies the corona formed by MRI turbulence is patchy. This mechanism to form an accretion disk corona may help explain the observed relation between the spectral index and luminosity from AGNs, and the soft X-ray excess from some AGNs. It may also be relevant to spectral state changes in X-ray binaries.

We perform a systematic study of the dynamics of dust particles in protoplanetary disks with embedded planets using global 2-D and 3-D inviscid hydrodynamic simulations. Lagrangian particles have been implemented into magnetohydrodynamic code Athena with cylindrical coordinates. We find two distinct outcomes depending on the mass of the embedded planet. In the presence of a low mass planet ($8 M_{\oplus}$), two narrow gaps start to open in the gas on each side of the planet where the density waves shock. These shallow gaps can dramatically affect particle drift speed and cause significant, roughly axisymmetric dust depletion. On the other hand, a more massive planet ($>0.1 M_{J}$) carves out a deeper gap with sharp edges, which are unstable to the vortex formation. Particles with a wide range of sizes ($0.02<\Omega t_{s}<20$) are trapped and settle to the midplane in the vortex, with the strongest concentration for particles with $\Omega t_{s}\sim 1$. The dust concentration is highly elongated in the $\phi$ direction, and can be as wide as 4 disk scale heights in the radial direction. Dust surface density inside the vortex can be increased by more than a factor of 10$^2$ in a very non-axisymmetric fashion. For very big particles ($\Omega t_{s}\gg 1$) we find strong eccentricity excitation, in particular around the planet and in the vicinity of the mean motion resonances, facilitating gap opening there. Our results imply that in weakly turbulent protoplanetary disk regions (e.g. the "dead zone") dust particles with a very wide range of sizes can be trapped at gap edges and inside vortices induced by planets with $M_{p}<M_{J}$, potentially accelerating planetesimal and planet formation there, and giving rise to distinctive features that can be probed by ALMA and EVLA.

Procedures and results on hardware level detector calibration in Super-Kamiokande (SK) are presented in this paper. In particular, we report improvements made in our calibration methods for the experimental phase IV in which new readout electronics have been operating since 2008. The topics are separated into two parts. The first part describes the determination of constants needed to interpret the digitized output of our electronics so that we can obtain physical numbers such as photon counts and their arrival times for each photomultiplier tube (PMT). In this context, we developed an in-situ procedure to determine high-voltage settings for PMTs in large detectors like SK, as well as a new method for measuring PMT quantum efficiency and gain in such a detector. The second part describes the modeling of the detector in our Monte Carlo simulation, including in particular the optical properties of its water target and their variability over time. Detailed studies on the water quality are also presented. As a result of this work, we achieved a precision sufficient for physics analysis over a wide energy range (from a few MeV to above a TeV). For example, the charge determination was understood at the 1% level, and the timing resolution was 2.1 nsec at the one-photoelectron charge level and 0.5 nsec at the 100-photoelectron charge level.

We describe Pegasus, a new hybrid-kinetic particle-in-cell code tailored for the study of astrophysical plasma dynamics. The code incorporates an energy-conserving particle integrator into a stable, second-order--accurate, three-stage predictor-predictor-corrector integration algorithm. The constrained transport method is used to enforce the divergence-free constraint on the magnetic field. A delta-f scheme is included to facilitate a reduced-noise study of systems in which only small departures from an initial distribution function are anticipated. The effects of rotation and shear are implemented through the shearing-sheet formalism with orbital advection. These algorithms are embedded within an architecture similar to that used in the popular astrophysical magnetohydrodynamics code Athena, one that is modular, well-documented, easy to use, and efficiently parallelized for use on thousands of processors. We present a series of tests in one, two, and three spatial dimensions that demonstrate the fidelity and versatility of the code.

We report on the first phase of our study of cloud irradiation. We study irradiation by means of numerical, two-dimensional time-dependent radiation-hydrodynamic simulations of a cloud irradiated by a strong radiation. We adopt a very simple treatment of the opacity, neglect photoionization and gravity, and instead focus on assessing the role of the type and magnitude of the opacity on the cloud evolution. Our main result is that even relatively dense clouds that are radiatively heated (i.e., with significant absorption opacity) do not move as a whole instead they undergo a very rapid and major evolution in its shape, size and physical properties. In particular, the cloud and its remnants become optical thin within less than one sound crossing time and before they can travel over a significant distance (a distance of a few radii of the initial cloud). We also found that a cloud can be accelerated as a whole under quite extreme conditions, e.g., the opacity must be dominated by scattering. However, the acceleration due to the radiation force is relatively small and unless the cloud is optically thin the cloud quickly changes its size and shape. We discuss implications for the modelling and interpetation broad line regions of active galactic nuclei.

We study the long-term thermal stability of radiation dominated disks in which the vertical structure is determined self-consistently by the balance of heating due to dissipation of MHD turbulence driven by the magneto-rotational instability (MRI), and cooling due to radiation emitted at the photosphere. The calculations adopt the local shearing box approximation, and utilize the recently developed radiation transfer module in the Athena MHD code based on a variable Eddington tensor rather than an assumed local closure. After saturation of the MRI, in many cases the disk maintains a steady vertical structure for many thermal times. However, in every case in which the box size in the horizontal directions is at least one pressure scale height, fluctuations associated with MRI turbulence and dynamo action in the disk eventually trigger a thermal runaway which causes the disk to either expand or contract until the calculation must be terminated. During runaway, the dependence of the heating and cooling rates on total pressure satisfy the simplest criterion for classical thermal instability. We identify several physical reasons why the thermal runaway observed in our simulations differ from the standard alpha disk model, for example the advection of radiation contributes a non-negligible fraction to the vertical energy flux at the largest radiation pressure, most of the dissipation does not happen in the disk mid-plane, and the change of dissipation scale height with mid-plane pressure is slower than the change of density scale height. We discuss how and why our results differ from those published previously. Such thermal runaway behavior might have important implications for interpreting temporal variability in observed systems, but fully global simulations are required to study the saturated state before detailed predictions can be made.

We compile and analyse long term (~10 year) submillimetre (1.3, 0.87, 0.43 mm, submm) wavelength light curves of the Galactic centre black hole, Sagittarius A*. The 0.87 and 0.43 mm data are taken from the literature, while the majority of the 1.3 mm light curve is from previously unpublished SMA and CARMA data. We use Monte Carlo simulations to show that on minute to few hour time-scales the variability is consistent with a red noise process with a 230 GHz power spectrum slope of 2.3+0.8-0.6 at 95% confidence. The light curve is de-correlated (white noise) on very long (month to year) times. In order to identify the transition time between red and white noise, we model the light curves as a stochastic damped random walk process. The models allow a quantitative estimate of this physical characteristic time-scale of 8-4+3 hours at 230 GHz at 95% confidence, with consistent results at 345 and 690 GHz. This corresponds to ~10 orbital times or ~1 inflow (viscous) time at R = 3 Rs, a typical radius producing the 230 GHz emission as measured by very long baseline interferometry and found in theoretical accretion flow and jet models. This time-scale is significantly shorter (longer) than those measured at radio (near-infrared, NIR) wavelengths, and is marginally inconsistent with the same variability mechanism operating in the submm and NIR for the expected t ~ R^3/2 scaling. It is crudely consistent with the analogous time-scale inferred in studies of quasar optical light curves after accounting for the difference in emission radius. We find evidence that the submm variability persists at least down to the ISCO, if not the event horizon. These results can be compared quantitatively with similar analyses at different wavebands to test for connections between the variability mechanisms, and with light curves from theoretical models of accreting black holes.

We carry out a series of local, vertically stratified shearing box simulations of protoplanetary disks that include ambipolar diffusion and a net vertical magnetic field. The ambipolar diffusion profiles we employ correspond to 30AU and 100AU in a minimum mass solar nebula (MMSN) disk model, which consists of a far-UV-ionized surface layer and low-ionization disk interior. These simulations serve as a follow up to Simon et al. (2013), in which we found that without a net vertical field, the turbulent stresses that result from the magnetorotational instability (MRI) are too weak to account for observed accretion rates. The simulations in this work show a very strong dependence of the accretion stresses on the strength of the background vertical field; as the field strength increases, the stress amplitude increases. For gas to magnetic pressure ratios of 1e4 and 1e5, we find accretion rates between 1e-8 and 1e-7 solar masses per year. These accretion rates agree with observational constraints, suggesting a vertical magnetic field strength between 60 and 200 microgauss at 30AU and 10 and 30 microgauss at 100AU in a MMSN disk. Furthermore, the stress has a non-negligible component due to a magnetic wind. For sufficiently strong vertical field strengths, MRI turbulence is quenched, and the flow becomes largely laminar, with accretion proceeding through large scale correlations in the radial and toroidal field components as well as through the magnetic wind. In all simulations, the presence of a low ionization region near the disk mid-plane, which we call the ambipolar damping zone, results in reduced stresses there.

We present near-IR spectra of a sample of T Tauri, Herbig Ae/Be, and FU Ori objects. Using the FSPEC instrument on the Bok 90-inch telescope, we obtained K-band spectra with a resolution of ~3500. Here we present spectra of the v=2->0 and v=3->1 bandheads of ro-vibrational transitions of carbon monoxide. We observed these spectra over multiple epochs spaced by a few days and approximately one month. Several of our targets show CO emission or absorption features. However we see little evidence of variability in these features across multiple epochs. We compare our results with previous observations, and discuss the physical implications of non-variable CO emission across the sampled timescales.

We perform local, vertically stratified shearing-box MHD simulations of protoplanetary disks (PPDs) at a fiducial radius of 1 AU that take into account the effects of both Ohmic resistivity and ambipolar diffusion (AD). The magnetic diffusion coefficients are evaluated self-consistently from a look-up table based on equilibrium chemistry. We first show that the inclusion of AD dramatically changes the conventional picture of layered accretion. Without net vertical magnetic field, the system evolves into a toroidal field dominated configuration with extremely weak turbulence in the far-UV ionization layer that is far too inefficient to drive rapid accretion. In the presence of a weak net vertical field (plasma beta~10^5 at midplane), we find that the MRI is completely suppressed, resulting in a fully laminar flow throughout the vertical extent of the disk. A strong magnetocentrifugal wind is launched that efficiently carries away disk angular momentum and easily accounts for the observed accretion rate in PPDs. Moreover, under a physical disk wind geometry, all the accretion flow proceeds through a strong current layer with thickness of ~0.3H that is offset from disk midplane with radial velocity of up to 0.4 times the sound speed. Both Ohmic resistivity and AD are essential for the suppression of the MRI and wind launching. The efficiency of wind transport increases with increasing net vertical magnetic flux and the penetration depth of the FUV ionization. Our laminar wind solution has important implications on planet formation and global evolution of PPDs.

The saturation level of the magneto-rotational instability (MRI) in a strongly radiation dominated accretion disk is studied using a new Godunov radiation MHD code in the unstratified shearing box approximation. Since vertical gravity is neglected in this work, our focus is on how the MRI saturates in the optically thick mid-plane of the disk. We confirm that turbulence generated by the MRI is very compressible in the radiation dominated regime, as found by previous calculations using the flux-limited diffusion approximation. We also find little difference in the saturation properties in calculations that use a larger horizontal domain (up to four times the vertical scale height in the radial direction). However, in strongly radiation pressure dominated disks (one in which the radiation energy density reaches 1% of the rest mass energy density of the gas), we find Maxwell stress from the MRI turbulence is larger than the value produced when radiation pressure is replaced with the same amount of gas pressure. At the same time, the ratio between Maxwell stress and Reynolds stress is increased by almost a factor of 8 compared with the gas pressure dominated case. We suggest that this effect is caused by radiation drag, which acts like bulk viscosity and changes the effective magnetic Prandtl number of the fluid. Radiation viscosity significantly exceeds both the microscopic plasma viscosity and resistivity, ensuring that radiation dominated systems occupy the high magnetic Prandtl number regime. Nevertheless, we find radiative shear viscosity is negligible compared to the Maxwell and Reynolds stress in the flow. This may have important implications for the structure of radiation dominated accretion disks.

In the light of the recent discovery of a neutron star with a mass accurately determined to be almost two solar masses, it has been suggested that hyperons cannot play a role in the equation of state of dense matter in $\beta$-equilibrium. We re-examine this issue in the most recent development of the quark-meson coupling model. Within a relativistic Hartree-Fock approach and including the full tensor structure at the vector-meson-baryon vertices, we find that not only must hyperons appear in matter at the densities relevant to such a massive star but that the maximum mass predicted is completely consistent with the observation.

We perform 3D vertically-stratified local shearing-box ideal MHD simulations of the magnetorotational instability (MRI) that include a net vertical magnetic flux, which is characterized by beta_0 (ratio of gas pressure to magnetic pressure of the net vertical field at midplane). We have considered beta_0=10^2, 10^3 and 10^4 and in the first two cases the most unstable linear MRI modes are well resolved in the simulations. We find that the behavior of the MRI turbulence strongly depends on beta_0: The radial transport of angular momentum increases with net vertical flux, achieving alpha=0.08 for beta_0=10^4 and alpha>1.0 for beta_0=100, where alpha is the Shakura-Sunyaev parameter. A critical value lies at beta_0=10^3: For beta_0>10^3, the disk consists of a gas pressure dominated midplane and a magnetically dominated corona. The turbulent strength increases with net flux, and angular momentum transport is dominated by turbulent fluctuations. The magnetic dynamo that leads to cyclic flips of large-scale fields still exists, but becomes more sporadic as net flux increases. For beta_0<10^3, the entire disk becomes magnetic dominated. The turbulent strength saturates, and the magnetic dynamo is quenched. Stronger large-scale fields are generated with increasing net flux, which dominates angular momentum transport. A strong outflow is launched from the disk by the magnetocentrifugal mechanism, and the mass flux increases linearly with net vertical flux and shows sign of saturation at beta_0=10^2. However, the outflow is unlikely to be directly connected to a global wind: for beta_0>10^3, the large-scale field has no permanent bending direction due to dynamo activities, while for beta_0<10^3, the outflows from the top and bottom sides of the disk bend towards opposite directions, inconsistent with a physical disk wind geometry. Global simulations are needed to address the fate of the outflow.

We study wakes and gap opening by low mass planets in gaseous protoplanetary disks threaded by net vertical magnetic fields which drive magnetohydrodynamical (MHD) turbulence through the magnetorotational instabilty (MRI), using three dimensional simulations in the unstratified local shearing box approximation. The wakes, which are excited by the planets, are damped by shocks similar to the wake damping in inviscid hydrodynamic (HD) disks. Angular momentum deposition by shock damping opens gaps in both MHD turbulent disks and inviscid HD disks even for low mass planets, in contradiction to the "thermal criterion" for gap opening. To test the "viscous criterion", we compared gap properties in MRI-turbulent disks to those in viscous HD disks having the same stress, and found that the same mass planet opens a significantly deeper and wider gap in net vertical flux MHD disks than in viscous HD disks. This difference arises due to the efficient magnetic field transport into the gap region in MRI disks, leading to a larger effective \alpha within the gap. Thus, across the gap, the Maxwell stress profile is smoother than the gap density profile, and a deeper gap is needed for the Maxwell stress gradient to balance the planetary torque density. We also confirmed the large excess torque close to the planet in MHD disks, and found that long-lived density features (termed zonal flows) produced by the MRI can affect planet migration. The comparison with previous results from net toroidal flux/zero flux MHD simulations indicates that the magnetic field geometry plays an important role in the gap opening process. Overall, our results suggest that gaps can be commonly produced by low mass planets in realistic protoplanetary disks, and caution the use of a constant \alpha-viscosity to model gaps in protoplanetary disks.

The microscopic composition and properties of matter at super-saturation densities have been the subject of intense investigation for decades. The scarcity of experimental and observational data has lead to the necessary reliance on theoretical models. However, there remains great uncertainty in these models, which, of necessity, have to go beyond the over-simple assumption that high density matter consists only of nucleons and leptons. Heavy strange baryons, mesons and quark matter in different forms and phases have to be included to fulfil basic requirements of fundamental laws of physics. In this review the latest developments in construction of the Equation of State (EoS) of high-density matter at zero and finite temperature assuming different composition of the matter are surveyed. Critical comparison of model EoS with available observational data on neutron stars, including gravitational masses, radii and cooling patterns is presented. The effect of changing rotational frequency on the composition of neutron stars during their lifetime is demonstrated. Compatibility of EoS of high-density, low temperature compact objects and low density, high temperature matter created in heavy-ion collisions is discussed.

We perform global unstratified 3D magnetohydrodynamic simulations of an astrophysical boundary layer (BL) -- an interface region between an accretion disk and a weakly magnetized accreting object such as a white dwarf -- with the goal of understanding the effects of magnetic field on the BL. We use cylindrical coordinates with an isothermal equation of state and investigate a number of initial field geometries including toroidal, vertical, and vertical with zero net flux. Our initial setup consists of a Keplerian disk attached to a non-rotating star. In a previous work, we found that in hydrodynamical simulations, sound waves excited by shear in the BL were able to efficiently transport angular momentum and drive mass accretion onto the star. Here we confirm that in MHD simulations, waves serve as an efficient means of angular momentum transport in the vicinity of the BL, despite the magnetorotational instability (MRI) operating in the disk. In particular, the angular momentum current due to waves is at times larger than the angular momentum current due to MRI. Our results suggest that angular momentum transport in the BL and its vicinity is a global phenomenon occurring through dissipation of waves and shocks. This point of view is quite different from the standard picture of transport by a local anomalous turbulent viscosity. In addition to angular momentum transport, we also study magnetic field amplification within the BL. We find that the field is indeed amplified in the BL, but only by a factor of a few and remains subthermal.

We use local numerical simulations to investigate the strength and nature of magnetohydrodynamic (MHD) turbulence in the outer regions of protoplanetary disks, where ambipolar diffusion is the dominant non-ideal MHD effect. The simulations include vertical stratification and assume zero net vertical magnetic flux. We employ a super time-stepping technique to ameliorate the Courant restriction on the diffusive time step. We find that in idealized stratified simulations, with a spatially constant ambipolar Elsasser number Am, turbulence driven by the magnetorotational instability (MRI) behaves in a similar manner as in prior unstratified calculations. Turbulence dies away for Am < 1, and becomes progressively more vigorous as ambipolar diffusion is decreased. Near-ideal MHD behavior is recovered for Am > 1000. In the intermediate regime (10 < Am < 1000) ambipolar diffusion leads to substantial increases in both the period of the MRI dynamo cycle and the characteristic scales of magnetic field structures. To quantify the impact of ambipolar physics on disk accretion, we run simulations at 30 AU and 100 AU that include a vertical Am profile based upon far ultraviolet (FUV) ionized disk models. These models develop a vertically layered structure analogous to the Ohmic dead zone that is present at smaller radii. We find that, although the levels of surface turbulence can be strong (and consistent with constraints on turbulent line widths at these radii), the inferred accretion rates are at least an order of magnitude smaller than those observed in T Tauri stars. This discrepancy is very likely due to the assumption of zero vertical magnetic field in our simulations and suggests that vertical magnetic fields are essential for MRI-driven accretion in the outer regions of protoplanetary disks.

The non-linear regime of Rayleigh-Taylor instability (RTI) in a radiation supported atmosphere, consisting of two uniform fluids with different densities, is studied numerically. We perform simulations using our recently developed numerical algorithm for multi-dimensional radiation hydrodynamics based on a variable Eddington tensor as implemented in Athena, focusing on the regime where scattering opacity greatly exceeds absorption opacity. We find that the radiation field can reduce the growth and mixing rate of RTI, but this reduction is only significant when radiation pressure significantly exceeds gas pressure. Small scale structures are also suppressed in this case. In the non-linear regime, dense fingers sink faster than rarefied bubbles can rise, leading to asymmetric structures about the interface. By comparing the calculations that use a variable Eddington tensor (VET) versus the Eddington approximation, we demonstrate that anisotropy in the radiation field can affect the non-linear development of RTI significantly. We also examine the disruption of a shell of cold gas being accelerated by strong radiation pressure, motivated by models of radiation driven outflows in ultraluminous infrared galaxies. We find that when the growth rate of RTI is smaller than acceleration time scale, the amount of gas that would be pushed away by the radiation field is reduced due to RTI.

The nature of angular momentum transport in the boundary layers of accretion disks has been one of the central and long-standing issues of accretion disk theory. In this work we demonstrate that acoustic waves excited by supersonic shear in the boundary layer serve as an efficient mechanism of mass, momentum and energy transport at the interface between the disk and the accreting object. We develop the theory of angular momentum transport by acoustic modes in the boundary layer, and support our findings with 3D hydrodynamical simulations, using an isothermal equation of state. Our first major result is the identification of three types of global modes in the boundary layer. We derive dispersion relations for each of these modes that accurately capture the pattern speeds observed in simulations to within a few percent. Second, we show that angular momentum transport in the boundary layer is intrinsically nonlocal, and is driven by radiation of angular momentum away from the boundary layer into both the star and the disk. The picture of angular momentum transport in the boundary layer by waves that can travel large distances before dissipating and redistributing angular momentum and energy to the disk and star is incompatible with the conventional notion of local transport by turbulent stresses. Our results have important implications for semianalytical models that describe the spectral emission from boundary layers.

We study the so called hybrid stars, which are hadronic stars that contain a core of deconfined quarks. For this purpose, we make use of an extended version of the SU(3) chiral model. Within this approach, the degrees of freedom change naturally from hadrons (baryon octet) to quarks (u, d, s) as the temperature and/or density increases. At zero temperature we are still able to reproduce massive stars, even with the inclusion of hyperons.

We propose a new method to conserve the total energy to round-off error in grid-based codes for hydrodynamic simulations with self-gravity. A formula for the energy flux due to the work done by the the self-gravitational force is given, so the change in total energy can be written in conservative form. Numerical experiments with the code Athena show that the total energy is indeed conserved with our new algorithm and the new algorithm is second order accurate. We have performed a set of tests that show the numerical errors in the traditional, non-conservative algorithm can affect the dynamics of the system. The new algorithm only requires one extra solution of the Poisson equation, as compared to the traditional algorithm which includes self-gravity as a source term. If the Poisson solver takes a negligible fraction of the total simulation time, such as when FFTs are used, the new algorithm is almost as efficient as the original method. This new algorithm is useful in Eulerian hydrodynamic simulations with self-gravity, especially when results are sensitive to small energy errors, as for radiation pressure dominated flow.

Disk accretion onto a weakly magnetized central object, e.g. a star, is inevitably accompanied by the formation of a boundary layer near the surface, in which matter slows down from the highly supersonic orbital velocity of the disk to the rotational velocity of the star. We perform high resolution 2D hydrodynamical simulations in the equatorial plane of an astrophysical boundary layer with the goal of exploring the dynamics of non-axisymmetric structures that form there. We generically find that the supersonic shear in the boundary layer excites non-axisymmetric quasi-stationary acoustic modes that are trapped between the surface of the star and a Lindblad resonance in the disk. These modes rotate in a prograde fashion, are stable for hundreds of orbital periods, and have a pattern speed that is less than and of order the rotational velocity at the inner edge of the disk. The origin of these intrinsically global modes is intimately related to the operation of a corotation amplifier in the system. Dissipation of acoustic modes in weak shocks provides a universal mechanism for angular momentum and mass transport even in purely hydrodynamic (i.e. non-magnetized) boundary layers. We discuss the possible implications of these trapped modes for explaining the variability seen in accreting compact objects.

We carry out local three dimensional (3D) hydrodynamic simulations of planet-disk interaction in stratified disks with varied thermodynamic properties. We find that whenever the Brunt-Vaisala frequency (N) in the disk is nonzero, the planet exerts a strong torque on the disk in the vicinity of the planet, with a reduction in the traditional "torque cutoff". In particular, this is true for adiabatic perturbations in disks with isothermal density structure, as should be typical for centrally irradiated protoplanetary disks. We identify this torque with buoyancy waves, which are excited (when N is non-zero) close to the planet, within one disk scale height from its orbit. These waves give rise to density perturbations with a characteristic 3D spatial pattern which is in close agreement with the linear dispersion relation for buoyancy waves. The torque due to these waves can amount to as much as several tens of per cent of the total planetary torque, which is not expected based on analytical calculations limited to axisymmetric or low-m modes. Buoyancy waves should be ubiquitous around planets in the inner, dense regions of protoplanetary disks, where they might possibly affect planet migration.

We present the result of an indirect search for high energy neutrinos from WIMP annihilation in the Sun using upward-going muon (upmu) events at Super-Kamiokande. Datasets from SKI-SKIII (3109.6 days) were used for the analysis. We looked for an excess of neutrino signal from the Sun as compared with the expected atmospheric neutrino background in three upmu categories: stopping, non-showering, and showering. No significant excess was observed. The 90% C.L. upper limits of upward-going muon flux induced by WIMPs of 100 GeV/c$^2$ were 6.4$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ and 4.0$\times10^{-15}$ cm$^{-2}$ sec$^{-1}$ for the soft and hard annihilation channels, respectively. These limits correspond to upper limits of 4.5$\times10^{-39}$ cm$^{-2}$ and 2.7$\times10^{-40}$ cm$^{-2}$ for spin-dependent WIMP-nucleon scattering cross sections in the soft and hard annihilation channels, respectively.

We present simulations of Keck Interferometer ASTRA and VLTI GRAVITY observations of mock star fields in orbit within ~50 milliarcseconds of Sgr A*. Dual-field phase referencing techniques, as implemented on ASTRA and planned for GRAVITY, will provide the sensitivity to observe Sgr A* with infrared interferometers. Our results show an improvement in the confusion noise limit over current astrometric surveys, opening a window to study stellar sources in the region. Since the Keck Interferometer has only a single baseline, the improvement in the confusion limit depends on source position angles. The GRAVITY instrument will yield a more compact and symmetric PSF, providing an improvement in confusion noise which will not depend as strongly on position angle. Our Keck results show the ability to characterize the star field as containing zero, few, or many bright stellar sources. We are also able to detect and track a source down to mK~18 through the least confused regions of our field of view at a precision of ~200 microarcseconds along the baseline direction. This level of precision improves with source brightness. Our GRAVITY results show the potential to detect and track multiple sources in the field. GRAVITY will perform ~10 microarcsecond astrometry on a mK=16.3 source and ~200 microarcsecond astrometry on a mK=18.8 source in six hours of monitoring a crowded field. Monitoring the orbits of several stars will provide the ability to distinguish between multiple post-Newtonian orbital effects, including those due to an extended mass distribution around Sgr A* and to low-order General Relativistic effects. Early characterizations of the field by ASTRA including the possibility of a precise source detection, could provide valuable information for future GRAVITY implementation and observation.

This paper presents a detailed assessment of the ability of the 240 Skyrme interaction parameter sets in the literature to satisfy a series of criteria derived from macroscopic properties of nuclear matter in the vicinity of nuclear saturation density at zero temperature and their density dependence, derived by the liquid drop model, experiments with giant resonances and heavy-ion collisions. The objective is to identify those parameterizations which best satisfy the current understanding of the physics of nuclear matter over a wide range of applications. Out of the 240 models, only 16 are shown to satisfy all these constraints. Additional, more microscopic, constraints on density dependence of the neutron and proton effective mass beta-equilibrium matter, Landau parameters of symmetric and pure neutron nuclear matter, and observational data on high- and low-mass cold neutron stars further reduce this number to 5, a very small group of recommended Skyrme parameterizations to be used in future applications of the Skyrme interaction of nuclear matter related observables. Full information on partial fulfillment of individual constraints by all Skyrme models considered is given. The results are discussed in terms of the physical interpretation of the Skyrme interaction and the validity of its use in mean-field models. Future work on application of the Skyrme forces, selected on the basis of variables of nuclear matter, in Hartree-Fock calculation of properties of finite nuclei, is outlined.

The symmetry energy contribution to the nuclear Equation of State (EoS) impacts various phenomena in nuclear astrophysics, nuclear structure, and nuclear reactions. Its determination is a key objective of contemporary nuclear physics with consequences for the understanding of dense matter within neutron stars. We examine the results of laboratory experiments that have provided initial constraints on the nuclear symmetry energy and its density dependence at and somewhat below normal nuclear matter density. Some of these constraints have been derived from properties of nuclei. Others have been derived from the nuclear response to electroweak and hadronic probes. We also examine the most frequently used theoretical models that predict the symmetry energy and its slope. By comparing existing constraints on the symmetry pressure to theories, we demonstrate how the contribution of the three-body force, an essential ingredient in neutron matter models, can be determined.

Barred galaxies are known to possess magnetic fields that may affect the properties of bar substructures such as dust lanes and nuclear rings. We use two-dimensional high-resolution magnetohydrodynamic (MHD) simulations to investigate the effects of magnetic fields on the formation and evolution of such substructures as well as on the mass inflow rates to the galaxy center. The gaseous medium is assumed to be infinitesimally-thin, isothermal, non-self-gravitating, and threaded by initially uniform, azimuthal magnetic fields. We find that there exists an outermost x1-orbit relative to which gaseous responses to an imposed stellar bar potential are completely different between inside and outside. Inside this orbit, gas is shocked into dust lanes and infalls to form a nuclear ring. Magnetic fields are compressed in dust lanes, reducing their peak density. Magnetic stress removes further angular momentum of the gas at the shocks, temporarily causing the dust lanes to bend into an 'L' shape and eventually leading to a smaller and more centrally distributed ring than in unmagnetized models. The mass inflow rates in magnetized models correspondingly become larger, by more than two orders of magnitude when the initial fields have an equipartition value with thermal energy, than in the unmagnetized counterparts. Outside the outermost x1-orbit, on the other hand, an MHD dynamo due to the combined action of the bar potential and background shear operates near the corotation and bar-end regions, efficiently amplifying magnetic fields. The amplified fields shape into trailing magnetic arms with strong fields and low density. The base of the magnetic arms has a thin layer in which magnetic fields with opposite polarity reconnect via a tearing-mode instability. This produces numerous magnetic islands with large density which propagate along the arms to turn the outer disk into a highly chaotic state.

GUT monopoles captured by the Sun's gravitation are expected to catalyze proton decays via the Callan-Rubakov process. In this scenario, protons, which initially decay into pions, will ultimately produce \nu_{e}, \nu_{\mu} and \bar{\nu}_{\mu}. After undergoing neutrino oscillation, all neutrino species appear when they arrive at the Earth, and can be detected by a 50,000 metric ton water Cherenkov detector, Super-Kamiokande (SK). A search for low energy neutrinos in the electron total energy range from 19 to 55 MeV was carried out with SK and gives a monopole flux limit of F_M(\sigma_0/1 mb) < 6.3 \times 10^{-24} (\beta_M/10^{-3})^2 cm^{-2} s^{-1} sr^{-1} at 90% C.L., where \beta_M is the monopole velocity in units of the speed of light and \sigma_0 is the catalysis cross section at \beta_M=1. The obtained limit is more than eight orders of magnitude more stringent than the current best cosmic-ray supermassive monopole flux limit, F_M < 1 \times 10^{-15} cm^{-2} s^{-1} sr^{-1} for \beta_M < 10^{-3} and also two orders of magnitude lower than the result of the Kamiokande experiment, which used a similar detection method.

The intracluster medium of galaxy clusters is a weakly collisional, high-beta plasma in which the transport of heat and momentum occurs primarily along magnetic-field lines. Anisotropic heat conduction allows convective instabilities to be driven by temperature gradients of either sign, the magnetothermal instability (MTI) in the outskirts of non-isothermal clusters and the heat-flux buoyancy-driven instability (HBI) in their cooling cores. We employ the Athena MHD code to investigate the nonlinear evolution of these instabilities, self-consistently including the effects of anisotropic viscosity (i.e. Braginskii pressure anisotropy), anisotropic conduction, and radiative cooling. We highlight the importance of the microscale instabilities that inevitably accompany and regulate the pressure anisotropies generated by the HBI and MTI. We find that, in all but the innermost regions of cool-core clusters, anisotropic viscosity significantly impairs the ability of the HBI to reorient magnetic-field lines orthogonal to the temperature gradient. Thus, while radio-mode feedback appears necessary in the central few tens of kpc, conduction may be capable of offsetting radiative losses throughout most of a cool core over a significant fraction of the Hubble time. Magnetically-aligned cold filaments are then able to form by local thermal instability. Viscous dissipation during the formation of a cold filament produces accompanying hot filaments, which can be searched for in deep Chandra observations of nearby cool-core clusters. In the case of the MTI, anisotropic viscosity maintains the coherence of magnetic-field lines over larger distances than in the inviscid case, providing a natural lower limit for the scale on which the field can fluctuate freely. In the nonlinear state, the magnetic field exhibits a folded structure in which the field-line curvature and field strength are anti-correlated.

We present a multimoment technique for signal classification and apply it to the detection of fast radio transients in incoherently dedispersed data. Specifically, we define a spectral modulation index in terms of the fractional variation in intensity across a spectrum. A signal whose intensity is distributed evenly across the entire band has a much lower modulation index than a spectrum with the same intensity localized in a single channel. We are interested in broadband pulses and use the modulation index to excise narrowband radio frequency interference (RFI) by applying a modulation index threshold above which candidate events are removed. The technique is tested both with simulations and using data from sources of known radio pulses (RRAT J1928+15 and giant pulses from the Crab pulsar). We find that our technique is effective at eliminating not only narrowband RFI but also spurious signals from bright, real pulses that are dedispersed at incorrect dispersion measures. The method is generalized to coherent dedispersion, image cubes, and astrophysical narrowband signals that are steady in time. We suggest that the modulation index, along with other statistics using higher-order moments, should be incorporated into signal detection pipelines to characterize and classify signals.

We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code which solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-LTE effects. The module is based on well-known and well-tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator-splitting and we describe this approach in detail. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself, and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate, and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad-hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g. LTE grey) problems can be comparable to or less than a single timestep of Athena's MHD integrators, and only few times more expensive than that for more general problems. (Abridged)

We describe a numerical algorithm to integrate the equations of radiation magnetohydrodynamics in multidimensions using Godunov methods. This algorithm solves the radiation moment equations in the mixed frame, without invoking any diffusion-like approximations. The moment equations are closed using a variable Eddington tensor whose components are calculated from a formal solution of the transfer equation at a large number of angles using the method of short characteristics. We use a comprehensive test suite to verify the algorithm, including convergence tests of radiation-modified linear acoustic and magnetosonic waves, the structure of radiation modified shocks, and two-dimensional tests of photon bubble instability and the ablation of dense clouds by an intense radiation field. These tests cover a very wide range of regimes, including both optically thick and thin flows, and ratios of the radiation to gas pressure of at least 10^{-4} to 10^{4}. Across most of the parameter space, we find the method is accurate. However, the tests also reveal there are regimes where the method needs improvement, for example when both the radiation pressure and absorption opacity are very large. We suggest modifications to the algorithm that will improve accuracy in this case. We discuss the advantages of this method over those based on flux-limited diffusion. In particular, we find the method is not only substantially more accurate, but often no more expensive than the diffusion approximation for our intended applications.

The inner regions of barred galaxies contain substructures such as off-axis shocks, nuclear rings, and nuclear spirals. These substructure may affect star formation, and control the activity of a central black hole (BH) by determining the mass inflow rate. We investigate the formation and properties of such substructures using high-resolution, grid-based hydrodynamic simulations. The gaseous medium is assumed to be infinitesimally-thin, isothermal, and non-self-gravitating. The stars and dark matter are represented by a static gravitational potential with four components: a stellar disk, the bulge, a central BH, and the bar. To investigate various galactic environments, we vary the gas sound speed c_s as well as the mass of the central BH M_BH. Once the flow has reached a quasi-steady state, off-axis shocks tend to move closer to the bar major axis as c_s increases. Nuclear rings shrink in size with increasing c_s, but are independent of M_BH, suggesting that ring position is not determined by the Lindblad resonances. Rings in low-c_s models are narrow since they are occupied largely by gas on x2-orbits and well decoupled from nuclear spirals, while they become broad because of large thermal perturbations in high-c_s models. Nuclear spirals persist only when either c_s is small or M_BH is large; they would otherwise be destroyed completely by the ring material on eccentric orbits. The shape and strength of nuclear spirals depend sensitively on c_s and M_BH such that they are leading if both c_s and M_BH are small, weak trailing if c_s is small and M_BH is large, and strong trailing if both c_s and M_BH are large. While the mass inflow rate toward the nucleus is quite small in low-c_s models because of the presence of a narrow nuclear ring, it becomes larger than 0.01 Msun/yr when c_s is large, providing a potential explanation of nuclear activity in Seyfert galaxies.

Thermal instability (TI) can strongly affect the structure and dynamics of the interstellar medium (ISM) in the Milky Way and other disk galaxies. Thermal conduction plays an important role in the TI by stabilizing small scales and limiting the size of the smallest condensates. In the magnetized ISM, however, heat is conducted anisotropically (primarily along magnetic field lines). We investigate the effects of anisotropic thermal conduction on the nonlinear regime of the TI by performing two-dimensional magnetohydrodynamic simulations. We present models with magnetic fields of different initial geometries and strengths, and compare them to hydrodynamic models with isotropic conduction. We find anisotropic conduction does not significantly alter the overall density and temperature statistics in the saturated state of the TI. However, it can strongly affect the shapes and sizes of cold clouds formed by the TI. For example, for uniform initial fields long filaments of cold gas are produced that are reminiscent of some observed HI clouds. For initially tangled fields, such filaments are not produced. We also show that anisotropic conduction suppresses turbulence generated by evaporative flows from the surfaces of cold blobs, which may have implications for mechanisms for driving turbulence in the ISM.

Density waves excited by planets embedded in protoplanetary disks play a central role in planetary migration and gap opening processes. We carry out 2D shearing sheet simulations to study the linear regime of wave evolution with the grid-based code Athena, and provide detailed comparisons with the theoretical predictions. Low mass planets (down to ~0.03 Earth mass at 1 AU) and high spatial resolution (256 grid points per scale height) are chosen to mitigate the effects of wave nonlinearity. To complement the existing numerical studies, we focus on the primary physical variables such as the spatial profile of the wave, torque density, and the angular momentum flux carried by the wave, instead of secondary quantities such as the planetary migration rate. Our results show percent level agreement with theory in both physical and Fourier space. New phenomena such as the change of the toque density sign far from the planet are discovered and discussed. Also, we explore the effect of the numerical algorithms, and find that a high order of accuracy, high resolution, and an accurate planetary potential are crucial to achieve good agreement with the theory. We find that the use of a too large time-step without properly resolving the dynamical time scale around the planet produces incorrect results, and may lead to spurious gap opening. Global simulations of planet migration and gap opening violating this requirement may be affected by spurious effects resulting in e.g. the incorrect planetary migration rate and gap opening mass.

We investigate numerically the propagation of density waves excited by a low-mass planet in a protoplanetary disk in the nonlinear regime, using 2D local shearing box simulations with the grid-based code Athena at high spatial resolution (256 grid points per scale height h). The nonlinear evolution results in the wave steepening into a shock, causing damping and angular momentum transfer to the disk. On long timescales this leads to spatial redistribution of the disk density, causing migration feedback and potentially resulting in gap opening. Previous numerical studies concentrated on exploring these secondary phenomena as probes of the nonlinear wave evolution. Here we focus on exploring the evolution of the basic wave properties, such as its density profile evolution, shock formation, post-shock wave behavior, and provide comparison with analytical theory. The generation of potential vorticity at the shock is computed analytically and is subsequently verified by simulations and used to pinpoint the shock location. We confirm the theoretical relation between the shocking length and the planet mass (including the effect of the equation of state), and the post-shock decay of the angular momentum flux carried by the wave. The post-shock evolution of the wave profile is explored, and we quantitatively confirm its convergence to the theoretically expected N-wave shape. The accuracy of various numerical algorithms used to compute the nonlinear wave evolution is also investigated: we find that higher order spatial reconstruction and high resolution are crucial for capturing the shock formation correctly.

We present a mechanism that may seed compact stellar objects with stable lumps of quark matter, or {\it strangelets}, through the self-annihilation of gravitationally accreted WIMPs. We show that dark matter particles with masses above a few GeV may provide enough energy in the nuclear medium for quark deconfinement and subsequent strangelet formation. If this happens this effect may then trigger a partial or full conversion of the star into a strange star. We set a new limit on the WIMP mass in the few-GeV range that seems to be consistent with recent indications in dark matter direct detection experiments.

Grid-based magnetohydrodynamic (MHD) simulations have proven invaluable for the study of astrophysical accretion disks. However, the fact that angular momentum transport in disks is mediated by MHD turbulence (with structure down to very small scales) raises the concern that the properties of the modeled accretion disks are affected by the finite numerical resolution of the simulation. By implementing an orbital advection algorithm into the Athena code in cylindrical geometry, we have performed a set of global (but unstratified) Newtonian disk simulations extending up to resolutions previously unattained. We study the convergence of these models as a function of spatial resolution and initial magnetic field geometry. The usual viscosity parameter ($\alpha$) or the ratio of thermal-to-magnetic pressure ($\beta$) are found to be poor diagnostics of convergence, whereas the average tilt angle of the magnetic field in the $(r,\phi)$-plane is a very good diagnostic of convergence. We suggest that this is related to the saturation of the MHD turbulence via parasitic modes of the magnetorotational instability. Even in the case of zero-net magnetic flux, we conclude that our highest resolution simulations (with 32-zones and 64-zones per vertical scale height) have achieved convergence. Our global simulations reach resolutions comparable to those used in local, shearing box models of MHD disk turbulence. We find that the saturation predictors derived from local simulations correspond well to the instantaneous correlations between local flux and stress found in our global simulations. However, the conservation of magnetic flux implicit in local models is not realized in our global disks. Thus, the magnetic connectivity of an accretion disk represents physics that is truly global and cannot be captured in any ab-initio local model.

We study the role of ambipolar diffusion (AD) on the non-linear evolution of the MRI in protoplanetary disks using the strong coupling limit, which applies when the electron recombination time is much shorter than the orbital time. The effect of AD in this limit is characterized by the dimensionless number Am, the frequency of which neutral particles collide with ions normalized to the orbital frequency. We perform three-dimensional unstratified shearing-box simulations of the MRI over a wide range of Am as well as different magnetic field strengths and geometries. The saturation level of the MRI turbulence depends on the magnetic geometry and increases with the net magnetic flux. There is an upper limit to the net flux for sustained turbulence, corresponding to the requirement that the most unstable vertical wavelength be less than the disk scale height. Correspondingly, at a given Am, there exists a maximum value of the turbulent stress alpha_max. For Am<1, the largest stress is associated with a field geometry that has both net vertical and toroidal flux. In this case, we confirm the results of linear analyses that show the fastest growing mode has a non-zero radial wave number with growth rate exceeding the pure vertical field case. We find there is a very tight correlation between the turbulent stress (alpha) and the plasma beta=P_gas/P_mag~1/(2alpha) at the saturated state of the MRI turbulence regardless of field geometry, and alpha_max rapidly decreases with decreasing Am. In particular, we quote alpha_max~0.007 for Am=1 and alpha_max~0.0006 for Am=0.1.

We describe a new Godunov algorithm for relativistic magnetohydrodynamics (RMHD) that combines a simple, unsplit second order accurate integrator with the constrained transport (CT) method for enforcing the solenoidal constraint on the magnetic field. A variety of approximate Riemann solvers are implemented to compute the fluxes of the conserved variables. The methods are tested with a comprehensive suite of multidimensional problems. These tests have helped us develop a hierarchy of correction steps that are applied when the integration algorithm predicts unphysical states due to errors in the fluxes, or errors in the inversion between conserved and primitive variables. Although used exceedingly rarely, these corrections dramatically improve the stability of the algorithm. We present preliminary results from the application of these algorithms to two problems in RMHD: the propagation of supersonic magnetized jets, and the amplification of magnetic field by turbulence driven by the relativistic Kelvin-Helmholtz instability (KHI). Both of these applications reveal important differences between the results computed with Riemann solvers that adopt different approximations for the fluxes. For example, we show that use of Riemann solvers which include both contact and rotational discontinuities can increase the strength of the magnetic field within the cocoon by a factor of ten in simulations of RMHD jets, and can increase the spectral resolution of three-dimensional RMHD turbulence driven by the KHI by a factor of 2. This increase in accuracy far outweighs the associated increase in computational cost. Our RMHD scheme is publicly available as part of the Athena code.

In a hot, dilute, magnetized accretion flow, the electron mean-free path can be much greater than the Larmor radius, thus thermal conduction is anisotropic and along magnetic field lines. In this case, if the temperature decreases outward, the flow may be subject to a buoyancy instability (the magnetothermal instability, or MTI). The MTI amplifies the magnetic field, and aligns field lines with the radial direction. If the accretion flow is differentially rotating, the magnetorotational instability (MRI) may also be present. Using two-dimensional, time-dependent magnetohydrodynamic simulations, we investigate the interaction between these two instabilities. We use global simulations that span over two orders of magnitude in radius, centered on the region around the Bondi radius where the infall time of gas is longer than the growth time of both the MTI and MRI. Significant amplification of the magnetic field is produced by both instabilities, although we find that the MTI primarily amplifies the radial component, and the MRI primarily the toroidal component, of the field, respectively. Most importantly, we find that if the MTI can amplify the magnetic energy by a factor $F_t$, and the MRI by a factor $F_r$, then when the MTI and MRI are both present, the magnetic energy can be amplified by a factor of $F_t \cdot F_r$. We therefore conclude that amplification of the magnetic energy by the MTI and MRI operates independently. We also find that the MTI contributes to the transport of angular momentum, because radial motions induced by the MTI increase the Maxwell (by amplifying the magnetic field) and Reynolds stresses. Finally, we find that thermal conduction decreases the slope of the radial temperature profile. The increased temperature near the Bondi radius decreases the mass accretion rate.

The latest observation of a Shapiro delay of the binary millisecond pulsar J1614-2230 by Demorest et al. Nature 467 1081 (2010) yielded the pulsar gravitational mass to be 1.97 +/- 0.04 solar mass, the heaviest observed pulsar to-date. This result produces a stringent constraint on Equation(s) of State (EoS) of high density neutron star matter. One of the main conclusions of Demorest et al. was that their result makes the presence of non-nucleonic components in the neutron star matter unlikely. We compare the result with our recent work and conclude that hyperons in high-density matter are fully consistent with the observation and that their presence is a necessary consequence of general physical laws.

In a thermally bistable medium, cold, dense gas is separated from warm, rareified gas by thin phase transition layers, or fronts, in which heating, radiative cooling, thermal conduction, and convection of material are balanced. We calculate the steady-state structure of such fronts in the presence of magnetic fields, including the processes of ion-neutral drift and ion-neutral frictional heating. We find that ambipolar diffusion efficiently transports the magnetic field across the fronts, leading to a flat magnetic field strength profile. The thermal profiles of such fronts are not significantly different from those of unmagnetized fronts. The near uniformity of the magnetic field strength across a front is consistent with the flat field strength-gas density relation that is observed in diffuse interstellar gas.

We show that self-annihilating neutralino WIMP dark matter accreted onto neutron stars may provide a mechanism to seed compact objects with long-lived lumps of strange quark matter, or strangelets, for WIMP masses above a few GeV. This effect may trigger a conversion of most of the star into a strange star. We use an energy estimate for the long-lived strangelet based on the Fermi gas model combined with the MIT bag model to set a new limit on the possible values of the WIMP mass that can be especially relevant for subdominant species of massive neutralinos.

(Abridged) We present local 2D and 3D hybrid numerical simulations of particles and gas in the midplane of protoplanetary disks (PPDs) using the Athena code. The particles are coupled to gas aerodynamically, with particle-to-gas feedback included. Magnetorotational turbulence is ignored as an approximation for the dead zone of PPDs, and we ignore particle self-gravity to study the precursor of planetesimal formation. Our simulations include a wide size distribution of particles, ranging from strongly coupled particles with dimensionless stopping time tau_s=Omega t_stop=1e-4 to marginally coupled ones with tau_s=1 (where Omega is the orbital frequency, t_stop is the particle friction time), and a wide range of solid abundances. Our main results are: 1. Particles with tau_s>=0.01 actively participate in the streaming instability, generate turbulence and maintain the height of the particle layer before Kelvin-Helmholtz instability is triggered. 2. Strong particle clumping as a consequence of the streaming instability occurs when a substantial fraction of the solids are large (tau_s>=0.01) and when height-integrated solid to gas mass ratio Z is super-solar. 3. The radial drift velocity is reduced relative to the conventional Nakagawa-Sekiya-Hayashi (NSH) model, especially at high Z. We derive a generalized NSH equilibrium solution for multiple particle species which fits our results very well. 4. Collision velocity between particles with tau_s>=0.01 is dominated by differential radial drift, and is strongly reduced at larger Z. 5. There exist two positive feedback loops with respect to the enrichment of local disk solid abundance and grain growth. All these effects promote planetesimal formation.

The Athena MHD code has been extended to integrates the motion of particles coupled with the gas via aerodynamic drag, in order to study the dynamics of gas and solids in protoplanetary disks and the formation of planetesimals. Our particle-gas hybrid scheme is based on a second order predictor-corrector method. Careful treatment of the momentum feedback on the gas guarantees exact conservation. The hybrid scheme is stable and convergent in most regimes relevant to protoplanetary disks. We describe a semi-implicit integrator generalized from the leap-frog approach. In the absence of drag force, it preserves the geometric properties of a particle orbit. We also present a fully-implicit integrator that is unconditionally stable for all regimes of particle-gas coupling. Using our hybrid code, we study the numerical convergence of the non-linear saturated state of the streaming instability. We find that gas flow properties are well converged with modest grid resolution (128 cells per pressure length \eta r for dimensionless stopping time tau_s=0.1), and equal number of particles and grid cells. On the other hand, particle clumping properties converge only at higher resolutions, and finer resolution leads to stronger clumping before convergence is reached. Finally, we find that measurement of particle transport properties resulted from the streaming instability may be subject to error of about 20%.

We describe the implementation of the shearing box approximation for the study of the dynamics of accretion disks in the Athena magnetohydrodynamics (MHD) code. Second-order Crank-Nicholson time differencing is used for the Coriolis and tidal gravity source terms that appear in the momentum equation for accuracy and stability. We show this approach conserves energy for epicyclic oscillations in hydrodynamic flows to round-off error. In the energy equation, the tidal gravity source terms are differenced as the gradient of an effective potential in a way which guarantees that total energy (including the gravitational potential energy) is also conserved to round-off error. We introduce an orbital advection algorithm for MHD based on constrained transport to preserve the divergence-free constraint on the magnetic field. This algorithm removes the orbital velocity from the time step constraint, and makes the truncation error more uniform in radial position. Modifications to the shearing box boundary conditions applied at the radial boundaries are necessary to conserve the total vertical magnetic flux. In principle similar corrections are also required to conserve mass, momentum and energy, however in practice we find the orbital advection method conserves these quantities to better than 0.03% over hundreds of orbits. The algorithms have been applied to studies of the nonlinear regime of the MRI in very wide (up to 32 scale heights) horizontal domains.

The streaming instability (SI) provides a promising mechanism for planetesimal formation because of its ability to concentrate solids into dense clumps. The degree of clumping strongly depends on the height-integrated solid to gas mass ratio Z in protoplanetary disks (PPDs). In this letter, we show that the magnitude of the radial pressure gradient (RPG) which drives the SI (characterized by $q\equiv{\eta}v_K/c_s$, where ${\eta}v_K$ is the reduction of Keplerian velocity due to the RPG and $c_s$ is the sound speed) also strongly affects clumping. We present local two-dimensional hybrid numerical simulations of aerodynamically coupled particles and gas in the midplane of PPDs. Magnetic fields and particle self-gravity are ignored. We explore three different RPG values appropriate for typical PPDs: $q=0.025, 0.05$ and 0.1. For each $q$ value, we consider four different particle size distributions ranging from sub millimeter to meter sizes and run simulations with solid abundance from Z=0.01 up to Z=0.07. We find that a small RPG strongly promotes particle clumping in that: 1) At fixed particle size distribution, the critical solid abundance $Z_{\rm crit}$ above which particle clumping occurs monotonically increases with $q$; 2) At fixed Z, strong clumping can occur for smaller particles when $q$ is smaller. Therefore, we expect planetesimals to form preferentially in regions of PPDs with a small RPG.

Recent observations by Chandra and XMM-Newton indicate there are complex structures at the cores of galaxy clusters, such as cavities and filaments. One plausible model for the formation of such structures is the interaction of radio jets with the intracluster medium (ICM). To investigate this idea, we use three-dimensional magnetohydrodynamic simulations including anisotropic (Braginskii) viscosity to study the effect of magnetic fields on the evolution and morphology of buoyant bubbles in the ICM. We investigate a range of different initial magnetic field geometries and strengths, and study the resulting x-ray surface brightness distribution for comparison to observed clusters. Magnetic tension forces and viscous transport along field lines tend to suppress instabilities parallel, but not perpendicular, to field lines. Thus, the evolution of the bubble depends strongly on the initial field geometry. We find toroidal field loops initially confined to the interior of the bubble are best able reproduce the observed cavity structures.

We examine the effects of density stratification on magnetohydrodynamic turbulence driven by the magnetorotational instability in local simulations that adopt the shearing box approximation. Our primary result is that, even in the absence of explicit dissipation, the addition of vertical gravity leads to convergence in the turbulent energy densities and stresses as the resolution increases, contrary to results for zero net flux, unstratified boxes. The ratio of total stress to midplane pressure has a mean of ~0.01, although there can be significant fluctuations on long (>~50 orbit) timescales. We find that the time averaged stresses are largely insensitive to both the radial or vertical aspect ratio of our simulation domain. For simulations with explicit dissipation, we find that stratification extends the range of Reynolds and magnetic Prandtl numbers for which turbulence is sustained. Confirming the results of previous studies, we find oscillations in the large scale toroidal field with periods of ~10 orbits and describe the dynamo process that underlies these cycles.

We measure the turbulent resistivity in the nonlinear regime of the MRI, and evaluate the turbulent magnetic Prandtl number. We perform a set of numerical simulations with the Eulerian finite volume codes Athena and Ramses in the framework of the shearing box model. We consider models including explicit dissipation coefficients and magnetic field topologies such that the net magnetic flux threading the box in both the vertical and azimuthal directions vanishes. We first demonstrate good agreement between the two codes by comparing the properties of the turbulent states in simulations having identical microscopic diffusion coefficients (viscosity and resistivity). We find the properties of the turbulence do not change when the box size is increased in the radial direction, provided it is elongated in the azimuthal direction. To measure the turbulent resistivity in the disk, we impose a fixed electromotive force on the flow and measure the amplitude of the saturated magnetic field that results. We obtain a turbulent resistivity that is in rough agreement with mean field theories like the Second Order Smoothing Approximation. The numerical value translates into a turbulent magnetic Prandtl number Pm_t of order unity. Pm_t appears to be an increasing function of the forcing we impose. It also becomes smaller as the box size is increased in the radial direction, in good agreement with previous results obtained in very large boxes. Our results are in general agreement with other recently published papers studying the same problem but using different methodology. Thus, our conclusion that Pm_t is of order unity appears robust.

We consider the scenario of a magnetic field orthogonal to a front separating two media of different temperatures and densities, such as cold and warm interstellar gas, in a 2-D plane-parallel geometry. A linear stability analysis is performed to assess the behavior of both evaporation and condensation fronts when subject to incompressible, corrugational perturbations with wavelengths larger than the thickness of the front. We discuss the behavior of fronts in both super-Alfvenic and sub-Alfvenic flows. Since the propagation speed of fronts is slow in the ISM, it is the sub-Alfvenic regime that is relevant, and magnetic fields are a significant influence on front dynamics. In this case we find that evaporation fronts, which are unstable in the hydrodynamic regime, are stabilized. Condensation fronts are unstable, but for parameters typical of the neutral ISM the growth rates are so slow that steady state fronts are effectively stable. However, the instability may become important if condensation proceeds at a sufficiently fast rate. This paper is the first in a series exploring the linear and nonlinear effects of magnetic field strength and orientation on the corrugational instability, with the ultimate goal of addressing outstanding questions about small-scale ISM structure.

This is a white paper submitted to the Stars and Stellar Evolution (SSE) Science Frontier Panel (SFP) of the NRC's Astronomy and Astrophysics 2010 Decadal Survey. The white paper is endorsed by the American Physical Society's (APS) Topical Group on Plasma Astrophysics (GPAP).

We present a set of numerical simulations addressing the effects of magnetic field strength and orientation on the flow-driven formation of molecular clouds. Fields perpendicular to the flows sweeping up the cloud can efficiently prevent the formation of massive clouds but permit the build-up of cold, diffuse filaments. Fields aligned with the flows lead to substantial clouds, whose degree of fragmentation and turbulence strongly depends on the background field strength. Adding a random field component leads to a "selection effect" for molecular cloud formation: high column densities are only reached at locations where the field component perpendicular to the flows is vanishing. Searching for signatures of colliding flows should focus on the diffuse, warm gas, since the cold gas phase making up the cloud will have lost the information about the original flow direction because the magnetic fields redistribute the kinetic energy of the inflows.

We present two-dimensional hydrodynamic models of thermally driven winds from highly irradiated, close-in extra-solar planets. We adopt a very simple treatment of the radiative heating processes at the base of the wind, and instead focus on the differences between the properties of outflows in multidimensions in comparison to spherically symmetric models computed with the same methods. For hot (T > 2 x 10^{4} K) or highly ionized gas, we find strong (supersonic) polar flows are formed above the planet surface which produce weak shocks and outflow on the night-side. In comparison to a spherically symmetric wind with the same parameters, the sonic surface on the day-side is much closer to the planet surface in multidimensions, and the total mass loss rate is reduced by almost a factor of four. We also compute the steady-state structure of interacting planetary and stellar winds. Both winds end in a termination shock, with a parabolic contact discontinuity which is draped over the planet separating the two shocked winds. The planetary wind termination shock and the sonic surface in the wind are well separated, so that the mass loss rate from the planet is essentially unaffected. However, the confinement of the planetary wind to the small volume bounded by the contact discontinuity greatly enhances the column density close to the planet, which might be important for the interpretation of observations of absorption lines formed by gas surrounding transiting planets.

We study the probability distribution function (PDF) of the mass density in simulations of supersonic turbulence with properties appropriate for molecular clouds. For this study we use Athena, a new higher-order Godunov code. We find there are surprisingly similar relationships between the mean of the time-averaged PDF and the turbulent Mach number for driven hydrodynamic and strong-field MHD turbulence. There is, however, a large scatter about these relations, indicating a high level of temporal and spatial variability in the PDF. Thus, the PDF of the mass density is unlikely to be a good measure of magnetic field strength. We also find the PDF of decaying MHD turbulence deviates from the mean-Mach relation found in the driven case. This implies that the instantaneous Mach number alone is not enough to determine the statistical properties of turbulence that is out of equilibrium. The scatter about the mean-Mach relation for driven turbulence, along with the large departure of decaying turbulence PDFs from those of driven turbulence, may illuminate one factor contributing to the large observed cloud-to-cloud variation in the star formation rate per solar mass.

We study energy dissipation and heating by supersonic MHD turbulence in molecular clouds using Athena, a new higher-order Godunov code. We analyze the dependence of the saturation amplitude, energy dissipation characteristics, power spectra, sonic scaling, and indicators of intermittency in the turbulence on factors such as the magnetic field strength, driving scale, energy injection rate, and numerical resolution. While convergence in the energies is reached at moderate resolutions, we find that the power spectra require much higher resolutions that are difficult to obtain. In a 1024^3 hydro run, we find a power law relationship between the velocity dispersion and the spatial scale on which it is measured, while for an MHD run at the same resolution we find no such power law. The time-variability and temperature intermittency in the turbulence both show a dependence on the driving scale, indicating that numerically driving turbulence by an arbitrary mechanism may not allow a realistic representation of these properties. We also note similar features in the power spectrum of the compressive component of velocity for supersonic MHD turbulence as in the velocity spectrum of an initially-spherical MHD blast wave, implying that the power law form does not rule out shocks, rather than a turbulent cascade, playing a significant role in the regulation of energy transfer between spatial scales.

We study the effects of anisotropic thermal conduction on magnetized spherical accretion flows using global axisymmetric MHD simulations. In low collisionality plasmas, the Bondi spherical accretion solution is unstable to the magnetothermal instability (MTI). The MTI grows rapidly at large radii where the inflow is subsonic. For a weak initial field, the MTI saturates by creating a primarily radial magnetic field, i.e., by aligning the field lines with the background temperature gradient. The saturation is quasilinear in the sense that the magnetic field is amplified by a factor of $\sim 10-30$ independent of the initial field strength (for weak fields). In the saturated state, the conductive heat flux is much larger than the convective heat flux, and is comparable to the field-free (Spitzer) value (since the field lines are largely radial). The MTI by itself does not appreciably change the accretion rate $\dot M$ relative to the Bondi rate $\dot M_B$. However, the radial field lines created by the MTI are amplified by flux freezing as the plasma flows in to small radii. Oppositely directed field lines are brought together by the converging inflow, leading to significant resistive heating. When the magnetic energy density is comparable to the gravitational potential energy density, the plasma is heated to roughly the virial temperature; the mean inflow is highly subsonic; most of the energy released by accretion is transported to large radii by thermal conduction; and the accretion rate $\dot M \ll \dot M_B$. The predominantly radial magnetic field created by the MTI at large radii in spherical accretion flows may account for the stable Faraday rotation measure towards Sgr A* in the Galactic Center.

The electron mean free path in the intracluster medium (ICM) of galaxy clusters is much larger than the gyroradius; thus, heat is transported anisotropically along magnetic field lines. We show that the intracluster medium is unstable to the magnetothermal instability (MTI) using MHD simulations with anisotropic thermal conduction. As a result of the MTI, we find that the temperature profile of the ICM can be substantially modified on timescales of several billion years while the magnetic field is amplified by dynamo action up to more than fifty times the original energy. We also show that the instability drives field lines to become preferentially radial leading to conduction that is a highly efficient fraction of the Spitzer conductivity. As such, we present the first self-consistent calculation of the effective thermal conductivity in the ICM.

A new code for astrophysical magnetohydrodynamics (MHD) is described. The code has been designed to be easily extensible for use with static and adaptive mesh refinement. It combines higher-order Godunov methods with the constrained transport (CT) technique to enforce the divergence-free constraint on the magnetic field. Discretization is based on cell-centered volume-averages for mass, momentum, and energy, and face-centered area-averages for the magnetic field. Novel features of the algorithm include (1) a consistent framework for computing the time- and edge-averaged electric fields used by CT to evolve the magnetic field from the time- and area-averaged Godunov fluxes, (2) the extension to MHD of spatial reconstruction schemes that involve a dimensionally-split time advance, and (3) the extension to MHD of two different dimensionally-unsplit integration methods. Implementation of the algorithm in both C and Fortran95 is detailed, including strategies for parallelization using domain decomposition. Results from a test suite which includes problems in one-, two-, and three-dimensions for both hydrodynamics and MHD are given, not only to demonstrate the fidelity of the algorithms, but also to enable comparisons to other methods. The source code is freely available for download on the web.

The magnetohydrodynamic evolution of a dense spherical cloud as it interacts with a strong planar shock is studied, as a model for shock interactions with density inhomogeneities in the interstellar medium. The cloud is assumed to be small enough that radiative cooling, thermal conduction, and self-gravity can be ignored. A variety of initial orientations (including parallel, perpendicular, and oblique to the incident shock normal) and strengths for the magnetic field are investigated. During the early stages of the interaction (less than twice the time taken for the transmitted shock to cross the interior of the cloud) the structure and dynamics of the shocked cloud is fairly insensitive to the magnetic field strength and orientation. However, at late times strong fields substantially alter the dynamics of the cloud, suppressing fragmentation and mixing by stabilizing the interface at the cloud surface. Even weak magnetic fields can drastically alter the evolution of the cloud compared to the hydrodynamic case. Weak fields of different geometries result in different distributions and amplifications of the magnetic energy density, which may affect the thermal and non-thermal x-ray emission expected from shocked clouds associated with, for example, supernovae remnants.

We perform numerical simulations of solid particle motion in a shearing box model of a protoplanetary disc. The accretion flow is turbulent due to the action of the magnetorotational instability. Aerodynamic drag on the particles is modelled using the Epstein law with the gas velocity interpolated to the particle position. The effect of the magnetohydrodynamic turbulence on particle velocity dispersions is quantified for solids of different stopping times t_s, or equivalently, different sizes. The anisotropy of the turbulence is reflected upon the dispersions of the particle velocity components, with the radial component larger than both the azimuthal and vertical components for particles larger than ~ 10 cm (assuming minimum-mass solar nebula conditions at 5 AU). The dispersion of the particle velocity magnitude, as well as that of the radial and azimuthal components, as functions of stopping time, agree with previous analytical results for isotropic turbulence. The relative speed between pairs of particles with the same value of t_s decays faster with decreasing separation than in the case of solids with different stopping time. Correlations in the particle number density introduce a non-uniform spatial distribution of solids in the 10 to 100 cm size range. Any clump of particles is disrupted by the turbulence in less than one tenth on an orbital period, and the maximally concentrated clumps are stable against self-gravitational collapse.

We present a single step, second-order accurate Godunov scheme for ideal MHD which is an extension of the method described by Gardiner & Stone (2005) to three dimensions. This algorithm combines the corner transport upwind (CTU) method of Colella for multidimensional integration, and the constrained transport (CT) algorithm for preserving the divergence-free constraint on the magnetic field. We describe the calculation of the PPM interface states for 3D ideal MHD which must include multidimensional ``MHD source terms'' and naturally respect the balance implicit in these terms by the ${\bf\nabla\cdot B}=0$ condition. We compare two different forms for the CTU integration algorithm which require either 6- or 12-solutions of the Riemann problem per cell per time-step, and present a detailed description of the 6-solve algorithm. Finally, we present solutions for test problems to demonstrate the accuracy and robustness of the algorithm.

The New Horizons ALICE instrument is a lightweight (4.4 kg), low-power (4.4 Watt) imaging spectrograph aboard the New Horizons mission to Pluto/Charon and the Kuiper Belt. Its primary job is to determine the relative abundances of various species in Pluto's atmosphere. ALICE will also be used to search for an atmosphere around Pluto's moon, Charon, as well as the Kuiper Belt Objects (KBOs) that New Horizons hopes to fly by after Pluto-Charon, and it will make UV surface reflectivity measurements of all of these bodies as well. The instrument incorporates an off-axis telescope feeding a Rowland-circle spectrograph with a 520-1870 angstroms spectral passband, a spectral point spread function of 3-6 angstroms FWHM, and an instantaneous spatial field-of-view that is 6 degrees long. Different input apertures that feed the telescope allow for both airglow and solar occultation observations during the mission. The focal plane detector is an imaging microchannel plate (MCP) double delay-line detector with dual solar-blind opaque photocathodes (KBr and CsI) and a focal surface that matches the instrument's 15-cm diameter Rowland-circle. In what follows, we describe the instrument in greater detail, including descriptions of its ground calibration and initial in flight performance.

The New Horizons instrument named Ralph is a visible/near infrared multi-spectral imager and a short wavelength infrared spectral imager. It is one of the core instruments on New Horizons, NASA's first mission to the Pluto/Charon system and the Kuiper Belt. Ralph combines panchromatic and color imaging capabilities with IR imaging spectroscopy. Its primary purpose is to map the surface geology and composition of these objects, but it will also be used for atmospheric studies and to map the surface temperature. It is a compact, low-mass (10.5 kg), power efficient (7.1 W peak), and robust instrument with good sensitivity and excellent imaging characteristics. Other than a door opened once in flight, it has no moving parts. These characteristics and its high degree of redundancy make Ralph ideally suited to this long-duration flyby reconnaissance mission.

These Faraday rotation calculations of hot, thick accretion flows are motivated by the measured steady rotation measure (RM) of $\approx -6 \times 10^5$ rad m$^{-2}$ from Sgr A*. In our numerical simulations, the quasi-steady state structure of the accretion flow, and the RM it produces, depends on the initial magnetic field. In spite of this dependence, we can draw several robust conclusions about Faraday rotation produced by geometrically thick accretion disks: i) the time averaged RM does not depend that sensitively on the viewing angle, but the stability of the RM can. Equatorial viewing angles show significant variability in RM (including sign reversals), while polar viewing angles are relatively stable if there is a large scale magnetic field threading the disk at large radii. ii) Most of the RM is produced at small radii for polar viewing angles while all radii contribute significantly near the midplane of the disk. Our simulations confirm previous analytic arguments that the accretion rate onto Sgr A* must satisfy $\dot M_{\rm in} \ll \dot M_{\rm Bondi} \sim 10^{-5} \mpy$ in order to not over-produce the measured RM. We argue that the steady RM $\approx -6 \times 10^5$ rad m$^{-2}$ from Sgr A* has two plausible explanations: 1) it is produced at $\sim 100$ Schwarzschild radii, requires $\dot{M}_{\rm in} \approx 3 \times 10^{-8} M_\odot$ yr$^{-1}$, and we view the flow at an angle of $\sim 30^\circ$ relative to the rotation axis of the disk; in our simulations, the variation in RM across a finite-sized source is sufficient to depolarize the emission below $\approx$ 100 GHz, consistent with observations. 2) Alternatively, the RM may be produced in the relatively spherical inflowing plasma near the circularization radius at $\sim 10^3-10^4$ Schwarzschild radii.

Local (shearing box) simulations of the nonlinear evolution of the magnetorotational instability in a collisionless plasma show that angular momentum transport by pressure anisotropy ($p_\perp \ne p_\parallel$, where the directions are defined with respect to the local magnetic field) is comparable to that due to the Maxwell and Reynolds stresses. Pressure anisotropy, which is effectively a large-scale viscosity, arises because of adiabatic invariants related to $p_\perp$ and $p_\parallel$ in a fluctuating magnetic field. In a collisionless plasma, the magnitude of the pressure anisotropy, and thus the viscosity, is determined by kinetic instabilities at the cyclotron frequency. Our simulations show that $\sim 50$ % of the gravitational potential energy is directly converted into heat at large scales by the viscous stress (the remaining energy is lost to grid-scale numerical dissipation of kinetic and magnetic energy). We show that electrons receive a significant fraction ($\sim [T_e/T_i]^{1/2}$) of this dissipated energy. Employing this heating by an anisotropic viscous stress in one dimensional models of radiatively inefficient accretion flows, we find that the radiative efficiency of the flow is greater than 0.5% for $\dot{M} \gtrsim 10^{-4} \dot{M}_{Edd}$. Thus a low accretion rate, rather than just a low radiative efficiency, is necessary to explain the low luminosity of many accreting black holes. For Sgr A* in the Galactic Center, our predicted radiative efficiencies imply an accretion rate of $\approx 3 \times 10^{-8} M_\odot {\rm yr^{-1}}$ and an electron temperature of $\approx 3 \times 10^{10}$ K at $\approx 10$ Schwarzschild radii; the latter is consistent with the brightness temperature inferred from VLBI observations.

We study the magnetic Rayleigh-Taylor instability in three dimensions, with focus on the nonlinear structure and evolution that results from different initial field configurations. We study strong fields in the sense that the critical wavelength l_c at which perturbations along the field are stable is a large fraction of the size of the computational domain. We consider magnetic fields which are initially parallel to the interface, but have a variety of configurations, including uniform everywhere, uniform in the light fluid only, and fields which change direction at the interface. Strong magnetic fields do not suppress instability, in fact by inhibiting secondary shear instabilities, they reduce mixing between the heavy and light fluid, and cause the rate of growth of bubbles and fingers to increase in comparison to hydrodynamics. Fields parallel to, but whose direction changes at, the interface produce long, isolated fingers separated by the critical wavelength l_c, which may be relevant to the morphology of the optical filaments in the Crab nebula.

First results from a fully self-consistent, temperature-dependent equation of state that spans the whole density range of neutron stars and supernova cores are presented. The equation of state (EoS) is calculated using a mean-field Hartree-Fock method in three dimensions (3D). The nuclear interaction is represented by the phenomenological Skyrme model in this work, but the EoS can be obtained in our framework for any suitable form of the nucleon-nucleon effective interaction. The scheme we employ naturally allows effects such as (i) neutron drip, which results in an external neutron gas, (ii) the variety of exotic nuclear shapes expected for extremely neutron heavy nuclei, and (iii) the subsequent dissolution of these nuclei into nuclear matter. In this way, the equation of state is calculated across phase transitions without recourse to interpolation techniques between density regimes described by different physical models. EoS tables are calculated in the wide range of densities, temperature and proton/neutron ratios on the ORNL NCCS XT3, using up to 2000 processors simultaneously.

We present a method for simulating the evolution of HII regions driven by point sources of ionizing radiation in magnetohydrodynamic media, implemented in the three-dimensional Athena MHD code. We compare simulations using our algorithm to analytic solutions and show that the method passes rigorous tests of accuracy and convergence. The tests reveal several conditions that an ionizing radiation-hydrodynamic code must satisfy to reproduce analytic solutions. As a demonstration of our new method, we present the first three-dimensional, global simulation of an HII region expanding into a magnetized gas. The simulation shows that magnetic fields suppress sweeping up of gas perpendicular to magnetic field lines, leading to small density contrasts and extremely weak shocks at the leading edge of the HII region's expanding shell.

We study the nonlinear evolution of the magnetic Rayleigh-Taylor instability using three-dimensional MHD simulations. We consider the idealized case of two inviscid, perfectly conducting fluids of constant density separated by a contact discontinuity perpendicular to the effective gravity g, with a uniform magnetic field B parallel to the interface. Modes parallel to the field with wavelengths smaller than l_c = [B B/(d_h - d_l) g] are suppressed (where d_h and d_l are the densities of the heavy and light fluids respectively), whereas modes perpendicular to B are unaffected. We study strong fields with l_c varying between 0.01 and 0.36 of the horizontal extent of the computational domain. Even a weak field produces tension forces on small scales that are significant enough to reduce shear (as measured by the distribution of the amplitude of vorticity), which in turn reduces the mixing between fluids, and increases the rate at which bubbles and finger are displaced from the interface compared to the purely hydrodynamic case. For strong fields, the highly anisotropic nature of unstable modes produces ropes and filaments. However, at late time flow along field lines produces large scale bubbles. The kinetic and magnetic energies transverse to gravity remain in rough equipartition and increase as t^4 at early times. The growth deviates from this form once the magnetic energy in the vertical field becomes larger than the energy in the initial field. We comment on the implications of our results to Z-pinch experiments, and a variety of astrophysical systems.

We perform a linear magnetohydrodynamic perturbation analysis for a stratified magnetized envelope where the diffusion of heat is mediated by charged particles that are confined to flow along magnetic field lines. We identify an instability, the ``coulomb bubble instability,'' which may be thought of as standard magnetosonic fast and slow waves, driven by the rapid diffusion of heat along the direction of the magnetic field. We calculate the growth rate and stability criteria for the coulomb bubble instability for various choices of equilibrium conditions. The coulomb bubble instability is most strongly driven for weakly magnetized atmospheres that are strongly convectively stable. We briefly discuss a possible application of astrophysical interest: diffusion of interstellar cosmic rays in the hot T ~ 10^6 K Galactic corona. We show that for commonly accepted values of the cosmic ray and gas pressure as well as its overall characteristic dimensions, the Galactic corona is in a marginal state of stability with respect to a cosmic ray coulomb bubble instability. The implication being that a cosmic ray coulomb bubble instability plays a role regulating both the pressure and transport properties of interstellar cosmic rays, while serving as a source of acoustic power above the galactic disk. (abridged)

Using fully three-dimensional hydrodynamic simulations, we investigate the effect of the Coriolis force on the hydrodynamic and observable properties of colliding wind binary systems. To make the calculations tractable, we assume adiabatic, constant velocity winds. The neglect of radiative driving, gravitational deceleration, and cooling limit the application of our models to real systems. However, these assumptions allow us to isolate the effect of the Coriolis force, and by simplifying the calculations, allow us to use a higher resolution (up to 640^3) and to conduct a larger survey of parameter space. We study the dynamics of collidng winds with equal mass loss rates and velocities emanating from equal mass stars on circular orbits, with a range of values for the ratio of the wind to orbital velocity. We also study the dynamics of winds from stars on elliptical orbits and with unequal strength winds. Orbital motion of the stars sweeps the shocked wind gas into an Archimedean spiral, with asymmetric shock strengths and therefore unequal postshock temperatures and densities in the leading and trailing edges of the spiral. We observe the Kelvin-Helmholtz instability at the contact surface between the shocked winds in systems with orbital motion even when the winds are identical. The change in shock strengths caused by orbital motion increases the volume of X-ray emitting post-shock gas with T > 0.59 keV by 63% for a typical system as the ratio of wind velocity to orbital velocity decreases to V_w/V_o = 2.5. This causes increased free-free emission from systems with shorter orbital periods and an altered time-dependence of the wind attenuation. We comment on the importance of the effects of orbital motion on the observable properties of colliding wind binaries.

In dilute astrophysical plasmas, thermal conduction is primarily along magnetic field lines, and therefore highly anisotropic. As a result, the usual convective stability criterion is modified from a condition on entropy to a condition on temperature. For small magnetic fields or small wavenumbers, instability occurs in any atmosphere where the temperature and pressure gradients point in the same direction. We refer to the resulting convective instability as the magnetothermal instability (MTI). We present fully three-dimensional simulations of the MTI and show that saturation results in an atmosphere with different vertical structure, dependent upon the boundary conditions. When the temperature at the boundary of the unstable layer is allowed to vary, the temperature gradient relaxes until the unstable region is almost isothermal. When the temperature at the boundary of the unstable region is fixed, the magnetic field is reoriented to an almost vertical geometry as a result of buoyant motions. This case exhibits more vigorous turbulence. In both cases the resulting saturated heat flux is almost one-half of the value expected if the conduction were purely isotropic. The action of the MTI results in dynamical processes that lead to significant transport perpendicular to the initial direction of the magnetic field. The resulting magnetoconvection in both cases amplifies the magnetic field until it is almost in equipartition with sustained subsonic turbulence. These results are relevant to understanding measurements of the temperature profiles of the intracluster medium of clusters of galaxies as well as the structure of radiatively inefficient accretion flows.

We carry out three-dimensional, high resolution (up to $1024^2\times 256$) hydrodynamic simulations of the evolution of vortices in vertically unstratified Keplerian disks using the shearing sheet approximation. The transient amplification of incompressible, linear amplitude leading waves (which has been proposed as a possible route to nonlinear hydrodynamical turbulence in disks) is used as one test of our algorithms; our methods accurately capture the predicted amplification, converges at second-order, and is free from aliasing. Waves expected to reach nonlinear amplitude at peak amplification become unstable to Kelvin-Helmholtz modes when $\mid W_{\rm max}\mid\gtrsim \Omega$ (where $W_{\rm max}$ is the local maximum of vorticity and $\Omega$ the angular velocity). We study the evolution of a power-law distribution of vorticity consistent with Kolmogorov turbulence; in two-dimensions long-lived vortices emerge and decay slowly, similar to previous studies. In three-dimensions, however, vortices are unstable to bending modes, leading to rapid decay. Only vortices with a length to width ratio smaller than one survive; in three-dimensions the residual kinetic energy and shear stress is at least one order of magnitude smaller than in two-dimensions. No evidence for sustained hydrodynamical turbulence and transport is observed in three-dimensions. Instead, at late times the residual transport is determined by the amplitude of slowly decaying, large-scale vortices (with horizontal extent comparable to the scale height of the disk), with additional contributions from nearly incompressible inertial waves possible. Evaluating the role that large-scale vortices play in astrophysical accretion disks will require understanding the mechanisms that generate and destroy them.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

In many magnetized, dilute astrophysical plasmas, thermal conduction occurs almost exclusively parallel to magnetic field lines. In this case, the usual stability criterion for convective stability, the Schwarzschild criterion, which depends on entropy gradients, is modified. In the magnetized long mean free path regime, instability occurs for small wavenumbers when (dP/dz)(dln T/dz) > 0, which we refer to as the Balbus criterion. We refer to the convective-type instability that results as the magnetothermal instability (MTI). We use the equations of MHD with anisotropic electron heat conduction to numerically simulate the linear growth and nonlinear saturation of the MTI in plane-parallel atmospheres that are unstable according to the Balbus criterion. The linear growth rates measured from the simulations are in excellent agreement with the weak field dispersion relation. The addition of isotropic conduction, e.g. radiation, or strong magnetic fields can damp the growth of the MTI and affect the nonlinear regime. The instability saturates when the atmosphere becomes isothermal as the source of free energy is exhausted. By maintaining a fixed temperature difference between the top and bottom boundaries of the simulation domain, sustained convective turbulence can be driven. MTI-stable layers introduced by isotropic conduction are used to prevent the formation of unresolved, thermal boundary layers. We find that the largest component of the time-averaged heat flux is due to advective motions as opposed to the actual thermal conduction itself. Finally, we explore the implications of this instability for a variety of astrophysical systems, such as neutron stars, the hot intracluster medium of galaxy clusters, and the structure of radiatively inefficient accretion flows.

We describe the design, performance and scientific objectives of the NASA-funded ALICE instrument aboard the ESA Rosetta asteroid flyby/comet rendezvous mission. ALICE is a lightweight, low-power, and low-cost imaging spectrograph optimized for cometary far-ultraviolet (FUV) spectroscopy. It will be the first UV spectrograph to study a comet at close range. It is designed to obtain spatially-resolved spectra of Rosetta mission targets in the 700-2050 A spectral band with a spectral resolution between 8 A and 12 A for extended sources that fill its ~0.05 deg x 6.0 deg field-of-view. ALICE employs an off-axis telescope feeding a 0.15-m normal incidence Rowland circle spectrograph with a concave holographic reflection grating. The imaging microchannel plate detector utilizes dual solar-blind opaque photocathodes (KBr and CsI) and employs a 2 D delay-line readout array. The instrument is controlled by an internal microprocessor. During the prime Rosetta mission, ALICE will characterize comet 67P/Churyumov-Gerasimenko's coma, its nucleus, and the nucleus/coma coupling; during cruise to the comet, ALICE will make observations of the mission's two asteroid flyby targets and of Mars, its moons, and of Earth's moon. ALICE has already successfully completed the in-flight commissioning phase and is operating normally in flight. It has been characterized in flight with stellar flux calibrations, observations of the Moon during the first Earth fly-by, and observations of comet Linear T7 in 2004 and comet 9P/Tempel 1 during the 2005 Deep Impact comet-collision observing campaign

Stellar atmosphere models of ionized accretion disks have generally neglected the contribution of magnetic fields to the vertical hydrostatic support, although magnetic fields are widely believed to play a critical role in the transport of angular momentum. Simulations of magnetorotational turbulence in a vertically stratified shearing box geometry show that magnetic pressure support can be dominant in the upper layers of the disk. We present calculations of accretion disk spectra that include this magnetic pressure support, as well as a vertical dissipation profile based on simulation. Magnetic pressure support generically produces a more vertically extended disk atmosphere with a larger density scale height. This acts to harden the spectrum compared to models that neglect magnetic pressure support. We estimate the significance of this effect on disk-integrated spectra by calculating an illustrative disk model for a stellar mass black hole, assuming that similar magnetic pressure support exists at all radii.

(shortened) We calculate the vertical structure of a local patch of an accretion disk in which heating by dissipation of MRI-driven MHD turbulence is balanced by radiative cooling. Heating, radiative transport, and cooling are computed self-consistently with the structure by solving the equations of radiation MHD in the shearing-box approximation. Using a fully 3-d and energy-conserving code, we compute the structure of this disk segment over a span of more than five cooling times. After a brief relaxation period, a statistically steady-state develops. Measuring height above the midplane in units of the scale-height H predicted by a Shakura-Sunyaev model, we find that the disk atmosphere stretches upward, with the photosphere rising to about 7H, in contrast to the approximately 3H predicted by conventional analytic models. This more extended structure, as well as fluctuations in the height of the photosphere, may lead to departures from Planckian form in the emergent spectra. Dissipation is distributed across the region within roughly 3H of the midplane, but is very weak at greater altitudes. Because fluctuations in the dissipation are particularly strong away from the midplane, the emergent radiation flux can track dissipation fluctuations with a lag that is only 0.1--0.2 times the mean cooling time of the disk. Long timescale asymmetries in the dissipation distribution can also cause significant asymmetry in the flux emerging from the top and bottom surfaces of the disk. Radiative diffusion dominates Poynting flux in the vertical energy flow throughout the disk.

We describe local shearing box simulations of turbulence driven by the magnetorotational instability (MRI) in a collisionless plasma. Collisionless effects may be important in radiatively inefficient accretion flows, such as near the black hole in the Galactic Center. The MHD version of ZEUS is modified to evolve an anisotropic pressure tensor. A fluid closure approximation is used to calculate heat conduction along magnetic field lines. The anisotropic pressure tensor provides a qualitatively new mechanism for transporting angular momentum in accretion flows (in addition to the Maxwell and Reynolds stresses). We estimate limits on the pressure anisotropy due to pitch angle scattering by kinetic instabilities. Such instabilities provide an effective ``collision'' rate in a collisionless plasma and lead to more MHD-like dynamics. We find that the MRI leads to efficient growth of the magnetic field in a collisionless plasma, with saturation amplitudes comparable to those in MHD. In the saturated state, the anisotropic stress is comparable to the Maxwell stress, implying that the rate of angular momentum transport may be moderately enhanced in a collisionless plasma.

In weakly magnetized, dilute plasmas in which thermal conduction along magnetic field lines is important, the usual convective stability criterion is modified. Instead of depending on entropy gradients, instability occurs for small wavenumbers when (dP/dz)(dln T/dz) > 0, which we refer to as the Balbus criterion. We refer to the convective instability that results in this regime as the magnetothermal instability (MTI). We use numerical MHD simulations which include anisotropic electron heat conduction to follow the growth and saturation of the MTI in two-dimensional, plane parallel atmospheres that are unstable according to the Balbus criterion. The linear growth rates measured in the simulations agree with the weak field dispersion relation. We investigate the effect of strong fields and isotropic conduction on the linear properties and nonlinear regime of the MTI. In the nonlinear regime, the instability saturates and convection decays away, when the atmosphere becomes isothermal. Sustained convective turbulence can be driven if there is a fixed temperature difference between the top and bottom edges of the simulation domain, and if isotropic conduction is used to create convectively stable layers that prevent the formation of unresolved, thermal boundary layers. The largest component of the time-averaged heat flux is due to advective motions. These results have implications for a variety of astrophysical systems, such as the temperature profile of hot gas in galaxy clusters, and the structure of radiatively inefficient accretion flows.

The double pulsar J0737--3039 has become an important astrophysical laboratory for testing fundamental physics. Here we demonstrate that the low measured mass of Pulsar B can be used to constrain the equation of state of neutron star matter {\em under the assumption} that it formed in an electron-capture supernova. We show that the observed orbital parameters as well as the likely evolutionary history of the system support such a hypothesis and discuss future refinements that will improve the constraints this test may provide.

In these proceedings we present recent efforts to understand the energetics of magnetohydrodynamic (MHD) turbulence driven by the magneto-rotational instability (MRI). These studies are carried out in the local (shearing box) approximation using the Athena simulation code. Athena is a higher order Godunov algorithm based on the piecewise parabolic method (PPM), the corner transport upwind (CTU) integration algorithm, and the constrained transport (CT) algorithm for evolving the magnetic field. This algorithm is particularly suited for these studies owing to the conservation properties of a Godunov scheme and the particular implementation of the shearing box source terms used here. We present a variety of calculations which may be compared directly to previously published results and discuss them in some detail. The only significant discrepancy found between the results presented here and in the published literature involves the turbulent heating rate. We observe the presence of recurrent channel solutions in calculations involving a mean vertical magnetic field and the associated time lag between the energy injection and thermalization rate. We also present the results of a shearing box calculation which includes an optically thin radiative term with a cooling rate selected to match the turbulent heating rate. Some properties which we find uniformly present in all of the calculations presented here are compressible fluctuations, spiral waves and weak shocks. It is found that these compressible modes dominate the temporal fluctuations in the probability distribution functions for most of the thermodynamic variables; only the specific entropy is relatively immune to their effects.

We describe a new numerical scheme for MHD which combines a higher order Godunov method (PPM) with Constrained Transport. The results from a selection of multidimensional test problems are presented. The complete test suite used to validate the method, as well as implementations of the algorithm in both F90 and C, are available from the web. A fully three-dimensional version of the algorithm has been developed, and is being applied to a variety of astrophysical problems including the decay of supersonic MHD turbulence, the nonlinear evolution of the MHD Rayleigh-Taylor instability, and the saturation of the magnetorotational instability in the shearing box. Our new simulations of the MRI represent the first time that a higher-order Godunov scheme has been applied to this problem, providing a quantitative check on the accuracy of previous results computed with ZEUS; the latter are found to be reliable.

We describe a single step, second-order accurate Godunov scheme for ideal MHD based on combining the piecewise parabolic method (PPM) for performing spatial reconstruction, the corner transport upwind (CTU) method of Colella for multidimensional integration, and the constrained transport (CT) algorithm for preserving the divergence-free constraint on the magnetic field. We adopt the most compact form of CT, which requires the field be represented by area-averages at cell faces. We demonstrate that the fluxes of the area-averaged field used by CT can be made consistent with the fluxes of the volume-averaged field returned by a Riemann solver if they obey certain simple relationships. We use these relationships to derive new algorithms for constructing the CT fluxes at grid cell corners which reduce exactly to the equivalent one-dimensional solver for plane-parallel, grid-aligned flow. We show that the PPM reconstruction algorithm must include multidimensional terms for MHD, and we describe a number of important extensions that must be made to CTU in order for it to be used for MHD with CT. We present the results of a variety of test problems to demonstrate the method is accurate and robust.

We calculate the radial diffusion coefficient for a passive contaminant in an accretion disc which is turbulent due to the action of the magnetorotational instability. Numerical MHD simulations are used to follow the evolution of a local patch of the disc using the shearing box formalism. A separate continuity equation for the mass fraction of contaminant is integrated along with the MHD system, and radial profiles of this fraction are obtained as a function of time. Solutions of a linear diffusion equation are fitted to the numerical measured profiles of the contaminant, treating the diffusion coefficient D as the fitting parameter. At early times, the value of D is found to vary, however once the contaminant is spread over scales comparable to the box size, it saturates at a steady value. The ratio of D to the transport coefficient of angular momentum due to shear stress is small. If D can be used as a proxy for the turbulent magnetic diffusivity, the effective magnetic Prandtl number P_eff=\nu/D (where \nu is the coefficient of "efective viscosity" associated with shear stress) would be large.

(Abridged) Ideal hydrodynamic models of the intracluster medium (ICM) in the core regions of galaxy clusters fail to explain both the observed temperature structure of this gas, and the observed morphology of radio-galaxy/ICM interactions. It has recently been suggested that, even in the presence of reasonable magnetic fields, thermal conduction in the ICM may be crucial for reproducing the temperature floor seen in many systems. If this is indeed correct, it raises the possibility that other transport processes may be important. With this motivation, we present a numerical investigation of the buoyant evolution of AGN-blown cavities in ICM that has a non-negligible shear viscosity. We use the ZEUS-MP code to follow the 3-d evolution of an initially static, hot bubble in a beta-model ICM atmosphere with varying degrees of shear viscosity. With no explicit viscosity, it is found that the combined action of Rayleigh-Taylor and Kelvin-Helmholtz instabilities shred the ICM cavity and one does not reproduce the intact and detached ``ghost cavities'' observed in systems such as Perseus-A. On the other hand, even a modest level of shear viscosity can be important in quenching the fluid instabilities and maintaining the integrity of the bubble. In particular, we show that the morphology of the NW ghost cavity found in Perseus-A can be reproduced, as can the flow pattern inferred from the morphology of H-alpha filaments. Finally, we discuss the possible relevance of ICM viscosity to the fact that many of the active ICM cavities are not bounded by strong shocks.

We present the results of indirect searches for Weakly Interacting Massive Particles (WIMPs) with 1679.6 live days of data from the Super-Kamiokande detector using neutrino-induced upward through-going muons. The search is performed by looking for an excess of high energy muon neutrinos from WIMP annihilations in the Sun, the core of the Earth, and the Galactic Center, as compared to the number expected from the atmospheric neutrino background. No statistically significant excess was seen. We calculate flux limits in various angular cones around each of the above celestial objects. We obtain conservative model-independent upper limits on WIMP-nucleon cross-section as a function of WIMP mass and compare these results with the corresponding results from direct dark matter detection experiments.

The saturation level of the magnetorotational instability (MRI) is investigated using three-dimensional MHD simulations. The shearing box approximation is adopted and the vertical component of gravity is ignored, so that the evolution of the MRI is followed in a small local part of the disk. We focus on the dependence of the saturation level of the stress on the gas pressure, which is a key assumption in the standard alpha disk model. From our numerical experiments it is found that there is a weak power-law relation between the saturation level of the Maxwell stress and the gas pressure in the nonlinear regime; the higher the gas pressure, the larger the stress. Although the power-law index depends slightly on the initial field geometry, the relationship between stress and gas pressure is independent of the initial field strength, and is unaffected by Ohmic dissipation if the magnetic Reynolds number is at least 10. The relationship is the same in adiabatic calculations, where pressure increases over time, and nearly-isothermal calculations, where pressure varies little with time. Our numerical results are qualitatively consistent with an idea that the saturation level of the MRI is determined by a balance between the growth of the MRI and the dissipation of the field through reconnection. The quantitative interpretation of the pressure-stress relation, however, may require advances in the theoretical understanding of non-steady magnetic reconnection.

Using local 3D MHD simulations, we investigate ways in which galactic turbulence associated with the magnetorotational instability (MRI) may influence the formation and properties of GMCs. Our disk models are vertically stratified and subject to uniform shear. Initial magnetic fields are weak and purely vertical. For simplicity, we adopt an isothermal equation of state. We find that MRI-driven turbulence develops rapidly, with the saturated-state Shakura & Sunyaev parameter alpha~(0.15-0.3) dominated by Maxwell stresses. Many of the dimensionless characteristics of the turbulence (e.g. the ratio of the Maxwell to Reynolds stresses) are similar to results from previous MRI studies of accretion disks, hence insensitive to the degree of vertical disk compression, shear rate, and the presence of self-gravity. The density-weighted velocity dispersions in non- or weakly self-gravitating disks are sigma_x ~ sigma_y ~ (0.4-0.6)c_s and sigma_z~(0.2-0.3)c_s, suggesting that MRI can contribute significantly to the observed level of galactic turbulence. The saturated-state magnetic field strength B ~ 2 \mu G is similar to typical galactic values. When self-gravity is strong enough, MRI-driven high-amplitude density perturbations are swing-amplified to form Jeans-mass (~10^7 Msun) bound clouds. Compared to previous unmagnetized or strongly-magnetized disk models, the threshold for nonlinear instability in the present models occurs for surface densities at least 50% lower, corresponding to the Toomre parameter Q_th~1.6. We present evidence that self-gravitating clouds like GMCs formed under conditions similar to our models can lose much of their original spin angular momenta by magnetic braking, preferentially via fields threading near-perpendicularly to their spin axes.

In the nuclear and particle astrophysics session of CIPANP 2003 we heard talks on a number of topics, focused for the most part into four broad areas. Here we outline the discussions of the standard cosmological model, dark matter searches, cosmic rays, and neutrino astrophysics. The robustness of theoretical and experimental programs in all of these areas is very encouraging, and we expect to have many questions answered, and new ones asked, in time for CIPANP 2006.

Observations reveal concentrations of molecular line emission on the sky, called ``clumps,'' in dense, star-forming molecular clouds. These clumps are believed to be the eventual sites of star formation. We study the three-dimensional analogs of clumps using a set of self-consistent, time-dependent numerical models of molecular clouds. The models follow the decay of initially supersonic turbulence in an isothermal, self-gravitating, magnetized fluid. We find the following. (1) Clumps are intrinsically triaxial. This explains the observed deficit of clumps with a projected axis ratio near unity, and the apparent prolateness of clumps. (2) Simulated clump axes are not strongly aligned with the mean magnetic field within clumps, nor with the large-scale mean fields. This is in agreement with observations. (3) The clump mass spectrum has a high-mass slope that is consistent with the Salpeter value. There is a low-mass break in the slope at $\sim 0.5 \msun$, although this may depend on model parameters including numerical resolution. (4) The typical specific spin angular momentum of clumps is $4 \times 10^{22} {\rm cm^2 s^{-1}}$. This is larger than the median specific angular momentum of binary stars. Scaling arguments suggest that higher resolution simulations may soon be able to resolve the scales at which the angular momentum of binary stars is determined.

We examine the small-scale dynamics of black hole accretion disks in which radiation pressure exceeds gas pressure. Local patches of disk are modeled by numerically integrating the equations of radiation MHD in the flux-limited diffusion approximation. The shearing-box approximation is used, and the vertical component of gravity is neglected. Magneto-rotational instability (MRI) leads to turbulence in which accretion stresses are due primarily to magnetic torques. When radiation is locked to gas over the length and time scales of fluctuations in the turbulence, the accretion stress, density contrast, and dissipation differ little from those in the corresponding calculations with radiation replaced by extra gas pressure. However, when radiation diffuses each orbit a distance that is comparable to the RMS vertical wavelength of the MRI, radiation pressure is less effective in resisting squeezing. Large density fluctuations occur, and radiation damping of compressive motions converts PdV work into photon energy. The accretion stress in calculations having a net vertical magnetic field is found to be independent of opacity over the range explored, and approximately proportional to the square of the net field. In calculations with zero net magnetic flux, the accretion stress depends on the portion of the total pressure that is effective in resisting compression. The stress is lower when radiation diffuses rapidly with respect to the gas. We show that radiation-supported Shakura-Sunyaev disks accreting via internal magnetic stresses are likely in their interiors to have radiation marginally coupled to turbulent gas motions.

We analyze the spectral properties of driven, supersonic compressible magnetohydrodynamic (MHD) turbulence obtained via high-resolution numerical experiments, for application to understanding the dynamics of giant molecular clouds. Via angle-averaged power spectra, we characterize the transfer of energy from the intermediate, driving scales down to smaller dissipative scales, and also present evidence for inverse cascades that achieve modal-equipartition levels on larger spatial scales. Investigating compressive versus shear modes separately, we evaluate their relative total power, and find that as the magnetic field strength decreases, (1) the shear fraction of the total kinetic power decreases, and (2) slopes of power-law fits over the inertial range steepen. To relate to previous work on incompressible MHD turbulence, we present qualitative and quantitative measures of the scale-dependent spectral anisotropy arising from the shear-Alfv\'{e}n cascade, and show how these vary with changing mean magnetic field strength. Finally, we propose a method for using anisotropy in velocity centroid maps as a diagnostic of the mean magnetic field strength in observed cloud cores.

The evolution of the magnetorotational instability (MRI) during the transition from outburst to quiescence in a dwarf nova disk is investigated using three-dimensional MHD simulations. The shearing box approximation is adopted for the analysis, so that the efficiency of angular momentum transport is studied in a small local patch of the disk: this is usually referred as to a one-zone model. To take account of the low ionization fraction of the disk, the induction equation includes both ohmic dissipation and the Hall effect. We induce a transition from outburst to quiescence by an instantaneous decrease of the temperature. The evolution of the MRI during the transition is found to be very sensitive to the temperature of the quiescent disk. As long as the temperature is higher than a critical value of about 2000 K, MHD turbulence and angular momentum transport is sustained by the MRI. However, MHD turbulence dies away within an orbital time if the temperature falls below this critical value. In this case, the stress drops off by more than 2 orders of magnitude, and is dominated by the Reynolds stress associated with the remnant motions from the outburst. The critical temperature depends slightly on the distance from the central star and the local density of the disk.

Using numerical MHD simulations, we have studied the evolution of the magnetorotational instability in stratified accretion disks in which the ionization fraction (and therefore resistivity) varies substantially with height. This model is appropriate to dense, cold disks around protostars or dwarf nova systems which are ionized by external irradiation of cosmic rays or high-energy photons. We find the growth and saturation of the MRI occurs only in the upper layers of the disk where the magnetic Reynolds number exceeds a critical value; in the midplane the disk remains queiscent. The vertical Poynting flux into the "dead", central zone is small, however velocity fluctuations in the dead zone driven by the turbulence in the active layers generate a significant Reynolds stress in the midplane. When normalized by the thermal pressure, the Reynolds stress in the midplane never drops below about 10% of the value of the Maxwell stress in the active layers, even though the Maxwell stress in the dead zone may be orders of magnitude smaller than this. Significant mass mixing occurs between the dead zone and active layers. Fluctuations in the magnetic energy in the active layers can drive vertical oscillations of the disk in models where the ratio of the column density in the dead zone to that in the active layers is <10. These results have important implications for the global evolution of a layered disk, in particular there may be residual mass inflow in the dead layer. We discuss the effects that dust in the disk may have on our results.

We use 3D MHD simulations to investigate nonlinear development of the Parker, magneto-Jeans (MJI), and swing mechanisms in galactic disks. The model disks are local, isothermal, and begin from a vertically-stratified equilibrium. We first construct axisymmetric equilibria and examine their stability. Finite disk thickness reduces the critical Toomre Q parameter below unity; we find Q_c \~ 0.75, 0.72, and 0.57 for \beta=\infty, 10, and 1 cases, respectively. We then pursue fully 3D models. In non-self-gravitating cases, the peak density enhancement from the `pure' Parker instability is less a factor of two. The dominant growing modes have radial wavelengths comparable to the disk scale height H, much shorter than the azimuthal wavelength (~10-20 H). Shearing disks, being more favorable to midplane-symmetric modes, have somewhat different late-time magnetic field profiles from nonshearing disks, but otherwise saturated states are similar. Late-time velocity fluctuations at 10% of the sound speed persist, but no characteristic structural signatures of Parker modes remain in the new equilibria. In self-gravitating cases, the development of density structure is qualitatively similar to our previous results from thin-disk simulations. The Parker instability, although it may help seed structure or tip the balance under marginal conditions, appears to play a secondary role. In shearing disks with Q less than a threshold level ~ 1, swing amplification can produce bound clouds of a few times the local Jeans mass. The most powerful cloud-condensing mechanism appears to be the MJI. Our simulations show that condensations of a local Jeans mass (~3\times 10^7 M_sun) grow very rapidly, supporting the idea that MJI is at least partly responsible for the formation of bound cloud complexes in spiral galaxies.

The nonlinear evolution of the magnetorotational instability (MRI) in weakly ionized accretion disks, including the effect of the Hall term and ohmic dissipation, is investigated using local three-dimensional MHD simulations and various initial magnetic field geometries. When the magnetic Reynolds number, Re_M \equiv v_A^2 / \eta \Omega (where v_A is the Alfven speed, \eta the magnetic diffusivity, and \Omega the angular frequency), is initially larger than a critical value Re_{M, crit}, the MRI evolves into MHD turbulence in which angular momentum is transported efficiently by the Maxwell stress. If Re_M < Re_{M, crit}, however, ohmic dissipation suppresses the MRI, and the stress is reduced by several orders of magnitude. The critical value is in the range of 1 - 30 depending on the initial field configuration. The Hall effect does not modify the critical magnetic Reynolds number by much, but enhances the saturation level of the Maxwell stress by a factor of a few. We show that the saturation level of the MRI is characterized by v_{Az}^2 / \eta \Omega, where v_{Az} is the Alfven speed in the nonlinear regime along the vertical component of the field. The condition for turbulence and significant transport is given by v_{Az}^2 / \eta \Omega \gtrsim 1, and this critical value is independent of the strength and geometry of the magnetic field or the size of the Hall term. If the magnetic field strength in an accretion disk can be estimated observationally, and the magnetic Reynolds number v_A^2 / \eta \Omega is larger than about 30, this would imply the MRI is operating in the disk.

The effect of the Hall term on the evolution of the magnetorotational instability (MRI) in weakly ionized accretion disks is investigated using local axisymmetric simulations. First, we show that the Hall term has important effects on the MRI when the temperature and density in the disk is below a few thousand K and between 10^13 and 10^18 cm^{-3} respectively. Such conditions can occur in the quiescent phase of dwarf nova disks, or in the inner part (inside 10 - 100 AU) of protoplanetary disks. When the Hall term is important, the properties of the MRI are dependent on the direction of the magnetic field with respect to the angular velocity vector \Omega. If the disk is threaded by a uniform vertical field oriented in the same sense as \Omega, the axisymmetric evolution of the MRI is an exponentially growing two-channel flow without saturation. When the field is oppositely directed to \Omega, however, small scale fluctuations prevent the nonlinear growth of the channel flow and the MRI evolves into MHD turbulence. These results are anticipated from the characteristics of the linear dispersion relation. In axisymmetry on a field with zero-net flux, the evolution of the MRI is independent of the size of the Hall term relative to the inductive term. The evolution in this case is determined mostly by the effect of ohmic dissipation.

We perform numerical simulations of magneto-rotational instability in a local patch of accretion disk in which radiation pressure exceeds gas pressure. Such conditions may occur in the central regions of disks surrounding compact objects in active galactic nuclei and Galactic X-ray sources. We assume axisymmetry, and neglect vertical stratification. The growth rates of the instability on initially uniform magnetic fields are consistent with the linear analysis of Blaes & Socrates 2001. As is the case when radiation effects are neglected, the non-linear development of the instability leads to transitory turbulence when the initial magnetic field has no net vertical flux. During the turbulent phase, angular momentum is transported outwards. The Maxwell stress is a few times the Reynolds stress, and their sum is about four times the mean pressure in the vertical component of the magnetic field. For magnetic pressure exceeding gas pressure, turbulent fluctuations in the field produce density contrasts about equal to the ratio of magnetic to gas pressure. These are many times larger than in the corresponding gas pressure dominated situation, and may have profound implications for the steady-state vertical structure of radiation-dominated disks. Diffusion of radiation from compressed regions damps turbulent motions, converting kinetic energy into photon energy.

We report here results of high-resolution hydrodynamical simulations of gas flows in barred galaxies, with a focus on gas dynamics in the central kiloparsec. In a single bar with an Inner Lindblad Resonance, we find either near-circular motion of gas in the nuclear ring, or a spiral shock extending towards the galaxy center, depending on the sound speed in the gas. From a simple model of a dynamically-possible doubly barred galaxy with resonant coupling, we infer that the secondary bar is likely to end well inside its corotation. Such a bar cannot create shocks in the gas flow, and therefore will not reveal itself in color maps through straight dust lanes: the gas flows induced by it are different from those caused by the rapidly rotating main bars. In particular, we find that secondary stellar bars are unlikely to increase the mass inflow rate into the galactic nucleus.

The standard equilibrium for radiation-dominated accretion disks has long been known to be viscously, thermally, and convectively unstable, but the nonlinear development of these instabilities---hence the actual state of such disks---has not yet been identified. By performing local two-dimensional hydrodynamic simulations of disks, we demonstrate that convective motions can release heat sufficiently rapidly as to substantially alter the vertical structure of the disk. If the dissipation rate within a vertical column is proportional to its mass, the disk settles into a new configuration thinner by a factor of two than the standard radiation-supported equilibrium. If, on the other hand, the vertically-integrated dissipation rate is proportional to the vertically-integrated total pressure, the disk is subject to the well-known thermal instability. Convection, however, biases the development of this instability toward collapse. The end result of such a collapse is a gas pressure-dominated equilibrium at the original column density.

We present two-dimensional hydrodynamical simulations of both jet- and wind-driven models for protostellar outflows in order to make detailed comparisons to the kinematics of observed molecular outflows. Comparing the different simulations with observations, we find that some outflows, e.g., HH 212, show features consistent with the jet-driven model, while others, e.g., VLA 05487, are consistent with the wind-driven model.

We analyze the dynamics of the shell produced when a bow shock from a collimated jet propagates into the surrounding medium. Under interstellar conditions, the shock is radiative, and a ballistic approximation for the shell flow is appropriate beyond the working surface where the jet impacts its surroundings. The solution is then determined by the ambient and jet densities and velocities and by the momentum impulse applied in the working surface. Using estimates for these impulses (confirmed by separate numerical simulations), we obtain solutions for the shell structure, and for the range of velocities in the shell at any point. We provide predictions for the position-velocity and mass-velocity relations expected for plane-of-sky bow-shock shells, and for the bulk shell properties. In a companion paper, we show that these analytic solutions are in excellent agreement with the results of direct numerical simulations. We argue that classical molecular (CO) outflows cannot be purely jet-driven, because the bow-shock shell solutions are much too elongated compared with observations. Finally, we suggest that the ``spur'' structures seen in position-velocity diagrams of observed molecular outflows are the manifestation of internal bow shocks which may be fit with our simple dynamical models.

We present a global magnetohydrodynamic (MHD) three dimensional simulation of a nonradiative accretion flow originating in a pressure supported torus. The evolution is controlled by the magnetorotational instability which produces turbulence. The flow forms a nearly Keplerian disk. The total pressure scale height in this disk is comparable to the vertical size of the initial torus. Gas pressure dominates only near the equator; magnetic pressure is more important in the surrounding atmosphere. A magnetically dominated bound outflow is driven from the disk. The accretion rate through the disk exceeds the final rate into the hole, and a hot torus forms inside 10 r_g. Hot gas, pushed up against the centrifugal barrier and confined by magnetic pressure, is ejected in a narrow, unbound, conical outflow. The dynamics are controlled by magnetic turbulence, not thermal convection, and a hydrodynamic alpha model is inadequate to describe the flow. The limitations of two dimensional MHD simulations are also discussed.

A module for the ZEUS-2D code is described which may be used to solve the equations of radiation hydrodynamics to order unity in v/c, in the flux-limited diffusion (FLD) approximation. In this approximation, the tensor Eddington factor f which closes the radiation moment equations is chosen to be an empirical function of radiation energy density. This is easier to implement and faster than full-transport techniques, in which f is computed by solving the transfer equation. However, FLD is less accurate when the flux has a component perpendicular to the gradient in radiation energy density, and in optically thin regions when the radiation field depends strongly on angle. The material component of the fluid is here assumed to be in local thermodynamic equilibrium. The energy equations are operator-split, with transport terms, radiation diffusion term, and other source terms evolved separately. Transport terms are applied using the same consistent transport algorithm as in ZEUS-2D. The radiation diffusion term is updated using an alternating-direction implicit method with convergence checking. Remaining source terms are advanced together implicitly using numerical root-finding. However when absorption opacity is zero, accuracy is improved by treating compression and expansion source terms using time-centered differencing. Results are discussed for test problems including radiation-damped linear waves, radiation fronts propagating in optically-thin media, subcritical and supercritical radiating shocks, and an optically-thick shock in which radiation dominates downstream pressure.

We present the results of axisymmetric, time-dependent magnetohydrodynamic simulations of accretion flows around black holes. The calculations begin from a rotationally supported thick torus which contains a weak poloidal field. Accretion is produced by growth and saturation of the magnetorotational instability (MRI) provided the wavelength of the fastest growing mode is less than the thickness of the torus. Using a computational grid that spans more than two decades in radius, we compare the time-averaged properties of the flow to previous hydrodynamical simulations. The net mass accretion rate is small compared to the mass inflow and outflow rates at large radii associated with turbulent eddies. Turbulence is driven by the MRI rather than convection. The two-dimensional structure of the time-averaged flow is significantly different compared to the hydrodynamical case. We discuss the limitations imposed on our results by the assumption of axisymmetry and the relatively small radial domain.

We use 3D numerical MHD simulations to follow the evolution of cold, turbulent, gaseous systems with parameters representing GMC conditions. We study three cloud simulations with varying mean magnetic fields, but identical initial velocity fields. We show that turbulent energy is reduced by a factor two after 0.4-0.8 flow crossing times (2-4 Myr), and that the magnetically supercritical cloud models collapse after ~6 Myr, while the subcritical cloud does not collapse. We compare density, velocity, and magnetic field structure in three sets of snapshots with matched Mach numbers. The volume and column densities are both log-normally distributed, with mean volume density a factor 3-6 times the unperturbed value, but mean column density only a factor 1.1-1.4 times the unperturbed value. We use a binning algorithm to investigate the dependence of kinetic quantities on spatial scale for regions of column density contrast (ROCs). The average velocity dispersion for the ROCs is only weakly correlated with scale, similar to the mean size-linewidth relation for clumps within GMCs. ROCs are often superpositions of spatially unconnected regions that cannot easily be separated using velocity information; the same difficulty may affect observed GMC clumps. We analyze magnetic field structure, and show that in the high density regime, total magnetic field strengths increase with density with logarithmic slope 1/3 -2/3. Mean line-of-sight magnetic field strengths vary widely across a projected cloud, and do not correlate with column density. We compute simulated interstellar polarization maps at varying orientations, and determine that the Chandrasekhar-Fermi formula multiplied by a factor ~0.5 yields a good estimate of the plane-of sky magnetic field strength provided the dispersion in polarization angles is < 25 degrees.

The effect of optically thin radiative cooling on the Kelvin-Helmholtz instability of three dimensional jets is investigated via linear stability theory and nonlinear hydrodynamical simulation. Two different cooling functions are considered: radiative cooling is found to have a significant effect on the stability of the jet in each case. The wavelengths and growth rates of unstable modes in the numerical simulations are found to be in good agreement with theoretical predictions. Disruption of the jet is found to be sensitive to the precessional frequency at the origin with lower frequencies leading to more rapid disruption. Strong nonlinear effects are observed as the result of the large number of normal modes in three dimensions which provide rich mode-mode interactions. These mode-mode interactions provide new mechanisms for the formation of knots in the flows. Significant structural features found in the numerical simulations appear similar to structures observed on protostellar jets.

We present the results of time-dependent hydrodynamic calculations of line- driven winds from accretion disks in AGN. To calculate the radiation force, we take into account radiation from the disk and the central engine. The gas temperature and ionization state in the wind are calculated self-consistently. We find that a disk accreting onto a 10^8 MSUN black hole at the rate of 1.8 MSUN/YR can launch a wind at about 10^16 cm from the central engine. The X-rays from the central object are attenuated by the disk atmosphere so they cannot prevent the local disk radiation from pushing matter away from the disk. For a reasonable X-ray opacity, the disk wind can be accelerated by the central UV radiation to velocities of up to 15000 km/s at a distance of about 10^17 cm from the central engine. The covering factor of the disk wind is $\sim 0.2$. The wind is unsteady and consists of an opaque, slow vertical flow near the disk that is bounded on the polar side by a high-velocity stream. A typical column density through the fast stream is a few 10^23 cm^-2. The fast stream contributes nearly 100% to the total wind mass loss rate of 0.5 MSUN/YR.

We present the results of a series of axisymmetric time-dependent magnetohydrodynamic (MHD) simulations of the propagation of cooling, overdense jets, motivated by the properties of outflows associated with young stellar objects. A variety of initial field strengths and configurations are explored for both steady and time-variable (pulsed) jets. Even apparently weak magnetic fields with strengths B < 60 micro-G in the pre-shocked jet beam can have a significant effect on the dynamics, for example by altering the density, width, and fragmentation of thin shells formed by cooling gas. A linear analysis predicts that axisymmetric pinch modes of the MHD Kelvin-Helmholtz instability should grow only slowly for the highly supermagnetosonic jets studied here; we find no evidence for them in our simulations. Some of our models appear unstable to current-driven pinch modes, however the resulting pressure and density variations induced in the jet beam are not large, making this mechanism an unlikely source of emission knots in the jet beam. In the case of pulsed jets, radial hoop stresses confine shocked jet material in the pulses to the axis, resulting in a higher density in the pulses in comparison to purely hydrodynamic models.

We present three-dimensional magnetohydrodynamic simulations of the nonlinear evolution of the magnetorotational instability (MRI) with a non-zero Ohmic resistivity. The properties of the saturated state depend on the initial magnetic field configuration. In simulations with an initial uniform vertical field, the MRI is able to support angular momentum transport even for large resistivities through the quasi-periodic generation of axisymmetric radial channel solutions rather than through the maintenance of anisotropic turbulence. Simulations with zero net flux show that the angular momentum transport and the amplitude of magnetic energy after saturation are significantly reduced by finite resistivity, even at levels where the linear modes are only slightly affected. This occurs at magnetic Reynolds numbers expected in low, cool states of dwarf novae, these results suggest that finite resistivity may account for the low and high angular momentum transport rates inferred for these systems.

We use three-dimensional magnetohydrodynamical (MHD) simulations to study the formation of a corona above an initially weakly magnetized, isothermal accretion disk. We also describe a modification to time-explicit numerical algorithms for MHD which enables us to evolve highly stratified disks for many orbital times. We find that MHD turbulence driven by the magnetorotational instability (MRI) produces strong amplification of weak fields within two scale heights of the disk midplane in a few orbital times. About 25 % of the magnetic energy generated by the MRI within two scale heights escapes due to buoyancy, producing a strongly magnetized corona above the disk. Most of the buoyantly rising magnetic energy is dissipated between 3 and 5 scale heights, suggesting the corona will also be hot. The average vertical disk structure consists of a weakly magnetized turbulent core below a strongly magnetized corona which is stable to the MRI. The largescale field structure in both the disk and corona is toroidal. The functional form of the stress is flat within two scale heights, but proportional to the density above two scale heights. For initially weak uniform vertical fields, we find the exponential growth of magnetic field via axisymmetric vertical modes of the MRI produces strongly buoyant sheets of magnetic energy which break the disk apart into horizontal channels. These channels rise several scale heights vertically before the onset of the Parker instability distorts the sheets and allows matter to flow back towards the midplane and reform a disk. We suggest this evolution may be relevant to the dynamical processes which disrupt the inner regions of a disk when it interacts with a strongly magnetized central object.

Two-dimensional (axially symmetric) numerical hydrodynamical calculations of accretion flows which cannot cool through emission of radiation are presented. The calculations begin from an equilibrium configuration consisting of a thick torus with constant specific angular momentum. Accretion is induced by the addition of a small anomalous azimuthal shear stress which is characterized by a function \nu. We study the flows generated as the amplitude and form of \nu are varied. A spherical polar grid which spans more than two orders of magnitude in radius is used to resolve the flow over a wide range of spatial scales. We find that convection in the inner regions produces significant outward mass motions that carry away both the energy liberated by, and a large fraction of the mass participating in, the accretion flow. Although the instantaneous structure of the flow is complex and dominated by convective eddies, long time averages of the dynamical variables show remarkable correspondence to certain steady-state solutions. Near the equatorial plane, the radial profiles of the time-averaged variables are power-laws with an index that depends on the radial scaling of the shear stress. We find that regardless of the adiabatic index of the gas, or the form or magnitude of the shear stress, the mass inflow rate is a strongly increasing function of radius, and is everywhere nearly exactly balanced by mass outflow. The net mass accretion rate through the disc is only a fraction of the rate at which mass is supplied to the inflow at large radii, and is given by the local, viscous accretion rate associated with the flow properties near the central object.

We describe an efficient method of calculating the radiation pressure due to spectral lines, including all the terms in the velocity gradient tensor. We apply this method to calculate the two-dimensional, time-dependent structure of winds from luminous disks. Qualitative features of our new models are very similar to those we calculated including only the dominant terms in the tensor (Proga, Stone & Drew 1998, hereafter PSD). In particular, we find that models which displayed unsteady behaviour in PSD are also unsteady with the new method, and gross properties of the winds, such as mass-loss rate and characteristic velocity are not changed by the more accurate approach. The largest change caused by the new method is in the disk-wind opening angle: winds driven only by the disk radiation are more polar with the new method while winds driven by the disk and central object radiation are typically more equatorial. In the closing discussion, we provide further insight into the way the geometry of the radiation field and consequent flow determines the time properties of the flow.

We review recent results from numerical simulations and related models of MHD turbulence in the interstellar medium (ISM) and in molecular clouds. We discuss the implications of turbulence for the processes of cloud formation and evolution, and the determination of clouds' physical properties. Numerical simulations of the turbulent ISM to date have included magnetic fields, self-gravity, parameterized heating and cooling, modeled star formation and other turbulent inputs. The structures which form reproduce well observed velocity-size scaling properties, while predicting the non-existence of a general density-size scaling law. Criteria for the formation of gravitationally-bound structures by turbulent compression are summarized. For flows with equations of state $P\propto \rho^\gamma$, the statistics of the density field depend on the exponent $\gamma$. Numerical simulations of both forced and decaying MHD compressible turbulence have shown that the decay rate is comparable to the non-magnetic case. For virialized clouds, the turbulent decay time is shorter than the gravitational free-fall time, so wholesale cloud collapse is only prevented by ongoing turbulent inputs and/or a strong mean magnetic field. Finally, perspectives for future work in this field are briefly discussed.

We present a numerical model for the extended steep-spectrum radio sources and the elongated X-ray structure in A3667 based on new 3-dimensional MHD/N-body simulations. The X-ray and optical analyses of A3667 indicate that it has undergone a recent subcluster merger event. We believe that the Mpc-scale radio sources identified in A3667 are also a consequence of the merger. Our previous numerical simulations show that mergers often produce large-scale shocks and turbulence capable of both magnetic field amplification and in-situ reacceleration of relativistic particles. Our model suggests that these radio structures, separated by ~2.6 h^{-1}_100 Mpc, are in fact causally linked via a slightly off-axis merger that occurred nearly in the plane of the sky approximately 1 Gyr ago with a subcluster having a total mass equal to ~20% of the primary cluster.

We present initial results from the first 3-dimensional numerical magnetohydrodynamical (MHD) simulations of magnetic field evolution in merging clusters of galaxies. Within the framework of idealized initial conditions similar to our previous work, we look at the gasdynamics and the magnetic field evolution during a major merger event in order to examine the suggestion that shocks and turbulence generated during a cluster/subcluster merger can produce magnetic field amplification and relativistic particle acceleration and, as such, may play a role in the formation and evolution of cluster-wide radio halos. The ICM, as represented by the equations of ideal MHD, is evolved self-consistently within a changing gravitational potential defined largely by the collisionless dark matter component represented by an N-body particle distribution. The MHD equations are solved by the Eulerian, finite-difference code, ZEUS. The particles are evolved by a standard particle-mesh (PM) code. We find significant evolution of the magnetic field structure and strength during two distinct epochs of the merger evolution.

The molecular component of the Galaxy is comprised of turbulent, magnetized clouds, many of which are self-gravitating and form stars. To understand how these clouds' evolution may depend on their level of turbulence, mean magnetization, and degree of self-gravity, we perform a survey of direct numerical MHD simulations in 2.5 dimensions. Two of our independent survey parameters allow us to model clouds which either meet or fail conditions for magneto-Jeans stability and magnetic criticality. Our third survey parameter allows us to initiate turbulence of either sub- or super-Alfv\'enic amplitude. We evaluate the times for each cloud model to become gravitationally bound, and measure each model's kinetic energy loss over the fluid flow crossing time. We compare the evolution of density and magnetic field structural morphology, and quantify the differences in the density contrast generated by internal stresses, for models of differing mean magnetization. We find that the values of (c_s/v_A)^2 and L/L_Jeans, but not the initial Mach number v/c_s, determine the time for cloud gravitational binding and collapse. We find, contrary to some previous expectations, less than a factor of two difference between turbulent decay times for models with varying magnetic field strength. In all models, we find turbulent amplification in the magnetic field strength, with the turbulent magnetic energy between 25-60% of the turbulent kinetic energy after one flow crossing time. We find that for non-self-gravitating stages of evolution, the mass-averaged density contrast magnitudes <log(\rho/\bar\rho)> are in the range 0.2-0.5; we note that only the more strongly-magnetized models appear consistent with estimates of clump/interclump density contrasts inferred in Galactic GMCs.

We report results of a three dimensional, high resolution (up to 512^3) numerical investigation of supersonic compressible magnetohydrodynamic turbulence. We consider both forced and decaying turbulence. The model parameters are appropriate to conditions found in Galactic molecular clouds. We find that the dissipation time of turbulence is of order the flow crossing time or smaller, even in the presence of strong magnetic fields. About half the dissipation occurs in shocks. Weak magnetic fields are amplified and tangled by the turbulence, while strong fields remain well ordered.

We carry out three-dimensional magnetohydrodynamical simulations of the magnetorotational (Balbus-Hawley) instability in weakly-ionized plasmas. We adopt a formulation in which the ions and neutrals each are treated as separate fluids coupled only through a collisional drag term. Ionization and recombination processes are not considered. The linear stability of the ion-neutral system has been previously considered by Blaes and Balbus (1994). Here we extend their results into the nonlinear regime by computing the evolution of Keplerian angular momentum distribution in the local shearing box approximation. We find significant turbulence and angular momentum transport when the collisional frequency is on order 100 times the orbital frequency. At higher collision rates, the two-fluid system studied here behaves much like the fully ionized systems studied previously. At lower collision rates the evolution of the instability is determined primarily by the properties of the ions, with the neutrals acting as a sink for the turbulence. Since in this regime saturation occurs when the magnetic field is superthermal with respect to the ion pressure, we find the amplitude of the magnetic energy and the corresponding angular momentum transport rate is proportional to the ion density. Our calculations show the ions and neutrals are essentially decoupled when the collision frequency is less than 0.01 times the orbital frequency; in this case the ion fluid behaves as in the single fluid simulations and the neutrals remain quiescent. We find that purely toroidal initial magnetic field configurations are unstable to the magnetorotational instability across the range of coupling frequencies.

Our understanding of the dynamical processes which control the structure and evolution of the interaction region between an accretion disk and the central star is reviewed. If the central star is unmagnetized, this interaction is in the form of a classical boundary layer. However, if the central star is strongly magnetized, it is expected that the inner disk is truncated by the stellar field, and the accretion flow follows the stellar field lines to the magnetic poles. Several outstanding questions remain regarding this picture. It is likely that numerical MHD simulations will prove essential to answering these questions. However, in order to model the MHD turbulence, angular momentum transport, and possible dynamo action associated with the dynamics of star-disk interaction, fully three-dimensional simulations spanning long dynamical times are required.

A radiation-driven disk wind model is proposed that offers great promise of explaining the extreme mass loss signatures of massive young stellar objects (the BN-type objects and more luminous Herbig Be stars). It is argued that the dense low-velocity winds associated with young late-O/early-B stars would be the consequence of continuing optically-thick accretion onto them. The launch of outflow from a Keplerian disk allows wind speeds of about 200 km/s that are substantially less than the escape speed from the stellar surface. The star itself is not required to be a rapid rotator. Disk irradiation is taken into account in the hydrodynamical calculation presented, and identified as an important issue both observationally and from the dynamical point of view.

We study the 2-D, time-dependent hydrodynamics of radiation-driven winds from accretion disks in which the radiation force is mediated by spectral lines. If the dominant contribution to the total radiation field comes from the disk, then we find the outflow is intrinsically unsteady and characterised by large amplitude velocity and density changes. Both infall and outflow can occur in different regions of the wind at the same time. On the other hand, if the total luminosity of the system is dominated by the central star, then the outflow is steady. In either case, we find the 2-D structure of the wind consists of a dense, slow outflow, typically confined to angles within about 45$^o$ of the equatorial plane, that is bounded on the polar side by a high -velocity, lower density stream. The flow geometry is controlled largely by the geometry of the radiation field. Global properties e.g., the total mass loss rate and terminal velocity depend more on the system luminosity and are insensitive to geometry. Matter is fed into the fast wind from within a few stellar radii of the central star. Our solutions agree qualitatively with the kinematics of outflows in CV systems inferred from spectroscopic observations. We predict that low luminosity systems may display unsteady behavior in wind-formed spectral lines. Our study also has application to winds from active galactic nuclei and from high mass YSOs.

A754 is a well-observed cluster of galaxies which exhibits a variety of morphological peculiarities. These include a bar of X-ray emission that is offset significantly from the galaxy distribution, an elongated X-ray surface brightness distribution extending between two distinct peaks in the galaxy distribution, and an extremely non-isothermal and asymmetric intracluster medium (ICM) temperature morphology. Using these observational constraints, we present a numerical Hydro/N-body model of A754 in which two clusters (2.5:1 mass ratio) have merged nearly in the plane of the sky less than 0.5 Gyrs ago with an impact parameter of ~120 kpc, and an impact velocity of ~2500 km/s (roughly the escape velocity of the primary cluster). Our models allow us to identify the origin of A754's peculiar X-ray and temperature morphologies, the underlying hydrodynamical processes that shape them, and their future evolution. We make detailed predictions for future high resolution X-ray spectroscopic observations (e.g. ASTRO-E). We discuss general properties of our models which will be characteristic of off-axis mergers. In particular, we find significant non-thermal pressure support within the central region which could bias cluster mass estimates. We find significant angular momentum imparted on the gas distribution in the cluster. We find that mixing of the subcluster gas components is an inefficient process, particularly at large radii. Finally, we find that subsequent dark matter core passages result in an extended relaxation timescale.

The nonlinear evolution of unstable C-type shocks in weakly ionized plasmas is studied by means of time-dependent magnetohydrodynamical simulations. This study is limited to shocks in magnetically dominated plasmas (in which the Alfven speed in the neutrals greatly exceeds the sound speed), and microphysical processes such as ionization and recombination are not followed. Both two-dimensional simulations of planar perpendicular and oblique C-type shocks, and fully three-dimensional simulation of a perpendicular shock are presented.

We have modeled the gas temperature structure in unstable C-type shocks and obtained predictions for the resultant CO and H2 rotational line emissions, using numerical simulations of the Wardle instability that were carried out by Stone (1997) and that have been described in a companion paper. We obtained results for the gas temperature distribution in -- and H2 and CO line emission from -- shocks of neutral Alfvenic Mach number 10 and velocity 20 or 40 km/sec in which the Wardle instability has saturated. Although the Wardle instability profoundly affects the density structure behind C-type shocks, most of the shock-excited molecular line emission is generated upstream of the region where the strongest effects of the instability are felt. Thus the Wardle instability has a relatively small effect upon the overall gas temperature distribution in -- and the emission line spectrum from -- C-type shocks, at least for the cases that we have considered.

Values of the Hubble constant reported to date which are based on measurement of the Sunyaev-Zeldovich (SZ) effect in clusters of galaxies are systematically lower than those derived by other methods (e.g., Cepheid variable stars, or the Tully-Fisher relation). We investigate the possibility that systematic errors may be introduced into the analysis by the generally adopted assumptions that observed clusters are in hydrostatic equilibrium, are spherically symmetric, and are isothermal. We construct self-consistent theoretical models of merging clusters of galaxies using hydrodynamical/N-body simulations. We then compute the magnitude of Ho derived from the SZ effect at different times and at different projection angles both from first principles, and by applying each of the standard assumptions used in the interpretation of observations. Our results indicate that the assumption of isothermality in the evolving clusters can result in Ho being underestimated by 10-30% depending on both epoch and projection angle. Moreover, use of the projected, emission-weighted temperature profile under the assumption of spherical symmetry does not significantly improve the situation except in the case of more extreme mergers (i.e., those involving relatively gas-rich subclusters).