Search

Collective cell migration plays an essential role in vertebrate development, yet the extent to which dynamically changing microenvironments influence this phenomenon remains unclear. Observations of the distribution of the extracellular matrix (ECM) component fibronectin during the migration of loosely connected neural crest cells (NCCs) lead us to hypothesize that NCC remodeling of an initially punctate ECM creates a scaffold for trailing cells, enabling them to form robust and coherent stream patterns. We evaluate this idea in a theoretical setting by developing an individual-based computational model that incorporates reciprocal interactions between NCCs and their ECM. ECM remodeling, haptotaxis, contact guidance, and cell-cell repulsion are sufficient for cells to establish streams in silico, however additional mechanisms, such as chemotaxis, are required to consistently guide cells along the correct target corridor. Further model investigations imply that contact guidance and differential cell-cell repulsion between leader and follower cells are key contributors to robust collective cell migration by preventing stream breakage. Global sensitivity analysis and simulated gain- and loss-of-function experiments suggest that long-distance migration without jamming is most likely to occur when leading cells specialize in creating ECM fibers, and trailing cells specialize in responding to environmental cues by upregulating mechanisms such as contact guidance., Comment: 66 pages, 19 figures (of which 14 are supplementary)

The classical, relativistic thin-disk theory of Novikov and Thorne (NT) predicts a maximum accretion efficiency of 40% for an optically thick, radiatively efficient accretion disk around a maximally spinning black hole (BH). However, when a strong magnetic field is introduced to numerical simulations of thin disks, large deviations in efficiency are observed, in part due to mass and energy carried by jets and winds launched by the disk or BH spin. The total efficiency of accretion can be significantly enhanced beyond that predicted by NT but it has remained unclear how the radiative component is affected. In order to study the effect of a dynamically relevant large-scale magnetic field on radiatively efficient accretion, we have performed numerical 3D general relativistic - radiative - magnetohydroynamic (GRRMHD) simulations of a disk with scale height to radius ratio of $H/R~0.1$ around a moderately spinning BH (a=0.5) using the code HARMRAD. Our simulations are fully global and allow us to measure the jet, wind, and radiative properties of a magnetically arrested disk (MAD) that is kept thin via self-consistent transport of energy by radiation using the M1 closure scheme. Our fiducial disk is MAD out to a radius of ~16R_g and the majority of the total ~13\% efficiency of the accretion flow is carried by a magnetically driven wind. We find that the radiative efficiency is slightly suppressed compared to NT, contrary to prior MAD GRMHD simulations with an ad hoc cooling function, but it is unclear how much of the radiation and thermal energy trapped in the outflows could ultimately escape., Comment: Accepted to be published in MNRAS, 15 pages and 13 figures

Magnetic fields are believed to drive accretion and relativistic jets in black hole accretion systems, but the magnetic-field structure that controls these phenomena remains uncertain. We perform general relativistic (GR) polarized radiative transfer of time-dependent three-dimensional GR magnetohydrodynamical (MHD) simulations to model thermal synchrotron emission from the Galactic Center source Sagittarius A$^\ast$ (Sgr A*). We compare our results to new polarimetry measurements by the Event Horizon Telescope (EHT) and show how polarization in the visibility (Fourier) domain distinguishes and constrains accretion flow models with different magnetic field structures. These include models with small-scale fields in disks driven by the magnetorotational instability (MRI) as well as models with large-scale ordered fields in magnetically-arrested disks (MAD). We also consider different electron temperature and jet mass-loading prescriptions that control the brightness of the disk, funnel-wall jet, and Blandford-Znajek-driven funnel jet. Our comparisons between the simulations and observations favor models with ordered magnetic fields near the black hole event horizon in Sgr A*, although both disk- and jet-dominated emission can satisfactorily explain most of the current EHT data. We show that stronger model constraints should be possible with upcoming circular polarization and higher frequency ($349~{\rm GHz}$) measurements., Comment: Submitted to ApJ, 20 pages, 19 figures, 1 table

