there has been much interest in superconductivity with very short coherence length .
this interest is due to its possible relevance to high temperature superconductors ( cuprates , doped bismuthates , iron - based systems , fullerenes ) and also to the several other exotic superconducting materials ( for a review , see @xcite and references therein ) .
it can also give relevant insight into behaviour of strongly bounded fermion pairs on the optical lattices .
the phase separations involving superconducting states have been evidenced in a  broad range of currently intensely investigated materials including iron - pnictides , cuprates and organic conductors ( see for example @xcite and references therein ) .    in our work we will study a  model which is a  simple generalization of the standard model of a local pair superconductor with on - site pairing ( i.  e. the model of hard core bosons on a  lattice @xcite ) to the case of finite pair binding energy .
the hamiltonian considered has the following form : @xmath4 where , , ; @xmath5 ( @xmath6 ) denotes the annihilation ( creation ) operator of an electron with spin at the site @xmath7 , which satisfy canonical anticommutation relations @xmath8 where @xmath9 is the kronecker delta . indicates the sum over nearest - neighbour sites @xmath7 and @xmath10 independently .
@xmath11 will denote the number of nearest neighbours .
@xmath3 is the chemical potential , depending on the concentration of electrons : @xmath12 with and @xmath13 is the total number of lattice sites .
there are two interaction parameters of the model : ( i )  the intersite pair hopping interaction @xmath1 , determining the electron pair mobility and responsible for the long - range superconducting order in the system and ( ii )  the on - site density - density interaction @xmath0 , which contributes ( together with @xmath1 ) to the pair binding energy by reducing ( ) or enhancing ( ) its value . to simplify our analysis
we do not include in hamiltonian ( [ row : ham1 ] ) the single electron hopping term ( @xmath14 ) as well as other inter - site interaction terms .
this assumption corresponds to the situation when the single particle mobility is much smaller than the pair mobility and can be neglected . in the presence of the @xmath14 term , model ( [ row : ham1 ] ) is called the penson - kolb - hubbard ( pkh ) model @xcite and the effects of will be pointed out in section  [ sec : conclusions ] .    model ( [ row : ham1 ] ) can be considered as a  relatively simple , effective model of a  superconductor with local electron pairing . moreover
, the knowledge of the limit can be used as starting point for a perturbation expansion in powers of the hopping @xmath14 and provides a benchmark for various approximate approaches analyzing the corresponding finite bandwidth models .    in terms of pseudospin charge operators hamiltonian ( [ row : ham1 ] ) takes the form of a special xy - model with a  single ion anisotropy in a  transverse field ( in the @xmath11-direction ) @xmath15 where and corresponds to the external magnetic field .
the operators satisfy standard commutation relations for momentum operators : @xmath16 = \pm \hat{\rho}^\pm_i \delta_{ij } , \qquad [ \hat{\rho}_i^{+ } , \hat{\rho}_j^{- } ] = 2\hat{\rho}_i^{z}\delta_{ij},\ ] ] what follows directly from  ( [ row : antycomutation ] ) .
however , one should stress that takes on four values , as for model with doubly degenerated @xmath17 value , whereas , as in the system .
condition ( [ row : condn1 ] ) for the electron concentration reads : @xmath18 it means that magnetic field @xmath19 is determined by a  fixed value of magnetization @xmath20 .
the interactions @xmath0 and @xmath1 will be treated as the effective ones and will be assumed to include all the possible contributions and renormalizations like those coming from the strong electron - phonon coupling or from the coupling between electrons and other electronic subsystems in solid or chemical complexes @xcite .
in such a general case arbitrary values and signs of @xmath0 and @xmath1 are important to consider .
it is notable that formally @xmath1 is one of the off - diagonal terms of the coulomb interaction  @xcite , describing a part of the so - called bond - charge interaction , and it sign of the _ coulomb - driven _
charge exchange is typically negative ( repulsive , ) .
however , the effective attractive of this form ( ) is also possible  @xcite and in particular it can originate from the coupling of electrons with intersite ( intermolecular ) vibrations via modulation of the hopping integral  @xcite , or from the on - site hybridization term in generalized periodic anderson model  @xcite .
the ferromagnetic xy order of pseudospins @xmath21 ( for ) corresponds to the ss phase ( @xmath22-pairing superconducting ) , whereas the antiferromagnetic xy order ( for )  to the @xmath23s phase ( @xmath23-pairing superconducting ) .
however , in the absence of the external field in the xy - direction conjugated with the order parameter ( , for and , for , @xmath24 is a  half of the smallest reciprocal lattice vector ) there is a symmetry between the ferromagnetic and antiferromagnetic cases ( except for an obvious redefinition of the order parameter ) for lattices consisting of two interpenetrating sublattices ( such as for example sc or bcc lattices ) , thus we restrict ourselves to the case in the following .
let us stress that the symmetry is broken by the single electron hopping term ( cf . section  [ sec : conclusions ] ) .    the first analysis of the phase diagram of model ( [ row : ham1 ] ) have been performed by bari  @xcite and ho and barry  @xcite using the variational method in order to examine the instability of the mott insulator to superconductivity for the special case of the half - filled band ( ) .
the effects of diagonal disorder on the critical temperature for and have been also determined  @xcite , arriving at a  satisfactory qualitative interpretation of quite a number of different experiments in amorphous superconductors .    in the analysis we have adopted the variational approach ( va ) which treats the on - site interaction @xmath0 exactly and the intersite interactions @xmath1 within the mean - field approximation ( mfa ) .
we have also derived exact ground state results for a  linear chain ( ) .
moreover , at several results obtained within the random phase approximation ( rpa ) and the spin - wave approximation ( swa ) ( for and lattices ) and within the systematic low density expansions ( lde ) ( for lattices ) are presented .    within the va
the phase diagrams of model ( [ row : ham1 ] ) have been investigated in  @xcite as a function of the electron concentration @xmath2 , but the stability conditions of states with phase separation have not been discussed .    in this paper , we investigate the thermodynamic properties of model hamiltonian ( [ row : ham1 ] ) for arbitrary chemical potential and arbitrary electron concentration .
we focus on the states with phase separation , which have not been analyzed till now .
arbitrary values of on - site interaction @xmath0 are considered and the range of pair hopping interaction @xmath1 is restricted to the nearest neighbours .
our investigation of the general case finds that , depending on the values of the interaction parameters and the electron concentration , the system can exhibit the superconducting ordered and nonordered homogeneous phases as well as the phase separation between them .
transitions between various states and phases can be continuous and discontinuous , what implies existence of tricritical points on the phase diagrams .
we present detailed results concerning the evolution of phase diagrams and thermodynamic properties as a function of interaction parameters , chemical potential @xmath3 and electron concentration @xmath2 . by
taking into account the phase separated states and presenting several new results for lattices of various dimensions derived beyond va ( rigorous solution for a  chain ( ) , rpa for sq and sc lattices and lde in ) our paper substantially extends and generalizes the results and conclusions of previous papers concerning the model considered  @xcite .
the paper is organized as follows . in section
[ sec : method ] we describe the method used in this work . there are also derived explicit formulas for the free energies of homogeneous phases and states with phase separation as well as equations determining the superconducting order parameter and the chemical potential in homogeneous phases obtained in va . in section
[ sec : gs ] we analyze the properties of the system at zero temperature and present exact ground state results for a  chain ( section  [ sec : gsexact ] ) , diagrams obtained in the va ( section  [ sec : gsmfa ] ) as well as results obtained in some approximations going beyond va for lattices ( section  [ sec : rpa ] and section  [ sec : lde3d ] ) .
section  [ sec : kt ] is devoted to the study of the finite temperature phase diagrams derived within the va as a function of @xmath3 ( section  [ sec : ktmi ] ) and as a function of @xmath2 ( section  [ sec : ktn ] ) .
some particular temperature dependencies of the thermodynamic parameters are discussed in section  [ sec : properties ] .
finally , section  [ sec : conclusions ] is devoted to conclusions and supplementary discussion .
the appendices report some relevant details : explicit expressions of site - dependent self - consistent va equations ( [ app : sitedependentequations ] ) , particle - hole transformations ( [ app : sym ] ) , and selected analytical va results ( [ app : anal ] ) .
the free energy of the system as well as the self - consistent equations for the average number of electrons in the system and the average value of the order parameter are derived within site - dependent va in appendix [ app : sitedependentequations ] . restricting our analysis to the nearest neighbours and assuming no spatial variations of the order parameter the grand potential per site obtained in va is given by : @xmath25 where @xmath26 and @xmath27 .
the free energy per site : is derived as @xmath28 the grand potential @xmath29 is minimized in the system considered for fixed @xmath3 , whereas the free energy @xmath30 takes the minimum when @xmath2 is fixed .
the condition for electron concentration ( [ row : condn1 ] ) and a  minimization of @xmath29 ( or @xmath30 ) with respect to the superconducting order parameter @xmath31 lead to the following self - consistent equations ( for homogeneous phases ) : @xmath32.\ ] ] we also introduce the concentration of locally paired electrons and the ratio . notice that @xmath33 is different from the condensate density ( a  fraction of pairs in the condensate ) @xmath34 .
these quantities differ even in the va at : and in the superconducting phase .    equations ( [ row : mfa1])([row : mfa2 ] ) are solved numerically for and we obtain @xmath31 and @xmath2 when @xmath3 is fixed or @xmath31 and @xmath3 when @xmath2 is fixed .
the superconducting phase ( ss ) is characterized by non - zero value of @xmath31 , whereas in the non - ordered ( normal ) phase ( no ) .
equation  ( [ row : mfa2 ] ) is only the necessary condition for an  extremum of ( [ row : grandpotential ] ) ( or ( [ row : freeenergy ] ) ) thus the solutions of can correspond to a  minimum or a  maximum ( or a  point of inflection ) of @xmath29 ( or @xmath30 ) . in addition , the number of minima can be larger than one , so it is very important to find the solution which corresponds to the global minimum of ( [ row : grandpotential ] ) ( or ( [ row : freeenergy ] ) ) .
we say that the solution of ( [ row : mfa1])([row : mfa2 ] ) corresponds to a  metastable phase if it corresponds to a  minimum of the free energy @xmath30 or grand potential @xmath29 and the stability condition @xmath35 is fulfilled .
otherwise , we say that the phase is unstable .
a  stable ( homogeneous ) phase is a  metastable phase with the lowest free energy ( among all metastable phases and the phase separated states described below ) .
phase separation ( ps ) is a state in which two domains with different electron concentration : @xmath36 and @xmath37 exist in the system ( coexistence of two homogeneous phases ) .
the free energies of the ps states are calculated from the expression : @xmath38 where @xmath39 are values of a free energy of two separating phases at @xmath40 corresponding to the lowest homogeneous solution for a  given phase , @xmath41 is a fraction of the system with a charge density @xmath36 , is a  fraction with density @xmath37 ( ) and .
the minimization of ( [ row : freeenergyps ] ) with respect to @xmath36 and @xmath37 yields the equality between the chemical potentials in both domains : @xmath42 ( chemical equilibrium ) and the following equation ( so - called maxwell s construction ) : @xmath43 which is equivalent with equality of grand potentials per site in domains : .
it implies that the transitions with a  discontinuous change of @xmath2 in the system considered for fixed @xmath3 can lead to occurrence of the regions of phase separation in the concentration range on the diagrams obtained as a function of @xmath2 . in these regions the homogeneous phases can be metastable as well as unstable , depending on the @xmath2-dependence of @xmath3 ( cf .
( [ row : stabcondition ] ) ) . in the ps states the chemical potential is independent of the electron concentration , i.e. . for more detailed discussion of maxwell s construction
we refer the reader to existing literature , e.  g.  @xcite .    in the model
considered only one type of ps states can occur , which is a coexistence of ss and no phases .    in the paper
we have used the following convention .
a  second ( first ) order transition is a  transition between homogeneous phases with a  ( dis-)continuous change of the order parameter at the transition temperature . a  transition between a  homogeneous phase and the ps state is symbolically named as a  `` third order '' transition .
at this transition a  size of one domain in the ps state decreases continuously to zero at the  transition temperature .
second order transitions are denoted by solid lines on phase diagrams , whereas dotted and dashed curves denote first order and `` third order '' transitions , respectively .
the phase diagrams obtained are symmetric with respect to half - filling because of the particle - hole symmetry of hamiltonian ( [ row : ham1 ] ) , so the diagrams will be presented only in the range and .
the transformation formulas of various thermodynamic parameters and quantities such as electron concentration , superconducting order parameter , chemical potential and free energy under the electron - hole transformation are collected in appendix [ app : sym ] .
the numerical computations have been done in c / c++ with source codes from  @xcite .
the minima have been found by using brent s method in one dimension and downhill simplex method in multidimensions .
we have also implemented newton - raphson method using derivative to solve non - linear equations .
let us consider the ground state ( gs ) of model ( [ row : ham1 ] ) .    in the absence of @xmath14 term , the parity of the number of particles at each site is a  good quantum number .
thus , one can decompose the eigenspace of @xmath44 into sectors specified by the parity of the occupation number at each site ( @xmath45 ) . within any sector ,
the system is a collection of _ even _ segments , within which pairs are confined ( the subspace with excluded singly occupied sites ) interleaved by _ odd
_ segments ( the subspace with excluded doubly occupied sites ) .    since , the interfaces between the _ odd _ and _ even _ regions will cost energy , so in the thermodynamic limit the ground state will be either a  single infinitely large _ even _ segment or a  single _ odd _ segment for fixed @xmath3 . only in the case when the grand potentials per site ( it is the chemical potential ) of the segments are equal there
is possibility of coexistence of both types of segments if
it appears as a  first order line on the gs phase diagram .
the situation is different when @xmath2 is fixed . in such case
the ps states can occur in definite ranges of @xmath2 .
the reason for this can be explained as the following .
let the energy of a single boundary ( interface ) between _ odd _ and _ even _ segments be .
if the number of intersite boundaries will increase as @xmath46 ( where ) the contribution from them to the total free energy disappears in the thermodynamic limit .    within the subspace excluding the double occupancy of sites ( by electrons for or holes for ) the product for any , i.  e. the @xmath1 term does not contribute .
effective hamiltonian ( [ row : ham1 ] ) for this case takes the form and the calculations of the gs grand potential and the gs energy are quite obvious . as a result in the limit from ( [ row : freeenergy ] ) and ( [ row : mfa1 ] ) one obtains for fixed @xmath2 in the no phase : @xmath47 this gs is infinitely degenerate for .
for ( ) , @xmath48 sites are occupied by electrons ( holes ) and the spins of all particles are independent of each other .
inclusion of the intersite spin - spin and density - density interactions can remove this degeneracy and generate various magnetic and charge orderings .
the magnetic moments are not suppressed in this state .
for there is no order pattern in the distribution of the particles , whereas for mott phase with one - electron - per - site configuration ( with random direction of electron spin per site ) occurs .    from ( [ row : grandpotential ] ) and ( [ row : mfa1 ] ) for fixed @xmath19 and one obtains : @xmath49 the average occupancy of sites can only be integer for fixed @xmath19 , what can be intuitively understand if one realizes that only and every site is equivalent ( there is no interactions between different sites in @xmath50 ) .
the state @xmath51 ( no  empty ) is stable with respect to @xmath52 ( no  mott ) if ( and ) , whereas if ( and ) .    within the subspace excluding the single occupancy of sites the eigenvalues of charge operators
are , for each site .
thus , from ( [ row : ham2 ] ) , the effective hamiltonian for this case takes the form : @xmath53 with auxiliary condition ( [ row : condn2 ] ) .
@xmath54 has the form of the xy quantum spin model ( ) with an effective external field in the @xmath11-direction , such that the average magnetization has a fixed value equal to ( cf .
( [ row : condn2 ] ) ) . the gs energy of hamiltonian ( [ row : ham3xy ] )
is derived as @xmath55 where .
the gs grand potential in this case has  the form @xmath56 where it is underlined that @xmath19 is an independent variable .
as it follows from previous analyzes  @xcite the gs of model ( [ row : ham3xy ] ) exhibits the long - range xy order for dimensional lattice in the whole concentration range .
for , only the short range xy correlations are developed @xcite , nevertheless even arbitrary small interchain interactions can induce the real long - range order also in this case .
the ss phase will be realized if for fixed @xmath2 ( if for fixed @xmath19 ) . for explicit calculation of the phase boundary between the no ( eqs .
( [ row : gsexe1 ] ) and ( [ row : gsexomega1 ] ) ) and the ss ( eqs .
( [ row : gsexe2 ] ) and ( [ row : gsexomega2 ] ) ) phases and determination of the gs properties of the ss phase for various dimensions and the electron concentration or the chemical potential one has to find the solution of model ( [ row : ham3xy ] ) at .
let us stress that if @xmath2 is fixed one has to consider also the ps state in which one domain is in the no phase and the other is in the ss phase and which energy is calculated from ( [ row : freeenergyps ] ) .
however , the gs energy of the ps ( ) state can be obtained in the most cases only numerically by minimization of the free energy given by ( [ row : freeenergyps ] ) . in the no domain ( mott phase with one - electron - per - site configuration ) , whereas concentration in the ss domain is dependent on @xmath57 ( and approach used ) .
+      the exact solutions of model ( [ row : ham3xy ] ) for arbitrary @xmath2 can be derived for and . for , by making use
the exact results for the ground state energy of xy model in a  transverse field @xcite one obtains : @xmath58 where ( because ) .
notice , that due to a  well known isomorphism between the planar ferromagnet and antiferromagnet for loose packed ( alternating ) lattices with nearest - neighbour @xcite , the results for @xmath59 are the same for both signs of @xmath1 .
thus , from ( [ row : gsexe2 ] ) and ( [ row : energyxy1d ] ) we have @xmath60 and the exact expressions for chemical potential @xmath3 and the superfluid density @xmath61 ( the short - range order parameter ) are : @xmath62 where @xmath63 denotes the nearest neighbour of the @xmath7-site . for ( ss case ) ferromagnetic correlations between xy components of neighbouring pseudospins are present , whereas for ( @xmath23s case ) antiferromagnetic correlations occur . from ( [ row : missexgs ] ) one obtains that the ss can exist in the range .
notice that stability condition ( [ row : stabcondition ] ) is always fulfilled in the ss phase .
the gs phase diagrams obtained are shown in figure  [ rys : gsdiagrams ] : as a function of @xmath19 ( panel ( a ) ) , and as a function of @xmath2 without ( panel ( b ) ) and with ( panel ( c ) ) consideration of the ps state . at half - filling
there is no phase separation and the no - mott phase is stable for .
the boundaries between different states corresponding to the exact results in are denoted by dotted lines ( if they are different from those obtained in va ) .
notice that for homogeneous phases in the ranges of the ps state occurrence .      for there are no exact solutions for model ( [ row : ham3xy ] ) available ( except of the lde for  see @xcite ) and the approximate approaches are necessary to find @xmath64 .
the simplest of these approaches is the va which at gives ( the limit of ( [ row : freeenergy])([row : mfa2 ] ) for ) @xmath65 equation ( [ row : missvags ] ) implies that the ss phase can exist in the range .
notice that , which means that all electrons in the system are locally paired in the ss phase at , whereas the finite value of the depletion ( ) at is due to the hard - core effect  @xcite .
the gs energy of the ps state can be obtained by minimization of ( [ row : freeenergyps ] ) . in the no domain ( mott phase ) , whereas concentration in the ss domain is dependent on @xmath57 and ( ) .
it corresponds to the first order boundary on the @xmath57  vs.  @xmath66 gs diagram determined by equation .
the gs phase diagrams obtained in va are shown in figure  [ rys : gsdiagrams ] .
the va phase boundaries are denoted by solid lines .
similarly as in section  [ sec : gsexact ] one finds that for homogeneous phases in the ranges of the ps state occurrence .    for hypercubic lattice of the dimension the va becomes the exact theory for model ( [ row : ham1 ] ) ( as well as for model ( [ row : ham3xy ] ) ) .
going beyond va for a  reliable approach is the self - consistent rpa  @xcite .
this approach has been proven to be a  very good approximation scheme in problems of quantum magnetism and it fully takes into account quantum fluctuations , which can be of crucial importance for the considered system for .
the gs energy of model ( [ row : ham3xy ] ) is given by ( [ row : gsexe2 ] ) , where .
@xmath67 and @xmath68 denote the space - fourier transform of the charge operators and the pair hopping integral , respectively , and the chemical potential is determined by    the energy @xmath69 is obtained by means of green s functions using the spectral theorem , and by applying the approximation for the longitudinal correlation functions @xcite . at this yields
@xmath70   \\ & - & ( 1+\cos^2\theta)\frac{r}{2n}\sum_{\vec{k } } i_{\vec{k } } \frac{a_{\vec{k}}}{e_{\vec{k } } } , \nonumber\end{aligned}\ ] ] where , , is the collective excitation spectrum , , @xmath71 is given by and the length of pseudospin @xmath72 is given as a  solution of the self - consistent equation @xmath73.\end{aligned}\ ] ] the result corresponding to the swa ( being the lowest order expansion of above equation ) is with @xmath74 determined by ( [ row : rpapsi ] ) .
for the ss phase the spectrum @xmath75 is gapless and linear in @xmath76 for small @xmath77 .
one gets with the sound velocity @xmath78 given by ( ) .    the chemical potential in the ss phase in the rpa is given by , which is the same as in va .
the quantum corrections to this va result for @xmath79 can be evaluated as in @xcite and one obtains : @xmath80,\ ] ] where , .
one should notice that rpa and swa give very similar result for the model ( [ row : ham3xy ] ) in the gs  @xcite .
the resulting phase boundaries for square ( sq ) lattice and simple cubic ( sc ) lattice are shown in figure  [ rys : gsdiagrams ] and denoted by dashed and dashed - dotted lines , respectively .      in @xcite several rigorous results concerning the gs properties of the model of hard - core charged bosons
have been derived for various lattices in using a  systematic lde based on the knowledge of the exact two - body scattering amplitude .
transforming the results of section  vi  b from @xcite into our problem one gets for the following exact expressions .
the gs energy of ss phase @xmath81 is given by @xmath82 , \nonumber\end{aligned}\ ] ] where is the exact scattering length for the considered case , is the watson integral for a  given lattice , , @xmath83 is the volume of the unit cell , @xmath84 is the lattice constant . in particular , for sc : , , ; for bcc : , , ; and for fcc : , , .    the chemical potential has a form @xmath85\ ] ] and the sound velocity @xmath78 is derived by @xmath86.\ ] ]    the phase diagrams obtained within lde are in very good agreement with rpa results ( for sc lattice ) in the range . the homogeneous ss phase boundary derived within lde
is moved little downwards in comparison with the rpa ones .
an analysis of lde results obtained for various cubic lattices in shows that the region of the ss phase stability shrinks with increasing number of the nearest neighbours in the lattice .     for increasing values of @xmath57 ( numbers above the curves ) .
solid and dotted lines indicate second order and first order transitions , respectively .
the dashed - dotted curve ( labeled as tcp ) is a  projection of the tricritical point line on the @xmath87  vs.  @xmath66 plane.,scaledwidth=49.0% ]     and @xmath57 .
solid and dotted lines indicate second order and first order transitions , respectively .
the tricritical point line is denoted by dashed - dotted line.,scaledwidth=49.0% ]      the general structure of the gs phase diagrams derived is similar for all lattice dimensionalities ( ) and approaches used .
there are two sequences of phase transitions with increasing @xmath57 for fixed @xmath2 :    * ss @xmath88 ps @xmath88 no ( for ) , * ss @xmath88 no ( mott ) ( for ) .    the effects of quantum fluctuations increase with decreasing lattice dimensionality and for a  linear chain their influence on phase diagram is the most prominent .
thus the region of the homogenous ss phase stability decreases with increasing @xmath89 . in the quantum fluctuations
are suppressed .
one should also notice that the transition ps - no at is independent of @xmath2 and dimensionality of the lattice .
the ps region expands with increasing dimensionality .
+         all phase transition boundaries , necessary to construct the complete phase diagram for fixed @xmath3 , have been obtained numerically by solving the set of equations ( [ row : mfa1])([row : mfa2 ] ) and comparing the grand potential ( [ row : grandpotential ] ) for the solutions found .
the resulting @xmath87  vs.  @xmath66 phase diagrams obtained for various fixed values of the on - site interaction @xmath57 are shown in figure  [ rys : psktvsmi ] .
the overall behaviour of the system as a function of @xmath19 is demonstrated in figure  [ rys : psktvsmi3d ]
.    for only the second order transitions occur with increasing temperature .
the transition temperature is maximal for , and it decreases monotonically with increasing @xmath57 and @xmath90 .
the most interesting is the range . in this range
the of the transition temperature becomes non - monotonic and with increasing @xmath57 the maximum of the transition temperature moves towards higher values of @xmath90 . in this range of the on - site interaction ,
the order of the transition changes from second order to first order with increasing of @xmath66 . a  _ tricritical point _ ( tcp ) is connected with the change of transition order .    for
only the no phase is stable at any .    as it will be showed in section  [ sec : ktn ] the first order phase boundaries are associated with the existence of the ps states in definite ranges of @xmath2 if the system is considered for fixed @xmath2 .
the transition temperatures have been determined by self - consistent solving of and and comparing free energies ( [ row : freeenergy ] ) and ( [ row : freeenergyps ] ) .
the concentration dependencies of ( i )  the second order transition temperatures and ( ii )  the `` third order '' transition temperatures can be also easily obtained from the diagrams as a function of @xmath19 , because of a  simple relation between the electron concentration and the chemical potential in ss phase . the transitions ( i ) and ( ii ) for fixed @xmath2 correspond to the second order and first order transitions for fixed @xmath3 , respectively
. however , there is no simple correspondence for the transitions ( cf .
( [ row : mino ] ) ) .
examples of the @xmath87  vs.  @xmath2 phase diagrams evaluated for various values of the on - site interaction @xmath57 are shown in figure  [ rys : psktvsn ] .
for the ps states do not occur and the obtained phase diagrams are the same as those derived in @xcite . the transition between homogeneous ss and no phases taking place with increasing temperature is second order for arbitrary @xmath2 . the transition temperature is maximal for , and it decreases monotonically with increasing @xmath57 and decreasing @xmath2 ( for )  cf .
figure  [ rys : psktvsn]a .    in the range with increasing @xmath57 the maximum of the transition temperature moves towards lower concentrations ( for ) . in a  definite range of @xmath57 and
@xmath2 the ps state is stable . for the ps state extends from the ground state ( cf .
figures  [ rys : psktvsn]d , e , f and figure  [ rys : gsdiagrams]c ) , whereas for it is stable only at finite temperatures ( cf .
figures  [ rys : psktvsn]b , c ) . the critical point for the phase separation ( denoted as @xmath91 , which is a  _ tricritical point _ , tcp ) lies on the end of the second order line .    between and only the no phase is stable at any and there are no transitions .                    for given value of @xmath2 the @xmath87  vs.  @xmath57 phase diagrams of the system
are shown in figure  [ rys : psktvsu ] .
all of the transition temperatures decrease with increasing @xmath57 . in the range the diagrams are rather similar and they consist of three regions (
exemplary diagrams are shown in figures  [ rys : psktvsu]a , b ) .
the regions of ss and no phases stability are separated by second order boundary or , for larger on - site repulsion , by the region of ps state occurrence . for half - filling ( , )
only homogeneous phases are present on the diagram ( figure  [ rys : psktvsu]c ) .
the @xmath91-point , connected with the change of the transition order , is located at and ( for ) .
the possible sequences of transitions with increasing temperatures and the transition orders of them are listed below :    * ss@xmath88no : second order for and second order or first order for ( and ) , * ps@xmath88no :  third order  , it can take place only for ( and ) , * ss@xmath88ps@xmath88no : both  third order  , it can take place only for ( and ) .
notice that increasing temperature usually suppresses heterogeneity in a  system what favours occurrence of homogeneous phases rather than ps states . in the model considered the ps state can develop at higher temperatures than the homogeneous ss phase .
this rather unusual behaviour can be connected with the interplay between the pair hopping interaction and the on - site repulsion .
in this section we present several representative dependencies of the thermodynamic characteristics for fixed model parameters ( figures  [ rys : homoprop][rys : seppropfunkn ] ) .    in particular for fixed @xmath19 ( figure  [ rys : homoprop ] ) , one can single out two limiting types of thermodynamic behaviour near transition temperature @xmath92 : ( i )  the _ local pair regime _ and ( ii )  the _ pair breaking regime_. in between , there is a  crossover between the two regimes ( see e.  g. the plots @xmath93 vs. @xmath87 in figure  [ rys : homoprop]b ) .    for large on - site attraction
the concentration of locally paired electrons exhibits no sharp feature as the temperature is lowered through @xmath92 .
the number of non - paired electrons at @xmath92 is negligible and the transition is to the state of dynamically disordered pairs ( only for and , the _ local pair regime _ ) .    in the second limit , for on - site repulsion
, @xmath94 has a sharp break at @xmath92 and a  substantial fraction of single particles can exists above @xmath92 .
we call this the _ pair breaking regime_. as temperature is lowered , the condensate growths both from a  condensation of pre - existing pairs and from binding and condensation of single particles . for small binding energies , if ( @xmath95 ) , there will be essentially no pre - formed pairs at @xmath92 .
notice that non - zero value of @xmath94 does not imply that local pairs are in coherent state and even significant values of @xmath93 are possible in the no phase . in the limit
@xmath94 increases to ( each of four states at a  given site can be occupied with equal probability ) . the condensate density ( which can be approximated as at least for , ) vanishes for , but
the doubly occupied states are still thermally excited above @xmath92 ( ) .    in figure
[ rys : homoprop]a the temperature dependencies of the superconducting order parameter @xmath96 are presented , where one can see clearly the discontinuous change of the order parameter ( lines c and g ) .
the other lines correspond to second order transitions .
finally , let us briefly summarize the behaviour of the specific heat at constant volume ( figure  [ rys : homoprop]c ) .
the no phase is characterized by the relatively broad maximum in @xmath97 connected with continuous changes in a  short - range electronic ordering ( in higher temperatures , not shown in figure  [ rys : homoprop ] ) .
the narrow peak in @xmath98 is associated with the first order transition , while the @xmath99-point behaviour is typical for the second order transition .
let us concentrate now on the thermodynamic properties of the system in the ps state and their changes at the `` third order '' transitions for fixed electron concentration @xmath2 ( figures  [ rys : seppropfunkkt ] and  [ rys : seppropfunkn ] ) .
as an example , we consider the temperature dependencies of : ( i )  the specific heat , ( ii )  the entropy , ( iii )  the concentration of paired electrons @xmath93 , ( iv )  the chemical potential @xmath66 , and ( v )  the fractions and ( sizes of domains ) with a  charge density @xmath36 and @xmath37 ( figure  [ rys : seppropfunkkt ] ) as well as the concentration dependencies of : ( a )  @xmath66 , ( b )  @xmath100 and ( c )  @xmath93 ( figure  [ rys : seppropfunkn ] ) .
thermodynamic characteristics for three different possible sequences of transitions with increasing temperature are shown in figure  [ rys : seppropfunkkt ] . the label ( a ) corresponds to the case of transition .
at such a  transition the specific heat @xmath97 exhibits a  finite jump at transition temperature , what is associated with the continuity of the entropy @xmath22 ( cf .
figure  [ rys : seppropfunkkt]a ) .
the lines labelled by ( b ) correspond to case . at the transition
there is also a  finite jump in @xmath97 , however the @xmath97 in the ps state is larger than that in the ss phase , contrary to the transition , where the jump has an opposite sign . finally , the lines ( c ) correspond to the case of second order ss@xmath88no transition .
the concentration ratio of locally paired electrons @xmath94 is presented in figure  [ rys : seppropfunkkt]b . on the inset of figure  [ rys : seppropfunkkt]b the chemical potential is shown as a function of temperature .
one can see that @xmath19 in the ss phase is independent of @xmath101 ( cf .
( [ row : miss ] ) ) . in figure  [ rys : seppropfunkkt]c there are shown dependencies of the fractions of the system @xmath102 and @xmath103 with charge density @xmath36 and @xmath37 , respectively .
for instance , in the case ( a ) the ss domain with lower concentration ( @xmath37 ) vanishes at `` third order '' transition .    the chemical potential @xmath66 and the free energy @xmath100 as a function of concentration for and are presented in figure  [ rys : seppropfunkn]a and figure  [ rys : seppropfunkn]b , respectively . for the ss phase is stable and its free energy is the lowest one . in the range
the ps state has the lowest energy ( solid line ) .
the solutions of corresponding to the lowest energy ones ( the homogeneous phases with the lowest energy ) are thermodynamically stable i.  e. @xmath104 in their existence ranges ( denoted by dashed lines in figures  [ rys : seppropfunkn]a ,  b ) . at point labelled as @xmath105 the transition between metastable phases occurs , eg . for the free energy @xmath106 of the ss phase
is lower than the energy @xmath107 of the no phase , whereas for .
the dotted lines in figures  [ rys : seppropfunkn]a ,  b correspond to the homogeneous phases with the highest energy .
one should notice that in the described case for and in the whole range of the ps state occurrence two homogeneous phases are metastable , whereas in the ranges of the homogeneous phases stability there is one metastable phase at the most .
the linear dependence is characteristic for the ss phase , whereas in the ps state @xmath66 is independent of @xmath2 .
this is consistent with square dependence of @xmath106 and linear dependence of @xmath108 as a function of @xmath2 .
the anomalous decrease of the chemical potential with doping in the lowest energy metastable homogeneous phases near @xmath105-point ( i.  e. an abrupt drop at @xmath105-point ) is a  signal that it is thermodynamically more convenient to phase separate into two subsystems with the different electron densities @xmath37 and @xmath36 ( figure  [ rys : seppropfunkn]a ) .
the density range for which the ps occurs can be readily found by the maxwell s construction described in section  [ sec : method ] .
the @xmath2-dependencies of the chemical potential @xmath66 , the free energy @xmath100 and the concentration of paired electrons are shown on the insets of figure  [ rys : seppropfunkn]a and figure  [ rys : seppropfunkn]b and in figure  [ rys : seppropfunkn]c , respectively , for at several fixed temperatures .
they are obtained for : ( a ) , ( b ) , ( c ) , and ( d ) .    at `` third order '' transitions ( ss - ps and ps - no ) the thermodynamic parameters such as @xmath19 , @xmath93 , @xmath22 , @xmath30 are continuous , whereas the @xmath97 exhibits a  finite jump as at the second order transitions .
however , one should notice that the order parameter for `` third order '' transitions is the concentration difference ( not @xmath96 , which is the order parameter in one domain ) and its change is discontinuous at transition temperature .
such transitions are present if the system is consider for fixed @xmath2 and they are associated with first order transitions at fixed @xmath19 .
we have studied a  simple model of a  superconductor with very short coherence length ( i.  e.  with the pair size being of the order of the radius of an effective lattice site ) and considered the situation where the single particle mobility is much smaller than the pair mobility and can be neglected .
we have evaluated the phase diagrams and thermodynamic characteristics of the model at fixed @xmath3 and at fixed @xmath2 and determined the ranges of the ps occurrence in the system .
one should stress that within the va the on - site @xmath0 term is treated exactly .
thus , the major conclusions of our paper concerning the evolution of the properties of the system with @xmath0 are reliable for arbitrary @xmath0 .
the mean - field used for the intersite term is best justified if the @xmath109 interactions are long - ranged or if the number of nearest neighbours is relatively large .
the derived va results are exact in the limit of infinite dimensions , where the mfa treatment of the intersite interaction @xmath1 term becomes the rigorous one . moreover , the va yields exact results ( in the thermodynamic limit ) for @xmath109 of infinite range ( for any @xmath110 ) regardless of the dimensionality of the system  @xcite .
the properties of the system are strongly dependent on the ratio of @xmath111 , which is related to the relative values of the pair binding energy @xmath112 and the effective pair mobility @xmath113 . for the case of a  single pair in a  system , i.  e. for , we have
, ( the pair effective mass is ) @xcite and increasing @xmath2 can modify @xmath112 .
there are two well defined limits of the model : ( i )  the _ local pair regime _ and ( ii )  the _ pair breaking limit_. in between there is a  crossover between the two regimes for fixed @xmath3 .
the _ local pair limit _ is realized for large pair binding energy ( i.  e. for ) .
in such a  case , the two interactions : the on - site interaction @xmath0 and the intersite pair hopping @xmath1 cooperate and transition temperature is determined by pair mobility ( i.  e. by @xmath1 ) . the transition with increasing temperature is second order and is to the no phase being a  state of dynamically disordered local pairs .    the opposite regime i.  e. the _ pair breaking limit _ is realized for substantial values of .
repulsive @xmath0 destabilizes on - site pairing and competes with @xmath1 .
as we have shown above , the increasing @xmath0 : ( i )  for fixed @xmath3 it changes first the nature of phase transitions from a  continuous to a discontinuous type and then it suppresses superconductivity for , ( ii )  it stabilizes the phase separation state ( ) in definite ranges of @xmath2 and temperature .
this behaviour is associated with the presence of the tricritical point on the phase diagrams . with increasing @xmath57 the maximum of the transition temperature moves towards lower concentrations ( for ) .
in definite range of @xmath57 and @xmath2 the ps state is stable . for the ps state extends from the ground state ,
whereas for it is stable only at finite temperatures .
finally , for large @xmath0 ( ) , the system remains in a  normal state at all temperatures for any chemical potential or concentration .
comparing the va results at with the other ones ( exact ones for chain , rpa ones for square lattice and cubic lattices ) , one can see that the general structure of the diagrams is similar .
it implies that the va can give qualitatively reasonable results also for lattices of finite dimensionality .
however , one should stress that in finite dimensions due to quantum fluctuations connected with @xmath109 interactions the regions of the homogeneous ss phase occurrence are extended in comparison with the va results .
moreover , in a  linear chain only the short - range order occurs at in contrary to the va and rpa ( swa ) in , where the long - range order is present .
one should emphasize that for the va overestimates transition temperatures and gives incorrect concentration dependence of them ( for , cf .
( [ row:2ndordertssno ] ) ) . for ideal bose gas with @xmath114 spectrum ( the and limit of the model considered ) the @xmath2-dependence of the temperature is , what can be obtained , for example in the rpa , even for finite values of @xmath57  @xcite .
the rpa predicts also the existence of tricritical point in finite temperatures  @xcite in , what confirms that the va can give qualitatively reasonable results also for real lattices in finite temperatures .
the rpa results for at @xcite are in agreement with the mermin - wagner theorem @xcite , which proofs the absence of the long - range order for at . in one
has the kosterlitz - thouless - type transition with increasing temperature .
the model considered can describe the systems in which the single electron mobility ( single electron hopping ) is much smaller than the pair mobility ( kinetic energy of pairs ) and it can be neglected in first approximation .
let us comment on the finite electron bandwidth effects .
the presence of the hopping term breaks a  symmetry between ( favouring ss ) and ( favouring @xmath23s ) cases .
the phase diagrams of the penson  kolb  hubbard model are more complicated than those obtained for hamiltonian ( [ row : ham1 ] ) @xcite , even in the ground state .
we can suppose that small but finite single electron hopping @xmath14 will not qualitatively alter the phase diagrams , at least for the case .
the main effect of @xmath14 ( for , ) is a  renormalization of the pair hopping term and an introduction of an effective intersite density - density repulsion . for and the charge density wave state can also occur @xcite . for and
both signs of @xmath1 the @xmath14 term generates magnetic correlations competing with superconducting ones and its effects essentially modify the phase diagrams and the properties of normal state ( particularly for larger values of @xmath14 ) . in such a  case
it is necessary to consider also various magnetic orderings .
moreover , several phase separation states involving superconducting , charge and ( or ) magnetic orderings could also be stable for .
for larger values of @xmath14 the studies of gs phase diagrams of the half - filled pkh model for @xmath89-dimensional lattices have been determined by means of the ( broken symmetry ) hartree - fock approximation ( hfa ) and by the slave - boson mean - field method ( sbmfa ) in @xcite . for the diagrams supported also by the results of the continuum - limit field theory and exact diagonalisation studies
@xcite are shown to consist of at least nine different phases including superconducting states , bond located antiferromagnetic and charge - density - wave states as well as mixed phases with coexisting site and bond orderings .
the stability range of the bond - type orderings shrinks with increasing lattice dimensionality and for the phase diagrams involve exclusively site - located orderings .
one should notice that the phase boundary at and determined by the hfa and sbmfa ( ) are in agreement with our va result , even for relatively large @xmath14 ( especially for ) @xcite . obviously , for the transition is to the antiferromagnetic mott phase .     between nearest neighbours for half - filling obtained within va ( solid lines ) .
the dashed lines denote the rigorous boundaries for 1d - chain .
all transitions between homogeneous phases are discontinuous . .,scaledwidth=49.0% ]    the intersite density - density interactions @xmath115 can introduce mixed phases and ( or ) phase separations between charge orderings and superconductivity @xcite .
the ground state phase diagram in presence of interaction @xmath116 between nearest neighbours ( added a  term to hamiltonian ( [ row : ham1 ] ) ) calculated for half - filling is presented in figure  [ rys : uwi ] ( cf . with @xcite , where the case has been analyzed ) . for sufficiently strong intersite repulsion
the homogeneous charge ordered ( co ) phase occurs ( with ) , whereas intersite attraction introduce the ps state in which two no phases coexists ( with and , i.  e. electron droplets , ps[no / no ] ) . for @xmath117 the regions of the ss ( or @xmath23s , depending on the sign of @xmath1 ) and co phases occurrence are extended , whereas a  region of the no ( mott ) phase is reduced .
the superconducting order ( in @xmath117 ) , similarly as in section  [ sec : gsexact ] , is a  short - range order whereas the co is a  long - range order ( cf .
also @xcite ) .
the electron concentration @xmath2 and chemical potential @xmath3 are ( thermodynamically ) conjugated variables in the bulk systems .
however , one can fit the concentration rather than the potential in a  controlled way experimentally .
in such a  case @xmath3 is a  dependent internal parameter , which is determined by the temperature , the value of @xmath2 , and other model parameters ( cf .
( [ row : condn1 ] ) ) .
thus the obtained phase diagrams as a  function of the concentration are quite important because in real systems @xmath2 can vary in a  large range .
although our model is ( in several aspects ) oversimplified , it can be useful in qualitative analysis of experimental data for real narrow - band materials and it can be used for better understanding of various systems mentioned below and in section  [ sec : intro ] . in particular
, our results predict the existence of the phase separation ( ) near the mott state and describe its possible evolution and phase transitions with increasing temperature and a  change of @xmath2 ( @xmath19 ) .
of course the ps instability is specific to the short - range nature of the interactions in the model . when the ( unscreened ) long - range coulomb interactions are included the large - scale ps of charged particles is prevented and only a  frustrated ps can occur ( mesoscale , nanoscale ) with the formation of various possible textures  @xcite .
the electron phase separation involving ss is shown experimentally in several systems .
in particular , for special cases of la@xmath118cuo@xmath119 and la@xmath120sr@xmath121cuo@xmath119 , muon and superconducting quantum interference measurements suggest that the electron inhomogeneities move beyond local variations to form fully phase separated regions @xcite .
organic compounds also exhibit the superconductor - insulator phase separations as a  result of the external pressure ( e.  g. quasi - one dimensional ( tmtsf)@xmath118pf@xmath122 @xcite , ( tmtsf)@xmath118reo@xmath123 @xcite ) and fast cooling rate through the glass - like structure transition ( e.  g. @xmath124-(et)@xmath118cu[n(co)@xmath118]br @xcite ) .
finally we mention the family of iron - pnictides , e.  g. in ba@xmath126k@xmath121fe@xmath127se@xmath118 , where mesoscopic ps between ss and insulating ( magnetic ) phases have been observed @xcite .
the detailed analysis of rpa solutions will be explored further in subsequent publication  @xcite .
it is of interest to analyze the effects of external magnetic field as well as the impact of density - density @xcite and magnetic @xcite interactions on the phase diagrams of model ( [ row : ham1 ] ) .
the authors wish to thank t.  kostyrko and g.  pawowski for helpful discussions and a careful reading of the manuscript .
the work ( k.  k. ) has been financed by national science center ( ncn , poland ) as a research project in years 2011 - 2013 , grant no .
dec-2011/01/n / st3/00413 .
k.  k. would like to thank the european commission and ministry of science and higher education ( poland ) for the partial financial support from european social fund  operational programme `` human capital ''
pokl.04.01.01 - 00 - 133/09 - 00  `` _ _ proinnowacyjne ksztacenie , kompetentna kadra , absolwenci przyszoci _ _ '' .
within the va the on - site interaction term is treated exactly and the intersite interactions are decoupled within the mfa ( site - dependent ) : @xmath128 a  variational hamiltonian has the following form : @xmath129,\end{aligned}\ ] ] where , and .
@xmath130 can be diagonalized easily and a  general expression for the grand potential @xmath131 in the grand canonical ensemble in the va is @xmath132\right\},\ ] ] where .
the average value of operator @xmath133 is defined as @xmath134}{\textrm{tr}\left[\exp(-\beta\hat{h}_{0})\right]}.\ ] ] @xmath135 means a  trace of operator @xmath136 calculated in the fock space .
the explicit formula for the grand potential obtained in the va has the following form @xmath137 where @xmath138 and .
the expression for the average number of electrons at @xmath7-site is given by @xmath139 and the average of the charge exchange operator at @xmath7-site is derived as @xmath140 so one has a  set of self - consistent equations to solve consisting of : ( i ) @xmath13 equations in form of ( [ row : nonsite ] ) and ( ii ) @xmath13 equations in form of ( [ row : rhoonsite ] ) ( two equations for each site from @xmath13 sites ) .    the free energy is obtained as , where , is the number of electrons in the system ( @xmath2 is defined by ( [ row : condn1 ] ) ) .
the double occupancy @xmath141 of the site @xmath7 is determined by the following equation : @xmath142=\frac{\partial \omega_i}{\partial u},\ ] ] where @xmath45 is defined by ( [ row : nonsite ] ) .
the solutions of the set can correspond to a  minimum , a  maximum or a  point of inflection of the grand potential ( [ row : sitedependentgrandpotential ] ) .
in other words , the solutions of the set fulfil the necessary condition for a  minimum of @xmath131 with respect to . to find the solutions corresponding to stable ( or metastable ) states of the system , one should find a  minimum of with respect to all @xmath143 .
one can show easily that in the case of neglected spatial variations of the superconducting order parameter ( i.  e.   for every site @xmath7 ) reduce to eqs .
( [ row : grandpotential ] ) and obtained in section  [ sec : method ] .
let us introduce the particle - hole transformation on the alternate lattice defined as @xmath144 where if and if ( @xmath145 , @xmath146 label two sublattices in the alternating lattice ) .
then the operators appearing in hamiltonian ( [ row : ham1 ] ) are transformed to @xmath147 whereas the thermodynamic parameters and averages : chemical potential , free energy , specific heat , superconducting - order parameter and double occupancy transform as follows @xmath148 respectively . from any characteristics of the considered model in the range we can find the behaviour of the system in the range . in particular
, the phase diagrams are symmetric with respect to half - filling .
from the set one can obtain the relation between the chemical potential and the electron concentration in both homogeneous phases : @xmath149      let us stress that @xmath92 determined from the above equation can be not only a  transition temperature between stable phases but also a  transition temperature between metastable phases ( e.  g. in the regions of ps states stability ) or unstable phases .
however , all second order transition temperatures which were obtained in section  [ sec : kt ] by comparison of free energies fulfill ( [ row:2ndordertssno ] ) . from ( [ row:2ndordertssno ] )
one can obtain @xmath92 as a function of @xmath2 ( or @xmath19  cf .
( [ row : miss ] ) ) . in the limiting cases it yields @xmath151


