You are an expert at summarizing long articles. Proceed to summarize the following text: 
both indirect evidence @xcite and later direct imaging @xcite , have shown that disks are quite common in young stellar systems . these disks are commonly thought @xcite to be sites for planet or brown dwarf formation . a large fraction of stars are formed in binary systems @xcite and in the same star formation regions as single stars . theory suggests that the most likely mechanisms responsible for forming jovian mass planets or low mass brown dwarfs are either gravitational collapse of large scale spiral structure or coagulation of small solid grains followed by later accretion of additional gas ( ` core accretion ' ) in the disks of forming stellar systems . evaluating the effectiveness of these mechanisms is important for understanding the origin of our own solar system as well as planetary systems in other mature single or multiple systems . the l1551  irs  5  system serves as a useful observational testbed for comparison to theoretical modeling because of its relative youth ( @xmath1  yr , * ? ? ? * ) and many previous detailed observations ( see e.g. * ? ? ? this system consists of an extended nebulosity some 2400@xmath21100  au in size with an inner core of 220@xmath276  au @xcite . two bipolar jets flow outward in each direction from the core in the plane perpendicular to its long axis . the core has been resolved into two sources with projected separation of about 50  au and inferred disk masses of @xmath3@xmath4 , each @xmath520 - 25  au in diameter @xcite . the total mass in the core has been estimated to contain 0.51.0@xmath4  of material @xcite , which produces @xmath530@xmath6  in luminous output @xcite . we present a numerical simulation of a binary star / disk@xmath7star / disk system using a two dimensional ( @xmath8 ) smoothed particle hydrodynamic ( sph ) code . the dimensions of the disks and semi - major axis of the binary are chosen to be similar to the inner core region of l1551  irs  5 . in the absence of strong constraints on the constituents of the binary ( e.g. the masses of the two stars ) , we choose to set up a binary system consisting of identical components , obtained by setting up a single system in isolation , then duplicating it exactly . we assume each star and disk have mass @xmath9@xmath4  and @xmath10@xmath4 , respectively . the disk radius is set to @xmath11  au which , for a semi - major axis of @xmath12  au , is comparable to the largest stable streamline @xcite . the mass and temperature of the disk are distributed according to @xmath13 and @xmath14 power laws respectively . the absolute scale of each power law is determined from the disk mass , the radial dimensions of the disk and the condition that the toomre stability parameter , @xmath15 , is no smaller than @xmath16=1.5 over the entire disk . this value ensures that the simulation begins in a state marginally stable against the growth of spiral structure , so that we do not accidently ` discover ' a collapsed object early in the evolution which in reality is an artifact of our initial condition . both density and temperature are free to vary in time and space , so the initial condition will not prevent spiral structure growth or fragmentation , if the evolution leads to such . the gas is set up on circular orbits around the star so that pressure and gravitational forces exactly balance centrifugal forces . radial motion is zero . the magnitudes of the pressure and self - gravitational forces are small compared to the stellar gravity , so the disk is nearly keplerian in character . approximately 60000 equal mass particles are set on a series of concentric rings around the star in a single , star / disk system , then duplicated , bringing the total number of particles to @xmath5120000 . the two stars and disks are offset equal distances in the @xmath17 and @xmath18 directions . we define the binary semi - major axis to be @xmath12  au , similar to l1551  irs  5 . only weak constraints on eccentricity exist in l1551  irs  5 , primarily consisting of the sizes of the observed disks : eccentricities larger than @xmath19 would lead to rapid roche lobe overflow . we set @xmath19 to be the initial value in this simulation . the system is at apoapse at time @xmath20 with the orbital velocities defined by approximating each star@xmath7disk system as a point mass , so that the orbit determination reduces to the solution of the two body problem . the disks are self gravitating and each star is modeled as a point mass free to move in response to gravitational forces from the rest of the system . the stellar gravitational forces are calculated using a plummer potential with a softening radius of 0.2  au , which also serves as an accretion radius , @xmath21 . sph particles with trajectories that pass closer than @xmath21 to a star are absorbed , and the star s mass and momentum increase accordingly .    the thermodynamic evolution is identical to that described in @xcite . thermal energy is added to the gas due to active hydrodynamic processes using an artificial viscosity scheme , which approximately models shocks and turbulence . this heating is roughly equivalent in magnitude to an alpha model with @xmath22 . thermal energy is removed from the disk gas by radiative cooling due to passive blackbody emission from the disk s photosphere surfaces . the blackbody temperature is calculated at each time step and for each sph particle . this treatment remedies a major shortcoming of previous models @xcite which used a ` locally isothermal ' or ` locally adiabatic ' approximation to show that relatively low mass disks can undergo fragmentation and/or collapse , despite earlier claims @xcite that a very massive disk is required . the system is evolved for eight binary orbits , or 2700  yr ( @xmath23=350  yr ) . figure [ fig : disk - pair ] shows the system shortly before and after the fourth periapse passage of the two components ( measured from the beginning of the simulation ) . before each periapse , the two disks are smooth and exhibit no visible spiral structure , although they are no longer perfectly ` round ' . during and after periapse , each disk develops strong , two armed spiral structures due to the mutual tidal interactions of the binary . the structures decay to a smooth condition like that in the top panel over the next @xmath24@xmath23 . the cycle repeats with little variation as the system again approaches periapse , and we expect that further evolution will be similar .    the spiral structures decay because internal heating in the disk increases the stability of the disks against spiral arm growth , as measured by the toomre stability ( fig . [ fig : temp - prof ] top ) of each disk . the minimum value of @xmath15 increases from its initial @xmath16=1.5 to @xmath16@xmath25 before periapse and @xmath16@xmath26 afterwards . these values are the same before and after each successive periapse passage and both are well above the @xmath27 values for which spiral structures are expected to grow . therefore , if we suppose that jovian planets form via gravitational collapse or fragmentation of spiral structure in disks , their formation will be unlikely in this system . the high stability is due to an increase in the disk temperature . in fact , the temperatures are high enough to cause some grain species to be vaporized . this is important because the core accretion model for planet formation requires that solid grains can coagulate and are not instead repeatedly returned to vapor state . water ice may be particularly important because it composes 40% by mass of the solid material in the disk and is among the most volatile grain species @xcite , vaporizing at @xmath28  k.    the bottom panel of figure [ fig : temp - prof ] shows the disk midplane temperatures obtained in this simulation . only in a region with @xmath29  au does the temperature reach low enough values that water ice can form , even temporarily . however in this region , the matter is most subject to shocks generated by the spiral structure produced by the binary interaction , which raise its temperature by as much as a factor three over the azimuth average . the spiral patterns co - rotate with the orbit of the binary and the orbital period at 15  au is @xmath582  yr , so material everywhere in the disk will have time to travel through the spiral arms several times before they decay . a simulation run with zero eccentricity also produces spiral structure . in this run the azimuth averaged temperatures were similar to fig . [ fig : temp - prof ] , but the temperature in the shocked regions are not as extreme . water ice would still be vaporized everywhere but more refractory species may not be .    in the outer disk , grains less than @xmath30  mm in size which pass through such a shock and which contain water ice will be vaporized on a timescale of @xmath31  s depending on the temperature @xcite and the remaining more inert species may dissagregate . in this region , passage through the warmest part of a spiral arm requires @xmath32  yr , so sufficient time exists to return grains of this size to the gas phase . grain growth may still occur between the spiral arms and when the spiral arms have decayed , but must begin with gaseous material each time , so growth of solid material into larger entities will be suppressed . temperatures at high altitudes are lower , but grains which form there will tend to sink to the midplane as they grow larger and also be destroyed . therefore jovian planet formation by the core accretion mechanism will occur much more slowly , if at all in this system . the weak link remaining in the argument against the core accretion mechanism is the lack of knowledge of the microphysics important for dust coagulation . for example , one could imagine that growth of silicate and iron grains is catalyzed by temporarily enhanced cross sections as mantles of more volatile material form on their surfaces and the gravitational torques produced by the binary interaction enhance mixing throughout the disks . even if this type of interaction takes place , the eventual formation of planet sized objects remains in doubt . as the rocky aggregates grow , conditions appropriate for dust coagulation ( e.g. ` perfect sticking ' ) break down and collisions between particles become increasingly disruptive due to the finite strength of the aggregates . the disruption of solid bodies depends strongly on the relative velocity of the impactor and target particles , with disruption occurring for velocities @xmath331 - 3  km / s for planetesimal sized targets ( @xmath34  km ) @xcite . on average in an accretion disk , the relative velocity of planetesimals will be proportional to their eccentricities , @xmath35 . for our model , a relative velocity of 1 km / s corresponds to an eccentricity of @xmath36 at 1  au , or @xmath37 at 10  au . we have seen that gravitational torques are strong enough to generate large amplitude spiral structure as they drive the gas onto eccentric orbits . gas eccentricities are quickly damped due to shock dissipation , but planetesimals are only weakly coupled to the gas and will have time to encounter other objects and collide or to increase their eccentricities still further as the evolution proceeds . if particle eccentricities can grow to @xmath38 , we expect that the growth of kilometer sized bodies will be suppressed in binary systems such as the one modeled here . however , a more detailed analysis must be done in order to constrain this possibility . the conclusions regarding planet formation will remain valid as long as the temperatures determined from the model are lower limits on the temperatures present in real systems . if the model produces temperatures which are too low , then real systems will be even more stable against spiral structure growth and fragmentation and less likely to produce large , coagulated grains . if they are too high , the model could inaccurately portray the disk as too stable . are the temperatures produced in the model too low ? we can constrain the temperatures in the model by comparing the radiated energy from observed systems ( here specifically to l1551  irs  5 ) , to that produced by the simulation in various wavelength bands . this comparison requires that we relate the luminous output to the temperature .    in regions where the optical depth is high , as it is in the accretion disks , the radiated emission can be approximated as a blackbody with a temperature of the disk photosphere . the disk s midplane and photosphere temperatures are then related to each other by a given rosseland opacity and the local vertical density / temperature profile . we determine such a profile as by - product of the cooling model in this work , under the assumption that the vertical structure is instantaneously adiabatic . other work @xcite has shown that the structure may instead be super adiabatic . if this is the case then the midplane temperatures will be higher , and our conclusion is stronger . the opposite case may be true instead : large relative heating can occur at high altitudes , even though the high altitude heating is small in an absolute sense @xcite . this means that vertical temperature structure may become distorted and an incorrect midplane temperature could be inferred . high altitude dynamical heating will play a role similar to high altitude passive heating from stellar photons . @xcite show that this process produces a temperature inversion high above the photosphere , but this region contributes negligibly to the radiated flux . therefore , we can rely on the modeled radiated output of the simulation to represent accurately the temperatures at the disk midplane , given the physical processes included in the calculation . our conclusions about the formation of planets will be confirmed if we find that the energy output from the disks is equal to or less than that observed .    to obtain a valid comparison between the observed and modeled fluxes , we must be certain that the flux from other parts of the system ( e.g. the circumbinary disk and envelope ) is not a significant contributor to the observed flux used in the comparison . we must also be certain that extinction between the source and the observer has not altered the emitted flux . we therefore require very high spatial resolution photometry at long wavelengths , not affected by extinction .    in figure [ fig : freq - cmp ] , the highest available resolution , long wavelength observations of l1551  irs  5  are plotted and compared to the flux densities produced from the simulation . for all wavelengths between 1.3  cm and 870  @xmath39 m , the observed fluxes exceed those obtained from the simulation by a factor of @xmath40 . the differences at 1.4  mm and 870  @xmath39 m are larger , a factor @xmath41 , however they do not resolve the binary and may contain some contaminating flux from the circumbinary environment . given these comparisons , we can conclude that the temperatures in the disk are conservatively estimated by the simulation ( too low ) and that the disks in l1551  irs  5are more stable against spiral arm growth and fragmentation and less likely to allow dust coagulation than in the simulation . our conclusion that planet formation is unlikely in the l1551  irs  5  system is secure . the inventory of physical processes considered in this calculation is not complete . the total luminosity will include direct contributions at the short wavelength end of the spectrum from not only the two circumstellar disks , but also the two stars . accretion of material from the inner disk edge ( 0.2  au for this simulation ) onto the stellar surface will also add to the total . at longer wavelengths , the circumbinary disk , the infalling envelope and radio emission from bipolar polar outflows will add to the total . of these processes , only outflows remove thermal energy from the disks , but such outflows are thought to originate very close to the star . they will not remove energy from the part of the disk important for planet formation . the rest will either add directly to the short wavelength spectrum or heat the cooler parts of the disk by absorption of short wavelength radiation , which subsequently reradiates at longer wavelengths ( see e.g. * ? ? ? * ) . including them in the calculation can only raise disk temperatures and strengthen our conclusions . estimates of the total luminosity ( not shown ) from all these sources are within a factor of two of the observed 30@xmath6  luminosity of l1551  irs  5 . we are encouraged by this agreement , and expect that the overall agreement between model and observations would become closer if radiative reprocessing were included . a number of questions remain . what is the distribution of planetesimal random velocities ? will lower mass disks be more susceptible to planet formation because of increased efficiency of radiative energy losses ? what happens for binary systems with different separations or with unequal stellar mass ratios ? how distant is distant enough , so that one component of the binary does not strongly influence the other ? these questions will be addressed more completely in a followup paper . i gratefully acknowledge david koerner for the generous release of the 1.3  cm fluxes before their publication . i thank brian pickett and willy kley for helpful conversations during the development of this paper and for brian s later comments as referee .
we show that planet formation via both gravitational collapse and core accretion is unlikely to occur in equal mass binary systems with moderate ( @xmath0  au ) semi - major axes . 
 internal thermal energy generation in the disks is sufficient to heat the gas everywhere so that spiral structures quickly decay rather than grow or fragment . this same heating will inhibit dust coagulation because the temperatures rise above the vaporization temperatures of many volatile materials . 
 we consider other processes not included in the model and conclude that our temperatures are conservatively estimated ( low ) , i.e. planet formation is less likely in real systems than in the model .