galaxy mergers lead to intense star formation , i.e. , starburst @xcite .
the gravitational perturbations and dynamical disturbances induced by the merging process remove angular momentum from the galactic interstellar medium ( ism ) , the ism thus settles into the gravitational center of the system leading to starburst @xcite .
the collective effects of supernovae and stellar winds from massive stars within the starburst region drive large - scale outflows , or `` superwinds '' ( e.g. , * ? ? ? * ) .
galaxy mergers and their successive processes , starbursts and superwinds , significantly alter the morphology , stellar population , gas content , and chemical abundance of the galaxy .
together with statistical investigations based on a large sample of mergers , detailed case studies of merging starburst galaxies provide us with important clues for understanding the above process and its impact on galaxy evolution .
ngc  6240 is a well studied merger / starburst galaxy in the local universe ( @xmath9 ; * ? ? ?
it has a highly disturbed morphology indicating that two spiral galaxies are in the process of colliding with each other .
ngc  6240 is often considered as an ultra - luminous infrared galaxy ( ulirg ) in the local universe , although its infrared luminosity ( @xmath10@xmath11 @xmath12 : * ? ? ?
* ) is slightly smaller than that of the criterion for ulirg ( @xmath13(8  1000 @xmath14m)@xmath15 @xmath12 ; * ? ? ?
this system includes double nuclei ( the northern and southern nuclei ; hereafter referred to as the n - nucleus and the s - nucleus ) whose separation is @xmath11 kpc .
these nuclei are hard x - ray point sources suggesting that these are the colliding active galactic nuclei ( agns ) of the progenitor galaxies @xcite .
the s - nucleus is @xmath13 times brighter than the n - nucleus in the 0.210 kev band @xcite .
it is well known that a bright optical emission line nebula is associated with ngc  6240 .
early studies by various authors ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) suggested that this optical nebula is primarily excited by shock - heating induced by a starburst superwind .
high spatial resolution _
observations have revealed that the central h@xmath0 nebula ( `` butterfly nebula '' ) is spatially coincident with the soft x - ray - emitting hot gas surrounding the double nuclei @xcite .
in addition , there is a highly extended ionized gas region surrounding ngc  6240 .
@xcite detected faint h@xmath0 + [ ] emission filaments extending @xmath1630 kpc from the center .
@xcite revealed a more extended optical emission line nebula approximately 80 kpc in diameter .
an extended hot gas nebula of ngc  6240 has also been found in soft x - rays @xcite .
a spatial correlation between the optical emission line nebula and the extended soft x - ray gas was identified by @xcite .
recently , @xcite showed the extent and detailed structure of the very extended ( its total size reaches @xmath1100 kpc ) soft x - ray - emitting gas .    although the central region of ngc  6240 has been studied many times , the detailed structure of the extended optical emission line gas is not yet understood .
here we present the results of deep optical narrow - band imaging of ngc  6240 using the subaru suprime - cam .
our deep observations are the first to reveal the detailed structure of the extended h@xmath0 emitting warm ionized gas surrounding this galaxy .
we assumed cosmological parameters of ( @xmath17 , @xmath18 , @xmath19 ) = ( 0.705 , 0.27 , 0.73 ) @xcite and that the luminosity distance of ngc  6240 is 107 mpc .
the linear scale at the galaxy is 492 pc arcsec@xmath3 under this assumption .
lccc n - a - l671 & 0.88 & 27.5 & 10 @xmath20 600 s + @xmath21 & 0.91 & 28.5 & 5 @xmath20 300 s + @xmath22 & 0.76 & 27.7 & 5 @xmath20 180 s + @xmath23 & 0.87 & 27.3 & 5 @xmath20 180 s +    [ tab : obs ]    narrow - band h@xmath0 + [ ] and broad - band @xmath21 , @xmath22 , and @xmath23 imaging observations of ngc  6240 were made using the suprime - cam @xcite attached to the subaru telescope on april 30 ( ut ) 2014 . a journal of the observations
is shown in table [ tab : obs ] .
we used the n - a - l671 narrow - band filter @xcite for h@xmath0 + [ ] imaging of the galaxy ( @xmath9 ) .
the central wavelength and full - width at half maximum ( fwhm ) of the transmission curve of the filter were 6712   and 120  respectively .
the total exposure time was 6,000 sec for n - a - l671 .
the sky was clear and photometric .
seeing was 0709 on the observing night ( table [ tab : obs ] ) .
data reduction was carried out in a standard manner .
overscan subtraction , flat - fielding , distortion correction , sky subtraction , and mosaicking were performed with custom - made suprime - cam data - reduction software , nekosoft @xcite . during the overscan subtraction procedure
, we also conducted cross - talk correction and gain correction between the ccds ( see * ? ? ?
we used sky flats , which had been created from a number of object frames taken during the same observing run , for the flat - fielding correction of the @xmath21 and @xmath22 band frames . for
the n - a - l671 and @xmath23 frames , dome flat data were used for flat - fielding .
we then extracted an area of 400  @xmath20 400  ( 197 kpc @xmath20 197 kpc ) centered on the s - nucleus of ngc  6240 from each flat - fielded frame .
after adjusting the positions and scales , the frames for each band were combined using a 3 sigma clipping method .
the flux calibration was performed using the 9th data release of the sloan digital sky survey ( sdss ) photometric catalog @xcite .
the details of the flux calibration are described in appendix a and also in @xcite .
the fwhm of the point - spread function ( psf ) was measured by fitting a two - dimensional ( 2d ) gaussian function to the brightness profiles of non - saturated star images for each band .
the limiting surface brightness was estimated by photometry using randomly - sampled 2 arcsec apertures on the blank regions .
the psf fwhms and limiting magnitudes of the final @xmath21 , @xmath22 , @xmath23 and n - a - l671 images are listed in table [ tab : obs ] .
the reduced n - a - l671 image is shown in the right panel of figure [ fig : r - h - alpha ] .    to obtain a net ( pure ) h@xmath0 + [ ] image of ngc  6240 , we subtracted the scaled @xmath22 frame from the n - a - l671 frame .
scaling was performed using several unsaturated field stars common to both frames . before subtraction , we convolved 2d double gaussian functions to the @xmath22 frame to adjust the psfs of both of the frames . after subtraction
, we multiplied the net h@xmath0 + [ ] image by a factor of 1.23 to correct the over - subtraction of the h@xmath0 + [ ] flux due to the contamination of the emission lines in the @xmath22 frame ( see appendix b ) .
in addition , we subtracted the net h@xmath0 + [ ] that was multiplied by a factor of 0.23 from the scaled @xmath22 frame to create the emission - line - free @xmath22 frame .
the emission - line - free @xmath22 image is shown in the left panel of figure [ fig : r - h - alpha ] .
the net h@xmath0 + [ ] image is shown in figure [ fig : h - alpha ] using three different contrasts .
the h@xmath0 + [ ] image thus created is affected by the positional variation of the emission line ratios in the nebula .
we assumed constant emission line ratios over the entire region of the h@xmath0 nebula in the above procedure .
however , it is well known that there is a tendency for the flux ratios of the forbidden lines to the h@xmath0 emission line to be large at larger distances from the center of ngc  6240 @xcite .
this causes the systematic underestimation of the net h@xmath0 + [ ] flux at the outer region of the nebula .
in addition , the combination of the recession velocity variation of the nebula and the transmission curve of the n - a - l671 filter affects the detected h@xmath0 + [ ] flux .
the error of the h@xmath0 + [ ] flux due to this effect was estimated to be @xmath10@xmath245% ; the variation of the transmission of the filter is @xmath105%  in the recession velocity range of the ionized gas nebula of ngc  6240 @xcite .
the morphology of ngc  6240 in optical broad bands ( e.g. @xmath22 band ) is complicated , but there are some characteristic features in it .
the bright part of the system is elongated as long as @xmath120 kpc with a position angle ( pa ) of @xmath130@xmath25 and the clear dust lane runs along this structure . in the southeastern side of this bright feature , a `` banana''-shaped tail
whose pa is @xmath1@xmath26 is connected to the center .
a `` s''-shaped extended faint tidal feature surrounds these bright parts of ngc  6240 ( figure [ fig : r - h - alpha ] ) .    from a morphological point of view
, we suspect that the elongated bright part would be the remnant of the disk of the primary galaxy of this colliding system .
we tried to extract the disk component from the central bright part by fitting an exponential disk to the @xmath23 band image .
we used galfit @xcite for profile fitting .
we added four curved s@xmath27rsic components to fit residuals of the exponential disk fitting . in the fitting procedure
, we excluded the circular region whose radius is 3at the center to avoid the influence of saturation and agn components .
in addition , we also excluded bright stars and blooming patterns around the galaxy .
cccc 12.67 & 4.7 kpc & 28@xmath28 & 62@xmath28 +    [ tab : expdisk ]    the best fit parameters of the exponential disk are listed in table [ tab : expdisk ] .
this exponential disk component dominates in the @xmath23 band brightness ; its brightness is more than half of the sdss @xmath29 band magnitude , 11.99@xmath240.05 @xcite .
the other s@xmath27rsic components are 23 magnitude fainter than the disk , thus do not contribute to the continuum light significantly .
@xcite performed a fabry - perot spectroscopy of the central region of ngc  6240 and found a rotating disk component whose pa and inclination are @xmath145@xmath28 and @xmath170@xmath28 , respectively .
these values are roughly consistent with the parameters of the exponential disk we fitted .
thus we conclude that there still is a large stellar disk , which is probably the remnant of the primary galaxy of the merger , in the center of the system .
we refer to this disk as the `` main stellar disk '' of ngc  6240 .
the net h@xmath0 + [ ] image obtained in this study revealed the rich and complex structure of the optical emission line nebula of ngc  6240 ( figures [ fig : h - alpha ] , [ fig : h - alpha - color ] , and [ fig : colorimg ] ) .
we refer to this nebula as the `` h@xmath0 nebula '' .
the h@xmath0 nebula has a diameter of @xmath190 kpc .
the central `` butterfly '' nebula ( e.g. * ? ? ?
* ) is surrounded by a relatively bright hourglass shaped region ( `` hourglass region '' ; light blue color in figure [ fig : h - alpha - color ] ; @xcite ) .
in addition to these well - known features , @xcite also detected a faint optical emission line region extending out to over 70 kpc @xmath20 80 kpc around ngc  6240 .
our deep suprime - cam observations are the first to reveal the detailed and complex structure of this extended component of the h@xmath0 nebula .
we applied an unsharp masking technique to the net h@xmath0 + [ ] image to enhance the filamentary structure of the h@xmath0 nebula .
we first masked bright stars and blooming patterns in the h@xmath0 + [ ] image .
second , we smoothed the masked image with a 35 @xmath20 35 pixels ( @xmath103.4 kpc @xmath20 3.4 kpc ) median filter , then subtracted the smoothed image from the masked image . the unsharp masked image is shown in figure [ fig : unsharp ] .
it is clear that a number of filaments are running radially from the galactic center .
a sketch of the structure of the h@xmath0 nebula is shown in figure [ fig : sketch ] .
we identified several characteristic features in this nebula . in the following subsections
, we describe these features in detail .
the most remarkable feature of the extended h@xmath0 nebula is a bright , wiggled filament ( `` w - filament '' ) in the western region .
it is prominent in the middle panel of figure [ fig : h - alpha ] ( see panel ( c ) of figure [ fig : detail ] for the detailed structure ) .
although the w - filament was previously visible in the narrow - band h@xmath0 images taken by @xcite and @xcite , its detailed structure and surrounding complex features were only revealed in our h@xmath0 + [ ] image .
the bright part of the w - filament begins at 18  kpc west from the s - nucleus and extends to the west - southwest .
the filament turns slightly west - northwest at @xmath125 kpc and reaches @xmath136 kpc from the s - nucleus .
there is a bright cloud to the east of the bright region of the w - filament ( @xmath112 kpc west from the s - nucleus ) .
we assume , from a morphological point of view , this could be a part of the w - filament ( figure [ fig : sketch ] ) .      on the northern side of the galaxy ,
there exists a complex of clouds and filaments ( `` n - complex '' : figure [ fig : sketch ] and panel ( a ) of figure [ fig : detail ] ) , which extends along the pa @xmath1310@xmath28 from @xmath125 kpc north to @xmath133 kpc north - northwest of the galaxy .
the unsharp masked image indicates that the n - complex is the tip of a large curved structure connected to the hourglass region ( figure [ fig : unsharp ] ) .
a clear bubble - like structure is visible in the southeastern side of the galaxy ( `` se - bubble '' : figure [ fig : sketch ] and panel ( b ) of figure [ fig : detail ] ) .
this faint bubble was unknown before this study .
the size and pa of the se - bubble are @xmath128 kpc @xmath20 38 kpc and @xmath1135@xmath28 , respectively . on the other side of the galaxy , the n - complex and w - filament appear to form a large broken bubble whose size is comparable to the se - bubble ( `` nw - bubble '' ) .
the pa of the nw - bubble axis is approximately @xmath3040@xmath28 .
the axes of these two bubbles are almost perpendicular to the major axis of the main body of ngc  6240 .
the opening angle of the se - bubble and the nw - bubble is @xmath170@xmath2890@xmath28 ( figures [ fig : unsharp ] and [ fig : sketch ] ) .
the western side of the h@xmath0 nebula is much brighter and has a more complex structure than the eastern side .
in particular , the w - filament and its surroundings have a rich structure consisting of many filaments and blobs .
there are several bright knots ( `` bk1''``bk3 '' : figure [ fig : sketch ] ) located in a loop - like filament extending from the central nebula towards the southwest ( `` w - loop1 '' : figure [ fig : sketch ] ) .
these knots are marginally resolved in our h@xmath0 + [ ] image ; the fwhms of the knots are @xmath101011 .
w - loop1 extends almost straight towards the southwest from the center , then turns to the north at @xmath123 kpc drawing a loop - like trajectory ( panel ( d ) of figure [ fig : detail ] ) .
one of the bright knots ( bk3 ) is located at the tip of this loop .
another loop ( `` w - loop2 '' ) is located at @xmath135 kpc from the center ( panel ( e ) of figure [ fig : detail ] ) .
the size and shape of w - loop2 are similar to those of w - loop1 , but its surface brightness is much lower than that of w - loop1 .
there are two streams extending along the north - south direction on the south side of the galaxy ( `` s - filaments '' and `` sw - filaments '' : figure [ fig : sketch ] ; see also figure [ fig : unsharp ] ) . the s - filaments consist of several filaments extending from the center towards the south . the sw - filaments consist of a number of filaments distributed from the region close to w - loop2 to @xmath159 kpc southwest from the s - nucleus
the extension of the s - filaments is similar to but slightly offset from that of the southern tidal feature which is the southern end of the s - shaped morphology of the optical continuum emission of ngc  6240 .
there is no clear optical continuum counterpart to the sw - filaments .
we performed multiple aperture photometry for the h@xmath0 emission of ngc  6240 .
we selected the region further out than 140  ( @xmath170 kpc ) from the s - nucleus for sky background estimation . before performing photometry ,
bright stars and artifacts were masked and interpolated by linear functions around the patterns .
we adopted circular apertures whose center was the brightest point of the h@xmath0 + [ ] image corresponding to the s - nucleus .
we excluded the very center ( @xmath31 ) region in this aperture photometry procedure , because the @xmath22 band frame that was subtracted from the n - a - l671 frame to create the net h@xmath0 + [ ] image was saturated at the central 1  region of the galaxy .
the contamination of the [ ] @xmath326548,6583 emission was corrected assuming @xmath336583 } / f_{\rm h\alpha } = 1 $ ] @xcite for the entire region of the nebula .
in addition , we assumed that the foreground extinction is the same as the @xmath34 band galactic extinction toward ngc  6240 ; @xmath35 = 0.165 magnitude @xcite .
the surface brightness profile of the h@xmath0 emission is given in figure [ fig : sb ] .
it is clear that the h@xmath0 surface brightness distribution indicates the core - halo structure of the nebula , which is similar to the case of the soft x - ray nebula @xcite , whereas the core is more compact in h@xmath0 .
the total h@xmath0 flux of ngc  6240 was derived by integrating the fluxes in the regions at which the h@xmath0 surface brightness is greater than 2@xmath36 of the sky background .
we assumed that the h@xmath0 surface brightness within the central 1  region which we excluded in the aperture photometry is the same as that at 1  away from the s - nucleus .
as mentioned above , we assumed that @xmath336583 } / f_{\rm h\alpha}$ ] @xmath37 for the whole region of the nebula .
this ratio tends to increase toward the outer part of the h@xmath0 nebula , reaching @xmath161.4 around the w - filament @xcite .
thus our method may underestimate the contribution of the [ ] emission in some parts of the nebula .
however , since most of the bright regions of the nebula have almost constant @xmath336583 } / f_{\rm h\alpha}$ ] @xmath38 @xcite , the variation of this ratio does not significantly affect the estimation of the total h@xmath0 flux .
the resultant total h@xmath0 flux @xmath39 is @xmath101.3 @xmath20 10@xmath40 erg  s@xmath3  @xmath41 .
the h@xmath0 luminosity @xmath42 is @xmath101.6 @xmath20 10@xmath43 erg  s@xmath3 .
cccccc 0 @xmath30 1.7 & @xmath44 & @xmath45 & @xmath46 & @xmath47 & @xmath48 + 1.7 @xmath30 3.4 & @xmath49 & @xmath50 & @xmath51 & @xmath52 & @xmath53 + 3.4 @xmath30 7.4 & @xmath54 & @xmath55 & @xmath56 & @xmath57 & @xmath58 + 7.4 @xmath30 11 & @xmath59 & @xmath60 & @xmath61 & @xmath62 & @xmath63 + 11 @xmath30 15 & @xmath64 & @xmath65 & @xmath66 & @xmath67 & @xmath68 + 15 @xmath30 20 & @xmath69 & @xmath70 & @xmath71 & @xmath72 & @xmath73 + 20 @xmath30 27 & @xmath74 & @xmath75 & @xmath76 & @xmath77 & @xmath78 + 27 @xmath30 39 & @xmath79 & @xmath80 & @xmath81 & @xmath82 & @xmath83 +    [ tab : eden ]    we derived the root mean square ( rms ) electron density @xmath84(rms ) of the h@xmath0 nebula assuming case b recombination .
we did not , however , apply an extinction correction to this calculation .
the surface brightness of h@xmath0 @xmath85 is defined as @xmath86 where @xmath87 and em are the photon production coefficient of h@xmath0 in case b recombination @xcite and the emission measure , respectively . assuming a fully - ionized pure hydrogen gas ( @xmath88 ; @xmath89 is proton density ) and a constant density distribution along the line of sight , the em can be written as @xmath90 , where @xmath91 is the line of sight length of the nebula . in the case of
spherically symmetric nebula , @xmath92 , where @xmath34 and @xmath93 are the radius of the nebula and the projected distance from the nebula center , respectively .
thus , the @xmath84(rms ) is given by @xmath94 the radial profile of the @xmath84(rms ) distribution is shown in figure [ fig : ne ] .
the mass of the h@xmath0 emitting gas @xmath95 is given by @xmath96 where @xmath97 and @xmath98 are the volume filling factor of the h@xmath0 emitting gas and the mass of a hydrogen atom , respectively .
we estimated the volume filling factor @xmath97 by the following procedure .
first , considering the close spatial correlation between the h@xmath0 nebula and the soft x - ray - emitting hot gas ( see section [ compx ] below ) , we assumed that the h@xmath0 gas and x - ray gas are locally coexisting , then we considered the following two cases for the coexistent state ; 1 ) a simple thermal pressure equilibrium case and 2 ) a ram pressure balance case .    in the case 1 ) , we calculated the local electron densities @xmath99 of the h@xmath0 emitting cloud assuming that the equation @xmath100 holds over each radial zone in which @xcite provided the averaged electron densities @xmath101 and temperatures @xmath102 of the hot gas individually . here , @xmath103 is the volume filling factor of the hot gas .
because the optical emission line filaments of ngc  6240 show clear shock excitation nature in its emission line ratios @xcite , we assumed that @xmath104 k ( typical temperature of optical shock excitation gas @xcite ) for all of the zones in this case .
the case 2 ) is the one that shock excitation is responsible for the h@xmath0 emission gas and x - ray emitting hot gas .
we calculated the @xmath99 assuming the ram pressure balance equation @xcite , @xmath105 , where @xmath106 and @xmath107 are the shock velocities at the shock front between the hot gas and warm gas , respectively .
the shock velocity required to produce x - ray emitting gas with a temperature of @xmath110@xmath108 k is @xmath11000 km  s@xmath3 for an adiabatic shock @xcite .
@xcite found that the averaged [ ] @xmath109/h@xmath0 and [ ] @xmath326717 + 6731/h@xmath0 emission line intensity ratios of the optical filaments of ngc  6240 are @xmath10.3 and @xmath10.8 , respectively .
these values are consistent with shock models with the shock velocity @xmath1100 km  s@xmath3 ( e.g. , * ? ? ?
hence , we assume that @xmath106 and @xmath107 are @xmath11000 km  s@xmath3 and @xmath1100 km  s@xmath3 , respectively .
we then derived @xmath97 by @xmath110 for the both cases . in the case 1 ) , @xmath97 is given by @xmath111 in the case 2 ) , we obtain @xmath112 where @xmath113
.    the volume filing factor @xmath103 is not well known , but the soft x - ray would be strongly enhanced at the upstream of the shocks in the wind @xcite and the @xmath103 of such shocked gas clouds is considerably smaller than 1 .
@xcite found that the typical @xmath103 of the x - ray gas in a superwind is of order 10@xmath114 . assuming that @xmath115 , @xmath104 k , and @xmath116 , we found that the resultant volume filling factors are orders of 10@xmath410@xmath5 for the both cases ( table [ tab : eden ] ) .
although these values are much smaller than the typical @xmath117 of regions or planetary nebulae in our galaxy ( @xmath110@xmath11410@xmath118 ; e.g. * ? ? ?
* ) , they are comparable to the @xmath117 derived for outflow / wind gas by the observations of nearby starburst galaxies and agns @xcite .    under the assumption of spherically symmetric nebula , we obtained @xmath119 @xmath7 using the @xmath97 derived above .
the spherical symmetry is , however , too simple assumption for the geometry of the emission line nebula .
hence , we next considered a more realistic geometry for the h@xmath0 nebula ; it has a filamentary structure and the emission line filaments occupy a fraction @xmath120 ( @xmath121 ) along the line of sight of the spherical nebula . in this case , the rms electron densities are increased by a factor of @xmath122 . since we think that the x - ray gas and h@xmath0 gas spatially coexist , the line - of - sight occupation factor @xmath120 is the same for the both gas .
thus the factor @xmath120 vanishes in the calculation of the volume filling factor ( see equations ( [ eq : fv1 ] ) and ( [ eq : fv2 ] ) ) . on the other hand ,
the total volume of the nebula is decreased by a factor of @xmath120 . as a result
, the @xmath95 alters by a factor of @xmath123 .
then we obtained @xmath124 @xmath7 for this partly occupied filamentary nebula case .
we found that the h@xmath0 nebula shows a very good spatial correlation with the highly extended soft x - ray - emitting hot gas nebula @xcite , both in large and small scales . a tight spatial correlation between the h@xmath0 and soft x - ray is well known for the butterfly nebula of ngc  6240 @xcite .
the overall coincidence between the optical emission line gas and the extended soft x - ray gas of ngc  6240 has been previously pointed out by some authors @xcite .
our deep observations , however , allow us to make a more detailed comparison between the hot and warm phase gases .
we overlaid the contours of the h@xmath0 + [ ] image on the @xmath1150 ks deep soft x - ray image to inspect the spatial correlation between the h@xmath0 nebula and the x - ray nebula in , more detail . ]
taken with _
@xcite in figure [ fig : x - alpha ] . the overall shape ( a `` diamond shape '' ) of the faint emission level of the soft x - ray image is well traced by the h@xmath0 + [ ] emission . on smaller scales
, there are a number of common features between the x - ray and h@xmath0 + [ ] , including the central butterfly nebula , the hourglass region , the n - complex , the se - bubble , the w - filament , and w - loop2 .
there is also an indication of the counterpart of s - filaments in the x - ray image , but the southern tips of the filaments are located outside the x - ray - emitting gas .
there seems to be no clear counterpart to the sw - filaments in the x - ray .
a very tight spatial correlation between the warm ionized gas and the hot gas is maintained from the central region to the edge of the nebula of ngc  6240 .
we additionally compared the h@xmath0 nebula image with ultraviolet ( nuv and fuv ) images of ngc  6240 taken with _
( figure [ fig : uv - h - alpha ] ) .
there also is a good spatial correlation between the h@xmath0 and the uv images , particularly in the western side of the galaxy .
the western uv emission of ngc  6240 is mostly overlapped with the western complex of filaments of the h@xmath0 nebula . the w - filament and n - complex
can be traced in the uv images ( figure [ fig : uv - h - alpha ] ) .
the bright knots ( bk1  bk3 ) are also visible in the uv images , and the s - filaments are well traced by the uv emission .
there also are faint uv counterparts for the sw - filaments .
we found that the overall morphology of the h@xmath0 nebula is dominated by the radially - extended filamentary structure ( figure [ fig : unsharp ] ) .
we also identified a large - scale bi - polar bubble ( se- and nw - bubble ) along with prominent loop structures ( w - loop1 and 2 ) .
these morphological characteristics indicate that most of the h@xmath0 nebula consists of radially - outflowing gas  the starburst superwind .    the close spatial coincidence between the h@xmath0 nebula and the extended soft x - ray gas supports this idea .
many common features , including the se - bubble , the n - complex , the w - filament , and w - loop2 , are observed in the extended ionized gas of ngc  6240 in both h@xmath0 and x - ray . a similar morphological correlation between the optical emission line region and
the soft x - ray gas has been widely observed in starburst galaxies @xcite .
such good spatial coincidence has been interpreted as an indication that the hot and warm phase gases are heated and excited via the same mechanism  probably from shocks in the outflowing gas  over the entire region of the nebula .
numerical simulations of superwinds have supported this scenario @xcite .
@xcite obtained the optical spectra of two of the bright knots ( bk1 and bk2 ) , and found that they are active star - forming regions .
although there are no spectral data for bk3 , we suggest that it is also a star - forming region because the continuum counterparts of the bright knots that show a blue color ( @xmath125 ( ab ) @xmath126 ) are clearly seen in our broad - band images ( figure [ fig : knots ] ) and @xmath23 images in figure [ fig : knots ] were not corrected for emission lines , we think that the most flux in these bands comes from the stellar continuum . in the @xmath21 band , there are several emission lines such as h@xmath127 , h@xmath128 and [ ] @xmath1295007 .
we estimated that the contribution of these lines to the @xmath21 band flux is less than 10% assuming the ordinary balmer decrement and @xmath130}$]/@xmath131 @xcite . in the @xmath23 band
, there is no strong emission lines . ] .
it appears that the bright knots are interacting with the outflowing filament , and the gas in the filament is intercepted by these knots , creating w - loop1 .
tail - like structures are also observed in the broad - band images , in particular , the tails extending to the southwest are clearly visible in the @xmath21 band image ( figure [ fig : knots ] ) .
indications of the tails are also marginally seen in the @xmath23 band image .
this feature is recognizable in a hubble space telescope ( hst ) public image of ngc  6240 .
the head - tail morphology of the knots suggests that the outflow gas and condensed gas clouds in the main stellar disk interact with each other and that its ram pressure and/or shock may induce active star formation .
the bright knots are also visible in the uv images supporting the view that these knots are the sites of active star formation . in all probability , the young stellar population observed in the tails of the bright knots was also formed by this process .
the s - filaments are well traced by the uv emission .
there also are faint uv counterparts for the sw - filaments . because these filaments have either a weak or no spatial correlation with the x - ray emission , we consider that the s- and sw - filaments are not closely related to the starburst winds .
these features are more likely to be star - forming tidal tails .
the uv emission associated with the optical filaments may partly be the precursor emission of a shock within the nebula .
the extended uv emission associated with a superwind has been observed in nearby starburst galaxies including m82 and ngc  253 ( e.g. * ? ? ? * ) . in these galaxies ,
most of the uv emission is thought to be from the stellar continuum scattered by dust entrained in the wind @xcite , because it is too bright to be provided by shock - heated or photoionized gas in the wind .
it is well - known that ngc  6240 is a dusty system , thus it is possible that the extended uv emission of ngc  6240 is mostly derived from light scattered by dust .    to estimate the contribution of dust scattering , we compared the uv intensity and h@xmath0 intensity at the w - filament .
we measured the uv and h@xmath0 intensities in a 195 @xmath20 105 rectangular region at the w - filament ( the center of the photometry region is located 54  west and 5  north from the s - nucleus ) , and determined that @xmath132 erg  s@xmath3  @xmath41 , @xmath133 erg  s@xmath3  @xmath41 , and @xmath134 erg  s@xmath3  @xmath41 ( assuming that @xmath1356583}/f_{\rm h\alpha } = 1.4 $ ] , @xcite ) .
the average h@xmath0 surface brightness in this region was @xmath136 erg  s@xmath3  @xmath41  arcsec@xmath114 .
both the @xmath137 and @xmath138 ratios  7.4 and 3.4 , respectively ,  are slightly smaller than , but almost consistent with , values in the outflow filaments in m82 and ngc  253 at the same h@xmath0 surface brightness level @xcite .
the ratios are , however , still larger than the values predicted by the shock model .
this indicates that the uv emission of ngc  6240 is too bright to be attributed to shock - heated or photoionized gas emission in the wind .
thus , some part of the uv emission is possibly light scattering by dust .
alternatively , this emission may reflect extended star - forming activity within the filaments . in this case
, gas compression by the interaction of the outflowing gas and the interstellar medium in the main stellar disk may have induced star formation .
@xcite found that the x - ray hot gas of ngc  6240 has a spatially - uniform sub - solar metal abundance whose abundance pattern is consistent with type - ii supernova enrichment .
they discussed the age of the nebula using the thermal velocity dispersion estimated from the temperature of the hot gas , and concluded that the whole ngc  6240 nebula could not be formed by the recent activity of the nuclear starburst or agns .
they suggested that extended continuous star formation from the early phase of the merging event would be responsible for the formation of the outer nebula . the extended uv emission and
its spatial coincidence with the h@xmath0 emission might support this scenario .
most of the h@xmath0 emission of the h@xmath0 nebula , however , probably originates from shocks in the outflow , because the extended ionized gas shows a large [ ] /h@xmath0 intensity ratio ( [ ] /h@xmath0 @xmath161 ; @xcite ) .
detailed polarimetry and spectroscopy of the extended nebula are necessary to qualitatively understand the contribution of dust - scattering and local star formation to the h@xmath0 nebula .
lccccccc butterfly & 4.5 & @xmath139 & 700 & 25 & @xmath140 & @xmath141 + hourglass & 13 & @xmath142 & 700 & 25 & @xmath143 & @xmath144 + outer & 45 & @xmath145 & 700 & @xmath146 & @xmath147 & @xmath148 +    [ tab : windparam ]    we discuss the ages and durations of the starburst superwinds of ngc  6240 . from a morphological point of view , we divided the h@xmath0 nebula into three regions ; the central butterfly nebula @xcite , the hourglass region @xcite , and the outer filaments .
the characteristic radii of the three regions are @xmath14.5 kpc , @xmath113 kpc , and @xmath145 kpc , respectively .
we assume that these three regions were formed in the past by different starburst / superwind activities , and estimate the ages of the regions and the durations of these activities .
first , we estimate the outflow velocity of the winds and derive the dynamical ages of the regions .
@xcite studied the kinematics of the optical ionized gas of the central region of ngc  6240 and found that the fwhm of the h@xmath0 line exceeds 1000 km  s@xmath3 within 45 kpc from the center ( see also * ? ? ?
* ; * ? ? ?
they also found that some parts of the extended region have radial velocities of @xmath1500 km  s@xmath3 .
we thus adopted a value of 500 km  s@xmath3 as the line - of - sight component of the outflow velocity of the gas .
assuming that the typical inclination of the outflow is 45@xmath28 , we estimate the outflow velocity of the h@xmath0 nebula to be @xmath149 km  s@xmath3 .
we also assume that the outflow velocity is constant over the whole nebula .
the ages of the three regions are thus @xmath18 myr , @xmath124 myr , and @xmath184 myr , for the butterfly nebula , the hourglass region , and the outer filaments , respectively .
next , we estimate the durations of the superwind activities that formed the three regions based on the standard superwind model .
the terminal velocity of a mass loading galactic wind is given by @xmath150 where @xmath151 and @xmath152 are the thermalization efficiency and mass loading factor of the wind , respectively @xcite .
substituting @xmath153 km  s@xmath3 into the equation ( [ eq : vinf ] ) , we find @xmath154 .
although the exact value of @xmath151 is unknown and is difficult to estimate , it would be expected to lie between 0.1 and 1 @xcite . on the other hand ,
the mass - loss rate from the starburst region is scaled by the star formation rate ( sfr ) and given by @xmath155 @xcite .
then the duration of the superwind @xmath156 is calculated by @xmath157 , leading to @xmath158 where @xmath159 is the total mass of the wind including the loaded mass .
the sfrs are estimated as follows .
since the butterfly nebula was formed most recent star formation activities have played important roles in its formation .
we found that there is a wide spread in the estimated values of the current sfr . a well - used equation for converting the far infrared luminosity to the sfr @xcite , @xmath160
@xmath7  yr@xmath3 @xmath161 @xmath162 , gives us an sfr of @xmath1120 @xmath7  yr@xmath3 .
this value is clearly overestimated , because a significant fraction of the bolometric luminosity of ngc  6240 is known to originate from the agns and extended star formation other than the current nuclear starburst @xcite .
@xcite found a sfr of @xmath125 @xmath7  yr@xmath3 for the nuclear starburst based on the @xmath163 band luminosity . on the other hand , @xcite found a sfr of 83 @xmath7  yr@xmath3 based on the 1.4 ghz radio continuum observation .
@xcite noted that the discrepancy between the above two values arises from the difference in the star formation histories they used ; @xcite assumed a constant star formation history over 20 myr , whereas @xcite adopted a merger - induced increasing star formation scenario .
the morphology of the h@xmath0 nebula suggests an intermittent star formation nature for the star - forming history of ngc  6240 .
in particular , the estimated age of the butterfly nebula is much shorter than the duration of the star formation assumed by @xcite .
we thus adopt a sfr of 25 @xmath7  yr@xmath3 for the butterfly nebula .
the age of the hourglass region estimated above ( @xmath124 myr ) is almost consistent with the starburst age ( @xmath11525 myr ) estimated by @xcite .
we thus interpret that the hourglass region was created by this starburst , and its sfr was 25 @xmath7  yr@xmath3 , which is the value estimated by @xcite from the central red stellar population .
we have no suitable observational basis for estimating the sfr for the outer filaments .
hence , for the outer filaments we tentatively assume a typical sfr from the local starburst galaxies , @xmath110 @xmath7  yr@xmath3 .    the h@xmath0 gas masses of the three regions were determined by integrating the h@xmath0 gas mass in three concentric regions ; @xmath93 @xmath164 4.5 kpc , 4.5 kpc @xmath165 @xmath93 @xmath164 13 kpc , and 13 kpc @xmath165 @xmath93 @xmath164 45 kpc .
the close morphological correlation between the h@xmath0 emission and soft x - ray of ngc  6240 indicates that these two - phase gases are spatially associated with each other .
therefore , the mass of the soft x - ray gas must be taken into account to derive the total mass of the superwinds @xmath159 .
@xcite estimated the total mass of the soft x - ray gas as @xmath110@xmath166 @xmath7 with the assumption that the volume filling factor @xmath103 is almost unity and a spherically symmetric geometry of the hot gas nebula . however , this assumption is not adequate for the superwind . as we described in section [ sec : ed ] , @xmath103 would be of order 10@xmath114 .
in addition , we introduced the line - of - sight occupation factor @xmath120 to represent the filamentary structure of the nebula ( see section [ sec : ed ] ) . using these factors
, we calculated the total wind mass @xmath159 .    the ages and durations of the superwinds that formed the three regions derived by the above calculations are listed in table [ tab : windparam ] .
the duration of the superwind is approximately one order of magnitude smaller than the age of each region .
this suggests that the starburst / superwind activity has an intermittent nature as suggested by numerical simulations of galaxy mergers @xcite .
it should be noted that the duration of the superwind for the hourglass region ( 4 myr ) is consistent with the starburst duration ( @xmath15 myr ) estimated by @xcite .    in the above discussion , we ignored the agn contribution to the wind evolution . currently , there are double active nuclei at the center of ngc  6240 and they play a non negligible role in energetics in the central region of the galaxy @xcite .
the agns are , however , deeply embedded and still not active enough to form the extended hot and warm gas nebula surrounding the galaxy @xcite .
thus the most part of the nebula was probably formed by the starburst activity @xcite . even in the vicinity of the nuclei ,
the influence of the agns to the surrounding excited gas region has not been clear yet @xcite .
therefore , we conclude that the starburst superwind is the prime mover in forming the extended nebula and the contribution of the agns would be negligible except for the central region ( a part of the butterfly nebula ) .
ngc  6240 is thought to be in the late stages of a galaxy merger @xcite .
its morphology indicates that two gaseous spiral galaxies first encountered each other @xmath11 gyr ago and their nuclei are now in the final stages of coalescence .
the h@xmath0 nebula was formed by intermittent starburst events which occurred within the last 100 myr . here
, we propose a possible scenario of the starburst history of ngc  6240 from the point of view of the formation of the h@xmath0 nebula by starburst superwinds .    a starburst which occurred
@xmath180 myr ago created powerful bi - polar galactic bubbles ; the nw- and se - bubbles . at that time
, the ism in the merging system might be more abundant in the northwestern and western sides of the merging galaxies than in the eastern and southeastern sides . a combination of the intermittent nature of the starburst activity and the lopsided ism distribution formed a rich structure of loops , filaments , or blobs in the western side of the system .
the extended active star formation as indicated by the uv images was also preferentially induced in the western side as the wind propagated outward from the central starburst .
the bright knots observed in the western side of the system were probably formed by collisions between the high - speed outflowing gas from the center and dense gas blobs .
the ram pressure of the outflow served to highly condense the gas blobs and initiated active star formation .
approximately 20 myr ago , a short duration starburst ( @xmath145 myr ) occurred at the center of the system .
this starburst and the subsequent superwind formed the central stellar population of ngc  6240 and the hourglass region of the h@xmath0 nebula .
active star formation ceased temporarily following the short duration burst described above .
several myr ago , the concentration of gas by the continuing approach of the two nuclei induced an intense starburst that continues to the present ; and also drove agn activity .
the sfr may be steadily increasing , and has currently reached @xmath125 @xmath7  yr@xmath3 @xcite .
this recent starburst and agn activity formed various features including the butterfly nebula @xcite , the central extended hard x - ray gas @xcite , the molecular outflow @xcite , and the intense h@xmath167 and [ ] emission @xcite .
we conducted unprecedented deep h@xmath0 narrow - band imaging observations of ngc  6240 using the subaru suprime - cam .
the 6,000-sec suprime - cam observation was the first to reveal the detailed structure of a highly - extended h@xmath0 emitting nebula ( the h@xmath0 nebula ) whose size reaches to @xmath190 kpc in diameter .
this nebula was probably formed by large - scale energetic outflows  or superwinds  driven by starburst events induced by a galaxy  galaxy merger .
in addition to the butterfly - shaped inner nebula , which has already been investigated by various authors , the outer nebula has a complex structure consisting of wiggled narrow filaments , loops , bubbles , and blobs .
we found a huge bi - polar bubble ( the se - bubble plus the nw - bubble ) , which extends along the minor axis of the main stellar disk of ngc  6240 . the total h@xmath0 luminosity is @xmath101.6 @xmath20 10@xmath43 erg  s@xmath3 .
we estimated the volume filling factor @xmath97 of the h@xmath0 nebula assuming a simple thermal pressure equilibrium and a ram pressure balance between the warm ionized gas and the soft x - ray gas .
the @xmath97 is of orders 10@xmath410@xmath5 for the both cases .
we calculated the mass of the h@xmath0 emitting gas to be @xmath95 of @xmath15 @xmath168 @xmath7 using the @xmath97 .
it is remarkable that the h@xmath0 extent and morphology are quite similar to the x - ray ( 0.71.1 kev ) emitting hot gas detected with a deep _ chandra _ observation @xcite .
such a close spatial coincidence between the optical emission line nebula and the hot gas is generally observed in starburst superwinds .
this strongly suggests that the same physical mechanism , shock - heating , plays an important role in the excitation of both phases of the gas at the outer regions of the h@xmath0 nebula .
the h@xmath0 nebula also has a close spatial correlation with the uv emission in the western region of the galaxy .
some bright star - forming knots are located in the western region .
these facts suggest , therefore , that star - forming activity was induced by the interaction of the outflow and the ism , and that the interaction was probably enhanced due to a lopsided gas distribution formed by the merger .
the morphology and surface brightness distribution of the h@xmath0 nebula suggest that there were probably three major episodes of starburst / superwind activities that formed the nebula .
the most extended part of the nebula was initiated @xmath180 myr ago .
multiple intermittent starbursts in this era created many outflow filaments and clouds extending @xmath140 kpc out from the center . at the same time , the outflowing gas interacted with the ism and induced extended star formation in the main stellar disk .
an intense short duration starburst then occurred @xmath120 myr ago , forming the hourglass region of the h@xmath0 nebula .
ngc  6240 entered a new active phase @xmath11 myr ago .
the galaxy is now in the final stages of the merger and its activity is increasing .
we thank the anonymous referee for the careful reading and constructive suggestions .
we are grateful to the staff of the subaru telescope for their kind help with the observations .
this work was in part carried out using the data obtained by a collaborative study on the cluster abell  1367 .
this study was in part carried out using the facilities at the astronomical data analysis center ( adac ) , national astronomical observatory of japan .
this research made use of nasa s astrophysics data system abstract service , nasa / ipac extragalactic database , galex gr6 database , and sdss dr9 database .
this work was financially supported in part by grant - in - aid for scientific research no.23244030 , no.15h02069 from the japan society for the promotion of science ( jsps ) , no.24103003 from the ministry of education , culture , sports , science and technology ( mext ) , and most grant 104 - 2112-m-001 - 034- .
ahn , c. p. , et al .
2012 , , 203 , 21    armus , l. , et al .  1990 , , 364 , 471    barnes , j. e.  1998 , in galaxies : interaction and induced star formation , saas - fee advanced course 26 , ed .
d. friedli , l.martinet , & d. pfenniger ( new york : springer ) , 275    beswick , r. j. , pedlar , a. , mundell , c. g. , & gallimore , j. f.  2001 , , 325 , 151    bland - hawthorn , j. , wilson , a.  s. , & tully , b. 1991 , , 371 , l19    blanton , m. r. , & roweis , s.  2007 , , 133 , 734    brown , m. j. i , et al .
2014 , , 212 , 18    cecil , g. , bland - hawthorn , j. , veilleux , s. , & filippenko , a. v.  2001 , , 555 , 338    cecil , g. , bland - hawthorn , j. , & veilleux , s.  2002 , , 576 , 745    cooper , j. l. , bicknell , g. v. , sutherland , r. s. , & bland - hawthorn , j.  2008 , , 674 , 157    cooper , j. l. , bicknell , g. v. , sutherland , r. s. , & bland - hawthorn , j.  2009 , , 703 , 330    downes , d. , solomon , p. m. , & radford , s. j. e.  1993 , , 414 , l13    elston , r. , & maloney , p.  1990
, , 357 , 91    engel , h. , et al .
2010 , , 524 , a56    feruglio , c. , et al .
2013 , , 549 , a51    genzel , r. , tacconi , l. j. , rigopoulou , d. , lutz , d. , & tecza , m.  2001 , , 563 , 527    gerssen , j. , et al .
2004 , , 127 , 75    heckman , t. m. , armus , l. , & miley , g. k.  1987 , , 92 , 276    heckman , t. m. , armus , l. , & miley , g. k.  1990 , , 74 , 833    hollenbach , d. , & mckee , c. f.  1979 , , 41 , 555    hoopes , g. g. , et al .
2005 , , 619 , l99    hopkins , p. f. , et al .
2006 , , 163 , 1    hopkins , p. f. , et al .
2008 , , 175 , 356    hopkins , p. f. , et al .  2013 , , 433 , 78    huo , z. y. , xia , x. y. , xue , s. j. , mao , s. , & deng , z. g.  2004 , , 611 , 208    hutton , s. , et al .
2014 , , 440 , 150    iono , d. , et al .
2007 , , 659 , 283    keel , w. c.  1990 , , 100 , 356    kennicutt , r. c.  1998 , , 498 , 541    komatsu , e. , et al .
2011 , , 192 , 18    komossa , s. , et al .
2003 , , 582 , l15    gerssen , j. , et al .
2004 , , 127 , 75    lira , p. , ward , m. j. , zezas , a. , & murray , s. s.  2002 , , 333 , 709    lonsdale , c. j. , farrah , d. , & smith , h. e.  2006 , astrophysics update 2 ( chichester : praxis ) , 285    lutz , d. , et al .  2003 , aap , 409 , 867    mihos , j. c. & hernquist , l.  1994 , , 437 , l47    mihos , j. c. & hernquist , l.  1996 , , 464 , 641    miyazaki , s. , et al .  2002 , , 54 , 833    mori , t. , et al .
2014 , , 66 , 93    nakanishi , k. , et al .
2005 , , 57 , 575    nardini , e. , et al .
2013 , , 765 , 141    ohyama , y. , et al .
2000 , , 52 , 563    ohyama , y. , yoshida , m. , & takata , t.  2003 , , 126 , 2291    osterbrock , d. e. & ferland , g. j.  2006 , astrophysics of gaseous nebulae and active galactic nuclei ( sausalito , ca : univ .
sci . )    peng , c.  y. , ho , l.  c. , impey , c.  d. , & rix , h.  -w .
2002 , , 124 , 266 .
peng , c.  y. , ho , l.  c. , impey , c.  d. , & rix , h.  -w .
2010 , , 139 , 2097 .
robinson , a. , et al .
1994 , , 291 , 351    rupke , d. s. , veilleux , s. , & sanders , d. b.  2005 , , 160 , 87    sanders , d. b. & mirabel , i. f.  1996 , ara&a , 34 , 749    schlafly , e. f. & finkbeiner , d. p.  2011
, , 737 , 103    scoville , n. , et al .  2015 , , 800 , 70    schweizer , f.  1998 , in galaxies : interaction and induced star formation , saas - fee advanced course 26 , ed .
d. friedli , l.martinet , & d. pfenniger ( new york : springer ) , 105    sharp , r. g. & bland - hawthorn , j.  2010 , , 711 , 818    shull , j. m. , & mckee , c. f.  1979 , , 227 , 122    smith , d. a. & wilson , a. s.  2001 , , 557 , 180    spitzer , l.  1978 , physical processes in the interstellar medium ( new york : wiley - interscience )    strickland , d. k. & stevens , i. r.  2000 , , 314 , 511    strickland , d. k. , heckman , t. m. , kimberly , a. w. , & dahlem , m.  2000 , , 120 , 2965    strickland , d. k. , heckman , t. m. , colbert , j. m. , hoopes , c. g. , & weaver , k. a.  2004 , , 151 , 193    sugai , h. , malkan , m. , ward , m. , davies , r. , & mclean , i.  1997 , , 481 , 186    tacconi , l. j. , et al .
1999 , , 524 , 732    tanaka , m. , hasegawa , t. , & gatley , i.  1991 , , 374 , 516    tecza , m. , et al .
2000 , , 537 , 178    teyssier , r. , chapon , d. , & bournaud , f.  2010 , , 720 , l149    tullmann , r. , pietsch , w. , rossa , j. , breitschwerdt , d. , & dettmar , r .- j .  2006 , , 448 , 43    van der werf , p. p. , et al .
1993 , , 405 , 522    veilleux , s. , shopbell , p. l. , rupke , d. s. , bland - hawthorn , j. , & cecil , g.  2003 , , 126 , 2185    veilleux , s. , et al .
2005 , 43 , 769    wang , j. , et al .
2014 , , 781 , 55    yagi , m. , et al .
2002 , , 123 , 66    yagi , m. , et al .
2007 , , 660 , 1209    yagi , m.  2012 , , 124 , 1347    yagi , m. , et al .
2013 , , 65 , 22    yoshida , m. & ohtani , h.  1993 , , 45 , 407    yoshida , m. , et al .
2008 , , 688 , 918
we adopted a color conversion from sdss to the suprime - cam filter system with hpk ccds as @xmath169 where the coefficients and the color range of the stars are given in table [ tab : conv ] .
the coefficients are estimated using the method given in @xcite , with the model quantum efficiency curve of the ccds .
we used stars of magnitude @xmath170 in sdss dr9 .
the k - correct @xcite v4 offset for sdss was applied ; @xmath171 0.012 , 0.010 and 0.028 for @xmath172 , @xmath93 , and @xmath29 , respectively .
the number of stars is approximately 60 - 200 , and the estimated rms error around the best fit ranged from 0.05 to 0.07 .
cccccccccccccc @xmath173 & @xmath174 & @xmath175 & @xmath176 & @xmath177 & 0.070 & @xmath178 & @xmath179 & 0.881 & ... & ... + @xmath180 & @xmath181 & @xmath182 & 0.006 & 0.312 & @xmath183 & @xmath184 & 1.601 & @xmath185 & @xmath186 & 9.992 + @xmath187 & @xmath181 & @xmath188 & 0.020 & 0.536 & 0.021 & @xmath189 & @xmath190 & ... & ... & ... + @xmath191 & @xmath181 & @xmath192 & @xmath193 & 0.089 & 0.019 & @xmath194 & ... & ... & ... & ... +    [ tab : conv ]
the use of the scaled @xmath22 frame as an off - band frame leads an over - subtraction of the continuum , because the h@xmath0 + [ ] emission is non - negligibly contaminated in the scaled @xmath22 frame .
the optical emission line nebula of ngc  6240 primarily consists of shock - heated gas ( e.g. * ? ? ?
* ; * ? ? ?
* ) , thus the [ ] @xmath1296300 and [ ] @xmath326717,6731 lines must also contribute towards contaminating the emission line flux in the @xmath22 frame .
assuming that the emission line intensity ratio @xmath195+[s ii ] } / f_{\rm h\alpha+[n ii]}$ ] is constant over the whole region of the nebula , we can write the h@xmath0 + [ ] flux in the net h@xmath0 + [ ] frame @xmath196}^\prime$ ] as @xmath197}^\prime = ( 1 - \beta ) f_{\rm h\alpha+[n ii ] } , \label{eq : f}\ ] ] where @xmath127 is the contamination factor in the @xmath22 frame , @xmath196}^\prime$ ] is the h@xmath0 + [ ] flux uncorrected for over - subtraction of line emission from the @xmath22 frame , and @xmath196}$ ] is that corrected for this effect . here
, we assumed that the filter transmissions are the same for all the emission lines , for simplicity .
we then derived the scaling factor @xmath198 in the following equation by subtracting the net h@xmath0 + [ ] frame multiplied by @xmath198 from the scaled @xmath22 frame so that the central butterfly nebula component disappears in the subtracted frame , @xmath199}^\prime = \beta\ ;   f_{\rm h\alpha+[n ii ] } \label{eq : c}.\ ] ] equations ( [ eq : f ] ) and ( [ eq : c ] ) give @xmath196 } = ( 1+c ) f_{\rm h\alpha+[n ii]}^\prime$ ] .
we found that the best scaling factor @xmath198 was 0.23 .
thus the correction factor for the h@xmath0 + [ ] flux of the net h@xmath0 + [ ] image is 1.23 .


