▪ Abstract We summarize the properties of FU Orionis variables, and show how accretion disk models simply explain many peculiarities of these objects. FU Ori systems demonstrate that disk accretion in early stellar evolution is highly episodic, varying from ∼ 10−7 yr−1 in the low (T Tauri) state to 10−4 yr−1 in the high (FU Ori) state. This variability in mass accretion is matched by a corresponding variability in mass ejection, with mass loss rates reaching ∼ 10−1 of the mass accretion rates in outburst. It appears that the FU Ori phenomenon is restricted to early phases of stellar evolution, probably with infall still occuring to the disk, which may help drive repetitive outbursts. Thermal instabilities are a promising way to produce FU Ori disk outbursts, although many uncertainties remain in the theory; triggering by interactions with companion stars on eccentric orbits may also play a role.

▪ Abstract Turbulence affects the structure and motions of nearly all temperature and density regimes in the interstellar gas. This two-part review summarizes the observations, theory, and simulations of interstellar turbulence and their implications for many fields of astrophysics. The first part begins with diagnostics for turbulence that have been applied to the cool interstellar medium and highlights their main results. The energy sources for interstellar turbulence are then summarized along with numerical estimates for their power input. Supernovae and superbubbles dominate the total power, but many other sources spanning a large range of scales, from swing-amplified gravitational instabilities to cosmic ray streaming, all contribute in some way. Turbulence theory is considered in detail, including the basic fluid equations, solenoidal and compressible modes, global inviscid quadratic invariants, scaling arguments for the power spectrum, phenomenological models for the scaling of higher-order structu...

▪ Abstract The inner few parsecs at the Galactic Center have come under intense scrutiny in recent years, in part due to the exciting broad-band observations of this region, but also because of the...

▪ Abstract More than four decades have elapsed since modeling of the core collapse supernova mechanism began in earnest. To date, the mechanism remains elusive, at least in detail, although significant progress has been made in understanding these multiscale, multiphysics events. One-, two-, and three-dimensional simulations of or relevant to core collapse supernovae have shown that (a) neutrino transport, (b) fluid instabilities, (c) rotation, and (d) magnetic fields, together with proper treatments of (e) the sub- and super- nuclear density stellar core equation of state, (f) the neutrino interactions, and (g) gravity are all important. The importance of these ingredients applies to both the explosion mechanism and to phenomena directly associated with the mechanism, such as neutron star kicks, supernova neutrino and gravitational wave emission, and supernova spectropolarimetry. Not surprisingly, current two- and three-dimensional models have yet to include (a)–(d) with sufficient realism. One-dimension...

▪ Abstract Extragalactic jets were discovered and initially studied by radio astronomers in connection with extended radio sources. At present, the combination of jets and disks is considered the crucial element in unification models for all active galactic nuclei (AGNs). The acceleration and propagation conditions of jets, together with the aspect ratio of the disk/jet geometry with respect to the observer, shape the morphologies of AGNs. However, these phenomenological models are very complex from the physical and mathematical point of view, as they involve different elements of the theories of gravitation, fluid dynamics, and electrodynamics in a highly nonlinear combination and in conditions not easily reproducible in laboratory plasma or fluid experiments. In the last ten years, theorists have attacked the subject with advanced analytical and numerical methods, and some important results have already been established that confirm the global scenario, although we are still far from a complete physical...

When surrounded by a transparent emission region, black holes are expected to reveal a dark shadow caused by gravitational light bending and photon capture at the event horizon. To image and study this phenomenon, we have assembled the Event Horizon Telescope, a global very long baseline interferometry array observing at a wavelength of 1.3 mm. This allows us to reconstruct event-horizon-scale images of the supermassive black hole candidate in the center of the giant elliptical galaxy M87. We have resolved the central compact radio source as an asymmetric bright emission ring with a diameter of 42 ± 3 μas, which is circular and encompasses a central depression in brightness with a flux ratio ≳10:1. The emission ring is recovered using different calibration and imaging schemes, with its diameter and width remaining stable over four different observations carried out in different days. Overall, the observed image is consistent with expectations for the shadow of a Kerr black hole as predicted by general relativity. The asymmetry in brightness in the ring can be explained in terms of relativistic beaming of the emission from a plasma rotating close to the speed of light around a black hole. We compare our images to an extensive library of ray-traced general-relativistic magnetohydrodynamic simulations of black holes and derive a central mass of M = (6.5 ± 0.7) × 109 M⊙. Our radio-wave observations thus provide powerful evidence for the presence of supermassive black holes in centers of galaxies and as the central engines of active galactic nuclei. They also present a new tool to explore gravity in its most extreme limit and on a mass scale that was so far not accessible.