in the last 20 years , the processes involved in the formation of solar type stars have been well characterized .
an evolutionary sequence has been developed in which a prestellar core ( i.e. * ? ? ? * ) begins to collapse under its own gravity and forms a class 0 object ( i.e. * ? ? ?
* ) with a large disk of material surrounding it .
this accreting protostar begins to produce an observable molecular outflow which releases the build up of angular momentum ( i.e. * ? ? ?
* ; * ? ? ?
* ) and evolves into the class 1 stage . later , as the protostar becomes more revealed , accretion slows and a pre - main - sequence star forms .
it then undergoes quasi - static contraction until core hydrogen burning commences and the star reaches the main sequence ( i.e. * ? ? ?
these processes have been well studied and characterized based in part on the high resolution observations possible towards the closest examples of low mass star forming regions .    the initial mass function ( imf ) of stars ( i.e. * ? ? ?
* ) shows that most of the mass contained in stars is within the lower mass stars , and that massive stars do not form nearly as often as their lower mass counterparts .
there are a number of nearby low mass star forming regions ( i.e. taurus , ophiucus , chameleon , and serpens are well studied nearby examples ) in which there is no massive star formation .
high mass stars do not , in general , form in isolation , but together with low mass stars ( i.e. orion , see * ? ? ?
the high densities required for massive star formation also mean that more of them form in the molecular ring .
this means that statistically , the closest examples of regions forming high mass stars are further away , and are also forming clusters of stars . both distance and clustering can cause confusion when attempting to determine the dynamics in a region
.    add to the large average distances and source confusion the higher energies involved in the formation of massive stars and the presence of hii regions surrounding the more evolved massive protostars , and it becomes clear why high mass star formation is not nearly as well understood as lower mass star formation . if we assume that massive stars form in a manner similar to their lower mass counterparts ( i.e. * ? ? ?
* ) , it is possible to suggest an evolutionary sequence starting from an infrared dark cloud ( i.e. * ? ? ?
* ) with a deeply embedded protostar which then begins to heat its surroundings and creates a hot core ( i.e. * ? ? ?
the protostar then begins burning hydrogen brightly enough to ionize its surroundings , forming a hypercompact hii region @xcite , which continues to expand into an ultracompact hii region @xcite as the protostar evolves towards and up the main sequence @xcite .
there are many observable phenomena , such as large scale infall and outflows , which appear to be common between both low and high mass star formation , which suggests that perhaps high mass stars may form via the accretion of matter onto a single object in a manner similar to low mass stars .
however , for high mass star formation , there are still a number of questions to be answered . for the large scale phenomena observable with a single dish telescope ,
we still do not know how large an area the outflowing gas shocks , how much energy is required to shock these regions , how large the infall region is , or how much mass from the protostellar core undergoes infall onto the central massive object .    here ,
we present follow - up observations to ( * ? ? ?
* hereafter paper 1 ) with which we address the questions posed above for seven high mass star forming regions mapped using the james clerk maxwell telescope . in paper 1
, we obtained single pointing observations towards 23 massive star forming regions previously known to have hii regions @xcite and molecular outflows as traced by cs @xcite or co @xcite .
we were able to show that , of our observed high mass star forming regions with ongoing outflow activity ( 14 sources , as traced by sio emission ) , half ( seven ) showed evidence for large scale infall motion as well ( as traced by hco@xmath0 line asymmetries ) .
we suggested that the sources with infall and outflow detections were still actively accreting material ( possibly in an ionized form , like that discussed in @xcite , while the other half may have already finished accreting ( like the source discussed in * ? ? ? * ) .
however , with single pointing observations , we could not determine the sizes of the infalling or shocked ( outflowing ) regions , and therefore could not discuss properties such as the extents of the infall and outflow regions , whether these phenomena were beam diluted , or how their sizes compared .
in this paper , we present maps of seven sources in the same tracers as paper 1 in order to understand better the nature of the large scale infall and outflow towards massive star forming regions .    for this paper
, we chose to map the sio ( j=8 - 7 ) and hco@xmath0/h@xmath1co@xmath0 ( j=4 - 3 ) emission towards seven massive star forming regions .
six sources were chose from paper 1 , and were chosen to be representative of the survey as a whole .
two sources ( g29.96 - 0.02 and g5.89 - 0.39 ) showed only outflow signatures , while two showed both infall and outflow signatures ( g10.47 + 0.03 and g19.61 - 0.23 ) , and one source ( k3 - 50a ) showed only a weak outflow signature .
it was unclear whether the emission towards the sixth source ( g20.08 - 0.14 ) was due to one cloud ( as suggested by the single peaked sio line profile ) or two ( as suggested by the double peaked hco@xmath0 and h@xmath1co@xmath0 line profiles ) .
the g20.08 region was previously known to have two ultracompact hii regions in close proximity ( i.e. * ? ? ?
* ) , although the two hii regions are separated by more than a single jcmt beam .
to this list of six sources , we added one more source , g10.6 - 0.4 .
this source was added to determine whether the ionized and molecular accretion observed on small scales in this region ( i.e. * ? ? ?
* ; * ? ? ?
* ) has a corresponding larger scale molecular infall signature .
g10.6 - 0.4 will act as a good test for determining whether our large scale results reflect the smaller scale dynamics , since the small scale structures are fairly well understood for this region . from this point forward ,
all source names based on galactic coordinates will be truncated to latitudes only .    in section [ sec : observations ] we present our observations and in section [ sec : results ] we determine outflow and infall properties based on these observations .
our results are then discussed in section [ sec : discussion ] , and we summarize in section [ sec : conclusions ]
observations of sio ( j=8 - 7 ) , hco@xmath0 and h@xmath1co@xmath0 ( j=4 - 3 ) were obtained at the jcmt towards g10.47 , g10.62 , g19.61 , g20.08 , g29.96 and k3 - 50a in april and august 2007 as project m06bc08 along with hco@xmath0 observations of g5.89 .
sio ( j=8 - 7 ) and h@xmath1co@xmath0 were observed simultaneously towards g5.89 - 0.39 as part of project m03ac13 ( see * ? ? ? * for a description of these observations ) and taken from the jcmt archive for this study .
the archive map is smaller than our new maps ; however , both the sio and h@xmath1co@xmath0 emission falls below the rms limits within the edges of this map .
the spectral resolution of the new datasets is 0.42 km s@xmath2 , the observations are spaced at 7@xmath3 to nyquist sample the 14@xmath3 full width of the jcmt beam at these frequencies , and the main beam efficiency was 0.62 .
map center positions and local standard of rest velocities are given for our six original sources in paper 1 .
for g10.6 , these properties are as follows : , , and . for reasons presented in section [ sec : g1047 ] , in this paper we adopt a distance of 5.8 kpc ( i.e. * ? ? ?
* ) for g10.37 , which is much closer than the value of 12 kpc ( i.e. * ? ? ?
* ) used in paper 1 .
the 2007 observations were obtained using harp - b and the acsis autororrelator system in either grid position switching or jiggle position switching mode . both methods provided similar noise limits in similar times
. the noise levels in each map are given in the respective figure captions and were generally consistent across each map , staying within approximately 30% of a mean value when receivers with bad noise levels were removed from the maps .
these observations were reduced using the joint astronomy center starlink software .
the 2003 observations were obtained with rxb in raster mapping mode using the das autocorrelator system and reduced using specx .
first order linear baselines were removed from each pixel within each map , and the reduced maps were exported to the miriad format for analysis .
sio and hco@xmath0 were chosen for this study because they trace outflow and infall , respectively , in star forming regions on all mass scales .
once each map was properly calibrated and exported into miriad , the peak brightness temperature and line full width at zero power ( fwzp ) were determined from the spectrum at the map pointing center .
the uncertainty in the integrated intensity for each line was determined using @xmath4t@xmath5n@xmath6 where t@xmath7 is the single channel rms uncertainty , @xmath8 is the velocity resolution of our observations , and n@xmath9 is the number of channels used to cover the full width of the line ( fwzp ) .
note that for sources with self absorbed hco@xmath0 profiles , the peak line temperature is an underestimate of the source brightness temperature .    in our maps ,
we find that the hco@xmath0 emission towards each source is extended , as is the h@xmath1co@xmath0 emission ; however , h@xmath1co@xmath0 does not reach the same spatial scales as the 77 times more abundant ( wilson & rood , 1994 ) hco@xmath0 . in only five of our seven sources
is the sio emission at least marginally resolved .
the deconvolved source size was determined using @xmath10 where the longest axis of the observed region full width at half maximum diameter ( @xmath11 ) was used to compare to the beam diameter ( @xmath12 ) . towards the two unresolved sources ( k3 - 50a and g19.61 ) , the deconvolved diameter of the emitting region is less than 50% of the jcmt primary beam , while the deconvolved diameter is greater than 70% of the jcmt primary beam towards our other sources .
our sio detection is weak in k3 - 50a ( the line peak was only detected at 3.5 - 4@xmath13 ) and was insufficient for detecting sio beyond the central pixel .
this was not the case for g19.61 ( the line peak was detected at 7@xmath13 ) , yet the sio emission in this source was also unresolved .
thus , if the shocked region in g19.61 were as extended as our jcmt beam , we would have detected the extended sio emission .
emission maps for each source are shown in figures [ fig : g5map]-[fig : k3 - 50map ] . in each case , hco@xmath0 emission spectra are shown for pixels with hco@xmath0 integrated intensities which are 10@xmath13 or greater .
for each tracer , the intensity scale for the hco@xmath0 spectra is constant for a given source , and the peak hco@xmath0 intensity for each source is listed in table [ tab : hco_prop ] .
the dashed contour in each plot represents the 3@xmath13 emission region for h@xmath1co@xmath0 and the solid contours represent sio emission , starting at 3@xmath13 and increasing in intervals of 3@xmath13 .
the 3@xmath13 rms uncertainties in the integrated intensities are given in each figure caption .
the star in each figure marks the position of the peak in the hco@xmath0 integrated intensity . in all cases except k3 - 50a and g29.96 ,
the positions of the hco@xmath0 peak is within the same beam as the sio peak .
for the two remaining sources , this offset may be due to pointing errors between the two maps .    using the fwzp from the central pixel ,
the integrated intensity for each species across each map was determined .
the rms uncertainty in the integrated intensity was determined using the same method as presented in paper 1 . for each tracer
, the line is broadest in the map center , and so using the same fwzp across the entire map captures all of the emission from that transition .
the integrated intensities averaged over each map in sio and hco@xmath0 are shown in tables [ tab : sio_prop ] and [ tab : hco_prop ] , respectively .
the integrated intensities in sio and hco@xmath0 were used to find the column density of gas emitting across the map in that tracer . for optically thin gas ,
the column density in the higher energy level is proportional to the integrated intensity of the line , with a scaling constant proportional to the frequency of the transition and the einstein coefficient for that transition ( see equation 1 of paper 1 , or * ? ? ?
assuming local thermodynamic equilibrium and an ambient temperature of 44 k ( as was done in paper 1 ) , we can use the partition function to determine the column density of optically thin gas in each tracer ( see for instance equation 14.38 of * ? ? ?
the average column densities of sio and hco@xmath0 are presented in tables [ tab : sio_prop ] and [ tab : hco_prop ] .
the values for the h@xmath1co@xmath0 column density are not presented here , but were found to be only 4 - 22 times less than those for hco@xmath0 . that we do not see a larger difference in the column densities , given the abundance ratio of 77 between the isotopologues ,
suggests that we are not seeing all of the hco@xmath0 emission because the line is optically thick .
thus , our derived column densities are lower limits to the true hco@xmath0 column densities due to optical depth effects .
this trend was also seen in paper 1 , and supports the claim that the double peaked hco@xmath0 profiles towards four of our sources are due to self absorption of the optically thick line , and not multiple clouds along the line of sight .      in an environment
which has recently been shocked , the abundance of sio can be enhanced up to a factor of 10@xmath14 ( from dark cloud abundances of approximately 10@xmath15 ) with respect to molecular hydrogen , due to silicon being liberated from dust grains and combining with oxygen to form sio ( i.e. * ? ? ?
* ; * ? ? ?
* ) . after approximately 10@xmath16 years
, sio depletes out of the gas phase again and the abundance drops ( i.e. * ? ? ?
because of the short enhancement lifetime , sio is generally only detectable in regions which have recently been shocked .
where we detect sio ( j=4 - 3 ) we suggest there is ongoing outflow activity .
we present the column density of sio towards each source in table [ tab : sio_prop ] .
we can divide this value by the abundance of sio with respect to h@xmath17 in order to determine the column density of shocked gas in each region .
these values are also presented in table [ tab : sio_prop ] .
we used the same methods as described in paper 1 to calculate the abundance of sio towards g10.6 and we have assumed a constant sio abundance across each map .
we note that the sio column densities listed in paper 1 are over estimated by a factor of 3.8 , and thus the logarithmic sio abundances listed in that paper are corrected in table [ tab : sio_prop ] .
because we have adopted a new distance for g10.47 , we have had in addition to modify our original abundance estimate from paper 1 .
the cs column density quoted in @xcite appears to be distance independent ; however , the cs abundance presented in @xcite is dependent on the assumed distance to the source ( see their equations 7 and 12 ) .
thus , our sio abundance has changed from -9.53 in paper 1 to -10.42 when the new distance and revised column density are taken into account .
the sio abundance determined for k3 - 50a in paper 1 is much closer to that found in dark clouds than in outflow regions .
the region in which sio is emitting is unresolved by our observations and highly beam diluted .
thus , for our analysis , we have assumed an sio abundance of -10.15 , the average abundance from our other six sources .
based on the abundance estimates presented in table [ tab : sio_prop ] , we have determined the total column of shocked material towards each source , from which we determined the mass of shocked gas because we know the area of the emitting region ( see below ) . for our sources which are unresolved in sio
, we have taken the radius of our primary beam to be the radius of the emitting region .    comparing the integrated intensities of our sio emission in the central pixel of each map , where our signal to noise ratio is highest , to the models presented in @xcite
, we suggest shock velocities of 20 - 22 km s@xmath2 , and ambient densities of @xmath18 @xmath19 .
the ambient densities suggested here are consistent with those of @xcite who determined densities and cs column densities for these regions using large velocity gradient models . for each of our sources ,
we adopt a shock velocity of @xmath20 km s@xmath2 .
because we have mapped the sio emission towards seven sources , we can determine the radius of the region in which sio is being emitted ( r@xmath21 ) .
the method used to determine the deconvolved source size was described in section [ sec : observations ] , and the radius of the emitting region is presented in table [ tab : sio_prop ] . from this radius and the shock velocity determined above , we can constrain the ages of our outflows from simple kinematic arguments . these approximate outflow
ages are presented in table [ tab : hco_prop ] and are used later to constrain the amount of infalling material .    in order for the outflow shocks to reach their observed sizes given the mass of shocked gas , the powering sources for these outflows must produce enough energy to push through the ambient medium .
if we assume a steady wind is driving a spherical bubble of gas outward , we can use the model of @xcite to determine the luminosity required to push the ambient gas .
equation 6 of @xcite can be re - expressed as :    @xmath22    where @xmath23 is the luminosity imparted by the wind , @xmath24 is the radius of the bubble , @xmath25 is the shock velocity , and related to the outflow age by @xmath26 , and @xmath27 is the ambient density , and related to the ambient number density by the mass of molecular hydrogen ( @xmath28 ) .
while much of this wind luminosity goes into accelerating the ambient medium , some of it is also transferred into thermal energy . from @xcite ,
the thermal energy can be expressed as @xmath29 , or , given the value for @xmath30 above , @xmath31 .
both the calculated wind luminosity and thermal energy required to power these wind driven shocks are listed in table [ tab : sio_prop ] .
these thermal energies are of approximately the same order of magnitude as the kinetic energies calculated for the co outflows in a number of massive star forming regions ( see for instance * ? ? ?
* ; * ? ? ?
* ) ; similarly the wind luminosities are comparable to the luminosities of their hii regions ( see for instance * ? ? ?
many authors suggest that an asymmetric self absorbed line profile of an optically thick line tracer in which the blue emission peak is brighter than the red , and the optically thin tracer has a single line peak , can be indicative of infall ( i.e. * ? ? ?
* ; * ? ? ?
, our optically thick line is hco@xmath0 ( j=4 - 3 ) and the optically thin isotopologue is h@xmath1co@xmath0 ( j=4 - 3 ) .    from paper 1
, we expected to observe infall signatures in g10.47 and g19.61 .
infall is detected towards these two sources , as well as towards g20.08 from our original study . in our single pointing observations ,
g20.08 was shown to have a double peaked line profile in both hco@xmath0 and h@xmath1co@xmath0 , which was suggested to be due to two overlapping sources . however , in the more extended map , we see that h@xmath1co@xmath0 becomes single peaked , suggesting that perhaps it is optically thick towards the center of the source and that mapping is required to resolve the ambiguity with respect to the chance super - position of line of sight clouds and self absorption of multiple isotopologues .
g10.6 was not in our original sample , but was added because it was known to have small scale infall motions .
it too is seen to have a large scale infall signature .    in paper 1
, we suggested that beam dilution may have been a contributing factor as to why we did not detect infall in a number of sources .
that we have detected extended infall in two distant ( @xmath32 6 kpc ) sources in this study suggests that if beam dilution of the signal is a problem , it only affects intrinsically smaller sources , and not intrinsically more distant ones , since our unresolved detections are for g19.61 and g20.08 at 4.5 and 4.1 kpc .
infall velocities were calculated for each self absorbed hco@xmath0 spectrum towards g10.47 , g10.6 , g19.61 , and g20.08 using the models of @xcite in the same way as described in paper 1 . from our rms uncertainties in line temperatures and the spectral resolution of our observations , we determined the uncertainties associated with our derived infall velocities .
we used these uncertainties to limit the regions in which we claim infall detections . only if the derived infall velocity is more than three times our uncertainty for that spectrum do we claim to detect infall in that pixel . towards g19.61 and g20.08 , robust detections were made in one and three pixels ( respectively ) suggesting that the infall in both of these sources is unresolved . towards g10.47 and g10.6 ( both at approximately 6 kpc ) , robust infall detections were made in seven nyquist sampled pixels per map , suggesting that the large scale infall in these regions is extended .    along the northern
, blue shifted , outflow lobe in g10.47 , the hco@xmath0 line profiles again become double peaked .
it is quite possible that there is a second instance of a collapsing region here since sio is marginally detected ( at the 3@xmath13 level ) in the middle of this northern extension .
the double peaked profiles do lead to infall velocities ( in three pixels ) which are greater than 3@xmath13 ; however , we do not include this region in our infall analysis for a number of reasons .
the infall signature is not continuous from the central region to this northern region , which suggests that there are two distinct regions of infall .
this discontinuity would confuse conclusions about the central infall region .
similarly , we are interested in the infall activity surrounding hii regions , where we know massive stars are forming .
neither the mid infrared nor radio maps of this region ( @xcite and @xcite , respectively ) extend far enough to the north to determine whether there are hii regions in this northern extension .
archival scuba maps of this region at the same resolution as our observations @xcite show emission which has the same extent and morphology as our hco@xmath0 observations .    for each region in which infall was observed
, we determined an average infall velocity by taking the mean of the infall velocities determined for each position in each map using the models of @xcite .
this averaging resulted in the velocities shown in table [ tab : hco_prop ] .
we chose to use this method over creating an averaged spectrum and determining infall velocities from it .
the latter method resulted in unreasonably high infall rates ( i.e. 5 km s@xmath2 in g10.47 ) , which give very large mass infall rates ( 0.22 m@xmath33 yr@xmath2 ) and total infall masses ( @xmath34 m@xmath33 ) assuming constant infall rates over the lifetime of the outflow . averaging the spectra and then determining the infall velocity increased the infall velocity by increasing the full width at half maximum ( fwhm ) of the line . since the infall velocity varies as the fwhm squared ( see equation 9 of * ? ? ? * ) , this broadening leads to an overestimate of the infall velocity .
thus , for our unresolved sources , the calculated infall velocities may be upper limits , since there could be unresolved velocity gradients within the jcmt beam .
higher resolution imaging of these sources could decrease our calculated infall velocities .    using the derived average infall velocities for each region , we determined mass infall rates for each sources using @xmath35 , where @xmath36 is the mass density
, @xmath37 is the infall velocity presented in table [ tab : hco_prop ] , and @xmath38 is the area of the infalling area ( which we assume to be circular ) . since the area in which we detected an infall signature was always similar to the area in which sio was detected , we used the radius of the sio emitting region as the radius of the infalling region .    from the radii of the sio emitting regions and a shock velocity of 21 km s@xmath2 , we were able to estimate outflow ages in section [ subsec : sio ] .
assuming infall and outflow are contemporaneous , these outflow ages can be used to constrain how much mass has undergone infall .
these values are also given in table [ tab : hco_prop ] , and represent the total mass of infalling molecular gas , not just the mass of the infalling hco@xmath0 .
it is interesting to note that the infall and outflow signatures in g19.61 are both unresolved , and the signatures in g20.08 are both marginally resolved .
this comparison suggests that the shocked region , as traced by sio , may have a similar physical size as the infall regions in these smaller sources . towards g10.6 and g10.47 , the regions over which we have robustly detected extended infall , our deconvolved sio emitting regions also appear to be on the same approximate spatial scales .
comparing infall mass and outflow mass , we note that our infall masses are lower than our outflow masses .
this is contrary to the results of @xcite ; however , their ratio was calculated for the protostellar jet , not the entrained molecular outflow .
this comparison to @xcite suggests that more of the ambient mass is entrained into the large scale molecular outflow ( produced by a protostellar jet ) than falls in onto the forming stars ( i.e. * ? ? ?
* ) .      the large spatial extent of the hco@xmath0 emission in our map of g10.47 places it among the largest molecular outflows due to star formation in our galaxy .
some herbig - haro objects have been observed emanating from nearby lower mass star forming regions at distances of up to 12 pc from their powering stars ( i.e. * ? ? ?
* ) ; however in those cases , the associated molecular outflow is not nearly as wide as the case of g10.47 , and those jets are not pushing through an ambient medium dense enough for massive stars to form ( i.e. 10@xmath39 @xmath19 ) . in paper 1 , we used a distance estimate of 12 kpc to g10.47 .
given the extent of hco@xmath0 gas ( @xmath40 ) in our current maps , this distance would give a linear scale of over 4 pc for a single outflow .
other authors have suggested a distance to this source of 5.8 kpc ( i.e. * ? ? ?
for the analysis presented here , we adopt this new distance , which changes the length of the outflow from 4 pc to 2 pc .
given that this source is along a line of sight reasonably close to both the galactic plane and galactic center , it is not unreasonable to suggest that what appears to be an extremely large outflow is instead the super - position of a number of clouds along the line of sight .
however , based on the arguments presented below , we suggest that the gas kinematics and morphology are due to a large outflow .    examining the spectra of our hco@xmath0 map , we note that while the hco@xmath0 emission is extended , it is not present in all pixels , suggesting that we have captured all of the emission in this velocity range .
figure [ fig : averaged_spectra ] shows the hco@xmath0 line profiles averaged over the blue and red shifted emission regions as well as the hco@xmath0 and h@xmath1co@xmath0 line profiles averaged over the central nine pixels of our map .
the red and blue spectra , while having very little overlap between their profiles , do both fit within the width of the source center line profile .
the velocities at line peak for the blue and red spectra are -64 and -70 km s@xmath2 respectively , and match the velocities of the two peaks of the self absorbed , averaged hco@xmath0 central profile .
the single peak of the optically thin h@xmath1co@xmath0 is at -67 km s@xmath2 as is the minimum of the hco@xmath0 self absorption feature .
this self - absorption lies directly between the peaks of the red and blue profiles , and is the same as the source rest velocity .
no other sources in the w31 region have rest velocities similar to the ones seen in the g10.47 region .
the double peak in the hco@xmath0 spectrum is due to self absorption , not multiple components , since there is a single peak in the h@xmath1co@xmath0 spectrum .
the same is true for the co isotopologues in this region .
the bima observations of @xcite show one peak in c@xmath41o , while the @xmath1co observations of @xcite with the iram interferometer show two peaks indicative of self absorption .
the first moment map of the hco@xmath0 emission ( grayscale in figure [ fig : g1047_mom1 ] ) in this region also suggests a single outflow .
the observed smooth transition from blue to red along the source would be more likely to be observed in one source with a velocity gradient than three separate sources with discontinuous cloud boundaries . the second moment map of this region ( contours in figure [ fig : g1047_mom1 ] ) suggest that the largest velocity dispersion is towards the center , and continuously decreases towards the emission edges ( outflow lobe ends ) , suggesting that most of the mass is centrally concentrated .
while we favor the outflow interpretation of this morphology , higher resolution , large scale observations of this source are necessary in order to resolve this ambiguity .
this region contains four individual uchii regions and what is likely a large scale outflow with a collimation factor ( @xmath42 ) of approximately 1.5 , which may be even more highly collimated at higher resolution .
objects like g10.47 , with large scale , fairly collimated outflow emanating from a cluster of hii regions pose interesting challenges for accretion models .
as stated previously , g10.6 was added to our source list because it was previously known to show evidence for small scale infall , both in molecular gas @xcite and ionized gas @xcite .
these small scale infall signatures appear to be consistent with a single accretion flow which changes from being molecular to ionized as it passes the ionization boundary of the hii region .
while we can not determine the orientation of the large scale infall from our observations , we note that our derived mass infall rate ( 0.03 m@xmath43 yr@xmath2 ) is comparable to , but slightly larger than , that of @xcite ( 0.02 @xmath44 yr@xmath2 ) .
this similarity may suggest that the large scale infall possibly maps to the small scale accretion flow ; in addition , despite resolving the infall region , we may still be over estimating the infall velocity due to velocity gradients in our larger jcmt beam ( see section [ sec : hco ] ) .
none of our basic conclusions drawn from our single pointing observations of sio ( i.e. column densities in the central pixel or outflow lifetimes ) have changed due to mapping save for the infall detection in g20.08 .
however , with respect to the outflowing gas , we have been able to place constraints on the energy required to shock these regions . with respect to infalling gas ,
we have placed better constraints on infall velocities and shown that the lack of infall signatures in some sources is not due to the signature being beam diluted out of observability at large distances .
we have found that one of our sources which was previously not thought to have an infall signature ( g20.08 ) in fact does appear to have large scale infall motions .
this discovery suggests that more accurate infall results are possible with maps of massive star forming environments than with single pointing observations .
returning to the results of paper 1 , we determined that seven out of the 14 sources with ongoing outflow signatures had infall signatures as well .
these numbers can now be revised to eight out of 15 outflow sources with infall signatures when g10.6 is added to the source list .
thus , we find that approximately half of the sources with ongoing outflows have infall signatures which we take to trace accretion .
because our sources were selected based on the presence of hii regions , ongoing infall suggests that there is a phase of massive star formation in which the star is still actively accreting , but is also burning hydrogen brightly enough to ionize its surroundings ( see for instance the discussion in * ? ? ?
we have found that just over half of the sources in our sample have large scale infall motions suggesting that the protostars powering the uchii regions and massive outflows are still actively accreting material .
where we have detected infall motions , they have been on quite large scales , suggesting that infall , as traced by hco@xmath0 , is in general observable with single dish maps at the distances to massive star forming regions .
maps are required in order to capture all of the infalling gas since the infall signature can be extended , even at 6 kpc . in mapping these regions ,
we have been able to determine how much mass is falling into the central regions of our sample sources , which we would not have been able to determine with our single pointing observations .
we find that when the infall is resolved , the infalling mass is greater than 30 m@xmath43 .
our derived mass infall rate towards g10.6 is similar to the small scale rate found by @xcite , suggesting that we may be observing the same physical process , but on different spatial scales .    in general
, we have found that sio is not very extended ; however , since we mapped these regions in sio , we were able to determine the radius of the emitting region , the mass of shocked gas and how much energy is required to shock the gas surrounding a massive protostar .
we find these energies and luminosities to be similar to the kinetic energy imparted to the outflowing gas and the luminosities of the hii regions themselves .
we also find that the region over which we detect an infall signature and the region in which we detect sio is quite similar , and find that in these regions , the total infall mass is less than the total outflow mass ( m@xmath45/m@xmath46 0.1 to 0.8 ) assuming the two processes are contemporaneous .
we would like to acknowledge the support of the national science and engineering research council of canada ( nserc ) .
p. d. k would also like to thank j. wadsley for helpful discussions during the preparation of this manuscript .
rccccccccccc g5.89 - 0.39 & 2.0 & 2.37 & 105 & -9.48 & 4.38 & 12.81 & 22.29 & 3.7 & 1.9 & 1.1 & 0.18 + g10.47 + 0.03 & 5.8 & 1.24 & 30 & -10.42 & 2.70 & 12.60 & 23.02 & 184 & 5.8 & 10.3 & 5.1 + g10.6 - 0.4 & 6.0 & 1.14 & 20 & -9.92 & 8.74 & 13.11 & 23.03 & 497 & 9.4 & 27.2 & 21.6 + g19.61 - 0.23 & 4.5 & 0.71 & 45 & -10.32 & 4.57 & 12.83 & 23.15 & 162 & @xmath474.7&@xmath476.8&@xmath472.7 + g20.08 - 0.14 & 4.1 & 0.66 & 30 & -11.06 & 2.68 & 12.60 & 23.66 & 374 & 4.0 & 4.8 & 1.6 + g29.96 - 0.02 & 9.0 & 0.63 & 40 & -9.70 & 5.14 & 12.88 & 22.58 & 124 & 7.9 & 19.2 & 12.8 + k3 - 50a & 8.6 & 0.48 & 20 & -10.15 & 3.15 & 12.67 & 22.82 & 276 & @xmath479.0&@xmath4724.7&@xmath4718.7 [ tab : sio_prop ]    rccccccccc g5.89 - 0.39 & 31.5 & 140 & 35.6 & 13.17 & & & 0.29 & & + g10.47 + 0.03 & 11.7 & 60 & 11.2 & 12.67 & 1.02@xmath480.1 & 8.89 & 0.88 & 78.0 & 0.42 + g10.6 - 0.4 & 31.2 & 50 & 28.4 & 13.08 & 1.19@xmath480.1 & 27.34 & 1.42 & 389 & 0.78 + g19.61 - 0.23 & 8.3 & 80 & 14.4 & 12.78 & 1.49@xmath480.6 & @xmath478.5&@xmath470.71&@xmath4760.3&@xmath470.37 + g20.08 - 0.14 & 5.9 & 55 & 9.2 & 12.59 & 1.35@xmath480.3 & 5.51 & 0.60 & 33.0 & 0.09 + g29.967 - 0.02 & 16.5 & 30 & 9.2 & 12.59&&&1.19 & & + k3 - 50a & 17.0 & 30 & 13.2 & 12.75&&&@xmath471.36 & & + [ tab : hco_prop ]              beuther , h. , & shepherd , d.  2005 , cores to clusters : star formation with next generation telescopes , 105 beuther , h. , schilke , p. , sridharan , t.  k. , menten , k.  m. , walmsley , c.  m. , & wyrowski , f.  2002 , a&a , 383 , 892 castor , j. , mccray , r. , & weaver , r.  1975 , apjl , 200 , l107 caselli , p. , hartquist , t.  w. , & havnes , o.  1997 , a&a , 322 , 296 di francesco , j. , johnstone , d. , kirk , h. , mackenzie , t. , & ledwosinska , e.  2008 , , 175 , 277 fuller , g.  a. , williams , s.  j. , & sridharan , t.  k.  2005 , a&a , 442 , 949 gibb , a.  g. , wyrowski , f. , & mundy , l.  g.  2004 , apj , 616 , 301 gregersen , e.  m. , evans , n.  j. , ii , zhou , s. , & choi , m.  1997 , apj , 484 , 256 hunter , t.  r.  1997 , ph.d .
thesis , keto , e. , & klaassen , p.  2008
, , 678 , l109 keto , e.  2007 , apj , 666 , 976        klaassen , p.
d. , & wilson , c.  d.  2007a , , 663 , 1092 klaassen , p.  d. , & wilson , c.  d.  2007b , arxiv e - prints , 710 , arxiv:0710.3785 klaassen , p.  d. , plume , r. , ouyed , r. , von benda - beckmann , a.  m. , & di francesco , j.  2006 , apj , 648 , 1079 kroupa , p.  2001 , mnras , 322 , 231 kurtz , s. , cesaroni , r. , churchwell , e. , hofner , p. , & walmsley , c.  m.  2000 , protostars and planets iv , 299 kurtz , s. , churchwell , e. , & wood , d.  o.  s.  1994 , apjs , 91 , 659 mardones , d. , myers , p.  c. , tafalla , m. , wilner , d.  j. , bachiller , r. , & garay , g.  1997 , apj , 489 , 719 mckee , c.  f. , & tan , j.  c.  2003 , , 585 , 850 myers , p.  c. , mardones , d. , tafalla , m. , williams , j.  p. , & wilner , d.  j.  1996 , apjl , 465 , l133 olmi , l. , cesaroni , r. , & walmsley , c.  m.  1996 , a&a , 307 , 599 pascucci , i. , apai , d. , henning , t. , stecklum , b. , & brandl , b.  2004 , a&a , 426 , 523 pillai , t. , wyrowski , f. , carey , s.  j. , & menten , k.  m.  2006 , , 450 , 569 pineau des forets , g. , flower , d.  r. , & chieze , j .-
1997 , herbig - haro flows and the birth of stars , 182 , 199 plume , r. , jaffe , d.  t. , evans , n.  j. , ii , martin - pintado , j. , & gomez - gonzalez , j.  1997 , apj , 476 , 730 riddick , f.  c. , roche , p.  f. , & lucas , p.  w.  2007 , , 381 , 1077 rholfs , k. , wilson , t.  l. 2004 , tools of radio astronomy ( 4th ed ; berlin ; springer )    schilke , p. , walmsley , c.  m. , pineau des forets , g. , & flower , d.  r.  1997 , a&a , 321 , 293 sewilo , m. , churchwell , e. , kurtz , s. , goss , w.  m. , & hofner , p.  2004 ,
apj , 605 , 285 shepherd , d.  s. , & churchwell , e.  1996 , , 472 , 225 shirley , y.  l. , evans , n.  j. , ii , young , k.  e. , knez , c. , & jaffe , d.  t.  2003 , apjs , 149 , 375 shu , f.  h. , adams , f.  c. , & lizano , s.  1987 , , 25 , 23 sollins , p.  k. , zhang , q. , keto , e. , & ho , p.  t.  p.  2005
, , 624 , l49 tielens , a. g. g. m. 2005 , the physics and chemistry of the interstellar medium ( cambridge university press ) wood , d.  o.  s. , & churchwell , e.  1989 , apjs , 69 , 831 wu , y. , wei , y. , zhao , m. , shi , y. , yu , w. , qin , s. , & huang , m.  2004 , a&a , 426 , 503 zinnecker , h. , & yorke , h.  w.  2007 , ara&a , 45 , 481