molecular clouds are observed through the emission of a number of molecular transitions that provide a wealth of information about their chemical composition , gas temperature and density , magnetic field strength , fractional ionization , structure and kinematics .
this information is essential to our understanding of the process of star formation .
the interpretation of molecular emission line maps is not always unique .
the main source of uncertainty is the absence of the third spatial dimension ( along the line of sight ) in the observational data .
statistical properties of the velocity and density distributions along the line of sight are difficult to disentangle . furthermore , the components of the gas velocity on the plane of the sky are unknown .    two dimensional images of molecular clouds are usually converted into three dimensional `` objects '' using the radial velocity instead of the third spatial dimension .
this method can be useful to separate individual mass condensations from each other , since it is conceivable that their relative velocity is larger than their internal one .
however , velocity blending or the lack of well defined condensations along the line of sight may cause significant uncertainties ( issa , maclaren & wolfendale 1990 ; adler & roberts 1992 ; ballesteros  paredes , vazquez  semadeni & scalo 1999 ; pichardo et al .
2000 ; ostriker , stone & gammie 2001 ; lazarian et al . 2001 ; ballesteros  paredes & mac low 2002 ) .    due to the difficulty of a direct interpretation of the observational data , a `` forward approach '' that starts from a rather general theoretical model and synthesizes
its observational properties can be more instructive .
different models may sometimes satisfy the same set of observational constraints , but they should also provide guidance for further observational studies that could help select the correct model .    ideally , numerical models to be compared with observed spectral line data cubes should be based on the numerical solutions of the mhd equations , in the regime of highly super  sonic turbulence , and on radiative transfer calculations . in some works ,
stochastic fields are used instead of the solution of the mhd equations and in most studies the radiative transfer calculation is omitted , in favor of density  weighted velocity profiles .
the first large synthetic spectral maps of molecular transitions computed by solving the non  lte radiative transfer through the density and velocity data cubes obtained as the numerical solution of the mhd equations were presented by padoan et al .
( 1998 ) , based on juvela s radiative transfer code ( juvela 1997 ) , and were used in a number of works ( e.g. padoan et al .
1999 , 2000 , 2001 ) .
another new radiative transfer code has also been used more recently to generate synthetic spectral maps from mhd simulations ( ossenkopf 2002 ) .
a number of statistical methods have been proposed to compare numerical models of turbulence with large spectral maps of molecular clouds ( see for example scalo 1984 ; kleiner & dickman 1985 , 1987 ; stutzki & gusten 1990 ; gill & henriksen 1990 ; houlahan & scalo 1992 ; hobson 1992 ; langer , wilson , & anderson 1993 ; williams , de  geus & blitz 1994 ; miesch & bally 1994 ; miesch & scalo 1995 ; lis et al .
1996 ; blitz & williams 1997 ; heyer & schloerb 1997 ; stutzki et al . 1998 ; miesch , scalo & bally 1999 ; falgarone et al . 1994 ; padoan et al . 1999 ; mac low & ossenkopf 2000 ; bensch , stutzki & ossenkopf 2001 ) .    in this work
we apply the spectral correlation function ( scf ) method , proposed by rosolowsky et al .
( 1999 ) and further developed in padoan , rosolowsky & goodman ( 2001 ) , to a number of observational and synthetic spectral maps .
we show that the slope and normalization of the scf of observational maps correlate with the spectral line width .
theoretical models of molecular clouds should therefore yield synthetic spectral maps reproducing such correlations , but not all of them can .    in the next section
we briefly define the scf , and in   3 we present the observational data used in this work .
the computation of the theoretical models and synthetic spectral maps is presented in   4 .
results from numerical models are compared with the observational data in   5 and are discussed in   6 .
conclusions are drawn in   7 .
the spectral correlation function ( scf ) measures the spatial correlation of spectral line profiles within a spectral map .
it is sensitive to the properties of both the gas mass distribution and the gas velocity field ( rosolowsky et al .
1999 ; padoan , rosolowsky & goodman 2001 ; padoan et al .
2001 ; ballesteros  paredes , vazquez  semadeni & goodman 2002 ) .
let @xmath2 be the antenna temperature as a function of velocity channel @xmath3 at map position @xmath4 .
the scf for spectra with spatial separation @xmath5 is : @xmath6 where the average is computed over all map positions @xmath4 .
@xmath7 is the scf uncorrected for the effects of noise , @xmath8 ^ 2 } { \sigma_vt(\vecr , v)^2+\sigma_vt(\vecr+\vecdr , v)^2 } } \right\rangle _ { \vecdr } , \label{2}\ ] ] where the average is limited to separation vectors @xmath9 with @xmath10 , and @xmath11 is the scf due to noise alone , @xmath12 and @xmath13 is the `` spectrum quality '' ( see discussion in padoan , rosolowsky & goodman 2001 ) .
@xmath13 is defined as the ratio of the rms signal within a velocity window @xmath14 and the rms noise , @xmath15 ( over all velocity channels ) , @xmath16 where @xmath17 is the width of the velocity channels .    in the present work
we compute the scf of both observational and synthetic spectral maps , obtained by computing the radiative transfer through the three  dimensional density and velocity fields of numerical simulations of super  sonic mhd turbulence .
the result is typically a power law for @xmath18 that extends up to a separation @xmath5 comparable to the map size , reflecting the self
similarity of super  sonic turbulence ( padoan , rosolowsky & goodman 2001 ) .
the power law behavior is sometimes interrupted at an intermediate scale , possibly suggesting the presence of a physical mechanism limiting the inertial range of turbulence .
an example of a scf that defines an intermediate scale is the scf of the hi survey of the large magellanic cloud ( lmc ) by kim et al .
( 1998 , 1999 ) .
padoan et al .
( 2001 ) have recently been able to map the gas disk thickness of the lmc , assuming it is related to the intermediate scale defined by the break in the scf power law .
the absolute value of @xmath18 at any @xmath19 and the slope of the @xmath18 power law for any given region depends on which molecular tracer is used ( padoan , rosolowsky & goodman 2001 ) .
transitions probing higher gas density produce more fragmented integrated intensity maps than transitions probing lower gas density , and their scf is therefore steeper . in order to compare the scf of observational and synthetic maps
it is therefore important to solve the radiative transfer through the model density and velocity fields accurately for the same molecular transition that is observed .    in this work
our aim is to compute the scf of observational data in order to provide constraints for theoretical models .
the best constraints come from computing the scf of spectral maps of a specific molecular transition over a large range of line width and linear size .
observationally , small scale and narrow line width objects are usually mapped out with high density tracers , while larger objects are instead usually probed with lower density tracers .
@xmath0co provides a good compromise , since it is the only molecule for which very large maps containing thousands of spectra have been obtained with a significant range of resolution . in this work
we have therefore chosen to use observational and synthetic maps of the j=10 line of @xmath0co .
we have used 11 @xmath0co maps . for each map ,
we have listed in table  1 the approximate size , the distance , the rms velocity computed as the standard deviation of the line profile averaged over the whole map , the telescope beam size , the spatial sampling , the width of the velocity channels and the spectral quality defined in the previous section .
smaller maps have been obtained from portions of the maps of the taurus , perseus and rosette molecular cloud complexes and the scf has been computed for each of them .
the position of these smaller maps within the molecular cloud complexes is shown in figures  1 and 2 .
they have been called t1 to t7 in taurus , p1 to p5 in perseus , r1b in the rosette molecular cloud map by blitz & stark ( 1986 ) and r1 and r2 in the rosette molecular cloud map by heyer et al .
( 2001 ) .
the scf of each map has been approximated with a power law , over the range of spatial separations where a power law fit is relevant .
for each power law fit we compute its slope , , and its absolute value at 1  pc , : @xmath20 the values of , , @xmath1 ( the line of sight rms velocity ) and the galactic coordinates of the center of each map are given in table  2 .
the scf of maps of molecular cloud complexes and some smaller regions are shown in figure  [ fig5 ] .
we solve the compressible mhd equations in a staggered mesh of @xmath21 computational cells , with volume centered mass density and thermal energy , face centered velocity and magnetic field components , edge centered electric currents and electric fields and with periodic boundary conditions .
the code uses shock and current sheet capturing techniques to ensure that magnetic and viscous dissipation at the smallest resolved scales provide the necessary dissipation paths for magnetic and kinetic energy .
a more detailed presentation of the numerical method can be found elsewhere .
( padoan & nordlund 1999 ) .    for the purpose of the present work we have computed numerical solutions of the mhd equations using an isothermal equation of state , and a random driving force . in all experiments , the initial density is uniform , and the initial velocity is random .
we generate the velocity field in fourier space , and we give power , with a normal distribution , only to the fourier components in the shell of wave - numbers @xmath22 . we perform a helmholtz decomposition , and use only the solenoidal component of the initial velocity . however , a compressional component of the velocity field develops almost immediately due to the flow compressibility .
the external driving force is generated on large scales in the same way as the velocity field .
the initial magnetic field is uniform , and is oriented parallel to the @xmath23 axis : @xmath24 .    because of the limited numerical resolution we have chosen not to model the collapse of turbulent density fluctuations .
gravity has therefore been neglected .
we have recently started to compute turbulent self  gravitating flows with a numerical mesh of 500@xmath25 cells .
results of the analysis of these larger simulations including self  gravity will be presented in future works .
we have run a number of mhd simulations in a @xmath21 computational mesh , with periodic boundary conditions .
the simulations are intended to describe the turbulent dynamics in the interior of molecular clouds .
the two most important numerical parameters in the models are the rms sonic and alfvnic mach numbers , @xmath26 and @xmath27 .
the rms sonic mach number is here defined as the ratio of the rms flow velocity and the speed of sound .
the alfvnic mach number is defined as the ratio of the rms flow velocity and the alfvn velocity , @xmath28 , where b is the volume  averaged magnetic field strength .
all the models used in this work have @xmath29 , except for model e that has @xmath30 , according to the suggestion that the dynamics of molecular clouds is essentially super  alfvnic ( padoan & nordlund 1999 ) .
our numerical simulations conserve magnetic flux , and so the volume averaged magnetic field is constant in time . as a consequence , also
the value of @xmath27 as defined above remains constant .
however , the value of @xmath31 grows with time ( until equilibrium is reached ) due to compression and stretching of magnetic field lines ( see padoan & nordlund 1999 ) .
if we define the alfvn velocity using the rms value of the magnetic field strength , instead of its volume average , then the typical alfvnic mach number in our super  alfvnic runs is @xmath32 , because of the formation of regions with large value of magnetic field strength ( mainly dense regions , as found in observations ) .
the sonic mach number of observed turbulent motions in molecular clouds is @xmath33 on the scale of several parsecs , and decreases toward smaller scale .
the turbulent velocity becomes comparable to the speed of sound only on very small scale , @xmath34  pc . in order to study the effect of the sonic mach number on the scf
, we have computed mhd models with different values of @xmath26 , @xmath35 , 5 , 2.5 , 1.25 and 0.625 .
each model has been run for approximately six dynamical times ( the dynamical time is here defined as the ratio of half the size of the computational box and the rms flow velocity ) , in order to achieve a statistically relaxed state , independent of the initial conditions .    the velocity and density fields from the final snapshot of each model have been used to compute @xmath0co @xmath36 spectra , solving the radiative transfer with a non  lte monte carlo code (   4.2 ) . while the mhd calculations are independent of the physical value of the average gas density , the size of the computational mesh ( or the column density ) and the kinetic temperature , these physical parameters are necessary inputs for the radiative transfer calculations .
the models are scaled to physical units assuming a value for i ) the kinetic temperature , @xmath37 , that determines the physical unit of velocity ( the numerical unit of velocity is the speed of sound ) ; ii ) the average gas density , @xmath38 ; iii ) the size of the computational box , @xmath39 .
for all models we have assumed @xmath40  k , typical of molecular clouds .
the dependence of observed average gas density and cloud size on the observed rms turbulent velocity ( or sonic mach number , assuming a constant value of @xmath37 ) is well  approximated by empirical larson type relations ( larson 1981 ) .
however , the size  velocity relation has a large intrinsic scatter ( falgarone , puget & perault 1992 ) , and both the size  velocity and density  size relations have been criticized by several authors ( loren 1989 ; kegel 1989 ; scalo 1990 ; issa , maclaren & wolfendale 1990 ; adler & roberts 1992 ; vazquez - semadeni , ballesteros - paredes & rodriguez 1997 ; ostriker , stone & gammie 2001 ; ballesteros  paredes & mac low 2002 ) . for these reasons , we scale the mhd models in four different ways . these four sets of models are all based on the same five mhd turbulence models and differ from each other only in the way they are rescaled to physical units when computing the radiative transfer .
models a1 to a5 and b1 to b5 have all the same value of the average density , @xmath41  @xmath42 .
models a1 to a5 have all the same size @xmath43  pc and column density @xmath44  @xmath45 ; models b1 to b5 have @xmath46  pc and @xmath47  @xmath45 .
models a1r to a5r and b1r to b5r are rescaled using the larson type relations : @xmath48 where a temperature @xmath40  k is assumed , and @xmath49 that is equivalent to a constant mean surface density .
models a1r to a5r have the same column density as models a1 to a5 , that is @xmath50  @xmath42 in equation ( [ larson2 ] ) ; they also have sizes @xmath5110 , 2.5 , 0.625 , 0.156 and 0.039 pc respectively , which implies @xmath52 in equation ( [ larson1 ] ) .
models b1r to b2r have the same column density as models b1 to b5 , that is @xmath53  @xmath42 in equation ( [ larson2 ] ) ; they have sizes @xmath5120 , 5 , 1.25 , 0.31 and 0.078 pc respectively , which implies @xmath54 in equation ( [ larson1 ] ) .
finally , the equipartition model ( model e ) has been computed only for one value of the rms sonic mach number , @xmath35 .
it is rescaled to the larson type relations only once , for a size of 10  pc and a column density of @xmath44  @xmath45 . for this model
we have computed spectral maps along 5 different directions , three orthogonal to the faces of the numerical mesh , as in the other experiments , and two along diagonal directions , at an angle of 54.7@xmath55 with the average magnetic field ( @xmath23 axis ) .
maps from diagonal directions sample lines of sight of different length at different map positions ( longer at the central position than near the corners ) . however , the number of computational cells along each line of sight is on the average even larger than in maps from orthogonal directions , since a diagonal line of sight often cuts through the computational cells away from their center ( close to their corners ) .
furthermore , the maps are computed only for a region of size equivalent to that of maps from orthogonal directions ( @xmath56 cells ) , eliminating the corners of the computational mesh . as a result ,
only a few percent of the spectra are generated from lines of sight sampling less than 50 computational cells .
a fraction of the the lines of sight close to the map edges are nevertheless shorter than the energy injection scale ( approximately half the size of the computational mesh ) .
this may introduce a bias toward smaller line width , since velocity differences are expected to grow with increasing distances .
this bias or its effect on the scf should be small , since our results seem to vary smoothly as a function of the angle between the line of sight and the direction of the average magnetic field .
the models a4 , a5 , b3 , b4 and b5 have velocity dispersion significantly smaller than found observationally at the scale of 5  pc .
they are not used here to test the validity of models with such low velocity dispersion , but rather to test the ability of the scf method to rule them out as poor description of molecular cloud turbulence .    in order to test the ability of the scf to rule out unphysical models ,
we have also computed two stochastic models , s2 and s4 . in both models
the density field is a random field with a log  normal probability distribution function , and a power law power spectrum with power law exponent equal to -1 ( the approximate value found in our mhd models ) .
the velocity field is generated as a gaussian field , also with power law power spectrum .
the power law exponent of the velocity field power spectrum is -2 ( close to the actual value in the mhd models ) in model s2 and -4 in model s4 .
for the purpose of computing the radiative transfer and the synthetic spectral maps , both models have been scaled to a physical size @xmath46  pc and a column density @xmath57  @xmath45 .
these two stochastic models are unphysical in the sense that they are not solutions of the fluid equations .
statistical properties such as the power spectrum and the probability density function of density and velocity may be similar to those of flows obtained by solving the fluid equations , but their phase correlations are unphysical .
this is in part illustrated by the fact that these stochastic models look clumpy , rather than filamentary as real clouds and mhd models .
furthermore , their velocity and density fields can not be self
consistent because they are computed independently of each other .
it is shown below that the scf method can indeed rule out these unphysical models .
the sonic rms mach number , @xmath26 , the average gas density , @xmath38 , and the physical size , @xmath39 , of the different models used for the radiative transfer computations are given in table  3 .
the radiative transfer calculations were carried out with a monte carlo program which is a generalization of the one - dimensional monte carlo method ( bernes 1979 ) into three dimensions .
the model cloud is divided into small , cubic cells in which physical properties are assumed to be constant .
the discretization allows the inclusion of arbitrary kinetic temperature and molecular abundance variations .
however , in the present calculations , the temperature and relative abundances are kept constant .
the 2.73  k cosmic background is used as the external radiation field .
there are important differences between our program and the normal monte carlo method , and some principles of the implementation are given below .
a detailed description is given elsewhere ( juvela 1997 ) .
in the basic monte carlo method radiation field is simulated with photon packages , each representing a number of real photons .
the packages are created at random velocities at random locations and sent toward random directions .
each package is followed through the cloud and interactions between photons and molecules are counted .
later this information is used to solve new estimates for the level populations of the molecules .    in our method
the radiative transfer is simulated along random lines going through the cloud .
initially , as a photon package enters the cloud it contains only background photons .
as the package goes through a cell in the cloud some photons emitted by this cell are added to the package and , in particular , the number of photons absorbed within the emitting cell is calculated explicitly .
this becomes important when cells are optically thick and , compared with normal monte carlo simulation , ensures more accurate estimation of the energy transfer between cells . in our program
each simulated photon package represents intensity of all simulated transitions and doppler shifts at the same time .
the number of individual photon packages is correspondingly smaller , and in the present case we use 240000 photon packages per iteration .
the lines are divided into 70 fixed velocity channels .
there is no noise associated with random sampling of doppler shifts .
the simulated velocity range was adjusted according to the velocity range found in the model clouds .
the channels are narrow compared with the total line widths and smaller than or equal to the smallest intrinsic line widths in the cells .
the velocity discretization is therefore not expected to affect the results of the calculations .
the density and velocity fields from the mhd simulations are sampled on a numerical mesh of @xmath21 cells . to speed up the radiative transfer calculations
the density and velocity fields were rebinned into a mesh of @xmath58 cells by linear interpolation .
the velocity dispersion between neighboring cells in the original @xmath21 data cube was used to approximate the turbulent line width within each cell of the new @xmath58 data cube .
this velocity dispersion should apply to a scale slightly larger than the size of the cells in the @xmath58 mesh .
however , this is approximately compensated by the fact that numerical dissipation in the mhd simulations decreases significantly the velocity dispersion on very small scale , below the actual turbulent inertial  range value at that scale .    on each iteration
new level populations are solved from the equilibrium equations and iterations are stopped when the relative change is below @xmath592.0@xmath6010@xmath61 in all cells . only the six lowest levels were tested for convergence .
the relative changes tend to be largest on the upper levels where the level populations become very small and , on the average , the convergence of the relevant first energy levels is much better than the quoted limit . the total number of energy levels included in the calculations was nine , a number clearly sufficient in case of excitation temperatures below 10k .
the collisional coefficients were taken from flower & launay ( 1985 ) and green & thaddeus ( 1976 ) .
the final level populations were used to calculate maps of 90@xmath6290 spectra toward three directions perpendicular to the faces of the mhd data cube . for the equipartition model e spectra
were calculated also along two diagonal directions . in these cases
the maps of 90@xmath6290 spectra do not extend over the whole projected cloud area .
each spectrum corresponds to the intensity calculated along one line of sight ( spectra are not convolved with a larger beam ) .
the spectra contain 60 velocity channels as in the monte carlo simulation .
the results were compared with spectra calculated assuming lte conditions .
the comparison showed that for typical physical conditions found in molecular clouds the lte assumption would be unsuitable ( padoan et al .
most observational and theoretical spectral maps yield a scf that can be approximated by a single power law within a range of spatial separations , often spanning over an order of magnitude . from each power law
fit we compute its slope , , and its value at 1  pc , , defined as in ( [ alphadef ] ) .
we also compute the value of the velocity dispersion , @xmath1 , from each map , measured as the standard deviation of the @xmath0co @xmath36 spectrum averaged over the entire map .
the values of ,  and @xmath1 computed from the observational maps are given in table  2 , while the values of the same quantities from the theoretical models are given in table  3 .
every model provides three sets of values , because spectral maps have been computed using three orthogonal directions for the line of sight .
each group of three sets of values can be interpreted as the same model cloud being `` observed '' from different directions , or as three different model clouds with comparable rms velocity . in the equipartition model , e , the rms velocities inferred from different directions are very different from each other , the largest rms velocity being found in the direction parallel to the mean magnetic field ( along the @xmath23 axis ) , and the lowest in the directions perpendicular to the magnetic field .
for this model we have also computed spectral maps from two more lines of sight , corresponding to diagonal directions across the computational box .
figure  [ fig6 ] shows the scf of the model a1r , with line of sight parallel to the direction of the mean magnetic field .
the scf of the equipartition model e is also shown for four lines of sight , two diagonal , one parallel to the direction of the mean magnetic field and one perpendicular to it to it .
the figure shows that the scf of model e is very sensitive to the line of sight , due to the large variations of the rms velocity in different directions relative to the mean magnetic field .
the velocity dispersion relative to the speed of sound , or the value of the sonic mach number , is the most important physical parameter characterizing the nature of the turbulence .
we therefore study the dependence of the scf on the turbulent velocity dispersion ( or the rms sonic mach number , ) and propose to use this dependence to test theoretical models against the observational data .    in figure  [ fig7 ] the slope of the scf
is plotted against the line of sight velocity dispersion .
the top panels show the models of constant size and constant average gas density ( models a1 to a5 -left , and b1 to b5 -right ) ; the bottom panels show the models rescaled according to larson type relations ( models a1r to a5r -left and b1r to b5r -right ) .
the observations indicate a strong correlation between  and @xmath1 , over an order of magnitude in @xmath1 .
a least square fit to the observational data gives : @xmath63 and for the super  sonic and super  alfvnic models rescaled with larson type relations ( a1r to a5r ) : @xmath64 consistent with the observational result ( the uncertainty in the exponent is the standard deviation from the least square fit ) .
the corresponding models not scaled with the larson type relations ( a1 to a5 ) are also indistinguishable from the observational result ( see figure  [ fig7 ] top left panel ) .
models of type b ( right panels of figure  [ fig7 ] ) have instead values of significantly smaller than the average ones from the observational data .
the stochastic models s2 and s4 are indistinguishable from each other ; they are also totally inconsistent with the empirical @xmath1 relation , which allows them to be ruled out as invalid by the scf method .
finally , the equipartition model e provides values that are consistent with the observations , and comparable to the super  alfvnic models , apart from a larger scatter of values between different lines of sight .
we interpret the increase of  with @xmath1 as a consequence of the increasing compressibility of the turbulent flow ( @xmath1 is roughly proportional to the rms sonic mach number of the flow because the temperature in all the models is @xmath40  k , and approximately the same in the observed regions ) .
the value of  is in general found to grow with increasing density contrast , probably due to the increasing concentration of the mass along the line of sight around one or few dense cores , which helps decorrelating the spectra from each other .
the value of  is plotted against @xmath1 in figure  [ fig8 ] .
the top panels show the models of constant size and the bottom panels the models scaled with the larson type relations , as in figure  [ fig7 ] .
the values of  and @xmath1 from the observational maps are weakly correlated , with  slightly increasing with increasing @xmath1 : @xmath65 a tight correlation is instead found in the models of constant size ( top panels of figure  [ fig8 ] ) , with  decreasing with increasing @xmath1 .
this inconsistency between the models and the observations is most likely due to the fact that molecular clouds of 5 to 20 pc of size ( as assumed by these models ) are never found with velocity dispersion as low as assumed in models a4 , a5 and b3 , b4 and b5 .
the bottom panels of figure  [ fig8 ] show that the inconsistency is in fact mostly resolved as soon as the model sizes are scaled according to the larson type relation . for the models
a1r to a5r we obtain : @xmath66 if models with @xmath67  km / s were not included ( justified by the absence of such low velocity dispersions in the observational sample ) , the slope of the least square fit would be @xmath68 , fully consistent with the observations .
the equipartition model yields values of  and @xmath1 consistent with the observations in all directions , but the one parallel to the mean magnetic field .
it could be concluded that either none of the observed objects has a significant component of the magnetic field along the line of sight , or that all of them have a magnetic field weaker than predicted by the equipartition model , consistent with the super  alfvnic models .    in figure  [ fig9 ]
we have plotted observations and models on the   plane .
the constant size models are again inconsistent with the observations , as is expected since the observational maps span a large range of scales .
when the models are scaled according with the larson type relations and the realistic average column density of @xmath69  @xmath45 ( myers & goodman 1988 ) , the observed scatter in the   plane is reproduced .
the trend of the absolute value of  to increase with for large values of both of them is also reproduced , between models with rms mach 5 and 10 ( a2r and a1r respectively ) ; however , models with rms mach of 20 or 30 would be necessary to fit the  and  values measured for the rosette molecular cloud , which can be appropriately resolved only with a numerical resolution in excess of @xmath70 computational cells .    while the stochastic models s2 and s4 are only marginally inconsistent with the observations in this plot , the line of sight parallel to the direction of the mean magnetic field and one of the two diagonal lines of sight in the equipartition model e are again inconsistent with the observational data .
the scf has been proposed as a statistical tool to test the validity of theoretical models describing the structure and dynamics of star forming clouds ( rosolowsky et al . 1999 ) . in a previous work we improved the scf method by studying its dependence on spatial and velocity resolution and on instrumental noise ( padoan , rosolowsky & goodman 2001 ) .
here we have applied that improved scf to a number of large @xmath0co maps of molecular cloud complexes and obtained empirical correlations that can be used to test theoretical models . of the theoretical models we have computed some
compare well with the empirical correlations and some do not , which shows that the scf can be used as an effective tool to rule out inappropriate or unphysical models .
the empirical correlations we have obtained relate the values of ,  and @xmath1 with each other .
the @xmath1 correlation rules out the unphysical stochastic models ( s2 and s4 ) . such models were found to produce spectral line profiles similar to observational ones by dubinski , narayan & phillips ( 1995 ) .
they have also been used as models of the density field in molecular clouds by stutzki et al .
( 1998 ) and to calibrate their principal component analysis by brunt & heyer ( 2002 ) .
the scf @xmath1 correlation shows that stochastic models are inappropriate to describe the structure of molecular cloud complexes .
models not scaled with larson type relations ( a1a5 , b1b5 ) and models with larger  than  average column density ( b1b5 , b1r
b5r ) have also been compared with the empirical scf correlations to show that incorrectly scaled models are readily ruled out by the scf method .
the @xmath1 and the
correlations do not favor the model with equipartition of kinetic and magnetic energies ( model e ) .
such model yields too small values of  or too large values of  compared with the observational data , when seen in the direction parallel to the average magnetic field . of the two diagonal directions , one is consistent with the data and the other is not .
a possible interpretation is that none of the observed regions has an average magnetic field oriented close to the direction of the line of sight .
the equipartition model starts to be inconsistent with the observational data when seen along the diagonal directions , at an angle of 54.7@xmath55 to the average magnetic field .
the line of sight should be within such an angle to the magnetic field in approximately 40% of the cases , assuming random orientation of the average magnetic field in the observed regions .
an alternative interpretation is that all the observed regions have an average magnetic field strength smaller than in the equipartition model , and consistent with super
alfvnic conditions .
the super  alfvnic models rescaled with larson type relations are in fact able to reproduce the empirical scf correlations . however , the total number of truly independent directions on the sky in the present observational sample is still small .
more regions should be studied to rule out the equipartition model based on the scf results .
the analysis of the mhd models could in principle give different results if self  gravity was taken into account .
however , the introduction of self
gravity is not expected to decrease the value of  and increase the value of , as necessary to make the equipartition model consistent with the observational   correlation .
the main effect of self
gravity is the collapse of the densest regions , increasing the density contrast beyond the level due to the turbulence alone .
this could slightly increase the value of  because we interpret the increase of  with @xmath1 in the mhd models as due to the increased compressibility of the turbulent flow .
an increase in the value of  is not expected because the local collapse of dense cores can not increase the correlation between spectra at large distances ( for a given value of an increase in  would correspond to an increase in the scf at large spatial separation ) . nevertheless , the effect of self
gravity should be tested by including it in the numerical solution of the mhd equations .
the numerical resolution should also be larger than in the present work to resolve the initial phase of the gravitational collapse of dense cores .
we have only recently started to compute self
gravitating flows in a numerical mesh of @xmath71 cells , and their analysis will be presented in future works .
padoan & nordlund ( 1997 , 1999 ) have proposed that the dynamics of molecular clouds on large scales is consistent with super  alfvnic turbulence and inconsistent with the equipartition model . in numerical simulations of super - alfvnic turbulence
the average magnetic energy grows with time , even if flux is conserved ( the average magnetic field is constant ) .
the magnetic field strength is increased locally mainly in regions of compression in super  sonic turbulence , and in part by stretching of field lines .
even if initial conditions are such that the turbulence is highly super - alfvnic , magnetic pressure is often larger than thermal pressure in the postshock gas , due to the amplification of the magnetic field components perpendicular to the shock direction .
equipartition of dynamic pressure , @xmath72 , and magnetic pressure , @xmath73 , is therefore achieved locally , but not necessarily over the whole flow .
for example , in super  sonic and super - alfvnic runs at a resolution of @xmath74 , the ratio of volume  average magnetic and dynamic pressures relaxes at a value @xmath75 , starting from initial conditions with @xmath76 ( padoan et al .
comparable values are found in the numerical experiments used in this work .
the amplification of the magnetic field by the turbulence therefore does not alter the super  alfvnic character of the flow .
the correlation between local magnetic field strength and gas density in super  sonic and super  alfvnic turbulence has a very large scatter , and a well defined upper envelope with @xmath77 , both consistent with the observational data ( padoan & nordlund 1997 , 1999 ; ostriker , stone & gammie 2001 ; passot & vazquez  semadeni 2002 ) .
the largest values of the magnetic field strength are generally found in dense cores , but some dense cores may have a relatively weak magnetic field .
however , dense cores assembled by turbulent shocks are not expected to have internal super  alfvnic turbulence , because of the dissipation of kinetic energy in the shocks and of the amplification of the magnetic field in the compressed gas .
observational evidence for an approximate equipartition of turbulent and magnetic energy in dense cores would therefore not be inconsistent with the super  alfvnic character of the large scale flow that assembles them .
the comparison between our theoretical models and the observational data could be improved if more regions with sub  sonic turbulence were available in the observational sample .
small velocity dispersion is found in small objects , according to larson s velocity  size relation , or to the power spectrum of turbulence .
the spatial resolution in single dish surveys is typically too low to sample a small object ( fraction of a parsec ) with a very large spectral maps ( several thousands of spectra ) .
the only exceptions in the observational sample used in this work are l1512 and l134a .
these two large maps of nearby clouds with very low velocity dispersion were obtained by falgarone et al .
( 1998 ) as part of their iram key project , focused on regions of relatively low column density at the edges of molecular cloud complexes .
maps of large regions with very large velocity dispersion are instead more easily obtained from observations than in numerical simulations . assuming a gas kinetic temperature of the order of 10  k
, a line of sight ( one dimensional ) velocity dispersion in excess of 2  km / s corresponds to a sonic rms mach number of the flow @xmath78 . in the present work we have not computed numerical flows with @xmath79 , since that would require a larger numerical resolution ( the density contrast grows linearly with the alfvnic mach number and therefore with the value of @xmath26 if the average magnetic field strength is not varied ) . for this reason the models do not reach the largest values of @xmath1 ,  and  obtained from the observations ( from the maps of the rosette molecular cloud complex ) .
the progression of models toward increasing values of @xmath26 suggests that a model with @xmath78 would likely fit the observed   values found in the rosette molecular cloud complex , where the observed velocity dispersion is in excess of 2  km / s .
this is illustrated in figure  [ fig10 ] .
the top panel of figure  [ fig10 ] shows the   for the observational data .
the shaded area shows the range of values covered by the theoretical models a1r to a5r .
the bottom panel shows the same plot for the model a1r to a5r .
each diagonal segment connects the values for the three directions of each model .
the values of  and @xmath39 of the models are also given in the plot .
the arrow marks the direction of increasing  suggesting that models with @xmath80 may fit the observations with the largest velocity dispersion .
in the present work we have computed the spectral correlation function ( scf ) of spectral maps of molecular cloud complexes and regions within them , observed in the j=10 transition of @xmath0co .
we have found that the scf is a power law over approximately an order of magnitude in spatial separation .
the power law slope of the scf , , its normalization , , and the spectral line width averaged over the whole map , @xmath1 , have been computed for all the observational maps .
we have obtained empirical correlations between these quantities and have proposed to use them to test the validity of theoretical models of molecular clouds .
theoretical models of spectral line maps have been generated by computing the radiative transfer through the numerical solutions ( density and velocity fields ) of the magneto - hydrodynamic ( mhd ) equations , for turbulent flows with different values of the rms sonic and alfvnic mach numbers , and also through stochastic density and velocity fields with different power spectra .
super - alfvnic mhd models rescaled according to larson type relations are in the best agreement with the empirical correlations .
unphysical stochastic models are instead ruled out .
mhd models with equipartition of magnetic and kinetic energy of turbulence do not reproduce the observational data when their average magnetic field is oriented approximately parallel to the line of sight .
finally , mhd models not rescaled according to larson type relations are also inconsistent with the observational data .
we can not exclude the possibility that different physical models for the dynamics of molecular clouds , or even unphysical models , that we have not tested here , would satisfy the empirical correlations found in this work .
however , we have shown that the scf method is able to rule out certain unphysical or incorrectly scaled models . reproducing these scf results
should be considered as a necessary ( but not sufficient ) condition for the validity of theoretical models describing the structure and the dynamics of molecular clouds .
models for which the scf or similar statistical tests can not be computed to allow a quantitative comparison with observed spectral maps can not be legitimately evaluated .
the comparison between theory and observations presented in this work requires significant computational resources .
numerical simulations of three dimensional turbulent flows must be run at large resolution and the radiative transfer has to be computed in three dimensions in order to generate synthetic spectral maps of the observed molecular transitions .
the type of models and the physical parameters investigated in this work are therefore limited to a few significant cases .
future work should investigate the scf of a larger variety of models , including different magnetic field intensities and flows with gravitationally collapsing cores .
we are grateful to eve ostriker and jim stone for helpful comments on our model
data comparison .
the referee report by enrique vazquez
semadeni has also contributed to improve this work .
this work was supported by an nsf galactic astronomy grant to ag .
the work of pp was partially performed while pp held a national research council associateship award at the jet propulsion laboratory , california institute of technology .
mj acknowledges the support of the academy of finland grant no .
1011055 .                                                                                * figure captions : * + * table [ t1 ] * main parameters of the observed spectral maps : approximate size , distance , rms velocity over the whole map , telescope beam , spatial sampling , velocity channel width , average spectrum quality and bibliographic reference .
+ * table [ t2 ] * spectral line width averaged over the whole map , @xmath1 , power law slope of the scf ,  and scf normalization , , galactic longitude , @xmath81 , and galactic latitude , @xmath82 , of the center of all the observed maps and selected regions within them . + * table [ t3 ] * first three columns from the left : model name , rms sonic mach number of the flow and physical size of the computational mesh .
following columns : line of sight velocity dispersion , scf slope and scf normalization , repeated for the three orthogonal directions for which synthetic spectral maps have been computed in each model .
values for the diagonal directions of model e are not given ( they are within the ranges of values covered by the other three directions parallel and perpendicular to the mean magnetic field ) . + *
figure [ fig1 ] : * top panel : velocity integrated intensity map of the perseus molecular cloud complex in the j=1 - 0 transition of @xmath0co ( padoan et al .
bottom panel : same as top panel , but for the taurus molecular cloud complex ( mizuno et al .
smaller regions within the maps where the scf has also been computed are highlighted .
+ * figure [ fig2 ] : * same as in figure  [ fig1 ] , but for the rosette molecular cloud complex .
top panel from heyer et al .
( 2001 ) ; bottom panel from blitz & stark ( 1986 ) . + * figure [ fig5 ] : * top left panel : the scf averaged over the entire map of the perseus , rosette and taurus molecular cloud complexes .
solid lines are least square fits to the power law sections of the scf .
the exponents @xmath83 of the power law fits are also given in the figure .
top right panel : scf of the whole map of the perseus molecular cloud complex and of smaller regions within the same map .
bottom left panel : scf of the taurus molecular cloud complex and of smaller regions within the same map .
bottom right panel : scf of pvceph and hh300 .
+ * figure [ fig6 ] : * the scf computed from mhd models .
asterisk symbols are for the super - alfvnic model a1r in the @xmath84 direction ( parallel to the mean magnetic field ) .
diamond symbols are for the equipartition model in the @xmath85 and @xmath84 direction ( perpendicular and parallel to the average magnetic field direction respectively ) and along two diagonal directions ( @xmath86 and @xmath87 ) . the slope of the scf increases with increasing rms velocity .
the scf is therefore weakly dependent on the direction of the line of sight for the super - alfvnic model , while it is much steeper in the direction parallel to the magnetic field ( larger rms velocity ) than in the perpendicular direction in the equipartition model . + * figure [ fig7 ] : * scf slope versus velocity dispersion .
the top panels show the models of constant size and constant average gas density as asterisks ( models a1 to a5 left , b1 to b5 right ) ; the bottom panels show the models rescaled according to larson type relations as asterisks ( models a1r to a5r left , b1r to b5r right ) .
observational values are shown as squares , the equipartition model as triangles and the stochastic models as diamonds .
+ * figure [ fig8 ] : * scf value at 1  pc versus velocity dispersion .
different panels show different models as in figure  [ fig7 ] .
symbols are also as in figure  [ fig7 ] .
+ * figure [ fig9 ] : * scf slope versus scf value at 1  pc . symbols and panels as in figure  [ fig7 ] .
+ * figure [ fig10 ] : * top panel : values of  and from the observations .
some of the symbols are labeled with the region name .
the shaded area shows the range of values covered by the models a1r to a5r .
bottom panel : same shaded area as in the top panel .
diagonal segments shows the range of values of  and  for the three directions of each model .
the rms sonic mach number of the corresponding model is given on the right hand side of each segment , while the value of the linear size is given on the left hand side .
the arrow indicates the progression of models toward larger values of sonic mach number , .
lcc|ccc|ccc|ccc & & & & x & & & y & & & z & + model & & 0 & & & & & & & & & +   + a1 & 10.0 & 5 & 1.13 & 0.30 & 0.37 & 1.36 & 0.32 & 0.34 & 1.24 & 0.35 & 0.31 + a2 & 5.0 & 5 & 0.56 & 0.25 & 0.46 & 0.71 & 0.30 & 0.39 & 0.67 & 0.28 & 0.42 + a3 & 2.5 & 5 & 0.30 & 0.19 & 0.55 & 0.38 & 0.24 & 0.48 & 0.36 & 0.21 & 0.52 + a4 & 1.2 & 5 & 0.17 & 0.13 & 0.66 & 0.17 & 0.15 & 0.64 & 0.18 & 0.16 & 0.61 + a5 & 0.6 & 5 & 0.12 & 0.13 & 0.68 & 0.11 & 0.10 & 0.71 & 0.12 & 0.14 & 0.66 + b1 & 10.0 & 20 & 1.21 & 0.26 & 0.64 & 1.42 & 0.27 & 0.61 & 1.33 & 0.27 & 0.59 + b2 & 5.0 & 20 & 0.61 & 0.22 & 0.70 & 0.77 & 0.26 & 0.65 & 0.72 & 0.24 & 0.67 + b3 & 2.5 & 20 & 0.33 & 0.18 & 0.73 & 0.41 & 0.23 & 0.68 & 0.39 & 0.20 & 0.71 + b4 & 1.2 & 20 & 0.18 & 0.14 & 0.77 & 0.19 & 0.15 & 0.77 & 0.21 & 0.17 & 0.75 + b5 & 0.6 & 20 & 0.13 & 0.13 & 0.79 & 0.13 & 0.11 & 0.81 & 0.14 & 0.14 & 0.78 +   + a1r & 10.0 & 10 & 1.13 & 0.30 & 0.44 & 1.36 & 0.32 & 0.42 & 1.23 & 0.35 & 0.39 + a2r & 5.0 & 2.5 & 0.56 & 0.23 & 0.41 & 0.72 & 0.30 & 0.32 & 0.66 & 0.26 & 0.37 + a3r & 2.5 & 0.62 & 0.30 & 0.18 & 0.41 & 0.37 & 0.21 & 0.34 & 0.35 & 0.20 & 0.37 + a4r & 1.2 & 0.16 & 0.16 & 0.12 & 0.45 & 0.15 & 0.13 & 0.42 & 0.17 & 0.15 & 0.38 + a5r & 0.6 & 0.04 & 0.10 & 0.11 & 0.41 & 0.10 & 0.08 & 0.49 & 0.11 & 0.12 & 0.39 + b1r & 10.0 & 20 & 1.21 & 0.26 & 0.64 & 1.42 & 0.27 & 0.61 & 1.33 & 0.27 & 0.59 + b2r & 5.0 & 5 & 0.60 & 0.22 & 0.55 & 0.76 & 0.23 & 0.49 & 0.69 & 0.20 & 0.53 + b3r & 2.5 & 1.25 & 0.33 & 0.15 & 0.52 & 0.40 & 0.18 & 0.45 & 0.37 & 0.16 & 0.49 + b4r & 1.2 & 0.31 & 0.18 & 0.12 & 0.50 & 0.18 & 0.13 & 0.47 & 0.19 & 0.14 & 0.44 + b5r & 0.6 & 0.08 & 0.12 & 0.10 & 0.46 & 0.12 & 0.09 & 0.52 & 0.13 & 0.12 & 0.42 +   + e & 10.0 & 5 & 0.73 & 0.21 & 0.46 & 0.81 & 0.25 & 0.43 & 1.90 & 0.49 & 0.22 +   + s2 & 10.0 & 20 & 1.37 & 0.15 & 0.52 & 1.29 & 0.15 & 0.53 & 1.31 & 0.14 & 0.53 + s4 & 10.0 & 20 & 1.23 & 0.16 & 0.47 & 1.17 & 0.15 & 0.49 & 1.17 & 0.17 & 0.46 +