entanglement is the key resource of the vast fields of quantum information processing , quantum computation , and quantum technology , for an introduction see e.g.  @xcite .
for example , applications of entangled states are those for quantum key distribution  @xcite , quantum dense coding  @xcite , and quantum teleportation  @xcite .
thus both , the identification and the quantification of entanglement , play a mayor role for future applications  @xcite .
the phenomenon entanglement is closely related to the superposition principle of quantum mechanics .
a pure separable state is represented by a product of states for both systems .
a general pure state is a superposition of factorizable states .
the minimal number of such global superpositions denotes the schmidt rank  @xcite .
a separable mixed quantum state is a convex combination of pure factorizable quantum states  @xcite .
the generalization of the schmidt rank to mixed quantum states delivers the schmidt number ( sn ) .
this generalization and the introduction of sn witnesses is given in  @xcite .
the experimental construction of states with a certain sn has been realized in  @xcite .
the sn of a mixed quantum state fulfills the axioms of an entanglement measure , cf .
more precisely , it is a convex roof measure as defined in  @xcite .
the identification of entanglement in terms of entanglement witnesses has been introduced in  @xcite . for a given witness an optimization
can be performed  @xcite .
recently , we proposed optimized , necessary and sufficient conditions for the detection of entanglement  @xcite .
note that the latter optimization and the optimization of an entanglement witness are inherently different .
our optimization procedure delivers so - called separability eigenvalue equations .
they resemble the well - known eigenvalue problem , but they include the factorization property of pure local quantum states .    in the present contribution we study the identification of the sn of a given quantum state .
we provide a method that deliver all optimized sn witnesses . the optimization yields equations which will be discussed in two forms .
these generalizations of the separability eigenvalue problem delivers similar equations for arbitrary mixed sn states .
our method will be compared with the well known spectral decomposition and the schmidt decomposition of quantum states .
properties of the solutions of these equations will be studied . with some fundamental examples ,
we generate some general classes of sn witnesses . we apply them to mixed quantum states .
the paper is structured as follows .
we motivate our method in sec .
[ sec : motiv ] . in sec .
[ sec : witness ] we reformulate the detection of the sn of a quantum state by witnesses in terms of arbitrary hermitian operators and optimized , necessary and sufficient conditions .
the optimization procedure will be discussed in sec .
[ sec : sneval ] , where we express the optimization problem in terms of a perturbed eigenvalue problem and discuss properties of these equations . in sec .
[ sec : sesol ] we solve these equations for a wide class of operators , including all one dimensional projectors and operators defined in continuous variables .
we apply the method in sec .
[ sec : ex ] to identify a sn greater than one and two , for the case of a phase diffused two - mode squeezed - vacuum state .
a summary and some conclusions are given in sec .
[ sec : sc ] .
let us consider the following experimental situation , cf . fig .
[ fig : exp ] .
we have a beam splitter with squeezed - vacuum states in both inputs .
these states have the same amount of squeezing , but in orthogonal quadratures .
the output is the two - mode squeezed - vacuum state @xmath0 , ( @xmath1 and @xmath2 ) , @xmath3 one output is disturbed in terms of a phase randomization , which equally randomizes the phase @xmath4 between @xmath5 and @xmath6 .
the scenario under consideration could be used for transferring one part of an entangled state through a noisy channel .
the sender keeps the other part of the state , e.g. in a delay line , such as an optical fiber .
the measured state is given by @xmath7 with the local unitary operation @xmath8 .
in general , this state is a mixed quantum state in continuous variable systems . due to the linear independence of the two - mode squeezed - vacuum states ,
the rank of the operator @xmath9 is , in general , also infinite .
, @xmath10 resemble the measurement of a , so far unknown , test operator @xmath11 .
the fluctuation of the distance between the beam splitter output and the detector @xmath12 yield the phase randomization @xmath13 .
the distance between the beam splitter output and the detector @xmath10 is fixed . , width=226 ]    the two - mode squeezed - vacuum state has been considered , for example , as a recource for quantum teleportation  @xcite , quantum dense coding  @xcite , and quantum memories  @xcite .
it has been shown that this state violates a continuous variable bell inequality  @xcite .
pertubations of the pure two - mode squeezed - vacuum state have been studied , such as phase and amplitude damping  @xcite , or noise due to the transmission in optical fibers  @xcite . here , we focus on the sn of the phase randomized two - mode squeezed - vacuum state @xmath14 .    in the case of zero phase randomization , we have an infinite sn for @xmath15 .
this state includes an entangled qubit ( @xmath16 ) , an entangled qutrit ( @xmath17 ) ,  , an entangled qudit ( @xmath18 ) , see eq  ( [ eq:1 ] ) .
the two - mode squeezed - vacuum state is a global superposition of infinitely many product states @xmath19 , which corresponds an infinite sn . in this case
, entanglement determines all the correlation between system @xmath20 and system  @xmath21 .
now let us consider the fully randomized state @xmath22 , @xmath23 this state is obviously separable .
thus , it does not include any entangled qudits , and the sn is one .
the subsystems @xmath20 and @xmath21 are only classically correlated .
some questions arise automatically . somewhere in between the two extreme cases , @xmath24 and @xmath22
, all the qudits disappear one after another , and the quantum correlations between @xmath20 and @xmath21 vanish .
we aim to answer questions such as : _ for which value @xmath25 the state becomes separable ? _ , or : _ which phase randomization @xmath25 still deliver a quantum state with , for example , an entangled qutrit ? _ with other words , we want to quantify the entanglement of the state @xmath14 under a certain phase randomization by the sn .
let us divide the problem into several sub - problems , which we solve in this manuscript :    [ [ conditions . ] ] conditions .
+ + + + + + + + + + +    we will derive general , necessary and sufficient conditions for the identification of a certain schmidt number @xmath26 ( sn-@xmath26 ) .
this method applies to all mixed quantum states .
these conditions are based on measuring a general observable @xmath11 , @xmath27 . for obtaining the sn of the state
, we relate @xmath27 to the maximal expectation value under all quantum states with a sn less or equal to @xmath26 , @xmath28 .    [ [ optimization . ] ] optimization .
+ + + + + + + + + + + + +    we derive equations in two different forms , to obtain the desired function @xmath29 .
these equations resemble the eigenvalue problem .
the global superposition property of sn-@xmath26 quantum states is encoded in these equations .
[ [ solution . ] ] solution .
+ + + + + + + + +    we analyze these new kinds of equations .
we find a close relation to the schmidt decomposition and the spectral decomposition of an operator .
afterward , we solve them for a large class of operators - including those for our considered example .
[ [ applying - the - method . ] ] applying the method .
+ + + + + + + + + + + + + + + + + + + +    we apply our method to the phase randomized two - mode squeezed - vacuum state .
we show how the squeezing and the phase randomization influence the quantum correlations between subsystems @xmath20 and @xmath21 .
we identify entanglement in general and the existence of a qutrit in the randomized state .
altogether , we derive a new mathematical method for the identification of the sn .
this method allows us to test up to which order entangled qudits are contained in an arbitrary bipartite quantum state .
these entangled qudits can be used for quantum information processing .
let us consider a bipartite quantum system which is given by compound hilbert space @xmath30 .
previously we have shown in , how entanglement in cantinuous variable systems can be identified in finite dimensional spaces  @xcite .
thus , it is not a restriction if we assume @xmath31 . in addition
let us denote the sets of linear operators acting on the hilbert space as @xmath32 , and the hermitian operators as @xmath33 .
statistical mixtures of pure factorizable states define the separable quantum states  @xcite .
following this idea , quantum states with a sn less or equal to @xmath26 are mixtures of pure states with a schmidt rank less or equal to @xmath26 .
pure states @xmath34 with a sn less or equal to @xmath26 , @xmath35 , are elements of the set @xmath36 and can be written as  @xcite @xmath37 with a local unitary operation @xmath38 , orthonormal @xmath39 , @xmath40 , and @xmath41 , and schmidt coefficients @xmath42 .
the mixed sn-@xmath26 states @xmath43 are elements of @xmath44 given as convex combinations of those pure ones , @xmath45 .
let us note that the maximal possible sn is @xmath46 , and the minimal one is @xmath47 for separable quantum states .
analogously to entanglement witnesses , a sn witnesses @xmath48 is a hermitian operator , with @xmath49 an operator fulfilling eq .
( [ eq : snwit ] ) is called sn-@xmath26 witness , @xmath50 .
the condition given in eq .
( [ eq : witwit ] ) will be considered at the end of our treatment . here
, let us only note that the set @xmath51 includes all sn-@xmath26 witnesses fulfilling eq .
( [ eq : witwit ] ) , and all positive semi - definite operators . here
, a sn witness @xmath48 is considered to be optimal , if there exits a sn-@xmath26 state @xmath52 with @xmath53 .      in this subsection
we aim to derive the conditions for the identification of sn-@xmath26 states .
it has been considered in  @xcite , that a quantum state has a sn greater than @xmath26 , if and only if @xmath54 the existence of such a witness is a consequence of the hahn - banach theorem . whereas , the general structure of sn-@xmath26 witnesses  elements of @xmath51  is not clear .
thus , it is an advantage to find conditions in terms of arbitrary hermitian operators with a well - known structure .    for this purpose
let us define the following function @xmath29 .
this function maps a hermitian operator to its maximal expectation value of all mixed and pure sn-@xmath26 states .
[ def : maxexp ] the maximal sn-@xmath26 expectation value is given by @xmath55    [ [ proof . ] ] proof .
+ + + + + +    this value exists due to the compactness of @xmath36 .
this value is the maximal expectation value of _ all _ sn-@xmath26 states :    * @xmath56 @xmath57 @xmath58 * @xmath59 , with @xmath60 and @xmath61 @xmath62 @xmath57 @xmath63    it follows @xmath64 .
@xmath65    in the following , the decomposition of all operators of the set @xmath51 in terms of arbitrary hermitian operators is considered .
we show that a certain form delivers a sn-@xmath26 witness , and all sn-@xmath26 witness can be written in this form .
this yields the identification of the sn .
let us consider an arbitrary hermitian operator @xmath66 , and a number @xmath67 .
we define the operator @xmath68 as @xmath69 with the bipartite identity operator @xmath70 .
this is a generalization of the construction of entanglement witnesses , cf .
. obviously holds for all sn-@xmath26 quantum states @xmath43 @xmath71 thus , @xmath68 is a sn-@xmath26 witness .    in the case
@xmath72 , we call this sn-@xmath26 witness optimized . for @xmath47 ,
our optimality condition is , in general , weaker then the definition in ref .
@xcite . to be precise
, our optimized witness , @xmath72 , is finer than the witness for @xmath73 .    the other way around we decompose all witnesses in this form .
let @xmath50 be a sn-@xmath26 witness .
using eq .
( [ eq : snwit ] ) multiplied with @xmath74 and the resulting choice @xmath75 , the operator @xmath76 delivers the decomposition , @xmath77 it follows that @xmath68 is a sn-@xmath26 witness , and any sn-@xmath26 witness can be written as given in eq .
( [ eq : snwitdec ] ) .
this means the set @xmath51 and the set which is defined by operators of the form in eq .
( [ eq : snwitdec ] ) are identical .    moreover , the decomposition in eq .
( [ eq : snwitdec ] )  and  ( [ eq : snpos ] ) proves us the fact that optimized witnesses are necessary and sufficient for the detection of a sn . summing up all these findings
, we can formulate optimized , necessary and sufficient conditions in terms of arbitrary operators .
[ theo : condition ] a quantum state @xmath78 has a sn greater than @xmath26 , if and only if there exists @xmath66 with @xmath79    [ [ proof.-1 ] ] proof .
+ + + + + +    the first part of the proof has been done above with the general structure of a witness in eq.([eq : snwitdec ] ) .
we use this general form of an optimized sn-@xmath26 witness @xmath80 .
in addition , we use the fact that a witness must exist : @xmath81 this condition means that the expectation value of @xmath78 exceeds the boundary given by all sn-@xmath26 states . the formulation of the theorem
could also be given by @xmath82 if we use the minimal value instead of the maximal value in definition  [ def : maxexp ] .
in conclusion , we have obtained sn-@xmath26 conditions in terms of arbitrary operators .
so far , only some special examples of sn-@xmath26 witnesses have been known , so that a general approach was unknown .
this shortcoming we have resolved in terms of the condition in theorem  [ theo : condition ] .
such a method is of some general interest , since the identification of the sn of a quantum state delivers the structure and the quantification of entanglement in the notion of global superpositions . in the following
, we will focus on the determination of @xmath29 , which is desired for a general characterization of entanglement through the sn .
now , we consider the function @xmath29 as given in definition  [ def : maxexp ] .
the value of @xmath28 denotes the maximal expectation value of a given operator @xmath11 under all sn-@xmath26 states .
first , we consider all extrema  not only the global maximum .
obviously the largest extrema is the desired value of @xmath29 .
the included optimization problem will lead us to equations which are of an algebraic nature .
this means we transform a optimization problem to an algebraic system of equations .
first of all , let us consider a weaker decomposition than the schmidt decomposition given in eq .
( [ eq : schmdec ] ) .
a pure quantum state @xmath34 has a sn less or equal to @xmath26 , if it can be written as @xmath83 with @xmath84 . here
, the product states @xmath85 are not necessarily orthonormal or linear independent
. some of them can also be the zero .
thus , @xmath34 might have a schmidt rank less than @xmath26 .
but in any case it follows that @xmath86 .    the optimization problem given in definition
[ def : maxexp ] is the following .
for a given @xmath66 , we want to find the maxima or minima  in general extrema  of @xmath87 . for a quantum state @xmath86
the expectation value is given as @xmath88 where @xmath89 denotes an optimal value .
we also make use of the normalization condition @xmath90    an optimization problem under a certain condition can be solved with the method of lagrangian multipliers @xmath91 .
this means for all @xmath92 holds @xmath93 using @xmath34 in the form of eq .
( [ ee : schmdecweak ] ) , we obtain from eqs .
( [ eq : lagrmult1 ] )  and  ( [ eq : lagrmult2 ] ) for all @xmath92 @xmath94 the lagrangian multiplier @xmath91 can be obtained by summing up eq .
( [ eq : lagrmult3 ] )  summation given by @xmath95  and using the condition @xmath96 in eq .
( [ eq : condnorm ] ) : @xmath97 this is already the algebraic conversion of the optimization problem .    further on ,
let us simplify and rewrite these equations .
therefore we may consider the following definitions @xmath98
\right)_{k , j=1,\dots , r},\\      & l_{\vec y}=\left ( { \rm tr}_b
l\left[\mathbb i\otimes|y_k\rangle\langle y_j|\right ]   \right)_{k , j=1,\dots , r}.\end{aligned}\ ] ] analogously to the operators @xmath99 and @xmath100 , we define the _ pseudo - metrics _ @xmath101 and @xmath102 . from eqs .
( [ eq : lagrmult3 ] )  and  ( [ eq : lagrmult4 ] ) we obtain the first set of equations for the sn-@xmath26 optimization problem , @xmath103 this form generalizes the separability eigenvalue equations , as derived in  @xcite .
in fact , for @xmath47 eqs .
( [ eq : sesnr1 ] )  and  ( [ eq : sesnr2 ] ) reduce to the separability eigenvalue equations .
[ def : rseprobl1 ] the @xmath26-schmidt eigenvalue problem is defined in form 1 as @xmath104 we use the abbreviation @xmath26-se for @xmath26-schmidt eigenvalue . the vector @xmath34 is the @xmath26-schmidt eigenvector ( @xmath26-se vector ) , and the value @xmath91 is the @xmath26-se for @xmath34 of @xmath11 .
the conversion of the optimization problem to its algebraic form is given in the following theorem which summarizes the above calculations .
the sn-@xmath26 state @xmath34 delivers an extremal expectation value @xmath105 , if and only if it solves the @xmath26-se problem defined in definition  [ def : rseprobl1 ] .
@xmath65    more sophisticated is the derivation of the second form of these equations .
we note that the sn-@xmath26 state in its weak form of eq .
( [ ee : schmdecweak ] ) must also have a _ strict _ schmidt decomposition , see eq .
( [ eq : schmdec ] ) .
obviously the optimization and the optimal value @xmath106 does not depend on the form of a decomposition of the pure state .
let us note that the @xmath26 in eq .
( [ ee : schmdecweak ] ) is not necessarily the schmidt rank @xmath107 of @xmath34 .
due to the independence of the decomposition and @xmath108 hence , we relabel formally @xmath109    now , let us reconsider which information is included in eqs .
( [ eq : lagrmult3 ] )  and  ( [ eq : lagrmult4 ] ) , and rewrite them as @xmath110 we use the notion @xmath111 and @xmath112 for the projection of the state @xmath34 on the @xmath40 and @xmath41 component , respectively .
note that @xmath42 .
thus , for each @xmath92 , we can multiply eq .
( [ eq : lagrmult11 ] ) with @xmath113 . equations  ( [ eq : lagrmult10 ] )  and  ( [ eq : lagrmult11 ] ) yield for @xmath114 @xmath115    on the other hand , let us consider the action of @xmath11 on the pure sn-@xmath26 state @xmath34 .
this delivers the ( not normalized ) result @xmath116 , @xmath117 with the decomposition of @xmath116 into a parallel component , @xmath118 , and an orthogonal one , @xmath119 . now on the subsystems @xmath20 and @xmath21 act @xmath40 and @xmath120 , respectively , @xmath121 we combine eqs .
( [ eq : lagrmult5 ] )    ( [ eq : lagrmult7 ] ) and obtain @xmath122 for all @xmath114 .
we conclude that @xmath119 is orthogonal to each @xmath123 for @xmath114 which leads to @xmath124 this means that @xmath119 is not only orthogonal to @xmath34 , but also to all linear combinations of @xmath125 of the schmidt decomposition of @xmath34 and @xmath114 .
we denote this property of eq .
( [ eq : vecchi ] ) as _ bi - orthogonality_.    from these considerations follows another algebraic conversion of the optimization problem as @xmath126 with a bi - orthogonal vector @xmath119 as it is given in eq .
( [ eq : vecchi ] ) .
we define the second algebraic form of the optimization problem for a state @xmath34 with maximally @xmath26 global superpositions of factorizable states .
[ def : rseprobl ] the @xmath26-schmidt eigenvalue problem is defined in form 2 as : @xmath126 with a bi - orthogonal perturbation @xmath119 .    due to the fact that definitions  [ def : rseprobl1 ]  and  [ def : rseprobl ] solve the same optimization problem , form 1 and form 2 of the @xmath26-se problem are equivalent .
the sn-@xmath26 state @xmath34 delivers an extremal expectation value @xmath105 , if and only if it solves the @xmath26-se problem defined in definition  [ def : rseprobl ] .
@xmath65      so far , we have derived two forms of the @xmath26-se problem . our initial intention was to determine the value of the function @xmath28 for a given linear operator @xmath11 .
the @xmath26-se @xmath91 is the extremal value of the function @xmath127 , @xmath128 under all extremal values the global maximum can be found as    @xmath129 .
@xmath65    in some cases @xmath11 can not be used for the identification of sn-@xmath26 states .
obviously holds that @xmath11 is not suitable in such a case , if and only if @xmath130 , cf .
( [ eq : snwit ] )  and  ( [ eq : witwit ] ) and @xmath131 denoting the maximal possible schmidt rank .
this means , no quantum state with an arbitrary sn exceeds the maximal expectation value of @xmath11 for sn-@xmath26 states .
thus , the test operator @xmath11 for the identification of the sn must have the following property .
the eigenspace of the largest eigenvalue does not include a vector with a schmidt rank less or equal to @xmath26 . for @xmath47
this is related to the range criterion for quantum states  @xcite .
it is clear that the bi - orthogonal perturbation @xmath119 has a schmidt decomposition as well . due to the bi - orthogonal property it can be given as @xmath132 now the @xmath26-se equation in definition  [ def : rseprobl ]
can be rewritten as @xmath133 thus , we can write a necessary condition for a solution @xmath34 of the @xmath26-se problem .    [ cor : commonsdec ] a @xmath26-se vector @xmath34 and the vector @xmath134 have a schmidt decomposition with the same local unitary operation @xmath38 , @xmath135    this condition is only necessary .
it would also be sufficient , if @xmath136 for @xmath137 .
but corollary  [ cor : commonsdec ] in its given form is already a quite strong restriction for finding @xmath26-se solutions .
we will see this later when we apply our method .
we consider the @xmath26-se problem for all possible @xmath138 .
in addition , let us consider the _
( ordinary ) _ eigenvalue problem of @xmath139 in linear algebra , @xmath140 with the eigenvector @xmath34 and the eigenvalue @xmath91 .
it is clear that in this case there is a trivial perturbation @xmath141 .
thus , eigenvalues and eigenvectors , @xmath91 and @xmath34 , are also solutions of the @xmath26-se problem for @xmath142 . the other way around , we conclude for @xmath143 that no bi - orthogonal perturbation can remain .
therefore , we obtain in this case that the @xmath131-se problem and the eigenvalue problem are identical .    now , let us consider local invertible operations .
these are operations which can be written as @xmath144 , with @xmath145 and @xmath146 invertible operations , @xmath147 and @xmath148 .
these operations can not change the amount of entanglement , given by the sn , for any quantum states  @xcite .
now we transform the @xmath26-se problem of the operator @xmath11 given in definition  [ def : rseprobl ] to the @xmath26-se of an operator @xmath149 @xmath150 this can be done by multiplying the @xmath26-se problem of @xmath11 with @xmath144 , @xmath151 we obtain that the @xmath26-se values remain , and the vectors are transformed by @xmath144 , @xmath152    [ cor : loctrans ] let @xmath149 be a locally transformed operator @xmath153 , @xmath91 is a @xmath26-se value and @xmath34 a is @xmath26-se vector of @xmath11 .
it follows that @xmath154 are @xmath26-se value and @xmath26-se vector of @xmath149 , respectively .
[ [ proof.-2 ] ] proof .
+ + + + + +    any local invertible map can be decomposed in terms of unitary and diagonal matrices , e.g .. by the singular value decomposition .
first let us consider the local unitary @xmath38 , @xmath155 .
we use @xmath156 with the @xmath26-se problem @xmath157 it follows from eq .
( [ eq:41 ] ) that @xmath158 in addition we consider local , diagonal , and invertible maps , @xmath159 and we conclude for the transformed @xmath26-se problem @xmath160 the schmidt coefficients change , @xmath161 , but the @xmath26-se value remains invariant . @xmath65    for the _ ( ordinary ) _ eigenvalue problem in linear algebra an invertible  in general global  transformation @xmath162 can be found for the diagonalization of the matrix . in general , a matrix @xmath11 can be transformed into another one by @xmath163 .
this transformation changes the eigenvectors in the same way , and the eigenvalues remain invariant . for the @xmath26-se problem
we have a related situation which , in addition , reflects the global superposition property of sn @xmath26 states .
let us consider a special case of local invertible operations , cf .
proof of corollary  [ cor : loctrans ] .
the same considerations are obviously true for @xmath164 , with @xmath165 and @xmath166 being unitary operations .
consequently , the @xmath26-se problem is invariant under local unitary , which is of importance for the quantification of entanglement , cf .
e.g.  @xcite .
so far we have considered general properties of the @xmath26-se problem .
we have compared our method with the eigenvalue problem in linear algebra .
now let us apply our method to some examples .
starting from low rank matrices , we explain how the @xmath26-se problem can also be solved for higher rank operators .
we consider an approach connecting the spectral decomposition with the @xmath26-se problem .
our sn-@xmath26 conditions will be tested for some examples in continuous variable systems .      for simplicity
let us here and in the following restrict to hilbert spaces with @xmath167 .
we consider a one dimensional projection @xmath168 , with @xmath169 , @xmath170 @xmath171 , and the normalization @xmath172 .
in addition , let us assume that the coefficients @xmath173 are ordered as @xmath174 .
it is also useful to define @xmath175 .
we aim to solve the @xmath26-se problem for @xmath11 .
first , there are a lot of trivial solutions given by orthogonal eigenvectors @xmath176 an example is @xmath177 for @xmath178 .
the non - trivial solutions , @xmath179 , are more interesting for obtaining @xmath28 . applying corollary  [ cor : commonsdec ]
we conclude that we only need to consider states @xmath34 with a decomposition as @xmath180 with @xmath181 .
the schmidt rank of @xmath34 is the number of coefficients with @xmath182 . for a schmidt rank less than @xmath131
some of the @xmath183 have to be zero .
thus , we obtain by applying the projection , @xmath184 , @xmath185 a closer look on eq .
( [ eq : coeffproj ] ) delivers all the possible solutions : @xmath186 with a constant number @xmath187 .
the bi - orthogonality property between @xmath34 and @xmath119 delivers that from @xmath182 follows @xmath188 .    for the case
@xmath178 we obtain all solutions as @xmath189 for @xmath190 and @xmath191 .
the maximal @xmath192-se vector and value are @xmath193    in conclusion , every @xmath34  with @xmath26 schmidt coefficients @xmath194  and the same schmidt decomposition as @xmath195 solves the @xmath26-se problem .
the @xmath26-se value is given as @xmath196 .
using the schmidt coefficients with the largest absolute value , we obtain @xmath197 .
for example , let us use @xmath198 and the projection operator @xmath11 constructed from @xmath199 for @xmath200 ( note that the index starts from 0 ) .
this is the two - mode squeezed - vacuum state . in this case
we obtain the maximal solution of the @xmath26-se problem as @xmath201 for the quantum state @xmath15 ( @xmath202 ) we can identify the entanglement with an arbitrary sn by @xmath203 this optimal violation , @xmath204 , does not depend on the amount of squeezing @xmath205 . even in the case of a minimal
squeezing , @xmath206 , the condition eq .
( [ eq : tmsvcond ] ) proves that the state @xmath207 has a sn larger than any @xmath26 .
another example could be @xmath208 and equally distributed schmidt coefficients , @xmath209 in this case we obtain the maximal solution of the @xmath26-se problem for @xmath168 as @xmath210 in conclusion , we can formulate the following necessary criteria .
a quantum state @xmath78 has a schmidt rank greater than @xmath26 if @xmath211 with @xmath212 ( @xmath137 ) being the @xmath26 largest schmidt coefficients of @xmath195 .
an arbitrary hermitian operator @xmath11 can be decomposed in terms of one - dimensional projectors . in the following ,
we generalize our results for one - dimensional projections to classes of higher rank operators .
it is of advantage to consider all one dimensional projections as a local transformation of a given state with maximal schmidt rank .
we start from pure vectors @xmath34 and
@xmath213 @xmath214 the linear map @xmath215 yields the desired property @xmath216 note that the singular value decomposition of @xmath217 resembles the schmidt decomposition of @xmath34  @xcite .
this rewriting delivers that the scalar product of two states @xmath34 and @xmath218 is the scalar product of the corresponding matrices @xmath219 a general hermitian operator @xmath66 can be decomposed as @xmath220 with @xmath221 .
an example would be the spectral decomposition of @xmath11 with the eigenvalues @xmath222 and the eigenvectors @xmath223 .
we apply @xmath11 on a given state @xmath224 which yields @xmath225\right ) ( m_k\otimes\mathbb i)|\phi\rangle\\      = & \left(\sum_k l_k { \rm tr}\left [ m_k^\dagger m_\psi\right ] m_k\right)\otimes\mathbb i|\phi\rangle.\end{aligned}\ ] ] from corollary  [ cor : commonsdec ] follows that @xmath217 and the resulting matrix @xmath226 m_k\right)$ ] have the same singular value ( schmidt ) decomposition , if @xmath227 is a @xmath26-se vector .
note that the schmidt rank of the state @xmath34 is the rank of the corresponding matrix , @xmath228 .
let us consider an example .
may the matrices @xmath229 have a decomposition as @xmath230 , with unitary operators @xmath165 and @xmath166 and a diagonal matrix @xmath231 with complex components .
we have seen , that unitary operations basically do not effect the @xmath26-se problem , cf .
corollary  [ cor : loctrans ] .
hence , we choose @xmath232 .
the effect of this example is that the mapping delivers @xmath233 d_k\right),\end{aligned}\ ] ] which is already given as a diagonal matrix .
in such a case it follows , that the vector @xmath234 is given in schmidt decomposition and we can treat this case analogously to those given in section  [ ss : odp ] for low rank operators .
we may choose a decomposition @xmath235 ( with @xmath236 and @xmath237 ) , and the operator @xmath11 as @xmath238 the new parameters @xmath239 are given by @xmath240 it is worth to note that this operator @xmath11 is symmetric with respect to exchanging the quantum systems .    for obtaining the nontrivial solutions , @xmath179 , we consider @xmath241 .
this yields @xmath242 due to the fact that some @xmath183 can be zero , we obtain the @xmath26-se solutions by neglecting some rows and the corresponding columns , cf .
( [ eq : schmcoeffsevec ] ) .
these rows are given by the schmidt coefficients of @xmath34 which are zero , @xmath243 .
this means we restrict to the rows and columns with the index @xmath244 , and solve the ( ordinary ) eigenvalue problem of the resulting operator , @xmath245    the simplest case is @xmath47 .
we have to stroke out @xmath246 rows and the corresponding columns .
it immediately yields @xmath247 the maximal expectation value for all possible choices of @xmath248 delivers the function @xmath249 as @xmath250    in the case @xmath178 , we have to consider all principal @xmath251 sub - matrices of @xmath239 with @xmath252 .
we calculate the maximal eigenvalue , and obtain the function @xmath253 as @xmath254 denoting the maximal eigenvalue of all principal @xmath251 sub - matrices .
analogously , the method can be applied to sn-@xmath26 states by finding the maximal eigenvalue of all @xmath255 principal sub - matrices of @xmath239 .
let us come back to our initially considered state @xmath14 , cf .
( [ eq : rhovarphi ] ) .
this state was generated by a local phase randomization of a two - mode squeezed - vacuum state , cf .
[ fig : exp ] . for witnessing the sn , we already considered the example of a projection defined by a pure two - mode squeezed - vacuum state in sec .
[ ss : odp ] .
this operator was suitable for the detection that the state @xmath207 , @xmath24 , has an infinite sn .
however , we already observed that a total randomization , @xmath22 , delivers a separable state @xmath256 , cf .
eq  ( [ eq : tmsvstate ] ) .
now , we apply our considerations to the general phase randomized ( mixed ) state @xmath14 to answer our initial questions .
namely , for which values of @xmath25 the entanglement survives ( sn @xmath257 ) , and which values of @xmath25 guarantee that the state contains an entangled qutrit ( sn @xmath258 ) ?
a phase randomization in a single mode state has been experimentally realized , and the non - classicality of such a state has been verified  @xcite . here , we consider a similar situation , for two modes , where the desired aspect of non - classicality is the number of global superpositions between these modes .
the method presented in the following also applies to a non - equally distributed phase randomization in both subsystems @xmath20 and @xmath21 , and an additional amplitude randomization of @xmath259 , as it occurs in the turbulent atmosphere  @xcite .    here , we chose a test operator @xmath11 of formally the same form as the state @xmath14 in eq .
( [ eq : rhovarphi ] ) , @xmath260 with @xmath1 and @xmath2 . for simplicity
, we have assumed an equally distributed randomization for the angle , @xmath261 .
therefore , the corresponding parameters @xmath239 read as @xmath262\right)}{\delta\varphi[m - n]}.\\      = & ( 1-\epsilon^2)\epsilon^{m+n}{\rm sinc}\left(\delta\varphi[m - n]\right).\end{aligned}\ ] ]    first , we may consider a fully randomized and equally distributed phase for the test operator , @xmath22 , @xmath263 the spectral decomposition of @xmath11 is given in factorizable states .
thus , @xmath264 , and the phase insensitive operator can not detect any entanglement .
we already observed this result for the corresponding state @xmath256 which is separable .
now we consider a given partially phase randomized operator , @xmath11 .
we obtain from eq .
( [ eq : sn1f1 ] ) for separable states ( sn @xmath47 ) a maximal expectation value of @xmath265 the maximal expectation value for two - qubit states ( sn @xmath178 ) , cf .
( [ eq : sn2f2 ] ) , yields @xmath266 ^ 2 + 4\epsilon^{2m+2n}{\rm sinc}^2\left(\delta\varphi[m - n]\right )      } \right ) \intertext{using $ n = m+k$ and $ k\geq1 $ we obtain }      \nonumber f_2(l)=&\left.\frac{1-\epsilon^2}{2}\max_{m\in\mathbb n}\epsilon^{2m}\max_{k\geq1}\right(1+\epsilon^{2k}\\      & \left.+\sqrt{[1-\epsilon^{2k}]^2 + 4\epsilon^{2k}{\rm sinc}^2\left(\delta\varphi[m - n]\right)}\right).\end{aligned}\ ] ] the maximum over @xmath267 yields @xmath268 .
the maximum over @xmath92 has to be calculated numerically .
note that @xmath269 , and therefore @xmath11 can detect entanglement .
now let us consider the state @xmath14 .
the expectation value @xmath270 of the observable @xmath11 for the state @xmath14 reads as @xmath271)\\      \nonumber = & ( 1-\epsilon^2)^2\sum_{m=0}^\infty\epsilon^{4m}\left(1 + 2\sum_{k=1}^\infty \epsilon^{2k}{\rm sinc}^2(\delta\varphi k)\right)\\      = & \frac{1-\epsilon^2}{1+\epsilon^2}\left(1 + 2\sum_{k=1}^\infty \epsilon^{2k}{\rm sinc}^2(\delta\varphi k)\right).\end{aligned}\ ] ] again , this can be calculated numerically .
we consider the case @xmath272 , which corresponds to a moderate squeezing of about @xmath273  @xcite .
the entanglement condition for the choice of an operator @xmath274 can be formulated as @xmath275 this means a positive number is a verification of a sn greater than @xmath26 .    .
in addition the state is manipulated by an equally distributed phase randomization with the interval of the length @xmath276 .
the functions are normalized to the maximal value .
the entanglement can be verified until @xmath277 ( blue solid function ) .
we also identified a greater schmidt rank than 2 for phase diffusion of @xmath278 ( red dashed function ) .
, width=321 ]    in fig  [ fig:2msvac ] , the identification of entanglement for angles of an equal distribution in phase for @xmath279 $ ] is given .
the values are obtained by numerical calculation for the matrix components @xmath239 with @xmath280 .
this means that the given positivity , i.e. the entanglement , can only increase for indices @xmath281 .
let us stress that we choose an arbitrary operator for the detection of entanglement and a sn greater than 2 .
there may exists operators which deliver entanglement also for randomization above the limitations of the chosen operator @xmath11 .
in addition , it is worth to note that a higher squeezing of the state delivers a higher sensitivity . for a realistic @xmath282 ( @xmath283 ) squeezing of the input states , cf .
@xcite , we can identify entanglement up to a phase diffusion of @xmath284 and a sn greater than 2 close to @xmath285 , cf . fig .
[ fig:2msvac2 ] .
here , the chosen test operator @xmath11 is given as @xmath286\right)|m , m\rangle\langle n , n|.\end{aligned}\ ] ] this means that the state @xmath14 is entangled even for a phase randomization up to to the total phase randomization .
the other interesting point is that the state surely contains a qutrit ( sn greater than 2 ) for a randomization of more then @xmath287 in the channel of the receiver .
we note that for such a phase randomization , the squeezing of the single - mode squeezed - vacuum state would vanish .    , and the test operator in eq .
( [ eq:82 ] ) . in the numerical calculation
, the entanglement can be verified for a phase randomization up to @xmath288 ( blue solid function ) .
we also identified a greater sn than 2 for phase diffusion up to @xmath289 ( red dashed function ) .
, width=321 ]
we have derived optimized , necessary and sufficient conditions for the detection of the schmidt number of an arbitrary bipartite quantum state .
these conditions have been formulated in terms of arbitrary hermitian test operators as measurable conditions .
we have shown that the optimization problem leads to the @xmath26-schmidt - eigenvalue equations .
we discussed the properties of the solutions and consequences of these equations in connection with entanglement and its quantification via the schmidt number .
for example , we have shown which operators can be used for the identification of the schmidt number , and we have considered the relation to the eigenvalue problem .
we have solved these equations for a wide class of operators , namely for all one dimensional projections and classes of higher rank operators in finite and infinite dimensional systems .
we have shown a direct identification of the schmidt number for all pure states with our condition .
we also applied our method to the identification of the schmidt number of mixed quantum states , including those with an infinite rank and infinite schmidt rank in continuous variable systems .
examples of test operators are considered , which identify the schmidt number for a broad class of mixed states . to illustrate the method , we have studied the example of a phase diffused two - mode squeezed - vacuum state for a squeezing of @xmath290 and @xmath282 .
we have identified entanglement ( schmidt number greater than one ) , and we could identify entangled pairs of qutrits in the state ( schmidt number greater than two ) for different values of the phase randomization .
moreover , the influence of the amount of squeezing and the randomization has been considered .
we are grateful to t. kiesel for useful discussions .
we are also grateful to the anonymous referee for helpful comments .
this work was supported by the deutsche forschungsgemeinschaft through sfb 652 .
52 m. a. nielsen and i. l. chuang , _ quantum computation and quantum information _
( cambridge university press , cambridge , 2000 ) .
r. horodecki , p. horodecki , m. horodecki , and k. horodecki , rev .
phys . * 81 * , 865 ( 2009 ) .
a. k. ekert , phys .
lett . * 67 * , 661 ( 1991 ) . c. h. bennett and s. j. wiesner , phys .
lett . * 69 * , 2881 ( 1992 ) . c. h. bennett , g. brassard , c. crepeau , r. jozsa , a. peres , and w. k. wootters , phys .
* 70 * , 1895 ( 1993 ) .
o. ghne and g. tth , physics reports * 474 * , 1 ( 2009 ) .
r. f. werner , phys .
a * 40 * , 4277 ( 1989 ) .
a. sanpera , d. bru , and m. lewenstein , phys .
a * 63 * , 050301(r ) ( 2001 ) .
b. m. terhal and p. horodecki , phys .
a * 61 * , 040301(r ) ( 2000 ) .
d. bru , j. i. cirac , p. horodecki , f. hulpke , b. kraus , m. lewenstein , and a. sanpera , j. mod . opt . * 49 * , 1399 - 1418 ( 2002 ) .
a. b. uren , r. k. erdmann , m. de la cruz - gutierrez , and i. a. walmsley , phys .
lett . * 97 * , 223602 ( 2006 ) .
v. vedral , m. b. plenio , m. a. rippin , and p. l. knight , phys .
lett . * 78 * , 2275 ( 1997 ) .
g. vidal , j. mod . opt . * 47 * , 355 ( 2000 ) .
v. vedral and m. b. plenio , phys .
a * 57 * , 1619 ( 1998 ) . c. h. bennett , d. p. divincenzo , j. a. smolin , and w. k. wootters , phys .
a * 54 * , 3824 ( 1996 ) .
a. uhlmann , open sys .
* 5 * , 209 ( 1998 ) .
m. horodecki , p. horodecki and r. horodecki , phys .
a * 223 * , 1 ( 1996 ) .
m. lewenstein , b. kraus , j. i. cirac and p. horodecki , phys .
a * 62 * , 052310 ( 2000 ) .
j. sperling and w. vogel , phys .
a * 79 * , 022318 ( 2009 ) .
p. t. cochrane , g. j. milburn , and w. j. munro , phys .
a * 62 * , 062307 ( 2000 ) .
j. clausen , t. opatrny , and d .-
welsch , phys .
. a * * 62 , 042308 ( 2000 ) .
j. mizuno , k. wakui , a. furusawa , and m. sasaki , phys .
a 71 , 012304 ( 2005 ) .
k. jensen , w. wasilewski , h. krauter , t. fernholz , b. m. nielsen , m. owari , m. b. plenio , a. serafini , m. m. wolf , and e. s. polzik , nature physics * 7 * , 1316 ( 2011 ) .
b . chen , j .- w .
pan , g. hou , and y .- d .
zhang , phys .
lett . * 88 * , 040406 ( 2002 ) .
t. hiroshima , phys .
a * 63 * , 022305 ( 2001 ) .
l. knll , s. scheel , e. schmidt , d .-
welsch , and a. v. chizhov , phys .
a * 59 * , 4716 ( 1999 ) .
j. sperling and w. vogel , phys .
a * 79 * , 042337 ( 2009 ) .
g. toth , rhys .
a * 71 * , 010301(r ) ( 2005 ) .
p. horodecki , phys .
lett . a * 232 * , 233 ( 1997 )
. j. m. leinaas , j. myrheim , and e. ovrum , phys .
a * 74 * , 012313 ( 2006 ) .
j. sperling and w. vogel , arxiv:0908.3974 [ quant - ph ] .
the value of @xmath259 resembles the parameter of the unitary squeezing operator , @xmath291 .
the variance of the squeezed quadrature is given by @xmath292 .
the squeezing in @xmath293 is given by @xmath294 .
+ w. vogel and d .-
welsch , _ quantum optics _ , 3rd ed .
( wiley - vch , warnheim , 2006 ) .
t. kiesel , w. vogel , b. hage , j. diguglielmo , a. samblowski , and r. schnabel , phys .
a * 79 * , 022122 ( 2009 ) .
a. a. semenov and w. vogel , phys .
a * 81 * , 023835 ( 2010 ) .
h. vahlbruch , m. mehmet , s. chelkowski , b. hage , a. franzen , n. lastzka , s. gossler , k. danzmann , and r. schnabel , phys .
lett . * 100 * , 033602 ( 2008 ) .