We utilize the multi-wavelength data of M33 to study the origin of turbulence in its interstellar medium. We find that the H i turbulent energy surface density inside 8 kpc is ∼1–3 × 1046 erg pc−2, and has no strong dependence on galactocentric radius because of the lack of variation in H i surface density and H i velocity dispersion. Then, we consider the energies injected by supernovae (SNe), the magneto-rotational instability (MRI), and the gravity-driven turbulence from accreted materials as the sources of turbulent energy. For a constant dissipation time of turbulence, the SNe energy can maintain turbulence inside ∼4 kpc radius (equivalent to ∼0.5 R25), while the MRI energy is always smaller than the turbulent energy within 8 kpc radius. However, when we let the dissipation time to be equal to the crossing time of turbulence across the H i scale height, the SNe energy is enough to maintain turbulence out to 7 kpc radius, and the sum of SNe and MRI energies is able to maintain turbulence out to 8 kpc radius. Due to lack of constraint in the mass accretion rate through the disk of M33, we cannot rule out the accretion driven turbulence as a possible source of energy. Furthermore, by resolving individual giant molecular clouds in M33, we also show that the SNe energy can maintain turbulence within individual molecular clouds with ∼1% of coupling efficiency. This result strengthens the proposition that stellar feedback is an important source of energy to maintain turbulence in nearby galaxies.

We use two-dimensional hydrodynamic simulations of viscous disks to examine whether dynamically interacting multiple giant planets can explain the large gaps (spanning over one order of magnitude in radius) inferred for the transitional and pre-transitional disks around T Tauri stars. In the absence of inner disk dust depletion, we find that it requires three to four giant planets to open up large enough gaps to be consistent with inferences from spectral energy distributions, because the gap width is limited by the tendency of the planets to be driven together into 2:1 resonances. With very strong tidal torques and/or rapid planetary accretion, fewer planets can also generate a large cavity interior to the locally formed gap(s) by preventing outer disk material from moving in. In these cases, however, the reduction of surface density produces a corresponding reduction in the inner disk accretion rate onto the star; this makes it difficult to explain the observed accretion rates of the pre-transitional/transitional disks. We find that even with four planets in disks, additional substantial dust depletion is required to explain observed disk gaps/holes. Substantial dust settling and growth, with consequent significant reductions in optical depths, is inferred for typical T Tauri disks in any case, and an earlier history of dust growth is consistent with the hypothesis that pre-transitional/transitional disks are explained by the presence of giant planets. We conclude that the depths and widths of gaps and disk accretion rates in pre-transitional/transitional disks cannot be reproduced by a planet-induced gap opening scenario alone. Significant dust depletion is also required within the gaps/holes. Order-of-magnitude estimates suggest that the mass of small dust particles (1 μm) relative to the gas must be depleted to 10 −5 to 10 −2 of the interstellar medium value, implying a very efficient mechanism of small dust removal or dust growth.

We study time evolution of sub-Keplerian transonic accretion flows onto black holes using a general relativistic numerical simulation code. We perform simulations in Schwarzschild spacetime. We first compare one-dimensional simulation results with theoretical results and validate the performance of our code. Next, we present results of axisymmetric, two-dimensional simulation of advective flows. We find that even in this case, for which no complete theoretical analysis is present in the literature, steady state shock formation is possible.

