nonthermal emission is commonly observed among the lowest - mass stars , traced by optical line ( e.g. , ca ii , h i sequence ) , x - ray , uv and radio emission . this emission can be both persistent ( quiescent ) and eruptive , with short - duration flares from m dwarfs occurring at a rate of roughly 3% @xcite . the incidence and strength of quiescent magnetic activity as traced by the h@xmath6 line reaches @xmath380% and @xmath7 @xmath8 @xmath04 , respectively , among nearby late - type m dwarfs @xcite , but both metrics decline precipitously for the cooler l and t dwarfs @xcite and sources far from the galactic plane @xcite .
similar declines are seen in x - ray emission , but surprisingly not at radio frequencies . assuming that this nonthermal emission arises from magnetic interaction , the decline with vertical scale height among m dwarfs is likely an age effect , as angular momentum loss over time results in weakened magnetic dynamos .
however , spin - down timescales exceed 5  gyr beyond spectral type m5 @xcite , and many late - type dwarfs are found to be rapid rotators ( p @xmath8 2 - 10  hr ; @xcite ) exhibiting kilogauss magnetic fields at their photospheres @xcite .
hence , the decline in magnetic emission with spectral type must arise from a different effect .
the favored cause is the decoupling of cooler , increasingly neutral photospheres from magnetic structures , which reduces magnetic stresses and the frequency of magnetic reconnection above the ( sub)stellar surface ( e.g. , ) .
deeper reconnection events may continue to power strong flaring bursts observed in a handful of weakly active or inactive late - m and l dwarfs ( e.g. , @xcite ) .
contrary to these trends , a very rare set of low - temperature `` hyperactive '' dwarfs exhibit unusually prodigious and persistent nonthermal emission .
one of the first examples of these to be identified was 2mass  j13153094@xmath02649513 ( hereafter 2mass  j1315@xmath02649 ; @xcite ) , a high proper motion l5e dwarf which has exhibited sustained but variable h@xmath6 emission on no fewer than seven epochs spanning nearly a decade . with @xmath7 @xmath8 @xmath04 ,
2mass  j1315@xmath02649 is as active as a mid - type m dwarf , but is 12 orders of magnitude brighter in h@xmath6 than comparably - classified l dwarfs @xcite .
in addition , h@xmath9 and na  i d lines have also been observed in emission .
since the photosphere of 2mass  j1315@xmath02649 is cool and likely to be highly neutral , the origin of its unexpected emission remains a mystery .
neither the kinematics nor spectral characteristics of 2mass  j1315@xmath02649 indicate youth , and the absence of mid - infrared excess argues against accretion from a protoplanetary or debris disk @xcite .
other hyperactive dwarfs , such as the l1e dwarf 2mass  j10224821 + 5825453 ( @xmath7 @xmath8 @xmath03.5 to @xmath02.7 ; @xcite ) and the t6.5e dwarf 2mass j1237392 + 652615 ( @xmath7 @xmath8 @xmath04.6 to @xmath04.2 ; @xcite ; hereafter 2mass  j1237 + 6526 ) also lack evidence of disk accretion and do not appear to be particularly young .
alternative mechanisms , such as acoustic heating , unusually strong magnetic fields , and roche lobe overflow from a substellar companion have been proposed , but none of these scenarios have been validated .    in this article , we report new observations of 2mass  j1315@xmath02649 that reveal both continued nonthermal optical line emission in several neutral atomic species , and the presence of a t dwarf companion at a projected separation of 7  au . in section  2 we describe our combined - light optical and near - infrared spectroscopic observations , the latter of which yields preliminary evidence for a brown dwarf companion . in section  3
we describe adaptive optics ( ao ) imaging and spectroscopic observations that confirm the presence of the companion and allow determination of its separation and classification . in section  4
we analyze the observed and inferred physical properties of the components , the latter based on comparison to evolutionary and atmospheric models .
we also perform a coevality test to examine the reliability of these models . in section  5
we discuss how the properties of 2mass  j1315@xmath02649 argue for a magnetic origin to its persistent emission , powered by a strong magnetic field retained by a relatively old and massive cool dwarf .
we summarize our results in section  6 .
moderate - resolution optical spectra of 2mass  j1315@xmath02649 were obtained on 2011 march 26 ( ut ) with the magellan echellette ( mage ; @xcite ) , mounted on the 6.5 m landon clay telescope at las campanas observatory .
conditions during the observations were clear with 0@xmath106 seeing .
two exposures totaling 3000  s were obtained at an average airmass of 1.003 using the 0@xmath107 slit aligned with the parallactic angle ; this setup provided 320010050   spectroscopy at a resolution @xmath11 @xmath8 4000 .
we also observed the spectrophotometric flux standard eg  274 @xcite on the same night for flux calibration .
thar lamps were obtained after each source observation for wavelength calibration , and internal quartz and dome flat field lamps were obtained during the night for pixel response calibration .
data were reduced using the mase reduction pipeline @xcite , following standard procedures for order tracing , flat field correction , wavelength calibration ( including heliocentric correction ) , optimal source extraction , order stitching , and flux calibration .
the red portion of the reduced spectrum is shown in figure  [ fig_optspec ] , flux - calibrated to an apparent @xmath12-band magnitude of 20.16@xmath10.14 as estimated from two micron all sky survey ( 2mass ) photometry and a mean @xmath13 = 4.97@xmath10.13 color for l5 dwarfs @xcite .
we confirm the characteristic mid - l dwarf features identified in previous studies , including strong feh and crh bands ; weak tio absorption ( relative to late - m and early - l dwarfs ) ; line absorption from na  i , rb  i and cs  i ; and the strongly pressure - broadened 7700   k  i doublet .
the overall spectral shape is well - matched to the l5 dwarf 2mass j15074769@xmath01627386 @xcite , consistent with previously reported classifications @xcite .
we confirm the absence of 6710   li  i absorption as reported by @xcite to an equivalent width ( ew ) limit of 0.5   .
this is well below measured ews for equivalently - classified l dwarfs ( e.g. , @xcite ) .
we also see no evidence of any peculiar spectral features associated with low surface gravities , such as enhanced vo absorption or weakened alkali lines @xcite
. line center measurements of the alkali lines indicate a heliocentric radial velocity of @xmath07@xmath19  km  s@xmath14 .    lccc + li  i ( 6710   ) & @xmath150.5 & & + rb  i ( 7802   ) & 4.9@xmath10.8 & & + rb  i ( 7950   ) & 5.4@xmath10.6 & & + na  i ( 8186   ) & 1.2@xmath10.3 & & + na  i ( 8197   ) & 2.5@xmath10.3 & & + cs  i ( 8523   ) & 4.2@xmath10.2 & & + cs  i ( 8946   ) & 2.9@xmath10.7 & & +   + h@xmath16 ( 4342   ) & & 4.0@xmath11.8 & -5.12@xmath10.20 + h@xmath9
( 4863   ) & & 13.1@xmath11.5 & -4.61@xmath10.08 + h@xmath6 ( 6565   ) & @xmath058@xmath14 & 34.9@xmath11.2 & -4.18@xmath10.06 + na  i ( 5892   ) & & 1.6@xmath10.2 & -5.50@xmath10.08 + na  i ( 5898   ) & & 2.1@xmath10.3 & -5.40@xmath10.08 + k  i ( 7667   ) & & 0.6@xmath10.1 & 5.92@xmath10.10 + k  i ( 7701   ) & & 0.9@xmath10.2 & 5.79@xmath10.13 +    the most striking feature in the optical spectrum of 2mass  j1315@xmath02649 is its pronounced h@xmath6 emission .
we measure an equivalent width of @xmath058@xmath14   in the mage data , in the middle of prior measurements that span @xmath024 to @xmath0160   .
we also detect h@xmath16 ( 4342   ) and h@xmath9 ( 4863   ) in emission , and confirm the presence of na  i emission as reported in .
in addition , we detect weak emission in the cores of the 7700   k
i doublets , but no emission from the 5877   he  i d3 ( ew @xmath17 18   ) or 85008660   ca  ii triplet lines ( ew @xmath17 0.3   ) . line profiles of detected emission features ( figure  [ fig_emission ] ) show no appreciable broadening at the 75  km  s@xmath14 velocity resolution of mage , nor do we detect any significant velocity shift between emission and absorption lines ( @xmath18@xmath19 = 4@xmath116  km  s@xmath14 ) .
as most of these lines are superimposed on an undetected continuum , we report in table  [ tab_emission ] line fluxes and relative line - to - bolometric luminosities , that latter computed using the bolometric correction / spectral type relations of @xcite .
our measurement of @xmath7 = @xmath04.18@xmath10.06 is similar to the first detection made by @xcite , and is again 12 orders of magnitude greater than measured quiescent or flaring fluxes for equivalently classified l dwarfs .
the balmer decrement @xmath20 = 2.7@xmath10.3 is roughly comparable to the mean values for non - flaring m dwarfs @xcite .
h@xmath6 and h@xmath9 emission contribute 61% and 23% of the total measured line flux of @xmath21 @xmath8 @xmath04 .
low resolution near - infrared spectra of 2mass  j1315@xmath02649 were obtained on 2009 june 30 ( ut ) with the 3 m nasa infrared telescope facility spex spectrograph @xcite .
observing conditions were clear and dry with 0@xmath107 seeing at @xmath22-band .
we used the prism - dispersed mode of spex with a 0@xmath105 slit ( aligned to the parallactic angle ) to obtain a continuous 0.72.5  @xmath23 spectrum with resolution @xmath11 @xmath24 .
a total of 8 exposures of 120  s each were obtained in two abba dither pairs , nodding along the slit , at an average airmass of 1.50 .
we also observed the a0  v star hd  125438 ( @xmath25 = 7.10 ) for flux calibration and telluric absorption correction , as well as internal flat field and argon arc lamps for pixel response and wavelength calibration .
data were reduced with the idl spextool package , version 3.4 @xcite , using standard settings ; see @xcite for details .
the reduced spectrum of 2mass  j1315@xmath02649 is shown in figure  [ fig_nirspec ] , compared to the l5 dwarf 2mass  j21373742 + 0808463 ( hereafter 2mass  j2137 + 0808 ; @xcite ) .
its near - infrared spectral morphology is consistent with its l5 optical type , with strong h@xmath4o and co absorption bands , feh absorption at 1.0 and 1.6  @xmath23 , and ( unresolved ) alkali line absorption in the 1.11.3  @xmath23 region .
the overall spectral shape is again consistent with a normal l5 field dwarf , with no evidence of an unusual surface gravity , metallicity or cloud content ( e.g. , @xcite )
.    there is , however , a subtle `` notch '' feature present at 1.62  @xmath23 that differs from the spectrum of 2mass  j2137 + 0808 ( see inset box of figure  [ fig_nirspec ] ) .
this feature has previously been noted in the combined - light spectra of l dwarf plus t dwarf binaries , arising from the overlap of feh absorption in the primary and ch@xmath26 absorption in the secondary ( e.g. @xcite ) . to characterize this feature
, we performed a spectral template fitting analysis similar to that described in @xcite , using 295 l2t8 spectral templates from the spex prism spectral libraries .
fluxes of these templates were scaled to the @xmath27/spectral type relation of @xcite , and all spectra were interpolated onto a common wavelength scale . restricting potential secondaries to have t spectral types , a total of 17889 binary templates were constructed and compared to the spectrum of 2mass  j1315@xmath02649 over the wavelength ranges 0.951.35
@xmath23 , 1.451.80  @xmath23 and 2.002.35  @xmath23 using the @xmath28 statistic .     +    figure  [ fig_binfit ] shows the four best - fitting binary templates from these comparisons , a combination of 2mass  j2137 + 0808 and either early- or late - type t dwarfs .
the addition of a t - type companion fills in the `` missing '' flux at 1.58  @xmath23 , creating the distinct notch feature at 1.62  @xmath23 ; and at 1.27  @xmath23 , producing a somewhat sharper @xmath29-band flux peak .
importantly , all of the binary templates shown in figure  [ fig_binfit ] provide statistically superior fits to the spectrum of 2mass  j1315@xmath02649 as compared to 2mass  j2137 + 0808 alone , based on the f - test statistic ( eqns .
26 in @xcite ) .
however , we can not precisely constrain the properties of the secondary from this analysis ; the average spectral type of all of the template fits weighted by the f - test probability distribution function ( f - pdf ) is t3@xmath14 .
to more accurately characterize this putative companion , high - resolution , near - infrared images of 2mass  j1315@xmath02649 were obtained with the 10 m keck ii laser guide star adaptive optics system ( lgsao ; @xcite ) and facility nirc2 near - infrared camera .
observations were conducted on two runs , 2010 march 24 and 2010 may 13 ( ut ) , both with clear skies and fair seeing ( @xmath171@xmath30 and 0@xmath105 , respectively ) .
we used the narrow camera with image scale 9.963@xmath10.011  mas  pixel@xmath14 @xcite covering a 10@xmath102@xmath3110@xmath102 field of view .
images were obtained through the @xmath29 , @xmath22 , and @xmath32 filters , using a three - point dither pattern that avoided the noisy , lower left quadrant of the focal plane array .
exposure times ranged from 30  s with 8  coadds to 120  s with 2 coadds per pointing position , with total integrations of 360  s to 720  s in a given filter .
the sodium lgs provided the wave front reference source for ao correction , while tip - tilt aberrations and quasi - static changes were measured by monitoring the @xmath33 = 12.8 field star usno - b1.0 0631 - 0348160 @xcite located @xmath34 = 38@xmath106 from 2mass  j1315@xmath02649 .
images were reduced using custom idl scripts , as described in @xcite .
@xmath311@xmath30 nirc2 @xmath29 , @xmath22 and @xmath32 images of the 2mass  j1315@xmath02649ab system .
images are linearly scaled to optimize visibility of the faint secondary ( southeast of the primary ) , and are oriented with north up and east to the left .
[ fig_nirc2],title="fig:",scaledwidth=30.0% ] @xmath311@xmath30 nirc2 @xmath29 , @xmath22 and @xmath32 images of the 2mass  j1315@xmath02649ab system .
images are linearly scaled to optimize visibility of the faint secondary ( southeast of the primary ) , and are oriented with north up and east to the left .
[ fig_nirc2],title="fig:",scaledwidth=30.0% ] @xmath311@xmath30 nirc2 @xmath29 , @xmath22 and @xmath32 images of the 2mass  j1315@xmath02649ab system .
images are linearly scaled to optimize visibility of the faint secondary ( southeast of the primary ) , and are oriented with north up and east to the left .
[ fig_nirc2],title="fig:",scaledwidth=30.0% ]    figure  [ fig_nirc2 ] displays the reduced nirc2 images from our may observations .
a faint source is clearly present southeast of 2mass  j1315@xmath02649 , and was visible during both imaging epochs .
relative astrometry for each epoch was measured on the individual frames using a centroiding algorithm , and these values were then averaged and multiplied by the pixel scale ( uncertainties include the standard deviations of the position measurements and 0.1% pixel scale uncertainty ) .
relative photometry was performed on the coadded mosaics through aperture photometry . as the point spread function ( psf ) of the brighter component contributes significant flux ( @xmath825% ) to the brightness of the fainter object , photometry for the latter
was extracted from a primary psf - subtracted image , constructed by rotating the frame 180@xmath35 about the centroid of the primary and differencing .
systematic errors in the primary subtraction were estimated by offsetting the rotation axis over a 5@xmath315-pixel grid centered on the original centroid , subtracting , and measuring aperture photometry on the secondary .
the final photometric values for each band and epoch were taken as the means and standard deviations of these 25 measurements .
lc + @xmath36 ( @xmath30 ) & 185@xmath16 + @xmath37 ( @xmath30 ) & @xmath0277@xmath13 + @xmath34 ( @xmath30 ) & 333@xmath17 + @xmath38 ( @xmath39 ) & 143.1@xmath11.1 + @xmath40 ( mag ) & 3.12@xmath10.05 + @xmath41 ( mag ) & 4.29@xmath10.14 + @xmath2 ( mag ) & 4.91@xmath10.18 +   + @xmath36 ( @xmath30 ) & 187@xmath12 + @xmath37 ( @xmath30 ) & @xmath0286@xmath17 + @xmath34 ( @xmath30 ) & 342@xmath17 + @xmath38 ( @xmath39 ) & 143.9@xmath10.5 + @xmath40 ( mag ) & 3.00@xmath10.03 + @xmath41 ( mag ) & 4.59@xmath10.07 + @xmath2 ( mag ) & 5.17@xmath10.12 +   + @xmath36 ( @xmath30 ) & 180@xmath15 + @xmath37 ( @xmath30 ) & @xmath0286@xmath16 + @xmath34 ( @xmath30 ) & 338@xmath14 + @xmath38 ( @xmath39 ) & 147.8@xmath11.2 +    results are listed in table  [ tab_astrometry ] .
separations in right ascension and declination are consistent between both epochs , and ( in conjunction with the osiris observations described below ) yield a mean separation of 336@xmath16  mas at position angle of @xmath42 , measured from primary to secondary .
the mean relative magnitudes are also statistically consistent between epochs ( to within 2@xmath43 ) , and indicate that the companion is both considerably fainter ( @xmath40 = 3.03@xmath10.03 , @xmath2 = 5.09@xmath10.10 ) and bluer in the near - infrared . to our knowledge , this is the second most extreme near - infrared flux ratio measured for an ultracool dwarf binary to date2649 is exceeded at @xmath29-band only by the young tw hydrae binary 2massw  j1207334 - 393254ab ( @xmath40 = 7.0@xmath10.2 , @xmath44 = 4.98@xmath10.14 ; ) ; and at @xmath45-band only by the old , widely - separated binary sdss j141624.08 + 134826.7ab ( @xmath40 = 4.31@xmath10.02 , @xmath44 = 6.85@xmath10.17 ; ) . ] .
we examine the physical association of the companion in section  4.3 ; hereafter , we refer to the two sources as 2mass  j1315@xmath02649a and b.      resolved @xmath22-band spectroscopy of 2mass  j1315@xmath02649ab was obtained using the keck ii oh - suppressing infrared integral field spectrograph ( osiris ; @xcite ) and lgsao in mostly clear conditions on 2010 may 19 ( ut ) .
we used the 35  mas - scale camera and h@xmath46 filter , providing 1.471.80  @xmath23 spectroscopy at an average resolution of 3800 and dispersion of 2.1    pixel@xmath14 over a 0@xmath1056@xmath312@xmath1024 field of view .
the instrument rotator was set at a position angle of @xmath045@xmath35 to accommodate both components in the rectangular field of view .
six exposures of 600  s each were obtained at an average airmass of 1.52 using a linear dither pattern with steps of 0@xmath104 along the long axis , and tip - tilt correction for lgsao operation was again provided by usno - b1.0 0631 - 0348160 .
these observations were followed by a 600  s exposure of a blank sky frame .
we also obtained three dithered 20  s exposures of the a0  v star hd  107120 ( @xmath25 = 9.90 ) in natural guide star ( ngs ) ao mode at an airmass of 1.56 for telluric absorption correction and flux calibration .
data were reduced with the osiris data reduction pipeline @xcite , version 2.3 .
we first subtracted the 600  s sky frame from each of the 2mass  j1315@xmath02649 images , and a median - combined dark frame from the calibrator images .
we then used the pipeline to adjust bias levels , remove detector artifacts and cosmic rays , extract and wavelength - calibrate the position - dependent spectra ( using the most current rectification files as of february 2011 ) , assemble 3d data cubes and correct for dispersion .
spectra for the 2mass  j1315@xmath02649 primary and hd  107120 were extracted directly from the data cube by aperture photometry in each image plane , using a 3-pixel ( 105  mas ) aperture and 10 - 20  pixel ( 350700  mas ) sky annulus . for the faint companion , light contamination from the primary was a concern , so we first performed a partial subtraction of the primary s radial brightness profile .
we sampled the profile over two position angle ranges 2040@xmath35 away from the separation axis , generating a mean profile as a function of wavelength and separation .
we then subtracted this profile @xmath125@xmath35 about the separation axis .
figure  [ fig_osiris_image ] displays mosaics of both the original and subtracted images , illustrating the reduced background achieved around the companion .
we measured aperture photometry for this component in each of the subtracted image planes , using a more restricted 1.5  pixel ( 52  mas ) aperture and 35  pixel ( 105 - 175  mas ) sky annulus .
the individual spectra for all three sources were scaled and combined using the _ xcombspec _ routine in spextool @xcite .
flux calibration and telluric correction of the 2mass  j1315@xmath02649ab spectra were performed using the _ xtellcor_general _ routine in spextool , assuming a 20  nm gaussian kernel for the a0v h  i lines @xcite .
figure  [ fig_osiris_spec ] displays the reduced spectra of the two components of 2mass  j1315@xmath02649 .
the spectrum of 2mass  j1315@xmath02649a has exceptionally high signal - to - noise ( s / n @xmath8 200 ) , and is similar to that of the combined - light spex spectrum , with strong h@xmath4o absorption wings shortward of 1.55  @xmath23 and longward of 1.7  @xmath23 , and weak feh absorption in the 1.571.64  @xmath23 region , all indicative of a mid - type l dwarf .
the notch feature , however , is no longer present .
the spectrum of 2mass  j1315@xmath02649b ( s / n @xmath8 25 at 1.6  @xmath23 ) is unambiguously that of a late - type t dwarf , with strong ch@xmath26 absorption at 1.6  @xmath23 .
component spectral types were determined by comparing the resolved osiris spectra to the spex prism spectral libraries templates , restricting the template sample to optically - classified l dwarfs and near - infrared - classified t dwarfs . following a @xmath28-fitting procedure similar to that described above over the 1.51.8  @xmath23 region
, we identified the l5 2mass  j21373742 + 0808463 @xcite and the t7 2mass  j07271824 + 1710012 @xcite as the best - fitting templates to 2mass  j1315@xmath02649a and b , respectively ( figure  [ fig_osiris_spec ] ) .
an f - test weighted average of all templates indicates mean classifications of l3.5@xmath12.5 and t7@xmath10.6 for the components .
the large uncertainty for the former is largely due to the broad diversity of near - infrared spectra exhibited by l dwarfs at a given optical spectral type , arising from variations in surface gravity , metallicity and cloud properties @xcite .
we therefore adopt the optical l5 classification for this component and the near - infrared t7 classification for the secondary .
component brightnesses on the mko system were determined by converting the combined - light 2mass @xmath47 and relative nirc2 @xmath47 magnitudes to mko @xmath48 .
these conversions were computed directly from the spex prism spectra of 2mass  j1315@xmath02649ab ( 2mass @xmath49 mko ) and the best - fit templates in figure  [ fig_osiris_spec ] ( nirc2 @xmath49 mko ) using the appropriate filter profiles and a kurucz model spectrum of vega ( see @xcite ) .
the resulting combined light and component magnitudes are listed in tables  [ tab_system ] and  [ tab_components ] , respectively .
both components appear to have relatively normal near - infrared colors for their spectral types @xcite .
lcl spectral type & l5 & 1,2 + est .
distance ( pc ) & 19@xmath13 & 1 + mko @xmath29 ( mag ) & 15.07@xmath10.05 & 1,3 + mko @xmath22 ( mag ) & 14.12@xmath10.04 & 1,3 + mko @xmath45 ( mag ) & 13.45@xmath10.04 & 1,3 + mko @xmath50 ( mag ) & 1.63@xmath10.07 & 1,3 + @xmath51 ( mas  yr@xmath14 ) & @xmath0682@xmath113 & 4 + @xmath52 ( mas  yr@xmath14 ) & @xmath0282@xmath114 & 4 + @xmath53 ( km  s@xmath14 ) & 65@xmath110 & 1,4 + @xmath19 ( km  s@xmath14 ) & @xmath07@xmath19 & 1 + u ( km  s@xmath14 ) & @xmath038@xmath18 & 1 + v ( km  s@xmath14 ) & @xmath041@xmath19 & 1 + w ( km  s@xmath14 ) & @xmath013@xmath15 & 1 + @xmath34 ( mas ) & 336@xmath13 & 1 + @xmath34 ( au ) & 6.6@xmath10.9 & 1 + @xmath38 ( @xmath35 ) & 336@xmath13 & 1 + @xmath40 ( mag ) & 3.03@xmath10.03 & 1 + @xmath41 ( mag ) & 4.53@xmath10.06 & 1 + @xmath2 ( mag ) & 5.09@xmath10.10 & 1 + age ( gyr ) & @xmath30.81.0 & 1 + est .  orbit period ( yr ) & 4560 ( 1595 ) & 1 +    lccc + nir spt & l3.5@xmath12.5 & t7@xmath10.6 & + mko @xmath29 ( mag ) & 15.14@xmath10.05 & 18.20@xmath10.06 & + mko @xmath22 ( mag ) & 14.14@xmath10.03 & 18.66@xmath10.07 & + mko @xmath45 ( mag ) & 13.45@xmath10.04 & 18.79@xmath10.11 & + mko @xmath50 ( mag ) & 1.69@xmath10.06 & @xmath00.59@xmath10.12 & + distance ( pc ) & 18@xmath14 & 36@xmath19 & 17@xmath19 + @xmath54 ( dex ) & -4.19@xmath10.16 & @xmath05.86@xmath10.16 & +   + ( k ) & 1760@xmath170 & 790@xmath170 & + @xmath55 ( cm  s@xmath56 ) & @xmath35.2 & 5.0@xmath10.5 & +   + mass ( m@xmath57 ) & 60@xmath16 & 16@xmath13 & 0.26@xmath10.03 + t@xmath58 ( k ) & 1720@xmath1150 & 630@xmath160 & + @xmath55 ( cm  s@xmath56 ) & 5.26@xmath10.04 & 4.57@xmath10.09 & + radius ( r@xmath57 ) & 0.90@xmath10.01 & 1.01@xmath10.02 & +   + mass ( m@xmath57 ) & 74@xmath12 & 28@xmath14 & 0.38@xmath10.05 + t@xmath58 ( k ) & 1770@xmath1140 & 670@xmath170 & + @xmath55 ( cm  s@xmath56 ) & 5.41@xmath10.01 & 4.93@xmath10.08 & + radius ( r@xmath57 ) & 0.84@xmath10.02 & 0.91@xmath10.02 & +   + mass ( m@xmath57 ) & 76@xmath12 & 37@xmath15 & 0.48@xmath10.06 + t@xmath58 ( k ) & 1790@xmath1140 & 690@xmath170 & + @xmath55 ( cm  s@xmath56 ) & 5.44@xmath10.02 & 5.10@xmath10.08 & + radius ( r@xmath57 ) & 0.83@xmath10.02 & 0.85@xmath10.02 & +   + mass ( m@xmath57 ) & 77@xmath11 & 48@xmath15 & 0.62@xmath10.06 + t@xmath58 ( k ) & 1790@xmath1140 & 720@xmath170 & + @xmath55 ( cm  s@xmath56 ) & 5.44@xmath10.02 & 5.28@xmath10.07 & + radius ( r@xmath57 ) & 0.83@xmath10.02 & 0.79@xmath10.02 & +    component distances were computed using the mko @xmath48 absolute magnitude / spectral type relations of @xcite ; we considered both the `` bright '' and `` faint '' relations .
propagating uncertainties in component spectral types and photometry , and scatter in the relations , through monte carlo simulation yields consistent distances for each component in all three bands and both relations , with mean values of 18@xmath14  pc and 35@xmath19  pc for 2mass  j1315@xmath02649a and b , respectively .
these distances are formally consistent with each other , differing at the 1.8@xmath43 level ; the latter has a larger uncertainty due to the larger photometric error for this component ( table  [ tab_components ] ) .
the error - weighted mean distance of 20@xmath13  pc matches the 21.7  pc estimate of @xcite . at this distance ,
the projected separation of the components is 6.6@xmath10.9  au , right at the peak of the separation distribution of resolved very low - mass field binaries @xcite .
component luminosities were computed from the individual mko @xmath48 magnitudes using bolometric correction ( bc ) relations as a function of spectral type , as quantified in @xcite .
apparent bolometric magnitudes were converted to absolute bolometric magnitudes by adopting a common distance of 20@xmath13  pc , and luminosities calculated assuming @xmath59 = 4.74 .
luminosities computed in each of the @xmath48 bands were again mutually consistent , resulting in @xmath54 = @xmath04.19@xmath10.16 and @xmath05.86@xmath10.16 for 2mass  j1315@xmath02649a and b , respectively .
the luminosity for 2mass  j1315@xmath02649a is similar to other l4.5l5.5 field dwarfs as compiled by @xcite and @xcite , while 2mass  j1315@xmath02649b is somewhat underluminous for its spectral type .
the similar distances and relative small projected separation of 2mass  j1315@xmath02649a and b indicates that these sources are cospatial ; we also find that their space motions are aligned .
the three epochs of nirc2 and osiris relative astrometry are consistent with each other in both right ascension and declination , and despite the short period between the observations the high proper motion of 2mass  j1315@xmath02649 ( @xmath51 = @xmath0682@xmath113  mas  yr@xmath14 , @xmath52 = @xmath0282@xmath114  mas  yr@xmath14 ; @xcite ) allows us to rule out either component as a ( non - moving ) background star at the 9@xmath43 level . furthermore , the magnitudes and position angles of the component proper motions are identical to within @xmath147  mas  yr@xmath14 ( 4.4  km  s@xmath14 at 20  pc ) and @xmath14@xmath35 .
the radial motions of the two components are also equivalent .
cross - correlation of the osiris spectra with zero - velocity spectral model templates from ( * ? ? ?
* see section  4.4 ) yield identical velocities to within the uncertainties ( @xmath18@xmath19 = 0@xmath19  km  s@xmath14 ) .    with a projected tangential velocity of @xmath53 = 70@xmath118  km  s@xmath14 , and adopting the radial velocity from the combined - light optical spectrum above ,
we find @xmath60 space velocities in the local standard of rest = 10  km  s@xmath14 , @xmath25 = 5.25  km  s@xmath14 and @xmath61 = 7.17  km  s@xmath14 in the lsr @xcite , where the directions of @xmath60 follow a right - handed coordinate system . ] ( lsr ) of @xmath62 = @xmath038@xmath18  km  s@xmath14 , @xmath25 = @xmath041@xmath19  km  s@xmath14 and @xmath61 = @xmath013@xmath15  km  s@xmath14 .
these values lie outside the 1@xmath43 velocity spheroid for the `` cold '' population of nearby l dwarfs reported by @xcite , but are within the velocity spheroid of local thick disk stars ( e.g. , ) .
both the kinematics and spectral properties of 2mass  j1315@xmath02649 are therefore consistent with an older dwarf system in the galactic disk population .      assuming that 2mass  j1315@xmath02649ab comprises a coeval system , insight into the physical properties of its components can be obtained by joint comparison to atmospheric and evolutionary models .
the spectral and kinematics analysis above suggests that 2mass  j1315@xmath02649ab is an older system ; a more quantitative constraint comes from the absence of li  i absorption in the combined light optical spectrum , which is dominated by the l5 primary .
the absence of this line sets a minimum mass of @xmath630.06  m@xmath64 for this component @xcite .
combined with its luminosity , the evolutionary models of @xcite , and @xcite indicate a minimum age for the primary ranging from 0.8 to 1.0  gyr .
this implies a minimum secondary mass of 0.013  m@xmath64 , around the deuterium - burning minimum mass limit .
+    figure  [ fig_evol ] displays the regions in t@xmath58/@xmath55 space occupied by the components as constrained by their luminosities , the minimum mass of 2mass  j1315@xmath02649a , and the evolutionary models listed above .
table  [ tab_components ] details specific physical properties for select ages based on the evolutionary models of . with no empirical upper limit on the age of this system ,
the mass of 2mass  j1315@xmath02649a could be above the hydrogen - burning mass limit ( m @xmath3 0.072  m@xmath64 for @xmath65 @xmath3 2  gyr ) , while the secondary must be substellar at any age .
the estimated mass ratio of the system , @xmath66 m@xmath4/m@xmath5 , is one of the smallest inferred for a very low - mass binary : @xmath630.3 ( @xmath630.6 ) for an age of 1  gyr ( 10  gyr ) .
this makes 2mass  j1315@xmath02649ab a unique system , given that @xmath6390% of resolved brown dwarf binaries identified to date have @xmath67 @xcite .
indeed , if low - mass binaries prefer to be in higher mass ratio systems , these values further support the hypothesis that 2mass  j1315@xmath02649 is quite old .
t@xmath58 constraints on the components15701930  k for the primary and 570790  k for the secondary ( 1@xmath43 ranges)are again consistent with comparably classified field dwarfs @xcite , with 2mass  j1315@xmath02649b being somewhat on the cool side for its spectral type . the lower mass limit for 2mass  j1315@xmath02649a tightly constrains its surface gravity to @xmath55 = 5.225.46  cm  s@xmath56 , while 2mass  j1315@xmath02649b has a broader constraint of 4.465.35  cm  s@xmath56 .
an independent assessment of the component atmospheric parameters was made by fitting the osiris spectra to the bt - settl models of @xcite .
these models are based on the phoenix code @xcite , and reflect an update to the original settl models of @xcite with a microturbulence velocity field determined from 2d hydrodynamic models and updated solar abundances from .
we followed the same fitting procedure described in ( * ? ? ?
* see also @xcite and @xcite ) , using a set of solar - metallicity models sampling t@xmath58 = 6002500  k ( 100  k steps ) and @xmath55 = 4.05.5  cm  s@xmath56 ( 0.5  cm  s@xmath56 steps ) .
model surface fluxes ( in @xmath68 units ) and the osiris spectra were smoothed to a common resolution of @xmath11 = 3500 using a gaussian kernel , and interpolated onto a common wavelength grid .
the data were then scaled to the appropriate @xmath22-band apparent magnitude .
data and models were compared over the 1.51.75  @xmath23 region using a @xmath28 statistic , with the degrees of freedom equal to the number of resolution elements sampled .
an optimal scaling factor @xmath69 was computed for each fit to minimize @xmath28 , where @xmath33 is the radius of the brown dwarf and @xmath70 its distance from the sun @xcite .
we further constrained our model selection by requiring that the model - inferred distance be within 2@xmath43 of the estimated spectrophotometric distance of the system ( 1326  pc ) .
means and uncertainties in the atmospheric parameters were determined using the f - pdf as a weighting factor , as above ; we also propagated sampling uncertainties of 50  k and 0.25  dex for t@xmath58 and @xmath55 , respectively .
note that these uncertainties quantify _ experimental _ errors ; they do not account for _ systematic _ errors that arise from the fidelity of the model fits , as discussed below .
figure  [ fig_osirisfit ] displays the best - fitting models for the osiris spectra . for wavelengths shortward of 1.55  @xmath23 ,
the models provide reasonably good fits to the forest of h@xmath4o lines present in both component spectra .
however , at longer wavelengths we see deviations in the primary arising from missing feh opacity and a premature downturn in fluxes longward of 1.65  @xmath23 .
the latter is symptomatic of overly blue spectral energy distributions across the near - infrared range exhibited by the models at these temperatures .
there are also deviations in model fits to the secondary near the 1.58  @xmath23 peak and within the 1.61.75  @xmath23 ch@xmath26 absorption system .
note that , quantitatively , the fit to the secondary is better ( @xmath71 = 1.52 ) than that to the primary ( @xmath71 = 38.4 ) , but this mainly stems from the high s / n data for the latter ; neither fit reproduces the observed spectrum with great fidelity .
nevertheless , visual inspection confirms that these are the best fits among the model sample , and the inferred mean parameters
t@xmath58 = 1760@xmath170  k and @xmath55 @xmath3 5.2  cm  s@xmath56 for the primary and t@xmath58 = 790@xmath170  k and @xmath55 = 5.2@xmath10.4  cm  s@xmath56 for the secondary  are roughly in line with estimates from the evolutionary model parameters above .      with independent determinations of luminosity and t@xmath58 for both components
, we can now examine whether the evolutionary and atmospheric models are consistent with each other assuming the system is coeval ; this is the so - called coevality test ( e.g. , @xcite ) .
prior studies of brown dwarf binaries have produced mixed results with respect to this test , with a few low - temperature systems showing evidence of @xmath6350100  k offsets ( both high and low ) between evolutionary and atmospheric models @xcite .
however , spectroscopic t@xmath58s have generally been taken from estimates of comparably classified field dwarfs , rather than spectroscopic fits to the binary components themselves .
such temperature estimation by proxy could result in systematic biases .
the single exception is the t1+t6 binary @xmath72 indi bc , for which resolved optical and near - infrared spectroscopy have enabled determinations of individual component luminosities and t@xmath58s , and comparison of these components on the hr diagram indicate that they are coeval with each other and with their stellar primary @xcite .
+    for 2mass  j1315@xmath02649ab , we find good agreement between atmospheric and evolutionary models for the primary but not for the secondary .
figure  [ fig_evol ] compares the t@xmath58 and @xmath55 constraints from our spectral models fits to those from the evolutionary model comparisons .
the atmospheric t@xmath58 and @xmath55 constraints for the primary overlap reasonably well for all four evolutionary model sets , although there is less agreement for the cloudy models of @xcite .
moreover , these values are consistent with the @xmath630.8 - 1.0  gyr minimum age of the system based on the absence of li  i in the optical spectrum . for the secondary , overlap in t@xmath58 and @xmath55 regions
is not as good , with essentially no overlap for the @xcite cloudy models .
this discrepancy can also be seen in figure  [ fig_hr ] , which compares the hr diagrams for all four evolutionary models to the luminosities and spectral model fit t@xmath58s for the 2mass  j1315@xmath02649 components .
both sit at the lower envelope of the evolutionary tracks , consistent with older ages ; however , 2mass  j1315@xmath02649b falls off the @xcite cloudy tracks entirely .
the sense of the deviations between the model comparisons of 2mass  j1315@xmath02649b is that the spectral model fit t@xmath58s are systematically higher than the evolutionary model t@xmath58s . as we do not have any other independent constraints on the system ( e.g. , age , metallicity or component masses )
, we can not _ a priori _ determine whether this mismatch is specifically attributable to errors in the spectral or evolutionary models .
however , based on the fits shown in figure  [ fig_osirisfit ] , we suspect the former given the poor match between the bt - settl models and spectral data around the 1.6  @xmath23 ch@xmath26 band .
this feature causes persistent problems in spectral model ftis due to incomplete ch@xmath26 opacities at t dwarf temperatures @xcite .
a decrease of just 100200  k in the derived secondary t@xmath58 would bring both components in precise alignment with evolutionary models in both t@xmath58/@xmath55 and hr diagrams , a shift previously suggested in prior low - mass binary analyses ( although not necessarily in the same direction ; @xcite ) .
it is also notable that spectral model fits for 2mass  j1315@xmath02649a are even worse than those for 2mass  j1315@xmath02649b .
it may be that the alignment of evolutionary and spectral model parameters for this component is merely an example of `` chance agreement '' . in any case , without more accurate fitting of l and t dwarf spectral data , such coevality tests are fundamentally inconclusive about the underlying accuracies of model - derived parameters .
we emphasize that the discrepancies noted here are only at the 1@xmath43 level , and should be verified through more precise constraints on the component luminosities and t@xmath58s .
these can best be accomplished through a parallax distance measurement of the system , as the distance uncertainty dominates the luminosity uncertainty .
in addition , resolved spectroscopy spanning the near - infrared ( and possibly optical ) range would provide a more robust test of the atmosphere models , and allow us to test different sets of models . unfortunately , mass measurements for this widely - separated system are probably not feasible in the near future .
a probability analysis of the possible orbits for this system.[foot_periodsim ] ] predicts likely periods of 4560  yr ( 15 - 95  yr at 90% confidence ) .
this rules out a `` fast '' astrometric orbit determination , and the relative orbital radial velocities ( @xmath730.5  km  s@xmath14 ) are comparable to current systematic uncertainties for isolated late - type dwarfs ( e.g. , @xcite ) . despite these challenges , the 2mass  j1315@xmath02649 system is an important benchmark for empirical tests of atmospheric and evolutionary models given its proximity , well - separated components , component types and relatively old age .
while the presence of a faint , substellar companion provides useful constraints on the physical properties of the 2mass  j1315@xmath02649 system , it also suggests that a binary interaction could be responsible for the unusual strength and persistence of its nonthermal emission . in a study of the comparably active t6.5e dwarf 2mass  j1237 + 6526 , @xcite proposed that accretion of material from a binary companion via roche lobe overflow could be a mechanism for sustained emission , as the inverted mass / radius relationship for brown dwarfs allows for sustained mass loss for @xmath74 , precisely the mass limit we find for 2mass  j1315@xmath02649ab .
however , the maximum separation for roche lobe overflow is only a few jupiter radii , well below the @xmath815,000  r@xmath57 projected separation for this pair .
moreover , the ballistic velocity of material impacting the surface of the primary , @xmath75  km  s@xmath14 @xmath8 500  km  s@xmath14
( where m@xmath76 and r@xmath76 are the mass and radius of the primary in jupiter units ) would have been readily detectable from line broadening over a broad range of viewing geometries .
combined with the absence of infrared excess associated with a circumstellar disk , we rule out accretion from the observed binary companion as the source of line emission .
the most likely explanation is that 2mass  j1315@xmath02649a possesses an unusually active chromosphere for its spectral type , perhaps reflecting an unusually strong magnetic field .
@xcite postulated an inverse relationship between age and activity among late - type m and l dwarfs , finding that stellar - mass ( and hence older ) l dwarfs were more likely to exhibit h@xmath6 emission than younger l dwarfs in their sample ( note that @xcite find marginal evidence for the opposite trend ) .
@xcite and infer a similar mass - dependence in the strength and persistence of magnetic fields on brown dwarfs based on the total energy flux available for field generation . in their model , objects above the hydrogen burning mass limit retain kilogauss fields for up to 10  gyr . if strong fields correlate with strong chromospheres in l dwarfs ( as they do for late - type m dwarfs ; @xcite ) , then the presence of nonthermal emission may align with kinematic , spectroscopic and binary mass ratio evidence that 2mass  j1315@xmath02649 is an old system , and 2mass  j1315@xmath02649a a relatively massive brown dwarf or low - mass star with a strong magnetic field .
a direct measurement of this component s magnetic field could be obtained through zeeman line broadening measurements in the 0.99  @xmath23 feh wing - ford band .
however , the high spectral resolution required for such a measurement ( @xmath11 @xmath77 30,000 ) makes it a challenge for this faint system ; indeed , no zeeman broadening measurements have been reported for an l dwarf to date .
an independent issue is how such a field could generate a persistent chromosphere in the presence of a highly neutral photosphere .
variability in the strength of h@xmath6 emission from 2mass  j1315@xmath02649 over the past decade may indicate microflaring as a viable source of heating , and a possible indicator of vigorous turbulent field generation @xcite .
alternately , @xcite have hypothesized that dust grain ionization and electron avalanche ( i.e. , lightning ) could locally increase the photospheric ionization fraction and magnetic field coupling in dusty l dwarfs , although early models appear to favor this mechanism in younger , lower - mass brown dwarfs and exoplanets .
a connection between magnetic activity and `` cloudiness '' has yet to be explored .    for completeness
, we note that because of its relatively wide separation , magnetic interaction between 2mass  j1315@xmath02649a and b is an unlikely source of emission in this system .
magnetospheric interactions over scales of up to @xmath6350  r@xmath78 have been implicated in outbursts from solar - mass t  tauri binaries . however , the estimated @xmath6330  r@xmath78 size of a typical l dwarf magnetic field
@xcite / r@xmath79 @xmath80 b@xmath81 . ]
is still several orders of magnitude smaller than the projected separation of the 2mass  j1315@xmath02649 binary .
in addition , the possibility that the observed nonthermal emission arises from the secondary , rather than the primary , is unlikely .
this scenario would imply a relative emission luminosity of @xmath7 @xmath8 @xmath02.5 for the secondary , an amplitude seen only during exceptionally large flare bursts from m and l dwarfs ( e.g. , @xcite ) . while other interactions could be postulated ( e.g. , magnetic wind interactions , kozai - like perturbations of an unseen third component )
, these scenarios are far more complex than the simple hypothesis of enhanced chromospheric emission powered by a strong magnetic field around a relatively old and massive cool dwarf . in other words , while the presence of a low - mass t dwarf secondary in this system is certainly intriguing , it appears to play no role in the observed nonthermal emission .
we have identified a t7 brown dwarf companion to the unusually active l5e dwarf 2mass  j1315@xmath02649 , a source that continues to exhibit strong h  i and alkali line emission despite its late spectral type . resolved imaging and spectroscopy
confirm the pair to be co - moving and co - spatial , and evolutionary models indicate an unusually low mass ratio as compared to other low - mass field binaries .
the spectral and kinematic properties of both components confirm prior indications that this system is relatively old ( @xmath65 @xmath3 0.81.0  gyr ) , and likely a member of the old galactic disk population .
joint comparison to atmospheric and evolutionary models ( a coevality test ) also supports an older age for this system , but reveals discrepancies in the case of the secondary ; we suspect these are due to continued shortcomings in the modeling cool brown dwarf spectra .
the age and separation of the system , coupled with the narrow emission lines and absence of mid - infrared excess , rule out accretion from or magnetic interaction with the t dwarf secondary as the emission source .
rather , we attribute it to an unusually strong magnetic field as predicted by energy flux scaling arguments for relatively old and massive low - mass dwarfs
. a direct test of this hypothesis can be achieved through zeeman line broadening measurements , although the source of chromospheric heating remains a separate issue . whether the larger sample of hyperactive cool dwarfs are also old , relatively massive , and possess strong magnetic fields remains to be determined , but they hint at a remarkable inversion of the standard age - activity relationship for low - mass stars .
the authors would like to thank paul sears and john rayner at irtf , al conrad , heather hersley and marc kassis at keck , and sergio vera and david osip at las campanas for their assistance in the observations reported here .
thanks also go to shelley wright for her help with the osiris data reduction , and our anonymous referee for her / his careful and timely critique .
this publication makes use of data from the two micron all sky survey , which is a joint project of the university of massachusetts and the infrared processing and analysis center , and funded by the national aeronautics and space administration and the national science foundation .
2mass data were obtained from the nasa / ipac infrared science archive , which is operated by the jet propulsion laboratory , california institute of technology , under contract with the national aeronautics and space administration .
this research has made use of the simbad database , operated at cds , strasbourg , france ; the very - low - mass binaries archive housed at http://www.vlmbinaries.org and maintained by nick siegler , chris gelino , and adam burgasser ; and the m , l , and t dwarf compendium housed at dwarfarchives.org and maintained by chris gelino , davy kirkpatrick , and adam burgasser .
the authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of mauna kea has always had within the indigenous hawaiian community .
we are most fortunate to have the opportunity to conduct observations from this mountain .                                                              , b. , leggett , s.  k. , lucas , p.  w. , pinfield , d.  j. , smart , r.  l. , day - jones , a.  c. , jones , h.  r.  a. , murray , d. , nickson , e. , tamura , m. , zhang , z. , lodieu , n. , tinney , c.  g. , & zapatero osorio , m.  r. 2010 , , 404 , 1952                                                                              , j.  d. , looper , d.  l. , burgasser , a.  j. , schurr , s.  d. , cutri , r.  m. , cushing , m.  c. , cruz , k.  l. , sweet , a.  c. , knapp , g.  r. , barman , t.  s. , bochanski , j.  j. , roellig , t.  l. , mclean , i.  s. , mcgovern , m.  r. , & rice , e.  l. 2010 , , 190 , 100          , a. , gasaway , t. , song , i. , iserlohe , c. , weiss , j. , larkin , j.  e. , barczys , m. , & lafreniere , d. 2004 , in presented at the society of photo - optical instrumentation engineers ( spie ) conference , vol .
5492 , society of photo - optical instrumentation engineers ( spie ) conference series , ed .
a.  f.  m.  moorwood & m.  iye , 14031410    , j. , barczys , m. , krabbe , a. , adkins , s. , aliado , t. , amico , p. , brims , g. , campbell , r. , canfield , j. , gasaway , t. , honey , a. , iserlohe , c. , johnson , c. , kress , e. , lafreniere , d. , lyke , j. , magnone , k. , magnone , n. , mcelwain , m. , moon , j. , quirrenbach , a. , skulason , g. , song , i. , spencer , m. , weiss , j. , & wright , s. 2006 , in presented at the society of photo - optical instrumentation engineers ( spie ) conference , vol .
6269 , society of photo - optical instrumentation engineers ( spie ) conference series                    , j.  l. , burles , s. , thompson , i.  b. , shectman , s.  a. , bigelow , b.  c. , burley , g. , birk , c. , estrada , j. , jones , p. , smith , m. , kowal , v. , castillo , j. , storts , r. , & ortiz , g. 2008 , in society of photo - optical instrumentation engineers ( spie ) conference series , vol .
7014 , society of photo - optical instrumentation engineers ( spie ) conference series                  , d.  g. , levine , s.  e. , canzian , b. , ables , h.  d. , bird , a.  r. , dahn , c.  c. , guetter , h.  h. , harris , h.  c. , henden , a.  a. , leggett , s.  k. , levison , h.  f. , luginbuhl , c.  b. , martini , j. , monet , a.  k.  b. , munn , j.  a. , pier , j.  r. , rhodes , a.  r. , riepe , b. , sell , s. , stone , r.  c. , vrba , f.  j. , walker , r.  l. , westerhout , g. , brucato , r.  j. , reid , i.  n. , schoening , w. , hartley , m. , read , m.  a. , & tritton , s.  b. 2003 , , 125 , 984                                                      , f.  j. , henden , a.  a. , luginbuhl , c.  b. , guetter , h.  h. , munn , j.  a. , canzian , b. , burgasser , a.  j. , kirkpatrick , j.  d. , fan , x. , geballe , t.  r. , golimowski , d.  a. , knapp , g.  r. , leggett , s.  k. , schneider , d.  p. , & brinkmann , j. 2004 , , 127 , 2948            , a.  a. , morgan , d.  p. , bochanski , j.  j. , andersen , j.  m. , bell , k.  j. , kowalski , a.  f. , davenport , j.  r.  a. , hawley , s.  l. , schmidt , s.  j. , bernat , d. , hilton , e.  j. , muirhead , p. , covey , k.  r. , rojas - ayala , b. , schlawin , e. , gooding , m. , schluns , k. , dhital , s. , pineda , j.  s. , & jones , d.  o. 2011 , , 141 , 97