One of the puzzles associated with tidal disruption event candidates (TDEs) is that there is a dichotomy between the color temperatures of ${\rm few}\times 10^4$~K for TDEs discovered with optical and UV telescopes, and the color temperatures of ${\rm few}\times 10^5 - 10^6$~K for TDEs discovered with X-ray satellites. Here we propose that high-temperature TDEs are produced when the tidal debris of a disrupted star self-intersects relatively close to the supermassive black hole, in contrast to the more distant self-intersection that leads to lower color temperatures. In particular, we note from simple ballistic considerations that greater apsidal precession in an orbit is the key to closer self-intersection. Thus larger values of $\beta$, the ratio of the tidal radius to the pericenter distance of the initial orbit, are more likely to lead to higher temperatures of more compact disks which are super-Eddington and geometrically and optically thick. For a given star and $\beta$, apsidal precession also increases for larger black hole masses, but larger black hole masses imply a lower temperature at the Eddington luminosity. Thus the expected dependence of the temperature on the mass of the black hole is non-monotonic. We find that in order to produce a soft X-ray temperature TDE, a deep plunging stellar orbit with $\beta> 3$ is needed and a black hole mass of $\lesssim 5\times 10^6 M_\odot$ is favored. Although observations of TDEs are comparatively scarce and are likely dominated by selection effects, it is encouraging that both expectations are consistent with current data., Comment: 10 pages, 4 figures, accepted by the Astrophysical Journal Letters

Accretion discs and black holes (BHs) have angular momenta that are generally misaligned with respect to each other, which can lead to warps in the discs and bends in any jets produced. We consider a disc that is misaligned at large radii and torqued by Lense-Thirring (LT) precession and a Blandford-Znajek (BZ) jet torque. We consider a variety of disc states that include radiatively inefficient thick discs, radiatively efficient thin discs, and super-Eddington accretion discs. The magnetic field strength of the BZ jet is chosen as either from standard equipartition arguments or from magnetically arrested disc (MAD) simulations. We show that standard thin accretion discs can reach spin-disc alignment out to large radii long before LT would play a role, as caused by the slow infall time that gives even a weak BZ jet time to align the disc. We show that geometrically thick radiatively inefficient discs and super-Eddington discs in the MAD state reach disc-spin alignment near the black hole when density profiles are shallow as in magnetohydrodynamical simulations, while the BZ jet aligns discs with steep density profiles (as in advection-dominated accretion flows) with the BH spin out to larger radii. Our results imply that the BZ jet torque should affect the cosmological evolution of BH spin magnitude and direction, BH spin measurements in active galactic nuclei and X-ray binaries, and the interpretations for Event Horizon Telescope observations of discs or jets in strong-field gravity regimes., Comment: 11 pages, 3 figures, submitted to MNRAS

High-frequency quasi-periodic oscillations (QPOs) appear in general-relativistic magnetohydrodynamic simulations of magnetically choked accretion flows around rapidly rotating black holes (BHs). We perform polarized radiative transfer calculations with ASTRORAY code to explore the manifestations of these QPOs for SgrA*. We construct a simulation-based model of a radiatively inefficient accretion flow and find model parameters by fitting the mean polarized source spectrum. The simulated QPOs have a total sub-mm flux amplitude up to 5% and a linearly polarized flux amplitude up to 2%. The oscillations reach high levels of significance 10-30\sigma and high quality factors Q~5. The oscillation period T~100M~35min corresponds to the rotation period of the BH magnetosphere that produces a trailing spiral in resolved disk images. The total flux signal is significant over noise for all tested frequencies 87GHz, 230GHz, and 857GHz and inclination angles 10deg, 37deg, and 80deg. The non-detection in the 230GHz Sub-Millimeter Array light curve is consistent with a low signal level and a low sampling rate. The presence of sub-mm QPOs in SgrA* will be better tested with the Atacama Large Millimeter Array., Comment: 7 pages, 5 figures, accepted to ApJ Letters, a movie of QPOs can be found at http://www.youtube.com/watch?v=NIMOFCtQDxQ

Recent advances in general relativistic magnetohydrodynamic modeling of jets offer unprecedented insights into the inner workings of accreting black holes that power the jets in active galactic nuclei (AGN) and other accretion systems. I will present the results of recent studies that determine spin-dependence of jet power and discuss the implications for the AGN radio loud/quiet dichotomy and recent observations of high jet power in a number of AGN., Comment: 8 pages, 4 figures. For an associated movie, see http://youtu.be/nRGCNaWST5Q . To appear in the proceedings of the conference "The Central Kiloparsec in Galactic Nuclei: Astronomy at High Angular Resolution 2011", open access Journal of Physics: Conference Series (JPCS), published by IOP Publishing