We study the x-ray intensity variations obtained from the time series of 155 light curves of x-ray binary systems collected by the instrument All Sky Monitor on board the satellite Rossi X-Ray Timing Explorer. These intensity distributions are adequately fitted by q-Gaussian distributions which maximize the Tsallis entropy and in turn satisfy a nonlinear Fokker–Planck equation, indicating their nonextensive and nonequilibrium behavior. From the values of the entropic index q obtained, we give a physical interpretation of the dynamics in x-ray binary systems based on the kinetic foundation of generalized thermostatistics. The present findings indicate that the binary systems display a nonextensive and turbulent behavior.

We study the properties of the turbulence driven by the magnetorotational instability in a stratified shearing box with outflow boundary conditions and an equation of state determined by self-consistent dissipation and radiation losses. A series of simulations with increasing resolution are performed within a fixed computational box. We achieve numerical convergence with respect to radial and azimuthal resolution. As vertical resolution is improved, the ratio of stress to pressure increases slowly, but the absolute levels of both the stress and the pressure increase noticeably. These results are in contrast with those of previous work on unstratified shearing boxes, in which improved resolution caused a diminution in the magnetic field strength. We argue that the persistence of strong magnetic field at higher resolution found in the stratified case is due to buoyancy. In addition, we find that the time-averaged vertical correlation length of the magnetic field near the disk midplane is 3 times larger than that found in previous unstratified simulations, decreasing slowly with improved vertical resolution. We further show that the undulatory Parker instability drives the magnetic field upwelling at several scale heights from the midplane that is characteristic of stratified magnetohydrodynamics-turbulent disks.

We study the possible origin of hydrodynamic turbulence in cold accretion disks such as those in star-forming systems and quiescent cataclysmic variables. As these systems are expected to have neutral gas, the turbulent viscosity is likely to be hydrodynamic in origin, not magnetohydrodynamic. Therefore, MRI will be sluggish or even absent in such disks. Although there are no exponentially growing eigenmodes in a hydrodynamic disk because of the non-normal nature of the eigenmodes, a large transient growth in the energy is still possible, which may enable the system to switch to a turbulent state. For a Keplerian disk, we estimate that the energy will grow by a factor of 1000 for a Reynolds number close to a million.

We study the effects of magnetic fields on the evolution of differentially rotating neutron stars, which can be formed in stellar core collapse or binary neutron star coalescence. Magnetic braking and the magnetorotational instability (MRI) both act on differentially rotating stars to redistribute angular momentum. Simulations of these stars are carried out in axisymmetry using our recently developed codes which integrate the coupled Einstein-Maxwell-MHD equations. We consider stars with two different equations of state (EOS), a gamma-law EOS with $\ensuremath{\Gamma}=2$, and a more realistic hybrid EOS, and we evolve them adiabatically. Our simulations show that the fate of the star depends on its mass and spin. For initial data, we consider three categories of differentially rotating, equilibrium configurations, which we label normal, hypermassive and ultraspinning. Normal configurations have rest masses below the maximum achievable with uniform rotation, and angular momentum below the maximum for uniform rotation at the same rest mass. Hypermassive stars have rest masses exceeding the mass limit for uniform rotation. Ultraspinning stars are not hypermassive, but have angular momentum exceeding the maximum for uniform rotation at the same rest mass. We show that a normal star will evolve to a uniformly rotating equilibrium configuration. An ultraspinning star evolves to an equilibrium state consisting of a nearly uniformly rotating central core, surrounded by a differentially rotating torus with constant angular velocity along magnetic field lines, so that differential rotation ceases to wind the magnetic field. In addition, the final state is stable against the MRI, although it has differential rotation. For a hypermassive neutron star, the MHD-driven angular momentum transport leads to catastrophic collapse of the core. The resulting rotating black hole is surrounded by a hot, massive, magnetized torus undergoing quasistationary accretion, and a magnetic field collimated along the spin axis\char22{}a promising candidate for the central engine of a short gamma-ray burst.