levitation melting ( elm ) was invented in the twenties of the last century @xcite whereas its usage started only at the beginning of the fifties @xcite when high - frequency power generators became available .
the basic principle of elm is simple : a conducting sample , usually metallic , weighting from several tens of grams up to several kilograms is placed in a coil fed by an ac current with a typical frequency ranging from about @xmath2 up to several @xmath3 .
the ac magnetic field induces eddy currents in the sample which , in turn , give rise to two effects . on one hand ,
eddy currents interact with those in the coil giving rise to a lorentz force that repels the sample from the coil . on the other hand ,
the ohmic dissipation due to the induced currents provides heating of the sample .
in such a way , the sample can be levitated and also melted , provided that the coil configuration and the current in it are properly adjusted .
elm is particularly useful for melting reactive metals with high melting points , for example such as _ ti _ , _ zr _ , _ v _ , _ ta _ , _ mo _ , which often react with the crucible material and so get polluted by it @xcite .
elm avoids contamination of the melt and allows one to carry out solidification from deeply under - cooled states which is of interest for certain material science applications @xcite .
besides , elm is a well - known method for measurements of material properties of liquid metals used on ground as well as in space @xcite .
the balance of gravity and electromagnetic forces is necessary but not sufficient for a successful levitation .
in addition , the sample has to be stable at least to perturbations of sufficiently small amplitude .
first , it means that the reaction force due to the displacement of the sample from its equilibrium position has to act against that displacement .
otherwise , the equilibrium will be statically unstable and the solid sample will fall out of the coil or touch it when it is slightly perturbed .
similarly , a molten sample may leak out of the coil in the result of surface folding and deformation @xcite .
the static stability alone may also be insufficient for a successful levitation because sometimes the sample exhibits overstability .
namely , there is a restoring force that makes the sample to execute oscillations with increasing amplitude so that the sample eventually hits the coil or leaves it @xcite .
similarly , sometimes the sample is observed to spin - up @xcite .
a purely electromagnetic theory trying to account for such dynamic instabilities has been proposed in our previous work and applied to some simple configurations of magnetic fields @xcite . in this study
, we apply our previous theory to analyze both the static and dynamic stability of a spherical sample in more realistic axisymmetric magnetic systems made of a set of coaxial circular current loops .
although the basic configuration is axisymmetric , the perturbation fields are , in general , three - dimensional that renders the problem mathematically more complicated . to calculate 3d fields for a spherical sample we use similar analytic techniques as in refs .
however , our approach differs from the previous ones by an original use of the gauge transformation for the 3d vector potential in order to satisfy boundary conditions for the induced current .
this allows us to carry out the analysis only in terms of the vector potential without considering the scalar electrostatic potential .
therefore our approach is considerably simpler compared to @xcite .
the paper is organized as follows .
the stability of small amplitude oscillations of arbitrary direction is considered in sec .
2 beginning with the governing equations and analytic solutions for the vector potential of a circular current loop in spherical harmonics . in this section
we also derive solutions for the axisymmetric base state and the perturbation field due to a small displacement .
section 3 presents governing equations and analytic solution for the spin - up instability . numerical results for both instabilities in several simple inductors are discussed in section 4 and summarized in section 5 .
consider a sphere with radius @xmath4 and conductivity @xmath5 moving at velocity @xmath6 in a magnetic field @xmath7 alternating with circular frequency @xmath8 .
the induced electric field follows from the first maxwell equation as @xmath9 where @xmath10 is the scalar potential of the electric field and @xmath11 is the vector potential of the magnetic field .
the density of the electric current induced in a moving medium is given by ohm s law @xmath12 the electric field induced by a translational solid body motion with @xmath6 being a spatially invariant vector can be represented as @xmath13 assuming that the frequency of the alternating magnetic field is sufficiently low to neglect the displacement current , the second maxwell equation leads to the following advection - diffusion equation @xmath14 where the gauge invariance of @xmath11 has been employed to define the scalar potential as @xmath15 boundary conditions at the surface follow from the continuity of the magnetic field @xmath16_{s}=\left[\partial_{n}\vec{a}_{0}\right]_{s}=0 \label{eq : bcnd-1}\ ] ] where @xmath17_{s}$ ] denotes a jump of the quantity @xmath18 across the boundary @xmath19 ; @xmath20 is the derivate normal to the boundary .
sketch to the formulation of the problem . ]
the magnetic field of an inductor consisting of a set of coaxial circular loops represents a superposition of the fields of the separate loops .
therefore , we first consider a single loop lying parallel to the _ @xmath21 _ plane and centered with respect to the
_ z_-axis of a cartesian coordinate system with a sphere at its origin , see fig .
[ fig : sketch ] .
approximation of a toroidal inductor of small cross - section by a circular loop is considered in more detail in appendix 1 .
the loop is supplied with an ac current with amplitude @xmath22 alternating harmonically with a circular frequency @xmath8 as @xmath23 .
henceforth we use the magnetic diffusion time @xmath24 and the radius of the sphere @xmath4 as time and length scales while the current density and the vector potential are scaled by @xmath25 and @xmath26 , respectively . in the following ,
we assume all quantities to be dimensionless with the same notation as for the dimensional counterparts used so far .
further we consider the vector potential generated by the loop @xmath27 with the axisymmetric amplitude @xmath28 where @xmath29 is the dimensionless ac frequency which will be the main parameter throughout this paper .
the amplitude of the current density of such a loop can be represented by using the dirac @xmath30 function : @xmath31 where @xmath32 and @xmath33 are the meridional angle and the spherical radius of the loop with respect to the center of the sphere as shown in fig .
[ fig : sketch ] .
further , it is advantageous to introduce a complex unit vector @xmath34 that allows us to represent the azimuthal unit vector as @xmath35 where @xmath36 and @xmath37 are @xmath38 and @xmath39 unit vectors , respectively , and the asterisk denotes the complex conjugate
. then the current density can be expanded in spherical harmonics as @xmath40 where @xmath41 are spherical harmonics as defined in @xcite .
the expansion coefficients are @xmath42 using the generating function of the spherical harmonics @xcite the vector potential of a circular loop is obtained as @xmath43 where @xmath44 is the cylindrical radius of the loop ; @xmath45 ; @xmath46 . in order to facilitate the following algebra , it is advantageous to introduce the functions @xmath47 and @xmath48 , which are the outer and the inner solutions of the laplace equation associated to the spherical harmonic @xmath41 .
basic properties of those functions are given in appendix 2 .
for the region occupied by the sphere ( @xmath49 ) the solution above takes the form @xmath50 where @xmath51 for @xmath52 is assumed .      for a sphere at rest
the vector potential is sought as @xmath53,$ ] where @xmath54 is complex . for the interior of the sphere , ( [ eq : vecpot ] ) takes the form @xmath55 while for the exterior we have @xmath56 whereas ( [ eq : gauge ] ) takes the form of the coulomb gauge @xmath57 . for the interior ,
the solution is @xmath58 where @xmath59 are unknown coefficients to be determined from the boundary conditions ; @xmath60 is the spherical bessel function of index @xmath0 @xcite . for the exterior ,
the solution is represented as a superposition of external and induced fields @xmath61 where the latter is sought as @xmath62 from the conditions ( [ eq : bcnd-1 ] ) applied at the surface of the sphere at @xmath63 we find @xmath64 some useful properties and an efficient algorithm for the calculation of @xmath65 and its derivative are given in appendix 3 .
the time - averaged total force on the sphere is @xmath66dv,\ ] ] where the integral is taken over the volume of the sphere @xmath67 .
taking into account that @xmath68 and @xmath69 where @xmath70 is a purely real quantity , we obtain @xmath71\times\vec{b}_{0}^{e}dv.$ ] according to the momentum conservation law the force exerted by a magnetic field on a body is opposite to that exerted by the body on the source of the magnetic field @xmath72dv$ ] where the integral is taken over the space outside the sphere @xmath73 the last integral can be taken straightforwardly because the external current ( [ eq : j0-sph ] ) is defined in terms of a dirac @xmath30 function .
taking into account that the induced magnetic field outside the sphere is @xmath74 where @xmath75 , we find @xmath76 f_{l-1}f_{l}\left(\frac{l^{2}-1}{4l^{2}-1}\right)^{1/2 } \label{eq : f0}\ ] ] where @xmath77 involves the normalised legendre function @xmath78 related to the spherical harmonics as @xmath79 @xcite . in case of @xmath80 current loops we have @xmath81 where @xmath82 is the dimensionless current carried by the _ n_-th loop ; @xmath83 , @xmath84 , and @xmath85 are defined analogously as for the single loop .
note that the solution ( [ eq : f0 ] ) alone could be obtained in a simpler way as in ref .
our following analysis , however , is aimed at non - axisymmetric solutions which require a more complicated algebra .      consider a small displacement @xmath86 of the center of the sphere from the position @xmath87 to @xmath88 in the frame of reference related to the sphere , this corresponds to a perturbation of the external magnetic field @xmath89 where @xmath90 @xmath91 and @xmath92 now the external field represents a superposition of the axisymmetric base field @xmath93 and a small but in general non - axisymmetric perturbation @xmath94 .
the perturbation of the inner field @xmath95 is governed by the same eqs .
( [ eq : a0-ex ] , [ eq : a0-in ] ) as the base field with the difference that now the external field is given by @xmath94 .
consequently , the perturbation of the inner field and the corresponding induced field outside the sphere can be obtained straightforwardly by replacing the coefficients @xmath96 by @xmath97 in the corresponding solutions ( [ eq : a0-insol ] , [ eq : a0-exsol ] ) for the base field .
although the solution of the perturbed field obtained in such a way satisfies ( [ eq : a0-ex ] , [ eq : a0-in ] ) and the boundary conditions , it turns out that the induced field outside does not satisfy the coulomb gauge unless the problem is axisymmetric , _
i.e. _ , the offset is along the symmetry axis .
thus , although the laplace equation is satisfied we formally have @xmath98 which implies that the current is leaking outside the sphere . in addition note that the continuity of the tangential components of induction following from the continuity of the vector potential and its normal derivative at the surface of the sphere ensures continuity of the radial , _
i.e. _ normal , component of the current because @xmath99 .
thus , the vector potential satisfying the laplace equation and the coloumb gauge outside the sphere and being continuous at the surface ensures that the current is closed in the sphere .
as shown above , the induced vector potential unambiguously follows from the imposed one .
thus , the induced vector potential can be modified only by changing the imposed one .
this can be done by employing the gauge invariance which allows us to add to the external vector potential @xmath94 a gradient of some gauge potential defined as @xmath100 in order to fulfill the coulomb gauge for the external vector potential , @xmath101 has to be a harmonic function which is reflected in the expression above .
so , we obtain a set of free coefficients @xmath102 which can be chosen to satisfy the coulomb gauge for the induced field .
this gauge transformation results in the replacement of the original coefficients @xmath97 by @xmath103 then , similarly to the base field , the perturbation of the induced field is obtained as @xmath104 from the coulomb gauge we find @xmath105 the corresponding induction is @xmath106.\end{aligned}\ ] ] according to our previous theory @xcite , small - amplitude oscillations of the sphere with @xmath107 give rise to a perturbation of the force on the sphere @xmath108 = -\vec{k}\cdot\vec{x}+\vec{\g}\cdot\frac{d\vec{x}}{dt}.\ ] ] where @xmath109 and @xmath110 are the effective electromagnetic stiffness and damping tensors , respectively .
further we refer to the elements of @xmath109 and @xmath111 as the spring constants and growth rates , respectively .
note that due to the axial symmetry both @xmath109 and @xmath111 are purely diagonal matrices with zero azimuthal components . at first
, the axial symmetry implies that a purely axial displacement gives rise to a similar reaction force .
second , if the displacement is purely radial the reaction force , on the one hand , has to change to the opposite together with the displacement because of linearity but , on the other hand , it has to rotate around the symmetry axis together with the latter .
thus , a purely radial displacement causes a purely radial reaction force and , therefore , both @xmath109 and @xmath111 are purely diagonal .
both diagonal elements of @xmath109 have to be positive for the static stability whereas the corresponding growth rates have to be negative for dynamic stability .
note that the growth rates can become positive and , thus , destabilizing when the dimensionless ac frequency exceeds a certain threshold depending on the configuration of the field @xcite .
both the spring constants and the growth rates can be calculated for the given direction of the displacement specified by @xmath112 by using the corresponding base and perturbation fields @xcite : @xmath113\times\vec{b}_{0}^{e } + \re\left[\vec{j}_{0}\right]\times\vec{b}_{1}^{e}\right]dv\\ \vec{\g}\cdot\vec{\epsilon}_{x } & = & \frac{1}{2}\partial_{\bar{\omega } } \int_{v}\im\left[\vec{j}_{1}\right]\times\vec{b}_{0}^{e}dv.\end{aligned}\ ] ] similarly to the integral force , it is advantageous to change the region of integration from the sphere to the current loop @xmath114+\vec{j}_{1}^{e}\times \re\left[\vec{b}_{0}^{i}\right]\right]dv,\\ \vec{\g}\cdot\vec{\epsilon}_{x } & = & -\frac{1}{2}\partial_{\bar{\omega } } \int_{\bar{v}}\vec{j}_{0}^{e}\times\im\left[\vec{b}_{1}^{i}\right]dv.\end{aligned}\ ] ] after some algebra we find @xmath115\frac{f_{l}^{2}}{2l+1}(l^{2}-1)\\   & & { \times}\left[\frac{1}{2}\vec{e}_{z}\times\vec{\epsilon}_{x}\times   \vec{e}_{z}(1-\frac{1}{l})+\vec{e}_{z}(\vec{\epsilon}_{x}\cdot\vec{e}_{z})\right ] .
\nonumber \label{eq : gamma}\end{aligned}\ ] ] taking into account that @xmath116 because @xmath117 where @xmath118 is defined via dirac @xmath30 functions we eventually obtain @xmath119\left\{\frac{f_{l}^{2}}{2l+1}(l^{2}-1)\right.\\   & & \left .
{ \times}\left[\frac{1}{2}\vec{e}_{z}\times\vec{\epsilon}_{x}\times \vec{e}_{z}(1-\frac{1}{l})+\vec{e}_{z}(\vec{\epsilon}_{x}\cdot \vec{e}_{z})\right]\right.\\   & & \left .
{ -}\:\frac{f_{l-1}f_{l+1}}{2l-1}\left ( \frac{l(l+2)(l^{2}-1)}{(2l-1)(2l+3)}\right)^{1/2}\right.\\   & & \left . { \times}\left[\frac{1}{2}\vec{e}_{z}\times\vec{\epsilon}_{x}\times\vec{e}_{z } -\vec{e}_{z}(\vec{\epsilon}_{x}\cdot\vec{e}_{z})\right]\right\}.\end{aligned}\ ] ]
now , we can use this series to calculate both the stiffness and the damping coefficients for a given coil defined by the coefficients @xmath120 according to ( [ eq : fl ] ) .
we consider a solid sphere , as in section 2 , rotating with angular velocity @xmath121 in an alternating magnetic field .
the density of the induced current is given by ohm s law for a moving medium @xmath122 where @xmath123 is the velocity of a solid - body rotation .
the electric field induced by such a rotation may be represented as @xmath124 then from the second maxwell equation with neglected displacement current , that corresponds to the quasistationary approximation assumed throughout this study , we obtain @xmath125 note that this equation governs the vector potential in the laboratory frame of reference where the body rotates while the source of the magnetic field is at rest .
further we proceed to the frame of reference rotating together with the sphere . in a rotating frame of reference , where the sphere is at rest while the source of the magnetic field rotates with velocity @xmath126
, the equation above takes the form as for a body at rest @xmath127 where @xmath128 is a radius vector , which is time - dependent in the rotating frame of reference with @xmath129 . in a rotating frame of reference
, the body rotation appears as an additional time - dependence of the external magnetic field modulating the time - dependence of the applied magnetic field .
the time - averaged total torque on the sphere is @xmath130\times\vec{b}_{0}^{e}dv.\ ] ] according to the conservation of the angular momentum the torque exerted by the magnetic field on a body is opposite to that exerted by the body on the source of the magnetic field @xmath131dv,\ ] ] where the integral is taken over the volume @xmath132 outside the sphere .
the last integral can be taken straightforwardly because @xmath133 is defined via dirac @xmath30 functions . because @xmath134 the integral above may be represented as @xmath135dv\ ] ]
further we consider a sphere rotating about an axis perpendicular to the symmetry axis of an axisymmetric external magnetic field of a circular loop .
the coordinate system related to the loop with the _ @xmath136 _ axis as symmetry axis is chosen so that the rotation takes place about the @xmath38-axis , as shown in fig .
[ fig : sketch ] .
it is convenient to analyze this rotation in a coordinate system turned about the @xmath39-axis by 90 degrees so that the _ @xmath137_-axis of the transformed coordinate system coincides with the original @xmath38-axis which is the axis of rotation and , respectively , the _
@xmath138_-axis is directed opposite to the _ @xmath136 _ axis of the original coordinate system .
then the basis vectors of both coordinate systems are related as @xmath139 @xmath140 and @xmath141 respectively , the complex unity vector is @xmath142 where @xmath143 to represent the vector potential of a circular loop in the turned coordinate system we use the addition theorem for legendre polynomials @xcite written as @xmath144 together with the relation @xmath145 which yields @xmath146 where @xmath147.\end{aligned}\ ] ] similarly , the current density of the loop can be rewritten as @xmath148 in this case , we define a gauge potential @xmath149 where @xmath150 @xmath102 are unknown coefficients to be determined subsequently from the coulomb gauge for the vector potential induced outside the sphere .
addition of @xmath151 to @xmath152 results in the replacement of the coefficients @xmath153 by @xmath154 now we can simply proceed to the frame of reference rotating together with the sphere by redefining the azimuthal angle as @xmath155 where @xmath156 is the azimuthal angle in the rotating frame of reference and @xmath157 is the dimensionless angular velocity of rotation . respectively , the complex unit vector and the spherical harmonics change in a rotating frame of reference as @xmath158 and @xmath159 in a rotating frame of reference , the vector potential takes the form @xmath160 further we consider the current alternating harmonically with the dimensionless frequency @xmath161 as @xmath162 and the vector potential alternating in the laboratory frame of reference , respectively , as @xmath163.\ ] ] in the laboratory frame of reference , the complex amplitude of the vector potential induced outside the sphere is obtained as @xmath164 where @xmath165 and @xmath166 . from the coulomb gauge @xmath167 we find @xmath168&{}={}&\\   { } = { } \left [   \frac{(\vec{e}_{\eta}'^{*}\cdot\vec{i}_{l}^{m-1})}{\sqrt{2}\left(n_{l}^{m-1}\right)^{2 } } -\frac{(\vec{e}_{\eta}'\cdot\vec{i}_{l}^{m+1})}{\sqrt{2}\left(n_{l}^{m+1}\right)^{2 } } -\frac{(\vec{e}_{z}'\cdot\vec{i}_{l}^{m})}{\left(n_{l}^{m}\right)^{2 } } \right ] & & \end{aligned}\ ] ] for the field induced outside the sphere we obtain @xmath169.\end{aligned}\ ] ] the gauge potential @xmath170 vanishes in the expression above and so it does not affect the total torque which is @xmath171\\   & & { \times}\left[(l - m+1)(l - m+2)\bar{p}_{l+1}^{m-1}(0)^{2}\right.\\    & & \left .
{ -}\:(l+m+1)(l+m+2)\bar{p}_{l+1}^{m+1}(0)^{2}\right].\end{aligned}\ ] ] other components of the torque , which are not parallel to the angular velocity of rotation , vanish because of @xmath172 . for slow rotations with @xmath173 we have @xmath174 that results in @xmath175 , where @xmath176f_{l}^{2}}{(2l+1)^{2}(l+1)}\sum_{m=1}^{l}m\\   & & { \times}\left[(l - m+1)(l - m+2)\bar{p}_{l+1}^{m-1}(0)^{2}\right.\nonumber\\   & & \left .
{ -}\:(l+m+1)(l+m+2)\bar{p}_{l+1}^{m+1}(0)^{2}\right].\nonumber\end{aligned}\ ] ] for the rest state of the sphere to be stable to rotational perturbations , the torque caused by these perturbations has to damp them that corresponds to @xmath177 in the opposite case , the rest state is unstable .
the state is marginally stable when @xmath178 which is the equation defining the critical frequency @xmath179 for the spin - up instability .
in the following we present some numerical results for an inductor consisting of two circular loops of equal radii placed coaxially at the distance @xmath180 from the mid - plane .
the current having the same amplitude may flow either in the same or in opposite directions in both loops .
subsequently , both these cases are referred to as symmetric and antisymmetric , respectively .
we start with a simple symmetric case when both loops coincide ( @xmath181 ) making up a single loop with doubled current amplitude . the forces along with axial and radial spring constants versus the distance of the sphere from such a loop of dimensionless radius @xmath182
are shown in fig .
[ fig : force ] for various dimensionless frequencies . as
seen in fig .
[ fig : force](_a _ ) , at large enough distances the force , which is repulsing , raises with decreasing distance between the loop and the sphere . at some distance
, the force attains a maximum and tends to zero as the sphere approaches the center of the loop .
such a reducing force implies that the position of the sphere is statically unstable to axial perturbations .
the position is statically stable to axial perturbations when the axial spring constants shown in fig .
[ fig : force](_b _ ) are positive that requires the sphere to be placed further away from the loop than the point of maximum force which is at @xmath183 .
a position being axially stable may become unstable radially when the distance from the loop becomes too large ( see fig .
[ fig : force]_b _ ) .
thus the range of statically stable positions for a single loop is rather limited ( @xmath184 ) and does not change significantly with the frequency of the magnetic field .
further let us turn briefly to the static stability of a sphere in an arrangement of two coaxial current loops of the same radius placed symmetrically at distance @xmath180 from the center of the sphere . for the sake of brevity
we restrict the static stability analysis to the perfect - conductor approximation @xmath185 axial and radial spring constants are shown in fig .
[ fig : finfz0 ] versus the distance @xmath180 at various radii of the loops for both symmetric and antisymmetric arrangements . in the symmetric arrangement , which is usually used for heating ,
the position of the sphere is statically unstable in radial direction when the loops are too much separated , and it becomes unstable in the axial direction when the loops are moved too tightly together .
overlapping of the stability ranges for radial and axial perturbations depends on the radius of the loops .
for @xmath182 the range of axial stability ends at @xmath186 where the range of radial stability begins .
thus , there is no overlapping of the stability ranges in this case that implies a static instability regardless of the distance between the loops .
a range of positions statically stable to both axial and radial perturbations is possible only for sufficiently small radii of the loops .
for instance at @xmath187 this range is approximately @xmath188 in contrast to the symmetric arrangement , an antisymmetric one , which is usually used for positioning of the sample , ensures a statically stable state to radial perturbations regardless of the distance between the loops ( see fig .
[ fig : finfz0]_b _ on right ) .
the position of the sphere is stable to axial perturbations for any radius of the loop provided the loops are separated by a distance larger than the radius of the sphere ( @xmath189 ) .
there might be an axial instability for smaller separations of the loops at @xmath190 ( see fig .
[ fig : finfz0]_a _ on right ) .
( 0,0)(a )    ( 0,0)(b )      ( 0,0)(a )    ( 0,0)(b )    ( 0,0)(c )     forces together with axial and radial spring constants versus the distance of the sphere from the mid - plane of both symmetric and antisymmetric arrangements of two loops are shown in fig .
[ fig : finfr2 ] for various distances between the loops in the perfect - conductor approximation ( @xmath191 ) .
it is seen that the position of the sphere in the middle - point ( @xmath192 ) of a symmetric arrangement ( fig .
[ fig : finfr2 ] on left ) is unstable either in radial or axial directions for almost all distances between the coils .
however , there is one distance between the coils , @xmath186 , at which both spring constants cross zero almost simultaneously .
thus , this distance is an optimal one for the static stability of the sphere in the middle - point .
note that the distance @xmath193 provides the most uniform magnetic field in the vicinity of the middle - point of the symmetric arrangement .
this is because the next to the leading order contribution for the magnetic field of symmetric arrangement vanishes at the aspect ratio @xmath194 , _
i.e_. , @xmath195 in ( [ eq : a0e ] ) .
similarly , the distance @xmath196 providing the most uniform gradient of the magnetic field in the vicinity of the middle - point of an antisymmetric arrangement is expected to ensure the widest range of statically stable positions for loops of large enough radius .
( 0,0)(a )    ( 0,0)(b )    ( 0,0)(c )    ( 0,0)(d )     the dynamic stability of small - amplitude radial and axial oscillations is characterized by the corresponding growth rates @xmath197 and @xmath198 defined by ( [ eq : gamma ] ) .
the threshold of linear stability is given by the frequency above which the corresponding growth rate becomes positive . at the threshold frequency
, the amplitude of oscillations turns from decaying to growing in time .
a characteristic feature of this instability is that the growth rate , becoming positive above the threshold frequency , reaches a maximum at a certain frequency and tends to zero at higher frequencies .
the frequency at which the growth rate attains a maximum is referred to as the most dangerous one .
threshold frequencies along with the most dangerous one and the corresponding maximal growth rates are plotted in fig .
[ fig : oscz0 ] versus the half - distance between the loops at various radii of the loops for both axial and radial oscillations of a sphere about the middle - point of both symmetric and antisymmetric arrangements of two loops .
the most dangerous frequency is always above the threshold one . as seen in figs .
[ fig : oscz0](_a _ ) and [ fig : oscz0](_b _ ) , at large enough distances between the loops , the threshold and the most dangerous frequencies tend to @xmath199 and @xmath200 in the symmetric arrangement , and to @xmath201 and @xmath202 in the antisymmetric one , respectively .
note that these critical frequencies for the oscillatory instability in spatially quadratic and linear magnetic fields coincide with those for the spin - up instability in linear and uniform magnetic fields , respectively , which will be considered later .
for each radius of a symmetric arrangement , there is a certain optimal distance between the loops that maximizes the critical frequency . for large enough radii the optimal distance between the loops is approximately equal to the radius , @xmath203 , which is optimal also for the static stability as considered above . with decreasing of the loop radius ,
the optimal distance becomes slightly larger ( see figs . [
fig : oscz0]_a _ and [ fig : oscz0]_b _ on left ) . in the antisymmetric case ,
the critical frequency changes weakly on reducing of the distance @xmath180 until the loops approach a distance @xmath204 to the mid - plane ( see figs .
[ fig : oscz0]_a _ and [ fig : oscz0]_b _ on right ) . further approaching of the loops results in the increase of the critical frequency which is the larger , the smaller the radius .
this raise is caused by the nonuniformity of the magnetic field in the vicinity of the loop .
although the critical frequency for axial oscillations is , in general , slightly lower than that for the radial oscillations , there is no significant difference between both .
the maximal growth rate for the symmetric case changes weakly with reducing distance between the loops until the distance becomes comparable to the radius of the sphere ( @xmath205 ) ( see figs .
[ fig : oscz0]_c _ and [ fig : oscz0]_d _ on left ) .
closer approaching of the loops causes a significant increase of the growth rate for both axial and radial oscillations .
again , the smaller the radius of the loops , the larger the increase of the growth rate . in the antisymmetric case
, the growth rate reaches a maximum as the loops approach to the mid - plane at a distance comparable to the radius of the sphere ( @xmath206 ) ( see figs .
[ fig : oscz0]_c _ and [ fig : oscz0]_d _ on right ) .
spin - up instability occurs when the frequency of the alternating magnetic field exceeds the threshold @xmath179 defined above .
the spin - up rate being proportional to @xmath207 becomes positive as the frequency raises over the threshold and attains a maximum at some higher frequency @xmath208 subsequently referred to as the most dangerous one .
these critical frequencies together with the corresponding maximal spin - up rates for a sphere at the middle - point between symmetric and antisymmetric arrangements of two loops are plotted in fig .
[ fig : rotz0 ] versus the half - distance between the loops of various radii . as seen in fig .
[ fig : rotz0](_a _ ) , at large enough distances between the loops the threshold and the most dangerous frequencies tend to @xmath201 and @xmath202 in the symmetric arrangement , and to @xmath199 and @xmath200 in the antisymmetric one , which correspond to the spatially uniform and linear fields , respectively .
in the symmetric case , the critical frequencies have a minimum at the distance between the loops equal to their radius : @xmath203 ( see fig . [
fig : rotz0]_a _ on left ) .
note that this is converse to the critical frequencies for the oscillatory instability considered previously which have a maximum at this distance .
this minimum is due to the maximal uniformity of the field in the vicinity of the middle - point achieved at this distance between the loops .
critical frequencies increase with the non - uniformity of the field which becomes particularly significant when the loops are approached closer than the radius of the sphere ( @xmath209 ) . in this case , the smaller the radius of the loops , the higher the critical frequencies .
the spin - up rate raises as the loops are approached and attains a maximum at @xmath209 . for the loops with large enough radius
the maximum of the spin - up rate is attained when both loops merge together forming a single one ( see fig .
[ fig : rotz0]_b _ on left ) . in the antisymmetric case ,
the critical frequencies begin to raise significantly when the loops approach to the mid - plane at a distance of the radius of the sphere : @xmath186 ( see fig .
[ fig : rotz0]_a _ on right ) . as seen in fig .
[ fig : rotz0]_b _ on right , the spin - up rate attains a maximum approximately at the same distance .
( 0,0)(a )    ( 0,0)(b )
this work presents an analytic stability analysis of a solid sphere levitated electromagnetically in an axisymmetric magnetic field induced by a set of coaxial circular loops .
stability is considered with respect to both oscillations of small amplitude and arbitrary direction and rotations perpendicular to the axis of symmetry of the field . both oscillations and rotations
are found to become growing in time when the dimensionless frequency of the field exceeds a certain critical threshold @xmath179 depending on the configuration of the field for each kind of instability .
the growth rates of both instabilities become positive as the frequency is raised over the corresponding threshold , attain a maximum at some higher frequency and tend to zero at high frequencies .
the frequency @xmath208 , at which the maximum is attained , is referred to as the most dangerous one for the corresponding instability .
critical frequencies and the corresponding maximal growth rates are calculated for arrangements of two loops of the same radius carrying the same total current which may flow either in the same or in opposite directions . these arrangements , corresponding to the usual heating and positioning fields , are referred to as the symmetric and antisymmetric ones .
it is found that critical frequencies for both oscillatory and rotational instabilities increase with the nonuniformity of the magnetic field .
the lowest dimensionless critical frequencies are @xmath201 and @xmath202 which are the same for both the spin - up instability in a uniform magnetic field and the oscillatory instability in a spatially linear field . this coincidence represents a particular case of a more general rule : the critical frequencies for the spin - up instability in a spherically harmonic field of degree @xmath0 coincide with the critical frequencies for oscillatory instability in a spherically harmonic field of degree @xmath1 .
thus , the critical frequencies @xmath199 and @xmath200 for the spin - up instability in a linear field coincide with the critical frequencies for the oscillatory instability in a spatially quadratic magnetic field which is the case at the mid - point of a symmetric arrangement of two loops except the case that the distance between the loops is equal to the radius of the loops . in this particular case ,
the quadratic term vanishes and the instability is dominated by the fourth - order term which has critical frequencies @xmath210 and @xmath211 .
thus , the distance @xmath203 is optimal in order to avoid the oscillatory instability of the sphere at the mid - point of the symmetric arrangement .
this distance , which ensures the most uniform magnetic field at the middle - point , yields a minimum of the critical frequency for the spin - up instability .
note that there is no similar optimal distance for the oscillatory instability in antisymmetric arrangement where the critical frequency can be increased only by the nonuniformity of the magnetic field when the loops are approached close to the surface of the sphere .
hence , the stabilization of a sphere by an optimal design of the inductor is rather limited and active means of stabilization , like an additional steady magnetic field , may be necessary .
note that the maximal dimensionless growth rates for both oscillatory and spin - up instabilities are small , typically @xmath212 . from the physical point of view
we have introduced the growth rate as a viscous - type friction coefficient at the velocity .
thus the physical dimension of @xmath213 is @xmath214 and , consequently , its scale is @xmath215 because we have used the magnetic diffusion time @xmath216 and the radius of the sample @xmath4 as time and length scales , respectively . in order to obtain the actual friction force we have to take the product of velocity , friction coefficient and the magnitude of the characteristic total electromagnetic force @xmath217 which usually is comparable to the gravity of the sample @xmath218 , where @xmath219 is the mass of the sample and @xmath220 the free fall acceleration . by comparing the negative effective friction force to the inertia
we obtain an estimate of the characteristic growth time of the instability @xmath221 which for a characteristic size @xmath222 and a conductivity @xmath223 of the sample leads to @xmath224 .
thus , although the magnetic field has a dynamically destabilizing effect , the development of the instability is expected to be very slow and , thus , hardly observable . in conclusion , note that the frequency of the magnetic field is the only parameter determining the threshold of dynamic instabilities as long as no external damping , for instance due to a surrounding gas , is taken into account . in the opposite case , the threshold of instability would depend not only on the frequency but also on the amplitude of the current .    in conclusion note that the dynamic instabilities resulting from the effect of motion of a conducting body in an ac magnetic field can be interpreted from an alternative physical point of view which can explain the weakness of this instability and also suggests other possible instability mechanisms . in the given ac magnetic field ,
the electromagnetic force on the spherical body at rest depends only on its position .
if the position changes , it takes some time for the electromagnetic force to relax to its time - averaged stationary value for the given position .
this delay is caused by the finite time of magnetic diffusion in the body which , taking place over the skin depth , is thus comparable to the ac oscillation period @xmath225 the electromagnetic force at instant @xmath226 on the body in motion may be represented as @xmath227 assuming a sufficiently high ac frequency , we can expand this representation in a power series of small @xmath228 yielding @xmath229 thus , the delay of the magnetic field results in a force component proportional to the velocity which is obviously analogous to a viscous friction force . since for a statically stable position
the effective electromagnetic reaction force @xmath230 is directed against the displacement @xmath86 , the effective friction force @xmath231 is directed with the motion .
consequently , the delay of the magnetic field results in an effective electromagnetic friction force with a negative and , thus , destabilizing friction coefficient .
it is important to note that the friction coefficient is small because it involves the delay time which is short and comparable to the ac period when the relaxation of the magnetic field is solely due to the electromagnetic diffusion over the skin depth .
the same arguments may be extended to the coupling of the magnetic field and the temperature distribution in the body caused by its temperature - dependent electrical conductivity . in this case , the delay of the magnetic field would be dominated by the thermal relaxation time which is much longer than the electromagnetic diffusion time and , thus , expected to cause stronger dynamic instabilities than the pure electromagnetic mechanism considered here .
in addition , the decrease of the electrical conductivity for increasing temperature , which is caused by the joule heating of the ac magnetic field , and the related reduction of the non - dimensional frequency @xmath161 might be important in levitation experiments .
consider a toroidal inductor of arbitrary cross - section carrying an axisymmetric and purely azimuthal current .
the magnetic field of such an inductor may be represented as a superposition of the fields of separate circular loops constituting the inductor , and the corresponding vector potential may be written as the integral over the inductor cross - section @xmath19 : @xmath232 where @xmath233 is the nondimensionalized vector potential of a single circular loop defined by the radius vector @xmath234 with @xmath235 and @xmath137 being the radius and axial position of the loop in cylindrical coordinates .
subsequently assume the skin effect to be negligible as if the inductor would be fed by a direct current driven by the gradient of an electrostatic potential @xmath236 an axisymmetric and purely azimuthal current distribution with density @xmath237 implies @xmath10 to depend solely on the azimuthal angle @xmath238 .
consequently , @xmath239 where @xmath240 is a constant that can be related to the total current in the inductor @xmath241 and its effective radius @xmath242 as @xmath243 further we assume the inductor cross - section to be located about some position @xmath244 which will be specified later , and approximate the field distribution by the following multipole type expansion : @xmath245\frac{ds'd\phi'}{\left|\vec{r}-\vec{r}_{0}\right| } , \label{eq : multexp}\end{aligned}\ ] ] where the operator @xmath246 acts on @xmath247 now , if we define @xmath52 as the position of the mass center of the cross - section @xmath248 the second dipole - like term cancels in the expression above which , thus , reduces to    @xmath249    the last expression defines the vector potential of a circular loop located at @xmath52 and carrying the effective current @xmath250 for example , in the case of an inductor represented by a torus with major and minor radii @xmath33 and @xmath251 the effective radius can easily be found from ( [ eq : effr ] ) as @xmath252 while the mass center of the cross - section coincides with its geometrical center .
note that the solution ( [ eq : a_approx ] ) is accurate up to the quadrapole - like term neglected in ( [ eq : multexp ] ) the magnitude of which relative to the remaining term may be estimated as @xmath253 where @xmath254 is the characteristic size of the cross - section and @xmath4 is the characteristic distance to the cross - section center .
obviously , the neglected term becomes significant only at distances comparable to @xmath255
similarly to ref . @xcite we also use a complex unity vector @xmath258 where @xmath36 and @xmath37 are the corresponding cartesian unit vectors , and employ the outer and inner solutions of the scalar laplace equation , @xmath259 and @xmath260 associated with the spherical harmonic @xmath261 @xcite .
these variables provide a simple algebra for calculation of the gradient operator @xmath262 used throughout this study : @xmath263 where the asterisk denotes the complex conjugate and @xmath75 .
the corresponding relations for the complex conjugate functions follow straightforwardly from the above relations and the property of spherical harmonics : @xmath264 another relation used for the calculation of the torque in sec .
3 is : @xmath265.\end{aligned}\ ] ] note that our definition of @xmath266 is slightly different from that used in @xcite leading to a bit more complicated algebra but simpler resulting expressions .
the coefficients @xmath268 and @xmath269 in ( [ eq : f0],[eq : gamma ] , [ eq : dm ] ) can efficiently be calculated using continued fractions @xcite . for this purpose
we rewrite @xmath270 where @xmath271 and @xmath272 is the spherical bessel function of index @xmath273 further , applying a recurrence relation for the spherical bessel function of index @xmath0 + 1 we obtain the continued fraction : @xmath274 allowing us to calculate @xmath275 provided that @xmath276 is known .
but asymptotic properties of bessel functions suggest that for large index @xmath0 @xmath277 .
thus for any @xmath136 we can choose sufficiently large @xmath278 and then truncate the fraction by approximating @xmath279 by the previous expression that allows us to calculate back the necessary @xmath280 note that the forward recurrence for @xmath275 is not practically applicable because it is numerically unstable similarly to its counterpart for bessel functions .
the other necessary quantity @xmath281 contained in ( [ eq : gamma],[eq : dm ] ) , can be calculated in a similar way by expressing it as @xmath282 .
for @xmath283 we have @xmath284 which follows from the corresponding asymptotics of the bessel functions @xcite .
i. egry , a. diefenbach , w. dreier and j. piller , `` containerless processing in space - thermophysical property measurements using electromagnetic levitation , '' int .
j. thermophys .
2 , pp . 569578 , 2001 .