Newborn black holes in collapsing massive stars can be accompanied by a fallback disk. The accretion rate is typically super-Eddington and strong disk outflows are expected. Such outflows could be directly observed in some failed explosions of compac
t (blue supergiants or Wolf-Rayet stars) progenitors, and may be more common than long-duration gamma-ray bursts. Using an analytical model, we show that the fallback disk outflows produce blue UV-optical transients with a peak bolometric luminosity of ~10^(42-43) erg s^-1 (peak R-band absolute AB magnitudes of -16 to -18) and an emission duration of ~ a few to ~ 10 days. The spectra are likely dominated intermediate mass elements, but will lack much radioactive nuclei and iron-group elements. The above properties are broadly consistent with some of the rapid blue transients detected by Pan-STARRS and PTF. This scenario can be distinguished from alternative models using radio observations within a few years after the optical peak.
We quantify the stellar abundances of neutron-rich r-process nuclei in cosmological zoom-in simulations of a Milky Way-mass galaxy from the Feedback In Realistic Environments project. The galaxy is enriched with r-process elements by binary neutron s
tar (NS) mergers and with iron and other metals by supernovae. These calculations include key hydrodynamic mixing processes not present in standard semi-analytic chemical evolution models, such as galactic winds and hydrodynamic flows associated with structure formation. We explore a range of models for the rate and delay time of NS mergers, intended to roughly bracket the wide range of models consistent with current observational constraints. We show that NS mergers can produce [r-process/Fe] abundance ratios and scatter that appear reasonably consistent with observational constraints. At low metallicity, [Fe/H]<-2, we predict there is a wide range of stellar r-process abundance ratios, with both supersolar and subsolar abundances. Low-metallicity stars or stars that are outliers in their r-process abundance ratios are, on average, formed at high redshift and located at large galactocentric radius. Because NS mergers are rare, our results are not fully converged with respect to resolution, particularly at low metallicity. However, the uncertain rate and delay time distribution of NS mergers introduces an uncertainty in the r-process abundances comparable to that due to finite numerical resolution. Overall, our results are consistent with NS mergers being the source of most of the r-process nuclei in the Universe.
It is typically assumed that radiation pressure driven winds are accelerated to an asymptotic velocity of V ~ v_esc, where v_esc is the escape velocity from the central source. We note that this is not the case for dusty shells and clouds. Instead, i
f the shell or cloud is initially optically-thick to the UV emission from the source of luminosity L, then there is a significant boost in V that reflects the integral of the momentum absorbed as it is accelerated. For shells reaching a generalized Eddington limit, we show that V ~ (4R_UV L/M_sh c)^1/2, in both point-mass and isothermal-sphere potentials, where R_UV is the radius where the shell becomes optically-thin to UV photons, and M_sh is the mass of the shell. The asymptotic velocity significantly exceeds v_esc for typical parameters, and can explain the ~1000-2000km/s outflows observed from rapidly star-forming galaxies and active galactic nuclei if the surrounding halo has low gas density. Similarly fast outflows from massive stars can be accelerated on few - 10^3 yr timescales. These results carry over to clouds that subtend only a small fraction of the solid angle from the source of radiation and that expand as a consequence of their internal sound speed. We further consider the dynamics of shells that sweep up a dense circumstellar or circumgalactic medium. We calculate the momentum ratio Mdot v/(L/c) in the shell limit and show that it can only significantly exceed ~2 if the effective optical depth of the shell to re-radiated FIR photons is much larger than unity. We discuss simple prescriptions for the properties of galactic outflows for use in large-scale cosmological simulations. We also briefly discuss applications to the dusty ejection episodes of massive stars, the disruption of giant molecular clouds, and AGN.
We present estimates of magnetic field strengths in the interstellar media of starburst galaxies derived from measurements of Zeeman splitting associated with OH megamasers. The results for eight galaxies with Zeeman detections suggest that the magne
tic energy density in the interstellar medium of starburst galaxies is comparable to their hydrostatic gas pressure, as in the Milky Way. We discuss the significant uncertainties in this conclusion, and possible measurements that could reduce these uncertainties. We also compare the Zeeman splitting derived magnetic field estimates to magnetic field strengths estimated using synchrotron fluxes and assuming that the magnetic field and cosmic rays have comparable energy densities, known as the minimum energy argument. We find that the minimum energy argument systematically underestimates magnetic fields in starburst galaxies, and that the conditions that would be required to produce agreement between the minimum energy estimate and the Zeeman derived estimate of interstellar medium magnetic fields are implausible. The conclusion that magnetic fields in starburst galaxies exceed the minimum energy magnetic fields is consistent with starburst galaxies adhering to the linearity of the FIR-radio correlation.
In the intracluster medium (ICM) of galaxy clusters, heat and momentum are transported almost entirely along (but not across) magnetic field lines. We perform the first fully self-consistent Braginskii-MHD simulations of galaxy clusters including bot
h of these effects. Specifically, we perform local and global simulations of the magnetothermal instability (MTI) and the heat-flux-driven buoyancy instability (HBI) and assess the effects of viscosity on their saturation and astrophysical implications. We find that viscosity has only a modest effect on the saturation of the MTI. As in previous calculations, we find that the MTI can generate nearly sonic turbulent velocities in the outer parts of galaxy clusters, although viscosity somewhat suppresses the magnetic field amplification. At smaller radii in cool-core clusters, viscosity can decrease the linear growth rates of the HBI. However, it has less of an effect on the HBIs nonlinear saturation, in part because three-dimensional interchange motions (magnetic flux tubes slipping past each other) are not damped by anisotropic viscosity. In global simulations of cool core clusters, we show that the HBI robustly inhibits radial thermal conduction and thus precipitates a cooling catastrophe. The effects of viscosity are, however, more important for higher entropy clusters. We argue that viscosity can contribute to the global transition of cluster cores from cool-core to non cool-core states: additional sources of intracluster turbulence, such as can be produced by AGN feedback or galactic wakes, suppress the HBI, heating the cluster core by thermal conduction; this makes the ICM more viscous, which slows the growth of the HBI, allowing further conductive heating of the cluster core and a transition to a non cool-core state.
We use three-dimensional MHD simulations with anisotropic thermal conduction to study turbulence due to the magnetothermal instability (MTI) in the intracluster medium (ICM) of galaxy clusters. The MTI grows on timescales of ~1 Gyr and is capable of
driving vigorous, sustained turbulence in the outer parts of galaxy clusters if the temperature gradient is maintained in spite of the rapid thermal conduction. If this is the case, turbulence due to the MTI can provide up to 5-30% of the pressure support beyond r_500 in galaxy clusters, an effect that is strongest for hot, massive clusters. The turbulence driven by the MTI is generally additive to other sources of turbulence in the ICM, such as that produced by structure formation. This new source of non-thermal pressure support reduces the observed Sunyaev-Zeldovich (SZ) signal and X-ray pressure gradient for a given cluster mass and introduces a cluster mass and temperature gradient-dependent bias in SZ and X-ray mass estimates of clusters. This additional physics may also need to be taken into account when estimating the matter power spectrum normalization, sigma-8, through simulation templates from the observed amplitude of the SZ power spectrum.
We study the excitation and damping of tides in close binary systems, accounting for the leading order nonlinear corrections to linear tidal theory. These nonlinear corrections include two distinct effects: three-mode nonlinear interactions and nonli
near excitation of modes by the time-varying gravitational potential of the companion. This paper presents the formalism for studying nonlinear tides and studies the nonlinear stability of the linear tidal flow. Although the formalism is applicable to binaries containing stars, planets, or compact objects, we focus on solar type stars with stellar or planetary companions. Our primary results include: (1) The linear tidal solution often used in studies of binary evolution is unstable over much of the parameter space in which it is employed. More specifically, resonantly excited gravity waves are unstable to parametric resonance for companion masses M > 10-100 M_Earth at orbital periods P = 1-10 days. The nearly static equilibrium tide is, however, parametrically stable except for solar binaries with P < 2-5 days. (2) For companion masses larger than a few Jupiter masses, the dynamical tide causes waves to grow so rapidly that they must be treated as traveling waves rather than standing waves. (3) We find a novel form of parametric instability in which a single parent wave excites a very large number of daughter waves (N = 10^3[P / 10 days]) and drives them as a single coherent unit with growth rates that are ~N times faster than the standard three wave parametric instability. (4) Independent of the parametric instability, tides excite a wide range of stellar p-modes and g-modes by nonlinear inhomogeneous forcing; this coupling appears particularly efficient at draining energy out of the dynamical tide and may be more important than either wave breaking or parametric resonance at determining the nonlinear dissipation of the dynamical tide.
Feedback from massive stars is believed to play a critical role in shaping the galaxy mass function, the structure of the interstellar medium (ISM), and the low efficiency of star formation, but the exact form of the feedback is uncertain. In this pa
per, the first in a series, we present and test a novel numerical implementation of stellar feedback resulting from momentum imparted to the ISM by radiation, supernovae, and stellar winds. We employ a realistic cooling function, and find that a large fraction of the gas cools to <100K, so that the ISM becomes highly inhomogeneous. Despite this, our simulated galaxies reach an approximate steady state, in which gas gravitationally collapses to form giant molecular clouds (GMCs), dense clumps, and stars; subsequently, stellar feedback disperses the GMCs, repopulating the diffuse ISM. This collapse and dispersal cycle is seen in models of SMC-like dwarfs, the Milky-Way, and z~2 clumpy disk analogues. The simulated global star formation efficiencies are consistent with the observed Kennicutt-Schmidt relation. Moreover, the star formation rates are nearly independent of the numerically imposed high-density star formation efficiency, density threshold, and density scaling. This is a consequence of the fact that, in our simulations, star formation is regulated by stellar feedback limiting the amount of very dense gas available for forming stars. In contrast, in simulations without stellar feedback, i.e. under the action of only gravity and gravitationally-induced turbulence, the ISM experiences runaway collapse to very high densities. In these simulations without feedback, the global star formation rates exceed observed galactic star formation rates by 1-2 orders of magnitude, demonstrating that stellar feedback is crucial to the regulation of star formation in galaxies.
The emission from Sgr A*, the supermassive black hole in the Galactic Center, shows order of magnitude variability (flares) a few times a day that is particularly prominent in the near-infrared (NIR) and X-rays. We present a time-dependent model for
these flares motivated by the hypothesis that dissipation of magnetic energy powers the flares. We show that episodic magnetic reconnection can occur near the last stable circular orbit in time-dependent magnetohydrodynamic simulations of black hole accretion - the timescales and energetics of these events are broadly consistent with the flares from Sgr A*. Motivated by these results, we present a spatially one-zone time-dependent model for the electron distribution function in flares, including energy loss due to synchrotron cooling and adiabatic expansion. Synchrotron emission from transiently accelerated particles can explain the NIR/X-ray lightcurves and spectra of a luminous flare observed 4 April 2007. A significant decrease in the magnetic field strength during the flare (coincident with the electron acceleration) is required to explain the simultaneity and symmetry of the simultaneous lightcurves. Our models predict that the NIR and X-ray spectral indices differ by 0.5 and that there is only modest variation in the spectral index during flares. We also explore implications of this model for longer wavelength (radio-submm) emission seemingly associated with X-ray and NIR flares; we argue that a few hour decrease in the submm emission is a more generic consequence of large-scale magnetic reconnection than delayed radio emission from adiabatic expansion.
We compile observations of the surface mass density profiles of dense stellar systems, including globular clusters in the Milky Way and nearby galaxies, massive star clusters in nearby starbursts, nuclear star clusters in dwarf spheroidals and late-t
ype disks, ultra-compact dwarfs, and galaxy spheroids spanning the range from low-mass cusp bulges and ellipticals to massive core ellipticals. We show that in all cases the maximum stellar surface density attained in the central regions of these systems is similar, Sigma_max ~ 10^11 M_sun/kpc^2 (~20 g/cm^2), despite the fact that the systems span 7 orders of magnitude in total stellar mass M_star, 5 in effective radius R_e, and have a wide range in effective surface density M_star/R_e^2. The surface density limit is reached on a wide variety of physical scales in different systems and is thus not a limit on three-dimensional stellar density. Given the very different formation mechanisms involved in these different classes of objects, we argue that a single piece of physics likely determines Sigma_max. The radiation fields and winds produced by massive stars can have a significant influence on the formation of both star clusters and galaxies, while neither supernovae nor black hole accretion are important in star cluster formation. We thus conclude that feedback from massive stars likely accounts for the observed Sigma_max, plausibly because star formation reaches an Eddington-like flux that regulates the growth of these diverse systems. This suggests that current models of galaxy formation, which focus on feedback from supernovae and active galactic nuclei, are missing a crucial ingredient.