We study the effects of high-energy particles (HEPs) on the accretion flows onto a supermassive black hole and luminosities of escaping particles such as protons, neutrons, gamma rays, and neutrinos. We formulate a one-dimensional model of the two-component accretion flow consisting of thermal particles and HEPs, supposing that some fraction of the released energy is converted to the acceleration of HEPs. The thermal component is governed by fluid dynamics while the HEPs obey the moment equations of the diffusion–convection equation. By solving the time evolution of these equations, we obtain advection-dominated flows as the steady state solutions. The effects of the HEPs on the flow structures turn out to be small even if the pressure of the HEPs dominates over the thermal pressure. For a model in which the escaping protons take away almost all the energy released, the HEPs have a large enough influence to make the flow have a Keplerian angular velocity at the inner region. We calculate the luminosities of the escaping particles for these steady solutions. The escaping particles can extract the energy from about to , where is the mass accretion rate. The luminosities of the escaping particles depend on parameters such as the injection Lorentz factors, the mass accretion rates, and the diffusion coefficients. We also discuss some implications on the relativistic jet production by the escaping particles.

We study the dust motion at the surface layer of protoplanetary disks. Dust grains in surface layer migrate outward due to angular momentum transport via gas-drag force induced by the stellar radiation pressure. In this study, we calculate mass flux of the outward motion of compact grains and porous dust aggregates by the radiation pressure. The radiation pressure force for porous dust aggregates is calculated using the T-Matrix Method for the Clusters of Spheres. First, we confirm that porous dust aggregates are forced by strong radiation pressure even if they grow to be larger aggregates in contrast to homogeneous and spherical compact grains to which efficiency of radiation pressure becomes lower when their sizes increase. In addition, we find that the outward mass flux of porous dust aggregates with monomer size of 0.1 $\mu$m is larger than that of compact grains by an order of magnitude at the disk radius of 1 AU, when their sizes are several microns. This implies that large compact grains like calcium-aluminum rich inclusions (CAIs) are hardly transported to outer region by stellar radiation pressure, whereas porous dust aggregates like chondritic-porous interplanetary dust particles (CP-IDPs) are efficiently transported to comet formation region. Crystalline silicates are possibly transported in porous dust aggregates by stellar radiation pressure from inner hot region to outer cold cometary region in the protosolar nebula.

We study the behavior of magnetorotational turbulence in shearing box simulations with a radial and azimuthal extent up to 10 scale heights. Maxwell and Reynolds stresses are found to increase by more than a factor of 2 when increasing the box size beyond two scale heights in the radial direction. Further increase of the box size has little or no effect on the statistical properties of the turbulence. An inverse cascade excites magnetic field structures at the largest scales of the box. The corresponding 10% variation in the Maxwell stress launches a zonal flow of alternating sub- and super-Keplerian velocity. This, in turn, generates a banded density structure in geostrophic balance between pressure and Coriolis forces. We present a simplified model for the appearance of zonal flows, in which stochastic forcing by the magnetic tension on short timescales creates zonal flow structures with lifetimes of several tens of orbits. We experiment with various improved shearing box algorithms to reduce the numerical diffusivity introduced by the supersonic shear flow. While a standard finite difference advection scheme shows signs of a suppression of turbulent activity near the edges of the box, this problem is eliminated by a new method where the Keplerian shear advection is advanced in time by interpolation in Fourier space.

We study model images and polarization maps of Sagittarius A* at 230 GHz. We post-process GRMHD simulations and perform a fully relativistic radiative transfer calculation of the emitted synchrotron radiation to obtain polarized images for a range of mass accretion rates and electron temperatures. At low accretion rates, the polarization map traces the underlying toroidal magnetic field geometry. At high accretion rates, we find that Faraday rotation internal to the emission region can depolarize and scramble the map. We measure the net linear polarization fraction and find that high accretion rate "jet-disc" models are heavily depolarized and are therefore disfavoured. We show how Event Horizon Telescope measurements of the polarized "correlation length" over the image provide a model-independent upper limit on the strength of these Faraday effects, and constrain plasma properties like the electron temperature and magnetic field strength.

