the physical mechanisms to set the mass ratios of binary stars is of fundamental astrophysical importance , since more than half of main - sequence , pre - main sequence , and even protostars having masses comparable to the sun are members of binary systems @xcite .
binary stars with the primary masses comparable to the solar mass exhibit a wide , essentially flat distribution of the mass ratios ( @xmath12 ) from @xmath130.1 to 1 @xcite .
the physical processes to reproduce such a flat distribution of the binary mass ratios are , however , still controversial .
protostellar binary systems , precursors of main - sequence binary systems , are often surrounded by disks of molecular gas and dusts , `` circumbinary disks '' ( hereafter cbds ) @xcite .
the radii of these cbds range from @xmath9100 to @xmath7500 au , and these cbds often show inner emission depressions and ringlike structures .
mass accretion from these cbds onto the central protostellar binaries has been considered as one of the key physical processes to set the binary mass ratios .
recent high - resolution observations of the cbds have indeed detected gas flows from cbds to the inner protostellar binaries , while the spatial resolutions are not high enough to disentangle the accretions onto the primaries and secondaries @xcite .
theoretically , bate & bonnell ( 1997 ) and bate ( 2000 ) have shown that from their smoothed particle hydrodynamic ( sph ) simulations the majority of circumbinary materials accrete onto the secondary .
high - resolution grid - based simulations by ochi et al .
( 2005 ) and hanawa et al .
( 2010 ) have shown the opposite results , that is , the primary accretes more than the secondary .
ochi et al .
( 2005 ) and hanawa et al .
( 2010 ) argued that high - resolution simulations to sufficiently resolve the gas motions near the l2 and l3 lagrangian points are essential to properly trace the mass accretion , since accreting gas flows into the binary through these lagrangian points . on the other hand
, latest sph simulations have demonstrated that mass accretion in cbds is a function of gas temperatures , and in the case of the high temperature the primary accretes more than the secondary , and vice versa @xcite .
young et al . ( 2015 ) argued that the difference between the above results can simply be explained by the difference of the adopted gas temperatures in the models .
it has been still controversial to construct a consistent , unified model of mass accretion from cbds onto protostellar binaries .
cbds in protostellar binary systems are considered to be embedded in larger - scale protostellar envelopes , which can replenish the central cbds with fresh materials @xcite .
since the observed masses of cbds around protostellar binaries ( @xmath90.1 @xmath14 ) are much smaller than the masses of the central protostellar binaries @xcite , further supply of materials from the protostellar envelopes to the cbds is essential to significantly change the mass and mass ratios of the protostellar binaries .
if a protostellar envelope keeps supplying materials to the central cbd , further growth of the mass and the mass ratio of the protostellar binary will be expected . on the other hand ,
if the most of the envelope materials are being dissipated , the present mass and mass ratio of the protostellar binary will be close to the final values .
such dissipations of protostellar envelopes through the interactions with the associated outflows or stellar winds have been investigated observationally @xcite .
thus , in addition to physics of accretion from cbds onto the protostellar binaries , mass replenishment from the surrounding protostellar envelopes , and the connection between protostellar envelopes and cbds , must be taken into account .    in the present paper ,
we focus on a protostellar envelope harboring the protostellar binary system l1551 ne .
l1551 ne is a young class i protostellar binary ( @xmath15=91 k , @xmath16=4.2 @xmath17 ; froebrich 2005 ) located to @xmath72@xmath185 northeast of another , brightest protostellar binary l1551 irs 5 in the l1551 region @xcite .
the protostellar binary consists of the northwestern source named `` source a '' , and the southeastern source `` source b '' , with the projected separation of @xmath770 au at a position angle of 300@xmath19 @xcite .
our previous sma observations of l1551 ne have identified a @xmath6300 au scale cbd in keplerian rotation with a central stellar mass of 0.8 @xmath14 @xcite , plus a possible outer infalling component @xcite .
our subsequent alma cycle 0 observation of l1551 ne @xcite , at a spatial resolution that is @xmath71.6 times higher ( in beam area ) and a sensitivity that is @xmath76 times better ( in brightness temperature ) than those attained in our previous sma observations , unveiled substructures of the cbd in the 0.9-mm dust - continuum emission .
the revealed substructures are consistent with our 3-d adaptive mesh refinement ( amr ) hydrodynamic simulation @xcite for the presence of two spiral arms driven by gravitational torques from the central binary system .
the alma data of l1551 ne in the c@xmath0o ( 32 ) line also exhibits the deviations from the keplerian motion in the cbd , consistent with our amr simulation that gravitational torques impart angular momentum along the spiral arms ( driving material outwards ) and extract angular momentum between the spiral arms ( driving infall ) . our theoretical model which reproduces the observed features of the cbd around l1551 ne predicts that secondary accretes more than the primary . since our series of the sma and alma observations of l1551
ne have revealed the cbd + protostellar binary system , the next question is the connection from the protostellar envelope to the cbd , and the mass replenishment from the envelope to the cbd .
we conducted single - dish mapping observations of the protostellar envelope surrounding the protostellar binary system of l1551 ne in the c@xmath0o ( @xmath1=32 ) , @xmath2co ( @xmath1=32 ) , cs ( @xmath1=76 ) , and the so ( @xmath3=7@xmath46@xmath5 ) lines with the atacama submillimeter telescope experiment ( aste ) 10-m telescope .
the aste mapping data in the c@xmath0o line are combined with our previous sma data in the same line , which have a larger field of view and a shorter spacing than those of the alma data and thus more suitable to study the connection between the extended protostellar envelope and the cbd . in this paper
we report these aste+sma results of the protostellar envelope in l1551 ne .
the structure of this paper is as follows . in section 2 , our new aste observations , as well as the process to combine
the aste mapping data with the sma data in the c@xmath0o line , are described . in section 3 ,
the results of the aste mapping observations and the combination of the aste + sma data are presented . in section 4 ,
gas motions and physical properties of the identified gas components are estimated . in section 5 , we discuss implications of the observed features in the context of growth of protostellar binaries , and concise summary is written in section 6 .
we performed observations of l1551 ne in the c@xmath0o ( @xmath1=32 ; 329.3305453 ghz ) , @xmath2co ( @xmath1=32 ; 330.587965 ghz ) , so ( @xmath3=7@xmath46@xmath5 ; 340.71416 ghz ) , and the cs ( @xmath1=76 ; 342.882857 ghz ) lines with the aste 10-m telescope on 2013 september 23 - 26 and october 4 .
remote observations were performed from the aste operation room of naoj at mitaka , japan , using the network observation system n - cosmos3 developed by naoj @xcite .
a cartridge - type , double - sideband 350 ghz receiver with an if frequency range of 4.5 to 7 ghz ( cats345 ) mounted on aste @xcite was used , and all the four lines were observed simultaneously , except that on september 23 only the so and cs lines were observed due to the mistake of the instrumental setting .
we used only the data with the dsb system noise temperature ranging @xmath7200 - 800 k , and the data with the higher noise temperatures are excluded .
the telescope pointing was checked every @xmath71.5 - 2 hours by five - point co ( @xmath1 = 32 ) observations of a late - type star , nml - tau , and was found to be better than @xmath72@xmath11 . as a standard source we also observed orion kl @xcite and l1551 irs 5 @xcite , and confirmed that the relative intensity was consistent within @xmath730@xmath20 with a main beam efficiency of 0.6 .
no further correction for the sideband ratio was performed .
hereafter we show the observed line intensities in the unit of @xmath21 .
the observations consisted of two parts ; one - point , deep observations toward the center of l1551 ne ( @xmath22,@xmath23)= ( 04@xmath2431@xmath2544@xmath2647 , 180832@xmath272 ) , which is the field center of our previous sma observations and matches approximately the positions of the protostellar binary @xcite , and mapping observations around l1551 ne .
both observations were conducted in the position - switching mode .
the so and cs line data toward the center of l1551 ne taken on september 23 , when the other line data were not taken due to the mistake of the instrumental setting , are included in the deep integrations .
the resultant total on - source integration times of the so and cs data toward the center of l1551 ne is 1620 sec , while those of the @xmath2co and c@xmath0o data 1220 sec .
the rms noise levels of the @xmath2co , c@xmath0o , so , and the cs spectra toward the center of l1551 ne are 0.055 k , 0.062 k , 0.025 k , and 0.029 k , respectively .
the mapping observations were conducted with a grid spacing of 10@xmath11 ( nyquist sampling ) , and a typical on - source integration time per point of @xmath780 seconds .
the map center is set to be the position of the deep integrations .
the spatial and spectral resolutions and the noise level of the mapping observations in the c@xmath0o line are summarized in table 1 .
we combined the aste image cube in the c@xmath0o line with the interferometric c@xmath0o data of l1551 ne taken with the sma in its subcompact and compact configurations @xcite , adopting the method described by takakuwa et al .
( 2007b ) .
details of the sma observations are described by takakuwa et al .
the conversion factor from @xmath21 ( k ) to s ( jy beam@xmath8 ) was derived to be 48.5 as @xmath28 where @xmath29 is the boltzmann constant , @xmath30 is the wavelength , and @xmath31 is the solid angle of the aste beam (= 23@xmath11 ) . the spatial and spectral resolutions and the noise level of the combined sma+aste image cube in the c@xmath0o line are summarized in table 1 .
figure [ astesp ] shows the observed aste spectra of the c@xmath0o ( @xmath1=32 ) , @xmath2co ( @xmath1=32 ) , cs ( @xmath1=76 ) , and the so ( @xmath3=7@xmath46@xmath5 ) lines toward the protobinary of l1551 ne .
the c@xmath0o spectrum appears to consist of narrow and wide components .
two - component gaussian fitting to the c@xmath0o spectrum shows that the peak brightness temperature , line width , and the centroid velocity of the narrow component are 2.4 k , 0.72 km s@xmath8 , and 6.68 km s@xmath8 , and those of the wide component 1.1 k , 2.24 km s@xmath8 , and 6.99 km s@xmath8 , respectively ( light lines in figure [ astesp ] ) .
previous cso observations also found that the c@xmath0o ( 32 ) spectrum toward l1551 ne can be decomposed of two gaussian components with similar line widths and centroid velocities @xcite , in particular , the centroid velocity of the broad component is similar to the centroid velocity of the cbd derived from our previous sma observations (= 6.9 km s@xmath8 ) @xcite .
as shown in the next subsection , the broad c@xmath0o component is present only toward the protobinary position .
these results indicate that the broader component in the aste c@xmath0o spectrum likely originates from the central cbd . on the other hand ,
the narrower c@xmath0o component traces a distinct component at a slightly blueshifted ( @xmath70.2 km s@xmath8 ) centroid velocity .
as shown below , this narrow component traces the extended envelope component around l1551 ne . hereafter in this paper , @xmath32 = 6.7 km s@xmath8 is adopted as the systemic velocity of the extended envelope component , while @xmath32 = 6.9 km s@xmath8 as the systemic velocity of the cbd in l1551 ne .
the @xmath2co ( @xmath1=32 ) spectrum shows an absorption dip at around the centroid velocity of the narrow component of the c@xmath0o spectrum , and a broad redshifted wing up to @xmath32 @xmath713 km s@xmath8 . while a part of the redshifted wing in the @xmath2co spectrum likely originates from the cbd as in the case of the c@xmath0o spectrum , mapping observations of the @xmath2co line demonstrate that the redshifted @xmath2co emission is most likely arising from the outflow component driven from l1551 irs 5 , as will be discussed below . with a deeper integration a much better cs ( @xmath1=76 ) spectrum of l1551 ne than our previous one @xcite
is obtained .
the cs spectrum shows a non - gaussian , flat - top spectral shape . whereas the so ( @xmath3=7@xmath46@xmath5 ) line is weak compared to the other lines , the peak intensity ( @xmath70.14 k ) is detected above 5.5@xmath33 .
figure [ astemom0 ] shows total integrated intensity maps of the four molecular lines observed with aste .
the c@xmath0o ( 32 ) emission exhibits condensed and diffuse components , the former of which is centered on the protobinary position .
two - component , two - dimensional gaussian fitting to the c@xmath0o total integrated map shows that the c@xmath0o emission can be decomposed of a central component with a size of @xmath75500 au @xmath34 2600 au ( p.a.=-46@xmath19 ) and a diffuse component with a size of @xmath721000 au .
emission peaks toward the protostellar position are also seen in the cs ( 76 ) and so ( 7@xmath46@xmath5 ) lines . on the other hand ,
the @xmath2co ( 32 ) peak emission does not coincide with the protostellar position but is offset at @xmath720@xmath11 west from the protostar .
the @xmath2co emission shows an elongated feature toward the southwest direction , which is also evident in the c@xmath0o and cs emission .
figure [ aste13co ] shows the aste velocity channel maps of the @xmath2co line at a velocity interval of 0.455 km s@xmath8 . in the blueshifted velocity range ( @xmath32=4.62  5.99 km s@xmath8 ) , an emission component to the northwest of the protostar and that centered on the protostellar position are seen . around the systemic velocity ( 6.45  6.90 km s@xmath8 )
the @xmath2co line shows less intensity contrast , suggesting that the @xmath2co emission traces the overall cloud component . in the redshifted velocity of 7.36  9.18 km s@xmath8 , two emission components , one located to @xmath720@xmath11 west of the protostar , and the other to @xmath740@xmath11 southwest ,
while the southwestern component diminishes at velocities higher than 9.63 km s@xmath8 , the western component is seen until @xmath712.36 km s@xmath8 .
the location and velocity of this high - velocity redshifted component is consistent with those of a redshifted outflow component seen in the @xmath35co ( 10 ; 21 ; 32 ) lines @xcite and the cs ( 21 ; 32 ) lines @xcite driven from the neighboring protostellar source , l1551 irs 5 .
thus , the high - velocity redshifted component seen in the @xmath2co line most likely traces the redshifted outflow component from l1551 irs 5 .
figure [ astec18o ] shows the aste velocity channel maps of the c@xmath0o line at the same velocity bin as that of the @xmath2co velocity channel maps .
the c@xmath0o emission component associated with the protostar is seen from @xmath32=5.53 km s@xmath8 .
the c@xmath0o emission component exhibits an elongated feature at 6.44  6.90 km s@xmath8 along the northwest to the southeast direction , which is approximately consistent with the major axis of the central disk .
while the c@xmath0o emission peak at 6.44 km s@xmath8 is located to the northeast of the protostar , the emission location is systematically shifted from northeast to southwest from the blueshifted to redshifted velocities ( 6.44 km s@xmath8 to 7.81 km s@xmath8 ) .
these results indicate that there is a velocity gradient in the c@xmath0o emission along the direction of the minor axis of the central disk .
these c@xmath0o emission components most likely trace the protostellar envelope surrounding l1551 ne ( hereafter `` env '' ) . in the redshifted velocities from 7.35 km s@xmath8 ,
another emission component to the southwest of the protostar becomes evident .
this component appears to trace the same southwestern gas component seen in the @xmath2co emission . in even higher redshifted velocities ( 8.72  10.54 km s@xmath8 ) ,
the c@xmath0o counterpart of the western redshifted component found in the @xmath2co emission is also seen , although the highest velocity end of the c@xmath0o counterpart is lower than that of the @xmath2co emission .
figure [ astecs ] shows the same velocity channel maps as those of figures [ aste13co ] and [ astec18o ] but for the cs ( 76 ) line .
while the cs emission traces a gas component associated with the protostar in the velocity range of 5.54  7.81 km s@xmath8 , the most prominent cs emission component is the redshifted ( 8.27  9.18 km s@xmath8 ) component to the southwest of the protostar .
this southwestern redshifted component is also seen in the @xmath2co and c@xmath0o emission .
previous bima and nma observations of the cs ( 21 ; 32 ) lines have also found the same cs component , where the cs abundance is enhanced by a factor of a few @xcite .
this cs component is considered to be shock - excited molecular gas of the outflow driven from l1551 irs 5 @xcite .
the detection of the submillimeter cs ( 76 ) emission with the present aste observations , which traces warm ( @xmath10 60 k ) molecular gas @xcite , further supports this interpretation .
hereafter we call this redshifted component to the southwest of the protostar `` red1 '' , and the highest - velocity redshifted component seen in the @xmath2co and c@xmath0o lines to the west of the protostar `` red2 ''
.    in figure [ colli ] , we compare spatial distributions of env as seen in the c@xmath0o emission , red1 in the cs emission , and red2 in the @xmath2co emission .
line profiles in the peak positions of red1 and red2 , as well as those 20@xmath11 east of the protostar , are also shown ( note that the line profiles of env are shown in figure [ astesp ] . ) .
the two outflow components of red1 and red2 are located to the western side of env . toward red1 , the @xmath2co line profile shows an intense redshifted peak plus a red wing up to @xmath32 @xmath712 km s@xmath8 .
the c@xmath0o line profile exhibits redshifted emission `` plateau '' up to @xmath32 @xmath710 km s@xmath8 , and the cs line profile shows the emission only in the redshifted velocity range .
similar @xmath2co redshifted wing and c@xmath0o emission plateau are also seen toward red2 , where a weak redshifted cs emission is also present . on the other hand , toward the 20@xmath11 east of the protostar , such a c@xmath0o emission plateau is not seen and the @xmath2co redshifted wing is much less prominent , and the submillimeter cs emission is not detected .
these results imply that the outflow streams of red1 and red2 are terminated at the western side of env .
one possible interpretation for this termination is the collision between the outflow components and env .    in figure
[ collipv ] , we show the position - velocity ( p - v ) diagrams of the @xmath2co , c@xmath0o , and the cs lines along the cuts passing through the protostellar position and the peak positions of red1 ( p.a .
= 45@xmath19 ) and red2 ( p.a .
= 90@xmath19 ) . in the c@xmath0o p
- vs the envelope component is seen as labeled . to the west of env
the @xmath2co , c@xmath0o , and the cs emission components arising from red1 , and the @xmath2co and c@xmath0o emission components from red2 , can be identified .
there appears velocity gradients both in red1 and red2 as delineated by dashed lines , and the velocities of red1 and red2 become closer to the systemic velocity as the positions become closer to env .
these results suggest that the flow velocities of red1 and red2 are decelerated as they become close to env , presumably due to the collision with env .
figure [ astepvs ] shows aste p - v diagrams of the c@xmath0o and @xmath2co lines along the major and minor axes of the central cbd passing through the protostellar position . in the c@xmath0o p - vs ,
the line broadening at the protostellar position arising from the cbd is clearly seen , and the velocity gradient originated from the keplerian rotation ( blue curves in figure [ astepvs ] ) can be identified along the major axis . as discussed in the c@xmath0o line profile toward the protostellar position ( figure [ astesp ] ) , the centroid velocity of the cbd component and that of the envelope component are offset by @xmath70.2 km s@xmath8 ( two vertical dashed lines in figure [ astepvs ] ) , and the c@xmath0o p - vs exhibit this offset clearly .
furthermore , along the major axis the envelope component as seen in the c@xmath0o emission does not show any noticeable velocity gradient , and the keplerian rotation of the central disk is not continuous to the outer envelope . on the other hand , along the minor axis
the c@xmath0o envelope shows a clear velocity gradient along the northeast ( blue ) to southwest ( red ) direction ( tilted black dashed line in figure [ astepvs ] ) .
since the associated nir jets are blueshifted to the southwest and redshifted to the northeast @xcite , the southwestern side of the disk and envelope must be far - side and the northeastern side near - side , assuming that the midplane of the disk and envelope is perpendicular to the axis of the jets .
correspondingly , the southwestern redshifted component would then be located on the far side and the northeastern blueshifted component on the near side , suggesting the outward radial motion on the plane . in the @xmath2co p - v diagram along the minor axis ,
both red1 and red2 and their decelerated velocity structures ( red and pink dashed lines ) are identified .
these results indicate that the envelope is being dissipated due to the interaction with the redshifted outflowing gas driven from l1551 irs 5 .
we will further discuss this dispersing motion of the envelope with the combined sma+aste images of the c@xmath0o emission below .
figure [ mom0s ] compares moment 0 maps of the aste , sma , and the combined sma+aste image cubes of the c@xmath0o ( 32 ) line in l1551 ne .
while the aste image exhibits an extended envelope structure as described in the last subsection , the sma image shows the central smaller - scale ( @xmath71000 au @xmath34 800 au ) structure with its position angle ( @xmath7167@xmath19 ) consistent with that of the keplerian cbd @xcite .
combining these two images shows that the central structure found by the sma observations is embedded in the more extended envelope .
figure [ combch ] shows velocity channel maps of the sma+aste image cube of the c@xmath0o ( 32 ) line . in the highly blueshifted ( @xmath32 = 4.2  6.1 km s@xmath8 ) and redshifted ( 7.5  9.2 km s@xmath8 ) velocities , compact c@xmath0o emission located to the north and south of the protostellar binary are seen , respectively .
these compact components trace the keplerian cbd as we already discussed in our previous papers @xcite .
the extended envelope component appears in the lower velocity range ( 6.1  8.2 km s@xmath8 ) . in the blueshifted velocity range of 6.3  6.6 km s@xmath8
, the extended c@xmath0o emission is located predominantly to the eastern side of the protostellar binary , whereas in the redshifted velocity range of 7.5  8.2 km s@xmath8 to the western side .
thus , there is a velocity gradient along the east ( blue ) to west ( red ) direction in the envelope , which is also identified in the aste only c@xmath0o image cube as discussed in the last subsection . in the lower velocity region
( 6.6  7.4 km s@xmath8 ) the envelope emission extends over the entire region .    to highlight these velocity structures in the sma+aste image cube , we integrated the velocity channel maps over the high - velocity , middle - velocity , and the low - velocity blueshifted and redshifted regions , and the resultant images are shown in figure [ comb3 ] . in the high - velocity region ,
compact blueshifted and redshifted c@xmath0o emission located to the north and south of the protostellar binary , which originates from the keplerian cbd , are seen . in the middle - velocity region , as well as
the central disk component , an extended envelope component , with the blueshifted emission located predominantly to the east and the redshifted emission to the west , are seen .
as we discussed in the last subsection , one interpretation of this velocity structure is the dispersing gas motion in the flattened envelope . in the low - velocity region
both the blueshifted and redshifted emission spread over almost the entire field of view .
we note , however , that the blueshifted and redshifted emission peaks in the central compact emission are located slightly to the west and east of the protostellar binary , respectively .
the sense of this velocity structure is opposite to that of the extended envelope seen in the middle velocity range .
this velocity structure is also seen in the sma only image cube and the centroid velocity of this velocity feature matches with that of the cbd . as discussed in our previous paper @xcite , one interpretation for this velocity structure is an infalling gas motion in the inner envelope toward the central cbd .
such coexistence of inner infall and outer dispersing motions has also been identified in an envelope surrounding a t - tauri star dg tauri @xcite .
figure [ combpv ] ( left ) shows the p - v diagrams of the combined sma+aste image cube of the c@xmath0o line along the major ( upper panel ) and minor axes ( lower ) of the cbd . for comparison , keplerian rotation curves derived from the sma data @xcite are drawn in black lines , and the centroid velocity of the cbd and the envelope in vertical dashed lines . in the p - v along the major axis , the central compact component , cbd , exhibits a velocity gradient consistent with the keplerian rotation . on the other hand , the extended envelope component does not show any clear velocity gradient along the major axis , and the offset of the centroid velocity of the envelope from that of the cbd is clearly seen .
these results indicate that the keplerian rotation of the central cbd is discontinuous to the outer envelope . in the p - v diagram along the minor axis
, the central component exhibits a wide velocity width but no velocity gradient , as expected from the keplerian rotation . on the other hand , along the minor axis the emission ridges in the envelope component
are located predominantly to the northeast in the blueshifted velocities , and to the southwest to the redshifted velocities , as seen in the middle velocity ranges of figure [ comb3 ] .
this velocity structure can be interpreted as the dispersing gas motion in the flattened envelope .
we will interpret these detected velocity structures of the combined sma+aste c@xmath0o image cube with our toy model below .
the combined sma+aste c@xmath0o image cube reveals a sign of an outer dispersing envelope motion , possible inner infalling motion , and keplerian rotation in the central cbd . to interpret these results quantitatively , we constructed a toy model of the envelope and disk assuming that the disk and envelope are geometrically thin and co - planar .
our model is essentially the same as that by takakuwa et al .
( 2013 ) , which reproduces the velocity structures of the sma only image cube , but now the outer dispersing envelope component is added to reproduce the velocity structure of the extended envelope component traced by aste . as a model moment 0 map , we adopted a combination of two 2-dimensional gaussians , which are derived from the gaussian fitting to the observed sma+aste moment 0 map ( figure [ mom0s ] bottom - right ) .
the velocity structures in the inner component traced by the sma are described as ,    @xmath36    where @xmath37 is the radius , @xmath38 the gravitational constant , @xmath39 the mass of the protostellar binary , @xmath40 and @xmath41 are the rotational and radial velocities , respectively , @xmath42 the outermost radius of the keplerian rotation , and @xmath43 an adopted constant infall velocity beyond @xmath42 . as discussed by takakuwa et al .
( 2013 ) , we adopt @xmath39 = 0.8 @xmath14 , @xmath42 = 300 au , and @xmath43 = -0.6 km s@xmath8 , and the disk inclination angle @xmath44 .
the velocity structure of the outer envelope traced by aste is expressed as , @xmath45 where @xmath46 and @xmath47 are the innermost radius of the dispersing envelope and the dispersing velocity , respectively .
@xmath46 = 700 au is adopted , inferred from the comparison between the sma only and sma+aste moment 0 maps .
@xmath47 can be roughly estimated from the velocity gradient seen in the combined p - v diagram along the minor axis as shown in figure [ combpv ] .
the velocity shift of the emission ridge of the envelope across the minor axis is @xmath70.5 km s@xmath8 , centered on the centroid velocity of 6.7 km s@xmath8 , and at an inclination of @xmath48 @xmath7 62@xmath19 the dispersing velocity along the flattened envelope is @xmath47 @xmath7(0.5 km s@xmath8 / 2 ) / sin 62@xmath19 @xmath70.3 km s@xmath8 .
note that the centroid velocity of the envelope component is offset by @xmath7 - 0.2 km s@xmath8 from that of the cbd , and in our model this offset of the centroid velocity is incorporated .
the internal velocity dispersion @xmath49 = 0.4 km s@xmath8 , which reproduces the velocity channel maps of the central cbd @xcite , is adopted throughout the entire components .
figure [ modelch ] shows the velocity channel maps of our toy model . in the high - velocity blueshifted ( 4.6  5.7 km s@xmath8 ) and redshifted ranges
( 7.7  9.2 km s@xmath8 ) the features of the keplerian cbd are seen . in the blueshifted range of 5.9  6.3 km
s@xmath8 the envelope component is located to the east of the protostellar binary while in the redshifted range of 7.2  7.5 km s@xmath8 to the west , reflecting the dispersing motion in the envelope . in the lower - velocity range ( 6.4  7.0 km s@xmath8 )
the envelope component is extended both in the eastern and western sides , but the east ( red ) - west ( blue ) velocity gradient arising from the infalling motion in the central component is evident .
these characteristics reproduces primary features of the observed sma+aste c@xmath0o velocity channel maps in l1551 ne ( figure [ combch ] , [ comb3 ] ) .
figure [ combpv ] ( right ) shows the p - v diagrams of our toy model along the major and minor axes .
the model p - vs can reproduce primary features of the observed p - vs , such as the high - velocity component arising from the keplerian disk , and the extended envelope component which shows a velocity gradient along the minor axis reflecting the dispersing motion , but no rotation along the major axis .
these results show that the observed sma+aste c@xmath0o image cube in l1551 ne can be understood with three distinct velocity components ; a central keplerian cbd at a radius of @xmath9300 au , a possible infalling region at 300 au @xmath50 @xmath37 @xmath50 500 au , and an outer dispersing envelope with a dispersing velocity of @xmath70.3 km s@xmath8 at @xmath37 @xmath10500 au .
our aste mapping observations of the l1551 ne region have revealed three distinct gas components ; @xmath51 , env , red1 , and red2 .
env is the component of the protostellar envelope surrounding the protostellar binary system l1551 ne , and exhibits the dispersing gas motion as discussed above .
red1 and 2 are most likely the redshifted outflow components driven from l1551 irs 5 , which are interacting with env .    to discuss the interaction between env and red1 and 2 quantitatively ,
physical parameters of these components were estimated as follows .
first , the moment 0 map of each component was made in the relevant tracers and the velocity ranges listed in table [ comphys ] . as shown above ,
env is most clearly identified in the c@xmath0o ( 32 ) emission , and red2 in the @xmath2co ( 32 ) emission .
we adopt the cs ( 76 ) emission to deduce the physical condition of red1 , because in the cs emission red1 can be identified as a distinct gas component most easily , and the submillimeter cs emission must trace higher - temperature gas associated with outflows .
we also note that it is not straightforward to separate red2 in the 3-dimensional space from the other lower - velocity components .
we here identified red2 as the @xmath2co emission component in the highest redshifted velocity range of 9.47  12.53 km s@xmath8 , since in this velocity range red2 is a distinct , highest redshifted component seen primarily in the @xmath2co emission . after making the moment 0 maps of these distinct components , two - dimensional gaussian fittings to the moment
0 maps were performed to derive the central positions , beam - deconvolved sizes along the major and minor axes , and the total flux . for env a single two - dimensional gaussian fitting to the moment 0
map did not provide a satisfactory fit .
thus , for env two - component , two - dimensional gaussian fitting was performed , and the position was defined as the peak position of the central gaussian component , while the size was defined as the size of the outer gaussian component .
the line widths of these components ( @xmath52 ) were derived from the fwhm values of the relevant spectra toward the central positions . projected flow velocities of red1 and red2 , @xmath53 , were defined as the mean velocity at the center of the components with respect to the velocity of env (= 6.7 km s@xmath8 ) .
masses of the components ( @xmath54 ) were derived from the total fluxes , on the assumption of the lte condition , optically - thin emission , and the molecular abundances of @xmath55 = 1.7 @xmath34 10@xmath56 @xcite , @xmath57 = 7.7 @xcite , and @xmath58 = 6.8 @xmath34 10@xmath59 @xcite . since the peak brightness temperature ( @xmath60 ) of the c@xmath0o emission toward the protobinary position exceeds @xmath613.4 k
, the excitation temperature ( @xmath62 ) must exceed @xmath109.0 k as @xmath63 must be higher than @xmath64 . on the other hand , to excite the submillimeter cs ( 76 ) emission a gas temperature as high as @xmath760 k is required @xcite .
thus , @xmath65 = 10 - 60 k is assumed as a probable range of @xmath65 .
virial masses of the components ( @xmath66 ) were derived as ; @xmath67 where @xmath68 in the above expressions @xmath69 denotes the geometrical mean of the sizes along the major and minor axes of the components , and @xmath70 the effective sound speed of gas .
the isotropic velocity dispersion inside the components is assumed . for red1 and 2 , the sizes ( @xmath74000 au ) are as small as the beam size ( @xmath73200 au ) , and thus the line widths derived from the central spectra represent the average velocity dispersions . for env ,
there is no discernible spatial variation of the c@xmath0o line width as shown in the p - v diagrams ( figures [ collipv ] , [ astepvs ] ) , except for that arising from the central disk .
thus , the assumption of the isotropic velocity dispersion inside these components is probably valid .
momenta ( @xmath71 ) , internal gas energies ( @xmath72 ) , and kinetic energies ( @xmath73 ) of the components were derived as ;    @xmath74    @xmath75    @xmath76    table [ comphys ] summarizes these derived physical parameters .
the size ( @xmath720000 au ) and lte mass ( @xmath71 - 3 @xmath14 ) of env are typical of protostellar envelopes @xcite .
the virial mass of env is comparable to or slightly larger than the lte mass plus the central protobinary mass . from the comparison between the virial and lte masses it is not straightforward to ascertain whether env is gravitational bound or not . on the other hand
, the identification of the dispersing gas motion in env implies that env is gravitationally unbound .
the outflow components of red1 and 2 are characterized with their compact sizes ( @xmath94000 au ) and wide line widths ( @xmath712 km s@xmath8 ) .
the lte masses of the outflow components are likely lower than the corresponding virial masses , although for red1 only the upper limit of the virial mass is obtained .
these results indicate that red1 and 2 are not gravitationally - bound gas condensations , consistent with our interpretation that red1 and 2 are outflowing gas components driven from l1551 irs 5 .
the summed flow momentum of red1 and 2 is comparable to that of the momentum of the dispersing motion in env .
furthermore , the net internal+kinetic energy of these outflow components is also comparable to the internal+kinetic gas energy of env .
these results imply that red1 and 2 can exert significant dynamical impacts on env , and can induce the observed dispersing gas motion of env through the interactions .
our new aste observations and combining the aste data with our previous sma data in l1551 ne have revealed that the protostellar envelope around the keplerian cbd is being dispersed .
the dispersing envelope does not show any rotating gas motion in contrast with the keplerian rotation in the cbd , and the centroid velocity of the envelope is @xmath70.2 km s@xmath8 blueshifted with respect to that of cbd .
these results indicate that the protostellar envelope is kinematically distinct from the central cbd .
the outflow components from l1551 irs 5 , as seen most prominently in the submillimeter cs ( 76 ) line ( red1 ) and the high - velocity redshifted @xmath2co ( 32 ) line ( red2 ) , appear to collide and interact with the protostellar envelope .
the net momentum , and kinetic and internal energy of these outflow components are comparable to those of env , suggesting that the outflow collisions can enforce significant dynamical impacts on the envelope . a natural interpretation of these results is that the protostellar envelope in l1551 ne is being disrupted through the interaction with the outflows driven from l1551 irs 5 .
figure [ scheme ] ( bottom ) shows a schematic , inferred configuration of the l1551 region viewed from the north , in parallel with the scuba 850 @xmath77 dust - continuum image @xcite as a representative projected image onto the sky plane ( top ) . in the scuba image a torus - like continuum feature
connecting the envelopes around l1551 irs 5 and ne is seen as delineated with a dashed green curve .
this feature is drawn schematically in the bottom as the northernmost feature .
since the envelope around l1551 ne is being collided with the redshifted outflow components driven from l1551 irs 5 ( red1 and 2 ) , the location of l1551 ne must be backward with respect to the location of l1551 irs 5 along the line of sight .
as already described the western side of the flattened envelope around l1551 ne is far - side and the eastern side near - side .
thus the blueshifted emission to the east and the redshifted emission to the west are interpreted as the dispersing gas motion in the flattened envelope .
our aste single - dish observations of the @xmath2co ( 32 ) , c@xmath0o ( 32 ) , cs ( 76 ) , and so ( 7@xmath46@xmath5 ) lines , as well as sma and alma observations of the @xmath2co ( 32 ) and c@xmath0o ( 32 ) lines @xcite , do not find any evidence for molecular outflows driven from l1551 ne itself .
on the other hand , previous nir observations of l1551 ne found collimated [ fe ii ] jets driven from source a as well as a reflection nebula with its apex close to the position of source b @xcite .
molecular outflows driven from l1551 ne have also been reported in the @xmath35co ( 32 ; 21 ; 10 ) lines @xcite , although it is not straightforward to disentangle the outflows from l1551 irs 5 and ne because of the close overlaps of the two outflows .
our aste observations show that interactions with the powerful l1551 irs 5 outflows have greater impact on the disruption of the envelope than the outflow driven from l1551 ne itself .
this is in contrast with previous findings of interactions between protostellar envelopes and outflows or stellar winds driven from the sources embedded in the envelopes @xcite .
yokogawa et al . ( 2003 ) have suggested that the interactions with the outflows driven from l1551 irs 5 triggered the star formation in l1551 ne .
we consider , however , that the formation of the protostellar binary system l1551 ne has initiated before the arrival of the irs 5 outflows , because at present the protostellar envelope is kinematically distinct and thus physically detached from the cbd and the protostellar binary system .
the similar binary masses ( @xmath70.8 @xmath14 in l1551 ne and @xmath70.5 @xmath14 in irs 5 ) , envelope masses ( @xmath71 @xmath14 ) , and the bolometric temperatures ( @xmath15 @xmath790 k ) indicate that both l1551 irs 5 and ne were born at a similar age .
after their birth both l1551 ne and irs 5 grow simultaneously . the angular separation between l1551 irs 5 and ne is @xmath72@xmath185 , and assuming the transverse velocity of the outflows driven from l1551 irs 5 is @xmath710 km s@xmath8 @xcite , the transverse time for the irs 5 outflow to reach the l1551 ne location is @xmath710@xmath78 yr .
this time scale is approximately consistent with the estimated age of l1551 ne of 0.65.0@xmath3410@xmath78 yr @xcite .
after the outflows from l1551 irs 5 have arrived at the protostellar envelope around l1551 ne , the outflows start disrupting the envelope through the interactions . the possible infalling motion in the inner envelope found with the sma may be a remnant gas motion before the initiation of the interactions with the irs 5 outflows .
such a remnant infalling gas motion embedded in the dispersing envelope has also been identified in a t - tauri star dg tauri @xcite .
the envelope dispersion via the interactions with the irs 5 outflows will impact the final fates of the protostellar binary system of l1551 ne significantly , as will be discussed in the next subsection .      binary stars with the primary masses comparable to that of the sun show a wide range of the binary mass ratios , and the distribution of the mass ratio ( @xmath12 ) is essentially flat from @xmath130.1 to 1 @xcite .
physical mechanisms to reproduce such a wide range of binary mass ratios have been controversial , and gas accretion in the cbds onto the protostellar binaries has been discussed as a mechanism to set the final binary mass ratios . smoothed particle hydrodynamic ( sph ) simulations of circumbinary accretion show that the majority of materials accrete onto the secondary , because the secondary orbits further from the center of mass of the binary and thus sweeps more materials in the cbd than the primary @xcite .
our 3-d adaptive mesh refinement ( amr ) hydrodynamic simulation @xcite to reproduce the alma results of the cbd in l1551 ne also shows that the secondary ( source b ) accretes more than the primary ( source a ) @xcite . on the other hand , high - resolution grid - based simulations by ochi et al .
( 2005 ) and hanawa et al .
( 2010 ) show that primaries accrete more than secondaries .
latest sph simulations demonstrate that flows from the secondary to the primary within the roche lobes are a sensitive function of gas temperatures , and in the case of the high temperature the primary accretes more than the secondary , and vice versa @xcite .
our aste observations of an archetypal protostellar binary l1551 ne imply that dispersion of protostellar envelopes , which replenish cbds with fresh materials , also needs to be taken into account as a physical mechanism to set the binary masses and the ratios .
the present masses of sources a and b in l1551 ne are estimated to be @xmath70.67 @xmath14 and @xmath70.13 @xmath14 , respectively , and the mass of the cbd @xmath70.026
@xmath14 @xcite .
the mass of the possible infalling component around the cbd is even smaller ( @xmath70.0023 @xmath14 ) @xcite .
thus , even if all the amount of the material in the cbd plus the infalling component is accreted onto source b , the mass ratio does not change much ( from 0.19 to 0.24 ) .
replenishment of the material to the cbd from the surrounding envelope is required to significantly change the mass and the mass ratio of the protostellar binary system .
our aste observations show that the envelope is being dispersed , however , and that further replenishment of materials from the envelope to the cbd will not be expected
. therefore , even though l1551 ne is a class i protostellar binary associated with the cbd and the protostellar envelope , the mass and the mass ratio have been already close to the final values . physical processes before the start of the envelope dispersion through the interaction with the outflow driven from l1551 irs 5 , such as the fragmentation of the initial pseudo - disk and subsequent accretion from the cbd and the envelope , must be responsible to set the mass and mass ratio of the binary system .
our aste observations of l1551 ne indicate that dispersion of natal envelopes surrounding the protostellar binary systems is one of the key physical phenomena to stop the growth of the binary and to determine the final masses and the mass ratios of the binary .
as described in the last subsection , interactions between protostellar envelopes and outflows , which causes the envelope dissipation , are ubiquitous .
systematic studies of envelope dispersion surrounding the protostellar binaries with various masses and mass ratios will be intriguing to understand the physical mechanism to set the final binary mass ratios and the origin of the wide range of the binary mass ratios .
we have conducted single - dish mapping observations of dense gas around an archetypal protostellar binary l1551 ne in the c@xmath0o ( @xmath1=32 ) , @xmath2co ( @xmath1=32 ) , cs ( @xmath1=76 ) , and the so ( @xmath3=7@xmath46@xmath5 ) lines with aste .
we have also combined the aste mapping data in the c@xmath0o line with the interferometric data taken with the sma .
the main results are summarized below :    * all the four molecular lines are detected toward the protostellar position , with multiple and/or non - gaussian spectral shapes .
in particular , the c@xmath0o spectrum consists of a narrow ( @xmath790.7 km s@xmath8 ) component with a central velocity of @xmath806.7 km s@xmath8 and a wide ( @xmath792.2 km s@xmath8 ) component with a central velocity of @xmath806.9 km s@xmath8 .
the wide component originates from the compact circumbinary disk ( cbd ) , while the narrow component from the extended protostellar envelope . *
the c@xmath0o map in l1551 ne primarily traces a @xmath720000-au scale protostellar envelope ( env ) , elongated along the northwest to southeast direction , approximately parallel to the major axis of the central cbd .
velocity channel maps of the c@xmath0o line show that the blueshifted ( @xmath806.26.7 km s@xmath8 ) emission is located to the northeast while the redshifted emission ( @xmath806.7 8.0 km s@xmath8 ) to the southwest ; @xmath51 , a velocity gradient of env along the minor axis , while there is no detectable velocity gradient along the major axis .
the observed velocity gradient in env can be interpreted as a dispersing gas motion .
velocity channel maps of the high - velocity redshifted ( @xmath818.0 km s@xmath8 ) @xmath2co emission show two @xmath94000-au scale redshifted components to the @xmath740@xmath11 southwest of the central protostar ( red1 ) and @xmath720@xmath11 west ( red2 ) .
the submillimeter cs ( 76 ) emission selectively traces red1 , suggesting high - temperature in this high - velocity component .
both red1 and red2 exhibit similar velocity gradients such that the velocities become closer to the systemic velocity of env as the positions become closer to that of env .
these high - velocity redshifted components most likely trace the outflow components driven from the neighboring protostellar source , l1551 irs 5 , which are colliding with env surrounding l1551 ne . *
the combined aste+sma map in the c@xmath0o emission shows that the @xmath6300 au scale cbd and @xmath6700 au scale infalling component found in our previous sma observations are surrounded by env with a slightly blueshifted velocity ( @xmath70.2 km s@xmath8 ) .
env does not show any rotational motion but dispersing gas motion with an outward velocity of @xmath70.3 km s@xmath8 in contrast to the inner gas components .
these results indicate that env is kinematically distinct from the inner gas components .
* the net momentum , kinetic and internal energies of the redshifted outflow components driven from l1551 irs 5 are comparable to those of env , suggesting that the outflow collisions can enforce significant dynamical impacts on env .
thus , the dispersing motion in env is likely caused by the interactions with the outflows driven from l1551 irs 5 .
* since env is being dispersed and the inner gas components do not have sufficient materials to alter the mass and mass ratio of the protostellar binary of l1551 ne , the current mass ( @xmath70.8 @xmath14 ) and the mass ratio ( @xmath70.19 ) are close to the final values , even though l1551 ne is a young class i protostellar binary .
our aste+sma observations of l1551 ne suggest that dispersion of natal protostellar envelopes which can replenish the cbds with new materials is one of the important physical mechanisms to set the final binary mass ratios .
the interactions with the outflows from l1551 irs 5 are unlikely to trigger the binary formation in l1551 ne as suggested by previous studies , but suppress the further growth of the protostellar binary system .
we would like to thank t. matsumoto , t. hanawa , j. lim , n. ohashi , and p. t. p. ho for their fruitful discussions .
we are grateful to s. ohashi for his support during our aste observations .
s.t . acknowledges a grant from the ministry of science and technology ( most ) of taiwan ( most 102 - 2119-m-001 - 012-my3 ) in support of this work .
andr , p. , ward - thompson , d. , & barsony , m. 2000 , in protostars and planets iv , ed . mannings , v. , boss , a. p. , & russell , s. s. ( tucson : univ .
arizona press ) , 59 arce , h. g. & sargent , a. i. 2006 , , 646 , 1070 bate , m. r. 1997 , , 285 , 16 bate , m. r. , & bonnell , i. a. 1997 , , 285 , 33 bate , m. r. 2000 , , 314 , 33 chen , x. , arce , h.  g. , zhang , q. @xmath82 2013 , , 768 , 110 chou , t .- l . , takakuwa , s. , yen , h .- w . , ohashi , n. , & ho , p. t. p. 2014
, , 796 , 70 crapsi , a. , caselli , p. , walmsley , c.  m. , @xmath82  2004 , , 420 , 957 dutrey , a. , di folco , e. , guilloteau , s. , @xmath82 2014 , , 514 , 600 froebrich , d. 2005 , , 156 , 169 fuente , a. , martn - pintado , j. , bachiller , r. , rodrguez - franco , a. , & palla , f. 2002 , , 387 , 977 fuller , g.  a. , & ladd , e.  f.  2002 , , 573 , 699 goodwin , s. p. 2013
, , 430 , l6 hanawa , t. , ochi , y. , & ando , k. 2010 , , 708 , 485 hayashi , m. & pyo , t .- s .
2009 , , 694 , 582 jrgensen , j.  k. , schier , f. l. , & van dishoeck , e. f. 2004 , , 416 , 603 kamazaki , t. , ezawa , h. , tatematsu , k. , @xmath82 2005 , astronomical data analysis software and systems xiv , asp conf .
ser . , 347 , ed .
p. shopbell , m. britton , & r. ebert , ( san francisco : asp ) , 533 kitamura , y. , kawabe , r. , & saito , m. 1996 , , 457 , 277 kohno , k. 2005 , the cool universe : observing cosmic dawn , asp conf .
344 , ed . c. lidman & d. alloin , ( san francisco : asp )
, 242 matsumoto , t. 2007 , , 59 , 905 momose , m. , ohashi , n. , kawabe , r. , hayashi , m. , & nakano , t. 1996 , , 470 , 1001 momose , m. , ohashi , n. , kawabe , r. , nakano , t. , & hayashi , m. 1998 , , 504 , 314 moriarty - schieven , g. h. , butner , h. m. & wannier , p. g. 1995 , , 445 , l55 moriarty - schieven , g. h. , johnstone , d. , bally , j. , & jenness , t. 2006 , , 645 , 357 motte , f. , & andr .
, , 365 , 440 myers , p. c. , evans , n. j. , ii & ohashi , n. 2000 , in protostars and planets iv , ed . mannings , v. , boss , a. p. , & russell , s. s. ( tucson : univ .
arizona press ) , 217 ochi , y. , sugimoto , k. , & hanawa , t. 2005 , , 623 , 922 plambeck , r. l. , & snell , r. l. 1995 , , 446 , 234 raghavan , d. , mcalister , h. a. , henry , t. j. @xmath82 2010 , , 190 , 1 reipurth , b. , yu , k. c. , heathcote , s. , bally , j. , & rodrguez , l. f. 2000 , , 120 , 1449 reipurth , b. , rodrguez , l. f. , anglada , g. , & bally , j. 2002 , , 124 , 1045 reipurth , b. , clarke , c. j. , boss , a. p. @xmath82 2014 , protostars and planets vi , 267 saito , m. , kawabe , r. , kitamura , y. , & sunada , k. 2001 , , 547 , 840 schilke , p. , groesbeck , t. d. , blake , g. a. , & phillips , t. g. 1997 , , 108 , 301 stojimirovi , i. , narayanan , g. , snell , r. l. , & bally , j. 2006 , , 649 , 280 takahashi , s. , saito , m. , takakuwa , s. , & kawabe , r. 2006 , , 651 , 933 takakuwa , s. , ohashi , n. , & hirano , n. 2003 , , 590 , 932 takakuwa , s. , ohashi , n. , ho , p. t. p. , @xmath82 2004
, , 616 , l15 takakuwa , s. , kamazaki , t. , saito , m. , yamaguchi , n. , & kohno , k. 2007 , , 59 , 1 takakuwa , s. , ohashi , n. , bourke , t.  l. , @xmath83 @xmath84  2007 , , 662 , 431 takakuwa , s. , & kamazaki , t. 2011 , , 63 , 921 takakuwa , s. , saito , m. , lim , j. , @xmath82 2012 , , 754 , 52 takakuwa , s. , saito , m. , lim , j. , & saigo , k. 2013 , , 776 , 51 takakuwa , s. , saito , m. , saigo , k. , @xmath82 2014 , , 796 , 1 tang , y .- w . , dutrey , a. , guilloteau , s. , @xmath82 2014 , , 793 , 10 tobin , j. j. , chandler , c. j. , wilner , d. j. , @xmath82 2013 , , 779 , 93 wilson , t.  l. , & rood , r.  1994 , , 32 , 191 wu , p .- f . , takakuwa , s. , & lim , j. 2009 , , 698 , 184 yokogawa , s. , kitamura , y. , momose , m. , & kawabe , r. 2003 , , 595 , 266 young , m. d. , baird , j. t. , & clarke , c. j. 2015 , , 447 , 2907      llllccccccccc env & c@xmath0o ( 32 ) & 04 31 44.4 & 18 08 36.4 & 5.31 - 8.15 & 22100@xmath3418200 & 0.72 & & 1.3 - 3.1 & 5.3 & 0.38 - 0.93 & 3.5 - 8.7 & 1.14 - 2.78 + red1 & cs ( 76 ) & 04 31 42.0 & 18 08 06.0 & 8.04 - 9.13 & @xmath504000 & 0.97 & 1.9 & 0.02 - 0.83 & @xmath501.9 & 0.039 - 1.574 & 0.10 - 4.23 & 0.73 - 29.74 + red2 & @xmath2co ( 32 ) & 04 31 42.8 & 18 08 35.0 & 9.47 - 12.53 & 4300@xmath344000 & 2.0 & 3.5 & 0.02 - 0.05 & 8.4 & 0.069 - 0.169 & 0.42 - 1.04 & 2.4 - 5.9 +