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: the twirl operation was introduced by bennett _ @xcite in the context of entanglement purification of mixed states . the first step in many of these protocols is the twirl operation which converts an arbitrary mixed state of a two qubit system , @xmath0 , into a werner singlet state @xmath1 @xcite , an incoherent mixture of @xmath2 parts of a maximally entangled singlet state , @xmath3 , and @xmath4 parts of the maximally mixed state , @xmath5 the key idea behind the twirl is that the singlet state is invariant under any bilateral unitary transformation of the two qubits ( that is , any operation where identical local unitaries are applied to the two qubits ) , whereas any other state will be affected . a randomly chosen bilateral rotation should serve to average any other state , converting an arbitrary mixed state to a werner singlet . ( note that the twirl sequence does not _ create _ new singlet , rather it _ preserves _ the singlet , while averaging out all other terms . ) this provides a good `` in principle '' definition of the twirl operation , but in practice it is important to consider what a random bilateral rotation really is , and how it can actually be implemented . early work in this field concentrated on reducing the infinite set of randomly chosen rotations down to a small finite set . originally it was suggested @xcite that the twirl could be achieved using a set of four bilateral rotations , although it was later shown that a set of twelve bilateral rotations are both necessary and sufficient @xcite . this approach is sensible for conventional quantum information processors , but with ensemble devices , such as nuclear magnetic resonance ( nmr ) quantum computers , it is useful to use a different approach . in such devices it is in fact easier to perform rotations by angles that vary continuously over the ensemble than by angles from some small carefully chosen set . in nmr devices , for example , the application of a strong magnetic field gradient @xcite performs a rotation around the @xmath6 axis with a rotation angle that depends strongly on the position of a given molecule within the spatial ensemble ; integrating over the position within the sample is then equivalent to applying a random rotation . in the language of nmr this is often referred to as a _ crush gradient _ pulse . for simplicity we will sometimes describe a crush gradient as a random rotation around the @xmath6-axis , as the effect on the averaged density matrix is the same , but when considering experimental implementations it is necessary to be more precise . the use of field gradients forms the basis of spatial averaging methods , used in nmr quantum computing to prepare pseudo - pure states @xcite . field gradients have also been used to average out error terms @xcite , to project qubits into the zeeman basis @xcite , and to simplify density matrices prior to partial state tomography @xcite . the alternative procedure , in which one rotation chosen from a small set is applied , requires several experiments to be performed , and in the context of nmr quantum computing is usually known as temporal averaging .    in subsequent sections we explore different methods for implementing the twirl , concentrating on those which are best suited to ensemble techniques . we begin by noting that any set of bilateral rotations which performs a twirl must correspond to a set of single qubit rotations which , when considered from the viewpoint of a single qubit subsystem , averages any state of a single qubit to the maximally mixed state . as we shall see , however , averaging a single qubit is a necessary but not a sufficient condition for a set of rotations to act as a twirl . averaging a single qubit gives some useful insight into the problem of twirling a two qubit state . a single qubit state can be represented by a ray on the bloch sphere and an arbitrary rotation will move this ray over the surface of the sphere . each rotation @xmath7 is defined by a rotation angle @xmath8 about an axis @xmath9 parameterized by a tilt angle @xmath10 and an azimuthal angle @xmath11 . an arbitrary single qubit state is completely defined by three real parameters , which are conventionally taken as the expectation values of the single spin pauli operators @xmath12 @xcite . for scrambling an arbitrary state to the maximally mixed state , we must choose a convenient set of operations , which scrambles each of these pauli operators individually . a tempting , but incorrect , approach is to apply random rotations about random axes , that is to take @xmath8 as being uniformly distributed between @xmath13 and @xmath14 , while the rotation axes defined by @xmath10 and @xmath11 are uniformly distributed over the sphere ; we call this set of rotations @xmath15 .    brute force integration shows that the continuous set of random rotations around random axes , @xmath15 , does not completely randomize the state of a single qubit , and so _ can not _ constitute the basis for a twirl operation ; instead it reduces the bloch vector to one third of its original length . in retrospect the reason for this behaviour is clear : random rotations around axes perpendicular to the original state will completely average it , while rotations around axes parallel or anti - parallel to the state will leave it unaffected . thus the overall effect of @xmath15 must be to scale down the state , rather than to average it completely . the significance of the scaling by one third is discussed below . a better definition of a random rotation is provided by considering the description of a rotation by means of its euler angles . while many different conventions for describing euler angles exist , the essential feature is that any rotation can be decomposed as a sequence of three rotations around two axes . for example any rotation can be achieved by the sequence of rotations @xmath16 where rotations are applied from left to right . a random distribution of euler rotations can be achieved by taking @xmath8 as uniformly distributed between @xmath13 and @xmath14 , with the rotation axes defined by @xmath10 and @xmath11 uniformly distributed over the sphere ; we call this set of rotations @xmath17 . it might seem that @xmath17 is the same as @xmath15 , but they are in fact quite different : in particular @xmath17 completely averages a single qubit , while @xmath15 does not . this is easily seen by noting that the first rotation will average @xmath18 and @xmath19 , while the second and third will average @xmath20 . as any state of a single qubit can be written as a linear combination of these basic matrices and the maximally mixed state , the first two operations will average _ any _ state . it is clear from the above that a random euler rotation will average any single qubit , but this process is perhaps excessive , as two rotations suffice . ( we will see below that this is not true when considering a true twirl applied to two qubits . ) the sequence of rotations @xmath21 where both @xmath10 and @xmath11 are now taken as uniformly distributed between @xmath13 and @xmath14 , will average any single qubit ; as before this is most simply seen by considering the result for @xmath18 , @xmath19 and @xmath20 . indeed this process can be simplified still further : as described by bennett _ @xcite is is possible to average a single qubit by randomly selecting from the four operations @xmath22 . since @xmath23 is equivalent to a @xmath24 rotation , and ( neglecting global phases ) @xmath25 it is clear that this operation is equivalent to applying either @xmath13 or @xmath26 at random , and then applying @xmath13 or @xmath27 at random ; a similar result has been described by hayden et al @xcite . thus instead of choosing the two angles in equation [ rotzy ] at random from uniform distributions , we can choose at random from two particular values . this process of replacing a continuous rotation by a small number of discrete values is an example of a quite general procedure . for rank @xmath28 and rank @xmath29 tensors , continuous rotations about a fixed axis , can in fact , always be replaced by discrete rotations : we can `` hop '' between the angles @xmath30 , with @xmath31 and @xmath32 , and the overall effect is the same as rotating continuously about the same axis @xcite . it is _ sometimes _ possible to use @xmath33 , but this can not be relied on in general . the rotation through @xmath30 is the member @xmath34 of the cyclic group @xmath35 @xcite ; in this way , the continuous rotations can be replaced by discrete group operations and the continuous integrals over the rotation angle @xmath8 , can be substituted by _ group theoretical averages _ , such as @xmath36 the simplest discrete group that can generally be used for averaging is the cyclic @xmath37 group , which is employed by nmr spectroscopists in the form of magic angle hopping @xcite , a discrete version of its continuous counterpart , magic angle spinning @xcite . the group theoretical framework has also been used in the context of bang - bang decoupling protocols @xcite , which seek to average out system - environment interactions by applying fast , discrete and periodic control impulses to _ only _ the system , with the control propagators faithfully representing the members of a discrete cyclic group @xmath38 . this approach hints at another method for performing complete averaging of a single qubit . we have already noted that random rotations around a random axis do not completely average a single qubit , but simply rescale its length by one third . this process is , as usual , equivalent to applying a rotation through an angle of @xmath39 , @xmath40 or @xmath41 around a random axis . clearly a @xmath39 rotation can have no effect , while the effect of the @xmath40 and @xmath41 rotations must be the same . from this it can be deduced that a @xmath40 rotation around a random axis will completely average a single qubit , and this is indeed the case . finally we turn to the issue of practical experimental implementations . for a conventional quantum information processor the obvious approach is to use the smallest discrete set of operations ; for the case of averaging a single qubit this is the set of four rotations described by bennett _ et al . _ with an ensemble processor , such as an nmr quantum computer , it is better to use a procedure based on continuous rotations , corresponding to spatial averaging . clearly the best approach is to apply random rotations around two orthogonal axes , such as @xmath6 and @xmath42 . a strong magnetic field gradient , denoted * g * , will effect a random rotation around the @xmath6 axis , and rotations around other axes can be achieved by combining gradients with single qubit gates ( radiofrequency pulses ) which can be treated as rotating the axis system . thus the sequence @xmath43 will completely average a single spin . note that this procedure is not completely equivalent to equation [ rotzy ] , as it should strictly speaking be followed by a @xmath44 pulse to rotate the axes back to their original positions , but as the maximally mixed state will not be affected by the rotation this final stage may be safely omitted . the discussion above assumes that the gradient pulses are instantaneous , so that their sole effect is average the spin state and evolution under any background hamiltonian can be ignored . in fact , gradient pulses take a finite time , but as the gradient and background hamiltonians commute they can be treated as an instantaneous pulse followed by a finite period of evolution under the background hamiltonian . for a single nmr spin this is not a problem , as the maximally mixed state does not evolve under the background hamiltonian ; for a two qubit system , however , it is necessary to consider this problem more carefully . we now turn to the problem of implementing a full twirl on two qubits . this is clearly related to the problem of averaging a single qubit , but is more complicated . in particular , any twirl procedure will average a single qubit , but not every averaging procedure can be converted to a twirl . for example , the sequence of four operations @xmath45 suggested by bennett _ @xcite does average a single qubit , but does not effect a twirl ; instead it is necessary to use a set of 12 operations @xcite corresponding to the rotational symmetry elements of the tetrahedral group @xmath46 @xcite . the group theoretical justification for using these rotations is discussed in @xcite .    turning to the other averaging methods discussed above , only the method based on random euler rotations ( @xmath17 ) will correctly implement a twirl ; the other methods will only partially average the state . in particular the sequence of two random rotations around orthogonal axes has already been studied @xcite , and shown to implement a partial twirl , in which the density matrix is reduced to bell diagonal form , but not to a werner state .    for an ensemble quantum information processor , it might seem that the best strategy for performing a full twirl would be to perform random euler rotations . unfortunately this approach can not be implemented using gradients or similar methods , as these generate rotations with a uniform distribution ( between 0 and @xmath14 ) _ for each rotation angle _ , while the euler rotation method requires that the angles @xmath10 and @xmath11 be uniformly distributed _ over a sphere_. thus it might seem that twirl operations can not be easily implemented by ensemble methods . it is , however , possible to do this by using sequential random rotations around three different axes . this procedure is a development of the partial twirl described above . consider again the sequence of rotations in equation [ one - qubit - crot-2axes ] . these two rotations constitute a partial twirl @xcite , and leave the density matrix bell diagonal with equal populations of @xmath47 and @xmath48 , where @xmath49 . the role of the third rotation is to scramble these states with the sole remaining undesired term , @xmath50 , leaving a werner singlet state . this can be achieved by rotating about an axis at an angle @xmath51 , commonly known in nmr studies as the `` magic angle '' @xcite . our final twirl sequence is then @xmath52 as was the case for the partial twirl @xcite it is necessary to choose the lengths of the gradient pulses as multiples of @xmath53 , where @xmath54 is the difference between the larmor frequencies of the two spins , to refocus evolution under the background hamiltonian , which occurs at this difference frequency . if desired this twirl can be implemented using discrete steps by replacing each random rotation with a three step averaging procedure , @xmath55 where @xmath56 indicates a set of @xmath6 rotations with rotation angles @xmath57 , where @xmath58 . this discrete twirl requires twenty - seven steps , but the number of steps can be reduced to eighteen , using the sequences @xmath59 or @xmath60 it is not possible to reduce the number of steps further while using sequential rotations around three axes . it is instructive to compare this discrete twirl with the twelve step version of bennett _ et al . their twelve rotations include the identity operation , @xmath61 rotations about the principal @xmath62 , @xmath42 and @xmath6 axes and rotations by @xmath63 around the four body diagonals , each of which is at the magic angle with the cardinal axes . ( these are the twelve rotations that leave a tetrahedron invariant . ) clearly our rotations are closely related to theirs ; the use of sequential rotations means that the number of steps required is larger , but the sequential approach may be simpler to implement in practice . finally , we demonstrate the experimental implementation of our ensemble twirl operation as described in equation [ 2-qubit-3-axes - gzonly ] . our implementation is based on a homonuclear system , comprising the two nuclei of cytosine dissolved in @xmath64 ; we use the product operator description @xcite for the nmr states and pulse sequences , and label the two spins as @xmath65 and @xmath66 . the hamiltonian of our system is then @xmath67 where @xmath68 and @xmath69 are the resonance offset frequencies of spins @xmath65 and @xmath66 , and @xmath70 is the @xmath71 spin  spin coupling constant ( assuming weak coupling ) , all measured in hertz . nmr experiments were performed on a varian unity inova spectrometer with the spectrometer frequency placed between the two cytosine resonances ; for our system @xmath72 , @xmath73 and @xmath74 . the experiment began with the preparation of a generic mixed state of the two qubits containing a contribution from the singlet state as well as a wide range of other terms . the preparation sequence ( which was applied to the thermal equilibrium state , @xmath75 ) involves shaped pulses selectively exciting the spins and separated by a delay @xmath76 , @xmath77 shaped pulses were implemented using strongly modulated composite pulses as described by fortunato _ the deviation density matrix of the prepared state @xmath0 has a singlet component proportional to @xmath78 clearly the amount of singlet in @xmath0 can be easily controlled through choice of the delay @xmath76 . note that as we are using a deviation density matrix description of a highly mixed state it is possible for the singlet component to be negative . we performed two sets of experiments on states drawn from this family of states . the first set uses a fixed initial state and demonstrates a stepwise progression through the twirl sequence , with the resulting state becoming closer to a werner state at each stage . the second set uses a range of initial states , with different singlet fractions , and shows that the twirl sequence works equally well over this range . .pulse sequences for the first set of experiments , which demonstrate the effects of the three stages of our twirl operation . see main text for details . [ cols="^ , < " , ]     the first set of experiments , listed in table [ experiment1 ] , comprises a stepwise progression through the twirl sequence , equation [ 2-qubit-3-axes - gzonly ] . each experiment uses the same initial state , with @xmath79 , a value chosen to maximize the singlet fraction . the three stages of the twirl correspond to the application of three different crush gradients , and can be characterized as follows : stage 0 , do nothing ; stage 1 , apply a single crush gradient ; stage 2 , apply the partial twirl sequence @xcite ; stage 3 , apply the full twirl sequence . two different measurements are performed to characterize the state after each step : simple direct acquisition ( which should show no signals for a werner singlet state ) and acquisition after a selective excitation pulse , which converts the singlet state into nmr observable terms . the selective pulse was implemented using the sequence @xmath80 - 90_{180},\ ] ] where @xmath81 is the frequency separation between the two resonances as before . this sequence is based on jump - and - return sequences described in previous work @xcite , and when applied to a singlet state , results in the observable nmr terms , @xmath82 corresponding to a pair of antiphase doublets with equal and opposite intensities . the resulting spectra are shown in figure [ graphics1 ] . the left hand column of this figure shows spectra from direct acquisition and the right hand column shows spectra obtained with the excitation pulse followed by acquisition . the four rows correspond to observation after 0 , 1 , 2 or all 3 stages of the twirl sequence . the initial state @xmath0 contains many different components , and the observed spectra ( top row ) are complicated , reflecting this fact . after the first stage of the twirl sequence ( the first field gradient ) all components which are directly observable by nmr are crushed ( averaged to zero ) , and so no signal is visible in the direct detection spectrum , but many other components remain , indicated by the variety of signals seen after excitation .    the second stage of the twirl ( the @xmath83 pulse and the second gradient ) removes most of these terms , producing a bell diagonal state with equal populations of @xmath47 and @xmath48 . as before there is no signal in the left hand spectrum , while the spectrum on the right contains two antiphase doublets with clearly different intensities . this intensity difference arises from the imbalance between the @xmath84 state and the two @xmath85 states @xcite .    after the third and final stage of the twirl ( the @xmath86 pulse and the last gradient ) these two antiphase doublets have equal and opposite intensities , characteristic of a werner singlet state . the final intensity of these doublets , compared with a calibration spectrum ( not shown ) acquired from the thermal state , is consistent with the twirl preserving the fraction of the singlet state as expected . , in steps of @xmath87 , which causes the fraction of singlet in the initial state to be modulated , with a period of ten steps . ] this last point is explored in more detail in the second experiment , in which the full twirl is applied to a range of initial states with different amounts of the singlet state . the size of the singlet component for our family of states is given in equation [ singlet - fraction - dirty ] , and shows two kinds of sinusoidal modulation with the variable delay @xmath76 : a fast variation , arising from the offset frequency @xmath68 , on top of a slow variation due to the coupling @xmath70 between the two qubits . if we choose @xmath76 to be close to @xmath88 , we are near the maximum of the slow @xmath70 modulation , and the effect of varying @xmath76 is dominated by a sinusoidal variation arising from @xmath68 .    in this way we can produce a range of density matrices @xmath0 with varying amounts of singlet , together with other terms , each of which can be twirled to produce a werner state . this state can then be observed using a selective excitation pulse prior to acquisition . in each case the expected spectrum is a pair of antiphase doublets , with the intensity of the signal showing a sinusoidal modulation at the frequency @xmath68 . the results of this experiment are shown in fig . [ graphics2 ] , which depicts the intensity variation in the righthand component of the lefthand doublet as @xmath76 is varied around a value of @xmath89 , with an increment between successive spectra of @xmath87 ; equivalent effects can be seen for the other three components of the nmr signal . as expected a sinusoidal modulation of the signal is clearly seen , and the observed modulation period of ten spectra is exactly as expected . the small out of phase signals observable near the zero - crossings of the sine wave can be ascribed to the effects of spin  lattice relaxation during the gradient pulses . we have described a variety of strategies for the practical implementation of twirl sequences on conventional and ensemble quantum computers , and have demonstrated an ensemble implementation on an nmr quantum computer . with a conventional quantum computer the implementation requiring the smallest number of different bilateral operations is the set of 12 rotations previously described @xcite , but our 18 step and 27 step averaging procedures require a smaller number of elementary operations and may be simpler to implement in practice . with an ensemble quantum computer , such as an nmr device , it can be simpler to replace the discrete averaging procedure by continuous averaging , exploiting the ensemble nature of the system . we have developed a scheme , involving the application of three successive crush gradients separated by rf pulses , which is well suited to nmr quantum computers , and have demonstrated that its experimental performance is consistent with theoretical expectations .
Abstract:

Answer: twirl operations , which convert impure singlet states into werner states , play an important role in many schemes for entanglement purification . in this paper we describe strategies for implementing twirl operations , with an emphasis on methods suitable for ensemble quantum information processors such as nuclear magnetic resonance ( nmr ) quantum computers . we implement our twirl operation on a general two - spin mixed state using liquid state nmr techniques , demonstrating that we can obtain the singlet werner state with high fidelity .


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: it is now 45 years since cocconi @xcite drew attention to the broadening effects of the geomagnetic field on the lateral development of electron - photon cascades in the atmosphere . allen @xcite discussed the interaction of the geomagnetic field and the cascade electrons / positrons as the origin of the radio frequency emission produced by large cosmic ray showers . et al _ @xcite suggested that the geomagnetic separation of muons could be used to estimate the height of origin of showers . porter @xcite and browning and turver @xcite made calculations for the production of erenkov radiation by electron - photon cascades in the atmosphere , with allowance for the effects of the geomagnetic field . ground based gamma ray astronomy exploits the magnifying effects of the atmosphere to enable the detection of very low fluxes of gamma rays from a variety of sources @xcite . the cascade of electrons initiated by a primary gamma ray produces , at an altitude of about 10 km , erenkov light which reaches the ground in a pool @xmath0 m in diameter . it is possible to detect gamma rays of energy @xmath3 300 gev via this erenkov light produced in the atmosphere with an effective collection area , for our mark 6 telescope , of @xmath4 m@xmath5 . bowden _ et al _ @xcite discussed the effect of the geomagnetic field on the performance of ground based gamma ray telescopes . the interaction of the field and the cascade electrons produces a broadening of the atmospheric erenkov light image resulting in a reduction in the density of light sampled by the telescope ; so the energy threshold for the telescope increases . this is observed in the higher count rate for a telescope detecting showers propagating along the lines of the field ( with no spreading ) than when observing cascades developing perpendicular to the field lines ( and being spread ) , all other factors being the same . typical differences in measured count rate were about 20% in these extreme cases . the possibility was noted , on the basis of simulations , that an associated rotation of the direction of the erenkov light images may occur for cascades developing under high magnetic fields and in unfavourable directions . this could be of importance in ground - based gamma ray studies since the orientation of the image in gamma ray initiated cascades is the key to rejection of @xmath3 99% of the charged cosmic ray background @xcite . lang _ et al _ @xcite showed that the measurements of tev gamma rays from the crab nebula using the imaging erenkov technique were not significantly affected when the magnetic field was @xmath6 0.35 g.    we report measurements made using a ground based gamma ray telescope of cascades developing under the influence of fields up to 0.55 g. observations with the mark 6 telescope operating in narrabri , australia which is discussed by armstrong _ et al _ @xcite are subject to such magnetic fields when observing objects to the south . the observational data demonstrate all the effects of the geomagnetic field on cascades which are predicted by simulations . the lateral distribution of electrons and positrons near the maximum of an electromagnetic shower in air is approximately axially symmetric . the presence of a transverse component of magnetic field will cause the electrons and positrons to separate and the lateral distribution to become wider along the radial direction perpendicular to the transverse field component . the order of magnitude of this effect may be estimated . a relativistic electron of energy @xmath7 ( mev ) will follow a circular path in a transverse magnetic field of strength @xmath8 ( g ) with a radius of curvature @xmath9 ( km ) of @xmath10 . the energy of the typical erenkov light - producing electron at shower maximum is @xmath11 100 mev ( the critical energy in air ) , which is at a height of @xmath11 10 km for a 250 gev primary gamma ray . the radius of curvature of these electrons for a maximum transverse magnetic field of 0.56 g is @xmath11 6 km . the typical length of such an electron s path in air at 10 km altitude is @xmath11 1 km , giving a lateral deflection of @xmath11 80 m. at a height of 10 km this corresponds to an angular deflection of @xmath12 , the same order as the coulomb scattering width . it should be possible to detect this as a broadening of the shower along one axis . the effect of the broadening will depend on the orientation of the image with respect to the direction of the magnetic field . the direction of the shower axis in the atmosphere is within @xmath111@xmath13 of the telescope - source line . the transverse magnetic field component therefore lies very close to the plane of the focal image of the shower . the parameters of an image which are used to discriminate between gamma rays and hadrons depend primarily on the second spatial moments . in a coordinate frame in which the x - axis is aligned with the projected field direction , these are @xmath14 , @xmath15 and @xmath16 . the orientation angle of the long axis of the image is . those shower images whose major axes are aligned with the projected direction of the magnetic field have . the moment @xmath17 will be expected to be increased by the magnetic field and the images will be widened only . for images at other orientations , and the result of an increase in @xmath18 will be to rotate their minor axes towards the projected field direction . the main parameter used to discriminate gamma rays is _ alpha _ , the angle contained within the long axis of the image and the radius vector from the source position to the centroid of the image . any rotation of the image may affect _ alpha _ and hence the sensitivity for gamma ray detection . this is discussed further in section  [ simulations ] . monte carlo simulations of gamma ray and hadron cascades have been made to investigate the strength of this geomagnetic effect . these simulations have been performed using a monte carlo code developed from that used in our earlier work @xcite which incorporates the responses of the university of durham telescopes . the results of the simulations have been validated against established simulation codes e.g. mocca @xcite . agreement has been found in all cases . the consequence of an increase in the lateral and angular spread of the shower electrons is a decrease in the density of erenkov light on the ground . this leads to a decrease in counting rate for fixed threshold , which was reported by us from measurements , using a non - imaging telescope , of counting rate against azimuth at fixed zenith angle @xcite . figure  [ fig : simtotal ] shows the simulated reduction in the counting rate of a gamma ray telescope of the type used by us for hadron primaries , relative to the zero - field case , caused by a transverse magnetic field of 0.5 g. this is consistent with our earlier measurements  a reduction of up to @xmath19 in count rate for gamma ray energies @xmath20 gev @xcite and , apparently , a greater reduction for @xmath21 gev cascades .    in our earlier work @xcite we suggested that the distribution of recorded photons in an imaging camera may be distorted . we show in figure [ fig : simimage ] the pattern of photons in individual simulated 1 tev gamma ray cascades developing under the effects of 0 and 0.56 g geomagnetic fields . we also show the predicted responses of our camera after allowing for optical distortion , pixellation of the camera , measurement noise and passing through our standard data analysis package . the simulations for the pairs of cascades developing under the influence of a field and no field have been initiated using the same random number seed ; this ensures , as far as is possible , that the early development of the cascade ( which dominates the resulting image and is nearly independent of field ) is identical for the two simulations . the expected effects on the shape of the image ( widening or lengthening and rotation , depending on the orientation of the image with respect to the magnetic field direction ) was apparent , both in the distribution of photons hitting the detector and in the processed image . the predicted width of a shower image was investigated as a function of magnetic field strength and angle for gamma rays of energy 500 gev . results are shown in figure  [ fig : simswid ] for transverse field strengths of 0.0 , 0.2 and 0.5 g. no allowance is made for the broadening effects of the defects in the mirror or the pixellation of the recording pmt camera . the effect of a 0.5 g field on the width ( minor axis ) of such undistorted gamma ray images is seen to be greatest when the image minor axis is perpendicular to the projected field direction , as expected . small transverse fields of 0.2 g have a limited effect on the image width . however , the effect on the pointing angle _ alpha _ become noticable and may vary both with the image s angle to the magnetic field and with its eccentricity ( defined as the ratio _ width_/_length _ ) . for example , an image with very low eccentricity may be widened only by a small amount but rotated through a large angle . histograms in _ alpha _ are shown in figure  [ fig : simsalpha ] for the simulated gamma rays , for which widths were shown in figure  [ fig : simswid ] , under the influence of fields of 0.0 , 0.2 and 0.5 g.    the effects of mirror quality and camera pixellation of a typical telescope may reduce this change in width but may not mask the change in orientation of the image . we show in figure  [ fig : tracks ] the celestial sphere , as viewed from the site of the mark 6 telescope in narrabri , australia , showing the loci of constant magnetic field strength . we indicate the tracks on the sky for the telescope during observations of a number of potential gamma ray sources which have been intensively studied . in fact , most of the sources studied at narrabri are in the azimuth range @xmath22 and in directions for which the magnetic field is @xmath3 0.35 g and in many cases @xmath3 0.5 g. the telescope sensitivity is therefore likely to be reduced by the broadening of the distribution of _ alpha _ for many of our observations , unless corrections are applied . all data reported here have been obtained with the mark 6 vhe gamma ray telescope which has been described in detail by armstrong _ et al _ @xcite . the telescope measures gamma rays in the energy range 300 gev  10 tev , with a 50% trigger probability at an energy of @xmath0 gev . it comprises an alt - azimuth mount with three 7 m diameter f/1.0 aluminium mirrors with an rms spread for a point source of @xmath23 . the pixel size for a @xmath24 diameter photomultiplier is @xmath25 . the focal plane of the central mirror contains a 109-pixel camera comprising a close - packed hexagonal array of 91 1-inch diameter photomultipliers surrounded by a ring of 18 2-inch diameter photomultipliers . the field of view of the imaging camera is @xmath26 wide . the mirrors on the left and right of the mount each contain a triggering detector of 19 2-inch photomultipliers in a close - packed hexagonal array . each photomultiplier is contained in a mumetal shell designed to give magnetic and electrostatic shielding . the effect of the geomagnetic field on the performance of the photomultipliers , with and without shielding , has been studied @xcite . for the 1-inch hamamatsu r1924 pmt the effect of a 0.5 g field on the gain was not measureable and less than 1% . for the 2-inch phillips xp3422 pmt the measured decrease in gain due to a 0.5 g field was 1% . shielded photomultipliers of both types show no measurable effects due to the magnetic field . the control of the attitude of the telescope is via dc servomotors driving onto gears mounted directly on the telescope structure . angles are sensed by absolute digital 14-bit shaft encoders with a resolution of @xmath27 . the calculated azimuth and zenith of a source is compared by a digital servomechanism ( employing 12-bit resolution ) with the shaft encoder outputs at 100 ms intervals . the error signals are passed via dacs to the dc motor amplifiers . these provide damping on acceleration and stabilise the movements of the telescope structure . thus although the telescope pointing is known to a resolution of @xmath28 , the source can be offset from the camera centre by up to @xmath29 by this mechanism . the attitude of the telescope is measured in two ways . the shaft encoder positions are recorded for each event to 14-bit accuracy . this gives a measurement of the pointing to @xmath30 within the telescope system . in addition , a coaxial optical ccd camera is mounted on the telescope . the output of this ccd camera is continuously monitored by microcomputer , which measures the position and brightness of a nominated guide star within the @xmath31 field . this information is integrated into the data stream on an event - by - event basis . guide stars of magnitude @xmath32 can be employed , providing absolute position sensing for the telescope to better than @xmath33 . false source analysis has been shown to be a useful method of demonstrating that @xmath34-ray like events originate from the source direction @xcite . conversely , such analyses may also be used to verify the steering performance of a telescope , if an established @xmath34-ray emitter is observed . however , we have only one data set for a gamma ray source containing data which are _ not _ subject to the effects of a strong geomagnetic field  see section [ pks2155alpha ] . most of the data considered here were taken during routine observations of potential gamma ray sources during 1996  1998 . a small amount of data was taken in dedicated measurements at fixed values of azimuth and zenith , corresponding to a range of values of the geomagnetic field . the only restriction imposed on data was that events should lie within @xmath35 of the centre of the camera ( to avoid edge effects ) and were large enough ( 5 times the triggering threshold ) to ensure that their shape was well measured . the investigation of the sensitivity of the count rate to the magnetic field which established the effect with the mark 3 and mark 4 telescopes @xcite has been repeated for the mark 6 telescope . we find that the count rate observed with the mark 6 telescope for cascades developing under minimum and maximum values of the geomagnetic field at a constant zenith angle of @xmath36 differs by about 15% , as expected . the widths of images recorded by the mark 6 telescope in directions parallel and perpendicular to the magnetic field have been investigated . for a large sample of cascades measured during observations of a range of sources the differences of the mean width in the parallel and perpendicular directions are shown in figure [ fig : distortion ] for a range of values of distorting field . the widths of images in directions perpendicular to the field are significantly greater than those in parallel directions for fields @xmath2 g. although the differences in widths are small , because of the effects of pixellation and noise which are common to all data , the measurements are free from systematic effects . this is because the orientation of images to the magnetic field depends on the orientation of the image in the camera . measurements of the width of the images parallel and perpendicular to the magnetic field are derived from events distributed throughout the same observation . we have attempted to measure the differences in widths of images in cascades measured at fixed zenith angles of @xmath36 for values of fields less than and greater than 0.35 g. such measurements are difficult and require large exposures ( @xmath37 hr ) to give significant differences . for example , with exposures @xmath38 hr we find for @xmath39 g , the mean width is @xmath40 ( sd @xmath41 ) and for @xmath42 g the mean width is @xmath43 ( sd @xmath44 ) . the small difference in means , significant at the @xmath45 level , demonstrates the difficulty in making such a small measurement for which the standard deviation is approximately similar for all @xmath46 .      in the absence of any magnetic field and detector triggering biases , the distribution of the orientation of the images in the camera of a telescope will be isotropic . if we define a coordinate system where the angle @xmath47 is the angle of the long axis of an image to the horizontal in the camera frame , then the distribution of @xmath47 should be flat between @xmath48 and @xmath49 . minor deviations from uniformity in angle might be expected near to threshold because of the increased effect of small changes in triggering probability and the six - fold symmetry of the hexagonal arrangement of close - packed detectors . we show in figure  [ fig : aysmallh](a ) the distribution in the angle @xmath47 for background hadronic events which are subject to small values of the transverse magnetic field ( 0.15 g ) . we show in figure  [ fig : aysmallh](b ) the distribution in @xmath47 for images due to hadron - induced cascades recorded in special observations at fixed azimuth and zenith angles which gave a transverse geomagnetic field of 0.52 g. all of these events were recorded in dedicated observations with the telescope at a fixed zenith angle ( @xmath36 ) . data at different values of the transverse magnetic field were obtained by varying the azimuth angle . all data presented here were recorded within a period of one hour , so potential variations due to changes in camera performance , atmospheric clarity , etc . were minimised . the data were processed following our normal procedures ( see e.g. @xcite ) . the requirement that images have a minimum brightness and fall within the camera ensures that the data are free of effects of variations in night sky brightness . the data for @xmath50 g show the expected distribution indicative of a near isotropic distribution of directions in the camera with a small peak . we note a strong anisotropy , with a peak containing twice the number of events resulting from the skewing or distortion of the image , for events subject to a @xmath51 g field . the maximum of the distortion occurs at @xmath52 , which is appropriate for those observations which were made at magnetic south . the magnitude of the anisotropy displayed in figure  [ fig : aysmallh](b ) should depend on the strength of the magnetic field . in figure [ fig : ayampvh ] we plot the amplitude of the peak of the distortion in the @xmath47 distribution as a function of the magnetic field strength . each point corresponds to the result for a 15 minute segment of data . note that for values of transverse component of the geomagnetic field less than 0.35 g there is no great distortion of the @xmath47 distribution , as suggested by the work of lang _ et al _ @xcite , but for values of field in excess of 0.4 g substantial distortion occurs . the position of the peak in the @xmath47 distribution depends on the angle between the projected magnetic field and the vertical direction in the camera  @xmath53 . we show in figure  [ fig : aydirvh ] the correlation between @xmath53 and the position of the peak in @xmath47 . the data demonstrate the expected relation between these angles . again , each point is based on a 15 min sample of data taken during routine observations of a range of potential gamma ray sources . we have evidence that the geomagnetic field should , and does , influence the lateral development of atmospheric cascades . the observed effects on the threshold of a telescope and on the shape of the image and the magnitude and phase of the distortion of the pointing of images produced by background cosmic ray protons are as expected . a number of potential gamma ray sources has been observed using the mark 6 telescope at narrabri . in several cases there is evidence for gamma ray emission @xcite . the evidence comes mainly from a comparison of the _ alpha _ distributions for data selected on the basis of image shape for the on - source and off - source scans . the difference between the _ alpha _ distributions should show an excess of events  the gamma ray candidates  at low values of _ alpha_. if the data were taken in geomagnetically unfavourable directions , as is the case for most of our data , it might be expected that the _ alpha _ distribution of the excess events would be wider than that for data taken in more favourable geomagnetic directions . consideration of data recorded from a gamma ray source which are subject to @xmath54 g may provide the only true indication of the _ alpha_-distribution for our telescope . observations of this object were made with the cascades recorded over a range of transverse geomagnetic field strengths between 0.25 and 0.5 g. the total data set contained 41 hrs of observation , as reported @xcite . the _ alpha _ plot for the difference between the on - source and off - source data taken at zenith angles @xmath55 is shown in figure  [ fig : alphaplot2](a ) . the significance of the excess at _ alpha _ @xmath56 is @xmath57 according to the maximum likelihood method @xcite . we are able to select a subset of data for which the strength of the projected field to which the cascades were exposed was @xmath6 0.35 g and for which minimal distorting effect would be expected . the _ alpha _ plot for this subset is shown in figure  [ fig : alphaplot2](b ) . the excess events all have _ alpha _ @xmath58 ( significance @xmath59 ) . the distribution is narrower than that of the total dataset and is typical of that expected for a @xmath60 pixel camera . ( in the absence of any other data for gamma rays detected with our telescope and _ not _ subject to the effects of the magnetic field , we assume that this is reasonable . )    the _ alpha _ plot for the majority of the events for which the field is @xmath3 0.35 g , is shown in figure  [ fig : alphaplot2](c ) . it is evident that the width of the peak is larger for these events recorded under the influence of higher transverse magnetic fields , with equal populations for values of _ alpha _ between @xmath61  @xmath62 ( @xmath63 ) and @xmath62  @xmath64 ( @xmath65 ) . it should be noted that the peak at @xmath66  @xmath67 is superimposed on an _ alpha_-plot with increasing frequency for _ alpha _ approaching @xmath68  see @xcite . the significance of this peak is therefore @xmath69 ( before allowing for the number of bins in the _ alpha_-plots ) .    on the basis of our only measurement of a gamma ray source not subject to the effects of the geomagnetic field , we have evidence that the excess events with @xmath70 _ alpha _ @xmath56 are confined to those measurements made with fields @xmath71 g. this conclusion is significant at the 2.5% level , on the basis of a @xmath72 contingency test of the populations of the @xmath73 bins in figures  [ fig : alphaplot2](b ) and [ fig : alphaplot2](c ) . we have demonstrated that , on the basis of simulations and measurement , the erenkov images from gamma rays and cosmic rays are broadened and rotated by the geomagnetic field . the broadening results in an increase in the telescope threshold and a reduction in counting rate . the rotation of the images away from the projected direction of the magnetic field in the image plane broadens the _ alpha _ distribution . these results suggest that the geomagnetic field can have an important effect on the operation of atmospheric erenkov telescopes in some directions and that , for detected and candidate sources , our narrabri site is particularly susceptible . these effects may be removed using appropriate correction techniques ; unlike noise , geomagnetic effects do not reduce the information contained in the image . we are grateful to the uk particle physics and astronomy research council for support of the project and the university of sydney for the lease of the narrabri site . the mark 6 telescope was designed and constructed with the assistance of the staff of the physics department , university of durham . 99                                  gibson a i , harrison a b , kirkman i w , lotts a p , macrae j h , orford k j , turver k e and walmsley m 1982 , _ proc . int . workshop on very high energy gamma ray astronomy _ , ed p v ramana murthy and t c weekes , ( bombay : tata institute ) p 97
Abstract:

Answer: atmospheric erenkov telescopes are used to detect electromagnetic showers from primary gamma rays of energy @xmath0 gev @xmath1 tev and to discriminate these from cascades due to hadrons using the erenkov images . the geomagnetic field affects the development of showers and is shown to diffuse and distort the images . when the component of the field normal to the shower axis is sufficiently large ( @xmath2 g ) the performance of gamma ray telescopes may be affected , although corrections should be possible .


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: the first accretion - powered millisecond pulsar , sax j1808.4 - 3658 , was discovered a decade ago @xcite . since then several other millisecond x - ray pulsars with spin periods between 1.7 and 5.4ms have been discovered ( e.g. xte j0929 - 314 , xte j1751 - 305 , xte 1807 - 294 and xte j1814 - 338 ; see @xcite for a review ) . these neutron stars have considerably weaker magnetic dipole moments , @xmath0tm@xmath1 , than the conventional x - ray pulsars . the accretion disc then extends much closer to the neutron star , where it reaches such a high temperature that the opacity is dominated by electron scattering and radiation pressure can grow stronger than the gas pressure . on the other hand the accretion rates in these systems can be so low that the alfvn radius might be located outside of the co - rotation radius of the neutron star . conventionally it has been thought that such systems will eject the mass that is transferred from the companion star through the propeller mechanism @xcite , though it has been pointed out by @xcite and @xcite that the velocity excess right outside the co - rotation radius is insufficient to eject all the matter . @xcite found that when the alfvn radius is only slightly larger than the co - rotation radius then the disc will cycle between one state in which the matter is accumulated outside the co - rotation radius and another state in which it is accreted onto the star . later @xcite constructed a steady accretion disc model , in which the entire disc is located outside of the co - rotation radius . the magnetic field of the neutron star is penetrating this disc and is mediating an exchange of angular momentum between the accretion disc and the neutron star as in the model by @xcite . since the entire accretion disc is rotating more slowly than the neutron star in this case , the disc brakes the rotation of the neutron star . the models of @xcite and @xcite have ignored the possibility of an internal dynamo in the accretion disc , and have considered only the magnetic field that is generated in the disc by the mismatch in the rotation of the neutron star and the disc . however the magnetic field due to a dynamo can strengthen the coupling between the accretion disc and the neutron star , and we @xcite have shown that the angular momentum exchange between the accretion disc and the neutron star can increase by more than an order of magnitude compared to the standard model by @xcite . furthermore the dynamo is insensitive to the rotation of the neutron star , and can therefore in principle provide a spin - up torque even at a low accretion rate .    a complementary approach to study the magnetic interaction between an accretion disc and a star is represented by the numerical simulations by @xcite . in these simulations it is possible to study more complicated magnetic field geometries than are accessible to semi - analytic studies ( see for instance @xcite ) , but the resolution is so low that it is not possible to study the appearance of small - scale turbulence and the resulting turbulent transport and dynamo action in the simulation . thus these simulations do not really address the viability of a dynamo in the disc . however there are other things one can learn from these simulations . they show how the matter is streaming via funnel flows from the accretion disc onto the magnetic poles of the neutron star . furthermore they show that angular momentum is also lost from the accretion disc through the jets that are accelerated away from the disc . @xcite make a distinction between a conical jet that transport most of the mass , and a less dense axial jet inside the conical jet . the latter component can dominate the energetics if the star is in the propeller regime . it is also noteworthy that these simulations confirm that some matter might still be accreted when the star is in the propeller regime . the purpose of this paper is to extend our previous study ( paper i ) of an accretion disc with an internal dynamo around a magnetised star to the discs around millisecond pulsars . the difference is that due to the weak magnetic field the disc extends to much smaller radii , at which it becomes so hot that a significant part of the opacity is due to electron scattering , and radiation pressure might dominate over the gas pressure . therefore in the spirit of @xcite we divide the disc into an outer region with gas pressure and free - free opacity , a middle region with gas pressure and electron scattering opacity , and an inner region with radiation pressure and electron scattering . in sect . 2 we derive a single ordinary differential equation , which is generally applicable for the radial structure of the accretion disc independently of the equation of state and the form of the opacity . we then introduce the equations of state and the opacity prescriptions for the different regions of the disc in sect . 3 and present the numerical solutions in sect . the properties of these solutions and their application to the millisecond x - ray pulsars are discussed in sect . 5 , and we finally summarise our results in sect . our approach is an extension of the standard model for a thin accretion disc ( @xcite ) , and it follows closely the method that we introduced in paper i , though we will now consider a general opacity law and equation of state . for this reason we only summarise the steps leading up to eq . ( 21 ) of paper i , after which we deviate from the treatment in paper i.    the surface density of the accretion disc is defined as @xmath3 where @xmath4 is the density , and @xmath5 the half thickness of the disc . we can then write the accretion rate as @xmath6 which is independent of the radial coordinate @xmath7 , and where @xmath8 is the radial velocity . a thin accretion disc is supported by the centrifugal force in the radial direction , so we get that @xmath9 where @xmath10 is the azimuthal velocity , @xmath11 the gravitational constant and @xmath12 is the mass of the accretor . this shows that the disc rotates in a keplerian fashion . the disc is in hydrostatic equilibrium in the vertical direction , which gives us that the pressure at the midplane of the disc is @xmath13 but the hydrostatic equilibrium can also be expressed as @xmath14 where @xmath15 is the isothermal speed of sound . this shows that the keplerian velocity is highly supersonic in a thin accretion disc . the conservation of angular momentum leads to the following height - integrated equation ( paper i ) @xmath16_{-h}^{h}r   + \frac{1}{r}\frac{\mbox{d}}{\mbox{d } r}\left[r^{3}\nu\sigma \frac{\mbox{d}}{\mbox{d } r }      \left(\frac{l}{r^{2}}\right)\right],\ ] ] where the specific angular momentum @xmath17 and @xmath18 is the viscosity . the magnetic term describes the exchange of angular momentum between the disc and the star via the magnetosphere . we assume that the vertical magnetic field is due to the dipolar field of the neutron star , so that its value in the equatorial plane is @xmath19 where @xmath20 is the magnetic dipole moment of the star . there are two sources for @xmath21 . the shear between the keplerian disc and the magnetosphere produces a field whose value in the upper half of the disc is @xmath22 where @xmath23 @xmath24 is the angular velocity of the star , and @xmath25 is a dimensionless parameter of the order of a few @xcite .    the magnetohydrodynamical turbulence in the accretion disc is also generating a magnetic field , @xmath26 , through a dynamo action ( see for instance @xcite ) . in order to estimate the strength of this magnetic field we start by making the assumptions that the turbulence is not affected by @xmath27 , which is justified because in general @xmath28-component of the turbulent maxwell stress tensor , @xmath29 , that transports the angular momentum , and we couple this to the pressure in the disc through the @xcite @xmath30-prescription @xmath31 note that a negative maxwell stress means that the angular momentum is transported outwards and that numerical simulations suggest that @xmath32 @xcite . @xcite suggested that than @xmath33 that was used by @xcite . ] @xmath34 where @xmath35 based on the numerical simulations by @xcite . we now assume that due to an inverse cascade this turbulent filed contains a large - scale field , @xmath36 @xcite , and we express the ratio of the large - scale field to the rms - value of the turbulent field as @xmath37 thus we can express the large - scale magnetic field , that couples to the magnetic field of the neutron star , as @xmath38 where @xmath39 , and a negative value describes a magnetic field which is pointing in the negative @xmath40-direction at the upper disc surface . we will in general be conservative and assume that @xmath41 , but the simulations by @xcite suggest that @xmath42 can be significantly larger . we can now re - express eq . ( [ ang_mom ] ) as @xmath43.\end{aligned}\ ] ] which can be formulated as an ordinary differential equation in @xmath44 ( paper i ) .      for a slow inflow of matter through an optically thick disc the local viscous dissipation @xmath45 is balanced by the radiative losses @xmath46 . this leads to the height - integrated equation @xmath47 where @xmath48 is the temperature at the midplane of the disc , @xmath49 the stefan - boltzmann constant , and the optical depth of the disc is given by : @xmath50 where @xmath51 is the sum of the electron scattering opacity @xmath52 and the free - free opacity , which is given by kramer s law @xmath53 with @xmath54 combining eqs . ( 14 ) and ( 15 ) , we get @xmath55      the total pressure is the sum of gas and radiation pressure @xmath56 where @xmath57 is boltzmann s constant , @xmath58 the mean molecular weight for a fully ionised gas , @xmath59 the mass of a proton , and @xmath60 the speed of light , but the pressure can also be expressed using eq . ( [ hydro - eq ] ) for hydrostatic equilibrium . we can thus express the scale height in terms of the total pressure @xmath61 the viscous stress tensor ( paper 1 , eq . 30 ) gives us @xmath62 but pressure and density are related according to eq . ( [ hydro - eq ] ) , so we can solve for the density of the gas @xmath63    combining eqs . ( [ hydro - eq ] ) , ( [ ang - mom1 ] ) and ( [ density ] ) and using @xmath64 we get @xmath65^{1/2}(gm)^{-1/4}r^{-3/2 } \nonumber\\ - \frac{4\mu^{2}\gamma}{3\mu_{0}\sqrt{gm } } r^{-9/2}\left[1 - \left(\frac{r}{r_{\rm c}}\right)^{3/2}\right],\end{aligned}\ ] ] where @xmath66 is the co - rotation radius , at which the keplerian angular velocity is the same as the stellar angular velocity . here @xmath67 is the spin period of the star and @xmath68 . as @xmath69 @xmath70 which we use as the outer boundary condition . we now introduce the dimensionless variable @xmath71 through @xmath72 and a dimensionless radial coordinate through @xmath73 where @xmath74 is the alfvn radius , which is given by setting the magnetic pressure equal to the ram pressure ( e.g. @xcite ) @xmath75 where @xmath76 represents the mass transfer rate in units of @xmath77 , and @xmath78 is the stellar magnetic moment in units of @xmath79tm@xmath1 . the alfvn radius should not be taken as the location of the inner edge of our discs , since this is a quantity that comes out of our solutions , but merely as a convenient length unit . furthermore we introduce the fastness parameter @xmath80 we can now write eq . ( [ diff_eq ] ) as : @xmath81^{1/2 } ( gm)^{-1/4}r_{a}^{-5/4}r^{-9/4 } \nonumber\\ - \frac{4\mu^{2}\gamma}{3\mu_{0}\sqrt{gm}\dot{m } } r_{a}^{-7/2}r^{-9/2 } \left[1 - \omega_{\rm s}r^{3/2}\right ] , % \left[1 - \left(\frac{r}{r_{c}}\right)^{3/2}\right],\end{aligned}\ ] ] in order to solve eq . ( [ diff_eq_2 ] ) numerically we need to determine @xmath5 for the different regions in the disc . with gas pressure and kramer s opacity we get @xmath82 and for the gas pressure and electron scattering opacity @xmath83 finally with radiation pressure and electron scattering the scale height is @xmath84 which is essentially independent of r. this comes from the fact that for radiation pressure @xmath85    we can now write down the self - similar structure of the different regions of the accretion disc in terms of the function @xmath86 . the toroidal magnetic field is given by the same expression @xmath87 everywhere in the disc , but for all other quantities we have to handle the different regions of the disc separately . in the outer disc where @xmath88 and @xmath89 we have that : @xmath90 @xmath91 @xmath92 @xmath93 @xmath94 @xmath95 @xmath96 @xmath97 the transition from the outer to the middle region occurs where @xmath98 . assuming that @xmath99 at this point we have that the transition radius is @xmath100 which is independent of @xmath101    in the middle disc , where @xmath88 and @xmath102 , we have that : @xmath103 @xmath104 @xmath105 @xmath106 @xmath107 @xmath108 @xmath109 @xmath110 the disc structure changes drastically in the inner disc where @xmath111 and @xmath112 , where we have that @xmath113 @xmath114 @xmath115 @xmath116 @xmath117 @xmath118 @xmath119 @xmath120 it is of particular significance here that @xmath121 . the transition radius between the middle and the inner region is estimated by equating the gas and radiation pressures in the middle region and approximating @xmath71 with @xmath122 @xmath123 we then see that the accretion disc right outside of the alfvn radius is dominated by radiation pressure only if @xmath124 this condition is not fulfilled for a conventional x - ray pulsar with a magnetic dipole moment of @xmath125tm@xmath1 @xcite , though it can be fulfilled for a millisecond x - ray pulsar . for our fiducial model we take a neutron star of @xmath126 with a magnetic moment of @xmath127 and a spin period of 4.8ms . we set the dimensionless parameters @xmath25 and @xmath128 to 1 and 10 , respectively , while @xmath129 @xmath25 , @xmath128 , and @xmath30 appear only in combinations with other parameters , and thus their values do not carry any special significance , but the parameters @xmath130 and @xmath131 appear on their own in the equation , and therefore influence the solution in unique ways . in particular @xmath131 influences the fastness parameter and the transition radii between the regions of the accretion disc . in this paper we consider three different accretion rates , @xmath132 , 0.12 and 1.5 . we list the corresponding fastness parameters and transition radii in tab . [ params ] . lllll + @xmath133 & @xmath134 [ m ] & @xmath135 & @xmath136 & @xmath137 +   + 0.012 & @xmath138 & 0.98 & + 0.12 & @xmath139 & 0.37 & 8.5 &  + 0.8 & @xmath140 & 0.31 & 52 & 6.0 + 1.5 & @xmath141 & 0.12 & 95 & 12 +   +      the accretion disc consists of only an outer region when @xmath143 , which corresponds to that @xmath144kgs@xmath145 for our fiducial model . assuming @xmath146 we get the solutions from the top to the bottom of fig . [ lambdafigout ] for @xmath147 and @xmath148 . we see here that @xmath149 at a finite @xmath150 for @xmath151,and@xmath152 , which we called case d in paper i , since @xmath153 at this radius ( fig . [ sigmafigout ] ) . in the other cases @xmath71 grows without a bound for decreasing @xmath150 , which we call case v in paper i , since then the viscous flux of angular momentum changes sign at a finite radius where @xmath154 ( fig . [ innerfigouter ] ) , which we take as the inner edge of the disc . inside of this radius the accretion flow is driven by the magnetic coupling to the neutron star . the disc has both an outer and a middle region when @xmath156 . we then use the analytical solution with @xmath157 for the outer region and solve eq . ( [ diff_eq_2 ] ) numerically for the middle region . for our calculations we take @xmath158kgs@xmath145 , which yields @xmath159 . the numerical solutions for the middle region are shown in fig . [ lambdafigmiddle ] for our fiducial neutron star with @xmath160 , and @xmath148 . we find a case d inner boundary for @xmath161 , while the other solutions have case v inner boundaries . the transition from case d to case v at the @xmath130 , that gives the smallest possible radius for the inner edge of the disc , which varies weakly with the accretion rate @xmath131 . the innermost region of the disc is dominated by radiation pressure and electron scattering when @xmath163 . we use two different accretion rates that fulfill this condition , @xmath164kgs@xmath145 and @xmath165kgs@xmath145 together with our fiducial neutron star in our calculations . we assume that @xmath157 in the outer and middle regions of the disc , and solve eq . ( [ diff_eq_2 ] ) for the inner region starting from @xmath166 our solutions for the inner disc region at the higher accretion rate and with @xmath167 , and @xmath168 are shown in fig . [ lambdafiginner ] . all these solutions have case v inner boundaries except the @xmath169 model . however because @xmath170 in this case we see that the density goes to infinity at the inner edge of the @xmath171 disc ( fig . [ sigmainn ] ) , at which point the flow stagnates ( fig . [ radialvelocityfiginn ] ) . in the case v solutions , on the other hand , we see that @xmath172 at a point which is inside the inner edge of the disc ( fig . [ innerfigin ] ) . in order for the neutron star to appear as an x - ray pulsar the accretion disc must be truncated by the stellar magnetic field ( tab . [ table1 ] ) at a radius @xmath173 outside of the neutron star . the inner edge of the disc can be either of type d ( @xmath174 ) or type v ( @xmath175 ) ( see paper i for an extensive discussion of these boundaries ) . in models without a dynamo ( for instance @xcite ) @xmath173 is comparable to or smaller than the alfvn radius , @xmath134 , but we find that the dynamo can make it significantly larger than @xmath134 . this opens up the possibility , as suggested by the referee , that a system can temporarily turn into an x - ray pulsar during the periods when the magnetic field generated by the dynamo is strong enough to make @xmath173 larger than the radius of the neutron star . sax j1748.9 - 2021 is such an intermittent millisecond x - ray pulsar @xcite , which is on for hundreds of seconds at a time .    in case @xmath173 is located in the outer or middle region of our disc solution we find that a type d boundary coincides with the point at which @xmath153 and @xmath176 though the accretion rate is finite . however a type d boundary in the inner part of the accretion disc behaves in the opposite way ; @xmath177 and @xmath178 . neither one of these behaviours is entirely physical and one must keep in mind that we have ignored the importance of the radial pressure gradient and the inertial term in eq . ( [ 3 ] ) ( see for instance @xcite ) , and we have not attempted to model how the plasma is entering the stellar magnetic field lines to form a funnel flow onto the magnetic poles of the neutron star . one further problem with the inner region of the disc solution is that it is known to be viscously and thermally unstable @xcite . this instability persists if the accretion disc is threaded by a stellar magnetic field @xcite , but so far it has not been investigated what happens if the accretion disc is sustaining its own magnetic field through a dynamo .      to understand the nature of the exchange of angular momentum between the accretion disc and the neutron star we multiply eq . ( [ ang - mom1 ] ) by @xmath179 and integrate it from @xmath173 , the inner radius of the disc , to a large radius @xmath180 . @xmath181 the left - hand - side is the difference between the angular momentum that is advected through the inner edge of the accretion disc and that which is fed into the disc at its outer edge and the right - hand - side describes the contribution of magnetic and viscous torques to the angular momentum balance . notice that the term @xmath182 describes the viscous tension at the outer edge of the disc , which will not be considered further in this paper . the material torque of the inner edge of the disc on the neutron star is @xmath183 it increases as the accretion rate increases ( tab.[table1 ] ) . then we have the magnetic torques on the neutron star , which we have divided into one part due to the shear @xmath184 \mbox{d}r    \end{aligned}\ ] ] and a second part due to the dynamo @xmath185 we split up the latter integral into one integral for each of the three different regions of the disc @xmath186 @xmath187 and @xmath188 finally we have the boundary terms describing the viscous stresses on the inner and the outer edges of the accretion disc . the viscous stress at the outer edge of the disc is responsible for removing the angular momentum that is transported outwards through the accretion disc , and does not affect the angular momentum evolution of the neutron star . the viscous stress at the inner edge of the disc on the other hand can result in a torque that transports angular momentum away from the neutron star @xmath189 this term is usually neglected in accretion disc theory since the standard accretion disc solution has a case d inner boundary . but it can transport angular momentum from the neutron star to the disc when we have a case v inner boundary . we report on all the components and the total sum of the torques on our fiducial neutron star in tab.[table1 ] . the main factor determining the total torque is the accretion rate . in paper i we found that the dominant contributor to the torque is the internal disc dynamo . this torque is still strong , but now @xmath190 and @xmath191 are comparable in strength , though these make opposing contributions to the net torque , while the contribution from @xmath192 is still small . there is therefore no simple dependence of @xmath193 on @xmath130 , though the dynamo is indirectly affecting @xmath190 and @xmath191 by causing the inner disc edge to move outwards compared to the @xcite model . @xmath191 is also highly sensitive to the nature of the inner disc boundary , and in particular it vanishes for a d inner boundary , which occurs at a sufficiently negative @xmath130 , in which case the net torque is determined by the balance between @xmath194 and @xmath190 .    [ cols="<,<,<,<,<,<,<,<,<,<",options="header " , ]      the largest spin up rate recorded so far for a millisecond x - ray pulsar is @xmath195hzs@xmath145 during the december 2004 outburst of igr j00291 + 5934 @xcite , but one should keep in mind that there are large uncertainties in the spin variations that have been reported for the millisecond x - ray pulsars . for instance @xcite reported @xmath196s between @xmath197 and @xmath198hzs@xmath145 for sax j1808.4 - 3658 , but @xcite noted that the measurements of this source are plagued by large variations in the pulse shape , and put an upper limit of @xmath199hzs@xmath145 on the spin variations at any instant . some of the timing noise in the millisecond x - ray pulsars can be explained by the motion of the hot spot on the neutron star @xcite , and may thus not correspond to real changes in the spin period of the neutron star . the spin variation corresponds to a torque @xmath200 where @xmath201 is the moment of inertia of the neutron star measured in @xmath202kgm@xmath203 , and @xmath204 is the spin derivative in units of @xmath205hzs@xmath145 . in the @xcite model , the spin variation of igr j00291 + 5934 would imply then an accretion rate of at least @xmath206kgs@xmath145 , which is uncomfortably large , but our tab . [ table1 ] shows that it is possible to enhance the torque severalfold by coupling the stellar magnetic field to the magnetic field generated by an internal dynamo in the disc , and through the fact that the inner edge of the accretion disc is then located outside the alfvn radius , which increases the amount of angular momentum that is advected onto the neutron star . thus it is possible that a lower accretion rate can produce the necessary torque , though we have not attempted to fit our model to this system .      a significant part of the luminosity is not generated in the accretion disc outside of @xmath173 , but rather it is released as the accreting matter falls onto the neutron star . we estimate that this luminosity is @xmath207 where we have assumed that both the potential energy difference between the inner edge of the accretion disc and the surface of the neutron star as well as the kinetic energy associated with the keplerian motion at the inner edge of the accretion disc is released at the surface of the neutron star . consequently @xmath208 . the neutron star can then contribute to the local heating in the accretion disc through its irradiation of the disc . assuming that the neutron star is a point source and the accretion disc is flaring , i.e. @xmath209 , the irradiation flux on the accretion disc is at most @xcite @xmath210 where @xmath211 is the disc albedo ( for a more complete treatment taking into account shadowing by the disc see @xcite ) . the local flux generated by the viscous heating in the accretion disc is on the other hand @xmath212 we thus find that @xmath213 and since @xmath214 we find that irradiation dominates for @xmath215 if @xmath216 . consequently irradiation dominates over viscous heating for @xmath217 km . thus the outer part of the accretion disc is hotter than is predicted by our model , which can be observed as an increased uv flux @xcite . thus we have underestimated the pressure in the outer disc , and since we have used the @xcite prescription @xmath218 to estimate @xmath26 , we have also underestimated @xmath26 and the torque from the outer disc on the neutron star . however , this torque is small in comparison to the total torque on the neutron star since @xmath219 is a rapidly decreasing function of @xmath7 , and thus the effect of the irradiation of the outer disc does not significantly change our results . irradiation is unlikely to significantly affect the inner disc , but , as pointed out above , other effects , such as radial advection and the transition from a disc flow to a funnel flow , that are not considered in this paper become important in the inner disc . we will investigate these effects in the future . the accretion rates that have been observed in millisecond x - ray pulsars cover several orders of magnitude , and the highest accretion rates are sufficient that the innermost part of the accretion disc is dominated by radiation pressure and electron scattering . for these reasons we have followed the approach by @xcite and divided our disc model in three regions , an outer region dominated by gas pressures and kramer s opacity , a middle region dominated by gas pressure and electron scattering , and an inner region dominated by radiation pressure and electron scattering . the large spin changes that have been observed in some of the accreting millisecond x - ray pulsars are difficult to explain in the model by @xcite , but the inclusion of a disc dynamo produces a significant enhancement of the torque between the accretion disc and the neutron star . this increase is not only due to the coupling between the magnetic fields of the neutron star and the accretion disc as such , but also to the fact that the accretion disc is truncated at a larger radius thus increasing the lever arm between the edge of the disc and the neutron star . sbt thanks the department of physics at the university of gothenburg for their hospitality and support during a part of this project . sbt is supported in part by the swedish institute ( si ) under their guest scholarship program . this research has made use of nasa s astrophysics data system . we thank an anonymous referee for comments that have improved the quality of the paper . we have also benefitted from discussions with dr . jonathan ferreira .
Abstract:

Answer: millisecond x - ray pulsars have weak magnetic dipole moments of @xmath0tm@xmath1 compared to ordinary x - ray pulsars with dipole moments of @xmath2tm@xmath1 . for this reason a surrounding accretion disc can extend closer to the neutron star , and thus reach a higher temperature , at which the opacity is dominated by electron scattering and radiation pressure is strong . we compute the self - similar structure of such a geometrically thin axisymmetric accretion disc with an internal dynamo . such models produce significantly stronger torques on the neutron star than models without dynamos , and can explain the strong spin variations in some millisecond x - ray pulsars . [ firstpage ]    accretion , accretion discs  magnetic fields  mhd stars : magnetic fields  stars : neutron  x - rays : binaries  pulsars : general