We study magnetohydrodynamic (MHD) effects arising in the collapse of magnetized, rotating, massive stellar cores to proto-neutron stars (PNSs). We perform axisymmetric numerical simulations in full general relativity with a hybrid equation of state. The formation and early evolution of a PNS are followed with a grid of $2500\ifmmode\times\else\texttimes\fi{}2500$ zones, which provides better resolution than in previous (Newtonian) studies. We confirm that significant differential rotation results even when the rotation of the progenitor is initially uniform. Consequently, the magnetic field is amplified both by magnetic winding and the magnetorotational instability (MRI). Even if the magnetic energy ${E}_{\mathrm{EM}}$ is much smaller than the rotational kinetic energy ${T}_{\mathrm{rot}}$ at the time of PNS formation, the ratio ${E}_{\mathrm{EM}}/{T}_{\mathrm{rot}}$ increases to 0.1\char21{}0.2 by the magnetic winding. Following PNS formation, MHD outflows lead to losses of rest mass, energy, and angular momentum from the system. The earliest outflow is produced primarily by the increasing magnetic stress caused by magnetic winding. The MRI amplifies the poloidal field and increases the magnetic stress, causing further angular momentum transport and helping to drive the outflow. After the magnetic field saturates, a nearly stationary, collimated magnetic field forms near the rotation axis and a Blandford-Payne\char21{}type outflow develops along the field lines. These outflows remove angular momentum from the PNS at a rate given by $\stackrel{\ifmmode \dot{}\else \textperiodcentered \fi{}}{J}\ensuremath{\sim}\ensuremath{\eta}{E}_{\mathrm{EM}}{C}_{B}$, where $\ensuremath{\eta}$ is a constant of order $\ensuremath{\sim}0.1$ and ${C}_{B}$ is a typical ratio of poloidal to toroidal field strength. As a result, the rotation period quickly increases for a strongly magnetized PNS until the degree of differential rotation decreases. Our simulations suggest that rapidly rotating, magnetized PNSs may not give rise to rapidly rotating neutron stars.

We study magnetic field steepening due to ambipolar diffusion in protoplanetary disk environments and draw the following conclusions. Current sheets are generated in magnetically active regions of the disk where the ionization fraction is high enough for the magnetorotational instability to operate. In late stages of solar nebula evolution, the surface density is expected to decrease and dust grains are expected to gravitationally settle to the midplane. If the local dust-to-gas mass ratio near the midplane is increased above cosmic abundances by factors ≳103, current sheets reach high enough temperatures to melt millimeter-sized dust grains and hence may provide the mechanism to form meteoritic chondrules. In addition, these current sheets possibly explain the near-infrared excesses observed in spectral energy distributions (SEDs) of young stellar objects. Direct imaging of protoplanetary disks via a nulling interferometer or, in the future, a multiband, adaptive optics coronagraph can test this hypothesis.

We study high-energy neutrino and cosmic-ray (CR) emission from the cores of low-luminosity active galactic nuclei (LLAGN). In LLAGN, the thermalization of particles is expected to be incomplete in radiatively inefficient accretion flows (RIAF), allowing the existence of non-thermal particles. In this work, assuming stochastic particle acceleration due to turbulence in RIAFs, we solve the Fokker–Planck equation and calculate spectra of escaping neutrinos and CRs. The RIAF in LLAGN can emit CR protons with ≳ 10 ?> PeV energies and TeV–PeV neutrinos generated via pp and/or p γ ?> reactions. We find that, if ∼1% of the accretion luminosity is carried away by non-thermal ions, the diffuse neutrino intensity from the cores of LLAGN may be as high as E &ngr; 2 Φ &ngr; ∼ 3 × 10 − 8 GeV cm − 2 s − 1 sr − 1 ?> , which can be compatible with the observed IceCube data. This result does not contradict either of the diffuse gamma-ray background observed by Fermi or observed diffuse CR flux. Our model suggests that, although very-high-energy gamma-rays may not escape, radio-quiet active galactic nuclei with RIAFs can emit GeV gamma-rays, which could be used for testing the model. We also calculate the neutron luminosity from RIAFs of LLAGN, and discuss a strong constraint on the model of jet mass loading mediated by neutrons from the diffuse neutrino observation.