We briefly summarize the method of simulating Sgr A* polarized sub-mm spectra from the accretion flow and fitting the observed spectrum. The dynamical flow model is based on three-dimensional general relativistic magneto hydrodynamic simulations. Fully self-consistent radiative transfer of polarized cyclo-synchrotron emission is performed. We compile a mean sub-mm spectrum of Sgr A* and fit it with the mean simulated spectra. We estimate the ranges of inclination angle theta=42-75deg, mass accretion rate Mdot=(1.4-7.0)*10^{-8}Msun/yr, and electron temperature Te=(3-4)*10^{10}K at 6M. We discuss multiple caveats in dynamical modeling, which must be resolved to make further progress., Comment: 2 pages, no figures, for Proceedings of IAU Symposium 290 "Feeding compact objects: Accretion on all scales", eds C.M. Zhang, T. Belloni, M. Mendez & S.N. Zhang

Unconfined relativistic outflows from rotating, magnetized compact objects are often well-modeled by assuming the field geometry is approximately a split-monopole at large radii. Earlier work has indicated that such an unconfined flow has an inefficient conversion of magnetic energy to kinetic energy. This has led to the conclusion that ideal magnetohydrodynamical (MHD) processes fail to explain observations of, e.g., the Crab pulsar wind at large radii where energy conversion appears efficient. In addition, as a model for astrophysical jets, the monopole field geometry has been abandoned in favor of externally confined jets since the latter appeared to be generically more efficient jet accelerators. We perform time-dependent axisymmetric relativistic MHD simulations in order to find steady state solutions for a wind from a compact object endowed with a monopole field geometry. Our simulations follow the outflow for 10 orders of magnitude in distance from the compact object, which is large enough to study both the initial "acceleration zone" of the magnetized wind as well as the asymptotic "coasting zone." We obtain the surprising result that acceleration is actually {\it efficient} in the polar region, which develops a jet despite not being confined by an external medium. Our models contain jets that have sufficient energy to account for moderately energetic long and short gamma-ray burst (GRB) events (~10^{51}--10^{52} erg), collimate into narrow opening angles (opening half-angle \theta_j \approx 0.03 rad), become matter-dominated at large radii (electromagnetic energy flux per unit matter energy flux \sigma, Comment: Accepted to ApJ, 22 pages, 10 figures, uses emulateapj.cls. Changes in v2: new section in the Appendix

Long-duration gamma-ray bursts (GRBs) require an engine capable of driving a jet of plasma to ultrarelativistic bulk Lorentz factors of up to several hundred and into narrow opening angles of a few degrees. We use global axisymmetric stationary solutions of magnetically-dominated (force-free) ultrarelativistic jets to test whether the popular magnetic-driving paradigm can generate the required Lorentz factors and opening angles. Our global solutions are obtained via time-dependent relativistic ideal magnetodynamical numerical simulations which follow the jet from the central engine to beyond six orders of magnitude in radius. Our model is primarily motivated by the collapsar model, in which a jet is produced by a spinning black hole or neutron star and then propagates through a massive stellar envelope. We find that the size of the presupernova progenitor star and the radial profile of pressure inside the star determine the terminal Lorentz factor and opening angle of the jet. At the radius where the jet breaks out of the star, our well-motivated fiducial model generates a Lorentz factor $\gamma\sim 400$ and a half-opening angle $\theta_j\sim 2^\circ$, consistent with observations of many long-duration GRBs. Other models with slightly different parameters give $\gamma$ in the range 100 to 5000 and $\theta_j$ from $0.1^\circ$ to $10^\circ$, thus reproducing the range of properties inferred for GRB jets. A potentially observable feature of some of our solutions is that the maximum Poynting flux in the jet is found at $\theta \sim \theta_j$ with the jet power concentrated in a hollow cone, while the maximum in the Lorentz factor occurs at an angle $\theta$ substantially smaller than $\theta_j$ also in a hollow cone. [abridged], Comment: 25 pages, 12 figures, Accepted to MNRAS. Changes in v2: new figure and 2 new sections in the Appendix

The formation and large-scale propagation of Poynting-dominated jets produced by accreting, rapidly rotating black hole systems are studied by numerically integrating the general relativistic magnetohydrodynamic equations of motion to follow the self-consistent interaction between accretion disks and black holes. This study extends previous similar work by studying jets till $t\approx 10^4GM/c^3$ out to $r\approx 10^4GM/c^2$, by which the jet is super- fast magnetosonic and moves at a lab-frame bulk Lorentz factor of $\Gamma\sim 10$ with a maximum terminal Lorentz factor of $\Gamma_\infty\lesssim 10^3$. The radial structure of the Poynting-dominated jet is piece-wise self-similar, and fits to flow quantities along the field line are provided. Beyond the \alf surface at $r\sim 10$--$100GM/c^2$, the jet becomes marginally unstable to (at least) current-driven instabilities. Such instabilities drive shocks in the jet that limit the efficiency of magnetic acceleration and collimation. These instabilities also induce jet substructure with $3\lesssim\Gamma\lesssim 15$. The jet is shown to only marginally satisfy the necessary and sufficient conditions for kink instability, so this may explain how astrophysical jets can extend to large distances without completely disrupting. At large distance, the jet angular structure is Gaussian-like (or uniform within the core with sharp exponential wings) with a half-opening angle of $\approx 5^\circ$ and there is an extended component out to $\approx 27^\circ$. Unlike in some hydrodynamic simulations, the environment is found to play a negligible role in jet structure, acceleration, and collimation as long as the ambient pressure of the surrounding medium is small compared to the magnetic pressure in the jet., Comment: 25 pages, 11 figures, accepted for publication in MNRAS Main Journal (astro-ph/0506369 with pair creation under revision to focus on GRBs)

In paper I, we showed that time-dependent general relativistic magnetohydrodynamic (GRMHD) numerical models of accretion disks, although being highly turbulent, have surprisingly simple electromagnetic properties. In particular, the toroidal current density in the disk takes the form $dI_��/dr \propto r^{-5/4}$. Guided by this simplicity, we use a time-dependent general relativistic force-free electrodynamics (GRFFE) code to study an idealized problem in which the accretion disk is replaced by an infinitely thin rotating equatorial current sheet. We consider both an $r^{-5/4}$ current profile and an $r^{-1}$ profile, the latter corresponding to the paraboloidal model of Blandford & Znajek (1977). The force-free magnetosphere we obtain with the $r^{-5/4}$ current sheet matches remarkably well to the Poynting-dominated jet seen in GRMHD numerical models. By comparing to the non-rotating force-free model studied in paper I, rotation is seen to lead to mild decollimation of the jet suggesting that hoop-stress forces nearly cancel centrifugal forces. In order to study the process that generates the corona and disk wind and destroys the ordered field in the corona in GRMHD numerical models, the force-free field with the $r^{-5/4}$ current distribution is embedded in an accretion disk and followed in a GRMHD simulation. Reconnection and magnetic stresses contribute to a magnetized, thermal wind without the aid of an ordered field threading the disk., 20 pages, 17 figures, accepted to MNRAS

We consider a self-similar force-free wind flowing out of an infinitely thin disk located in the equatorial plane. On the disk plane, we assume that the magnetic stream function $P$ scales as $P\propto R^\nu$, where $R$ is the cylindrical radius. We also assume that the azimuthal velocity in the disk is constant: $v_\phi = Mc$, where $M, Comment: 21 pages, 15 figures, accepted to MNRAS

The force-free limit of magnetohydrodynamics (MHD) is often a reasonable approximation to model black hole and neutron star magnetospheres. We describe a general relativistic force-free (GRFFE) formulation that allows general relativistic magnetohydrodynamic (GRMHD) codes to directly evolve the GRFFE equations of motion. Established, accurate, and well-tested conservative GRMHD codes can simply add a new inversion piece of code to their existing code, while continuing to use all the already-developed facilities present in their GRMHD code. We show how to enforce the $\mathbf{E}\cdot\mathbf{B}=0$ constraint and energy conservation, and we introduce a simplified general model of the dissipation of the electric field to enforce the $B^2-E^2>0$ constraint. We also introduce a simplified yet general method to resolve current sheets, without much reconnection, over many dynamical times. This formulation is incorporated into an existing GRMHD code (HARM), which is demonstrated to give accurate and robust GRFFE results for Minkowski and black hole space-times., 14 pages, 4 figures, accepted for publication in MNRAS Main Journal

Some active galactic nuclei, microquasars, and gamma ray bursts may be powered by the electromagnetic braking of a rapidly rotating black hole. We investigate this possibility via axisymmetric numerical simulations of a black hole surrounded by a magnetized plasma. The plasma is described by the equations of general relativistic magnetohydrodynamics, and the effects of radiation are neglected. The evolution is followed for $2000 G M/c^3$, and the computational domain extends from inside the event horizon to typically $40 G M/c^2$. We compare our results to two analytic steady state models, including the force-free magnetosphere of Blandford & Znajek. Along the way we present a self-contained rederivation of the Blandford-Znajek model in Kerr-Schild (horizon penetrating) coordinates. We find that (1) low density polar regions of the numerical models agree well with the Blandford-Znajek model; (2) many of our models have an outward Poynting flux on the horizon in the Kerr-Schild frame; (3) none of our models have a net outward energy flux on the horizon; and (4) one of our models, in which the initial disk has net magnetic flux, shows a net outward angular momentum flux on the horizon. We conclude with a discussion of the limitations of our model, astrophysical implications, and problems to be addressed by future numerical experiments., 45 pages, 11 Postscript figures, 4 tables, uses subfigure.sty, ApJ, in press

We consider astrophysical maser radiation that is created in the presence of mildly supersonic, magnetohydrodynamic (MHD) turbulence. The focus is on the OH masers for which the magnetic field is strong enough that the separations of the Zeeman components are greater than the spectral linebreadths. A longstanding puzzle has been the absence of the Zeeman $\pi $ components and the high circular polarization in the observed spectra of these masers. We first argue that the elongation of eddies along the field that has recently been recognized in MHD turbulence will enhance the optical depth parallel to the magnetic field in comparison with that perpendicular to the magnetic field. We then simulate maser emission with a numerical model of MHD turbulence to demonstrate quantitatively how the intensities of the linearly polarized $\pi $ components are suppressed and the intensities of the nearly circularly polarized $\sigma $ components are enhanced. This effect is also generic in the sense that most spectral lines in MHD turbulence with Mach number $M\sim 1$ should have larger optical depths parallel to the magnetic field than perpendicular. The effect is reduced considerably when $M, Comment: accepted for Astrophysical Journal

We consider a subset of the physical processes that determine the spin j = a/M of astrophysical black holes. These include: (1) Initial conditions. Recent models suggest that the collapse of supermassive stars are likely to produce black holes with j ~ 0.7. (2) Major mergers. The outcome of a nearly equal mass black hole-black hole merger is not yet known, but we review the current best guesses and analytic bounds. (3) Minor mergers. We recover the result of Blandford & Hughes that accretion of small companions with isotropically distributed orbital angular momenta results in spindown, with j ~ M^{-7/3}. (4) Accretion. We present new results from fully relativistic magnetohydrodynamic accretion simulations. These show that, at least for one sequence of flow models, spin equilibrium (dj/dt = 0) is reached for j ~ 0.9, far less than the canonical value 0.998 of Thorne that was derived in the absence of MHD effects. This equilibrium value may not apply to all accretion flows, particularly thin disks. Nevertheless, it opens the possibility that black holes that have grown primarily through accretion are not maximally rotating., 22 pp, 4 figures, accepted to ApJ

We report on a numerical study of viscous fluid accretion onto a black hole. The flow is axisymmetric and uses a pseudo-Newtonian potential to model relativistic effects near the event horizon. The numerical method is a variant of the ZEUS code. As a test of our numerical scheme, we are able to reproduce results from earlier, similar work by Igumenshchev and Abramowicz and Stone et al. We consider models in which mass is injected onto the grid as well as models in which an initial equilibrium torus is accreted. In each model we measure three ``eigenvalues'' of the flow: the accretion rate of mass, angular momentum, and energy. We find that the eigenvalues are sensitive to r_{in}, the location of the inner radial boundary. Only when the flow is always supersonic on the inner boundary are the eigenvalues insensitive to small changes in r_{in}. We also report on the sensitivity of the results to other numerical parameters., Comment: 14 pages, 4 figures, 2 tables, to appear in v573 n2 pt1 ApJ July 10, 2002

We study the motion of electrons in a single miniband of a semiconductor superlattice driven by THz electric field polarized along the growth direction. We work in the semiclassical balance-equation model, including different elastic and inelastic scattering rates, and incorporating the self-consistent electric field generated by electron motion. We explore regions of complex dynamics, which can include chaotic behaviour and symmetry-breaking. We estimate the magnitudes of dc current and dc voltage that spontaneously appear in regions of broken-symmetry for parameters characteristic of modern semiconductor superlattices. This work complements PRL 80(1998)2669 [ cond-mat/9709026 ]., Comment: 4 pages, 3 figures, RevTEX, EPSF