Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: magnetic fields are now generally acknowledged to be essential to accretion @xcite . amplified exponentially from weak seed fields by the magneto - rotational instability ( mri ) , they are stretched in a consistent direction so as to produce a net stress @xmath0 capable of carrying enough angular momentum outward through the disk to support sizable accretion rates . because the nonlinear state of the mri is fully - developed magnetohydrodynamic ( mhd ) turbulence , it is often thought ( or at least hoped ) that its properties have little to do with the character of the seed magnetic field brought to the disk with the accretion flow .    however , there are certain properties of the magnetic field that are preserved despite the development of nonlinear turbulence . because magnetic fields are divergence - free and the plasma in accretion flows has very little electrical resistivity , the global topology of the field persists , no matter what happens . it is therefore natural to ask whether any interesting properties of the accretion flow may depend on the global field topology . for example , could the field topology alter the relation between surface density and accretion rate ? or how the stress varies through the marginally stable region ? in addition to swallowing matter , accreting black holes often expel relativistic jets . although the mechanisms by which these jets are launched have long been the subject of investigation and speculation , the current consensus focuses on a combination of large - scale magnetic fields and the rotation of either the central hole , the accretion disk , or both . two specific mhd models have received the greatest attention : the blandford - znajek mechanism @xcite and the blandford - payne mechanism @xcite . in the blandford - payne model a large - scale vertical magnetic field is anchored in the disk , rotating with the orbital frequency . above and below the disk magnetic tension dominates , and the field lines force the matter to rotate with this same frequency . if the fieldlines are angled outward sufficiently with respect to the disk , there can be a net outward force , accelerating the matter along the rotating fieldlines . the blandford - znajek model also depends upon rotating field lines , but the source of that rotation ( and the source of the power in the jet ) is the black hole rather than the disk : any field lines connecting the immediate vicinity of the black hole to infinity are forced to rotate by frame - dragging , so that poynting flux carries off the energy .    given the strong dependence of both the blandford - znajek and the blandford - payne mechanisms on the nature of the magnetic field in the disk and near the black hole , it is equally natural to inquire whether the topology of the field embedded in the accretion flow has any influence on their operation . it is the goal of this paper to begin to investigate these questions , both as they relate to accretion disk dynamics and as they relate to jet properties . unfortunately , it is difficult even to scratch the surface by purely analytic techniques ( although there have been some attempts using highly - simplified models : * ? ? ? * ; * ? ? ? on the other hand , direct numerical simulations provide a means by which aspects of these problems can be investigated in detail . @xcite made a first exploration of the dependence of accretion properties on field topology in a highly - simplified approximation : newtonian dynamics in a pseudo - newtonian potential . they found that a purely toroidal initial field led to a systematically smaller accretion rate than one whose initial field was dipolar . more recently , fully general relativistic mhd ( grmhd ) simulations have become possible @xcite . these works have focused primarily on a single initial field structure : a set of nested dipole field loops contained entirely within the initial gas torus . other field geometries have been examined , but only in rather limited fashion . for example , @xcite simulated the evolution of a torus with an initial toroidal field around a schwarzschild black hole and noted the lack of any resulting jet , but did not consider why no jet formed . @xcite , using axisymmetric simulations , compared the electromagnetic luminosity of the standard dipole jet with that resulting from a few alternate magnetic field topologies . in particular , they studied cases in which the initial field was arranged in loops of alternating sign arranged either vertically or horizontally within the initial torus . the principal effect they reported was a factor of 23 reduction in the ratio of the time - averaged jet poynting flux to the mean mass accretion rate . they also noted that an initial field that is purely vertical led to significant enhancements in both the mean mass accretion rate and the time - averaged jet poynting flux . axisymmetric simulations of tori threaded by net vertical field led @xcite to concur that such a field can greatly enhance jet power . on the basis of these results , @xcite proposed that black holes inevitably become threaded by an organized vertical field if the initial disk contains a large - scale poloidal field with few parity changes . the intriguing character of these results motivates a more systematic study , one that employs three - dimensional simulations , considers the accretion flow as well as the jet , and investigates time - dependence . the question immediately arises : how can one best attempt to describe the very complex mixture of topologies likely to be found in real accretion flows ? the magnetic field near a black hole can be thought of as simply what the accretion process itself has self - consistently brought there from whatever source supplies it . in order to gain understanding of what , in realistic circumstances , must be a truly complex combination of structures , in this paper we do not try to construct a single  realistic " field geometry ; instead , we study the effects due to individual members of a  basis set " of field topologies . our hope is that by understanding the action of pure types , we can gain some understanding of how realistic mixtures behave .    to be specific , we will begin with a carefully chosen set of different kinds of closed field loops that we expect to be generic . closed field loops have the advantage of requiring many fewer free parameters to describe them fully within a simulation than do field geometries with net flux . net flux geometries may well exist in nature and have important effects ; we will consider them in the second stage of this project . we begin with the conjecture that , at a qualitative level , the different sorts of closed field loop behavior can all be represented by a combination of those due to : dipolar loops , or fieldlines initially entirely poloidal and which cross the equatorial plane at well - separated radii ; toroidal loops , which encircle the black hole ; and quadrupolar loops , which ( again , in the initial state ) are poloidal , but close without crossing the equatorial plane . for each of these cases , we have run a lengthy three - dimensional grmhd simulation whose results we hope characterize the accretion flow structure created by that field topology . we have also run a series of high resolution axisymmetric simulations of each of the initially poloidal topologies ( i.e. dipole and quadrupole ) . in addition , in order to understand better the implications of the finite lifetimes of closed loops , we have performed an axisymmetric simulation whose initial magnetic field structure was a series of dipolar loops , each contained within a span of radii narrow enough that the difference in inflow times across a loop is comparable to the equilibration time for the magnetic field in the jet . the rest of this work is structured as follows . in  [ numerics ] we briefly review the numerical scheme employed to solve the equations of grmhd and give a detailed description of the initial conditions used to generate each topology . in  [ flowstructure ] we examine the late time accretion flow structure that evolves from each of the initial field topologies , paying attention to each of the structures ( main disk body , magnetized corona ) that were identified by @xcite as arising from the evolution of the dipole field topology . in  [ jets ] we examine the properties of the third structure identified by @xcite , the unbound outflow , and examine how the initial evolution of the field during the linear growth phase of the mri gives rise to the components of the accretion flow at late times . finally in  [ conclusion ] we summarize our results and highlight their importance for jets in astrophysical settings . in this work we undertake a study of both two- and three - dimensional simulations to investigate the influence of magnetic field topology on both the evolved accretion flow structure and relativistic jet formation . in doing so we continue a program of black hole accretion disk simulations begun in @xcite , @xcite , @xcite , @xcite and @xcite . three - dimensional simulations , such as those presented in the previous work , are the only means by which one can gain insight into the evolved accretion flow properties due to the fundamental restriction posed by the anti - dynamo theorem @xcite . further , toroidal field models can not be studied in axisymmetry , since the toroidal mri requires nonzero azimuthal wavevectors . on the other hand , because jet formation occurs promptly during the evolution of the flow @xcite , restrictions due to the anti - dynamo theorem are less critical for jet studies . in this case the higher resolution afforded by the two - dimensional simulations is particularly valuable . the simulation code we use is described in @xcite . this code solves the equations of ideal non - radiative mhd in the static kerr metric of a rotating black hole using boyer - lindquist coordinates . values are expressed in gravitational units @xmath1 with line element @xmath2 and signature @xmath3 . since we are focusing in this paper on the influence of the magnetic field topology , we fix the black hole spin to @xmath4 for all of the new simulations reported here . the relativistic fluid at each grid zone is described by its density @xmath5 , specific internal energy @xmath6 , @xmath7-velocity @xmath8 , and isotropic pressure @xmath9 . the relativistic enthalpy is @xmath10 . the pressure is related to @xmath5 and @xmath6 via the equation of state for an ideal gas , @xmath11 . the magnetic field is described by two sets of variables . the first is the constrained transport magnetic field @xmath12 f_{jk}$ ] , where @xmath13 $ ] is the completely anti - symmetric symbol , and @xmath14 are the spatial components of the electromagnetic field strength tensor . from the constrained transport ( ct ) magnetic field components , we derive the magnetic field four - vector , @xmath15 , and the magnetic field scalar , @xmath16 . the electromagnetic component of the stress - energy tensor is @xmath17 .    in all of the simulations , the initial condition for the matter is an isolated , hydrostatic torus orbiting near the black hole . the initial state of the three - dimensional simulations was chosen so as to be consistent with that used in earlier simulations in this series ( see * ? ? ? the matter s adiabatic index , @xmath18 was @xmath19 and the angular momentum distribution was slightly sub - keplerian , with a specific angular momentum at the inner edge of the torus ( located at @xmath20 m ) @xmath21 ( we define @xmath22 ) . for this choice of @xmath23 , the pressure maximum is at @xmath24 .    in the two - dimensional simulations , our fiducial torus was chosen to have a constant angular momentum distribution with @xmath25 , @xmath18 was 4/3 , and the inner edge was again located at @xmath20 m . this torus is significantly thicker than the nearly - keplerian torus used for the three - dimensional simulations ; its density scale height at the pressure maximum , measured along surfaces of constant radius , is a factor @xmath26 greater than that of the three - dimensional torus . this thicker configuration makes it easier to resolve the multiple field loop model , where several field loops must fit into a restricted space . this specific combination of @xmath27 and @xmath18 was also necessary to manipulate the shape of the torus so that there was an appropriate balance in magnetic flux between the different components of the multiple dipole loop simulation . in the quadrupolar field case , in order to study how the evolution of the field depends on its detailed geometry , we examined what happened as @xmath27 varied from 4.9 to 5.2 in steps of 0.05 . the initial torus is overlaid with a variety of magnetic field configurations , all with zero net flux . the initial poloidal magnetic field is determined from the curl of the four - vector potential , i.e. , @xmath28 with only @xmath29 . for a single dipolar loop , we set @xmath30 for a quadrupole , the azimuthal component of the vector potential is @xmath31 in both of these cases @xmath32 is set to zero where the density is less than @xmath33 ; this device confines the initial field to lie well inside the torus . in all the poloidal simulations , the initial field intensity ( i.e. , the parameter @xmath34 ) is determined by setting the volume - averaged ratio of gas to magnetic pressure @xmath35 . the vector potential for @xmath36 dipolar loops ( used only for the two - dimensional simulations ) is more complicated . it is given by @xmath37,\ ] ] where the total vector potential is the superposition of the ( non - relativistic ) vector potentials from @xmath36 current loops ( see * ? ? ? 182 , eqn . 5.37 ) . here @xmath38 is the current in the @xmath39th loop , located at radius @xmath40 in the equatorial plane . @xmath41 and @xmath42 are the complete elliptic integrals of the first and second kinds , and @xmath43 . the @xmath44 are chosen so that the difference between the inflow times of adjacent loops is approximately equal to one jet equilibration time , @xmath45@xmath46 . the signs of the @xmath47 alternate with increasing radius so that the senses of the dipolar field loops alternate . the magnitudes of the currents were chosen with the goal of ensuring that the evolution of the two radially innermost loops was governed by the field local to those loops . the purpose of this simulation was to study the interaction of the two inner loops as they accrete ; in that context , the only function of the two outer loops is to supply the outer part of the torus with magnetic field , so that it could accept the angular momentum brought to it from the true accretion flow . thus , the most appropriate outer field structure for this particular numerical experiment is one that has the least influence on dynamics inside the pressure maximum . after some experimentation , we found that this goal could be achieved only by giving the two outer loops relatively small currents . when that was done , the outer two loops moved outward with the outer portion of the disk and were quickly destroyed ( i.e. forced to reconnect ) as the second loop from the inside expanded toward them . we make no pretense that this choice is a fair representation of generic four - loop evolution , but it does support our study of what happens when two complete loops of opposite sense are brought into the inner disk one after the other .    for a purely toroidal initial field , all that is necessary is to set @xmath48 and determine @xmath49 by a choice of volume - averaged @xmath50 . here we use @xmath51 , as , in the purely toroidal case , accretion can occur only once mri - induced turbulence on the large scales has grown to sufficiently large amplitude , and long wavelength modes grow comparatively slowly . beginning with a comparatively strong field ensured that accretion activity started soon after the simulation began . after saturation , the field ( in relative terms ) actually weakened somewhat : the volume - averaged @xmath50 grew to @xmath52 ( see section [ diskbody ] ) . as before , the initial magnetic field is set to zero for densities below @xmath33 . we take the dipole simulation kdpg ( originally presented in * ? ? ? * where further discussion of our boundary conditions and grids may be found ) as our fiducial three - dimensional run . the three - dimensional quadrupole simulation is designated qdpa and the three - dimensional toroidal simulation is tdpa . kdpg and qdpa were run for @xmath53 in time , which corresponds to approximately @xmath54 orbits at the initial pressure maximum . the toroidal field model , tdpa , evolves initially at a slower rate ; it was run to @xmath55 in time . for each simulation , the time step @xmath56 was determined by the minimal light crossing time for a zone on the spatial grid and remained constant for the entire simulation @xcite . each of the three - dimensional simulations used @xmath57 @xmath58 grid zones . the radial grid extends from an inner boundary located at @xmath59 , which lies just outside the black hole event horizon , to an outer boundary located at @xmath60 . for the two poloidal topologies , the radial grid was graded using a hyperbolic cosine distribution ; for the toroidal case a logarithmic distribution was used . an outflow condition is applied at both the inner and outer radial boundary . the @xmath61-grid spans the range @xmath62 using an exponential distribution that concentrates zones near the equator . a reflecting boundary condition is used along the conical cutout surrounding the coordinate axis . the @xmath63-grid spans the quarter plane , @xmath64 , with periodic boundary conditions applied in @xmath63 . the use of this restricted angular domain significantly reduces the computational requirements of the simulation . the two - dimensional simulations were performed on a @xmath65 @xmath66 grid in which the radial cells are graded logarithmically , and the polar angle cells are concentrated toward the equator in the same way as for the three - dimensional simulations . the @xmath61 grid , however , goes much closer to the axis than in the three - dimensional simulations , with a cut - out of only @xmath67  radians . in addition to allowing higher poloidal resolution at a reasonable computational cost , two - dimensional simulations have the additional advantage that field - lines can be easily visualized as level surfaces of the vector potential component @xmath32 . using this approach , we illustrate the magnetic field topologies for the three different poloidal configurations in figure [ topologies ] . we begin by examining the average properties of the accretion flows in simulations kdpg , qdpa and tdpa . @xcite described the accretion flow structure in terms of three distinct regions : the main disk body , the coronal envelope , and the unbound outflow . in this section , we study the first two , paying special attention to contrasts ( or lack of contrast ) when the magnetic topology changes . we will discuss the unbound outflow in the following section . much of the analysis will be in terms of various time- and space - averaged quantities . the radial profile of the time - averaged shell integral of @xmath68 is defined as @xmath69 similarly , the shell average of a quantity @xmath68 at a radius @xmath70 is defined to be : @xmath71 the angular profile of the time - average of @xmath68 at radius @xmath70 is defined by @xmath72 lastly , the time- and volume - average of @xmath68 is : @xmath73 in these equations @xmath74 is the time over which the integral is computed , and @xmath75 is the usual metric determinant . typically @xmath76 ; for kdpg and qdpg this is the last @xmath77 of the full evolution , while for tdpa we choose a @xmath77 window in the middle of the simulation after the accretion flow is established . the spatial extent of the shell integration is the full @xmath61 and @xmath63 computational domain . during a given simulation , various shell integrals and radial fluxes are computed and stored every @xmath78 in time . this data can then be integrated over time to obtain quantities such as the total or time - averaged jet outflow or accretion rate . it is also useful to divide the shell and volume integrals into two parts , one for bound and one for unbound flow . for simplicity we define a particular zone to be  unbound " if @xmath79 . unbound outflow can further be defined as those unbound cells with @xmath80 . in these simulations only the outflow near the axis ( the jet outflow ) is unbound ; the coronal backflow from the disk itself remains bound . the initial evolution of the accretion disks in kdpg and qdpa is qualitatively similar . both field configurations begin with considerable radial field within the torus . this is sheared out , generating toroidal field which , by @xmath81 , is sufficiently strong that the resulting poloidal gradient in @xmath82 begins to drive the inner edge of the torus ( initially located at @xmath83 ) inward . the inner edge of the disk arrives at the black hole at @xmath84 . within the disk body , the mri generates the turbulence that will determine the subsequent evolution of the disk , and by @xmath85 a statistically stationary turbulent accretion flow has been established inside the radius of the inner edge of the initial torus . the toroidal field model tdpa evolves more slowly than the two poloidal field cases , consistent with the results from earlier toroidal field pseudo - newtonian simulations @xcite . as discussed at length in @xcite , this behavior stems both from the absence of an initial radial field ( which means that there is no toroidal field amplification due to shear ) and from the fact that long - wavelength modes  which are the most effective in driving accretion  grow relatively slowly . inflow can begin only when the mri has produced turbulence of sufficient amplitude , which occurs by @xmath86 , corresponding to about 5 orbits at the radius of the torus pressure maximum . the accretion rate into the hole increases until about @xmath87 , after which it shows large fluctuations without an overall trend . figure  [ bndavgplt1 ] shows time - averaged , shell - integrated radial profiles of a number of quantities relevant to the accretion flow : accretion rate @xmath88 , surface density @xmath89 , the net accreted angular momentum per unit rest mass , @xmath90 , @xmath91 , the electromagnetic poynting flux @xmath92 and the em angular momentum flux @xmath93 . the subscripts fl and em denote the fluid and electromagnetic contributions to the stress - energy tensor , respectively . all quantities are computed in the coordinate frame after the turbulent accretion flow is established , and in all cases the volume integral was restricted to cells where the matter was bound . the poloidal field simulations were averaged over time @xmath94@xmath95 while the toroidal field data were averaged over @xmath96@xmath97 . these six plots may be divided into two groups : those with little dependence on initial field topology ( surface density and accreted angular momentum per unit mass ) , and those with a stronger dependence ( accretion rate , magnetic field strength , poynting flux , and electromagnetic angular momentum flux ) . although the radial distribution of the surface density in all three cases is fairly similar , there is a slight distinction between the two poloidal models on the one hand and the toroidal on the other . the surface density in the toroidal model rises somewhat more steeply with radius , with lower density near the black hole and higher values outside of @xmath98 . for all three topologies , the specific angular momentum carried into the hole is close to the angular momentum of the marginally stable orbit . for tdpa , the value is essentially equal to that of the isco , while both qdpa and kdpa have smaller values , indicative of continued stress operating near or inside the isco .    in all three cases , the time - averaged accretion rate is constant with radius out to @xmath9910 ; in this sense , all three have achieved a quasi - steady state within @xmath100 , where @xmath101 is the radial coordinate of the innermost stable circular orbit , @xmath102 ( the rise in @xmath103 and @xmath104 at large radius in tdpa is likely an artifact of the stronger magnetic field with which that simulation was begun ; equilibration takes place on timescales proportional to the local orbital period , which is , of course , @xmath105 , so inflow equilibration is slower at large radius ) . however , the actual value of the accretion rate is significantly different in the three simulations . the time - averaged mass accretion rate is highest in qdpa , exceeding that of kdpg by about @xmath106 , while @xmath107 is lower in tdpa by more than a factor of 2 compared to qdpa . accretion is , of course , driven by angular momentum transport , which is largely due to magnetic torques . to support these torques , the shell - integrated electromagnetic angular momentum flux must increase outward , as the net torque on the matter is the divergence of the em angular momentum flux . the dependence of the electromagnetic angular momentum flux on initial field topology is therefore similar to that of the mass accretion rate : it is consistently greater for kdpg and qdpa than for tdpa . the most significant difference between the kdpg and qdpa curves is that in the region near the isco , the electromagnetic angular momentum flux in kdpg becomes almost constant with decreasing radius ; this effect has previously been seen in other simulations that begin with dipole magnetic field loops @xcite . in the other two cases , the em angular momentum flux continues to decline inward . the mean magnetic field strength and poynting flux show a very similar dependence on initial field topology : qdpa and kdpg closely resemble each other , tdpa is somewhat weaker . in the toroidal case , the disk body average magnetic energy density is weaker by a factor of several than in the other two simulations , and decreases inward toward the horizon rather more rapidly . correspondingly , the time- and volume - averaged @xmath50 parameter was 11 in qdpa , 16 in kdpg , and 30 in tdpa ( because this is a volume - average , the corona is weighted heavily in this statistic ) . the electromagnetic energy flux shows the greatest variation between models . @xmath108 in kdpg is consistently larger than the value in qdpa : by a factor of two at @xmath109 and a factor of six at @xmath110 . the value in tdpa is smaller still by a factor of @xmath111 at these same radii . all these relations can be understood in a straightforward fashion : there is a nearly - constant ratio between the em angular momentum flux ( and therefore the accretion rate ) , the poynting flux , and the magnetic field intensity . in figure [ ratio ] , we plot the ratios of @xmath103 , @xmath108 and@xmath93 to @xmath91 . the inter - simulation contrast in the means of these quantities falls by a factor of several from the plot of their absolute values ( fig . [ bndavgplt1 ] ) to the plot of their ratios to @xmath82 ( fig .  [ ratio ] ) . thus , we see that their primary dependence is on the magnetic field strength . only in the ratio of the electromagnetic energy flux to magnetic field strength do we see any global variation between the field topologies : @xmath108 in kdpg is consistently a factor of two above qdpa and tdpa . this contrast arises from the strong poynting flux in kdpg associated with the slower - moving , bound portion of the funnel - wall outflow ( as can be seen from figure 10 of * ? ? ? * see  [ jets ] ) . the behavior of the fluid - frame stress can be seen in figure  [ maxwell ] , which depicts its time - averaged vertically - integrated value as a function of radius . to compute this quantity , we project the stress tensor onto a local tetrad system in each cell and multiply by the fluid - frame cell volume , which is computed by a similar projection @xcite . after integrating over the bound matter on that spherical shell , we normalize the result to the surface area that this finite - thickness shell occupies in the equatorial plane . we present the result in this form in order to compare with the conventional representation of the vertically - integrated stress given by @xcite , which is presented in these terms . the novikov - thorne fluid - frame stress , unlike ours , is _ assumed _ to be zero at and within the isco . in the figure , the novikov - thorne stress is normalized by the time - averaged value of the accretion rate in the simulation .    outside the isco , the radial profile of the maxwell stress is similar to the prediction of the novikov - thorne model , but in no case is there a roll - over in the value near the isco . as has been noted previously for the initial dipole field topology @xcite , the @xmath112 stress does not go to zero at the isco , but is continuous down to the event horizon . this is true for all three field topologies , at least in a time - averaged sense , although the amplitude of the stress declines as one goes from dipole to quadrupole to toroidal field topology . in the dipole case , the logarithmic gradient of the stress as a function of radius actually steepens as the accretion flow approaches the horizon .    as can be seen from figure  [ bndavgplt1 ] , this extra stress does reduce the specific angular momentum carried into the hole by the accretion flow . at the isco the accretion flow in the toroidal field simulation has @xmath113 ( where @xmath114 is the specific angular momentum of the isco ) and the accretion flows in the dipole and quadrupole simulations have @xmath115 ( a change of @xmath116 compared to @xmath117 ) . below the isco , @xmath104 in the quadrupole and toroidal cases remains constant with decreasing @xmath70 . in the dipole simulation , @xmath104 decreases to @xmath118 at @xmath119 , a decrease of @xmath120 compared to @xmath117 . the fall off in the dipole case below the isco is presumably due to the steepening of the maxwell stress in this region .    as has been remarked in earlier papers @xcite , @xcite , when @xmath121 , the magnitude of the maxwell stress ( i.e. , the electromagnetic angular momentum flux ) is smaller than the angular momentum flux associated directly with matter ( @xmath122 ) in the inner disk . it follows ( as shown in figure  [ ratio ] ) , that the mean accreted angular momentum per unit rest - mass is only slightly below the novikov - thorne value ( i.e. , @xmath123 at the isco ) . nonetheless , this continuity in the stress is important for other reasons : because the dissipation is rather more centrally - concentrated than the stress , this additional stress can be significant in the disk heat budget ; and the relative importance of the additional stress increases sharply for higher black hole spins . thus , we see that the initial magnetic field geometry makes relatively little qualitative difference to the character of the accretion flow . the sense of its ( weak ) influence is that the simulation beginning with a purely toroidal magnetic field results in a saturated field amplitude a factor of two smaller than those resulting from poloidal seed fields . the poynting flux , electromagnetic angular momentum flux , and accretion rate scale in direct proportion to the magnetic field energy density .    in their study based on axisymmetric simulations , @xcite also found that certain disk properties are largely independent of the initial field configuration . they found that the shell - averages of magnetic field quantities within the bound portion of the accretion flow ( strictly speaking , where @xmath124 ) follow power - law distributions that are the same regardless of the initial field topology . for example , they found that the field amplitude @xmath125 from the isco out to about @xmath126 . integrating the same magnetic field quantities ( magnetic field strength and the lab - frame projections of the ct - fields ) over the bound portion of the flow in our simulations reveals that these quantities do scale roughly as power - laws in radius , but with shallower gradients than as determined by @xcite . in addition , we see a slow trend with field topology , as well as significant fluctuations in time : @xmath127 in simulation kdpg ( where @xmath128 denotes 1 std . deviation about the time - average ) , @xmath129 in simulation qdpa , and @xmath130 in simulation tdpa . the contrast between our results and @xcite could result from a number of considerations : the initial torus configurations differed in their radial pressure and angular momentum profiles , and our simulations are fully three - dimensional , rather than axisymmetric . these scalings may or may not be universal , but they appear to have little to do with the jet . as we will shortly show in  [ jets ] , jet properties differ strongly from one initial magnetic topology to another ; it must be , then , that jets depend on other properties in addition to the current distribution in the disk . @xcite noted that the main disk body is surrounded by a low density region where the magnetic and thermal energies are in approximate equipartition . this region was designated as the coronal envelope . in contrast to the main disk body , where poloidal fluid velocities are dominated by turbulence , the poloidal motions in the coronal envelope are smooth on small scales and create an outgoing , but bound , backflow .    figure [ beta ] displays the spatial distribution of the time- and azimuthally - averaged @xmath50 parameter overlaid with contours of the gas density for each of the three different topologies . kdpg shows the familiar strongly - magnetized axial funnel and a corona with relatively large regions of @xmath131 . qdpa shows a lower level of magnetization in the funnel region , while tdpa has much larger regions with @xmath132 overall . to quantify these impressions , we adopt a more precise definition of the corona than that given in @xcite . we define it as those regions of _ bound _ material that lie outside of @xmath133 density - scale heights from the mid - plane , that is those cells that satisfy both @xmath134 and @xmath135 , where @xmath136 is the equatorial density and we measure along curves of constant @xmath70 . using this definition , we integrate over the corona and over time to compute the average total volume and mass fraction in the corona , the volume - averaged @xmath50 parameter , and the mean fluid - frame thermal and magnetic energy per unit mass . these data are given in table  [ corona ] . the coronal regions for all three simulations take up half or more of the total simulation volume but contain only @xmath137 of the mass . the volume assigned to the corona for kdpg is a smaller percentage compared to the other models , but kdpg also contains an unbound jet along the funnel axis , which is excluded from the corona by definition . within the coronal volume , the average thermal and magnetic energies per mass are similar , with volume averaged @xmath50 values of @xmath1384 . figure [ angcor ] plots the angular profiles of density @xmath139 , gas pressure @xmath140 , magnetic pressure @xmath141 and @xmath50 @xmath142 at @xmath126 . all three magnetic topology cases have very similar density profiles in the disk body , but appear to differ in the corona . although our averaging period of @xmath77 is 30 orbital periods at the radius shown ( @xmath143 ) , it is not long enough to erase completely fluctuations in the coronal density profile : observe that the quadrupolar case is the most extended in the  southern hemisphere " corona while the toroidal case is the most extended in the  northern hemisphere " . consistent with the general behavior we have already emphasized , the gas and magnetic pressure profiles in the disk body are very similar in the dipolar and quadrupolar cases , but both pressures are rather smaller in the toroidal field case . nonetheless , all three cases share one important property : the magnetic pressure is almost flat for several scale - heights around the midplane . this finding is consistent with results from shearing - box simulations with much better vertical resolution @xcite . in the corona , the magnetic pressure drops rapidly with height in qdpa and tdpa , but remains nearly constant in kdpg . in all three cases , the magnetic and gas pressures track each other closely , so that the @xmath50 profiles in the three runs are very similar . figure [ angtrptrt ] plots the angular profiles of the electromagnetic contributions to @xmath144 and @xmath145 ( @xmath146 and @xmath147 respectively ) at two different radii , @xmath110 and @xmath109 . at the larger radius , it should come as no surprise that the profile of magnetic pressure shown at that radius in figure  [ angcor ] accurately predicts the profile of these two electromagnetic quantities . at the smaller radius , as we have already shown ( fig . [ bndavgplt1 ] ) , the magnetic field is considerably more intense in the dipole case than in the other two , so both the electromagnetic angular momentum flux and energy flux follow suit .    to summarize this section , we have found that , just as for the main disk body , the character of disk coronae is only weakly dependent upon the initial field topology . there are modest quantitative contrasts of the same sort as seen in the disk body  dipolar and quadrupolar initial fields both create somewhat stronger magnetic fields on average than does initial toroidal field  but that is all . the formation of unbound outflows , or jets , is one of the most striking features of the kd models @xcite , all of which were initialized with single - loop dipolar magnetic fields . these outflows generically have two components : a poynting flux jet inside the axial funnel formed by the centrifugal barrier , and a sheath of much denser unbound gas moving outward at sub - relativistic speed along the funnel s outer edge . as described in @xcite , energy for the poynting flux comes from the rotation of the black hole , in a manner closely related to the classical blandford - znajek mechanism . in the case of the schwarzschild hole , the absence of black hole rotation means that the magnetic field in the funnel , although relatively strong , is mostly radial and carries no poynting flux . the funnel wall jet , on the other hand , appears to be powered by the significant gas and magnetic pressure near the isco . the gas retains enough angular momentum to be excluded from the funnel and is squeezed outward along the centrifugal barrier . in contrast , the three - dimensional quadrupole and toroidal field simulations have substantially weaker unbound outflows , if any . figure [ ubndavgplt1 ] shows the time - averaged shell - integrals of the mass outflow rate @xmath148 , magnetic field strength@xmath91 ; along with the electromagnetic energy flux , @xmath92 and angular momentum flux @xmath149 in unbound material between @xmath150 m and @xmath151 m . the values in kdpg are 10100 times larger than those in qdpa , whose outflow quantities are , in turn , an order of magnitude larger than those in tdpa . for example , the unbound mass outflow rate at @xmath152 in kdpg is @xmath153 times that of qdpa , while the mass outflow rate in tdpa is 10 times smaller than in qdpa . the strength of the outflow is directly related to the strength of the funnel field . in kdpg the time - averaged funnel magnetic energy density is 20@xmath154 times greater than in the quadrupole model , while the toroidal model is weaker by another factor of 10 . the mean poynting flux in kdpg is therefore @xmath155 times the poynting flux for qdpa . in other words , unlike the qualitative similarity of the accretion flows produced by these different field topologies , the outflows differ dramatically . to understand better the origin of this sharp contrast , in the next section we focus more closely on exactly how field is introduced into the funnel , and how its magnitude is controlled . we have seen that the strength of the poynting flux jet is very sensitive to the initial magnetic field topology . because the establishment of large - scale field within the funnel appears to be critical to driving a jet with substantial poynting flux , we devote the remainder of this section to a close examination of just how different magnetic topologies lead to different funnel magnetic fields , leaning heavily on high - resolution 2-d axisymmetric simulations . we begin by analyzing the axisymmetric simulation whose initial condition has already been shown to produce a significant jet , the dipole . figure  [ dfpevln ] shows a time series of color plots of the @xmath50 parameter overlaid with white contours depicting poloidal field lines . as the radial field is sheared , a toroidal field ( visible in fig . [ dfpevln ] through its effect on @xmath50 ) is created and amplified , leading to poloidal magnetic pressure gradients @xcite . because the initial radial field has opposite sign on opposite sides of the equatorial plane , the toroidal field created by the shear also changes sign across the equator . the inner edge of the torus begins to move inward toward the event horizon , dragging the magnetic field with it . as the gas moves inward , its rotation velocity increases , further strengthening the toroidal field . as the magnetic pressure increases , gas is forced toward the equator where the field is comparatively weak , bringing the regions of oppositely - directed toroidal field closer together ( see , e.g. , the panel illustrating @xmath156 ) . as a result , a strong current sheet forms along the equatorial plane . once the inner edge of the accretion flow arrives at the black hole ( @xmath157 ) , field lines attach to the event horizon , and gas rapidly drains off them . field lines slide toward the poles as they expand outward at relativistic speed , driven by the contrast between the high magnetic pressure near the equatorial plane and the much lower pressure higher up along the axis . this flow moves outward so swiftly because it is almost entirely unencumbered by the inertia of matter . as discussed in @xcite , a magnetic tower forms , similar to the one predicted on analytic grounds @xcite and seen in other simulations @xcite however , in contrast to these models , tightly - wrapped field lines persist only briefly ; the large vertical gradient in magnetic pressure propels the upper loops so rapidly outward that in the steady - state the field lines are better described as loosely helical than as a tightly - wrapped coil . the radial field inside the funnel has the same sign as the vertical component of the initial dipole loop inside the initial pressure maximum ; these field lines close outside of @xmath151 m , reentering the problem volume throughout the corona . a time history of the unbound poynting flux at @xmath152 , normalized by the instantaneous ( coordinate frame ) magnetic field energy within the bound accretion flow , is shown in figure  [ 2dtrtdu ] for both the three - dimensional simulation kdpg and the high - resolution axisymmetric model . in both cases , once the jet is established ( a process taking @xmath158 ) , the poynting flux jet persists . it is frequently perturbed , and the amplitude undergoes fluctuations , but its basic properties persist throughout the simulation . excluding initial transients ( i.e. @xmath159 ) , the jet histories of the two simulations are very similar to one another , despite the differences in dimensionality , ( poloidal ) resolution , torus configuration and equation of state : in both cases , the ratio of jet em luminosity to energy density in magnetic energy is @xmath1600.2 per time - unit . this similarity demonstrates that , for this topology at least , the properties of the funnel field depend primarily on the field strength within the initial accretion flow , and that funnel field creation is a primarily axisymmetric process .    although the initial magnetic field had structure only on the scale of the disk thickness , the _ essentially laminar _ dynamics of initial inflow lead quickly to the formation of truly large - scale field . on the other hand , because the accretion flow takes longer to be established than the jet does , the magnetic field in the funnel acts in effect as a permanent boundary condition for the accretion flow . in other words , these simulations mimic the long - term behavior of an accretion flow in which there is large - scale net magnetic flux threading the event horizon and filling a cone around the axis , but no large - scale flux piercing the disk proper . as both the three - dimensional and two - dimensional simulations demonstrate , it is entirely possible for long - term quasi - stationary accretion to coexist with a strong poynting flux jet whose foundation is magnetic flux trapped long before . we next consider the evolution of an initial quadrupolar field in high resolution axisymmetry . figure [ qfpevln ] shows a time series of plots illustrating @xmath50 and the poloidal field line structure for this case . initially , the loops of weak poloidal magnetic field enclose two regions , located above and below the pressure maximum ( see figure [ topologies ] , center panel ) . in a similar fashion to the dipole case , these loops are sheared out , creating ever growing toroidal field and producing strong poloidal magnetic pressure gradients . in contrast to the dipole case , however , two current sheets are formed , one above the equatorial plane , one below . the magnetic pressure forces gas toward these current sheets , where the magnetic pressure has a local minimum . as the inner edge of the torus moves radially inward toward the black hole , the current sheets maintain their structure . when the field reaches the black hole , the current sheets are located above and below the equator . as with the dipole , gas drains off the field lines , which anchor themselves in the event horizon , and the field lines between the current sheets and the axis slide toward the poles , where the density and pressure are very low . these field lines expand outward relativistically , filling the funnel region and establishing a large - scale quadrupole field . while this field is in place it supports a poynting - flux jet , but the power in the jet is significantly less than that in the axisymmetric dipolar simulation . averaged from @xmath161 to @xmath162 , the em flux in the dipolar case is twice as strong as in the quadrupolar simulation when measured at @xmath126 , and 10 times stronger measured at @xmath152 ( we will discuss later why the contrast was more than an order of magnitude greater in the 3-d simulations ) . figure  [ trtducomp ] compares the dipole jet poynting flux with that of the quadrupole simulation as a function of time at @xmath152 , normalized by the total magnetic field energy in the disk . in terms of this normalized power , the jet in the quadrupole simulation is a factor of 210 weaker than that associated with dipole model ; in other words , much of the actual power contrast between the dipolar and quadrupolar cases can be explained in terms of a reduction in this  magnetic efficiency " . unlike the situation for accretion properties , contrasts in jet luminosity from one topology to another can not be described solely in terms of contrasts in mean magnetic field strength . @xcite reported that replacing the dipole with a quadrupole field led to a jet that is weaker by a factor of 23 , but for several reasons one should expect only qualitative , not quantitative agreement with our results . they measured the jet luminosity at @xmath163 and their initial pressure maximum was much closer to the black hole ( @xmath164 ) , so that the radial dynamic range of their genuine accretion flow was considerably smaller . their disk was also geometrically thicker , partly because their code solved a total energy equation , and partly because their initial state was hotter . simple geometry accounts for the qualitative contrast between the jet produced by an initially quadrupolar magnetic geometry and that driven by an initially dipolar field . in the quadrupolar case , relatively small amplitude vertical oscillations in the inner region of the accretion flow can lead to reconnection between magnetic field in the funnel and magnetic field in the disk because vertical field of both signs is available in the same hemisphere . the result is to break the magnetic connection between the black hole horizon and large radius . put another way , the presence of closed field loops on the same side of the equator offers the opportunity for easy collapse of these loops to zero . in the dipole case , however , the sign of the vertical magnetic field passing through the equator is the same everywhere inside @xmath165 . the net result of this geometric contrast is that , while quadrupolar field configurations can pump magnetic field into the outflow funnel and support a jet , they can also equally easily destroy such a jet by reconnection with oppositely - directed field . jets in this case are consequently very unsteady , and exhibit strong variability in the poynting flux , as shown in figure  [ trtducomp ] . this variability is a major factor in depressing the time - averaged em luminosity of jets driven by quadrupolar fields .    because the code simulates the equations of ideal mhd , reconnection is numerical and occurs at the grid scale . as such , the reconnection rate is resolution dependent ( lower resolution promotes reconnection ) . we can not therefore make specific predictions regarding the variability characteristics of jets in nature , as this is determined ( in our simulations at least ) by competition between processes that are physical in character ( e.g. , orbital dynamics , magnetic forces ) and purely numerical ( the reconnection rate ) .    the rate of reconnection is also determined by details of the field geometry , in addition to its basic topology . to investigate this sort of dependence , we have conducted a series of lower resolution ( @xmath166 ) axisymmetric quadrupolar simulations in which we varied the size of the initial torus , and hence the scale of the initial field loops . we kept the inner edge of the torus fixed at @xmath83 , but reduced the angular momentum @xmath27 in order to reduce the vertical thickness of the torus , and consequently reduce the vertical extent of the quadrupolar loops . simulations were carried out for a variety of values of @xmath27 , ranging from @xmath167 to @xmath168 in steps of @xmath169 .    for the purposes of our discussion here , we focus on the @xmath170 torus , which serves to illustrate the impact of vertical scale height on the evolution of the quadrupole topology . specifically , setting @xmath171 results in a torus which is initially a factor of two thinner ( in terms of the density scale height , measured along surfaces of constant radius at the location of the pressure maximum ) than the @xmath25 torus previously discussed . the vertical thickness of the resulting quadrupolar field loops is then very similar to that of loops in the sub - keplerian torus used in qdpa . the initial evolution of this system is shown in figure [ qnpevln ] . the evolution proceeds as before : current sheets form above and below the equatorial plane as the inner edge of the torus moves toward the black hole . this time , however , the two current sheets are much closer to the equatorial plane ( contrast fig . [ qnpevln ] with fig . [ qfpevln ] ) . around @xmath172 the current sheets buckle , apparently due to an instability , and the coherent poloidal field structure is destroyed . a similar phenomenon was observed in qdpa , which had a similar scale - height : the current sheets destabilized far from the black hole , destroying the coherent poloidal field structure . as a result , no large scale poloidal field was dragged into the near horizon region by the torus , and the resulting field inside the funnel region is substantially reduced . we note that in three - dimensions the presence of non - axisymmetric mri modes are likely to enhance the rate of reconnection , further decreasing the likelihood that an organized quadrupole configuration will survive to create a significant funnel field . these additional simulations show that the strength of poynting - dominated jets depends on field geometry in two senses : it is topology - dependent and shape - dependent . quadrupolar magnetic topology leads to easier reconnection , and therefore weaker funnel fields and jets , than does dipolar topology , and geometrically thinner quadrupole loops make reconnection still more rapid . we have shown how large dipolar loops can create strong jets , but what happens if the oppositely - directed vertical field on the far side of the loop is accreted ? this did not occur in kdpg , in part because it did not run long enough and in part because the finite mass in our torus meant that the matter outside the pressure maximum moves outward over time , not inward . in a real accretion disk , however , continuing mass re - supply will force the far sides of even large dipolar loops to accrete .    to explore this sort of event , we ran a two - dimensional simulation whose initial state contained a series of four narrow dipolar loops ( described in  [ numerics ] ) . the initial evolution of this topology is shown in figures [ mfpevln1 ] and [ mfpevln2 ] . the first of these figures shows the descent of the innermost radial loop into the near - horizon region and the establishment of a poynting flux jet . the second figure shows the descent of the second ( moving out in radius ) of these loops and its interaction with the field structures established by the first loop .    over the timescale considered , the evolution of this model is determined by the two innermost loops , both initially inside the pressure maximum ( see fig . [ topologies ] ) . as for a single dipole loop , the first loop forms a funnel dipole field by @xmath173 . however , its shape differs significantly from the single loop because the second loop begins its descent toward the inner disk even before the inner edge of the first loop nears the event horizon ( fig . [ mfpevln1 ] ) . as the second loop moves inward , it displaces the field lines of the first loop , forcing the midplane portions of the outer half of the first loop inward , and the off - plane portions of the first loop to higher altitude at their original radius . the result is that the first loop is very nearly transformed into a pair of loops , one on each side of the midplane . after arriving in the inner disk , the field of the second loop expands rapidly toward the funnel region . it pushes on the almost closed sections of the first loop , triggering rapid reconnection where its oppositely - directed segments pass very close to one another at small radius ( fig . [ mfpevln2 ] ) . at this point , the transformation of the first loop into a disconnected quadrupolar loop pair is complete , and the loops fly off to large distance , emptying the funnel of magnetic field . the funnel is then refilled with flux from the second loop . these dramatic events are mirrored in the very large fluctuations in normalized poynting luminosity shown in figure  [ trtducomp ] . jumps in output even greater in relative terms than those displayed by the quadrupolar simulation are frequent , and can be associated with these episodes of funnel field establishment and annihilation . note that while the multiple - loop jet is highly variable , similar to the quadrupolar jet , the instantaneous strength of the poynting - flux jet ( when it is strong ) is similar to that of model kdpg when normalized for disk magnetic field strength . in this sense , the  efficiency " of these systems in using magnetic field generated in the disk for jets appears to be a function of poloidal field character , declining as the relative importance of higher multipole moments increases .    in interpreting this simulation , it is important to bear in mind that , unlike any of the others reported in this paper , its accretion flow never achieves anything resembling a quasi - stationary state . the multiple loop structure also requires high resolution to prevent reconnection from occurring at too rapid a rate . it is for this reason that , unlike the large dipolar loop and quadrupolar caess , we report no 3-d simulation results for this field geometry . thus , we can be less confident than in the other cases that the specific phenomenology we observe may be associated with real systems in a state of continuing accretion . what it _ does _ demonstrate , however , is that accretion of a series of closed dipolar loops can lead to a sequence of field establishment and destruction rather similar to what was seen in the quadrupolar case , and through mechanisms that also resemble those seen when the initial field is quadrupolar . it also supports our qualitative inference that maintenance of a consistent sign of vertical field is essential to the long - term support of a poynting - dominated jet . the lifetime for an episode of strong outflow is roughly the inflow time for the outer anchor point of the loop .      in the poloidal cases , magnetic field is injected into the funnel by strong toroidal field pressure in the plunging region . one might expect that this process could act equally efficiently when the magnetic field is toroidal from the start . if the toroidal field has a consistent sign , the radial component created when it expands outward would also have a consistent sign , satisfying the coherent vertical field criterion we have inferred from study of the poloidal simulations . if so , the toroidal field case should also drive a powerful jet . however , it does nothing of the sort . figures  [ tdpevln1 ]  [ tdpevln2 ] are a series of snapshots of the @xmath50 parameter and the radial field structure in the evolving torus of simulation tdpa . in contrast to the poloidal loop simulations , there is no initial radial field to generate toroidal field through shear . instead , the evolution is governed entirely by the development of the nonaxisymmetric mri . the fastest - growing modes for the mri in this context have large @xmath174 and azimuthal wavelength @xmath175 , which is much shorter than the disk thickness when the field is weak @xcite . the magnetic structures that develop are then predominantly on short lengthscales . note the rapid spatial variation in sign in the radial field in figure  [ tdpevln2 ] . with time , the radial field features become somewhat more radially elongated , but field lines do not connect extended regions ; no large - scale coherent poloidal field is generated from the initial toroidal field . magnetic field amplification is less than in the poloidal field initial conditions , and @xmath50 never becomes less than 1 near the hole . matter tends to remain with the field rather than drain either into the hole or toward a localized region such as the equator , as seen in the dipole case . any vertical expansion of the field would therefore have to overcome the inertia of the matter and , as can be seen from figure [ tdpevln2 ] , this does not occur . the absence of a funnel - field in the toroidal topology follows directly from the specific properties of the toroidal field mri : the poloidal fields created by the mri are typically small scale . no radially - extended poloidal field connects the near horizon region to the disk . what radial field there is is accreted faster than it can tap into the differential rotation and generate strong near - hole fields . the main disk - body in the toroidal case remains gas pressure - dominated at all times , and there are no strong magnetic pressure gradients to drive an expansion of the field into the funnel region . the absence of a funnel field means there is no poynting flux jet . it appears , then , that creation of a strong funnel field , and therefore a strong poynting jet when the black hole spins rapidly , depends essentially on the existence of some coherent poloidal field in the initial state . purely toroidal field initial conditions can lead to short coherence - length poloidal loops with rapidly alternating sign , but no large - scale poloidal field components . through a series of numerical studies , we have examined the influence of the initial field topology on the properties of magnetized black hole accretion flows . in particular , we have performed full three - dimensional simulations of three different initial field topologies : a large - scale dipolar loop , a pair of quadrupolar loops separated by the equatorial plane , and a purely toroidal field . apart from the imposed initial field topologies , the initial conditions used in each of these simulations were identical : a @xmath176 hydrostatic torus with a slightly sub - keplerian distribution of angular momentum around a rotating ( @xmath177 ) black hole . the two poloidal field topologies were evolved to @xmath178 m , while the toroidal case was evolved to @xmath179 m . @xcite described the late - time structure of the resulting accretion flow in terms of distinct components : the main disk body , the coronal envelope , and an unbound outflow which itself could be separated into a matter - dominated funnel - wall jet and a poynting - flux jet in the funnel interior . we have found that the late - time averaged properties of the main disk body exhibit only modest quantitative contrasts attributable to different initial field topologies : the dipolar and quadrupolar cases are roughly similar , and both create somewhat stronger magnetic fields than the toroidal case . so long as there is some weak magnetic field present , the mri generates turbulence and drives accretion , although the character of the purely toroidal mri is somewhat different from that produced by poloidal fields . because the maxwell stress controls the surface density required to support a given accretion rate , and the cascade of energy to smaller scales in the mhd turbulence leads ultimately to disk heating , we expect that the radiation emitted by the disk will likewise be only weakly dependent on magnetic topology . it is possible that different magnetic topologies could produce different vertical profiles of dissipation ( which these simulations are not designed to study ) , but the similarity in magnetic pressure profiles ( fig . [ angcor ] ) suggests that even those differences are likely to be minor .    in regard to the accretion flow , we found only two dependences on field topology worth noting : for fixed outer torus mass , the initial toroidal field simulation gives a somewhat lower accretion rate . in addition , the magnetic field in the plunging region is , in relative terms , somewhat stronger in the large dipolar loop case than for quadrupolar or toroidal fields , so that stresses near the isco are somewhat stronger when the field topology is dipolar . in other words , the short - wavelength behavior of the mri - driven mhd turbulence in the disk produces near - universal behavior , independent of magnetic field structure . jets are quite different . _ laminar _ flow dynamics can , given the right conditions ( e.g. , the right field topology ) , inflate field with no pre - existing structure on scales larger than a disk scale - height into a truly global field . simulations that begin with a large dipole loop in the initial torus produce a strong poynting - flux jet in the axial funnel and a comparably strong matter - dominated funnel - wall outflow . the inner field lines of the dipolar loop move quickly inward to fill the funnel , and remain there until the opposite side of the loop is accreted , an event which may not happen until quite a long time later . thus , this case mimics a situation in which net magnetic flux has been trapped against the black hole event horizon by previous accretion . when this trapped flux creates a coherent poloidal field in the funnel , the black hole s spin forces nearby field lines into rotation and drives a large outward - directed poynting flux . this sequence of events is not as easily achieved with any of the other field geometries we explored . unlike the dipole configuration , the quadrupole features a pair of current sheets located above and below the equator . these current sheets are subject to instabilities that can lead to reconnection in the magnetic field and loss of extended field coherence . even when the quadrupole is successful in establishing a funnel field , the jet may live only briefly because fluctuations in the accretion flow can generate reconnection events . by the same token , quadrupolar field can also restore field to an empty funnel . consequently , any jet based on quadrupolar field tends to be unstable and episodic , with the timescale and amplitude of changes set both by details of the field geometry and the rate of reconnection . unfortunately , it is very hard to obtain any information about the former in real systems , and in our simulations the latter is controlled by numerical effects and resolution , so it is impossible for us to make predictive statements about the timescale or amplitude of jet variability in this case . all of the models reported on above used a black hole with a spin @xmath180 . to test our conclusions against other black hole spins we also ran two three - dimensional initially quadrupolar simulations with black hole spins @xmath181 and , @xmath182 . the effect of a wide range of black hole spins ( including a counter - rotating hole ) on dipolar models was studied by @xcite . no poynting flux jet is present for a schwarzschild hole regardless of the initial field topology . the energy to drive the poynting flux originates with the spin of the hole . however , the funnel is still filled with radial field , which is considerably stronger in the dipole case compared to the quadrupole . on the other hand , the funnel field strength increases with black hole spin for a given field topology . for the rapidly - spinning hole , the poynting flux jet is stronger than in the @xmath121 simulations for both topologies , but as we have discussed , it is considerably weaker when the initial field is quadrupolar than when it is dipolar . specifically , the ratio between the time - averaged poynting luminosity with dipolar field to the poynting luminosity with quadrupolar field is @xmath183 for @xmath184 and @xmath185 for @xmath121 . thus , the accretion dynamics control what sort of field is brought to the vicinity of the black hole and fill the funnel , and the black hole spin determines how energetic the resulting jet is . the purely toroidal field topology does not generate any jet at all . the only way for toroidal field to create any poloidal field is by initiating non - axisymmetric perturbations of alternating sign . the lengthscales of these perturbations must be relatively short , and the portions of opposite sign must balance in magnitude . fields of this sort can do very little to fill the funnel region because they are even more subject to reconnection than is the quadrupolar field .    when the initial field is a single dipolar loop that wraps around the initial pressure maximum of the simulation , all the vertical field brought down to the funnel by the initial accretion has the same direction . in real disks containing closed dipolar loops , oppositely directed vertical field will , sooner or later , arrive in the inner disk , and it may cancel the funnel field . to investigate this possibility , we ran a high - resolution axisymmetric simulation that began with four narrow dipolar loops within the initial torus . its subsequent evolution qualitatively supports this conjecture . the innermost loop produces a funnel field as in the single dipole loop model . when the second loop , oriented oppositely from the first one , descends , it displaces the field established by the first loop . the funnel field established by the first loop reconnects , detaches from the event horizon , and is rapidly expelled from the grid . the second loop then expands poloidally and reestablishes the funnel field . in other words , the accretion of dipolar field loops leads to jet creation , but its structure is highly - variable : the narrower the dipolar loops present in the accretion flow , the more variability we can expect in the jet .    when a simulation is initialized with ( possibly multiple ) dipolar field loop(s ) that cross the equatorial plane , the strength of the resulting poynting - flux jet depends on the strength of the magnetic field within the disk itself and on the spin of the hole . poynting - flux jets that are formed from quadrupolar fields are generically much weaker than their dipolar counterparts , even when magnetic field strength in the disk body is taken into account . how much weaker depends on other factors as well , including disk thickness , location of the inner edge of the initial torus , and numerical resolution . in the limited range of numerical experiments reported here , we have seen suppression factors ranging from several to several hundred .    real accretion flows are likely to have a mixture , possibly time - varying , of all these field topologies . real jets from black holes will exhibit behavior combining the sorts of events we have seen in this paper : stability over longer timescales when the field is mostly made of large dipolar loops , greater variability and occasional jet quenching when the field has a larger contribution from higher multipoles , jet elimination when the field is mostly toroidal . it is possible that in addition to the sort of coherent funnel field spontaneously created by large closed - loop poloidal structures , there is a large - scale field imposed by external boundary conditions . we plan to investigate the impact of such an externally - imposed field in future work . this study is only one small step in the process of exploring the mechanisms to create large scale field configurations within accretion flows . here the large scale field formation process we have observed results from an interaction between poloidal field components in the accretion flow and the boundary represented by the black hole . the mechanism operates in axisymmetry as well as in three dimensions , so it is not a dynamo if that word is defined as in the antidynamo theorem . many authors have conjectured that an inverse cascade within the three dimensional turbulent disk itself might generate large scale fields ; we see no evidence in these simulations of any disk field with structure on scales more than a few scale - heights , but we caution that the resolution limitations of three dimensional global simulations make this observation less than definitive . the results we have reported also stimulate speculations about whether the contrasts in jet phenomenology from one object , or class of objects , to another is due in part to differences in the mixture of different field topologies , or the interaction of different topologies in the flow . for example , could the jet outbursts sometimes seen in galactic black hole binaries when they cross the  jet line " @xcite be due to rapid field annihilation and loop expulsion of the sort observed in the multiple dipolar loop simulation ? perhaps the radio - weakness of most agn ( see e.g. * ? ? ? * ) is due to the predominance of small - scale multipolar field in their accretion flows as a result of accreting from a turbulent interstellar medium ? as we learn more about how it acts , we may find that magnetic topology can , in fact , be inferred from phenomenology . surface normalized to the total magnetic field strength within the disk for simulation kdpg ( black lines ) and the high resolution axisymmetric simulation of the dipole ( magenta line).,scaledwidth=50.0% ]         surface normalized by the total magnetic field strength within the disk . the kdpg dipole simulation is shown in black , the quadrupole simulation in blue , and the multiple loop simulation is green.,scaledwidth=50.0% ]
Abstract:

Answer: because the magnetorotational instability is capable of exponentially amplifying weak preexisting magnetic fields , it might be hoped that the character of the magnetic field in accretion disks would be independent of the nature of the seed field . however , the divergence - free nature of magnetic fields in highly conducting fluids ensures that their large - scale topology is preserved , no matter how greatly the field intensity is changed . by performing global two - dimensional and three - dimensional general relativistic magnetohydrodynamic disk simulations with several different topologies for the initial magnetic field , we explore the degree to which the character of the flows around black holes depends on the initial topology . we find that while the qualitative properties of the accretion flow are nearly independent of field topology , jet - launching is very sensitive to it : a sense of vertical field consistent for at least an inner disk inflow time is essential to the support of strong jets .